Genome-wide Association Study of Axial Length in Population-based Cohorts in Japan

Population-specific putative causal variants shape quantitative traits

Human genetic variants are associated with many traits through largely unknown mechanisms. Here, combining approximately 260,000 Japanese study participants, a Japanese-specific genotype reference panel and statistical fine-mapping, we identified 4,423 significant loci across 63 quantitative traits, among which 601 were new, and 9,406 putatively causal variants. New associations included Japanese-specific coding, splicing and noncoding variants, exemplified by a damaging missense variant rs730881101 in TNNT2 associated with lower heart function and increased risk for heart failure (P = 1.4 × 10-15 and odds ratio = 4.5, 95% confidence interval = 3.1-6.5). Putative causal noncoding variants were supported by state-of-art in silico functional assays and had comparable effect sizes to coding variants. A plausible example of new mechanisms of causal variants is an enrichment of causal variants in 3' untranslated regions (UTRs), including the Japanese-specific rs13306436 in IL6 associated with pro-inflammatory traits and protection against tuberculosis. We experimentally showed that transcripts with rs13306436 are resistant to mRNA degradation by regnase-1, an RNA-binding protein. Our study provides a list of fine-mapped causal variants to be tested for functionality and underscores the importance of sequencing, genotyping and association efforts in diverse populations.

Associations of genetic variants for refractive error and axial length in adults with ocular endophenotypes in children: a cross-sectional and longitudinal study

AIMS

To investigate the associations of genetic variants previously linked to axial length (AL) and spherical equivalent refraction (SE) in adults with refractive error and related endophenotypes in children, at baseline and 3-year follow-up.

METHODS

15 candidate single-nucleotide polymorphisms (SNPs), selected from previous Genome-Wide Association Studies and meta-analyses, were genotyped in 2819 Chinese children, who had undergone baseline and 3-year follow-up cycloplegic refraction, ocular biometry and ocular health examinations. Linear regression analyses were conducted to assess the associations of the SNPs with baseline measurements and longitudinal changes in SE, spherical power (SPH), AL, corneal radius of curvature (CR) and AL/CR ratio.

RESULTS

SNPs ZMAT4 rs7829127, ZMAT4 rs16890057, TOX rs7837791, GRIA4 rs11601239 and RDH5 rs3138142 were associated with SE (β=0.233, p=4.21×10-4; β=0.221, p=7.87×10-4; β=0.106, p=0.0076; β=0.084, p=0.041; β=0.14, p=0.013, respectively) and SPH (β=0.24, p=2.3×10-4; β=0.232, p=3.8×10-4; β=0.088, p=0.025; β=0.086, p=0.034; β=0.14, p=0.012, respectively). Among them, ZMAT4 rs7829127 and rs16890057, were also associated with AL (β=-0.128, p=5.6×10-4; β=-0.128, p=5.21×10-4) and AL/CR ratio (β=-0.014, p=0.0028; β=-0.014, p=0.0034), whereas TOX rs7837791 was associated with AL (β=-0.062, p=0.0058) and GRIA4 11 601 239 with AL/CR ratio (β=-0.0058, p=0.049). Additionally, CD55 rs1652333 and RDH5 rs3138142 were associated with 3-year longitudinal changes in AL (β=0.062, p=0.018; β=-0.079, p=0.029) and CR (β=0.014, p=0.027; β=-0.018, p=0.035).

CONCLUSION

Among SNPs previously associated with AL and SE in adults, variants in ZMAT4, TOX and GRIA4 were associated with AL, SE, SPH, and/or AL/CR ratio, while variants in RDH5 and CD55 showed associations with AL and CR changes in children.

Progress Report of the Tohoku Medical Megabank Community-based Cohort Study: Study Profile of the Repeated Center-based Survey During Second Period in Miyagi Prefecture

BACKGROUND

The purpose of this study was to report the basic profile of the Miyagi Prefecture part of a repeated center-based survey during the second period of the Tohoku Medical Megabank Community-Based Cohort Study (TMM CommCohort Study), as well as the participants' characteristics based on their participation type in the baseline survey.

METHODS

The second period survey, conducted from June 2017 to March 2021, included participants of the TMM CommCohort Study (May 2013 to March 2016). In addition to the questionnaire, blood, urine, and physiological function tests were performed during the second period survey. There were three main ways of participation in the baseline survey: Type 1, Type 1 additional, or Type 2 survey. The second period survey was conducted in the same manner as the Type 2 survey, which was based on the community support center (CSC).

