Original Research ARTICLE
Immunofibrogenic Gene Expression Patterns in Tanzanian Children with Ocular Chlamydia trachomatis Infection, Active Trachoma and Scarring: Baseline Results of a 4-Year Longitudinal Study
- 1Clinical Research Department, Faculty of Infectious and Tropical Diseases, London School of Hygiene and Tropical Medicine, London, United Kingdom
- 2Kilimanjaro Christian Medical Centre, Moshi, Tanzania
- 3Department of Infectious Disease Epidemiology, Faculty of Epidemiology and Population Health, London School of Hygiene and Tropical Medicine, London, United Kingdom
- 4Support Services (Statistics), Medical Research Council Unit The Gambia, Fajara, Gambia
Trachoma, caused by Chlamydia trachomatis, is the world's leading infectious cause of blindness and remains a significant public health problem. Much of trachomatous disease pathology is thought to be caused indirectly by host cellular and immune responses, however the immune response during active trachoma and how this initiates progressive scarring is not clearly understood. Defining protective vs. pathogenic immune response to C. trachomatis is important for vaccine design and evaluation. This study reports the baseline results of a longitudinal cohort of Tanzanian children, who were monitored for 4 years in order to determine the immunofibrogenic and infectious correlates of progressive scarring trachoma. In this cohort baseline, 506 children aged 6–10 years were assessed for clinical signs, infection status and the expression of 91 genes of interest prior to mass azithromycin administration for trachoma control. C. trachomatis was detected using droplet digital PCR and gene expression was measured using quantitative real-time PCR. The prevalence of follicles, papillary inflammation and scarring were 33.6, 31.6, and 28.5%, respectively. C. trachomatis was detected in 78/506 (15.4%) individuals, 62/78 of whom also had follicles. C. trachomatis infection was associated with a strong upregulation of IFNG and IL22, the enrichment of Th1 and NK cell pathways and Th17 cell-associated cytokines. In individuals with inflammation in the absence of infection the IFNG/IL22 and NK cell response was reduced, however, pro-inflammatory, growth and matrix factors remained upregulated and mucins were downregulated. Our data suggest that, strong IFNG/IL22 responses, probably related to Th1 and NK cell involvement, is important for clearance of C. trachomatis and that the residual pro-inflammatory and pro-fibrotic phenotype that persists after infection might contribute to pathological scarring. Interestingly, females appear more susceptible to developing papillary inflammation and scarring than males, even at this young age, despite comparable levels of C. trachomatis infection. Females also had increased expression of a number of IFNγ pathway related genes relative to males, suggesting that overexpression of this pathway in response to infection might contribute to more severe scarring. Longitudinal investigation of these factors will reveal their relative contributions to protection from C. trachomatis infection and development of scarring complications.
Trachoma, a Neglected Tropical Disease caused by the bacterium Chlamydia trachomatis, remains the leading infectious cause of blindness worldwide. It is characterized by repeated conjunctival infection initiated early in childhood, triggering inflammation (trachomatous inflammation—intense and trachomatous inflammation—follicular) that drives scarring and trichiasis (in-turned eyelashes) and eventually blinding corneal opacification (Hu et al., 2013a; Taylor et al., 2014; Ramadhani et al., 2016a). The burden of this disease is high; current estimates indicate 200 million people live in trachoma endemic areas in 42 countries (World Health Organization, 2016). Approximately 1.9 million people are visually impaired or irreversibly blind from trachoma (Bourne et al., 2013). To meet this public health challenge, the WHO-led Global Alliance for the Elimination of Trachoma recommends the implementation of the SAFE strategy which tackles the disease at different stages: Surgery to correct trichiasis, Antibiotics to treat chlamydial infection and Facial cleanliness and Environmental improvements to suppress transmission of infection (Taylor et al., 2014).
Our understanding of this disease process is only partial (Hu et al., 2013a). There are few detailed long-term longitudinal studies that investigate the risk factors for and pathophysiology of scarring trachoma (Dawson et al., 1990; West et al., 2001; Wolle et al., 2009a,b; Burton et al., 2015). These studies have consistently found incident and progressive scarring to be strongly associated with clinically apparent conjunctival inflammation (Ramadhani et al., 2016a). However, only one study has prospectively examined the relationship between C. trachomatis infection and the development of incident scarring in younger people; this found that the detection of infection and/or severe inflammation on multiple occasions was associated with increased risk of new scarring (Wolle et al., 2009b). The relationship between chlamydial infection and progression of previously established scarring in adults has been prospectively examined in two longitudinal cohorts; neither found an association (Burton et al., 2015). This raises the possibility that progressive scarring may not be entirely dependent on continual re-exposure to C. trachomatis. It is possible that repeated chlamydial infection results in long-term physiological changes in conjunctival tissue responsiveness, such that other pro-inflammatory factors may stimulate on-going fibrotic responses (Kechagia et al., 2016). For trachoma control programmes, this raises the possibility that scarring may continue to progress after ocular chlamydial infection has been eliminated from the population. The implication of this is that services for managing incident trichiasis may be needed for many years.
Previously, we and others have explored the immunological correlates of the different stages of trachoma mostly through cross-sectional studies using a range of methodologies, including gene expression and protein analysis, immunohistochemistry and genome wide association studies (Bobo et al., 1996; Burton et al., 2004, 2011b, 2015; Skwor et al., 2008; Holland et al., 2010; Natividad et al., 2010; Derrick et al., 2013, 2016a,b; Hu et al., 2013a, 2016; Roberts et al., 2015). These studies have shown an increase in pro-inflammatory and matrix factors (IL1B (interleukin 1 beta), TNF (tumor necrosis factor), S100A7 (psoriasin), IL17A, IFNG (interferon gamma), perforin, IL12, IL10, CXCL5, CTGF (connective tissue growth factor) and MMP9 (matrix metalloproteinase 9) in individuals with active trachoma and/or chlamydial infection, indicating the involvement of a type 1 T helper (Th1) cell response, NK cell cytotoxicity and potentially Th17 cells. In individuals with trachomatous scarring and/or trichiasis, factors involved in innate pro-inflammatory responses and matrix remodeling (IL1B, CXCL5, S100A7, CTGF, MMP7, and MMP9) were upregulated.
To understand the immunofibrogenic correlates of progressive conjunctival scarring we conducted a long-term cohort study. Here we present the baseline findings for a large panel of factors, to define those to be examined in the prospective study.
This study was reviewed and approved by the Tanzanian National Institute for Medical Research Ethics Committee, the Kilimanjaro Christian Medical Centre Ethics Committee, and the London School of Hygiene and Tropical Medicine Ethics Committee. It adhered to the tenets of the Declaration of Helsinki.
The study was conducted in three adjacent trachoma endemic villages in Kilimanjaro and Arusha regions, Northern Tanzania. The villages are relatively remote, geographically neighbors and have similar patterns of life and traditions. This area is predominantly inhabited by people of the Maasai tribe. Pastoral activities are the main occupation. The area is dry for much of the year, except for the rainy season (February to May). Water supply is therefore limited, and largely depends on a long-distance water pipe scheme from Mount Kilimanjaro. Family units are organized in Bomas, with living huts arranged in a circle around a central animal enclosure, which is often characterized by a high density of flies.
In January 2012, we recruited a cohort of children aged 6–10 years from these villages. The cohort has subsequently been followed-up every 3 months for 4 years to investigate the pathogenesis of conjunctival scarring. All children, aged 6–10 years, who were normally resident in one of the three villages, were eligible for inclusion. We chose this restricted age group as we considered that they were more likely to show evidence of incident or progressive conjunctival scarring during the 4 years of the study. At the outset community meetings were held to introduce the study. Each household was then visited to meet the parents or legal guardians of children eligible for enrolment. A field worker explained the nature of the study in detail in either Swahili or Maasai language. There was an opportunity to discuss and ask questions. Finally, if the parent or guardian agreed to allow the child to be enrolled into the study this was documented on a consent form in Kiswahili, and witnessed by a third person.
Clinical Assessments and Sample Collection
The left eye of each child was examined by an ophthalmic nurse experienced in grading trachoma. The examinations were all performed under standardized conditions, using x2.5 loupes and a bright touch. The conjunctiva was anesthetized with a drop of preservative-free proxymetacaine hydrochloride 0.5%w/v (Minims®, Chauvin Pharmaceuticals Ltd, Surrey, UK). The eyelid was everted, examined and photographed (Nikon D90 with 105 mm Macro lens). Two conjunctival swab samples (Dacron polyester, Puritan Medical Products Company, Maine, USA) were collected for C. trachomatis detection and gene expression analysis. The swabs were passed across the upper tarsal conjunctiva four times, with a quarter turn between each pass. The first swab was placed directly into a tube containing RNAlater (Thermo Fisher Scientific, Massachusetts, USA) and the second into a dry tube. The samples were placed into a cool box. Later the same day the dry swab samples were stored directly at −80°C and the RNAlater samples kept at 4–8°C overnight and then stored at −80°C. Air control swabs were collected after every 50 samples by passing a swab 10 cm from a participant's everted eye, these were labeled and processed identically to participant samples.
Clinical signs were graded using the 1981 WHO Detailed Trachoma Grading System (FPC) (Dawson et al., 1981). This sub-divides the features into several four-point severity scales: follicles (F), papillary inflammation (P), and conjunctival scarring (C). This system corresponds to the WHO Simplified Trachoma Grading System in the following way: Trachomatous inflammation-Follicular (TF) is equivalent to F2/F3 and Trachomatous inflammation-Intense (TI) is equivalent to P3 (Thylefors et al., 1987). For the purpose of this study, we consider that both P2 and P3 represent clinically significant papillary inflammation, and refer to this as “TP” (Burton et al., 2015). We followed the widely used definition of “Clinically Active Trachoma”: TF and/or TI. The digital photographs were graded by an ophthalmologist experienced in trachoma assessment, using the FPC system, supplemented by a previously described system for fine grading of conjunctival scarring, that quantifies the extent of the conjunctiva involved (Hu et al., 2011). The ophthalmologist was masked to the infection status.
Extensive public health education about trachoma was provided to the community through village level meetings and during the house-to-house visits, including the importance of face washing for children and environmental improvements. Trichiasis surgery was provided free of charge within the community. Subsequently, all residents of the three villages have been offered three rounds of annual mass antibiotic treatment with oral azithromycin (and topical tetracycline ointment for infants under 6 months and pregnant women), during the course of the longitudinal study.
C. trachomatis Detection
To detect C. trachomatis, we extracted DNA from swab samples stored in dry tubes using the PowerSoil DNA Isolation Kit (Mo Bio Laboratories, California, USA), according to the manufacturer's instructions. The cells attached to the swabs were initially disrupted by bead beating to release their contents. C. trachomatis DNA was detected using a previously described droplet digital PCR assay (Roberts et al., 2013). All samples were tested for the C. trachomatis plasmid and the human gene RPP30. Samples in which C. trachomatis plasmid was detected were retested for C. trachomatis omcB (a single-copy gene on the chlamydial chromosome). Five microliters of template DNA was added per reaction. PCR reaction conditions were as follows: 95°C for 10 min, then 40 cycles of 95°C for 10 s and 60°C for 30 s and finally 98°C for 12 min. Droplets were then examined for fluorescence on a QX100 Droplet Reader (Bio-Rad, UK), providing a quantitative result.
Human Gene Expression Analysis
Total RNA was extracted from the swab stored in RNALater using a DNA/RNA Purification Kit (Norgen Biotek Corp, Canada), following the manufacturer's instructions. RNA was reverse transcribed using the SuperScript® VILO® cDNA Synthesis Kit (Life Technologies). Quantitative real-time PCR was performed to measure the relative abundance of a panel of human transcripts, using the TaqMan® Microfluidic 384-well Low-Density Array (TLDA) (Life Technologies) on a ViiA7 real-time PCR instrument (Thermo Fisher Scientific, Massachusetts, USA). These cards have separate channels for eight samples, each containing 48 wells, pre-printed with the assay primers and probes. We used two different TLDA designs and measured the expression of 91 different genes of interest for each sample. The selected genes are listed in Table 5. We measured expression of three different reference genes: HPRT1, GAPDH, and RPLP0. HPRT1 was selected as the most suitable for normalization as it was expressed at a relatively similar level to the majority of other transcripts of interest, whereas GAPDH and RPLP0 were very highly expressed. The choice of transcripts tested was informed by previous gene expression studies including transcriptome analysis experiments on samples from The Gambia, Ethiopia and Tanzania and by a genome-wide association study from The Gambia (Holland et al., 2010; Natividad et al., 2010; Burton et al., 2011b,c, 2015; Roberts et al., 2015).
Data were managed in Access and transferred to STATA v14 for analysis. The ophthalmologist's grading of the digital photographs was used for analysis. The ddPCR results were initially exported to R for analysis; the critical cut-off value for designating a positive result for each target was 0.2 plasmid copies/μl of eluted DNA, 0.2 omcB copies/μl and >0.3 RPP30 copies/μl. Samples were designated C. trachomatis positive if plasmid and human RPP30 were detected.
The associations between each of the clinical signs (follicular inflammation, papillary inflammation and scarring) and C. trachomatis infection were assessed using univariable logistic regression with infection as the outcome variable. The associations between scarring and each of the other clinical signs, infection, sex and age were investigated with logistic regression, using presence of scarring as the outcome variable. Initially, univariable logistic regression was performed using each exposure in turn, and then a multivariable analysis was performed in order to provide an unbiased estimate of the association between exposure and outcome adjusting for all clinical features.
The gene expression data were normalized relative to the expression of HPRT1 in the same sample, to adjust for variable concentrations of cDNA. This was done by the ΔCT method. Distributions of ΔCT values were plotted to assess them for normality. Fold change differences in gene expression between different phenotypic groups (Infection, follicular inflammation, papillary inflammation and scarring), vs. individuals without those clinical signs or infection, were calculated using the ΔΔCT method. Linear regression was used to compare the two groups in each comparison, with ΔCT values as the dependent variable and each phenotype group in turn as the independent variable, adjusting for age and sex. To take account of multiple comparisons we used a false discovery rate (FDR) of 5% (Benjamini and Hochberg, 1995). The linear regression analyses were repeated, but rather than taking the four phenotypic groups individually, they were all included as exposures in a multivariable model with age and sex, providing estimates of the fold-changes associated with each phenotype.
A heatmap of ΔCT values for each gene of each participant was produced, to indicate visually whether the expression of individual genes was consistently of higher or lower expression in groups of individuals with clinically active Trachoma (TF and/or TI) and/or infection. A principal component analysis (PCA) was performed on the gene expression data and clinical disease and infection status were overlaid on plots in order to visually indicate any associations between the first two components and infection, active trachoma and scarring. These associations were tested formally using a logistic regression analysis, with the principal components as the exposure variables. For the heatmap and PCA, only individuals with complete expression data could be included. Therefore, so as not to lose too many individuals, any target with >5% of observations missing was excluded from these analyses, as well as the lasso regression and co-expression analyses described below. Any individuals with missing observations on the remaining targets were excluded from these analyses also.
For each of infection, active trachoma and scarring, a multivariable logistic regression was performed using this subset of gene expression levels as exposure variables, with the aim of retaining those expressions most strongly associated with the outcome of interest, adjusting for all other expression levels. Due to the large number of exposure variables, a penalized logistic regression was performed using the lasso technique (Tibshirani, 1997). The list of genes most strongly associated with each outcome was entered into ConsensusPathDB (http://cpdb.molgen.mpg.de/) gene set over-representation analysis, with all genes tested entered as background. Enriched pathway-based sets with a minimum overlap between pathway and input gene list of two and a p < 0.05 were identified.
A network graph based on the specimen-to-specimen Pearson correlation was generated using miru (https://kajeka.com/miru). The overall expression correlation matrix and graph were constructed from the raw cycle threshold expression values. The filtered dataset including only individuals and genes with complete expression data was used. Pearson correlation coefficients (r) ≥0.85 were retained and used as cut-offs in network construction. Nodes in the graph are individual mRNA transcripts linked by an edge if r was ≥0.85. The graph was then clustered using a Markov Clustering algorithm with an inflation value of 2.2. The partitioned clusters of expression contain sets of transcripts that exhibit a very high degree of co-expression across the sample set. The co-expression clusters or modules were then investigated for enrichment at the pathway level using ConsensusPathDB as described above, using all genes tested as background. A single value for each cluster of each individual was defined using the first principal component and this value was used to investigate differential expression associated with disease and infection phenotypes.
A total of 506 children between the ages of 6 and 10 years were recruited and examined. Their demographic and clinical characteristics are described in Table 1. There were almost equal numbers of male and female children, with a mean age of 7.5 years. Most (97.4%) of the children were of the Maasai ethnic group.
Trachomatous inflammation-Follicular (TF, F2/F3) was present in 170 (33.6%) children (Table 1). However, a further 135 (26.7%) had evidence of mild follicular conjunctivitis (F1). Significant conjunctival papillary inflammation (TP, P2/P3) was observed in 160 (31.7%) children, of which 64 (12.7%) had intense papillary inflammation (TI, P3).
C. trachomatis Infection and Active Trachoma
C. trachomatis plasmid was detected by ddPCR in 78 (15.4%) individuals. There was a strong association between the presence of TF and the detection of C. trachomatis plasmid; 62/78 (79.5%) of C. trachomatis positive individuals had TF (Table 2). Of the 16 infected children who did not have TF, 11 (68.8%) had signs of mild follicular conjunctivitis (F1). However, only 62/170 (36.5%) individuals with TF had detectable infection. There was a similar strong association between the presence of TP and chlamydial infection (Table 2).
Table 2. C. trachomatis plasmid detection and clinical signs; both the simplified and detailed WHO “FPC” grading system.
Infection load and the concentration of RPP30 endogenous control DNA were quantified. The concentration of plasmid DNA ranged from 0.22 to 3023 copies/μl, omcB ranged from 0.2 to 742 copies/μl and RPP30 ranged from 0.32 to 1,081 copies/μl of eluted DNA. All samples had detectable RPP30. There was a marked positive trend in both the proportion infected and the load of infection with increasing F-Score and increasing P-Score (Table 2). There was a consistent ratio of plasmid to omcB copies/μl: geometric mean ratio of 4.57 to 1 (95% CI 3.91–5.34) and a median of 4.79 to 1 (95%CI 4.01–5.21). The correlation between plasmid and omcB copies/μl is illustrated in Figure 1. The ratio did not vary with either the follicular or papillary clinical severity scores (data not shown).
Figure 1. The relationship between paired measures of copies of Chlamydia trachomatis plasmid/μl vs. omcB/μl in the same sample. Correlation R-squared = 0.989.
Conjunctival scarring was relatively frequent, being found in 144 (28.5%) individuals; this was mostly mild scarring (Table 1). There was a significant association between the detection of C. trachomatis and scarring [p = 0.002, OR = 2.23 (95%CI 1.4–3.7)]. The load of infection also increased with scarring severity (Table 2). There were univariate associations between conjunctival scarring and TF, TP (P2/P3), C. trachomatis detection and female sex (Table 3). However, in a multivariable model only TP, female sex and age were independently associated with scarring.
Table 3. Univariable and multivariable associations between conjunctival scarring and other clinical features, C. trachomatis infection, sex and age.
There was strong evidence (p = 0.0047) that even at this young age females are proportionately more likely to have signs of established scarring than males (Table 4). Similarly, there was some evidence (p = 0.048) that females had slightly increased odds of TP. No evidence was found of an association between sex and the odds of TF or chlamydial infection (p = 0.12 and 0.25, respectively) (Table 4).
Table 4. The relationship between sex and (i) clinical signs, (ii) Chlamydia trachomatis infection and (iii) infection load.
Conjunctival Gene Expression
The expression of 91 target genes was measured relative to that of HPRT1. Out of the 506 participants, samples from 12 individuals failed the qPCR for all target genes, leaving 494 individuals for gene expression analysis. The number of participant samples with detectable expression for each target is shown in Supplementary Table 1 (N). Multivariable linear regression analysis was performed using data for all 91 genes from 494 individuals.
Four sets of comparisons were performed: (1) infected v non-infected, (2) TF v no TF, (3) TP v no TP, and (4) TS v no TS. The relative Fold Change (FC) between these paired groups, calculated by the ΔΔCT method, are presented in Supplementary Table 1 along with the p-values for the linear regression (adjusted for age and sex) for each comparison.
As clinical signs and infection may be highly correlated in their association with the expression of specific targets we tested for independent associations with each target's expression using multivariable linear regression models, Table 5. More genes were significantly differentially expressed in individuals with C. trachomatis infection (55/91) relative to individuals with TP (33/91), TF (11/91), and TS (17/91) in these multivariable models (P < 0.05). C. trachomatis Infection was associated with increased expression of multiple cytokines and chemokines (CCL2, CXCL13, CCL18, CSF2, FOXP3, IFNG, IDO1, IL1B, IL6, IL8, IL10, IL12B, IL17A, IL19, IL21, IL22, IL23A), cell cycle components (CDC25C, TTK, TYMS), matrix modifiers (MMP7, MMP9, MMP12, TGFB1), NK cell markers (CD247, NCR1, NCAM1) and intracellular signaling molecules/regulators (CD274, IKZF1, RHOH, SAMSN1, SERPINB3, SERPINB4, STAT1, STAT4, SOCS1, TBX21), Table 5. Infection was associated with particularly marked increases in expression of IFNG (FC: 7.12, p-value: 1.32E-42) and IL22 (FC: 6.49, p-value: 1.6E-23). There was reduced expression of several factors including mucins (MUC1, MUC4, MUC5AC, MUC7) and SPARCL1.
Table 5. Multivariable linear regression models for conjunctival gene expression in the presence of clinical signs, C. trachomatis plasmid, female sex and age.
Clinical signs of active trachoma were associated with increased expression of multiple cytokines, chemokines, antimicrobial peptides, matrix modifiers, following a similar pattern to infection (Table 5). TP was associated with more substantial increases in a wider range of factors than TF (Table 5): antimicrobial peptides (DEFB4A, S100A7), cytokine/chemokines (CXCL5, CXCL13, CCL18, CCL2, IL1B, IL8, IL10, IL17A, IL19, IL23A, PTGS2), matrix modifiers (MMP7, MMP9, MMP12) and intracellular signaling molecules/regulators (CD274, SAMSN1, SERPINB3, SERPINB4, SOCS3). Several cytokines and chemokines (CCL18, CXCL13, IL19, and S100A7) had a FC >2 in children with TP. Conjunctival scarring was associated with modestly increased expression (>1.2 FC) and borderline significance of several chemokines and cytokines (CCL2, CCL20, CXCL5, CSF3, IL1B, IL8, IL10, IL21, IL33, FGF2) in the multivariable linear regression models (Table 5). SPARCL1 was significantly downregulated in TF and TP but not TS.
Only 108 individuals had detectable gene expression for all 91 targets, so in order to retain a reasonable number of individuals in the PCA and lasso analyses, which require complete records for all individuals, targets which were missing >5% of observations were excluded. This resulted in the exclusion of eight transcripts (USP6, IL13, CTNND2, FGF2, SERPINB4, IL22, PREX2v1, PREX2v2) for which there was no detectable expression in >5% individuals. In individuals with detectable expression of these eight targets, raw CT values were generally very high (CT > 33). These eight targets were excluded from the heatmap, PCA, lasso regression and co-expression network analyses, retaining a complete expression dataset for 83 genes. Among the 494 individuals, 36 were excluded as they had missing data among the 83 included targets, resulting in a final dataset of 458 individuals and 83 genes. Excluding further genes offered minimal gains in sample size. A heatmap of ΔCT values separated by clinical and infection phenotype is presented in Figure 2. Infected individuals with or without active trachoma had a visibly distinct pattern of gene expression, whereas the ΔCT values of normal healthy controls and individuals with active trachoma in the absence of infection appeared more similar overall.
Figure 2. Heatmap of the expression of each of the 83 retained genes in 458 individuals. Each row represents an individual, with individuals grouped as per their disease and infection status and the relevant group indicated on the right: Group (A) Infected, Active Trachoma; (B) Infected, no Active Trachoma; (C) Uninfected, Active Trachoma; (D) Uninfected, no Active Trachoma. ΔCT values of genes are represented in the columns, with genes that tend to be expressed together located next to one another. Red means that an individual has a higher expression of that gene, blue, a lower expression, relative to that individual's expression of the endogenous control gene HPRT1.
PCA analysis was performed for 458 individuals that had complete gene expression data for 83 targets and active trachoma/infection and scarring status were overlaid on plots, Supplementary Figures 1A,B. Principal component 1 (PC1) was strongly associated with infection and active trachoma; for a one unit increase in PC1 there was an estimated 35% decrease in the odds of infection (95% CI 27–41%, p < 0.001) and an estimated 23% decrease in the odds of active trachoma (95% CI 9–28%, p < 0.001). PC1 loadings showed that as genes in the chemokines/cytokines and regulator/signaling pathways groups were up-regulated, PC1 became smaller, suggesting that the odds of infection/active trachoma was greater among individuals where these genes are up-regulated. The expression of mucins, CDH1, SPARCL1, and S100A4 showed the reverse relationship, as when these genes were down-regulated the PC1 was smaller. Principal component 2 (PC2) was associated with increased odds of active trachoma alone, with a unit increase in PC2 being associated with an estimated 15% increase in the odds (95% CI 7–24%, p < 0.001), however there was no evidence of an association between PC2 and infection. There was evidence of a weaker association between PC1 and scarring, with an estimated 10% decrease in the odds of scarring (95% CI 6–14%, p < 0.001) per unit increase in PC1, Supplementary Figure 1B. There was no evidence of an association between scarring and PC2.
Three lasso logistic regressions were performed, first using infection (with or without active disease) as the outcome and the 83 targets, age and sex of 458 individuals as exposures. Then, excluding all infected individuals, active trachoma (uninfected only) was used as the outcome against uninfected individuals without active trachoma and the same 83 targets, age and sex as the exposures of interest. Finally, the regression was performed using scarring as the outcome and the same set of exposures in all 458 individuals. The transcripts most strongly associated with each of the three outcomes, and biological pathways with members that were over-represented in each relative to the background (83 targets) are listed in Table 6. Genes most strongly associated with infection were associated with Th1 cell development and NK cell pathways, whilst genes most associated with active trachoma in the absence of C. trachomatis infection were enriched for members of leukocyte transendothelial migration, TNF receptor, MMP activation, collagen formation and collagen fibrils assembly pathways. Genes most strongly associated with the TS phenotype were enriched for members of immunoregulatory interactions between lymphoid and non-lymphoid cells.
Table 6. Lasso logistic regression of gene expression by infection, active disease (in only uninfected individuals) and scarring status.
Networks of co-expression, independent of differential expression, were explored in the complete expression data of 83 transcripts and 458 individuals using Miru. The undirected graph contained 51 nodes connected by 121 edges. Markov clustering partitioned the network into 5 separate clusters of co-expressed genes that accounted for 38/51 (74.5%) of the transcripts in the original network. Cluster 1 contained 11 co-expressed transcripts, cluster 2 contained 10 transcripts, there were 6 transcripts in each of cluster 3 and 4, and 5 transcripts in cluster 5 (Supplementary Table 2A), the remaining nodes were not connected in the graph. The most connected gene transcript or hub in each module was MyD88/REL (cluster 1), IKZF1/VIM (cluster 2), STAT1 (cluster 3), SRGN (cluster 4), and IL10/IL17A/IL21 (cluster 5).
Pathway enrichment analysis was performed for each cluster relative to the background of 83 transcripts (Supplementary Table 2B). ΔCT values for the transcripts in each cluster were collapsed into a single value for each of the 458 individuals, represented by the first principal component. Multivariable linear regression was then performed for each cluster, using the first principal component as the outcome variable and TF, TP, TS, infection, age, and sex as independent variables (Supplementary Table 2C). Cluster 1 was enriched for cell cycle and apoptosis pathways, however the combined expression value was not differentially expressed in any condition. Cluster 2 was significantly associated with infection and TF and was marginally enriched for the retinoblastoma pathway. Cluster 3 was highly enriched for Th1/2 cell differentiation, IFNγ, IL-12, CD8, and NK cell signaling pathways and was strongly associated with infection and age and marginally associated with sex, but interestingly not with TP or TF. Conversely, cluster 4 was enriched for TNF, inflammasome, TLR and NFkB signaling pathways and was associated with TP, TS, and age but not infection. Cluster 5 was enriched for the allograft rejection pathway and was associated with infection, TP, TS and age.
Clinical Disease and C. trachomatis Infection
In this study, conducted in communities prior to azithromycin MDA, there was a relatively high TF (34%) and moderate C. trachomatis (15.4%) prevalence in 6–10 year olds. There was a strong relationship between infection and clinical signs. The large majority (79.5%) of individuals with C. trachomatis detected by ddPCR had TF (F2/F3). Moreover, most individuals (68.8%) with detectable chlamydial infection but without TF had a mild follicular conjunctivitis (F1). There were few cases of infection in the absence of clinical signs. This is consistent with earlier studies which found a stronger correlation as the underlying prevalence of infection increased (Ramadhani et al., 2016b). However, in common with many other studies reporting the relationship between disease and infection, a minority (36%) of individuals with TF had detectable infection, probably due to the shorter duration of infection episodes relative to disease (Grassly et al., 2008; Burton et al., 2011a,c; Lee et al., 2014). With increasing severity of both TF and TP scores, there was an increasing proportion with infection and increasing loads of infection. Consistent with earlier studies, increasing load appeared to be more closely related to increasing TP rather than TF (Burton et al., 2003; Solomon et al., 2003; Michel et al., 2011; Derrick et al., 2016a).
In the present cross-sectional study conjunctival scarring was associated with TP, female sex and increasing age but not C. trachomatis infection, which is consistent with several earlier cohort studies (Dawson et al., 1990; West et al., 2001; Wolle et al., 2009b; Burton et al., 2015; Hu et al., 2016). However, longitudinal data supporting an association between TF or C. trachomatis infection and progressive scarring are limited (Ramadhani et al., 2016a). We found females had more TP and TS than males, whereas C. trachomatis infection and TF prevalence were not significantly different between the sexes (Table 4). This might suggest that there are additional determinants of TP and TS development beyond the initial C. trachomatis infection. It is generally believed that the difference in the proportion of males and females developing scarring sequelae of trachoma is attributable to a greater life-time exposure to repeated C. trachomatis infection among females (Taylor et al., 2014). While that may well be the case, the data from this study also suggest that even at a relatively young age, females appear more susceptible to developing TP and TS than males, despite comparable levels of C. trachomatis infection. This finding is consistent with the observation that in general, females generate stronger immune responses than males, making them more susceptible to diseases resulting from immune-mediated pathology (Klein and Flanagan, 2016).
Among infected individuals, the C. trachomatis plasmid:omcB ratio (4.57:1) was very similar to the ratio reported from a range of ocular and genital C. trachomatis isolates (mean 4.0:1) and from a population-based trachoma study in Guinea Bissau (5.3:1) (Pickett et al., 2005; Last et al., 2014). In common with the study from Guinea Bissau, we did not find evidence of variation in the plasmid copy number with increasing disease severity scores.
Host Gene Expression
IFNG was strongly upregulated in association with C. trachomatis infection, consistent with previous findings in trachoma endemic populations (Burton et al., 2004; Holland et al., 2006; Natividad et al., 2010). Several NK cell markers were also upregulated in association with C. trachomatis infection: NCR1, CD56 (NCAM1), and CD247. IFNγ, NCAM1, and NCR1 were among the genes most strongly associated with increased odds of C. trachomatis infection, and Th1 and NK cell pathways were enriched in this gene set (Table 6). A strong IFNG response and transcriptional signatures associated with NK cell activation and cytotoxicity were closely associated with C. trachomatis infection and active trachoma in a West African population (Natividad et al., 2010), and infiltrates suggestive of NK cells were identified in scarred conjunctival biopsy tissue from individuals with trachomatous trichiasis from Tanzania (Hu et al., 2016). In peripheral blood from adults with trichiasis and controls that were stimulated ex-vivo with C. trachomatis antigens, NK cells and to a lesser extent CD8+ T cells were the major sources of IFNγ with <50% of the IFNγ produced by CD4+ T cells (Gall et al., 2011). In a non-human primate live-attenuated ocular chlamydial vaccine study, depletion of CD8+ cells abrogated protective immunity, however the contribution of CD4+ and NK cells was not investigated (Olivares-Zavaleta et al., 2014). IFNγ is a crucial component of the Th1 response to restrict C. trachomatis survival via tryptophan depletion, macrophage activation and promoting leukocyte adhesion, migration and activation (Beatty et al., 1994; Mosser and Edwards, 2008; Redgrove and McLaughlin, 2014). Interestingly, whilst IFNG was strongly upregulated in individuals with C. trachomatis infection, there was much less upregulation in individuals with disease but not infection, suggesting that the IFNG response is quickly down-regulated once C. trachomatis has been cleared. Supporting this observation, cluster 3, which was enriched for Th1, NK, IL-12, and IFNγ pathways, was strongly associated with infection but not with TF or TP. Cluster 3 expression was also marginally associated with female sex, and IFNG, IDO1, IL12B, and IL19 were each upregulated in females relative to males (Table 5). Both the number of activated CD4+ cells and their relative expression of pro-inflammatory genes including IFNG have previously been shown to be higher in females than males (Hewagama et al., 2009; Sankaran-Walters et al., 2013). Our data suggest that increased expression of IFNG and related pathways in females may contribute to the observed differences in tissue responsiveness (TP and TS) between the sexes in the context of equal infection exposure. In support of this hypothesis, a longitudinal study conducted in The Gambia that visited participants every 2 weeks found that the likelihood of developing infection or active disease and the duration of infection were independent of sex, however males had a reduced duration of active disease relative to females (Grassly et al., 2008). Whilst Th1 responses and IFNG are essential for controlling chlamydial infection, excessive responses are thought to cause collateral tissue damage and contribute to chlamydial pathology (Van Voorhis et al., 1997; Rank et al., 2000). It is therefore possible that a female predisposition to overproduce Th1 responses and IFNG may tip the balance from protection toward pathological tissue damage. However, further research is required to clarify this hypothesis and a direct causal relationship cannot be determined from this cross-sectional analysis.
The cytokine genes IL1β, IL6, IL17, IL21, IL22, IL23, and TGFβ1 were strongly expressed in individuals with C. trachomatis infection. Th17 cells are generated in the presence of TGFβ1, IL-1β, IL-6, IL-21, and IL-23 and classically produce IL-17, IL-22, and IL-21 (Manel et al., 2008; Zielinski et al., 2012). IL21 expression was also upregulated in individuals with TF and TS. Individuals with TP significantly upregulated IL17 and IL23 expression and IL21 at marginal significance. The combined expression of cluster 5, containing IL10, IL17A, IL21, CXCL13, and MMP12, was significantly associated with C. trachomatis infection, TF, TP, and TS. Th17 cells have important roles in mucosal immunity against bacteria and fungi, however they are also associated with the inflammatory pathology of diseases such as Crohn's disease, psoriasis and rheumatoid arthritis (Tesmer et al., 2008). Epithelial cells stimulated by IL-17 and IL-22 upregulate chemokines, cytokines, antimicrobial peptides and growth factors such as IL-6, IL-8, S100A7, and GM-CSF (CSF2) (Wolk et al., 2006; Eyerich et al., 2010); these transcripts were upregulated in individuals with C. trachomatis infection (S100A7 at marginal significance). IL8 and S100A7 were also upregulated in individuals with TP. IL-8 and S100A7, amongst other factors induced by Th17 cells, are chemotactic for neutrophils, an influx of which is associated with C. trachomatis infection and contributes to pathology in animal models of urogenital and ocular infection (Frazer et al., 2011; Lacy et al., 2011). IL-22, which was highly upregulated with infection, is also important in the maintenance of epithelial health and barrier function (Radaeva et al., 2004; Rutz et al., 2013). The relative contributions of IL-17 and associated responses to protection from C. trachomatis infection and maintenance of barrier function vs. pathological tissue damage remain unclear. In murine models of genital C. trachomatis infection higher levels of IL-17/Th17 cells have been associated with greater pathology (Lu et al., 2012; Vicetti Miguel et al., 2016), whereas in IL-23 knock-out mice which had greatly reduced levels of IL-17 and IL-22 there was no difference in infection burden or oviduct pathology compared to wild type mice (Frazer et al., 2013). Longitudinal investigation of these factors is expected to further our understanding of the role of Th17 cell associated cytokines in human chlamydial pathogenesis.
STAT1 and STAT4, which promote Th1 cell development (Nishikomori et al., 2002; Mikhak et al., 2006; Zhu et al., 2010), were upregulated in response to C. trachomatis infection. STAT1 and STAT4 have previously been shown to be associated with active trachoma and ocular C. trachomatis infection, and STAT1 was upregulated in cervical epithelial cells infected with C. trachomatis in vitro which was found to inhibit chlamydial growth (Lad et al., 2005; Natividad et al., 2010). STAT3, which is required for Th17 development (Zhu et al., 2010; Yang et al., 2011), was not differentially regulated, contrasting the upregulation of Th17 cytokines. SOCS1, IL10, and FOXP3 were upregulated in individuals with C. trachomatis infection and SOCS3 and IL10 were upregulated in individuals with TS and TP, indicative of inflammatory regulation. Upregulation of IL10 expression has previously been reported in individuals with active trachoma and a genetic polymorphism in IL10 was associated with trachomatous scarring (Mozzato-Chamay et al., 2000; Burton et al., 2004; Natividad et al., 2005; Faal et al., 2006; Skwor et al., 2008). However, the role of Tregs and IL-10 in trachoma remains unclear; dampening inflammation might result in less tissue damage, but conversely it might also prolong survival of C. trachomatis (Zhang et al., 2009).
MMP9, MMP12, TGFβ1, and PDGF were significantly upregulated in individuals with C. trachomatis infection. MMP7, MMP9, and MMP12 were upregulated in individuals with TP, MMP9 was upregulated in TF and FGF2 and PDGF were upregulated in individuals with TS. MMP7 was downregulated in infected individuals. MMP7 has consistently been found to be upregulated in scarring trachoma and trichiasis alongside MMP9 and MMP12 but there is less evidence for a role of MMP7 in active trachoma (Holland et al., 2010; Burton et al., 2011b, 2015; Hu et al., 2012). Degradation of the extracellular matrix and basement membrane by MMP9 and MMP12, in part derived from macrophage and neutrophil degranulation (Shapiro et al., 1993; Kjeldsen et al., 1994), facilitates connective tissue remodeling and leukocyte migration into the epithelium (Ozdemir et al., 1999; Zeng et al., 1999; Misko et al., 2002). Whilst MMPs are an essential part of the wound healing process (Nwomeh et al., 1998; Wolf et al., 2017), MMP9 overexpression has been linked to aggressive contractile scarring in proliferative vitroretinopathy and failure of rat corneal cells to re-epithelialize following corneal injury (Fini et al., 1996; Kon et al., 1998).
SPARCL1 was down-regulated in individuals with C. trachomatis infection, TF and TP, consistent with previous studies (Hu et al., 2012; Burton et al., 2015). SPARCL1 codes for a secreted glycoprotein which regulates extracellular matrix (ECM) and cell-matrix adhesion (Girard and Springer, 1996) SPARCL1 inhibits cell proliferation and migration and reduced SPARCL1 expression has been linked to cancer metastases (Nelson et al., 1998; Sullivan and Sage, 2004). In a murine corneal injury model, knock-out of SPARCL1 led to excessive accumulation of collagen type IV from myofibroblasts and irregular fibrotic ECM formation (Chaurasia et al., 2013). Addition of exogenous SPARCL1 rescued the wild type phenotype with a cessation of collagen IV production, replacement with collagen I from keratinocytes and restoration of normal ECM architecture. Consistent with this, immunohistochemistry and histology studies using tissue from individuals with scarring trachoma have demonstrated a reduction in collagens I and III, an increase in collagens IV and V and progressive disorganization of collagen fibers in scarring trachoma relative to healthy conjunctival tissue (Abu El-Asrar et al., 1998; Hu et al., 2013b; Derrick et al., 2016b). Genes that were most strongly associated with active trachoma in the absence of C. trachomatis infection were enriched for pathways of leukocyte transendothelial migration, TNF receptors, MMP activation, collagen formation and assembly (Table 6). Active trachoma signs often persist after the clearance of C. trachomatis infection (Grassly et al., 2008; Burton et al., 2011a,c; Lee et al., 2014) and inflammation is linked to progressive scarring (Ramadhani et al., 2016a). We found that TP but not infection was associated with TS in the multivariable model (Table 3). These results suggest that the inflammatory phenotype following the clearance of C. trachomatis infection is characterized by collagen formation and assembly, MMP-mediated matrix remodeling and leukocyte influx, and that these pathways might be causative factors in the scarring process.
We found significantly reduced expression of MUC1, MUC4, MUC5AC, and MUC7 in individuals with C. trachomatis infection, and MUC5AC expression was reduced in those with TP. MUC1 and MUC4 are primarily expressed by apical cells of the stratified epithelia, MUC7 by lacrimal gland epithelia and MUC5AC by goblet cells (Gipson and Argueso, 2003). MUC7 expression was also significantly reduced in females relative to males. MUC5AC and MUC7 down-regulation has previously been reported in individuals with trachomatous trichiasis, however expression of MUC1 and MUC4 was increased in this later stage of disease (Burton et al., 2011b). IL-22 stimulates the upregulation of MUC4 and MUC1, therefore, one might expect a corresponding increase in mucin expression in infected individuals (Turner et al., 2013). This inconsistency and the broad reduction in mucin expression during infection may reflect a loss of epithelial homeostasis and a potential reduction in barrier function. This in turn could result in increased access of bacteria to the epithelial cell surface, promoting inflammation. In patients with ulcerative colitis, reduction of goblet cells, mucin production and microbial diversity were identified as causative factors in disease etiology (Alipour et al., 2016). Microbial dysbiosis can itself be caused by inflammation, creating a positive feedback loop (Winter et al., 2013). Decreased diversity has previously been found in the conjunctival microbiome of individuals with scarring trachoma relative to healthy controls (Zhou et al., 2014). Furthermore, penetration of the mucin barrier and contact between epithelial cells and bacteria was demonstrated in ulcerative colitis patients with acute inflammation, whereas the mucin barrier of healthy controls was impenetrable (Johansson et al., 2014). It is tempting to speculate that reduced barrier function resulting from inflammation and loss of mucins causes microbial dysbiosis and contact between bacteria and epithelial cells, leading to persistent conjunctival inflammation, which in turn contributes to trachomatous scarring.
Consistent with our previous studies in various clinical stages of trachoma (Burton et al., 2011b,c, 2012, 2015; Hu et al., 2012), the antimicrobial peptide S100A7 was upregulated in individuals with TP and was marginally associated with C. trachomatis infection. S100A7 is secreted by epithelial cells in response to microbial challenge and expression can also be stimulated by IL-22 (Wolk et al., 2006). It is chemotactic for neutrophils and interferes with pathogen membrane permeability and enzyme function (Ganz, 2003; Guilhelmelli et al., 2013). The upregulation of S100A7 in various clinical stages of trachoma could support the loss of epithelial barrier function hypothesis described above and could contribute to chronic inflammation.
There was differential expression of several cell cycle genes (CD53, CDK13, CDC25C, PHYH, TYMS, TTK, TP53) in individuals with C. trachomatis infection and/or TP (Table 5). These expression changes could reflect regulation of the cell cycle in response to infection or inflammation, such as the proliferation or contraction of infiltrating inflammatory cells and cells effecting the wound healing response. Studies have reported that progression of the cell cycle is slow in cells infected with C. trachomatis (Kun et al., 2013; Elwell et al., 2016) and C. trachomatis is known to interfere with host cell apoptosis as a survival mechanism (Fischer et al., 2004; Ying et al., 2007; Siegl et al., 2014). Cluster 1 was enriched for genes in apoptosis pathways, however the combined expression of this cluster was not associated with clinical phenotypes or infection.
Our results, illustrated in a graphical summary (Figure 3), suggest that IFNG and IL22 were acutely upregulated in response to C. trachomatis infection and expression appeared to be reduced following the clearance of infection, despite the persistence of inflammation. This might suggest that the sources of these cytokines are important in the clearance of infection and that the inflammatory phenotype differs once infection has been cleared, shifting away from a Th1/NK cell dominated response. The upregulation of Th17 cell-associated cytokines in infection and TP phenotypes implies the involvement of Th17 cells, although it remains unclear whether they contribute more to protection or pathology. Epithelial cells appear to be highly active, with possible loss of epithelial homeostasis and barrier function resulting in loss of mucins and a pro-inflammatory anti-microbial phenotype, perhaps exacerbating inflammation. The upregulation of growth factors and MMPs, derived in part from infiltrating macrophages and neutrophils, and the down-regulation of SPARCL1 are likely to contribute to matrix remodeling, collagen deposition and conjunctival scarring in the post-infection inflammatory phase. Females were more susceptible to TP and TS and although this was not associated with a significant difference in infection prevalence, there was some evidence of increased IFNG related pathway activity in females. A caveat of this study is that gene expression may not correlate directly to protein expression and this will require further work. Despite this, we anticipate that longitudinal investigation of these transcriptional signatures and relation to infection and ongoing clinical signs will enhance our understanding of their relative contributions to protection and pathology.
Figure 3. Graphical summary illustrating the host immune pathways hypothesized to be associated with different trachoma phenotypes. Genes most strongly associated with each phenotype are shown in italics with the direction and strength of expression illustrated by arrows. Pathways that were most enriched for each phenotype are shown. Ct EBs, C. trachomatis elementary bodies; N, epithelial cell nucleus; NK, Natural killer cells; Th1, Th1 T cells; Th17, Th17 T cells; PMN, neutrophil; Mϕ, macrophage; red cross-hatching, fibrosis. The mucus layer is shown in green, the epithelial cell layer in blue, and the stroma in beige.
Substantial contributions to the conception: MB, MH, DCM, and RB. Design of the work: MB, AR, MH, PM, TM and CR. Analysis: AR, MB, DM, MH, TD, and DJ. Interpretation of data for the work: AR, MB, MH, and TD. Drafting the work: AR and MB. Revising it critically for important intellectual content: AR, MB, MH, TD, DCM, RB, DM, CR, PM, TM, and DJ. Final approval of the version to be published: AR, MB, MH, TD, DCM, RB, DM, CR, PM, TM, and DJ. Agreement to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved: AR, MB, MH, TD, DCM, RB, DM, CR, PM, TM, and DJ.
This work is supported by the Wellcome Trust grant number 098481/Z/12/Z.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We express our gratitude to the field team and to the study participants.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fcimb.2017.00406/full#supplementary-material
Abu El-Asrar, A. M., Geboes, K., Al-Kharashi, S. A., Al-Mosallam, A. A., Tabbara, K. F., Al-Rajhi, A. A., et al. (1998). An immunohistochemical study of collagens in trachoma and vernal keratoconjunctivitis. Eye 12(Pt 6), 1001–1006. doi: 10.1038/eye.1998.257
Alipour, M., Zaidi, D., Valcheva, R., Jovel, J., Martinez, I., Sergi, C., et al. (2016). Mucosal barrier depletion and loss of bacterial diversity are primary abnormalities in paediatric ulcerative colitis. J. Crohns. Colitis 10, 462–471. doi: 10.1093/ecco-jcc/jjv223
Beatty, W. L., Belanger, T. A., Desai, A. A., Morrison, R. P., and Byrne, G. I. (1994). Tryptophan depletion as a mechanism of gamma interferon-mediated chlamydial persistence. Infect. Immun. 62, 3705–3711.
Bobo, L., Novak, N., Mkocha, H., Vitale, S., West, S., and Quinn, T. C. (1996). Evidence for a predominant proinflammatory conjunctival cytokine response in individuals with trachoma. Infect. Immun. 64, 3273–3279.
Bourne, R. R., Stevens, G. A., White, R. A., Smith, J. L., Flaxman, S. R., Price, H., et al. (2013). Causes of vision loss worldwide, 1990-2010: a systematic analysis. Lancet Glob. Health 1, E339–E349. doi: 10.1016/S2214-109X(13)70113-X
Burton, M. J., Bailey, R. L., Jeffries, D., Mabey, D. C., and Holland, M. J. (2004). Cytokine and fibrogenic gene expression in the conjunctivas of subjects from a gambian community where trachoma is endemic. Infect. Immun. 72, 7352–7356. doi: 10.1128/IAI.72.12.7352-7356.2004
Burton, M. J., Holland, M. J., Faal, N., Aryee, E. A., Alexander, N. D., Bah, M., et al. (2003). Which members of a community need antibiotics to control trachoma? conjunctival Chlamydia trachomatis infection load in gambian villages. Invest. Ophthalmol. Vis. Sci. 44, 4215–4222. doi: 10.1167/iovs.03-0107
Burton, M. J., Hu, V. H., Massae, P., Burr, S. E., Chevallier, C., Afwamba, I. A., et al. (2011a). What is causing active trachoma? the role of non-chlamydial bacterial pathogens in a low prevalence setting. Invest. Ophthalmol. Vis. Sci. 52, 6012–6017. doi: 10.1167/iovs.11-7326
Burton, M. J., Rajak, S. N., Hu, V. H., Ramadhani, A., Habtamu, E., Massae, P., et al. (2015). Pathogenesis of progressive scarring trachoma in ethiopia and tanzania and its implications for disease control: two cohort studies. PLoS Negl. Trop. Dis. 9:E0003763. doi: 10.1371/journal.pntd.0003763
Burton, M. J., Rajak, S. N., Ramadhani, A., Weiss, H. A., Habtamu, E., Abera, B., et al. (2012). Post-operative recurrent trachomatous trichiasis is associated with increased conjunctival expression of S100a7 (Psoriasin). PLoS Negl. Trop. Dis. 6:E1985. doi: 10.1371/journal.pntd.0001985
Burton, M. J., Ramadhani, A., Weiss, H. A., Hu, V., Massae, P., Burr, S. E., et al. (2011c). Active trachoma is associated with increased conjunctival expression of Il17a and profibrotic cytokines. Infect. Immun. 79, 4977–4983. doi: 10.1128/IAI.05718-11
Dawson, C. R., Marx, R., Daghfous, T., Juster, R., and Schachter, J. (1990). What clinical signs are critical in evaluating the intervention in trachoma? in Chlamydial Infections, ed W. R. Bowie (Cambridge: Cambridge University Press), 271–278.
Derrick, T., Last, A. R., Burr, S. E., Roberts, C. H., Nabicassa, M., Cassama, E., et al. (2016a). Inverse relationship between microrna-155 And -184 expression with increasing conjunctival inflammation during ocular Chlamydia trachomatis infection. BMC Infect. Dis. 16, 60. doi: 10.1186/s12879-016-1367-8
Derrick, T., Luthert, P. J., Jama, H., Hu, V. H., Massae, P., Essex, D., et al. (2016b). Increased epithelial expression of ctgf and S100a7 with elevated subepithelial expression Of Il-1beta in trachomatous trichiasis. PLoS Negl. Trop. Dis. 10:E0004752. doi: 10.1371/journal.pntd.0004752
Derrick, T., Roberts, C., Rajasekhar, M., Burr, S. E., Joof, H., Makalo, P., et al. (2013). Conjunctival microrna expression in inflammatory trachomatous scarring. PLoS Negl. Trop. Dis. 7:E2117. doi: 10.1371/journal.pntd.0002117
Faal, N., Bailey, R. L., Jeffries, D., Joof, H., Sarr, I., Laye, M., et al. (2006). Conjunctival Foxp3 expression in trachoma: do regulatory t cells have a role in human ocular Chlamydia trachomatis infection? PLoS Med. 3:E266. doi: 10.1371/journal.pmed.0030266
Fini, M. E., Parks, W. C., Rinehart, W. B., Girard, M. T., Matsubara, M., Cook, J. R., et al. (1996). Role of matrix metalloproteinases in failure to re-epithelialize after corneal injury. Am. J. Pathol. 149, 1287–1302.
Fischer, S. F., Vier, J., Kirschnek, S., Klos, A., Hess, S., Ying, S., et al. (2004). Chlamydia inhibit host cell apoptosis by degradation of proapoptotic Bh3-only proteins. J. Exp. Med. 200, 905–916. doi: 10.1084/jem.20040402
Frazer, L. C., O'connell, C. M., Andrews, C. W. Jr., Zurenski, M. A., and Darville, T. (2011). Enhanced neutrophil longevity and recruitment contribute to the severity of oviduct pathology during Chlamydia muridarum infection. Infect. Immun. 79, 4029–4041. doi: 10.1128/IAI.05535-11
Frazer, L. C., Scurlock, A. M., Zurenski, M. A., Riley, M. M., Mintus, M., Pociask, D. A., et al. (2013). Il-23 Induces Il-22 and Il-17 production in response to Chlamydia muridarum genital tract infection, but the absence of these cytokines does not influence disease pathogenesis. Am. J. Reprod. Immunol. 70, 472–484. doi: 10.1111/aji.12171
Gall, A., Horowitz, A., Joof, H., Natividad, A., Tetteh, K., Riley, E., et al. (2011). Systemic effector and regulatory immune responses to chlamydial antigens in trachomatous trichiasis. Front. Microbiol. 2:10. doi: 10.3389/fmicb.2011.00010
Girard, J. P., and Springer, T. A. (1996). Modulation of endothelial cell adhesion by hevin, an acidic protein associated with high endothelial venules. J. Biol. Chem. 271, 4511–4517. doi: 10.1074/jbc.271.8.4511
Grassly, N. C., Ward, M. E., Ferris, S., Mabey, D. C., and Bailey, R. L. (2008). The natural history of trachoma infection and disease in a gambian cohort with frequent follow-up. PLoS Negl. Trop. Dis. 2:E341. doi: 10.1371/journal.pntd.0000341
Guilhelmelli, F., Vilela, N., Albuquerque, P., Derengowski Lda, S., Silva-Pereira, I., and Kyaw, C. M. (2013). Antibiotic development challenges: the various mechanisms of action of antimicrobial peptides and of bacterial resistance. Front. Microbiol. 4:353. doi: 10.3389/fmicb.2013.00353
Hewagama, A., Patel, D., Yarlagadda, S., Strickland, F. M., and Richardson, B. C. (2009). Stronger inflammatory/cytotoxic t-cell response in women identified by microarray analysis. Genes Immun. 10, 509–516. doi: 10.1038/gene.2009.12
Holland, M. J., Faal, N., Sarr, I., Joof, H., Laye, M., Cameron, E., et al. (2006). The Frequency Of Chlamydia trachomatis major outer membrane protein-specific Cd8+ T lymphocytes in active trachoma is associated with current ocular infection. Infect. Immun. 74, 1565–1572. doi: 10.1128/IAI.74.3.1565-1572.2006
Holland, M. J., Jeffries, D., Pattison, M., Korr, G., Gall, A., Joof, H., et al. (2010). Pathway-focused arrays reveal increased matrix metalloproteinase-7 (matrilysin) transcription in trachomatous trichiasis. Invest. Ophthalmol. Vis. Sci. 51, 3893–3902. doi: 10.1167/iovs.09-5054
Hu, V. H., Holland, M. J., and Burton, M. J. (2013a). Trachoma: protective and pathogenic ocular immune responses to Chlamydia trachomatis. PLoS Negl. Trop. Dis. 7:E2020. doi: 10.1371/journal.pntd.0002020
Hu, V. H., Holland, M. J., Cree, I. A., Pullin, J., Weiss, H. A., Massae, P., et al. (2013b). In vivo confocal microscopy and histopathology of the conjunctiva in trachomatous scarring and normal tissue: a systematic comparison. Br. J. Ophthalmol. 97, 1333–1337. doi: 10.1136/bjophthalmol-2013-303126
Hu, V. H., Luthert, P. J., Derrick, T., Pullin, J., Weiss, H. A., Massae, P., et al. (2016). Immunohistochemical analysis of scarring trachoma indicates infiltration by natural killer and undefined Cd45 negative cells. PLoS Negl. Trop. Dis. 10:E0004734. doi: 10.1371/journal.pntd.0004734
Hu, V. H., Massae, P., Weiss, H. A., Chevallier, C., Onyango, J. J., Afwamba, I. A., et al. (2011). Bacterial infection in scarring trachoma. Invest. Ophthalmol. Vis. Sci. 52, 2181–2186. doi: 10.1167/iovs.10-5829
Hu, V. H., Weiss, H. A., Ramadhani, A. M., Tolbert, S. B., Massae, P., Mabey, D. C., et al. (2012). Innate immune responses and modified extracellular matrix regulation characterize bacterial infection and cellular/connective tissue changes in scarring trachoma. Infect. Immun. 80, 121–130. doi: 10.1128/IAI.05965-11
Johansson, M. E., Gustafsson, J. K., Holmen-Larsson, J., Jabbar, K. S., Xia, L., Xu, H., et al. (2014). Bacteria penetrate the normally impenetrable inner colon mucus layer in both murine colitis models and patients with ulcerative colitis. Gut 63, 281–291. doi: 10.1136/gutjnl-2012-303207
Kechagia, J. Z., Ezra, D. G., Burton, M. J., and Bailly, M. (2016). Fibroblasts profiling in scarring trachoma identifies Il-6 as a functional component of a fibroblast-macrophage pro-fibrotic and pro-inflammatory feedback loop. Sci. Rep. 6:28261. doi: 10.1038/srep28261
Kon, C. H., Occleston, N. L., Charteris, D., Daniels, J., Aylward, G. W., and Khaw, P. T. (1998). A prospective study of matrix metalloproteinases in proliferative vitreoretinopathy. Invest. Ophthalmol. Vis. Sci. 39, 1524–1529.
Lacy, H. M., Bowlin, A. K., Hennings, L., Scurlock, A. M., Nagarajan, U. M., and Rank, R. G. (2011). Essential role for neutrophils in pathogenesis and adaptive immunity in Chlamydia caviae ocular infections. Infect. Immun. 79, 1889–1897. doi: 10.1128/IAI.01257-10
Lad, S. P., Fukuda, E. Y., Li, J., De La Maza, L. M., and Li, E. (2005). Up-regulation of the jak/stat1 signal pathway during Chlamydia trachomatis infection. J. Immunol. 174, 7186–7193. doi: 10.4049/jimmunol.174.11.7186
Last, A. R., Roberts, C., Cassama, E., Nabicassa, M., Molina-Gonzalez, S., Burr, S. E., et al. (2014). Plasmid copy number and disease severity in naturally occurring ocular Chlamydia trachomatis infection. J. Clin. Microbiol. 52, 324–327. doi: 10.1128/JCM.02618-13
Lee, J. S., Munoz, B. E., Mkocha, H., Gaydos, C. A., Quinn, T. C., and West, S. K. (2014). The effect of multiple rounds of mass drug administration on the association between ocular Chlamydia trachomatis infection and follicular trachoma in preschool-aged children. PLoS Negl. Trop. Dis. 8:E2761. doi: 10.1371/journal.pntd.0002761
Lu, C., Zeng, H., Li, Z., Lei, L., Yeh, I. T., Wu, Y., et al. (2012). Protective immunity against mouse upper genital tract pathology correlates with high ifngamma but low Il-17 T cell and anti-secretion protein antibody responses induced by replicating chlamydial organisms in the airway. Vaccine 30, 475–485. doi: 10.1016/j.vaccine.2011.10.059
Manel, N., Unutmaz, D., and Littman, D. R. (2008). The differentiation of human TH-17 cells requires transforming growth factor-beta and induction of the nuclear receptor rorgammat. Nat. Immunol. 9, 641–649. doi: 10.1038/ni.1610
Michel, C. E., Roper, K. G., Divena, M. A., Lee, H. H., and Taylor, H. R. (2011). Correlation of clinical trachoma and infection in aboriginal communities. PLoS Negl. Trop. Dis. 5:E986. doi: 10.1371/journal.pntd.0000986
Mikhak, Z., Fleming, C. M., Medoff, B. D., Thomas, S. Y., Tager, A. M., Campanella, G. S., et al. (2006). Stat1 in peripheral tissue differentially regulates homing of antigen-specific Th1 and Th2 cells. J. Immunol. 176, 4959–4967. doi: 10.4049/jimmunol.176.8.4959
Misko, A., Ferguson, T., and Notterpek, L. (2002). Matrix Metalloproteinase mediated degradation of basement membrane proteins in trembler J neuropathy nerves. J. Neurochem. 83, 885–894. doi: 10.1046/j.1471-4159.2002.01200.x
Mozzato-Chamay, N., Mahdi, O. S., Jallow, O., Mabey, D. C., Bailey, R. L., and Conway, D. J. (2000). Polymorphisms in candidate genes and risk of scarring trachoma in a Chlamydia trachomatis–endemic population. J. Infect. Dis. 182, 1545–1548. doi: 10.1086/315891
Natividad, A., Freeman, T. C., Jeffries, D., Burton, M. J., Mabey, D. C., Bailey, R. L., et al. (2010). Human conjunctival transcriptome analysis reveals the prominence of innate defense in Chlamydia trachomatis infection. Infect. Immun. 78, 4895–4911. doi: 10.1128/IAI.00844-10
Natividad, A., Wilson, J., Koch, O., Holland, M. J., Rockett, K., Faal, N., et al. (2005). Risk of trachomatous scarring and trichiasis in gambians varies with snp haplotypes at the interferon-gamma and interleukin-10 loci. Genes Immun. 6, 332–340. doi: 10.1038/sj.gene.6364182
Nelson, P. S., Plymate, S. R., Wang, K., True, L. D., Ware, J. L., Gan, L., et al. (1998). Hevin, an antiadhesive extracellular matrix protein, is down-regulated in metastatic prostate adenocarcinoma. Cancer Res. 58, 232–236.
Nishikomori, R., Usui, T., Wu, C. Y., Morinobu, A., O'shea, J. J., and Strober, W. (2002). Activated Stat4 has an essential role in Th1 differentiation and proliferation that is independent of its role in the maintenance of Il-12r beta 2 chain expression and signaling. J. Immunol. 169, 4388–4398. doi: 10.4049/jimmunol.169.8.4388
Nwomeh, B. C., Liang, H. X., Diegelmann, R. F., Cohen, I. K., and Yager, D. R. (1998). Dynamics of the matrix metalloproteinases mmp-1 and Mmp-8 in acute open human dermal wounds. Wound Repair Regen. 6, 127–134. doi: 10.1046/j.1524-475X.1998.60206.x
Olivares-Zavaleta, N., Whitmire, W. M., Kari, L., Sturdevant, G. L., and Caldwell, H. D. (2014). Cd8+ T cells define an unexpected role in live-attenuated vaccine protective immunity against Chlamydia trachomatis infection in macaques. J. Immunol. 192, 4648–4654. doi: 10.4049/jimmunol.1400120
Ozdemir, E., Kakehi, Y., Okuno, H., and Yoshida, O. (1999). Role of matrix metalloproteinase-9 in the basement membrane destruction of superficial urothelial carcinomas. J. Urol. 161, 1359–1363. doi: 10.1016/S0022-5347(01)61684-7
Pickett, M. A., Everson, J. S., Pead, P. J., and Clarke, I. N. (2005). The plasmids of Chlamydia trachomatis and chlamydophila pneumoniae (N16): accurate determination of copy number and the paradoxical effect of plasmid-curing agents. Microbiology 151, 893–903. doi: 10.1099/mic.0.27625-0
Radaeva, S., Sun, R., Pan, H. N., Hong, F., and Gao, B. (2004). Interleukin 22 (Il-22) plays a protective role in T cell-mediated murine hepatitis: Il-22 is a survival factor for hepatocytes via stat3 activation. Hepatology 39, 1332–1342. doi: 10.1002/hep.20184
Ramadhani, A. M., Derrick, T., Holland, M. J., and Burton, M. J. (2016a). Blinding trachoma: systematic review of rates and risk factors for progressive disease. PLoS Negl. Trop. Dis. 10:E0004859. doi: 10.1371/journal.pntd.0004859
Ramadhani, A. M., Derrick, T., Macleod, D., Holland, M. J., and Burton, M. J. (2016b). The relationship between active trachoma and ocular Chlamydia trachomatis infection before and after mass antibiotic treatment. PLoS Negl. Trop. Dis. 10:E0005080. doi: 10.1371/journal.pntd.0005080
Rank, R. G., Bowlin, A. K., and Kelly, K. A. (2000). Characterization of lymphocyte response in the female genital tract during ascending chlamydial genital infection in the guinea pig model. Infect. Immun. 68, 5293–5298. doi: 10.1128/IAI.68.9.5293-5298.2000
Redgrove, K. A., and McLaughlin, E. A. (2014). The role of the immune response in Chlamydia trachomatis infection of the male genital tract: a double-edged sword. Front. Immunol. 5:534. doi: 10.3389/fimmu.2014.00534
Roberts, C. H., Franklin, C. S., Makalo, P., Joof, H., Sarr, I., Mahdi, O. S., et al. (2015). Conjunctival fibrosis and the innate barriers to Chlamydia trachomatis intracellular infection: a genome wide association study. Sci. Rep. 5:17447. doi: 10.1038/srep17447
Roberts, C. H., Last, A., Molina-Gonzalez, S., Cassama, E., Butcher, R., Nabicassa, M., et al. (2013). Development and evaluation of a next-generation digital pcr diagnostic assay for ocular Chlamydia trachomatis infections. J. Clin. Microbiol. 51, 2195–2203. doi: 10.1128/JCM.00622-13
Sankaran-Walters, S., Macal, M., Grishina, I., Nagy, L., Goulart, L., Coolidge, K., et al. (2013). Sex differences matter in the gut: effect on mucosal immune activation and inflammation. Biol. Sex Differ. 4:10. doi: 10.1186/2042-6410-4-10
Shapiro, S. D., Kobayashi, D. K., and Ley, T. J. (1993). Cloning and characterization of a unique elastolytic metalloproteinase produced by human alveolar macrophages. J. Biol. Chem. 268, 23824–23829.
Siegl, C., Prusty, B. K., Karunakaran, K., Wischhusen, J., and Rudel, T. (2014). Tumor Suppressor P53 alters host cell metabolism to limit Chlamydia trachomatis infection. Cell Rep. 9, 918–929. doi: 10.1016/j.celrep.2014.10.004
Skwor, T. A., Atik, B., Kandel, R. P., Adhikari, H. K., Sharma, B., and Dean, D. (2008). Role of secreted conjunctival mucosal cytokine and chemokine proteins in different stages of trachomatous disease. PLoS Negl. Trop. Dis. 2:E264. doi: 10.1371/journal.pntd.0000264
Solomon, A. W., Holland, M. J., Burton, M. J., West, S. K., Alexander, N. D., Aguirre, A., et al. (2003). Strategies for control of trachoma: observational study with quantitative Pcr. Lancet 362, 198–204. doi: 10.1016/S0140-6736(03)13909-8
Sullivan, M. M., and Sage, E. H. (2004). Hevin/Sc1, A matricellular glycoprotein and potential tumor-suppressor of the Sparc/Bm-40/osteonectin family. Int. J. Biochem. Cell Biol. 36, 991–996. doi: 10.1016/j.biocel.2004.01.017
Turner, J. E., Stockinger, B., and Helmby, H. (2013). Il-22 Mediates Goblet cell hyperplasia and worm expulsion in intestinal helminth infection. PLoS Pathog. 9:E1003698. doi: 10.1371/journal.ppat.1003698
Van Voorhis, W. C., Barrett, L. K., Sweeney, Y. T., Kuo, C. C., and Patton, D. L. (1997). Repeated Chlamydia trachomatis infection of Macaca nemestrina fallopian tubes produces a Th1-like cytokine response associated with fibrosis and scarring. Infect. Immun. 65, 2175–2182.
Vicetti Miguel, R. D., Quispe Calla, N. E., Pavelko, S. D., and Cherpes, T. L. (2016). Intravaginal Chlamydia trachomatis challenge infection elicits Th1 and Th17 immune responses in mice that promote pathogen clearance and genital tract damage. PLoS ONE 11:E0162445. doi: 10.1371/journal.pone.0162445
West, S. K., Munoz, B., Mkocha, H., Hsieh, Y. H., and Lynch, M. C. (2001). Progression of active trachoma to scarring in a cohort of tanzanian children. Ophthalmic Epidemiol. 8, 137–144. doi: 10.1076/opep.18.104.22.16858
Wolf, M., Maltseva, I., Clay, S. M., Pan, P., Gajjala, A., and Chan, M. F. (2017). Effects of Mmp12 on cell motility and inflammation during corneal epithelial repair. Exp. Eye Res. 160, 11–20. doi: 10.1016/j.exer.2017.04.007
Wolk, K., Witte, E., Wallace, E., Docke, W. D., Kunz, S., Asadullah, K., et al. (2006). Il-22 regulates the expression of genes responsible for antimicrobial defense, cellular differentiation, and mobility in keratinocytes: a potential role in psoriasis. Eur. J. Immunol. 36, 1309–1323. doi: 10.1002/eji.200535503
Wolle, M. A., Munoz, B. E., Mkocha, H., and West, S. K. (2009b). Constant ocular infection with Chlamydia trachomatis predicts risk of scarring in children in tanzania. Ophthalmology 116, 243–247. doi: 10.1016/j.ophtha.2008.09.011
Wolle, M. A., Munoz, B., Mkocha, H., and West, S. K. (2009a). Age, sex, and cohort effects in a longitudinal study of trachomatous scarring. Invest. Ophthalmol. Vis. Sci. 50, 592–596. doi: 10.1167/iovs.08-2414
Yang, X. P., Ghoreschi, K., Steward-Tharp, S. M., Rodriguez-Canales, J., Zhu, J., Grainger, J. R., et al. (2011). Opposing regulation of the locus encoding Il-17 through direct, reciprocal actions of stat3 and stat5. Nat. Immunol. 12, 247–254. doi: 10.1038/ni.1995
Zeng, Z. S., Cohen, A. M., and Guillem, J. G. (1999). Loss of basement membrane type iv collagen is associated with increased expression of metalloproteinases 2 and 9 (Mmp-2 and Mmp-9) during human colorectal tumorigenesis. Carcinogenesis 20, 749–755. doi: 10.1093/carcin/20.5.749
Zhang, H., Liu, D., Sun, S., Liu, J., Wu, S., Mao, Y., et al. (2009). Selective expression of calcium-binding protein S100a7 in lung cancer. Zhongguo. Fei. Ai. Za. Zhi. 12, 54–58. doi: 10.3779/j.issn.1009-3419.2009.01.009
Zhou, Y., Holland, M. J., Makalo, P., Joof, H., Roberts, C. H., Mabey, D. C., et al. (2014). The conjunctival microbiome in health and trachomatous disease: a case control study. Genome Med. 6:99. doi: 10.1186/s13073-014-0099-x
Keywords: trachoma, Chlamydia trachomatis, gene expression, mass azithromycin administration, longitudinal study, active trachoma, conjunctival scarring
Citation: Ramadhani AM, Derrick T, Macleod D, Massae P, Mtuy T, Jeffries D, Roberts CH, Bailey RL, Mabey DCW, Holland MJ and Burton MJ (2017) Immunofibrogenic Gene Expression Patterns in Tanzanian Children with Ocular Chlamydia trachomatis Infection, Active Trachoma and Scarring: Baseline Results of a 4-Year Longitudinal Study. Front. Cell. Infect. Microbiol. 7:406. doi: 10.3389/fcimb.2017.00406
Received: 28 June 2017; Accepted: 31 August 2017;
Published: 15 September 2017.
Edited by:Rey Carabeo, Washington State University, United States
Reviewed by:Ashlesh Murthy, Midwestern University, United States
Xiaojing Zheng, University of North Carolina at Chapel Hill, United States
Peter Timms, University of the Sunshine Coast, Australia
Copyright © 2017 Ramadhani, Derrick, Macleod, Massae, Mtuy, Jeffries, Roberts, Bailey, Mabey, Holland and Burton. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Matthew J. Burton, firstname.lastname@example.org