Bioinformatic Prioritization and Functional Annotation of GWAS-Based Candidate Genes for Primary Open-Angle Glaucoma

Utilizing systems genetics to enhance understanding into molecular targets of skin cancer

Despite progress made with immune checkpoint inhibitors and targeted therapies, skin cancer remains a significant public health concern in the United States. The intricacies of the disease, encompassing genetics, immune responses, and external factors, call for a comprehensive approach. Techniques in systems genetics, including transcriptional correlation analysis, functional pathway enrichment analysis, and protein-protein interaction network analysis, prove valuable in deciphering intricate molecular mechanisms and identifying potential diagnostic and therapeutic targets for skin cancer. Recent studies demonstrate the efficacy of these techniques in uncovering molecular processes and pinpointing diagnostic markers for various skin cancer types, highlighting the potential of systems genetics in advancing innovative therapies. While certain limitations exist, such as generalizability and contextualization of external factors, the ongoing progress in AI technologies provides hope in overcoming these challenges. By providing protocols and a practical example involving Braf, we aim to inspire early-career experimental dermatologists to adopt these tools and seamlessly integrate these techniques into their skin cancer research, positioning them at the forefront of innovative approaches in combating this devastating disease.

A genome-wide association meta-analysis implicates Hedgehog and Notch signaling in Dupuytren's disease

Dupuytren's disease (DD) is a highly heritable fibrotic disorder of the hand with incompletely understood etiology. A number of genetic loci, including Wnt signaling members, have been previously identified. Our overall aim was to identify novel genetic loci, to prioritize genes within the loci for functional studies, and to assess genetic correlation with associated disorders. We performed a meta-analysis of six DD genome-wide association studies from three European countries and extensive bioinformatic follow-up analyses. Leveraging 11,320 cases and 47,023 controls, we identified 85 genome-wide significant single nucleotide polymorphisms in 56 loci, of which 11 were novel, explaining 13.3-38.1% of disease variance. Gene prioritization implicated the Hedgehog and Notch signaling pathways. We also identified a significant genetic correlation with frozen shoulder. The pathways identified highlight the potential for new therapeutic targets and provide a basis for additional mechanistic studies for a common disorder that can severely impact hand function.

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.

Understanding the complex genetics and molecular mechanisms underlying glaucoma

Glaucoma is the leading cause of irreversible blindness worldwide. Currently the only effective treatment for glaucoma is to reduce the intraocular pressure, which can halt the progression of the disease. Highlighting the importance of identifying individuals at risk of developing glaucoma and those with early-stage glaucoma will help patients receive treatment before sight loss. However, some cases of glaucoma do not have raised intraocular pressure. In fact, glaucoma is caused by a variety of different mechanisms and has a wide range of different subtypes. Understanding other risk factors, the underlying mechanisms, and the pathology of glaucoma might lead to novel treatments and treatment of underlying diseases. In this review we present the latest research into glaucoma including the genetics and molecular basis of the disease.

Computational methods in glaucoma research: Current status and future outlook

Advancements in computational techniques have transformed glaucoma research, providing a deeper understanding of genetics, disease mechanisms, and potential therapeutic targets. Systems genetics integrates genomic and clinical data, aiding in identifying drug targets, comprehending disease mechanisms, and personalizing treatment strategies for glaucoma. Molecular dynamics simulations offer valuable molecular-level insights into glaucoma-related biomolecule behavior and drug interactions, guiding experimental studies and drug discovery efforts. Artificial intelligence (AI) technologies hold promise in revolutionizing glaucoma research, enhancing disease diagnosis, target identification, and drug candidate selection. The generalized protocols for systems genetics, MD simulations, and AI model development are included as a guide for glaucoma researchers. These computational methods, however, are not separate and work harmoniously together to discover novel ways to combat glaucoma. Ongoing research and progresses in genomics technologies, MD simulations, and AI methodologies project computational methods to become an integral part of glaucoma research in the future.

Regional Gene Expression in the Retina, Optic Nerve Head, and Optic Nerve of Mice with Optic Nerve Crush and Experimental Glaucoma

