Comprehensive Integration of Multiple Single-Cell Transcriptomic Data Sets Defines Distinct Cell Populations and Their Phenotypic Changes in Murine Atherosclerosis.

BACKGROUND: The application of single-cell transcriptomic (single-cell RNA sequencing) analysis to the study of atherosclerosis has provided unique insights into the molecular and genetic mechanisms that mediate disease risk and pathophysiology. However, nonstandardized methodologies and relatively high costs associated with the technique have limited the size and replication of existing data sets and created disparate or contradictory findings that have fostered misunderstanding and controversy. METHODS: To address these uncertainties, we have performed a conservative integration of multiple published single-cell RNA sequencing data sets into a single meta-analysis, performed extended analysis of native resident vascular cells, and used in situ hybridization to map the disease anatomic location of the identified cluster cells. To investigate the transdifferentiation of smooth muscle cells to macrophage phenotype, we have developed a classifying algorithm based on the quantification of reporter transgene expression. RESULTS: The reporter gene expression tool indicates that within the experimental limits of the examined studies, transdifferentiation of smooth muscle cell to the macrophage lineage is extremely rare. Validated transition smooth muscle cell phenotypes were defined by clustering, and the location of these cells was mapped to lesion anatomy with in situ hybridization. We have also characterized 5 endothelial cell phenotypes and linked these cellular species to different vascular structures and functions. Finally, we have identified a transcriptomically unique cellular phenotype that constitutes the aortic valve. CONCLUSIONS: Taken together, these analyses resolve a number of outstanding issues related to differing results reported with vascular disease single-cell RNA sequencing studies, and significantly extend our understanding of the role of resident vascular cells in anatomy and disease.

Early clinical outcomes and molecular smooth muscle cell phenotyping using a prophylactic aortic arch replacement strategy in Loeys-Dietz syndrome.

OBJECTIVES: Patients with Loeys-Dietz syndrome demonstrate a heightened risk of distal thoracic aortic events after valve-sparing aortic root replacement. This study assesses the clinical risks and hemodynamic consequences of a prophylactic aortic arch replacement strategy in Loeys-Dietz syndrome and characterizes smooth muscle cell phenotype in Loeys-Dietz syndrome aneurysmal and normal-sized downstream aorta. METHODS: Patients with genetically confirmed Loeys-Dietz syndrome (n = 8) underwent prophylactic aortic arch replacement during valve-sparing aortic root replacement. Four-dimensional flow magnetic resonance imaging studies were performed in 4 patients with Loeys-Dietz syndrome (valve-sparing aortic root replacement + arch) and compared with patients with contemporary Marfan syndrome (valve-sparing aortic root replacement only, n = 5) and control patients (without aortopathy, n = 5). Aortic tissues from 4 patients with Loeys-Dietz syndrome and 2 organ donors were processed for anatomically segmented single-cell RNA sequencing and histologic assessment. RESULTS: Patients with Loeys-Dietz syndrome valve-sparing aortic root replacement + arch had no deaths, major morbidity, or aortic events in a median of 2 years follow-up. Four-dimensional magnetic resonance imaging demonstrated altered flow parameters in patients with postoperative aortopathy relative to controls, but no clear deleterious changes due to arch replacement. Integrated analysis of aortic single-cell RNA sequencing data (>49,000 cells) identified a continuum of abnormal smooth muscle cell phenotypic modulation in Loeys-Dietz syndrome defined by reduced contractility and enriched extracellular matrix synthesis, adhesion receptors, and transforming growth factor-beta signaling. These modulated smooth muscle cells populated the Loeys-Dietz syndrome tunica media with gradually reduced density from the overtly aneurysmal root to the nondilated arch. CONCLUSIONS: Patients with Loeys-Dietz syndrome demonstrated excellent surgical outcomes without overt downstream flow or shear stress disturbances after concomitant valve-sparing aortic root replacement + arch operations. Abnormal smooth muscle cell-mediated aortic remodeling occurs within the normal diameter, clinically at-risk Loeys-Dietz syndrome arch segment. These initial clinical and pathophysiologic findings support concomitant arch replacement in Loeys-Dietz syndrome.

Single-nucleus chromatin accessibility profiling highlights regulatory mechanisms of coronary artery disease risk.

