Deep-transcriptome and ribonome sequencing redefines the molecular networks of pluripotency and the extracellular space in human embryonic stem cells

Transcriptional Profiling of Stem Cells: Moving from Descriptive to Predictive Paradigms

Transcriptional profiling is a powerful tool commonly used to benchmark stem cells and their differentiated progeny. As the wealth of stem cell data builds in public repositories, we highlight common data traps, and review approaches to combine and mine this data for new cell classification and cell prediction tools. We touch on future trends for stem cell profiling, such as single-cell profiling, long-read sequencing, and improved methods for measuring molecular modifications on chromatin and RNA that bring new challenges and opportunities for stem cell analysis.

Cardiomyogenic differentiation is fine-tuned by differential mRNA association with polysomes

BACKGROUND

Cardiac cell fate specification occurs through progressive steps, and its gene expression regulation features are still being defined. There has been an increasing interest in understanding the coordination between transcription and post-transcriptional regulation during the differentiation processes. Here, we took advantage of the polysome profiling technique to isolate and high-throughput sequence ribosome-free and polysome-bound RNAs during cardiomyogenesis.

RESULTS

We showed that polysome-bound RNAs exhibit the cardiomyogenic commitment gene expression and that mesoderm-to-cardiac progenitor stages are strongly regulated. Additionally, we compared ribosome-free and polysome-bound RNAs and found that the post-transcriptional regulation vastly contributes to cardiac phenotype determination, including RNA recruitment to and dissociation from ribosomes. Moreover, we found that protein synthesis is decreased in cardiomyocytes compared to human embryonic stem-cells (hESCs), possibly due to the down-regulation of translation-related genes.

CONCLUSIONS

Our data provided a powerful tool to investigate genes potentially controlled by post-transcriptional mechanisms during the cardiac differentiation of hESC. This work could prospect fundamental tools to develop new therapy and research approaches.

Polysome profiling followed by RNA-seq of cardiac differentiation stages in hESCs

The regulation of gene expression acts at numerous complementary levels to control and refine protein abundance. The analysis of mRNAs associated with polysomes, called polysome profiling, has been used to investigate the post-transcriptional mechanisms that are involved in different biological processes. Pluripotent stem cells are able to differentiate into a variety of cell lineages, and the cell commitment progression is carefully orchestrated. Genome-wide expression profiling has provided the possibility to investigate transcriptional changes during cardiomyogenesis; however, a more accurate study regarding post-transcriptional regulation is required. In the present work, we isolated and high-throughput sequenced ribosome-free and polysome-bound RNAs from NKX2-5^eGFP/w HES3 undifferentiated pluripotent stem cells at the subsequent differentiation stages of cardiomyogenesis: embryoid body aggregation, mesoderm, cardiac progenitor and cardiomyocyte. The expression of developmental markers was followed by flow cytometry, and quality analyses were performed as technical controls to ensure high quality data. Our dataset provides valuable information about hESC cardiac differentiation and can be used to investigate genes potentially controlled by post-transcriptional mechanisms.

Downregulation of the protein synthesis machinery is a major regulatory event during early adipogenic differentiation of human adipose-derived stromal cells

Commitment of adult stem cells involves the activation of specific gene networks regulated from transcription to protein synthesis. Here, we used ribosome profiling to identify mRNAs regulated at the translational level, through both differential association to polysomes and modulation of their translational rates. We observed that translational regulation during the differentiation of human adipose-derived stromal cells (hASCs, also known as adipose-derived mesenchymal stem cells), a subset of which are stem cells, to adipocytes was a major regulatory event. hASCs showed a significant reduction of whole protein synthesis after adipogenic induction and a downregulation of the expression and translational efficiency of ribosomal proteins. Additionally, focal adhesion and cytoskeletal proteins were downregulated at the translational level. This negative regulation of the essential biological functions of hASCs resulted in a reduction in cell size and the potential of hASCs to migrate. We analyzed whether the inactivation of key translation initiation factors was involved in this observed major repression of translation. We showed that there was an increase in the hypo phosphorylated forms of 4E-BP1, a negative regulator of translation, during early adipogenesis. Our results showed that extensive translational regulation occurred during the early stage of the adipogenic differentiation of hASCs.