A major risk factor for glaucomatous optic neuropathy is the level of intraocular pressure (IOP), which can lead to retinal ganglion cell axon injury and cell death. The optic nerve has a rostral unmyelinated portion at the optic nerve head followed by a caudal myelinated region. The unmyelinated region is differentially susceptible to IOP-induced damage in rodent models and human glaucoma. While several studies have analyzed gene expression changes in the mouse optic nerve following optic nerve injury, few were designed to consider the regional gene expression differences that exist between these distinct areas. We performed bulk RNA-sequencing on the retina and separately micro-dissected unmyelinated and myelinated optic nerve regions from naïve C57BL/6 mice, mice after optic nerve crush, and mice with microbead-induced experimental glaucoma (total = 36). Gene expression patterns in the naïve unmyelinated optic nerve showed significant enrichment of the Wnt, Hippo, PI3K-Akt, and transforming growth factor β pathways, as well as extracellular matrix-receptor and cell membrane signaling pathways, compared to the myelinated optic nerve and retina. Gene expression changes induced by both injuries were more extensive in the myelinated optic nerve than the unmyelinated region, and greater after nerve crush than glaucoma. Changes present three and fourteen days after injury largely subsided by six weeks. Gene markers of reactive astrocytes did not consistently differ between injury states. Overall, the transcriptomic phenotype of the mouse unmyelinated optic nerve was significantly different from immediately adjacent tissues, likely dominated by expression in astrocytes, whose junctional complexes are inherently important in responding to IOP elevation.

On the interpretation of transcriptome-wide association studies

Transcriptome-wide association studies (TWAS) aim to detect relationships between gene expression and a phenotype, and are commonly used for secondary analysis of genome-wide association study (GWAS) results. Results from TWAS analyses are often interpreted as indicating a genetic relationship between gene expression and a phenotype, but this interpretation is not consistent with the null hypothesis that is evaluated in the traditional TWAS framework. In this study we provide a mathematical outline of this TWAS framework, and elucidate what interpretations are warranted given the null hypothesis it actually tests. We then use both simulations and real data analysis to assess the implications of misinterpreting TWAS results as indicative of a genetic relationship between gene expression and the phenotype. Our simulation results show considerably inflated type 1 error rates for TWAS when interpreted this way, with 41% of significant TWAS associations detected in the real data analysis found to have insufficient statistical evidence to infer such a relationship. This demonstrates that in current implementations, TWAS cannot reliably be used to investigate genetic relationships between gene expression and a phenotype, but that local genetic correlation analysis can serve as a potential alternative.

Elevated Intraocular Pressure and Glaucomatous Optic Neuropathy: Genes to Disease Mechanisms, Therapeutic Drugs, and Gene Therapies

This review article focuses on the pathogenesis of and genetic defects linked with chronic ocular hypertension (cOHT) and glaucoma. The latter ocular disease constitutes a group of ocular degenerative diseases whose hallmark features are damage to the optic nerve, apoptotic demise of retinal ganglion cells, disturbances within the brain regions involved in visual perception and considerable visual impairment that can lead to blindness. Even though a number of pharmaceuticals, surgical and device-based treatments already exist addressing cOHT associated with the most prevalent of the glaucoma types, primary open-angle glaucoma (POAG), they can be improved upon in terms of superior efficacy with reduced side-effects and with longer duration of activity. The linkage of disease pathology to certain genes via genome-wide associated studies are illuminating new approaches to finding novel treatment options for the aforementioned ocular disorders. Gene replacement, gene editing via CRISPR-Cas9, and the use of optogenetic technologies may replace traditional drug-based therapies and/or they may augment existing therapeutics for the treatment of cOHT and POAG in the future.

Regional Gene Expression in the Retina, Optic Nerve Head, and Optic Nerve of Mice with Experimental Glaucoma and Optic Nerve Crush

