Transcription factor protein interactomes reveal genetic determinants in heart disease.

Congenital heart disease (CHD) is present in 1% of live births, yet identification of causal mutations remains challenging. We hypothesized that genetic determinants for CHDs may lie in the protein interactomes of transcription factors whose mutations cause CHDs. Defining the interactomes of two transcription factors haplo-insufficient in CHD, GATA4 and TBX5, within human cardiac progenitors, and integrating the results with nearly 9,000 exomes from proband-parent trios revealed an enrichment of de novo missense variants associated with CHD within the interactomes. Scoring variants of interactome members based on residue, gene, and proband features identified likely CHD-causing genes, including the epigenetic reader GLYR1. GLYR1 and GATA4 widely co-occupied and co-activated cardiac developmental genes, and the identified GLYR1 missense variant disrupted interaction with GATA4, impairing in vitro and in vivo function in mice. This integrative proteomic and genetic approach provides a framework for prioritizing and interrogating genetic variants in heart disease.

Integrative Analyses of Human Reprogramming Reveal Dynamic Nature of Induced Pluripotency.

Induced pluripotency is a promising avenue for disease modeling and therapy, but the molecular principles underlying this process, particularly in human cells, remain poorly understood due to donor-to-donor variability and intercellular heterogeneity. Here, we constructed and characterized a clonal, inducible human reprogramming system that provides a reliable source of cells at any stage of the process. This system enabled integrative transcriptional and epigenomic analysis across the human reprogramming timeline at high resolution. We observed distinct waves of gene network activation, including the ordered re-activation of broad developmental regulators followed by early embryonic patterning genes and culminating in the emergence of a signature reminiscent of pre-implantation stages. Moreover, complementary functional analyses allowed us to identify and validate novel regulators of the reprogramming process. Altogether, this study sheds light on the molecular underpinnings of induced pluripotency in human cells and provides a robust cell platform for further studies. PAPERCLIP.

Transcriptional and epigenetic dynamics during specification of human embryonic stem cells.

Differentiation of human embryonic stem cells (hESCs) provides a unique opportunity to study the regulatory mechanisms that facilitate cellular transitions in a human context. To that end, we performed comprehensive transcriptional and epigenetic profiling of populations derived through directed differentiation of hESCs representing each of the three embryonic germ layers. Integration of whole-genome bisulfite sequencing, chromatin immunoprecipitation sequencing, and RNA sequencing reveals unique events associated with specification toward each lineage. Lineage-specific dynamic alterations in DNA methylation and H3K4me1 are evident at putative distal regulatory elements that are frequently bound by pluripotency factors in the undifferentiated hESCs. In addition, we identified germ-layer-specific H3K27me3 enrichment at sites exhibiting high DNA methylation in the undifferentiated state. A better understanding of these initial specification events will facilitate identification of deficiencies in current approaches, leading to more faithful differentiation strategies as well as providing insights into the rewiring of human regulatory programs during cellular transitions.

A transcriptional switch governs fibroblast activation in heart disease.

In diseased organs, stress-activated signalling cascades alter chromatin, thereby triggering maladaptive cell state transitions. Fibroblast activation is a common stress response in tissues that worsens lung, liver, kidney and heart disease, yet its mechanistic basis remains unclear1,2. Pharmacological inhibition of bromodomain and extra-terminal domain (BET) proteins alleviates cardiac dysfunction3-7, providing a tool to interrogate and modulate cardiac cell states as a potential therapeutic approach. Here we use single-cell epigenomic analyses of hearts dynamically exposed to BET inhibitors to reveal a reversible transcriptional switch that underlies the activation of fibroblasts. Resident cardiac fibroblasts demonstrated robust toggling between the quiescent and activated state in a manner directly correlating with BET inhibitor exposure and cardiac function. Single-cell chromatin accessibility revealed previously undescribed DNA elements, the accessibility of which dynamically correlated with cardiac performance. Among the most dynamic elements was an enhancer that regulated the transcription factor MEOX1, which was specifically expressed in activated fibroblasts, occupied putative regulatory elements of a broad fibrotic gene program and was required for TGFß-induced fibroblast activation. Selective CRISPR inhibition of the single most dynamic cis-element within the enhancer blocked TGFß-induced Meox1 activation. We identify MEOX1 as a central regulator of fibroblast activation associated with cardiac dysfunction and demonstrate its upregulation after activation of human lung, liver and kidney fibroblasts. The plasticity and specificity of BET-dependent regulation of MEOX1 in tissue fibroblasts provide previously unknown trans- and cis-targets for treating fibrotic disease.

