Human C-to-U Coding RNA Editing Is Largely Nonadaptive

Can transcriptome size and off-target effects explain the contrasting evolution of mitochondrial vs nuclear RNA editing?

Mitochondrial RNA editing has evolved independently in numerous eukaryotic lineages, where it generally restores conserved sequences and functional reading frames in mRNA transcripts derived from altered or disrupted mitochondrial protein-coding genes. In contrast to this "restorative" RNA editing in mitochondria, most editing of nuclear mRNAs introduces novel sequence variants and diversifies the proteome. This Perspective addresses the hypothesis that these completely opposite effects of mitochondrial vs. nuclear RNA editing arise from the enormous difference in gene number between the respective genomes. Because mitochondria produce a much smaller transcriptome, they likely create less opportunity for off-target editing, which has been supported by recent experimental work expressing mitochondrial RNA editing machinery in foreign contexts. In addition, there is recent evidence that the size and complexity of RNA targets may slow the kinetics and reduce efficiency of on-target RNA editing. These findings suggest that efficient targeting and a low risk of off-target editing have facilitated the repeated emergence of disrupted mitochondrial genes and associated restorative RNA editing systems via (potentially non-adaptive) evolutionary pathways that are not feasible in larger nuclear transcriptomes due to lack of precision.

Constraints on the optimization of gene product diversity

RNA and proteins can have diverse isoforms due to post-transcriptional and post-translational modifications. A fundamental question is whether these isoforms are mostly beneficial or the result of noisy molecular processes. To assess the plausibility of these explanations, we developed mathematical models depicting different regulatory architectures and investigated isoform evolution under multiple population genetic regimes. We found that factors beyond selection, such as effective population size and the number of cis-acting loci, significantly influence evolutionary outcomes. We found that sub-optimal phenotypes are more likely to evolve when populations are small and/or when the number of cis-loci is large. We also discovered that opposing selection on cis- and trans-acting loci can constrain adaptation, leading to a non-monotonic relationship between effective population size and optimization. More generally, our models provide a quantitative framework for developing statistical tests to analyze empirical data; as a demonstration of this, we analyzed A-to-I RNA editing levels in coleoids and found these to be largely consistent with non-adaptive explanations.

BIRC3 RNA Editing Modulates Lipopolysaccharide-Induced Liver Inflammation: Potential Implications for Animal Health

Animals and humans are frequently infected by bacteria or exposed to bacterial derivatives in contaminated food, drinking water, or air, which significantly impacts their health. Among these bacterial sources, LPS (lipopolysaccharide) is the primary culprit. While it is widely known that LPS can cause liver inflammation and damage in animals, few studies have investigated this mechanism from the perspective of RNA editing. In this study, we administered LPS to mice via gavage to induce a liver injury model. We then used RNA editing omics approaches (RE-seq) to analyze RNA editing events potentially leading to liver inflammation following LPS administration, aiming to reveal the crucial role of RNA editing in LPS-induced processes. At the RNA editing level, we observed significant differences between the LPS group and the control (CON) group. Specifically, we identified 354 differentially edited genes, with 192 upregulated and 162 downregulated. These differentially edited genes were significantly enriched in pathways related to apoptosis, mTOR signaling, oxidative stress, and Nf-Kappa B signaling. By further integrating gene expression profiles and using a nine-quadrant analysis, we identified an important gene, Birc3, which showed significantly higher editing and expression levels in the LPS group. This gene is directly linked to liver inflammation and damage. The RNA editing of Birc3 represents a significant potential mechanism underlying LPS-induced liver damage, providing a novel approach for addressing animal and human health issues.

WTAP-Mediated N6-Methyladenosine Modification Promotes Gastric Cancer Progression by Regulating MAP2K6 Expression

Wilms tumor 1 associated protein (WTAP) is a key RNA N6-methyladenosine (m6A) methylase, which is involved in gastric cancer (GC) development, but its pathogenic mechanism is not clear. This study aims to thoroughly explore the underlying molecular mechanism of WTAP-mediated m6A modification in GC pathogenesis. qRT-PCR and immunohistochemistry showed that significantly elevated WTAP expression in GC tissues and is related to advanced age, poorly differentiation, lymph node metastasis and high TNM stage. Overexpression and knockdown of WTAP could promote or inhibit the proliferation, migration and invasion of GC cells in vitro, furthermore, suppression of WTAP expression impeded the growth of xenograft tumors in vivo. Utilizing RNA sequencing, methylated RNA immunoprecipitation (MeRIP) sequencing and bioinformatics analysis, we identified MAP2K6 as direct downstream target of WTAP with m6A modification in GC. The interaction between WTAP and MAP2K6 was confirmed by MeRIP-qPCR, luciferase reporter assay, Co-IF and bioinformatics prediction. Immunofluorescence and rescue studies were performed to verify WTAP-mediated m6A modification promotes the proliferation, migration, and invasion of GC cells by positively regulating the target gene MAP2K6. This underscores the potential therapeutic significance of targeting the WTAP-MAP2K6 axis in combating GC occurrence and progression.

Research progress on N6-methyladenosine RNA modification in osteosarcoma: functions, mechanisms, and potential clinical applications

Osteosarcoma (OS) is the most commonly diagnosed primary malignant bone tumor in children and adolescents. Despite significant advancements in therapeutic strategies against OS over the past few decades, the prognosis for this disease remains poor, largely due to its high invasiveness and challenges associated with its treatment. N6-methyladenosine (m6A) modification is one of the most abundant epigenetic modifications of RNAs, and many studies have highlighted its crucial role in OS. This article provides a comprehensive summary and introduction to m6A regulators, including methyltransferases, demethylases, and binding proteins. The article emphasizes how regulated m6A modifications can either promote or inhibit OS. It also delves into the mechanisms by which m6A-modified messenger RNAs (mRNAs) and noncoding RNAs (ncRNAs) participate in signaling pathways such as the Wnt/β-catenin, PI3K/AKT, and STAT3 pathways, and discusses these mechanisms in detail. Given the abnormal expression of m6A regulators in OS, the article also explores their potential applications as biomarkers or therapeutic targets in clinical settings. It is anticipated that this review will provide new insights into the diagnosis and treatment of OS.

Most m5C Modifications in Mammalian mRNAs are Nonadaptive

5-Methylation (m5C) on mRNA molecules is a prevalent internal posttranscriptional modification in eukaryotes. Although m5C modification has been reported to regulate some biological processes, whether most mRNA m5C modifications are functional is unknown. To address this question, we analyzed the genome-wide evolutionary characteristics of m5C modifications in protein-coding genes of humans and mice. Our analysis of RNA sequencing data from 13 tissues of both species revealed that (i) the occurrence of m5C modification is exceedingly low, (ii) the fraction of m5Cs decreases with the amount of Cs across genes or tissues, (iii) m5C modifications are mostly unshared between species, and (iv) m5C sites and motifs do not exhibit greater evolutionary conservation. Additionally, we estimate that a large fraction of the observed mRNA m5C modifications may be deleterious. Together, these observations suggest that most m5C modifications in mammalian mRNAs are nonadaptive, which has important implications for understanding the biological significance of m5C and other posttranscriptional modifications.