RESULTS

In Miyagi Prefecture, 29,383 (57.7%) of 50,967 participants participated in the second period survey. The participation rate among individuals who had visited the CSC was approximately 80%. Although some factors differed depending on the participation type in the baseline survey, the second period survey respondents in the Type 1 and Type 2 survey groups at baseline had similar traits.

CONCLUSION

The second period survey of the TMM CommCohort Study provided detailed follow-up information. Following up on the health conditions of the participants will clarify the long-term effects of disasters and contribute to personalized prevention.

Genetic Risk, Healthy Lifestyle Adherence, and Risk of Developing Diabetes in the Japanese Population

AIM

This study examined the relationship between genetic risk, healthy lifestyle, and risk of developing diabetes.

METHODS

This prospective cohort study included 11,014 diabetes-free individuals ≥ 20 years old from the Tohoku Medical Megabank Community-based cohort study. Lifestyle scores, including the body mass index, smoking, physical activity, and gamma-glutamyl transferase (marker of alcohol consumption), were assigned, and participants were categorized into ideal, intermediate, and poor lifestyles. A polygenic risk score (PRS) was constructed based on the type 2 diabetes loci from the BioBank Japan study. A multiple logistic regression model was used to estimate the association between genetic risk, healthy lifestyle, and diabetes incidence and to calculate the area under the receiver operating characteristic curve (AUROC).

RESULT

Of the 11,014 adults included (67.8% women; mean age [standard deviation], 59.1 [11.3] years old), 297 (2.7%) developed diabetes during a mean 4.3 (0.8) years of follow-up. Genetic and lifestyle score is independently associated with the development of diabetes. Compared with the low genetic risk and ideal lifestyle groups, the odds ratio was 3.31 for the low genetic risk and poor lifestyle group. When the PRS was integrated into a model including the lifestyle and family history, the AUROC significantly improved to 0.719 (95% confidence interval [95% CI]: 0.692-0.747) compared to a model including only the lifestyle and family history (0.703 [95% CI, 0.674-0.732]).

CONCLUSION

Our findings indicate that adherence to a healthy lifestyle is important for preventing diabetes, regardless of genetic risk. In addition, genetic risk might provide information beyond lifestyle and family history to stratify individuals at high risk of developing diabetes.

Bi-Gaussian analysis reveals distinct education-related alterations in spherical equivalent and axial length-results from the Gutenberg Health Study

PURPOSE

The aim of this study is to investigate the distribution of spherical equivalent and axial length in the general population and to analyze the influence of education on spherical equivalent with a focus on ocular biometric parameters.

METHODS

The Gutenberg Health Study is a population-based cohort study in Mainz, Germany. Participants underwent comprehensive ophthalmologic examinations as part of the 5-year follow-up examination in 2012-2017 including genotyping. The spherical equivalent and axial length distributions were modeled with gaussian mixture models. Regression analysis (on person-individual level) was performed to analyze associations between biometric parameters and educational factors. Mendelian randomization analysis explored the causal effect between spherical equivalent, axial length, and education. Additionally, effect mediation analysis examined the link between spherical equivalent and education.

RESULTS

A total of 8532 study participants were included (median age: 57 years, 49% female). The distribution of spherical equivalent and axial length follows a bi-Gaussian function, partially explained by the length of education (i.e., < 11 years education vs. 11-20 years). Mendelian randomization indicated an effect of education on refractive error using a genetic risk score of education as an instrument variable (- 0.35 diopters per SD increase in the instrument, 95% CI, - 0.64-0.05, p = 0.02) and an effect of education on axial length (0.63 mm per SD increase in the instrument, 95% CI, 0.22-1.04, p = 0.003). Spherical equivalent, axial length and anterior chamber depth were associated with length of education in regression analyses. Mediation analysis revealed that the association between spherical equivalent and education is mainly driven (70%) by alteration in axial length.

CONCLUSIONS

The distribution of axial length and spherical equivalent is represented by subgroups of the population (bi-Gaussian). This distribution can be partially explained by length of education. The impact of education on spherical equivalent is mainly driven by alteration in axial length.

Novel loci for ocular axial length identified through extreme-phenotype genome-wide association study in Chinese populations

PURPOSE

To investigate genetic loci associated with ocular axial length (AL) in the Chinese population.

METHODS