Single-cell analysis of cardiogenesis reveals basis for organ-level developmental defects.

Organogenesis involves integration of diverse cell types; dysregulation of cell-type-specific gene networks results in birth defects, which affect 5% of live births. Congenital heart defects are the most common malformations, and result from disruption of discrete subsets of cardiac progenitor cells1, but the transcriptional changes in individual progenitors that lead to organ-level defects remain unknown. Here we used single-cell RNA sequencing to interrogate early cardiac progenitor cells as they become specified during normal and abnormal cardiogenesis, revealing how dysregulation of specific cellular subpopulations has catastrophic consequences. A network-based computational method for single-cell RNA-sequencing analysis that predicts lineage-specifying transcription factors2,3 identified Hand2 as a specifier of outflow tract cells but not right ventricular cells, despite the failure of right ventricular formation in Hand2-null mice4. Temporal single-cell-transcriptome analysis of Hand2-null embryos revealed failure of outflow tract myocardium specification, whereas right ventricular myocardium was specified but failed to properly differentiate and migrate. Loss of Hand2 also led to dysregulation of retinoic acid signalling and disruption of anterior-posterior patterning of cardiac progenitors. This work reveals transcriptional determinants that specify fate and differentiation in individual cardiac progenitor cells, and exposes mechanisms of disrupted cardiac development at single-cell resolution, providing a framework for investigating congenital heart defects.

Transcription Factor GATA4 Regulates Cell Type-Specific Splicing Through Direct Interaction With RNA in Human Induced Pluripotent Stem Cell-Derived Cardiac Progenitors.

BACKGROUND: GATA4 (GATA-binding protein 4), a zinc finger-containing, DNA-binding transcription factor, is essential for normal cardiac development and homeostasis in mice and humans, and mutations in this gene have been reported in human heart defects. Defects in alternative splicing are associated with many heart diseases, yet relatively little is known about how cell type- or cell state-specific alternative splicing is achieved in the heart. Here, we show that GATA4 regulates cell type-specific splicing through direct interaction with RNA and the spliceosome in human induced pluripotent stem cell-derived cardiac progenitors. METHODS: We leveraged a combination of unbiased approaches including affinity purification of GATA4 and mass spectrometry, enhanced cross-linking with immunoprecipitation, electrophoretic mobility shift assays, in vitro splicing assays, and unbiased transcriptomic analysis to uncover GATA4's novel function as a splicing regulator in human induced pluripotent stem cell-derived cardiac progenitors. RESULTS: We found that GATA4 interacts with many members of the spliceosome complex in human induced pluripotent stem cell-derived cardiac progenitors. Enhanced cross-linking with immunoprecipitation demonstrated that GATA4 also directly binds to a large number of mRNAs through defined RNA motifs in a sequence-specific manner. In vitro splicing assays indicated that GATA4 regulates alternative splicing through direct RNA binding, resulting in functionally distinct protein products. Correspondingly, knockdown of GATA4 in human induced pluripotent stem cell-derived cardiac progenitors resulted in differential alternative splicing of genes involved in cytoskeleton organization and calcium ion import, with functional consequences associated with the protein isoforms. CONCLUSIONS: This study shows that in addition to its well described transcriptional function, GATA4 interacts with members of the spliceosome complex and regulates cell type-specific alternative splicing via sequence-specific interactions with RNA. Several genes that have splicing regulated by GATA4 have functional consequences and many are associated with dilated cardiomyopathy, suggesting a novel role for GATA4 in achieving the necessary cardiac proteome in normal and stress-responsive conditions.

