Spermatogenesis is a continual process that occurs in the testes, in which diploid spermatogonial stem cells (SSCs) differentiate and generate haploid spermatozoa. This highly efficient and intricate process is orchestrated at multiple levels. N⁶-methyladenosine (m⁶A), an epigenetic modification prevalent in mRNAs, is implicated in the transcriptional regulation during spermatogenesis. However, the dynamics of m⁶A modification in non-rodent mammalian species remains unclear. Here, we systematically investigated the profile and role of m⁶A during spermatogenesis in pigs. By analyzing the transcriptomic distribution of m⁶A in spermatogonia, spermatocytes, and round spermatids, we identified a globally conserved m⁶A pattern between porcine and murine genes with spermatogenic function. We found that m⁶A was enriched in a group of genes that specifically encode the metabolic enzymes and regulators. In addition, transcriptomes in porcine male germ cells could be subjected to the m⁶A modification. Our data show that m⁶A plays the regulatory roles during spermatogenesis in pigs, which is similar to that in mice. Illustrations of this point are three genes (SETDB1, FOXO1, and FOXO3) that are crucial to the determination of the fate of SSCs. To the best of our knowledge, this study for the first time uncovers the expression profile and role of m⁶A during spermatogenesis in large animals and provides insights into the intricate transcriptional regulation underlying the lifelong male fertility in non-rodent mammalian species.
Animals ; Male ; Mice ; Cell Differentiation/genetics* ; Mammals/metabolism* ; Methylation ; RNA, Messenger/metabolism* ; Spermatogenesis/genetics* ; Spermatozoa/metabolism* ; Swine/genetics* ; Testis/metabolism* ; Transcriptome ; RNA Methylation/genetics*

Serine/arginine-rich splicing factor 7 (SRSF7), a known splicing factor, has been revealed to play oncogenic roles in multiple cancers. However, the mechanisms underlying its oncogenic roles have not been well addressed. Here, based on N⁶-methyladenosine (m⁶A) co-methylation network analysis across diverse cell lines, we find that the gene expression of SRSF7 is positively correlated with glioblastoma (GBM) cell-specific m⁶A methylation. We then indicate that SRSF7 is a novel m⁶A regulator, which specifically facilitates the m⁶A methylation near its binding sites on the mRNAs involved in cell proliferation and migration, through recruiting the methyltransferase complex. Moreover, SRSF7 promotes the proliferation and migration of GBM cells largely dependent on the presence of the m⁶A methyltransferase. The two m⁶A sites on the mRNA for PDZ-binding kinase (PBK) are regulated by SRSF7 and partially mediate the effects of SRSF7 in GBM cells through recognition by insulin-like growth factor 2 mRNA-binding protein 2 (IGF2BP2). Together, our discovery reveals a novel role of SRSF7 in regulating m⁶A and validates the presence and functional importance of temporal- and spatial-specific regulation of m⁶A mediated by RNA-binding proteins (RBPs).
Humans ; Cell Line, Tumor ; Cell Proliferation ; Gene Expression Regulation, Neoplastic ; Glioblastoma/genetics* ; Methyltransferases/metabolism* ; RNA Splicing Factors/metabolism* ; RNA, Messenger/genetics* ; RNA-Binding Proteins/metabolism* ; Serine-Arginine Splicing Factors/metabolism* ; RNA Methylation/genetics*

