Hand2 Selectively Reorganizes Chromatin Accessibility to Induce Pacemaker-like Transcriptional Reprogramming

Active repression of cell fate plasticity by PROX1 safeguards hepatocyte identity and prevents liver tumorigenesis

Cell fate plasticity enables development, yet unlocked plasticity is a cancer hallmark. While transcription master regulators induce lineage-specific genes to restrict plasticity, it remains unclear whether plasticity is actively suppressed by lineage-specific repressors. Here we computationally predict so-called safeguard repressors for 18 cell types that block phenotypic plasticity lifelong. We validated hepatocyte-specific candidates using reprogramming, revealing that prospero homeobox protein 1 (PROX1) enhanced hepatocyte identity by direct repression of alternative fate master regulators. In mice, Prox1 was required for efficient hepatocyte regeneration after injury and was sufficient to prevent liver tumorigenesis. In line with patient data, Prox1 depletion caused hepatocyte fate loss in vivo and enabled the transition of hepatocellular carcinoma to cholangiocarcinoma. Conversely, overexpression promoted cholangiocarcinoma to hepatocellular carcinoma transdifferentiation. Our findings provide evidence for PROX1 as a hepatocyte-specific safeguard and support a model where cell-type-specific repressors actively suppress plasticity throughout life to safeguard lineage identity and thus prevent disease.

Single-cell sequencing delineates T-cell clonality and pathogenesis of the parapsoriasis disease group

BACKGROUND

Mycosis fungoides (MF), the most common cutaneous T-cell lymphoma, is often underdiagnosed in early stages because of similarities with benign dermatoses such as atopic dermatitis (AD). Furthermore, the delineation from what is called "parapsoriasis en plaque", a disease that can appear either in a small- or large-plaque form, is still controversial.

OBJECTIVE

We sought to characterize the parapsoriasis disease spectrum.

METHODS

We performed single-cell RNA sequencing of skin biopsies from patients within the parapsoriasis-to-early-stage MF spectrum, stratified for small and large plaques, and compared them to AD, psoriasis, and healthy control skin.

RESULTS

Six of 8 large-plaque lesions harbored either an expanded alpha/beta or gamma/delta T-cell clone with downregulation of CD7 expression, consistent with a diagnosis of early-stage MF. In contrast, 6 of 7 small-plaque lesions were polyclonal in nature, thereby lacking a lymphomatous phenotype, and also revealed a less inflammatory microenvironment than early-stage MF or AD. Of note, polyclonal small- and large-plaque lesions characteristically harbored a population of NPY+ innate lymphoid cells and displayed a stromal signature of complement upregulation and antimicrobial hyperresponsiveness in fibroblasts and sweat gland cells, respectively. These conditions were clearly distinct from AD or psoriasis, which uniquely harbored CD3+CRTH2+ IL-13 expressing "TH2A" cells, or strong type 17 inflammation, respectively.

CONCLUSION

These data position polyclonal small- and large-plaque parapsoriasis lesions as a separate disease entity that characteristically harbors a so far undescribed innate lymphoid cell population. We thus propose a new term, "polyclonal parapsoriasis en plaque", for this kind of lesion because they can be clearly differentiated from early- and advanced-stage MF, psoriasis, and AD on several cellular and molecular levels.

Chromatin accessibility: biological functions, molecular mechanisms and therapeutic application

The dynamic regulation of chromatin accessibility is one of the prominent characteristics of eukaryotic genome. The inaccessible regions are mainly located in heterochromatin, which is multilevel compressed and access restricted. The remaining accessible loci are generally located in the euchromatin, which have less nucleosome occupancy and higher regulatory activity. The opening of chromatin is the most important prerequisite for DNA transcription, replication, and damage repair, which is regulated by genetic, epigenetic, environmental, and other factors, playing a vital role in multiple biological progresses. Currently, based on the susceptibility difference of occupied or free DNA to enzymatic cleavage, solubility, methylation, and transposition, there are many methods to detect chromatin accessibility both in bulk and single-cell level. Through combining with high-throughput sequencing, the genome-wide chromatin accessibility landscape of many tissues and cells types also have been constructed. The chromatin accessibility feature is distinct in different tissues and biological states. Research on the regulation network of chromatin accessibility is crucial for uncovering the secret of various biological processes. In this review, we comprehensively introduced the major functions and mechanisms of chromatin accessibility variation in different physiological and pathological processes, meanwhile, the targeted therapies based on chromatin dynamics regulation are also summarized.

