Cardiometabolic risk loci share downstream cis- and trans-gene regulation across tissues and diseases

Integrating Genes Affecting Coronary Artery Disease in Functional Networks by Multi-OMICs Approach

Coronary artery disease (CAD) and myocardial infarction (MI) remain among the leading causes of mortality worldwide, urgently demanding a better understanding of disease etiology, and more efficient therapeutic strategies. Genetic predisposition as well as the environment and lifestyle are thought to contribute to disease risk. It is likely that non-linear and complex interactions occur between these multiple factors, involving simultaneous pathological changes in diverse cell types, tissues, and organs, at multiple molecular levels. Recent technological advances have exponentially expanded the breadth of available -omics data, from genome, epigenome, transcriptome, proteome, metabolome to even the microbiome. Integration of multiple layers of information across several -omics domains, i.e., the so-called multi-omics approach, currently holds the promise as a path toward precision medicine. Indeed, a more meaningful interpretation of genotype-phenotype relationships and the development of successful therapeutics tailored to individual patients are urgently needed. In this review, we will summarize recent findings and applications of integrative multi-omics in elucidating the etiology of CAD/MI; with a special focus on established disease susceptibility loci sequentially identified in genome-wide association studies (GWAS) over the last 10 years. Moreover, in addition to the autosomal genome, we will also consider the genetic variation in our “second genome”—the mitochondrial genome. Finally, we will summarize the current challenges in the field and point to future research directions required in order to successfully and effectively apply these approaches for precision medicine.

Multiscale Analysis of Independent Alzheimer's Cohorts Finds Disruption of Molecular, Genetic, and Clinical Networks by Human Herpesvirus

Investigators have long suspected that pathogenic microbes might contribute to the onset and progression of Alzheimer's disease (AD) although definitive evidence has not been presented. Whether such findings represent a causal contribution, or reflect opportunistic passengers of neurodegeneration, is also difficult to resolve. We constructed multiscale networks of the late-onset AD-associated virome, integrating genomic, transcriptomic, proteomic, and histopathological data across four brain regions from human post-mortem tissue. We observed increased human herpesvirus 6A (HHV-6A) and human herpesvirus 7 (HHV-7) from subjects with AD compared with controls. These results were replicated in two additional, independent and geographically dispersed cohorts. We observed regulatory relationships linking viral abundance and modulators of APP metabolism, including induction of APBB2, APPBP2, BIN1, BACE1, CLU, PICALM, and PSEN1 by HHV-6A. This study elucidates networks linking molecular, clinical, and neuropathological features with viral activity and is consistent with viral activity constituting a general feature of AD.

CD90 Identifies Adventitial Mesenchymal Progenitor Cells in Adult Human Medium- and Large-Sized Arteries

Mesenchymal stem cells (MSCs) reportedly exist in a vascular niche occupying the outer adventitial layer. However, these cells have not been well characterized in vivo in medium- and large-sized arteries in humans, and their potential pathological role is unknown. To address this, healthy and diseased arterial tissues were obtained as surplus surgical specimens and freshly processed. We identified that CD90 marks a rare adventitial population that co-expresses MSC markers including PDGFRα, CD44, CD73, and CD105. However, unlike CD90, these additional markers were widely expressed by other cells. Human adventitial CD90+ cells fulfilled standard MSC criteria, including plastic adherence, spindle morphology, passage ability, colony formation, and differentiation into adipocytes, osteoblasts, and chondrocytes. Phenotypic and transcriptomic profiling, as well as adoptive transfer experiments, revealed a potential role in vascular disease pathogenesis, with the transcriptomic disease signature of these cells being represented in an aortic regulatory gene network that is operative in atherosclerosis.

•

We identify, in situ and in vivo, adventitial CD90+ MSCs in human arteries

•

Human adventitial CD90+ cells fulfill all criteria for an MSC population

•

Other markers, such as CD44 and PDGFRα, were non-specific for adventitial MSCs

•

The CD90+ MSC transcriptomic signature suggests a major role in vascular disease

Using Gene Expression to Annotate Cardiovascular GWAS Loci

Genetic variants at hundreds of loci associated with cardiovascular phenotypes have been identified by genome wide association studies. Most of these variants are located in intronic or intergenic regions rendering the functional and mechanistic follow up difficult. These non-protein-coding regions harbor regulatory sequences. Thus the study of genetic variants associated with transcription—so called expression quantitative trait loci—has emerged as a promising approach to identify regulatory sequence variants. The genes and pathways they control constitute candidate causal drivers at cardiovascular risk loci. This review provides an overview of the expression quantitative trait loci resources available for cardiovascular genetics research and the most commonly used approaches for candidate gene identification.

Translating GWAS Findings to Novel Therapeutic Targets for Coronary Artery Disease

The success of genome-wide association studies (GWAS) has significantly advanced our understanding of the etiology of coronary artery disease (CAD) and opens new opportunities to reinvigorate the stalling CAD drug development. However, there exists remarkable disconnection between the CAD GWAS findings and commercialized drugs. While this could implicate major untapped translational and therapeutic potentials in CAD GWAS, it also brings forward extensive technical challenges. In this review we summarize the motivation to leverage GWAS for drug discovery, outline the critical bottlenecks in the field, and highlight several promising strategies such as functional genomics and network-based approaches to enhance the translational value of CAD GWAS findings in driving novel therapeutics

Exploring the phenotypic consequences of tissue specific gene expression variation inferred from GWAS summary statistics

Scalable, integrative methods to understand mechanisms that link genetic variants with phenotypes are needed. Here we derive a mathematical expression to compute PrediXcan (a gene mapping approach) results using summary data (S-PrediXcan) and show its accuracy and general robustness to misspecified reference sets. We apply this framework to 44 GTEx tissues and 100+ phenotypes from GWAS and meta-analysis studies, creating a growing public catalog of associations that seeks to capture the effects of gene expression variation on human phenotypes. Replication in an independent cohort is shown. Most of the associations are tissue specific, suggesting context specificity of the trait etiology. Colocalized significant associations in unexpected tissues underscore the need for an agnostic scanning of multiple contexts to improve our ability to detect causal regulatory mechanisms. Monogenic disease genes are enriched among significant associations for related traits, suggesting that smaller alterations of these genes may cause a spectrum of milder phenotypes.