Multilayered Control of Alternative Splicing Regulatory Networks by Transcription Factors

Networks of coordinated alternative splicing (AS) events play critical roles in development and disease. However, a comprehensive knowledge of the factors that regulate these networks is lacking. We describe a high-throughput system for systematically linking trans-acting factors to endogenous RNA regulatory events. Using this system, we identify hundreds of factors associated with diverse regulatory layers that positively or negatively control AS events linked to cell fate. Remarkably, more than one-third of the regulators are transcription factors. Further analyses of the zinc finger protein Zfp871 and BTB/POZ domain transcription factor Nacc1, which regulate neural and stem cell AS programs, respectively, reveal roles in controlling the expression of specific splicing regulators. Surprisingly, these proteins also appear to regulate target AS programs via binding RNA. Our results thus uncover a large "missing cache" of splicing regulators among annotated transcription factors, some of which dually regulate AS through direct and indirect mechanisms.

New Monoclonal Antibodies to Defined Cell Surface Proteins on Human Pluripotent Stem Cells

The study and application of human pluripotent stem cells (hPSCs) will be enhanced by the availability of well‐characterized monoclonal antibodies (mAbs) detecting cell‐surface epitopes. Here, we report generation of seven new mAbs that detect cell surface proteins present on live and fixed human ES cells (hESCs) and human iPS cells (hiPSCs), confirming our previous prediction that these proteins were present on the cell surface of hPSCs. The mAbs all show a high correlation with POU5F1 (OCT4) expression and other hPSC surface markers (TRA‐160 and SSEA‐4) in hPSC cultures and detect rare OCT4 positive cells in differentiated cell cultures. These mAbs are immunoreactive to cell surface protein epitopes on both primed and naive state hPSCs, providing useful research tools to investigate the cellular mechanisms underlying human pluripotency and states of cellular reprogramming. In addition, we report that subsets of the seven new mAbs are also immunoreactive to human bone marrow‐derived mesenchymal stem cells (MSCs), normal human breast subsets and both normal and tumorigenic colorectal cell populations. The mAbs reported here should accelerate the investigation of the nature of pluripotency, and enable development of robust cell separation and tracing technologies to enrich or deplete for hPSCs and other human stem and somatic cell types. Stem Cells2017;35:626–640

Estimating mRNA lengths from Plasmodium falciparum genes by Virtual Northern RNA-seq analysis

Accurate gene models are essential for understanding parasite biology. However, transcript structure information is lacking for most parasite genes. Here, we describe "Virtual Northern" analysis of the malaria parasite Plasmodium falciparum to address this issue. RNA-seq libraries were made from size-fractionated RNA. Transcript sizes for 3052 genes were inferred from the read counts in each library. The data show that for almost half of the transcripts, the combined untranslated regions are more than twice the length of the open reading frame. Furthermore, we identified novel polycistronic, or gene overlapping, transcripts that suggest revisions to current gene models are needed.

Musashi2 sustains the mixed-lineage leukemia-driven stem cell regulatory program

