3'UTR Diversity: Expanding Repertoire of RNA Alterations in Human mRNAs

Reovirus infection induces transcriptome-wide unique A-to-I editing changes in the murine fibroblasts

The conversion of Adenosine (A) to Inosine (I), by Adenosine Deaminases Acting on RNA or ADARs, is an essential post-transcriptional modification that contributes to proteome diversity and regulation in metazoans including humans. In addition to its transcriptome-regulating role, ADARs also play a major part in immune response to viral infection, where an interferon response activates interferon-stimulated genes, such as ADARp150, in turn dynamically regulating host-virus interactions. A previous report has shown that infection from reoviruses, despite strong activation of ADARp150, does not influence the editing of some of the major known editing targets, while likely editing others, suggesting a potentially nuanced editing pattern that may depend on different factors. However, the results were based on a handful of selected editing sites and did not cover the entire transcriptome. Thus, to determine whether and how reovirus infection specifically affects host ADAR editing patterns, we analyzed a publicly available deep-sequenced RNA-seq dataset, from murine fibroblasts infected with wild-type and mutant reovirus strains that allowed us to examine changes in editing patterns on a transcriptome-wide scale. To the best of our knowledge, this is the first transcriptome-wide report on host editing changes after reovirus infection. Our results demonstrate that reovirus infection induces unique nuanced editing changes in the host, including introducing sites uniquely edited in infected samples. Genes with edited sites are overrepresented in pathways related to immune regulation, cellular signaling, metabolism, and growth. Moreover, a shift in editing targets has also been observed, where the same genes are edited in infection and control conditions but at different sites, or where the editing rate is increased for some and decreased for other differential targets, supporting the hypothesis of dynamic and condition-specific editing by ADARs.

Revealing miRNAs patterns by employing matrix representations and energy analysis

MicroRNAs (miRNAs) are small, non-coding RNA molecules that regulate gene expression. Despite their relatively short length (about 21 nucleotides), they can regulate thousands of transcripts within a cell. Due to their low complementarity to targets, studying their activity and binding region preferences (3'UTR, 5'UTR, or CDS) is challenging. In this paper, we analyzed a set of human miRNAs to uncover their general patterns. We began with a sequence logo to verify conservation at specific positions. To discover long-range correlations, we employed chaos game representation (CGR) and genomatrix, methods that enable both graphical and analytical analysis of sequence sets and are well-established in bioinformatics. Our results showed that miRNAs exhibit strongly non-random and characteristic patterns. To incorporate physicochemical properties into the analysis, we applied the electron-ion interaction potential (EIIP) parameter. An important part of our study was to validate the division of miRNAs into two parts-seed and puzzle. The seed region is responsible for target binding, while the puzzle region likely interacts with the RISC complex. We estimated duplex binding energy within the 3'UTR, 5'UTR, and CDS regions using the miRanda tool. Based on the median energy distribution, we divided the miRNAs into two subsets, reflecting different patterns in chaos game representation. Interestingly, one subset displayed significant similarity to conserved and highly confidential miRNAs. Our results confirm the low complementarity of miRNA/mRNA interactions and support the functional division of miRNA structure. Additionally, we present findings related to the localization of transcript target sites, which form the basis for further analyses.

Variability in SOD1-associated amyotrophic lateral sclerosis: geographic patterns, clinical heterogeneity, molecular alterations, and therapeutic implications

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease characterized by progressive loss of motor neurons, resulting in global health burden and limited post-diagnosis life expectancy. Although primarily sporadic, familial ALS (fALS) cases suggest a genetic basis. This review focuses on SOD1, the first gene found to be associated with fALS, which has been more recently confirmed by genome sequencing. While informative, databases such as ALSoD and STRENGTH exhibit regional biases. Through a systematic global examination of SOD1 mutations from 1993 to 2023, we found different geographic distributions and clinical presentations. Even though different SOD1 variants are expressed at different protein levels and have different half-lives and dismutase activities, these alterations lead to loss of function that is not consistently correlated with disease severity. Gain of function of toxic aggregates of SOD1 resulting from mutated SOD1 has emerged as one of the key contributors to ALS. Therapeutic interventions specifically targeting toxic gain of function of mutant SOD1, including RNA interference and antibodies, show promise, but a cure remains elusive. This review provides a comprehensive perspective on SOD1-associated ALS and describes molecular features and the complex genetic landscape of SOD1, highlighting its importance in determining diverse clinical manifestations observed in ALS patients and emphasizing the need for personalized therapeutic strategies.