Dissecting neural differentiation regulatory networks through epigenetic footprinting.

Models derived from human pluripotent stem cells that accurately recapitulate neural development in vitro and allow for the generation of specific neuronal subtypes are of major interest to the stem cell and biomedical community. Notch signalling, particularly through the Notch effector HES5, is a major pathway critical for the onset and maintenance of neural progenitor cells in the embryonic and adult nervous system. Here we report the transcriptional and epigenomic analysis of six consecutive neural progenitor cell stages derived from a HES5::eGFP reporter human embryonic stem cell line. Using this system, we aimed to model cell-fate decisions including specification, expansion and patterning during the ontogeny of cortical neural stem and progenitor cells. In order to dissect regulatory mechanisms that orchestrate the stage-specific differentiation process, we developed a computational framework to infer key regulators of each cell-state transition based on the progressive remodelling of the epigenetic landscape and then validated these through a pooled short hairpin RNA screen. We were also able to refine our previous observations on epigenetic priming at transcription factor binding sites and suggest here that they are mediated by combinations of core and stage-specific factors. Taken together, we demonstrate the utility of our system and outline a general framework, not limited to the context of the neural lineage, to dissect regulatory circuits of differentiation.

Differentiation of V2a interneurons from human pluripotent stem cells.

The spinal cord consists of multiple neuronal cell types that are critical to motor control and arise from distinct progenitor domains in the developing neural tube. Excitatory V2a interneurons in particular are an integral component of central pattern generators that control respiration and locomotion; however, the lack of a robust source of human V2a interneurons limits the ability to molecularly profile these cells and examine their therapeutic potential to treat spinal cord injury (SCI). Here, we report the directed differentiation of CHX10+ V2a interneurons from human pluripotent stem cells (hPSCs). Signaling pathways (retinoic acid, sonic hedgehog, and Notch) that pattern the neural tube were sequentially perturbed to identify an optimized combination of small molecules that yielded ∼25% CHX10+ cells in four hPSC lines. Differentiated cultures expressed much higher levels of V2a phenotypic markers (CHX10 and SOX14) than other neural lineage markers. Over time, CHX10+ cells expressed neuronal markers [neurofilament, NeuN, and vesicular glutamate transporter 2 (VGlut2)], and cultures exhibited increased action potential frequency. Single-cell RNAseq analysis confirmed CHX10+ cells within the differentiated population, which consisted primarily of neurons with some glial and neural progenitor cells. At 2 wk after transplantation into the spinal cord of mice, hPSC-derived V2a cultures survived at the site of injection, coexpressed NeuN and VGlut2, extended neurites >5 mm, and formed putative synapses with host neurons. These results provide a description of V2a interneurons differentiated from hPSCs that may be used to model central nervous system development and serve as a potential cell therapy for SCI.

Threonine phosphorylation post-translationally regulates protein secretion in Pseudomonas aeruginosa.

Secreted proteins are crucial to the arsenal of bacterial pathogens. Although optimal activity of these proteins is likely to require precise regulation of release, the signalling events that trigger secretion are poorly understood. Here, we identify a threonine phosphorylation event that post-translationally regulates the Hcp secretion island-I-encoded type VI secretion system of Pseudomonas aeruginosa (H-T6SS). We show that a serine-threonine kinase, PpkA, is required for assembly of the H-T6SS and for secretion of Hcp1. PpkA activity is antagonized by PppA, a Ser-Thr phosphatase. These proteins exhibit reciprocal effects on the H-T6SS by acting on an FHA domain-containing protein, termed Fha1. Colocalization experiments with the T6S AAA+ family protein, ClpV1, indicate that Fha1 is a core scaffolding protein of the H-T6SS. Mutations affecting this H-T6S regulatory pathway provide a molecular explanation for the variation in Hcp1 secretion among clinical P. aeruginosa isolates. This mechanism of triggering secretion may be general, as many T6SSs contain orthologues of these proteins. Post-translational regulation of protein secretion by Thr phosphorylation is unprecedented in bacteria, and is likely to reflect the requirement for T6S to respond rapidly and reversibly to its environment.