A gene desert required for regulatory control of pleiotropic Shox2 expression and embryonic survival

Approximately a quarter of the human genome consists of gene deserts, large regions devoid of genes often located adjacent to developmental genes and thought to contribute to their regulation. However, defining the regulatory functions embedded within these deserts is challenging due to their large size. Here, we explore the cis-regulatory architecture of a gene desert flanking the Shox2 gene, which encodes a transcription factor indispensable for proximal limb, craniofacial, and cardiac pacemaker development. We identify the gene desert as a regulatory hub containing more than 15 distinct enhancers recapitulating anatomical subdomains of Shox2 expression. Ablation of the gene desert leads to embryonic lethality due to Shox2 depletion in the cardiac sinus venosus, caused in part by the loss of a specific distal enhancer. The gene desert is also required for stylopod morphogenesis, mediated via distributed proximal limb enhancers. In summary, our study establishes a multi-layered role of the Shox2 gene desert in orchestrating pleiotropic developmental expression through modular arrangement and coordinated dynamics of tissue-specific enhancers.

Cadherin-5 facilitated the differentiation of human induced pluripotent stem cells into sinoatrial node-like pacemaker cells by regulating β-catenin

Our study was conducted to investigate whether cadherin-5 (CDH5), a vascular endothelial cell adhesion glycoprotein, could facilitate the differentiation of human induced pluripotent stem cells (hiPSCs) into sinoatrial node-like pacemaker cells (SANLPCs), following previous findings of silk-fibroin hydrogel-induced direct conversion of quiescent cardiomyocytes into pacemaker cells in rats through the activation of CDH5. In this study, the differentiating hiPSCs were treated with CDH5 (40 ng/mL) between Day 5 and 7 during cardiomyocytes differentiation. The findings in the present study demonstrated that CDH5 stimulated the expression of pacemaker-specific markers while suppressing markers associated with working cardiomyocytes, resulting in an increased proportion of SANLPCs among hiPSCs-derived cardiomyocytes (hiPSC-CMs) population. Moreover, CDH5 induced typical electrophysiological characteristics resembling cardiac pacemaker cells in hiPSC-CMs. Further mechanistic investigations revealed that the enriched differentiation of hiPSCs into SANLPCs induced by CDH5 was partially reversed by iCRT14, an inhibitor of β-catenin. Therefore, based on the aforementioned findings, it could be inferred that the regulation of β-catenin by CDH5 played a crucial role in promoting the enriched differentiation of hiPSCs into SANLPCs, which presents a novel avenue for the construction of biological pacemakers in forthcoming research.

Can we stop one heart from breaking: triumphs and challenges in cardiac reprogramming

Ischemic cardiac injury causes irreversible muscle loss and scarring, but recent years have seen dramatic advances in cardiac reprogramming, the field focused on regenerating cardiac muscle. With SARS-CoV2 increasing the age-adjusted cardiovascular disease mortality rate, it is worth evaluating the state of this field. Here, we summarize novel innovations in reprogramming strategies, insights into their mechanisms, and technologies for factor delivery. We also propose a broad model of reprogramming to suggest directions for future research. Poet Emily Dickinson wrote, "If I can stop one heart from breaking, I shall not live in vain." Today, researchers studying cardiac reprogramming view this line as a call to action to translate this revolutionary approach into life-saving treatments for patients with cardiovascular diseases.

Hand2 represses non-cardiac cell fates through chromatin remodeling at cis- regulatory elements