A Strategy for Discovery of Endocrine Interactions with Application to Whole-Body Metabolism

Inter-tissue communication via secreted proteins has been established as a vital mechanism for proper physiologic homeostasis. Here, we report a bioinformatics framework using a mouse reference population, the Hybrid Mouse Diversity Panel (HMDP), which integrates global multi-tissue expression data and publicly available resources to identify and functionally annotate novel circuits of tissue-tissue communication. We validate this method by showing that we can identify known as well as novel endocrine factors responsible for communication between tissues. We further show the utility of this approach by identification and mechanistic characterization of two new endocrine factors. Adipose-derived Lipocalin-5 is shown to enhance skeletal muscle mitochondrial function, and liver-secreted Notum promotes browning of white adipose tissue, also known as "beiging." We demonstrate the general applicability of the method by providing in vivo evidence for three additional novel molecules mediating tissue-tissue interactions.

Advances in Transcriptomics: Investigating Cardiovascular Disease at Unprecedented Resolution

Whole-genome transcriptional profiling has become a standard genomic approach to investigate biological processes. RNA sequencing (RNAseq) in particular has witnessed myriad applications in genetics and various biomedical fields. RNAseq involves a relatively simple experimental protocol of RNA extraction and cDNA library preparation and, because of decreasing next-generation sequencing cost and lower computational burden for data processing, has obtained a central role in the modern biology. The recent application of RNAseq methodology to single-cell transcriptional profiling has enabled the more precise characterization of cell lineage and cell state genetic profiles. The development of bioinformatic and statistical tools has provided for differential gene expression analysis, RNA isoform analysis, haplotype-specific analysis of gene expression (allele-specific expression), and analysis of expression quantitative trait loci. We give an overview of these and recent developments in RNAseq methodology with emphasis on quality control, read mapping, feature counting, differential gene expression, allele-specific expression and expression quantitative trait loci analysis, and fusion transcript detection. We describe utilization of RNAseq as a diagnostic tool in Mendelian diseases, complex phenotypes, and cancer and give an overview of long read RNAseq technology. Furthermore, we discuss in detail the recent revolution in single-cell transcriptomics that is reshaping modern biology.

Multiancestry genome-wide association study of 520,000 subjects identifies 32 loci associated with stroke and stroke subtypes

Stroke has multiple etiologies, but the underlying genes and pathways are largely unknown. We conducted a multiancestry genome-wide-association meta-analysis in 521,612 individuals (67,162 cases and 454,450 controls) and discovered 22 new stroke risk loci, bringing the total to 32. We further found shared genetic variation with related vascular traits, including blood pressure, cardiac traits, and venous thromboembolism, at individual loci (n = 18), and using genetic risk scores and linkage-disequilibrium-score regression. Several loci exhibited distinct association and pleiotropy patterns for etiological stroke subtypes. Eleven new susceptibility loci indicate mechanisms not previously implicated in stroke pathophysiology, with prioritization of risk variants and genes accomplished through bioinformatics analyses using extensive functional datasets. Stroke risk loci were significantly enriched in drug targets for antithrombotic therapy.

The Body-wide Transcriptome Landscape of Disease Models

Virtually all diseases affect multiple organs. However, our knowledge of the body-wide effects remains limited. Here, we report the body-wide transcriptome landscape across 13–23 organs of mouse models of myocardial infarction, diabetes, kidney diseases, cancer, and pre-mature aging. Using such datasets, we find (1) differential gene expression in diverse organs across all models; (2) skin as a disease-sensor organ represented by disease-specific activities of putative gene-expression network; (3) a bone-skin cross talk mediated by a bone-derived hormone, FGF23, in response to dysregulated phosphate homeostasis, a known risk-factor for kidney diseases; (4) candidates for the signature activities of many more putative inter-organ cross talk for diseases; and (5) a cross-species map illustrating organ-to-organ and model-to-disease relationships between human and mouse. These findings demonstrate the usefulness and the potential of such body-wide datasets encompassing mouse models of diverse disease types as a resource in biological and medical sciences. Furthermore, the findings described herein could be exploited for designing disease diagnosis and treatment.

•

Body-wide multi-organ transcriptome datasets encompassing diverse disease models

•

Skin is a disease-sensor organ, and FGF23 mediates a bone-skin cross talk in diseases

•

Diverse putative inter-organ cross talk selectively associates with diseases

•

A cross-species map illustrating the mouse-human relationships

Genetic alterations in the NO-cGMP pathway and cardiovascular risk

In the past ten years, several chromosomal loci have been identified by genome-wide association studies to influence the risk of coronary artery disease (CAD) and its risk factors. The GUCY1A3 gene encoding the α₁ subunit of the soluble guanylyl cyclase (sGC) resides at one of these loci and has been strongly associated with blood pressure and CAD risk. More recently, further genes in the pathway encoding the endothelial nitric oxide synthase, the phosphodiesterases 3A and 5A, and the inositol 1,4,5-trisphosphate receptor I-associated protein (IRAG), i.e., NOS3, PDE3A, PDE5A, and MRVI1, respectively, were likewise identified as CAD risk genes. In this review, we highlight the genetic findings linking variants in NO-cGMP signaling and cardiovascular disease, discuss the potential underlying mechanisms which might propagate the development of atherosclerosis, and speculate about therapeutic implications.

Genetic Drivers of Pancreatic Islet Function

