Search

Abstract Adipogenesis is a process in which fat-specific progenitor cells (preadipocytes) differentiate into adipocytes that carry out the key metabolic functions of the adipose tissue, including glucose uptake, energy storage, and adipokine secretion. Several cell lines are routinely used to study the molecular regulation of adipogenesis, in particular the immortalized mouse 3T3-L1 line and the primary human Simpson-Golabi-Behmel syndrome (SGBS) line. However, the cell-to-cell variability of transcriptional changes prior to and during adipogenesis in these models is not well understood. Here, we present a single-cell RNA-Sequencing (scRNA-Seq) dataset collected before and during adipogenic differentiation of 3T3-L1 and SGBS cells. To minimize the effects of experimental variation, we mixed 3T3-L1 and SGBS cells and used computational analysis to demultiplex transcriptomes of mouse and human cells. In both models, adipogenesis results in the appearance of three cell clusters, corresponding to preadipocytes, early and mature adipocytes. These data provide a groundwork for comparative studies on these widely used in vitro models of human and mouse adipogenesis, and on cell-to-cell variability during this process.

Summary: Cardiac lymphatic vessels play important roles in fluid homeostasis, inflammation, disease, and regeneration of the heart. The developing cardiac lymphatics in human fetal hearts are closely associated with coronary arteries, similar to those in zebrafish hearts. We identify a population of cardiac lymphatic endothelial cells (LECs) that reside in the epicardium. Single-nuclei multiomic analysis of the human fetal heart reveals the plasticity and heterogeneity of the cardiac endothelium. Furthermore, we find that VEGFC is highly expressed in arterial endothelial cells and epicardium-derived cells, providing a molecular basis for the arterial association of cardiac lymphatic development. Using a cell-type-specific integrative analysis, we identify a population of cardiac lymphatic endothelial cells marked by the PROX1 and the lymphangiocrine RELN and enriched in binding motifs of erythroblast transformation specific (ETS) variant (ETV) transcription factors. We report the in vivo molecular characterization of human cardiac lymphatics and provide a valuable resource to understand fetal heart development.

Genes encode risk for coronary disease, identifying how they function is critical to developing new therapies. In work reported the authors have identified one culprit gene, PDGFD, and studied how it functions to promote disease risk.

Abstract Background Identification of causal genes for polygenic human diseases has been extremely challenging, and our understanding of how physiological and pharmacological stimuli modulate genetic risk at disease-associated loci is limited. Specifically, insulin resistance (IR), a common feature of cardiometabolic disease, including type 2 diabetes, obesity, and dyslipidemia, lacks well-powered genome-wide association studies (GWAS), and therefore, few associated loci and causal genes have been identified. Methods Here, we perform and integrate linkage disequilibrium (LD)-adjusted colocalization analyses across nine cardiometabolic traits (fasting insulin, fasting glucose, insulin sensitivity, insulin sensitivity index, type 2 diabetes, triglycerides, high-density lipoprotein, body mass index, and waist-hip ratio) combined with expression and splicing quantitative trait loci (eQTLs and sQTLs) from five metabolically relevant human tissues (subcutaneous and visceral adipose, skeletal muscle, liver, and pancreas). To elucidate the upstream regulators and functional mechanisms for these genes, we integrate their transcriptional responses to 21 relevant physiological and pharmacological perturbations in human adipocytes, hepatocytes, and skeletal muscle cells and map their protein-protein interactions. Results We identify 470 colocalized loci and prioritize 207 loci with a single colocalized gene. Patterns of shared colocalizations across traits and tissues highlight different potential roles for colocalized genes in cardiometabolic disease and distinguish several genes involved in pancreatic β-cell function from others with a more direct role in skeletal muscle, liver, and adipose tissues. At the loci with a single colocalized gene, 42 of these genes were regulated by insulin and 35 by glucose in perturbation experiments, including 17 regulated by both. Other metabolic perturbations regulated the expression of 30 more genes not regulated by glucose or insulin, pointing to other potential upstream regulators of candidate causal genes. Conclusions Our use of transcriptional responses under metabolic perturbations to contextualize genetic associations from our custom colocalization approach provides a list of likely causal genes and their upstream regulators in the context of IR-associated cardiometabolic risk.

Abstract 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.

Genetic variants in the FAM13A locus have been associated with anthropometric and glycemic traits. Here, using fine-mapping, in vitro knockdown studies in pre-adipocytes and in vivo knockout in mice, the authors show that FAM13A is involved in regulating fat distribution and metabolic traits.

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. : Clinical Relevance: Endovascular intervention causes an injury to the arterial wall that subsequently induces a healing response to restore the vessel wall homeostasis. Complications after vascular interventions related to defective or excessive healing response, such as thrombosis or restenosis, are common and result in increased morbidity, suffering of the patient, need for repeated interventions, and possibly death. Thus, there is a need for better understanding of the underlying molecular mechanisms during vascular injury and healing response to identify and to assess the risk of complications in patients. Using an experimental model of vascular injury, this study demonstrates the full landscape of dynamic transcriptional changes in the resolution of vascular injury, accompanied also by systemic variations in plasma lipid levels and reaching homeostasis several weeks after injury. These results can guide the development of new strategies and molecular targets for modulation of the intimal response on endovascular interventions.