A genome-wide association study meta-analysis was conducted in totalling 2644 Chinese individuals from 3 cohorts: the Guangzhou cohort (GZ, 537 high myopes and 151 hyperopes), Wenzhou cohort (334 high myopes and 6 hyperopes) and Guangzhou Twin Eye Study (1051 participants with normally distributed AL). Functional mapping was performed to annotate the significant signals, possible tissues and cell types by integrating available multiomics data. Logistic regression models using AL-associated SNPs were constructed to predict three AL status in GZ.

RESULTS

Two novel loci (1q25.2 FAM163A and 7p22.2 SDK1) showed genome-wide significant associations with AL, together explaining 29.63% of AL variance in GZ. The two lead SNPs improved the prediction accuracy for AL status, especially for hyperopes. The frequencies of AL decreasing (less myopic) alleles of the two SNPs were lowest in East Asians as compared with other populations (rs17370084: f EAS=0.03, f EUR=0.24, f AFR=0.05; rs73046501: f EAS=0.06, f EUR=0.07, f AFR=0.20), which was in line with the global distribution of myopia. The cerebral cortex and gamma-aminobutyric acidergic interneurons showed possible functional involvement in myopia development, and the galactose metabolic pathways were significantly enriched.

CONCLUSION

Our study identified two population-specific novel loci for AL, expanding our understanding of the genetic basis of AL and providing evidence for a role of the nervous system and glucose metabolism in myopia pathogenesis.

A Case of Strabismus Fixus With Extremely Long Axial Length Results in Improvement of Limitation of Abduction by the Yokoyama Procedure

Acquired strabismus in high myopia is typically fixed in the positions of adduction and depression, with restrictions in both abduction and elevation. As a treatment for myopic strabismus fixus, the Yokoyama procedure is effective. We report a case of strabismus fixus with a long axial length (34 mm), in which abduction limitation was improved by the Yokoyama procedure with medial rectus recession. A 68-year-old woman was referred for strabismus fixus in her right eye. Her right eye was fixed in the positions of adduction and depression, with restrictions in both abduction and elevation. The axial length of her right eye was extremely long 33.97mm. Magnetic resonance imaging (MRI) showed that the posterior eyeball of her right eye had dislocated out of the superotemporal muscle cone, and she was diagnosed with strabismus fixus with high myopia. She underwent the Yokoyama procedure in her right eye, and medial rectus recession was performed at the same time because abduction limitations remained at the end of the Yokoyama procedure. After surgery, there was a small residual esotropia, but abduction beyond the midline was possible, and the patient's satisfaction was high. A combination of the Yokoyama procedure and medial rectus recession for a patient with myopic strabismus fixus with long axial length resulted in good improvement of ocular misalignment and limitation of abduction.

Association of physiological factors with grip and leg extension strength: tohoku medical megabank community-based cohort study

BACKGROUND

Upper and lower extremity muscle strength can be used to predict health outcomes. However, the difference between the relation of upper extremity muscle and of lower extremity muscle with physiological factors is unclear. This study aimed to evaluate the association between physiological data and muscle strength, measured using grip and leg extension strength, among Japanese adults.

METHODS

We conducted a cross-sectional study of 2,861 men and 6,717 women aged ≥ 20 years living in Miyagi Prefecture, Japan. Grip strength was measured using a dynamometer. Leg extension strength was measured using a hydraulic isokinetic leg press machine. Anthropometry and physiological data, including blood pressure, calcaneal ultrasound bone status, pulmonary function, carotid echography, and blood information, were assessed. We used a general linear model adjusted for age, body composition, and smoking status to evaluate the association between muscle strength and physiological factors.

RESULTS

Grip and leg extension strength were positively associated with bone area ratio, vital capacity, forced vital capacity, forced expiratory volume in one second, and estimated glomerular filtration rate, and negatively associated with waist circumference and percentage body fat mass in both the sexes. Diastolic blood pressure was positively associated with grip strength in both the sexes and leg extension strength in men, but not women. High-density lipoprotein cholesterol and red blood cell counts were positively associated with grip and leg extension strength in women, but not men. In both the sexes, pulse rate, total cholesterol, and uric acid were consistently associated with only leg extension strength, but not grip strength. In women, glycated hemoglobin demonstrated negative and positive associations with grip and leg extension strength, respectively.

CONCLUSIONS

Grip and leg extension strength demonstrated similar associations with anthropometry, pulmonary function, and estimated glomerular filtration rate, but the associations with the other factors were not always consistent.