The majority of gene loci that have been associated with type 2 diabetes play a role in pancreatic islet function. To evaluate the role of islet gene expression in the etiology of diabetes, we sensitized a genetically diverse mouse population with a Western diet high in fat (45% kcal) and sucrose (34%) and carried out genome-wide association mapping of diabetes-related phenotypes. We quantified mRNA abundance in the islets and identified 18,820 expression QTL. We applied mediation analysis to identify candidate causal driver genes at loci that affect the abundance of numerous transcripts. These include two genes previously associated with monogenic diabetes (PDX1 and HNF4A), as well as three genes with nominal association with diabetes-related traits in humans (FAM83E, IL6ST, and SAT2). We grouped transcripts into gene modules and mapped regulatory loci for modules enriched with transcripts specific for α-cells, and another specific for δ-cells. However, no single module enriched for β-cell-specific transcripts, suggesting heterogeneity of gene expression patterns within the β-cell population. A module enriched in transcripts associated with branched-chain amino acid metabolism was the most strongly correlated with physiological traits that reflect insulin resistance. Although the mice in this study were not overtly diabetic, the analysis of pancreatic islet gene expression under dietary-induced stress enabled us to identify correlated variation in groups of genes that are functionally linked to diabetes-associated physiological traits. Our analysis suggests an expected degree of concordance between diabetes-associated loci in the mouse and those found in human populations, and demonstrates how the mouse can provide evidence to support nominal associations found in human genome-wide association mapping.

Global analysis of A-to-I RNA editing reveals association with common disease variants

RNA editing modifies transcripts and may alter their regulation or function. In humans, the most common modification is adenosine to inosine (A-to-I). We examined the global characteristics of RNA editing in 4,301 human tissue samples. More than 1.6 million A-to-I edits were identified in 62% of all protein-coding transcripts. mRNA recoding was extremely rare; only 11 novel recoding sites were uncovered. Thirty single nucleotide polymorphisms from genome-wide association studies were associated with RNA editing; one that influences type 2 diabetes (rs2028299) was associated with editing in ARPIN. Twenty-five genes, including LRP11 and PLIN5, had editing sites that were associated with plasma lipid levels. Our findings provide new insights into the genetic regulation of RNA editing and establish a rich catalogue for further exploration of this process.

Genetic study links components of the autonomous nervous system to heart-rate profile during exercise

Heart rate (HR) responds to exercise by increasing during exercise and recovering after exercise. As such, HR is an important predictor of mortality that researchers believe is modulated by the autonomic nervous system. However, the mechanistic basis underlying inter-individual differences has yet to be explained. Here, we perform a large-scale genome-wide analysis of HR increase and HR recovery in 58,818 UK Biobank individuals. Twenty-five independent SNPs in 23 loci are identified to be associated (p < 8.3 × 10⁻⁹) with HR increase or HR recovery. A total of 36 candidate causal genes are prioritized that are enriched for pathways related to neuron biology. No evidence is found of a causal relationship with mortality or cardiovascular diseases. However, a nominal association with parental lifespan requires further study. In conclusion, the findings provide new biological and clinical insight into the mechanistic underpinnings of HR response to exercise. The results also underscore the role of the autonomous nervous system in HR recovery.

Response of the heart rate (HR) to exercise is associated with cardiac fitness and risk of cardiac death. Here, in a genome-wide association study, Verweij et al. identify 23 loci for HR increase during exercise or HR recovery, and highlight pleiotropy with blood pressure by polygenic risk score analysis.

Network analysis of coronary artery disease risk genes elucidates disease mechanisms and druggable targets

Genome-wide association studies (GWAS) have identified over two hundred chromosomal loci that modulate risk of coronary artery disease (CAD). The genes affected by variants at these loci are largely unknown and an untapped resource to improve our understanding of CAD pathophysiology and identify potential therapeutic targets. Here, we prioritized 68 genes as the most likely causal genes at genome-wide significant loci identified by GWAS of CAD and examined their regulatory roles in 286 metabolic and vascular tissue gene-protein sub-networks (“modules”). The modules and genes within were scored for CAD druggability potential. The scoring enriched for targets of cardiometabolic drugs currently in clinical use and in-depth analysis of the top-scoring modules validated established and revealed novel target tissues, biological processes, and druggable targets. This study provides an unprecedented resource of tissue-defined gene–protein interactions directly affected by genetic variance in CAD risk loci.

HDL Cholesterol Metabolism and the Risk of CHD: New Insights from Human Genetics

PURPOSE OF REVIEW

Elevated high-density lipoprotein cholesterol levels in the blood (HDL-C) represent one of the strongest epidemiological surrogates for protection against coronary heart disease (CHD), but recent human genetic and pharmacological intervention studies have raised controversy about the causality of this relationship. Here, we review recent discoveries from human genome studies using new analytic tools as well as relevant animal studies that have both addressed, and in some cases, fueled this controversy.

RECENT FINDINGS

Methodologic developments in genotyping and sequencing, such as genome-wide association studies (GWAS), exome sequencing, and exome array genotyping, have been applied to the study of HDL-C and risk of CHD in large, multi-ethnic populations. Some of these efforts focused on population-wide variation in common variants have uncovered new polymorphisms at novel loci associated with HDL-C and, in some cases, CHD risk. Other efforts have discovered loss-of-function variants for the first time in genes previously implicated in HDL metabolism through common variant studies or animal models. These studies have allowed the genetic relationship between these pathways, HDL-C and CHD to be explored in humans for the first time through analysis tools such as Mendelian randomization. We explore these discoveries for selected key HDL-C genes CETP, LCAT, LIPG, SCARB1, and novel loci implicated from GWAS including GALNT2, KLF14, and TTC39B. Recent human genetics findings have identified new nodes regulating HDL metabolism while reshaping our current understanding of known candidate genes to HDL and CHD risk through the study of critical variants across model systems. Despite their effect on HDL-C, variants in many of the reviewed genes were found to lack any association with CHD. These data collectively indicate that HDL-C concentration, which represents a static picture of a very dynamic and heterogeneous metabolic milieu, is unlikely to be itself causally protective against CHD. In this context, human genetics represent an extremely valuable tool to further explore the biological mechanisms regulating HDL metabolism and investigate what role, if any, HDL plays in the pathogenesis of CHD.

Genetic regulation of adipose tissue transcript expression is involved in modulating serum triglyceride and HDL-cholesterol

