The TLR9 agonist MGN1703 triggers a potent type I interferon response in the sigmoid colon.

Toll-like receptor 9 (TLR9) agonists are being developed for treatment of colorectal and other cancers, yet the impact of these drugs on human intestines remains unknown. This, together with the fact that there are additional potential indications for TLR9 agonist therapy (e.g., autoimmune and infectious diseases), led us to investigate the impact of MGN1703 (Lefitolimod) on intestinal homeostasis and viral persistence in HIV-positive individuals. Colonic sigmoid biopsies were collected (baseline and week four) from 11 HIV+ individuals on suppressive antiretroviral therapy, who received MGN1703 (60 mg s.c.) twice weekly for 4 weeks in a single-arm, phase 1b/2a study. Within sigmoid mucosa, global transcriptomic analyses revealed 248 modulated genes (false discovery rate<0.05) including many type I interferon (IFN)-stimulated genes. MGN1703 increased the frequencies of cells exhibiting MX1 (P=0.001) and ISG15 (P=0.014) protein expression. No changes were observed in neutrophil infiltration (myeloperoxidase; P=0.97). No systematic effect on fecal microbiota structure was observed (analysis of similarity Global R=-0.105; P=0.929). TLR9 expression at baseline was inversely proportional to the change in integrated HIV DNA during MGN1703 treatment (P=0.020). In conclusion, MGN1703 induced a potent type I IFN response, without a concomitant general inflammatory response, in the intestines.

Introduction

The intracellular innate immune surveillance molecule Toll-like receptor 9 (TLR9) is a pattern recognition receptor found in plasmacytoid dendritic cells (pDC), B cells and in some mucosal epithelial cells [1-4]. DNA containing non-methylated cytosine-guanine sequence motifs (CG- or CpG-motifs) bind to and activate TLR9 in either of two types of endosomes, depending on the structural context and CG-motif sequence environment. In one type of endosome, TLR9 signal transduction is via activation of nuclear factor kappa B (NFκB) which preferentially leads to transcription of proinflammatory cytokine genes. In the other endosome type, signal transduction is through interferon regulatory factor 7 (IRF7) which preferentially drives transcription of type I interferon genes [5, 6]. Paracrine and autocrine cytokine signaling then leads to activation of innate and adaptive immune responses [7, 8]. To clinically utilize TLR9’s vital immune surveillance activity, families of synthetic DNA molecules were developed as TLR9 ligands (agonists and antagonists). TLR9 agonists have been and are being tested as vaccine adjuvants and as therapeutics for cancers – including colorectal carcinoma, infectious diseases, inflammatory bowel diseases (IBD), asthmas and allergies [e.g. ClinicalTrials.gov ID: NCT02077868; NCT02200081; NCT02254772; and NCT02668770] [9-14]. However, the effects of different TLR9 agonists within the intestines have only been described in mice, sheep and fish [15-17] while the physiological effects of any TLR9 agonist in human intestines were essentially unknown prior to this study.

In the context of IBD, TLR9 agonists could play a role in modulating aberrant immune responses that lead to inflammation, for example by inducing tolerance towards commensal bacteria [9-12]. In the context of infectious diseases, TLR9 agonists may enhance immune surveillance activities that facilitate reductions in pathogen burden [18]. Specifically in the case of HIV eradication, therapeutic use of TLR9 agonists will likely be undertaken in a combination strategy within HIV-positive individuals on suppressive antiretroviral therapy (ART). Therefore, we enrolled such persons and performed a comprehensive assessment of the impact on intestinal homeostasis and inflammation by the TLR9 agonist and immune surveillance reactivator MGN1703 (lefitolimod).

MGN1703 is synthetic TLR9 agonist of the double-stem loop immunomodulators or dSLIM family where multiple CG-motifs are embedded into a stable, covalently-closed double stem-loop DNA molecule harboring only natural phosphodiester bonds between nucleotides [19-21]. MGN1703’s pre-clinically-determined efficacy parameters are maintained in clinical studies without dose-limiting toxicities [19-21] and this drug is currently in phase III testing for treating metastatic colorectal cancer [21]. The clinical proof-of-concept for MGN1703 was achieved in a clinical phase II study evaluating switch-maintenance treatment following standard chemotherapy of patients with metastatic colorectal carcinoma [21, 22].

Currently, there is a renewed interest in developing effective and accessible therapies that can lead to eradication of HIV infection or that facilitate sustained viral remission in the absence of ART [23]. The main barrier to HIV eradication is the ability of the virus to establish latent infection in long-lived CD4+ T cells that persist in the body despite ART [24]. These latently infected cells are the likely source of plasma viral rebound after interruption of ART and there is compelling evidence that the intestines are a major reservoir for HIV persistence [25, 26]. To eradicate persistent HIV, robust immune surveillance is required [27]. We previously showed (ex vivo and in vivo) that the cells of HIV-positive individuals responded appropriately to TLR9 agonist stimulation [18, 28, 29]. Thus, TLR9 agonist-mediated enhancement of the activation status and killing capacity of CD8+ T and NK cells may augment the reduction or elimination of persistent viruses.