Patient Factors Influencing Intraocular Penetration of Brimonidine-Related Eye Drops in Adults: A Post Hoc Pooled Analysis

INTRODUCTION

The factors related to the ocular penetration of drugs after the administration of eye drops in humans have not been examined in detail. Therefore, this study assessed the influence of patient factors on the intraocular penetration of eye drops.

METHODS

A pooled analysis was performed on the data of 42 participants from three studies to evaluate the ocular pharmacokinetics in humans after the topical application of brimonidine-related eye drops. The patients were scheduled for vitrectomy and received brimonidine-related eye drops (0.1% brimonidine tartrate ophthalmic solution, 0.1% brimonidine tartrate and 0.5% timolol fixed-combination ophthalmic solution, or 0.1% brimonidine tartrate and 1% brinzolamide fixed-combination suspension) twice daily for 1 week. We analyzed the effects of patient factors (sex, the presence or absence of lens, age, corneal thickness, corneal endothelial cell density, tear secretion, eye axial length, height, weight and body mass index [BMI]) on brimonidine, timolol and brinzolamide concentrations in the aqueous and vitreous humor after topical application.

RESULTS

The drug concentrations in the aqueous and vitreous humor were not significantly different, regardless of sex or the presence or absence of lens. Age correlated positively with brimonidine (r = 0.3948, p = 0.012) and brinzolamide (r = 0.6809, p = 0.030) concentrations in the aqueous humor; the correlation with timolol showed a trend towards significance (r = 0.6425, p = 0.086). Corneal thickness, corneal endothelial cell density, tear secretion, eye axial length, height and BMI did not correlate with the drug concentrations in the aqueous or vitreous humor. Timolol concentration in the vitreous humor was negatively correlated with weight (r =  - 0.8333, p = 0.010).

CONCLUSION

The findings of this study emphasize the necessity of considering individual differences in ocular pharmacokinetics during drug therapy (formulation design of the eye drops and dose regimen).

Epitranscriptomic investigation of myopia-associated RNA editing in the retina

Myopia is one of the most common causes of vision loss globally and is significantly affected by epigenetics. Adenosine-to-inosine (A-to-I RNA) editing is an epigenetic process involved in neurological disorders, yet its role in myopia remains undetermined. We performed a transcriptome-wide analysis of A-to-I RNA editing in the retina of form-deprivation myopia mice. Our study identified 91 A-to-I RNA editing sites in 84 genes associated with myopia. Notably, at least 27 (32.1%) of these genes with myopia-associated RNA editing showed existing evidence to be associated with myopia or related ocular phenotypes in humans or animal models, such as very low-density lipoprotein receptor (Vldlr) in retinal neovascularization and hypoxia-induced factor 1 alpha (Hif1a). Moreover, functional enrichment showed that RNA editing enriched in FDM was primarily involved in response to fungicides, a potentially druggable process for myopia prevention, and epigenetic regulation. In contrast, RNA editing enriched in controls was mostly involved in post-embryonic eye morphogenesis. Our results demonstrate altered A-to-I RNA editing associated with myopia in an experimental mouse model and warrant further study on its role in myopia development.

A multiethnic genome-wide analysis of 19,420 individuals identifies novel loci associated with axial length and shared genetic influences with refractive error and myopia

Introduction: Long axial length (AL) is a risk factor for myopia. Although family studies indicate that AL has an important genetic component with heritability estimates up to 0.94, there have been few reports of AL-associated loci. Methods: Here, we conducted a multiethnic genome-wide association study (GWAS) of AL in 19,420 adults of European, Latino, Asian, and African ancestry from the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort, with replication in a subset of the Consortium for Refractive Error and Myopia (CREAM) cohorts of European or Asian ancestry. We further examined the effect of the identified loci on the mean spherical equivalent (MSE) within the GERA cohort. We also performed genome-wide genetic correlation analyses to quantify the genetic overlap between AL and MSE or myopia risk in the GERA European ancestry sample. Results: Our multiethnic GWA analysis of AL identified a total of 16 genomic loci, of which 5 are novel. We found that all AL-associated loci were significantly associated with MSE after Bonferroni correction. We also found that AL was genetically correlated with MSE (rg = -0.83; SE, 0.04; p = 1.95 × 10-89) and myopia (rg = 0.80; SE, 0.05; p = 2.84 × 10-55). Finally, we estimated the array heritability for AL in the GERA European ancestry sample using LD score regression, and found an overall heritability estimate of 0.37 (s.e. = 0.04). Discussion: In this large and multiethnic study, we identified novel loci, associated with AL at a genome-wide significance level, increasing substantially our understanding of the etiology of AL variation. Our results also demonstrate an association between AL-associated loci and MSE and a shared genetic basis between AL and myopia risk.