Abstract Rationale Vascular calcification is a prominent feature of late‐stage diabetes, renal and cardiovascular disease (CVD), and has been linked to adverse events. Recent studies in patients reported that plasma levels of osteomodulin (OMD), a proteoglycan involved in bone mineralisation, associate with diabetes and CVD. We hypothesised that OMD could be implicated in these diseases via vascular calcification as a common underlying factor and aimed to investigate its role in this context. Methods and results In patients with chronic kidney disease, plasma OMD levels correlated with markers of inflammation and bone turnover, with the protein present in calcified arterial media. Plasma OMD also associated with cardiac calcification and the protein was detected in calcified valve leaflets by immunohistochemistry. In patients with carotid atherosclerosis, circulating OMD was increased in association with plaque calcification as assessed by computed tomography. Transcriptomic and proteomic data showed that OMD was upregulated in atherosclerotic compared to control arteries, particularly in calcified plaques, where OMD expression correlated positively with markers of smooth muscle cells (SMCs), osteoblasts and glycoproteins. Immunostaining confirmed that OMD was abundantly present in calcified plaques, localised to extracellular matrix and regions rich in α‐SMA+ cells. In vivo, OMD was enriched in SMCs around calcified nodules in aortic media of nephrectomised rats and in plaques from ApoE−/− mice on warfarin. In vitro experiments revealed that OMD mRNA was upregulated in SMCs stimulated with IFNγ, BMP2, TGFβ1, phosphate and β‐glycerophosphate, and by administration of recombinant human OMD protein (rhOMD). Mechanistically, addition of rhOMD repressed the calcification process of SMCs treated with phosphate by maintaining their contractile phenotype along with enriched matrix organisation, thereby attenuating SMC osteoblastic transformation. Mechanistically, the role of OMD is exerted likely through its link with SMAD3 and TGFB1 signalling, and interplay with BMP2 in vascular tissues. Conclusion We report a consistent association of both circulating and tissue OMD levels with cardiovascular calcification, highlighting the potential of OMD as a clinical biomarker. OMD was localised in medial and intimal α‐SMA+ regions of calcified cardiovascular tissues, induced by pro‐inflammatory and pro‐osteogenic stimuli, while the presence of OMD in extracellular environment attenuated SMC calcification.

Summary: 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. : Genetics; Human Genetics; Specialized Functions of Cells; Transcriptomics Subject Areas: Genetics, Human Genetics, Specialized Functions of Cells, Transcriptomics

Abstract 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.

Abstract Background Human-induced pluripotent stem cells (hiPSCs) show a great promise as a renewable source of cells with broad biomedical applications. Since insulin has been used in the maintenance of hiPSCs, in this study we explored the role of insulin in culture of these cells. Methods We report conditions for insulin starvation and stimulation of hiPSCs. Crystal violet staining was used to study the adhesion and proliferation of hiPSCs. Apoptosis and cell cycle assays were performed through flow cytometry. Protein arrays were used to confirm phosphorylation targets, and mRNA sequencing was used to evaluate the effect of transcriptome. Results Insulin improved the seeding and proliferation of hiPSCs. We also observed an altered cell cycle profile and increase in apoptosis in hiPSCs in the absence of insulin. Furthermore, we confirmed phosphorylation of key components of insulin signaling pathway in the presence of insulin and demonstrated the significant effect of insulin on regulation of the mRNA transcriptome of hiPSCs. Conclusion Insulin is a major regulator of seeding, proliferation, phosphorylation and mRNA transcriptome in hiPSCs. Collectively, our work furthers our understanding of human pluripotency and paves the way for future studies that use hiPSCs for modeling genetic ailments affecting insulin signaling pathways.

Abstract 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.

Objectives and Aims: Vascular smooth muscle cells (VSMCs) are key constituents of both normal arteries and atherosclerotic plaques. They have an ability to adapt to changes in the local environment by undergoing phenotypic modulation. An improved understanding of the mechanisms that regulate VSMC phenotypic changes may provide insights that suggest new therapeutic targets in treatment of cardiovascular disease (CVD). The amino-acid glutamate has been associated with CVD risk and VSMCs metabolism in experimental models, and glutamate receptors regulate VSMC biology and promote pulmonary vascular remodeling. However, glutamate-signaling in human atherosclerosis has not been explored.Methods and Results: We identified glutamate receptors and glutamate metabolism-related enzymes in VSMCs from human atherosclerotic lesions, as determined by single cell RNA sequencing and microarray analysis. Expression of the receptor subunits glutamate receptor, ionotropic, α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic (AMPA)-type subunit 1 (GRIA1) and 2 (GRIA2) was restricted to cells of mesenchymal origin, primarily VSMCs, as confirmed by immunostaining. In a rat model of arterial injury and repair, changes of GRIA1 and GRIA2 mRNA level were most pronounced at time points associated with VSMC proliferation, migration, and phenotypic modulation. In vitro, human carotid artery SMCs expressed GRIA1, and selective AMPA-type receptor blocking inhibited expression of typical contractile markers and promoted pathways associated with VSMC phenotypic modulation. In our biobank of human carotid endarterectomies, low expression of AMPA-type receptor subunits was associated with higher content of inflammatory cells and a higher frequency of adverse clinical events such as stroke.Conclusion: AMPA-type glutamate receptors are expressed in VSMCs and are associated with phenotypic modulation. Patients suffering from adverse clinical events showed significantly lower mRNA level of GRIA1 and GRIA2 in their atherosclerotic lesions compared to asymptomatic patients. These results warrant further mapping of neurotransmitter signaling in the pathogenesis of human atherosclerosis.