Developmental studies have revealed the importance of the transcription factor Hand2 in cardiac development. Hand2 promotes cardiac progenitor differentiation and epithelial maturation, while repressing other tissue types. The mechanisms underlying the promotion of cardiac fates are far better understood than those underlying the repression of alternative fates. Here, we assess Hand2-dependent changes in gene expression and chromatin remodeling in cardiac progenitors of zebrafish embryos. Cell-type specific transcriptome analysis shows a dual function for Hand2 in activation of cardiac differentiation genes and repression of pronephric pathways. We identify functional cis- regulatory elements whose chromatin accessibility are increased in hand2 mutant cells. These regulatory elements associate with non-cardiac gene expression, and drive reporter gene expression in tissues associated with Hand2-repressed genes. We find that functional Hand2 is sufficient to reduce non-cardiac reporter expression in cardiac lineages. Taken together, our data support a model of Hand2-dependent coordination of transcriptional programs, not only through transcriptional activation of cardiac and epithelial maturation genes, but also through repressive chromatin remodeling at the DNA regulatory elements of non-cardiac genes.

Global chromatin landscapes identify candidate noncoding modifiers of cardiac rhythm

Comprehensive cis-regulatory landscapes are essential for accurate enhancer prediction and disease variant mapping. Although cis-regulatory element (CRE) resources exist for most tissues and organs, many rare - yet functionally important - cell types remain overlooked. Despite representing only a small fraction of the heart's cellular biomass, the cardiac conduction system (CCS) unfailingly coordinates every life-sustaining heartbeat. To globally profile the mouse CCS cis-regulatory landscape, we genetically tagged CCS component-specific nuclei for comprehensive assay for transposase-accessible chromatin-sequencing (ATAC-Seq) analysis. Thus, we established a global CCS-enriched CRE database, referred to as CCS-ATAC, as a key resource for studying CCS-wide and component-specific regulatory functions. Using transcription factor (TF) motifs to construct CCS component-specific gene regulatory networks (GRNs), we identified and independently confirmed several specific TF sub-networks. Highlighting the functional importance of CCS-ATAC, we also validated numerous CCS-enriched enhancer elements and suggested gene targets based on CCS single-cell RNA-Seq data. Furthermore, we leveraged CCS-ATAC to improve annotation of existing human variants related to cardiac rhythm and nominated a potential enhancer-target pair that was dysregulated by a specific SNP. Collectively, our results established a CCS-regulatory compendium, identified novel CCS enhancer elements, and illuminated potential functional associations between human genomic variants and CCS component-specific CREs.

An atrial fibrillation-associated regulatory region modulates cardiac Tbx5 levels and arrhythmia susceptibility

Heart development and rhythm control are highly Tbx5 dosage-sensitive. TBX5 haploinsufficiency causes congenital conduction disorders, whereas increased expression levels of TBX5 in human heart samples has been associated with atrial fibrillation (AF). We deleted the conserved mouse orthologues of two independent AF-associated genomic regions in the Tbx5 locus, one intronic (RE(int)) and one downstream (RE(down)) of Tbx5. In both lines, we observed a modest (30%) increase of Tbx5 in the postnatal atria. To gain insight into the effects of slight dosage increase in vivo, we investigated the atrial transcriptional, epigenetic and electrophysiological properties of both lines. Increased atrial Tbx5 expression was associated with induction of genes involved in development, ion transport and conduction, with increased susceptibility to atrial arrhythmias, and increased action potential duration of atrial cardiomyocytes. We identified an AF-associated variant in the human RE(int) that increases its transcriptional activity. Expression of the AF-associated transcription factor Prrx1 was induced in Tbx5^RE(int)KO cardiomyocytes. We found that some of the transcriptional and functional changes in the atria caused by increased Tbx5 expression were normalized when reducing cardiac Prrx1 expression in Tbx5^RE(int)KO mice, indicating an interaction between these two AF genes. We conclude that modest increases in expression of dose-dependent transcription factors, caused by common regulatory variants, significantly impact on the cardiac gene regulatory network and disease susceptibility.