In our preclinical evaluation of MGN1703 ex vivo using peripheral blood mononuclear cells isolated from HIV infected individuals, we found that MGN1703 incubation results in NK cell activation and increased NK cell inhibition of HIV spreading in a culture of autologous CD4+ T cells [18]. As we examined the systemic effects evaluating MGN1703 effects in our clinical trial, we found that treatment in HIV infected individuals has a dual potential in HIV eradication efforts as MGN1703 increased HIV-1 transcription and enhanced cytotoxic NK cell activation in vivo [30]. Here we report a sub-study of this trial where we define MGN1703-treatment outcomes related to: (i) TLR9 agonist-induced intestinal immune activity; (ii) maintenance of intestinal homeostasis during TLR9 agonist treatment; (iii) stability of the gut microbiota composition during TLR9 agonist treatment; and (iv) TLR9 agonist-associated changes in HIV persistence in the intestine.

Results

MGN1703 induced a potent type I interferon response in the intestines

The clinical trial design is depicted in (Figure 1a). Participants received subcutaneous injections of MGN1703 (60 mg) two times per week for four weeks while maintaining ART. We previously reported that MGN1703 dosing led to activation of peripheral blood pDC, cytotoxic NK cells and CD8+ T cells as well as substantial increases in plasma interferon-α2 levels in this cohort [30]. For the current study, endoscopic collection of colonic sigmoid biopsies were performed at baseline and during the last week of MGN1703 administration. We used quantitative IHC on intestinal biopsy tissues to detect the canonical type I interferon-induced proteins MX1 and ISG15 [31]. We observed an increase in the numbers of MX1 positive cell profiles from baseline to the fourth week of treatment (Figure 1b) in both the intestinal lamina propria (p=0.001) (Figure 1c) and epithelium (p=0.003) (Figure 1d). The number of ISG15 positive cell profiles in lamina propria also increased during MGN1703 treatment (p=0.014) (Figure 1e) while no significant changes in the number of ISG15-producing cell profiles in the epithelial region were observed (Figure 1f). These protein-detection observations were consistent with the global transcriptomic analyses performed using RNASeq on intestinal tissue mononuclear cells isolated from these individuals. We observed hierarchical clustering of samples in a treatment-dependent manner (Slatkin-Maddison p=0.0167; Figure 2a), and segregation of the response was also observed with principal component analyses (Figure 2b and Supplemental Figure 1). These findings indicate that there was a potent, dominant effect of MGN1703 administration on the intestinal transcriptome that supersedes inter-individual variability in gene expression. There were 248 genes whose expression was significantly modulated by MGN1703 treatment (fold change >2; FDR<0.05) (Figure 2c and Table S1). In addition to upregulating MX1 and ISG15 expression, MGN1703 treatment led to significant increases in the expression levels of multiple additional interferon-stimulated genes (e.g. IFI6, OAS and IFIT1), antiviral genes specifically associated with the control of HIV-1 in vivo (IFITM1, MX2, TRIM22), as well as interferon pathway signaling molecules (e.g. STAT1 and IRF7) (Figures 2c and 3a).

Figure 1

MGN1703 induced an increase in MX1 and ISG15 positive cell profiles in the sigmoid colon

(a) The timing of MGN1703 dosing and sample collection are illustrated. (b) Representative images for MX1 IHC (brown staining) from biopsies collected prior to dosing (baseline) and during the fourth dosing week (treatment) are presented. Scale bars=100 µm. (c–f) Plots show the number of cell profiles expressing MX1 (c–d) or ISG15 (e–f) at baseline and during MGN1703 treatment (n=11 in all graphs). Gray circles indicate results for each individual. Paired data are connected with gray lines. Black squares represent the mean value for all individuals at the indicated time point. Black lines connect longitudinal mean data points. (c and e) The numbers of positive cell profiles per mm2 within the lamina propria regions are illustrated. (d and f) The numbers of positive cell profiles per mm2 within the epithelial regions are illustrated. Statistics: Wilcoxon matched-pairs signed rank test.

Figure 2

MGN1703 significantly modulated the expression of 248 genes in intestinal mononuclear cells, including many genes involved in type I interferon responses

(a) Ward link correlation dendrogram for all genes assayed reveals that gene expression patterns cluster with respect to treatment status, rather than by individual (Slatkin-Maddison p=0.0167). (b) Principal component analysis scores plot for PC1 vs. PC4; the two principal components which contributed the maximum variance in the analysis. Plots depicting PC2 and PC3 are presented in Supplemental Figure 1. Samples collected at baseline are represented by red circles and samples collected during treatment are represented by blue circles. (c) RNASeq heat map was generated using standardized Z-scores of the genes that were significantly regulated (up or down) in intestinal mononuclear cells after adjusting for a FDR <0.05. Individual participants are depicted in the columns and one gene is depicted per row. The dendrogram reveals the relatedness of individual gene expression patterns. A full list of significantly regulated genes can be found in Table S3. (n=10 in all graphs)

