Adaptive strategies and ecological roles of phages in habitats under physicochemical stress.

Bacteriophages (phages) play a vital role in ecosystem functions by influencing the composition, genetic exchange, metabolism, and environmental adaptation of microbial communities. With recent advances in sequencing technologies and bioinformatics, our understanding of the ecology and evolution of phages in stressful environments has substantially expanded. Here, we review the impact of physicochemical environmental stress on the physiological state and community dynamics of phages, the adaptive strategies that phages employ to cope with environmental stress, and the ecological effects of phage-host interactions in stressful environments. Specifically, we highlight the contributions of phages to the adaptive evolution and functioning of microbiomes and suggest that phages and their hosts can maintain a mutualistic relationship in response to environmental stress. In addition, we discuss the ecological consequences caused by phages in stressful environments, encompassing biogeochemical cycling. Overall, this review advances an understanding of phage ecology in stressful environments, which could inform phage-based strategies to improve microbiome performance and ecosystem resilience and resistance in natural and engineering systems.

COBRA improves the completeness and contiguity of viral genomes assembled from metagenomes.

Viruses are often studied using metagenome-assembled sequences, but genome incompleteness hampers comprehensive and accurate analyses. Contig Overlap Based Re-Assembly (COBRA) resolves assembly breakpoints based on the de Bruijn graph and joins contigs. Here we benchmarked COBRA using ocean and soil viral datasets. COBRA accurately joined the assembled sequences and achieved notably higher genome accuracy than binning tools. From 231 published freshwater metagenomes, we obtained 7,334 bacteriophage clusters, ~83% of which represent new phage species. Notably, ~70% of these were circular, compared with 34% before COBRA analyses. We expanded sampling of huge phages (≥200 kbp), the largest of which was curated to completion (717 kbp). Improved phage genomes from Rotsee Lake provided context for metatranscriptomic data and indicated the in situ activity of huge phages, whiB-encoding phages and cysC- and cysH-encoding phages. COBRA improves viral genome assembly contiguity and completeness, thus the accuracy and reliability of analyses of gene content, diversity and evolution.

CasPEDIA Database: a functional classification system for class 2 CRISPR-Cas enzymes.