A major risk factor for glaucomatous optic neuropathy is the level of intraocular pressure (IOP), which can lead to retinal ganglion cell axon injury and cell death. The optic nerve has a rostral unmyelinated portion at the optic nerve head followed by a caudal myelinated region. The unmyelinated region is differentially susceptible to IOP-induced damage in rodent models and in human glaucoma. While several studies have analyzed gene expression changes in the mouse optic nerve following optic nerve injury, few were designed to consider the regional gene expression differences that exist between these distinct areas. We performed bulk RNA-sequencing on the retina and on separately micro-dissected unmyelinated and myelinated optic nerve regions from naïve C57BL/6 mice, mice after optic nerve crush, and mice with microbead-induced experimental glaucoma (total = 36). Gene expression patterns in the naïve unmyelinated optic nerve showed significant enrichment of the Wnt, Hippo, PI3K-Akt, and transforming growth factor β pathways, as well as extracellular matrix-receptor and cell membrane signaling pathways, compared to the myelinated optic nerve and retina. Gene expression changes induced by both injuries were more extensive in the myelinated optic nerve than the unmyelinated region, and greater after nerve crush than glaucoma. Changes three and fourteen days after injury largely subsided by six weeks. Gene markers of reactive astrocytes did not consistently differ between injury states. Overall, the transcriptomic phenotype of the mouse unmyelinated optic nerve was significantly different from immediately adjacent tissues, likely dominated by expression in astrocytes, whose junctional complexes are inherently important in responding to IOP elevation.

All-trans Retinoic Acids Synergistically and Beneficially Affect In Vitro Glaucomatous Trabecular Meshwork (TM) Models Using 2D and 3D Cell Cultures of Human TM Cells

We report herein on the effects of all-trans retinoic acid (ATRA) on two-dimensional (2D) and three-dimensional (3D) cultures of human trabecular meshwork (HTM) cells that were treated with transforming growth factor β2 (TGF-β2). In the presence of 5 ng/mL TGF-β2, the effects of ATRA on the following were observed: (1) the barrier function of the 2D HTM monolayers, as determined by trans-endothelial electrical resistance (TEER) and fluorescein isothiocyanate (FITC) dextran permeability measurements; (2) a Seahorse cellular bio-metabolism analysis; (3) physical properties, including the size and stiffness, of 3D spheroids; (4) the gene expression of extracellular matrix (ECM) molecules, ECM modulators including tissue inhibitor of metalloproteinases (TIMPs), matrix metalloproteinases (MMPs), tight junction (TJ)-related molecules, and endoplasmic reticulum (ER)-stress-related factors. ATRA significantly inhibited the TGF-β2-induced increase in the TEER values and FITC dextran permeability of the 2D monolayers, while an ATRA monotreatment induced similar effects as TGF-β2. A real-time metabolic analysis revealed that ATRA significantly inhibited the TGF-β2-induced shift in metabolic reserve from mitochondrial oxidative phosphorylation to glycolysis in 2D HTM cells, whereas ATRA alone did not induce significant metabolic changes. In contrast, ATRA induced the formation of substantially downsized and softer 3D spheroids in the absence and presence of TGF-β2. The different effects induced by ATRA toward 2D and 3D HTM cells were also supported by the qPCR analysis of several proteins as above. The findings reported here indicate that ATRA may induce synergistic and beneficial effects on TGF-β2-treated 2D- and 3D-cultured HTM cells; those effects varied significantly between the 2D and 3D cultures.

Large-Scale Multi-Omics Studies Provide New Insights into Blood Pressure Regulation

Recent genome-wide association studies uncovered part of blood pressure’s heritability. However, there is still a vast gap between genetics and biology that needs to be bridged. Here, we followed up blood pressure genome-wide summary statistics of over 750,000 individuals, leveraging comprehensive epigenomic and transcriptomic data from blood with a follow-up in cardiovascular tissues to prioritise likely causal genes and underlying blood pressure mechanisms. We first prioritised genes based on coding consequences, multilayer molecular associations, blood pressure-associated expression levels, and coregulation evidence. Next, we followed up the prioritised genes in multilayer studies of genomics, epigenomics, and transcriptomics, functional enrichment, and their potential suitability as drug targets. Our analyses yielded 1880 likely causal genes for blood pressure, tens of which are targets of the available licensed drugs. We identified 34 novel genes for blood pressure, supported by more than one source of biological evidence. Twenty-eight (82%) of these new genes were successfully replicated by transcriptome-wide association analyses in a large independent cohort (n = ~220,000). We also found a substantial mediating role for epigenetic regulation of the prioritised genes. Our results provide new insights into genetic regulation of blood pressure in terms of likely causal genes and involved biological pathways offering opportunities for future translation into clinical practice.