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.
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.
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.
Recent approaches regarding hESCs differentiation made it possible to mimic features of developmental biology and address the key events that regulate early lineage commitment . Cardiomyocytes have been derived from hESCs as an in vitro model to study cardiomyogenesis and as an attempt to produce clinically relevant cell populations [2,3,4,5]. Modeling congenital abnormalities of the heart or testing the cardiac toxicity of new drugs strengthen a particular interest in deriving cardiovascular lineages from pluripotent stem cells . Additionally, in vitro cardiomyocyte generation and promotion of an endogenous regenerative capacity offers new therapeutic strategies to replace heart tissue damaged by age or disease [6, 7]. Activation of a specific genetic program is essential to drive cells into the cell type of interest and contribute to these challenging goals. Cardiac tissue formation is controlled by sequential gene regulatory steps that define specialized cell fates , although understanding of molecular signatures of intermediate differentiation states of the cardiomyogenic lineage are still lacking.
Gene expression regulation acts at both transcriptional and post-transcriptional levels. There has been an increasing interest in understanding how these regulatory mechanisms interact and coordinate along the differentiation processes . Large-scale studies help to uncover the complex and dynamic multi-layered regulation involving chromatin modifications, transcriptional, and post-transcriptional networks that control the reconfiguration of ESC gene expression program when cells are stimulated to differentiate [10, 11]. Genome-wide expression profiling, using microarrays or RNA-seq analysis, has provided the possibility to investigate transcriptional changes during cardiomyogenic differentiation [12,13,14,15,16,17]. However, most of these studies were centered at total transcriptome analysis. Since the cellular abundance of proteins is also controlled at the translation level , the most-reliable gene expression patterns cannot be based solely on total RNA analysis.
The role of translational control as a key regulatory node in gene expression during development and disease is still being established. Actively translated mRNAs are predicted to be associated with many ribosomes, which form large complexes called polysomes . These complexes define the translatome, the subpopulation of mRNAs that are engaged with the translation machinery. Methods to infer the translatome have helped to expand our knowledge of protein synthesis control and revealed post-transcriptional mechanisms involved in cell fate commitment [20,21,22]. A classical technique called polysome profiling has been used to assess the mRNAs bound to ribosomes through ultracentrifugation in a sucrose gradient, and further identify them by high-throughput methods [23,24,25]. Previously, we used polysome profiling to study the fate of human adipose stem cells (hASCs) and their commitment to adipogenesis [25, 26]. We showed that 60% of the genes, which were differentially expressed after 72 h of differentiation induction, were controlled by post-transcriptional regulation.
Here, we used polysome profiling to assess the ribosome-free and polysome-bound mRNAs during the developmental steps of cardiomyogenic commitment and investigate genes potentially controlled by post-transcriptional mechanisms. Our results showed that 60–80% of differentially expressed genes (DEG) showed some degree of post-transcriptional regulation, suggesting its crucial role in shaping of the developmental pathways. Additionally, committed cells showed decreased protein synthesis rates compared to hESC.
Polysome profiling during hESCs cardiomyogenic differentiation
The hES-NKX2–5eGFP/w reporter human embryonic stem cell (hESC) line  was used to derive cardiomyocytes using a developmentally staged protocol [2, 28] that includes the induction of a cardiac mesoderm population on days 3 and 4 and a NKX2–5+/cTNT+ population by day 15 (Fig. 1a and b). Cardiomyogenesis progression was followed by flow cytometry using CD56 as a mesoderm marker  on day 4 (cut off < 40%) and NKX2–5/eGFP expression on day 9 as a cardiac progenitor marker (cut off < 50%) (Fig. 1b). Beating clusters were observed after 10 days of differentiation (Fig. 1c and Online Additional file 1: Video S1), yielding a population of 50–60% cTNT+ cardiomyocytes on day 15 (D15) (Fig. 1b). At day 20, cTNI immunostaining showed the striations characteristic of sarcomere structures (Fig. 1c).
To investigate the differential association of mRNAs with polysomes and, therefore, post-transcriptional changes in gene expression during cardiac differentiation, we first performed polysome profiling on days D0, D1, D4, D9 and D15, which represent pluripotency, embryoid body (EB) aggregation, cardiac mesoderm, cardiac progenitor and cardiomyocyte stages, respectively (Fig. 1a). After 10 min of cycloheximide treatment, active ribosomes got arrested with associated RNAs and we analyzed cellular extracts by ultracentrifugation in a sucrose gradient (Fig. 1d). Differential density throughout the gradient allowed the isolation of ribosome-free (fractions 1–3) and polysome-bound (fractions 9–22) RNAs. Based on the polysome profile, pooled ribosome-free and polysome-bound RNA fractions were sequenced using the Illumina platform, yielding nearly 30 million reads for each sample. Approximately 70–80% of the reads were mapped onto the reference genome (GRCh38), and more than 17,000 genes were detected in each type of fraction (Additional file 2: Figure S1). As a control, D15 cells were also treated with puromycin to disassemble the polysomes and cardiomyocyte markers were evaluated by qPCR (Additional file 2: Figure S2).
Correspondence analysis (COA) showed that samples were grouped according to the type of RNA fraction (ribosome-free vs. polysome-bound) (Fig. 1e) and according to the day of differentiation (D0, D1, D4, D9 and D15) (Fig. 1f and g). Polysome-bound samples showed more distinct groups relative to the day of differentiation, indicating high similarity between translated genes in experimental replicates (Fig. 1g). On the other hand, ribosome-free samples showed more dispersion and less similarity between the replicates.
The translatome delineates cardiomyogenic gene expression
During heart development, temporal gene expression changes occur to define each step of cardiogenic commitment [6, 30]. RPKM values (reads per kilobase per million mapped reads) for polysome-bound RNA fraction showed the expression levels of known lineage-specific genes throughout cardiac differentiation (Fig. 2a) which were confirmed by qPCR (Additional file 2: Figures S2 and S3). Pluripotency marker genes, including SOX2, POU5F1 (OCT4) and NANOG, were expressed at higher levels on D0 and D1 and down-regulated at following time-points. The mesoderm marker genes T and EOMES, and early cardiac gene MESP1 were highly specifically expressed on D4. Expression of cardiac-related genes such as GATA4, NKX2–5 and TBX5 as well as cardiomyocyte-specific sarcomeric genes TNNI, TNNT, MYL7 and MYH6 were observed from D9 to D15. In contrast, the endoderm (PECAM1 and PDX1) and ectoderm (PAX6 and FOXP2) marker genes did not change their polysome association during cardiomyogenic differentiation (Fig. 2a).
Comparisons between each differentiation time-point and the preceding time-point, considering an overall FDR of ≤0.05 and − 2 ≥ logFC ≥2, identified differentially expressed genes (DEGs) in polysome-bound RNA fractions (Fig. 2b, data available in Additional file 3). Aggregation of embryoid bodies during the first 24 h of differentiation induced differential expression of 288 genes. Mesoderm commitment from D1 to D4 showed 1264 DEGs, and cardiac progenitor progression to D9 showed 1582 DEGs. The final step of differentiation analysis on D15 showed 743 DEGs compared to D9. In general, the majority of DEGs were up-regulated, except for D15 compared to D9, where the number of down-regulated genes was slightly higher than up-regulated genes (Fig. 2b). Similar numbers of genes and patterns of up- and down-regulated genes were shown in the ribosome-free samples (Additional file 2: Figure S4A and Additional file 4). Regarding the protein-coding and non-coding genes in polysome-bound samples, approximately 20% of DEGs in each time-point analysis are annotated as non-coding RNAs, in which 39% correspond to “lincRNA”, 22% to “antisense” and 18% to “processed pseudogene” RNAs (Additional file 2: Figure S4B).
Gene Ontology (GO) analysis of polysome-bound up-regulated DEGs (FDR ≤ 0.05 and logFC ≥2) revealed developmental and cardiac-related “biological process” (BP) during cardiac differentiation (Fig. 2c). BP terms p-values (−log10) of DEGs for each time-point compared to the preceding one are represented in Fig. 2c heatmap. “Cardiac muscle tissue morphogenesis” and “regulation of muscle contraction” are highly enriched on D9 (compared to D4) and only slightly enriched on D15 (compared to D9), indicating that on D9, most cardiac characteristics are already committed. This pattern is similar for “muscle system process” and “muscle contraction”. “Extracellular matrix organization” seems to have important roles in two distinct phases: mesoderm commitment (D4) and cardiomyocyte final differentiation (D15) (Fig. 2c and Additional file 2: Figure S5).
To assess more information about variations in gene expression during cardiac differentiation, we performed gene clustering using logCPM (counts per million mapped reads). Distinct pattern expression groups were shown: genes with decreased expression during the differentiation, called pluripotency-related cluster and enriched in early developmental BP terms, such as “anterior/posterior axis specification” and “BMP signaling pathway”; and genes with increased expression during cardiac differentiation, called cardiac-related cluster and enriched in lineage specific commitment terms, such as “muscle tissue development” and “regulation of heart contraction” (Fig. 2d).
Strong gene expression coordination is observed during mesoderm-to-cardiac progenitor commitment
When considering polysome-bound RNAs, the largest gene expression variation showing 1582 DEGs occurred during the D4 to D9 shift, which represents mesoderm-to-cardiac progenitor commitment (Fig. 2b). GO analysis revealed that some of D4 up-regulated BP terms were also enriched in the D9 down-regulated analysis (Fig. 3a and Additional file 2: Figure S5), such as “pattern specification process” and “embryonic morphogenesis”. These findings suggest a crucial gene expression regulation at this stage. Comparisons between D4 up-regulated and D9 down-regulated genes showed 217 in common, which are related to “mesoderm development” and “embryonic pattern specification” BP terms (Fig. 3b).
Moreover, an expression gene clustering pattern showed a distinct group of genes highly and specifically expressed on D4 (Fig. 3c). Those genes are called mesoderm-related genes and are enriched on developmental processes, such as “pattern specification process”, “regionalization” and “somitogenesis”. Altogether, these results indicate that the developmental progress of mesoderm-to-cardiac progenitor is carefully regulated and can be assessed by polysome-bound RNA analysis.
Cardiac commitment is intensely tuned by differential mRNA association with polysomes
Gene expression fluctuations could be a consequence of coordination or lack thereof between transcription and translation rate changes. To assess the post-transcriptional regulation during the cardiomyogenic differentiation, we performed the two step-analysis: (1) DEGs were identified through the comparison between each differentiation time-point and the preceding time-point, considering an overall FDR ≤ 0.05, − 2 ≥ logFC ≥2 and RPKM > 1 on ribosome-free or polysome-bound samples (Additional files 3 and 4, respectively); and (2) ribosome-free DEGs and polysome-bound DEGs were compared and classified according to the following categories (Additional file 5). Genes that were up- or down-regulated in both fractions were labeled “up-coordinated” or “down-coordinated”. Transcripts that were up-regulated in ribosome-free fraction but might be neutralized by post-transcriptional mechanisms were labeled “up-buffered”, or in the opposite case, as “down-buffered”. Moreover, genes showing that their polysome association was increased or decreased, were labeled “up-loaded” or “down-loaded”, respectively (Fig. 4a).
We used the coordinated, buffered and loaded classification of DEGs and included one more label category considering the final gene expression as coordinated, post-transcriptional positive or negative regulation (Fig. 4a). Therefore, when compared to the preceding time-point, genes up- or down-coordinated were labeled being under coordinated regulation, genes up-buffered and down-loaded under post-transcriptional negative regulation and genes down-buffered and up-loaded under post-transcriptional positive regulation. Approximately 60–80% of DEGs showed one or another kind of post-transcriptional regulation, positive or negative, suggesting a crucial role of this level of gene expression control (Fig. 4b). Interestingly, during the initial steps of differentiation (D0-D1 and D1-D4), there was a prevalence of up- (10/19 and 142/271) or down-loaded (60/138 and 190/351) genes, suggesting a strong post-transcriptional regulation at these stages. During cardiac progenitor commitment on D9, numbers of coordinated, buffered and loaded genes were similar, either on positive or negative regulation. On the other hand, between D9 and D15 most genes were classified as buffered, in particular, on positive regulation (234/370), indicating that transcriptional variations might be controlled by post-transcriptional mechanisms. These results corroborate our previous findings showing that most of cardiac characteristics are already committed on D9 (Fig. 2).
Considering that Gene Ontology-annotated genes are usually protein-coding, we filtered our data and performed GO analysis using only protein-coding genes. Analysis of coordinately regulated genes showed well-established pathways during hESC differentiation, such as the Reactome pathway terms “POU5F1 (OCT4), SOX2, NANOG repress genes related to differentiation”, “Transcriptional regulation of pluripotent stem cells” and “Developmental Biology”, which were upregulated on D1 and D4 and down-regulated later (D9 and D15) (Additional file 2: Figure S6). Additionally, the cardiac-related pathway “Muscle contraction” and cardiac-specific pathway “Cardiac conduction” were coordinately up-regulated at D9 and D15 time-points (Additional file 2: Figure S6).
Genes classified as buffered or loaded showed enriched pathway terms with a diversity of biological processes (Fig. 5). The Reactome pathways “Developmental Biology” and “NCAM signaling for neurite out-growth” were up-loaded on D4 and down-loaded on D9, once more suggesting the critical regulation at this stage. For instance, developmental-related genes, such as NOTUM, CER1 and SOX17, appeared as up-loaded on D4 and down-loaded on D9 (Additional file 5 and Additional file 2: Figure S7A), indicating the polysomal loading regulation. The “M phase” and “Mitotic Metaphase and Anaphase” terms were shown as down-loaded on D15, while “Cyclin A/B1 associated events during G2/M transition” as down-buffered, indicating the fine adjustment of cell cycle during differentiation. The E2F Transcription Factor 1 plays a crucial role in the control of cell cycle  and was shown as involved in myoblast proliferation and differentiation through the auto-regulation loop with miR-20a-5p and miR-20b-5p . E2F1 gene appeared less associated to polysomes on D9 (FDR < 0.05, logFC − 1.73, not included on down-loaded group), and down-buffered on D15, illustrating an initial polysomal dissociation followed by transcriptional down regulation of this gene (Additional files 3 and 4).
We have previously shown that “Extracellular matrix organization” seems to have an important role during differentiation (Fig. 2c), and it is probably also being post-transcriptionally regulated once it showed as enriched for up-loaded genes on D4, D9 and D15 (Fig. 5). Between the D15 up-loaded genes are LUM, COL6A3 and COL3A1 (Additional file 5). COL3A1 was already shown as post-transcriptionally regulated by the interaction of the heterogeneous nuclear ribonucleoprotein (hnRNP) A1 mRNA-binding protein with its 3′-UTR . Interestingly, the cardiac-specific NKX2–5 gene was shown as up-buffered on D4, suggesting that its RNA was transcribed but not translated yet. On D9, this gene was shown as up-coordinated, suggesting, in this case, its transcription and translation (Additional file 2 and Additional file 5: Figure S7B). In addition, other crucial cardiac transcription factors, such as MEF2A and TBX5 appeared as up-loaded on D9, indicating their polysome-bounding increase on this stage (Additional file 5 and Additional file 2: Figure S7C).
To further characterize how recruitment to and dissociation from ribosomes influence gene expression, we performed polysome/ribosome-free ratio analysis using RPKM values. The results showed genes affected by polysome recruitment (FDR ≤ 0.05, logFC ≥2) or dissociation (FDR ≤ 0.05, − 2 ≥ logFC) during cardiomyocyte differentiation (Additional file 6). Some of them are also DEGs in polysome-bound or ribosome-free fractions, but most are not differentially expressed (Fig. 6a and Additional file 2: Figure S8). GO analysis for these non-differentially expressed and ratio-variated genes showed a variety of BP terms (Fig. 6b and Additional file 2: Figure S8). For instance, the development-related pathways JUN, Wnt and Notch were strongly regulated between D1, D4 and D9, were recruited from D1 to D4 and dissociated from D4 to D9 (Fig. 6b). The polysome/ribosome-free ratio of the JUN, Wnt and Notch pathway genes MTCH1, GALNT11, NCLN and TMEM237 were plotted to visualize the variations on D1, D4 and D9 as an example (Fig. 6c).
mRNA loading into polysomes fine-tunes crucial processes during hESC cardiomyogenesis
To better understand the changes in translation between hESC and cardiomyocytes, we performed Gene Ontology (GO) analysis with DEGs of D0 (hESC) vs. D15 (cardiomyocytes), considering FDR ≤ 0.05, − 1 ≥ logFC ≥1 and ribosome-free and polysome-bound data combined (Additional file 7). Genes down-regulated on D15 were enriched in Biological Process (BP) terms such as “rRNA processing”, “tRNA aminoacylation for protein translation” and “cytoplasmic translation”, which were grouped as “RNA-related terms” (Fig. 7a). Combining the genes annotated in this group (143 genes down-regulated on D15) and analyzing their post-transcriptional regulation classification, 44 (30.7%) of them were down-coordinated and 72 (50.3%) were down-loaded (Fig. 7b and Additional file 7). In addition, many ribosomal proteins were shown down-regulated on D15, predominantly down-loaded (90.9%) (Fig. 7c). Other translation machinery proteins were also down-regulated after cardiomyocyte commitment, for instance, the initiation factors EIF5AL1 and EIF4E1B, and elongation factor EEF1E1 (Additional file 7). On the other hand, the cardiac elongation factor EEF1A2  was up-regulated on D15 (Additional file 7).
Moreover, to explore if non-differentially expressed genes could have been affected by the variation in polysome occupancy (recruitment vs. dissociation), we also compared the polysome/ribosome-free RPKM ratio between hESCs (D0) and cardiomyocytes (D15) (Additional file 6). Interestingly, genes related to translation processes also showed polysomal dissociation on D15 when compared to D0, illustrated by the decreased polysome/ribosome-free ratio (Fig. 7d and e). Among them, there are translation initiation factors (EIF4A3, EIF4E, EIF4B), ribosomal proteins (RPL6, RPL14) and RNA helicases (DDX52), whose ratios are plotted in Fig. 7f (Additional file 6). Some of these observations were confirmed by qPCR (Additional file 2: Figure S9).
To further confirm that the down-regulation of translation-related genes after cardiomyogenic differentiation could affect protein synthesis, we performed a protein synthesis quantification assay. Cells on D0 (hESC) and D15 (cardiomyocytes) were treated with O-propargyl-puromycin (OPP) which is incorporated into newly translated proteins and then fluorescently labeled. Quantification of fluorescence intensity showed a decrease in protein synthesis after cardiac commitment compared to undifferentiated cells (Fig. 7g and h). Taken together, these findings suggest a translation adjustment during hESC-to-cardiomyocyte differentiation.
Interestingly, other crucial cellular processes also showed regulation by post-transcriptional mechanisms. For instance, on D1 down-buffered and D15 up-buffered GO analysis, the Reactome pathway terms related to cellular metabolism were enriched, suggesting an important post-transcriptional regulation of this process during cardiomyogenic differentiation. Metabolic properties differ between cardiomyocytes and hESCs , thus, to understand the regulation of metabolic genes, we grouped the genes annotated as the Reactome pathway terms related to cellular metabolism and called them “Metabolism-related genes” (Additional file 2: Figure S10A and Additional file 7). Comparing the ribosome-free and polysome-bound data, 494 genes were up-regulated on D15, among which 243 (49.2%) were up-loaded, 46 (9.3%) were up-buffered and 205 (41.5%) were up coordinated (Additional file 2: Figure S10B). This likely reflects the post-transcriptional contribution to cardiomyocyte metabolic remodeling demonstrated by metabolic gene recruitment to polysomes.
In this report, we differentiated hESC to cardiomyocytes and analyzed distinct time-points during this process to assess the temporal cell fate transition. The polysome profiling approach followed by ribosome-free and polysome-bound RNA-seq allowed us to evaluate gene regulation during the cardiogenic commitment. Polysome profiling analysis has been used as a robust method to assess the association of ribosomes with mRNAs, providing information about their translational status [23,24,25]. Here, we show that polysome-bound RNAs reflect the cardiac commitment phenotype, illustrated by down-regulation of pluripotency core regulatory circuitry (OCT4, SOX2 and NANOG)  followed by up-regulation of cardiomyogenesis-related genes. The recapitulation of developmental steps is a powerful strategy to control a specific cell fate , where the first step is the transition into one of the three embryonic germ layers. The heart originates from the mesoderm emerging from the primitive streak . Mesodermal development genes were found to be strongly regulated on D4, which represents the cardiac mesoderm commitment time-point. The massive number of DEGs between D4 and D9, added to the up-regulation of developmental pathways and pattern specification genes on D4 followed by their down-regulation at progenitor specification stage (D9) are consistent with the complexity of multiple mesodermal lineage choices, recently mapped by Loh et al. (2016) .
Moreover, we showed that 60–80% of DEG along cardiomyogenic differentiation were under some degree of post-transcriptional regulation. Uncoupling between the transcriptome and translatome changes, characterized by the independency of the machineries responsible for mRNA availability and engagement in translation, has been observed in most cell types . Here, we observed that many metabolic and cellular processes are regulated exclusively at the level of polysomal association, probably because of post-transcriptional regulatory mechanisms in cellular development [25, 39, 40]. Surprisingly, genes with coordinated regulation (additive regulation) represent a minority of DEGs. In addition to that, changes in the abundance of mRNAs do not directly influence their polysomal engagement, represented by the buffered expression genes. The mRNA access to the translation machinery is governed by the interaction of translation factors with its sequence and structural motifs and is mediated by a complex network of trans-acting proteins and regulatory RNAs [41, 42]. While mRNA levels in the cytoplasm can randomly fluctuate due to transcriptional bursts, the final gene expression levels are adjusted at the time of protein synthesis . McManus et al. (2013)  showed that buffering effects on gene expression (translation efficiency opposite to mRNA abundance) were common and would reduce divergent interspecies expression at the protein level. Translational regulation has already been reported in mouse mesoderm commitment  and hESC neuronal differentiation . Here, we show unprecedented data about post-transcriptional regulation during human cardiomyogenic commitment.
We also found that the translational activity is reduced during differentiation; this became more evident when comparing hESC on D0 to cardiomyocytes on D15. Protein synthesis rate in cardiac tissue decreases from fetal to adult development, and in adult heart is lower than in other tissues [47, 48]. Chorghade et al. (2017) showed that the translation initiation RNA-binding protein PABPC1 is post-transcriptionally down-regulated in the adult heart tissue, which possibly decreases the translational capacity of the heart. Our data showed that the D15 down-regulated genes were enriched in RNA processing and translation genes, which were mostly down-regulated only in the polysome-bound fraction. Translation-related genes also showed polysomal dissociation on D15 when compared to D0, illustrated by the decreased polysome/ribosome-free ratio, further suggesting a strong post-transcriptional regulation and corroborating previous observations [38, 49, 50].
Protein synthesis is an indispensable process for maintaining homeostasis in the cell, and aberrations in its regulation may contribute to a wide range of diseases . We have previously investigated the translational regulation during cell commitment to adipogenesis in hASCs. Our results showed the translational control as a key mechanism regulating the early steps of adipogenic differentiation, with a significant reduction in protein synthesis and a lower translational efficiency of ribosomal proteins . In vitro differentiation of other cell types also showed decreasing translation rate, due to, e.g., eIF2a phosphorylation or rRNA biogenesis [53,54,55,56,57]. Blair et al. (2017) reported the translational down-regulation of translation-related genes after hESC neuronal differentiation (including translation initiation factors and ribosomal proteins) , indicating that this mechanism doesn’t seem to be cardiac-specific. In contrast, Blanco et al. (2016) showed that adult stem cells have lower protein synthesis rates than committed cells, using skin as a model. In normal skin, the RNA methyltransferase NSUN2 expression is restricted to committed hair follicle populations. They have demonstrated that loss of NSUN2 causes hypomethylation of tRNAs, accumulation of 5′ tRNA fragments, which repress cap-dependent protein translation [58,59,60]. This reduced translation in adult stem cells is probably related to the biology of these specific cell types. Adult stem and progenitor cells are metabolically quiescent and low, if at all, proliferative. After stimulation by injury, they enter an activation state where they proliferate and show increased translational rates . These particular characteristics make the comparison of the regulatory pathways and mechanisms underlying translational regulation in adult and embryonic stem cells not possible as these cells are in different stages of differentiation and in completely different biological environments. Our model recapitulates embryonic differentiation, where highly proliferative stem cells commit and differentiate into cell types with defined gene expression patterns and, in most cases, lower proliferative rates.
Noncontractile stem cells require less energy than beating cardiomyocytes, therefore, a transition in energetic infrastructure is necessary to support the increased energetic needs during cardiac differentiation [62, 63]. This switch in energy metabolism is associated with increased mitochondrial maturation and oxygen consumption and reduced glycolysis, due to a metabolic transcriptome remodeling [35, 62]. We showed that the metabolic reconfiguration is also a consequence of post-transcriptional level gene regulation. Almost 50% of metabolism-related genes regulated on D15 were increased only in polysome-bound fraction, characterizing the higher recruitment of these genes to the translational machinery. Increasing evidence has been showing the role of translational control in regulating metabolic function and that its defect is implicated in the pathogenesis of metabolic disorders [64, 65]. Additionally, cardiomyocyte mRNA interactome revealed that metabolic enzymes can act as RBPs , and energy metabolism proteins were found associated with ribosomes in ESCs, for instance, controlling the translation of mRNAs . All these lines of evidence suggest a close relationship between cellular metabolism and translational control which deserves further exploration.
Polysome profiling followed by ribosome-free and polysome-bound RNA-seq of hESC cardiogenic commitment allowed us to evaluate gene expression regulation, revealing a massive control of developmental steps. Post-transcriptional regulation might be acting on 60–80% of DEG along cardiomyogenic differentiation, contributing to the final phenotype. Cardiomyocytes translational activity was reduced when compared to undifferentiated cells showing that translation regulation is also controlled during differentiation. Understanding the biological processes that trigger differentiation is crucial for an improvement of in vitro differentiation strategies and discovery of endogenous regenerative pathways. Our results helped to uncover transcriptional and post-transcriptional regulation involved in cardiomyogenic commitment gene expression, collected unprecedented data regarding polysome-bound RNAs and provided advances in comprehension of differentiation processes. Altogether, this work could prospect fundamental tools to develop new therapy and research approaches.
Cell culture and cardiomyocyte differentiation
hES-NKX2–5eGFP/w cell line was a gift from Dr. David A. Elliott (Monash Immunology and Stem Cell Laboratories, Monash University, Australia) . Cells were maintained on irradiated MEFs (mouse embryonic fibroblasts) using DMEM-F12 (Gibco™) with 20% KSR (Gibco™), 0.1 mM non-essential amino acids, 55 μM β-mercaptoethanol, 100 μg/mL penicillin-streptomycin, 2 mM L-glutamine and 10 ng/mL bFGF. The cardiomyocyte differentiation protocol was adapted from the previously published one . Seventy-two hr. before differentiation, 7 × 105 cells/well were passaged onto a 6-well plate coated with Growth Factor Reduced BD® Matrigel Matrix for feeder depletion. Then, cells were dissociated with collagenase I (1 mg/mL) for 20 min followed by trypsin-EDTA (0.05%) for approximately 30 s and scraped to form small clusters. After washing, the clusters were cultured in 6-well Ultralow Attachment Corning plates using StemPro-34 (Invitrogen), supplemented with 100 μg/ml penicillin-streptomycin, 2 mM L-glutamine, 150 μg/mL transferrin, 50 μg/mL ascorbic acid and 0.45 mM monothioglycerol (MTG) and kept in humid incubator at 37 °C, 5% CO2 and 5% O2. For embryoid bodies (EBs) aggregation (D0-D1), the basal medium was supplemented with 1 ng/mL BMP4. On D1, the medium was replaced with supplementation of 10 ng/mL BMP4, 6 ng/mL Activin A and 5 ng/mL bFGF; on D4 with 10 μg/mL VEGF and 10 μM XAV 939, and starting on D8, the medium was replaced every two or three days supplemented with 10 μg/mL VEGF and 1 ng/mL BMP4. Three independent differentiation assays were used as experimental replicates.
EBs were dissociated on D4 using trypsin-EDTA (0.05%) and incubated with the surface marker PE-conjugated anti-CD56 (1:25 in 0.5% PBS/BSA, BD cat. 347,747) and 1 μg/μL DAPI. On D9, cells were disaggregated with trypsin-EDTA (0.05%) for 5 min and resuspended in PBS to evaluate eGFP expression. On D15, EBs were disaggregated using 1 mg/mL collagenase I for 16 h and trypsin-EDTA (0.05%) for 5 min, fixed with 4% paraformaldehyde, permeabilized with 0.5% Triton X-100 and incubated with anti-troponin T antibody (1:100 in 0.5% PBS/BSA, cardiac isoform Ab-1, Thermo Scientific™, cat. #MS-295-P0) followed by Pacific Blue-conjugated anti-mouse antibody (11000). Analyses were carried out using a FACSCanto II flow cytometer and FlowJo software.
Immunofuorescence and fluorescent microscopy
On D15, EBs were disaggregated using 1 mg/mL collagenase I for 16 h and trypsin-EDTA (0.05%) for 5 min and plated on Matrigel-coated wells. After 2–5 days, cardiomyocytes were fixed with 4% paraformaldehyde, permeabilized and blocked with 0.5% Triton X-100 and 1% PBS-BSA. Overnight incubation with anti-troponin I antibody (1:100 in 0.5% PBS/BSA, Santa Cruz Biotechnology, cat.: sc-15,368) followed by Alexa 546-conjugated anti-rabbit IgG (1:800, Invitrogen) and 1 μg/μL DAPI. EBs or fixed cardiomyocytes were visualized using a Leica DMI6000B optical microscope and images and videos acquired by LAS AF software.
Polysome profile and RNA isolation
At the indicated time points, monolayer hESC on D0 or differentiating EBs were treated with 0.1 mg/mL cycloheximide (Sigma-Aldrich) for 10 min at 37 °C, disaggregated with trypsin-EDTA (0.05%) for 10 min and washed twice with PBS. Cells were resuspended in polysome lysis buffer (15 mM Tris HCl, pH 7.4, 15 mM MgCl2, 300 mM NaCl, 1% Triton X-100, 40 U/μL RNAse Out, 24 U/mL DNAse and 100 μg/mL cycloheximide), incubated for 10 min on ice and centrifuged at 12000 x g for 10 min at 4 °C. For the puromycin control, cells were treated with 2 mg/mL puromycin for 1 h before disaggregated and the buffers were used without cycloheximide. The supernatants were loaded onto 10 to 50% sucrose gradients (prepared with BioComp model 108 Gradient Master) and centrifuged at 150000 x g (SW40 rotor, HIMAC CP80WX HITACHI) for 160 min at 4 °C. Different sucrose gradient fractions were separated using ISCO gradient fractionation system (ISCO Model 160 Gradient Former Foxy Jr. Fraction Collector), connected to a UV detector, which monitored the absorbance at 275 nm to record the polysome profile. RNA from ribosome-free and pooled polysomal fractions was isolated using the Direct-zol RNA MiniPrep (Zymo Research), following the manufacturer’s instructions.
High-throughput sequencing and data analysis
For cDNA library preparation, 200 to 500 ng of ribosome-free or 2 μg of polysome-bound RNA were used to perform three independent sample replicates. The cDNA libraries were prepared using the TruSeq Stranded mRNA Sample Preparation kit (Illumina, Inc.), and RNA-seq was carried out in an Illumina HiSeq platform. Mapping and counting of sequencing data was performed with the Rsubread package  against the new version of the human genome GRCh38. Mapping parameters were set for unique mapping of the reads (the rest were default). For some calculations, CPM values (counts per million) were determined, where each sample was normalized to one million reads to account for the library size (Additional file 8).
To assess the quality of the experiment and reproducibility of results, we performed a correspondence analysis (COA), a dimension reduction method of the matrix of counts. In COA, it is possible to simultaneously visualize samples and genes, revealing associations between them.
Differential expression analysis was done using the Bioconductor R package edgeR . Several comparisons were performed for both ribosome-free and polysome-bound RNA fractions - each sample against the preceding time-point: D0 vs D1, D1 vs D4, D4 vs D9 and D9 vs D15; and D0 vs D15. For these analyses, we retained only those genes with at least one count per million in at least three samples. After a normalization procedure using three recommended methods (estimateGLMCommonDisp, estimateGLMTrendedDisp, estimateGLMTagwiseDisp), differential expression analysis for all 10 comparisons was performed using the generalized linear mixed model (glmFit and glmLRT). Correction for multiple testing was performed with FDR and RPKM values for each sample were also determined to compare the expression between the samples (Additional file 3, Additional file 4 and Additional file 7).
Clustering was performed with the k-means algorithm on the table of read log-counts of each gene. Biological replicates were averaged (mean of 3 replicates) and only differentially expressed genes in either condition (D0 vsD1, D1 vs D4, D4vs D9, D9 vs D15) were considered. Different numbers of k were tested by measuring the within sum of squares. After visual inspection of within sum of squares distribution, we decided on k = 9.
Polysome/ribosome-free ratio was determined by dividing RPKM values derived from the polysomal fraction by the ones derived from the ribosome-free experiment (three biological replicates were averaged).
Protein synthesis was measured using Click- iT® Plus OPP Protein Synthesis Assay (Molecular Probes, Grand Island, NY). Cells on D0 (hESC) or D15 (cardiomyocytes) were seeded in matrigel-coated 96-well plates, and after 48 h, the staining and detection was performed following the manufacturer’s instructions. The quantitative analysis was performed using an Operetta HTS imaging system (PerkinElmer, Waltham MA, USA). Images of 25 fields per well were evaluated with Harmony 3.5.2 software (PerkinElmer). Fluorescence intensities were measured at the cytoplasm regions around the nucleus. For each condition, 1400 cells were randomly chosen for intensity analysis.
Counts per million mapped reads
Differentially expressed genes
False discovery rate
Human adipose stem cells
Human embryonic stem-cell
Reads per kilobase per million mapped reads
Murry CE, Keller G. Differentiation of embryonic stem cells to clinically relevant populations: lessons from embryonic development. Cell. 2008;132:661–80.
Kattman SJ, Witty AD, Gagliardi M, Dubois NC, Niapour M, Hotta A, et al. Stage-specific optimization of activin/nodal and BMP signaling promotes cardiac differentiation of mouse and human pluripotent stem cell lines. Cell Stem Cell. 2011;8:228–40. https://doi.org/10.1016/j.stem.2010.12.008.
Kehat I, Kenyagin-Karsenti D, Snir M, Segev H, Amit M, Gepstein A, et al. Human embryonic stem cells can differentiate into myocytes with structural and functional properties of cardiomyocytes. J Clin Invest. 2001;108:407–14. https://doi.org/10.1172/JCI200112131.
Laflamme MA, Chen KY, Naumova AV, Muskheli V, Fugate JA, Dupras SK, et al. Cardiomyocytes derived from human embryonic stem cells in pro-survival factors enhance function of infarcted rat hearts. Nat Biotechnol. 2007;25:1015–24.
Burridge PW, Keller G, Gold JD, Wu JC. Production of de novo cardiomyocytes: human pluripotent stem cell differentiation and direct reprogramming. Cell Stem Cell. 2012;10:16–28. https://doi.org/10.1016/j.stem.2011.12.013.
den Hartogh SC, Wolstencroft K, Mummery CL, Passier R. A comprehensive gene expression analysis at sequential stages of in vitro cardiac differentiation from isolated MESP1-expressing-mesoderm progenitors. Sci Rep. 2016;6. https://doi.org/10.1038/srep19386.
Li Y, Lin B, Yang L. Comparative transcriptomic analysis of multiple cardiovascular fates from embryonic stem cells predicts novel regulators in human Cardiogenesis. Sci Rep. 2015;5. https://doi.org/10.1038/srep09758.
Liu Q, Jiang C, Xu J, Zhao MT, Van Bortle K, Cheng X, et al. Genome-wide temporal profiling of transcriptome and open chromatin of early cardiomyocyte differentiation derived from hiPSCs and hESCs. Circ Res. 2017;121:376–91.
Tompkins JD, Jung M, Chen CY, Lin Z, Ye J, Godatha S, et al. Mapping Human Pluripotent-to-Cardiomyocyte Differentiation: Methylomes, Transcriptomes, and Exon DNA Methylation “Memories”. EBioMedicine. 2016;4:74–85. https://doi.org/10.1016/j.ebiom.2016.01.021.
King HA, Gerber AP. Translatome profiling: Methods for genome-scale analysis of mRNA translation. Brief Funct Genomics. 2014;15:22–31.
Kolle G, Shepherd JL, Gardiner B, Kassahn KS, Cloonan N, Wood DLA, et al. Deep-transcriptome and ribonome sequencing redefines the molecular networks of pluripotency and the extracellular space in human embryonic stem cells. Genome Res. 2014;21:2014–25.
Sampath P, Pritchard DK, Pabon L, Reinecke H, Schwartz SM, Morris DR, et al. A hierarchical network controls protein translation during murine embryonic stem cell self-renewal and differentiation. Cell Stem Cell. 2008;2:448–60.
Zych J, Spangenberg L, Stimamiglio MA, Abud APR, Shigunov P, Marchini FK, et al. Polysome profiling shows the identity of human adipose-derived stromal/stem cells in detail and clearly distinguishes them from dermal fibroblasts. Stem Cells Dev. 2014;23:2791–802. https://doi.org/10.1089/scd.2013.0496.
Elliott DA, Braam SR, Koutsis K, Ng ES, Jenny R, Lagerqvist EL, et al. NKX2-5eGFP/w hESCs for isolation of human cardiac progenitors and cardiomyocytes. Nat Methods. 2011;8:1037–40. https://doi.org/10.1038/nmeth.1740.
Yang L, Soonpaa MH, Adler ED, Roepke TK, Kattman SJ, Kennedy M, et al. Human cardiovascular progenitor cells develop from a KDR+ embryonic-stem-cell-derived population. Nature. 2008;453:524–8. https://doi.org/10.1038/nature06894.
Thiele BJ, Doller A, Kähne T, Pregla R, Hetzer R, Regitz-Zagrosek V. RNA-binding proteins heterogeneous nuclear ribonucleoprotein A1, E1, and K are involved in post-transcriptional control of collagen I and III synthesis. Circ Res. 2004;95:1058–66.
Svobodová K, Horák P, Stratil A, Bartenschlager H, Van Poucke M, Chalupová P, et al. Porcine EEF1A1 and EEF1A2 genes: genomic structure, polymorphism, mapping and expression. Mol Biol Rep. 2015;42:1257–64.
Loh KMM, Chen A, Koh PWW, Deng TZZ, Sinha R, Tsai JMM, et al. Mapping the pairwise choices leading from pluripotency to human bone, heart, and other mesoderm cell types. Cell. 2016;166:451–68. https://doi.org/10.1016/j.cell.2016.06.011.
Tebaldi T, Re A, Viero G, Pegoretti I, Passerini A, Blanzieri E, et al. Widespread uncoupling between transcriptome and translatome variations after a stimulus in mammalian cells. BMC Genomics. 2012;13:220. https://doi.org/10.1186/1471-2164-13-220.
Fujii K, Shi Z, Zhulyn O, Denans N, Barna M. Pervasive translational regulation of the cell signalling circuitry underlies mammalian development. Nat Commun. 2017;8:1–13. https://doi.org/10.1038/ncomms14443.
Chorghade S, Seimetz J, Emmons R, Yang J, Bresson SM, De Lisio M, et al. Poly(a) tail length regulates PABPC1 expression to tune translation in the heart. Elife. 2017;6:e24139. https://doi.org/10.7554/eLife.24139.
Marcon BH, Holetz FB, Eastman G, Origa-Alves AC, Amorós MA, de Aguiar AM, et al. Downregulation of the protein synthesis machinery is a major regulatory event during early adipogenic differentiation of human adipose-derived stromal cells. Stem Cell Res. 2017;25:191–201. https://doi.org/10.1016/j.scr.2017.10.027.
You KT, Park J, Kim VN. Role of the small subunit processome in the maintenance of pluripotent stem cells. Genes Dev. 2012;29:2004–9.
Blanco S, Bandiera R, Popis M, Hussain S, Lombard P, Aleksic J, et al. Stem cell function and stress response are controlled by protein synthesis. Nature. 2016;534:335–40. https://doi.org/10.1038/nature18282.
Rodgers JT, King KY, Brett JO, Cromie MJ, Charville GW, Maguire KK, et al. MTORC1 controls the adaptive transition of quiescent stem cells from G0to GAlert. Nature. 2014;510:393–6. https://doi.org/10.1038/nature13255.
Lopaschuk GD, Jaswal JS. Energy metabolic phenotype of the cardiomyocyte during development, differentiation, and postnatal maturation. J Cardiovasc Pharmacol. 2010;56:130–40. https://doi.org/10.1097/FJC.0b013e3181e74a14.
Adeli K. Translational control mechanisms in metabolic regulation: critical role of RNA binding proteins, microRNAs, and cytoplasmic RNA granules. AJP Endocrinol Metab. 2011;301:E1051–64. https://doi.org/10.1152/ajpendo.00399.2011.
The authors would like to thank the Program for Technological Development in Tools for Health-RPT-FIOCRUZ for the use of flow cytometry and microscopy facilities at Carlos Chagas Institute – Fiocruz/PR, the National Cancer Institute (Rio de Janeiro, Brazil) for the use of sequencing facility and Dr. David Elliot from Monash University (Australia) for kindly providing the cells.
Financial support from CNPq: INCT-REGENERA grant 465656/2014–5. ITP and AWR received fellowships from CAPES, and BD from CNPq.
Availability of data and materials
The datasets generated and/or analysed during the current study are available in the Short Read Archive, NCBI (https://www.ncbi.nlm.nih.gov/sra/SRP150416). This site serves as a landing page for the study: description of the project, metadata and raw sequencing files can be found there. Counts data and RPKM can be found in Additional file 8. One tab corresponds to the read counts of each sample, and the other, to RPKM values. Each column of each file is labeled with the sample condition, e.g., 1D4P corresponds to biological replicate 1, at day 4 and polysomal RNA fraction.
Authors and Affiliations
Basic Stem-cell Biology Laboratory, Instituto Carlos Chagas - FIOCRUZ-PR, Rua Professor Algacyr Munhoz Mader, 3775, Curitiba, PR, 81.350-010, Brazil
Isabela Tiemy Pereira, Anny Waloski Robert, Marco Augusto Stimamiglio & Bruno Dallagiovanna
Bioinformatics Unit, Institut Pasteur de Montevideo, Mataojo 2020, 11400, Montevideo, Uruguay
ITP and AWR performed the experiments. ITP and BD wrote the manuscript. LS performed the bioinformatic analyses and wrote the manuscript. RAH performed the data mapping. MAS coordinated the cardiac differentiation experiments. HN coordinated the bioinformatic analyses. BD conceived and coordinated the work. All authors read and approved the final manuscript.
Figure S1. Summary of high-throughput ribosome-free and polysome-bound RNA-seq of distinct cardiomyogenic differentiation time-points (n = 3). Figure S2. (A) Polysome profiling of D15 cells treated with cycloheximide or puromycin. Ribosome-free (fractions 1–3), monosome (fractions 5–7), light polysomes (fractions 9–15) and heavy polysomes (fractions 16–22) fractions were pooled and isolated. (B) Cardiomyocyte markers evaluated by qPCR on distinct polysome fractions. Figure S3. Polysome-bound validation of developmental markers expression using qPCR. Figure S4. (A) Number of differentially expressed genes on each differentiation time-point, compared to previous time-point (FDR < 0.05, − 2 > logFC> 2) on ribosome-free samples. Number of protein-coding and non-coding genes are also shown (bottom panel). (B) Non-coding genes categories of DEGs (all time-points combined, each time-point against its previous for analysis) (FDR < 0.05, − 2 > logFC> 2) on polysome-bound samples. Figure S5. Gene Ontology EnrichR BP enriched terms for up (FDR < 0.05, logFC> 2) and down (FDR < 0.05, logFC<− 2) polysome-bound regulated genes. Figure S6. Genes coordinately regulated are under control of transcriptional and post-transcriptional regulation during cardiomyogenic differentiation. Figure S7. qPCR validation of DEG during cardiomyogenic differentiation. Figure S8. Non-DEGs showed differences on polysome recruitment and dissociation on D0 vs. D1 and D9 vs. D15. Figure S9. RNA related-genes validation by qPCR. Figure S10. Cardiomyocytes (D15) showed up-regulation of cellular metabolism genes. (DOCX 2575 kb)
Counts and RPKM values of all mapped genes. (XLSX 30119 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
Pereira, I.T., Spangenberg, L., Robert, A.W. et al. Cardiomyogenic differentiation is fine-tuned by differential mRNA association with polysomes.
BMC Genomics20, 219 (2019). https://doi.org/10.1186/s12864-019-5550-3