Transcriptome-wide analysis of N6-methyladenosine uncovers its regulatory role in gene expression in the lepidopteran Bombyx mori

Targeted DNA N6-methyladenine editing by dCas9 fused to METTL4 in the lepidopteran model insect Bombyx mori

We have established a novel CRISPR-dCas9-METTL4 epigenome editing tool that can methylate target regions to achieve site-specific DNA 6mA methylation in both hypermethylated and hypomethylated genes. Targeted methylation on genes by dCas9-METTL4 results in misexpression, allowing for the functional investigation of target genes of interest in silkworm.

Juvenile hormone regulates silk gene expression by m6A RNA methylation

Juvenile hormone (JH) is an indispensable insect hormone that is critical in regulating insect development and physiology. N6-methyladenosine (m6A) is the most abundant modification of RNA that regulates RNA fate in eukaryotic organisms. However, the relationship between m6A and JH remains largely unknown. Here, we found that the application of a Juvenile hormone analog (JHA) extended the larval period of Bombyx mori and increased the weight and thickness of the cocoon. Interestingly, global transcriptional patterns revealed that m6A-related genes are specifically regulated by JHA in the posterior silk gland (PSG) that synthesizes the major component of cocoon silk. By transcriptome and m6A sequencing data conjointly, we discovered that JHA significantly regulated the m6A modification in the PSG of B. mori and many m6A-containing genes are related to nucleic acid binding, nucleus, and nucleobase-containing compound metabolism. Notably, 547 genes were significantly regulated by JHA at both the m6A modification and expression levels, especially 16 silk-associated genes, including sericin2, seroin1, Serine protease inhibitors 4 (BmSPI4), Serine protease inhibitors 5 (BmSPI5), and LIM domain-binding protein 2 (Ldb). Among them, 11 silk associated genes were significantly affected by METTL3 knockdown, validating that these genes are targets of m6A modification. Furthermore, we confirm that JHA directly regulates the expression of BmSPI4 and BmSPI5 through m6A modification of CDS regions. These results demonstrate the essential role of m6A methylation regulated by JH in PSG, and elucidate a novel mechanism by which JH affects silk gland development via m6A methylation. This study uncovers that m6A modification is a critical factor mediating the effect of JH in insects.

Transcriptome-wide analysis of mRNA N6 -methyladenosine modification in the embryonic development of Spodoptera frugiperda

N6 -methyladenosine (m6 A) RNA is the most abundant modification of mRNA, and has been demonstrated in regulating various post-transcriptional processes. Many studies have shown that m6 A methylation plays key roles in sex determination, neuronal functions, and embryonic development in Drosophila and mammals. Here, we analyzed transcriptome-wide profile of m6 A modification in the embryonic development of the destructive agricultural pest Spodoptera frugiperda. We found that the 2 key mRNA m6 A methyltransferases SfrMETTL3 and SfrMETTL14 have high homologies with other insects and mammals, suggesting that SfrMETTL3 and SfrMETTL14 may have conserved function among different species. From methylated RNA immunoprecipitation sequencing analysis, we obtained 46 869 m6 A peaks representing 8 587 transcripts in the 2-h embryos after oviposition, and 41 389 m6 A peaks representing 9 230 transcripts in the 24-h embryos. In addition, 5 995 m6 A peaks were differentially expressed including 3 752 upregulated and 2243 downregulated peaks. Functional analysis with Gene Ontology and Kyoto Encyclopedia of Genes and Genomes suggested that differentially expressed m6 A peak-modified genes were enriched in cell and organ development between the 2- and 24-h embryos. By conjoint analysis of methylated RNA immunoprecipitation-seq and RNA-seq data, we found that RNA m6 A methylation may regulate the transcriptional levels of genes related to tissue and organ development from 2- to 24-h embryos. Our study reveals the role of RNA m6 A epigenetic regulation in the embryonic development of S. frugiperda, and provides new insights for the embryonic development of insects.

Enthralling genetic regulatory mechanisms meddling insecticide resistance development in insects: role of transcriptional and post-transcriptional events