Figure 3

MGN1703 did not induce inflammation in the sigmoid colon

(a) Key type I interferon response and inflammatory parameters from the RNASeq were assessed in an equivalence evaluation together with MPO data from the IHC analyses. In this forest plot representation, an upregulation of 2-fold from baseline was set as the threshold for concluding that upregulation had been induced. Confidence intervals in the forest plot are at a 10% significance level, such that the comparison to the upper equivalence limit is performed at an alpha of 0.05 (n=10 for gene expression parameters and 11 for the MPO in LP protein parameter). The median of all 16 type I “median ratio” entries is 5.5 +/− 3.0 SD. The median of all 11 type II “median ratio” entries is 0.9 +/− 0.3 SD. (b–e) Plots show the number of cell profiles expressing CXCL10 (b–c), MPO (d), or IL-21 (e) protein assessed by IHC analyses at baseline and during MGN1703 treatment (n=11 in all graphs). Gray circles indicate results for each individual. Paired data are connected with gray lines. Black squares represent the mean value for all individuals at the indicated time point. Black lines connect longitudinal mean data points. (b, d and e) The numbers of positive cell profiles per mm2 within the lamina propria regions are illustrated. (c) The numbers of positive cell profiles per mm2 within the epithelial regions are illustrated. Statistics: Wilcoxon matched-pairs signed rank test.

MGN1703 treatment was not associated with an induction of type II interferon responses or markers of inflammation in the intestines

Secondary effects of IFN-α production can include the activation of a range of effector cells [e.g. IFN-α can activate NK cells to secrete type II interferon (i.e. IFN-γ) [32]]. In our study, significant modulation of IFN-γ and type II interferon response genes [e.g. RELA (NF-κBp65), IL6, CXCL8 (IL-8) and TNF (TNF-α)] was not detected in RNASeq analyses of intestinal mononuclear cells (Figure 3a). Consistent with no change in the levels of IFN-γ mRNA in the intestinal tissue mononuclear cells, we did not detect significant changes in the numbers of cell profiles positive for the IFN-γ-induced protein CXCL10 using IHC in either the intestinal lamina propria (Figure 3b) or epithelium (Figure 3c). Nor was MGN1703 treatment associated with changes in intestinal T cell activation (Supplemental Figure 2a), T cell memory subset proportions (Supplemental Figure 2b) or expression of genes known to be involved in T cell activation (Supplemental Figure 2c). Furthermore, there was no change in the number of cell profiles in the intestinal lamina propria positive for myeloperoxidase (MPO) (Figure 3a and 3d), an enzyme primarily produced by neutrophils [33] and used as an indicator of general inflammation in the intestinal mucosa [34]. There was, however, an MGN1703-associated increase in the number of cell profiles in the lamina propria that expressed IL-21 (p=0.019, Figure 3e). Together these data indicate that MGN1703 treatment did not induce a type II interferon response or inflammation in the intestines.

MGN1703 treatment was associated with subtle changes in the intestinal microbiota

Clinical interventions (e.g., antibiotic treatment for ulcerative colitis) may decrease intestinal inflammation, but often at the expense of reduced microbial diversity in the intestines [35-37]. Here, we performed a comprehensive analysis of the effects of MGN1703 on the gastrointestinal (GI) microbiota to determine whether MGN1703 therapy altered microbial diversity. Ribosomal RNA gene amplicon libraries from fecal samples were sequenced and all datasets were rarefied to 14,500 reads for calculations of both beta and alpha diversity. No effect of treatment was observed across all fecal samples using analysis of similarity (ANOSIM; R=-0.105; p=0.929; 999 permutations), while participant ID was strongly correlated with observed community structure (R=0.935; p=0.001; 999 permutations). This observation was supported by hierarchical clustering of baseline and treatment fecal microbiota structure (Figure 4a). Bray-Curtis similarity values of paired samples (baseline and treatment) ranged from 72.6 to 91.3. Across all pairs of samples no individual taxon, at any taxonomic level, was significantly different between baseline and treatment samples (Figure 4b; Table S2). However, small, but statistically significant, increases between baseline and treatment fecal samples were detected in several alpha diversity indices, including evenness (p<0.03, taxonomic levels of class and order), Shannon index (p<0.04, taxonomic levels of class, order and family), and Simpson’s index (p<0.05, taxonomic level of family) (Figure 4c–d, Table S3). The number of bacterial species identified in the fecal samples trended upwards, but did not change significantly with MGN1703 treatment (Figure 4e). Finally, the ratio of Firmicutes-to-Bacteroidetes sequences was slightly, but not significantly, lowered in treatment samples relative to baseline (Figure 4f). These findings indicate that MGN1703 impacts alpha diversity, but had no effect on overall microbial community structure.