Genetic analysis and multimodal imaging identify novel mtDNA 12148T>C leading to multisystem dysfunction with tissue-specific heteroplasmy.

Patients with mitochondrial disorders present with clinically diverse symptoms, largely driven by heterogeneous mutations in mitochondrial-encoded and nuclear-encoded mitochondrial genes. These mutations ultimately lead to complex biochemical disorders with a myriad of clinical manifestations, often accumulating during childhood on into adulthood, contributing to life-altering and sometimes fatal events. It is therefore important to diagnose and characterize the associated disorders for each mitochondrial mutation as early as possible since medical management might be able to improve the quality and longevity of life in mitochondrial disease patients. Here we identify a novel mitochondrial variant in a mitochondrial transfer RNA for histidine (mt-tRNA-his) [m.12148T>C], that is associated with the development of ocular, aural, neurological, renal, and muscular dysfunctions. We provide a detailed account of a family harboring this mutation, as well as the molecular underpinnings contributing to cellular and mitochondrial dysfunction. In conclusion, this investigation provides clinical, biochemical, and morphological evidence of the pathogenicity of m.12148T>C. We highlight the importance of multiple tissue testing and in vitro disease modeling in diagnosing mitochondrial disease.

The multi-lineage transcription factor ISL1 controls cardiomyocyte cell fate through interaction with NKX2.5.

Congenital heart disease often arises from perturbations of transcription factors (TFs) that guide cardiac development. ISLET1 (ISL1) is a TF that influences early cardiac cell fate, as well as differentiation of other cell types including motor neuron progenitors (MNPs) and pancreatic islet cells. While lineage specificity of ISL1 function is likely achieved through combinatorial interactions, its essential cardiac interacting partners are unknown. By assaying ISL1 genomic occupancy in human induced pluripotent stem cell-derived cardiac progenitors (CPs) or MNPs and leveraging the deep learning approach BPNet, we identified motifs of other TFs that predicted ISL1 occupancy in each lineage, with NKX2.5 and GATA motifs being most closely associated to ISL1 in CPs. Experimentally, nearly two-thirds of ISL1-bound loci were co-occupied by NKX2.5 and/or GATA4. Removal of NKX2.5 from CPs led to widespread ISL1 redistribution, and overexpression of NKX2.5 in MNPs led to ISL1 occupancy of CP-specific loci. These results reveal how ISL1 guides lineage choices through a combinatorial code that dictates genomic occupancy and transcription.

SARS-CoV-2 Susceptibility and ACE2 Gene Variations Within Diverse Ethnic Backgrounds.

There is considerable variability in the susceptibility and progression for COVID-19 and it appears to be strongly correlated with age, gender, ethnicity and pre-existing health conditions. However, to our knowledge, cohort studies of COVID-19 in clinically vulnerable groups are lacking. Host genetics has also emerged as a major risk factor for COVID-19, and variation in the ACE2 receptor, which facilitates entry of the SARS-CoV-2 virus into the cell, has become a major focus of attention. Thus, we interrogated an ethnically diverse cohort of National Health Service (NHS) patients in the United Kingdom (United Kingdom) to assess the association between variants in the ACE2 locus and COVID-19 risk. We analysed whole-genome sequencing (WGS) data of 1,837 cases who were tested positive for SARS-CoV-2, and 37,207 controls who were not tested, from the UK's 100,000 Genomes Project (100KGP) for the presence of ACE2 coding variants and extract expression quantitative trait loci (eQTLs). We identified a splice site variant (rs2285666) associated with increased ACE2 expression with an overrepresentation in SARS-CoV-2 positive patients relative to 100KGP controls (p = 0.015), and in hospitalised European patients relative to outpatients in intra-ethnic comparisons (p = 0.029). We also compared the prevalence of 288 eQTLs, of which 23 were enriched in SARS-CoV-2 positive patients. The eQTL rs12006793 had the largest effect size (d = 0.91), which decreases ACE2 expression and is more prevalent in controls, thus potentially reducing the risk of COVID-19. We identified three novel nonsynonymous variants predicted to alter ACE2 function, and showed that three variants (p.K26R, p. H378R, p. Y515N) alter receptor affinity for the viral Spike (S) protein. Variant p. N720D, more prevalent in the European population (p < 0.001), potentially increases viral entry by affecting the ACE2-TMPRSS2 complex. The spectrum of genetic variants in ACE2 may inform risk stratification of COVID-19 patients and could partially explain the differences in disease susceptibility and severity among different ethnic groups.

