Search

Leukocytes and resident cells in the arterial wall contribute to atherosclerosis, especially at sites of disturbed blood flow. Expression of endothelial Tie1 receptor tyrosine kinase is enhanced at these sites, and attenuation of its expression reduces atherosclerotic burden and decreases inflammation. However, Tie2 tyrosine kinase function in atherosclerosis is unknown. Here we provide genetic evidence from humans and from an atherosclerotic mouse model to show that TIE2 is associated with protection from coronary artery disease. We show that deletion of Tie2, or both Tie2 and Tie1, in the arterial endothelium promotes atherosclerosis by increasing Foxo1 nuclear localization, endothelial adhesion molecule expression and accumulation of immune cells. We also show that Tie2 is expressed in a subset of aortic fibroblasts, and its silencing in these cells increases expression of inflammation-related genes. Our findings indicate that unlike Tie1, the Tie2 receptor functions as the dominant endothelial angiopoietin receptor that protects from atherosclerosis.

BackgroundNon-alcoholic fatty liver disease (NAFLD) is a fast-growing, underdiagnosed, epidemic. We hypothesise that obesity-related inflammation compromises adipose tissue functions, preventing efficient fat storage, and thus driving ectopic fat accumulation into the liver.MethodsTo identify adipose-based mechanisms and potential serum biomarker candidates (SBCs) for NAFLD, we utilise dual-tissue RNA-sequencing (RNA-seq) data in adipose tissue and liver, paired with histology-based NAFLD diagnosis, from the same individuals in a cohort of obese individuals. We first scan for genes that are differentially expressed (DE) for NAFLD in obese individuals' subcutaneous adipose tissue but not in their liver; encode proteins secreted to serum; and show preferential adipose expression. Then the identified genes are filtered to key adipose-origin NAFLD genes by best subset analysis, knockdown experiments during human preadipocyte differentiation, recombinant protein treatment experiments in human liver HepG2 cells, and genetic analysis.FindingsWe discover a set of genes, including 10 SBCs, that may modulate NAFLD pathogenesis by impacting adipose tissue function. Based on best subset analysis, we further follow-up on two SBCs CCDC80 and SOD3 by knockdown in human preadipocytes and subsequent differentiation experiments, which show that they modulate crucial adipogenesis genes, LPL, SREBPF1, and LEP. We also show that treatment of the liver HepG2 cells with the CCDC80 and SOD3 recombinant proteins impacts genes related to steatosis and lipid processing, including PPARA, NFE2L2, and RNF128. Finally, utilizing the adipose NAFLD DE gene cis-regulatory variants associated with serum triglycerides (TGs) in extensive genome-wide association studies (GWASs), we demonstrate a unidirectional effect of serum TGs on NAFLD with Mendelian Randomization (MR) analysis. We also demonstrate that a single SNP regulating one of the SBC genes, rs2845885, produces a significant MR result by itself. This supports the conclusion that genetically regulated adipose expression of the NAFLD DE genes may contribute to NAFLD through changes in serum TG levels.InterpretationOur results from the dual-tissue transcriptomics screening improve the understanding of obesity-related NAFLD by providing a targeted set of 10 adipose tissue-active genes as new serum biomarker candidates for the currently grossly underdiagnosed fatty liver disease.FundingThe work was supported by NIH grants R01HG010505 and R01DK132775. The Genotype-Tissue Expression (GTEx) Project was supported by the Common Fund of the Office of the Director of the National Institutes of Health, and by NCI, NHGRI, NHLBI, NIDA, NIMH, and NINDS. The KOBS study (J. P.) was supported by the Finnish Diabetes Research Foundation, Kuopio University Hospital Project grant (EVO/VTR grants 2005-2019), and the Academy of Finland grant (Contract no. 138006). This study was funded by the European Research Council under the European Union's Horizon 2020 research and innovation program (Grant No. 802825 to M. U. K.). K. H. P. was funded by the Academy of Finland (grant numbers 272376, 266286, 314383, and 335443), the Finnish Medical Foundation, Gyllenberg Foundation, Novo Nordisk Foundation (grant numbers NNF10OC1013354, NNF17OC0027232, and NNF20OC0060547), Finnish Diabetes Research Foundation, Finnish Foundation for Cardiovascular Research, University of Helsinki, and Helsinki University Hospital and Government Research Funds. I. S. was funded by the Instrumentarium Science Foundation. Personal grants to U. T. A. were received from the Matti and Vappu Maukonen Foundation, Ella och Georg Ehrnrooths Stiftelse and the Finnish Foundation for Cardiovascular Research.