Insecticide resistance in insects severely threatens both human health and agriculture, making insecticides less compelling and valuable, leading to frequent pest management failures, rising input costs, lowering crop yields, and disastrous public health. Insecticide resistance results from multiple factors, mainly indiscriminate insecticide usage and mounted selection pressure on insect populations. Insects respond to insecticide stress at the cellular level by modest yet significant genetic propagations. Transcriptional, co-transcriptional, and post-transcriptional regulatory signals of cells in organisms regulate the intricate processes in gene expressions churning the genetic information in transcriptional units into proteins and non-coding transcripts. Upregulation of detoxification enzymes, notably cytochrome P450s (CYPs), glutathione S-transferases (GSTs), esterases [carboxyl choline esterase (CCE), carboxyl esterase (CarE)] and ATP Binding Cassettes (ABC) at the transcriptional level, modification of target sites, decreased penetration, or higher excretion of insecticides are the noted insect physiological responses. The transcriptional regulatory pathways such as AhR/ARNT, Nuclear receptors, CncC/Keap1, MAPK/CREB, and GPCR/cAMP/PKA were found to regulate the detoxification genes at the transcriptional level. Post-transcriptional changes of non-coding RNAs (ncRNAs) such as microRNAs (miRNA), long non-coding RNAs (lncRNA), and epitranscriptomics, including RNA methylation, are reported in resistant insects. Additionally, genetic modifications such as mutations in the target sites and copy number variations (CNV) are also influencing insecticide resistance. Therefore, these cellular intricacies may decrease insecticide sensitivity, altering the concentrations or activities of proteins involved in insecticide interactions or detoxification. The cellular episodes at the transcriptional and post-transcriptional levels pertinent to insecticide resistance responses in insects are extensively covered in this review. An overview of molecular mechanisms underlying these biological rhythms allows for developing alternative pest control methods to focus on insect vulnerabilities, employing reverse genetics approaches like RNA interference (RNAi) technology to silence particular resistance-related genes for sustained insect management.

m6A Profile Dynamics Indicates Regulation of Oyster Development by m6A-RNA Epitranscriptomes

The N6-methylation of RNA adenosines (N6-methyladenosine, m6A) 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 m6A-RNA methylomes covering the entire development of the oyster from oocytes to juveniles. Oyster RNA classes display specific m6A 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 m6A 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 m6A level of TE transcripts is also regulated and peaks during the gastrulation. We demonstrate that m6A-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.

Dynamic N6-methyladenosine RNA methylation landscapes reveal epi-transcriptomic modulation induced by ammonia nitrogen exposure in the Pacific whiteleg shrimp Litopenaeus vannamei

Despite the versatility of RNA m⁶A methylation in regulating various biological processes, its involvement in the physiological response to ammonia nitrogen toxicity in decapod crustaceans like shrimp remains enigmatic. Here, we provided the first characterization of dynamic RNA m⁶A methylation landscapes induced by toxic ammonia exposure in the Pacific whiteleg shrimp Litopenaeus vannamei. The global m⁶A methylation level showed significant decrease following ammonia exposure, and most of the m⁶A methyltransferases and m⁶A binding proteins were significantly repressed. Distinct from many well-studied model organisms, m⁶A methylated peaks in the transcriptome of L. vannamei were enriched not only near the termination codon and in the 3' untranslated region (UTR), but also around the start codon and in the 5' UTR. Upon ammonia exposure, 11,430 m⁶A peaks corresponding to 6113 genes were hypo-methylated, and 5660 m⁶A peaks from 3912 genes were hyper-methylated. The differentially methylated genes showing significant changes in expression were over-represented by genes associated with metabolism, cellular immune defense and apoptotic signaling pathways. Notably, the m⁶A-modified ammonia-responsive genes encompassed a subset of genes related to glutamine synthesis, purine conversion and urea production, implying that m⁶A methylation may modulate shrimp ammonia stress responses partly through these ammonia metabolic processes.

RNA methyltransferase BmMettl3 and BmMettl14 in silkworm (Bombyx mori) and the regulation of silkworm embryonic development

N6-methyladenosine (m6A) is a ubiquitous reversible epigenetic RNA modification that plays an important role in regulating many biological processes, especially embryonic development. However, regulation of m6A methylation during silkworm embryonic development and diapause remains to be investigated. In this study, we analyzed the phylogeny of subunits of methyltransferases BmMettl3 and BmMettl14, and detected the expression patterns of BmMettl3 and BmMettl14 in different tissues and at different developmental stages in silkworm. To investigate the function of m6A on the development of silkworm embryo, we analyzed the m6A/A ratio in diapause and diapause termination eggs. The results showed that BmMettl3 and BmMettl14 were highly expressed in gonads and eggs. Moreover, the expression of BmMettl3 and BmMettl14 and the m6A/A ratio were significantly increased in diapause termination eggs compared with diapause eggs in the early stage of silkworm embryonic development. Furthermore, in BmN cell cycle experiments, the percentage of cells in the S phase increased when lacking BmMettl3 or BmMettl14. This work contributes to understanding the role of m6A methylation during insect embryogenesis and gametogenesis. It also provides a research orientation to further analyze the role of m6A methylation in diapause initiation and termination during insect embryonic development.

