Whole-exome sequencing study identifies rare variants and genes associated with intraocular pressure and glaucoma

Integration of variant annotations using deep set networks boosts rare variant association testing

Rare genetic variants can have strong effects on phenotypes, yet accounting for rare variants in genetic analyses is statistically challenging due to the limited number of allele carriers and the burden of multiple testing. While rich variant annotations promise to enable well-powered rare variant association tests, methods integrating variant annotations in a data-driven manner are lacking. Here we propose deep rare variant association testing (DeepRVAT), a model based on set neural networks that learns a trait-agnostic gene impairment score from rare variant annotations and phenotypes, enabling both gene discovery and trait prediction. On 34 quantitative and 63 binary traits, using whole-exome-sequencing data from UK Biobank, we find that DeepRVAT yields substantial gains in gene discoveries and improved detection of individuals at high genetic risk. Finally, we demonstrate how DeepRVAT enables calibrated and computationally efficient rare variant tests at biobank scale, aiding the discovery of genetic risk factors for human disease traits.

ACAD10 and ACAD11 allow entry of 4-hydroxy fatty acids into β-oxidation

Hydroxylated fatty acids are important intermediates in lipid metabolism and signaling. Surprisingly, the metabolism of 4-hydroxy fatty acids remains largely unexplored. We found that both ACAD10 and ACAD11 unite two enzymatic activities to introduce these metabolites into mitochondrial and peroxisomal β-oxidation, respectively. First, they phosphorylate 4-hydroxyacyl-CoAs via a kinase domain, followed by an elimination of the phosphate to form enoyl-CoAs catalyzed by an acyl-CoA dehydrogenase (ACAD) domain. Studies in knockout cell lines revealed that ACAD10 preferentially metabolizes shorter chain 4-hydroxy fatty acids than ACAD11 (i.e. 6 carbons versus 10 carbons). Yet, recombinant proteins showed comparable activity on the corresponding 4-hydroxyacyl-CoAs. This suggests that the localization of ACAD10 and ACAD11 to mitochondria and peroxisomes, respectively, might influence their physiological substrate spectrum. Interestingly, we observed that ACAD10 is cleaved internally during its maturation generating a C-terminal part consisting of the ACAD domain, and an N-terminal part comprising the kinase domain and a haloacid dehalogenase (HAD) domain. HAD domains often exhibit phosphatase activity, but negligible activity was observed in the case of ACAD10. Yet, inactivation of a presumptive key residue in this domain significantly increased the kinase activity, suggesting that this domain might have acquired a regulatory function to prevent accumulation of the phospho-hydroxyacyl-CoA intermediate. Taken together, our work reveals that 4-hydroxy fatty acids enter mitochondrial and peroxisomal fatty acid β-oxidation via two enzymes with an overlapping substrate repertoire.

Increasing the diagnostic yield of childhood glaucoma cases recruited into the 100,000 Genomes Project

Childhood glaucoma (CG) encompasses a heterogeneous group of genetic eye disorders that is responsible for approximately 5% of childhood blindness worldwide. Understanding the molecular aetiology is key to improving diagnosis, prognosis and unlocking the potential for optimising clinical management. In this study, we investigated 86 CG cases from 78 unrelated families of diverse ethnic backgrounds, recruited into the Genomics England 100,000 Genomes Project (GE100KGP) rare disease cohort, to improve the genetic diagnostic yield. Using the Genomics England/Genomic Medicine Centres (GE/GMC) diagnostic pipeline, 13 unrelated families were solved (13/78, 17%). Further interrogation using an expanded gene panel yielded a molecular diagnosis in 7 more unrelated families (7/78, 9%). This analysis effectively raises the total number of solved CG families in the GE100KGP to 26% (20/78 families). Twenty-five percent (5/20) of the solved families had primary congenital glaucoma (PCG), while 75% (15/20) had secondary CG; 53% of this group had non-acquired ocular anomalies (including iris hypoplasia, megalocornea, ectopia pupillae, retinal dystrophy, and refractive errors) and 47% had non-acquired systemic diseases such as cardiac abnormalities, hearing impairment, and developmental delay. CYP1B1 was the most frequently implicated gene, accounting for 55% (11/20) of the solved families. We identified two novel likely pathogenic variants in the TEK gene, in addition to one novel pathogenic copy number variant (CNV) in FOXC1. Variants that passed undetected in the GE100KGP diagnostic pipeline were likely due to limitations of the tiering process, the use of smaller gene panels during analysis, and the prioritisation of coding SNVs and indels over larger structural variants, CNVs, and non-coding variants.

