ORIGINAL RESEARCH article
Active Trachoma Cases in the Solomon Islands Have Varied Polymicrobial Community Structures but Do Not Associate with Individual Non-Chlamydial Pathogens of the Eye
- 1Clinical Research Department, London School of Hygiene & Tropical Medicine, London, United Kingdom
- 2Eye Department, Ministry of Health and Medical Services, Honiara, Solomon Islands
- 3Primary Care Department, Lata Hospital, Santa Cruz Island, Solomon Islands
- 4Bellona Rural Health Centre, Bellona, Solomon Islands
- 5The Fred Hollows Foundation, Carlton, VIC, Australia
Background: Several non-chlamydial microbial pathogens are associated with clinical signs of active trachoma in trachoma-endemic communities with a low prevalence of ocular Chlamydia trachomatis (Ct) infection. In the Solomon Islands, the prevalence of Ct among children is low despite the prevalence of active trachoma being moderate. Therefore, we set out to investigate whether active trachoma was associated with a common non-chlamydial infection or with a dominant polymicrobial community dysbiosis in the Solomon Islands.
Methods: We studied DNA from conjunctival swabs collected from 257 Solomon Islanders with active trachoma and matched controls. Droplet digital PCR was used to test for pathogens suspected to be able to induce follicular conjunctivitis. Polymicrobial community diversity and composition were studied by sequencing of hypervariable regions of the 16S ribosomal ribonucleic acid gene in a subset of 54 cases and 53 controls.
Results: Although Ct was associated with active trachoma, the number of infections was low (cases, 3.9%; controls, 0.4%). Estimated prevalence (cases and controls, respectively) of each non-chlamydial infection was as follows: Staphylococcus aureus: 1.9 and 1.9%, Adenoviridae: 1.2 and 1.2%, coagulase-negative Staphylococcus: 5.8 and 4.3%, Haemophilus influenzae: 7.4 and 11.7%, Moraxella catarrhalis: 2.3 and 4.7%, and Streptococcus pneumoniae: 7.0 and 6.2%. There was no statistically significant association between the clinical signs of trachoma and the presence or load of any of the non-Ct infections that were assayed. Interindividual variations in the conjunctival microbiome were characterized by differences in the levels of Corynebacterium, Propionibacterium, Helicobacter, and Paracoccus, but diversity and relative abundance of these specific genera did not differ significantly between cases and controls.
Discussion: It is unlikely that the prevalent trachoma-like follicular conjunctivitis in this region of the Solomon Islands has a dominant bacterial etiology. Before implementing community-wide azithromycin distribution for trachoma, policy makers should consider that clinical signs of trachoma can be observed in the absence of any detectable azithromycin-susceptible organism.
Trachoma, caused by Chlamydia trachomatis (Ct), is the leading infectious cause of preventable blindness (1) and is targeted for elimination as a public health problem by 2020 through the SAFE strategy (Surgery, Antibiotics, Facial cleanliness and Environmental improvement). The decision to implement community-wide trachoma control interventions, which include mass drug administration (MDA) with azithromycin, is based on population prevalence estimates of one clinical sign of active trachoma [trachomatous inflammation–follicular (TF)] in the 1- to 9-year-old age group (2). In the Solomon Islands, the prevalence of TF is sufficient to warrant MDA, but the prevalence of trachomatous trichiasis (TT) suggests that trachoma may not pose a significant public health problem. A recent trachoma survey in a treatment-naive population of the Solomon Islands estimated 26% of 1 to 9 year olds to have TF, but, surprisingly, conjunctival Ct infection was detected in only 1.3% of that age group (3). While Ct infection is not always detectable in TF cases (4, 5), infection prevalence in the Solomon Islands is far lower than is seen in other countries with similar TF prevalence (3).
Chlamydia trachomatis is not the only pathogen to associate with signs of TF. A number of differential diagnoses for follicular conjunctivitis are described (6), some of which are rare and distinguishable by patient history (such as Parinaud’s oculoglandular syndrome). Other causes of conjunctivitis, which can present with follicles, are Streptococcus pneumoniae, Haemophilus influenzae (7), and adenovirus (8). Non-ocular Chlamydia serotypes (9) and non-Ct species (10) have also been suggested as possible causes of TF. A number of studies have examined conjunctival microbiology in trachoma-endemic areas. A greater diversity of pathogens can be cultured from the conjunctivae of people with TF compared to counterparts without TF (7). Among bacterial species isolated during a study in Tanzanian children, S. pneumoniae, H. influenzae B, and H. influenzae non-type B associated more strongly than Ct with clinical signs of active trachoma (7). In The Gambia, S. pneumoniae and H. influenzae type B infection correlated closely with signs of active trachoma in communities that had previously received MDA, whereas Moraxella catarrhalis and Staphylococcus aureus did not (11). Common non-chlamydial pathogens have associated with trachomatous scarring (TS) (12) and recurrence of TT after surgery (13) in some, but not all (13, 14), studies.
Although traditionally thought to be “sterile”, several studies have now described polymicrobial communities colonizing the conjunctival epithelium (15, 16). Staphylococcus, Streptococcus, Haemophilus, and Moraxella genera can be readily cultured from swabs taken from the inferior fornix (7, 11, 12). Attempts to detect bacteria from the conjunctiva have shown community diversity and composition to vary significantly between individuals. It is unclear whether a “core” microbiota persists at the conjunctiva, but profiles closely related to that of the skin have been reported (17). Corynebacterium, Propionibacterium, Staphylococcus, and Streptococcus have been consistently dominant, whereas other genera such as Acinetobacter, Brevundimonas, Pseudomonas, Bradyrhizobium, Sphingomonas, Bacillus, Simonsiella, and Elizabethkingia have been identified more sporadically (18–20). Conjunctival polymicrobial community composition is known to vary with age and season (21) and appears to be responsive to external stimuli such as regular contact lens wear (22). The conjunctival microbiome varies significantly between people with TS and those without; however, significant associations between microbiome and follicular conjunctival inflammation or microbiome and other diseases have yet to be described (21, 23).
We hypothesized that clinical signs of active trachoma in the Solomon Islands where ocular Ct is uncommon could be explained by a common non-chlamydial infection or by a dominant polymicrobial community dysbiosis.
Materials and Methods
A parent or guardian gave written informed consent for each child to take part in the survey in accordance with the Declaration of Helsinki. The protocol was approved by the London School of Hygiene & Tropical Medicine (6319/6360) and the Solomon Islands National Health Research (HRC13/18) Ethics Committees.
The samples tested during this study were a subset of specimens from a population-based trachoma prevalence survey of 3,674 people (1,135 children aged 1–9 years) in 32 clusters from Temotu and Rennell & Bellona Provinces, Solomon Islands; data from that survey are presented elsewhere (3). At the time of the survey, no trachoma interventions had been implemented in the Solomon Islands. Sterile swabs were passed three times over the right conjunctiva of every child before storage in 300 µL RNALater® (Life Technologies). Swabs were kept cool in the field for up to 48 h then frozen until processing (3). Inclusion criteria for cases (n = 257) were age 1–9 years, detectable human DNA in the conjunctival swab, and a field grade of TF or trachomatous inflammation—intense (TI) in the right eye. An equal number (n = 257) of age-, gender- and island-matched controls without TF or TI were selected using the “e1071” R package. A random subset of cases (n = 54) and controls (n = 53) were further characterized using V1–V3 16S ribosomal RNA (rRNA) gene sequencing. Field controls (clean swabs passed within 20 cm of seated participants, without touching them, and then treated identically to subjects’ specimens), extraction controls (extraction carried out in the absence of a specimen and PCR carried out on the eluate), and PCR controls (PCR carried out in the absence of any added material) were used to determine background levels of microbial contamination in the 16S rRNA gene sequencing protocol.
Droplet Digital PCR
DNA was extracted from conjunctival swabs using the Qiagen AllPrep DNA/RNA kit. We chose not to use mechanical lysis on these low biomass specimens due to the reported lower yields compared to chemical lysis (24). All specimens have previously been tested for Ct and Homo sapiens ribonuclease P protein subunit p30 (acting as an endogenous control) using a validated assay (25); the community prevalence of ocular Ct infection is described elsewhere (3).
Duplex ddPCR assays were developed for (1) Adenoviridae (26) and S. pneumoniae (2, 27), H. influenzae (28) and M. catarrhalis (29, 30), and S. aureus and coagulase-negative Staphylococcus (31, 32) based on the published assays. Each assay contained 10 µL 2× ddPCR Supermix for probes (Bio-rad, Hemel Hempstead, UK), primers and probes at custom concentrations (Table 1), and 8 µL of template DNA. Thermal cycling was performed for 10 mins at 95°C, 40 cycles of 15 seconds at 95°C then 1 min at 60°C, followed by 12 minutes at 98°C. A positive result for a clinical specimen was defined as >95% confidence of non-zero target load, as described previously (25). Assay performance was assessed by repeat testing of PCR product dilution series’ and cultured pathogen material before rolling out to clinical samples. The reproducibility, linearity, and limits of detection of each assay were considered appropriate (Table S1 in Supplementary Material) and were also similar to published performance of the Ct assay used (25).
Table 1. Oligonucleotides used in this study, with assay concentrations derived from in vitro optimization.
16S rRNA Gene Sequencing
An approximately 530-bp region of the 16S rRNA gene (variable regions 1–3) was amplified using forward (modified 27F; 5′-[adaptor]-AGAGTTTGGATCCTGGCTCAG-3′) and custom barcoded reverse primers (534R; 5′-[adaptor]-[barcode]-AGTCAGTCAGCCATTACCGCGGCTGCTGG-3′). Each 15 µL reaction contained 7.5 µL 2× Phusion High Fidelity Master Mix (New England Biosciences, MA, USA), 0.45 µL DMSO, 0.1 µM primers, and 5.55 µL DNA. The thermal cycling conditions were as follows: 30 seconds at 98°C followed by 31 cycles of 10 seconds at 98°C then 30 seconds at 62°C followed by 7 minutes at 72°C. Amplicons were cleaned using 0.6 v/v AMPure XP beads (Beckman Coulter, CA, USA), quantified using Qubit (Thermo Fisher Scientific, MA, USA), and then pooled at equimolar concentrations. The 4-nM sequencing library was mixed 0.75 v/v with a 4 nM Phi-X control library (Illumina, CA, USA). 3 µL of 100 µM custom read primers were mixed into wells 12 (Read 1; CTACACTATGGTAATTGTAGAGTTTGGATCCTGGCTCAG), 13 (Index; CCAGCAGCCGCGGTAATGGCTGACTGACT), and 14 (Read 2; AGTCAGTCAGCCATTACCGCGGCTGCTGG) of the MiSeq reagents cartridge before 2 × 300 bp paired-end sequencing with the 600-cycle MiSeq v3 sequencing reagent kit on the MiSeq platform using a standard Illumina protocol. 96 uniquely barcoded specimens were run in multiplex on each MiSeq run.
ddPCR data were analyzed using R version 3.2.2 (34). Binomial univariate regression was used to test the relationship between each individual infection and active trachoma. The most accurate final multivariate model, determined by the lowest Akaike Information Criterion value, was determined by stepwise removal of variables from a binomial multivariate regression model incorporating all tested pathogens. The chance of the Ct association occurring due to chance was assessed by counting the number of significant results obtained from randomly reordering the Ct data 1,000 times.
Raw 16S amplicon sequences were directly assigned to genera using Illumina BaseSpace “16S Metagenomics” app version 126.96.36.199, which uses algorithms from Wang et al. (35). Two clinical specimens yielded fewer than 1,000 reads in total and were removed from the analysis. Amplicons were sequenced from six no-template control (NTC) specimens (two “field” controls, two “extraction” controls, and two “PCR” controls, defined above) to identify any contaminants endogenous to the collection, PCR, or sequencing process. Any genera represented by more than 600 reads in the 6 combined NTCs (average 100 reads per NTC) were eliminated from the clinical specimen analysis. This resulted in the removal of 21 genera, which accounted for 95.1% of the reads from NTCs and 88.3% of the reads from clinical specimens. The genera removed included common contaminants of reagent kits (36) such as Pelomonas, as well as some previously described conjunctival microbiota constituents such as Brevundimonas, Staphylococcus, and Streptococcus (18, 21).
Between-group differences in the Shannon Diversity Index (a measure of diversity, which increases with increasing species abundance and evenness) and Inverse Simpson Index (another measure of diversity where 1 is infinite diversity and 0 is no diversity) were compared using a t-test. Discriminant analysis of principal components, a multivariate clustering method to analyze highly dimensional datasets, was performed using the “adegenet” package in R (37), and cross-validation was used to determine the optimal number of principal components (PCs) to include in a discriminant function aimed at separating cases from controls on the basis of their polymicrobial community structures.
Specimen Set Demographics
A total of 257 cases and 257 controls (n = 514) were tested. All 257 cases had TF, and one also had TI. Case status was defined by clinical signs in the right eye, but 236/257 (91.8%) and 8/257 (3.1%) of the cases and controls, respectively, had TF in the left eye. Both case and control groups were 38% female. Mean age was 5.6 years (cases) and 5.5 years (controls; Student’s t-test p = 0.76). There was no significant difference between cases and controls in terms of the clusters represented (Kolmogorov–Smirnov test p = 0.97), and all 32 clusters were represented. The case specimens had higher loads of human DNA than the controls (18,030 versus 9,354 copies/swab; p = 0.00003). There were no significant differences in age, gender, or location within the subset selected for 16S-amplicon sequencing.
Quantitative PCR Tests for Ocular Pathogens
In this study, 19.8% of children had evidence of infection with at least one of the targeted organisms (Table 2). We considered the prevalence of Ct in this sample set to be too low (96.1% of TF cases were Ct infection negative) to account for the level of active disease. The prevalence of both Adenoviridae (1.2%) and S. aureus (1.9%) were very low. There was no association between active trachoma and infection with H. influenzae, S. pneumoniae, S. aureus, coagulase-negative Staphylococcus, or M. catarrhalis. Ct infection was associated with active trachoma (logistic regression p = 0.026; odds ratio, 10.4). The association between Ct and active trachoma was still significant in a multivariate analysis (p = 0.025). The permuted p value for the association between Ct and active trachoma was highly significant (p = 0.001). Active disease was neither associated with infection with at least one pathogen (any pathogen, p = 0.185) nor the number of concurrent infections in the same eye (number of concurrent pathogens, p = 0.207). Infection loads of those who tested positive are shown in Figure 1. Despite some numerical differences between groups, none of the differences were statistically significant.
Figure 1. Target copies per swab of each pathogen identified in conjunctival swabs collected from children with and without active trachoma in the Solomon Islands. Numbers show p values for logistic regression comparison between active trachoma case and control groups for each pathogen. None of the differences observed were statistically significant. TF, trachomatous inflammation–follicular; TI, trachomatous inflammation—intense.
16S rRNA Gene Sequencing
The total number of genera identified across all clinical specimens and NTCs was 659; however, 125/659 (19%) had a cumulative total of 10 reads or fewer. The median number of reads per clinical specimen was 55,104 (interquartile range [IQR]: 40,705–89,541), and the median percentage of reads mapped to genus level was 51.9%. The median number of genera identified was 151 (IQR: 117–212) per clinical specimen. After the removal of presumed contaminants, approximately 40% of cleaned reads were assigned to genera that each constituted less than 1% of the total polymicrobial community. 34 genera were represented by at least 1% of remaining reads in clinical samples.
The diversity of bacterial communities was in general low; the mean Inverse Simpsons Diversity index over all specimens was 0.061 (IQR: 0.046–0.092). The median Inverse Simpson’s Diversity index in cases was 0.061 (IQR: 0.048–0.095) versus controls 0.060 (IQR: 0.044–0.078) (Student’s t-test p = 0.56). The median Shannon Diversity Index was 3.38 in cases and 3.37 in controls (Student’s t-test p = 0.96). There was no significant difference in alpha-diversity between cases and controls by either measure of diversity. In cases of active trachoma, the most dominant bacterial genera were Corynebacterium (12.0% of total reads), Propionibacterium (6.2%), and Helicobacter (4.8%). In controls, the most dominant genera were Corynebacterium (13.9%), Paracoccus (5.2%), Propionibacterium (4.7%), and Neisseria (4.1%; Figure S1 in Supplementary Material). The genus level membership of the bacterial community varied between cases and controls. Of 21 genera represented by at least 1% of reads in specimens from those with active trachoma, 10 were not found in controls, including Helicobacter, Mesoplasma, Brachybacterium, and Haemophilus. Of 23 genera found in controls, 12 were not found in cases, including Neisseria, Prevotella, Rhodococcus, and Porphyromonas.
PC analysis (Figure 2; Figure S2 in Supplementary Material) revealed that Corynebacterium (PC1), Paracoccus (PC2), Propionibacterium (PC2 and PC3), and Helicobacter (PC3) are major contributors to the variation in the conjunctival microbiomes of children in the Solomon Islands. While individual cases or controls have distinctive profiles dominated by these genera, the majority of cases and controls are indistinguishable (Figure 2) using the first three PCs (explaining 45% of total variation). We condensed PCs 1–20 into a single discriminant function, which explained 87% of variation in the data. In this analysis, the genus level polymicrobial community structure is significantly different between cases and controls (logistic regression p = 0.0000062; Figure 3). That discriminant function was dominated by a number of genera such as Curvibacterium and Mesoplasma (Figure S3 in Supplementary Material). Cross-validation to predict group membership using discriminant functions of between 10 and 90 PCs had a low success rate (median 53.3% success; range: 50.1–62.3%). This suggests that although differences do exist between the polymicrobial communities in cases and controls, they are too subtle and varied to be predictive of phenotype, at least when working with this number of specimens.
Figure 2. (A) First and second and (B) second and third principal components describing variation between the 16S sequences identified in Solomon Island children with and without active trachoma. Dark blue spots indicate controls. Light blue spots indicate cases. Red arrows show principal component loadings. PC, principal component; TF, trachomatous inflammation–follicular; TI, trachomatous inflammation—intense.
Figure 3. Discriminant analysis of the association of 20 combined principal components with active trachoma, showing a significant discrimination of phenotype groups (p = 0.000006). TF, trachomatous inflammation–follicular; TI, trachomatous inflammation—intense.
Follicular conjunctivitis meeting criteria for the trachoma phenotype TF is highly prevalent in the Solomon Islands, but the prevalence of Ct infection is curiously low. Although Ct was the only pathogen to associate with TF in this study, the prevalence was much lower than may be expected of a population with this level of TF (3). Both S. pneumoniae and H. influenzae were detected at moderate prevalence in our population, while we found very few cases of Adenoviridae and M. catarrhalis infection. Comparator data from other populations are scarce, but it is clear that staphylococcal species were detected in the conjunctivae in the Solomon Islands at a substantially lower prevalence (7%) than in children in Tanzania (14.8% S. epidermidis) (7), Sierra Leone (20% S. aureus and 29% coagulase-negative Staphylococcus) (38), or The Gambia (post-MDA, 14.7% S. aureus) (11). By ddPCR testing for a number of bacteria and viruses that have previously been linked to TF, we can discount the possibility that any of them can account either singly or en masse for the high TF prevalence.
Through 16S amplicon sequencing and community profiling, we have shown that there is apparently no dominant bacterial genus associated with this disease. The dominant features of variation in conjunctival bacterial communities in Solomon Islands children (Corynebacterium and Propionibacterium) have consistently been identified in other studies. Other genera (e.g., Paracoccus, Helicobacter, Haemophilus) have not previously been identified as being important contributors to the conjunctival microbiome in 16S studies, while some reported in other studies (e.g., Simonsiella, Pseudomonas) were not found in these specimens (18–21). Many studies have sequenced the V3–4 region of the 16S gene, whereas we targeted V1–3. V region choice has been shown to have an impact on the genera identified at other mucosal sites (39), and while no data have yet emerged on how this affects profiles from the eye, it is possible this could account for some of the differences. Consistent with previous studies (36), we found a background of genera that amplified in NTCs that also appeared in clinical specimens. The microbiological biomass at the conjunctiva is known to be very low, even compared with other nearby sites such as skin or oral mucosa (19). When true resident bacteria are scarce, 16S rRNA gene PCR readily amplifies reagent contaminants. We took stringent measures to exclude reagent contaminants, but this resulted in the removal of some important genera such as Staphylococcus and Streptococcus. Among the biggest contributors to between-conjunctiva variation were Corynebacterium and Propionibacterium, but these genera did not differentiate TF cases from controls (Figure 2). Previous investigations of the role of the conjunctival microbiome in ocular disease have also not shown significant differences. In The Gambia, there was an increased abundance of Haemophilus in cases of TF, compared to controls, although this difference was not significant (21). People with ocular manifestations of rosacea, Sjögren’s syndrome, and healthy controls did not have significant differences in diversity or relative abundance of key phyla identified across all three groups (23). However, other disease associations have been identified. For example, the abundance of Corynebacterium and Streptococcus varied significantly between individuals with and without conjunctival scarring (21).
The microbiota of these children were heterogeneous and had subtle variations that appeared to reach significant association with TF when the full community structure was considered in comparative statistical models. It remains possible (though perhaps unlikely) that a multitude of factors, operating at the level of single bacterial species, polymicrobial communities, or viral species, all contribute to the presentation of the phenotype. It is perhaps more likely that as yet unidentified viral or allergic causes could explain the prevalence of TF in the Solomon Islands.
There are some limitations to our study. First, although the ddPCR assays are based on the validated assays and appear to be accurate (Table S1 in Supplementary Material), they have not been formally evaluated in this format. Second, by not using mechanical lysis in our extraction process, the absolute prevalence of some difficult-to-lyse Gram-positive genera (e.g., Staphylococcus) may have been underestimated. This would not impair the comparison of cases to controls, between which protocols were consistent. Third, 16S amplicon sequencing of low biomass samples is known to result in amplification of reagent and environmental contaminants (36). This study focused on gross differences in community structure between those with and without disease, and these should be independent of contaminating reads. We also employed stringent quality control measures to ensure data were not confused by artifacts of the sequencing process.
Given the current international commitment to trachoma elimination, further characterization of the role of non-Ct stimuli in TF is important. We might expect those with TF in the absence of Ct to be at lower risk of progression to TS and TT than those with repeated Ct infection and inflammation (40); however, such individuals have not been studied longitudinally to assess outcomes, and therefore, it is unclear how these cases should be managed. Furthermore, years after MDA has been administered, TF persists at levels above the threshold for elimination (> 5%) in communities where Ct infection is not readily detectable (41, 42). The public health risk of trachoma in these communities is also unclear.
Public health scenarios such as the one in the Solomon Islands will become increasingly common as trachoma elimination programmes reduce the global prevalence of Ct, and the positive predictive value of TF for ocular Ct infection declines as a result (43). MDA might be inappropriately delivered when clinical signs of trachoma are the only indicator used for programmatic decision-making. While MDA is effective for the treatment of trachoma (44), mass antibiotic exposure may increase macrolide resistance in non-chlamydial bacteria (45) and theoretically could make the population more susceptible to later infections by preventing accumulation of acquired immunity to Ct (46). Therefore, decisions to undertake such a program should be carefully considered. The case for using tests for Ct infection during trachoma surveys is strengthened by our data. The value of nucleic acid amplification tests for detecting non-chlamydial infection remains questionable, as multiple infectious agents might be important, and these may differ between populations.
A parent or guardian gave written informed consent for each child to take part in the survey in accordance with the Declaration of Helsinki. The protocol was approved by the London School of Hygiene & Tropical Medicine (6319/6360) and the Solomon Islands National Health Research (HRC13/18) Ethics Committees.
Conceived and designed the study: RB, OS, RM, AS, DM, and ChR. Performed the fieldwork: RB, OS, KJ, EK, LS, and CR. Performed the experiments: RB, MH, JH, CP, and ChR. Analyzed the data: RB and ChR. Wrote the manuscript: RB and ChR. Revised and approved the manuscript: RB, OS, KJ, EK, LS, CR, JH, CP, MH, RM, AS, DM, and ChR.
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 thank Victoria Miari and Emma Cobb (Pathogen Molecular Biology Department, London School of Hygiene & Tropical Medicine, UK) for providing bacterial culture isolates for PCR assay validation; and David Nelson (Indiana University, Indianapolis, USA) for the provision of 16S sequencing primers and protocols. We would also like to thank the residents of the communities who took part in this study.
Fieldwork was jointly funded by the United Kingdom’s Department for International Development Global Trachoma Mapping Project grant (ARIES: 203145) to Sightsavers, and by the Fred Hollows Foundation, Australia (1041). Laboratory costs were funded by the Fred Hollows Foundation, Australia (1041). RB and AS were funded by a Wellcome Trust Intermediate Fellowship (098521). OS, KJ, EK, LS, and CR were employees of the Solomon Islands Ministry of Health and Medical Services, but fees for conducting the survey were paid by the study grants. MH was supported by the Wellcome Trust (093368/Z/10/Z). ChR was supported by the Wellcome Trust Institutional Strategic Support Fund (105609/Z/14/Z).
The Supplementary Material for this article can be found online at http://www.frontiersin.org/articles/10.3389/fmed.2017.00251/full#supplementary-material.
1. Bourne RA, Stevens GA, White RA, Smith JL, Flaxman SR, Price H, et al. Causes of vision loss worldwide, 1990-2010: a systematic analysis. Lancet Glob Health (2013) 1:e339–49. doi:10.1016/S2214-109X(13)70113-X
3. Butcher RMR, Sokana O, Jack K, Macleod CK, Marks ME, Kalae E, et al. Low prevalence of conjunctival infection with Chlamydia trachomatis in a treatment-naïve trachoma-endemic region of the Solomon islands. PLoS Negl Trop Dis (2016) 10:e0004863. doi:10.1371/journal.pntd.0004863
4. Thein J, Zhao PN, Liu HS, Xu JJ, Jha HC, Miao YH, et al. Does clinical diagnosis indicate ocular chlamydial infection in areas with a low prevalence of trachoma? Ophthalmic Epidemiol (2002) 9:263–9. doi:10.1076/opep.9.4.263.1508
5. Grassly NC, Ward ME, Ferris S, Mabey DC, Bailey RL. The natural history of trachoma infection and disease in a Gambian cohort with frequent follow-up. PLoS Negl Trop Dis (2008) 2:e341. doi:10.1371/journal.pntd.0000341
7. Burton MJ, Hu VH, Massae P, Burr SE, Chevallier C, Afwamba IA, et al. What is causing active trachoma? The role of nonchlamydial bacterial pathogens in a low prevalence setting. Invest Ophthalmol Vis Sci (2011) 52:6012–7. doi:10.1167/iovs.11-7326
9. Andersson P, Harris SR, Seth Smith HM, Hadfield J, O’Neill C, Cutcliffe LT, et al. Chlamydia trachomatis from Australian aboriginal people with trachoma are polyphyletic composed of multiple distinctive lineages. Nat Commun (2016) 7:10688. doi:10.1038/ncomms10688
10. Lietman T, Brooks D, Moncada J, Schachter J, Dawson C, Dean D. Chronic follicular conjunctivitis associated with Chlamydia psittaci or Chlamydia pneumoniae. Clin Infect Dis (1998) 26:1335–40. doi:10.1086/516373
11. Burr SE, Hart JD, Edwards T, Baldeh I, Bojang E, Harding-Esch EM, et al. Association between ocular bacterial carriage and follicular trachoma following mass azithromycin distribution in The Gambia. PLoS Negl Trop Dis (2013) 7:e2347. doi:10.1371/journal.pntd.0002347
13. Burton MJ, Adegbola RA, Kinteh F, Ikumapayi UN, Foster A, Mabey DC, et al. Bacterial infection and trachoma in the Gambia: a case control study. Invest Ophthalmol Vis Sci (2007) 48:4440–4. doi:10.1167/iovs.07-0315
14. Cox JT, Kasubi MJ, Mu BE, Zambrano AI, Greene GS, Mkocha H, et al. Trachomatous scarring and infection with non–Chlamydia trachomatis bacteria in wmoen in Kongwa, Tanzania. Invest Ophthalmol Vis Sci (2017) 58:3249–53. doi:10.1167/iovs.17-21519
19. Doan T, Akileswaran L, Andersen D, Johnson B, Ko N, Shrestha A, et al. Paucibacterial microbiome and resident DNA virome of the healthy conjunctiva. Invest Opthalmol Vis Sci (2016) 57:5116. doi:10.1167/iovs.16-19803
21. Zhou Y, Holland MJ, Makalo P, Joof H, Roberts CH, Mabey D, et al. The conjunctival microbiome in health and trachomatous disease: a case control study. Genome Med (2014) 6:99. doi:10.1186/s13073-014-0099-x
24. Abusleme L, Hong B-Y, Dupuy AK, Strausbaugh LD, Diaz PI. Influence of DNA extraction on oral microbial profiles obtained via 16S rRNA gene sequencing. J Oral Microbiol (2014) 6:23990. doi:10.3402/jom.v6.23990
25. Roberts CH, Last A, Molina-Gonzalez S, Cassama E, Butcher R, Nabicassa M, et al. Development and evaluation of a next-generation digital PCR diagnostic assay for ocular Chlamydia trachomatis infections. J Clin Microbiol (2013) 51:2195–203. doi:10.1128/JCM.00622-13
27. Carvalho M, da GS, Tondella ML, McCaustland K, Weidlich L, McGee L, et al. Evaluation and improvement of real-time PCR assays targeting lytA, ply, and psaA genes for detection of pneumococcal DNA. J Clin Microbiol (2007) 45:2460–6. doi:10.1128/JCM.02498-06
28. Meyler KL, Meehan M, Bennett D, Cunney R, Cafferkey M. Development of a diagnostic real-time polymerase chain reaction assay for the detection of invasive Haemophilus influenzae in clinical samples. Diagn Microbiol Infect Dis (2012) 74:356–62. doi:10.1016/j.diagmicrobio.2012.08.018
29. Greiner O, Day PJ, Altwegg M, Nadal D. Quantitative detection of Moraxella catarrhalis in nasopharyngeal secretions by real-time PCR. J Clin Microbiol (2003) 41:1386–90. doi:10.1128/JCM.41.4.1386-1390.2003
30. Kais M, Spindler C, Kalin M, Ortqvist A, Giske CG. Quantitative detection of Streptococcus pneumoniae, Haemophilus influenzae, and Moraxella catarrhalis in lower respiratory tract samples by real-time PCR. Diagn Microbiol Infect Dis (2006) 55:169–78. doi:10.1016/j.diagmicrobio.2006.01.007
31. Nakagawa S, Taneike I, Mimura D, Iwakura N, Nakayama T, Emura T, et al. Gene sequences and specific detection for Panton-Valentine leukocidin. Biochem Biophys Res Commun (2005) 328:995–1002. doi:10.1016/j.bbrc.2005.01.054
32. Okolie CE, Wooldridge KG, Turner DP, Cockayne A, James R. Development of a new pentaplex real-time PCR assay for the identification of poly-microbial specimens containing Staphylococcus aureus and other staphylococci, with simultaneous detection of staphylococcal virulence and methicillin resistance markers. Mol Cell Probes (2015) 29:144–50. doi:10.1016/j.mcp.2015.03.002
33. Last AR, Roberts CH, Cassama E, Nabicassa M, Molina-Gonzalez S, Burr SE, et al. Plasmid copy number and disease severity in naturally occurring ocular Chlamydia trachomatis infection. J Clin Microbiol (2013) 52:324. doi:10.1128/JCM.02618-13
35. Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol (2007) 73:5261–7. doi:10.1128/AEM.00062-07
36. Salter SJ, Cox MJ, Turek EM, Calus ST, Cookson WO, Moffatt MF, et al. Reagent and laboratory contamination can critically impact sequence-based microbiome analyses. BMC Biol (2014) 12:87. doi:10.1186/s12915-014-0087-z
37. Jombart T, Devillard S, Balloux F, Falush D, Stephens M, Pritchard J, et al. Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genet (2010) 11:94. doi:10.1186/1471-2156-11-94
38. Capriotti JA, Pelletier JS, Shah M, Caivano DM, Ritterband DC. Normal ocular flora in healthy eyes from a rural population in Sierra Leone. Int Ophthalmol (2009) 29:81–4. doi:10.1007/s10792-008-9196-4
39. Walker AW, Martin JC, Scott P, Parkhill J, Flint HJ, Scott KP. 16S rRNA gene-based profiling of the human infant gut microbiota is strongly influenced by sample processing and PCR primer choice. Microbiome (2015) 3:26. doi:10.1186/s40168-015-0087-4
40. Wolle MA, Muñoz BE, Mkocha H, West SK. Constant ocular infection with Chlamydia trachomatis predicts risk of scarring in children in Tanzania. Ophthalmology (2009) 116:243–7. doi:10.1016/j.ophtha.2008.09.011
41. Martin DL, Bid R, Sandi F, Goodhew EB, Massae PA, Lasway A, et al. Serology for trachoma surveillance after cessation of mass drug administration. PLoS Negl Trop Dis (2015) 9:e0003555. doi:10.1371/journal.pntd.0003555
42. Bird M, Dawson CR, Schachter JS, Miao Y, Shama A, Osman A, et al. Does the diagnosis of trachoma adequately identify ocular chlamydial infection in trachoma-endemic areas? J Infect Dis (2003) 187:1669–73. doi:10.1086/374743
43. Ramadhani AM, Derrick T, Macleod D, Holland MJ, Burton MJ. The relationship between active trachoma and ocular Chlamydia trachomatis infection before and after mass antibiotic treatment. PLoS Negl Trop Dis (2016) 10(10):e0005080. doi:10.1371/journal.pntd.0005080
44. Solomon AW, Holland MJ, Alexander NDE, Massae PA, Aguirre A, Natividad-Sancho A, et al. Mass treatment with single-dose azithromycin for trachoma. N Engl J Med (2004) 351:1962–71. doi:10.1056/NEJMoa040979
45. Seidman JC, Coles CL, Silbergeld EK, Levens J, Mkocha H, Johnson LB, et al. Increased carriage of macrolide-resistant fecal E. coli following mass distribution of azithromycin for trachoma control. Int J Epidemiol (2014) 43:1105–13. doi:10.1093/ije/dyu062
Keywords: droplet digital PCR, Solomon Islands, 16S rRNA gene sequencing, trachoma, non-chlamydial infections
Citation: Butcher RMR, Sokana O, Jack K, Kalae E, Sui L, Russell C, Houghton J, Palmer C, Holland MJ, Le Mesurier RT, Solomon AW, Mabey DCW and Roberts Ch (2018) Active Trachoma Cases in the Solomon Islands Have Varied Polymicrobial Community Structures but Do Not Associate with Individual Non-Chlamydial Pathogens of the Eye. Front. Med. 4:251. doi: 10.3389/fmed.2017.00251
Received: 27 September 2017; Accepted: 19 December 2017;
Published: 23 January 2018
Edited by:Olivier Vandenberg, Université libre de Bruxelles, Belgium
Reviewed by:Aleksandra Barac, University of Belgrade, Serbia
Mariane B. Melo, Massachusetts Institute of Technology, United States
Copyright: © 2018 Butcher, Sokana, Jack, Kalae, Sui, Russell, Houghton, Palmer, Holland, Le Mesurier, Solomon, Mabey and Roberts. 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: Robert M. R. Butcher, firstname.lastname@example.org