Dyslipidemia is a major contributor to the increased cardiovascular disease and mortality associated with obesity and type 2 diabetes. We hypothesized that variation in expression of adipose tissue transcripts is associated with serum lipid concentrations in African Americans (AAs), and common genetic variants regulate expression levels of these transcripts. Fasting serum lipid levels, genome-wide transcript expression profiles of subcutaneous adipose tissue, and genome-wide SNP genotypes were analyzed in a cohort of non-diabetic AAs (N=250). Serum triglyceride (TRIG) and high density lipoprotein-cholesterol (HDL-C) levels were associated (FDR<0.01) with expression level of 1021 and 1875 adipose tissue transcripts, respectively, but none associated with total cholesterol or LDL-C levels. Serum HDL-C-associated transcripts were enriched for salient biological pathways, including branched-chain amino acid degradation, and oxidative phosphorylation. Genes in immuno-inflammatory pathways were activated among individuals with higher serum TRIG levels. We identified significant cis-regulatory SNPs (cis-eSNPs) for 449 serum lipid-associated transcripts in adipose tissue. The cis-eSNPs of 12 genes were nominally associated (p<0.001) with serum lipid level in genome wide association studies in Global Lipids Genetics Consortium (GLGC) cohorts. Allelic effect direction of cis-eSNPs on expression of MARCH2, BEST1 and TMEM258 matched with effect direction of these SNP alleles on serum TRIG or HDL-C levels in GLGC cohorts. These data suggest that expressions of serum lipid-associated transcripts in adipose tissue are dependent on common cis-eSNPs in African Americans. Thus, genetically-mediated transcriptional regulation in adipose tissue may play a role in reducing HDL-C and increasing TRIG in serum.

Network-directed cis-mediator analysis of normal prostate tissue expression profiles reveals downstream regulatory associations of prostate cancer susceptibility loci

Large-scale genome-wide association studies have identified multiple single-nucleotide polymorphisms associated with risk of prostate cancer. Many of these genetic variants are presumed to be regulatory in nature; however, follow-up expression quantitative trait loci (eQTL) association studies have to-date been restricted largely to cis-acting associations due to study limitations. While trans-eQTL scans suffer from high testing dimensionality, recent evidence indicates most trans-eQTL associations are mediated by cis-regulated genes, such as transcription factors. Leveraging a data-driven gene co-expression network, we conducted a comprehensive cis-mediator analysis using RNA-Seq data from 471 normal prostate tissue samples to identify downstream regulatory associations of previously identified prostate cancer risk variants. We discovered multiple trans-eQTL associations that were significantly mediated by cis-regulated transcripts, four of which involved risk locus 17q12, proximal transcription factor HNF1B, and target trans-genes with known HNF response elements (MIA2, SRC, SEMA6A, KIF12). We additionally identified evidence of cis-acting down-regulation of MSMB via rs10993994 corresponding to reduced co-expression of NDRG1. The majority of these cis-mediator relationships demonstrated trans-eQTL replicability in 87 prostate tissue samples from the Gene-Tissue Expression Project. These findings provide further biological context to known risk loci and outline new hypotheses for investigation into the etiology of prostate cancer.

Physical Interactions and Expression Quantitative Traits Loci Identify Regulatory Connections for Obesity and Type 2 Diabetes Associated SNPs

The mechanisms that underlie the association between obesity and type 2 diabetes are not fully understood. Here, we investigated the role of the 3D genome organization in the pathogeneses of obesity and type-2 diabetes. We interpreted the combined and differential impacts of 196 diabetes and 390 obesity associated single nucleotide polymorphisms (SNPs) by integrating data on the genes with which they physically interact (as captured by Hi-C) and the functional [i.e., expression quantitative trait loci (eQTL)] outcomes associated with these interactions. We identified 861 spatially regulated genes (e.g., AP3S2, ELP5, SVIP, IRS1, FADS2, WFS1, RBM6, HORMAD1, PYROXD2), which are enriched in tissues (e.g., adipose, skeletal muscle, pancreas) and biological processes and canonical pathways (e.g., lipid metabolism, leptin, and glucose-insulin signaling pathways) that are important for the pathogenesis of type 2 diabetes and obesity. Our discovery-based approach also identifies enrichment for eQTL SNP-gene interactions in tissues that are not classically associated with diabetes or obesity. We propose that the combinatorial action of active obesity and diabetes spatial eQTL SNPs on their gene pairs within different tissues reduces the ability of these tissues to contribute to the maintenance of a healthy energy metabolism.

Open chromatin profiling of human postmortem brain infers functional roles for non-coding schizophrenia loci

Open chromatin provides access to DNA-binding proteins for the correct spatiotemporal regulation of gene expression. Mapping chromatin accessibility has been widely used to identify the location of cis regulatory elements (CREs) including promoters and enhancers. CREs show tissue- and cell-type specificity and disease-associated variants are often enriched for CREs in the tissues and cells that pertain to a given disease. To better understand the role of CREs in neuropsychiatric disorders we applied the Assay for Transposase Accessible Chromatin followed by sequencing (ATAC-seq) to neuronal and non-neuronal nuclei isolated from frozen postmortem human brain by fluorescence-activated nuclear sorting (FANS). Most of the identified open chromatin regions (OCRs) are differentially accessible between neurons and non-neurons, and show enrichment with known cell type markers, promoters and enhancers. Relative to those of non-neurons, neuronal OCRs are more evolutionarily conserved and are enriched in distal regulatory elements. Transcription factor (TF) footprinting analysis identifies differences in the regulome between neuronal and non-neuronal cells and ascribes putative functional roles to a number of non-coding schizophrenia (SCZ) risk variants. Among the identified variants is a Single Nucleotide Polymorphism (SNP) proximal to the gene encoding SNX19. In vitro experiments reveal that this SNP leads to an increase in transcriptional activity. As elevated expression of SNX19 has been associated with SCZ, our data provide evidence that the identified SNP contributes to disease. These results represent the first analysis of OCRs and TF-binding sites in distinct populations of postmortem human brain cells and further our understanding of the regulome and the impact of neuropsychiatric disease-associated genetic risk variants.