Sex differences in the relationship between axial length and dry eye in elderly patients

Purpose

The aim of this study was to explore the association between myopia and dry eye (DE)-related ocular parameters.

Methods

We recruited a total of 460 patients (mean age, 73.6 years; 40.2% men) and performed DE-related, axial length (AL) and retinal examinations. Statistical analysis revealed a significant sex difference in AL, strip meniscometry value, corneal staining score, corneal endothelial cell density, ganglion cell complex (GCC) thickness, and full macular thickness. AL was strongly age- and sex-dependent, so subsequent analyses were stratified by sex.

Results

Among DE-related parameters, strip meniscometry value (ß = -0.167, p = 0.033) and corneal endothelial cell density (ß = -0.139, p = 0.023) were correlated with AL in women but not in men. Regarding retinal parameters, GCC thickness and full macular thickness were correlated with AL in women but not in men.

Conclusion

The current results suggest a relationship between tear production and AL in elderly women and support the hypothesis that there may be a common upstream factor including the parasympathetic nervous system in the association between tear production and AL or DE and myopia.

A new polygenic score for refractive error improves detection of children at risk of high myopia but not the prediction of those at risk of myopic macular degeneration

BACKGROUND

High myopia (HM), defined as a spherical equivalent refractive error (SER) ≤ -6.00 diopters (D), is a leading cause of sight impairment, through myopic macular degeneration (MMD). We aimed to derive an improved polygenic score (PGS) for predicting children at risk of HM and to test if a PGS is predictive of MMD after accounting for SER.

METHODS

The PGS was derived from genome-wide association studies in participants of UK Biobank, CREAM Consortium, and Genetic Epidemiology Research on Adult Health and Aging. MMD severity was quantified by a deep learning algorithm. Prediction of HM was quantified as the area under the receiver operating curve (AUROC). Prediction of severe MMD was assessed by logistic regression.

FINDINGS

In independent samples of European, African, South Asian and East Asian ancestry, the PGS explained 19% (95% confidence interval 17-21%), 2% (1-3%), 8% (7-10%) and 6% (3-9%) of the variation in SER, respectively. The AUROC for HM in these samples was 0.78 (0.75-0.81), 0.58 (0.53-0.64), 0.71 (0.69-0.74) and 0.67 (0.62-0.72), respectively. The PGS was not associated with the risk of MMD after accounting for SER: OR = 1.07 (0.92-1.24).

INTERPRETATION

Performance of the PGS approached the level required for clinical utility in Europeans but not in other ancestries. A PGS for refractive error was not predictive of MMD risk once SER was accounted for.

FUNDING

Supported by the Welsh Government and Fight for Sight (24WG201).

Design and Progress of Child Health Assessments at Community Support Centers in the Birth and Three-Generation Cohort Study of the Tohoku Medical Megabank Project

The Tohoku Medical Megabank Project (TMM) has been conducting a birth and three-generation cohort study (the BirThree Cohort Study). We recruited 73,529 pregnant women and their family members for this cohort study, which included 23,143 newborns and 9,459 of their siblings. We designed and are in the process of conducting three-step health assessments for each newborn at approximately ages of 5, 10 and 16. These health assessments are administered at seven community support centers. Trained genome medical research coordinators conduct physical examinations of and collect biological specimens from each participant. The Sendai Children's Health Square has been established as the headquarters for these child health assessments and is utilized to accumulate knowledge that can facilitate the proper practice of child health assessments. We designed all the relevant health assessments facilities to allow parents and their children to participate in the health assessments concomitantly. Our centers serve as places where child participants and their parents can feel at ease as a result of the implementation of safety measures and child hospitality measures. The TMM BirThree Cohort Study is in the process of conducting strategically detailed health assessments and genome analysis, which can facilitate studies concerning the gene-environment interactions relevant to noncommunicable diseases. Through these operations, our study allows for a significant depth of data to be collected in terms of the number of biospecimens under study and the comprehensiveness of both basic and clinical data alongside relevant family information.