Unraveling the evolutionary origin of the complex Nuclear Receptor Element (cNRE), a cis-regulatory module required for preferential expression in the atrial chamber.

Cardiac function requires appropriate proteins in each chamber. Atria requires slow myosin to act as reservoirs, while ventricles demand fast myosin for swift pumping. Myosins are thus under chamber-biased cis-regulation, with myosin gene expression imbalances leading to congenital heart dysfunction. To identify regulatory inputs leading to cardiac chamber-biased expression, we computationally and molecularly dissected the quail Slow Myosin Heavy Chain III (SMyHC III) promoter that drives preferential expression to the atria. We show that SMyHC III gene states are orchestrated by a complex Nuclear Receptor Element (cNRE) of 32 base pairs. Using transgenesis in zebrafish and mice, we demonstrate that preferential atrial expression is achieved by a combinatorial regulatory input composed of atrial activation motifs and ventricular repression motifs. Using comparative genomics, we show that the cNRE might have emerged from an endogenous viral element through infection of an ancestral host germline, revealing an evolutionary pathway to cardiac chamber-specific expression.

Systematic molecular profiling of acute leukemia cancer stem cells allows identification of druggable targets.

Acute myeloid leukemia (AML) is one of the most prevalent and acute blood cancers with a poor prognosis and low overall survival rate, especially in the elderly. Although several new AML markers and drug targets have been recently identified, the rate of long-term cancer eradication has not improved significantly due to the presence and drug resistance of AML cancer stem cells (CSCs). Here we develop a novel computational pipeline to analyze the transcriptomic profiles of AML cancer (stem) cells and identify novel candidate AML CSC markers and drug targets. In our novel pipeline we apply a top-down meta-analysis strategy to integrate The Cancer Genome Atlas data with CSC datasets to infer cell stemness features. As a result, a set of genes termed the "AML key CSC genes" along with all the available drugs/compounds that could target them were identified. Overall, our novel computational pipeline could retrieve known cancer drugs (Carfilzomib) and predicted novel drugs such as Zonisamide, Amitriptyline, and their targets amongst the top ranked drugs and drug targets for targeting AML. Additionally, the pipeline applied in this study could be used for the identification of CSC-specific markers, drivers and their respective targeting drugs in other cancer types.

Human pluripotent stem cell-derived macrophages host Mycobacterium abscessus infection.

Human macrophages are a natural host of many mycobacterium species, including Mycobacterium abscessus (M. abscessus), an emerging pathogen affecting immunocompromised and cystic fibrosis patients with few available treatments. The search for an effective treatment is hindered by the lack of a tractable in vitro intracellular infection model. Here, we established a reliable model for M. abscessus infection using human pluripotent stem cell-derived macrophages (hPSC-macrophages). hPSC differentiation permitted reproducible generation of functional macrophages that were highly susceptible to M. abscessus infection. Electron microscopy demonstrated that M. abscessus was present in the hPSC-macrophage vacuoles. RNA sequencing analysis revealed a time-dependent host cell response, with differing gene and protein expression patterns post-infection. Engineered tdTOMATO-expressing hPSC-macrophages with GFP-expressing mycobacteria enabled rapid image-based high-throughput analysis of intracellular infection and quantitative assessment of antibiotic efficacy. Our study describes the first to our knowledge hPSC-based model for M. abscessus infection, representing a novel and accessible system for studying pathogen-host interaction and drug discovery.

CD90 Marks a Mesenchymal Program in Human Thymic Epithelial Cells In Vitro and In Vivo.

Thymic epithelium is critical for the structural integrity of the thymus and for T cell development. Within the fully formed thymus, large numbers of hematopoietic cells shape the thymic epithelium into a scaffold-like structure which bears little similarity to classical epithelial layers, such as those observed in the skin, intestine or pancreas. Here, we show that human thymic epithelial cells (TECs) possess an epithelial identity that also incorporates the expression of mesenchymal cell associated genes, whose expression levels vary between medullary and cortical TECs (m/cTECs). Using pluripotent stem cell (PSC) differentiation systems, we identified a unique population of cells that co-expressed the master TEC transcription factor FOXN1, as well as the epithelial associated marker EPCAM and the mesenchymal associated gene CD90. Using the same serum free culture conditions, we also observed co-expression of EPCAM and CD90 on cultured TECs derived from neonatal human thymus in vitro. Single cell RNA-sequencing revealed these cultured TECs possessed an immature mTEC phenotype and expressed epithelial and mesenchymal associated genes, such as EPCAM, CLDN4, CD90 and COL1A1. Importantly, flow cytometry and single cell RNA-sequencing analysis further confirmed the presence of an EPCAM+CD90+ population in the CD45- fraction of neonatal human thymic stromal cells in vivo. Using the human thymus cell atlas, we found that cTECs displayed more pronounced mesenchymal characteristics than mTECs during embryonic development. Collectively, these results suggest human TECs possess a hybrid gene expression program comprising both epithelial and mesenchymal elements, and provide a basis for the further exploration of thymus development from primary tissues and from the in vitro differentiation of PSCs.

