Genomics, Proteomics & Bioinformatics最新文献

Circular RNAs (circRNAs) are endogenous non-coding RNAs with covalently closed structures, which have important functions in plants. However, their biogenesis, degradation, and function upon treatment with gibberellins (GAs) and auxins (1-naphthaleneacetic acid, NAA) remain unknown. Here, we systematically identified and characterized the expression patterns, evolutionary conservation, genomic features, and internal structures of circRNAs using RNase R-treated libraries from moso bamboo (Phyllostachys edulis) seedlings. Moreover, we investigated the biogenesis of circRNAs dependent on both cis- and trans-regulation. We explored the function of circRNAs, including their roles in regulating microRNA (miRNA)-related genes and modulating the alternative splicing of their linear counterparts. Importantly, we developed a customized degradome sequencing approach to detect miRNA-mediated cleavage of circRNAs. Finally, we presented a comprehensive view of the participation of circRNAs in the regulation of hormone metabolism upon treatment of bamboo seedlings with GA and NAA. Collectively, our study provides insights into the biogenesis, function, and miRNA-mediated degradation of circRNAs in moso bamboo.

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.

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).

The epitranscriptomic mark N⁶-methyladenosine (m⁶A), which is the predominant internal modification in RNA, is important for plant responses to diverse stresses. Multiple environmental stresses caused by the tea-withering process can greatly influence the accumulation of specialized metabolites and the formation of tea flavor. However, the effects of the m⁶A-mediated regulatory mechanism on flavor-related metabolic pathways in tea leaves remain relatively uncharacterized. We performed an integrated RNA methylome and transcriptome analysis to explore the m⁶A-mediated regulatory mechanism and its effects on flavonoid and terpenoid metabolism in tea (Camellia sinensis) leaves under solar-withering conditions. Dynamic changes in global m⁶A level in tea leaves were mainly controlled by two m⁶A erasers (CsALKBH4A and CsALKBH4B) during solar-withering treatments. Differentially methylated peak-associated genes following solar-withering treatments with different shading rates were assigned to terpenoid biosynthesis and spliceosome pathways. Further analyses indicated that CsALKBH4-driven RNA demethylation can directly affect the accumulation of volatile terpenoids by mediating the stability and abundance of terpenoid biosynthesis-related transcripts and also indirectly influence the flavonoid, catechin, and theaflavin contents by triggering alternative splicing-mediated regulation. Our findings revealed a novel layer of epitranscriptomic gene regulation in tea flavor-related metabolic pathways and established a link between the m⁶A-mediated regulatory mechanism and the formation of tea flavor under solar-withering conditions.

The N⁶-methylation of RNA adenosines (N⁶-methyladenosine, m⁶A) is an important regulator of gene expression with critical implications in vertebrate and insect development. However, the developmental significance of epitranscriptomes in lophotrochozoan organisms remains unknown. Using methylated RNA immunoprecipitation sequencing (MeRIP-seq), we generated transcriptome-wide m⁶A-RNA methylomes covering the entire development of the oyster from oocytes to juveniles. Oyster RNA classes display specific m⁶A signatures, with messenger RNAs (mRNAs) and long non-coding RNAs (lncRNAs) exhibiting distinct profiles and being highly methylated compared to transposable element (TE) transcripts. Epitranscriptomes are dynamic and correspond to the chronological steps of development (cleavage, gastrulation, organogenesis, and metamorphosis), with minimal mRNA and lncRNA methylation at the morula stage followed by a global increase. mRNA m⁶A levels are correlated with transcript levels, and shifts in methylation profiles correspond to expression kinetics. Differentially methylated transcripts cluster according to embryo-larval stages and bear the corresponding developmental functions (cell division, signal transduction, morphogenesis, and cell differentiation). The m⁶A level of TE transcripts is also regulated and peaks during the gastrulation. We demonstrate that m⁶A-RNA methylomes are dynamic and associated with gene expression regulation during oyster development. The putative epitranscriptome implication in the cleavage, maternal-to-zygotic transition, and cell differentiation in a lophotrochozoan model brings new insights into the control and evolution of developmental processes.