Insulin resistance (IR) precedes the development of type 2 diabetes (T2D) and increases cardiovascular disease risk. Although genome wide association studies (GWAS) have uncovered new loci associated with T2D, their contribution to explain the mechanisms leading to decreased insulin sensitivity has been very limited. Thus, new approaches are necessary to explore the genetic architecture of insulin resistance. To that end, we generated an iPSC library across the spectrum of insulin sensitivity in humans. RNA-seq based analysis of 310 induced pluripotent stem cell (iPSC) clones derived from 100 individuals allowed us to identify differentially expressed genes between insulin resistant and sensitive iPSC lines. Analysis of the co-expression architecture uncovered several insulin sensitivity-relevant gene sub-networks, and predictive network modeling identified a set of key driver genes that regulate these co-expression modules. Functional validation in human adipocytes and skeletal muscle cells (SKMCs) confirmed the relevance of the key driver candidate genes for insulin responsiveness.

Previous studies have shown an association between elevated atrial NADPH-dependent oxidative stress and decreased plasma apelin in patients with atrial fibrillation (AF), though the basis for this relationship is unclear. In the current study, RT-PCR and immunofluorescence studies of human right atrial appendages (RAAs) showed expression of the apelin receptor, APJ, and reduced apelin content in the atria, but not in plasma, of patients with AF versus normal sinus rhythm. Disruption of the apelin gene in mice increased (2.4-fold) NADPH-stimulated superoxide levels and slowed atrial conduction velocities in optical mapping of a Langendorff-perfused isolated heart model, suggesting that apelin levels may influence AF vulnerability. Indeed, in mice with increased AF vulnerability (induced by chronic intense exercise), apelin administration reduced the incidence and duration of induced atrial arrhythmias in association with prolonged atrial refractory periods. Moreover, apelin decreased AF induction in isolated atria from exercised mice while accelerating conduction velocity and increasing action potential durations. At the cellular level, these changes were associated with increased atrial cardiomyocyte sodium currents. These findings support the conclusion that reduced atrial apelin is maladaptive in fibrillating human atrial myocardium and that increasing apelin bioavailability may be a worthwhile therapeutic strategy for treating and preventing AF.

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.

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.

