ForestPMPlot: A Flexible Tool for Visualizing Heterogeneity Between Studies in Meta-analysis.

Meta-analysis has become a popular tool for genetic association studies to combine different genetic studies. A key challenge in meta-analysis is heterogeneity, or the differences in effect sizes between studies. Heterogeneity complicates the interpretation of meta-analyses. In this paper, we describe ForestPMPlot, a flexible visualization tool for analyzing studies included in a meta-analysis. The main feature of the tool is visualizing the differences in the effect sizes of the studies to understand why the studies exhibit heterogeneity for a particular phenotype and locus pair under different conditions. We show the application of this tool to interpret a meta-analysis of 17 mouse studies, and to interpret a multi-tissue eQTL study

High-Resolution Mutation Mapping Reveals Parallel Experimental Evolution in Yeast

Understanding the genetic basis of evolutionary adaptation is limited by our ability to efficiently identify the genomic locations of adaptive mutations. Here we describe a method that can quickly and precisely map the genetic basis of naturally and experimentally evolved complex traits using linkage analysis. A yeast strain that expresses the evolved trait is crossed to a distinct strain background and DNA from a large pool of progeny that express the trait of interest is hybridized to oligonucleotide microarrays that detect thousands of polymorphisms between the two strains. Adaptive mutations are detected by linkage to the polymorphisms from the evolved parent. We successfully tested our method by mapping five known genes to a precision of 0.2–24 kb (0.1–10 cM), and developed computer simulations to test the effect of different factors on mapping precision. We then applied this method to four yeast strains that had independently adapted to a fluctuating glucose–galactose environment. All four strains had acquired one or more missense mutations in GAL80, the repressor of the galactose utilization pathway. When transferred into the ancestral strain, the gal80 mutations conferred the fitness advantage that the evolved strains show in the transition from glucose to galactose. Our results show an example of parallel adaptation caused by mutations in the same gene

Machine learning-derived baseline visual field patterns predict future glaucoma onset in the Ocular Hypertension Treatment Study

PURPOSE: The Ocular Hypertension Treatment Study (OHTS) identified risk factors for primary open-angle glaucoma (POAG) in patients with ocular hypertension, including pattern standard deviation (PSD). Archetypal analysis, an unsupervised machine learning method, may offer a more interpretable approach to risk stratification by identifying patterns in baseline visual fields (VFs). METHODS: There were 3272 eyes available in the OHTS. Archetypal analysis was applied using 24-2 baseline VFs, and model selection was performed with cross-validation. Decomposition coefficients for archetypes (ATs) were calculated. A penalized Cox proportional hazards model was implemented to select discriminative ATs. The AT model was compared to the OHTS model. Associations were identified between ATs with both POAG onset and VF progression, defined by mean deviation change per year. RESULTS: We selected 8494 baseline VFs. Optimal AT count was 19. The highest prevalence ATs were AT9, AT11, and AT7. The AT-based prediction model had a C-index of 0.75 for POAG onset. Multivariable models demonstrated that a one-interquartile range increase in the AT5 (hazard ratio [HR] = 1.14; 95% confidence interval [CI], 1.04-1.25), AT8 (HR = 1.22; 95% CI, 1.09-1.37), AT15 (HR = 1.26; 95% CI, 1.12-1.41), and AT17 (HR = 1.17; 95% CI, 1.03-1.31) coefficients conferred increased risk of POAG onset. AT5, AT10, and AT14 were significantly associated with rapid VF progression. In a subgroup analysis by high-risk ATs (\u3e95th percentile or \u3c75th percentile coefficients), PSD lost significance as a predictor of POAG in the low-risk group. CONCLUSIONS: Baseline VFs, prior to detectable glaucomatous damage, contain occult patterns representing early changes that may increase the risk of POAG onset and VF progression in patients with ocular hypertension. The relationship between PSD and POAG is modified by the presence of high-risk patterns at baseline. An AT-based prediction model for POAG may provide more interpretable glaucoma-specific information in a clinical setting