Leukemia stem cells (LSCs) are found in most aggressive myeloid diseases and contribute to therapeutic resistance. Leukemia cells exhibit a dysregulated developmental program as the result of genetic and epigenetic alterations. Overexpression of the RNA-binding protein Musashi2 (MSI2) has been previously shown to predict poor survival in leukemia. Here, we demonstrated that conditional deletion of Msi2 in the hematopoietic compartment results in delayed leukemogenesis, reduced disease burden, and a loss of LSC function in a murine leukemia model. Gene expression profiling of these Msi2-deficient animals revealed a loss of the hematopoietic/leukemic stem cell self-renewal program and an increase in the differentiation program. In acute myeloid leukemia patients, the presence of a gene signature that was similar to that observed in Msi2-deficent murine LSCs correlated with improved survival. We determined that MSI2 directly maintains the mixed-lineage leukemia (MLL) self-renewal program by interacting with and retaining efficient translation of Hoxa9, Myc, and Ikzf2 mRNAs. Moreover, depletion of MLL target Ikzf2 in LSCs reduced colony formation, decreased proliferation, and increased apoptosis. Our data provide evidence that MSI2 controls efficient translation of the oncogenic LSC self-renewal program and suggest MSI2 as a potential therapeutic target for myeloid leukemia.

Alternative polyadenylation in the nervous system: to what lengths will 3' UTR extensions take us?

Alternative cleavage and polyadenylation (APA) can diversify coding and non-coding regions, but has particular impact on increasing 3' UTR diversity. Through the gain or loss of regulatory elements such as RNA binding protein and microRNA sites, APA can influence transcript stability, localization, and translational efficiency. Strikingly, the central nervous systems of invertebrate and vertebrate species express a broad range of transcript isoforms bearing extended 3' UTRs. The molecular mechanism that permits proximal 3' end bypass in neurons is mysterious, and only beginning to be elucidated. This landscape of neural 3' UTR extensions, many reaching unprecedented lengths, may help service the unique post-transcriptional regulatory needs of neurons. A combination of approaches, including transcriptome-wide profiling, genetic screening to identify APA factors, biochemical dissection of alternative 3' end formation, and manipulation of individual neural APA targets, will be necessary to gain fuller perspectives on the mechanism and biology of neural-specific 3' UTR lengthening.

Gene expression variability as a unifying element of the pluripotency network

Heterogeneity is a hallmark of stem cell populations, in part due to the molecular differences between cells undergoing self-renewal and those poised to differentiate. We examined phenotypic and molecular heterogeneity in pluripotent stem cell populations, using public gene expression data sets. A high degree of concordance was observed between global gene expression variability and the reported heterogeneity of different human pluripotent lines. Network analysis demonstrated that low-variability genes were the most highly connected, suggesting that these are the most stable elements of the gene regulatory network and are under the highest regulatory constraints. Known drivers of pluripotency were among these, with lowest expression variability of POU5F1 in cells with the highest capacity for self-renewal. Variability of gene expression provides a reliable measure of phenotypic and molecular heterogeneity and predicts those genes with the highest degree of regulatory constraint within the pluripotency network.

•

Gene expression variability is highly concordant with population heterogeneity

•

Genes within the pluripotency network have distinct variability profiles

•

Expression variability is a network property important for pluripotency

Role of cell–cell adhesion complexes in embryonic stem cell biology

Pluripotent embryonic stem cells (ESCs) can self-renew or differentiate into any cell type within an organism. Here, we focus on the roles of cadherins and catenins – their cytoplasmic scaffold proteins – in the fate, maintenance and differentiation of mammalian ESCs. E-cadherin is a master stem cell regulator that is required for both mouse ESC (mESC) maintenance and differentiation. E-cadherin interacts with key components of the naive stemness pathway and ablating it prevents stem cells from forming well-differentiated teratomas or contributing to chimeric animals. In addition, depleting E-cadherin converts naive mouse ESCs into primed epiblast-like stem cells (EpiSCs). In line with this, a mesenchymal-to-epithelial transition (MET) occurs during reprogramming of somatic cells towards induced pluripotent stem cells (iPSCs), leading to downregulation of N-cadherin and acquisition of high E-cadherin levels. β-catenin exerts a dual function; it acts in cadherin-based adhesion and in WNT signaling and, although WNT signaling is important for stemness, the adhesive function of β-catenin might be crucial for maintaining the naive state of stem cells. In addition, evidence is rising that other junctional proteins are also important in ESC biology. Thus, precisely regulated levels and activities of several junctional proteins, in particular E-cadherin, safeguard naive pluripotency and are a prerequisite for complete somatic cell reprogramming.

