Molecular mechanisms of coronary disease revealed using quantitative trait loci for TCF21 binding, chromatin accessibility, and chromosomal looping.

BACKGROUND: To investigate the epigenetic and transcriptional mechanisms of coronary artery disease (CAD) risk, as well as the functional regulation of chromatin structure and function, we create a catalog of genetic variants associated with three stages of transcriptional cis-regulation in primary human coronary artery vascular smooth muscle cells (HCASMCs). RESULTS: We use a pooling approach with HCASMC lines to map regulatory variants that mediate binding of the CAD-associated transcription factor TCF21 with ChIPseq studies (bQTLs), variants that regulate chromatin accessibility with ATACseq studies (caQTLs), and chromosomal looping with Hi-C methods (clQTLs). We examine the overlap of these QTLs and their relationship to smooth muscle-specific genes and transcription factors. Further, we use multiple analyses to show that these QTLs are highly associated with CAD GWAS loci and correlate to lead SNPs where they show allelic effects. By utilizing genome editing, we verify that identified functional variants can regulate both chromatin accessibility and chromosomal looping, providing new insights into functional mechanisms regulating chromatin state and chromosomal structure. Finally, we directly link the disease-associated TGFB1-SMAD3 pathway to the CAD-associated FN1 gene through a response QTL that modulates both chromatin accessibility and chromosomal looping. CONCLUSIONS: Together, these studies represent the most thorough mapping of multiple QTL types in a highly disease-relevant primary cultured cell type and provide novel insights into their functional overlap and mechanisms that underlie these genomic features and their relationship to disease risk.

Environment-Sensing Aryl Hydrocarbon Receptor Inhibits the Chondrogenic Fate of Modulated Smooth Muscle Cells in Atherosclerotic Lesions.

BACKGROUND: Smooth muscle cells (SMC) play a critical role in atherosclerosis. The Aryl hydrocarbon receptor (AHR) is an environment-sensing transcription factor that contributes to vascular development, and has been implicated in coronary artery disease risk. We hypothesized that AHR can affect atherosclerosis by regulating phenotypic modulation of SMC. METHODS: We combined RNA-sequencing, chromatin immunoprecipitation followed by sequencing, assay for transposase-accessible chromatin using sequencing, and in vitro assays in human coronary artery SMCs, with single-cell RNA-sequencing, histology, and RNAscope in an SMC-specific lineage-tracing Ahr knockout mouse model of atherosclerosis to better understand the role of AHR in vascular disease. RESULTS: Genomic studies coupled with functional assays in cultured human coronary artery SMCs revealed that AHR modulates the human coronary artery SMC phenotype and suppresses ossification in these cells. Lineage-tracing and activity-tracing studies in the mouse aortic sinus showed that the Ahr pathway is active in modulated SMCs in the atherosclerotic lesion cap. Furthermore, single-cell RNA-sequencing studies of the SMC-specific Ahr knockout mice showed a significant increase in the proportion of modulated SMCs expressing chondrocyte markers such as Col2a1 and Alpl, which localized to the lesion neointima. These cells, which we term "chondromyocytes," were also identified in the neointima of human coronary arteries. In histological analyses, these changes manifested as larger lesion size, increased lineage-traced SMC participation in the lesion, decreased lineage-traced SMCs in the lesion cap, and increased alkaline phosphatase activity in lesions in the Ahr knockout in comparison with wild-type mice. We propose that AHR is likely protective based on these data and inference from human genetic analyses. CONCLUSIONS: Overall, we conclude that AHR promotes the maintenance of lesion cap integrity and diminishes the disease-related SMC-to-chondromyocyte transition in atherosclerotic tissues.

Genomic profiling of human vascular cells identifies TWIST1 as a causal gene for common vascular diseases.

