Search

Atrial fibrillation (AF) is the most common sustained arrhythmia in humans, yet the molecular basis of AF remains incompletely understood. To determine the cell type-specific transcriptional changes underlying AF, we perform single-nucleus RNA-seq (snRNA-seq) on left atrial (LA) samples from patients with AF and controls. From more than 175,000 nuclei we find that only cardiomyocytes (CMs) and macrophages (MΦs) have a significant number of differentially expressed genes in patients with AF. Attractin Like 1 (ATRNL1) was overexpressed in CMs among patients with AF and localized to the intercalated disks. Further, in both knockdown and overexpression experiments we identify a potent role for ATRNL1 in cell stress response, and in the modulation of the cardiac action potential. Finally, we detect an unexpected expression pattern for a leading AF candidate gene, KCNN3. In sum, we uncover a role for ATRNL1 which may serve as potential therapeutic target for this common arrhythmia., Competing Interests: Competing interests C.A.K. is an employee of Bayer US LLC (a subsidiary of Bayer AG) and may own stock in Bayer. The Precision Cardiology Laboratory is joint effort between the Broad Institute and Bayer AG. P.T.E. has received sponsored research support from Bayer AG, Novo-Nordisk, Bristol Myers Squibb and Pfizer; he has also served on advisory boards or consulted for Bayer AG. C.R. is a full-time employee at GSK as of July 2024. All remaining authors declare no competing interests., (© 2024. The Author(s).)

Background: Many specialized cells in adult organs acquire a state of cell cycle arrest and quiescence through unknown mechanisms. Our limited understanding of mammalian cell cycle arrest is derived primarily from cell culture models. Adult mammalian cardiomyocytes, a classic example of cell cycle arrested cells, exit the cell cycle postnatally and remain in an arrested state for the life of the organism. Cardiomyocytes can be induced to re-enter the cell cycle by YAP5SA, an active form of the Hippo signaling pathway effector YAP., Methods: We performed clonal analyses to determine the cell kinetics of YAP5SA cardiomyocytes. We also performed single-cell RNA sequencing, marker gene analysis, and functional studies to examine how YAP5SA cardiomyocytes progress through the cell cycle., Results: We discovered that YAP5SA-expressing cardiomyocytes divided efficiently, with >20% of YAP5SA cardiomyocyte clones containing ≥2 cardiomyocytes. YAP5SA cardiomyocytes re-entered cell cycle at the G1/S transition and had an S phase lasting ≈48 hours. Sarcomere disassembly is required for cardiomyocyte progression from S to G2 phase and the induction of mitotic rounding. Although oscillatory Cdk expression was induced in YAP5SA cardiomyocytes, these cells inefficiently progressed through G2 phase. This is improved by inhibiting P21 function, implicating checkpoint activity as an additional barrier to YAP5SA-induced cardiomyocyte division., Conclusions: Our data reveal that YAP5SA overcomes the mechanically constrained myocardial microenvironment to induce mitotic rounding with cardiomyocyte division, thus providing new insights into the in vivo mechanisms that maintain cell cycle quiescence in adult mammals.