Figure 4

MGN1703 treatment did not radically alter the intestinal flora

(a) Group average dendrogram illustrates relationships between samples at the taxonomic level of family. Individual participants clustered at baseline and treatment with regard to diversity in their microbiota compositions. (b) The rarefied number of sequences (14,500) of microbial taxa (family) at baseline and post-treatment is shown. (c–f) Plots show the Shannon diversity (family) index values (c), the Simpson diversity (family) index values (d), the number of bacterial species identified (e), and the fermicutes over bacterioides ratios (f) at baseline and during MGN1703 treatment (n=8 in all graphs). Gray circles indicate results for each individual. Paired data are connected with gray lines. Black squares represent the mean value for all individuals at the indicated time point. Black lines connect longitudinal mean data points. Statistics for (c–f): paired t tests.

In the intestines, baseline TLR9 and IFNAR1 gene expression is inversely proportional to changes in copies of integrated HIV-1 DNA during MGN1703 treatment

The number of integrated HIV-1 DNA copies reflects the magnitude of the persistent HIV-1 reservoir [38]. As in peripheral blood CD4+ T cells [30], we did not detect significant changes in the copies of HIV-1 DNA in the intestinal mononuclear cells associated with MGN1703 dosing in this cohort (p=0.23; Wilcoxon matched-pairs signed rank test). However, we did observe that some individuals exhibited up to 10-fold reduction in integrated HIV-1 DNA copies with MGN1703 dosing. To begin elucidating host factors that contribute to the robust decrease in intestinal viral DNA observed in some participants, we analyzed MGN1703-associated changes in copies of integrated HIV-1 DNA together with expression patterns of TLR9, IFNAR1, and IL17RE. The TLR9 gene codes for the target of MGN1703. The IFNAR1 (IFN-α receptor alpha chain) gene codes for a key molecule in the IFN-α receptor. The IL17RE (IL-17 receptor E) gene encodes the IL17C receptor and IL17C is important in innate epithelial immune responses [39-41]. We observed that baseline gene expression levels of TLR9 are inversely proportional to the fold-change in integrated HIV-1 DNA copies in the intestines (p=0.020; Figure 5a). Furthermore, when the data were stratified according to having a greater than or less than 50% reduction in integrated HIV-1 DNA copies (stratification point chosen because the data naturally diverged with 5 participants above and 5 participants below this value), we observed a significant difference in the number of TLR9 transcripts present at baseline between the two groups (p=0.008; Figure 5b). We also observed a trend for higher baseline levels of IFNAR1 transcripts to associate with greater reductions in integrated HIV-1 DNA (p=0.087; Figure 5c) and we detected a significant difference in the baseline IFNAR1 expression when the reduction in integrated HIV-1 DNA copies was greater than 50% versus less than 50% (p=0.024; Figure 5d). Next, we found that the MGN1703-induced change in IL17RE gene expression is inversely proportional to the fold-change in integrated HIV-1 DNA copies in the intestines (p=0.035; Figure 5e). Finally, we observed a significant difference in the MGN1703-induced reduction in IL17RE transcription when the fold-change in integrated HIV-1 DNA copies was greater than 50% versus less than 50% (p=0.032; Figure 5f). Together, these intestinal data highlight that higher TLR9 and IFNAR1 gene expression prior to MGN1703 treatment is associated with a greater reduction in integrated HIV-1 DNA levels following treatment.

Figure 5

MGN1703-induced reductions in intestinal mononuclear cell HIV-1 DNA copies is associated with higher baseline TLR9 and IFNAR1 expression

(a, c and e) Bi-variant plots show the fold change in the copies of integrated HIV-1 DNA (left axis) versus the baseline number of TLR9 transcripts (a), the baseline number of IFNAR1 (IFN-α receptor alpha chain) transcripts (c) or the MGN1703-induced fold-change in transcription of IL17RE (IL-17 receptor E) (e). Spearman correlation r and p values are shown. (b, d and f) Data for the baseline number of TLR9 transcripts (b), the baseline number of IFNAR1 transcripts (d) or the MGN1703-induced fold-change in transcription of the IL17RE gene (f) were stratified based on whether the reduction in integrated HIV-1 DNA copies was greater or less than 50%. The HIV-1 DNA values were generated from isolated mononuclear cells. Statistics: Wilcoxon matched-pairs signed rank test.

Discussion