Background: Sporadic venous malformation (VM) and angiomatosis of soft tissue (AST) are benign, congenital vascular anomalies affecting venous vasculature. Depending on the size and location of the lesion, symptoms vary from motility disturbances to pain and disfigurement. Due to high recurrence of the lesions more effective therapies are needed. Methods: As targeting stromal cells has been an emerging concept in anti-angiogenic therapies, here, by using VM/AST patient samples, RNA-sequencing, cell culture techniques and a xenograft mouse model, we investigated the crosstalk of endothelial cells (EC) and fibroblasts and its effect on vascular lesion growth. Results: We report, for the first time, expression and secretion of transforming growth factor A (TGFA) in ECs or intervascular stromal cells in AST and VM lesions. TGFA induced secretion of VEGF-A paracrinally, and regulated EC proliferation. Oncogenic PIK3CA variant in p.H1047R, a common somatic mutation found in these lesions, increased TGFA expression, enrichment of hallmark hypoxia, and in a mouse xenograft model, lesion size and vascularization. Treatment with afatinib, a pan-ErbB tyrosine-kinase inhibitor, decreased vascularization and lesion size in mouse xenograft model with ECs expressing oncogenic PIK3CA p.H1047R variant and fibroblasts. Conclusions: Based on the data, we suggest that targeting of both intervascular stromal cells and ECs is a potential treatment strategy for vascular lesions having a fibrous component. Funding: Academy of Finland, Ella and Georg Ehnrooth foundation, the ERC grants, Sigrid Jusélius Foundation, Finnish Foundation for Cardiovascular Research, Jane and Aatos Erkko Foundation, and Department of Musculosceletal and Plastic Surgery, Helsinki University Hospital.

Coronary artery disease (CAD) is a pandemic disease where up to half of the risk is explained by genetic factors. Advanced insights into the genetic basis of CAD require deeper understanding of the contributions of different cell types, molecular pathways, and genes to disease heritability. Here, we investigate the biological diversity of atherosclerosis-associated cell states and interrogate their contribution to the genetic risk of CAD by using single-cell and bulk RNA sequencing (RNA-seq) of mouse and human lesions. We identified 12 disease-associated cell states that we characterized further by gene set functional profiling, ligand-receptor prediction, and transcription factor inference. Importantly, Vcam1+ smooth muscle cell state genes contributed most to SNP-based heritability of CAD. In line with this, genetic variants near smooth muscle cell state genes and regulatory elements explained the largest fraction of CAD-risk variance between individuals. Using this information for variant prioritization, we derived a hybrid polygenic risk score (PRS) that demonstrated improved performance over a classical PRS. Our results provide insights into the biological mechanisms associated with CAD risk, which could make a promising contribution to precision medicine and tailored therapeutic interventions in the future. publishedVersion

Nuclear factor erythroid-2 related factor 2 (Nrf2) is a transcription factor that regulates protection against a wide variety of toxic insults to cells, including cytotoxic cancer chemotherapeutic drugs. Many lung cancer cells harbor a mutation in either Nrf2 or its inhibitor Keap1 resulting in permanent activation of Nrf2 and chemoresistance. In this study, we sought to examine whether this attribute could be exploited in cancer suicide gene therapy by using a lentiviral (LV) vector expressing herpes simplex virus thymidine kinase (HSV-TK/GCV) under the regulation of antioxidant response element (ARE), a cis-acting enhancer sequence that binds Nrf2. In human lung adenocarcinoma cells in which Nrf2 is constitutively overexpressed, ARE activity was found to be high under basal conditions. In this setting, ARE-HSV-TK was more effective than a vector in which HSV-TK expression was driven by a constitutively active promoter. In a mouse xenograft model of lung cancer, suicide gene therapy with LV-ARE-TK/GCV was effective compared with LV-PGK-TK/GCV in reducing tumor size. We conclude that ARE-regulated HSV-TK/GCV therapy offers a promising approach for suicide cancer gene therapy in cells with high constitutive ARE activity, permitting a greater degree of therapeutic targeting to those cells. Cancer Res; 72(23); 6227–35. ©2012 AACR.

PDF file - 57K

PDF file - 324K, The effect of ganciclovir and doxorubicin on ARE activity

Brain functionality resides on finely tuned regulation of gene expression by networks of non-coding RNAs (ncRNAs) such as the one composed by the circular RNA ciRS-7, the microRNA miR-7 and the long non-coding RNA Cyrano. However, very little is known on how this network regulates stress responses in neurodegeneration. Here we describe ischemia induced alterations in the ncRNA network bothin vitroandin vivoand in transgenic mice lacking ciRS-7 or miR-7. Our data show that cortical neurons downregulate ciRS-7 and Cyrano and upregulate miR-7 expression upon ischemic insults. Mice lacking ciRS-7 show reduced lesion size and motor impairment, whilst the absence of miR-7 alone leads to an increase in the ischemia induced neuronal death. Moreover, miR-7 levels in pyramidal excitatory neurons regulate dendrite morphology and glutamatergic signaling suggesting a potential molecular link to thein vivophenotype. Our data reveal a new endogenous mechanism by which ciRS-7 and miR-7 regulate the outcome of ischemic stroke and shed new light into the pathophysiological roles of intracellular networks of non-coding RNAs in the brain.

