Analysis of genetic networks regulating refractive eye development in collaborative cross progenitor strain mice reveals new genes and pathways underlying human myopia

Abstract: Background: Population studies suggest that genetic factors play an important role in refractive error development; however, the precise role of genetic background and the composition of the signaling pathways underlying refractive eye development remain poorly understood. Methods: Here, we analyzed normal refractive development and susceptibility to form-deprivation myopia in the eight progenitor mouse strains of the Collaborative Cross (CC). We used RNA-seq to analyze gene expression in the retinae of these mice and reconstruct genetic networks and signaling pathways underlying refractive eye development. We also utilized genome-wide gene-based association analysis to identify mouse genes and pathways associated with myopia in humans. Results: Genetic background strongly influenced both baseline refractive development and susceptibility to environmentally-induced myopia. Baseline refractive errors ranged from − 21.2 diopters (D) in 129S1/svlmj mice to + 22.0 D in CAST/EiJ mice and represented a continuous distribution typical of a quantitative genetic trait. The extent of induced form-deprivation myopia ranged from − 5.6 D in NZO/HILtJ mice to − 20.0 D in CAST/EiJ mice and also followed a continuous distribution. Whole-genome (RNA-seq) gene expression profiling in retinae from CC progenitor strains identified genes whose expression level correlated with either baseline refractive error or susceptibility to myopia. Expression levels of 2,302 genes correlated with the baseline refractive state of the eye, whereas 1,917 genes correlated with susceptibility to induced myopia. Genome-wide gene-based association analysis in the CREAM and UK Biobank human cohorts revealed that 985 of the above genes were associated with myopia in humans, including 847 genes which were implicated in the development of human myopia for the first time. Although the gene sets controlling baseline refractive development and those regulating susceptibility to myopia overlapped, these two processes appeared to be controlled by largely distinct sets of genes. Conclusions: Comparison with data for other animal models of myopia revealed that the genes identified in this study comprise a well-defined set of retinal signaling pathways, which are highly conserved across different vertebrate species. These results identify major signaling pathways involved in refractive eye development and provide attractive targets for the development of anti-myopia drugs

Genome-wide association studies for corneal and refractive astigmatism in UK Biobank demonstrate a shared role for myopia susceptibility loci.

Previous studies have suggested that naturally occurring genetic variation contributes to the risk of astigmatism. The purpose of this investigation was to identify genetic markers associated with corneal and refractive astigmatism in a large-scale European ancestry cohort (UK Biobank) who underwent keratometry and autorefraction at an assessment centre. Genome-wide association studies for corneal and refractive astigmatism were performed in individuals of European ancestry (N = 86,335 and 88,005 respectively), with the mean corneal astigmatism or refractive astigmatism in fellow eyes analysed as a quantitative trait (dependent variable). Genetic correlation between the two traits was calculated using LD Score regression. Gene-based and gene-set tests were carried out using MAGMA. Single marker-based association tests for corneal astigmatism identified four genome-wide significant loci (P < 5 × 10-8) near the genes ZC3H11B (1q41), LINC00340 (6p22.3), HERC2/OCA2 (15q13.1) and NPLOC4/TSPAN10 (17q25.3). Three of these loci also demonstrated genome-wide significant association with refractive astigmatism: LINC00340, HERC2/OCA2 and NPLOC4/TSPAN10. The genetic correlation between corneal and refractive astigmatism was 0.85 (standard error = 0.068, P = 1.37 × 10-35). Here, we have undertaken the largest genome-wide association studies for corneal and refractive astigmatism to date and identified four novel loci for corneal astigmatism, two of which were also novel loci for refractive astigmatism. These loci have previously demonstrated association with axial length (ZC3H11B), myopia (NPLOC4), spherical equivalent refractive error (LINC00340) and eye colour (HERC2). The shared role of these novel candidate genes for astigmatism lends further support to the shared genetic susceptibility of myopia and astigmatism

Time Outdoors at Specific Ages During Early Childhood and the Risk of Incident Myopia.