Network-based screen in iPSC-derived cells reveals therapeutic candidate for heart valve disease.

Mapping the gene-regulatory networks dysregulated in human disease would allow the design of network-correcting therapies that treat the core disease mechanism. However, small molecules are traditionally screened for their effects on one to several outputs at most, biasing discovery and limiting the likelihood of true disease-modifying drug candidates. Here, we developed a machine-learning approach to identify small molecules that broadly correct gene networks dysregulated in a human induced pluripotent stem cell (iPSC) disease model of a common form of heart disease involving the aortic valve (AV). Gene network correction by the most efficacious therapeutic candidate, XCT790, generalized to patient-derived primary AV cells and was sufficient to prevent and treat AV disease in vivo in a mouse model. This strategy, made feasible by human iPSC technology, network analysis, and machine learning, may represent an effective path for drug discovery.

Oligogenic inheritance of a human heart disease involving a genetic modifier.

Complex genetic mechanisms are thought to underlie many human diseases, yet experimental proof of this model has been elusive. Here, we show that a human cardiac anomaly can be caused by a combination of rare, inherited heterozygous mutations. Whole-exome sequencing of a nuclear family revealed that three offspring with childhood-onset cardiomyopathy had inherited three missense single-nucleotide variants in the MKL2, MYH7, and NKX2-5 genes. The MYH7 and MKL2 variants were inherited from the affected, asymptomatic father and the rare NKX2-5 variant (minor allele frequency, 0.0012) from the unaffected mother. We used CRISPR-Cas9 to generate mice encoding the orthologous variants and found that compound heterozygosity for all three variants recapitulated the human disease phenotype. Analysis of murine hearts and human induced pluripotent stem cell-derived cardiomyocytes provided histologic and molecular evidence for the NKX2-5 variant's contribution as a genetic modifier.

Context-Specific Transcription Factor Functions Regulate Epigenomic and Transcriptional Dynamics during Cardiac Reprogramming.

Ectopic expression of combinations of transcription factors (TFs) can drive direct lineage conversion, thereby reprogramming a somatic cell's identity. To determine the molecular mechanisms by which Gata4, Mef2c, and Tbx5 (GMT) induce conversion from a cardiac fibroblast toward an induced cardiomyocyte, we performed comprehensive transcriptomic, DNA-occupancy, and epigenomic interrogation throughout the reprogramming process. Integration of these datasets identified new TFs involved in cardiac reprogramming and revealed context-specific roles for GMT, including the ability of Mef2c and Tbx5 to independently promote chromatin remodeling at previously inaccessible sites. We also find evidence for cooperative facilitation and refinement of each TF's binding profile in a combinatorial setting. A reporter assay employing newly defined regulatory elements confirmed that binding of a single TF can be sufficient for gene activation, suggesting that co-binding events do not necessarily reflect synergy. These results shed light on fundamental mechanisms by which combinations of TFs direct lineage conversion.

Genetic determinants and epigenetic effects of pioneer-factor occupancy.