Background: Coronary artery disease (CAD) is the leading cause of death worldwide. Recent meta-analyses of genome-wide association studies have identified over 175 loci associated with CAD. The majority of these loci are in noncoding regions and are predicted to regulate gene expression. Given that vascular smooth muscle cells (SMCs) play critical roles in the development and progression of CAD, we aimed to identify the subset of the CAD loci associated with the regulation of transcription in distinct SMC phenotypes. Methods: We measured gene expression in SMCs isolated from the ascending aortas of 151 heart transplant donors of various genetic ancestries in quiescent or proliferative conditions and calculated the association of their expression and splicing with ~6.3 million imputed single-nucleotide polymorphism markers across the genome. Results: We identified 4910 expression and 4412 splicing quantitative trait loci (sQTLs) representing regions of the genome associated with transcript abundance and splicing. A total of 3660 expression quantitative trait loci (eQTLs) had not been observed in the publicly available Genotype-Tissue Expression dataset. Further, 29 and 880 eQTLs were SMC-specific and sex-biased, respectively. We made these results available for public query on a user-friendly website. To identify the effector transcript(s) regulated by CAD loci, we used 4 distinct colocalization approaches. We identified 84 eQTL and 164 sQTL that colocalized with CAD loci, highlighting the importance of genetic regulation of mRNA splicing as a molecular mechanism for CAD genetic risk. Notably, 20% and 35% of the eQTLs were unique to quiescent or proliferative SMCs, respectively. One CAD locus colocalized with a sex-specific eQTL ( TERF2IP ), and another locus colocalized with SMC-specific eQTL ( ALKBH8 ). The most significantly associated CAD locus, 9p21, was an sQTL for the long noncoding RNA CDKN2B-AS1 , also known as ANRIL , in proliferative SMCs. Conclusions: Collectively, our results provide evidence for the molecular mechanisms of genetic susceptibility to CAD in distinct SMC phenotypes.

Background: Sporadic venous malformation (VM) and angiomatosis of soft tissue (AST) are benign, congenital vascular anomalies affecting venous vasculature. Depending on the size and location of the lesion, symptoms vary from motility disturbances to pain and disfigurement. Due to high recurrence of the lesions more effective therapies are needed. Methods: As targeting stromal cells has been an emerging concept in anti-angiogenic therapies, here, by using VM/AST patient samples, RNA-sequencing, cell culture techniques and a xenograft mouse model, we investigated the crosstalk of endothelial cells (EC) and fibroblasts and its effect on vascular lesion growth. Results: We report, for the first time, expression and secretion of transforming growth factor A (TGFA) in ECs or intervascular stromal cells in AST and VM lesions. TGFA induced secretion of VEGF-A paracrinally, and regulated EC proliferation. Oncogenic PIK3CA variant in p.H1047R, a common somatic mutation found in these lesions, increased TGFA expression, enrichment of hallmark hypoxia, and in a mouse xenograft model, lesion size and vascularization. Treatment with afatinib, a pan-ErbB tyrosine-kinase inhibitor, decreased vascularization and lesion size in mouse xenograft model with ECs expressing oncogenic PIK3CA p.H1047R variant and fibroblasts. Conclusions: Based on the data, we suggest that targeting of both intervascular stromal cells and ECs is a potential treatment strategy for vascular lesions having a fibrous component. Funding: Academy of Finland, Ella and Georg Ehnrooth foundation, the ERC grants, Sigrid Jusélius Foundation, Finnish Foundation for Cardiovascular Research, Jane and Aatos Erkko Foundation, and Department of Musculosceletal and Plastic Surgery, Helsinki University Hospital.

Supplemental Digital Content is available in the text., Objective: Atherosclerosis is the underlying cause of most cardiovascular diseases. The main cell types associated with disease progression in the vascular wall are endothelial cells, smooth muscle cells, and macrophages. Although their role in atherogenesis has been extensively described, molecular mechanisms underlying gene expression changes remain unknown. The objective of this study was to characterize microRNA (miRNA)-related regulatory mechanisms taking place in the aorta during atherosclerosis. Approach and Results: We analyzed the miRNA expression changes in primary human aortic endothelial cells and human umbilical vein endothelial cells, human aortic smooth muscle cells, and macrophages (CD14+) under various proatherogenic stimuli by integrating GRO-seq, miRNA-seq, and RNA-seq data. Despite the highly cell-type-specific expression of multi-variant primary miRNAs, the majority of mature miRNAs were found to be common to all cell types and dominated by 2 to 5 abundant miRNA species. We demonstrate that transcription contributes significantly to the mature miRNA levels although this is dependent on miRNA stability. An analysis of miRNA effects in relation to target mRNA pools highlighted pathways and targets through which miRNAs could affect atherogenesis in a cell-type-dependent manner. Finally, we validate miR-100-5p as a cell-type specific regulator of inflammatory and HIPPO-YAP/TAZ-pathways. Conclusions: This integrative approach allowed us to characterize miRNA dynamics in response to a proatherogenic stimulus and identify potential mechanisms by which miRNAs affect atherogenesis in a cell-type-specific manner.