Rare genetic variation in VE-PTP is associated with central serous chorioretinopathy, venous dysfunction and glaucoma

Central serous chorioretinopathy (CSC) is a fluid maculopathy whose etiology is not well understood. Abnormal choroidal veins in CSC patients have been shown to have similarities with varicose veins. To identify potential mechanisms, we analyzed genotype data from 1,477 CSC patients and 455,449 controls in FinnGen. We identified an association for a low-frequency (AF=0.5%) missense variant (rs113791087) in the gene encoding vascular endothelial protein tyrosine phosphatase (VE-PTP) (OR=2.85, P=4.5×10-9). This was confirmed in a meta-analysis of 2,452 CSC patients and 865,767 controls from 4 studies (OR=3.06, P=7.4×10-15). Rs113791087 was associated with a 56% higher prevalence of retinal abnormalities (35.3% vs 22.6%, P=8.0×10-4) in 708 UK Biobank participants and, surprisingly, with varicose veins (OR=1.31, P=2.3×10-11) and glaucoma (OR=0.82, P=6.9×10-9). Predicted loss-of-function variants in VEPTP, though rare in number, were associated with CSC in All of Us (OR=17.10, P=0.018). These findings highlight the significance of VE-PTP in diverse ocular and systemic vascular diseases.

40 Hz light flickering promotes sleep through cortical adenosine signaling

Flickering light stimulation has emerged as a promising non-invasive neuromodulation strategy to alleviate neuropsychiatric disorders. However, the lack of a neurochemical underpinning has hampered its therapeutic development. Here, we demonstrate that light flickering triggered an immediate and sustained increase (up to 3 h after flickering) in extracellular adenosine levels in the primary visual cortex (V1) and other brain regions, as a function of light frequency and intensity, with maximal effects observed at 40 Hz frequency and 4000 lux. We uncovered cortical (glutamatergic and GABAergic) neurons, rather than astrocytes, as the cellular source, the intracellular adenosine generation from AMPK-associated energy metabolism pathways (but not SAM-transmethylation or salvage purine pathways), and adenosine efflux mediated by equilibrative nucleoside transporter-2 (ENT2) as the molecular pathway responsible for extracellular adenosine generation. Importantly, 40 Hz (but not 20 and 80 Hz) light flickering for 30 min enhanced non-rapid eye movement (non-REM) and REM sleep for 2-3 h in mice. This somnogenic effect was abolished by ablation of V1 (but not superior colliculus) neurons and by genetic deletion of the gene encoding ENT2 (but not ENT1), but recaptured by chemogenetic inhibition of V1 neurons and by focal infusion of adenosine into V1 in a dose-dependent manner. Lastly, 40 Hz light flickering for 30 min also promoted sleep in children with insomnia by decreasing sleep onset latency, increasing total sleep time, and reducing waking after sleep onset. Collectively, our findings establish the ENT2-mediated adenosine signaling in V1 as the neurochemical basis for 40 Hz flickering-induced sleep and unravel a novel and non-invasive treatment for insomnia, a condition that affects 20% of the world population.

Ubiquitin proteasome system and glaucoma: A survey of genetics and molecular biology studies supporting a link with pathogenic and therapeutic relevance