Stress granules (SGs) are cytoplasmic ribonucleoprotein assemblies formed under stress conditions and are related to various biological processes and human diseases. Previous studies have reported the regulatory role of some proteins and linear RNAs in SG assembly. However, the relationship between circular RNAs (circRNAs) and SGs has not been discovered. Here, we screened both linear RNAs and circRNAs in SGs using improved total RNA sequencing of purified SG cores in mammalian cells and identified circular transcripts specifically localized in SGs. circRNAs with higher SG-related RNA-binding protein (RBP) binding abilities are more likely to be enriched in SGs. Furthermore, some SG-enriched circRNAs are differentially expressed in hepatocellular carcinoma (HCC) and adjacent tissues. These results suggest the regulatory role of circRNAs in SG formation and provide insights into the biological function of circRNAs and SGs in HCC.
Animals ; Humans ; RNA, Circular/metabolism* ; Carcinoma, Hepatocellular/metabolism* ; Stress Granules ; Cytoplasmic Granules/metabolism* ; Liver Neoplasms/metabolism* ; RNA/metabolism* ; Stress, Physiological/genetics* ; Mammals/genetics*

The unicellular green alga Chlamydomonas reinhardtii (hereafter Chlamydomonas) possesses both plant and animal attributes, and it is an ideal model organism for studying fundamental processes such as photosynthesis, sexual reproduction, and life cycle. N⁶-methyladenosine (m⁶A) is the most prevalent mRNA modification, and it plays important roles during sexual reproduction in animals and plants. However, the pattern and function of m⁶A modification during the sexual reproduction of Chlamydomonas remain unknown. Here, we performed transcriptome and methylated RNA immunoprecipitation sequencing (MeRIP-seq) analyses on six samples from different stages during sexual reproduction of the Chlamydomonas life cycle. The results show that m⁶A modification frequently occurs at the main motif of DRAC (D = G/A/U, R = A/G) in Chlamydomonas mRNAs. Moreover, m⁶A peaks in Chlamydomonas mRNAs are mainly enriched in the 3' untranslated regions (3'UTRs) and negatively correlated with the abundance of transcripts at each stage. In particular, there is a significant negative correlation between the expression levels and the m⁶A levels of genes involved in the microtubule-associated pathway, indicating that m⁶A modification influences the sexual reproduction and the life cycle of Chlamydomonas by regulating microtubule-based movement. In summary, our findings are the first to demonstrate the distribution and the functions of m⁶A modification in Chlamydomonas mRNAs and provide new evolutionary insights into m⁶A modification in the process of sexual reproduction in other plant organisms.
Animals ; Chlamydomonas reinhardtii/metabolism* ; Reproduction/genetics* ; Life Cycle Stages/genetics* ; Transcriptome ; Plants/genetics*

After implantation, complex and highly specialized molecular events render functionally distinct organ formation, whereas how the epigenome shapes organ-specific development remains to be fully elucidated. Here, nano-hmC-Seal, RNA bisulfite sequencing (RNA-BisSeq), and RNA sequencing (RNA-Seq) were performed, and the first multilayer landscapes of DNA 5-hydroxymethylcytosine (5hmC) and RNA 5-methylcytosine (m⁵C) epigenomes were obtained in the heart, kidney, liver, and lung of the human foetuses at 13-28 weeks with 123 samples in total. We identified 70,091 and 503 organ- and stage-specific differentially hydroxymethylated regions (DhMRs) and m⁵C-modified mRNAs, respectively. The key transcription factors (TFs), T-box transcription factor 20 (TBX20), paired box 8 (PAX8), krueppel-like factor 1 (KLF1), transcription factor 21 (TCF21), and CCAAT enhancer binding protein beta (CEBPB), specifically contribute to the formation of distinct organs at different stages. Additionally, 5hmC-enriched Alu elements may participate in the regulation of expression of TF-targeted genes. Our integrated studies reveal a putative essential link between DNA modification and RNA methylation, and illustrate the epigenetic maps during human foetal organogenesis, which provide a foundation for for an in-depth understanding of the epigenetic mechanisms underlying early development and birth defects.
Humans ; DNA Methylation ; RNA ; Epigenesis, Genetic ; DNA/genetics* ; Basic Helix-Loop-Helix Transcription Factors/genetics*