The m6A and immune regulatory gene signature predicts the prognosis and correlates with immune infiltration of head and neck squamous cell carcinoma

Recent investigations have underscored the epigenetic modulation of the immune response; however, the interplay between RNA N6-methyladenosine (m6A) modification and immunomodulation in head and neck squamous cell carcinoma (HNSC) remains relatively unexplored. To bridge this knowledge gap, we undertook an extensive examination of the potential contributions of m6A modification and immunomodulation in HNSC. We amalgamated and deduplicated 27 m6A -related genes (m6AGs) and 1342 immune regulation-related genes (IMRGs), resulting in a comprehensive dataset encompassing 1358 genes. This dataset was scrutinized for m6A modification and immunomodulatory patterns within HNSC specimens. Employing Cox regression analysis and the Least Absolute Shrinkage and Selection Operator (LASSO) technique, we developed a prognostic risk model for m6A regulator-mediated methylation modification and immunomodulation-related differentially expressed genes (m6A&IMRDEGs). Our differential expression analysis delineated 29 m6A&IMRDEGs, and Weighted Gene Co-expression Network Analysis (WGCNA) elucidated two module genes (IL11 and MMP13) subjected to correlation analysis. The prognostic prediction models revealed that the clinical predictive efficacy peaked for 1-year forecasts, followed sequentially by 3-year and 5-year predictions. The risk scores derived from the model adeptly categorized HNSC patients into high- and low-risk cohorts, with the high-risk group exhibiting a more unfavorable prognosis. Protein-Protein Interaction (PPI) analysis identified 7 hub genes implicated in m6A and immune regulation, namely BPIFB1, BPIFB2, GP2, MUC5B, MUC7, PIP, and SCGB3A1. Furthermore, we noted marked disparities in the expression profiles of 18 immune cell types between the high- and low-risk groups. Our results substantiate that the clustering subpopulations and risk models associated with m6A and immune regulatory genes portend a poor prognosis in HNSC. The risk score emerges as a potent prognostic biomarker and predictive metric for HNSC patients. A thorough assessment of m6A and immune regulatory genes in HNSC will augment our comprehension of the tumor immune microenvironment and facilitate the advancement of HNSC therapeutics.

METTL3-mediated m6A modification promotes chemoresistance of intrahepatic cholangiocarcinoma by up-regulating NRF2 to inhibit ferroptosis in cisplatin-resistant cells

This study explores the relationship between m6A modification and ferroptosis in intrahepatic cholangiocarcinoma (ICC) and its impact on cisplatin resistance. We established cisplatin-resistant cells. CCK-8 and Transwell assays were conducted to evaluate the effects of METTL3 on drug resistance, migration, and invasion. RT-qPCR and Western blotting were used to measure target gene expression and the effects of overexpression and suppression. RIP, luciferase reporter assay, and other experiments were utilized to investigate the interaction between METTL3 and NRF2. Additionally, rescue experiments were performed to confirm the role of the METTL3/NRF2 axis in tumor drug resistance. METTL3 was found to be highly expressed in cisplatin-resistant cells, enhancing m6A modification levels, stabilizing NRF2 mRNA, and increasing NRF2 protein expression to inhibit ferroptosis. These findings indicate that the METTL3/NRF2 axis inhibits ferroptosis in cisplatin-resistant cells, thereby promoting chemotherapy resistance in ICC. This provides a potential direction for future research and treatment of ICC.

Regulatory effect of N6-methyladenosine on tumor angiogenesis

Previous studies have demonstrated that genetic alterations governing epigenetic processes frequently drive tumor development and that modifications in RNA may contribute to these alterations. In the 1970s, researchers discovered that N6-methyladenosine (m6A) is the most prevalent form of RNA modification in advanced eukaryotic messenger RNA (mRNA) and noncoding RNA (ncRNA). This modification is involved in nearly all stages of the RNA life cycle. M6A modification is regulated by enzymes known as m6A methyltransferases (writers) and demethylases (erasers). Numerous studies have indicated that m6A modification can impact cancer progression by regulating cancer-related biological functions. Tumor angiogenesis, an important and unregulated process, plays a pivotal role in tumor initiation, growth, and metastasis. The interaction between m6A and ncRNAs is widely recognized as a significant factor in proliferation and angiogenesis. Therefore, this article provides a comprehensive review of the regulatory mechanisms underlying m6A RNA modifications and ncRNAs in tumor angiogenesis, as well as the latest advancements in molecular targeted therapy. The aim of this study is to offer novel insights for clinical tumor therapy.

METTL3-Mediated Maturation of miR-99a-5p Promotes Cell Migration and Invasion in Oral Squamous Cell Carcinoma by Targeting ZBTB7A

METTL3 is an important methyltransferase in N(6)-methyladenosine (m6A) modification. Recently, METTL3 mediates methylation of pri-microRNA (miRNA) to accelerate miRNA maturation, regulating tumor development. This study explored whether METTL3 mediated miR-99a-5p to influence oral squamous cell carcinoma (OSCC) cell metastasis. MiR-99a-5p, ZBTB7A, and MATTL3 expression was measured using quantitative real-time PCR. Biological behaviors were assessed using cell counting kit-8, flow cytometry, Transwell assay, as well as western blot. Luciferase reporter assay evaluated the interaction between miR-99a-5p and ZBTB7A. METTL3-regulated pri-miR-99a-5p processing was determined by RNA binding protein immunoprecipitation (RIP) and methylated RNA immunoprecipitation (MeRIP) assays. The consequences clarified that miR-99a-5p was upregulated in OSCC cells. Downregulation of miR-99a-5p suppressed cellular viability, migration, invasion, and epithelial-mesenchymal transition (EMT), and induced apoptosis. ZBTB7A acted as a miR-99a-5p target and reversed the effects on cellular behaviors induced by miR-99a-5p inhibitor. m6A content and METTL3 expression were increased in OSCC cells. METTL3 promoted the m6A modification of pri-miR-99a-5p and thereby facilitated miR-99a-5p processing. Moreover, knockdown of METTL3 inhibited OSCC metastasis by downregulating miR-99a-5p. Taken together, METTL3 promoted miR-99a-5p maturation in an m6A-dependent manner, which further targets ZBTB7A to accelerate the progression of OSCC. These findings suggest potential targets for OSCC therapy.

Genetic and selective constraints on the optimization of gene product diversity

RNA and protein expressed from the same gene can have diverse isoforms due to various post-transcriptional and post-translational modifications. For the vast majority of alternative isoforms, It is unknown whether they are adaptive or simply biological noise. As we cannot experimentally probe the function of each isoform, we can ask whether the distribution of isoforms across genes and across species is consistent with expectations from different evolutionary processes. However, there is currently no theoretical framework that can generate such predictions. To address this, we developed a mathematical model where isoform abundances are determined collectively by cis-acting loci, trans-acting factors, gene expression levels, and isoform decay rates to predict isoform abundance distributions across species and genes in the face of mutation, genetic drift, and selection. We found that factors beyond selection, such as effective population size and the number of cis-acting loci, significantly influence evolutionary outcomes. Notably, suboptimal phenotypes are more likely to evolve when the population is small and/or when the number of cis-loci is large. We also explored scenarios where modification processes have both beneficial and detrimental effects, revealing a non-monotonic relationship between effective population size and optimization, demonstrating how opposing selection pressures on cis- and trans-acting loci can constrain the optimization of gene product diversity. As a demonstration of the power of our theory, we compared the expected distribution of A-to-I RNA editing levels in coleoids and found this to be largely consistent with non-adaptive explanations.