Evolution of m6A-related genes in insects and the function of METTL3 in silkworm embryonic development

N6-methyladenosine (m6A) plays a key role in many biological processes. However, the function and evolutionary relationship of m6A-related genes in insects remain largely unknown. Here we analysed the phylogeny of m6A-related genes among 207 insect species and found that m6A-related genes are evolutionarily conserved in insects. Subcellular localization experiments of m6A-related proteins in BmN cells confirmed that BmYTHDF3 was localized in the cytoplasm, BmMETTL3, BmMETTL14, and BmYTHDC were localized in the nucleus, and FL2D was localized to both the nucleus and cytoplasm. We examined the expression patterns of m6A-related genes during the embryonic development of Bombyx mori. To elucidate the function of BmMETTL3 during the embryonic stage, RNA sequencing was performed to measure changes in gene expression in silkworm eggs after BmMETTL3 knockdown, as well as in BmN cells overexpressing BmMETTL3. The global transcriptional pattern showed that knockdown of BmMETTL3 affected multiple cellular processes, including oxidoreductase activity, transcription regulator activity, and the cation binding. In addition, transcriptomic data revealed that many observed DEGs were associated with fundamental metabolic processes, including carbon metabolism, purine metabolism, amino acid biosynthesis, and the citrate cycle. Interestingly, we found that knockdown of BmMETTL3 significantly affected Wnt and Toll/Imd pathways in embryos. Taken together, these results suggest that BmMETTL3 plays an essential role in the embryonic development of B. mori, and deepen our understanding of the function of m6A-related genes in insects.

Potential role of N6-adenine DNA methylation in alternative splicing and endosymbiosis in Paramecium bursaria

N⁶-adenine DNA methylation (6mA), a rediscovered epigenetic mark in eukaryotic organisms, diversifies in abundance, distribution, and function across species, necessitating its study in more taxa. Paramecium bursaria is a typical model organism with endosymbiotic algae of the species Chlorella variabilis. This consortium therefore serves as a valuable system to investigate the functional role of 6mA in endosymbiosis, as well as the evolutionary importance of 6mA among eukaryotes. In this study, we report the first genome-wide, base pair-resolution map of 6mA in P. bursaria and identify its methyltransferase PbAMT1. Functionally, 6mA exhibits a bimodal distribution at the 5' end of RNA polymerase II-transcribed genes and possibly participates in transcription by facilitating alternative splicing. Evolutionarily, 6mA co-evolves with gene age and likely serves as a reverse mark of endosymbiosis-related genes. Our results offer new insights for the functional diversification of 6mA in eukaryotes as an important epigenetic mark.

Epitranscriptomic mRNA modifications governing plant stress responses: underlying mechanism and potential application

Plants inevitably encounter environmental adversities, including abiotic and biotic stresses, which significantly impede plant growth and reduce crop yield. Thus, fine-tuning the fate and function of stress-responsive RNAs is indispensable for plant survival under such adverse conditions. Recently, post-transcriptional RNA modifications have been studied as a potent route to regulate plant gene expression under stress. Among over 160 mRNA modifications identified to date, N⁶ -methyladenosine (m⁶ A) in mRNAs is notable because of its multifaceted roles in plant development and stress response. Recent transcriptome-wide mapping has revealed the distribution and patterns of m⁶ A in diverse stress-responsive mRNAs in plants, building a foundation for elucidating the molecular link between m⁶ A and stress response. Moreover, the identification and characterization of m⁶ A writers, readers and erasers in Arabidopsis and other model crops have offered insights into the biological roles of m⁶ A in plant abiotic stress responses. Here, we review the recent progress of research on mRNA modifications, particularly m⁶ A, and their dynamics, distribution, regulation and biological functions in plant stress responses. Further, we posit potential strategies for breeding stress-tolerant crops by engineering mRNA modifications and propose the future direction of research on RNA modifications to gain a much deeper understanding of plant stress biology.

Identification and functional analysis of m6A in the mammary gland tissues of dairy goats at the early and peak lactation stages

