Plasma cortisol-linked gene networks in hepatic and adipose tissues implicate corticosteroid-binding globulin in modulating tissue glucocorticoid action and cardiovascular risk

Genome-wide association meta-analysis (GWAMA) by the Cortisol Network (CORNET) consortium identified genetic variants spanning the SERPINA6/SERPINA1 locus on chromosome 14 associated with morning plasma cortisol, cardiovascular disease (CVD), and SERPINA6 mRNA expression encoding corticosteroid-binding globulin (CBG) in the liver. These and other findings indicate that higher plasma cortisol levels are causally associated with CVD; however, the mechanisms by which variations in CBG lead to CVD are undetermined. Using genomic and transcriptomic data from The Stockholm Tartu Atherosclerosis Reverse Networks Engineering Task (STARNET) study, we identified plasma cortisol-linked single-nucleotide polymorphisms (SNPs) that are trans-associated with genes from seven different vascular and metabolic tissues, finding the highest representation of trans-genes in the liver, subcutaneous fat, and visceral abdominal fat, [false discovery rate (FDR) = 15%]. We identified a subset of cortisol-associated trans-genes that are putatively regulated by the glucocorticoid receptor (GR), the primary transcription factor activated by cortisol. Using causal inference, we identified GR-regulated trans-genes that are responsible for the regulation of tissue-specific gene networks. Cis-expression Quantitative Trait Loci (eQTLs) were used as genetic instruments for identification of pairwise causal relationships from which gene networks could be reconstructed. Gene networks were identified in the liver, subcutaneous fat, and visceral abdominal fat, including a high confidence gene network specific to subcutaneous adipose (FDR = 10%) under the regulation of the interferon regulatory transcription factor, IRF2. These data identify a plausible pathway through which variation in the liver CBG production perturbs cortisol-regulated gene networks in peripheral tissues and thereby promote CVD

On reliable discovery of molecular signatures

<p>Abstract</p> <p>Background</p> <p>Molecular signatures are sets of genes, proteins, genetic variants or other variables that can be used as markers for a particular phenotype. Reliable signature discovery methods could yield valuable insight into cell biology and mechanisms of human disease. However, it is currently not clear how to control error rates such as the false discovery rate (FDR) in signature discovery. Moreover, signatures for cancer gene expression have been shown to be unstable, that is, difficult to replicate in independent studies, casting doubts on their reliability.</p> <p>Results</p> <p>We demonstrate that with modern prediction methods, signatures that yield accurate predictions may still have a high FDR. Further, we show that even signatures with low FDR may fail to replicate in independent studies due to limited statistical power. Thus, neither stability nor predictive accuracy are relevant when FDR control is the primary goal. We therefore develop a general statistical hypothesis testing framework that for the first time provides FDR control for signature discovery. Our method is demonstrated to be correct in simulation studies. When applied to five cancer data sets, the method was able to discover molecular signatures with 5% FDR in three cases, while two data sets yielded no significant findings.</p> <p>Conclusion</p> <p>Our approach enables reliable discovery of molecular signatures from genome-wide data with current sample sizes. The statistical framework developed herein is potentially applicable to a wide range of prediction problems in bioinformatics.</p

Correction to: Integrative analysis of loss-of-function variants in clinical and genomic data reveals novel genes associated with cardiovascular traits