SummaryThe NRF2 pathway is frequently activated in various cancer types, yet a comprehensive analysis of its effects across different malignancies is currently lacking. We developed a robust NRF2 activity metric and utilized it to conduct a pan-cancer wide analysis of oncogenic NRF2 signaling. We identified a distinct immunoevasive phenotype where high NRF2 activity is associated with low interferon-gamma (IFNγ), HLA-I expression and T-cell infiltration spanning non-small cell lung cancer (NSCLC) and squamous malignancies of head and neck area, cervix and esophagus. In squamous cell cancers, NRF2 overactive tumors comprise a molecular phenotype with SOX2/TP63 amplification, TP53 mutation and CDKN2A loss. These immune-cold NRF2 hyperactive diseases are associated with upregulation of immunomodulatory NAMPT, WNT5A, SPP1, SLC7A11 and SLC2A1 that represent candidate NRF2 target genes, suggesting direct modulation of the tumor immune milieu. Based on single-cell mRNA data, coupled with a priori information on intercellular ligand-receptor interactions, cancer cells of this subtype exhibit decreased expression of IFNγ responsive ligands, and increased expression of immunosuppressive ligands NAMPT, SPP1 and WNT5A that mediate signaling in intercellular crosstalk. As we observed differential cytokine mRNA expression with IFNγ treatment in NSCLC adenocarcinoma subtype, we explored the cytokine secretome in vitro. We found that secreted neutrophil chemoattractants interleukin-8 (CXCL8) and ENA-78 (CXCL5) are elevated in NRF2 overactive cells, suggesting contribution of immunosuppressive neutrophils in NRF2 driven immune escape. Importantly, as overactive NRF2 is associated with immune-cold characteristics, our results highlight the utility of NRF2 pathway activation as a putative biomarker for stratifying immune-checkpoint blockade responders and non-responders across NSCLC and squamous cancers.Abstract Figure

The BMP/TGFβ-Smad, Notch and VEGF signaling guides formation of endothelial tip and stalk cells. However, the crosstalk of bone morphogenetic proteins (BMPs) and vascular endothelial growth factor receptor 2 (VEGFR2) signaling has remained largely unknown. We demonstrate that BMP family members regulate VEGFR2 and Notch signaling, and act via TAZ-Hippo signaling pathway. BMPs were found to be regulated after VEGF gene transfer in C57/Bl6 mice and in a porcine myocardial ischemia model. BMPs 2/4/6 were identified as endothelium-specific targets of VEGF. BMP2 modulated VEGF-mediated endothelial sprouting via Delta like Canonical Notch Ligand 4 (DLL4). BMP6 modulated VEGF signaling by regulating VEGFR2 expression and acted via Hippo signaling effector TAZ, known to regulate cell survival/proliferation, and to be dysregulated in cancer. In a matrigel plug assay in nude mice BMP6 was further demonstrated to induce angiogenesis. BMP6 is the first member of BMP family found to directly regulate both Hippo signaling and neovessel formation. It may thus serve as a target in pro/anti-angiogenic therapies. Electronic supplementary material The online version of this article (10.1007/s10456-020-09748-4) contains supplementary material, which is available to authorized users.

Olfaction is orchestrated by olfactory mucosal cells located in the upper nasal cavity. Olfactory dysfunction manifests early in several neurodegenerative disorders including Alzheimer’s disease, however, disease-related alterations to the olfactory mucosal cells remain poorly described. The aim of this study was to evaluate the olfactory mucosa differences between cognitively healthy individuals and Alzheimer’s disease patients. We report increased amyloid-beta secretion in Alzheimer’s disease olfactory mucosal cells and detail cell-type-specific gene expression patterns, unveiling 240 differentially expressed disease-associated genes compared to the cognitively healthy controls, and five distinct cell populations. Overall, alterations of RNA and protein metabolism, inflammatory processes, and signal transduction were observed in multiple cell populations, suggesting their role in Alzheimer’s disease-related olfactory mucosa pathophysiology. Furthermore, the single-cell RNA-sequencing proposed alterations in gene expression of mitochondrially located genes in AD OM cells, which were verified by functional assays, demonstrating altered mitochondrial respiration and a reduction of ATP production. Our results reveal disease-related changes of olfactory mucosal cells in Alzheimer’s disease and demonstrate the utility of single-cell RNA sequencing data for investigating molecular and cellular mechanisms associated with the disease.