The RNA-binding protein SRSF3 has an essential role in megakaryocyte maturation and platelet production.

RNA processing is increasingly recognized as a critical control point in the regulation of different hematopoietic lineages including megakaryocytes responsible for the production of platelets. Platelets are anucleate cytoplasts that contain a rich repertoire of RNAs encoding proteins with essential platelet functions derived from the parent megakaryocyte. It is largely unknown how RNA binding proteins contribute to the development and functions of megakaryocytes and platelets. We show that serine-arginine-rich splicing factor 3 (SRSF3) is essential for megakaryocyte maturation and generation of functional platelets. Megakaryocyte-specific deletion of Srsf3 in mice led to macrothrombocytopenia characterized by megakaryocyte maturation arrest, dramatically reduced platelet counts, and abnormally large functionally compromised platelets. SRSF3 deficient megakaryocytes failed to reprogram their transcriptome during maturation and to load platelets with RNAs required for normal platelet function. SRSF3 depletion led to nuclear accumulation of megakaryocyte mRNAs, demonstrating that SRSF3 deploys similar RNA regulatory mechanisms in megakaryocytes as in other cell types. Our study further suggests that SRSF3 plays a role in sorting cytoplasmic megakaryocyte RNAs into platelets and demonstrates how SRSF3-mediated RNA processing forms a central part of megakaryocyte gene regulation. Understanding SRSF3 functions in megakaryocytes and platelets provides key insights into normal thrombopoiesis and platelet pathologies as SRSF3 RNA targets in megakaryocytes are associated with platelet diseases.

Navigating the non-coding genome in heart development and Congenital Heart Disease.

Congenital Heart Disease (CHD) is characterised by a wide range of cardiac defects, from mild to life-threatening, which occur in babies worldwide. To date, there is no cure to CHD, however, progress in surgery has reduced its mortality allowing children affected by CHD to reach adulthood. In an effort to understand its genetic basis, several studies involving whole-genome sequencing (WGS) of patients with CHD have been undertaken and generated a great wealth of information. The majority of putative causative mutations identified in WGS studies fall into the non-coding part of the genome. Unfortunately, due to the lack of understanding of the function of these non-coding mutations, it is challenging to establish a causal link between the non-coding mutation and the disease. Thus, here we review the state-of-the-art approaches to interpret non-coding mutations in the context of CHD and address the following questions: What are the non-coding sequences important for cardiac function? Which technologies are used to identify them? Which resources are available to analyse them? What mutations are expected in these non-coding sequences? Learning from developmental process, what is their expected role in CHD?

Control of cardiac jelly dynamics by NOTCH1 and NRG1 defines the building plan for trabeculation.

In vertebrate hearts, the ventricular trabecular myocardium develops as a sponge-like network of cardiomyocytes that is critical for contraction and conduction, ventricular septation, papillary muscle formation and wall thickening through the process of compaction 1 . Defective trabeculation leads to embryonic lethality2-4 or non-compaction cardiomyopathy (NCC) 5 . There are divergent views on when and how trabeculation is initiated in different species. In zebrafish, trabecular cardiomyocytes extrude from compact myocardium 6 , whereas in chicks, chamber wall thickening occurs before overt trabeculation 7 . In mice, the onset of trabeculation has not been described, but is proposed to begin at embryonic day 9.0, when cardiomyocytes form radially oriented ribs 2 . Endocardium-myocardium communication is essential for trabeculation, and numerous signalling pathways have been identified, including Notch2,8 and Neuregulin (NRG) 4 . Late disruption of the Notch pathway causes NCC 5 . Whereas it has been shown that mutations in the extracellular matrix (ECM) genes Has2 and Vcan prevent the formation of trabeculae in mice9,10 and the matrix metalloprotease ADAMTS1 promotes trabecular termination 3 , the pathways involved in ECM dynamics and the molecular regulation of trabeculation during its early phases remain unexplored. Here we present a model of trabeculation in mice that integrates dynamic endocardial and myocardial cell behaviours and ECM remodelling, and reveal new epistatic relationships between the involved signalling pathways. NOTCH1 signalling promotes ECM degradation during the formation of endocardial projections that are critical for individualization of trabecular units, whereas NRG1 promotes myocardial ECM synthesis, which is necessary for trabecular rearrangement and growth. These systems interconnect through NRG1 control of Vegfa, but act antagonistically to establish trabecular architecture. These insights enabled the prediction of persistent ECM and cardiomyocyte growth in a mouse NCC model, providing new insights into the pathophysiology of congenital heart disease.

