ABSTRACT
BACKGROUND: While our understanding of the single-cell gene expression patterns underlying the transformation of vascular cell types during the progression of atherosclerosis is rapidly improving, the clinical and pathophysiological relevance of these changes remains poorly understood. METHODS: Single-cell RNA sequencing data generated with SmartSeq2 (≈8000 genes/cell) in 16â 588 single cells isolated during atherosclerosis progression in Ldlr-/-Apob100/100 mice with human-like plasma lipoproteins and from humans with asymptomatic and symptomatic carotid plaques was clustered into multiple subtypes. For clinical and pathophysiological context, the advanced-stage and symptomatic subtype clusters were integrated with 135 tissue-specific (atherosclerotic aortic wall, mammary artery, liver, skeletal muscle, and visceral and subcutaneous, fat) gene-regulatory networks (GRNs) inferred from 600 coronary artery disease patients in the STARNET (Stockholm-Tartu Atherosclerosis Reverse Network Engineering Task) study. RESULTS: Advanced stages of atherosclerosis progression and symptomatic carotid plaques were largely characterized by 3 smooth muscle cells (SMCs), and 3 macrophage subtype clusters with extracellular matrix organization/osteogenic (SMC), and M1-type proinflammatory/Trem2-high lipid-associated (macrophage) phenotypes. Integrative analysis of these 6 clusters with STARNET revealed significant enrichments of 3 arterial wall GRNs: GRN33 (macrophage), GRN39 (SMC), and GRN122 (macrophage) with major contributions to coronary artery disease heritability and strong associations with clinical scores of coronary atherosclerosis severity. The presence and pathophysiological relevance of GRN39 were verified in 5 independent RNAseq data sets obtained from the human coronary and aortic artery, and primary SMCs and by targeting its top-key drivers, FRZB and ALCAM in cultured human coronary artery SMCs. CONCLUSIONS: By identifying and integrating the most gene-rich single-cell subclusters of atherosclerosis to date with a coronary artery disease framework of GRNs, GRN39 was identified and independently validated as being critical for the transformation of contractile SMCs into an osteogenic phenotype promoting advanced, symptomatic atherosclerosis.
Subject(s)
Atherosclerosis , Gene Regulatory Networks , Single-Cell Analysis , Humans , Animals , Atherosclerosis/genetics , Atherosclerosis/metabolism , Atherosclerosis/pathology , Mice , Myocytes, Smooth Muscle/metabolism , Myocytes, Smooth Muscle/pathology , Male , Plaque, Atherosclerotic , Disease Progression , Female , Macrophages/metabolism , Macrophages/pathology , Mice, Knockout , Receptors, LDL/genetics , Receptors, LDL/metabolism , Mice, Inbred C57BL , Muscle, Smooth, Vascular/metabolism , Muscle, Smooth, Vascular/pathologyABSTRACT
Genome wide association studies (GWAS) have identified thousands of single nucleotide polymorphisms (SNPs) associated with the risk of common disorders. However, since the large majority of these risk SNPs reside outside gene-coding regions, GWAS generally provide no information about causal mechanisms regarding the specific gene(s) that are affected or the tissue(s) in which these candidate gene(s) exert their effect. The 'gold standard' method for understanding causal genes and their mechanisms of action are laborious basic science studies often involving sophisticated knockin or knockout mouse lines, however, these types of studies are impractical as a high-throughput means to understand the many risk variants that cause complex diseases like coronary artery disease (CAD). As a solution, we developed a streamlined, data-driven informatics pipeline to gain mechanistic insights on complex genetic loci. The pipeline begins by understanding the SNPs in a given locus in terms of their relative location and linkage disequilibrium relationships, and then identifies nearby expression quantitative trait loci (eQTLs) to determine their relative independence and the likely tissues that mediate their disease-causal effects. The pipeline then seeks to understand associations with other disease-relevant genes, disease sub-phenotypes, potential causality (Mendelian randomization), and the regulatory and functional involvement of these genes in gene regulatory co-expression networks (GRNs). Here, we applied this pipeline to understand a cluster of SNPs associated with CAD within and immediately adjacent to the gene encoding HDAC9. Our pipeline demonstrated, and validated, that this locus is causal for CAD by modulation of TWIST1 expression levels in the arterial wall, and by also governing a GRN related to metabolic function in skeletal muscle. Our results reconciled numerous prior studies, and also provided clear evidence that this locus does not govern HDAC9 expression, structure or function. This pipeline should be considered as a powerful and efficient way to understand GWAS risk loci in a manner that better reflects the highly complex nature of genetic risk associated with common disorders.
Subject(s)
Coronary Artery Disease , Genome-Wide Association Study , Twist-Related Protein 1/metabolism , Animals , Coronary Artery Disease/genetics , Genetic Predisposition to Disease , Genome-Wide Association Study/methods , Histone Deacetylases/metabolism , Linkage Disequilibrium , Mice , Polymorphism, Single Nucleotide , Quantitative Trait Loci/genetics , Repressor Proteins/metabolismABSTRACT
Atherosclerosis is the process underlying heart attack and stroke. Despite decades of research, its pathogenesis remains unclear. Dogma suggests that atherosclerotic plaques expand primarily via the accumulation of cholesterol and inflammatory cells. However, recent evidence suggests that a substantial portion of the plaque may arise from a subset of "dedifferentiated" vascular smooth muscle cells (SMCs) which proliferate in a clonal fashion. Herein we use multicolor lineage-tracing models to confirm that the mature SMC can give rise to a hyperproliferative cell which appears to promote inflammation via elaboration of complement-dependent anaphylatoxins. Despite being extensively opsonized with prophagocytic complement fragments, we find that this cell also escapes immune surveillance by neighboring macrophages, thereby exacerbating its relative survival advantage. Mechanistic studies indicate this phenomenon results from a generalized opsonin-sensing defect acquired by macrophages during polarization. This defect coincides with the noncanonical up-regulation of so-called don't eat me molecules on inflamed phagocytes, which reduces their capacity for programmed cell removal (PrCR). Knockdown or knockout of the key antiphagocytic molecule CD47 restores the ability of macrophages to sense and clear opsonized targets in vitro, allowing for potent and targeted suppression of clonal SMC expansion in the plaque in vivo. Because integrated clinical and genomic analyses indicate that similar pathways are active in humans with cardiovascular disease, these studies suggest that the clonally expanding SMC may represent a translational target for treating atherosclerosis.
Subject(s)
Atherosclerosis/metabolism , Cloning, Molecular , Complement Activation , Myocytes, Smooth Muscle/metabolism , Phagocytosis/physiology , Animals , CD47 Antigen/metabolism , Cell Lineage , Cell Proliferation , Complement C3/genetics , Complement C3/metabolism , Female , Humans , Inflammation , Macrophages/metabolism , Male , Mice, Knockout, ApoE , Myocytes, Smooth Muscle/cytology , Plaque, Atherosclerotic/metabolism , Sequence Analysis, RNA , Up-RegulationABSTRACT
BACKGROUND: Coronary artery disease (CAD) is a multifactorial condition with both genetic and exogenous causes. The contribution of tissue-specific functional networks to the development of atherosclerosis remains largely unclear. The aim of this study was to identify and characterize central regulators and networks leading to atherosclerosis. METHODS: Based on several hundred genes known to affect atherosclerosis risk in mouse (as demonstrated in knockout models) and human (as shown by genome-wide association studies), liver gene regulatory networks were modeled. The hierarchical order and regulatory directions of genes within the network were based on Bayesian prediction models, as well as experimental studies including chromatin immunoprecipitation DNA-sequencing, chromatin immunoprecipitation mass spectrometry, overexpression, small interfering RNA knockdown in mouse and human liver cells, and knockout mouse experiments. Bioinformatics and correlation analyses were used to clarify associations between central genes and CAD phenotypes in both human and mouse. RESULTS: The transcription factor MAFF (MAF basic leucine zipper transcription factor F) interacted as a key driver of a liver network with 3 human genes at CAD genome-wide association studies loci and 11 atherosclerotic murine genes. Most importantly, expression levels of the low-density lipoprotein receptor (LDLR) gene correlated with MAFF in 600 CAD patients undergoing bypass surgery (STARNET [Stockholm-Tartu Atherosclerosis Reverse Network Engineering Task]) and a hybrid mouse diversity panel involving 105 different inbred mouse strains. Molecular mechanisms of MAFF were tested in noninflammatory conditions and showed positive correlation between MAFF and LDLR in vitro and in vivo. Interestingly, after lipopolysaccharide stimulation (inflammatory conditions), an inverse correlation between MAFF and LDLR in vitro and in vivo was observed. Chromatin immunoprecipitation mass spectrometry revealed that the human CAD genome-wide association studies candidate BACH1 (BTB domain and CNC homolog 1) assists MAFF in the presence of lipopolysaccharide stimulation with respective heterodimers binding at the MAF recognition element of the LDLR promoter to transcriptionally downregulate LDLR expression. CONCLUSIONS: The transcription factor MAFF was identified as a novel central regulator of an atherosclerosis/CAD-relevant liver network. MAFF triggered context-specific expression of LDLR and other genes known to affect CAD risk. Our results suggest that MAFF is a missing link between inflammation, lipid and lipoprotein metabolism, and a possible treatment target.
Subject(s)
Atherosclerosis/metabolism , Cholesterol/metabolism , DNA-Binding Proteins/metabolism , Inflammation/metabolism , MafF Transcription Factor/metabolism , Nuclear Proteins/metabolism , Animals , Disease Models, Animal , Humans , Male , Mice , Mice, KnockoutABSTRACT
The majority of risk loci identified by genome-wide association studies (GWAS) are in non-coding regions, hampering their functional interpretation. Instead, transcriptome-wide association studies (TWAS) identify gene-trait associations, which can be used to prioritize candidate genes in disease-relevant tissue(s). Here, we aimed to systematically identify susceptibility genes for coronary artery disease (CAD) by TWAS. We trained prediction models of nine CAD-relevant tissues using EpiXcan based on two genetics-of-gene-expression panels, the Stockholm-Tartu Atherosclerosis Reverse Network Engineering Task (STARNET) and the Genotype-Tissue Expression (GTEx). Based on these prediction models, we imputed gene expression of respective tissues from individual-level genotype data on 37,997 CAD cases and 42,854 controls for the subsequent gene-trait association analysis. Transcriptome-wide significant association (i.e. P < 3.85e-6) was observed for 114 genes. Of these, 96 resided within previously identified GWAS risk loci and 18 were novel. Stepwise analyses were performed to study their plausibility, biological function, and pathogenicity in CAD, including analyses for colocalization, damaging mutations, pathway enrichment, phenome-wide associations with human data and expression-traits correlations using mouse data. Finally, CRISPR/Cas9-based gene knockdown of two newly identified TWAS genes, RGS19 and KPTN, in a human hepatocyte cell line resulted in reduced secretion of APOB100 and lipids in the cell culture medium. Our CAD TWAS work (i) prioritized candidate causal genes at known GWAS loci, (ii) identified 18 novel genes to be associated with CAD, and iii) suggested potential tissues and pathways of action for these TWAS CAD genes.
Subject(s)
Coronary Artery Disease , Genome-Wide Association Study , Animals , Coronary Artery Disease/genetics , Genetic Predisposition to Disease , Genome-Wide Association Study/methods , Mice , Polymorphism, Single Nucleotide , TranscriptomeABSTRACT
The stress hormone cortisol modulates fuel metabolism, cardiovascular homoeostasis, mood, inflammation and cognition. The CORtisol NETwork (CORNET) consortium previously identified a single locus associated with morning plasma cortisol. Identifying additional genetic variants that explain more of the variance in cortisol could provide new insights into cortisol biology and provide statistical power to test the causative role of cortisol in common diseases. The CORNET consortium extended its genome-wide association meta-analysis for morning plasma cortisol from 12,597 to 25,314 subjects and from ~2.2 M to ~7 M SNPs, in 17 population-based cohorts of European ancestries. We confirmed the genetic association with SERPINA6/SERPINA1. This locus contains genes encoding corticosteroid binding globulin (CBG) and α1-antitrypsin. Expression quantitative trait loci (eQTL) analyses undertaken in the STARNET cohort of 600 individuals showed that specific genetic variants within the SERPINA6/SERPINA1 locus influence expression of SERPINA6 rather than SERPINA1 in the liver. Moreover, trans-eQTL analysis demonstrated effects on adipose tissue gene expression, suggesting that variations in CBG levels have an effect on delivery of cortisol to peripheral tissues. Two-sample Mendelian randomisation analyses provided evidence that each genetically-determined standard deviation (SD) increase in morning plasma cortisol was associated with increased odds of chronic ischaemic heart disease (0.32, 95% CI 0.06-0.59) and myocardial infarction (0.21, 95% CI 0.00-0.43) in UK Biobank and similarly in CARDIoGRAMplusC4D. These findings reveal a causative pathway for CBG in determining cortisol action in peripheral tissues and thereby contributing to the aetiology of cardiovascular disease.
Subject(s)
Cardiovascular Diseases/genetics , Myocardial Infarction/genetics , Transcortin/genetics , alpha 1-Antitrypsin/genetics , Adrenal Cortex Hormones/blood , Adult , Biological Specimen Banks , Cardiovascular Diseases/blood , Cardiovascular Diseases/epidemiology , Cardiovascular Diseases/pathology , Female , Gene Expression Regulation , Genetic Predisposition to Disease , Genome-Wide Association Study , Humans , Liver/metabolism , Liver/pathology , Male , Mendelian Randomization Analysis , Middle Aged , Myocardial Infarction/blood , Myocardial Infarction/pathology , Polymorphism, Single Nucleotide/genetics , Quantitative Trait Loci/genetics , United KingdomABSTRACT
The associations between diseases/traits and copy number variants (CNVs) have not been systematically investigated in genome-wide association studies (GWASs), primarily due to a lack of robust and accurate tools for CNV genotyping. Herein, we propose a novel ensemble learning framework, ensembleCNV, to detect and genotype CNVs using single nucleotide polymorphism (SNP) array data. EnsembleCNV (a) identifies and eliminates batch effects at raw data level; (b) assembles individual CNV calls into CNV regions (CNVRs) from multiple existing callers with complementary strengths by a heuristic algorithm; (c) re-genotypes each CNVR with local likelihood model adjusted by global information across multiple CNVRs; (d) refines CNVR boundaries by local correlation structure in copy number intensities; (e) provides direct CNV genotyping accompanied with confidence score, directly accessible for downstream quality control and association analysis. Benchmarked on two large datasets, ensembleCNV outperformed competing methods and achieved a high call rate (93.3%) and reproducibility (98.6%), while concurrently achieving high sensitivity by capturing 85% of common CNVs documented in the 1000 Genomes Project. Given this CNV call rate and accuracy, which are comparable to SNP genotyping, we suggest ensembleCNV holds significant promise for performing genome-wide CNV association studies and investigating how CNVs predispose to human diseases.
Subject(s)
DNA Copy Number Variations/genetics , Genotyping Techniques/methods , Machine Learning , Polymorphism, Single Nucleotide/genetics , Datasets as Topic , Genome, Human/genetics , Humans , Quality ControlABSTRACT
GWAS identified variants associated with birth weight (BW), childhood obesity (CO) and childhood BMI (CBMI), and placenta is a critical organ for fetal development and postnatal health. We examined the role of placental transcriptome and eQTLs in mediating the genetic causes for BW, CO and CBMI, and applied integrative analysis (Colocalization and MetaXcan). GWAS loci associated with BW, CO, and CBMI were substantially enriched for placenta eQTLs (6.76, 4.83 and 2.26 folds, respectively). Importantly, compared to eQTLs of adult tissues, only placental eQTLs contribute significantly to both anthropometry outcomes at birth (BW) and childhood phenotypes (CO/CBMI). Eight, six and one transcripts colocalized with BW, CO and CBMI risk loci, respectively. Our study reveals that placental transcription in utero likely plays a key role in determining postnatal body size, and as such may hold new possibilities for therapeutic interventions to prevent childhood obesity.
Subject(s)
Birth Weight/genetics , Pediatric Obesity/genetics , Placenta/metabolism , Transcriptome , Body Mass Index , Case-Control Studies , Child , Child, Preschool , Female , Fetal Development/genetics , Gene Expression Regulation , Genome-Wide Association Study , Humans , Infant, Newborn , Male , Pediatric Obesity/pathology , Polymorphism, Single Nucleotide , Pregnancy , Quantitative Trait Loci , Risk FactorsABSTRACT
Recent genome-wide association studies (GWAS) have identified multiple new loci which appear to alter coronary artery disease (CAD) risk via arterial wall-specific mechanisms. One of the annotated genes encodes LMOD1 (Leiomodin 1), a member of the actin filament nucleator family that is highly enriched in smooth muscle-containing tissues such as the artery wall. However, it is still unknown whether LMOD1 is the causal gene at this locus and also how the associated variants alter LMOD1 expression/function and CAD risk. Using epigenomic profiling we recently identified a non-coding regulatory variant, rs34091558, which is in tight linkage disequilibrium (LD) with the lead CAD GWAS variant, rs2820315. Herein we demonstrate through expression quantitative trait loci (eQTL) and statistical fine-mapping in GTEx, STARNET, and human coronary artery smooth muscle cell (HCASMC) datasets, rs34091558 is the top regulatory variant for LMOD1 in vascular tissues. Position weight matrix (PWM) analyses identify the protective allele rs34091558-TA to form a conserved Forkhead box O3 (FOXO3) binding motif, which is disrupted by the risk allele rs34091558-A. FOXO3 chromatin immunoprecipitation and reporter assays show reduced FOXO3 binding and LMOD1 transcriptional activity by the risk allele, consistent with effects of FOXO3 downregulation on LMOD1. LMOD1 knockdown results in increased proliferation and migration and decreased cell contraction in HCASMC, and immunostaining in atherosclerotic lesions in the SMC lineage tracing reporter mouse support a key role for LMOD1 in maintaining the differentiated SMC phenotype. These results provide compelling functional evidence that genetic variation is associated with dysregulated LMOD1 expression/function in SMCs, together contributing to the heritable risk for CAD.
Subject(s)
Autoantigens/genetics , Coronary Artery Disease/genetics , Cytoskeletal Proteins/genetics , Myocytes, Smooth Muscle/metabolism , Alleles , Animals , Autoantigens/metabolism , Becaplermin/metabolism , Binding Sites/genetics , Cells, Cultured , Chromosome Mapping , Coronary Artery Disease/etiology , Coronary Artery Disease/metabolism , Coronary Vessels/metabolism , Cytoskeletal Proteins/antagonists & inhibitors , Cytoskeletal Proteins/deficiency , Cytoskeletal Proteins/metabolism , Disease Models, Animal , Forkhead Box Protein O3/metabolism , Gene Knockdown Techniques , Genome-Wide Association Study , Humans , Linkage Disequilibrium , Male , Mice , Mice, Transgenic , Models, Cardiovascular , Muscle Proteins/deficiency , Muscle Proteins/genetics , Protein Binding , Quantitative Trait Loci , Risk FactorsABSTRACT
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.
Subject(s)
Gene Expression Regulation , Genetic Predisposition to Disease , Genetic Variation , Quantitative Trait Loci/genetics , Quantitative Trait, Heritable , Educational Status , Gene Regulatory Networks , Genetic Pleiotropy , Genome-Wide Association Study , Humans , Schizophrenia/geneticsABSTRACT
OBJECTIVE: Using a multi-tissue, genome-wide gene expression approach, we recently identified a gene module linked to the extent of human atherosclerosis. This atherosclerosis module was enriched with inherited risk for coronary and carotid artery disease (CAD) and overlapped with genes in the transendothelial migration of leukocyte (TEML) pathway. Among the atherosclerosis module genes, the transcription cofactor Lim domain binding 2 (LDB2) was the most connected in a CAD vascular wall regulatory gene network. Here, we used human genomics and atherosclerosis-prone mice to evaluate the possible role of LDB2 in TEML and atherosclerosis. APPROACH AND RESULTS: mRNA profiles generated from blood macrophages in patients with CAD were used to infer transcription factor regulatory gene networks; Ldlr(-/-)Apob(100/100) mice were used to study the effects of Ldb2 deficiency on TEML activity and atherogenesis. LDB2 was the most connected gene in a transcription factor regulatory network inferred from TEML and atherosclerosis module genes in CAD macrophages. In Ldlr(-/-)Apob(100/100) mice, loss of Ldb2 increased atherosclerotic lesion size ≈2-fold and decreased plaque stability. The exacerbated atherosclerosis was caused by increased TEML activity, as demonstrated in air-pouch and retinal vasculature models in vivo, by ex vivo perfusion of primary leukocytes, and by leukocyte migration in vitro. In THP1 cells, migration was increased by overexpression and decreased by small interfering RNA inhibition of LDB2. A functional LDB2 variant (rs10939673) was associated with the risk and extent of CAD across several cohorts. CONCLUSIONS: As a key driver of the TEML pathway in CAD macrophages, LDB2 is a novel candidate to target CAD by inhibiting the overall activity of TEML.
Subject(s)
Atherosclerosis/physiopathology , Carotid Artery Diseases/pathology , Chemotaxis, Leukocyte/physiology , Coronary Artery Disease/pathology , LIM Domain Proteins/physiology , Transcription Factors/physiology , Transendothelial and Transepithelial Migration/physiology , Animals , Apolipoprotein B-100/genetics , Carotid Artery Diseases/genetics , Cell Line, Tumor , Chemokine CCL2/pharmacology , Coronary Artery Disease/genetics , Gene Expression Profiling , Gene Expression Regulation , Genome-Wide Association Study , Humans , LIM Domain Proteins/deficiency , LIM Domain Proteins/genetics , Macrophages/metabolism , Mice , Mice, Knockout , RNA, Messenger/biosynthesis , Transcription Factors/deficiency , Transcription Factors/genetics , Transendothelial and Transepithelial Migration/geneticsABSTRACT
The accumulation of omics and biobank resources allows for a genome-wide understanding of the shared pathologic mechanisms between diseases and for strategies to identify drugs that could be repurposed as novel treatments. Here, we present a computational protocol, implemented as a Snakemake workflow, to identify shared transcriptional processes and screen compounds that could result in mutual benefit. This protocol also includes a description of a pharmacovigilance study designed to validate the effect of compounds using electronic health records. For complete details on the use and execution of this protocol, please refer to Gao et al.1 and Baylis et al.2.
Subject(s)
WorkflowABSTRACT
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.
Subject(s)
Cardiovascular Diseases , Glucocorticoids , Transcortin , Humans , Adipose Tissue , Cardiovascular Diseases/genetics , Gene Regulatory Networks , Genome-Wide Association Study , Heart Disease Risk Factors , Hydrocortisone , Liver , Receptors, Glucocorticoid/genetics , Risk Factors , Transcortin/geneticsABSTRACT
Senescent cells are a major contributor to age-dependent cardiovascular tissue dysfunction, but knowledge of their in vivo cell markers and tissue context is lacking. To reveal tissue-relevant senescence biology, we integrate the transcriptomes of 10 experimental senescence cell models with a 224 multi-tissue gene co-expression network based on RNA-seq data of seven tissues biopsies from â¼600 coronary artery disease (CAD) patients. We identify 56 senescence-associated modules, many enriched in CAD GWAS genes and correlated with cardiometabolic traits-which supports universality of senescence gene programs across tissues and in CAD. Cross-tissue network analyses reveal 86 candidate senescence-associated secretory phenotype (SASP) factors, including COL6A3. Experimental knockdown of COL6A3 induces transcriptional changes that overlap the majority of the experimental senescence models, with cell-cycle arrest linked to modulation of DREAM complex-targeted genes. We provide a transcriptomic resource for cellular senescence and identify candidate biomarkers, SASP factors, and potential drivers of senescence in human tissues.
Subject(s)
Cellular Senescence , Transcriptome , Humans , Transcriptome/genetics , Cellular Senescence/genetics , Phenotype , Biomarkers , Collagen , Collagen Type VI/geneticsABSTRACT
Cystatin C (CyC), a secreted cysteine protease inhibitor, has unclear biological functions. Many patients exhibit elevated plasma CyC levels, particularly during glucocorticoid (GC) treatment. This study links GCs with CyC's systemic regulation by utilizing genome-wide association and structural equation modeling to determine CyC production genetics in the UK Biobank. Both CyC production and a polygenic score (PGS) capturing predisposition to CyC production were associated with increased all-cause and cancer-specific mortality. We found that the GC receptor directly targets CyC, leading to GC-responsive CyC secretion in macrophages and cancer cells. CyC-knockout tumors displayed significantly reduced growth and diminished recruitment of TREM2+ macrophages, which have been connected to cancer immunotherapy failure. Furthermore, the CyC-production PGS predicted checkpoint immunotherapy failure in 685 patients with metastatic cancer from combined clinical trial cohorts. In conclusion, CyC may act as a GC effector pathway via TREM2+ macrophage recruitment and may be a potential target for combination cancer immunotherapy.
ABSTRACT
Background: Leukocyte progenitors derived from clonal hematopoiesis of undetermined potential (CHIP) are associated with increased cardiovascular events. However, the prevalence and functional relevance of CHIP in coronary artery disease (CAD) are unclear, and cells affected by CHIP have not been detected in human atherosclerotic plaques. Methods: CHIP mutations in blood and tissues were identified by targeted deep-DNA-sequencing (DNAseq: coverage >3,000) and whole-genome-sequencing (WGS: coverage >35). CHIP-mutated leukocytes were visualized in human atherosclerotic plaques by mutaFISH™. Functional relevance of CHIP mutations was studied by RNAseq. Results: DNAseq of whole blood from 540 deceased CAD patients of the Munich cardIovaScular StudIes biObaNk (MISSION) identified 253 (46.9%) CHIP mutation carriers (mean age 78.3 years). DNAseq on myocardium, atherosclerotic coronary and carotid arteries detected identical CHIP mutations in 18 out of 25 mutation carriers in tissue DNA. MutaFISH™ visualized individual macrophages carrying DNMT3A CHIP mutations in human atherosclerotic plaques. Studying monocyte-derived macrophages from Stockholm-Tartu Atherosclerosis Reverse Networks Engineering Task (STARNET; n=941) by WGS revealed CHIP mutations in 14.2% (mean age 67.1 years). RNAseq of these macrophages revealed that expression patterns in CHIP mutation carriers differed substantially from those of non-carriers. Moreover, patterns were different depending on the underlying mutations, e.g. those carrying TET2 mutations predominantly displayed upregulated inflammatory signaling whereas ASXL1 mutations showed stronger effects on metabolic pathways. Conclusions: Deep-DNA-sequencing reveals a high prevalence of CHIP mutations in whole blood of CAD patients. CHIP-affected leukocytes invade plaques in human coronary arteries. RNAseq data obtained from macrophages of CHIP-affected patients suggest that pro-atherosclerotic signaling differs depending on the underlying mutations. Further studies are necessary to understand whether specific pathways affected by CHIP mutations may be targeted for personalized treatment.
ABSTRACT
BACKGROUND: Ischemic preconditioning induces tolerance against ischemia-reperfusion injury prior a sustained ischemic insult. In experimental studies, exposure to hyperoxia for a limited time before ischemia induces a low-grade systemic oxidative stress and evokes an (ischemic) preconditioning-like effect of the myocardium. We hypothesised that pre-treatment by hyperoxia favours enchanced myocardial protection described by decreased release of cTn T in the 1st postoperative morning and reduces the release of inflammatory cytokines. METHODS: Forty patients with stable coronary artery disease underwent coronary artery bypass grafting with cardiopulmonary bypass. They were ventilated with 40 or >96% oxygen for 60 minutes followed by by 33 (18-59) min normoxia before cardioplegia. RESULTS: In the 1st postoperative morning concentrations of cTnT did not differ between groups ((0.44 (0.26-0.55) ng/mL in control and 0.45 (0.37-0.71) ng/mL in hyperoxia group). Sixty minutes after declamping the aorta, ratios of IL-10/IL-6 (0.73 in controls and 1.47 in hyperoxia, p = 0.03) and IL-10/TNF-α (2.91 and 8.81, resp., p = 0.015) were significantly drifted towards anti-inflammatory, whereas interleukins 6, 8and TNF-α and interferon-γ showed marked postoperative rise, but no intergroup differences were found. CONCLUSIONS: Pre-treatment by 60 minutes of hyperoxia did not reduce postoperative leak of cTn T in patients undergoing coronary artery bypass surgery. In the hyperoxia group higher release of anti-inflammatory IL-10 caused drifting of IL-10/IL-6 and IL-10/TNF-α towards anti-inflammatory.
Subject(s)
Coronary Artery Bypass , Myocardium/pathology , Oxygen/administration & dosage , Aged , Female , Humans , Inflammation/metabolism , Interleukin-10/metabolism , Interleukin-6/metabolism , Male , Middle Aged , Myocardium/metabolism , Postoperative Period , Tumor Necrosis Factor-alpha/metabolismABSTRACT
BACKGROUND: Hundreds of candidate genes have been associated with coronary artery disease (CAD) through genome-wide association studies. However, a systematic way to understand the causal mechanism(s) of these genes, and a means to prioritize them for further study, has been lacking. This represents a major roadblock for developing novel disease- and gene-specific therapies for patients with CAD. Recently, powerful integrative genomics analyses pipelines have emerged to identify and prioritize candidate causal genes by integrating tissue/cell-specific gene expression data with genome-wide association study data sets. METHODS: We aimed to develop a comprehensive integrative genomics analyses pipeline for CAD and to provide a prioritized list of causal CAD genes. To this end, we leveraged several complimentary informatics approaches to integrate summary statistics from CAD genome-wide association studies (from UK Biobank and CARDIoGRAMplusC4D) with transcriptomic and expression quantitative trait loci data from 9 cardiometabolic tissue/cell types in the STARNET study (Stockholm-Tartu Atherosclerosis Reverse Network Engineering Task). RESULTS: We identified 162 unique candidate causal CAD genes, which exerted their effect from between one and up to 7 disease-relevant tissues/cell types, including the arterial wall, blood, liver, skeletal muscle, adipose, foam cells, and macrophages. When their causal effect was ranked, the top candidate causal CAD genes were CDKN2B (associated with the 9p21.3 risk locus) and PHACTR1; both exerting their causal effect in the arterial wall. A majority of candidate causal genes were represented in cross-tissue gene regulatory co-expression networks that are involved with CAD, with 22/162 being key drivers in those networks. CONCLUSIONS: We identified and prioritized candidate causal CAD genes, also localizing their tissue(s) of causal effect. These results should serve as a resource and facilitate targeted studies to identify the functional impact of top causal CAD genes.