Here we present the first examination of the in vivo effects of TLR9 agonist therapy in the human colonic sigmoid mucosa and feces. The TLR9 agonist used herein is MGN1703 which is well-tolerated in patients and is currently in phase III testing for the treatment of metastatic colorectal cancer (ClinicalTrials.gov ID: NCT02077868) [21]. Moreover, this drug is being investigated clinically as treatment for small cell lung cancer and melanoma (ClinicalTrials.gov ID number: NCT02200081 and NCT02668770, respectively). There is also interest in using TLR9 agonists to treat IBD like ulcerative colitis [9-12] and infectious diseases, including HIV, that cause significant gut pathologies [18]. Given these clinical possibilities, we performed this study investigating the impact of MGN1703 on intestinal homeostasis and HIV persistence in HIV-positive individuals on ART. While untreated HIV infection is associated with disrupted barrier function and intestinal pathologies [42-44], the immune function restoration that occurs during ART [45, 46] supports the generalization of our findings to other indications including malignancy and IBD.

MGN1703 activates pDC to secrete IFN-α [19] which plays a critical role in the immune response to viral infections by inducing the production of proteins with described antiviral activity (e.g. MX1 and ISG15) [47, 48]. Accordingly, we observed significant increases in cell profiles positive for MX1 and ISG15 in the colon. These protein-based observations were confirmed by gene expression data demonstrating that the MX1 and ISG15 genes were significantly upregulated during MGN1703 therapy. Further evidence of a strong type I interferon response in the intestine during this treatment was the upregulated expression of many additional interferon response genes including: IFIT3, IFIT5, IFI1, IFI44, OAS1, OAS2, and OAS3. In summary, protein and gene expression data confirm that MGN1703 treatment induced a robust type I interferon response in the sigmoid colon of these individuals that may contribute to clearance of the viral reservoir and durable virologic suppression following cessation of ART.

Given the potency of the induced type I interferon response and the known inflammatory properties of TLR agonists, it is conceivable that a type II interferon response and inflammatory cascade could also be triggered by MGN1703 [13, 14]. Such an outcome could have a profound negative impact on intestinal homeostasis. Notably, we did not detect any indications of a type II interferon or inflammatory response to MGN1703 – either macroscopically during biopsy collection or via the microscopic and genetic analyses employed. The type I interferon cascade can activate monocytes to secrete CXCL10 and induce NK cells to secrete type II interferon (IFN-γ); a potent inducer of CXCL10 production [49]. While an increase in CXCL10 protein was observed in plasma during MGN1703 treatment [30], we did not find any cohort wide changes in the number of CXCL10 positive cell profiles in the sigmoid colon as would be expected during an inflammatory response [50, 51]. We also quantitated the number of infiltrating neutrophils at baseline versus on treatment. In our equivalence evaluation, we found no evidence of inflammatory neutrophil infiltration even though we set a low threshold (2-fold increase over baseline) for concluding that inflammation occurred (5-to 7-fold increases in infiltrating neutrophils have been reported in animal models of intestinal inflammation [52, 53]). Moreover, the gene expression levels of major inflammatory mediators [e.g. IFNG, RELA, IL6, CXCL8 and TNF] remained within the same 2-fold bound. Given that IL-21 combined with ART is known to reduce inflammation in vivo [54], the increased number of cell profiles producing IL-21 protein could, at least in part, account for the absence of neutrophil infiltration and inflammation in the intestines during MGN1703 treatment. Together, these observations lead us to conclude that MGN1703 did not induce an inflammatory response that could disrupt intestinal homeostasis.

The stability of the microbiota is another major factor that contributes to the maintenance of intestinal homeostasis [55-57]. While HIV-positive individuals have lower microbiota diversity [58], there were no major alterations in microbiota diversity associated with MGN1703 treatment. In fact, only subtle increases in multiple diversity parameters were observed. We also observed a nonsignificant trend in this cohort towards a lower firmicutes to bacteroidetes ratio following MGN1703 dosing. Given these microbiota results and our data revealing a lack of MGN1703-associated intestinal inflammation, we conclude that MGN1703 dosing minimally impacts intestinal homeostasis and might lead to modest improvements in the overall GI microflora composition.

In the context of a successful HIV eradication strategy, MGN1703 will likely be combined with other interventions. However, it is essential that we have a comprehensive understanding of the individual impact of this drug on HIV-1 reservoirs to fully interpret outcomes in future trials testing combinations of HIV eradication interventions (e.g. to detect synergy). Here we observed that increased transcription of either TLR9 or IFNAR1 at baseline was associated with a reduction in copies of integrated HIV-1 DNA in intestinal mononuclear cells. Thus, a greater number of targets for MGN1703 (i.e. TLR9) and for IFN-α (i.e. IFN-α receptor) at the time of dosing is associated with improved efficacy in eliminating HIV-1 DNA-containing cells in the intestines during the course of treatment. Further downstream in the IFN-α signaling cascade is the action of IL-17C through IL-17 receptor E which leads to antimicrobial peptide production to help maintain intestinal epithelial barrier integrity [39]. When examining changes in IL17RE gene expression levels following MGN1703 dosing, we found that this parameter was inversely proportional to the reduction in integrated HIV-1 DNA copies in the intestines. These data indicate that each participant’s expression levels for key genes directly related to TLR9 function (e.g. TLR and IFN-α receptor alpha chain) substantially impacts the capacity of MGN1703 to reduce the intestinal HIV reservoir.