Identification of Novel Type 2 Diabetes Candidate Genes Involved in the Crosstalk between the Mitochondrial and the Insulin Signaling Systems

Type 2 Diabetes (T2D) is a highly prevalent chronic metabolic disease with strong co-morbidity with obesity and cardiovascular diseases. There is growing evidence supporting the notion that a crosstalk between mitochondria and the insulin signaling cascade could be involved in the etiology of T2D and insulin resistance. In this study we investigated the molecular basis of this crosstalk by using systems biology approaches. We combined, filtered, and interrogated different types of functional interaction data, such as direct protein–protein interactions, co-expression analyses, and metabolic and signaling dependencies. As a result, we constructed the mitochondria-insulin (MITIN) network, which highlights 286 genes as candidate functional linkers between these two systems. The results of internal gene expression analysis of three independent experimental models of mitochondria and insulin signaling perturbations further support the connecting roles of these genes. In addition, we further assessed whether these genes are involved in the etiology of T2D using the genome-wide association study meta-analysis from the DIAGRAM consortium, involving 8,130 T2D cases and 38,987 controls. We found modest enrichment of genes associated with T2D amongst our linker genes (p = 0.0549), including three already validated T2D SNPs and 15 additional SNPs, which, when combined, were collectively associated to increased fasting glucose levels according to MAGIC genome wide meta-analysis (p = 8.12×10−5). This study highlights the potential of combining systems biology, experimental, and genome-wide association data mining for identifying novel genes and related variants that increase vulnerability to complex diseases

Identification of novel type 2 diabetes candidate genes involved in the crosstalk between the mitochondrial and the insulin signaling systems

Type 2 Diabetes (T2D) is a highly prevalent chronic metabolic disease with strong co-morbidity with obesity and cardiovascular diseases. There is growing evidence supporting the notion that a crosstalk between mitochondria and the insulin signaling cascade could be involved in the etiology of T2D and insulin resistance. In this study we investigated the molecular basis of this crosstalk by using systems biology approaches. We combined, filtered, and interrogated different types of functional interaction data, such as direct protein-protein interactions, co-expression analyses, and metabolic and signaling dependencies. As a result, we constructed the mitochondria-insulin (MITIN) network, which highlights 286 genes as candidate functional linkers between these two systems. The results of internal gene expression analysis of three independent experimental models of mitochondria and insulin signaling perturbations further support the connecting roles of these genes. In addition, we further assessed whether these genes are involved in the etiology of T2D using the genome-wide association study meta-analysis from the DIAGRAM consortium, involving 8,130 T2D cases and 38,987 controls. We found modest enrichment of genes associated with T2D amongst our linker genes (p = 0.0549), including three already validated T2D SNPs and 15 additional SNPs, which, when combined, were collectively associated to increased fasting glucose levels according to MAGIC genome wide meta-analysis (p = 8.12×10(-5)). This study highlights the potential of combining systems biology, experimental, and genome-wide association data mining for identifying novel genes and related variants that increase vulnerability to complex diseases

Genetic Analysis in UK Biobank Links Insulin Resistance and Transendothelial Migration Pathways to Coronary Artery Disease

UK Biobank is among the world’s largest repositories for phenotypic and genotypic information in individuals of European ancestry1. We performed a genome-wide association study in UK Biobank testing ~9 million DNA sequence variants for association with coronary artery disease (4,831 cases; 115,455 controls) and carried out meta-analysis with previously published results. We identified fifteen novel loci, bringing the total number of coronary artery disease-associated loci to 95. Phenome-wide association scanning revealed that CCDC92 likely affects coronary artery disease through insulin resistance pathways whereas experimental analysis suggests that ARHGEF26 impacts the transendothelial migration of leukocytes