Germline Stem Cell Activity Is Sustained by SALL4-Dependent Silencing of Distinct Tumor Suppressor Genes.

Sustained spermatogenesis in adult males and fertility recovery following germ cell depletion are dependent on undifferentiated spermatogonia. We previously demonstrated a key role for the transcription factor SALL4 in spermatogonial differentiation. However, whether SALL4 has broader roles within spermatogonia remains unclear despite its ability to co-regulate genes with PLZF, a transcription factor required for undifferentiated cell maintenance. Through development of inducible knockout models, we show that short-term integrity of differentiating but not undifferentiated populations requires SALL4. However, SALL4 loss was associated with long-term functional decline of undifferentiated spermatogonia and disrupted stem cell-driven regeneration. Mechanistically, SALL4 associated with the NuRD co-repressor and repressed expression of the tumor suppressor genes Foxl1 and Dusp4. Aberrant Foxl1 activation inhibited undifferentiated cell growth and survival, while DUSP4 suppressed self-renewal pathways. We therefore uncover an essential role for SALL4 in maintenance of undifferentiated spermatogonial activity and identify regulatory pathways critical for germline stem cell function.

Multicellular Transcriptional Analysis of Mammalian Heart Regeneration.

BACKGROUND: The inability of the adult mammalian heart to regenerate following injury represents a major barrier in cardiovascular medicine. In contrast, the neonatal mammalian heart retains a transient capacity for regeneration, which is lost shortly after birth. Defining the molecular mechanisms that govern regenerative capacity in the neonatal period remains a central goal in cardiac biology. Here, we assemble a transcriptomic framework of multiple cardiac cell populations during postnatal development and following injury, which enables comparative analyses of the regenerative (neonatal) versus nonregenerative (adult) state for the first time. METHODS: Cardiomyocytes, fibroblasts, leukocytes, and endothelial cells from infarcted and noninfarcted neonatal (P1) and adult (P56) mouse hearts were isolated by enzymatic dissociation and fluorescence-activated cell sorting at day 3 following surgery. RNA sequencing was performed on these cell populations to generate the transcriptome of the major cardiac cell populations during cardiac development, repair, and regeneration. To complement our transcriptomic data, we also surveyed the epigenetic landscape of cardiomyocytes during postnatal maturation by performing deep sequencing of accessible chromatin regions by using the Assay for Transposase-Accessible Chromatin from purified mouse cardiomyocyte nuclei (P1, P14, and P56). RESULTS: Profiling of cardiomyocyte and nonmyocyte transcriptional programs uncovered several injury-responsive genes across regenerative and nonregenerative time points. However, the majority of transcriptional changes in all cardiac cell types resulted from developmental maturation from neonatal stages to adulthood rather than activation of a distinct regeneration-specific gene program. Furthermore, adult leukocytes and fibroblasts were characterized by the expression of a proliferative gene expression network following infarction, which mirrored the neonatal state. In contrast, cardiomyocytes failed to reactivate the neonatal proliferative network following infarction, which was associated with loss of chromatin accessibility around cell cycle genes during postnatal maturation. CONCLUSIONS: This work provides a comprehensive framework and transcriptional resource of multiple cardiac cell populations during cardiac development, repair, and regeneration. Our findings define a regulatory program underpinning the neonatal regenerative state and identify alterations in the chromatin landscape that could limit reinduction of the regenerative program in adult cardiomyocytes.

Combinatorial Ranking of Gene Sets to Predict Disease Relapse: The Retinoic Acid Pathway in Early Prostate Cancer.