GTDrift: a resource for exploring the interplay between genetic drift, genomic and transcriptomic characteristics in eukaryotes

We present GTDrift, a comprehensive data resource that enables explorations of genomic and transcriptomic characteristics alongside proxies of the intensity of genetic drift in individual species. This resource encompasses data for 1506 eukaryotic species, including 1413 animals and 93 green plants, and is organized in three components. The first two components contain approximations of the effective population size, which serve as indicators of the extent of random genetic drift within each species. In the first component, we meticulously investigated public databases to assemble data on life history traits such as longevity, adult body length and body mass for a set of 979 species. The second component includes estimations of the ratio between the rate of non-synonymous substitutions and the rate of synonymous substitutions (dN/dS) in protein-coding sequences for 1324 species. This ratio provides an estimate of the efficiency of natural selection in purging deleterious substitutions. Additionally, we present polymorphism-derived N e estimates for 66 species. The third component encompasses various genomic and transcriptomic characteristics. With this component, we aim to facilitate comparative transcriptomics analyses across species, by providing easy-to-use processed data for more than 16 000 RNA-seq samples across 491 species. These data include intron-centered alternative splicing frequencies, gene expression levels and sequencing depth statistics for each species, obtained with a homogeneous analysis protocol. To enable cross-species comparisons, we provide orthology predictions for conserved single-copy genes based on BUSCO gene sets. To illustrate the possible uses of this database, we identify the most frequently used introns for each gene and we assess how the sequencing depth available for each species affects our power to identify major and minor splice variants.

Random genetic drift sets an upper limit on mRNA splicing accuracy in metazoans

Most eukaryotic genes undergo alternative splicing (AS), but the overall functional significance of this process remains a controversial issue. It has been noticed that the complexity of organisms (assayed by the number of distinct cell types) correlates positively with their genome-wide AS rate. This has been interpreted as evidence that AS plays an important role in adaptive evolution by increasing the functional repertoires of genomes. However, this observation also fits with a totally opposite interpretation: given that 'complex' organisms tend to have small effective population sizes (Ne), they are expected to be more affected by genetic drift, and hence more prone to accumulate deleterious mutations that decrease splicing accuracy. Thus, according to this 'drift barrier' theory, the elevated AS rate in complex organisms might simply result from a higher splicing error rate. To test this hypothesis, we analyzed 3496 transcriptome sequencing samples to quantify AS in 53 metazoan species spanning a wide range of Ne values. Our results show a negative correlation between Ne proxies and the genome-wide AS rates among species, consistent with the drift barrier hypothesis. This pattern is dominated by low abundance isoforms, which represent the vast majority of the splice variant repertoire. We show that these low abundance isoforms are depleted in functional AS events, and most likely correspond to errors. Conversely, the AS rate of abundant isoforms, which are relatively enriched in functional AS events, tends to be lower in more complex species. All these observations are consistent with the hypothesis that variation in AS rates across metazoans reflects the limits set by drift on the capacity of selection to prevent gene expression errors.

Role of N6-methyladenosine methylation in glioma: recent insights and future directions

Glioma is the most pervasive intracranial tumor in the central nervous system (CNS), with glioblastoma (GBM) being the most malignant type having a highly heterogeneous cancer cell population. There is a significantly high mortality rate in GBM patients. Molecular biomarkers related to GBM malignancy may have prognostic values in predicting survival outcomes and therapeutic responses, especially in patients with high-grade gliomas. In particular, N6-methyladenine (m6A) mRNA modification is the most abundant form of post-transcriptional RNA modification in mammals and is involved in regulating mRNA translation and degradation. Cumulative findings indicate that m6A methylation plays a crucial part in neurogenesis and glioma pathogenesis. In this review, we summarize recent advances regarding the functional significance of m6A modification and its regulatory factors in glioma occurrence and progression. Significant advancement of m6A methylation-associated regulators as potential therapeutic targets is also discussed.

Identification and Interpretation of A-to-I RNA Editing Events in Insect Transcriptomes

Adenosine-to-inosine (A-to-I) RNA editing is the most prevalent RNA modification in the nervous systems of metazoans. To study the biological significance of RNA editing, we first have to accurately identify these editing events from the transcriptome. The genome-wide identification of RNA editing sites remains a challenging task. In this review, we will first introduce the occurrence, regulation, and importance of A-to-I RNA editing and then describe the established bioinformatic procedures and difficulties in the accurate identification of these sit esespecially in small sized non-model insects. In brief, (1) to obtain an accurate profile of RNA editing sites, a transcriptome coupled with the DNA resequencing of a matched sample is favorable; (2) the single-cell sequencing technique is ready to be applied to RNA editing studies, but there are a few limitations to overcome; (3) during mapping and variant calling steps, various issues, like mapping and base quality, soft-clipping, and the positions of mismatches on reads, should be carefully considered; (4) Sanger sequencing of both RNA and the matched DNA is the best verification of RNA editing sites, but other auxiliary evidence, like the nonsynonymous-to-synonymous ratio or the linkage information, is also helpful for judging the reliability of editing sites. We have systematically reviewed the understanding of the biological significance of RNA editing and summarized the methodology for identifying such editing events. We also raised several promising aspects and challenges in this field. With insightful perspectives on both scientific and technical issues, our review will benefit the researchers in the broader RNA editing community.

Leveraging molecular quantitative trait loci to comprehend complex diseases/traits from the omics perspective

Comprehending the molecular basis of quantitative genetic variation is a principal goal for complex diseases or traits. Molecular quantitative trait loci (molQTLs) have made it possible to investigate the effects of genetic variants hiding behind large-scale omics data. A deeper understanding of molQTL is urgently required in light of the multi-dimensionalization of omics data to more fully elucidate the pertinent biological mechanisms. Herein, we reviewed molQTLs with the corresponding resource from the omics perspective and further discussed the integrative strategy of GWAS-molQTL to infer their causal effects. Subsequently, we described the opportunities and challenges encountered by molQTL. The case studies showed that molQTL is essential for complex diseases and traits, whether single- or multi-omics QTLs. Overall, we highlighted the functional significance of genetic variants to employ the discovery of molQTL in complex diseases and traits.

Differential adaptive RNA editing signals between insects and plants revealed by a new measurement termed haplotype diversity

BACKGROUND

C-to-U RNA editing in plants is believed to confer its evolutionary adaptiveness by reversing unfavorable DNA mutations. This "restorative hypothesis" has not yet been tested genome-wide. In contrast, A-to-I RNA editing in insects like Drosophila and honeybee is already known to benefit the host by increasing proteomic diversity in a spatial-temporal manner (namely "diversifying hypothesis").

