- 1Ophthalmology and Optometry Department, Affiliated Eye Hospital of Shandong University of Traditional Chinese Medicine, Jinan, China
- 2Medical College of Optometry and Ophthalmology, Shandong University of Traditional Chinese Medicine, Jinan, China
- 3Shandong Key Laboratory of Integrated Traditional Chinese and Western Medicine for Prevention and Therapy of Ocular Diseases, Jinan, China
Introduction: The increasing prevalence of myopia worldwide poses significant public health concerns. Accumulating evidence suggests a potential link between ocular diseases and the gut microbiota (GM); however, whether the GM directly contributes to myopia development remains to be established.
Methods: This study investigated the potential causal link between the GM and myopia through bidirectional Mendelian randomization (MR) analysis, further validated by experiments conducted on a form-deprivation myopia (FDM) guinea pig model. Bidirectional two-sample MR analysis was performed using genome-wide association study summary statistics comprising data on 196 GM taxa from the MiBioGen consortium and myopiaassociated data from the FinnGen database. Instrumental variables were carefully selected according to predetermined standards. Subsequently, MR estimates were calculated using the inverse variance weighted, MR-Egger regression, and weighted median approaches, along with supplementary sensitivity evaluations. Concurrently, FDM was experimentally induced in guinea pigs, and fecal samples were subjected to comprehensive full-length 16S rRNA gene sequencing analysis.
Results: MR analysis identified five bacterial taxa linked to the risk of myopia. Specifically, higher Bifidobacterium abundance was associated with lower myopia risk (odds ratio = 0.834, 95% confidence interval = 0.705–0.986, p < 0.05). Animal experiments validated the MR findings, demonstrating a significant enrichment of Bifidobacteria in control animals.
Discussion: Conclusively, supplementation with Bifidobacteria is a potential strategy for reducing the risk of myopia. Future research should focus on developing and testing Bifidobacterium-based interventions to validate their effectiveness in controlling myopia.
1 Introduction
Estimates suggest that almost half of the world’s population will develop myopia by 2050, imposing a significant burden on the global healthcare systems (1). Myopia development is influenced by both hereditary and environmental influences (2). For instance, prolonged near work (especially at short viewing distances) and higher levels of education are associated with an increased risk of myopia, whereas spending more time outdoors tends to be protective (3). Notably, high myopia substantially increases the likelihood of vision-threatening eye disorders such as retinal detachment, glaucoma, and cataract formation (4, 5). Uncorrected refractive errors and pathological myopia are major contributors to visual impairment and, in some cases, blindness (6, 7). Therefore, preventing or delaying the onset of myopia is crucial.
Recent research has revealed a potential connection between myopia and gut microbiota (GM). The GM constitutes a complex ecosystem comprising approximately 4 × 1013 symbiotic organisms, including bacteria, protozoa, fungi, archaea, and viruses (8, 9). Over the past decade, mounting evidence has identified gut microbial dysbiosis as a crucial factor underlying the initiation and progression of numerous diseases (10). Recent findings indicate GM involvement in various ocular pathologies, including optic neuritis (11), age-related macular degeneration (12), and diabetic retinopathy (13). Nevertheless, the relationship between GM and myopia remains incompletely characterized and warrants further investigation. Li et al. (14) demonstrated significant alterations in GM and metabolite profiles in myopic mice, with dysbiosis and plasma metabolomic changes suggesting gut barrier disruption. Subsequent investigations have explored GM-associated metabolic pathways; Wu et al. (15) characterized Ruminococcus albus and glutamate metabolic pathways as critical factors in form-deprivation myopia (FDM). In a cohort analysis of 52 adults, Omar et al. (16) observed that myopic participants displayed increased abundances of multiple gut bacterial genera, including Bacteroides, Faecalibacterium, Dorea, Roseburia, and Blautia, relative to non-myopic controls. While these investigations have identified associations, their limited sample sizes and absence of causal inference preclude establishing a definitive causal relationship between GM and myopia.
Mendelian randomization (MR) approaches have gained prominence in recent years for identifying causal risk factors underlying myopia (17–19). Notably, Mountjoy et al. employed bidirectional MR analysis to establish that extended educational duration contributes to increasing myopia prevalence (19). Additionally, bidirectional MR revealed a strong bidirectional genetic causal relationship between myopia and primary open-angle glaucoma, primarily mediated through intraocular pressure (17). MR leverages genetic variations as “natural randomized experiments” to replicate the effects of randomized controlled trials. This approach provides a robust framework for evaluating the causal relationships between modifiable exposures or risk factors and clinically relevant outcomes (20). Through application of Mendelian segregation and independent assortment principles, MR minimizes confounding biases that frequently compromise traditional observational epidemiology, enabling more precise delineation of causal pathways in complex phenotypic traits (21).
This investigation employed bidirectional two-sample MR analysis to comprehensively evaluate potential causal relationships between diverse GM taxa and myopia. To explore whether myopia might reciprocally influence GM composition, we performed reverse MR analysis. For validation of MR findings, we conducted full-length 16S rRNA gene sequencing on fecal samples from form-deprived myopic guinea pigs using third-generation sequencing technology. Full-length 16S rRNA gene sequencing provides enhanced species-level taxonomic resolution compared to traditional short-read 16S approaches (22). These advanced techniques enable the accurate identification of microbial taxa at the species and strain levels, facilitating the precise characterization of myopia-associated microbiota structures. This investigation, to our knowledge, represents the first integration of bidirectional MR with full-length 16S rRNA sequencing to elucidate the gut microbial taxa contributing to the onset and progression of myopia. Our findings provide a foundational framework for future investigations exploring gut microbiota modulation as a therapeutic strategy for myopia management.
2 Materials and methods
2.1 Study design
Using a two-sample bidirectional MR framework, this study examined potential causal relationships between GM and myopia onset by leveraging summary-level data from an independent genome-wide association study. For the MR methodology to yield valid inferences, it was necessary to adhere to three foundational assumptions: (1) the chosen genetic instruments must demonstrate robust statistical associations with the exposure variable; (2) these instrumental variants should exert effects on the outcome solely through their influence on the exposure of interest; and (3) genetic proxies must remain unassociated with confounding variables that could otherwise distort the observed exposure-outcome association. This analytical approach ensures a rigorous evaluation of causality while minimizing bias from pleiotropic pathways or residual confounding (23). Guinea pigs were selected to validate the in vivo experiments for MR analysis (Figure 1).