Coronary artery disease (CAD) is the leading cause of death worldwide. Recent meta-analyses of genome-wide association studies (GWAS) have identified over 175 loci associated with CAD. The majority of these loci are in non-coding regions and are predicted to regulate gene expression. Given that vascular smooth muscle cells (SMCs) play critical roles in the development and progression of CAD, we hypothesized that a subset of the CAD GWAS risk loci are associated with the regulation of transcription in distinct SMC phenotypes. Here, we measured gene expression in SMCs isolated from the ascending aortas of 151 ethnically diverse heart transplant donors in quiescent or proliferative conditions and calculated the association of their expression and splicing with ∼6.3 million imputed single nucleotide polymorphism (SNP) markers across the genome. We identified 4,910 expression and 4,412 splice quantitative trait loci (sQTL) that represent regions of the genome associated with transcript abundance and splicing. 3,660 of the eQTLs had not been observed in the publicly available Genotype-Tissue Expression dataset. Further, 29 and 880 of the eQTLs were SMC- and sex-specific, respectively. To identify the effector transcript(s) regulated by CAD GWAS loci, we used four distinct colocalization approaches and identified 84 eQTL and 164 sQTLs that colocalized with CAD loci, highlighting the importance of genetic regulation of mRNA splicing as a molecular mechanism for CAD genetic risk. Notably, 20% and 35% of the eQTLs were unique to quiescent or proliferative SMCs, respectively. Two CAD loci colocalized with a SMC sex-specific eQTL (AL160313.1 and TERF2IP) and another locus colocalized with SMC-specific eQTL (ALKBH8). Also, 27% and 37% of the sQTLs were unique to quiescent or proliferative SMCs, respectively. The most significantly associated CAD locus, 9p21, was an sQTL for the long non-coding RNA CDKN2B-AS1, also known as ANRIL, in proliferative SMCs. Collectively, these results provide evidence for the molecular mechanisms of genetic susceptibility to CAD in distinct SMC phenotypes.

MicroRNAs (miRNAs) are small RNA molecules that act as regulators of gene expression through targeted mRNA degradation. They are involved in many biological and pathophysiological processes and are widely studied as potential biomarkers and therapeutics agents for human diseases, including cardiovascular disorders. Recently discovered isoforms of miRNAs (isomiRs) exist in high quantities and are very diverse. Despite having few differences with their corresponding reference miRNAs, they display specific functions and expression profiles, across tissues and conditions. However, they are still overlooked and understudied, as we lack a comprehensive view on their condition-specific regulation and impact on differential expression analysis. Here, we show that isomiRs can have major effects on differential expression analysis results, as their expression is independent of their host miRNA genes or reference sequences. We present two miRNA-seq datasets from human umbilical vein endothelial cells, and assess isomiR expression in response to senescence and compartment-specificity (nuclear/cytosolic) under hypoxia. We compare three different methods for miRNA analysis, including isomiR-specific analysis, and show that ignoring isomiRs induces major biases in differential expression. Moreover, isomiR analysis permits higher resolution of complex signal dissection, such as the impact of hypoxia on compartment localization, and differential isomiR type enrichments between compartments. Finally, we show important distribution differences across conditions, independently of global miRNA expression signals. Our results raise concerns over the quasi exclusive use of miRNA reference sequences in miRNA-seq processing and experimental assays. We hope that our work will guide future isomiR expression studies, which will correct some biases introduced by golden standard analysis, improving the resolution of such assays and the biological significance of their downstream studies.