METHODS

We profiled the RNA editomes of multiple tissues of Arabidopsis thaliana, Drosophila melanogaster, and Apis melifera. We unprecedentedly defined the haplotype diversity (HD) of RNA molecules based on nonsynonymous editing events (recoding sites).

RESULTS

Signals of adaptation is confirmed in Arabidopsis by observing higher frequencies and levels at nonsynonymous editing sites over synonymous sites. Compared to A-to-I recoding sites in Drosophila, the C-to-U recoding sites in Arabidopsis show significantly lower HD, presumably due to the stronger linkage between C-to-U events.

CONCLUSIONS

C-to-U RNA editing in Arabidopsis is adaptive but it is not designed for diversifying the proteome like A-to-I editing in Drosophila. Instead, C-to-U recoding sites resemble DNA mutations. Our observation supports the restorative hypothesis of plant C-to-U editing which claims that editing is used for fixing unfavorable genomic sequences.

The Repertoire of RNA Modifications Orchestrates a Plethora of Cellular Responses

Although a plethora of DNA modifications have been extensively investigated in the last decade, recent breakthroughs in molecular biology, including high throughput sequencing techniques, have enabled the identification of post-transcriptional marks that decorate RNAs; hence, epitranscriptomics has arisen. This recent scientific field aims to decode the regulatory layer of the transcriptome and set the ground for the detection of modifications in ribose nucleotides. Until now, more than 170 RNA modifications have been reported in diverse types of RNA that contribute to various biological processes, such as RNA biogenesis, stability, and transcriptional and translational accuracy. However, dysfunctions in the RNA-modifying enzymes that regulate their dynamic level can lead to human diseases and cancer. The present review aims to highlight the epitranscriptomic landscape in human RNAs and match the catalytic proteins with the deposition or deletion of a specific mark. In the current review, the most abundant RNA modifications, such as N6-methyladenosine (m⁶A), N5-methylcytosine (m⁵C), pseudouridine (Ψ) and inosine (I), are thoroughly described, their functional and regulatory roles are discussed and their contributions to cellular homeostasis are stated. Ultimately, the involvement of the RNA modifications and their writers, erasers, and readers in human diseases and cancer is also discussed.

Lactobacillus for ribosome peptide editing cancer

This study reviews newly discovered insect peptide point mutations as new possible cancer research targets. To interpret newly discovered peptide point mutations in insects as new possible cancer research targets, we focused on the numerous peptide changes found in the 'CSP' family on the sex pheromone gland of the female silkworm moth Bombyx mori. We predict that the Bombyx peptide modifications will have a significant effect on cancer CUP (cancers of unknown primary) therapy and that bacterial peptide editing techniques, specifically Lactobacillus combined to CRISPR, will be used to regulate ribosomes and treat cancer in humans.

METTL16 regulates m6A methylation on chronic hepatitis B associated gene HLA-DPB1 involved in liver fibrosis

The role of genetic factors in the occurrence and progression of CHB (CHB) is still not fully explored. In recent years, genome-wide association studies on CHB patients have demonstrated that a large number of CHB-associated single nucleotide polymorphisms exist in the gene intron, which may regulate expression at the transcriptional level. Modification of RNA m6A methylation is one of the key mechanisms regulating gene expression. Here we show that METTL16, an m6A regulator involved in mRNA intron splicing, is differentially expressed in CHB the tissue of patients who has definite diagnosis of mild and severe fibrosis. At the same time, there are also significant differences in the expression of CHB-associated genes such as HLA-DPA1 and HLA-DPB1. The expression of HLA-DPB1 is related to METTL16. Furthermore, analyses of RNA binding of METTL16 and HLA-DPB1 show that the silencing of METTL16 in astrocytes downregulates m6A and expression of HLA-DPB1. In conclusion, METTL16 participates in the progression of CHB fibrosis by regulating the m6A level and expression of HLA-DPB1.

Gene product diversity: adaptive or not?

One gene does not equal one RNA or protein. The genomic revolution has revealed numerous different RNA and protein molecules that can be produced from one gene, such as circular RNAs generated by back-splicing, proteins with residues mismatching the genomic encoding because of RNA editing, and proteins extended in the C terminus via stop codon readthrough in translation. Are these diverse products results of exquisite gene regulations or imprecise biological processes? While there are cases where the gene product diversity appears beneficial, genome-scale patterns suggest that much of this diversity arises from nonadaptive, molecular errors. This finding has important implications for studying the functions of diverse gene products and for understanding the fundamental properties and evolution of cellular life.

Building integrative functional maps of gene regulation

Every cell in the human body inherits a copy of the same genetic information. The three billion base pairs of DNA in the human genome, and the roughly 50 000 coding and non-coding genes they contain, must thus encode all the complexity of human development and cell and tissue type diversity. Differences in gene regulation, or the modulation of gene expression, enable individual cells to interpret the genome differently to carry out their specific functions. Here we discuss recent and ongoing efforts to build gene regulatory maps, which aim to characterize the regulatory roles of all sequences in a genome. Many researchers and consortia have identified such regulatory elements using functional assays and evolutionary analyses; we discuss the results, strengths and shortcomings of their approaches. We also discuss new techniques the field can leverage and emerging challenges it will face while striving to build gene regulatory maps of ever-increasing resolution and comprehensiveness.

C-to-U RNA Editing: A Site Directed RNA Editing Tool for Restoration of Genetic Code

The restoration of genetic code by editing mutated genes is a potential method for the treatment of genetic diseases/disorders. Genetic disorders are caused by the point mutations of thymine (T) to cytidine (C) or guanosine (G) to adenine (A), for which gene editing (editing of mutated genes) is a promising therapeutic technique. In C-to-Uridine (U) RNA editing, it converts the base C-to-U in RNA molecules and leads to nonsynonymous changes when occurring in coding regions; however, for G-to-A mutations, A-to-I editing occurs. Editing of C-to-U is not as physiologically common as that of A-to-I editing. Although hundreds to thousands of coding sites have been found to be C-to-U edited or editable in humans, the biological significance of this phenomenon remains elusive. In this review, we have tried to provide detailed information on physiological and artificial approaches for C-to-U RNA editing.

m6A Topological Transition Coupled to Developmental Regulation of Gene Expression During Mammalian Tissue Development

N6-methyladenosine (m6A) is the most prevalent internal modification and reversible epitranscriptomic mark in messenger RNAs (mRNAs) and plays essential roles in a variety of biological processes. However, the dynamic distribution patterns of m6A and their significance during mammalian tissue development are poorly understood. Here, we found that based on m6A distribution patterns, protein-coding genes were classified into five groups with significantly distinct biological features and functions. Strikingly, comparison of the m6A methylomes of multiple mammalian tissues between fetal and adult stages revealed dynamic m6A topological transition during mammalian tissue development, and identified large numbers of genes with significant m6A loss in 5′UTRs or m6A gain around stop codons. The genes with m6A loss in 5′UTRs were highly enriched in developmental stage-specific genes, and their m6A topological transitions were strongly associated with gene expression regulation during tissue development. The genes with m6A gain around the stop codons were associated with tissue-specific functions. Our findings revealed the existence of different m6A topologies among protein-coding genes that were associated with distinct characteristics. More importantly, these genes with m6A topological transitions were crucial for tissue development via regulation of gene expression, suggesting the importance of dynamic m6A topological transitions during mammalian tissue development.