CRISPR-Cas enzymes enable RNA-guided bacterial immunity and are widely used for biotechnological applications including genome editing. In particular, the Class 2 CRISPR-associated enzymes (Cas9, Cas12 and Cas13 families), have been deployed for numerous research, clinical and agricultural applications. However, the immense genetic and biochemical diversity of these proteins in the public domain poses a barrier for researchers seeking to leverage their activities. We present CasPEDIA (http://caspedia.org), the Cas Protein Effector Database of Information and Assessment, a curated encyclopedia that integrates enzymatic classification for hundreds of different Cas enzymes across 27 phylogenetic groups spanning the Cas9, Cas12 and Cas13 families, as well as evolutionarily related IscB and TnpB proteins. All enzymes in CasPEDIA were annotated with a standard workflow based on their primary nuclease activity, target requirements and guide-RNA design constraints. Our functional classification scheme, CasID, is described alongside current phylogenetic classification, allowing users to search related orthologs by enzymatic function and sequence similarity. CasPEDIA is a comprehensive data portal that summarizes and contextualizes enzymatic properties of widely used Cas enzymes, equipping users with valuable resources to foster biotechnological development. CasPEDIA complements phylogenetic Cas nomenclature and enables researchers to leverage the multi-faceted nucleic-acid targeting rules of diverse Class 2 Cas enzymes.

Infant gut DNA bacteriophage strain persistence during the first 3 years of life.

Bacteriophages are key components of gut microbiomes, yet the phage colonization process in the infant gut remains uncertain. Here, we establish a large phage sequence database and use strain-resolved analyses to investigate DNA phage succession in infants throughout the first 3 years of life. Analysis of 819 fecal metagenomes collected from 28 full-term and 24 preterm infants and their mothers revealed that early-life phageome richness increases over time and reaches adult-like complexity by age 3. Approximately 9% of early phage colonizers, which are mostly maternally transmitted and infect Bacteroides, persist for 3 years and are more prevalent in full-term than in preterm infants. Although rare, phages with stop codon reassignment are more likely to persist than non-recoded phages and generally display an increase in in-frame reassigned stop codons over 3 years. Overall, maternal seeding, stop codon reassignment, host CRISPR-Cas locus prevalence, and diverse phage populations contribute to stable viral colonization.

Long-Term Axial Length Shortening in Myopic Orthokeratology: Incident Probability, Time Course, and Influencing Factors.

Purpose: Long-term axial length (AL) shortening in myopia is uncommon but noteworthy. Current understanding on the condition is limited due to difficulties in case collection. The study reported percentage, probability, and time course of long-term AL shortening in myopic orthokeratology based on a large database. Methods: This study reviewed 142,091 medical records from 29,825 subjects in a single-hospital orthokeratology database that were collected over 10 years. Long-term AL shortening was defined as a change in AL of -0.1 mm or less at any follow-up beyond 1 year. Incident probability was calculated based on multivariate logistic regression. Time course was estimated using mixed-effect regression model. Results: A total of 10,093 subjects (mean initial age, 11.70 ± 2.52 years; 58.8% female) with 80,778 visits were included. The number of subjects experienced long-term AL shortening was 1,662 (16.47%; 95% confidence interval, 15.75%-17.21%). Initial age showed significant impact on the incident occurrence (OR, 1.37; 95% confidence interval, 1.34-1.40; P < 0.001). The estimated probability of AL shortening was approximately 2% for subjects with initial age of 6 years and 50% for those aged 18. Among the 1662 AL shortening cases, the median magnitude of the maximum AL reduction was 0.19 mm. The shortening process mostly occurred within the initial 2 years. Subject characteristics had limited associations with the shortening rate. Conclusions: Long-term AL shortening is possible in subjects receiving myopic orthokeratology. Although age notably affect the incident probability, the time course seems to not vary significantly.

Host translation machinery is not a barrier to phages that interact with both CPR and non-CPR bacteria.

IMPORTANCE: Here, we profiled putative phages of Saccharibacteria, which are of particular importance as Saccharibacteria influence some human oral diseases. We additionally profiled putative phages of Gracilibacteria and Absconditabacteria, two Candidate Phyla Radiation (CPR) lineages of interest given their use of an alternative genetic code. Among the phages identified in this study, some are targeted by spacers from both CPR and non-CPR bacteria and others by both bacteria that use the standard genetic code as well as bacteria that use an alternative genetic code. These findings represent new insights into possible phage replication strategies and have relevance for phage therapies that seek to manipulate microbiomes containing CPR bacteria.

Eukaryotic RNA-guided endonucleases evolved from a unique clade of bacterial enzymes.

RNA-guided endonucleases form the crux of diverse biological processes and technologies, including adaptive immunity, transposition, and genome editing. Some of these enzymes are components of insertion sequences (IS) in the IS200/IS605 and IS607 transposon families. Both IS families encode a TnpA transposase and a TnpB nuclease, an RNA-guided enzyme ancestral to CRISPR-Cas12s. In eukaryotes, TnpB homologs occur as two distinct types, Fanzor1s and Fanzor2s. We analyzed the evolutionary relationships between prokaryotic TnpBs and eukaryotic Fanzors, which revealed that both Fanzor1s and Fanzor2s stem from a single lineage of IS607 TnpBs with unusual active site arrangement. The widespread nature of Fanzors implies that the properties of this particular lineage of IS607 TnpBs were particularly suited to adaptation in eukaryotes. Biochemical analysis of an IS607 TnpB and Fanzor1s revealed common strategies employed by TnpBs and Fanzors to co-evolve with their cognate transposases. Collectively, our results provide a new model of sequential evolution from IS607 TnpBs to Fanzor2s, and Fanzor2s to Fanzor1s that details how genes of prokaryotic origin evolve to give rise to new protein families in eukaryotes.

Deep-branching evolutionary intermediates reveal structural origins of form I rubisco.

The enzyme rubisco (ribulose-1,5-bisphosphate carboxylase/oxygenase) catalyzes the majority of biological carbon fixation on Earth. Although the vast majority of rubiscos across the tree of life assemble as homo-oligomers, the globally predominant form I enzyme-found in plants, algae, and cyanobacteria-forms a unique hetero-oligomeric complex. The recent discovery of a homo-oligomeric sister group to form I rubisco (named form I') has filled a key gap in our understanding of the enigmatic origins of the form I clade. However, to elucidate the series of molecular events leading to the evolution of form I rubisco, we must examine more distantly related sibling clades to contextualize the molecular features distinguishing form I and form I' rubiscos. Here, we present a comparative structural study retracing the evolutionary history of rubisco that reveals a complex structural trajectory leading to the ultimate hetero-oligomerization of the form I clade. We structurally characterize the oligomeric states of deep-branching form Iα and I'' rubiscos recently discovered from metagenomes, which represent key evolutionary intermediates preceding the form I clade. We further solve the structure of form I'' rubisco, revealing the molecular determinants that likely primed the enzyme core for the transition from a homo-oligomer to a hetero-oligomer. Our findings yield new insight into the evolutionary trajectory underpinning the adoption and entrenchment of the prevalent assembly of form I rubisco, providing additional context when viewing the enzyme family through the broader lens of protein evolution.

Myopic progression in school-aged children with moderate intermittent exotropia.

Objective: It is still controversial whether intermittent exotropia (IXT) affects myopic progression during the critical period of visual development. This study retrospectively analyzed the long-term myopic changes and the impact of IXT surgery on myopic progression in school-aged children with moderate IXT. Methods: The medical records of 65 children from 5 to 13 years old with or without IXT between 2015 and 2021 were retrospectively reviewed. Patients whose spherical equivalent refraction (SER) were less than -3.00 diopter (D) were included and divided into three groups: IXT surgery group (Group A), which comprised 22 IXT patients who received IXT surgery, IXT observation group (Group B), which comprises 19 IXT patients who only received long-term observational follow-up; and normal control group (Group C), which comprised 24 normal controls without IXT. The main outcome measurement was the rate of myopic progression, which was defined as the mean myopic shift in SER per year. Results: The 3- and 5-years long-term follow-up rates of myopic progression were -0.47 ± 0.28 D per year and -0.48 ± 0.23 D per year respectively in Group B, and those were significantly slower than that in Group C (-0.73 ± 0.32 D per year and -0.76 ± 0.19 D per year respectively). However, there was no significant difference in the rate of myopic progression between Group A and B or between Group A and C. Conclusion: Moderate IXT may have lower rate of myopic progression in school-aged children. Whether IXT surgery influence the rate of myopic progression still needs further study.

Eukaryotic RNA-guided endonucleases evolved from a unique clade of bacterial enzymes.

RNA-guided endonucleases form the crux of diverse biological processes and technologies, including adaptive immunity, transposition, and genome editing. Some of these enzymes are components of insertion sequences (IS) in the IS200/IS605 and IS607 transposon families. Both IS families encode a TnpA transposase and TnpB nuclease, an RNA-guided enzyme ancestral to CRISPR-Cas12. In eukaryotes and their viruses, TnpB homologs occur as two distinct types, Fanzor1 and Fanzor2. We analyzed the evolutionary relationships between prokaryotic TnpBs and eukaryotic Fanzors, revealing that a clade of IS607 TnpBs with unusual active site arrangement found primarily in Cyanobacteriota likely gave rise to both types of Fanzors. The wide-spread nature of Fanzors imply that the properties of this particular group of IS607 TnpBs were particularly suited to adaptation and evolution in eukaryotes and their viruses. Experimental characterization of a prokaryotic IS607 TnpB and virally encoded Fanzor1s uncovered features that may have fostered coevolution between TnpBs/Fanzors and their cognate transposases. Our results provide insight into the evolutionary origins of a ubiquitous family of RNA-guided proteins that shows remarkable conservation across domains of life.

Inference and reconstruction of the heimdallarchaeial ancestry of eukaryotes.

In the ongoing debates about eukaryogenesis-the series of evolutionary events leading to the emergence of the eukaryotic cell from prokaryotic ancestors-members of the Asgard archaea play a key part as the closest archaeal relatives of eukaryotes1. However, the nature and phylogenetic identity of the last common ancestor of Asgard archaea and eukaryotes remain unresolved2-4. Here we analyse distinct phylogenetic marker datasets of an expanded genomic sampling of Asgard archaea and evaluate competing evolutionary scenarios using state-of-the-art phylogenomic approaches. We find that eukaryotes are placed, with high confidence, as a well-nested clade within Asgard archaea and as a sister lineage to Hodarchaeales, a newly proposed order within Heimdallarchaeia. Using sophisticated gene tree and species tree reconciliation approaches, we show that analogous to the evolution of eukaryotic genomes, genome evolution in Asgard archaea involved significantly more gene duplication and fewer gene loss events compared with other archaea. Finally, we infer that the last common ancestor of Asgard archaea was probably a thermophilic chemolithotroph and that the lineage from which eukaryotes evolved adapted to mesophilic conditions and acquired the genetic potential to support a heterotrophic lifestyle. Our work provides key insights into the prokaryote-to-eukaryote transition and a platform for better understanding the emergence of cellular complexity in eukaryotic cells.

O2 partitioning of sulfur oxidizing bacteria drives acidity and thiosulfate distributions in mining waters.

The acidification of water in mining areas is a global environmental issue primarily catalyzed by sulfur-oxidizing bacteria (SOB). Little is known about microbial sulfur cycling in circumneutral pH mine tailing impoundment waters. Here we investigate biological sulfur oxidation over four years in a mine tailings impoundment water cap, integrating aqueous sulfur geochemistry, genome-resolved metagenomics and metatranscriptomics. The microbial community is consistently dominated by neutrophilic, chemolithoautotrophic SOB (relative abundances of ~76% in 2015, ~55% in 2016/2017 and ~60% in 2018). Results reveal two SOB strategies alternately dominate across the four years, influencing acid generation and sulfur speciation. Under oxic conditions, novel Halothiobacillus drive lower pH conditions (as low as 4.3) and lower [S2O32-] via the complete Sox pathway coupled to O2. Under anoxic conditions, Thiobacillus spp. dominate in activity, via the incomplete Sox and rDSR pathways coupled to NO3-, resulting in higher [S2O32-] and no net significant acidity generation. This study provides genomic evidence explaining acidity generation and thiosulfate accumulation patterns in a circumneutral mine tailing impoundment and has significant environmental applications in preventing the discharge of sulfur compounds that can impact downstream environments. These insights illuminate opportunities for in situ biotreatment of reduced sulfur compounds and prediction of acidification events using gene-based monitoring and in situ RNA detection.

Alum Addition Triggers Hypoxia in an Engineered Pit Lake.

Here, we examine the geobiological response to a whole-lake alum (aluminum sulfate) treatment (2016) of Base Mine Lake (BML), the first pilot-scale pit lake established in the Alberta oil sands region. The rationale for trialing this management amendment was based on its successful use to reduce internal phosphorus loading to eutrophying lakes. Modest increases in water cap epilimnetic oxygen concentrations, associated with increased Secchi depths and chlorophyll-a concentrations, were co-incident with anoxic waters immediately above the fluid fine tailings (FFT) layer post alum. Decreased water cap nitrate and detectable sulfide concentrations, as well as increased hypolimnetic phospholipid fatty acid abundances, signaled greater anaerobic heterotrophic activity. Shifts in microbial community to groups associated with greater organic carbon degradation (i.e., SAR11-LD12 subclade) and the SRB group Desulfuromonodales emerged post alum and the loss of specialist groups associated with carbon-limited, ammonia-rich restricted niches (i.e., MBAE14) also occurred. Alum treatment resulted in additional oxygen consumption associated with increased autochthonous carbon production, watercap anoxia and sulfide generation, which further exacerbate oxygen consumption associated with on-going FFT mobilized reductants. The results illustrate the importance of understanding the broader biogeochemical implications of adaptive management interventions to avoid unanticipated outcomes that pose greater risks and improve tailings reclamation for oil sands operations and, more broadly, the global mining sector.

Long-Term Incubation of Lake Water Enables Genomic Sampling of Consortia Involving Planctomycetes and Candidate Phyla Radiation Bacteria.

Microbial communities in lakes can profoundly impact biogeochemical processes through their individual activities and collective interactions. However, the complexity of these communities poses challenges, particularly for studying rare organisms such as Candidate Phyla Radiation bacteria (CPR) and enigmatic entities such as aster-like nanoparticles (ALNs). Here, a reactor was inoculated with water from Lake Fargette, France, and maintained under dark conditions at 4°C for 31 months and enriched for ALNs, diverse Planctomycetes, and CPR bacteria. We reconstructed draft genomes and predicted metabolic traits for 12 diverse Planctomycetes and 9 CPR bacteria, some of which are likely representatives of undescribed families or genera. One CPR genome representing the little-studied lineage "Candidatus Peribacter" was curated to completion (1.239 Mbp) and unexpectedly encodes the full gluconeogenesis pathway. Metatranscriptomic data indicate that some planctomycetes and CPR bacteria were active under the culture conditions, accounting for ∼30% and ∼1% of RNA reads mapping to the genome set, respectively. We also reconstructed genomes and obtained transmission electron microscope images for numerous viruses, including one with a >300-kbp genome and several predicted to infect Planctomycetes. Together, our analyses suggest that freshwater Planctomycetes are central players in a subsystem that includes ALNs, symbiotic CPR bacteria, and viruses. IMPORTANCE Laboratory incubations of natural microbial communities can aid in the study of member organisms and their networks of interaction. This is particularly important for understudied lineages for which key elements of basic biology are still emerging. Using genomics and microscopy, we found that members of the bacterial lineage Planctomycetes may be central players in a subset of a freshwater lake microbiome that includes other bacteria, archaea, viruses, and mysterious entities, called aster-like nanoparticles (ALNs), whose origin is unknown. Our results help constrain the possible origins of ALNs and provide insight into possible interactions within a complex lake ecosystem.

Detection and enumeration of Lak megaphages in microbiome samples by endpoint and quantitative PCR.

Lak megaphages are prevalent across diverse gut microbiomes and may potentially impact animal and human health through lysis of Prevotella. Given their large genome size (up to 660 kbp), Lak megaphages are difficult to culture, and their identification relies on molecular techniques. Here, we present optimized protocols for identifying Lak phages in various microbiome samples, including procedures for DNA extraction, followed by detection and quantification of genes encoding Lak structural proteins using diagnostic endpoint and SYBR green-based quantitative PCR, respectively. For complete details on the use and execution of this protocol, please refer to Crisci et al., (2021).

Phage-encoded ribosomal protein S21 expression is linked to late-stage phage replication.

The ribosomal protein S21 (bS21) gene has been detected in diverse viruses with a large range of genome sizes, yet its in situ expression and potential significance have not been investigated. Here, we report five closely related clades of bacteriophages (phages) represented by 47 genomes (8 curated to completion and up to 331 kbp in length) that encode a bS21 gene. The bS21 gene is on the reverse strand within a conserved region that encodes the large terminase, major capsid protein, prohead protease, portal vertex proteins, and some hypothetical proteins. Based on CRISPR spacer targeting, the predominance of bacterial taxonomic affiliations of phage genes with those from Bacteroidetes, and the high sequence similarity of the phage bS21 genes and those from Bacteroidetes classes of Flavobacteriia, Cytophagia and Saprospiria, these phages are predicted to infect diverse Bacteroidetes species that inhabit a range of depths in freshwater lakes. Thus, bS21 phages have the potential to impact microbial community composition and carbon turnover in lake ecosystems. The transcriptionally active bS21-encoding phages were likely in the late stage of replication when collected, as core structural genes and bS21 were highly expressed. Thus, our analyses suggest that the phage bS21, which is involved in translation initiation, substitutes into the Bacteroidetes ribosomes and selects preferentially for phage transcripts during the late-stage replication when large-scale phage protein production is required for assembly of phage particles.

Closely related Lak megaphages replicate in the microbiomes of diverse animals.

Lak phages with alternatively coded ∼540 kbp genomes were recently reported to replicate in Prevotella in microbiomes of humans that consume a non-Western diet, baboons, and pigs. Here, we explore Lak phage diversity and broader distribution using diagnostic polymerase chain reaction and genome-resolved metagenomics. Lak phages were detected in 13 animal types, including reptiles, and are particularly prevalent in pigs. Tracking Lak through the pig gastrointestinal tract revealed significant enrichment in the hindgut compared to the foregut. We reconstructed 34 new Lak genomes, including six curated complete genomes, all of which are alternatively coded. An anomalously large (∼660 kbp) complete genome reconstructed for the most deeply branched Lak from a horse microbiome is also alternatively coded. From the Lak genomes, we identified proteins associated with specific animal species; notably, most have no functional predictions. The presence of closely related Lak phages in diverse animals indicates facile distribution coupled to host-specific adaptation.

Brockarchaeota, a novel archaeal phylum with unique and versatile carbon cycling pathways.

Geothermal environments, such as hot springs and hydrothermal vents, are hotspots for carbon cycling and contain many poorly described microbial taxa. Here, we reconstructed 15 archaeal metagenome-assembled genomes (MAGs) from terrestrial hot spring sediments in China and deep-sea hydrothermal vent sediments in Guaymas Basin, Gulf of California. Phylogenetic analyses of these MAGs indicate that they form a distinct group within the TACK superphylum, and thus we propose their classification as a new phylum, 'Brockarchaeota', named after Thomas Brock for his seminal research in hot springs. Based on the MAG sequence information, we infer that some Brockarchaeota are uniquely capable of mediating non-methanogenic anaerobic methylotrophy, via the tetrahydrofolate methyl branch of the Wood-Ljungdahl pathway and reductive glycine pathway. The hydrothermal vent genotypes appear to be obligate fermenters of plant-derived polysaccharides that rely mostly on substrate-level phosphorylation, as they seem to lack most respiratory complexes. In contrast, hot spring lineages have alternate pathways to increase their ATP yield, including anaerobic methylotrophy of methanol and trimethylamine, and potentially use geothermally derived mercury, arsenic, or hydrogen. Their broad distribution and their apparent anaerobic metabolic versatility indicate that Brockarchaeota may occupy previously overlooked roles in anaerobic carbon cycling.

Association of refraction and ocular biometry in highly myopic eyes.

Clinical relevance: Ocular biometry is key to understanding the determinants of ocular development and pathology changes, especially for the thriving myopic population in Asia. Investigating biometric data in highly myopic eyes within a wide age spectrum is therefore of high importance.Background: To report the magnitude of change in spherical equivalent for each unit of change in the ocular biometry parameters in a highly myopic population in China.Methods: Highly myopic patients aged 7 to 70 years were recruited from the Zhongshan Ophthalmic Center, China. Each patient had a cycloplegic refraction and a measurement of ocular biometry.Results: Data from 823 right eyes were available for analysis, with a mean age of 22.7 years and a median spherical equivalent of -8.88 D. Axial length and lens thickness increased with age, while anterior chamber depth (ACD) decreased in older subjects. There was a significant trend of increasing axial length, lens thickness, vitreous chamber depth (VCD) and decreasing ACD and calculated lens power over spherical equivalent quartiles (all p < 0.001). The univariate linear regression models showed that 1-D change in refraction equalled to a 0.33- to 0.34-mm increase of axial length, and a 0.32 to 0.33-mm increase of VCD in highly myopic eyes. Among the three components of axial length, lens thickness was associated with myopia shift in the groups of 7-18 years and 19-39 years (both p < 0.001), and VCD was significant in all groups (all p < 0.001), while ACD was not significant in any age group.Conclusion: The associations between refraction and axial length were consistent in children, young adults and the elderly with high myopia. Lens thickening with a higher degree of myopia appeared at a very early age, and vitreous chamber depth remained to be a prominent factor of refraction change in highly myopic eyes throughout seven to 70 years of age.

Development and validation of a prediction model for axial length elongation in myopic children treated with overnight orthokeratology.

PURPOSE: To develop and validate a standardized prediction model aiming at 1-year axial length elongation and to guide the orthokeratology lens practice. METHODS: This retrospective study was based on medical records of myopic children treated with orthokeratology. Individuals aged 8-15 years (n = 1261) were included and divided into the primary cohort (n = 757) and validation cohort (n = 504). Feature selection was primarily performed to sort out influential predictors by high-throughput extraction. Then, the prediction model was developed using multivariable linear regression analysis completed by backward stepwise selection. Finally, the validation of the prediction model was performed by evaluation metrics (mean-square error, root-mean-square error, mean absolute error and R ad 2 ). RESULTS: No significant difference was found between primary and validation cohort (all p > 0.05). After the feature selection, the crude model was adjusted by demographic information in multivariable linear regression analysis, and five final predictors were identified (all p < 0.01). The interaction effect of age with 1-month change of zone-3 mm flat K was detected (p < 0.01); hence, two final prediction models were developed based on two age subgroups. The validation proved an acceptable performance. CONCLUSION: An effective multivariable prediction model aiming at 1-year axial length elongation was developed and validated. It can potentially help clinicians to predict orthokeratology efficacy and make valid adjustments. The influential variables revealed in this model can also provide designers directions to optimize the design of lens to improve the efficacy of myopia control.