Adenosine-to-inosine (A-to-I) RNA editing, constituting nearly 90% of all RNA editing events in humans, has been reported to contribute to the tumorigenesis in diverse cancers. However, the comprehensive map for functional A-to-I RNA editing events in cancers is still insufficient. To fill this gap, we systematically and intensively analyzed multiple tumorigenic mechanisms of A-to-I RNA editing events in samples across 33 cancer types from The Cancer Genome Atlas. For individual candidate among ∼ 1,500,000 quantified RNA editing events, we performed diverse types of downstream functional annotations. Finally, we identified 24,236 potentially functional A-to-I RNA editing events, including the cases in APOL1, IGFBP3, GRIA2, BLCAP, and miR-589-3p. These events might play crucial roles in the scenarios of tumorigenesis, due to their tumor-related editing frequencies or probable effects on altered expression profiles, protein functions, splicing patterns, and microRNA regulations of tumor genes. Our functional A-to-I RNA editing events (https://ccsm.uth.edu/CAeditome/) will help better understand the cancer pathology from the A-to-I RNA editing aspect.

The genetic information coded in DNA leads to trait innovation via a gene regulatory network (GRN) in development. Here, we developed a conserved non-coding element interpretation method to integrate multi-omics data into gene regulatory network (CNEReg) to investigate the ruminant multi-chambered stomach innovation. We generated paired expression and chromatin accessibility data during rumen and esophagus development in sheep, and revealed 1601 active ruminant-specific conserved non-coding elements (active-RSCNEs). To interpret the function of these active-RSCNEs, we defined toolkit transcription factors (TTFs) and modeled their regulation on rumen-specific genes via batteries of active-RSCNEs during development. Our developmental GRN revealed 18 TTFs and 313 active-RSCNEs regulating 7 rumen functional modules. Notably, 6 TTFs (OTX1, SOX21, HOXC8, SOX2, TP63, and PPARG), as well as 16 active-RSCNEs, functionally distinguished the rumen from the esophagus. Our study provides a systematic approach to understanding how gene regulation evolves and shapes complex traits by putting evo-devo concepts into practice with developmental multi-omics data.

The Chinese crested (CC) duck is a unique indigenous waterfowl breed, which has a crest cushion that affects its survival rate. Therefore, the CC duck is an ideal model to investigate the genetic compensation response to maintain genetic stability. In the present study, we first generated a chromosome-level genome of CC ducks. Comparative genomics revealed that genes related to tissue repair, immune function, and tumors were under strong positive selection, indicating that these adaptive changes might enhance cancer resistance and immune response to maintain the genetic stability of CC ducks. We also assembled a Chinese spot-billed (Csp-b) duck genome, and detected the structural variations (SVs) in the genome assemblies of three ducks (i.e., CC duck, Csp-b duck, and Peking duck). Functional analysis revealed that several SVs were related to the immune system of CC ducks, further strongly suggesting that genetic compensation in the anti-tumor and immune systems supports the survival of CC ducks. Moreover, we confirmed that the CC duck originated from the mallard ducks. Finally, we revealed the physiological and genetic basis of crest traits and identified a causative mutation in TAS2R40 that leads to crest formation. Overall, the findings of this study provide new insights into the role of genetic compensation in adaptive evolution.

Alternative polyadenylation (APA) contributes to transcriptome complexity and gene expression regulation and has been implicated in various cellular processes and diseases. Single-cell RNA sequencing (scRNA-seq) has enabled the profiling of APA at the single-cell level; however, the spatial information of cells is not preserved in scRNA-seq. Alternatively, spatial transcriptomics (ST) technologies provide opportunities to decipher the spatial context of the transcriptomic landscape. Pioneering studies have revealed potential spatially variable genes and/or splice isoforms; however, the pattern of APA usage in spatial contexts remains unappreciated. In this study, we developed a toolkit called stAPAminer for mining spatial patterns of APA from spatially barcoded ST data. APA sites were identified and quantified from the ST data. In particular, an imputation model based on the k-nearest neighbors algorithm was designed to recover APA signals, and then APA genes with spatial patterns of APA usage variation were identified. By analyzing well-established ST data of the mouse olfactory bulb (MOB), we presented a detailed view of spatial APA usage across morphological layers of the MOB. We compiled a comprehensive list of genes with spatial APA dynamics and obtained several major spatial expression patterns that represent spatial APA dynamics in different morphological layers. By extending this analysis to two additional replicates of the MOB ST data, we observed that the spatial APA patterns of several genes were reproducible among replicates. stAPAminer employs the power of ST to explore the transcriptional atlas of spatial APA patterns with spatial resolution. This toolkit is available at https://github.com/BMILAB/stAPAminer and https://ngdc.cncb.ac.cn/biocode/tools/BT007320.

Over the past 20 years, tremendous advances in sequencing technologies and computational algorithms have spurred plant genomic research into a thriving era with hundreds of genomes decoded already, ranging from those of nonvascular plants to those of flowering plants. However, complex plant genome assembly is still challenging and remains difficult to fully resolve with conventional sequencing and assembly methods due to high heterozygosity, highly repetitive sequences, or high ploidy characteristics of complex genomes. Herein, we summarize the challenges of and advances in complex plant genome assembly, including feasible experimental strategies, upgrades to sequencing technology, existing assembly methods, and different phasing algorithms. Moreover, we list actual cases of complex genome projects for readers to refer to and draw upon to solve future problems related to complex genomes. Finally, we expect that the accurate, gapless, telomere-to-telomere, and fully phased assembly of complex plant genomes could soon become routine.