Coronary artery disease (CAD) is a complex inflammatory disease involving genetic influences across cell types. Genome-wide association studies have identified over 200 loci associated with CAD, where the majority of risk variants reside in noncoding DNA sequences impacting cis-regulatory elements. Here, we applied single-nucleus assay for transposase-accessible chromatin with sequencing to profile 28,316 nuclei across coronary artery segments from 41 patients with varying stages of CAD, which revealed 14 distinct cellular clusters. We mapped ~320,000 accessible sites across all cells, identified cell-type-specific elements and transcription factors, and prioritized functional CAD risk variants. We identified elements in smooth muscle cell transition states (for example, fibromyocytes) and functional variants predicted to alter smooth muscle cell- and macrophage-specific regulation of MRAS (3q22) and LIPA (10q23), respectively. We further nominated key driver transcription factors such as PRDM16 and TBX2. Together, this single-nucleus atlas provides a critical step towards interpreting regulatory mechanisms across the continuum of CAD risk.

ZEB2 Shapes the Epigenetic Landscape of Atherosclerosis.

BACKGROUND: Smooth muscle cells (SMCs) transition into a number of different phenotypes during atherosclerosis, including those that resemble fibroblasts and chondrocytes, and make up the majority of cells in the atherosclerotic plaque. To better understand the epigenetic and transcriptional mechanisms that mediate these cell state changes, and how they relate to risk for coronary artery disease (CAD), we have investigated the causality and function of transcription factors at genome-wide associated loci. METHODS: We used CRISPR-Cas 9 genome and epigenome editing to identify the causal gene and cells for a complex CAD genome-wide association study signal at 2q22.3. Single-cell epigenetic and transcriptomic profiling in murine models and human coronary artery smooth muscle cells were used to understand the cellular and molecular mechanism by which this CAD risk gene exerts its function. RESULTS: CRISPR-Cas 9 genome and epigenome editing showed that the complex CAD genetic signals within a genomic region at 2q22.3 lie within smooth muscle long-distance enhancers for ZEB2, a transcription factor extensively studied in the context of epithelial mesenchymal transition in development of cancer. Zeb2 regulates SMC phenotypic transition through chromatin remodeling that obviates accessibility and disrupts both Notch and transforming growth factor ß signaling, thus altering the epigenetic trajectory of SMC transitions. SMC-specific loss of Zeb2 resulted in an inability of transitioning SMCs to turn off contractile programing and take on a fibroblast-like phenotype, but accelerated the formation of chondromyocytes, mirroring features of high-risk atherosclerotic plaques in human coronary arteries. CONCLUSIONS: These studies identify ZEB2 as a new CAD genome-wide association study gene that affects features of plaque vulnerability through direct effects on the epigenome, providing a new therapeutic approach to target vascular disease.

Genomic integrity of human induced pluripotent stem cells across nine studies in the NHLBI NextGen program.

Human induced pluripotent stem cell (hiPSC) lines have previously been generated through the NHLBI sponsored NextGen program at nine individual study sites. Here, we examined the structural integrity of 506 hiPSC lines as determined by copy number variations (CNVs). We observed that 149 hiPSC lines acquired 258 CNVs relative to donor DNA. We identified six recurrent regions of CNVs on chromosomes 1, 2, 3, 16 and 20 that overlapped with cancer associated genes. Furthermore, the genes mapping to regions of acquired CNVs show an enrichment in cancer related biological processes (IL6 production) and signaling cascades (JNK cascade & NFκB cascade). The genomic region of instability on chr20 (chr20q11.2) includes transcriptomic signatures for cancer associated genes such as ID1, BCL2L1, TPX2, PDRG1 and HCK. Of these HCK shows statistically significant differential expression between carrier and non-carrier hiPSC lines. Overall, while a low level of genomic instability was observed in the NextGen generated hiPSC lines, the observation of structural instability in regions with known cancer associated genes substantiates the importance of systematic evaluation of genetic variations in hiPSCs before using them as disease/research models.

PCSK6 Is a Key Protease in the Control of Smooth Muscle Cell Function in Vascular Remodeling.