Brief guide to RNA sequencing analysis for nonexperts in bioinformatics

Transcriptome analysis is widely used for current biological research but remains challenging for many experimental scientists. Here, we present a brief but broad guideline for transcriptome analysis, focusing on RNA sequencing, by providing the list of publicly available datasets, tools, and R packages for practical transcriptome analysis. This work will be useful for biologists to perform key transcriptomic analysis with minimum expertise in bioinformatics.

Meta-analysis of the transcriptome identifies aberrant RNA processing as common feature of aging in multiple species

Aging is accompanied by the gradual deregulation of the transcriptome. However, whether age-dependent changes in the transcriptome are evolutionarily conserved or diverged remains largely unexplored. Here, we performed a meta-analysis examining the age-dependent changes in the transcriptome using publicly available datasets of 11 representative metazoans, ranging from Caenorhabditis elegans to humans. To identify the transcriptomic changes associated with aging, we analyzed various aspects of the transcriptome, including genome composition, RNA processing, and functional consequences. The use of introns and novel splice sites tended to increase with age, particularly in the brain. In addition, our analysis suggests that the age-dependent accumulation of premature termination codon-containing transcripts is a common feature of aging across multiple animal species. Using C. elegans as a test model, we showed that several splicing factors that are evolutionarily conserved and age-dependently downregulated were required to maintain a normal lifespan. Thus, aberrant RNA processing appears to be associated with aging and a short lifespan in various species.

Large-scale RNA-seq mining reveals ciclopirox triggers TDP-43 cryptic exons

Nuclear clearance and cytoplasmic aggregation of TDP-43 in neurons, initially identified in ALS-FTD, are hallmark pathological features observed across a spectrum of neurodegenerative diseases. We previously found that TDP-43 loss-of-function leads to the transcriptome-wide inclusion of deleterious cryptic exons in brains and biofluids post-mortem as well as during the presymptomatic stage of ALS-FTD, but upstream mechanisms that lead to TDP-43 dysregulation remain unclear. Here, we developed a web-based resource (SnapMine) to determine the levels of TDP-43 cryptic exon inclusion across hundreds of thousands of publicly available RNA sequencing datasets. We established cryptic exon inclusion across a variety of human cells and tissues to provide ground truth references for future studies on TDP-43 dysregulation. We then explored studies that were entirely unrelated to TDP-43 or neurodegeneration and found that ciclopirox olamine (CPX), an FDA-approved antifungal, can trigger the inclusion of TDP-43-associated cryptic exons in a variety of mouse and human primary cells. CPX induction of cryptic exon occurs via heavy metal toxicity and oxidative stress, suggesting that similar vulnerabilities could play a role in neurodegeneration. Our work demonstrates how diverse datasets can be linked through common biological features and underscores that public archives of sequencing data represent a vastly underutilized resource with tremendous potential for uncovering novel insights into complex biological mechanisms and diseases.

Utilizing non-coding RNA-mediated regulation of ATP binding cassette (ABC) transporters to overcome multidrug resistance to cancer chemotherapy

Multidrug resistance (MDR) is one of the primary factors that produces treatment failure in patients receiving cancer chemotherapy. MDR is a complex multifactorial phenomenon, characterized by a decrease or abrogation of the efficacy of a wide spectrum of anticancer drugs that are structurally and mechanistically distinct. The overexpression of the ATP-binding cassette (ABC) transporters, notably ABCG2 and ABCB1, are one of the primary mediators of MDR in cancer cells, which promotes the efflux of certain chemotherapeutic drugs from cancer cells, thereby decreasing or abolishing their therapeutic efficacy. A number of studies have suggested that non-coding RNAs (ncRNAs), particularly microRNAs (miRNAs), long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs), play a pivotal role in mediating the upregulation of ABC transporters in certain MDR cancer cells. This review will provide updated information about the induction of ABC transporters due to the aberrant regulation of ncRNAs in cancer cells. We will also discuss the measurement and biological profile of circulating ncRNAs in various body fluids as potential biomarkers for predicting the response of cancer patients to chemotherapy. Sequence variations, such as alternative polyadenylation of mRNA and single nucleotide polymorphism (SNPs) at miRNA target sites, which may indicate the interaction of miRNA-mediated gene regulation with genetic variations to modulate the MDR phenotype, will be reviewed. Finally, we will highlight novel strategies that could be used to modulate ncRNAs and circumvent ABC transporter-mediated MDR.