Genome-wide association studies have identified multiple novel genomic loci associated with vascular diseases. Many of these loci are common non-coding variants that affect the expression of disease-relevant genes within coronary vascular cells. To identify such genes on a genome-wide level, we performed deep transcriptomic analysis of genotyped primary human coronary artery smooth muscle cells (HCASMCs) and coronary endothelial cells (HCAECs) from the same subjects, including splicing Quantitative Trait Loci (sQTL), allele-specific expression (ASE), and colocalization analyses. We identified sQTLs for TARS2, YAP1, CFDP1, and STAT6 in HCASMCs and HCAECs, and 233 ASE genes, a subset of which are also GTEx eGenes in arterial tissues. Colocalization of GWAS association signals for coronary artery disease (CAD), migraine, stroke and abdominal aortic aneurysm with GTEx eGenes in aorta, coronary artery and tibial artery discovered novel candidate risk genes for these diseases. At the CAD and stroke locus tagged by rs2107595 we demonstrate colocalization with expression of the proximal gene TWIST1. We show that disrupting the rs2107595 locus alters TWIST1 expression and that the risk allele has increased binding of the NOTCH signaling protein RBPJ. Finally, we provide data that TWIST1 expression influences vascular SMC phenotypes, including proliferation and calcification, as a potential mechanism supporting a role for TWIST1 in CAD.

Coronary Disease-Associated Gene TCF21 Inhibits Smooth Muscle Cell Differentiation by Blocking the Myocardin-Serum Response Factor Pathway.

RATIONALE: The gene encoding TCF21 (transcription factor 21) has been linked to coronary artery disease risk by human genome-wide association studies in multiple racial ethnic groups. In murine models, Tcf21 is required for phenotypic modulation of smooth muscle cells (SMCs) in atherosclerotic tissues and promotes a fibroblast phenotype in these cells. In humans, TCF21 expression inhibits risk for coronary artery disease. The molecular mechanism by which TCF21 regulates SMC phenotype is not known. OBJECTIVE: To better understand how TCF21 affects the SMC phenotype, we sought to investigate the possible mechanisms by which it regulates the lineage determining MYOCD (myocardin)-SRF (serum response factor) pathway. METHODS AND RESULTS: Modulation of TCF21 expression in human coronary artery SMC revealed that TCF21 suppresses a broad range of SMC markers, as well as key SMC transcription factors MYOCD and SRF, at the RNA and protein level. We conducted chromatin immunoprecipitation-sequencing to map SRF-binding sites in human coronary artery SMC, showing that binding is colocalized in the genome with TCF21, including at a novel enhancer in the SRF gene, and at the MYOCD gene promoter. In vitro genome editing indicated that the SRF enhancer CArG box regulates transcription of the SRF gene, and mutation of this conserved motif in the orthologous mouse SRF enhancer revealed decreased SRF expression in aorta and heart tissues. Direct TCF21 binding and transcriptional inhibition at colocalized sites were established by reporter gene transfection assays. Chromatin immunoprecipitation and protein coimmunoprecipitation studies provided evidence that TCF21 blocks MYOCD and SRF association by direct TCF21-MYOCD interaction. CONCLUSIONS: These data indicate that TCF21 antagonizes the MYOCD-SRF pathway through multiple mechanisms, further establishing a role for this coronary artery disease-associated gene in fundamental SMC processes and indicating the importance of smooth muscle response to vascular stress and phenotypic modulation of this cell type in coronary artery disease risk.

CRISPR-Cas9-mediated knockout of SPRY2 in human hepatocytes leads to increased glucose uptake and lipid droplet accumulation.