Erratum for Integrative analysis of loss-of-function variants in clinical and genomic data reveals novel genes associated with cardiovascular traits. [BMC Med Genomics. 2019

Detecting multivariate differentially expressed genes

<p>Abstract</p> <p>Background</p> <p>Gene expression is governed by complex networks, and differences in expression patterns between distinct biological conditions may therefore be complex and multivariate in nature. Yet, current statistical methods for detecting differential expression merely consider the univariate difference in expression level of each gene in isolation, thus potentially neglecting many genes of biological importance.</p> <p>Results</p> <p>We have developed a novel algorithm for detecting multivariate expression patterns, named Recursive Independence Test (RIT). This algorithm generalizes differential expression testing to more complex expression patterns, while still including genes found by the univariate approach. We prove that RIT is consistent and controls error rates for small sample sizes. Simulation studies confirm that RIT offers more power than univariate differential expression analysis when multivariate effects are present. We apply RIT to gene expression data sets from diabetes and cancer studies, revealing several putative disease genes that were not detected by univariate differential expression analysis.</p> <p>Conclusion</p> <p>The proposed RIT algorithm increases the power of gene expression analysis by considering multivariate effects while retaining error rate control, and may be useful when conventional differential expression tests yield few findings.</p

Histone deacetylase 9 promotes endothelial to mesenchymal transition and an unfavorable atherosclerotic plaque phenotype

Endothelial-mesenchymal transition (EndMT) is associated with various cardiovascular diseases and in particular with atherosclerosis and plaque instability. However, the molecular pathways that govern EndMT are poorly defined. Specifically, the role of epigenetic factors and histone deacetylases (HDACs) in controlling EndMT and the atherosclerotic plaque phenotype remains unclear. Here, we identified histone deacetylation, specifically that mediated by HDAC9 (a class IIa HDAC), as playing an important role in both EndMT and atherosclerosis. Using in vitro models, we found class IIa HDAC inhibition sustained the expression of endothelial proteins and mitigated the increase in mesenchymal proteins, effectively blocking EndMT. Similarly, ex vivo genetic knockout of Hdac9 in endothelial cells prevented EndMT and preserved a more endothelial-like phenotype. In vivo, atherosclerosis-prone mice with endothelial-specific Hdac9 knockout showed reduced EndMT and significantly reduced plaque area. Furthermore, these mice displayed a more favorable plaque phenotype, with reduced plaque lipid content and increased fibrous cap thickness. Together, these findings indicate that HDAC9 contributes to vascular pathology by promoting EndMT. Our study provides evidence for a pathological link among EndMT, HDAC9, and atherosclerosis and suggests that targeting of HDAC9 may be beneficial for plaque stabilization or slowing the progression of atherosclerotic disease

Transcriptomic-based clustering of human atherosclerotic plaques identifies subgroups with different underlying biology and clinical presentation

Histopathological studies have revealed key processes of atherosclerotic plaque thrombosis. However, the diversity and complexity of lesion types highlight the need for improved sub-phenotyping. Here we analyze the gene expression profiles of 654 advanced human carotid plaques. The unsupervised, transcriptome-driven clustering revealed five dominant plaque types. These plaque phenotypes were associated with clinical presentation and showed differences in cellular compositions. Validation in coronary segments showed that the molecular signature of these plaques was linked to coronary ischemia. One of the plaque types with the most severe clinical symptoms pointed to both inflammatory and fibrotic cell lineages. Further, we did a preliminary analysis of potential circulating biomarkers that mark the different plaques phenotypes. In conclusion, the definition of the plaque at risk for a thrombotic event can be fine-tuned by in-depth transcriptomic-based phenotyping. These differential plaque phenotypes prove clinically relevant for both carotid and coronary artery plaques and point to distinct underlying biology of symptomatic lesions

Female Gene Networks Are Expressed in Myofibroblast-Like Smooth Muscle Cells in Vulnerable Atherosclerotic Plaques

BACKGROUND: Women presenting with coronary artery disease more often present with fibrous atherosclerotic plaques, which are currently understudied. Phenotypically modulated smooth muscle cells (SMCs) contribute to atherosclerosis in women. How these phenotypically modulated SMCs shape female versus male plaques is unknown. METHODS: Gene regulatory networks were created using RNAseq gene expression data from human carotid atherosclerotic plaques. The networks were prioritized based on sex bias, relevance for smooth muscle biology, and coronary artery disease genetic enrichment. Network expression was linked to histologically determined plaque phenotypes. In addition, their expression in plaque cell types was studied at single-cell resolution using single-cell RNAseq. Finally, their relevance for disease progression was studied in female and male Apoe -/- mice fed a Western diet for 18 and 30 weeks. RESULTS: Here, we identify multiple sex-stratified gene regulatory networks from human carotid atherosclerotic plaques. Prioritization of the female networks identified 2 main SMC gene regulatory networks in late-stage atherosclerosis. Single-cell RNA sequencing mapped these female networks to 2 SMC phenotypes: a phenotypically modulated myofibroblast-like SMC network and a contractile SMC network. The myofibroblast-like network was mostly expressed in plaques that were vulnerable in women. Finally, the mice ortholog of key driver gene MFGE8 (milk fat globule EGF and factor V/VIII domain containing) showed retained expression in advanced plaques from female mice but was downregulated in male mice during atherosclerosis progression. CONCLUSIONS: Female atherosclerosis is characterized by gene regulatory networks that are active in fibrous vulnerable plaques rich in myofibroblast-like SMCs

Carotid Plaque Age Is a Feature of Plaque Stability Inversely Related to Levels of Plasma Insulin

C-declination curve (a result of the atomic bomb tests in the 1950s and 1960s) to determine the average biological age of carotid plaques.C content by accelerator mass spectrometry. The average plaque age (i.e. formation time) was 9.6±3.3 years. All but two plaques had formed within 5–15 years before surgery. Plaque age was not associated with the chronological ages of the patients but was inversely related to plasma insulin levels (p = 0.0014). Most plaques were echo-lucent rather than echo-rich (2.24±0.97, range 1–5). However, plaques in the lowest tercile of plaque age (most recently formed) were characterized by further instability with a higher content of lipids and macrophages (67.8±12.4 vs. 50.4±6.2, p = 0.00005; 57.6±26.1 vs. 39.8±25.7, p<0.0005, respectively), less collagen (45.3±6.1 vs. 51.1±9.8, p<0.05), and fewer smooth muscle cells (130±31 vs. 141±21, p<0.05) than plaques in the highest tercile. Microarray analysis of plaques in the lowest tercile also showed increased activity of genes involved in immune responses and oxidative phosphorylation.C, can improve our understanding of carotid plaque stability and therefore risk for clinical complications. Our results also suggest that levels of plasma insulin might be involved in determining carotid plaque age

Multi-ancestry genetic analysis of gene regulation in coronary arteries prioritizes disease risk loci

Genome-wide association studies (GWASs) have identified hundreds of risk loci for coronary artery disease (CAD). However, non-European populations are underrepresented in GWASs, and the causal gene-regulatory mechanisms of these risk loci during atherosclerosis remain unclear. We incorporated local ancestry and haplotypes to identify quantitative trait loci for expression (eQTLs) and splicing (sQTLs) in coronary arteries from 138 ancestrally diverse Americans. Of 2,132 eQTL-associated genes (eGenes), 47% were previously unreported in coronary artery; 19% exhibited cell-type-specific expression. Colocalization revealed subgroups of eGenes unique to CAD and blood pressure GWAS. Fine-mapping highlighted additional eGenes, including TBX20 and IL5. We also identified sQTLs for 1,690 genes, among which TOR1AIP1 and ULK3 sQTLs demonstrated the importance of evaluating splicing to accurately identify disease-relevant isoform expression. Our work provides a patient-derived coronary artery eQTL resource and exemplifies the need for diverse study populations and multifaceted approaches to characterize gene regulation in disease processes