N6-methyladenosine (m⁶A) is the most common reversible epigenetic RNA modification in the mRNA of all higher eukaryotic organisms and plays an important role in the regulation of gene expression and cell function. In this study, m⁶A-modified methylated RNA immunoprecipitation sequencing (MeRIP-seq) and transcriptome sequencing (RNA-seq) were used to identify the key genes with m⁶A modification during mammary gland development and lactation in dairy goats. The results showed that m⁶A methylation occurred at 3,927 loci, which were significantly enriched in the 3′ untranslated region (3′UTR) and the termination codon region. In the early stage and peak stage of lactation, m⁶A methylation occurred extensively in mammary tissues, and a total of 725 differentially expressed m⁶A-modified genes were obtained, all negatively correlated with mRNA expression. In addition, Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that different methylated genes were mainly involved in the growth and apoptosis of mammary epithelial cells through signaling pathways, such as the mitogen-activated protein kinase (MAPK) and phospholipase D pathways, and then affected the development and lactation of mammary gland. All in all, we identified and analyzed the methylation events related to the development and lactation regulation of mammary gland at the early and peak lactation stages, and provided a theoretical basis to reveal the physiological regulatory system of mammary gland development and lactation in dairy goats.

Transcriptional expression of m6A and m5C RNA methyltransferase genes in the brain and fat body of honey bee adult workers

Honey bee (Apis mellifera) adult workers change behaviors and nutrition according to age progression. Young workers, such as nurses, perform in-hive tasks and consume protein-rich pollen, while older workers (foragers) leave the colony to search for food, and consume carbohydrate-rich nectar. These environmentally stimulated events involve transcriptional and DNA epigenetic marks alterations in worker tissues. However, post-transcriptional RNA modifications (epitranscriptomics) are still poorly explored in bees. We investigated the transcriptional profiles of m6A and m5C RNA methyltransferases in the brain and fat body of adult workers of 1) different ages and performing different tasks [nurses of 8 days-old (N-8D) and foragers of 29 days-old (F-29D), sampled from wild-type colonies], and 2) same-aged young workers caged in an incubator and treated with a pollen-rich [PR] or a pollen-deprived [PD] diet for 8 days. In the brain, METTL3, DNMT2, NOP2, NSUN2, NSUN5, and NSUN7 genes increased expression during adulthood (from N-8D to F-29D), while the opposite pattern was observed in the fat body for METTL3, DNMT2, and NSUN2 genes. Regarding diet treatments, high expression levels were observed in the brains of the pollen-deprived group (DNMT2, NOP2, and NSUN2 genes) and the fat bodies of the pollen-rich group (NOP2, NSUN4, and NSUN5 genes) compared to the brains of the PR group and the fat bodies of the PD group, respectively. Our data indicate that RNA epigenetics may be an important regulatory layer in the development of adult workers, presenting tissue-specific signatures of RNA methyltransferases expression in response to age, behavior, and diet content.

Attenuated Duck Hepatitis A Virus Infection Is Associated With High mRNA Maintenance in Duckling Liver via m6A Modification

Host translation is generally modulated by viral infection, including duck hepatitis A virus (DHAV) infection. Previously, we reported that cellular protein synthesis in a cell model of duck embryo fibroblasts is significantly inhibited by DHAV infection but not viral proteins, suggesting that an important viral-host interaction occurs at the translational level. In this study, we aim to further understand the impact of DHAV virulence on cellular N6-methyladenosine (m6A) modification, which is essential to a wide variety of RNA biological processes, such as mRNA stabilization and translation. Using m6A antibody-based immunoprecipitation, m6A-seq, and LC–MS/MS, we observed that m6A-modified mRNA exists in both virulent and attenuated DHAV-infected duckling livers. Importantly, m6A levels in mRNA were much higher in attenuated DHAV-infected livers compared with virulent DHAV-infected livers, suggesting virulence-dependent regulation of m6A modification. Analysis of modification motifs indicated that GAAGAAG is the most enriched motif. Combined m6A-seq and RNA-seq data analysis indicated a generally positive correlation between m6A and mRNA expression levels in DHAV-infected duckling livers. GO analysis of genes with decreased or increased m6A levels showed that these genes were enriched in various terms, including oxidation–reduction processes and antiviral immune responses. Collectively, our work reveals DHAV virulence-dependent coordination between m6A modification and mRNA expression in duckling livers.

N6-methyladenosine modification of the Aedes aegypti transcriptome and its alteration upon dengue virus infection in Aag2 cell line