BACKGROUND: The prevalence of obesity and its comorbidities, including type 2 diabetes mellitus (T2DM), is dramatically increasing throughout the world; however, the underlying aetiology is incompletely understood. Genome-wide association studies (GWAS) have identified hundreds of genec susceptibility loci for obesity and T2DM, although the causal genes and mechanisms are largely unknown. SPRY2 is a candidate gene identified in GWAS of body fat percentage and T2DM, and has recently been linked to insulin production in pancreatic ß-cells. In the present study, we aimed to further understand SPRY2 via functional characterisation in HepG2 cells, an in vitro model of human hepatocytes widely used to investigate T2DM and insulin resistance. METHODS: CRISPR-Cas9 genome editing was used to target SPRY2 in HepG2 cells, and the functional consequences of SPRY2 knockout (KO) and overexpression subsequently assessed using glucose uptake and lipid droplet assays, measurement of protein kinase phosphorylation and RNA sequencing. RESULTS: The major functional consequence of SPRY2 KO was a significant increase in glucose uptake, along with elevated lipid droplet accumulation. These changes were attenuated, but not reversed, in cells overexpressing SPRY2. Phosphorylation of protein kinases across key signalling pathways (including Akt and mitogen activated protein kinases) was not altered after SPRY2 KO. Transcriptome profiling in SPRY2 KO and mock (control) cells revealed a number of differentially expressed genes related to cholesterol biosynthesis, cell cycle regulation and cellular signalling pathways. Phospholipase A2 group IIA (PLA2G2A) mRNA level was subsequently validated as significantly upregulated following SPRY2 KO, highlighting this as a potential mediator downstream of SPRY2. CONCLUSION: These findings suggest a role for SPRY2 in glucose and lipid metabolism in hepatocytes and contribute to clarifying the function of this gene in the context of metabolic diseases.

Detailed Functional Characterization of a Waist-Hip Ratio Locus in 7p15.2 Defines an Enhancer Controlling Adipocyte Differentiation.

We combined CAGE sequencing in human adipocytes during differentiation with data from genome-wide association studies to identify an enhancer in the SNX10 locus on chromosome 7, presumably involved in body fat distribution. Using reporter assays and CRISPR-Cas9 gene editing in human cell lines, we characterized the role of the enhancer in adipogenesis. The enhancer was active during adipogenesis and responded strongly to insulin and isoprenaline. The allele associated with increased waist-hip ratio in human genetic studies was associated with higher enhancer activity. Mutations of the enhancer resulted in less adipocyte differentiation. RNA sequencing of cells with disrupted enhancer showed reduced expression of established adipocyte markers, such as ADIPOQ and LPL, and identified CHI3L1 on chromosome 1 as a potential gene involved in adipocyte differentiation. In conclusion, we identified and characterized an enhancer in the SNX10 locus and outlined its plausible mechanisms of action and downstream targets.

Atheroprotective roles of smooth muscle cell phenotypic modulation and the TCF21 disease gene as revealed by single-cell analysis.

In response to various stimuli, vascular smooth muscle cells (SMCs) can de-differentiate, proliferate and migrate in a process known as phenotypic modulation. However, the phenotype of modulated SMCs in vivo during atherosclerosis and the influence of this process on coronary artery disease (CAD) risk have not been clearly established. Using single-cell RNA sequencing, we comprehensively characterized the transcriptomic phenotype of modulated SMCs in vivo in atherosclerotic lesions of both mouse and human arteries and found that these cells transform into unique fibroblast-like cells, termed 'fibromyocytes', rather than into a classical macrophage phenotype. SMC-specific knockout of TCF21-a causal CAD gene-markedly inhibited SMC phenotypic modulation in mice, leading to the presence of fewer fibromyocytes within lesions as well as within the protective fibrous cap of the lesions. Moreover, TCF21 expression was strongly associated with SMC phenotypic modulation in diseased human coronary arteries, and higher levels of TCF21 expression were associated with decreased CAD risk in human CAD-relevant tissues. These results establish a protective role for both TCF21 and SMC phenotypic modulation in this disease.

TCF21 and AP-1 interact through epigenetic modifications to regulate coronary artery disease gene expression.