PURPOSE: Time outdoors during childhood is negatively associated with incident myopia. Consequently, additional time outdoors has been suggested as a public health intervention to reduce the prevalence of myopia. We investigated whether there were specific ages during early childhood when the time outdoors versus incident myopia association was strongest. METHODS: Children participating in the Avon Longitudinal Study of Parents and Children (ALSPAC) were studied from age 2 to 15 years. Parentally reported time outdoors and time spent reading were assessed longitudinally in early childhood (ages 2, 3, 4, 5, 7, and 9 years). Noncycloplegic autorefraction was carried out longitudinally in later childhood (ages 10, 11, 12, and 15 years). Information was available for 2833 participants. Cox proportional hazards regression was used to test for association between time outdoors and incident myopia. RESULTS: From 3 years of age onward, greater time outdoors was associated with a reduced risk of incident myopia. The hazard ratio for myopia changed progressively from 0.90 (95% CI 0.83-0.98, P = 0.012) at age 3 years, to 0.86 (95% CI 0.78-0.93, P = 0.001) at age 9 years, for each additional SD of time spent outdoors per day. These associations were independent of two major risk factors for myopia: time reading and number of myopic parents. CONCLUSIONS: Additional time spent outdoors across the 3 to 9 years age range was associated with a reduced incidence of myopia between ages 10 and 15 years. There was a trend for the association to increase toward the older end of the 3 to 9 years range.Supported by the UK Medical Research Council. The Wellcome Trust (Grant 102215/2/13/2) and the University of Bristol provide core support for ALSPAC. This research was specifically funded by National Institute for Health Research Career Development Fellowship CDF-2009-02-35 (CW) and a Wellcome Trust Institutional Strategic Support Fund Populations Pilot Award (Grant 508353/509506)

Analysis of genetic networks regulating refractive eye development in collaborative cross progenitor strain mice reveals new genes and pathways underlying human myopia

Background Population studies suggest that genetic factors play an important role in refractive error development; however, the precise role of genetic background and the composition of the signaling pathways underlying refractive eye development remain poorly understood. Methods Here, we analyzed normal refractive development and susceptibility to form-deprivation myopia in the eight progenitor mouse strains of the Collaborative Cross (CC). We used RNA-seq to analyze gene expression in the retinae of these mice and reconstruct genetic networks and signaling pathways underlying refractive eye development. We also utilized genome-wide gene-based association analysis to identify mouse genes and pathways associated with myopia in humans. Results Genetic background strongly influenced both baseline refractive development and susceptibility to environmentally-induced myopia. Baseline refractive errors ranged from − 21.2 diopters (D) in 129S1/svlmj mice to + 22.0 D in CAST/EiJ mice and represented a continuous distribution typical of a quantitative genetic trait. The extent of induced form-deprivation myopia ranged from − 5.6 D in NZO/HILtJ mice to − 20.0 D in CAST/EiJ mice and also followed a continuous distribution. Whole-genome (RNA-seq) gene expression profiling in retinae from CC progenitor strains identified genes whose expression level correlated with either baseline refractive error or susceptibility to myopia. Expression levels of 2,302 genes correlated with the baseline refractive state of the eye, whereas 1,917 genes correlated with susceptibility to induced myopia. Genome-wide gene-based association analysis in the CREAM and UK Biobank human cohorts revealed that 985 of the above genes were associated with myopia in humans, including 847 genes which were implicated in the development of human myopia for the first time. Although the gene sets controlling baseline refractive development and those regulating susceptibility to myopia overlapped, these two processes appeared to be controlled by largely distinct sets of genes. Conclusions Comparison with data for other animal models of myopia revealed that the genes identified in this study comprise a well-defined set of retinal signaling pathways, which are highly conserved across different vertebrate species. These results identify major signaling pathways involved in refractive eye development and provide attractive targets for the development of anti-myopia drugs

Genome-wide Association Study for Vitamin D Levels Reveals 69 Independent Loci.

We aimed to increase our understanding of the genetic determinants of vitamin D levels by undertaking a large-scale genome-wide association study (GWAS) of serum 25 hydroxyvitamin D (25OHD). To do so, we used imputed genotypes from 401,460 white British UK Biobank participants with available 25OHD levels, retaining single-nucleotide polymorphisms (SNPs) with minor allele frequency (MAF) > 0.1% and imputation quality score > 0.3. We performed a linear mixed model GWAS on standardized log-transformed 25OHD, adjusting for age, sex, season of measurement, and vitamin D supplementation. These results were combined with those from a previous GWAS including 42,274 Europeans. In silico functional follow-up of the GWAS results was undertaken to identify enrichment in gene sets, pathways, and expression in tissues, and to investigate the partitioned heritability of 25OHD and its shared heritability with other traits. Using this approach, the SNP heritability of 25OHD was estimated to 16.1%. 138 conditionally independent SNPs were detected (p value < 6.6 × 10-9) among which 53 had MAF < 5%. Single variant association signals mapped to 69 distinct loci, among which 63 were previously unreported. We identified enrichment in hepatic and lipid metabolism gene pathways and enriched expression of the 25OHD genes in liver, skin, and gastrointestinal tissues. We observed partially shared heritability between 25OHD and socio-economic traits, a feature which may be mediated through time spent outdoors. Therefore, through a large 25OHD GWAS, we identified 63 loci that underline the contribution of genes outside the vitamin D canonical metabolic pathway to the genetic architecture of 25OHD

A commonly occurring genetic variant within the NPLOC4-TSPAN10-PDE6G gene cluster is associated with the risk of strabismus.