Increased glutamine metabolism by macrophages is associated with development of atherosclerotic lesions. Positron emission tomography/computed tomography (PET/CT) with a glutamine analog (2S,4

Vascular Smooth Muscle Cells (SMCs) exhibit a significant degree of phenotypic plasticity in atherosclerotic lesions. Previously, this phenotypic plasticity, and specifically the re-differentiation of SMCs into macrophages, has been modeled in vitro using cholesterol treatment. We performed a meta-analysis of five different murine SMC lineage tracing single-cell RNA sequencing (scRNA-seq) experiments performed by four different research groups; we identified 24 different cell types and find considerable differences in the proportion of recovered normal ‘contractile’ SMCs and SMC-derived macrophages between experiments. In addition, we performed bulk RNA-seq on cultured murine SMCs with and without cholesterol treatment for 0, 24, or 48 hours. We find that cholesterol treatment regulates transcripts associated with the induction of the unfolded protein response, reduction of cholesterol biosynthesis, and reduction in cell proliferation. Additionally, chromatin immunoprecipitation followed by high-throughput sequencing for H3K27ac performed on in vitro SMCs treated and untreated for 0, 24, or 48 hours reveals that AP-1 and ATF4 binding motifs are found more frequently at cholesterol treatment specific H3K27ac chromatin sites. Finally, we compared transcriptomic data between in vivo vascular cell types identified in our meta-analysis and in vitro SMCs. We additionally included scRNA-seq data and microarray data collected for in vitro SMCs and bulk RNA-seq data collected for M1, M2, and oxLDL treated macrophages. We find that compared to the phenotypic plasticity observed in vivo, cholesterol exhibits modest effects on the transcriptomes of in vitro SMCs and does not reproduce the changes observed in vivo, while in vitro macrophages, especially oxLDL treated macrophages, more closely approximate their in vivo counterparts. We conclude that differences between SMC lineage-tracing scRNA-seq experiments obscures the degree to which SMC-derived macrophages are present in murine atherosclerotic plaques, that cholesterol has modest effects on in vitro SMC transcriptomes, and that the effects of cholesterol do not reproduce the phenotypic plasticity observed in vivo.

MicroRNAs are well characterized in their role in silencing gene expression by targeting 3´-UTR of mRNAs in cytoplasm. However, recent studies have shown that miRNAs have a role in the regulation of genes in the nucleus, where they are abundantly located. We show here that in mouse endothelial cell line (C166), nuclear microRNA miR-466c participates in the regulation of vascular endothelial growth factor a (Vegfa) gene expression in hypoxia. Upregulation of Vegfa expression in response to hypoxia was significantly compromised after removal of miR-466c with CRISPR-Cas9 genomic deletion. We identified a promoter-associated long non-coding RNA on mouse Vegfa promoter and show that miR-466c directly binds to this transcript to modulate Vegfa expression. Collectively, these observations suggest that miR-466c regulates Vegfa gene transcription in the nucleus by targeting the promoter, and expands on our understanding of the role of miRNAs well beyond their canonical role.

Functional consequences of genetic variation in the non-coding human genome are difficult to ascertain despite demonstrated associations to common, complex disease traits. To elucidate properties of functional non-coding SNPs with effects in human endothelial cells (EC), we utilized molecular Quantitative Trait Locus (molQTL) analysis for transcription factor binding, chromatin accessibility, and H3K27 acetylation to nominate a set of likely functional non-coding SNPs. Together with information from genome-wide association studies for vascular disease traits, we tested the ability of 34,344 variants to perturb enhancer function in ECs using the highly multiplexed STARR-seq assay. Of these, 5,592 variants validated, whose enriched attributes included: 1) mutations to TF binding motifs for ETS or AP1 that are regulators of EC state, 2) location in accessible and H3K27ac-marked EC chromatin, and 3) molQTLs associations whereby alleles associate with differences in chromatin accessibility and TF binding across genetically diverse ECs. Next, using pro-inflammatory IL1B as an activator of cell state, we observed robust evidence (>50%) of context-specific SNP effects, underscoring the prevalence of non-coding gene-by-environment (GxE) effects. Lastly, using these cumulative data, we fine-mapped vascular disease loci and highlight evidence suggesting mechanisms by which non-coding SNPs at two loci affect risk for Pulse Pressure/Large Artery Stroke, and Abdominal Aortic Aneurysm through respective effects on transcriptional regulation of POU4F1 and LDAH. Together, we highlight the attributes and context dependence of functional non-coding SNPs, and provide new mechanisms underlying vascular disease risk.

Functional consequences of genetic variation in the noncoding human genome are difficult to ascertain despite demonstrated associations to common, complex disease traits. To elucidate properties of functional noncoding SNPs with effects in human endothelial cells (ECs), we utilized our previous molecular quantitative trait locus (molQTL) analysis for transcription factor binding, chromatin accessibility, and H3K27 acetylation to nominate a set of likely functional noncoding SNPs. Together with information from genome-wide association studies (GWASs) for vascular disease traits, we tested the ability of 34,344 variants to perturb enhancer function in ECs using the highly multiplexed STARR-seq assay. Of these, 5711 variants validated, whose enriched attributes included: (1) mutations to TF binding motifs for ETS or AP-1 that are regulators of the EC state; (2) location in accessible and H3K27ac-marked EC chromatin; and (3) molQTL associations whereby alleles associate with differences in chromatin accessibility and TF binding across genetically diverse ECs. Next, using pro-inflammatory IL1B as an activator of cell state, we observed robust evidence (>50%) of context-specific SNP effects, underscoring the prevalence of noncoding gene-by-environment (GxE) effects. Lastly, using these cumulative data, we fine-mapped vascular disease loci and highlighted evidence suggesting mechanisms by which noncoding SNPs at two loci affect risk for pulse pressure/large artery stroke and abdominal aortic aneurysm through respective effects on transcriptional regulation of POU4F1 and LDAH. Together, we highlight the attributes and context dependence of functional noncoding SNPs and provide new mechanisms underlying vascular disease risk.