Potential limitations of our trial include the fact that the study was powered to detect differences in NK cell activation in peripheral blood [30] and that the study participants had the option of not consenting to the intestinal biopsy analyses for this sub-study but still take part in the clinical trial. Thus, the secondary endpoint analyses reported here may be unable to detect MGN1703-induced effects. A second potential limitation of this study is the single arm design. However, we established a longitudinal study where each individual’s baseline samples served as participant-matched control samples for measuring MGN1703-induced effects. This design accounts for the fact that “normal” human immune parameters vary greatly between individuals but remain relatively stable over time within a given person in the absence of an intervention [59]. Other potential limitations of the study are that NK cell activation in the intestines was not determined and that gene expression data were obtained using total intestinal mononuclear cell populations. The latter precludes ascertaining cell-type specific expression patterns, but does permit a broader characterization of the impact of MGN1703 in the intestines. The analyses reported here are based upon stool samples and sigmoid colon biopsies and it is possible that MGN1703 effects vary in other regions of the GI tract. Thus, future studies characterizing other regions of the GI tract are needed. Regarding microbiota analysis, the data reported here are for luminal [stool] microbiota composition and the TLR9 agonist might have a more robust effect on mucosal-associated microbiota, thus future studies need to not only interrogate mucosal-associated microbiota but also assess microbiota function.

In our trial we found that repetitive TLR9 stimulation induced systemic IFN-α production, without causing the severe side effects commonly observed in exogenous PEG-IFN-α treatment [30]. Furthermore, the potent type I interferon response, without concomitant inflammation, in the intestines of MGN1703 treated individuals did not perturb intestinal homeostasis but was associated with notable reductions in the tissue viral reservoir in some individuals. Nevertheless, a full understanding of the role(s) of IFN-α in HIV disease and treatment is lacking [60]. Multiple human and non-human primate studies have associated multiple benefits with exogenous IFN-α treatment of HIV/SIV infection [47, 48, 61–64]. More recent reports in humanized mice suggest that blocking the IFN-α receptor during HIV infection could be beneficial [65, 66], however, this strategy needs to be further explored in HIV-infected patients to fully evaluate such approaches.

In conclusion, we report the in vivo effects of the TLR9 agonist MGN1703 in the human colonic sigmoid mucosa and fecal samples. We demonstrate that intestinal homeostasis was maintained during treatment. These analyses have identified a tissue-specific phenotype whereby MGN1703 treatment did improve immune surveillance via a type I interferon response without inducing inflammation in the sigmoid colon. These data also demonstrate that baseline gene expression levels for TLR9 and IFNAR1 may significantly impact the downstream efficacy of a TLR9 agonist in the intestines. In an HIV treatment context, these findings strongly support the continued development of MGN1703 as an immune surveillance reactivator in HIV eradication approaches. These data also suggest that screening for TLR9 and IFNAR1 expression may aid in predicting drug efficacy regardless of the indication being treated (e.g. malignancies; infections; and inflammatory disorders).

Materials and Methods

Study Design and Participants

Between April and December 2015 we conducted the TEACH study (Toll-like receptor 9 enhancement of antiviral immunity in chronic HIV-1 infection: an investigator-initiated, single-arm, phase Ib/IIa clinical trial) (ClinicalTrials.gov ID number: NCT02443935) [30]. The study was approved by the National Health Ethics Committee, Denmark (case no. 1-10-72-10-15) and the Danish Medical and Health Authorities (case no. 2015014125; April 7, 2015) prior to enrollment. Each participant provided written informed consent before any study procedure. Fifteen HIV-infected individuals (13 males, 2 females) who had received ART for a minimum of 12 months and exhibited CD4+ T cell counts above 350 cells per µL and virological suppression (<50 copies/mL, at least two measurements per year) received subcutaneous injections of 60 mg MGN1703 (Lefitolimod; Mologen AG, Berlin, Germany) two times per week for four weeks while maintaining ART (Figure 1a) [30]. As part of the study described in [30], eleven individuals (9 males, 2 females) consented to the optional endoscopic collection of colonic sigmoid biopsies at baseline and during the last week of MGN1703 administration (Figure 1a). The sigmoid biopsies from the fourth dosing-week were collected 24 hours after the last injection of MGN1703. Eight of the 11 intestinal sub-study participants also supplied fecal samples for microbiota analyses at baseline and within two weeks of the last study drug dose.

Microbiome Analyses