BACKGROUND: Genome-wide association studies have identified over 160 loci that are associated with coronary artery disease. As with other complex human diseases, risk in coronary disease loci is determined primarily by altered expression of the causal gene, due to variation in binding of transcription factors and chromatin-modifying proteins that directly regulate the transcriptional apparatus. We have previously identified a coronary disease network downstream of the disease-associated transcription factor TCF21, and in work reported here extends these studies to investigate the mechanisms by which it interacts with the AP-1 transcription complex to regulate local epigenetic effects in these downstream coronary disease loci. METHODS: Genomic studies, including chromatin immunoprecipitation sequencing, RNA sequencing, and protein-protein interaction studies, were performed in human coronary artery smooth muscle cells. RESULTS: We show here that TCF21 and JUN regulate expression of two presumptive causal coronary disease genes, SMAD3 and CDKN2B-AS1, in part by interactions with histone deacetylases and acetyltransferases. Genome-wide TCF21 and JUN binding is jointly localized and particularly enriched in coronary disease loci where they broadly modulate H3K27Ac and chromatin state changes linked to disease-related processes in vascular cells. Heterozygosity at coronary disease causal variation, or genome editing of these variants, is associated with decreased binding of both JUN and TCF21 and loss of expression in cis, supporting a transcriptional mechanism for disease risk. CONCLUSIONS: These data show that the known chromatin remodeling and pioneer functions of AP-1 are a pervasive aspect of epigenetic control of transcription, and thus, the risk in coronary disease-associated loci, and that interaction of AP-1 with TCF21 to control epigenetic features, contributes to the genetic risk in loci where they co-localize.

Functional regulatory mechanism of smooth muscle cell-restricted LMOD1 coronary artery disease locus.

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.

Coronary artery disease genes SMAD3 and TCF21 promote opposing interactive genetic programs that regulate smooth muscle cell differentiation and disease risk.

Although numerous genetic loci have been associated with coronary artery disease (CAD) with genome wide association studies, efforts are needed to identify the causal genes in these loci and link them into fundamental signaling pathways. Recent studies have investigated the disease mechanism of CAD associated gene SMAD3, a central transcription factor (TF) in the TGFß pathway, investigating its role in smooth muscle biology. In vitro studies in human coronary artery smooth muscle cells (HCASMC) revealed that SMAD3 modulates cellular phenotype, promoting expression of differentiation marker genes while inhibiting proliferation. RNA sequencing and chromatin immunoprecipitation sequencing studies in HCASMC identified downstream genes that reside in pathways which mediate vascular development and atherosclerosis processes in this cell type. HCASMC phenotype, and gene expression patterns promoted by SMAD3 were noted to have opposing direction of effect compared to another CAD associated TF, TCF21. At sites of SMAD3 and TCF21 colocalization on DNA, SMAD3 binding was inversely correlated with TCF21 binding, due in part to TCF21 locally blocking chromatin accessibility at the SMAD3 binding site. Further, TCF21 was able to directly inhibit SMAD3 activation of gene expression in transfection reporter gene studies. In contrast to TCF21 which is protective toward CAD, SMAD3 expression in HCASMC was shown to be directly correlated with disease risk. We propose that the pro-differentiation action of SMAD3 inhibits dedifferentiation that is required for HCASMC to expand and stabilize disease plaque as they respond to vascular stresses, counteracting the protective dedifferentiating activity of TCF21 and promoting disease risk.

Genetic Regulatory Mechanisms of Smooth Muscle Cells Map to Coronary Artery Disease Risk Loci.