Ischemic heart disease is the single most common cause of death worldwide with an annual death rate of over 9 million people. Genome-wide association studies have uncovered over 200 genetic loci underlying the disease, providing a deeper understanding of the causal mechanisms leading to it. However, in order to understand ischemic heart disease at the cellular and molecular level, it is necessary to identify the cell-type-specific circuits enabling dissection of driver variants, genes, and signaling pathways in normal and diseased tissues. Here, we provide the first detailed single-cell dissection of the cell types and disease-associated gene expression changes in the living human heart, using cardiac biopsies collected during open-heart surgery from control, ischemic heart disease, and ischemic and non-ischemic heart failure patients. We identify 84 cell types/states, grouped in 12 major cell types. We define markers for each cell type, providing the first extensive reference set for the live human heart. These major cell types include cardiovascular cells (cardiomyocytes, endothelial cells, fibroblasts), rarer cell types (B lymphocytes, neurons, Schwann cells), and rich populations of previously understudied layer-specific epicardial and endocardial cells. In addition, we reveal substantial differences in disease-associated gene expression at the cell subtype level, revealing t arterial pericytes as having a central role in the pathogenesis of ischemic heart disease and heart failure. Our results demonstrate the importance of high-resolution cellular subtype mapping in gaining mechanistic insight into human cardiovascular disease.

Background: Vascular smooth muscle cells (VSMC) exhibit phenotypic plasticity in atherosclerotic plaques, and among other approaches, has been modeled in vitro by cholesterol loading. Methods: Meta-analysis of scRNA-seq data from VSMC lineage traced cells across five experiments of murine atherosclerosis was performed. In vivo expression profiles were compared to three in vitro datasets of VSMCs loaded with cholesterol and three datasets of polarized macrophages. Results: We identified 24 cell clusters in the meta-analysis of single cells from mouse atherosclerotic lesions with notable heterogeneity across studies, especially for macrophage populations. Trajectory analysis of VSMC lineage positive cells revealed several possible paths of state transitions with one traversing from contractile VSMC to macrophages by way of a proliferative cell cluster. Transcriptome comparisons between in vivo and in vitro states underscored that data from three in vitro cholesterol-treated VSMC experiments did not mirror cell state transitions observed in vivo. However, all in vitro macrophage profiles analyzed (M1, M2, and oxLDL) were more similar to in vivo profiles of macrophages than in vitro VSMCs were to in vivo profiles of VSMCs. oxLDL loaded macrophages showed the most similarity to in vivo states. In contrast to the in vitro data, comparison between mouse and human in vivo data showed many similarities. Conclusions: Identification of the sources of variation across single cell datasets in atherosclerosis will be an important step towards understanding VSMC fate transitions in vivo. Also, we conclude that cholesterol-loading in vitro is insufficient to model the VSMC cell state transitions observed in vivo, which underscores the need to develop better cell models. Mouse models, however, appear to reproduce a number of the features of VSMCs in human plaques

Vascular endothelial growth factors (VEGFs) are best known as key regulators of angiogenesis and lymphangiogenesis. Although VEGFs have been promising therapeutic targets for various cardiovascular diseases, their regulatory landscape in endothelial cells remains elusive. Several studies have highlighted the involvement of noncoding RNAs (ncRNAs) in the modulation of VEGF expression. In this study, we investigated the role of two classes of ncRNAs, long ncRNAs (lncRNAs) and enhancer RNAs (eRNAs), in the transcriptional regulation of VEGFA and VEGFC. By integrating genome-wide global run-on sequencing (GRO-Seq) and chromosome conformation capture (Hi-C) data, we identified putative lncRNAs and eRNAs associated with VEGFA and VEGFC genes in endothelial cells. A subset of the identified putative enhancers demonstrated regulatory activity in a reporter assay. Importantly, we demonstrate that deletion of enhancers and lncRNAs by CRISPR/Cas9 promoted significant changes in VEGFA and VEGFC expression. Transcriptome sequencing (RNA-Seq) data from lncRNA deletions showed downstream factors implicated in VEGFA- and VEGFC-linked pathways, such as angiogenesis and lymphangiogenesis, suggesting functional roles for these lncRNAs. Our study uncovers novel lncRNAs and eRNAs regulating VEGFA and VEGFC that can be targeted to modulate the expression of these important molecules in endothelial cells.

The global obesity epidemic is driving the concomitant rise in nonalcoholic fatty liver disease (NAFLD). To identify new genes involved in central liver functions, we examined liver RNA‐sequence data from 259 patients who underwent morbidly obese bariatric surgery. Of these patients, 84 had normal liver histology, 40 simple steatosis, 43 nonalcoholic steatohepatitis, and the remaining 92 patients had varying degrees of NAFLD based on liver histology. We discovered oligodendrocyte maturation‐associated long intergenic noncoding RNA (OLMALINC), a long intervening noncoding RNA (lincRNA) in a human liver co‐expression network (n = 75 genes) that was strongly associated with statin use and serum triglycerides (TGs). OLMALINC liver expression was highly correlated with the expression of known cholesterol biosynthesis genes and stearoyl‐coenzyme A desaturase (SCD). SCD is the rate‐limiting enzyme in monounsaturated fatty acids and a key TG gene that is known to be up‐regulated in liver steatosis and NAFLD and resides adjacent to OLMALINC on the human chromosome 10q24.31. Next, we functionally demonstrated that OLMALINC regulates SCD as an enhancer‐RNA (eRNA), thus describing the first lincRNA that functions as an eRNA to regulate lipid metabolism. Specifically, we show that OLMALINC promotes liver expression of SCD in cis through regional chromosomal DNA–DNA looping interactions. Conclusion: The primate‐specific lincRNA OLMALINC is a novel epigenetic regulator of the key TG and NAFLD gene SCD.