Deep transcriptome profiling of mammalian stem cells supports a regulatory role for retrotransposons in pluripotency maintenance

The importance of microRNAs and long noncoding RNAs in the regulation of pluripotency has been documented; however, the noncoding components of stem cell gene networks remain largely unknown. Here we investigate the role of noncoding RNAs in the pluripotent state, with particular emphasis on nuclear and retrotransposon-derived transcripts. We have performed deep profiling of the nuclear and cytoplasmic transcriptomes of human and mouse stem cells, identifying a class of previously undetected stem cell-specific transcripts. We show that long terminal repeat (LTR)-derived transcripts contribute extensively to the complexity of the stem cell nuclear transcriptome. Some LTR-derived transcripts are associated with enhancer regions and are likely to be involved in the maintenance of pluripotency.

PUNCH-P for global translatome profiling: Methodology, insights and comparison to other techniques

Regulation of mRNA translation is a major modulator of gene expression, allowing cells to fine tune protein levels during growth and differentiation and in response to physiological signals and environmental changes. Mass-spectrometry and RNA-sequencing methods now enable global profiling of the translatome, but these still involve significant analytical and economical limitations. We developed a novel system-wide proteomic approach for direct monitoring of translation, termed PUromycin-associated Nascent CHain Proteomics (PUNCH-P), which is based on the recovery of ribosome-nascent chain complexes from cells or tissues followed by incorporation of biotinylated puromycin into newly-synthesized proteins. Biotinylated proteins are then purified by streptavidin and analyzed by mass-spectrometry. Here we present an overview of PUNCH-P, describe other methodologies for global translatome profiling (pSILAC, BONCAT, TRAP/Ribo-tag, Ribo-seq) and provide conceptual comparisons between these methods. We also show how PUNCH-P data can be combined with mRNA measurements to determine relative translation efficiency for specific mRNAs.

Role of alternative polyadenylation during adipogenic differentiation: an in silico approach

Post-transcriptional regulation of stem cell differentiation is far from being completely understood. Changes in protein levels are not fully correlated with corresponding changes in mRNAs; the observed differences might be partially explained by post-transcriptional regulation mechanisms, such as alternative polyadenylation. This would involve changes in protein binding, transcript usage, miRNAs and other non-coding RNAs. In the present work we analyzed the distribution of alternative transcripts during adipogenic differentiation and the potential role of miRNAs in post-transcriptional regulation. Our in silico analysis suggests a modest, consistent, bias in 3'UTR lengths during differentiation enabling a fine-tuned transcript regulation via small non-coding RNAs. Including these effects in the analyses partially accounts for the observed discrepancies in relative abundance of protein and mRNA.

Transcriptome dynamics during human erythroid differentiation and development

To explore the mechanisms controlling erythroid differentiation and development, we analyzed the genome-wide transcription dynamics occurring during the differentiation of human embryonic stem cells (HESCs) into the erythroid lineage and development of embryonic to adult erythropoiesis using high throughput sequencing technology. HESCs and erythroid cells at three developmental stages: ESER (embryonic), FLER (fetal), and PBER (adult) were analyzed. Our findings revealed that the number of expressed genes decreased during differentiation, whereas the total expression intensity increased. At each of the three transitions (HESCs-ESERs, ESERs-FLERs, and FLERs-PBERs), many differentially expressed genes were observed, which were involved in maintaining pluripotency, early erythroid specification, rapid cell growth, and cell-cell adhesion and interaction. We also discovered dynamic networks and their central nodes in each transition. Our study provides a fundamental basis for further investigation of erythroid differentiation and development, and has implications in using ESERs for transfusion product in clinical settings.

Polysome profiling shows extensive posttranscriptional regulation during human adipocyte stem cell differentiation into adipocytes