Coronary artery disease (CAD) is the leading cause of death globally. Genome-wide association studies (GWASs) have identified more than 95 independent loci that influence CAD risk, most of which reside in non-coding regions of the genome. To interpret these loci, we generated transcriptome and whole-genome datasets using human coronary artery smooth muscle cells (HCASMCs) from 52 unrelated donors, as well as epigenomic datasets using ATAC-seq on a subset of 8 donors. Through systematic comparison with publicly available datasets from GTEx and ENCODE projects, we identified transcriptomic, epigenetic, and genetic regulatory mechanisms specific to HCASMCs. We assessed the relevance of HCASMCs to CAD risk using transcriptomic and epigenomic level analyses. By jointly modeling eQTL and GWAS datasets, we identified five genes (SIPA1, TCF21, SMAD3, FES, and PDGFRA) that may modulate CAD risk through HCASMCs, all of which have relevant functional roles in vascular remodeling. Comparison with GTEx data suggests that SIPA1 and PDGFRA influence CAD risk predominantly through HCASMCs, while other annotated genes may have multiple cell and tissue targets. Together, these results provide tissue-specific and mechanistic insights into the regulation of a critical vascular cell type associated with CAD in human populations.

Advances in Transcriptomics: Investigating Cardiovascular Disease at Unprecedented Resolution.

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

The role of polycarbonate monomer bisphenol-A in insulin resistance.

Bisphenol A (BPA) is a synthetic unit of polycarbonate polymers and epoxy resins, the types of plastics that could be found in essentially every human population and incorporated into almost every aspect of the modern human society. BPA polymers appear in a wide range of products, from liquid storages (plastic bottles, can and glass linings, water pipes and tanks) and food storages (plastics wraps and containers), to medical and dental devices. BPA polymers could be hydrolyzed spontaneously or in a photo- or temperature-catalyzed process, providing widespread environmental distribution and chronic exposure to the BPA monomer in contemporary human populations. Bisphenol A is also a xenoestrogen, an endocrine-disrupting chemical (EDC) that interferes with the endocrine system mimicking the effects of an estrogen and could potentially keep our endocrine system in a constant perturbation that parallels endocrine disruption arising during pregnancy, such as insulin resistance (IR). Gestational insulin resistance represents a natural biological phenomenon of higher insulin resistance in peripheral tissues of the pregnant females, when nutrients are increasingly being directed to the embryo instead of being stored in peripheral tissues. Gestational diabetes mellitus may appear in healthy non-diabetic females, due to gestational insulin resistance that leads to increased blood sugar levels and hyperinsulinemia (increased insulin production from the pancreatic beta cells). The hypothesis states that unnoticed and constant exposure to this environmental chemical might potentially lead to the formation of chronic low-level endocrine disruptive state that resembles gestational insulin resistance, which might contribute to the development of diabetes. The increasing body of evidence supports the major premises of this hypothesis, as exemplified by the numerous publications examining the association of BPA and insulin resistance, both epidemiological and mechanistic. However, to what extent BPA might contribute to the development of diabetes in the modern societies still remains unknown. In this review, I discuss the chemical properties of BPA and the sources of BPA contamination found in the environment and in human tissues. I provide an overview of mechanisms for the proposed role of bisphenol A in insulin resistance and diabetes, as well as other related diseases, such as cardiovascular diseases. I describe the transmission of BPA effects to the offspring and postulate that gender related differences might originate from differences in liver enzyme levels, such as UDP-glucuronosyltransferase, which is involved in BPA processing and its elimination from the organism. I discuss the molecular mechanisms of BPA action through nuclear and membrane-bound ER receptors, non-monotonic dose response, epigenetic modifications of the DNA and propose that chronic exposure to weak binders, such as BPA, may mimic the effects of strong binders, such as estrogens.

TCF21 and the environmental sensor aryl-hydrocarbon receptor cooperate to activate a pro-inflammatory gene expression program in coronary artery smooth muscle cells.