The N⁶-methyladenosine (m⁶A) modification of RNA has been reported to affect viral infections. Studies have confirmed the role of m⁶A in replication of several vector-borne flaviviruses, including dengue virus (DENV), in mammalian cells. Here, we explored the role of m⁶A in DENV replication in the mosquito Aedes aegypti Aag2 cell line. We first determined the presence of m⁶A on the RNAs from mosquito cells and using methylated RNA immunoprecipitation and sequencing (MeRIP-Seq) identified m⁶A modification of the mosquito transcriptome and those that changed upon DENV infection. Depletion of m⁶A methyltransferases and the m⁶A binding protein YTHDF3 RNAs decreased the replication of DENV. In particular, we found that the Ae. aegypti ubiquitin carrier protein 9 (Ubc9) is m⁶A modified and its expression increases after DENV infection. Silencing of the gene and ectopic expression of Ubc9 led to reduced and increased DENV replication, respectively. The abundance of Ubc9 mRNA and its stability were reduced with the inhibition of m⁶A modification, implying that m⁶A modification of Ubc9 might enhance expression of the gene. We also show that the genome of DENV is m⁶A modified at five sites in mosquito cells. Altogether, this work reveals the involvement of m⁶A modification in Ae. aegypti-DENV interaction.

Analysis of m⁶A RNA modifications in the mosquito transcriptome and their changes upon dengue virus infection provides insight into the role of epigenetics in regulating viral replication in mosquitoes.

The m6A methylation regulates gonadal sex differentiation in chicken embryo

Background

As a ubiquitous reversible epigenetic RNA modification, N6-methyladenosine (m6A) plays crucial regulatory roles in multiple biological pathways. However, its functional mechanisms in sex determination and differentiation during gonadal development of chicken embryos are not clear. Therefore, we established a transcriptome-wide m6A map in the female and male chicken left gonads of embryonic day 7 (E7) by methylated RNA immunoprecipitation sequencing (MeRIP-seq) to offer insight into the landscape of m6A methylation and investigate the post-transcriptional modification underlying gonadal differentiation.

Results

The chicken embryonic gonadal transcriptome was extensively methylated. We found 15,191 and 16,111 m6A peaks in the female and male left gonads, respectively, which were mainly enriched in the coding sequence (CDS) and stop codon. Among these m6A peaks, we identified that 1013 and 751 were hypermethylated in females and males, respectively. These differential peaks covered 281 and 327 genes, such as BMP2, SMAD2, SOX9 and CYP19A1, which were primarily associated with development, morphogenesis and sex differentiation by functional enrichment. Further analysis revealed that the m6A methylation level was positively correlated with gene expression abundance. Furthermore, we found that YTHDC2 could regulate the expression of sex-related genes, especially HEMGN and SOX9, in male mesonephros/gonad mingle cells, which was verified by in vitro experiments, suggesting a regulatory role of m6A methylation in chicken gonad differentiation.

Conclusions

This work provided a comprehensive m6A methylation profile of chicken embryonic gonads and revealed YTHDC2 as a key regulator responsible for sex differentiation. Our results contribute to a better understanding of epigenetic factors involved in chicken sex determination and differentiation and to promoting the future development of sex manipulation in poultry industry.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-022-00710-6.

Profiling of Transcriptome-Wide N6-Methyladenosine (m6A) Modifications and Identifying m6A Associated Regulation in Sperm Tail Formation in Anopheles sinensis

Recent discoveries of reversible N6-methyladenosine (m6A) methylation on messenger RNA (mRNA) and mapping of m6A methylomes in many species have revealed potential regulatory functions of this RNA modification by m6A players—writers, readers, and erasers. Here, we first profile transcriptome-wide m6A in female and male Anopheles sinensis and reveal that m6A is also a highly conserved modification of mRNA in mosquitoes. Distinct from mammals and yeast but similar to Arabidopsis thaliana, m6A in An. sinensis is enriched not only around the stop codon and within 3′-untranslated regions but also around the start codon and 5′-UTR. Gene ontology analysis indicates the unique distribution pattern of m6A in An. sinensis is associated with mosquito sex-specific pathways such as tRNA wobble uridine modification and phospholipid-binding in females, and peptidoglycan catabolic process, exosome and signal recognition particle, endoplasmic reticulum targeting, and RNA helicase activity in males. The positive correlation between m6A deposition and mRNA abundance indicates that m6A can play a role in regulating gene expression in mosquitoes. Furthermore, many spermatogenesis-associated genes, especially those related to mature sperm flagellum formation, are positively modulated by m6A methylation. A transcriptional regulatory network of m6A in An. sinensis is first profiled in the present study, especially in spermatogenesis, which may provide a new clue for the control of this disease-transmitting vector.

Methyltransferase-like (METTL)14-mediated N6-methyladenosine modification modulates retinal pigment epithelial (RPE) activity by regulating the methylation of microtubule-associated protein (MAP)2