The HASTER lncRNA promoter is a cis-acting transcriptional stabilizer of HNF1A

The biological purpose of long non-coding RNAs (lncRNAs) is poorly understood. Haploinsufficient mutations in HNF1A homeobox A (HNF1A), encoding a homeodomain transcription factor, cause diabetes mellitus. Here, we examine HASTER, the promoter of an lncRNA antisense to HNF1A. Using mouse and human models, we show that HASTER maintains cell-specific physiological HNF1A concentrations through positive and negative feedback loops. Pancreatic β cells from Haster mutant mice consequently showed variegated HNF1A silencing or overexpression, resulting in hyperglycaemia. HASTER-dependent negative feedback was essential to prevent HNF1A binding to inappropriate genomic regions. We demonstrate that the HASTER promoter DNA, rather than the lncRNA, modulates HNF1A promoter-enhancer interactions in cis and thereby regulates HNF1A transcription. Our studies expose a cis-regulatory element that is unlike classic enhancers or silencers, it stabilizes the transcription of its target gene and ensures the fidelity of a cell-specific transcription factor program. They also show that disruption of a mammalian lncRNA promoter can cause diabetes mellitus.

Generation of a recombinant version of a biologically active cell-permeant human HAND2 transcription factor from E. coli

Transcription factor HAND2 has a significant role in vascularization, angiogenesis, and cardiac neural crest development. It is one of the key cardiac factors crucial for the enhanced derivation of functional and mature myocytes from non-myocyte cells. Here, we report the generation of the recombinant human HAND2 fusion protein from the heterologous system. First, we cloned the full-length human HAND2 gene (only protein-coding sequence) after codon optimization along with the fusion tags (for cell penetration, nuclear translocation, and affinity purification) into the expression vector. We then transformed and expressed it in Escherichia coli strain, BL21(DE3). Next, the effect (in terms of expression) of tagging fusion tags with this recombinant protein at two different terminals was also investigated. Using affinity chromatography, we established the one-step homogeneous purification of recombinant human HAND2 fusion protein; and through circular dichroism spectroscopy, we established that this purified protein had retained its secondary structure. We then showed that this purified human protein could transduce the human cells and translocate to its nucleus. The generated recombinant HAND2 fusion protein showed angiogenic potential in the ex vivo chicken embryo model. Following transduction in MEF2C overexpressing cardiomyoblast cells, this purified recombinant protein synergistically activated the α-MHC promoter and induced GFP expression in the α-MHC-eGFP reporter assay. Prospectively, the purified bioactive recombinant HAND2 protein can potentially be a safe and effective molecular tool in the direct cardiac reprogramming process and other biological applications.

A review of protocols for human iPSC culture, cardiac differentiation, subtype-specification, maturation, and direct reprogramming

The methods for the culture and cardiomyocyte differentiation of human embryonic stem cells, and later human induced pluripotent stem cells (hiPSC), have moved from a complex and uncontrolled systems to simplified and relatively robust protocols, using the knowledge and cues gathered at each step. HiPSC-derived cardiomyocytes have proven to be a useful tool in human disease modelling, drug discovery, developmental biology, and regenerative medicine. In this protocol review, we will highlight the evolution of protocols associated with hPSC culture, cardiomyocyte differentiation, sub-type specification, and cardiomyocyte maturation. We also discuss protocols for somatic cell direct reprogramming to cardiomyocyte-like cells.

Neonatal Deletion of Hand1 and Hand2 within Murine Cardiac Conduction System Reveals a Novel Role for HAND2 in Rhythm Homeostasis