N⁶-methyladenine (m⁶A) is the most abundant RNA modification in mammalian messenger RNAs (mRNAs), which participates in and regulates many important biological activities, such as tissue development and stem cell differentiation. Due to an improved understanding of m⁶A, researchers have discovered that the biological function of m⁶A can be linked to many stages of mRNA metabolism and that m⁶A can regulate a variety of complex biological processes. In addition to its location on mammalian mRNAs, m⁶A has been identified on viral transcripts. m⁶A also plays important roles in the life cycle of many viruses and in viral replication in host cells. In this review, we briefly introduce the detection methods of m⁶A, the m⁶A-related proteins, and the functions of m⁶A. We also summarize the effects of m⁶A-related proteins on viral replication and infection. We hope that this review provides researchers with some insights for elucidating the complex mechanisms of the epitranscriptome related to viruses, and provides information for further study of the mechanisms of other modified nucleobases acting on processes such as viral replication. We also anticipate that this review can stimulate collaborative research from different fields, such as chemistry, biology, and medicine, and promote the development of antiviral drugs and vaccines.
Animals ; Viruses/genetics* ; RNA, Messenger/metabolism* ; Adenosine/metabolism* ; Cell Differentiation ; Mammals/metabolism*

Although the function of tRNAs in the translational process is well established, it remains controversial whether tRNA abundance is tightly associated with translational efficiency (TE) in mammals. Moreover, how critically the expression of tRNAs contributes to the establishment of tissue-specific proteomes in mammals has not been well addressed. Here, we measured both tRNA expression using demethylase-tRNA sequencing (DM-tRNA-seq) and TE of mRNAs using ribosome-tagging sequencing (RiboTag-seq) in the brain, heart, and testis of mice. Remarkable variation in the expression of tRNA isodecoders was observed among different tissues. When the statistical effect of isodecoder-grouping on reducing variations is considered through permutating the anticodons, we observed an expected reduction in the variation of anticodon expression across all samples, an unexpected smaller variation of anticodon usage bias, and an unexpected larger variation of tRNA isotype expression at amino acid level. Regardless of whether or not they share the same anticodons, the isodecoders encoding the same amino acids are co-expressed across different tissues. Based on the expression of tRNAs and the TE of mRNAs, we find that the tRNA adaptation index (tAI) and TE are significantly correlated in the same tissues but not between tissues; and tRNA expression and the amino acid composition of translating peptides are positively correlated in the same tissues but not between tissues. We therefore hypothesize that the tissue-specific expression of tRNAs might be due to post-transcriptional mechanisms. This study provides a resource for tRNA and translation studies, as well as novel insights into the dynamics of tRNAs and their roles in translational regulation.
Animals ; Mice ; Anticodon/genetics* ; Transcriptome ; Protein Biosynthesis ; RNA, Transfer/chemistry* ; Amino Acids/metabolism* ; Mammals/metabolism*

As the most pervasive epigenetic marker present on mRNAs and long non-coding RNAs (lncRNAs), N⁶-methyladenosine (m⁶A) RNA methylation has been shown to participate in essential biological processes. Recent studies have revealed the distinct patterns of m⁶A methylome across human tissues, and a major challenge remains in elucidating the tissue-specific presence and circuitry of m⁶A methylation. We present here a comprehensive online platform, m6A-TSHub, for unveiling the context-specific m⁶A methylation and genetic mutations that potentially regulate m⁶A epigenetic mark. m6A-TSHub consists of four core components, including (1) m6A-TSDB, a comprehensive database of 184,554 functionally annotated m⁶A sites derived from 23 human tissues and 499,369 m⁶A sites from 25 tumor conditions, respectively; (2) m6A-TSFinder, a web server for high-accuracy prediction of m⁶A methylation sites within a specific tissue from RNA sequences, which was constructed using multi-instance deep neural networks with gated attention; (3) m6A-TSVar, a web server for assessing the impact of genetic variants on tissue-specific m⁶A RNA modifications; and (4) m6A-CAVar, a database of 587,983 The Cancer Genome Atlas (TCGA) cancer mutations (derived from 27 cancer types) that were predicted to affect m⁶A modifications in the primary tissue of cancers. The database should make a useful resource for studying the m⁶A methylome and the genetic factors of epitranscriptome disturbance in a specific tissue (or cancer type). m6A-TSHub is accessible at www.xjtlu.edu.cn/biologicalsciences/m6ats.
Humans ; Methylation ; RNA, Messenger/metabolism* ; Neoplasms/genetics* ; Mutation