BACKGROUND: Quantitative high-throughput data deposited in consortia such as International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA) present opportunities and challenges for computational analyses. METHODS: We present a computational strategy to systematically rank and investigate a large number (210-220) of clinically testable gene sets, using combinatorial gene subset generation and disease-free survival (DFS) analyses. This approach integrates protein-protein interaction networks, gene expression, DNA methylation, and copy number data, in association with DFS profiles from patient clinical records. RESULTS: As a case study, we applied this pipeline to systematically analyze the role of ALDH1A2 in prostate cancer (PCa). We have previously found this gene to have multiple roles in disease and homeostasis, and here we investigate the role of the associated ALDH1A2 gene/protein networks in PCa, using our methodology in combination with PCa patient clinical profiles from ICGC and TCGA databases. Relationships between gene signatures and relapse were analyzed using Kaplan-Meier (KM) log-rank analysis and multivariable Cox regression. Relative expression versus pooled mean from diploid population was used for z-statistics calculation. Gene/protein interaction network analyses generated 11 core genes associated with ALDH1A2; combinatorial ranking of the power set of these core genes identified two gene sets (out of 211 - 1 = 2,047 combinations) with significant correlation with disease relapse (KM log rank p < 0.05). For the more significant of these two sets, referred to as the optimal gene set (OGS), patients have median survival 62.7 months with OGS alterations compared to >150 months without OGS alterations (p = 0.0248, hazard ratio = 2.213, 95% confidence interval = 1.1-4.098). Two genes comprising OGS (CYP26A1 and RDH10) are strongly associated with ALDH1A2 in the retinoic acid (RA) pathways, suggesting a major role of RA signaling in early PCa progression. Our pipeline complements human expertise in the search for prognostic biomarkers in large-scale datasets.

Tyrosine phosphorylation of LRP6 by Src and Fer inhibits Wnt/ß-catenin signalling.

Low-density lipoprotein receptor-related proteins 5 and 6 (LRP5/6) function as transmembrane receptors to transduce Wnt signals. A key mechanism for signalling is Wnt-induced serine/threonine phosphorylation at conserved PPPSPxS motifs in the LRP6 cytoplasmic domain, which promotes pathway activation. Conserved tyrosine residues are positioned close to all PPPSPxS motifs, which suggests they have a functional significance. Using a cell culture-based cDNA expression screen, we identified the non-receptor tyrosine kinases Src and Fer as novel LRP6 modifiers. Both Src and Fer associate with LRP6 and phosphorylate LRP6 directly. In contrast to the known PPPSPxS Ser/Thr kinases, tyrosine phosphorylation by Src and Fer negatively regulates LRP6-Wnt signalling. Epistatically, they function upstream of ß-catenin to inhibit signalling and in agreement with a negative role in regulating LRP6, MEF cells lacking these kinases show enhanced Wnt signalling. Wnt3a treatment of cells enhances tyrosine phosphorylation of endogenous LRP6 and, mechanistically, Src reduces cell surface LRP6 levels and disrupts LRP6 signalosome formation. Interestingly, CK1γ inhibits Fer-induced LRP6 phosphorylation, suggesting a mechanism whereby CK1γ acts to de-represses inhibitory LRP6 tyrosine phosphorylation. We propose that LRP6 tyrosine phosphorylation by Src and Fer serves a negative regulatory function to prevent over-activation of Wnt signalling at the level of the Wnt receptor, LRP6.

Haematopoietic stem cell induction by somite-derived endothelial cells controlled by meox1.

Haematopoietic stem cells (HSCs) are self-renewing stem cells capable of replenishing all blood lineages. In all vertebrate embryos that have been studied, definitive HSCs are generated initially within the dorsal aorta (DA) of the embryonic vasculature by a series of poorly understood inductive events. Previous studies have identified that signalling relayed from adjacent somites coordinates HSC induction, but the nature of this signal has remained elusive. Here we reveal that somite specification of HSCs occurs via the deployment of a specific endothelial precursor population, which arises within a sub-compartment of the zebrafish somite that we have defined as the endotome. Endothelial cells of the endotome are specified within the nascent somite by the activity of the homeobox gene meox1. Specified endotomal cells consequently migrate and colonize the DA, where they induce HSC formation through the deployment of chemokine signalling activated in these cells during endotome formation. Loss of meox1 activity expands the endotome at the expense of a second somitic cell type, the muscle precursors of the dermomyotomal equivalent in zebrafish, the external cell layer. The resulting increase in endotome-derived cells that migrate to colonize the DA generates a dramatic increase in chemokine-dependent HSC induction. This study reveals the molecular basis for a novel somite lineage restriction mechanism and defines a new paradigm in induction of definitive HSCs.