Comprehensive analyses indicated the association between m6A related long non-coding RNAs and various pathways in glioma

BACKGROUND

Glioma is one of the most malignant brain tumors and diseases. N6-methyladenosine modification (m6A) is the most abundant and prevalent internal chemical modification of mRNA and long non-coding RNAs (lncRNAs) in eukaryotes. Nevertheless, the correlated pathways and clinical utilization of m6A-related lncRNAs have not been fully evaluated in glioma.

METHODS

Public RNA-sequencing and clinical annotation data were retrieved from TCGA, CGGA and GEO database. Differential expression analysis and univariate Cox regression analysis were performed to identify the m6A-related and differentially expressed lncRNAs with prognostic function (m6A-DELPF). The consensus clustering was performed to identify the expression pattern of m6A-DELPF. LASSO Cox regression analysis was performed to construct the lncRNA-based signature. The CIBERSORT and ESTIMATE algorithms were performed to analyze immune infiltration and tumor microenvironment, respectively. Immunotherapy sensitivity analysis was performed using data from TCIA. The small molecule drugs prediction analysis was performed using The Connectivity Map (CMap) database and STITCH database. A competing endogenous RNAs (ceRNA) network was constructed based on miRcode, miRDB, miRTarBase, TargetScan database.

RESULTS

Two clusters (cluster1 and cluster2) were identified after unsupervised cluster analysis based on m6A-DELPF. Additionally, a 15-gene prognostic signature namely m6A-DELPFS was constructed. Analyses of epithelial-mesenchymal-transition score, tumor microenvironment, immune infiltration, clinical characterization analysis, and putative drug prediction were performed to confirm the clinical utility and efficacy of m6A-DELPFS. The potential mechanisms including tumor immune microenvironment of m6A-DELPF influence the initiation and progression of glioma. A clinically accessible nomogram was also constructed based on the m6A-DELPF and other survival-relevant clinical parameters. Two miRNAs and 114 mRNAs were identified as the downstream of seven m6A-related lncRNAs in a ceRNA network.

CONCLUSION

Our present research confirmed the clinical value of m6A related lncRNAs and their high correlation with tumor immunity, tumor microenvironment, tumor mutation burden and drug sensitivity in glioma.

Analysis of the expression patterns and clinical relevance of m6A regulators in 33 cancer types

Background: The role of N6-methyladenine (m6A) RNA methylation in a variety of biological processes is gradually being revealed. Methods: Here, we systematically describe the correlation between the expression pattern of m6A RNA methylation regulatory factors and clinical phenotype, immunity, drug sensitivity, stem cells and prognosis in more than 10,000 samples of 33 types of cancer. Results: The results show that there are significant differences in the expression of 20 m6A RNA methylation regulatory factors in different cancers, and there was a significant correlation with the analysis indicators. Conclusion: In this study, the m6A RNA methylation regulatory factor was found not only to potentially assist in stratifying the prognosis but also to predict or improve the sensitivity of clinical drug therapy.

Fast evolution of SARS-CoV-2 driven by deamination systems in hosts

As an RNA virus, the fast evolution of SARS-CoV-2 is driven by the extensive RNA deamination by the host cells.

EditPredict: Prediction of RNA editable sites with convolutional neural network

RNA editing exerts critical impacts on numerous biological processes. While millions of RNA editings have been identified in humans, much more are expected to be discovered. In this work, we constructed Convolutional Neural Network (CNN) models to predict human RNA editing events in both Alu regions and non-Alu regions. With a validation dataset resulting from CRISPR/Cas9 knockout of the ADAR1 enzyme, the validation accuracies reached 99.5% and 93.6% for Alu and non-Alu regions, respectively. We ported our CNN models in a web service named EditPredict. EditPredict not only works on reference genome sequences but can also take into consideration single nucleotide variants in personal genomes. In addition to the human genome, EditPredict tackles other model organisms including bumblebee, fruitfly, mouse, and squid genomes. EditPredict can be used stand-alone to predict novel RNA editing and it can be used to assist in filtering for candidate RNA editing detected from RNA-Seq data.

Large-scale omics data reveal the cooperation of mutation-circRNA-miRNA-target gene network in liver cancer oncogenesis

Aims: Clarifying the initial trigger of the differentially expressed genes in cancers helps researchers understand the cellular system as a whole network. Materials & methods: We retrieve the transcriptome and translatome of tumor and normal tissues from ten liver cancer patients and define differentially expressed genes and tumor-specific mutations. We associate the oncogenesis with the mutations by target prediction and experimental verification. Results: Upregulated genes have tumor-specific mutations in 3'UTRs that abolish the binding of miRNAs. For downregulated genes, their corresponding miRNAs are mutually targeted by two circRNAs, with mutations in base-pairing regions. Transfection experiments support the oncogenic role of these mutations. Conclusions: The tumor-specific mutations serve as the initial trigger of liver cancer. The mutation-circRNA-miRNA-target gene chain is completed.

Comparative genomic analysis of a naturally born serpentized pig reveals putative mutations related to limb and bone development

BACKGROUND

It is believed that natural selection acts on the phenotypical changes caused by mutations. Phenotypically, from fishes to amphibians to reptiles, the emergence of limbs greatly facilitates the landing of ancient vertebrates, but the causal mutations and evolutionary trajectory of this process remain unclear.

RESULTS

We serendipitously obtained a pig of limbless phenotype. Mutations specific to this handicapped pig were identified using genome re-sequencing and comparative genomic analysis. We narrowed down the causal mutations to particular chromosomes and even several candidate genes and sites, such like a mutation-containing codon in gene BMP7 (bone morphogenetic protein) which was conserved in mammals but variable in lower vertebrates.

CONCLUSIONS

We parsed the limbless-related mutations in the light of evolution. The limbless pig shows phenocopy of the clades before legs were evolved. Our findings might help deduce the emergence of limbs during vertebrate evolution and should be appealing to the broad community of human genetics and evolutionary biology.

Mammalian circular RNAs result largely from splicing errors

Ubiquitous in eukaryotes, circular RNAs (circRNAs) comprise a large class of mostly non-coding RNAs produced by back-splicing. Although some circRNAs have demonstrated biochemical activities, whether most circRNAs are functional is unknown. Here, we test the hypothesis that circRNA production primarily results from splicing error and so is deleterious instead of beneficial. In support of the error hypothesis, our analysis of RNA sequencing data from 11 shared tissues of humans, macaques, and mice finds that (1) back-splicing is much rarer than linear-splicing, (2) the rate of back-splicing diminishes with the splicing amount, (3) the overall prevalence of back-splicing in a species declines with its effective population size, and (4) circRNAs are overall evolutionarily unconserved. We estimate that more than 97% of the observed circRNA production is deleterious. We identify a small number of functional circRNA candidates, and the genome-wide trend strongly suggests that circRNAs are largely non-functional products of splicing errors.