The expression of METTL14 is significantly reduced in patients with retinitis pigmentosa (RP). To clarify the significance of the N6-methyladenosine (m⁶A) RNA modification in RP, we examined phagocytosis, apoptosis, and cell cycle distribution in a human RPE cell line, ARPE-19, following lentivirus-mediated knockdown of METTL14. Differentially expressed genes and changes in m⁶A level were evaluated by RNA sequencing (RNA-seq) and methylated RNA immunoprecipitation sequencing (MeRIP-seq), respectively. The results showed that phagocytosis and proliferation were decreased whereas apoptosis was increased in RPE cells by METTL14 silencing. We found that METTL14 directly regulated m⁶A level and the expression of MAP2, as determined by RNA-seq, MeRIP-seq, MeRIP quantitative PCR, and the RNA pull-down assay. Additionally, MAP2 could bind to neuronal differentiation (NEUROD)1, a pathogenic gene in RPE-associated diseases. A family member of the YTH domain, (YTHDF)2 was recognized as an m⁶A reader of MAP2 mRNA. MAP2 overexpression had the same effects as METTL14 knockdown in RPE cells. Thus, METTL14 regulates the expression of MAP2 via the modification of m⁶A, resulting in the dysregulation of NEUROD1 and pathologic changes in RPE cells. These findings suggest that therapeutic strategies targeting the m⁶A modification of MAP2 or the METTL14/YTHDF2/MAP2/NEUROD1 signaling axis may be effective in the treatment of RPE-associated ocular diseases.

Identification of an m6A RNA Methylation Regulator Risk Score Model for Prediction of Clinical Prognosis in Astrocytoma

Astrocytoma (AS) is the most ubiquitous primary malignancy of the central nervous system (CNS). The vital involvement of the N6-methyladenosine (m6A) RNA modification in the growth of multiple human tumors is known. This study entailed probing m6A regulators with AS prognosis to construct a risk prediction model (RS) for potential clinical use. A total of 579 AS patients' (of the Chinese Glioma Genome Atlas,CGGA) data and the expression of 12 published m6A-related genes were included in this study. Cox and selection operator (LASSO) regression analyses for independent prognostic factors and multifactor Cox analysis established an R.S. model to predict the AS patient prognosis. This was subject to verification employing 331 samples from the TCGA data set followed by gene ontology and pathway enrichment study with gene set enrichment analysis (GSEA). The R.S. constructed with three m6A genes inclusive of WTAP, RBM15, and YTHDF2 emerged as independent prognostic factors in AS patients with vital involvement in the advancement and development of the malignancy. In a nutshell, this work reported an m6A-related gene risk model to predict the prognosis of AS patients to pave the way for discerning diagnostic and prognostic biomarkers. Further corroboration employing relevant wet-lab assays of this model is warranted.

Active genic machinery for epigenetic RNA modifications in bees

Epitranscriptomics is an emerging field of investigation dedicated to the study of post-transcriptional RNA modifications. RNA methylations regulate RNA metabolism and processing, including changes in response to environmental cues. Although RNA modifications are conserved from bacteria to eukaryotes, there is little evidence of an epitranscriptomic pathway in insects. Here we identified genes related to RNA m⁶ A (N6-methyladenine) and m⁵ C (5-methylcytosine) methylation machinery in seven bee genomes (Apis mellifera, Melipona quadrifasciata, Frieseomelitta varia, Eufriesea mexicana, Bombus terrestris, Megachile rotundata and Dufourea novaeangliae). In A. mellifera, we validated the expression of methyltransferase genes and found that the global levels of m⁶ A and m⁵ C measured in the fat body and brain of adult workers differ significantly. Also, m⁶ A levels were differed significantly mainly between the fourth larval instar of queens and workers. Moreover, we found a conserved m⁵ C site in the honeybee 28S rRNA. Taken together, we confirm the existence of epitranscriptomic machinery acting in bees and open avenues for future investigations on RNA epigenetics in a wide spectrum of hymenopteran species.

The dynamics of N6-methyladenine RNA modification in interactions between rice and plant viruses

BACKGROUND

N6-methyladenosine (m6A) is the most common RNA modification in eukaryotes and has been implicated as a novel epigenetic marker that is involved in various biological processes. The pattern and functional dissection of m6A in the regulation of several major human viral diseases have already been reported. However, the patterns and functions of m6A distribution in plant disease bursting remain largely unknown.

RESULTS