RATIONALE: PCSKs (Proprotein convertase subtilisins/kexins) are a protease family with unknown functions in vasculature. Previously, we demonstrated PCSK6 upregulation in human atherosclerotic plaques associated with smooth muscle cells (SMCs), inflammation, extracellular matrix remodeling, and mitogens. OBJECTIVE: Here, we applied a systems biology approach to gain deeper insights into the PCSK6 role in normal and diseased vessel wall. METHODS AND RESULTS: Genetic analyses revealed association of intronic PCSK6 variant rs1531817 with maximum internal carotid intima-media thickness progression in high-cardiovascular risk subjects. This variant was linked with PCSK6 mRNA expression in healthy aortas and plaques but also with overall plaque SMA+ cell content and pericyte fraction. Increased PCSK6 expression was found in several independent human cohorts comparing atherosclerotic lesions versus healthy arteries, using transcriptomic and proteomic datasets. By immunohistochemistry, PCSK6 was localized to fibrous cap SMA+ cells and neovessels in plaques. In human, rat, and mouse intimal hyperplasia, PCSK6 was expressed by proliferating SMA+ cells and upregulated after 5 days in rat carotid balloon injury model, with positive correlation to PDGFB (platelet-derived growth factor subunit B) and MMP (matrix metalloprotease) 2/MMP14. Here, PCSK6 was shown to colocalize and cointeract with MMP2/MMP14 by in situ proximity ligation assay. Microarrays of carotid arteries from Pcsk6-/- versus control mice revealed suppression of contractile SMC markers, extracellular matrix remodeling enzymes, and cytokines/receptors. Pcsk6-/- mice showed reduced intimal hyperplasia response upon carotid ligation in vivo, accompanied by decreased MMP14 activation and impaired SMC outgrowth from aortic rings ex vivo. PCSK6 silencing in human SMCs in vitro leads to downregulation of contractile markers and increase in MMP2 expression. Conversely, PCSK6 overexpression increased PDGFBB (platelet-derived growth factor BB)-induced cell proliferation and particularly migration. CONCLUSIONS: PCSK6 is a novel protease that induces SMC migration in response to PDGFB, mechanistically via modulation of contractile markers and MMP14 activation. This study establishes PCSK6 as a key regulator of SMC function in vascular remodeling. Visual Overview: An online visual overview is available for this article.

Transcriptomic profiling of experimental arterial injury reveals new mechanisms and temporal dynamics in vascular healing response.

OBJECTIVE: Endovascular interventions cause arterial injury and induce a healing response to restore vessel wall homeostasis. Complications of defective or excessive healing are common and result in increased morbidity and repeated interventions. Experimental models of intimal hyperplasia are vital for understanding the vascular healing mechanisms and resolving the clinical problems of restenosis, vein graft stenosis, and dialysis access failure. Our aim was to systematically investigate the transcriptional, histologic, and systemic reaction to vascular injury during a prolonged time. METHODS: Balloon injury of the left common carotid artery was performed in male rats. Animals (n = 69) were euthanized before or after injury, either directly or after 2 hours, 20 hours, 2 days, 5 days, 2 weeks, 6 weeks, and 12 weeks. Both injured and contralateral arteries were subjected to microarray profiling, followed by bioinformatic exploration, histologic characterization of the biopsy specimens, and plasma lipid analyses. RESULTS: Immune activation and coagulation were key mechanisms in the early response, followed by cytokine release, tissue remodeling, and smooth muscle cell modulation several days after injury, with reacquisition of contractile features in later phases. Novel pathways related to clonal expansion, inflammatory transformation, and chondro-osteogenic differentiation were identified and immunolocalized to neointimal smooth muscle cells. Analysis of uninjured arteries revealed a systemic component of the reaction after local injury, underlined by altered endothelial signaling, changes in overall tissue bioenergy metabolism, and plasma high-density lipoprotein levels. CONCLUSIONS: We demonstrate that vascular injury induces dynamic transcriptional landscape and metabolic changes identifiable as early, intermediate, and late response phases, reaching homeostasis after several weeks. This study provides a temporal "roadmap" of vascular healing as a publicly available resource for the research community.

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.

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.

Opportunities and challenges for transcriptome-wide association studies.