N6-Methyladenosine Modification and Its Regulation of Respiratory Viruses

N6-methyladenosine (m⁶A) is a ubiquitous RNA modification in eukaryotes. It plays important roles in the translocation, stabilization and translation of mRNA. Many recent studies have shown that the dysregulation of m⁶A modification is connected with diseases caused by pathogenic viruses, and studies on the role of m⁶A in virus-host interactions have shown that m⁶A plays a wide range of regulatory roles in the life cycle of viruses. Respiratory viruses are common pathogens that can impose a large disease burden on young children and elderly people. Here, we review the effects of m⁶A modification on respiratory virus replication and life cycle and host immunity against viruses.

The role of the HIF-1α/ALYREF/PKM2 axis in glycolysis and tumorigenesis of bladder cancer

BACKGROUND

As a rate-limiting enzyme of glycolysis, pyruvate kinase muscle isozyme M2 (PKM2) participates in tumor metabolism and growth. The regulatory network of PKM2 in cancer is complex and has not been fully studied in bladder cancer. The 5-methylcytidine (m5C) modification in PKM2 mRNA might participate in the pathogenesis of bladder cancer and need to be further clarified. This study aimed to investigate the biological function and regulatory mechanism of PKM2 in bladder cancer.

METHODS

The expression of PKM2 and Aly/REF export factor (ALYREF) was measured by Western blotting, qRT-PCR, and immunohistochemistry. The bioprocesses of bladder cancer cells were demonstrated by a series of experiments in vitro and in vivo. RNA immunoprecipitation, RNA-sequencing, and dual-luciferase reporter assays were conducted to explore the potential regulatory mechanisms of PKM2 in bladder cancer.

RESULTS

In bladder cancer, we first demonstrated that ALYREF stabilized PKM2 mRNA and bound to its m5C sites in 3'-untranslated regions. Overexpression of ALYREF promoted bladder cancer cell proliferation by PKM2-mediated glycolysis. Furthermore, high expression of PKM2 and ALYREF predicted poor survival in bladder cancer patients. Finally, we found that hypoxia-inducible factor-1alpha (HIF-1α) indirectly up-regulated the expression of PKM2 by activating ALYREF in addition to activating its transcription directly.

CONCLUSIONS

The m5C modification in PKM2 mRNA in the HIF-1α/ALYREF/PKM2 axis may promote the glucose metabolism of bladder cancer, providing a new promising therapeutic target for bladder cancer.

Effective Population Size Predicts Local Rates but Not Local Mitigation of Read-through Errors

In correctly predicting that selection efficiency is positively correlated with the effective population size (Ne), the nearly neutral theory provides a coherent understanding of between-species variation in numerous genomic parameters, including heritable error (germline mutation) rates. Does the same theory also explain variation in phenotypic error rates and in abundance of error mitigation mechanisms? Translational read-through provides a model to investigate both issues as it is common, mostly nonadaptive, and has good proxy for rate (TAA being the least leaky stop codon) and potential error mitigation via "fail-safe" 3' additional stop codons (ASCs). Prior theory of translational read-through has suggested that when population sizes are high, weak selection for local mitigation can be effective thus predicting a positive correlation between ASC enrichment and Ne. Contra to prediction, we find that ASC enrichment is not correlated with Ne. ASC enrichment, although highly phylogenetically patchy, is, however, more common both in unicellular species and in genes expressed in unicellular modes in multicellular species. By contrast, Ne does positively correlate with TAA enrichment. These results imply that local phenotypic error rates, not local mitigation rates, are consistent with a drift barrier/nearly neutral model.

Retrieving the deleterious mutations before extinction: genome-wide comparison of shared derived mutations in liver cancer and normal population

STUDY PURPOSE

Deleterious mutations would be rapidly purged from natural populations along with the extinction of their carriers. The currently observed mutations in existing species are mostly neutral. The inaccessibility of deleterious mutations impedes the functional studies on how these mutations affect the fitness at individual level.

STUDY DESIGN

The connection between the deleterious genotype and the non-adaptive phenotype could be bridged by sequencing the genome before extinction. Although this approach is no longer feasible for evolutionary biologists, it is feasible for cancer biologists by profiling the mutations in tumour samples which are so deleterious that the carriers hardly live.

RESULTS

By comparing the derived mutation profile between normal populations and patients with liver cancer, we found that the shared mutations, which are highly deleterious, are suppressed to low allele frequencies in normal populations and tissues, but show remarkably high frequency in tumours. The density of shared mutations is negatively correlated with gene conservation and expression levels.

CONCLUSIONS

Deleterious mutations are suppressed in functionally important genes as well as in normal populations. This work deepened our understanding on how natural selection act on deleterious mutations by analogising the cancer evolution to species evolution, which are essentially the same molecular process but at different time scales.

Synonymous mutations that regulate translation speed might play a non-negligible role in liver cancer development

Background

Synonymous mutations do not change the protein sequences. Automatically, they have been regarded as neutral events and are ignored in the mutation-based cancer studies. However, synonymous mutations will change the codon optimality, resulting in altered translational velocity.

Methods

We fully utilized the transcriptome and translatome of liver cancer and normal tissue from ten patients. We profiled the mutation spectrum and examined the effect of synonymous mutations on translational velocity.

Results

Synonymous mutations that increase the codon optimality significantly enhanced the translational velocity, and were enriched in oncogenes. Meanwhile, synonymous mutations decreasing codon optimality slowed down translation, and were enriched in tumor suppressor genes. These synonymous mutations significantly contributed to the translational changes in tumor samples compared to normal samples.

Conclusions

Synonymous mutations might play a role in liver cancer development by altering codon optimality and translational velocity. Synonymous mutations should no longer be ignored in the genome-wide studies.

Readthrough Errors Purge Deleterious Cryptic Sequences, Facilitating the Birth of Coding Sequences

De novo protein-coding innovations sometimes emerge from ancestrally noncoding DNA, despite the expectation that translating random sequences is overwhelmingly likely to be deleterious. The "preadapting selection" hypothesis claims that emergence is facilitated by prior, low-level translation of noncoding sequences via molecular errors. It predicts that selection on polypeptides translated only in error is strong enough to matter and is strongest when erroneous expression is high. To test this hypothesis, we examined noncoding sequences located downstream of stop codons (i.e., those potentially translated by readthrough errors) in Saccharomyces cerevisiae genes. We identified a class of "fragile" proteins under strong selection to reduce readthrough, which are unlikely substrates for co-option. Among the remainder, sequences showing evidence of readthrough translation, as assessed by ribosome profiling, encoded C-terminal extensions with higher intrinsic structural disorder, supporting the preadapting selection hypothesis. The cryptic sequences beyond the stop codon, rather than spillover effects from the regular C-termini, are primarily responsible for the higher disorder. Results are robust to controlling for the fact that stronger selection also reduces the length of C-terminal extensions. These findings indicate that selection acts on 3' UTRs in Saccharomyces cerevisiae to purge potentially deleterious variants of cryptic polypeptides, acting more strongly in genes that experience more readthrough errors.