The cardiac conduction system, a network of specialized cells, is required for the functioning of the heart. The basic helix loop helix factors Hand1 and Hand2 are required for cardiac morphogenesis and have been implicated in cardiac conduction system development and maintenance. Here we use embryonic and post-natal specific Cre lines to interrogate the role of Hand1 and Hand2 in the function of the murine cardiac conduction system. Results demonstrate that loss of HAND1 in the post-natal conduction system does not result in any change in electrocardiogram parameters or within the ventricular conduction system as determined by optical voltage mapping. Deletion of Hand2 within the post-natal conduction system results in sex-dependent reduction in PR interval duration in these mice, suggesting a novel role for HAND2 in regulating the atrioventricular conduction. Surprisingly, results show that loss of both HAND factors within the post-natal conduction system does not cause any consistent changes in cardiac conduction system function. Deletion of Hand2 in the embryonic left ventricle results in inconsistent prolongation of PR interval and susceptibility to atrial arrhythmias. Thus, these results suggest a novel role for HAND2 in homeostasis of the murine cardiac conduction system and that HAND1 loss potentially rescues the shortened HAND2 PR phenotype.

Patient-specific TBX5-G125R Variant Induces Profound Transcriptional Deregulation and Atrial Dysfunction

Background: The pathogenic missense variant p.G125R in TBX5 causes Holt-Oram syndrome (HOS; hand-heart syndrome) and early onset of atrial fibrillation. Revealing how an altered key developmental transcription factor modulates cardiac physiology in vivo will provide unique insights into the mechanisms underlying atrial fibrillation in these patients. Methods: We analyzed electrocardiograms (ECGs) of an extended family pedigree of HOS patients. Next, we introduced the TBX5-p.G125R variant in the mouse genome (Tbx5^G125R) and performed electrophysiological analyses (ECG, optical mapping, patch clamp, intracellular calcium measurements), transcriptomics (single nuclei and tissue RNA sequencing) and epigenetic profiling (ATAC-sequencing, H3K27ac CUT&RUN-sequencing). Results: We discovered high incidence of atrial extra systoles and atrioventricular conduction disturbances in HOS patients. Tbx5^G125R/+ mice were morphologically unaffected and displayed variable RR intervals, atrial extra systoles and susceptibility to atrial fibrillation, reminiscent of TBX5-p.G125R patients. Atrial conduction velocity was not affected but systolic and diastolic intracellular calcium concentrations were decreased and action potentials prolonged in isolated cardiomyocytes of Tbx5^G125R/+ mice compared to controls. Transcriptional profiling of atria revealed most profound transcriptional changes in cardiomyocytes versus other cell types, and identified over a thousand coding and non-coding transcripts that were differentially expressed. Epigenetic profiling uncovered thousands of TBX5-p.G125R sensitive putative regulatory elements (including enhancers) that gained accessibility in atrial cardiomyocytes. The majority of sites with increased accessibility were occupied by Tbx5. The small group of sites with reduced accessibility was enriched for DNA binding motifs of members of the SP- and KLF families of transcription factors. These data show that Tbx5-p.G125R induces changes in regulatory element activity, altered transcriptional regulation and changed cardiomyocyte behavior, possibly caused by altered DNA binding and cooperativity properties. Conclusions: Our data reveal how a disease-causing missense variant in TBX5 induces profound changes in the atrial transcriptional regulatory network and epigenetic state in vivo, leading to arrhythmia reminiscent of those seen in human TBX5-p.G125R variant carriers.

Improving Cardiac Reprogramming for Heart Regeneration in Translational Medicine

Direct reprogramming of fibroblasts into CM-like cells has emerged as an attractive strategy to generate induced CMs (iCMs) in heart regeneration. However, low conversion rate, poor purity, and the lack of precise conversion of iCMs are still present as significant challenges. In this review, we summarize the recent development in understanding the molecular mechanisms of cardiac reprogramming with various strategies to achieve more efficient iCMs. reprogramming. Specifically, we focus on the identified critical roles of transcriptional regulation, epigenetic modification, signaling pathways from the cellular microenvironment, and cell cycling regulation in cardiac reprogramming. We also discuss the progress in delivery system optimization and cardiac reprogramming in human cells related to preclinical applications. We anticipate that this will translate cardiac reprogramming-based heart therapy into clinical applications. In addition to optimizing the cardiogenesis related transcriptional regulation and signaling pathways, an important strategy is to modulate the pathological microenvironment associated with heart injury, including inflammation, pro-fibrotic signaling pathways, and the mechanical properties of the damaged myocardium. We are optimistic that cardiac reprogramming will provide a powerful therapy in heart regenerative medicine.

