{'fraction_non_alphanumeric': 0.029916746280878943, 'fraction_numerical': 0.01965333697284018, 'mean_word_length': 5.437532946758039, 'pattern_counts': {'":': 0, '<': 5, '<?xml version=': 0, '>': 2, 'https://': 0, 'lorem ipsum': 0, 'www.': 2, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
The cardiac translational landscape reveals that micropeptides are new players involved in cardiomyocyte hypertrophy

Hypertrophic growth of cardiomyocytes is one of the major compensatory responses in the heart after physiological or pathological stimulation. Protein synthesis enhancement, which is mediated by the translation of messenger RNAs, is one of the main features of cardiomyocyte hypertrophy. Although the transcriptome shift caused by cardiac hypertrophy induced by different stimuli has been extensively investigated, translatome dynamics in this cellular process has been less studied. Here, we generated a nucleotide-resolution translatome as well as transcriptome data from isolated primary cardiomyocytes undergoing hypertrophy. More than 10,000 open reading frames (ORFs) were detected from the deep sequencing of ribosome-protected fragments (Ribo-seq), which orchestrated the shift of the translatome in hypertrophied cardiomyocytes. Our data suggest that rather than increase the translational rate of ribosomes, the increased efficiency of protein synthesis in cardiomyocyte hypertrophy was attributable to an increased quantity of ribosomes. In addition, more than 100 uncharacterized short ORFs (sORFs) were detected in long noncoding RNA genes from Ribo-seq with potential of micropeptide coding. In a random test of 15 candidates, the coding potential of 11 sORFs was experimentally supported. Three micropeptides were identified to regulate cardiomyocyte hypertrophy by modulating the activities of oxidative phosphorylation, the calcium signaling pathway, and the mitogen-activated protein kinase (MAPK) pathway. Our study provides a genome-wide overview of the translational controls behind cardiomyocyte hypertrophy and demonstrates an unrecognized role of micropeptides in cardiomyocyte biology.

Hypertrophic growth of cardiomyocytes is one of the major compensatory responses in the heart after physiological or pathological stimulation. Protein synthesis enhancement, which is mediated by the translation of messenger RNAs, is one of the main features of cardiomyocyte hypertrophy. Although the transcriptome shift caused by cardiac hypertrophy induced by different stimuli has been extensively investigated, translatome dynamics in this cellular process has been less studied. Here, we generated a nucleotide-resolution translatome as well as transcriptome data from isolated primary cardiomyocytes undergoing hypertrophy. More than 10,000 open reading frames (ORFs) were detected from the deep sequencing of ribosome-protected fragments (Ribo-seq), which orchestrated the shift of the translatome in hypertrophied cardiomyocytes. Our data suggest that rather than increase the translational rate of ribosomes, the increased efficiency of protein synthesis in cardiomyocyte hypertrophy was attributable to an increased quantity of ribosomes. In addition, more than 100 uncharacterized short ORFs (sORFs) were detected in long noncoding RNA genes from Ribo-seq with potential of micropeptide coding. In a random test of 15 candidates, the coding potential of 11 sORFs was experimentally supported. Three micropeptides were identified to regulate cardiomyocyte hypertrophy by modulating the activities of oxidative phosphorylation, the calcium signaling pathway, and the mitogenactivated protein kinase (MAPK) pathway. Our study provides a genome-wide overview of the translational controls behind cardiomyocyte hypertrophy and demonstrates an unrecognized role of micropeptides in cardiomyocyte biology.

INTRODUCTION
The heart undergoes remodeling in response to various physiological and pathological stimuli. 1 Hypertrophic growth of cardiomyocytes, which account for 70% of the volume of the heart, is one of the major cellular responses during cardiac remodeling. 2 Increased protein synthesis and sarcomere assembly occur in this process. Both phosphatidylinositol 3-kinase (PI3K)/AKT and mitogen-activated protein kinase (MAPK) signaling pathways have been widely reported as one of the main signaling pathways transducing external stimuli to a shift in the transcriptome in cardiac hypertrophy. For example, phenylephrine (PE) activates the PI3K/AKT and MAPK signal pathways through the a1-receptor, which induces hypertrophic growth of cardiomyocytes. 3,4 Studies have demonstrated that PE induces protein synthesis in cardiomyocytes through the downstream targets of mammalian target of rapamycin (mTOR), S6K1, and 4EBP1; 5 however, the mechanism of this induction at the translation level has yet to be determined. Whether the increased efficiency of protein synthesis results from an increase in the quantity or the translational rate of ribosomes is unknown.
The development of new technologies for massively parallel sequencing provides new tools for genome-wide investigations of regulation by various DNA and RNA modifications. The shift of the cardiac transcriptome under different conditions has been intensively investigated by our group 6 and others 7,8 in the past decade; yet, transcriptome data deviate from mass spectrum data, which reflect the protein levels, to some extent. Therefore, deep sequencing of actively translated transcripts might provide a more accurate and complete measure of gene expression in addition to conventional RNA sequencing (RNA-seq) analyses. Translating ribosome affinity purification sequencing (TRAP-seq) was previously applied to study the actively translating transcripts associated with polyribosomes in the heart. 9 Recently, the cardiac translatome has been analyzed with a combination of deep sequencing of ribosome-protected fragments (Ribo-seq) and RNA-seq. 10,11 Similar to the cardiac transcriptome analysis, extension of the analysis of the translational landscape in cardiac remodeling under different stimuli may help to obtain a full picture of how the gene program is shifted in cardiac disease.
Micropeptides refer to those translational products containing 100 amino acids (aa) or less derived from short open reading frames (sORFs) in the transcripts. 12 Since the sORFs are too short, most of them are considered as false-positive hits during ORF scanning for genomic annotation. 13 Therefore, genes producing, for example, sORF-containing transcripts, are usually mis-annotated as noncoding RNA genes. Accumulating evidence has demonstrated that a large number of micropeptides are present in mammalian cells with unknown function in most cases. 10,14,15 The biological function of these "tiny" protein molecules needs to be explored because it could help distinguish whether the regulatory function of the DNA locus comes from the "noncoding" transcript or if they are encoded micropeptides. Recently, micropeptides have been shown to regulate striated muscle development and function. 16 Myomixer, a micropeptide enriched in developing and regenerating skeletal muscle, was reported to control the critical step in myofiber formation during muscle development. 17 Furthermore, dwarf open reading frame (DWORF), a sarcoplasmic reticulum-located micropeptide, activates sarcoplasmic endoreticulum calcium ATPase (SERCA), which is a target candidate for gene therapy of heart failure, 18,19 and prevents pathological remodeling and Ca 2+ dysregulation in a model of heart failure. 20 In this study, we performed Ribo-seq and RNA-seq with hypertrophic cardiomyocytes induced by PE. Our transcriptional and translational profiling data, respectively, showed that the PI3K/AKT and MAPK signaling pathways were upregulated translationally and transcriptionally. Our data demonstrated that the increased efficiency of protein synthesis in PE-induced cardiomyocyte hypertrophy was mainly attributable to the increased quantity of ribosomes, as evidenced by the increased translation of ribosomal protein genes. In addition, hundreds of sORF candidates were identified in this genome-wide translational analysis with a portion of the micropeptides coded by transcripts from the complementary strand of protein-coding genes. Three micropeptides were further demonstrated to regulate the hypertrophic growth of cardiomyocytes through modulating various signaling pathways as assessed by a functional assay.

Genome-wide translatomic and transcriptomic profiling in hypertrophic cardiomyocytes
Enhanced protein synthesis is one of the main features of hypertrophic growth of cardiomyocytes under stimuli. 21,22 To study how translation is regulated at a genome-wide scale during the process of hypertrophy, primary neonatal cardiomyocytes were isolated and stimulated with PE for 24 h. Hypertrophic growth was confirmed by measuring the size of cardiomyocytes (Figures S1A and S1B) and detecting the gene expression of markers, such as natriuretic peptide A (NPPA), NPPB, and ACTA1 ( Figure S1C). Cells were then harvested and prepared for Ribo-seq, which detects ribosome-protected fragments for measuring ribosome occupancy of protein-coding genes, and RNA-seq, which allows the study of changes in the transcriptome during cardiomyocyte hypertrophy. In total, we obtained $185 M clean reads from Riboseq ($90 M from control and $95 M from PE-treated cells; Table  S1) and $820 M clean reads from RNA-seq ($400 M from control and $420 M from PE-treated cells; Table S1). The majority of the detected ribosome-protected fragments showed a length distribution from 26 nt to 32 nt ( Figure 1A). Uniquely mapped reads from Riboseq were analyzed with Ribocode, 23 and we found that over 80% of the main ribosome-protected fragments (from 28 nt to 30 nt) were qualified reads ( Figure 1B), which displayed the 3-nt codon movement characteristic of actively translating ribosomes ( Figure 1C). Further analyses with Ribocode found that the majority of qualified reads (81.4%) mapped to ORFs of 13,762 annotated protein-coding genes. 6.5% qualified reads mapped to uORFs, a sORF localizing upstream of the main ORF of protein-coding genes. Interestingly, 3.9% of qualified reads mapped to long noncoding RNA (lncRNA) genes ( Figure 1D; Table  S2), indicating that these genes carried coding potential, which was studied below.
To identify the differentially expressed (DE) genes in transcription and translation, we analyzed RNA-seq and Ribo-seq data from PE-treated and control cardiomyocytes. As a result, we detected 1,201 transcriptionally DE genes (545 upregulated and 656 downregulated with |Log2(fold change)| > 1 and false discovery rate [FDR] < 0.05) ( Figure 2A;  Table S4). The increase of expressed and translating marker mRNAs, such as NPPA, NPPB, and ACTA1, were further confirmed in RNA-seq and Ribo-seq analyses (Figures 2AÀ2C; Tables S3 and S4. Next, we asked which signaling pathways were regulated during PEinduced cardiomyocyte hypertrophic growth. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that translationally dysregulated genes were enriched in several pathways known to regulate cardiac hypertrophy, such as "MAPK signaling pathway" and "PI3K/AKT signaling pathway" (Figures 2D and 2E). We performed a similar KEGG analysis for transcriptionally dysregulated genes. Similar KEGG terms were enriched with MAPK signaling pathway and PI3K/AKT signaling pathway at the top of the list (Figure 2F). Furthermore, enriched ontology clusters analysis of the transcriptionally dysregulated genes showed that terms related to muscle system process, contraction, and development were among the most highly enriched, which were expected in the process of cardiomyocyte hypertrophic growth ( Figure 2G). The similarity of top KEGG lists of the dysregulated gene enrichment between Ribo-seq and RNA-seq provided us a hint that the rate of utilization of the ribosome for translation of the majority of mRNAs did not significantly change in PE-induced cardiomyocyte hypertrophy.
Increased protein synthesis is mainly attributable to an enhanced amount of ribosomal proteins during PE-induced cardiomyocyte hypertrophy Increased protein synthesis during cardiomyocyte hypertrophy may result from an increased number of ribosomes (with stable translational rate) and/or an increased translational rate from these ribosomes. Analysis with a less stringent criteria (FDR < 0.05) provided us a broad view of the translatome shift during this cellular process. We found that "ribosome" is the first hit of the upregulated KEGG terms ( Figure 3A), indicating that the synthesis of ribosomal proteins was enhanced in hypertrophic cardiomyocytes. In fact, 47 out of 53 significantly altered (FDR < 0.05) ribosomal protein genes were upregulated when we looked at the gene list of Ribo-seq (Figure 3B), indicating an increased ribosomal number and enhanced translational capacity. Increases in the levels of two ribosomal proteins, Rpl35 and Rps23, in PE-treated cardiomyocytes were also confirmed by western blotting ( Figure 3C). Next, we asked whether increased translational rates play a role in PE-induced cardiomyocyte growth. Ribo-seq data were incorporated together with RNAseq data for analysis of gene translational efficiency (TE) using Ribodiff. 24 It turned out that only 109 protein-coding genes were significantly altered in TE (FDR < 0.05; |Log2(fold change)| > 1) (Table S5). Among these, most were genes related to muscle contraction ( Figure 3D); however, most TE-altered genes showed a buffered effect in transcription versus translation (the tendency of dysregulation is opposite between transcript expression and TE) ( Figure 3E). Then, we aligned the changes of translation, transcription, and TE for dysregulated genes in the MAPK signaling pathway, PI3K/AKT signaling pathway, and ribosome. No clear pattern of gene TE in these pathways was observed ( Figures 3F and S2), although dysregulated ribosome genes have a tendency to shift upon TE enhancement. All of these results suggest an increased number of ribosomes but not an increased translational rate, contributed to the enhanced protein synthesis in PE-induced cardiomyocyte hypertrophy.
The translational activity of uORFs does not affect the TE of the main ORF in PE-induced cardiomyocyte hypertrophy The uORF was reported to regulate the translational activity of the main ORF of protein-coding genes. 25 In our study, we identified The proportion of different categories of ORFs detected in our study. Annotated, annotated ORFs for proteins; uORF, upstream ORF; dORF, downstream ORF; overlapped, uORF/dORF overlapped with main ORF; novel_PCGs, novel ORFs from protein-coding genes; novel_-NonPCGs, novel ORFs from nonprotein-coding genes. www.moleculartherapy.org 1,834 uORFs in 1,491 annotated protein-coding genes (Table S6). The length of these uORFs is around 20À50 aa (60À150 nucleotides) (Figure 4A). Among 13,692 detectable ORFs for proteins identified in this study, 12,201 of them (89.11%) did not have an uORF; 1,312 of them (9.58%) possessed one uORF, whereas 179 of them (1.31%) possessed two or more uORFs ( Figure 4B). Those uORF-containing proteincoding genes were further subjected to KEGG analysis. Interestingly, terms of PI3K/AKT signaling pathway and MAPK signaling pathway were shown at the top of the list ( Figure 4C). Next, we tried to address whether uORFs regulate the translational activity of the main ORFs by checking the correlation of TE between the main protein-coding ORF and its preceding uORF (uORF-ORF pair). Multiple uORF-ORF pairs were applied if multiple preceding uORFs for one main protein-coding ORF were identified. No clear distribution pattern was found when we plotted TE ratios of all detected uORF-ORF pairs in Ribo-seq ( Figure 4D). Similar negative results were obtained when we carefully checked uORF-ORF pairs associated with PI3K/AKT signaling pathway or MAPK signaling pathway ( Figure S3). Therefore, the TE of primary ORFs was largely independent from the frequency of uORF translation in PE-induced cardiomyocyte hypertrophy.

Identification of small ORF-coded micropeptides from primary cardiomyocytes
Emerging evidence shows that many micropeptides, which are 100 aa or less, are encoded by sORFs in transcripts originally annotated as noncoding RNAs. In order to identify novel sORFs on a genomewide basis, we carefully examined reads from Ribo-seq, which displayed the 3-nt codon movement pattern, and matched these to noncoding transcripts, which indicated that actively translating ORFs were present in these noncoding transcripts. As a result, 265 novel ORF candidates were detected. Among them, 126 sORF candidates with 100 aa or less were identified (Table S7). These sORF candidates were matched to 103 annotated lncRNAs (8.67%, 103 out of 1,188). We found that 27 sORF candidates overlapped with sORFs detected in adult rat heart previously ( Figure S4). 10 To further evaluate the coding capability, 15 sORF candidates were randomly selected for further testing (Table 1). We built expression constructs containing a 3XFLAG epitope tag in-frame with the C terminus of sORFs. A previously reported sORF, DWORF, 16 was used as a positive control. As expected, a 3XFLAG-tagged micropeptide was detected as ectopically expressed in transfected HEK293 cells. It turned out that 11 out of 15 constructs yielded 3XFLAG-tagged fusion micropeptides ranging from 10 to 20 kDa in HEK293 cells, as demonstrated by western blotting (Table 1; Figure 5A), and the length of 5 0 UTR seemed to have a limited effect on the translation of sORF ( Figure S5). More importantly, the expression of these fusion micropeptides was abolished when the ATG codon was deleted ( Figure 5A) or a single nucleotide was inserted immediately after ATG to cause the frameshift of sORFs in mutant constructs ( Figure S6). Many micropeptides were reported to be localized within mitochrondia. 10 To further explore subcellular localizations of these identified micropeptides, we detected the cellular distribution of these ectopic-expressing 3XFLAGtagged fusion micropeptides in H9C2 cells, a cardiomyocyte-like cell line, with immunostaining. As shown in Figures 5B and S7, most micropeptides showed a predominant nuclear or cytoplasmic localization. Only RNO-sORF7, which is encoded by transcript ENSRNOG00000054516, showed a mitochondrial-specific localization ( Figure 5C). Analyses combined with bioinformatics searches and comparisons to recent data from humans 10 found that most of the validated micropeptides were not conserved in human and mouse, except sORF9 and sORF11 (Table 1; Figure 5D).
To further demonstrate the existence of sORFs in vivo, two antibodies against sORF9 and sORF11 were developed to detect the endogenous expression of these micropeptides, respectively. As shown in Figure 5E, anti-sORF9 antibodies detected endogenous sORF9 micropeptides with molecular weights around 12 kDa and 15 kDa from H9C2 cardiomyocytes in a western blotting analysis. The detected signals were significantly decreased in two mixed H9C2 cell populations in which sORF9 DNA was edited with the CRISPR-Cas9 system using two independent guide RNAs (gRNAs), respectively ( Figure S8). Similarly, anti-sORF11 antibodies detected endogenous sORF11 micropeptides with molecular weights around 17 kDa and 28 kDa from H9C2 cardiomyocytes ( Figure 5F). A small interfering RNA (siRNA) targeting the transcripts of sORF11 was applied to knock down the expression of sORF11. As expected, both the expression of RNA and protein of sORF11 were significantly decreased ( Figure 5F). To further test the conservation of sORF11, we edited sORF11 DNA with the CRISPR-Cas9 system in mouse NIH 3T3 cells. Since the editing induced frameshift mutations in most of the cells (11 out of 12 clones showed a frameshift mutation in a clonal analysis of DNA; Figure S9), a band with a molecular weight around 17 kDa detected in wild-type NIH 3T3 cells was undetected in CRISPR-Cas9-edited cells in a western blotting analysis ( Figure 5G). We noticed that smeared bands larger than predicted size were detected in western blotting with anti-sORF9 and anti-sORF11 antibodies (15 kDa bands for sORF9 and 28 kDa bands for sORF11), which could possibly be caused by protein modifications of these micropeptides.
Interestingly, some transcripts, such as ENSRNOG00000059100 and ENSRNOG00000049537, were found to possess two sORFs  ( Figure 5H). In our analyses, we also noticed that $20% (21 out of 103) of micropeptide-coding lncRNAs were derived from the complementary strand and likely shared a bi-directional promotor driving expression of a protein-coding gene ( Figure 5I), similar to Hand2-UPPERHAND gene pairs. 26,27 Consistent with a previous report, 10 sORFs within UPPERHAND were detected in our Ribo-seq from rat cardiomyocytes but were negative upon further validation (Table  1; ENSRNOG00000052518). However, our successful validation of the coding capability of some sORFs from these lncRNAs ( Figure 5I; Table 1) suggests that partial of these complementary strand-derived lncRNAs could have micropeptide-coding capability.

Micropeptides modulate PE-induced cardiomyocyte hypertrophy through key signaling pathways
The micropeptide DWORF has been reported to benefit the failing heart by regulating Ca 2+ dynamics. 20 However, whether micropeptides participate in the regulation of cardiac hypertrophy has not been studied. To address this question, we first examined the expression of transcripts harboring sORFs validated above in PE-induced cardiomyocyte hypertrophy. The expression of these sORF-harboring transcripts was analyzed from RNA-seq data (Table S8). A more careful examination with quantitative reverse-transcriptase PCR (qRT-PCR) was applied to determine their dysregulation. The result showed that the expression of 5 out of 9 sORF-harbored transcripts was repressed by PE treatment, with transcripts ENSRNOG00000058926, ENSRNOG00000057834, and ENSRNOG00000048986 being significantly upregulated ( Figure 6A). As shown above (Table S5), only $100 ORFs were significantly altered in TE genome-wide during PE treatment. All of these 11 sORFs were found being translated in a stable To further test whether some of these micropeptides modulate the hypertrophic growth of cardiomyocytes, adenovirus was applied to overexpress all validated micropeptides from Figure 5A in primary cardiomyocytes (Figure S10). As a result, 3 micropeptides-RNO-sORF6, RNO-sORF7, and RNO-sORF8-were found to alter hypertrophic growth of cardiomyocytes induced by PE, but not in basal conditions, as determined by the measurement of cellsurface area of a-actinin-positive cardiomyocytes. Cytoplasmic RNO-sORF6 and RNO-sORF8 promoted cardiac hypertrophy. On the contrary, the mitochondrial RNO-sORF7 exerted an inhibitive effect during this process (Figures 6B and 6C). The expression of hypertrophic molecular markers, such as NPPA, NPPB, and ACTA1, further demonstrated the regulatory functions of these micropeptides in cardiac hypertrophy ( Figure 6D), consistent with the size measurement of cardiomyocytes described above. Furthermore, the regulatory function of these micropeptides were abolished when sORFs with ATG start codon deletion were overexpressed ( Figure S11), indicating that micropeptides, but not RNA transcripts, regulate PE-induced cardiomyocyte hypertrophy.
Finally, RNA-seq experiments were performed to explore how these micropeptides drive the shift of the transcriptome during the process of hypertrophy. In RNO-sORF8-overexpressing cardiomyocytes, the KEGG terms of "oxidative phosphorylation" and "citrate cycle (tricarboxylic acid [TCA] cycle)" were enriched in upregulated genes, along with the increase of "hypertrophic cardiomyopathy"-related genes, indicating that the micropeptide RNO-sORF8 induced hypertrophic growth of cardiomyocytes by regulating mitochondrial function (Figure 7A). Interestingly, mitochondrial RNO-sORF7 downregulated genes in the "calcium signaling pathway" and MAPK signaling pathway, which are central in the gene regulatory network for cardiac remodeling and repressed in PE-induced cardiomyocyte hypertrophy ( Figure 7B). In addition, genes related to the MAPK signaling pathway were upregulated, whereas those related to "TNF signaling pathway" and "cell cycle" were downregulated during the promotion of cardiomyocyte hypertrophic growth of micropeptide RNO-sORF6, www.moleculartherapy.org indicating that multiple signaling pathways were coordinated by RNO-sORF6 ( Figure 7C). Dysregulated genes, induced by overexpression of sORFs during hypertrophy in each Gene Ontology (GO) term, were further confirmed by qPCR ( Figures 7DÀ7F). Since the activity of extracellular signal-regulated kinase (ERK) in the MAPK signaling pathway, which was shown to be dysregulated in hypertrophic phenotypes modulated by sORF6 and sORF7, is critical to cardiac hypertrophy, the phosphorylation of ERK in sORF-overexpressing cardiomyocytes was examined by western blotting. Consistently, phosphorylation of ERK was increased in sORF6-overexpressing cardiomyocytes ( Figure 7G) but decreased in sORF7-overexpressing cardiomyocytes during PE-induced hypertrophic growth ( Figure 7H). The phosphorylation level of ERK remained unchanged in sORF8overexpressing cardiomyocytes during PE treatment ( Figure S12). Together, our data suggest that the micropeptides sORF6 and sORF7 modulate cardiac hypertrophy, at least in part, through the ERK signaling pathway.

DISCUSSION
Most biological processes depend on dynamic protein expression. The translation from mRNAs to proteins often undergoes post-transcriptional modifications. Therefore, the dynamics of the transcriptome reflect this process in a less accurate manner. The development of Ribo-seq provided us a chance to study the translatome in a genome-wide and nucleotide-resolution fashion. Combined with RNA-seq and Ribo-seq data, we found that the translational rate from ribosomes was altered for only a small portion of protein-coding genes ($100 genes) during PE-induced cardiomyocyte hypertrophy. Among these TE-altered genes, most of them are contraction-related genes, and their decreased translational rate buffered the increase in mRNA transcription. In addition, PI3K/ AKT and MAPK signaling pathways, which play critical roles in cardiac hypertrophy, were shown transcriptionally controlled with stable translational rates in PE-treated cardiomyocytes. Furthermore, the increase in translation of ribosomal protein genes provided indirect evidence supporting the idea that increased ribosomal numbers (instead of increased translational rates from ribosomes) were responsible for the enhanced protein synthesis during PEinduced cardiomyocyte hypertrophy. More direct evidence with a combination of RNA-seq data and genome-wide proteinomic data detecting newly generated peptides in hypertrophic cardiomyocyte should further prove this point in the future.
Recently, several groups have reported a shift in the translatome during cardiac remodeling. 11,28 In the present study, we provided an integrative analysis using both transcriptome and translatome data of cardiomyocyte hypertrophy induced by PE. In recent studies, the translation landscape has been assessed in the failing human heart. 10 The TE of "extracellular matrix" genes was found to be significantly upregulated, which reflects the activation of cardiac fibroblasts and the manifestation of fibrosis. Combined with the strategy of Ribo-tag, labeling ribosomes in vivo using a genetic approach, another recent report showed that transcripts related to translation, protein quality control, and metabolism are translationally controlled in cardiomyocytes when the heart undergoes pressure overload. 11 Together with these recent studies, the translational landscape of PE-induced cardiomyocyte hypertrophy presented here shows a comprehensive picture of how translation controls cardiac remodeling in different situations.
More and more studies have demonstrated that large numbers of micropeptides encoded by sORFs, a group of previously overlooked molecules, are presented in eukaryotic cells. 29,30 The function of most of these micropeptides has yet to be explored. Interestingly, we found that a portion of micropeptides encoded by transcripts derived from the complementary strand of protein-coding genes might share a bi-directional promotor. Several studies suggest that most of these complementary stranded transcripts are noncoding RNAs. 31,32 Some of these noncoding transcripts, such as UPPER-HAND 26,27 and NEXN-AS1, 33 have been shown to play important roles in cardiovascular development and pathology. The potential micropeptides derived from sORFs in these functional lncRNAs need to be carefully determined. If it is confirmative of micropeptides coding, whether the molecular functions of these transcripts are mediated by the micropeptides or RNA transcripts should be addressed.
Although the functions of some micropeptides have begun to be investigated, 17,20,34,35 the majority of them remained uncharacterized. The current study identified three micropeptides-RNO-sORF6, RNO-sORF7, and RNO-sORF8-that possess a function of modulating cardiac hypertrophy. The regulatory functions of sORF6 and sORF7 were shown to be mediated, at least in part, by the ERK signaling pathway. Our data showed that sORF7 dominantly locates in mitochondria and mediates the downregulation of the calcium signaling pathway and MAPK signaling pathway in PE-induced cardiomyocyte hypertrophy, which is consistent with a previous report that repressing calcium signaling by chelating intracellular Ca 2+ abolished Ang II-induced activation of MAPKs, including ERK. 36 Recently, a mitochondrial micropeptide, mitoregulin, has been shown to regulate intracellular Ca 2+ dynamics. 34 Given that many micropeptides were found to locate in mitochondria in the heart, 10 it is interesting to investigate whether a majority of these mitochondrial micropeptides share a similar signaling cascade in affecting cardiac physiology and pathology. Although detailed mechanisms warrant investigation in detail in the future, our study indicates that micropeptides are a group of new players in regulating cardiac remodeling. Unlike ORFs of conventional protein-coding genes, the majority of which are evolutionarily conserved, a large portion of sORFs identified in this study are not conserved among humans, mice, and rats (except sORF9 and sORF11). Similar observations were also found in the human heart. 10 Although our study of rat-specific sORF-coding micropeptides may not directly contribute to the identification of novel regulators of human cardiac hypertrophy and failure, it suggests that it is promising to search for functional sORF-coding micropeptides in the human heart, as this may reveal novel therapeutic targets.

Cardiomyocyte isolation, culture, and transfection
Neonatal rat ventricular cardiomyocytes (NRVCs) were prepared as previously described. 37 Briefly, 2-day-old neonatal Sprague-Dawley rats were decapitated using sterile scissors, and the chest was opened to extract the heart. NRVCs were isolated by repeated enzymatic dissociation. All isolated cells were pre-plated for 1 h to remove noncardiomyocytes. Nonadherent cells were then plated on 0.5% gelatincoated plates and cultured in DMEM high-glucose medium (Gibco), supplemented with sodium pyruvate (1Â), GlutaMAX (1Â), and 10% fetal bovine serum with antibiotics. Cardiomyocytes were transduced with adenovirus at a multiplicity of infection (MOI) of 20 for 24 h prior to treatment with hypertrophic agent PE (20 mM). Cells were harvested 24 h after PE treatment for RNA isolation, 48 h for detection of ribosomal proteins or immunochemistry, or 1 h for detection of phosphorylated protein. All animal procedures were approved by the Medical Ethics Committee of The First Affiliated Hospital, Sun Yet-sen University.

Ribo-seq
Ribo-seq was performed according to a previous report 38 with minor modifications. A total of 5 Â 10 6 NRVCs were prepared for cell lysis. The cell medium was aspirated, and cells were washed twice with icecold PBS containing 100 mg/mL cycloheximide (MP Biomedicals). After thorough removal of the PBS, the plates were immersed in liquid nitrogen for 10 s and transferred to dry ice. We then dripped   resulting supernatant was used for library preparation. The prepared supernatant was digested with RNase I (Life Technologies; AM2294) for 45 min at room temperature. SUPERase$In RNase Inhibitor (Life Technologies; AM2694) was then added to stop the reaction. Monosomes generated from the digestion were purified with MicroSpin S-400 HR columns (GE Healthcare; 27-5140-01). The RNAs were then isolated from the flowthrough with RNAzol. To enrich the RFPs (ribosome footprints), the RNAs were subjected to the Ribo-Zero Gold rRNA Removal Kit (Illumina; MRZG12324) and then separated with a 17% 7.5 M urea-PAGE gel. RNAs, with sizes ranging from 26 nt to 30 nt, were isolated and purified by PAGE. Strand-specific ribosome-profiling libraries were prepared as described previously. 39 Libraries were sequenced on the HiSeq X Ten platform (Illumina).

Micropeptide expression vector and adenoviral construction
The sORF sequences with an extra 15-nt sequence upstream of the ATG or with predicted full-length 5 0 UTR sequence was cloned from NRVC cDNA and inserted into pCFH vector 3. Frameshift mutant vectors, in which an adenine (A) was inserted behind the ATG translation start codon of the sORF, or ATG translation start codon deletion mutant vectors were constructed according to mutagenesis, a protocol previously described. 40 To yield adenoviral constructs, 3 Â FLAG-tagged sORF cDNAs or sORF cDNAs with deletion of the ATG translation start codon were inserted into the adenoviral empty vector harboring the murine cytomegalovirus (mCMV) promoter. The prepared adenovirus vectors were then transfected in 293AD cells with PEI (1 mg/mL) to produce packed adenovirus. 293AD cells were subsequently transduced with adenovirus for viral amplification. Adenovirus was collected from 293AD cells and purified with the ViraTrap Adenovirus Purification Miniprep Kit, according to the manufacturer's instructions (Biomiga; V1160).

Knockdown of micropeptide expression in cells with siRNA and gRNA
To knock down the expression of micropeptide in H9C2 cardiomyocyte with siRNA, 50 nM targeting siRNA or control siRNA was transfected into cardiomyocytes by using the Lipofectamine RNAiMAX transfection reagent. 6 h later, medium with transfection reagent was removed. Cells were then cultured with full medium for 48 h before harvest for protein extraction. To knock down the expression of micropeptides in cells with CRISPR-Cas 9 system, NIH 3T3 cells and H9C2 cardiomyocytes were transfected with the PX459 vector containing the targeting gRNA to create a frameshift mutation in sORF for 1 day, followed by another 3-day selection with puromycin (3.2 mg/mL). The mixed cell population was then cultured in full medium for protein extraction.

Quantitative real-time PCR and western blotting analysis
Total RNA was isolated from cells, and 2 mg RNA was used for cDNA synthesis. 0.1 mL cDNA was used as a template for quantitative real-time PCR. Primers for qRT-PCR were listed in Table S9. Target gene expression levels were normalized against the expression of glyceraldehyde 3-phosphate dehydrogenase (GAPDH) or 18S rRNA.
The raw data and processed data of RNA-seq and Ribo-seq in this study were deposited in Gene Expression Omnibus (GEO) database of the National Center for Biotechnology Information (NCBI) (series GEO: GSE155914).

Conservation analysis of micropeptide
To analyze the conservation of micropeptide, we first used the BLAT tool (with default parameters) in the University of California, Santa Cruz (UCSC), genome database to search the aa sequences of micropeptides against mouse and human genomes. Hits with over 80% length of an aa sequence matched, and over 80% identical were considered conserved micropeptides. In addition, we also used the BLASTP tool (with default parameters) in the NCBI database to search the aa sequences of micropeptides against nonredundant protein sequence databases. Similarly, hits from genomes other than the rat with over 80% length of aa sequence matched, and over 80% identical were considered conserved micropeptides.

Statistics
Values are reported as mean ± SD unless indicated otherwise. An analysis of variance (ANOVA) followed by Dunnett's post hoc testing was used to evaluate the statistical significance for multiple-group comparisons. In addition, the 2-tailed Mann-Whitney U test was used for 2-group comparisons. Values of p <0.05 were considered statistically significant.