Statins are effective drugs used to prevent and treat cardiovascular diseases but their effects in the absence of low density lipoprotein receptor (LDLR) and on the risk of diabetes are not yet well characterized. The aim of this study was to clarify systemic and pleiotropic effects of rosuvastatin on cardiovascular and diabetic phenotypes. IGF-II/LDLR-/-ApoB100/100 hypercholesterolemic prediabetic mice were used to test the effects of rosuvastatin on plasma glucose, insulin, lipids, atherosclerosis and liver steatosis. To get a more comprehensive view about changes in gene expression RNA-sequencing was done from the liver. Rosuvastatin significantly reduced plasma cholesterol in hypercholesterolemic mice in the absence of LDLR but had no effects on atherosclerosis at aortic sinus level or in coronary arteries. Rosuvastatin also significantly reduced liver steatosis without any harmful effects on glucose or insulin metabolism. RNA-sequencing showed relatively specific effects of rosuvastatin on genes involved in cholesterol metabolism together with a significant anti-inflammatory gene expression profile in the liver. In addition, significant changes were found in the expression of Perilipin 4 and 5 which are involved in lipid droplet formation in the liver. For the first time it could be shown that Tribbles proteins are affected by rosuvastatin treatment in the hyperlipidemic mice. Rosuvastatin had several positive effects on hypercholesterolemic mice showing early signs of diabetes, many of which are unrelated to cholesterol and lipoprotein metabolism. These results increase our understanding about the systemic and pleiotropic effects of rosuvastatin in the absence of LDLR expression.

Simple Summary Proteasome inhibitors are currently used in the treatment of certain blood cancers, and clinical trials to treat solid tumors, including liver cancer, have also been conducted. However, different malignancies are not equally susceptible to proteasome inhibitors, and resistance to the drug may develop during the therapy. Here, we characterize the molecular mechanisms underlying the resilience of liver cancer cells to the proteasome inhibitor bortezomib. The results demonstrate that the activity of the eIF2α–ATF4 stress response pathway affects the viability of cells treated with bortezomib. We found that the pseudokinase TRIB3, an endogenous regulator of ATF4 and a gene highly expressed in liver cancer, resides predominantly at the same chromatin sites as ATF4 and constrains ATF4 activity. The survival of bortezomib-exposed hepatoma cells proved sensitive to TRIB3 overexpression and inactivation. Thus, TRIB3 is a novel factor contributing to bortezomib resistance of liver cancer cells. Abstract The proteasome is an appealing target for anticancer therapy and the proteasome inhibitor bortezomib has been approved for the treatment of several types of malignancies. However, the molecular mechanisms underlying cancer cell resistance to bortezomib remain poorly understood. In the current article, we investigate how modulation of the eIF2α–ATF4 stress pathway affects hepatoma cell response to bortezomib. Transcriptome profiling revealed that many ATF4 transcriptional target genes are among the most upregulated genes in bortezomib-treated HepG2 human hepatoma cells. While pharmacological enhancement of the eIF2α–ATF4 pathway activity results in the elevation of the activities of all branches of the unfolded protein response (UPR) and sensitizes cells to bortezomib toxicity, the suppression of ATF4 induction delays bortezomib-induced cell death. The pseudokinase TRIB3, an inhibitor of ATF4, is expressed at a high basal level in hepatoma cells and is strongly upregulated in response to bortezomib. To map genome-wide chromatin binding loci of TRIB3 protein, we fused a Flag tag to endogenous TRIB3 in HepG2 cells and performed ChIP-Seq. The results demonstrate that TRIB3 predominantly colocalizes with ATF4 on chromatin and binds to genomic regions containing the C/EBP–ATF motif. Bortezomib treatment leads to a robust enrichment of TRIB3 binding near genes induced by bortezomib and involved in the ER stress response and cell death. Disruption of TRIB3 increases C/EBP–ATF-driven transcription, augments ER stress and cell death upon exposure to bortezomib, while TRIB3 overexpression enhances cell survival. Thus, TRIB3, colocalizing with ATF4 and limiting its transcriptional activity, functions as a factor increasing resistance to bortezomib, while pharmacological over-activation of eIF2α–ATF4 can overcome the endogenous restraint mechanisms and sensitize cells to bortezomib.