Abstract C-reactive protein (CRP) encoded by CRP gene is a reflection of systemic inflammation. Many studies associated CRP level with diabetes and glucose levels, but the association of CRP gene with these traits is unclear. We conducted a cross-sectional study consisting of 945 siblings from 330 families collected by the Stanford Asian Pacific Program in Hypertension and Insulin Resistance (SAPPHIRe) to investigate associations between CRP polymorphisms, circulating CRP, diabetes, and glucose levels. Five single-nucleotide polymorphisms were analyzed: rs3093059, rs2794521, rs1417938, rs1800947, and rs1205. The generalized estimating equation approach was used to deal with correlated data within families. CRP level was positively correlated with diabetes prevalence and levels of fasting and 2-hour glucose (each P

Abstract Background Fasting glucose and fasting insulin are glycemic traits closely related to diabetes, and understanding the role of genetic factors in these traits can help reveal the etiology of type 2 diabetes. Although single nucleotide polymorphisms (SNPs) in several candidate genes have been found to be associated with fasting glucose and fasting insulin, copy number variations (CNVs), which have been reported to be associated with several complex traits, have not been reported for association with these two traits. We aimed to identify CNVs associated with fasting glucose and fasting insulin. Results We conducted a genome-wide CNV association analysis for fasting plasma glucose (FPG) and fasting plasma insulin (FPI) using a family-based genome-wide association study sample from a Han Chinese population in Taiwan. A family-based CNV association test was developed in this study to identify common CNVs (i.e., CNVs with frequencies ≥ 5%), and a generalized estimating equation approach was used to test the associations between the traits and counts of global rare CNVs (i.e., CNVs with frequencies

Lysyl oxidase-like 2 (LOXL2) is an enzyme that promotes scaffolding of extracellular matrix proteins. Here the authors show that LOXL2 is crucial for pressure-overload induced cardiac fibrosis, and that antibody-mediated inhibition or genetic disruption ofLoxl2in mice shows therapeutic potential for treatment of cardiac fibrosis.

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.

Coronary heart disease is the leading cause of death worldwide with multiple environmental and genetic risk factors. Here the authors integrate genomic, epigenomic and transcriptomic mapping to elucidate causal variation and mechanisms of known genetic associations.

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.

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.

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.

[This corrects the article DOI: 10.1371/journal.pgen.1005378.].

Congenital heart disease (CHD) has a complex genetic etiology, and recent studies suggest that high penetrance de novo mutations may account for only a small fraction of disease. In a multi-institutional cohort surveyed by exome sequencing, combining analysis of 987 individuals (discovery cohort of 59 affected trios and 59 control trios, and a replication cohort of 100 affected singletons and 533 unaffected singletons) we observe variation at novel and known loci related to a specific cardiac malformation the atrioventricular septal defect (AVSD). In a primary analysis, by combining developmental coexpression networks with inheritance modeling, we identify a de novo mutation in the DNA binding domain of NR1D2 (p.R175W). We show that p.R175W changes the transcriptional activity of Nr1d2 using an in vitro transactivation model in HUVEC cells. Finally, we demonstrate previously unrecognized cardiovascular malformations in the Nr1d2tm1-Dgen knockout mouse. In secondary analyses we map genetic variation to protein-interaction networks suggesting a role for two collagen genes in AVSD, which we corroborate by burden testing in a second replication cohort of 100 AVSDs and 533 controls (p = 8.37e-08). Finally, we apply a rare-disease inheritance model to identify variation in genes previously associated with CHD (ZFPM2, NSD1, NOTCH1, VCAN, and MYH6), cardiac malformations in mouse models (ADAM17, CHRD, IFT140, PTPRJ, RYR1 and ATE1), and hypomorphic alleles of genes causing syndromic CHD (EHMT1, SRCAP, BBS2, NOTCH2, and KMT2D) in 14 of 59 trios, greatly exceeding variation in control trios without CHD (p = 9.60e-06). In total, 32% of trios carried at least one putatively disease-associated variant across 19 loci,suggesting that inherited and de novo variation across a heterogeneous group of loci may contribute to disease risk.

By means of a combination of genome-wide and follow-up studies, recent large-scale association studies of populations of European descent have now identified over 46 loci associated with coronary artery disease (CAD). As part of the TAICHI Consortium, we have collected and genotyped 8556 subjects from Taiwan, comprising 5423 controls and 3133 cases with coronary artery disease, for 9087 CAD SNPs using the CardioMetaboChip. We applied penalized logistic regression to ascertain the top SNPs that contribute together to CAD susceptibility in Taiwan. We observed that the 9p21 locus contributes to CAD at the level of genome-wide significance (rs1537372, with the presence of C, the major allele, the effect estimate is -0.216, standard error 0.033, p value 5.8x10-10). In contrast to a previous report, we propose that the 9p21 locus is a single genetic contribution to CAD in Taiwan because: 1) the penalized logistic regression and the follow-up conditional analysis suggested that rs1537372 accounts for all of the CAD association in 9p21, and 2) the high linkage disequilibrium observed for all associated SNPs in 9p21. We also observed evidence for the following loci at a false discovery rate >5%: SH2B3, ADAMTS7, PHACTR1, GGCX, HTRA1, COL4A1, and LARP6-LRRC49. We also took advantage of the fact that penalized methods are an efficient approach to search for gene-by-gene interactions, and observed that two-way interactions between the PHACTR1 and ADAMTS7 loci and between the SH2B3 and COL4A1 loci contribute to CAD risk. Both the similarities and differences between the significance of these loci when compared with significance of loci in studies of populations of European descent underscore the fact that further genetic association of studies in additional populations will provide clues to identify the genetic architecture of CAD across all populations worldwide.

