Bacterial DNA was recovered from all eyes, nasal mucosa, and pharyngeal mucosa of all patients (100 eyes of 50 patients, 50 nasal mucosa, 50 pharyngeal mucosa) using 16S sequencing. The maximum total number of de novo raw OTUs was 776, which was decreased to 745 after filtering. The total number of reads was 70,060,846, with a combined number of sequences (total filtered reads + strain hit reads) of 32,515,209 and 31,356,744 sequences after removal of OTUs unclassified at the kingdom level. Filtering reduced the final number of OTUs to 551. OTUs from the sequences that passed sample quality check were used for downstream analyses. One hundred percent of sequences were classified at the kingdom level, 99.97% at the phylum and class levels, 99.95% of sequences were classified at order, 99.45% at family, 92.28% at genus, 1.413% at species, and 1.664% at strain levels (
Supplementary Fig. S1). Rarefaction curves reached saturation, which is indicative that samples were sequenced to sufficient depth and composition from samples from all conditions were likely captured. The top five most abundant phyla (in order of abundance) were:
Firmicutes,
Proteobacteria,
Actinobacteria,
Cyanobacteria, and
Bacteroidetes. The top seven most abundant families (in order of abundance) were:
Staphylococcaceae, Streptococcaceae, Corynebacteriaceae, Moraxellaceae, Enterobacteriaceae, Oceanospirillaceae, and
Bacillaceae.
No differences were noted between right and left eyes (
Supplementary Fig. S2). The OTU richness for the right eye was 80.6 (SD 30.8) alpha diversity units/condition and the left eye was 79 (SD 26.2) (
P = 0.49); similarly, the Shannon diversity was 0.632 (SD 0.877) in the right eye compared to 0.585 (SD 0.722) in the left eye. The proportional abundance for the top three bacterial classes at the family level was similar between the eyes (right versus left eye):
Staphylococcaceae (57.8 vs. 55.1),
Streptococcaceae (12.7 vs. 21), and
Corynebacteriaceae (8.76 vs. 3.73).
Similar to culture dependent methods, older infants and children (>6 months old) tended to have greater diversity of the OSM than young infants (<6 months old). Of the 264 OTUs tested, 75 were significantly different. Forty-two of the 75 significantly different OTU were enriched in children older than 6 months. The phylum breakdown for these significantly different OTUs in the older children included eight enriched in the Actinobacteria phylum, 18 enriched the Firmicutes phylum, and 16 enriched in the Proteobacteria phylum. Of note, older children also had 21 reduced OTUs in the Firmicutes phylum and six reduced OTUs in the Proteobacteria phylum when compared to young infants. Of the 75 significant OTUs, 18 were able to be identified at the strain level of which 14 were enriched in the over 6-month age class. They included three strains from Corynebacterium (C. tuberculostearicum, C. bovis, and C. mastitidis) and two strains from Psychromonas (P. arctica and P. ingrahamii). The observed alpha-diversity estimate was lower in children >6 months old compared to young infants (84 ± 19.4 vs. 93.1 ± 15.3); however, the Shannon diversity was higher in older children (0.986 ± 0.903 vs. 0.458 ± 0.59). The percent proportional abundance of the top three families also differed between older (>6 months) and younger (<6 months) groups: Staphylococcaceae (44.1 ± 45.9 vs. 58.2 ± 45.5), Streptococcaceae (18.7 ± 29.5 vs. 23.5 ± 40.7), and Moraxellaceae (3.49 ± 7.19 vs. 12.5 ± 26.6); however, these differences were not statistically significant. Older children did have statistically significant differences in the relative abundance of the Oceanospirillaceae (7.32 vs. 0), Listeriaceae (4.42 vs. 0), Psychomonadaceae (2.57 vs. 0.002), and Leuconostocaceae (2.07 vs. 0) (P < 0.001 for all, Kruskal-Wallis rank sum test). There were significant beta diversity differences between the two age classes (P = 0.05), regardless of factors such as sex or race.
When comparing OSM of the eye to that of the nose and the throat, the eye microbiome was more similar to the nose than the throat. Throat samples had higher OTU richness than eye samples and differed significantly (
Fig. 1 and
Supplementary Fig. S3). Sixty-six of the 137 significantly different OTUs were enriched, with greater
Bacteroidetes, Firmicutes, Fusobacteria, and
Proteobacteria in the throat. Thirty-three of the 137 significant OTUs were able to be identified at the strain level. There were three strains from
Staphylococcus (
S. condimenti, S. devriesei, and
S. aureus) and three strains from
Corynebacterium (
C. tuberculostearicum, C. propinquum, and
C. bovis) enriched in the eye samples. There were three strains from
Streptococcus (
S. HKU30, S. uberis, and
S. infantis) and five strains from
Prevotella (
P. veroralis, P. nanceiensis, P. pallens, P. sp. oral taxon 299, and
P. salivae) enriched in throat samples. Weighted ordination using abundance reveals most samples separated according to mucosal type (
P = 0.001). Wilcoxon signed rank test performed on alpha diversity measure of the sampled mucosa and revealed no significant differences between the OTU richness of the right versus left eye (
P = 0.36) or the eye versus nose (
P = 0.28); however, there were differences between the eyes and throat (
P = 0.03) and the nose and throat (
P = 0.004). Similar findings were elicited when analyzing Shannon diversity (right eye versus left eye
P = 0.94, eye versus nose
P = 0.14, eye versus throat
P = 0.03, and nose versus throat
P = 0.032). The throat had significantly higher relative abundance of
Streptococcaceae than eye samples (throat 71.6 ± 35.2 versus eye 19.9 ± 33.1) and nose samples (throat 71.6 ± 35.2 versus nose 30.2 ± 36.8). Eye samples had significantly higher relative abundance of
Staphylococcaceae, Moraxellaceae, Corynebacteriaceae, Oceanospirillaceae, and
Listeriaceae than throat samples. Nasal samples had significantly higher relative abundance of
Staphylococcaceae, Corynebacteriaceae, and
Veillonellaceae than throat samples.