New Blood Pressure-Associated Loci Identified in Meta-Analyses of 475 000 Individuals

BACKGROUND

Genome-wide association studies have recently identified >400 loci that harbor DNA sequence variants that influence blood pressure (BP). Our earlier studies identified and validated 56 single nucleotide variants (SNVs) associated with BP from meta-analyses of exome chip genotype data. An additional 100 variants yielded suggestive evidence of association.

METHODS AND RESULTS

Here, we augment the sample with 140 886 European individuals from the UK Biobank, in whom 77 of the 100 suggestive SNVs were available for association analysis with systolic BP or diastolic BP or pulse pressure. We performed 2 meta-analyses, one in individuals of European, South Asian, African, and Hispanic descent (pan-ancestry, ≈475 000), and the other in the subset of individuals of European descent (≈423 000). Twenty-one SNVs were genome-wide significant (P<5×10-8) for BP, of which 4 are new BP loci: rs9678851 (missense, SLC4A1AP), rs7437940 (AFAP1), rs13303 (missense, STAB1), and rs1055144 (7p15.2). In addition, we identified a potentially independent novel BP-associated SNV, rs3416322 (missense, SYNPO2L) at a known locus, uncorrelated with the previously reported SNVs. Two SNVs are associated with expression levels of nearby genes, and SNVs at 3 loci are associated with other traits. One SNV with a minor allele frequency <0.01, (rs3025380 at DBH) was genome-wide significant.

CONCLUSIONS

We report 4 novel loci associated with BP regulation, and 1 independent variant at an established BP locus. This analysis highlights several candidate genes with variation that alter protein function or gene expression for potential follow-up.

Enhancer connectome in primary human cells identifies target genes of disease-associated DNA elements

The challenge of linking intergenic mutations to target genes has limited molecular understanding of human diseases. Here we show that H3K27ac HiChIP generates high-resolution contact maps of active enhancers and target genes in rare primary human T cell subtypes and coronary artery smooth muscle cells. Differentiation of naive T cells into T helper 17 cells or regulatory T cells creates subtype-specific enhancer-promoter interactions, specifically at regions of shared DNA accessibility. These data provide a principled means of assigning molecular functions to autoimmune and cardiovascular disease risk variants, linking hundreds of noncoding variants to putative gene targets. Target genes identified with HiChIP are further supported by CRISPR interference and activation at linked enhancers, by the presence of expression quantitative trait loci, and by allele-specific enhancer loops in patient-derived primary cells. The majority of disease-associated enhancers contact genes beyond the nearest gene in the linear genome, leading to a fourfold increase in the number of potential target genes for autoimmune and cardiovascular diseases.

Exploring regulation in tissues with eQTL networks

Characterizing the collective regulatory impact of genetic variants on complex phenotypes is a major challenge in developing a genotype to phenotype map. Using expression quantitative trait locus (eQTL) analyses, we constructed bipartite networks in which edges represent significant associations between genetic variants and gene expression levels and found that the network structure informs regulatory function. We show, in 13 tissues, that these eQTL networks are organized into dense, highly modular communities grouping genes often involved in coherent biological processes. We find communities representing shared processes across tissues, as well as communities associated with tissue-specific processes that coalesce around variants in tissue-specific active chromatin regions. Node centrality is also highly informative, with the global and community hubs differing in regulatory potential and likelihood of being disease associated.

Carbonyl reductase 1 catalyzes 20β-reduction of glucocorticoids, modulating receptor activation and metabolic complications of obesity

Carbonyl Reductase 1 (CBR1) is a ubiquitously expressed cytosolic enzyme important in exogenous drug metabolism but the physiological function of which is unknown. Here, we describe a role for CBR1 in metabolism of glucocorticoids. CBR1 catalyzes the NADPH- dependent production of 20β-dihydrocortisol (20β-DHF) from cortisol. CBR1 provides the major route of cortisol metabolism in horses and is up-regulated in adipose tissue in obesity in horses, humans and mice. We demonstrate that 20β-DHF is a weak endogenous agonist of the human glucocorticoid receptor (GR). Pharmacological inhibition of CBR1 in diet-induced obesity in mice results in more marked glucose intolerance with evidence for enhanced hepatic GR signaling. These findings suggest that CBR1 generating 20β-dihydrocortisol is a novel pathway modulating GR activation and providing enzymatic protection against excessive GR activation in obesity.

Expression quantitative trait loci (eQTLs) in human placentas suggest developmental origins of complex diseases

Epidemiologic studies support that at least part of the risk of chronic diseases in childhood and even adulthood may have an in utero origin, and the placenta is a key organ that plays a pivotal role in fetal growth and development. The transcriptomes of 159 human placenta tissues were profiled by genome-wide RNA sequencing (Illumina High-Seq 2500), and linked to fetal genotypes assessed by a high density single nucleotide polymorphism (SNP) genotyping array (Illumina MegaEx). Expression quantitative trait loci (eQTLs) across all annotated transcripts were mapped and examined for enrichment for disease susceptibility loci annotated in the genome-wide association studies (GWAS) catalog. We discovered 3218 cis- and 35 trans-eQTLs at ≤10% false discovery rate in human placentas. Among the 16 439 known disease loci of genome-wide significance, 835 were placental eSNPs (enrichment fold = 1.68, P = 7.41e-42). Stronger effect sizes were observed between GWAS SNPs and gene expression in placentas than what has been reported in other tissues, such as the correlation between asthma risk allele, rs7216389-T and Gasdermin-B (GSDMB) in placenta (r2=27%) versus lung (r2=6%). Finally, our results suggest the placental eQTLs may mediate the function of GWAS loci on postnatal disease susceptibility. Results suggest that transcripts in placenta are under tight genetic control, and that placental gene networks may influence postnatal risk of multiple human diseases lending support for the Developmental Origins of Health and Disease.

Shared genetic regulatory networks for cardiovascular disease and type 2 diabetes in multiple populations of diverse ethnicities in the United States