Figure 1. Overview of the research workflow. FDM, form-deprivation myopia; NC, normal control; SNPs, single nucleotide polymorphisms. (A) Illustrates the analytical workflow of Mendelian randomization; (B) depicts the experimental workflow of 16S rRNA sequencing.
2.2 Data sources
GM data encompassing 196 taxa were obtained from the MiBioGen consortium, which performed meta-analysis on data of 18,340 individuals (24). Myopia outcome data were acquired from the Finnish FinnGen Biobank. FinnGen identified myopia cases using ICD-10 code H52.1 (nearsightedness) within health registries (https://r9.risteys.finngen.fi/endpoints/H7_MYOPIA). The FinnGen dataset comprised 4,106 myopia cases and 394,028 controls of European ancestry with over 1.0 million genotyped SNPs. All database populations consisted of individuals of European descent. Data were fully anonymized, with no individual-level information accessible for this investigation.
2.3 IVs selection
SNPs associated with the GM were selected with a combined threshold (p < 1 × 10−5) and linkage disequilibrium analysis (PLINK v1.9) with a distance >10,000 kb and r2 < 0.001, ensuring the robustness of Instrumental variables (IVs) validated by F-statistic >10. For myopia-associated IVs in the reverse MR analysis, a stricter threshold (p < 5 × 10−6) was applied (Supplementary Tables S1, S2).
2.4 MR analysis and sensitivity analysis
The inverse variance weighted (IVW) method served as our primary analytical approach for evaluating potential associations between GM and myopia. We performed additional analyses using MR-Egger regression and weighted median models to assess the robustness of causation estimates. Consistent directional results across different MR analytical models strengthened the credibility of causal inferences. We employed the MR-Egger intercept test to detect horizontal pleiotropy (25) and utilized MR-PRESSO to identify outliers and adjust for horizontal pleiotropy, with statistical significance set at p < 0.05 (26), indicating potential direct instrumental variable effects on outcomes independent of exposure. We applied Cochran’s Q test to assess heterogeneity among genetic variants (27), with p < 0.05 indicating significant heterogeneity. SNPs demonstrating high heterogeneity were excluded from analysis. We performed leave-one-out analysis to assess potential bias in MR estimations attributable to individual SNPs (28). Statistical analyses were conducted using R software (version 4.2.1, R Foundation for Statistical Computing, Vienna, Austria) with the two-sample MR package (version 0.5.6) and MR PRESSO package (version 1.0).
2.5 Myopia model and groups
We obtained 18 clinically healthy male tricolor guinea pigs from Jinan Jinfeng Laboratory Animal Co. and housed them in a dedicated, ventilated facility under controlled conditions: 12-h light–dark cycles, temperature maintained at 21–25°C, and relative humidity at 40–70%. Animals received food and water ad libitum. We randomly allocated animals into FDM and normal control (NC) groups (n = 9 each). The FDM model was selected for its reliable induction of significant axial elongation in guinea pigs within short timeframes. Xiao et al. (29) demonstrated no significant differences between FDM and lens-induced myopia models regarding final refraction or axial length, noting guinea pigs’ high sensitivity to FDM. This approach is extensively utilized in myopia research. In the FDM group, we induced myopia by occluding the right eye with translucent latex for 4 weeks, confirming myopia development through increased axial length (AL) and decreased spherical equivalent. NC animals received no experimental intervention. We excluded guinea pigs presenting with ocular disease or pre-existing myopia. Statistical analyses compared treated eyes from FDM animals with eyes from normal control animals, rather than employing paired within-animal comparisons. All procedures received approval from the Ethical Review Committee for Laboratory Animals of Shandong Eye Disease Control Institute (AEHSDUTCM-LAEC2024001) and complied with the ARVO Statement for the Use of Animals in Ophthalmic and Vision Research. Reporting followed ARRIVE guidelines to ensure transparency and completeness.
2.6 Refractive status and axial length measurement
We assessed guinea pig refractive errors using a streak retinoscope (YZ24; Six Six Vision Technology Co., Ltd., China). Prior to ocular measurements, we administered one drop of tropicamide/phenylephrine eye medication (Santen Pharmaceutical, Japan) to the conjunctival sac at 5-min intervals, repeating three times. Animals were subsequently placed in darkness for a 30-min adaptation period. An experienced optometrist, operating under dark conditions, determined equivalent refractive errors by averaging vertical and horizontal meridian readings and calculating the mean from three consecutive measurements. We measured axial length using A-scan ultrasound biometry (Cinescan, Quantel Medical, France) following device calibration specifically for guinea pig eyes. Sound-propagation velocities were configured at 1557 m s−1 for the anterior chamber, 1723 m s−1 for the crystalline lens, and 1,540 m s−1 for the vitreous cavity. Prior to each recording, we applied topical anesthesia to the corneal surface using one drop of 0.5% proparacaine hydrochloride (Santen Pharmaceutical, Japan) administered to the conjunctival sac. Final AL values represented the mean of 10 repeated measurements. No animals underwent euthanasia during or following the experiment. All guinea pigs continued under standard husbandry conditions post-experiment, in compliance with ARRIVE guidelines and ethical approval requirements.
2.7 DNA extraction, library preparation, and full-length 16S rRNA sequencing
We collected fecal samples from guinea pigs in both NC and FDM cohorts (nine animals per group). Samples underwent immediate snap-freezing and storage at −80°C to preserve molecular integrity. We extracted genomic DNA using a TGuide S96 magnetic-bead soil/feces kit (DP812; Tiangen Biochemical, Beijing, China) and evaluated quality using a femto-pulse capillary system. We amplified the bacterial 16S rRNA locus using universal primers 27F (5′-AGRGTTTGATYNTGGCTCAG-3′) and 1492R (5′-TASGGHTACCTTGTTASGACTT-3′). Polymerase chain reactions (PCRs) were performed in 20 μL volumes containing approximately 10 ng template DNA, 1.5 μL of each primer, 10 μL KOD One™ PCR Master Mix (KMM-101; Bellink Biotechnology, Beijing), and nuclease-free water. We visualized amplicons through agarose gel electrophoresis or LabChip GX Touch analyzer, then purified and converted them into SMRTbell libraries using PacBio Prep Kit v3.0, incorporating damage repair, end polishing, and adapter ligation steps. Libraries underwent cleaning with AMPure PB magnetic beads, quantification on a Qubit fluorometer, and polymerase pairing using the Revio kit. Sequencing was performed for 24 h on a PacBio Revio platform (Annoroad Gene Technology, Beijing, China) following manufacturer’s protocols.
2.8 Data analysis of gut microbiome
The effective CCS sequences were processed for clustering and denoising using USEARCH1 (30). Taxonomic annotation of the operational taxonomic units (OTU) sequences was performed using the naive Bayesian classifier in the SILVA database2 (31). Subsequently, the QIIME2 software3 (32) was used to analyze and compare the biodiversity across individual samples and their subgroups. Normalized data were utilized for both alpha and beta diversity analyses. Community richness and evenness were quantified with the Feature, Chao1, Shannon, and Simpson metrics, whereas β-diversity dissimilarities were derived from binary Jaccard, Bray–Curtis, and (weighted and unweighted) UniFrac distance matrices. The ordination of the resulting distance space was visualized using principal coordinate analysis and non-metric multidimensional scaling. Differentially abundant taxa were identified using the linear discriminant analysis effect size (LEfSe). Species-level abundance tables were further interrogated using Spearman’s rank correlations, retaining network edges where ρ > 0.10 and p < 0.05. Finally, functional potentials were inferred using PICRUSt2 against reference databases, such as KEGG and COG.
2.9 Statistical analysis
Continuous variables are presented as mean ± standard deviation. Comparisons between two groups were conducted using Student’s t-test or the Mann–Whitney U test. Differences between categorical variables were assessed using the chi-squared test. Spearman’s rank correlation analysis was used to evaluate the relationships between microbial taxa. All statistical analyses were performed using R (version 4.2.1, R Foundation for Statistical Computing) and SPSS (version 26, IBM Corp., Armonk, N.Y., United States).
3 Results
3.1 Effect of genetically determined GM on myopia
Through application of the inverse-variance–weighted MR model, we identified causal associations involving five gut bacterial taxa, with two demonstrating positive relationships with increased myopia susceptibility: family Oxalobacteraceae (odds ratio [OR] = 1.145, 95% confidence interval [CI] 1.008–1.301, p < 0.05) and genus Butyrivibrio (OR = 1.104, 95% CI = 1.001–1.218, p < 0.05). Conversely, genus Bifidobacterium (OR = 0.834, 95% CI = 0.705–0.986, p < 0.05), genus Erysipelotrichaceae UCG003 (OR = 0.823, 95% CI = 0.686–0.986, p < 0.05), and genus FamilyXIIIAD3011group (OR = 0.799, 95% CI = 0.643–0.994, p < 0.05) demonstrated protective effects against myopia (Figure 2; Supplementary Tables S3, S4).

Figure 2. Associations between the GM and myopia risk identified through IVW. IVW, inverse variance weighted.
Inverse MR analysis examining the relationship between myopia and GM species revealed that myopia onset was associated with reduced abundance of genus Collinsella (OR = 0.941, 95% CI = 0.887–0.999, p < 0.05), genus Coprococcus 2 (OR = 0.921, 95% CI = 0.853–0.993, p < 0.05), genus Lachnospiraceae NC2004 group (OR = 0.909, 95% CI = 0.830–0.995, p < 0.05), and genus Odoribacter (OR = 0.906, 95% CI = 0.855–0.961, p < 0.05). Myopia onset was associated with increased abundance of genus Escherichia/Shigella (OR = 1.068, 95% CI = 1.001–1.140, p < 0.05; Figure 2; Supplementary Tables S5, S6).
3.2 Sensitivity analysis
To strengthen the inferred causal relationships between the GM and myopia, sensitivity analyses were performed, including a leave-one-out analysis, Q-tests, and Egger’s intercept evaluation. The results indicated no significant heterogeneity or directional horizontal pleiotropy across the IVW and MR-Egger models (p > 0.05; Supplementary Tables S7–S10; Supplementary Figures S1–S10).
3.3 GM in animal experiments profiles
Table 1 illustrates the changes in SE and AL between groups before and after modeling. At baseline, no statistically significant differences existed in SE or AL between groups (all p > 0.05). However, following 4 weeks of form deprivation, the FDM group exhibited a pronounced myopic shift in both refraction (FDM: −4.50 ± 0.43D vs. NC: −0.31 ± 1.23D, p < 0.0001) and AL (FDM: 8.52 ± 0.04 mm vs. NC: 8.25 ± 0.09 mm, p < 0.0001) compared to values in the NC group.
Following data filtration, we obtained 834,367 high-quality sequences across 18 samples, identifying 14,283 OTUs. Rarefaction curves for richness and diversity metrics plateaued for both sample groups, confirming adequate sequencing depth. Alpha diversity analysis revealed no significant differences between the NC and FDM groups across the Feature, Chao, Shannon, and Simpson indices (Figure 3A). However, PERMANOVA analysis based on Bray–Curtis distance matrices demonstrated significant β-diversity differences between the microbial communities of the two groups (p < 0.05; Figure 3B).

Figure 3. (A) Differences in alpha diversity between NC and FDM based on the ob-served Feature, Chao1, Shannon, and Simpson indices. (B) Beta diversity differences between the NC and FDM were estimated by Principle coordinates analysis (PCoA).
3.4 Myopia-related changes in the composition of gut microflora
Taxonomic dependency analysis revealed that both NC and FDM groups harbored 13 phyla, with Firmicutes, Bacteroidota, and Verrucomicrobiota representing the most abundant taxa. Firmicutes dominated both groups, comprising 49.75 and 46.04% of the GM in the FDM and NC groups, respectively. The FDM group exhibited a higher proportion of Bacteroidota compared to the NC group (34.07% vs. 31.82%; Supplementary Figure S11A). Using the Wilcoxon rank-sum test (Mann–Whitney U test), we identified significant differences in the relative abundances of 26 genera between the NC and FDM groups. Supplementary Figure S11B displays the relative abundance heat map.
To elucidate the differences in taxonomic abundance between the FDM and NC groups, we generated a cladogram using the LEfSe methodology (Figure 4A). This analysis revealed 32 OTUs with significant differences (LDA score>3). Unclassified Candidatus Saccharimonas and Pediococcus pentosaceus demonstrated higher abundance in the NC group compared to the FDM group. Conversely, the Eubacterium coprostanoligenes group and unclassified Eubacterium coprostanoligenes group exhibited significantly elevated abundance in the FDM group relative to the NC group (Figure 4B). Importantly, genus Bifidobacterium and species Ileibacterium valens showed significant enrichment in the NC group (LDA score>3).

Figure 4. Gut‑microbiota composition and predicted functions in form‑deprivation myopia (FDM) and normal control (NC) guinea pigs. (A) Relative abundance of bacterial phyla in the FDM and NC groups. (B) Heatmap displaying the 26 operational taxonomic units (OTUs) that differed significantly between groups. (C) Predicted KEGG functions: box plot showing enrichment of signal‑transduction pathways in FDM. (D) Predicted COG categories: box plots showing decreased nucleoside transport and metabolism and increased inorganic‑ion transport and metabolism in FDM.
3.5 Microbial network structure and functional differences
Network analysis unveiled distinct structural differences between NC and FDM groups, despite equivalent node and edge numbers in both networks. The NC group demonstrated superior modularity (0.6396 vs. 0.5327), clustering coefficient (0.467 vs. 0.406), and betweenness centralization (0.259 vs. 0.140). Conversely, the FDM group displayed reduced average path length and increased network diameter (Supplementary Figures S12, S13).
Predictive functional analysis revealed pronounced functional differences between the NC and FDM groups. KEGG analysis demonstrated significant enrichment of differentially expressed GM in signal transduction pathways (p-adjusted = 0.0314), indicating upregulation of these pathways under myopic conditions (Figure 4C). COG analysis revealed significant reduction in nucleoside transport and metabolism within the FDM group (p = 0.0018), while inorganic ion transport and metabolism showed significant elevation (p = 0.024; Figure 4D). Although the NC group exhibited marginally higher relative frequency of translation, ribosomal structure, and biogenesis, these differences lacked statistical significance (p = 0.123).
4 Discussion
This investigation represents the first application of bidirectional two-sample MR combined with 16S rRNA sequencing to examine the potential relationship between GM and myopia risk. We specifically investigated whether particular gut microbial taxa might modulate myopia risk, thereby advancing understanding of the gut-ocular axis. Our MR analysis revealed that family Oxalobacteraceae and genus Butyrivibrio correlated positively with increased myopia risk, while genus Bifidobacterium, genus Erysipelotrichaceae UCG003, and genus Family XIII AD3011 exhibited protective effects against myopia. Inverse MR findings demonstrated associations between myopia onset and decreased abundance of specific bacterial species, including genus Collinsella, genus Coprococcus 2, genus Lachnospiraceae NC2004 group, and genus Odoribacter. Additionally, our analysis identified increased genus Escherichia/Shigella abundance associated with myopia onset. Full-length 16S sequencing of guinea pig fecal samples revealed significant genus Bifidobacterium and species Ileibacterium valens enrichment in the NC group, corroborating our MR findings.
Both MR analysis and 16S rRNA sequencing substantiated the association between Bifidobacterium and decreased myopia risk. While the ORs are small (OR = 0.834), our instruments’ high F- statistics imply that weak bias is unlikely (33). Prior investigations have suggest that Bifidobacterium can ameliorate ocular diseases through inhibition of oxidative stress-induced inflammatory responses (34, 35). Bifidobacterium administration to mice with compromised retinal ganglion cells reduced microglial number and activation while enhancing Müller cell activation. This process coincided with inflammatory cytokine downregulation and neurotrophic factor upregulation (35). As a probiotic organism, Bifidobacterium generates short-chain fatty acids (SCFAs), including acetic acid (36). SCFAs have been reported to modulate inflammatory factor expression (e.g., tumor necrosis factor tumor necrosis factor (TNF)-α, interleukin interleukin (IL)-6, IL-10) through G protein-coupled receptor activation and other receptor pathways (37, 38) or via histone deacetylase inhibition (39). One investigation revealed that proinflammatory cytokines IL-6 and TNF-α synergistically facilitate axial elongation within the inflammatory microenvironment characteristic of myopic eyes (40). IL-6 triggers proinflammatory signaling through soluble IL-6 receptor binding and subsequent JAK/MAPK and NF-κB pathway activation, resulting in oxidative stress, endothelial dysfunction, and retinal inflammation (41–44). Thus, we speculate that SCFAs produced by Bifidobacterium could help down regulate inflammatory factors such as IL-6 and TNF-α, potentially contributing to slower myopic axial elongation. We note, however, that this mechanism remains hypothetical and was not directly shown by our data. Furthermore, Bifidobacterium demonstrates substantial antioxidant capabilities. Reactive oxygen species (ROS), primary mediators of oxidative stress, can be induced by mechanical stretching (45), a process intimately associated with AL elongation during myopia progression (46, 47). Zeng et al. (48) demonstrated that Bifidobacterium NCU-08 attenuated neuronal apoptosis and ROS generation through PI3K/AKT signaling pathway upregulation, consequently enhancing antioxidant enzyme activities, including glutathione peroxidase, superoxide dismutase, and catalase.
Recent investigations have yielded contradictory findings regarding Bifidobacterium’s role in myopia. Li et al. (14) documented enrichment in myopic mice, while Omar et al. (16) reported elevated abundance in myopic humans. These inconsistencies likely stem from species-specific variations in microbial ecology, including functional divergence among murine, human, and guinea pig models, as well as methodological constraints. Cross-sectional human studies cannot adequately control for potential confounding variables, particularly dietary factors. In contrast, our MR analysis combined with animal experimentation, provides compelling evidence supporting Bifidobacterium’s causal protective role. Crucially, Bifidobacterium encompasses diverse subspecies exhibiting strain-specific effects. Protective strains identified in guinea pig models may differ from those present in human myopes. Future investigations should emphasize strain-level characterization within human cohorts and examine probiotic interventions to elucidate context-dependent effects.
Our MR analysis identified Erysipelotrichaceae UCG003, a genus within the Erysipelotrichaceae family, as a protective factor against myopia. Animal experiments revealed significant enrichment of Ileibacterium valens, another Erysipelotrichaceae family member, in the NC group. Erysipelotrichaceae, classified within the phylum Firmicutes, comprise Gram-positive, slender, filamentous, rod-shaped bacteria (49). Erysipelotrichaceae abundance varies considerably across hosts, demonstrating elevation in inflammatory mouse models but reduction in human inflammatory bowel disease (IBD) patients (50, 51). Palm et al. (52) found no significant abundance differences for specific immunogenic Erysipelotrichaceae species between IBD patients and healthy controls; however, subsequent germ-free mouse infection experiments demonstrated that colonization with IgA-coated Erysipelotrichaceae-containing bacterial communities induced more severe colitis responses. Zhao et al. documented increased Erysipelotricaceae incertae sedis levels in mice receiving high-dose gentamicin treatment (53). These findings suggest that Erysipelotrichaceae play a crucial role in intestinal immune defense, with immunoregulatory properties that could be protective in reducing ocular inflammation, which is a key factor in the progression of myopia. Nonetheless, given the preliminary nature of the current research on Erysipelotrichaceae, these results should be interpreted with caution. Future studies should investigate the specific roles of different genera and species within the Erysipelotrichaceae family to clarify their distinct effects on myopia and validate our findings.
The NC group’s microbial network displayed enhanced modularity and cohesion, suggesting a more stable microecological environment potentially better suited to maintaining host-microbe equilibrium. KEGG analysis demonstrated significant GM enrichment in signal transduction pathways. Additionally, COG analysis revealed marked nucleoside transport and metabolism reduction coupled with substantially elevated inorganic ion transport and metabolism in the FDM group. Previous investigations have established the retina’s substantial energy demands (54), and decreased nucleoside transport and metabolism may restrict retinal cell energy supply. Ganglion cells demonstrate particular vulnerability to energy deficiency due to their exceptionally active metabolism (55). Our study observed an increase in inorganic ion transport in the FDM group. Wu et al. (56) similarly documented potassium flux alterations in ciliary muscles of guinea pigs with lens-induced myopia, identifying myopia-induced potassium influx. Additionally, Zhong et al. (57) reported increased Piezo1 expression—a mechanically sensitive ion channel—in FDM guinea pig retinas. Notably, Piezo1 channel function inhibition slows myopia progression and partially reduces retinal ROS levels in FDM guinea pigs. Collectively, these findings suggest that enhanced inorganic ion transport, particularly through mechanisms including inward potassium flux and Piezo1 ion channels, may constitute a critical factor in myopia development. Modulation of these pathways may provide promising therapeutic approaches for addressing retinal oxidative stress and decelerating myopia progression.
Unlike previous investigations, our study employed MR analysis to effectively eliminate confounding factors and reverse causation bias. This approach circumvents limitations inherent in traditional randomized controlled trials, including ethical constraints, extended observation periods, and substantial costs. Our findings underwent rigorous robustness checks to enhance their credibility. Distinguishing our work from previous MR studies, we utilized MR as an exploratory tool for identifying GM potentially influencing myopia, subsequently validating these findings through animal experimentation to ensure comprehensive and reliable results.
Despite employing rigorous analytical methods, the inherent observational nature of Mendelian randomization (MR) and microbiome sequencing studies means residual confounding factors—such as dietary patterns, antibiotic use and lifestyle behaviors—cannot be fully excluded (58, 59). Moreover, both gut-microbiome composition and the distribution of myopia risk factors vary markedly across ethnic and geographic groups, limiting the external validity of our findings beyond the predominantly European ancestry examined. Our microbiome analysis predominantly identifies Bifidobacterium at the genus level, along with a limited number of abundant species. This approach does not allow for the capture of fine-scale strain differentiation or the functional gene content of the microbiota. Given that Bifidobacterium is a heterogeneous genus consisting of numerous species with overlapping 16S rRNA gene sequences but distinct metabolic functions, our findings may not fully reflect the nuanced biological roles of specific strains. Finally, MR estimates represent the lifelong effect of genetically proxied variation. Because myopia typically manifests during childhood and adolescence, the impact of the gut microbiota during critical developmental windows may differ from the average lifetime effect inferred in adult MR analyses. Findings should therefore be interpreted cautiously when extrapolating to age-specific microbiome–host interactions.
Bifidobacterium species demonstrate widespread use as oral probiotics with established long-term safety profiles. A clinical trial administering approximately 3 × 109 CFU/day of B. infantis or B. bifidum to infants for 8 weeks reported no serious adverse events (60). Consequently, Bifidobacterium-based interventions are generally well tolerated. When planning myopia prevention trials, researchers must select appropriate strains, doses, and delivery methods to ensure the most effective outcomes. Oral delivery remains most feasible (e.g., capsules or fortified foods). Nevertheless, optimal dosing, duration, and timing (particularly in pediatric populations) require further investigation.
5 Conclusion
In summary, our investigation provides evidence that Bifidobacterium might exert protective effects against myopia. These findings suggest that Bifidobacterium supplementation could be a promising strategy for myopia risk reduction. Future research should prioritize the development and evaluation of Bifidobacterium-based interventions to establish their efficacy in myopia control.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.
Ethics statement
Ethical approval was not required for the study involving humans in accordance with the local legislation and institutional requirements. Written informed consent to participate in this study was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and the institutional requirements. The animal study was approved by the Ethical Review Committee for Laboratory Animals of the Shandong Eye Disease Control Institute (AEHSDUTCM-LAEC2024001) and conformed to the ARVO Statement for the Use of Animals in Ophthalmic and Vision Research. Reporting adhered to ARRIVE guidelines to ensure transparency and completeness. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
GT: Conceptualization, Writing – review & editing, Writing – original draft, Methodology. YH: Investigation, Writing – original draft. XG: Writing – original draft, Investigation. XW: Visualization, Validation, Writing – review & editing. JL: Resources, Writing – review & editing, Investigation. JZ: Writing – review & editing, Software. JW: Writing – review & editing, Validation. JS: Project administration, Writing – review & editing, Writing – original draft. HB: Resources, Conceptualization, Project administration, Writing – review & editing, Supervision.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the Key Research and Development Program of Shandong Province (Grant No. 2021LCZX09), the Natural Science Foundation of Shandong Province (Grant Nos. ZR2021LZY045), the Medical and Health Science and Technology Project of Shandong Province (Grant No. 202307021729), the Science and Technology Project of Traditional Chinese Medicine of Shandong Province (Grant No. M-2023075), the Scientific Research Foundation of Shandong University of Traditional Chinese Medicine (General Program, Grant No. KYZK2024M06), and the Clinical Research Special Project of Shandong University of Traditional Chinese Medicine for 2024 (Grant No. LCKY202431).
Acknowledgments
We would like to express our sincere gratitude to the MiBioGen consortium project and the FinnGen database for making all the data needed.
Conflict of interest
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.
Generative AI statement
The author(s) declare that no Gen AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2025.1634120/full#supplementary-material
Footnotes
References
1. Holden, BA, Fricke, TR, Wilson, DA, Jong, M, Naidoo, KS, Sankaridurg, P, et al. Global prevalence of myopia and high myopia and temporal trends from 2000 through 2050. Ophthalmology. (2016) 123:1036–42. doi: 10.1016/j.ophtha.2016.01.006
2. Zhao, X, Lu, X, Yu, L, Zhang, Y, Li, J, Liu, Y, et al. Prevalence of myopia and associated risk factors among key schools in Xi’an, China. BMC Ophthalmol. (2022) 22:519. doi: 10.1186/s12886-022-02735-x
3. Martínez-Albert, N, Bueno-Gimeno, I, and Gené-Sampedro, A. Risk factors for myopia: a review. J Clin Med. (2023) 12:6062. doi: 10.3390/jcm12186062
4. Morgan, IG, French, AN, Ashby, RS, Guo, X, Ding, X, He, M, et al. The epidemics of myopia: Aetiology and prevention. Prog Retin Eye Res. (2018) 62:134–49. doi: 10.1016/j.preteyeres.2017.09.004
5. Haarman, AEG, Enthoven, CA, Tideman, JWL, Tedja, MS, Verhoeven, VJM, and Klaver, CCW. The complications of myopia: a review and Meta-analysis. Invest Ophthalmol Vis Sci. (2020) 61:49. doi: 10.1167/iovs.61.4.49
6. Flaxman, SR, Bourne, RRA, Resnikoff, S, Ackland, P, Braithwaite, T, Cicinelli, MV, et al. Global causes of blindness and distance vision impairment 1990-2020: a systematic review and meta-analysis. Lancet Glob Health. (2017) 5:e1221–34. doi: 10.1016/S2214-109X(17)30393-5
7. Saw, S-M, Matsumura, S, and Hoang, QV. Prevention and Management of Myopia and Myopic Pathology. Invest Ophthalmol Vis Sci. (2019) 60:488–99. doi: 10.1167/iovs.18-25221
8. Chen, Y, Zhou, J, and Wang, L. Role and mechanism of gut microbiota in human disease. Front Cell Infect Microbiol. (2021) 11:625913. doi: 10.3389/fcimb.2021.625913
9. Martino, C, Dilmore, AH, Burcham, ZM, Metcalf, JL, Jeste, D, and Knight, R. Microbiota succession throughout life from the cradle to the grave. Nat Rev Microbiol. (2022) 20:707–20. doi: 10.1038/s41579-022-00768-z
10. Wang, R, Tang, R, Li, B, Ma, X, Schnabl, B, and Tilg, H. Gut microbiome, liver immunology, and liver diseases. Cell Mol Immunol. (2021) 18:4–17. doi: 10.1038/s41423-020-00592-6
11. Liu, Y, Fan, H, Shao, Y, Zhang, J, Zuo, Z, Wang, J, et al. Gut microbiota dysbiosis associated with different types of demyelinating optic neuritis in patients. Mult Scler Relat Disord. (2023) 72:104619. doi: 10.1016/j.msard.2023.104619
12. Mao, D, Tao, B, Sheng, S, Jin, H, Chen, W, Gao, H, et al. Causal effects of gut microbiota on age-related macular degeneration: a Mendelian randomization study. Invest Ophthalmol Vis Sci. (2023) 64:32. doi: 10.1167/iovs.64.12.32
13. Huang, Y, Wang, Z, Ye, B, Ma, JH, Ji, S, Sheng, W, et al. Sodium butyrate ameliorates diabetic retinopathy in mice via the regulation of gut microbiota and related short-chain fatty acids. J Transl Med. (2023) 21:451. doi: 10.1186/s12967-023-04259-4
14. Li, H, Liu, S, Zhang, K, Zhu, X, Dai, J, and Lu, Y. Gut microbiome and plasma metabolome alterations in myopic mice. Front Microbiol. (2023) 14:1251243. doi: 10.3389/fmicb.2023.1251243
15. Wu, Y, Fan, H, Feng, Y, Yang, J, Cen, X, and Li, W. Unveiling the gut microbiota and metabolite profiles in guinea pigs with form deprivation myopia through 16S rRNA gene sequencing and untargeted metabolomics. Heliyon. (2024) 10:e30491. doi: 10.1016/j.heliyon.2024.e30491
16. Omar, WEW, Singh, G, McBain, AJ, Cruickshank, F, and Radhakrishnan, H. Gut microbiota profiles in Myopes and Nonmyopes. Invest Ophthalmol Vis Sci. (2024) 65:2. doi: 10.1167/iovs.65.5.2
17. Chong, RS, Li, H, Cheong, AJY, Fan, Q, Koh, V, Raghavan, L, et al. Mendelian randomization implicates bidirectional association between myopia and primary open-angle Glaucoma or intraocular pressure. Ophthalmology. (2023) 130:394–403. doi: 10.1016/j.ophtha.2022.11.030
18. Qin, Y, Lei, C, Lin, T, Han, X, and Wang, D. Identification of potential drug targets for myopia through Mendelian randomization. Invest Ophthalmol Vis Sci. (2024) 65:13. doi: 10.1167/iovs.65.10.13
19. Mountjoy, E, Davies, NM, Plotnikov, D, Smith, GD, Rodriguez, S, Williams, CE, et al. Education and myopia: assessing the direction of causality by mendelian randomisation. BMJ. (2018) 362:k2932. doi: 10.1136/bmj.k2932
20. Sekula, P, Del Greco, MF, Pattaro, C, and Köttgen, A. Mendelian randomization as an approach to assess causality using observational data. J Am Soc Nephrol. (2016) 27:3253–65. doi: 10.1681/ASN.2016010098
21. Davies, NM, Holmes, MV, and Davey Smith, G. Reading Mendelian randomisation studies: a guide, glossary, and checklist for clinicians. BMJ. (2018) 362:k601. doi: 10.1136/bmj.k601
22. Jeong, J, Yun, K, Mun, S, Chung, W-H, Choi, S-Y, Nam, Y, et al. The effect of taxonomic classification by full-length 16S rRNA sequencing with a synthetic long-read technology. Sci Rep. (2021) 11:1727. doi: 10.1038/s41598-020-80826-9
23. Boef, AGC, Dekkers, OM, and le Cessie, S. Mendelian randomization studies: a review of the approaches used and the quality of reporting. Int J Epidemiol. (2015) 44:496–511. doi: 10.1093/ije/dyv071
24. Kurilshikov, A, Medina-Gomez, C, Bacigalupe, R, Radjabzadeh, D, Wang, J, Demirkan, A, et al. Large-scale association analyses identify host factors influencing human gut microbiome composition. Nat Genet. (2021) 53:156–65. doi: 10.1038/s41588-020-00763-1
25. Verbanck, M, Chen, C-Y, Neale, B, and Do, R. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet. (2018) 50:693–8. doi: 10.1038/s41588-018-0099-7
26. Greco, MFD, Minelli, C, Sheehan, NA, and Thompson, JR. Detecting pleiotropy in Mendelian randomisation studies with summary data and a continuous outcome. Stat Med. (2015) 34:2926–40. doi: 10.1002/sim.6522
27. Bowden, J, Del Greco, MF, Minelli, C, Zhao, Q, Lawlor, DA, Sheehan, NA, et al. Improving the accuracy of two-sample summary-data Mendelian randomization: moving beyond the NOME assumption. Int J Epidemiol. (2019) 48:728–42. doi: 10.1093/ije/dyy258
28. Burgess, S, and Thompson, SG. Interpreting findings from Mendelian randomization using the MR-egger method. Eur J Epidemiol. (2017) 32:377–89. doi: 10.1007/s10654-017-0255-x
29. Xiao, H, Fan, Z-Y, Tian, X-D, and Xu, Y-C. Comparison of form-deprived myopia and lens-induced myopia in guinea pigs. Int J Ophthalmol. (2014) 7:245–50. doi: 10.3980/j.issn.2222-3959.2014.02.10
30. Prodan, A, Tremaroli, V, Brolin, H, Zwinderman, AH, Nieuwdorp, M, and Levin, E. Comparing bioinformatic pipelines for microbial 16S rRNA amplicon sequencing. PLoS One. (2020) 15:e0227434. doi: 10.1371/journal.pone.0227434
31. Glöckner, FO, Yilmaz, P, Quast, C, Gerken, J, Beccati, A, Ciuprina, A, et al. 25 years of serving the community with ribosomal RNA gene reference databases and tools. J Biotechnol. (2017) 261:169–76. doi: 10.1016/j.jbiotec.2017.06.1198
32. Bolyen, E, Rideout, JR, Dillon, MR, Bokulich, NA, Abnet, CC, Al-Ghalith, GA, et al. Author correction: reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. (2019) 37:1091. doi: 10.1038/s41587-019-0252-6
33. Burgess, S, and Thompson, SGCRP CHD Genetics Collaboration. Avoiding bias from weak instruments in Mendelian randomization studies. Int J Epidemiol. (2011) 40:755–64. doi: 10.1093/ije/dyr036
34. Yun, S-W, Son, Y-H, Lee, D-Y, Shin, Y-J, Han, MJ, and Kim, D-H. Lactobacillus plantarum and Bifidobacterium bifidum alleviate dry eye in mice with exorbital lacrimal gland excision by modulating gut inflammation and microbiota. Food Funct. (2021) 12:2489–97. doi: 10.1039/d0fo02984j
35. Zhao, X, Xu, M, Zhao, Z, Wang, Y, Liu, Y, Zhang, T, et al. Bifidobacterium promotes retinal ganglion cell survival by regulating the balance of retinal glial cells. CNS Neurosci Ther. (2023) 29:146–60. doi: 10.1111/cns.14165
36. Markowiak-Kopeć, P, and Śliżewska, K. The effect of probiotics on the production of short-chain fatty acids by human intestinal microbiome. Nutrients. (2020) 12:1107. doi: 10.3390/nu12041107
37. Yao, Y, Cai, X, Fei, W, Ye, Y, Zhao, M, and Zheng, C. The role of short-chain fatty acids in immunity, inflammation and metabolism. Crit Rev Food Sci Nutr. (2022) 62:1–12. doi: 10.1080/10408398.2020.1854675
38. Ang, Z, and Ding, JL. GPR41 and GPR43 in obesity and inflammation - protective or causative? Front Immunol. (2016) 7:28. doi: 10.3389/fimmu.2016.00028
39. Wang, J, and Friedman, E. Downregulation of p53 by sustained JNK activation during apoptosis. Mol Carcinog. (2000) 29:179–88. doi: 10.1002/1098-2744(200011)29:3<179::aid-mc7>3.0.co;2-k
40. Ku, H, Chen, JJ-Y, Hu, M, Tien, P-T, Lin, H-J, Xu, G, et al. Myopia development in tree shrew is associated with chronic inflammatory reactions. Curr Issues Mol Biol. (2022) 44:4303–13. doi: 10.3390/cimb44090296
41. Masjedi, A, Hashemi, V, Hojjat-Farsangi, M, Ghalamfarsa, G, Azizi, G, Yousefi, M, et al. The significant role of interleukin-6 and its signaling pathway in the immunopathogenesis and treatment of breast cancer. Biomed Pharmacother. (2018) 108:1415–24. doi: 10.1016/j.biopha.2018.09.177
42. Kaur, S, Bansal, Y, Kumar, R, and Bansal, G. A panoramic review of IL-6: structure, pathophysiological roles and inhibitors. Bioorg Med Chem. (2020) 28:115327. doi: 10.1016/j.bmc.2020.115327
43. Zahir-Jouzdani, F, Atyabi, F, and Mojtabavi, N. Interleukin-6 participation in pathology of ocular diseases. Pathophysiology. (2017) 24:123–31. doi: 10.1016/j.pathophys.2017.05.005
44. Hunter, CA, and Jones, SA. IL-6 as a keystone cytokine in health and disease. Nat Immunol. (2015) 16:448–57. doi: 10.1038/ni.3153
45. Li, X, Römer, G, Kerindongo, RP, Hermanides, J, Albrecht, M, Hollmann, MW, et al. Sodium glucose co-transporter 2 inhibitors ameliorate endothelium barrier dysfunction induced by cyclic stretch through inhibition of reactive oxygen species. Int J Mol Sci. (2021) 22:6044. doi: 10.3390/ijms22116044
46. Yu, F-J, Lam, TC, Sze, AY-H, Li, K-K, Chun, RK-M, and Shan, S-W. Alteration of retinal metabolism and oxidative stress may implicate myopic eye growth: evidence from discovery and targeted proteomics in an animal model. J Proteome. (2020) 221:103684. doi: 10.1016/j.jprot.2020.103684
47. Jonas, JB, Ohno-Matsui, K, and Panda-Jonas, S. Myopia: anatomic changes and consequences for its etiology. Asia Pac J Ophthalmol (Phila). (2019) 8:355–9. doi: 10.1097/01.APO.0000578944.25956.8b
48. Zeng, Q, Qi, Z, He, X, Luo, C, Wen, J, Wei, J, et al. Bifidobacterium pseudocatenulatum NCU-08 ameliorated senescence via modulation of the AMPK/Sirt1 signaling pathway and gut microbiota in mice. Food Funct. (2024) 15:4095–108. doi: 10.1039/d3fo04575g
49. Verbarg, S, Rheims, H, Emus, S, Frühling, A, Kroppenstedt, RM, Stackebrandt, E, et al. Erysipelothrix inopinata sp. nov., isolated in the course of sterile filtration of vegetable peptone broth, and description of Erysipelotrichaceae fam. Nov. Int J Syst Evol Microbiol. (2004) 54:221–5. doi: 10.1099/ijs.0.02898-0
50. Zschaler, J, Schlorke, D, and Arnhold, J. Differences in innate immune response between man and mouse. Crit Rev Immunol. (2014) 34:433–54. doi: 10.1615/CritRevImmunol.2014011600
51. Nguyen, TLA, Vieira-Silva, S, Liston, A, and Raes, J. How informative is the mouse for human gut microbiota research? Dis Model Mech. (2015) 8:1–16. doi: 10.1242/dmm.017400
52. Palm, NW, de Zoete, MR, Cullen, TW, Barry, NA, Stefanowski, J, Hao, L, et al. Immunoglobulin a coating identifies colitogenic bacteria in inflammatory bowel disease. Cell. (2014) 158:1000–10. doi: 10.1016/j.cell.2014.08.006
53. Zhao, Y, Wu, J, Li, JV, Zhou, N-Y, Tang, H, and Wang, Y. Gut microbiota composition modifies fecal metabolic profiles in mice. J Proteome Res. (2013) 12:2987–99. doi: 10.1021/pr400263n
54. Ng, SK, Wood, JPM, Chidlow, G, Han, G, Kittipassorn, T, Peet, DJ, et al. Cancer-like metabolism of the mammalian retina. Clin Experiment Ophthalmol. (2015) 43:367–76. doi: 10.1111/ceo.12462
55. Ito, YA, and Di Polo, A. Mitochondrial dynamics, transport, and quality control: a bottleneck for retinal ganglion cell viability in optic neuropathies. Mitochondrion. (2017) 36:186–92. doi: 10.1016/j.mito.2017.08.014
56. Wu, S, Guo, D, Wei, H, Yin, X, Zhang, L, Guo, B, et al. Disrupted potassium ion homeostasis in ciliary muscle in negative lens-induced myopia in Guinea pigs. Arch Biochem Biophys. (2020) 688:108403. doi: 10.1016/j.abb.2020.108403
57. Zhong, W, Lan, C, Gu, Z, Tan, Q, Xiang, X, Zhou, H, et al. The mechanosensitive Piezo1 channel mediates Mechanochemical transmission in myopic eyes. Invest Ophthalmol Vis Sci. (2023) 64:1. doi: 10.1167/iovs.64.7.1
58. Dutheil, F, Oueslati, T, Delamarre, L, Castanon, J, Maurin, C, Chiambaretta, F, et al. Myopia and near work: a systematic review and Meta-analysis. Int J Environ Res Public Health. (2023) 20:875. doi: 10.3390/ijerph20010875
59. Rose, KA, Morgan, IG, Ip, J, Kifley, A, Huynh, S, Smith, W, et al. Outdoor activity reduces the prevalence of myopia in children. Ophthalmology. (2008) 115:1279–85. doi: 10.1016/j.ophtha.2007.12.019
Keywords: myopia, gut microbiota, Mendelian randomization, full-length 16S rRNA gene sequencing, Bifidobacterium
Citation: Tang G, Han Y, Gong X, Wang X, Li J, Zhang J, Wang J, Song J and Bi H (2025) Potential protective role of Bifidobacteria in myopia prevention: evidence from full-length 16S rRNA sequencing and bidirectional Mendelian randomization analysis. Front. Med. 12:1634120. doi: 10.3389/fmed.2025.1634120
Edited by:
Weihua Yang, Southern Medical University, ChinaReviewed by:
Miaoran Zhang, Jilin University, ChinaTainan Lin, Fujian Provincial Governmental Hospital, China
Copyright © 2025 Tang, Han, Gong, Wang, Li, Zhang, Wang, Song and Bi. 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) and the copyright owner(s) 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: Hongsheng Bi, aG9uZ3NoZW5nYmkxQDE2My5jb20=; Jike Song, ZWR1c2prQDE2My5jb20=