Total DNA was extracted from human fecal samples (FastDNA bead-beating Spin Kit for Soil, MP Biomedicals, Solon, OH) collected at baseline and within two weeks of the last study drug dose. Fecal DNA was processed using high-throughput Illumina amplicon sequencing, testing primers described previously [67], and implementing a modified two-step targeted amplicon sequencing (“TAS”) approach [68]. Negative controls were used with each set of amplifications, which indicated no contamination. Raw sequence data (FASTQ files) were deposited in the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) under project PRJNA361500.

Raw FASTQ files for each sample were processed to merge reads, remove low quality data and chimeras, and perform annotation using the Greengenes 13_8 reference database, as previously described [69]. Each sample’s sequence set was rarefied to 14,500 sequences [70] and data were clustered into operational taxonomic units (OTU) at 97% similarity. Alpha diversity indices, Bray-Curtis similarity dendograms, and analysis of similarity (ANOSIM) measurements were calculated within the software package Primer7 [71]. ANOSIM was performed at the taxonomic level of family, using square-root transformed data. Differences in relative abundance of individual taxa (>1% of dataset), between baseline and treatment samples, were assessed for significance using the Kruskal-Wallis test controlling for false-discovery rate (FDR), implemented within the software package QIIME [72]. In SPSS (v.22), paired t-tests were used to analyze differences for parametric data satisfying normality test assumptions, via histograms, skewness/kurtosis, and Shaprio-Wilk test; Wilcoxon signed ranked paired t-test was used to analyze nonparametric data. Graphs were created using GraphPad Prism (v5.00) software.

Sigmoid Biopsies

Sigmoidoscopy was performed as previously described [73]. In brief, macroscopically normal intestinal tissue samples were taken approximately 35 cm from the anal verge. Two biopsies were randomly selected to be fixed in 4% PFA (16 hours at room temperature) and embedded in paraffin blocks. The remaining (~16 total) biopsies were collected in ice-cold Dulbecco’s phosphate buffered saline (PBS) and immediately placed on ice for transportation to the laboratory, where the tissues were immediately digested [enzyme cocktail in RPMI 1640 media: DNase I (1 kunitz unit/ml); Elastase (0.064 units/ml); (hyaluronidase I-S) (4–10 units/ml); and collagenase II-S (0.25 mg/ml)] at 37°C for 60 minutes with constant agitation. Epithelial cells and intraepithelial lymphocytes were not removed prior to digestion to ensure a more complete analysis of the effects of MGN1703 in the intestines. After digestion, intestinal tissue mononuclear cells were collected at the interface of a 40–70% Percoll density separation centrifugation. Viable cells were counted using trypan blue exclusion and cells were allocated to RNA sequencing (RNASeq) analysis, quantitative PCR for viral DNA, and flow cytometry analyses.

RNASeq Analysis of intestinal mononuclear cells