Transcriptional and Epigenetic Landscape of Cardiac Pacemaker Cells: Insights Into Cellular Specialization in the Sinoatrial Node

Cardiac pacemaker cells differentiate and functionally specialize early in embryonic development through activation of critical gene regulatory networks. In general, cellular specification and differentiation require that combinations of cell type-specific transcriptional regulators activate expression of key effector genes by binding to DNA regulatory elements including enhancers and promoters. However, because genomic DNA is tightly packaged by histones that must be covalently modified in order to render DNA regulatory elements and promoters accessible for transcription, the process of development and differentiation is intimately connected to the epigenetic regulation of chromatin accessibility. Although the difficulty of obtaining sufficient quantities of pure populations of pacemaker cells has limited progress in this field, the advent of low-input genomic technologies has the potential to catalyze a rapid growth of knowledge in this important area. The goal of this review is to outline the key transcriptional networks that control pacemaker cell development, with particular attention to our emerging understanding of how chromatin accessibility is modified and regulated during pacemaker cell differentiation. In addition, we will discuss the relevance of these findings to adult sinus node function, sinus node diseases, and origins of genetic variation in heart rhythm. Lastly, we will outline the current challenges facing this field and promising directions for future investigation.

Molecular Profiling of the Cardiac Conduction System: the Dawn of a New Era

PURPOSE OF REVIEW

Recent technological advances have led to an increased ability to define the gene expression profile of the cardiac conduction system (CCS). Here, we review the most salient studies to emerge in recent years and discuss existing gaps in our knowledge as well as future areas of investigation.

RECENT FINDINGS

Molecular profiling of the CCS spans several decades. However, the advent of high-throughput sequencing strategies has allowed for the discovery of unique transcriptional programs of the many diverse CCS cell types. The CCS, a diverse structure with significant inter- and intra-component cellular heterogeneity, is essential to the normal function of the heart. Progress in transcriptomic profiling has improved the resolution and depth of characterization of these unique and clinically relevant CCS cell types. Future studies leveraging this big data will play a crucial role in improving our understanding of CCS development and function as well as translating these findings into tangible translational tools for the improved detection, prevention, and treatment of cardiac arrhythmias.

A Variant Noncoding Region Regulates Prrx1 and Predisposes to Atrial Arrhythmias

[Figure: see text].

Variant Intronic Enhancer Controls SCN10A-short Expression and Heart Conduction

BACKGROUND

Genetic variants in SCN10A, encoding the neuronal voltage-gated sodium channel Na_V1.8, are strongly associated with atrial fibrillation, Brugada syndrome, cardiac conduction velocities, and heart rate. The cardiac function of SCN10A has not been resolved, however, and diverging mechanisms have been proposed. Here, we investigated the cardiac expression of SCN10A and the function of a variant-sensitive intronic enhancer previously linked to the regulation of SCN5A, encoding the major essential cardiac sodium channel Na_V1.5.

METHODS

The expression of SCN10A was investigated in mouse and human hearts. With the use of CRISPR/Cas9 genome editing, the mouse intronic enhancer was disrupted, and mutant mice were characterized by transcriptomic and electrophysiological analyses. The association of genetic variants at SCN5A-SCN10A enhancer regions and gene expression were evaluated by genome-wide association studies single-nucleotide polymorphism mapping and expression quantitative trait loci analysis.

RESULTS