Both environmental factors and genetic loci have been associated with coronary artery disease (CAD), however gene-gene and gene-environment interactions that might identify molecular mechanisms of risk are not easily studied by human genetic approaches. We have previously identified the transcription factor TCF21 as the causal CAD gene at 6q23.2 and characterized its downstream transcriptional network that is enriched for CAD GWAS genes. Here we investigate the hypothesis that TCF21 interacts with a downstream target gene, the aryl hydrocarbon receptor (AHR), a ligand-activated transcription factor that mediates the cellular response to environmental contaminants, including dioxin and polycyclic aromatic hydrocarbons (e.g., tobacco smoke). Perturbation of TCF21 expression in human coronary artery smooth muscle cells (HCASMC) revealed that TCF21 promotes expression of AHR, its heterodimerization partner ARNT, and cooperates with these factors to upregulate a number of inflammatory downstream disease related genes including IL1A, MMP1, and CYP1A1. TCF21 was shown to bind in AHR, ARNT and downstream target gene loci, and co-localization was noted for AHR-ARNT and TCF21 binding sites genome-wide in regions of HCASMC open chromatin. These regions of co-localization were found to be enriched for GWAS signals associated with cardio-metabolic as well as chronic inflammatory disease phenotypes. Finally, we show that similar to TCF21, AHR gene expression is increased in atherosclerotic lesions in mice in vivo using laser capture microdissection, and AHR protein is localized in human carotid atherosclerotic lesions where it is associated with protein kinases with a critical role in innate immune response. These data suggest that TCF21 can cooperate with AHR to activate an inflammatory gene expression program that is exacerbated by environmental stimuli, and may contribute to the overall risk for CAD.

Genetics and Genomics of Coronary Artery Disease.

Coronary artery disease (or coronary heart disease), is the leading cause of mortality in many of the developing as well as the developed countries of the world. Cholesterol-enriched plaques in the heart's blood vessels combined with inflammation lead to the lesion expansion, narrowing of blood vessels, reduced blood flow, and may subsequently cause lesion rupture and a heart attack. Even though several environmental risk factors have been established, such as high LDL-cholesterol, diabetes, and high blood pressure, the underlying genetic composition may substantially modify the disease risk; hence, genome composition and gene-environment interactions may be critical for disease progression. Ongoing scientific efforts have seen substantial advancements related to the fields of genetics and genomics, with the major breakthroughs yet to come. As genomics is the most rapidly advancing field in the life sciences, it is important to present a comprehensive overview of current efforts. Here, we present a summary of various genetic and genomics assays and approaches applied to coronary artery disease research.

Integrative functional genomics identifies regulatory mechanisms at coronary artery disease loci.

Coronary artery disease (CAD) is the leading cause of mortality and morbidity, driven by both genetic and environmental risk factors. Meta-analyses of genome-wide association studies have identified >150 loci associated with CAD and myocardial infarction susceptibility in humans. A majority of these variants reside in non-coding regions and are co-inherited with hundreds of candidate regulatory variants, presenting a challenge to elucidate their functions. Herein, we use integrative genomic, epigenomic and transcriptomic profiling of perturbed human coronary artery smooth muscle cells and tissues to begin to identify causal regulatory variation and mechanisms responsible for CAD associations. Using these genome-wide maps, we prioritize 64 candidate variants and perform allele-specific binding and expression analyses at seven top candidate loci: 9p21.3, SMAD3, PDGFD, IL6R, BMP1, CCDC97/TGFB1 and LMOD1. We validate our findings in expression quantitative trait loci cohorts, which together reveal new links between CAD associations and regulatory function in the appropriate disease context.

Coronary Artery Disease Associated Transcription Factor TCF21 Regulates Smooth Muscle Precursor Cells That Contribute to the Fibrous Cap.