Cardiovascular diseases (CVD) and type 2 diabetes (T2D) are closely interrelated complex diseases likely sharing overlapping pathogenesis driven by aberrant activities in gene networks. However, the molecular circuitries underlying the pathogenic commonalities remain poorly understood. We sought to identify the shared gene networks and their key intervening drivers for both CVD and T2D by conducting a comprehensive integrative analysis driven by five multi-ethnic genome-wide association studies (GWAS) for CVD and T2D, expression quantitative trait loci (eQTLs), ENCODE, and tissue-specific gene network models (both co-expression and graphical models) from CVD and T2D relevant tissues. We identified pathways regulating the metabolism of lipids, glucose, and branched-chain amino acids, along with those governing oxidation, extracellular matrix, immune response, and neuronal system as shared pathogenic processes for both diseases. Further, we uncovered 15 key drivers including HMGCR, CAV1, IGF1 and PCOLCE, whose network neighbors collectively account for approximately 35% of known GWAS hits for CVD and 22% for T2D. Finally, we cross-validated the regulatory role of the top key drivers using in vitro siRNA knockdown, in vivo gene knockout, and two Hybrid Mouse Diversity Panels each comprised of >100 strains. Findings from this in-depth assessment of genetic and functional data from multiple human cohorts provide strong support that common sets of tissue-specific molecular networks drive the pathogenesis of both CVD and T2D across ethnicities and help prioritize new therapeutic avenues for both CVD and T2D.

Association analyses based on false discovery rate implicate new loci for coronary artery disease

Genome-wide association studies (GWAS) in coronary artery disease (CAD) had identified 66 loci at 'genome-wide significance' (P < 5 × 10-8) at the time of this analysis, but a much larger number of putative loci at a false discovery rate (FDR) of 5% (refs. 1,2,3,4). Here we leverage an interim release of UK Biobank (UKBB) data to evaluate the validity of the FDR approach. We tested a CAD phenotype inclusive of angina (SOFT; ncases = 10,801) as well as a stricter definition without angina (HARD; ncases = 6,482) and selected cases with the former phenotype to conduct a meta-analysis using the two most recent CAD GWAS. This approach identified 13 new loci at genome-wide significance, 12 of which were on our previous list of loci meeting the 5% FDR threshold, thus providing strong support that the remaining loci identified by FDR represent genuine signals. The 304 independent variants associated at 5% FDR in this study explain 21.2% of CAD heritability and identify 243 loci that implicate pathways in blood vessel morphogenesis as well as lipid metabolism, nitric oxide signaling and inflammation.

Efficient and accurate causal inference with hidden confounders from genome-transcriptome variation data

Mapping gene expression as a quantitative trait using whole genome-sequencing and transcriptome analysis allows to discover the functional consequences of genetic variation. We developed a novel method and ultra-fast software Findr for higly accurate causal inference between gene expression traits using cis-regulatory DNA variations as causal anchors, which improves current methods by taking into consideration hidden confounders and weak regulations. Findr outperformed existing methods on the DREAM5 Systems Genetics challenge and on the prediction of microRNA and transcription factor targets in human lymphoblastoid cells, while being nearly a million times faster. Findr is publicly available at https://github.com/lingfeiwang/findr.

Understanding how genetic variation between individuals determines variation in observable traits or disease risk is one of the core aims of genetics. It is known that genetic variation often affects gene regulatory DNA elements and directly causes variation in expression of nearby genes. This effect in turn cascades down to other genes via the complex pathways and gene interaction networks that ultimately govern how cells operate in an ever changing environment. In theory, when genetic variation and gene expression levels are measured simultaneously in a large number of individuals, the causal effects of genes on each other can be inferred using statistical models similar to those used in randomized controlled trials. We developed a novel method and ultra-fast software Findr which, unlike existing methods, takes into account the complex but unknown network context when predicting causality between specific gene pairs. Findr’s predictions have a significantly higher overlap with known gene networks compared to existing methods, using both simulated and real data. Findr is also nearly a million times faster, and hence the only software in its class that can handle modern datasets where the expression levels of ten-thousands of genes are simultaneously measured in hundreds to thousands of individuals.

The Genetic Architecture of Coronary Artery Disease: Current Knowledge and Future Opportunities

PURPOSE OF REVIEW

We provide an overview of our current understanding of the genetic architecture of coronary artery disease (CAD) and discuss areas of research that provide excellent opportunities for further exploration.

RECENT FINDINGS

Large-scale studies in human populations, coupled with rapid advances in genetic technologies over the last decade, have clearly established the association of common genetic variation with risk of CAD. However, the effect sizes of the susceptibility alleles are for the most part modest and collectively explain only a small fraction of the overall heritability. By comparison, evidence that rare variants make a substantial contribution to risk of CAD has been somewhat disappointing thus far, suggesting that other biological mechanisms have yet to be discovered. Emerging data suggests that novel pathways involved in the development of CAD can be identified through complementary and integrative systems genetics strategies in mice or humans. There is also convincing evidence that gut bacteria play a previously unrecognized role in the development of CAD, particularly through metabolism of certain dietary nutrients that lead to proatherogenic metabolites in the circulation. A major effort is now underway to functionally understand the newly discovered genetic and biological associations for CAD, which could lead to the development of potentially novel therapeutic strategies. Other important areas of investigation for understanding the pathophysiology of CAD, including epistatic interactions between genes or with either sex and environmental factors, have not been studied on a broad scope and represent additional opportunities for future studies.

A data-driven modeling approach to identify disease-specific multi-organ networks driving physiological dysregulation