We found that cardiomyocytes of the atria, sinoatrial node, and ventricular conduction system express a short transcript comprising the last 7 exons of the gene (Scn10a-short). Transcription occurs from an intronic enhancer-promoter complex, whereas full-length Scn10a transcript was undetectable in the human and mouse heart. Expression quantitative trait loci analysis revealed that the genetic variants in linkage disequilibrium with genetic variant rs6801957 in the intronic enhancer associate with SCN10A transcript levels in the heart. Genetic modification of the enhancer in the mouse genome led to reduced cardiac Scn10a-short expression in atria and ventricles, reduced cardiac sodium current in atrial cardiomyocytes, atrial conduction slowing and arrhythmia, whereas the expression of Scn5a, the presumed enhancer target gene, remained unaffected. In patch-clamp transfection experiments, expression of Scn10a-short-encoded Na_V1.8-short increased Na_V1.5-mediated sodium current. We propose that noncoding genetic variation modulates transcriptional regulation of Scn10a-short in cardiomyocytes that impacts Na_V1.5-mediated sodium current and heart rhythm.

CONCLUSIONS

Genetic variants in and around SCN10A modulate enhancer function and expression of a cardiac-specific SCN10A-short transcript. We propose that noncoding genetic variation modulates transcriptional regulation of a functional C-terminal portion of Na_V1.8 in cardiomyocytes that impacts on Na_V1.5 function, cardiac conduction velocities, and arrhythmia susceptibility.

Genome-Wide Analysis Identifies an Essential Human TBX3 Pacemaker Enhancer

RATIONALE

The development and function of the pacemaker cardiomyocytes of the sinoatrial node (SAN), the leading pacemaker of the heart, are tightly controlled by a conserved network of transcription factors, including TBX3 (T-box transcription factor 3), ISL1 (ISL LIM homeobox 1), and SHOX2 (short stature homeobox 2). Yet, the regulatory DNA elements (REs) controlling target gene expression in the SAN pacemaker cells have remained undefined.

OBJECTIVE

Identification of the regulatory landscape of human SAN-like pacemaker cells and functional assessment of SAN-specific REs potentially involved in pacemaker cell gene regulation.

METHODS AND RESULTS

We performed Assay for Transposase-Accessible Chromatin using sequencing on human pluripotent stem cell-derived SAN-like pacemaker cells and ventricle-like cells and identified thousands of putative REs specific for either human cell type. We validated pacemaker cell-specific elements in the SHOX2 and TBX3 loci. CRISPR-mediated homozygous deletion of the mouse ortholog of a noncoding region with candidate pacemaker-specific REs in the SHOX2 locus resulted in selective loss of Shox2 expression from the developing SAN and embryonic lethality. Putative pacemaker-specific REs were identified up to 1 Mbp upstream of TBX3 in a region close to MED13L harboring variants associated with heart rate recovery after exercise. The orthologous region was deleted in mice, which resulted in selective loss of expression of Tbx3 from the SAN and (cardiac) ganglia and in neonatal lethality. Expression of Tbx3 was maintained in other tissues including the atrioventricular conduction system, lungs, and liver. Heterozygous adult mice showed increased SAN recovery times after pacing. The human REs harboring the associated variants robustly drove expression in the SAN of transgenic mouse embryos.

CONCLUSIONS

We provided a genome-wide collection of candidate human pacemaker-specific REs, including the loci of SHOX2, TBX3, and ISL1, and identified a link between human genetic variants influencing heart rate recovery after exercise and a variant RE with highly conserved function, driving SAN expression of TBX3.

The genetic structure and adaptation of Andean highlanders and Amazonians are influenced by the interplay between geography and culture