We analyse the high-quality m6A methylomes in rice plants infected with two devastating viruses. We find that the m6A methylation is mainly associated with genes that are not actively expressed in virus-infected rice plants. We also detect different m6A peak distributions on the same gene, which may contribute to different antiviral modes between rice stripe virus or rice black-stripe dwarf virus infection. Interestingly, we observe increased levels of m6A methylation in rice plant response to virus infection. Several antiviral pathway-related genes, such as RNA silencing-, resistance-, and fundamental antiviral phytohormone metabolic-related genes, are also m6A methylated. The level of m6A methylation is tightly associated with its relative expression levels.

CONCLUSIONS

We revealed the dynamics of m6A modification during the interaction between rice and viruses, which may act as a main regulatory strategy in gene expression. Our investigations highlight the significance of m6A modifications in interactions between plant and viruses, especially in regulating the expression of genes involved in key pathways.

Epitranscriptomic regulation of insecticide resistance

N 6-methyladenosine (m6A) is the most prevalent messenger RNA modification in eukaryotes and an important posttranscriptional regulator of gene expression. However, the biological roles of m6A in most insects remain largely unknown. Here, we show that m6A regulates a cytochrome P450 gene (CYP4C64) in the global whitefly pest, Bemisia tabaci, leading to insecticide resistance. Investigation of the regulation of CYP4C64, which confers resistance to the insecticide thiamethoxam, revealed a mutation in the 5' untranslated region of this gene in resistant B. tabaci strains that introduces a predicted m6A site. We provide several lines of evidence that mRNA methylation of the adenine at this position, in combination with modified expression of m6A writers, acts to increase expression of CYP4C64 and resistance. Collectively, these results provide an example of the epitranscriptomic regulation of the xenobiotic response in insects and implicate the m6A regulatory axis in the development of insecticide resistance.

Identification and functional annotation of m6A methylation modification in granulosa cells during antral follicle development in pigs

The N⁶-methyladenosine (m6A) derivative has the capacity for ubiquitous epigenetic modification of messenger RNA (mRNA) that regulates gene expression through post-transcriptional mRNA modifications. Findings with mapping of m6A methylomes have indicated there are potential functions of this derivative in different cell types of several species. A profile of m6A methylomes and potential functions in granulosa cells of pigs during antral follicle development, however, has not yet occurred. In the present study, there was profiling of an epitranscriptome-wide map of m6A methylation in granulosa cells of pigs derived from small and large follicles using methylated RNA immunoprecipitation techniques, next-generation sequencing and further annotation of the potential functions of m6A utilizing bioinformatic analyses procedures. The m6A modification is abundant in granulosa cells of pigs, and there are dynamic changes in m6A methylomes during the developmental transition from small (< 3 mm) to large (> 5 mm) sized follicles. In particular, there was a prevalence of 7289 and 6882 m6A in granulosa cells from follicles of two different sizes. There was an increased prevalence of m6A in close proximity to the 5' or 3'-untranslated coding regions and a shared conserved consensus motif. Results from further analysis indicated there was significant enrichment of differentially expressed m6A methylated genes in several signaling pathways associated with steroidogenesis, granulosa cell proliferation and follicular development. When considered as a whole, these results indicate there are differential m6A modifications in granulosa cells of pigs during follicle development that are potentially associated with steroidogenesis and folliculogenesis.

SGID: a comprehensive and interactive database of the silkworm

Although the domestic silkworm (Bombyx mori) is an important model and economic animal, there is a lack of comprehensive database for this organism. Here, we developed the silkworm genome informatics database (SGID). It aims to bring together all silkworm-related biological data and provide an interactive platform for gene inquiry and analysis. The function annotation in SGID is thorough and covers 98% of the silkworm genes. The annotation details include function description, Gene Ontology, Kyoto Encyclopedia of Genes and Genomes pathway, subcellular location, transmembrane topology, protein secondary/tertiary structure, homologous group and transcription factor. SGID provides genome-scale visualization of population genetics test results based on high-depth resequencing data of 158 silkworm samples. It also provides interactive analysis tools of transcriptomic and epigenomic data from 79 NCBI BioProjects. SGID will be extremely useful to silkworm research in the future.

The emerging molecular mechanism of m6A modulators in tumorigenesis and cancer progression

N⁶-methyladenosine (m⁶A) is the most abundant RNA modification; m⁶A modifications are installed by methyltransferases, removed by demethylases and recognized by reader proteins. M⁶A plays crucial roles in a variety of biological processes by regulating target RNA translation, splicing, nuclear export, and decay. Since the establishment of methylated RNA immunoprecipitation-sequencing methodology, over three hundred articles about m⁶A modulators, including "writers", "erasers" and "readers", have been reported in the last four years. In addition, an increasing number of molecular mechanisms underlying m⁶A RNA methylation in human cancers have been comprehensively clarified. The recently emerged molecular mechanisms of m⁶A modulators in cancer cell proliferation, cell cycle progression, migration and invasion, apoptosis, and autophagy remain to be summarized. Hence, this review specifically summarizes these recent advances in the understanding of m⁶A molecular mechanisms in tumorigenesis and cancer progression. In addition, we discuss the prospect of using an m⁶A methylation modulator as a new diagnostic biomarker and therapeutic target for human cancers.