Multiple physiological systems interact throughout the development of a complex disease. Knowledge of the dynamics and connectivity of interactions across physiological systems could facilitate the prevention or mitigation of organ damage underlying complex diseases, many of which are currently refractory to available therapeutics (e.g., hypertension). We studied the regulatory interactions operating within and across organs throughout disease development by integrating in vivo analysis of gene expression dynamics with a reverse engineering approach to infer data-driven dynamic network models of multi-organ gene regulatory influences. We obtained experimental data on the expression of 22 genes across five organs, over a time span that encompassed the development of autonomic nervous system dysfunction and hypertension. We pursued a unique approach for identification of continuous-time models that jointly described the dynamics and structure of multi-organ networks by estimating a sparse subset of ∼12,000 possible gene regulatory interactions. Our analyses revealed that an autonomic dysfunction-specific multi-organ sequence of gene expression activation patterns was associated with a distinct gene regulatory network. We analyzed the model structures for adaptation motifs, and identified disease-specific network motifs involving genes that exhibited aberrant temporal dynamics. Bioinformatic analyses identified disease-specific single nucleotide variants within or near transcription factor binding sites upstream of key genes implicated in maintaining physiological homeostasis. Our approach illustrates a novel framework for investigating the pathogenesis through model-based analysis of multi-organ system dynamics and network properties. Our results yielded novel candidate molecular targets driving the development of cardiovascular disease, metabolic syndrome, and immune dysfunction.

Complex diseases such as hypertension often involve maladaptive autonomic nervous system control over the cardiovascular, renal, hepatic, immune, and endocrine systems. We studied the pathogenesis of physiological homeostasis by examining the temporal dynamics of gene expression levels from multiple organs in an animal model of autonomic dysfunction characterized by cardiovascular disease, metabolic dysregulation, and immune system aberrations. We employed a data-driven modeling approach to jointly predict continuous gene expression dynamics and gene regulatory interactions across organs in the disease and control phenotypes. We combined our analyses of multi-organ gene regulatory network dynamics and connectivity with bioinformatic analyses of genetic mutations that could regulate gene expression. Our multi-organ modeling approach to investigate the mechanisms of complex disease pathogenesis revealed novel candidates for therapeutic interventions against the development and progression of complex diseases involving autonomic nervous system dysfunction.

Fifteen new risk loci for coronary artery disease highlight arterial-wall-specific mechanisms

Coronary artery disease (CAD) is a leading cause of morbidity and mortality worldwide. Although 58 genomic regions have been associated with CAD thus far, most of the heritability is unexplained, indicating that additional susceptibility loci await identification. An efficient discovery strategy may be larger-scale evaluation of promising associations suggested by genome-wide association studies (GWAS). Hence, we genotyped 56,309 participants using a targeted gene array derived from earlier GWAS results and performed meta-analysis of results with 194,427 participants previously genotyped, totaling 88,192 CAD cases and 162,544 controls. We identified 25 new SNP-CAD associations (P < 5 × 10-8, in fixed-effects meta-analysis) from 15 genomic regions, including SNPs in or near genes involved in cellular adhesion, leukocyte migration and atherosclerosis (PECAM1, rs1867624), coagulation and inflammation (PROCR, rs867186 (p.Ser219Gly)) and vascular smooth muscle cell differentiation (LMOD1, rs2820315). Correlation of these regions with cell-type-specific gene expression and plasma protein levels sheds light on potential disease mechanisms.

Enabling Precision Cardiology Through Multiscale Biology and Systems Medicine

The traditional paradigm of cardiovascular disease research derives insight from large-scale, broadly inclusive clinical studies of well-characterized pathologies. These insights are then put into practice according to standardized clinical guidelines. However, stagnation in the development of new cardiovascular therapies and variability in therapeutic response implies that this paradigm is insufficient for reducing the cardiovascular disease burden. In this state-of-the-art review, we examine 3 interconnected ideas we put forth as key concepts for enabling a transition to precision cardiology: 1) precision characterization of cardiovascular disease with machine learning methods; 2) the application of network models of disease to embrace disease complexity; and 3) using insights from the previous 2 ideas to enable pharmacology and polypharmacology systems for more precise drug-to-patient matching and patient-disease stratification. We conclude by exploring the challenges of applying a precision approach to cardiology, which arise from a deficit of the required resources and infrastructure, and emerging evidence for the clinical effectiveness of this nascent approach.

Identification of 15 novel risk loci for coronary artery disease and genetic risk of recurrent events, atrial fibrillation and heart failure

Coronary artery disease (CAD) is the major cause of morbidity and mortality in the world. Identification of novel genetic determinants may provide new opportunities for developing innovative strategies to predict, prevent and treat CAD. Therefore, we meta-analyzed independent genetic variants passing P <× 10^-5 in CARDIoGRAMplusC4D with novel data made available by UK Biobank. Of the 161 genetic variants studied, 71 reached genome wide significance (p < 5 × 10^-8) including 15 novel loci. These novel loci include multiple genes that are involved in angiogenesis (TGFB1, ITGB5, CDH13 and RHOA) and 2 independent variants in the TGFB1 locus. We also identified SGEF as a candidate gene in one of the novel CAD loci. SGEF was previously suggested as a therapeutic target based on mouse studies. The genetic risk score of CAD predicted recurrent CAD events and cardiovascular mortality. We also identified significant genetic correlations between CAD and other cardiovascular conditions, including heart failure and atrial fibrillation. In conclusion, we substantially increased the number of loci convincingly associated with CAD and provide additional biological and clinical insights.

Large-Scale Identification of Common Trait and Disease Variants Affecting Gene Expression

Genome-wide association studies (GWASs) have identified a multitude of genetic loci involved with traits and diseases. However, it is often unclear which genes are affected in such loci and whether the associated genetic variants lead to increased or decreased gene function. To mitigate this, we integrated associations of common genetic variants in 57 GWASs with 24 studies of expression quantitative trait loci (eQTLs) from a broad range of tissues by using a Mendelian randomization approach. We discovered a total of 3,484 instances of gene-trait-associated changes in expression at a false-discovery rate < 0.05. These genes were often not closest to the genetic variant and were primarily identified in eQTLs derived from pathophysiologically relevant tissues. For instance, genes with expression changes associated with lipid traits were mostly identified in the liver, and those associated with cardiovascular disease were identified in arterial tissue. The affected genes additionally point to biological processes implicated in the interrogated traits, such as the interleukin-27 pathway in rheumatoid arthritis. Further, comparing trait-associated gene expression changes across traits suggests that pleiotropy is a widespread phenomenon and points to specific instances of both agonistic and antagonistic pleiotropy. For instance, expression of SNX19 and ABCB9 is positively correlated with both the risk of schizophrenia and educational attainment. To facilitate interpretation, we provide this lexicon of how common trait-associated genetic variants alter gene expression in various tissues as the online database GWAS2Genes.