Hidden features of NAD-RNA epitranscriptome in Drosophila life cycle

Nicotinamide adenine dinucleotide (NAD), a nucleotide-containing metabolite, can be incorporated into the RNA 5'-terminus to result in NAD-capped RNA (NAD-RNA). Since NAD has been heightened as one of the most essential metabolites in cells, its linkage to RNA represents a critical but poorly studied modification at the epitranscriptomic level. Here, we design a highly sensitive method, DO-seq, to capture NAD-RNAs. Using Drosophila, we identify thousands of previously unexplored NAD-RNAs and their dynamics in the fly life cycle, from embryo to adult. We show the evidence that chromosomal clustering might be the structural basis by which co-expression can couple with NAD capping on physically and functionally linked genes. Furthermore, we note that NAD capping of cuticle genes inversely correlates with their gene expression. Combined, we propose NAD-RNA epitranscriptome as a hidden layer of regulation that underlies biological processes. DO-seq empowers the identification of NAD-capped RNAs, facilitating functional investigation into this modification.

Synaptopodin-2 Isoforms Have Specific Binding Partners and Display Distinct, Muscle Cell Type-Specific Expression Patterns

Synaptopodin-2 (SYNPO2) is a protein associated with the Z-disc in striated muscle cells. It interacts with α-actinin and filamin C, playing a role in Z-disc maintenance under stress by chaperone-assisted selective autophagy (CASA). In smooth muscle cells, SYNPO2 is a component of dense bodies. Furthermore, it has been proposed to play a role in tumor cell proliferation and metastasis in many different kinds of cancers. Alternative transcription start sites and alternative splicing predict the expression of six putative SYNPO2 isoforms differing by extended amino- and/or carboxy-termini. Our analyses at mRNA and protein levels revealed differential expression of SYNPO2 isoforms in cardiac, skeletal and smooth muscle cells. We identified synemin, an intermediate filament protein, as a novel binding partner of the PDZ-domain in the amino-terminal extension of the isoforms mainly expressed in cardiac and smooth muscle cells, and demonstrated colocalization of SYNPO2 and synemin in both cell types. A carboxy-terminal extension, mainly expressed in smooth muscle cells, is sufficient for association with dense bodies and interacts with α-actinin. SYNPO2 therefore represents an additional and novel link between intermediate filaments and the Z-discs in cardiomyocytes and dense bodies in smooth muscle cells, respectively. In pathological skeletal muscle samples, we identified SYNPO2 in the central and intermediate zones of target fibers of patients with neurogenic muscular atrophy, and in nemaline bodies. Our findings help to understand distinct functions of individual SYNPO2 isoforms in different muscle tissues, but also in tumor pathology.

Synonymous codon usage regulates translation initiation

Nonoptimal synonymous codons repress gene expression, but the underlying mechanisms are poorly understood. We and others have previously shown that nonoptimal codons slow translation elongation speeds and thereby trigger messenger RNA (mRNA) degradation. Nevertheless, transcript levels are often insufficient to explain protein levels, suggesting additional mechanisms by which codon usage regulates gene expression. Using reporters in human and Drosophila cells, we find that transcript levels account for less than half of the variation in protein abundance due to codon usage. This discrepancy is explained by translational differences whereby nonoptimal codons repress translation initiation. Nonoptimal transcripts are also less bound by the translation initiation factors eIF4E and eIF4G1, providing a mechanistic explanation for their reduced initiation rates. Importantly, translational repression can occur without mRNA decay and deadenylation, and it does not depend on the known nonoptimality sensor, CNOT3. Our results reveal a potent mechanism of regulation by codon usage where nonoptimal codons repress further rounds of translation.

Functional Consequences of Shifting Transcript Boundaries in Glucose Starvation