Western South America was one of the worldwide cradles of civilization. The well-known Inca Empire was the tip of the iceberg of an evolutionary process that started 11,000 to 14,000 years ago. Genetic data from 18 Peruvian populations reveal the following: 1) The between-population homogenization of the central southern Andes and its differentiation with respect to Amazonian populations of similar latitudes do not extend northward. Instead, longitudinal gene flow between the northern coast of Peru, Andes, and Amazonia accompanied cultural and socioeconomic interactions revealed by archeology. This pattern recapitulates the environmental and cultural differentiation between the fertile north, where altitudes are lower, and the arid south, where the Andes are higher, acting as a genetic barrier between the sharply different environments of the Andes and Amazonia. 2) The genetic homogenization between the populations of the arid Andes is not only due to migrations during the Inca Empire or the subsequent colonial period. It started at least during the earlier expansion of the Wari Empire (600 to 1,000 years before present). 3) This demographic history allowed for cases of positive natural selection in the high and arid Andes vs. the low Amazon tropical forest: in the Andes, a putative enhancer in HAND2-AS1 (heart and neural crest derivatives expressed 2 antisense RNA1, a noncoding gene related to cardiovascular function) and rs269868-C/Ser1067 in DUOX2 (dual oxidase 2, related to thyroid function and innate immunity) genes and, in the Amazon, the gene encoding for the CD45 protein, essential for antigen recognition by T and B lymphocytes in viral-host interaction.

Development of the Cardiac Conduction System

The cardiac conduction system initiates and propagates each heartbeat. Specialized conducting cells are a well-conserved phenomenon across vertebrate evolution, although mammalian and avian species harbor specific components unique to organisms with four-chamber hearts. Early histological studies in mammals provided evidence for a dominant pacemaker within the right atrium and clarified the existence of the specialized muscular axis responsible for atrioventricular conduction. Building on these seminal observations, contemporary genetic techniques in a multitude of model organisms has characterized the developmental ontogeny, gene regulatory networks, and functional importance of individual anatomical compartments within the cardiac conduction system. This review describes in detail the transcriptional and regulatory networks that act during cardiac conduction system development and homeostasis with a particular emphasis on networks implicated in human electrical variation by large genome-wide association studies. We conclude with a discussion of the clinical implications of these studies and describe some future directions.

ATAC-Seq Reveals an Isl1 Enhancer That Regulates Sinoatrial Node Development and Function

RATIONALE

Cardiac pacemaker cells (PCs) in the sinoatrial node (SAN) have a distinct gene expression program that allows them to fire automatically and initiate the heartbeat. Although critical SAN transcription factors, including Isl1 (Islet-1), Tbx3 (T-box transcription factor 3), and Shox2 (short-stature homeobox protein 2), have been identified, the cis-regulatory architecture that governs PC-specific gene expression is not understood, and discrete enhancers required for gene regulation in the SAN have not been identified.

OBJECTIVE

To define the epigenetic profile of PCs using comparative ATAC-seq (assay for transposase-accessible chromatin with sequencing) and to identify novel enhancers involved in SAN gene regulation, development, and function.

METHODS AND RESULTS

We used ATAC-seq on sorted neonatal mouse SAN to compare regions of accessible chromatin in PCs and right atrial cardiomyocytes. PC-enriched assay for transposase-accessible chromatin peaks, representing candidate SAN regulatory elements, were located near established SAN genes and were enriched for distinct sets of TF (transcription factor) binding sites. Among several novel SAN enhancers that were experimentally validated using transgenic mice, we identified a 2.9-kb regulatory element at the Isl1 locus that was active specifically in the cardiac inflow at embryonic day 8.5 and throughout later SAN development and maturation. Deletion of this enhancer from the genome of mice resulted in SAN hypoplasia and sinus arrhythmias. The mouse SAN enhancer also directed reporter activity to the inflow tract in developing zebrafish hearts, demonstrating deep conservation of its upstream regulatory network. Finally, single nucleotide polymorphisms in the human genome that occur near the region syntenic to the mouse enhancer exhibit significant associations with resting heart rate in human populations.

CONCLUSIONS

(1) PCs have distinct regions of accessible chromatin that correlate with their gene expression profile and contain novel SAN enhancers, (2) cis-regulation of Isl1 specifically in the SAN depends upon a conserved SAN enhancer that regulates PC development and SAN function, and (3) a corresponding human ISL1 enhancer may regulate human SAN function.