Transcription factors (TFs) direct developmental transitions by binding to target DNA sequences, influencing gene expression and establishing complex gene-regultory networks. To systematically determine the molecular components that enable or constrain TF activity, we investigated the genomic occupancy of FOXA2, GATA4 and OCT4 in several cell types. Despite their classification as pioneer factors, all three TFs exhibit cell-type-specific binding, even when supraphysiologically and ectopically expressed. However, FOXA2 and GATA4 can be distinguished by low enrichment at loci that are highly occupied by these factors in alternative cell types. We find that expression of additional cofactors increases enrichment at a subset of these sites. Finally, FOXA2 occupancy and changes to DNA accessibility can occur in G1-arrested cells, but subsequent loss of DNA methylation requires DNA replication.

Heart disease modelling adds a Notch to its belt.

The heart is a complex organ, consisting of multiple cell types that coordinately regulate blood flow. Reciprocal Notch pathway signalling in endocardial and myocardial cells is now shown to promote maturation of the ventricular chambers. These insights reveal mechanisms that, when disrupted, can lead to cardiomyopathies.

Transcriptional and Chromatin Dynamics of Muscle Regeneration after Severe Trauma.

Following injury, adult skeletal muscle undergoes a well-coordinated sequence of molecular and physiological events to promote repair and regeneration. However, a thorough understanding of the in vivo epigenomic and transcriptional mechanisms that control these reparative events is lacking. To address this, we monitored the in vivo dynamics of three histone modifications and coding and noncoding RNA expression throughout the regenerative process in a mouse model of traumatic muscle injury. We first illustrate how both coding and noncoding RNAs in tissues and sorted satellite cells are modified and regulated during various stages after trauma. Next, we use chromatin immunoprecipitation followed by sequencing to evaluate the chromatin state of cis-regulatory elements (promoters and enhancers) and view how these elements evolve and influence various muscle repair and regeneration transcriptional programs. These results provide a comprehensive view of the central factors that regulate muscle regeneration and underscore the multiple levels through which both transcriptional and epigenetic patterns are regulated to enact appropriate repair and regeneration.

A qPCR ScoreCard quantifies the differentiation potential of human pluripotent stem cells.

Research on human pluripotent stem cells has been hampered by the lack of a standardized, quantitative, scalable assay of pluripotency. We previously described an assay called ScoreCard that used gene expression signatures to quantify differentiation efficiency. Here we report an improved version of the assay based on qPCR that enables faster, more quantitative assessment of functional pluripotency. We provide an in-depth characterization of the revised signature panel (commercially available as the TaqMan hPSC Scorecard Assay) through embryoid body and directed differentiation experiments as well as a detailed comparison to the teratoma assay. We further show that the improved ScoreCard enables a wider range of applications, such as screening of small molecules, genetic perturbations and assessment of culture conditions. Our approach can be extended beyond stem cell applications to characterize and assess the utility of other cell types and lineages.

In vivo Monitoring of Transcriptional Dynamics After Lower-Limb Muscle Injury Enables Quantitative Classification of Healing.

Traumatic lower-limb musculoskeletal injuries are pervasive amongst athletes and the military and typically an individual returns to activity prior to fully healing, increasing a predisposition for additional injuries and chronic pain. Monitoring healing progression after a musculoskeletal injury typically involves different types of imaging but these approaches suffer from several disadvantages. Isolating and profiling transcripts from the injured site would abrogate these shortcomings and provide enumerative insights into the regenerative potential of an individual's muscle after injury. In this study, a traumatic injury was administered to a mouse model and healing progression was examined from 3 hours to 1 month using high-throughput RNA-Sequencing (RNA-Seq). Comprehensive dissection of the genome-wide datasets revealed the injured site to be a dynamic, heterogeneous environment composed of multiple cell types and thousands of genes undergoing significant expression changes in highly regulated networks. Four independent approaches were used to determine the set of genes, isoforms, and genetic pathways most characteristic of different time points post-injury and two novel approaches were developed to classify injured tissues at different time points. These results highlight the possibility to quantitatively track healing progression in situ via transcript profiling using high- throughput sequencing.