Glucose is a major source of carbon and essential for the survival of many organisms, ranging from yeast to human. A sudden 60-fold reduction of glucose in exponentially growing fission yeast induces transcriptome-wide changes in gene expression. This regulation is multilayered, and the boundaries of transcripts are known to vary, with functional consequences at the protein level. By combining direct RNA sequencing with 5'-CAGE and short-read sequencing, we accurately defined the 5'- and 3'-ends of transcripts that are both poly(A) tailed and 5'-capped in glucose starvation, followed by proteome analysis. Our results confirm previous experimentally validated loci with alternative isoforms and reveal several transcriptome-wide patterns. First, we show that sense-antisense gene pairs are more strongly anticorrelated when a time lag is taken into account. Second, we show that the glucose starvation response initially elicits a shortening of 3'-UTRs and poly(A) tails, followed by a shortening of the 5'-UTRs at later time points. These result in domain gains and losses in proteins involved in the stress response. Finally, the relatively poor overlap both between differentially expressed genes (DEGs), differential transcript usage events (DTUs), and differentially detected proteins (DDPs) highlight the need for further study on post-transcriptional regulation mechanisms in glucose starvation.

Functional mutation, splice, distribution, and divergence analysis of impactful genes associated with heart failure and other cardiovascular diseases

Cardiovascular disease (CVD) is caused by a multitude of complex and largely heritable conditions. Identifying key genes and understanding their susceptibility to CVD in the human genome can assist in early diagnosis and personalized treatment of the relevant patients. Heart failure (HF) is among those CVD phenotypes that has a high rate of mortality. In this study, we investigated genes primarily associated with HF and other CVDs. Achieving the goals of this study, we built a cohort of thirty-five consented patients, and sequenced their serum-based samples. We have generated and processed whole genome sequence (WGS) data, and performed functional mutation, splice, variant distribution, and divergence analysis to understand the relationships between each mutation type and its impact. Our variant and prevalence analysis found FLNA, CST3, LGALS3, and HBA1 linked to many enrichment pathways. Functional mutation analysis uncovered ACE, MME, LGALS3, NR3C2, PIK3C2A, CALD1, TEK, and TRPV1 to be notable and potentially significant genes. We discovered intron, 5' Flank, 3' UTR, and 3' Flank mutations to be the most common among HF and other CVD genes. Missense mutations were less common among HF and other CVD genes but had more of a functional impact. We reported HBA1, FADD, NPPC, ADRB2, ADBR1, MYH6, and PLN to be consequential based on our divergence analysis.

METTL3 enhances dentinogenesis differentiation of dental pulp stem cells via increasing GDF6 and STC1 mRNA stability

BACKGROUND

The dentinogenesis differentiation of dental pulp stem cells (DPSCs) is controlled by the spatio-temporal expression of differentiation related genes. RNA N6-methyladenosine (m⁶A) methylation, one of the most abundant internal epigenetic modification in mRNA, influences various events in RNA processing, stem cell pluripotency and differentiation. Methyltransferase like 3 (METTL3), one of the essential regulators, involves in the process of dentin formation and root development, while mechanism of METTL3-mediated RNA m⁶A methylation in DPSC dentinogenesis differentiation is still unclear.

METHODS

Immunofluorescence staining and MeRIP-seq were performed to establish m⁶A modification profile in dentinogenesis differentiation. Lentivirus were used to knockdown or overexpression of METTL3. The dentinogenesis differentiation was analyzed by alkaline phosphatase, alizarin red staining and real time RT-PCR. RNA stability assay was determined by actinomycin D. A direct pulp capping model was established with rat molars to reveal the role of METTL3 in tertiary dentin formation.

RESULTS

Dynamic characteristics of RNA m⁶A methylation in dentinogenesis differentiation were demonstrated by MeRIP-seq. Methyltransferases (METTL3 and METTL14) and demethylases (FTO and ALKBH5) were gradually up-regulated during dentinogenesis process. Methyltransferase METTL3 was selected for further study. Knockdown of METTL3 impaired the DPSCs dentinogenesis differentiation, and overexpression of METTL3 promoted the differentiation. METTL3-mediated m⁶A regulated the mRNA stabiliy of GDF6 and STC1. Furthermore, overexpression of METTL3 promoted tertiary dentin formation in direct pulp capping model.

CONCLUSION

The modification of m⁶A showed dynamic characteristics during DPSCs dentinogenesis differentiation. METTL3-mediated m⁶A regulated in dentinogenesis differentiation through affecting the mRNA stability of GDF6 and STC1. METTL3 overexpression promoted tertiary dentin formation in vitro, suggesting its promising application in vital pulp therapy (VPT).