Adipocyte stem cells (hASCs) can proliferate and self-renew and, due to their multipotent nature, they can differentiate into several tissue-specific lineages, making them ideal candidates for use in cell therapy. Most attempts to determine the mRNA profile of self-renewing or differentiating stem cells have made use of total RNA for gene expression analysis. Several lines of evidence suggest that self-renewal and differentiation are also dependent on the control of protein synthesis by posttranscriptional mechanisms. We used adipogenic differentiation as a model, to investigate the extent to which posttranscriptional regulation controlled gene expression in hASCs. We focused on the initial steps of differentiation and isolated both the total mRNA fraction and the subpopulation of mRNAs associated with translating ribosomes. We observed that adipogenesis is committed in the first days of induction and three days appears as the minimum time of induction necessary for efficient differentiation. RNA-seq analysis showed that a significant percentage of regulated mRNAs were posttranscriptionally controlled. Part of this regulation involves massive changes in transcript untranslated regions (UTR) length, with differential extension/reduction of the 3'UTR after induction. A slight correlation can be observed between the expression levels of differentially expressed genes and the 3'UTR length. When we considered association to polysomes, this correlation values increased. Changes in the half lives were related to the extension of the 3'UTR, with longer UTRs mainly stabilizing the transcripts. Thus, changes in the length of these extensions may be associated with changes in the ability to associate with polysomes or in half-life.

Translation regulation gets its 'omics' moment

The fate of cellular RNA is largely determined by complex networks of protein-RNA interactions through ribonucleoprotein (RNP) complexes. Despite their relatively short half-life, transcripts associate with many different proteins that process, modify, translate, and degrade the RNA. Following biogenesis some mRNPs are immediately directed to translation and produce proteins, but many are diverted and regulated by processes including miRNA-mediated mechanisms, transport and localization, as well as turnover. Because of this complex interplay estimates of steady-state expression by methods such as RNAseq alone cannot capture critical aspects of cellular fate, environmental response, tumorigenesis, or gene expression regulation. More selective and integrative tools are needed to measure protein-RNA complexes and the regulatory processes involved. One focus area is measurements of the transcriptome associated with ribosomes and translation. These so-called polysome or ribosome profiling techniques can evaluate translation efficiency as well as the interplay between translation initiation, elongation, and termination-subject areas not well understood at a systems biology level. Ribosome profiling is a highly promising technique that provides mRNA positional information of ribosome occupancy, potentially bridging the gap between gene expression (i.e., RNAseq and microarray analysis) and protein quantification (i.e., mass spectrometry). In combination with methods such as RNA immunoprecipitation, miRNA profiling, or proteomics, we obtain a fresh view of global post-transcriptional and translational gene regulation. In addition, these techniques also provide new insight into new regulatory elements, such as alternative open reading frames, and translation regulation under different conditions.

Non-canonical TAF complexes regulate active promoters in human embryonic stem cells

The general transcription factor TFIID comprises the TATA-box-binding protein (TBP) and approximately 14 TBP-associated factors (TAFs). Here we find, unexpectedly, that undifferentiated human embryonic stem cells (hESCs) contain only six TAFs (TAFs 2, 3, 5, 6, 7 and 11), whereas following differentiation all TAFs are expressed. Directed and global chromatin immunoprecipitation analyses reveal an unprecedented promoter occupancy pattern: most active genes are bound by only TAFs 3 and 5 along with TBP, whereas the remaining active genes are bound by TBP and all six hESC TAFs. Consistent with these results, hESCs contain a previously undescribed complex comprising TAFs 2, 6, 7, 11 and TBP. Altering the composition of hESC TAFs, either by depleting TAFs that are present or ectopically expressing TAFs that are absent, results in misregulated expression of pluripotency genes and induction of differentiation. Thus, the selective expression and use of TAFs underlies the ability of hESCs to self-renew.DOI:http://dx.doi.org/10.7554/eLife.00068.001.