Recent genome wide association studies have identified a number of genes that contribute to the risk for coronary heart disease. One such gene, TCF21, encodes a basic-helix-loop-helix transcription factor believed to serve a critical role in the development of epicardial progenitor cells that give rise to coronary artery smooth muscle cells (SMC) and cardiac fibroblasts. Using reporter gene and immunolocalization studies with mouse and human tissues we have found that vascular TCF21 expression in the adult is restricted primarily to adventitial cells associated with coronary arteries and also medial SMC in the proximal aorta of mouse. Genome wide RNA-Seq studies in human coronary artery SMC (HCASMC) with siRNA knockdown found a number of putative TCF21 downstream pathways identified by enrichment of terms related to CAD, including "vascular disease," "disorder of artery," and "occlusion of artery," as well as disease-related cellular functions including "cellular movement" and "cellular growth and proliferation." In vitro studies in HCASMC demonstrated that TCF21 expression promotes proliferation and migration and inhibits SMC lineage marker expression. Detailed in situ expression studies with reporter gene and lineage tracing revealed that vascular wall cells expressing Tcf21 before disease initiation migrate into vascular lesions of ApoE-/- and Ldlr-/- mice. While Tcf21 lineage traced cells are distributed throughout the early lesions, in mature lesions they contribute to the formation of a subcapsular layer of cells, and others become associated with the fibrous cap. The lineage traced fibrous cap cells activate expression of SMC markers and growth factor receptor genes. Taken together, these data suggest that TCF21 may have a role regulating the differentiation state of SMC precursor cells that migrate into vascular lesions and contribute to the fibrous cap and more broadly, in view of the association of this gene with human CAD, provide evidence that these processes may be a mechanism for CAD risk attributable to the vascular wall.

Characterization of TCF21 Downstream Target Regions Identifies a Transcriptional Network Linking Multiple Independent Coronary Artery Disease Loci.

To functionally link coronary artery disease (CAD) causal genes identified by genome wide association studies (GWAS), and to investigate the cellular and molecular mechanisms of atherosclerosis, we have used chromatin immunoprecipitation sequencing (ChIP-Seq) with the CAD associated transcription factor TCF21 in human coronary artery smooth muscle cells (HCASMC). Analysis of identified TCF21 target genes for enrichment of molecular and cellular annotation terms identified processes relevant to CAD pathophysiology, including "growth factor binding," "matrix interaction," and "smooth muscle contraction." We characterized the canonical binding sequence for TCF21 as CAGCTG, identified AP-1 binding sites in TCF21 peaks, and by conducting ChIP-Seq for JUN and JUND in HCASMC confirmed that there is significant overlap between TCF21 and AP-1 binding loci in this cell type. Expression quantitative trait variation mapped to target genes of TCF21 was significantly enriched among variants with low P-values in the GWAS analyses, suggesting a possible functional interaction between TCF21 binding and causal variants in other CAD disease loci. Separate enrichment analyses found over-representation of TCF21 target genes among CAD associated genes, and linkage disequilibrium between TCF21 peak variation and that found in GWAS loci, consistent with the hypothesis that TCF21 may affect disease risk through interaction with other disease associated loci. Interestingly, enrichment for TCF21 target genes was also found among other genome wide association phenotypes, including height and inflammatory bowel disease, suggesting a functional profile important for basic cellular processes in non-vascular tissues. Thus, data and analyses presented here suggest that study of GWAS transcription factors may be a highly useful approach to identifying disease gene interactions and thus pathways that may be relevant to complex disease etiology.

Molecular characterization of a human matrix attachment region epigenetic regulator.

Matrix attachment regions (MAR) generally act as epigenetic regulatory sequences that increase gene expression, and they were proposed to partition chromosomes into loop-forming domains. However, their molecular mode of action remains poorly understood. Here, we assessed the possible contribution of the AT-rich core and adjacent transcription factor binding motifs to the transcription augmenting and anti-silencing effects of human MAR 1-68. Either flanking sequences together with the AT-rich core were required to obtain the full MAR effects. Shortened MAR derivatives retaining full MAR activity were constructed from combinations of the AT-rich sequence and multimerized transcription factor binding motifs, implying that both transcription factors and the AT-rich microsatellite sequence are required to mediate the MAR effect. Genomic analysis indicated that MAR AT-rich cores may be depleted of histones and enriched in RNA polymerase II, providing a molecular interpretation of their chromatin domain insulator and transcriptional augmentation activities.