Genome-wide association studies (GWAS) have identified more than 100 genetic variants contributing to BMI, a measure of body size, or waist-to-hip ratio (adjusted for BMI, WHRadjBMI), a measure of body shape. Body size and shape change as people grow older and these changes differ substantially between men and women. To systematically screen for age- and/or sex-specific effects of genetic variants on BMI and WHRadjBMI, we performed meta-analyses of 114 studies (up to 320,485 individuals of European descent) with genome-wide chip and/or Metabochip data by the Genetic Investigation of Anthropometric Traits (GIANT) Consortium. Each study tested the association of up to ~2.8M SNPs with BMI and WHRadjBMI in four strata (men ≤50y, men >50y, women ≤50y, women >50y) and summary statistics were combined in stratum-specific meta-analyses. We then screened for variants that showed age-specific effects (G x AGE), sex-specific effects (G x SEX) or age-specific effects that differed between men and women (G x AGE x SEX). For BMI, we identified 15 loci (11 previously established for main effects, four novel) that showed significant (FDR

High throughput sequencing has facilitated a precipitous drop in the cost of genomic sequencing, prompting predictions of a revolution in medicine via genetic personalization of diagnostic and therapeutic strategies. There are significant barriers to realizing this goal that are related to the difficult task of interpreting personal genetic variation. A comprehensive, widely accessible application for interpretation of whole genome sequence data is needed. Here, we present a series of methods for identification of genetic variants and genotypes with clinical associations, phasing genetic data and using Mendelian inheritance for quality control, and providing predictive genetic information about risk for rare disease phenotypes and response to pharmacological therapy in single individuals and father-mother-child trios. We demonstrate application of these methods for disease and drug response prognostication in whole genome sequence data from twelve unrelated adults, and for disease gene discovery in one father-mother-child trio with apparently simplex congenital ventricular arrhythmia. In doing so we identify clinically actionable inherited disease risk and drug response genotypes in pre-symptomatic individuals. We also nominate a new candidate gene in congenital arrhythmia, ATP2B4, and provide experimental evidence of a regulatory role for variants discovered using this framework.

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.

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.

Growth arrest-specific 6 (Gas6), a vitamin K-dependent protein, has been implicated in systemic inflammation, obesity, and insulin resistance (IR). Data from recent studies suggest that polymorphisms in the Gas6 gene are associated with cardiovascular disorders and type 2 diabetes (T2D). However, the association of Gas6 gene variants with obesity, IR, and T2D development has not been explored.Four common single nucleotide polymorphisms (SNPs) in the Gas6 gene were genotyped in 984 participants from the Stanford Asia-Pacific Program for Hypertension and Insulin Resistance (SAPPHIRe) family cohort. An insulin suppression test was performed to determine IR based on steady-state plasma glucose (SSPG). Associations between IR indices and obesity, and SNP genotypes, based on previously-reported data for this cohort (Phase I), were analyzed. In the present follow-up study (Phase II), the effects of gene variants of Gas6 on the progression to T2D were explored in individuals who were free of T2D in Phase I. The mean follow-up period for Phase II was 5.7 years.The mean age of the study population in Phase I was 49.5 years and 16.7% of individuals developed T2D during follow-up. After adjusting for covariates, three SNPs (rs8191973, rs8197974, and rs7323932) were found to be associated with SSPG levels (p = 0.007, p = 0.03, and p = 0.011, respectively). This association remained significant after multiple testing and showed a significant interaction with physical activity for SNP rs8191973. However, no other significant correlations were observed between Gas6 polymorphisms and other indices of IR or obesity. A specific haplotype, AACG (from rs8191974, rs7323932, rs7331124, and rs8191973), was positively associated with SSPG levels (p = 0.0098). None of the polymorphisms were associated with an increased risk of T2D development.Our results suggest that Gas6 gene variants are associated with IR, although their effects on subsequent progression to T2D were minimal in this prospective Asian cohort.

Both LPL and HL are synthesized in parenchymal cells, are secreted, and bind to endothelial cells. To learn where endothelial lipase (EL) is synthesized in adult animals, the localization of EL in mouse and rat liver was studied by immunohistochemical analysis. Furthermore, to test whether EL could play a role in atherogenesis, the expression of EL in the aorta and liver of apolipoprotein E knockout (EKO) mice was determined. EL in both mouse and rat liver was colocalized with vascular endothelial cells but not with hepatocytes. In contrast, HL was present in both hepatocytes and endothelial cells. By in situ hybridization, EL mRNA was present only in endothelial cells in liver sections. EL was also present at low levels in aorta of normal mice. We fed EKO mice and wild-type mice a variety of diets and determined EL expression in liver and aorta. EKO mice showed significant expression of EL in aorta. EL expression was lower in the liver of EKO mice than in normal mice. Cholesterol feeding decreased EL in liver of both types of mice. In the aorta, EL was higher in EKO than in wild-type mice, and cholesterol feeding had no effect.Together, these data suggest that EL may be upregulated at the site of atherosclerotic lesions and thus could supply lipids to the area.

Endothelial lipase (EL) is a newly described member of the triglyceride lipase gene family. It has a considerable molecular homology with lipoprotein lipase (LPL) (44%) and hepatic lipase (HL) (41%). Unlike LPL and HL, this enzyme is synthesized by endothelial cells and functions at the site where it is synthesized. Furthermore, its tissue distribution is different from that of LPL and HL. As a lipase, EL has primarily phospholipase A1 activity. Animals that overexpress EL showed reduced HDL cholesterol levels. Conversely, animals that are deficient in EL showed a marked elevation in HDL cholesterol levels, suggesting that it plays a physiologic role in HDL metabolism. Unlike LPL and HL, EL is located in the vascular endothelial cells and its expression is highly regulated by cytokines and physical forces, suggesting that it may play a role in the development of atherosclerosis. However, there is only a limited amount of information available about this enzyme.Some of our unpublished data in addition to previously published data support the possibility that the enzyme plays a role in the formation of atherosclerotic lesion.

The majority of the heritability of coronary artery disease (CAD) remains unexplained, despite recent successes of genome-wide association studies (GWAS) in identifying novel susceptibility loci. Integrating functional genomic data from a variety of sources with a large-scale meta-analysis of CAD GWAS may facilitate the identification of novel biological processes and genes involved in CAD, as well as clarify the causal relationships of established processes. Towards this end, we integrated 14 GWAS from the CARDIoGRAM Consortium and two additional GWAS from the Ottawa Heart Institute (25,491 cases and 66,819 controls) with 1) genetics of gene expression studies of CAD-relevant tissues in humans, 2) metabolic and signaling pathways from public databases, and 3) data-driven, tissue-specific gene networks from a multitude of human and mouse experiments. We not only detected CAD-associated gene networks of lipid metabolism, coagulation, immunity, and additional networks with no clear functional annotation, but also revealed key driver genes for each CAD network based on the topology of the gene regulatory networks. In particular, we found a gene network involved in antigen processing to be strongly associated with CAD. The key driver genes of this network included glyoxalase I (GLO1) and peptidylprolyl isomerase I (PPIL1), which we verified as regulatory by siRNA experiments in human aortic endothelial cells. Our results suggest genetic influences on a diverse set of both known and novel biological processes that contribute to CAD risk. The key driver genes for these networks highlight potential novel targets for further mechanistic studies and therapeutic interventions.