N 6-Methyladenosine Level in Silkworm Midgut/Ovary Cell Line Is Associated With Bombyx mori Nucleopolyhedrovirus Infection

Bombyx mori nucleopolyhedrovirus (BmNPV) is one of the most serious pathogens in sericulture and causes huge economic loss annually. The roles of N6-methyladenosine (m6A) modification in silkworms following BmNPV infection are currently unclear. Here, methylated RNA immunoprecipitation with next-generation sequencing were applied to investigate the m6A profiles in silkworm midgut following BmNPV infection. A total of 9144 and 7384 m6A peaks were identified from the BmNPV-infected (TEST) and uninfected silkworm midguts (CON), respectively, which were distributed predominantly near stop codons. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of common m6A peaks in nuclear genes revealed that these m6A-related transcripts were associated with crucial signaling pathways. Comparative transcriptome analysis showed that 1221 differential expressed m6A peaks were identified between TEST and CON, indicating that m6A modification is regulated following BmNPV infection. GO and KEGG pathway analysis of the differentially expressed m6A peaks showed their association with signal transduction, translation, and degradation. To understand further the effect of the m6A machinery on virus infection, expression levels of m6A-related genes were altered in silencing and overexpression experiments. Expression of viral structural protein VP39 was increased in BmN cells by siRNA-mediated depletion of methyltransferase-like (METTL) enzyme genes (BmMETTL3, BmMETTL14) and cytoplasmic YTH-domain family 3 (BmYTHDF3), while the reverse results were found after overexpression of the m6A-related enzymes in BmN cells. Overall, m6A modification might be a novel epigenetic mechanism that regulation BmNPV infection and interference with this mechanism may provide a novel antiviral strategy for preventing BmNPV disease.

Heat Shock Cognate 70 Functions as A Chaperone for the Stability of Kinetochore Protein CENP-N in Holocentric Insect Silkworms

The centromere, in which kinetochore proteins are assembled, plays an important role in the accurate congression and segregation of chromosomes during cell mitosis. Although the function of the centromere and kinetochore is conserved from monocentric to holocentric, the DNA sequences of the centromere and components of the kinetochore are varied among different species. Given the lack of core centromere protein A (CENP-A) and CENP-C in the lepidopteran silkworm Bombyx mori, which possesses holocentric chromosomes, here we investigated the role of CENP-N, another important member of the centromere protein family essential for kinetochore assembly. For the first time, cellular localization and RNA interference against CENP-N have confirmed its kinetochore function in silkworms. To gain further insights into the regulation of CENP-N in the centromere, we analyzed the affinity-purified complex of CENP-N by mass spectrometry and identified 142 interacting proteins. Among these factors, we found that the chaperone protein heat shock cognate 70 (HSC70) is able to regulate the stability of CENP-N by prohibiting ubiquitin-proteasome pathway, indicating that HSC70 could control cell cycle-regulated degradation of CENP-N at centromeres. Altogether, the present work will provide a novel clue to understand the regulatory mechanism for the kinetochore activity of CENP-N during the cell cycle.

Epigenetic Methylations on N6-Adenine and N6-Adenosine with the same Input but Different Output

Epigenetic modifications on individual bases in DNA and RNA can encode inheritable genetic information beyond the canonical bases. Among the nucleic acid modifications, DNA N6-methadenine (6mA) and RNA N6-methyladenosine (m6A) have recently been well-studied due to the technological development of detection strategies and the functional identification of modification enzymes. The current findings demonstrate a wide spectrum of 6mA and m6A distributions from prokaryotes to eukaryotes and critical roles in multiple cellular processes. It is interesting that the processes of modification in which the methyl group is added to adenine and adenosine are the same, but the outcomes of these modifications in terms of their physiological impacts in organisms are quite different. In this review, we summarize the latest progress in the study of enzymes involved in the 6mA and m6A methylation machinery, including methyltransferases and demethylases, and their functions in various biological pathways. In particular, we focus on the mechanisms by which 6mA and m6A regulate the expression of target genes, and we highlight the future challenges in epigenetic regulation.