Rare protein-altering variants in ANGPTL7 lower intraocular pressure and protect against glaucoma

Rare genetic variation in PTPRB is associated with central serous chorioretinopathy, varicose veins and glaucoma

Central serous chorioretinopathy is an eye disease characterized by fluid buildup under the central retina whose etiology is not well understood. Abnormal choroidal veins in central serous chorioretinopathy patients have been shown to have similarities with varicose veins. To identify potential mechanisms, we analyzed genotype data from 1,477 patients and 455,449 controls in FinnGen. We identified an association for a low-frequency (allele frequency = 0.5%) missense variant (rs113791087) in PTPRB, the gene encoding vascular endothelial protein tyrosine phosphatase (odds ratio=2.85, P = 4.5 × 10-9). This was confirmed in a meta-analysis of 2,452 patients and 865,767 controls from 4 studies (odds ratio=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 increased risk of varicose veins (odds ratio=1.31, P = 2.3 × 10-11) and reduced risk of glaucoma (odds ratio=0.82, P = 6.9 × 10-9). Predicted loss-of-function variants in PTPRB, though rare in number, were associated with central serous chorioretinopathy in All of Us (odds ratio=17.09, P = 0.018). These findings highlight the significance of vascular endothelial protein tyrosine phosphatase in diverse ocular and systemic veno-vascular diseases.

Efficient storage and regression computation for population-scale genome sequencing studies

MOTIVATION

The growing availability of large-scale population biobanks has the potential to significantly advance our understanding of human health and disease. However, the massive computational and storage demands of whole genome sequencing (WGS) data pose serious challenges, particularly for underfunded institutions or researchers in developing countries. This disparity in resources can limit equitable access to cutting-edge genetic research.

RESULTS

We present novel algorithms and regression methods that dramatically reduce both computation time and storage requirements for WGS studies, with particular attention to rare variant representation. By integrating these approaches into PLINK 2.0, we demonstrate substantial gains in efficiency without compromising analytical accuracy. In an exome-wide association analysis of 19.4 million variants for the body mass index phenotype in 125 077 individuals (AllofUs project data), we reduced runtime from 695.35 min (11.5 h) on a single machine to 1.57 min with 30 GB of memory and 50 threads (or 8.67 min with 4 threads). Additionally, the framework supports multi-phenotype analyses, further enhancing its flexibility.

AVAILABILITY AND IMPLEMENTATION

Our optimized methods are fully integrated into PLINK 2.0 and can be accessed at: https://www.cog-genomics.org/plink/2.0/.

Efficient storage and regression computation for population-scale genome sequencing studies

In the era of big data in human genetics, large-scale biobanks aggregating genetic data from diverse populations have emerged as important for advancing our understanding of human health and disease. However, the computational and storage demands of whole genome sequencing (WGS) studies pose significant challenges, especially for researchers from underfunded institutions or developing countries, creating a disparity in research capabilities. We introduce new approaches that significantly enhance computational efficiency and reduce data storage requirements for WGS studies. By developing algorithms for compressed storage of genetic data, focusing particularly on optimizing the representation of rare variants, and designing regression methods tailored for the scale and complexity of WGS data, we significantly lower computational and storage costs. We integrate our approach into PLINK 2.0. The implementation demonstrates considerable reductions in storage space and computational time without compromising analytical accuracy, as evidenced by the application to the AllofUs project data. We optimized the runtime of an exome-wide association analysis involving 19.4 million variants and the body mass index phenotype of 125,077 individuals, reducing it from 695.35 minutes (approximately 11.5 hours) on a single machine to just 1.57 minutes using 30 GB of memory and 50 threads (or 8.67 minutes with 4 threads). Additionally, we extended this approach to support multi-phenotype analyses. We anticipate that our approach will enable researchers across the globe to unlock the potential of population biobanks, accelerating the pace of discoveries that can improve our understanding of human health and disease.

Corneal Stromal Stem Cell-Derived Extracellular Vesicles Attenuate ANGPTL7 Expression in the Human Trabecular Meshwork

Purpose

Regulating intraocular pressure (IOP), mainly via the trabecular meshwork (TM), is critical in developing glaucoma. Whereas current treatments aim to lower IOP, directly targeting the dysfunctional TM tissue for therapeutic intervention has proven challenging. In our study, we utilized Dexamethasone (Dex)-treated TM cells as a model to investigate how extracellular vesicles (EVs) from immortalized corneal stromal stem cells (imCSSCs) could influence ANGPTL7 and MYOC genes expression within TM cells.

Methods

Human TM cell lines were isolated and cultured from donor corneoscleral rims. EVs were purified from imCSSC conditioned media (CM) using size exclusion chromatography and characterized by nanoparticle tracking analysis, transmission electron microscopy (TEM), and ExoView technology. TM cells were treated with either Dex alone or with EVs for 5 days. Quantitative polymerase chain reaction (PCR) was carried out to quantify the mRNA level of MYOC and ANGPTL7.

Results

A notable increase in the expression levels of MYOC and ANGPTL7 genes was observed compared with untreated TM cells (control). Furthermore, upon comparing Dex-treated TM cells with those receiving both Dex and EV treatments, a statistically significant reduction in ANGPTL7 expression (P < 0.05) was detected.

Conclusions

The present study demonstrates that imCSSCs-derived EVs can effectively decrease the expression of ANGPLT7, a gene associated with fibrosis and implicated in the abnormal elevation of IOP in patients with glaucoma.

Translational Relevance

Our study shows that imCSSC-derived EVs can specifically target ANGPTL7 expression, making them a promising preclinical therapy for glaucoma.

Hypometric genetics: Improved power in genetic discovery by incorporating quality control flags

Balancing the tradeoff between quantity and quality of phenotypic data is critical in omics studies. Measurements below the limit of quantification (BLQ) are often tagged in quality control fields, but these flags are currently underutilized in human genetics studies. Extreme phenotype sampling is advantageous for mapping rare variant effects. We hypothesize that genetic drivers, along with environmental and technical factors, contribute to the presence of BLQ flags. Here, we introduce "hypometric genetics" (hMG) analysis and uncover a genetic basis for BLQ flags, indicating an additional source of genetic signal for genetic discovery, especially from phenotypic extremes. Applying our hMG approach to n = 227,469 UK Biobank individuals with metabolomic profiles, we reveal more than 5% heritability for BLQ flags and report biologically relevant associations, for example, at APOC3, APOA5, and PDE3B loci. For common variants, polygenic scores trained only for BLQ flags predict the corresponding quantitative traits with 91% accuracy, validating the genetic basis. For rare coding variant associations, we find an asymmetric 65.4% higher enrichment of metabolite-lowering associations for BLQ flags, highlighting the impact of putative loss-of-function variants with large effects on phenotypic extremes. Joint analysis of binarized BLQ flags and the corresponding quantitative metabolite measurements improves power in Bayesian rare variant aggregation tests, resulting in an average of 181% more prioritized genes. Our approach is broadly applicable to omics profiling. Overall, our results underscore the benefit of integrating quality control flags and quantitative measurements and highlight the advantage of joint analysis of population-based samples and phenotypic extremes in human genetics studies.

A deep catalogue of protein-coding variation in 983,578 individuals

Rare coding variants that substantially affect function provide insights into the biology of a gene1-3. However, ascertaining the frequency of such variants requires large sample sizes4-8. Here we present a catalogue of human protein-coding variation, derived from exome sequencing of 983,578 individuals across diverse populations. In total, 23% of the Regeneron Genetics Center Million Exome (RGC-ME) data come from individuals of African, East Asian, Indigenous American, Middle Eastern and South Asian ancestry. The catalogue includes more than 10.4 million missense and 1.1 million predicted loss-of-function (pLOF) variants. We identify individuals with rare biallelic pLOF variants in 4,848 genes, 1,751 of which have not been previously reported. From precise quantitative estimates of selection against heterozygous loss of function (LOF), we identify 3,988 LOF-intolerant genes, including 86 that were previously assessed as tolerant and 1,153 that lack established disease annotation. We also define regions of missense depletion at high resolution. Notably, 1,482 genes have regions that are depleted of missense variants despite being tolerant of pLOF variants. Finally, we estimate that 3% of individuals have a clinically actionable genetic variant, and that 11,773 variants reported in ClinVar with unknown significance are likely to be deleterious cryptic splice sites. To facilitate variant interpretation and genetics-informed precision medicine, we make this resource of coding variation from the RGC-ME dataset publicly accessible through a variant allele frequency browser.

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.

ANGPTL7 and Its Role in IOP and Glaucoma

Purpose

Loss-of-function variants in the ANGPTL7 gene are associated with protection from glaucoma and reduced intraocular pressure (IOP). We investigated the role of ANGPTL7 in IOP homeostasis and its potential as a target for glaucoma therapeutics.

Methods

IOP, outflow facility, and outflow tissue morphology of Angptl7 knockout (KO) mice were assessed with and without dexamethasone (Dex). ANGPTL7 was quantified in conditioned media from human trabecular meshwork cells in response to Dex, in effluent from perfused human donor eyes, and in aqueous humor from human patients treated with steroids. Antibodies to ANGPTL7 were generated and tested in three-dimensional (3D) culture of outflow cells and perfused human donor eyes. Rabbits were injected intravitreally with a neutralizing antibody targeting ANGPTL7, and IOP was measured.

Results

IOP was significantly elevated, but outflow facility and outflow tissue morphology were not different between Angptl7 KO mice and littermates. When challenged with Dex, IOP increased in wild-type but not Angptl7 KO mice. In human samples, increased ANGPTL7 was seen in the aqueous humor of patients treated with steroids, regardless of glaucoma status. Using 3D culture, recombinant ANGPTL7 decreased, and ANGPTL7-blocking antibodies increased hydraulic conductivity. Significantly, outflow facility increased in human eyes treated ex vivo with ANGPTL7-blocking antibodies, and IOP decreased for 21 days in rabbits after a single injection of blocking antibodies.

Conclusions

Using multiple models, we have demonstrated that excess ANGPTL7 increases outflow resistance and IOP and that neutralizing ANGPTL7 has beneficial effects in both naïve and steroid-induced hypertensive eyes, thus motivating the development of ANGPTL7-targeting therapeutics for the treatment of glaucoma.

Novel ancestry-specific primary open-angle glaucoma loci and shared biology with vascular mechanisms and cell proliferation

Primary open-angle glaucoma (POAG), a leading cause of irreversible blindness globally, shows disparity in prevalence and manifestations across ancestries. We perform meta-analysis across 15 biobanks (of the Global Biobank Meta-analysis Initiative) (n = 1,487,441: cases = 26,848) and merge with previous multi-ancestry studies, with the combined dataset representing the largest and most diverse POAG study to date (n = 1,478,037: cases = 46,325) and identify 17 novel significant loci, 5 of which were ancestry specific. Gene-enrichment and transcriptome-wide association analyses implicate vascular and cancer genes, a fifth of which are primary ciliary related. We perform an extensive statistical analysis of SIX6 and CDKN2B-AS1 loci in human GTEx data and across large electronic health records showing interaction between SIX6 gene and causal variants in the chr9p21.3 locus, with expression effect on CDKN2A/B. Our results suggest that some POAG risk variants may be ancestry specific, sex specific, or both, and support the contribution of genes involved in programmed cell death in POAG pathogenesis.

Novel mutation in EFEMP1 identified from two Chinese POAG families differentially activated endoplasmic reticulum stress markers and induced glaucoma in mouse

Primary open-angle glaucoma (POAG) is the most common type of glaucoma. Using whole-exome sequencing, we identified two independent families diagnosed as POAG from the China with a novel EFEMP1 variant (Exon3, c.175A>C p.Met59Leu); Three previously reported variants c.1160G>A p.R387Q, c.1189T>C p.Y397H, and c.1429C>T p.R477C in EFEPM1 from 55 sporadic POAG individuals were also identified. The variant c.175A>C p.Met59Leu co-segregated with the disease phenotype within the families. Immunoprecipitation and western blot assays showed that all three EFEMP1 mutants (p.Met59Leu, pArg140Trp, pArg345Trp) increased intracellular protein aggregations, and pMet59Leu and pArg140Arg also enhanced their extracellular proteins secretion, compared to WT in HEK293T. The differential regulations to endoplasmic reticulum (ER) stress markers ATF4, GPR78/94, and CHOP, and differential phosphorylation activations to CREB at Ser133, AKT at Ser473, p44/42 at Thr202/Tyr204, and STAT3 at Tyr705, were also detected among the mutants and WT. Finally, we revealed a significant increment of intraocular pressure and obvious reduction of RGC cells at the sixth week following intravitreal injection of adenovirus 5 (Ad5) expressing in pMet59Leu compared to WT and GFP controls. Together, variant c.175A>C p.Met59Leu in EFEMP1 is pathogenic and different mutants in EFEMP1 triggered distinct signaling pathways, explaining the reason of mutation-dependent disease phenotypes of EFEMP1.

Gene therapies and gene product-based drug candidates for normalizing and preserving tissue functions in animal models of ocular hypertension and glaucoma

More than 76 million people worldwide are afflicted with the neurodegenerative eye diseases described and grouped together as glaucoma. A common feature amongst the many forms of glaucoma is chronically elevated intraocular pressure (IOP) within the anterior chamber of the eye that physically damages the retina, optic nerve and parts of the brain connected with visual perception. The mediators of the contusing raised IOP responsible for such damage and loss of vision include locally released inflammatory agents, tissue remodeling enzymes and infiltrating immune cells which damage the retinal ganglion cell (RGC) axons and eventually kill a significant number of the RGCs. Additional culprits include genetic defects of the patient that involve aberrations in receptors, enzymes and/or endogenous ligands and possible over- or under-production of the latter. Other genetic abnormalities may include issues with signal transduction machinery within key cells of critical tissues in the front (e.g. trabecular meshwork [TM] and Schlemm's canal [SC]) and back of the eye (e.g. retinal ganglion cells and their axons). Genome-wide associated studies (GWAS) coupled with next generation sequencing have provided powerful linkage of certain gene defects and polymorphic variants to the onset and progression of diseases of the tissues involved in fluid dynamics in the TM and SC, and many retinal elements (lamina cribosa, optic nerve head) at the back of the eye which cause ocular hypertension (OHT) and glaucomatous optic neuropathy (GON), respectively. Despite the availability of some drugs, fluid drainage microshunts and full surgical techniques to lower and control intraocular pressure, the major modifiable biomarker of open-angle and other forms of glaucoma, their side-effect profiles, less than optimum effectiveness and short duration of action present opportunities to clinically manage the glaucomas with next generation of treatments with high therapeutic indices, including gene therapies. Thus, identification, characterization and deployment of genetic data coupled with traditional drug discovery and novel gene replacement, gene editing and genetic engineering technologies may provide some solutions to the aforementioned problems. These aspects will be discussed in this article.

Angiopoietin-Like Proteins: Cardiovascular Biology and Therapeutic Targeting for the Prevention of Cardiovascular Diseases

Despite the best pharmacologic tools available, cardiovascular diseases (CVDs) remain a major cause of morbidity and mortality in developed countries. After 2 decades of research, new therapeutic targets, such as angiopoietin-like proteins (ANGPTLs), are emerging. ANGPTLs belong to a family of 8 members, from ANGPTL1 to ANGPTL8; they have structural homology with angiopoietins and are secreted in the circulation. ANGPTLs display a multitude of physiological and pathologic functions; they contribute to inflammation, angiogenesis, cell death, senescence, hematopoiesis, and play a role in repair, maintenance, and tissue homeostasis. ANGPTLs-particularly the triad ANGPTL3, 4, and 8-have an established role in lipid metabolism through the regulation of triacylglycerol trafficking according to the nutritional status. Some ANGPTLs also contribute to glucose metabolism. Therefore, dysregulation in ANGPTL expression associated with abnormal circulating levels are linked to a plethora of CVD and metabolic disorders including atherosclerosis, heart diseases, diabetes, but also obesity and cancers. Because ANGPTLs bind to different receptors according to the cell type, antagonists are therapeutically inadequate. Recently, direct inhibitors of ANGPTLs, mainly ANGPTL3, have been developed, and specific monoclonal antibodies and antisense oligonucleotides are currently being tested in clinical trials. The aim of the current review is to provide an up-to-date preclinical and clinical overview on the function of the 8 members of the ANGPTL family in the cardiovascular system, their contribution to CVD, and the therapeutic potential of manipulating some of them.

Power of inclusion: Enhancing polygenic prediction with admixed individuals

Admixed individuals offer unique opportunities for addressing limited transferability in polygenic scores (PGSs), given the substantial trans-ancestry genetic correlation in many complex traits. However, they are rarely considered in PGS training, given the challenges in representing ancestry-matched linkage-disequilibrium reference panels for admixed individuals. Here we present inclusive PGS (iPGS), which captures ancestry-shared genetic effects by finding the exact solution for penalized regression on individual-level data and is thus naturally applicable to admixed individuals. We validate our approach in a simulation study across 33 configurations with varying heritability, polygenicity, and ancestry composition in the training set. When iPGS is applied to n = 237,055 ancestry-diverse individuals in the UK Biobank, it shows the greatest improvements in Africans by 48.9% on average across 60 quantitative traits and up to 50-fold improvements for some traits (neutrophil count, R2 = 0.058) over the baseline model trained on the same number of European individuals. When we allowed iPGS to use n = 284,661 individuals, we observed an average improvement of 60.8% for African, 11.6% for South Asian, 7.3% for non-British White, 4.8% for White British, and 17.8% for the other individuals. We further developed iPGS+refit to jointly model the ancestry-shared and -dependent genetic effects when heterogeneous genetic associations were present. For neutrophil count, for example, iPGS+refit showed the highest predictive performance in the African group (R2 = 0.115), which exceeds the best predictive performance for the White British group (R2 = 0.090 in the iPGS model), even though only 1.49% of individuals used in the iPGS training are of African ancestry. Our results indicate the power of including diverse individuals for developing more equitable PGS models.

Integrated Analysis of Transcriptome and Proteome of the Human Cornea and Aqueous Humor Reveal Novel Biomarkers for Corneal Endothelial Cell Dysfunction

Earlier studies have reported that elevated protein levels in the aqueous humor (AH) are associated with corneal endothelial cell dysfunction (CECD), but the details of the underlying mechanism as well as specific biomarkers for CECD remain elusive. In the present study, we aimed to identify protein markers in AH directly associated with changes to corneal endothelial cells (CECs), as AH can be easily obtained for analysis. We carried out an in-depth proteomic analysis of patient-derived AH as well as transcriptomic analysis of CECs from the same patients with bullous keratopathy (BK) resulting from CECD. We first determined differentially expressed genes (DEGs) and differentially expressed proteins (DEPs) from CECs and AH in CECD, respectively. By combining transcriptomic and proteomic analyses, 13 shared upregulated markers and 22 shared downregulated markers were observed between DEGs and DEPs. Among these 35 candidates from biomarker profiling, three upregulated markers were finally verified via data-independent acquisition (DIA) proteomic analysis using additional individual AH samples, namely metallopeptidase inhibitor 1 (TIMP1), Fc fragment of IgG binding protein (FCGBP), and angiopoietin-related protein 7 (ANGPTL7). Furthermore, we confirmed these AH biomarkers for CECD using individual immunoassay validation. Conclusively, our findings may provide valuable insights into the disease process and identify biofluid markers for the assessment of CEC function during BK development.

From target discovery to clinical drug development with human genetics

The substantial investments in human genetics and genomics made over the past three decades were anticipated to result in many innovative therapies. Here we investigate the extent to which these expectations have been met, excluding cancer treatments. In our search, we identified 40 germline genetic observations that led directly to new targets and subsequently to novel approved therapies for 36 rare and 4 common conditions. The median time between genetic target discovery and drug approval was 25 years. Most of the genetically driven therapies for rare diseases compensate for disease-causing loss-of-function mutations. The therapies approved for common conditions are all inhibitors designed to pharmacologically mimic the natural, disease-protective effects of rare loss-of-function variants. Large biobank-based genetic studies have the power to identify and validate a large number of new drug targets. Genetics can also assist in the clinical development phase of drugs-for example, by selecting individuals who are most likely to respond to investigational therapies. This approach to drug development requires investments into large, diverse cohorts of deeply phenotyped individuals with appropriate consent for genetically assisted trials. A robust framework that facilitates responsible, sustainable benefit sharing will be required to capture the full potential of human genetics and genomics and bring effective and safe innovative therapies to patients quickly.

Rare protective variants and glaucoma-relevant cell stressors modulate Angiopoietin-like 7 expression

Rare missense and nonsense variants in the Angiopoietin-like 7 (ANGPTL7) gene confer protection from primary open-angle glaucoma (POAG), though the functional mechanism remains uncharacterized. Interestingly, a larger variant effect size strongly correlates with in silico predictions of increased protein instability (r = -0.98), suggesting that protective variants lower ANGPTL7 protein levels. Here, we show that missense and nonsense variants cause aggregation of mutant ANGPTL7 protein in the endoplasmic reticulum (ER) and decreased levels of secreted protein in human trabecular meshwork (TM) cells; a lower secreted:intracellular protein ratio strongly correlates with variant effects on intraocular pressure (r = 0.81). Importantly, accumulation of mutant protein in the ER does not increase expression of ER stress proteins in TM cells (P > 0.05 for all variants tested). Cyclic mechanical stress, a glaucoma-relevant physiologic stressor, also significantly lowers ANGPTL7 expression in primary cultures of human Schlemm's canal (SC) cells (-2.4-fold-change, P = 0.01). Collectively, these data suggest that the protective effects of ANGPTL7 variants in POAG stem from lower levels of secreted protein, which may modulate responses to physiologic and pathologic ocular cell stressors. Downregulation of ANGPTL7 expression may therefore serve as a viable preventative and therapeutic strategy for this common, blinding disease.

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.

Metastasis from the tumor interior and necrotic core formation are regulated by breast cancer-derived angiopoietin-like 7

Necrosis in the tumor interior is a common feature of aggressive cancers that is associated with poor clinical prognosis and the development of metastasis. How the necrotic core promotes metastasis remains unclear. Here, we report that emergence of necrosis inside the tumor is correlated temporally with increased tumor dissemination in a rat breast cancer model and in human breast cancer patients. By performing spatially focused transcriptional profiling, we identified angiopoietin-like 7 (Angptl7) as a tumor-specific factor localized to the perinecrotic zone. Functional studies showed that Angptl7 loss normalizes central necrosis, perinecrotic dilated vessels, metastasis, and reduces circulating tumor cell counts to nearly zero. Mechanistically, Angptl7 promotes vascular permeability and supports vascular remodeling in the perinecrotic zone. Taken together, these findings show that breast tumors actively produce factors controlling central necrosis formation and metastatic dissemination from the tumor core.

Using human genetics to improve safety assessment of therapeutics

Human genetics research has discovered thousands of proteins associated with complex and rare diseases. Genome-wide association studies (GWAS) and studies of Mendelian disease have resulted in an increased understanding of the role of gene function and regulation in human conditions. Although the application of human genetics has been explored primarily as a method to identify potential drug targets and support their relevance to disease in humans, there is increasing interest in using genetic data to identify potential safety liabilities of modulating a given target. Human genetic variants can be used as a model to anticipate the effect of lifelong modulation of therapeutic targets and identify the potential risk for on-target adverse events. This approach is particularly useful for non-clinical safety evaluation of novel therapeutics that lack pharmacologically relevant animal models and can contribute to the intrinsic safety profile of a drug target. This Review illustrates applications of human genetics to safety studies during drug discovery and development, including assessing the potential for on- and off-target associated adverse events, carcinogenicity risk assessment, and guiding translational safety study designs and monitoring strategies. A summary of available human genetic resources and recommended best practices is provided. The challenges and future perspectives of translating human genetic information to identify risks for potential drug effects in preclinical and clinical development are discussed.

Research Note: Expression levels of collagen-related genes in PSE conditions and white striping features of broiler pectoralis major muscle

The objective of this study was to compare the meat quality characteristics and collagen-related gene expression levels in the broiler pectoralis major (PM) muscle among the meat quality groups, including normal, pale, soft, and exudative (PSE), and white striping (WS) groups. The group was classified by their WS degree (moderate or severe striping) and quality traits including pH and lightness values at 24 h postmortem (normal group: pH_{24 h} ≥5.7, 48≤L*≤53, without WS features; PSE group: pH_{24 h} <5.7, L* >53, without WS features; WS group: pH_{24 h} ≥5.7, 48≤L*≤53, with moderate or severe striping). The WS group revealed no differences in all measured meat quality traits compared to the normal group (P > 0.05). PM muscles exhibiting PSE conditions without WS indicated lower pH_{15 min} and pH_{24 h} values (P < 0.05). Whereas, lower lightness and cooking loss values were observed in the normal and WS groups compared to the PSE group (P < 0.05). No significant difference was observed in the level of type I collagen among the groups (P > 0.05), whereas a higher type III collagen level was observed in the WS group than in the other groups (P < 0.05). Additionally, the WS group showed a higher type IV collagen level compared to the normal group (P < 0.05) and a level not different from that of the PSE group (P > 0.05). In contrast, the expression levels of matrix metalloproteinase (MMP) 2, involved in type IV collagen degradation, and angiopoietin-like protein 7, associated with collagen accumulation, were higher in the WS group compared to the normal group (P < 0.05). However, no difference was detected in the MMP1 level among the all groups (P > 0.05). These results suggest that the occurrence of WS features in broiler PM muscle, unlike PSE and normal conditions, can be influenced by the expression levels of collagen-related genes associated with abnormalities in extracellular matrix components.

FinnGen provides genetic insights from a well-phenotyped isolated population

Population isolates such as those in Finland benefit genetic research because deleterious alleles are often concentrated on a small number of low-frequency variants (0.1% ≤ minor allele frequency < 5%). These variants survived the founding bottleneck rather than being distributed over a large number of ultrarare variants. Although this effect is well established in Mendelian genetics, its value in common disease genetics is less explored1,2. FinnGen aims to study the genome and national health register data of 500,000 Finnish individuals. Given the relatively high median age of participants (63 years) and the substantial fraction of hospital-based recruitment, FinnGen is enriched for disease end points. Here we analyse data from 224,737 participants from FinnGen and study 15 diseases that have previously been investigated in large genome-wide association studies (GWASs). We also include meta-analyses of biobank data from Estonia and the United Kingdom. We identified 30 new associations, primarily low-frequency variants, enriched in the Finnish population. A GWAS of 1,932 diseases also identified 2,733 genome-wide significant associations (893 phenome-wide significant (PWS), P < 2.6 × 10-11) at 2,496 (771 PWS) independent loci with 807 (247 PWS) end points. Among these, fine-mapping implicated 148 (73 PWS) coding variants associated with 83 (42 PWS) end points. Moreover, 91 (47 PWS) had an allele frequency of <5% in non-Finnish European individuals, of which 62 (32 PWS) were enriched by more than twofold in Finland. These findings demonstrate the power of bottlenecked populations to find entry points into the biology of common diseases through low-frequency, high impact variants.

ANGPTL7, a therapeutic target for increased intraocular pressure and glaucoma

Glaucoma is a leading cause of blindness. Current glaucoma medications work by lowering intraocular pressure (IOP), a risk factor for glaucoma, but most treatments do not directly target the pathological changes leading to increased IOP, which can manifest as medication resistance as disease progresses. To identify physiological modulators of IOP, we performed genome- and exome-wide association analysis in >129,000 individuals with IOP measurements and extended these findings to an analysis of glaucoma risk. We report the identification and functional characterization of rare coding variants (including loss-of-function variants) in ANGPTL7 associated with reduction in IOP and glaucoma protection. We validated the human genetics findings in mice by establishing that Angptl7 knockout mice have lower (~2 mmHg) basal IOP compared to wild-type, with a trend towards lower IOP also in heterozygotes. Conversely, increasing murine Angptl7 levels via injection into mouse eyes increases the IOP. We also show that acute Angptl7 silencing in adult mice lowers the IOP (~2-4 mmHg), reproducing the observations in knockout mice. Collectively, our data suggest that ANGPTL7 is important for IOP homeostasis and is amenable to therapeutic modulation to help maintain a healthy IOP that can prevent onset or slow the progression of glaucoma.

Systematic single-variant and gene-based association testing of thousands of phenotypes in 394,841 UK Biobank exomes

Genome-wide association studies have successfully discovered thousands of common variants associated with human diseases and traits, but the landscape of rare variations in human disease has not been explored at scale. Exome-sequencing studies of population biobanks provide an opportunity to systematically evaluate the impact of rare coding variations across a wide range of phenotypes to discover genes and allelic series relevant to human health and disease. Here, we present results from systematic association analyses of 4,529 phenotypes using single-variant and gene tests of 394,841 individuals in the UK Biobank with exome-sequence data. We find that the discovery of genetic associations is tightly linked to frequency and is correlated with metrics of deleteriousness and natural selection. We highlight biological findings elucidated by these data and release the dataset as a public resource alongside the Genebass browser for rapidly exploring rare-variant association results.

Short and long-term effect of dexamethasone on the transcriptome profile of primary human trabecular meshwork cells in vitro

In the quest of identifying newer molecular targets for the management of glucocorticoid-induced ocular hypertension (GC-OHT) and glaucoma (GCG), several microarray studies have attempted to investigate the genome-wide transcriptome profiling of primary human trabecular meshwork (TM) cells in response to dexamethasone (DEX). However, no studies are reported so far to demonstrate the temporal changes in the expression of genes in the cultured human TM cells in response to DEX treatment. Therefore, in the present study, the time-dependent changes in the genome-wide expression of genes in primary human TM cells after short (16 hours: 16 h) and long exposure (7 days: 7 d) of DEX was investigated using RNA sequencing. There were 199 (118 up-regulated; 81 down-regulated) and 525 (119 up-regulated; 406 down-regulated) DEGs in 16 h and 7 d treatment groups respectively. The unique genes identified in 16 h and 7 d treatment groups were 152 and 478 respectively. This study found a distinct gene signature and pathways between two treatment regimes. Longer exposure of DEX treatment showed a dys-regulation of Wnt and Rap1 signaling and so highlighted potential therapeutic targets for pharmacological management of GC-OHT/glaucoma.

The genetic basis for adult onset glaucoma: Recent advances and future directions

Glaucoma, a diverse group of eye disorders that results in the degeneration of retinal ganglion cells, is the world's leading cause of irreversible blindness. Apart from age and ancestry, the major risk factor for glaucoma is increased intraocular pressure (IOP). In primary open-angle glaucoma (POAG), the anterior chamber angle is open but there is resistance to aqueous outflow. In primary angle-closure glaucoma (PACG), crowding of the anterior chamber angle due to anatomical alterations impede aqueous drainage through the angle. In exfoliation syndrome and exfoliation glaucoma, deposition of white flaky material throughout the anterior chamber directly interfere with aqueous outflow. Observational studies have established that there is a strong hereditable component for glaucoma onset and progression. Indeed, a succession of genome wide association studies (GWAS) that were centered upon single nucleotide polymorphisms (SNP) have yielded more than a hundred genetic markers associated with glaucoma risk. However, a shortcoming of GWAS studies is the difficulty in identifying the actual effector genes responsible for disease pathogenesis. Building on the foundation laid by GWAS studies, research groups have recently begun to perform whole exome-sequencing to evaluate the contribution of protein-changing, coding sequence genetic variants to glaucoma risk. The adoption of this technology in both large population-based studies as well as family studies are revealing the presence of novel, protein-changing genetic variants that could enrich our understanding of the pathogenesis of glaucoma. This review will cover recent advances in the genetics of primary open-angle glaucoma, primary angle-closure glaucoma and exfoliation glaucoma, which collectively make up the vast majority of all glaucoma cases in the world today. We will discuss how recent advances in research methodology have uncovered new risk genes, and how follow up biological investigations could be undertaken in order to define how the risk encoded by a genetic sequence variant comes into play in patients. We will also hypothesise how data arising from characterising these genetic variants could be utilized to predict glaucoma risk and the manner in which new therapeutic strategies might be informed.

NEDD9 Mediates the FAK/Src Signaling Pathway to Promote the Adhesion of Human Trabecular Meshwork Cells after Dexamethasone Treatment

PURPOSE

The differential gene expression of neural precursor cell expressed developmentally downregulated 9 (NEDD9) in human trabecular meshwork (HTM) cells after dexamethasone (Dex) treatment was confirmed through gene expression profiling. However, the regulatory mechanism of NEDD9 expression in HTM cells remains unknown. In this study, we investigated NEDD9 expression in HTM cells and gained a better understanding of glucocorticoid-induced glaucoma (GIG) pathophysiology.

METHODS

The Gene Expression Omnibus database and GEO2R tool were used to identify differentially expressed genes in the GSE37474 and GSE124114 datasets, and NEDD9 gene expression was found to be upregulated. Human corneoscleral segments and HTM cells were treated with 100 nM Dex or an equal volume of ethanol (0.01%) for 7 days. NEDD9 expression in TM tissues was evaluated by immunohistochemistry, and NEED9 expression in HTM cells was confirmed by RT-qPCR and western blotting. HTM cell adhesive behaviors were assessed with a cell adhesion detection kit. NEDD9 expression was knocked down with short hairpin RNA in HTM cells, and FAK/Src signaling pathway activation was found to be regulated by NEDD9.

RESULTS

After 7 days of HTM cell Dex treatment, NEDD9 expression was upregulated to be approximately twice that of control. FAK, Src, phospho-FAK, and phospho-Src expression in Dex-treated HTM cells was markedly increased. Downregulation of NEDD9 expression reduced HTM cell adhesion to the surface of culture wells and simultaneously led to a reduction in FAK, Src, phospho-FAK and phospho-Src expression.

CONCLUSIONS

NEDD9 expression is upregulated in HTM cells after Dex treatment and promotes HTM cell adhesion. These findings underscore the contribution of NEDD9 overexpression to altered HTM cell adhesion during glucocorticoid therapy and may play a key role in GIG pathological progression. Considering the similarity between GIG and primary open-angle glaucoma (POAG), our findings suggest that targeting NEDD9 may be a new therapeutic strategy for POAG patients.

Cox regression is robust to inaccurate EHR-extracted event time: an application to EHR-based GWAS

MOTIVATION

Logistic regression models are used in genomic studies to analyze the genetic data linked to electronic health records (EHRs), and do not take full usage of the time-to-event information available in EHRs. Previous work has shown that Cox regression, which can account for left truncation and right censoring in EHRs, increased the power to detect genotype-phenotype associations compared to logistic regression. We extend this to evaluate the relative performance of Cox regression and various logistic regression models in the presence of positive errors in event time (delayed event time), relating to recorded event time accuracy.

RESULTS

One Cox model and three logistic regression models were considered under different scenarios of delayed event time. Extensive simulations and a genomic study application were used to evaluate the impact of delayed event time. While logistic regression does not model the time-to-event directly, various logistic regression models used in the literature were more sensitive to delayed event time than Cox regression. Results highlighted the importance to identify and exclude the patients diagnosed before entry time. Cox regression had similar or modest improvement in statistical power over various logistic regression models at controlled type I error. This was supported by the empirical data, where the Cox models steadily had the highest sensitivity to detect known genotype-phenotype associations under all scenarios of delayed event time.

AVAILABILITY AND IMPLEMENTATION

Access to individual-level EHR and genotype data is restricted by the IRB. Simulation code and R script for data process are at: https://github.com/QingxiaCindyChen/CoxRobustEHR.git.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The genetics of glaucoma: Disease associations, personalised risk assessment and therapeutic opportunities-A review

Glaucoma refers to a heterogenous group of disorders characterised by progressive loss of retinal ganglion cells and associated visual field loss. Both early-onset and adult-onset forms of the disease have a strong genetic component. Here, we summarise the known genetic associations for various forms of glaucoma and the possible functional roles for these genes in disease pathogenesis. We also discuss efforts to translate genetic knowledge into clinical practice, including gene-based tests for disease diagnosis and risk-stratification as well as gene-based therapies.

Significant sparse polygenic risk scores across 813 traits in UK Biobank

We present a systematic assessment of polygenic risk score (PRS) prediction across more than 1,500 traits using genetic and phenotype data in the UK Biobank. We report 813 sparse PRS models with significant (p < 2.5 x 10-5) incremental predictive performance when compared against the covariate-only model that considers age, sex, types of genotyping arrays, and the principal component loadings of genotypes. We report a significant correlation between the number of genetic variants selected in the sparse PRS model and the incremental predictive performance (Spearman's ⍴ = 0.61, p = 2.2 x 10-59 for quantitative traits, ⍴ = 0.21, p = 9.6 x 10-4 for binary traits). The sparse PRS model trained on European individuals showed limited transferability when evaluated on non-European individuals in the UK Biobank. We provide the PRS model weights on the Global Biobank Engine (https://biobankengine.stanford.edu/prs).

ANGPTL7 is transcriptionally regulated by SP1 and modulates glucocorticoid-induced cross-linked actin networks in trabecular meshwork cells via the RhoA/ROCK pathway

Glaucoma is one of the leading causes of worldwide irreversible blindness. Lowering elevated intraocular pressure (IOP) is currently the only effective approach for controlling the progress of glaucoma. Angiopoietin-like 7 (ANGPTL7) takes a key part in elevated outflow resistance of aqueous humor in dysfunctional trabecular meshwork (TM), along with the formation of cross-linked actin networks (CLANs), leading to high IOP. In this study, we explored the role of the ANGPTL7 signaling pathway in CLAN formation. We detected the expression of ANGPTL7 in cultured primary TM cells treated with dexamethasone (DEX) and ethanol as a control using qRT-PCR and western blotting. Actin filaments were revealed by phalloidin staining. ANGPTL7 short hairpin RNA (shRNA) was applied to TM cells to examine the effect of ANGPTL7 on DEX-induced CLAN formation. Western blotting was used to assess the effect of ANGPTL7 on the RhoA/Rho-associated kinase (Rho-kinase/ROCK) signaling pathway. Bioinformatics, dual-luciferase reporter assays, and chromatin immunoprecipitation were employed to identify the transcription factors of ANGPTL7. Transcription factor specificity protein 1 (SP1) overexpression and silencing were performed to determine their roles in the modulation of ANGPTL7 expression. We found DEX-induced ANGPTL7 expression and stress fiber rearrangement in TM cells. ANGPTL7 knockdown effectively inhibited the formation of CLANs. Moreover, it was involved in the regulation of the RhoA/ROCK signaling pathway, further affecting DEX-induced CLAN formation. SP1 was identified as a transcription factor of ANGPTL7 which regulated ANGPTL7 level to mediate CLAN formation through the RhoA/ROCK signaling pathway. This study contributes to revealing the molecular mechanisms of ANGPTL7 in CLAN formation, which is involved in TM dysfunction and glaucoma pathogenesis.

EFEMP1 rare variants cause familial juvenile-onset open-angle glaucoma

Juvenile open-angle glaucoma (JOAG) is a severe type of glaucoma with onset before age 40 and dominant inheritance. Using exome sequencing we identified 3 independent families from the Philippines with novel EFEMP1 variants (c.238A>T, p.Asn80Tyr; c.1480T>C, p.Ter494Glnext*29; and c.1429C>T, p.Arg477Cys) co-segregating with disease. Affected variant carriers (N = 34) exhibited severe disease with average age of onset of 16 years and with 76% developing blindness. To investigate functional effects, we transfected COS7 cells with vectors expressing the three novel EFEMP1 variants and showed that all three variants found in JOAG patients caused significant intracellular protein aggregation and retention compared to wild type and also compared to EFEMP1 variants associated with other ocular phenotypes including an early-onset form of macular degeneration, Malattia Leventinese/Doyne's Honeycomb retinal dystrophy. These results suggest that rare EFEMP1 coding variants can cause JOAG through a mechanism involving protein aggregation and retention, and that the extent of intracellular retention correlates with disease phenotype. This is the first report of EFEMP1 variants causing JOAG, expanding the EFEMP1 disease spectrum. Our results suggest that EFEMP1 mutations appear to be a relatively common cause of JOAG in Filipino families, an ethnically diverse population.

Association of the MYOC p.(Gln368Ter) Variant With Glaucoma in a Finnish Population

Importance

The c.1102C>T, p.(Gln368Ter) variant in the myocilin (MYOC) gene is a known risk allele for glaucoma. It is the most common MYOC risk variant for glaucoma among individuals of European ancestry, and its prevalence is highest in Finland. Furthermore, exfoliation syndrome has high prevalence in Scandinavia, making the Finnish population ideal to study the association of the variant with different types of glaucoma.

Objectives

To examine the association and penetrance of MYOC p.(Gln368Ter) (rs74315329) variant with different types of glaucoma in a Finnish population.

Design, Setting, and Participants

This genetic association study included individuals of Finnish ancestry in the FinnGen project. The participants were collected from Finnish biobanks, and the disease end points were defined using nationwide registries. The MYOC c.1102C>T variant was either directly genotyped or imputed with microarrays. Recruitment of samples to FinnGen was initiated in 2017, and data analysis was performed between December 2019 and May 2020.

Main Outcomes and Measures

The main outcomes were odds ratios (ORs) and penetrance with different types of glaucoma and in different age groups.

Results

A total of 218 792 individuals were included in this study (mean [SD] age 52.4 [17.5] years; 123 579 women [56.5%]), including 8591 (3.9%) with glaucoma, 3412 (1.6%) with primary open-angle glaucoma, 1515 (0.7%) with exfoliation glaucoma, 892 (0.4%) with normal-tension glaucoma, and 4766 (2.2%) with suspected glaucoma. The minor allele frequency of MYOC p.(Gln368Ter) was 0.28%. Individuals with the heterozygous variant had higher odds of primary open-angle glaucoma (OR, 3.36; 95% CI, 2.55-4.37), overall glaucoma (OR, 2.58; 95% CI, 2.12-3.13), suspected glaucoma (OR, 2.53; 95% CI, 1.93-3.26), exfoliation glaucoma (OR, 2.61; 95% CI, 1.60-4.02), and undergoing glaucoma-related operations (OR, 5.45; 95% CI, 2.95-9.28). The penetrance of heterozygous MYOC p.(Gln368Ter) was 5.2% in individuals with primary open-angle glaucoma, 9.6% in individuals with glaucoma, 5.4% in individuals with suspected glaucoma, and 1.9% in individuals with exfoliation glaucoma. There was no significant association with normal-tension glaucoma (OR, 1.69; 95% CI, 0.72-3.35).

Conclusions and Relevance

This genetic association study found that the MYOC p.(Gln368Ter) variant was associated with exfoliation glaucoma. The association with normal-tension glaucoma could not be replicated. These findings suggest that MYOC p.(Gln368Ter) was associated with open-angle glaucoma and exfoliation glaucoma in a Finnish population.

Extraocular, periocular, and intraocular routes for sustained drug delivery for glaucoma

Although once daily anti-glaucoma drug therapy is a current clinical reality, most therapies require multiple dosing and there is an unmet need to develop convenient, safe, and effective sustained release drug delivery systems for long-term treatment to improve patient adherence and outcomes. One of the first sustained release drug delivery systems was approved for the reduction of intraocular pressure in glaucoma patients. It is a polymeric reservoir-type insert delivery system, Ocusert™, placed under the eyelid and on the ocular surface for zero-order drug release over one week. The insert, marketed in two strengths, released pilocarpine on the eye surface. While many clinicians appreciated this drug product, it was eventually discontinued. No similar sustained release non-invasive drug delivery system has made it to the market to date for treating glaucoma. Drug delivery systems under development include punctal plugs, ring-type systems, contact lenses, implants, microspheres, nanospheres, gels, and other depot systems placed in the extraocular, periocular, or intraocular regions including intracameral, supraciliary, and intravitreal spaces. This article discusses the advantages and disadvantages of the various routes of administration and delivery systems for sustained glaucoma therapy. It also provides the reader with some examples and discussion of drug delivery systems that could potentially be applied for glaucoma treatment. Interestingly, one intracamerally injected implant, Durysta™, was approved recently for sustained intraocular pressure reduction. However, long-term acceptance of such devices has yet to be established. The ultimate success of the delivery system will depend on efficacy relative to eye drop dosing, safety, reimbursement options, and patient acceptance. Cautious development efforts are warranted considering prior failed approaches for sustained glaucoma drug delivery. Neuroprotective approaches for glaucoma therapy including cell, gene, protein, and drug-combination therapies, mostly administered intravitreally, are also rapidly progressing towards assessment in humans.

The GGLEAM Study: Understanding Glaucoma in the Ohio Amish

Glaucoma leads to millions of cases of visual impairment and blindness around the world. Its susceptibility is shaped by both environmental and genetic risk factors. Although over 120 risk loci have been identified for glaucoma, a large portion of its heritability is still unexplained. Here we describe the foundation of the Genetics of GLaucoma Evaluation in the AMish (GGLEAM) study to investigate the genetic architecture of glaucoma in the Ohio Amish, which exhibits lower genetic and environmental heterogeneity compared to the general population. To date, we have enrolled 81 Amish individuals in our study from Holmes County, Ohio. As a part of our enrollment process, 62 GGLEAM study participants (42 glaucoma-affected and 20 unaffected individuals) received comprehensive eye examinations and glaucoma evaluations. Using the data from the Anabaptist Genealogy Database, we found that 80 of the GGLEAM study participants were related to one another through a large, multigenerational pedigree containing 1586 people. We plan to integrate the health and kinship data obtained for the GGLEAM study to interrogate glaucoma genetics and pathophysiology in this unique population.