Strabismus refers to an abnormal alignment of the eyes leading to the loss of central binocular vision. Concomitant strabismus occurs when the angle of deviation is constant in all positions of gaze and often manifests in early childhood when it is considered to be a neurodevelopmental disorder of the visual system. As such, it is inherited as a complex genetic trait, affecting 2-4% of the population. A genome-wide association study (GWAS) for self-reported strabismus (1345 cases and 65,349 controls from UK Biobank) revealed a single genome-wide significant locus on chromosome 17q25. Approximately 20 variants across the NPLOC4-TSPAN10-PDE6G gene cluster and in almost perfect linkage disequilibrium (LD) were most strongly associated (lead variant: rs75078292, OR = 1.26, p = 2.24E-08). A recessive model provided a better fit to the data than an additive model. Association with strabismus was independent of refractive error, and the degree of association with strabismus was minimally attenuated after adjustment for amblyopia. The association with strabismus was replicated in an independent cohort of clinician-diagnosed children aged 7 years old (116 cases and 5084 controls; OR = 1.85, p = 0.009). The associated variants included 2 strong candidate causal variants predicted to have functional effects: rs6420484, which substitutes tyrosine for a conserved cysteine (C177Y) in the TSPAN10 gene, and a 4-bp deletion variant, rs397693108, predicted to cause a frameshift in TSPAN10. The population-attributable risk for the locus was approximately 8.4%, indicating an important role in conferring susceptibility to strabismus

Time spent outdoors partly accounts for the effect of education on myopia

Purpose: The purpose of this study was to investigate if education contributes to the risk of myopia because educational activities typically occur indoors or because of other factors, such as prolonged near viewing. Methods: This was a two-sample Mendelian randomization study. Participants were from the UK Biobank, Avon Longitudinal Study of Parents and Children, and Generation R. Genetic variants associated with years spent in education or time spent outdoors were used as instrumental variables. The main outcome measures were: (1) spherical equivalent refractive error attained by adulthood, and (2) risk of an early age-of-onset of spectacle wear (EAOSW), defined as an age-of-onset of 15 years or below. Results: Time spent outdoors was found to have a small genetic component (heritability 9.8%) that tracked from childhood to adulthood. A polygenic score for time outdoors was associated with children's time outdoors; a polygenic score for years spent in education was inversely associated with children's time outdoors. Accounting for the relationship between time spent outdoors and myopia in a multivariable Mendelian randomization analysis reduced the size of the causal effect of more years in education on myopia to −0.17 diopters (D) per additional year of formal education (95% confidence interval [CI] = −0.32 to −0.01) compared with the estimate from a univariable Mendelian randomization analysis of −0.27 D per year (95% CI = −0.41 to −0.13). Comparable results were obtained for the outcome EAOSW. Conclusions: Accounting for the effects of time outdoors reduced the estimated causal effect of education on myopia by 40%. These results suggest about half of the relationship between education and myopia may be mediated by children not being outdoors during schooling

Quantitative Evidence for Revising the Definition of Primary Graft Dysfunction after Lung Transplant

RATIONALE: Primary graft dysfunction (PGD) is a form of acute lung injury that occurs after lung transplantation. The definition of PGD was standardized in 2005. Since that time, clinical practice has evolved, and this definition is increasingly used as a primary endpoint for clinical trials; therefore, validation is warranted. OBJECTIVES: We sought to determine whether refinements to the 2005 consensus definition could further improve construct validity. METHODS: Data from the Lung Transplant Outcomes Group multicenter cohort were used to compare variations on the PGD definition, including alternate oxygenation thresholds, inclusion of additional severity groups, and effects of procedure type and mechanical ventilation. Convergent and divergent validity were compared for mortality prediction and concurrent lung injury biomarker discrimination. MEASUREMENTS AND MAIN RESULTS: A total of 1,179 subjects from 10 centers were enrolled from 2007 to 2012. Median length of follow-up was 4 years (interquartile range = 2.4-5.9). No mortality differences were noted between no PGD (grade 0) and mild PGD (grade 1). Significantly better mortality discrimination was evident for all definitions using later time points (48, 72, or 48-72 hours; P < 0.001). Biomarker divergent discrimination was superior when collapsing grades 0 and 1. Additional severity grades, use of mechanical ventilation, and transplant procedure type had minimal or no effect on mortality or biomarker discrimination. CONCLUSIONS: The PGD consensus definition can be simplified by combining lower PGD grades. Construct validity of grading was present regardless of transplant procedure type or use of mechanical ventilation. Additional severity categories had minimal impact on mortality or biomarker discrimination

The relationship between plasma lipid peroxidation products and primary graft dysfunction after lung transplantation is modified by donor smoking and reperfusion hyperoxia