Genetic regulation of gene expression in the epileptic human hippocampus

Epilepsy is a serious and common neurological disorder. Expression quantitative loci (eQTL) analysis is a vital aid for the identification and interpretation of disease-risk loci. Many eQTLs operate in a tissue- and condition-specific manner. We have performed the first genome-wide cis-eQTL analysis of human hippocampal tissue to include not only normal (n = 22) but also epileptic (n = 22) samples. We demonstrate that disease-associated variants from an epilepsy GWAS meta-analysis and a febrile seizures (FS) GWAS are significantly more enriched with epilepsy-eQTLs than with normal hippocampal eQTLs from two larger independent published studies. In contrast, GWAS meta-analyses of two other brain diseases associated with hippocampal pathology (Alzheimer's disease and schizophrenia) are more enriched with normal hippocampal eQTLs than with epilepsy-eQTLs. These observations suggest that an eQTL analysis that includes disease-affected brain tissue is advantageous for detecting additional risk SNPs for the afflicting and closely related disorders, but not for distinct diseases affecting the same brain regions. We also show that epilepsy eQTLs are enriched within epilepsy-causing genes: an epilepsy cis-gene is significantly more likely to be a causal gene for a Mendelian epilepsy syndrome than to be a causal gene for another Mendelian disorder. Epilepsy cis-genes, compared to normal hippocampal cis-genes, are more enriched within epilepsy-causing genes. Hence, we utilize the epilepsy eQTL data for the functional interpretation of epilepsy disease-risk variants and, thereby, highlight novel potential causal genes for sporadic epilepsy. In conclusion, an epilepsy-eQTL analysis is superior to normal hippocampal tissue eQTL analyses for identifying the variants and genes underlying epilepsy.

Multidimensional Integrative Genomics Approaches to Dissecting Cardiovascular Disease

Elucidating the mechanisms of complex diseases such as cardiovascular disease (CVD) remains a significant challenge due to multidimensional alterations at molecular, cellular, tissue, and organ levels. To better understand CVD and offer insights into the underlying mechanisms and potential therapeutic strategies, data from multiple omics types (genomics, epigenomics, transcriptomics, metabolomics, proteomics, microbiomics) from both humans and model organisms have become available. However, individual omics data types capture only a fraction of the molecular mechanisms. To address this challenge, there have been numerous efforts to develop integrative genomics methods that can leverage multidimensional information from diverse data types to derive comprehensive molecular insights. In this review, we summarize recent methodological advances in multidimensional omics integration, exemplify their applications in cardiovascular research, and pinpoint challenges and future directions in this incipient field.

Decreased STARD10 Expression Is Associated with Defective Insulin Secretion in Humans and Mice

Genetic variants near ARAP1 (CENTD2) and STARD10 influence type 2 diabetes (T2D) risk. The risk alleles impair glucose-induced insulin secretion and, paradoxically but characteristically, are associated with decreased proinsulin:insulin ratios, indicating improved proinsulin conversion. Neither the identity of the causal variants nor the gene(s) through which risk is conferred have been firmly established. Whereas ARAP1 encodes a GTPase activating protein, STARD10 is a member of the steroidogenic acute regulatory protein (StAR)-related lipid transfer protein family. By integrating genetic fine-mapping and epigenomic annotation data and performing promoter-reporter and chromatin conformational capture (3C) studies in β cell lines, we localize the causal variant(s) at this locus to a 5 kb region that overlaps a stretch-enhancer active in islets. This region contains several highly correlated T2D-risk variants, including the rs140130268 indel. Expression QTL analysis of islet transcriptomes from three independent subject groups demonstrated that T2D-risk allele carriers displayed reduced levels of STARD10 mRNA, with no concomitant change in ARAP1 mRNA levels. Correspondingly, β-cell-selective deletion of StarD10 in mice led to impaired glucose-stimulated Ca²⁺ dynamics and insulin secretion and recapitulated the pattern of improved proinsulin processing observed at the human GWAS signal. Conversely, overexpression of StarD10 in the adult β cell improved glucose tolerance in high fat-fed animals. In contrast, manipulation of Arap1 in β cells had no impact on insulin secretion or proinsulin conversion in mice. This convergence of human and murine data provides compelling evidence that the T2D risk associated with variation at this locus is mediated through reduction in STARD10 expression in the β cell.

EXPLORING THE REPRODUCIBILITY OF PROBABILISTIC CAUSAL MOLECULAR NETWORK MODELS

Network reconstruction algorithms are increasingly being employed in biomedical and life sciences research to integrate large-scale, high-dimensional data informing on living systems. One particular class of probabilistic causal networks being applied to model the complexity and causal structure of biological data is Bayesian networks (BNs). BNs provide an elegant mathematical framework for not only inferring causal relationships among many different molecular and higher order phenotypes, but also for incorporating highly diverse priors that provide an efficient path for incorporating existing knowledge. While significant methodological developments have broadly enabled the application of BNs to generate and validate meaningful biological hypotheses, the reproducibility of BNs in this context has not been systematically explored. In this study, we aim to determine the criteria for generating reproducible BNs in the context of transcription-based regulatory networks. We utilize two unique tissues from independent datasets, whole blood from the GTEx Consortium and liver from the Stockholm-Tartu Atherosclerosis Reverse Network Engineering Team (STARNET) study. We evaluated the reproducibility of the BNs by creating networks on data subsampled at different levels from each cohort and comparing these networks to the BNs constructed using the complete data. To help validate our results, we used simulated networks at varying sample sizes. Our study indicates that reproducibility of BNs in biological research is an issue worthy of further consideration, especially in light of the many publications that now employ findings from such constructs without appropriate attention paid to reproducibility. We find that while edge-to-edge reproducibility is strongly dependent on sample size, identification of more highly connected key driver nodes in BNs can be carried out with high confidence across a range of sample sizes.