Genome-wide association studies (GWAS) have identified chromosomal loci that affect risk of coronary heart disease (CHD) independent of classical risk factors. One such association signal has been identified at 6q23.2 in both Caucasians and East Asians. The lead CHD-associated polymorphism in this region, rs12190287, resides in the 3' untranslated region (3'-UTR) of TCF21, a basic-helix-loop-helix transcription factor, and is predicted to alter the seed binding sequence for miR-224. Allelic imbalance studies in circulating leukocytes and human coronary artery smooth muscle cells (HCASMC) showed significant imbalance of the TCF21 transcript that correlated with genotype at rs12190287, consistent with this variant contributing to allele-specific expression differences. 3' UTR reporter gene transfection studies in HCASMC showed that the disease-associated C allele has reduced expression compared to the protective G allele. Kinetic analyses in vitro revealed faster RNA-RNA complex formation and greater binding of miR-224 with the TCF21 C allelic transcript. In addition, in vitro probing with Pb2+ and RNase T1 revealed structural differences between the TCF21 variants in proximity of the rs12190287 variant, which are predicted to provide greater access to the C allele for miR-224 binding. miR-224 and TCF21 expression levels were anti-correlated in HCASMC, and miR-224 modulates the transcriptional response of TCF21 to transforming growth factor-β (TGF-β) and platelet derived growth factor (PDGF) signaling in an allele-specific manner. Lastly, miR-224 and TCF21 were localized in human coronary artery lesions and anti-correlated during atherosclerosis. Together, these data suggest that miR-224 interaction with the TCF21 transcript contributes to allelic imbalance of this gene, thus partly explaining the genetic risk for coronary heart disease associated at 6q23.2. These studies implicating rs12190287 in the miRNA-dependent regulation of TCF21, in conjunction with previous studies showing that this variant modulates transcriptional regulation through activator protein 1 (AP-1), suggests a unique bimodal level of complexity previously unreported for disease-associated variants.

In this paper we try to define insulin resistance (IR) precisely for a group of Chinese women. Our definition deliberately does not depend upon body mass index (BMI) or age, although in other studies, with particular random effects models quite different from models used here, BMI accounts for a large part of the variability in IR. We accomplish our goal through application of Gauss mixture vector quantization (GMVQ), a technique for clustering that was developed for application to lossy data compression. Defining data come from measurements that play major roles in medical practice. A precise statement of what the data are is in Section 1. Their family structures are described in detail. They concern levels of lipids and the results of an oral glucose tolerance test (OGTT). We apply GMVQ to residuals obtained from regressions of outcomes of an OGTT and lipids on functions of age and BMI that are inferred from the data. A bootstrap procedure developed for our family data supplemented by insights from other approaches leads us to believe that two clusters are appropriate for defining IR precisely. One cluster consists of women who are IR, and the other of women who seem not to be. Genes and other features are used to predict cluster membership. We argue that prediction with "main effects" is not satisfactory, but prediction that includes interactions may be.

Genome-wide association studies (GWAS) have identified ~100 loci associated with blood lipid levels, but much of the trait heritability remains unexplained, and at most loci the identities of the trait-influencing variants remain unknown. We conducted a trans-ethnic fine-mapping study at 18, 22, and 18 GWAS loci on the Metabochip for their association with triglycerides (TG), high-density lipoprotein cholesterol (HDL-C), and low-density lipoprotein cholesterol (LDL-C), respectively, in individuals of African American (n = 6,832), East Asian (n = 9,449), and European (n = 10,829) ancestry. We aimed to identify the variants with strongest association at each locus, identify additional and population-specific signals, refine association signals, and assess the relative significance of previously described functional variants. Among the 58 loci, 33 exhibited evidence of association at P

Tumor necrosis factor (TNF) is a powerful pro-inflammatory cytokine and immuno-regulatory molecule, and modulates susceptibility to musculoskeletal diseases. Several meta-analyses and replicated association studies have implicated the minor 'A' variant within the TNF promoter single nucleotide polymorphism (SNP) rs361525 (-238A/G) as a risk allele in joint related disorders, including psoriatic and juvenile idiopathic arthritis, and osteolysis after joint arthroplasty. Here we characterized the effect of this variant on TNF promoter function. A transcriptional reporter, encoding the -238A variant of the TNF promoter, resulted in 2.2 to 2.8 times greater transcriptional activation versus the 'G' variant in murine macrophages when stimulated with pro-inflammatory stimuli. Bioinformatic analysis predicted a putative binding site for thyroid hormone receptor (TR) for the -238A but not the -238G allele. Overexpression of TR-α induced promoter expression 1.8-fold in the presence of the 'A' allele only. TR-α expression both potentiated and sensitized the -238A response to LPS or a titanium particulate stimulus, whilst siRNA knockdown of either THRA or THRB impaired transcriptional activation for the -238A variant only. This effect was independent of receptor-ligand binding of triiodothyronine. Immunohistochemical analysis of osteolysis interface membranes from patients undergoing revision surgery confirmed expression of TR-α within osteoclast nuclei at the resorption surface. The 'A' allele at rs361525 confers increased transcriptional activation of the TNF promoter and influences susceptibility to several arthritic conditions. This effect is modulated, at least in part, by binding of TR, which both sensitizes and potentiates transcriptional activation of the 'A' variant independent of its endogenous ligand.