RNA modifications affect many biological processes and physiological diseases. The 5-methylcytosine (m⁵C) modification regulates the progression of multiple tumors. However, its characteristics and functions in hepatocellular carcinoma (HCC) remain largely unknown. Here, we found that HCC tissues had a higher m⁵C methylation level than the adjacent normal tissues. Transcriptome analysis revealed that the hypermethylated genes mainly participated in the phosphokinase signaling pathways, such as the Ras and PI3K-Akt pathways. The m⁵C methyltransferase NSUN2 was highly expressed in HCC tissues. Interestingly, the expression of many genes was positively correlated with the expression of NSUN2, including GRB2, RNF115, AATF, ADAM15, RTN3, and HDGF. Real-time PCR assays further revealed that the expression of the mRNAs of GRB2, RNF115, and AATF decreased significantly with the down-regulation of NSUN2 expression in HCC cells. Furthermore, NSUN2 could regulate the cellular sensitivity of HCC cells to sorafenib via modulating the Ras signaling pathway. Moreover, knocking down NSUN2 caused cell cycle arrest. Taken together, our study demonstrates the vital role of NSUN2 in the progression of HCC.
Humans ; ADAM Proteins ; Apoptosis Regulatory Proteins ; Carcinoma, Hepatocellular/metabolism* ; Liver Neoplasms/metabolism* ; Membrane Proteins ; Methyltransferases/metabolism* ; Phosphatidylinositol 3-Kinases/metabolism* ; Repressor Proteins ; RNA, Messenger/metabolism* ; Sorafenib ; RNA Methylation/genetics*

Epidermal growth factor receptor-tyrosine kinase inhibitors (EGFR-TKIs) positively affect the initial control of non-small cell lung cancer (NSCLC). Rapidly acquired resistance to EGFR-TKIs is a major hurdle in successful treatment. However, the mechanisms that control the resistance of EGFR-TKIs remain largely unknown. RNA structures have widespread and crucial functions in many biological regulations; however, the functions of RNA structures in regulating cancer drug resistance remain unclear. Here, the psoralen analysis of RNA interactions and structures (PARIS) method is used to establish the higher-order RNA structure maps of EGFR-TKIs-resistant and -sensitive cells of NSCLC. Our results show that RNA structural regions are enriched in untranslated regions (UTRs) and correlate with translation efficiency (TE). Moreover, yrdC N⁶-threonylcarbamoyltransferase domain containing (YRDC) promotes resistance to EGFR-TKIs. RNA structure formation in YRDC 3' UTR suppresses embryonic lethal abnormal vision-like 1 (ELAVL1) binding, leading to EGFR-TKI sensitivity by impairing YRDC translation. A potential therapeutic strategy for cancer treatment is provided using antisense oligonucleotide (ASO) to perturb the interaction between RNA and protein. Our study reveals an unprecedented mechanism through which the RNA structure switch modulates EGFR-TKI resistance by controlling YRDC mRNA translation in an ELAVL1-dependent manner.
Humans ; Carcinoma, Non-Small-Cell Lung/genetics* ; Cell Line, Tumor ; Drug Resistance, Neoplasm/genetics* ; ErbB Receptors/metabolism* ; GTP-Binding Proteins/therapeutic use* ; Lung Neoplasms/genetics* ; Mutation ; Protein Kinase Inhibitors/therapeutic use* ; RNA ; RNA-Binding Proteins/genetics* ; Tyrosine Kinase Inhibitors/therapeutic use*