Fibrotic diseases affect multiple organs and are associated with morbidity and mortality. To examine organ-specific and shared biologic mechanisms that underlie fibrosis in different organs, we developed machine learning models to quantify T1 time, a marker of interstitial fibrosis, in the liver, pancreas, heart and kidney among 43,881 UK Biobank participants who underwent magnetic resonance imaging. In phenome-wide association analyses, we demonstrate the association of increased organ-specific T1 time, reflecting increased interstitial fibrosis, with prevalent diseases across multiple organ systems. In genome-wide association analyses, we identified 27, 18, 11 and 10 independent genetic loci associated with liver, pancreas, myocardial and renal cortex T1 time, respectively. There was a modest genetic correlation between the examined organs. Several loci overlapped across the examined organs implicating genes involved in a myriad of biologic pathways including metal ion transport (SLC39A8, HFE and TMPRSS6), glucose metabolism (PCK2), blood group antigens (ABO and FUT2), immune function (BANK1 and PPP3CA), inflammation (NFKB1) and mitosis (CENPE). Finally, we found that an increasing number of organs with T1 time falling in the top quintile was associated with increased mortality in the population. Individuals with a high burden of fibrosis in ≥3 organs had a 3-fold increase in mortality compared to those with a low burden of fibrosis across all examined organs in multivariable-adjusted analysis (hazard ratio = 3.31, 95% confidence interval 1.77-6.19; P = 1.78 × 10 -4 ). By leveraging machine learning to quantify T1 time across multiple organs at scale, we uncovered new organ-specific and shared biologic pathways underlying fibrosis that may provide therapeutic targets., (© 2024. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Protein-specific Chromatin Conformation Capture (3C)-based technologies have become essential for identifying distal genomic interactions with critical roles in gene regulation. The standard techniques include Chromatin Interaction Analysis by Paired-End Tag (ChIA-PET), in situ Hi-C followed by chromatin immunoprecipitation (HiChIP) also known as PLAC-seq. To identify chromatin interactions from these data, a variety of computational methods have emerged. Although these state-of-art methods address many issues with loop calling, only few methods can fit different data types simultaneously, and the accuracy as well as the efficiency these approaches remains limited. Here we have generated a pipeline, MMCT-Loop, which ensures the accurate identification of strong loops as well as dynamic or weak loops through a mixed model. MMCT-Loop outperforms existing methods in accuracy, and the detected loops show higher activation functionality. To highlight the utility of MMCT-Loop, we applied it to conformational data derived from neural stem cell (NSCs) and uncovered several previously unidentified regulatory regions for key master regulators of stem cell identity. MMCT-Loop is an accurate and efficient loop caller for targeted conformation capture data, which supports raw data or pre-processed valid pairs as input, the output interactions are formatted and easily uploaded to a genome browser for visualization., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Background: ZFHX3 (zinc finger homeobox 3), a gene that encodes a large transcription factor, is at the second-most significantly associated locus with atrial fibrillation (AF), but its function in the heart is unknown. This study aims to identify causative genetic variation related to AF at the ZFHX3 locus and examine the impact of Zfhx3 loss on cardiac function in mice., Methods: CRISPR-Cas9 genome editing, chromatin immunoprecipitation, and luciferase assays in pluripotent stem cell-derived cardiomyocytes were used to identify causative genetic variation related to AF at the ZFHX3 locus. Cardiac function was assessed by echocardiography, magnetic resonance imaging, electrophysiology studies, calcium imaging, and RNA sequencing in mice with heterozygous and homozygous cardiomyocyte-restricted Zfhx3 loss ( Zfhx3 Het and knockout, respectively). Human cardiac single-nucleus ATAC (assay for transposase-accessible chromatin)-sequencing data was analyzed to determine which genes in atrial cardiomyocytes are directly regulated by ZFHX3 ., Results: We found single-nucleotide polymorphism (SNP) rs12931021 modulates an enhancer regulating ZFHX3 expression, and the AF risk allele is associated with decreased ZFHX3 transcription. We observed a gene-dose response in AF susceptibility with Zfhx3 knockout mice having higher incidence, frequency, and burden of AF than Zfhx3 Het and wild-type mice, with alterations in conduction velocity, atrial action potential duration, calcium handling and the development of atrial enlargement and thrombus, and dilated cardiomyopathy. Zfhx3 loss results in atrial-specific differential effects on genes and signaling pathways involved in cardiac pathophysiology and AF., Conclusions: Our findings implicate ZFHX3 as the causative gene at the 16q22 locus for AF, and cardiac abnormalities caused by loss of cardiac Zfhx3 are due to atrial-specific dysregulation of pathways involved in AF susceptibility. Together, these data reveal a novel and important role for Zfhx3 in the control of cardiac genes and signaling pathways essential for normal atrial function., Competing Interests: Disclosures P.T. Ellinor receives sponsored research support from Bayer AG, IBM Research, Bristol Myers Squibb, and Pfizer; he has also served on advisory boards or consulted for Bayer AG, MyoKardia, and Novartis. The other authors report no conflicts.

Mapping gene networks requires large amounts of transcriptomic data to learn the connections between genes, which impedes discoveries in settings with limited data, including rare diseases and diseases affecting clinically inaccessible tissues. Recently, transfer learning has revolutionized fields such as natural language understanding 1,2 and computer vision 3 by leveraging deep learning models pretrained on large-scale general datasets that can then be fine-tuned towards a vast array of downstream tasks with limited task-specific data. Here, we developed a context-aware, attention-based deep learning model, Geneformer, pretrained on a large-scale corpus of about 30 million single-cell transcriptomes to enable context-specific predictions in settings with limited data in network biology. During pretraining, Geneformer gained a fundamental understanding of network dynamics, encoding network hierarchy in the attention weights of the model in a completely self-supervised manner. Fine-tuning towards a diverse panel of downstream tasks relevant to chromatin and network dynamics using limited task-specific data demonstrated that Geneformer consistently boosted predictive accuracy. Applied to disease modelling with limited patient data, Geneformer identified candidate therapeutic targets for cardiomyopathy. Overall, Geneformer represents a pretrained deep learning model from which fine-tuning towards a broad range of downstream applications can be pursued to accelerate discovery of key network regulators and candidate therapeutic targets., (© 2023. The Author(s), under exclusive licence to Springer Nature Limited.)

Myocardial interstitial fibrosis is associated with cardiovascular disease and adverse prognosis. Here, to investigate the biological pathways that underlie fibrosis in the human heart, we developed a machine learning model to measure native myocardial T1 time, a marker of myocardial fibrosis, in 41,505 UK Biobank participants who underwent cardiac magnetic resonance imaging. Greater T1 time was associated with diabetes mellitus, renal disease, aortic stenosis, cardiomyopathy, heart failure, atrial fibrillation, conduction disease and rheumatoid arthritis. Genome-wide association analysis identified 11 independent loci associated with T1 time. The identified loci implicated genes involved in glucose transport (SLC2A12), iron homeostasis (HFE, TMPRSS6), tissue repair (ADAMTSL1, VEGFC), oxidative stress (SOD2), cardiac hypertrophy (MYH7B) and calcium signaling (CAMK2D). Using a transforming growth factor β1-mediated cardiac fibroblast activation assay, we found that 9 of the 11 loci consisted of genes that exhibited temporal changes in expression or open chromatin conformation supporting their biological relevance to myofibroblast cell state acquisition. By harnessing machine learning to perform large-scale quantification of myocardial interstitial fibrosis using cardiac imaging, we validate associations between cardiac fibrosis and disease, and identify new biologically relevant pathways underlying fibrosis., (© 2023. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Ischemic cardiomyopathy (ICM) is the leading cause of heart failure worldwide, yet the cellular and molecular signature of this disease is largely unclear. Using single-nucleus RNA sequencing (snRNA-seq) and integrated computational analyses, we profile the transcriptomes of over 99,000 human cardiac nuclei from the non-infarct region of the left ventricle of 7 ICM transplant recipients and 8 non-failing (NF) controls. We find the cellular composition of the ischemic heart is significantly altered, with decreased cardiomyocytes and increased proportions of lymphatic, angiogenic, and arterial endothelial cells in patients with ICM. We show that there is increased LAMININ signaling from endothelial cells to other cell types in ICM compared with NF. Finally, we find that the transcriptional changes that occur in ICM are similar to those in hypertrophic and dilated cardiomyopathies and that the mining of these combined datasets can identify druggable genes that could be used to target end-stage heart failure., Competing Interests: Declaration of interests C.A.K. is an employee of Bayer US LLC (a subsidiary of Bayer AG) and may own stock in Bayer. S.H., S.B., and C.A.K. were full-time employees of Bayer when this work was performed. P.T.E. has received sponsored research support from Bayer AG, IBM Health, Bristol Myers Squibb, and Pfizer; he has also served on advisory boards or consulted for Bayer AG, MyoKardia, and Novartis. K.B.M. has research grant funding from Amgen, USA and has also served on advisory boards for MyoKardia, Bristol-Myers Squibb, and Pfizer., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Neural circuit plasticity and sensory response dynamics depend on forming new synaptic connections. Despite recent advances toward understanding the consequences of circuit plasticity, the mechanisms driving circuit plasticity are unknown. Adult-born neurons within the olfactory bulb have proven to be a powerful model for studying circuit plasticity, providing a broad and accessible avenue into neuron development, migration, and circuit integration. We and others have shown that efficient adult-born neuron circuit integration hinges on presynaptic activity in the form of diverse signaling peptides. Here, we demonstrate a novel oxytocin-dependent mechanism of adult-born neuron synaptic maturation and circuit integration. We reveal spatial and temporal enrichment of oxytocin receptor expression within adult-born neurons in the murine olfactory bulb, with oxytocin receptor expression peaking during activity-dependent integration. Using viral labeling, confocal microscopy, and cell type-specific RNA-seq, we demonstrate that oxytocin receptor signaling promotes synaptic maturation of newly integrating adult-born neurons by regulating their morphological development and expression of mature synaptic AMPARs and other structural proteins., (© 2022 Pekarek et al.; Published by Cold Spring Harbor Laboratory Press.)

Heart failure encompasses a heterogeneous set of clinical features that converge on impaired cardiac contractile function 1,2 and presents a growing public health concern. Previous work has highlighted changes in both transcription and protein expression in failing hearts 3,4 , but may overlook molecular changes in less prevalent cell types. Here we identify extensive molecular alterations in failing hearts at single-cell resolution by performing single-nucleus RNA sequencing of nearly 600,000 nuclei in left ventricle samples from 11 hearts with dilated cardiomyopathy and 15 hearts with hypertrophic cardiomyopathy as well as 16 non-failing hearts. The transcriptional profiles of dilated or hypertrophic cardiomyopathy hearts broadly converged at the tissue and cell-type level. Further, a subset of hearts from patients with cardiomyopathy harbour a unique population of activated fibroblasts that is almost entirely absent from non-failing samples. We performed a CRISPR-knockout screen in primary human cardiac fibroblasts to evaluate this fibrotic cell state transition; knockout of genes associated with fibroblast transition resulted in a reduction of myofibroblast cell-state transition upon TGFβ1 stimulation for a subset of genes. Our results provide insights into the transcriptional diversity of the human heart in health and disease as well as new potential therapeutic targets and biomarkers for heart failure., (© 2022. The Author(s), under exclusive licence to Springer Nature Limited.)

The heart, the first organ to develop in the embryo, undergoes complex morphogenesis that when defective results in congenital heart disease (CHD). With current therapies, more than 90% of patients with CHD survive into adulthood, but many suffer premature death from heart failure and non-cardiac causes 1 . Here, to gain insight into this disease progression, we performed single-nucleus RNA sequencing on 157,273 nuclei from control hearts and hearts from patients with CHD, including those with hypoplastic left heart syndrome (HLHS) and tetralogy of Fallot, two common forms of cyanotic CHD lesions, as well as dilated and hypertrophic cardiomyopathies. We observed CHD-specific cell states in cardiomyocytes, which showed evidence of insulin resistance and increased expression of genes associated with FOXO signalling and CRIM1. Cardiac fibroblasts in HLHS were enriched in a low-Hippo and high-YAP cell state characteristic of activated cardiac fibroblasts. Imaging mass cytometry uncovered a spatially resolved perivascular microenvironment consistent with an immunodeficient state in CHD. Peripheral immune cell profiling suggested deficient monocytic immunity in CHD, in agreement with the predilection in CHD to infection and cancer 2 . Our comprehensive phenotyping of CHD provides a roadmap towards future personalized treatments for CHD., (© 2022. The Author(s), under exclusive licence to Springer Nature Limited.)

Chromatin conformation capture (3C)-based technologies have enabled the accurate detection of topological genomic interactions, and the adoption of ChIP techniques to 3C-based protocols makes it possible to identify long-range interactions. To analyze these large and complex datasets, computational methods are undergoing rapid and expansive evolution. Thus, a thorough evaluation of these analytical pipelines is necessary to identify which commonly used algorithms and processing pipelines need to be improved. Here we present a comprehensive benchmark framework, Bacon, to evaluate the performance of several computational methods. Finally, we provide practical recommendations for users working with HiChIP and/or ChIA-PET analyses., (© 2022. The Author(s).)

Combined methylmalonic acidemia and homocystinuria (cblC) is the most common inborn error of intracellular cobalamin metabolism and due to mutations in Methylmalonic Aciduria type C and Homocystinuria (MMACHC). Recently, mutations in the transcriptional regulators HCFC1 and RONIN (THAP11) were shown to result in cellular phenocopies of cblC. Since HCFC1/RONIN jointly regulate MMACHC, patients with mutations in these factors suffer from reduced MMACHC expression and exhibit a cblC-like disease. However, additional de-regulated genes and the resulting pathophysiology is unknown. Therefore, we have generated mouse models of this disease. In addition to exhibiting loss of Mmachc, metabolic perturbations, and developmental defects previously observed in cblC, we uncovered reduced expression of target genes that encode ribosome protein subunits. We also identified specific phenotypes that we ascribe to deregulation of ribosome biogenesis impacting normal translation during development. These findings identify HCFC1/RONIN as transcriptional regulators of ribosome biogenesis during development and their mutation results in complex syndromes exhibiting aspects of both cblC and ribosomopathies., (© 2022. The Author(s).)

The cardiac conduction system is a network of heterogeneous cell population that initiates and propagates electric excitations in the myocardium. Purkinje fibers, a network of specialized myocardial cells, comprise the distal end of the conduction system in the ventricles. The developmental origins of Purkinje fibers and their roles during cardiac physiology and arrhythmia have been reported. However, it is not clear if they play a role during ischemic injury and heart regeneration. Here we introduce a novel tamoxifen-inducible Cre allele that specifically labels a broad range of components in the cardiac conduction system while excludes other cardiac cell types and vital organs. Using this new allele, we investigated the cellular and molecular response of Purkinje fibers to myocardial injury. In a neonatal mouse myocardial infarction model, we observed significant increase in Purkinje cell number in regenerating myocardium. RNA-Seq analysis using laser-captured Purkinje fibers showed a unique transcriptomic response to myocardial infarction. Our finds suggest a novel role of cardiac Purkinje fibers in heart injury., Competing Interests: Declaration of competing interest None., (Copyright © 2021 Elsevier Inc. All rights reserved.)

Histone variants contribute to the complexity of the chromatin landscape and play an integral role in defining DNA domains and regulating gene expression. The histone H3 variant H3.3 is incorporated into genic elements independent of DNA replication by its chaperone HIRA. Here we demonstrate that Hira is required for the self-renewal of adult hematopoietic stem cells (HSCs) and to restrain erythroid differentiation. Deletion of Hira led to rapid depletion of HSCs while differentiated hematopoietic cells remained largely unaffected. Depletion of HSCs after Hira deletion was accompanied by increased expression of bivalent and erythroid genes, which was exacerbated upon cell division and paralleled increased erythroid differentiation. Assessing H3.3 occupancy identified a subset of polycomb-repressed chromatin in HSCs that depends on HIRA to maintain the inaccessible, H3.3-occupied state for gene repression. HIRA-dependent H3.3 incorporation thus defines distinct repressive chromatin that represses erythroid differentiation of HSCs., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

The adult mammalian heart's potential for regeneration is very inefficient. Importantly, adult mammalian cardiomyocytes (CMs) are characterized as a cell population with very limited mitotic potential. Conversely, the neonatal mouse heart possesses a brief, yet robust, regenerative capacity within the first week of life. Cell type-specific enrichment procedures are essential for characterizing the full spectrum of epigenomic landscapes and gene regulatory networks deployed by mammalian CMs. In this chapter, we describe a protocol useful for purifying CM nuclei from mammalian cardiac tissue. Furthermore, we detail a low-input procedure suitable for the parallel genome-wide profiling of chromatin accessibility, histone modifications, and transcription factor-binding sites. The CM nuclei purified using this process are suitable for multi-omic profiling approaches.

Transcriptional enhancers commonly work over long genomic distances to precisely regulate spatiotemporal gene expression patterns. Dissecting the promoters physically contacted by these distal regulatory elements is essential for understanding developmental processes as well as the role of disease-associated risk variants. Modern proximity-ligation assays, like HiChIP and ChIA-PET, facilitate the accurate identification of long-range contacts between enhancers and promoters. However, these assays are technically challenging, expensive, and time-consuming, making it difficult to investigate enhancer topologies, especially in uncharacterized cell types. To overcome these shortcomings, we therefore designed LoopPredictor, an ensemble machine learning model, to predict genome topology for cell types which lack long-range contact maps. To enrich for functional enhancer-promoter loops over common structural genomic contacts, we trained LoopPredictor with both H3K27ac and YY1 HiChIP data. Moreover, the integration of several related multi-omics features facilitated identifying and annotating the predicted loops. LoopPredictor is able to efficiently identify cell type-specific enhancer-mediated loops, and promoter-promoter interactions, with a modest feature input requirement. Comparable to experimentally generated H3K27ac HiChIP data, we found that LoopPredictor was able to identify functional enhancer loops. Furthermore, to explore the cross-species prediction capability of LoopPredictor, we fed mouse multi-omics features into a model trained on human data and found that the predicted enhancer loops outputs were highly conserved. LoopPredictor enables the dissection of cell type-specific long-range gene regulation and can accelerate the identification of distal disease-associated risk variants., (© 2020 Tang et al.; Published by Cold Spring Harbor Laboratory Press.)

Genome-wide association studies found that increased risk for atrial fibrillation (AF), the most common human heart arrhythmia, is associated with noncoding sequence variants located in proximity to PITX2 Cardiomyocyte-specific epigenomic and comparative genomics uncovered 2 AF-associated enhancers neighboring PITX2 with varying conservation in mice. Chromosome conformation capture experiments in mice revealed that the Pitx2c promoter directly contacted the AF-associated enhancer regions. CRISPR/Cas9-mediated deletion of a 20-kb topologically engaged enhancer led to reduced Pitx2c transcription and AF predisposition. Allele-specific chromatin immunoprecipitation sequencing on hybrid heterozygous enhancer knockout mice revealed that long-range interaction of an AF-associated region with the Pitx2c promoter was required for maintenance of the Pitx2 c promoter chromatin state. Long-range looping was mediated by CCCTC-binding factor (CTCF), since genetic disruption of the intronic CTCF-binding site caused reduced Pitx2c expression, AF predisposition, and diminished active chromatin marks on Pitx2 AF risk variants located at 4q25 reside in genomic regions possessing long-range transcriptional regulatory functions directed at PITX2 ., Competing Interests: Competing interest statement: P.T.E. is supported by a grant from Bayer AG to the Broad Institute focused on the genetics and therapeutics of cardiovascular disease and has consulted for Bayer AG, Novartis, and Quest Diagnostics. J.F.M. is a founder and owns shares in Yap Therapeutics.

Cardiac fibroblasts (CFs) respond to injury by transitioning through multiple cell states, including resting CFs, activated CFs, and myofibroblasts. We report here that Hippo signaling cell-autonomously regulates CF fate transitions and proliferation, and non-cell-autonomously regulates both myeloid and CF activation in the heart. Conditional deletion of Hippo pathway kinases, Lats1 and Lats2 , in uninjured CFs initiated a self-perpetuating fibrotic response in the adult heart that was exacerbated by myocardial infarction (MI). Single cell transcriptomics showed that uninjured Lats1/2 mutant CFs spontaneously transitioned to a myofibroblast cell state. Through gene regulatory network reconstruction, we found that Hippo-deficient myofibroblasts deployed a network of transcriptional regulators of endoplasmic reticulum (ER) stress, and the unfolded protein response (UPR) consistent with elevated secretory activity. We observed an expansion of myeloid cell heterogeneity in uninjured Lats1/2 CKO hearts with similarity to cells recovered from control hearts post-MI. Integrated genome-wide analysis of Yap chromatin occupancy revealed that Yap directly activates myofibroblast cell identity genes, the proto-oncogene Myc , and an array of genes encoding pro-inflammatory factors through enhancer-promoter looping. Our data indicate that Lats1/2 maintain the resting CF cell state through restricting the Yap-induced injury response., (© 2019 Xiao et al.; Published by Cold Spring Harbor Laboratory Press.)

Disease-associated genetic variants that lie in non-coding regions found by genome-wide association studies are thought to alter the functionality of transcription regulatory elements and target gene expression. To uncover causal genetic variants, variant regulatory elements and their target genes, here we cross-reference human transcriptomic, epigenomic and chromatin conformation datasets. Of 104 genetic variant regions associated with atrial fibrillation candidate target genes are prioritized. We optimize EMERGE enhancer prediction and use accessible chromatin profiles of human atrial cardiomyocytes to more accurately predict cardiac regulatory elements and identify hundreds of sub-threshold variants that co-localize with regulatory elements. Removal of mouse homologues of atrial fibrillation-associated regions in vivo uncovers a distal regulatory region involved in Gja1 (Cx43) expression. Our analyses provide a shortlist of genes likely affected by atrial fibrillation-associated variants and provide variant regulatory elements in each region that link genetic variation and target gene regulation, helping to focus future investigations.

Tet-mediated DNA demethylation plays an important role in shaping the epigenetic landscape and chromatin accessibility to control gene expression. While several studies demonstrated pivotal roles of Tet in regulating embryonic development, little is known about their functions in heart development. Here we analyze DNA methylation and hydroxymethylation dynamics during early cardiac development in both human and mice. We find that cardiac-specific deletion of Tet2 and Tet3 in mice (Tet2/3-DKO) leads to ventricular non-compaction cardiomyopathy (NCC) with embryonic lethality. Single-cell RNA-seq analyses reveal a reduction in cardiomyocyte numbers and transcriptional reprogramming in cardiac tissues upon Tet2/3 depletion. Impaired DNA demethylation and reduced chromatin accessibility in Tet2/3-DKO mice further compromised Ying-yang1 (YY1) binding to its genomic targets, and perturbed high-order chromatin organization at key genes involved in heart development. Our studies provide evidence of the physiological role of Tet in regulating DNA methylation dynamics and chromatin organization during early heart development.

Background: Surviving cells in the postinfarction border zone are subjected to intense fluctuations of their microenvironment. Recently, border zone cardiomyocytes have been specifically implicated in cardiac regeneration. Here, we defined their unique transcriptional and regulatory properties, and comprehensively validated new molecular markers, including Nppb, encoding B-type natriuretic peptide, after infarction., Methods: Transgenic reporter mice were used to identify the Nppb-positive border zone after myocardial infarction. Transcriptome analysis of remote, border, and infarct zones and of purified cardiomyocyte nuclei was performed using RNA-sequencing. Top candidate genes displaying border zone spatial specificity were histologically validated in ischemic human hearts. Mice in which Nppb was deleted by genome editing were subjected to myocardial infarction. Chromatin accessibility landscapes of border zone and control cardiomyocyte nuclei were assessed by using assay for transposase-accessible chromatin using sequencing., Results: We identified the border zone as a spatially confined region transcriptionally distinct from the remote myocardium. The transcriptional response of the border zone was much stronger than that of the remote ventricular wall, involving acute downregulation of mitochondrial oxidative phosphorylation, fatty acid metabolism, calcium handling, and sarcomere function, and the activation of a stress-response program. Analysis of infarcted human hearts revealed that the transcriptionally discrete border zone is conserved in humans, and led to the identification of novel conserved border zone markers including NPPB, ANKRD1, DES, UCHL1, JUN, and FOXP1. Homozygous Nppb mutant mice developed acute and lethal heart failure after myocardial infarction, indicating that B-type natriuretic peptide is required to preserve postinfarct heart function. Assay for transposase-accessible chromatin using sequencing revealed thousands of cardiomyocyte lineage-specific MEF2-occupied regulatory elements that lost accessibility in the border zone. Putative injury-responsive enhancers that gained accessibility were highly associated with AP-1 (activator protein 1) binding sites. Nuclear c-Jun, a component of AP-1, was observed specifically in border zone cardiomyocytes., Conclusions: Cardiomyocytes in a discrete zone bordering the infarct switch from a MEF2-driven homeostatic lineage-specific to an AP-1-driven injury-induced gene expression program. This program is conserved between mouse and human, and includes Nppb expression, which is required to prevent acute heart failure after infarction.

Tissue injury induces changes in cellular identity, but the underlying molecular mechanisms remain obscure. Here, we show that upon damage in a mouse model, epidermal cells at the wound edge convert to an embryonic-like state, altering particularly the cytoskeletal/extracellular matrix (ECM) components and differentiation program. We show that SOX11 and its closest relative SOX4 dictate embryonic epidermal state, regulating genes involved in epidermal development as well as cytoskeletal/ECM organization. Correspondingly, postnatal induction of SOX11 represses epidermal terminal differentiation while deficiency of Sox11 and Sox4 accelerates differentiation and dramatically impairs cell motility and re-epithelialization. Amongst the embryonic genes reactivated at the wound edge, we identify fascin actin-bundling protein 1 (FSCN1) as a critical direct target of SOX11 and SOX4 regulating cell migration. Our study identifies the reactivated embryonic gene program during wound repair and demonstrates that SOX11 and SOX4 play a central role in this process.

Oncogenic mutations confer on cells the ability to propagate indefinitely, but whether oncogenes alter the cell fate of these cells is unknown. Here, we show that the transcriptional regulator PRDM16s causes oncogenic fate conversion by transforming cells fated to form platelets and erythrocytes into myeloid leukemia stem cells (LSCs). Prdm16s expression in megakaryocyte-erythroid progenitors (MEPs), which normally lack the potential to generate granulomonocytic cells, caused AML by converting MEPs into LSCs. Prdm16s blocked megakaryocytic/erythroid potential by interacting with super enhancers and activating myeloid master regulators, including PU.1. A CRISPR dropout screen confirmed that PU.1 is required for Prdm16s -induced leukemia. Ablating PU.1 attenuated leukemogenesis and reinstated the megakaryocytic/erythroid potential of leukemic MEPs in mouse models and human AML with PRDM16 rearrangement. Thus, oncogenic PRDM16 s expression gives MEPs an LSC fate by activating myeloid gene regulatory networks., (© 2019 by The American Society of Hematology.)

The Pitx2 gene encodes a homeobox transcription factor that is required for mammalian development. Disruption of PITX2 expression in humans causes congenital heart diseases and is associated with atrial fibrillation; however, the cellular and molecular processes dictated by Pitx2 during cardiac ontogeny remain unclear. To characterize the role of Pitx2 during murine heart development we sequenced over 75,000 single cardiac cell transcriptomes between two key developmental timepoints in control and Pitx2 null embryos. We found that cardiac cell composition was dramatically altered in mutants at both E10.5 and E13.5. Interestingly, the differentiation dynamics of both anterior and posterior second heart field-derived progenitor cells were disrupted in Pitx2 mutants. We also uncovered evidence for defects in left-right asymmetry within atrial cardiomyocyte populations. Furthermore, we were able to detail defects in cardiac outflow tract and valve development associated with Pitx2 Our findings offer insight into Pitx2 function and provide a compilation of gene expression signatures for further detailing the complexities of heart development that will serve as the foundation for future studies of cardiac morphogenesis, congenital heart disease and arrhythmogenesis., Competing Interests: Competing interestsThe authors declare no competing or financial interests., (© 2019. Published by The Company of Biologists Ltd.)

In response to retinal damage, the Müller glial cells (MGs) of the zebrafish retina have the ability to undergo a cellular reprogramming event in which they enter the cell cycle and divide asymmetrically, thereby producing multipotent retinal progenitors capable of regenerating lost retinal neurons. However, mammalian MGs do not exhibit such a proliferative and regenerative ability. Here, we identify Hippo pathway-mediated repression of the transcription cofactor YAP as a core regulatory mechanism that normally blocks mammalian MG proliferation and cellular reprogramming. MG-specific deletion of Hippo pathway components Lats1 and Lats2, as well as transgenic expression of a Hippo non-responsive form of YAP (YAP5SA), resulted in dramatic Cyclin D1 upregulation, loss of adult MG identity, and attainment of a highly proliferative, progenitor-like cellular state. Our results reveal that mammalian MGs may have latent regenerative capacity that can be stimulated by repressing Hippo signaling., (Copyright © 2019 The Author(s). Published by Elsevier Inc. All rights reserved.)

The mammalian cochlea loses its ability to regenerate new hair cells prior to the onset of hearing. In contrast, the adult vestibular system can produce new hair cells in response to damage, or by reprogramming of supporting cells with the hair cell transcription factor Atoh1. We used RNA-seq and ATAC-seq to probe the transcriptional and epigenetic responses of utricle supporting cells to damage and Atoh1 transduction. We show that the regenerative response of the utricle correlates with a more accessible chromatin structure in utricle supporting cells compared to their cochlear counterparts. We also provide evidence that Atoh1 transduction of supporting cells is able to promote increased transcriptional accessibility of some hair cell genes. Our study offers a possible explanation for regenerative differences between sensory organs of the inner ear, but shows that additional factors to Atoh1 may be required for optimal reprogramming of hair cell fate., Competing Interests: HJ, MH, LT, KS, WC, HZ, HY, JL, CZ, JM, NS, AG No competing interests declared, (© 2019, Jen et al.)

Specialized adult somatic cells, such as cardiomyocytes (CMs), are highly differentiated with poor renewal capacity, an integral reason underlying organ failure in disease and aging. Among the least renewable cells in the human body, CMs renew approximately 1% annually. Consistent with poor CM turnover, heart failure is the leading cause of death. Here, we show that an active version of the Hippo pathway effector YAP, termed YAP5SA, partially reprograms adult mouse CMs to a more fetal and proliferative state. One week after induction, 19% of CMs that enter S-phase do so twice, CM number increases by 40%, and YAP5SA lineage CMs couple to pre-existing CMs. Genomic studies showed that YAP5SA increases chromatin accessibility and expression of fetal genes, partially reprogramming long-lived somatic cells in vivo to a primitive, fetal-like, and proliferative state., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Cellular heterogeneity within the mammalian brain poses a challenge toward understanding its complex functions. Within the olfactory bulb, odor information is processed by subtypes of inhibitory interneurons whose heterogeneity and functionality are influenced by ongoing adult neurogenesis. To investigate this cellular heterogeneity and better understand adult-born neuron development, we utilized single-cell RNA sequencing and computational modeling to reveal diverse and transcriptionally distinct neuronal and nonneuronal cell types. We also analyzed molecular changes during adult-born interneuron maturation and uncovered developmental programs within their gene expression profiles. Finally, we identified that distinct neuronal subtypes are differentially affected by sensory experience. Together, these data provide a transcriptome-based foundation for investigating subtype-specific neuronal function in the olfactory bulb (OB), charting the molecular profiles that arise during the maturation and integration of adult-born neurons and how they dynamically change in an activity-dependent manner., (Copyright © 2018 The Author(s). Published by Elsevier Inc. All rights reserved.)

Loss of the paired-like homeodomain transcription factor 2 ( Pitx2 ) in cardiomyocytes predisposes mice to atrial fibrillation and compromises neonatal regenerative capacity. In addition, Pitx2 gain-of-function protects mature cardiomyocytes from ischemic injury and promotes heart repair. Here, we characterized the long-term myocardial phenotype following myocardial infarction (MI) in Pitx2 conditional-knockout ( Pitx2 CKO) mice. We found adipose-like tissue in Pitx2 CKO hearts 60 days after MI induced surgically at postnatal day 2 but not at day 8. Molecular and cellular analyses showed the onset of adipogenic signaling in mutant hearts after MI. Lineage tracing experiments showed a non-cardiomyocyte origin of the de novo adipose-like tissue. Interestingly, we found that Pitx2 promotes mitochondrial function through its gene regulatory network, and that the knockdown of a key mitochondrial Pitx2 target gene, Cox7c , also leads to the accumulation of myocardial fat tissue. Single-nuclei RNA-seq revealed that Pitx2 -deficient hearts were oxidatively stressed. Our findings reveal a role for Pitx2 in maintaining proper cardiac cellular composition during heart regeneration via the maintenance of proper mitochondrial structure and function., Competing Interests: Competing interestsThe authors declare no competing or financial interests., (© 2018. Published by The Company of Biologists Ltd.)

Decoding the molecular composition of individual Ngn3 + endocrine progenitors (EPs) during pancreatic morphogenesis could provide insight into the mechanisms regulating hormonal cell fate. Here, we identify population markers and extensive cellular diversity including four EP subtypes reflecting EP maturation using high-resolution single-cell RNA-sequencing of the e14.5 and e16.5 mouse pancreas. While e14.5 and e16.5 EPs are constantly born and share select genes, these EPs are overall transcriptionally distinct concomitant with changes in the underlying epithelium. As a consequence, e16.5 EPs are not the same as e14.5 EPs: e16.5 EPs have a higher propensity to form beta cells. Analysis of e14.5 and e16.5 EP chromatin states reveals temporal shifts, with enrichment of beta cell motifs in accessible regions at later stages. Finally, we provide transcriptional maps outlining the route progenitors take as they make cell fate decisions, which can be applied to advance the in vitro generation of beta cells.

During development, progenitors progress through transition states. The cardiac epicardium contains progenitors of essential non-cardiomyocytes. The Hippo pathway, a kinase cascade that inhibits the Yap transcriptional co-factor, controls organ size in developing hearts. Here, we investigated Hippo kinases Lats1 and Lats2 in epicardial diversification. Epicardial-specific deletion of Lats1/2 was embryonic lethal, and mutant embryos had defective coronary vasculature remodeling. Single-cell RNA sequencing revealed that Lats1/2 mutant cells failed to activate fibroblast differentiation but remained in an intermediate cell state with both epicardial and fibroblast characteristics. Lats1/2 mutant cells displayed an arrested developmental trajectory with persistence of epicardial markers and expanded expression of Yap targets Dhrs3, an inhibitor of retinoic acid synthesis, and Dpp4, a protease that modulates extracellular matrix (ECM) composition. Genetic and pharmacologic manipulation revealed that Yap inhibits fibroblast differentiation, prolonging a subepicardial-like cell state, and promotes expression of matricellular factors, such as Dpp4, that define ECM characteristics., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Mammalian organs vary widely in regenerative capacity. Poorly regenerative organs, such as the heart are particularly vulnerable to organ failure. Once established, heart failure commonly results in mortality. The Hippo pathway, a kinase cascade that prevents adult cardiomyocyte proliferation and regeneration, is upregulated in human heart failure. Here we show that deletion of the Hippo pathway component Salvador (Salv) in mouse hearts with established ischaemic heart failure after myocardial infarction induces a reparative genetic program with increased scar border vascularity, reduced fibrosis, and recovery of pumping function compared with controls. Using translating ribosomal affinity purification, we isolate cardiomyocyte-specific translating messenger RNA. Hippo-deficient cardiomyocytes have increased expression of proliferative genes and stress response genes, such as the mitochondrial quality control gene, Park2. Genetic studies indicate that Park2 is essential for heart repair, suggesting a requirement for mitochondrial quality control in regenerating myocardium. Gene therapy with a virus encoding Salv short hairpin RNA improves heart function when delivered at the time of infarct or after ischaemic heart failure following myocardial infarction was established. Our findings indicate that the failing heart has a previously unrecognized reparative capacity involving more than cardiomyocyte renewal.

Diseases of lipid metabolism are a major cause of human morbidity, but no animal model entirely recapitulates human lipoprotein metabolism. Here we develop a xenograft mouse model using hepatocytes from a patient with familial hypercholesterolaemia caused by loss-of-function mutations in the low-density lipoprotein receptor (LDLR). Like familial hypercholesterolaemia patients, our familial hypercholesterolaemia liver chimeric mice develop hypercholesterolaemia and a 'humanized' serum profile, including expression of the emerging drug targets cholesteryl ester transfer protein and apolipoprotein (a), for which no genes exist in mice. We go on to replace the missing LDLR in familial hypercholesterolaemia liver chimeric mice using an adeno-associated virus 9-based gene therapy and restore normal lipoprotein profiles after administration of a single dose. Our study marks the first time a human metabolic disease is induced in an experimental animal model by human hepatocyte transplantation and treated by gene therapy. Such xenograft platforms offer the ability to validate human experimental therapies and may foster their rapid translation into the clinic.