Glaucoma represents a group of progressive neurodegenerative diseases characterized by the loss of retinal ganglion cells (RGCs) and their axons with subsequent visual field impairment. The disease develops through largely uncharacterized molecular mechanisms, that are likely to occur in different localized cell types, either in the anterior (e.g., trabecular meshwork cells) or posterior (e.g., Muller glia, retinal ganglion cells) segments of the eye. Genomic and preclinical studies suggest that glaucoma pathogenesis may develop through altered ubiquitin (Ub) signaling. Ubiquitin conjugation, referred to as ubiquitylation, is a major post-synthetic modification catalyzed by E1-E2-E3 enzymes, that profoundly regulates the turnover, trafficking and biological activity of the targeted protein. The development of new technologies, including proteomics workflows, allows the biology of ubiquitin signaling to be described in health and disease. This post-translational modification is emerging as a key role player in neurodegeneration, gaining relevance for novel therapeutic options, such as in the case of Proteolysis Targeting Chimeras technology. Although scientific evidence supports a link between Ub and glaucoma, their relationship is still not well-understood. Therefore, this review provides a detailed research-oriented discussion on current evidence of Ub signaling in glaucoma. A review of genomic and genetic data is provided followed by an in-depth discussion of experimental data on ASB10, parkin and optineurin, which are proteins that play a key role in Ub signaling and have been associated with glaucoma.

10 Years of GWAS in intraocular pressure

Intraocular pressure (IOP) is the only modifiable risk factor for glaucoma, the leading cause of irreversible blindness worldwide. In this review, we summarize the findings of genome-wide association studies (GWASs) of IOP published in the past 10 years and prior to December 2022. Over 190 genetic loci and candidate genes associated with IOP have been uncovered through GWASs, although most of these studies were conducted in subjects of European and Asian ancestries. We also discuss how these common variants have been used to derive polygenic risk scores for predicting IOP and glaucoma, and to infer causal relationship with other traits and conditions through Mendelian randomization. Additionally, we summarize the findings from a recent large-scale exome-wide association study (ExWAS) that identified rare variants associated with IOP in 40 novel genes, six of which are drug targets for clinical treatment or are being evaluated in clinical trials. Finally, we discuss the need for future genetic studies of IOP to include individuals from understudied populations, including Latinos and Africans, in order to fully characterize the genetic architecture of IOP.

Causal Association between Iritis or Uveitis and Glaucoma: A Two-Sample Mendelian Randomisation Study

Recent studies have suggested an association between iritis or uveitis and glaucoma. This study investigated the causal relationship between glaucoma and iritis and uveitis as exposures in a multi-ethnic population. Single-nucleotide polymorphisms associated with exposures to iritis and uveitis from the genome-wide association study (GWAS) data of Biobank Japan (BBJ) and the meta-analysis data from BBJ and UK Biobank (UKB) were used as instrumental variables (IVs). The GWAS dataset for glaucoma was extracted from the meta-analysis data (n = 240,302) of Genetic Epidemiology Research in Adult Health and Aging and UKB. The casual estimates were assessed with a two-sample Mendelian randomisation (MR) test using the inverse-variance-weighted (IVW) method, weighted median method, MR–Egger method, and MR-Pleiotropy Residual Sum and Outlier test. The IVW method revealed a significant causal association between iritis and glaucoma using IVs (p < 5.0 × 10−8) from the East Asian population (n = 2) (odds ratio [OR] = 1.01, p = 0.017), a significant association between iritis exposures (p < 5.0 × 10−8) in the multi-ethnic population (n = 11) (OR = 1.04, p = 0.001), and a significant causal association between uveitis exposures (n = 10 with p < 5.0 × 10−8) and glaucoma in the multi-ethnic population (OR = 1.04, p = 0.001). Iritis and uveitis had causal effects on glaucoma risk based on IVs from the multi-ethnic population. These findings imply that the current classifications of uveitic glaucoma and open-angle glaucoma overlap, indicating the need for further investigating these complex relationships.