Mammalian Alternative Translation Initiation Is Mostly Nonadaptive

Alternative translation initiation (ATLI) refers to the existence of multiple translation initiation sites per gene and is a widespread phenomenon in eukaryotes. ATLI is commonly assumed to be advantageous through creating proteome diversity or regulating protein synthesis. We here propose an alternative hypothesis that ATLI arises primarily from nonadaptive initiation errors presumably due to the limited ability of ribosomes to distinguish sequence motifs truly signaling translation initiation from similar sequences. Our hypothesis, but not the adaptive hypothesis, predicts a series of global patterns of ATLI, all of which are confirmed at the genomic scale by quantitative translation initiation sequencing in multiple human and mouse cell lines and tissues. Similarly, although many codons differing from AUG by one nucleotide can serve as start codons, our analysis suggests that using non-AUG start codons is mostly disadvantageous. These and other findings strongly suggest that ATLI predominantly results from molecular error, requiring a major revision of our understanding of the precision and regulation of translation initiation.

Characterization of Bicistronic Transcription in Budding Yeast

Bicistronic transcripts (operon-like transcripts) have occasionally been reported in eukaryotes, including unicellular yeasts, plants, and humans, despite the fact that they lack trans-splice mechanisms. However, the characteristics of eukaryotic bicistronic transcripts are poorly understood, except for those in nematodes. Here, we describe the genomic, transcriptomic, and ribosome profiling features of bicistronic transcripts in unicellular yeasts. By comparing the expression level of bicistronic transcripts with their monocistronic equivalents, we identify two main categories of bicistronic transcripts: highly and lowly expressed. These two categories exhibit quite different features. First, highly expressed bicistronic transcripts have higher conservation within and between strains and shorter intergenic spacers with higher GC content and less stable secondary structure. Second, genes in highly expressed bicistronic transcripts have lower translation efficiency, with the second gene showing statistically significant lower translation efficiency than the first. Finally, the genes found in these highly expressed bicistronic transcripts tend to be younger, with more recent origins. Together, these results suggest that bicistronic transcripts in yeast are heterogeneous. We further propose that at least some highly expressed bicistronic transcripts appear to play a role in modulating monocistronic translation.IMPORTANCE Operons, where a single mRNA transcript encodes multiple adjacent proteins, are a widespread feature of bacteria and archaea. In contrast, the genes of eukaryotes are generally considered monocistronic. However, a number of studies have revealed the presence of bicistronic transcripts in eukaryotes, including humans. The basic features of these transcripts are largely unknown in eukaryotes, especially in organisms lacking trans-splice mechanisms. Our analyses characterize bicistronic transcripts in one such eukaryotic group, yeasts. We show that highly expressed bicistronic transcripts have unusual features compared to lowly expressed bicistronic transcripts, with several features influencing translational modulation.

Specificity of RNA Folding and Its Association with Evolutionarily Adaptive mRNA Secondary Structures

The secondary structure is a fundamental feature of both noncoding and messenger RNAs. However, our understanding of the secondary structure of mRNA, especially that of the coding regions, remains elusive, likely due to translation and the lack of RNA-binding proteins that sustain the consensus structure, such as those that bind to noncoding RNA. Indeed, mRNA has recently been found to adopt diverse alternative structures, the overall functional significance of which remains untested. We hereby approached this problem by estimating the folding specificity, i.e., the probability that a fragment of RNA folds back to the same partner once refolded. We showed that the folding specificity of mRNA is lower than that of noncoding RNA and exhibits moderate evolutionary conservation. Notably, we found that specific rather than alternative folding is likely evolutionarily adaptive since specific folding is frequently associated with functionally important genes or sites within a gene. Additional analysis in combination with ribosome density suggests the ability to modulate ribosome movement as one potential functional advantage provided by specific folding. Our findings revealed a novel facet of the RNA structurome with important functional and evolutionary implications and indicated a potential method for distinguishing the mRNA secondary structures maintained by natural selection from molecular noise.

Mutations in m6A consensus motifs are suppressed in the m6A modified genes in human cancer cells

N⁶-methyladenosine (m⁶A) is the most prevalent type of RNA modification. METTL3 in the methyltransferase complex is the core enzyme responsible for methylation. METTL3 selectively catalyzes the adenosines centered in the RRAC motif. Functional studies established that m⁶A could enhance the translation efficiency (TE) of modified genes by recruiting reader protein YTHDF1 and other initiation factors. We downloaded the m⁶A peaks in HeLa cells from a previous study and defined the m⁶A modified genes and sites. Ancestral mutations in the genic region fixed in the HeLa cell samples were defined using their mRNA-Seq data and the alignment between human and mouse genomes. Furthermore, in the small interfering (si)-METTL3 sample, the calculated TE foldchange of all genes was compared to that in the negative control. The TE of m⁶A genes was globally down-regulated in si-METTL3 versus control compared to the non-m⁶A genes. In m⁶A modified genes, RRAC motif mutations were suppressed compared to mutations in non-motif regions or non-m⁶A genes. Among the m⁶A genes, a fraction RRAC motif mutations negatively correlated with the TE foldchange (si-METTL3 versus control). The TE of m⁶A modified genes was enhanced in HeLa cells. RRAC motif mutations could potentially prevent methylation of adenosines and consequently abolish the enhanced translation. Such mutations in the RRAC motif might be deleterious. Accordingly, we observed lower fractions of mutations in RRAC motifs than in other regions. This prevention of mutations in the RRAC motif could be a strategy adopted by cancer cells to maintain the elevated translation of particular genes.

Mechanisms of RNA N6-Methyladenosine in Hepatocellular Carcinoma: From the Perspectives of Etiology

N6-Methyladenosine (m⁶A) is the most common RNA internal modification in eukaryotic cells. Its regulatory effects at the post-transcriptional level on both messenger RNAs (mRNAs) and noncoding RNAs have been widely studied; these include alternative splicing, stability, translation efficiency, nucleus export, and degradation. m⁶A modification is implicated in a series of physiological and pathological activities, such as embryonic stem cell differentiation, immunoregulation, adipogenesis, and cancer development. Recently, the significance of m⁶A methylation has been identified in both viral hepatitis and non-alcohol fatty liver disease (NAFLD), which are major risk factors in the development of hepatocellular carcinoma (HCC). Given the high incidence and mortality rate of HCC worldwide, it is of great importance to elucidate the mechanisms underlying HCC initiation and progression. m⁶A as an emerging research focus has great potential to facilitate the understanding of HCC, particularly from an etiological perspective. Thus, in this review, we summarize recent progress in understanding m⁶A modification related to viral hepatitis, NAFLD, and HCC, including their mechanisms and clinical applications.

A Census and Categorization Method of Epitranscriptomic Marks