Transcriptome-wide association studies (TWAS) integrate genome-wide association studies (GWAS) and gene expression datasets to identify gene-trait associations. In this Perspective, we explore properties of TWAS as a potential approach to prioritize causal genes at GWAS loci, by using simulations and case studies of literature-curated candidate causal genes for schizophrenia, low-density-lipoprotein cholesterol and Crohn's disease. We explore risk loci where TWAS accurately prioritizes the likely causal gene as well as loci where TWAS prioritizes multiple genes, some likely to be non-causal, owing to sharing of expression quantitative trait loci (eQTL). TWAS is especially prone to spurious prioritization with expression data from non-trait-related tissues or cell types, owing to substantial cross-cell-type variation in expression levels and eQTL strengths. Nonetheless, TWAS prioritizes candidate causal genes more accurately than simple baselines. We suggest best practices for causal-gene prioritization with TWAS and discuss future opportunities for improvement. Our results showcase the strengths and limitations of using eQTL datasets to determine causal genes at GWAS loci.

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.

Large-Scale Single-Cell RNA-Seq Reveals Molecular Signatures of Heterogeneous Populations of Human Induced Pluripotent Stem Cell-Derived Endothelial Cells.

RATIONALE: Human-induced pluripotent stem cell-derived endothelial cells (iPSC-ECs) have risen as a useful tool in cardiovascular research, offering a wide gamut of translational and clinical applications. However, inefficiency of the currently available iPSC-EC differentiation protocol and underlying heterogeneity of derived iPSC-ECs remain as major limitations of iPSC-EC technology. OBJECTIVE: Here, we performed droplet-based single-cell RNA sequencing (scRNA-seq) of the human iPSCs after iPSC-EC differentiation. Droplet-based scRNA-seq enables analysis of thousands of cells in parallel, allowing comprehensive analysis of transcriptional heterogeneity. METHODS AND RESULTS: Bona fide iPSC-EC cluster was identified by scRNA-seq, which expressed high levels of endothelial-specific genes. iPSC-ECs, sorted by CD144 antibody-conjugated magnetic sorting, exhibited standard endothelial morphology and function including tube formation, response to inflammatory signals, and production of NO. Nonendothelial cell populations resulting from the differentiation protocol were identified, which included immature cardiomyocytes, hepatic-like cells, and vascular smooth muscle cells. Furthermore, scRNA-seq analysis of purified iPSC-ECs revealed transcriptional heterogeneity with 4 major subpopulations, marked by robust enrichment of CLDN5, APLNR, GJA5, and ESM1 genes, respectively. CONCLUSIONS: Massively parallel, droplet-based scRNA-seq allowed meticulous analysis of thousands of human iPSCs subjected to iPSC-EC differentiation. Results showed inefficiency of the differentiation technique, which can be improved with further studies based on identification of molecular signatures that inhibit expansion of nonendothelial cell types. Subtypes of bona fide human iPSC-ECs were also identified, allowing us to sort for iPSC-ECs with specific biological function and identity.

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.

Alternative Progenitor Cells Compensate to Rebuild the Coronary Vasculature in Elabela- and Apj-Deficient Hearts.

Organogenesis during embryonic development occurs through the differentiation of progenitor cells. This process is extraordinarily accurate, but the mechanisms ensuring high fidelity are poorly understood. Coronary vessels of the mouse heart derive from at least two progenitor pools, the sinus venosus and endocardium. We find that the ELABELA (ELA)-APJ signaling axis is only required for sinus venosus-derived progenitors. Because they do not depend on ELA-APJ, endocardial progenitors are able to expand and compensate for faulty sinus venosus development in Apj mutants, leading to normal adult heart function. An upregulation of endocardial SOX17 accompanied compensation in Apj mutants, which was also seen in Ccbe1 knockouts, indicating that the endocardium is activated in multiple cases where sinus venosus angiogenesis is stunted. Our data demonstrate that by diversifying their responsivity to growth cues, distinct coronary progenitor pools are able to compensate for each other during coronary development, thereby providing robustness to organ development.

Endothelial APLNR regulates tissue fatty acid uptake and is essential for apelin's glucose-lowering effects.