Intestinal tissue mononuclear cells were lysed in RLT Plus Buffer (Allprep isolation kit, Qiagen #80204) and lysates were stored at −80°C. Total RNA and total DNA were extracted simultaneously from the cell lysates using the Allprep DNA/RNA/miRNA Universal Kit (Qiagen) with on-column DNase treatment (Qiagen RNase-Free DNase Set) during the RNA extraction. RNA samples were quantified using Qubit 2.0 Fluorometer (Life Technologies, Carlsbad, CA, USA) and RNA integrity was checked with 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA). RNA sequencing libraries were prepared using the NEBNext Ultra RNA Library Prep Kit for Illumina following manufacturer’s recommendations (NEB, Ipswich, MA, USA). Briefly, mRNA was first enriched with Oligo d(T) beads. Enriched mRNAs were fragmented for 15 minutes at 94°C. First strand and second strand cDNA were subsequently synthesized. cDNA fragments were end repaired and adenylated at 3’ends, and a universal adapter was ligated to cDNA fragments, followed by index addition and library enrichment with limited cycle PCR. Sequencing libraries were validated using a DNA Chip on the Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA), and quantified by using Qubit 2.0 Fluorometer (Invitrogen, Carlsbad, CA) and quantitative PCR (Applied Biosystems, Carlsbad, CA, USA).

The sequencing libraries were multiplexed and clustered onto a flowcell. After clustering, the flowcell was loaded on the Illumina HiSeq 2500 instrument according to manufacturer’s instructions. The samples were sequenced using a 1x50bp Single Read (SR) configuration. Raw sequence data files (21,638 genes) were deposited in the NCBI SRA under project # SRP096929. Image analysis and base calling were conducted by the HiSeq Control Software (HCS) on the HiSeq 2500 instrument. Raw sequence data (.bcl files) generated from the Illumina HiSeq 2500 were converted into fastq files and de-multiplexed using the Illumina bcl2fastq v 1.8.4 program. One mismatch was allowed for index sequence identification. The sequencing reactions were conducted at GENEWIZ, LLC (South Plainfield, NJ, USA).

Quantification of Integrated HIV-1 DNA

Integrated HIV-1 DNA copies were quantified using two-step PCR reactions, as previously described [74]. Integrated HIV-1 DNA was pre-amplified with two Alu primers and a primer specific for the HIV LTR region, in addition to primers specific for the CD3 gene to determine cell count. Nested qPCR was then used to amplify HIV and CD3 sequences from the first round of amplification. Specimens were assayed with up to 500 ng cellular DNA in triplicate reaction wells and copy number was determined by extrapolation against a 5-point standard curve (3–30,000 copies), performed in triplicate using extracted DNA from ACH-2 cells. Given that baseline measurements of integrated HIV-1 DNA differ between individuals and given that baseline measures for each individual represent the control samples in this longitudinal study, integrated HIV-1 DNA data are presented as fold-change (treatment/baseline).

Cell Surface Flow Cytometry Analyses

For flow cytometric analyses, Fc-receptors were blocked with Human Trustain FcX blocking buffer (Biolegend, USA) for 5 minutes at room temperature, and cells were then stained with CD4 PE-Cy7 (SK3); CD8 PerCP-Cy5.5 (RPA-T8); CCR7 FITC (G043H7); CD45RA BV421 (HI100); HLA-DR PE (L243); CD38 BV605 (HIT2) (all from Biolegend, USA). Antibodies were mixed in BD HorizonTM Brilliant Stain Buffer (BD Biosciences, USA). Cells were incubated in the dark for 20 minutes at 4°C, then washed twice in PBS containing 2% FBS with centrifugation at 300 g for 5 minutes at 5°C. All sample data were acquired on a FACSVerse (BD) and analyses were performed using FlowJo software (v. 10.0.8r1) (Treestar, Inc., OR, USA).

Immunohistochemistry (IHC)

For IHC [75], 5 µm tissue sections were cut, mounted, dewaxed and rehydrated and then treated as described in Table S4. Epitope retrieval (15 min at 98–100°C) was followed by blocking (TBS+0.5% casein; 10 min), rinsing [TBS+0.05% Tween-20 (TBS-Tw)] and an incubation with primary antibody (60 min). Samples were washed in TBS-Tw and blocked using 1.5% H2O2 in TBS-Tw (5 min) (peroxidase block was not performed for MPO-IHC). Horseradish peroxidase-conjugated secondary antibodies (Table S4) were detected using ImmPact-DAB substrate (Vector Labs). All slides were washed in distilled water, counterstained with hematoxylin, mounted with Clear Mount (Invitrogen), dried for 60 min at 40°C and mounted with Permount (Electron Microscopy Sciences).

Quantification of IHC signal

A cell profile was counted as positive on blinded sections via visual inspection (magnification = 200x) when IHC-signal was associated with a nucleus and when the nucleus was sampled by randomly selected, non-overlaying 2-dimensional counting frame (25.600µm2) [73]. This was performed using compound light microscopy with a robotic stage that randomly selected non-overlaying areas for quantification (Visiopharm newCAST, version 5.0.3.1247). Sections representing the same individual and time point were cut from two separate blocks (biological replicates), and the average of the counted numbers in these replicates are the values presented in the figures. Validation of the IHC quantification strategy was performed according to standard methods before unblinding (Table S5) [73]. Quantitative IHC data are stratified into positive cell profiles per mm2 lamina propria or epithelium (insufficient events for IL-21 and MPO in the epithelium precluded conclusions for these parameters). Paraffin embedding and preparation cause tissue shrinkage [76, 77]. Thus, the area measurements reported here will not be the same as the true parameters present in vivo.

Statistics

Unless noted, all statistical tests were performed using an alpha of 0.05. Specific statistical tests (e.g. Spearman correlation) are indicated in the appropriate figure legend or table footnote. Equivalence evaluation and forest plot generation were performed using R statistical package [78]. For RNASeq analyses: RNAseq reads were mapped using Tophat to the reference genome hg19 [79, 80]. Gene expression quantification in fragments per kilobase of transcript per million mapped reads was calculated using Cufflinks suite including Cufflinks, Cuffmerge, Cuffquant and Cuffnorm [81, 82]. FDR were computed using the Benjamini-Hochberg procedure to adjust for multiple comparisons in the RNAseq data [83]. The clustering dendrogram depicting relatedness between gene expression profiles was generated using hierarchical clustering with complete linkage using R. The maximum number of unique tree permutations (184,756) were generated to determine statistical significance, applying a standard p<0.05 cutoff (significance was achieved when less than 5% of randomly-generated trees exhibited the same or greater level of segregation between treatment groups than the clustering dendrogram-depicted relatedness in actual expression data). The Slatkin-Maddison test was used to evaluate MGN1703 treatment-associated segregation in gene expression patterns (HyPhy_v2.2.4). Principal component analysis was performed using R on gene expression data for all 21,638 genes assessed. Each gene across samples was standardized with zero mean and unit variance.