BACKGROUND: Donor smoking history and higher fraction of inspired oxygen (FIO2) at reperfusion are associated with primary graft dysfunction (PGD) after lung transplantation. We hypothesized that oxidative injury biomarkers would be elevated in PGD, with higher levels associated with donor exposure to cigarette smoke and recipient hyperoxia at reperfusion. METHODS: We performed a nested case-control study of 72 lung transplant recipients from the Lung Transplant Outcomes Group cohort. Using mass spectroscopy, F2-isoprostanes and isofurans were measured in plasma collected after transplantation. Cases were defined in 2 ways: grade 3 PGD present at day 2 or day 3 after reperfusion (severe PGD) or any grade 3 PGD (any PGD). RESULTS: There were 31 severe PGD cases with 41 controls and 35 any PGD cases with 37 controls. Plasma F2-isoprostane levels were higher in severe PGD cases compared with controls (28.6 pg/ml vs 19.8 pg/ml, p = 0.03). Plasma F2-isoprostane levels were higher in severe PGD cases compared with controls (29.6 pg/ml vs 19.0 pg/ml, p = 0.03) among patients reperfused with FIO2 >40%. Among recipients of lungs from donors with smoke exposure, plasma F2-isoprostane (38.2 pg/ml vs 22.5 pg/ml, p = 0.046) and isofuran (66.9 pg/ml vs 34.6 pg/ml, p = 0.046) levels were higher in severe PGD compared with control subjects. CONCLUSIONS: Plasma levels of lipid peroxidation products are higher in patients with severe PGD, in recipients of lungs from donors with smoke exposure, and in recipients exposed to higher Fio2 at reperfusion. Oxidative injury is an important mechanism of PGD and may be magnified by donor exposure to cigarette smoke and hyperoxia at reperfusion

The association between circulating 25-hydroxyvitamin D metabolites and type 2 diabetes in European populations: A meta-analysis and Mendelian randomisation analysis

Funder: NIHR Biomedical Research Centre Cambridge: Nutrition, Diet, and Lifestyle Research Theme; Grant(s): IS-BRC-1215-20014Background: Prior research suggested a differential association of 25-hydroxyvitamin D (25(OH)D) metabolites with type 2 diabetes (T2D), with total 25(OH)D and 25(OH)D3 inversely associated with T2D, but the epimeric form (C3-epi-25(OH)D3) positively associated with T2D. Whether or not these observational associations are causal remains uncertain. We aimed to examine the potential causality of these associations using Mendelian randomisation (MR) analysis. Methods and findings: We performed a meta-analysis of genome-wide association studies for total 25(OH)D (N = 120,618), 25(OH)D3 (N = 40,562), and C3-epi-25(OH)D3 (N = 40,562) in participants of European descent (European Prospective Investigation into Cancer and Nutrition [EPIC]–InterAct study, EPIC-Norfolk study, EPIC-CVD study, Ely study, and the SUNLIGHT consortium). We identified genetic variants for MR analysis to investigate the causal association of the 25(OH)D metabolites with T2D (including 80,983 T2D cases and 842,909 non-cases). We also estimated the observational association of 25(OH)D metabolites with T2D by performing random effects meta-analysis of results from previous studies and results from the EPIC-InterAct study. We identified 10 genetic loci associated with total 25(OH)D, 7 loci associated with 25(OH)D3 and 3 loci associated with C3-epi-25(OH)D3. Based on the meta-analysis of observational studies, each 1–standard deviation (SD) higher level of 25(OH)D was associated with a 20% lower risk of T2D (relative risk [RR]: 0.80; 95% CI 0.77, 0.84; p < 0.001), but a genetically predicted 1-SD increase in 25(OH)D was not significantly associated with T2D (odds ratio [OR]: 0.96; 95% CI 0.89, 1.03; p = 0.23); this result was consistent across sensitivity analyses. In EPIC-InterAct, 25(OH)D3 (per 1-SD) was associated with a lower risk of T2D (RR: 0.81; 95% CI 0.77, 0.86; p < 0.001), while C3-epi-25(OH)D3 (above versus below lower limit of quantification) was positively associated with T2D (RR: 1.12; 95% CI 1.03, 1.22; p = 0.006), but neither 25(OH)D3 (OR: 0.97; 95% CI 0.93, 1.01; p = 0.14) nor C3-epi-25(OH)D3 (OR: 0.98; 95% CI 0.93, 1.04; p = 0.53) was causally associated with T2D risk in the MR analysis. Main limitations include the lack of a non-linear MR analysis and of the generalisability of the current findings from European populations to other populations of different ethnicities. Conclusions: Our study found discordant associations of biochemically measured and genetically predicted differences in blood 25(OH)D with T2D risk. The findings based on MR analysis in a large sample of European ancestry do not support a causal association of total 25(OH)D or 25(OH)D metabolites with T2D and argue against the use of vitamin D supplementation for the prevention of T2D