Integrating genetic regulation and single-cell expression with GWAS prioritizes causal genes and cell types for glaucoma

Primary open-angle glaucoma (POAG), characterized by retinal ganglion cell death, is a leading cause of irreversible blindness worldwide; however, the molecular and cellular causes are not well understood. Elevated intraocular pressure (IOP) is a major risk factor, but many patients have normal IOP. Colocalization and Mendelian randomization analysis of >240 POAG and IOP genome-wide association study (GWAS) loci and of overlapping expression and splicing quantitative trait loci (e/QTLs and sQTLs) in 49 GTEx tissues and retina prioritizesd causal genes for 60% of loci. These genes awere enriched in pathways implicated in extracellular matrix organization, cell adhesion, and vascular development. Analysis of single-nucleus RNA-seq of glaucoma-relevant eye tissues revealesd that the colocalizing genes and genome-wide POAG and IOP associations awere enriched in specific cell types in the aqueous outflow pathways, retina, optic nerve head, peripapillary sclera, and choroid. This study nominatesd IOP-dependent and independent regulatory mechanisms, genes, and cell types that may contribute to POAG pathogenesis

Large-scale multitrait genome-wide association analyses identify hundreds of glaucoma risk loci

Glaucoma, a leading cause of irreversible blindness, is a highly heritable human disease. Previous genome-wide association studies have identified over 100 loci for the most common form, primary open-angle glaucoma. Two key glaucoma-associated traits also show high heritability: intraocular pressure and optic nerve head excavation damage quantified as the vertical cup-to-disc ratio. Here, since much of glaucoma heritability remains unexplained, we conducted a large-scale multitrait genome-wide association study in participants of European ancestry combining primary open-angle glaucoma and its two associated traits (total sample size over 600,000) to substantially improve genetic discovery power (263 loci). We further increased our power by then employing a multiancestry approach, which increased the number of independent risk loci to 312, with the vast majority replicating in a large independent cohort from 23andMe, Inc. (total sample size over 2.8 million; 296 loci replicated at P < 0.05, 240 after Bonferroni correction). Leveraging multiomics datasets, we identified many potential druggable genes, including neuro-protection targets likely to act via the optic nerve, a key advance for glaucoma because all existing drugs only target intraocular pressure. We further used Mendelian randomization and genetic correlation-based approaches to identify novel links to other complex traits, including immune-related diseases such as multiple sclerosis and systemic lupus erythematosus

Integrative Genomics Reveals Novel Molecular Pathways and Gene Networks for Coronary Artery Disease

The majority of the heritability of coronary artery disease (CAD) remains unexplained, despite recent successes of genome-wide association studies (GWAS) in identifying novel susceptibility loci. Integrating functional genomic data from a variety of sources with a large-scale meta-analysis of CAD GWAS may facilitate the identification of novel biological processes and genes involved in CAD, as well as clarify the causal relationships of established processes. Towards this end, we integrated 14 GWAS from the CARDIoGRAM Consortium and two additional GWAS from the Ottawa Heart Institute (25,491 cases and 66,819 controls) with 1) genetics of gene expression studies of CAD-relevant tissues in humans, 2) metabolic and signaling pathways from public databases, and 3) data-driven, tissue-specific gene networks from a multitude of human and mouse experiments. We not only detected CAD-associated gene networks of lipid metabolism, coagulation, immunity, and additional networks with no clear functional annotation, but also revealed key driver genes for each CAD network based on the topology of the gene regulatory networks. In particular, we found a gene network involved in antigen processing to be strongly associated with CAD. The key driver genes of this network included glyoxalase I (GLO1) and peptidylprolyl isomerase I (PPIL1), which we verified as regulatory by siRNA experiments in human aortic endothelial cells. Our results suggest genetic influences on a diverse set of both known and novel biological processes that contribute to CAD risk. The key driver genes for these networks highlight potential novel targets for further mechanistic studies and therapeutic interventions