Coronary heart disease (CHD) is the leading cause of mortality in both developed and developing countries worldwide. Genome-wide association studies (GWAS) have now identified 46 independent susceptibility loci for CHD, however, the biological and disease-relevant mechanisms for these associations remain elusive. The large-scale meta-analysis of GWAS recently identified in Caucasians a CHD-associated locus at chromosome 6q23.2, a region containing the transcription factor TCF21 gene. TCF21 (Capsulin/Pod1/Epicardin) is a member of the basic-helix-loop-helix (bHLH) transcription factor family, and regulates cell fate decisions and differentiation in the developing coronary vasculature. Herein, we characterize a cis-regulatory mechanism by which the lead polymorphism rs12190287 disrupts an atypical activator protein 1 (AP-1) element, as demonstrated by allele-specific transcriptional regulation, transcription factor binding, and chromatin organization, leading to altered TCF21 expression. Further, this element is shown to mediate signaling through platelet-derived growth factor receptor beta (PDGFR-β) and Wilms tumor 1 (WT1) pathways. A second disease allele identified in East Asians also appears to disrupt an AP-1-like element. Thus, both disease-related growth factor and embryonic signaling pathways may regulate CHD risk through two independent alleles at TCF21.

Circulating levels of adiponectin, a hormone produced predominantly by adipocytes, are highly heritable and are inversely associated with type 2 diabetes mellitus (T2D) and other metabolic traits. We conducted a meta-analysis of genome-wide association studies in 39,883 individuals of European ancestry to identify genes associated with metabolic disease. We identified 8 novel loci associated with adiponectin levels and confirmed 2 previously reported loci (P = 4.5×10(-8)-1.2×10(-43)). Using a novel method to combine data across ethnicities (N = 4,232 African Americans, N = 1,776 Asians, and N = 29,347 Europeans), we identified two additional novel loci. Expression analyses of 436 human adipocyte samples revealed that mRNA levels of 18 genes at candidate regions were associated with adiponectin concentrations after accounting for multiple testing (p

Apelin is a peptide ligand for an orphan G-protein coupled receptor (APJ receptor) and serves as a critical gradient for migration of mesodermal cells fated to contribute to the myocardial lineage. The present study was designed to establish a robust cardiac differentiation protocol, specifically, to evaluate the effect of apelin on directed differentiation of mouse and human embryonic stem cells (mESCs and hESCs) into cardiac lineage. Different concentrations of apelin (50, 100, 500 nM) were evaluated to determine its differentiation potential. The optimized dose of apelin was then combined with mesodermal differentiation factors, including BMP-4, activin-A, and bFGF, in a developmentally specific temporal sequence to examine the synergistic effects on cardiac differentiation. Cellular, molecular, and physiologic characteristics of the apelin-induced contractile embryoid bodies (EBs) were analyzed. It was found that 100 nM apelin resulted in highest percentage of contractile EB for mESCs while 500 nM had the highest effects on hESCs. Functionally, the contractile frequency of mESCs-derived EBs (mEBs) responded appropriately to increasing concentration of isoprenaline and diltiazem. Positive phenotype of cardiac specific markers was confirmed in the apelin-treated groups. The protocol, consisting of apelin and mesodermal differentiation factors, induced contractility in significantly higher percentage of hESC-derived EBs (hEBs), up-regulated cardiac-specific genes and cell surface markers, and increased the contractile force. In conclusion, we have demonstrated that the treatment of apelin enhanced cardiac differentiation of mouse and human ESCs and exhibited synergistic effects with mesodermal differentiation factors.

Phenotypic alteration of vascular smooth muscle cells (SMC) in response to injury or inflammation is an essential component of vascular disease. Evidence suggests that this process is dependent on epigenetic regulatory processes. P300, a histone acetyltransferase (HAT), activates crucial muscle-specific promoters in terminal (non-SMC) myocyte differentiation, and may be essential to SMC modulation as well.We performed a subanalysis examining transcriptional time-course microarray data obtained using the A404 model of SMC differentiation. Numerous chromatin remodeling genes (up to 62% of such genes on our array platform) showed significant regulation during differentiation. Members of several chromatin-remodeling families demonstrated involvement, including factors instrumental in histone modification, chromatin assembly-disassembly and DNA silencing, suggesting complex, multi-level systemic epigenetic regulation. Further, trichostatin A, a histone deacetylase inhibitor, accelerated expression of SMC differentiation markers in this model. Ontology analysis indicated a high degree of p300 involvement in SMC differentiation, with 60.7% of the known p300 interactome showing significant expression changes. Knockdown of p300 expression accelerated SMC differentiation in A404 cells and human SMCs, while inhibition of p300 HAT activity blunted SMC differentiation. The results suggest a central but complex role for p300 in SMC phenotypic modulation.Our results support the hypothesis that chromatin remodeling is important for SMC phenotypic switching, and detail wide-ranging involvement of several epigenetic modification families. Additionally, the transcriptional coactivator p300 may be partially degraded during SMC differentiation, leaving an activated subpopulation with increased HAT activity and SMC differentiation-gene specificity.

Human induced pluripotent stem cells (hiPSCs) generated by de-differentiation of adult somatic cells offer potential solutions for the ethical issues surrounding human embryonic stem cells (hESCs), as well as their immunologic rejection after cellular transplantation. However, although hiPSCs have been described as "embryonic stem cell-like", these cells have a distinct gene expression pattern compared to hESCs, making incomplete reprogramming a potential pitfall. It is unclear to what degree the difference in tissue of origin may contribute to these gene expression differences. To answer these important questions, a careful transcriptional profiling analysis is necessary to investigate the exact reprogramming state of hiPSCs, as well as analysis of the impression, if any, of the tissue of origin on the resulting hiPSCs. In this study, we compare the gene profiles of hiPSCs derived from fetal fibroblasts, neonatal fibroblasts, adipose stem cells, and keratinocytes to their corresponding donor cells and hESCs. Our analysis elucidates the overall degree of reprogramming within each hiPSC line, as well as the "distance" between each hiPSC line and its donor cell. We further identify genes that have a similar mode of regulation in hiPSCs and their corresponding donor cells compared to hESCs, allowing us to specify core sets of donor genes that continue to be expressed in each hiPSC line. We report that residual gene expression of the donor cell type contributes significantly to the differences among hiPSCs and hESCs, and adds to the incompleteness in reprogramming. Specifically, our analysis reveals that fetal fibroblast-derived hiPSCs are closer to hESCs, followed by adipose, neonatal fibroblast, and keratinocyte-derived hiPSCs.

Single nucleotide polymorphisms (SNPs) in the peroxisome proliferator-activated receptor γ (PPARG) gene have been associated with cardiovascular risk factors, particularly obesity and diabetes. We assessed the relationship between 4 PPARG SNPs (C-681G, C-689T, Pro12Ala, and C1431T) and coronary heart disease (CHD) in the PRIME (249 cases/494 controls, only men) and ADVANCE (1,076 cases/805 controls, men or women) studies. In PRIME, homozygote individuals for the minor allele of the PPARG C-689T, Pro12Ala, and C1431T SNPs tended to have a higher risk of CHD than homozygote individuals for the frequent allele (adjusted OR [95% CI] = 3.43 [0.96–12.27], P=.058, 3.41 [0.95–12.22], P=.060 and 5.10 [0.99–26.37], P=.050, resp.). No such association could be detected in ADVANCE. Haplotype distributions were similar in cases and control in both studies. A meta-analysis on the Pro12Ala SNP, based on our data and 11 other published association studies (6,898 CHD cases/11,287 controls), revealed that there was no evidence for a significant association under the dominant model (OR=0.99 [0.92–1.07], P=.82). However, there was a borderline association under the recessive model (OR=1.29 [0.99–1.67], P=.06) that became significant when considering men only (OR=1.73 [1.20–2.48], P=.003). In conclusion, the PPARG Ala12Ala genotype might be associated with a higher CHD risk in men but further confirmation studies are needed.

BACKGROUND: Recently, the type 2 diabetes medication, rosiglitazone, has come under scrutiny for possibly increasing the risk of cardiac disease and death. To investigate the effects of rosiglitazone on the diabetic heart, we performed cardiac transcriptional profiling and imaging studies of a murine model of type 2 diabetes, the C57BL/KLS-lepr(db)/lepr(db) (db/db) mouse. METHODS AND FINDINGS: We compared cardiac gene expression profiles from three groups: untreated db/db mice, db/db mice after rosiglitazone treatment, and non-diabetic db/+ mice. Prior to sacrifice, we also performed cardiac magnetic resonance (CMR) and echocardiography. As expected, overall the db/db gene expression signature was markedly different from control, but to our surprise was not significantly reversed with rosiglitazone. In particular, we have uncovered a number of rosiglitazone modulated genes and pathways that may play a role in the pathophysiology of the increase in cardiac mortality as seen in several recent meta-analyses. Specifically, the cumulative upregulation of (1) a matrix metalloproteinase gene that has previously been implicated in plaque rupture, (2) potassium channel genes involved in membrane potential maintenance and action potential generation, and (3) sphingolipid and ceramide metabolism-related genes, together give cause for concern over rosiglitazone's safety. Lastly, in vivo imaging studies revealed minimal differences between rosiglitazone-treated and untreated db/db mouse hearts, indicating that rosiglitazone's effects on gene expression in the heart do not immediately turn into detectable gross functional changes. CONCLUSIONS: This study maps the genomic expression patterns in the hearts of the db/db murine model of diabetes and illustrates the impact of rosiglitazone on these patterns. The db/db gene expression signature was markedly different from control, and was not reversed with rosiglitazone. A smaller number of unique and interesting changes in gene expression were noted with rosiglitazone treatment. Further study of these genes and molecular pathways will provide important insights into the cardiac decompensation associated with both diabetes and rosiglitazone treatment.