In the past few years, thorough investigation of chemical modifications operated in the cells on ribonucleic acid (RNA) molecules is gaining momentum. This new field of research has been dubbed “epitranscriptomics”, in analogy to best-known epigenomics, to stress the potential of ensembles of RNA modifications to constitute a post-transcriptional regulatory layer of gene expression orchestrated by writer, reader, and eraser RNA-binding proteins (RBPs). In fact, epitranscriptomics aims at identifying and characterizing all functionally relevant changes involving both non-substitutional chemical modifications and editing events made to the transcriptome. Indeed, several types of RNA modifications that impact gene expression have been reported so far in different species of cellular RNAs, including ribosomal RNAs, transfer RNAs, small nuclear RNAs, messenger RNAs, and long non-coding RNAs. Supporting functional relevance of this largely unknown regulatory mechanism, several human diseases have been associated directly to RNA modifications or to RBPs that may play as effectors of epitranscriptomic marks. However, an exhaustive epitranscriptome’s characterization, aimed to systematically classify all RNA modifications and clarify rules, actors, and outcomes of this promising regulatory code, is currently not available, mainly hampered by lack of suitable detecting technologies. This is an unfortunate limitation that, thanks to an unprecedented pace of technological advancements especially in the sequencing technology field, is likely to be overcome soon. Here, we review the current knowledge on epitranscriptomic marks and propose a categorization method based on the reference ribonucleotide and its rounds of modifications (“stages”) until reaching the given modified form. We believe that this classification scheme can be useful to coherently organize the expanding number of discovered RNA modifications.

High Transcriptional Error Rates Vary as a Function of Gene Expression Level

Errors in gene transcription can be costly, and organisms have evolved to prevent their occurrence or mitigate their costs. The simplest interpretation of the drift barrier hypothesis suggests that species with larger population sizes would have lower transcriptional error rates. However, Escherichia coli seems to have a higher transcriptional error rate than species with lower effective population sizes, for example Saccharomyces cerevisiae. This could be explained if selection in E. coli were strong enough to maintain adaptations that mitigate the consequences of transcriptional errors through robustness, on a gene by gene basis, obviating the need for low transcriptional error rates and associated costs of global proofreading. Here, we note that if selection is powerful enough to evolve local robustness, selection should also be powerful enough to locally reduce error rates. We therefore predict that transcriptional error rates will be lower in highly abundant proteins on which selection is strongest. However, we only expect this result when error rates are high enough to significantly impact fitness. As expected, we find such a relationship between expression and transcriptional error rate for non-C→U errors in E. coli (especially G→A), but not in S. cerevisiae. We do not find this pattern for C→U changes in E. coli, presumably because most deamination events occurred during sample preparation, but do for C→U changes in S. cerevisiae, supporting the interpretation that C→U error rates estimated with an improved protocol, and which occur at rates comparable with E. coli non-C→U errors, are biological.

The significance of m6A RNA methylation regulators in predicting the prognosis and clinical course of HBV-related hepatocellular carcinoma

Background

Hepatocarcinogenesis is reportedly correlated with abnormal m6A modifications; however, it is unknown whether m6A RNA methylation regulators facilitate the occurrence of hepatitis B virus (HBV)-related hepatocellular carcinoma (HCC). Thus, we constructed an m6A-related model that may enhance HBV-related HCC prognosis.

Methods

Gene signatures of HNRNPA2B1 and RBM15 were generated by univariate and Lasso Cox regression analyses using the gene set and clinical information from The Cancer Genome Atlas (TCGA) database. High-risk and low-risk groups were confirmed based on the gene signature model. Furthermore, we validated the predictive roles of the two genes for overall survival (OS) in the GSE14520 dataset. The relative expression of 22 paired mRNAs was measured using quantitative real-time polymerase chain reaction (qRT-PCR) analysis to determine whether the two genes had a predictive role in our Guilin cohort.

Results

The differences in OS between the high-risk and low-risk groups were statistically significant in the TCGA (p = 0.003) and GSE14520 (p = 0.045) datasets, but not in the Guilin cohort, owing to differences in clinical information among the three cohorts (mainly the TNM stage and survival state). Stratified analysis of TNM stages showed that the two-gene signature acted as a prognostic indicator of HBV-related HCC patients in the early TNM stage; both TCGA and GSE14520 cohorts showed statistical significance. Moreover, multivariate Cox regression analysis indicated that the two-gene signature was an independent factor for predicting prognosis (HR = 1.087, 95% CI: 1.007–1.172). Correlation analysis between the gene signature and clinical features revealed that the risk stratification was significantly correlated with grade and survival state. Finally, Gene Set Enrichment Analysis (GSEA) revealed that the KEGG pathways associated with the cell cycle, DNA replication, the spliceosome, repair, and metabolism-related processes were all significantly enriched in the high-risk group. Among the enriched genes, the expression levels of the replication protein RPA1 and the pre-mRNA splicing factor SF3B1 were significantly upregulated in the high-risk group. These results might help in elucidating the underlying molecular mechanisms of HBV-related HCC.

Conclusions

Our data may provide new predictive signatures and potential therapeutic targets to identify and treat HBV-related HCC patients in the early disease stage.

Non-canonical RNA-DNA differences and other human genomic features are enriched within very short tandem repeats

Very short tandem repeats bear substantial genetic, evolutional, and pathological significance in genome analyses. Here, we compiled a census of tandem mono-nucleotide/di-nucleotide/tri-nucleotide repeats (MNRs/DNRs/TNRs) in GRCh38, which we term "polytracts" in general. Of the human genome, 144.4 million nucleotides (4.7%) are occupied by polytracts, and 0.47 million single nucleotides are identified as polytract hinges, i.e., break-points of tandem polytracts. Preliminary exploration of the census suggested polytract hinge sites and boundaries of AAC polytracts may bear a higher mapping error rate than other polytract regions. Further, we revealed landscapes of polytract enrichment with respect to nearly a hundred genomic features. We found MNRs, DNRs, and TNRs displayed noticeable difference in terms of locational enrichment for miscellaneous genomic features, especially RNA editing events. Non-canonical and C-to-U RNA-editing events are enriched inside and/or adjacent to MNRs, while all categories of RNA-editing events are under-represented in DNRs. A-to-I RNA-editing events are generally under-represented in polytracts. The selective enrichment of non-canonical RNA-editing events within MNR adjacency provides a negative evidence against their authenticity. To enable similar locational enrichment analyses in relation to polytracts, we developed a software Polytrap which can handle 11 reference genomes. Additionally, we compiled polytracts of four model organisms into a Track Hub which can be integrated into USCS Genome Browser as an official track for convenient visualization of polytracts.

The Epigenetic Drug Discovery Landscape for Metabolic-associated Fatty Liver Disease

Despite decades of research, effective therapies for metabolic (dysfunction)-associated fatty liver disease (MAFLD) are lacking. An increasing body of evidence suggests that epigenetic dysregulation is frequent in MAFLD, and orchestrates many aspects of its development and progression. Furthermore, the high plasticity of epigenetic modifications in response to environmental cues renders epigenetics a novel area for therapeutic drug discovery. Over recent years, several epigenetics-based drugs and diagnostic biomarkers have entered clinical development and/or obtained regulatory approval. Here, we review recent advances in our understanding of epigenetic regulation and programming during MAFLD, including DNA methylation, histone modifications, chromatin remodelling, transcriptional control, and noncoding (nc)RNAs. We also discuss the potential translational implications and challenges of epigenetics in the context of MAFLD.