Treatment of type 2 diabetes mellitus continues to pose an important clinical challenge, with most existing therapies lacking demonstrable ability to improve cardiovascular outcomes. The atheroprotective peptide apelin (APLN) enhances glucose utilization and improves insulin sensitivity. However, the mechanism of these effects remains poorly defined. We demonstrate that the expression of APLNR (APJ/AGTRL1), the only known receptor for apelin, is predominantly restricted to the endothelial cells (ECs) of multiple adult metabolic organs, including skeletal muscle and adipose tissue. Conditional endothelial-specific deletion of Aplnr (AplnrECKO ) resulted in markedly impaired glucose utilization and abrogation of apelin-induced glucose lowering. Furthermore, we identified inactivation of Forkhead box protein O1 (FOXO1) and inhibition of endothelial expression of fatty acid (FA) binding protein 4 (FABP4) as key downstream signaling targets of apelin/APLNR signaling. Both the Apln-/- and AplnrECKO mice demonstrated increased endothelial FABP4 expression and excess tissue FA accumulation, whereas concurrent endothelial Foxo1 deletion or pharmacologic FABP4 inhibition rescued the excess FA accumulation phenotype of the Apln-/- mice. The impaired glucose utilization in the AplnrECKO mice was associated with excess FA accumulation in the skeletal muscle. Treatment of these mice with an FABP4 inhibitor abrogated these metabolic phenotypes. These findings provide mechanistic insights that could greatly expand the therapeutic repertoire for type 2 diabetes and related metabolic disorders.

Induced Pluripotent Stem Cell-Derived Endothelial Cells in Insulin Resistance and Metabolic Syndrome.

Insulin resistance leads to a number of metabolic and cellular abnormalities including endothelial dysfunction that increase the risk of vascular disease. Although it has been particularly challenging to study the genetic determinants that predispose to abnormal function of the endothelium in insulin-resistant states, the possibility of deriving endothelial cells from induced pluripotent stem cells generated from individuals with detailed clinical phenotyping, including accurate measurements of insulin resistance accompanied by multilevel omic data (eg, genetic and genomic characterization), has opened new avenues to study this relationship. Unfortunately, several technical barriers have hampered these efforts. In the present review, we summarize the current status of induced pluripotent stem cell-derived endothelial cells for modeling endothelial dysfunction associated with insulin resistance and discuss the challenges to overcoming these limitations.

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.

Nat1 Deficiency Is Associated with Mitochondrial Dysfunction and Exercise Intolerance in Mice.

We recently identified human N-acetyltransferase 2 (NAT2) as an insulin resistance (IR) gene. Here, we examine the cellular mechanism linking NAT2 to IR and find that Nat1 (mouse ortholog of NAT2) is co-regulated with key mitochondrial genes. RNAi-mediated silencing of Nat1 led to mitochondrial dysfunction characterized by increased intracellular reactive oxygen species and mitochondrial fragmentation as well as decreased mitochondrial membrane potential, biogenesis, mass, cellular respiration, and ATP generation. These effects were consistent in 3T3-L1 adipocytes, C2C12 myoblasts, and in tissues from Nat1-deficient mice, including white adipose tissue, heart, and skeletal muscle. Nat1-deficient mice had changes in plasma metabolites and lipids consistent with a decreased ability to utilize fats for energy and a decrease in basal metabolic rate and exercise capacity without altered thermogenesis. Collectively, our results suggest that Nat1 deficiency results in mitochondrial dysfunction, which may constitute a mechanistic link between this gene and IR.

CD47-blocking antibodies restore phagocytosis and prevent atherosclerosis.

Atherosclerosis is the disease process that underlies heart attack and stroke. Advanced lesions at risk of rupture are characterized by the pathological accumulation of diseased vascular cells and apoptotic cellular debris. Why these cells are not cleared remains unknown. Here we show that atherogenesis is associated with upregulation of CD47, a key anti-phagocytic molecule that is known to render malignant cells resistant to programmed cell removal, or 'efferocytosis'. We find that administration of CD47-blocking antibodies reverses this defect in efferocytosis, normalizes the clearance of diseased vascular tissue, and ameliorates atherosclerosis in multiple mouse models. Mechanistic studies implicate the pro-atherosclerotic factor TNF-α as a fundamental driver of impaired programmed cell removal, explaining why this process is compromised in vascular disease. Similar to recent observations in cancer, impaired efferocytosis appears to play a pathogenic role in cardiovascular disease, but is not a fixed defect and may represent a novel therapeutic target.