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

Coupling of alternative splicing and alternative polyadenylation

RNA splicing and 3'-cleavage and polyadenylation (CPA) are essential processes for the maturation of RNA. There have been extensive independent studies of these regulated processing events, including alternative splicing (AS) and alternative polyadenylation (APA). However, growing evidence suggests potential crosstalk between splicing and 3'-end processing in regulating AS or APA. Here, we first provide a brief overview of the molecular machines involved in splicing and 3'-end processing events, and then review recent studies on the functions and mechanisms of the crosstalk between the two processes. On one hand, 3'-end processing can affect splicing, as 3'-end processing factors and CPA-generated polyA tail promote the splicing of the last intron. Beyond that, 3'-end processing factors can also influence the splicing of internal and terminal exons. Those 3'-end processing factors can also interact with different RNA-binding proteins (RBPs) to exert their effects on AS. The length of 3' untranslated region (3' UTR) can affect the splicing of upstream exons. On the other hand, splicing and CPA may compete within introns in generating different products. Furthermore, splicing within the 3' UTR is a significant factor contributing to 3' UTR diversity. Splicing also influences 3'-end processing through the actions of certain splicing factors. Interestingly, some classical RBPs play dual roles in both splicing and 3'-end processing. Finally, we discuss how long-read sequencing technologies aid in understanding the coordination of AS-APA events and envision that these findings may potentially promote the development of new strategies for disease diagnosis and treatment.

The low-density lipoprotein receptor: Emerging post-transcriptional regulatory mechanisms

Cholesterol is a vital component of cellular membranes and is an essential molecule in mammalian physiology. Yet dysregulation of hepatic cholesterol metabolism and an increase in plasma cholesterol is linked to development of atherosclerotic cardiovascular disease. Maintaining tight regulation of cholesterol homeostasis is therefore essential, elegantly highlighted by the control of hepatic low-density lipoprotein receptor (LDLR) abundance and associated lipoprotein clearance. The LDLR was discovered in the 1970's in the seminal work of Brown and Goldstein. This was followed by the development of statins, which promote hepatic clearance of LDL via the LDLR pathway. The discovery two decades ago of Proprotein Convertase Subtilisin-Kexin Type 9 (PCSK9), a secreted protein that binds to the LDLR ectodomain and promotes its degradation, and the clinical development of PCSK9 inhibitors has ushered an effort to uncover additional mechanisms that govern the function and abundance of the LDLR. In recent years this has led to the identification of novel post-transcriptional and post-translational mechanisms that govern the LDLR. This review focuses on these emerging regulatory mechanisms and specifically discusses: (1) Regulation of the LDLR mRNA by RNA-binding proteins and microRNAs, (2) Ubiquitin-dependent degradation of the LDLR protein by the E3 ubiquitin ligases inducible degrader of the LDLR (IDOL) and GOLIATH (RNF130), (3) Control of the LDLR pathway by the asialoglycoprotein receptor 1 (ASGR1), and (4) The role of LDLR ectodomain shedding mediated by membrane-type 1 matrix metalloprotease (MT1-MMP), Bone morphogenetic protein 1 (BMP1), and γ-secretase. Understanding the contribution of these emerging mechanisms to regulation of the LDLR is important for the development of novel LDLR-focused lipid-lowering strategies.

Transcriptome analysis of two isolates of the tomato pathogen Cladosporium fulvum, uncovers genome-wide patterns of alternative splicing during a host infection cycle

Alternative splicing (AS) is a key element of eukaryotic gene expression that increases transcript and proteome diversity in cells, thereby altering their responses to external stimuli and stresses. While AS has been intensively researched in plants and animals, its frequency, conservation, and putative impact on virulence, are relatively still understudied in plant pathogenic fungi. Here, we profiled the AS events occurring in genes of Cladosporium fulvum isolates Race 5 and Race 4, during nearly a complete compatible infection cycle on their tomato host. Our studies revealed extensive heterogeneity in the transcript isoforms assembled from different isolates, infections, and infection timepoints, as over 80% of the transcript isoforms were singletons that were detected in only a single sample. Despite that, nearly 40% of the protein-coding genes in each isolate were predicted to be recurrently AS across the disparate infection timepoints, infections, and the two isolates. Of these, 37.5% were common to both isolates and 59% resulted in multiple protein isoforms, thereby putatively increasing proteome diversity in the pathogen by 31% during infections. An enrichment analysis showed that AS mostly affected genes likely to be involved in the transport of nutrients, regulation of gene expression, and monooxygenase activity, suggesting a role for AS in finetuning adaptation of C. fulvum on its tomato host during infections. Tracing the location of the AS genes on the fungal chromosomes showed that they were mostly located in repeat-rich regions of the core chromosomes, indicating a causal connection between gene location on the genome and propensity to AS. Finally, multiple cases of differential isoform usage in AS genes of C. fulvum were identified, suggesting that modulation of AS at different infection stages may be another way by which pathogens refine infections on their hosts.

Generation of Transcript Length Variants and Reprogramming of mRNA Splicing During Atherosclerosis Progression in ApoE-Deficient Mice

Background. Variant 3'UTRs provide mRNAs with different binding sites for miRNAs or RNA-binding proteins (RBPs) allowing the establishment of new regulatory environments. Regulation of 3'UTR length impacts on the control of gene expression by regulating accessibility of miRNAs or RBPs to homologous sequences in mRNAs. Objective. Studying the dynamics of mRNA length variations in atherosclerosis (ATS) progression and reversion in ApoE-deficient mice exposed to a high-fat diet and treated with an αCD40-specific siRNA or with a sequence-scrambled siRNA as control. Methods. We gathered microarray mRNA expression data from the aortas of mice after 2 or 16 weeks of treatments, and used these data in a Bioinformatics analysis. Results. Here, we report the lengthening of the 5'UTR/3'UTRs and the shortening of the CDS in downregulated mRNAs during ATS progression. Furthermore, treatment with the αCD40-specific siRNA resulted in the partial reversion of the 3'UTR lengthening. Exon analysis showed that these length variations were actually due to changes in the number of exons embedded in mRNAs, and the further examination of transcripts co-expressed at weeks 2 and 16 in mice treated with the control siRNA revealed a process of mRNA isoform switching in which transcript variants differed in the patterns of alternative splicing or activated latent/cryptic splice sites. Conclusion. We document length variations in the 5'UTR/3'UTR and CDS of mRNAs downregulated during atherosclerosis progression and suggest a role for mRNA splicing reprogramming and transcript isoform switching in the generation of disease-related mRNA sequence diversity and variability.

Site-blocking antisense oligonucleotides as a mechanism to fine-tune MeCP2 expression

Rett syndrome (RTT) is a neurodevelopmental disorder caused by loss-of-function mutations in the methyl-CpG-binding protein 2 (MECP2) gene. Despite its severe phenotypes, studies in mouse models suggest that restoring MeCP2 levels can reverse RTT symptomology. Nevertheless, traditional gene therapy approaches are hindered by MeCP2's narrow therapeutic window, complicating the safe delivery of viral constructs without overshooting the threshold for toxicity. The 3' untranslated region (3' UTR) plays a key role in gene regulation, where factors like miRNAs bind to pre-mRNA and fine-tune expression. Given that each miRNA's contribution is modest, blocking miRNA binding may represent a potential therapeutic strategy for diseases with high dosage sensitivity, like RTT. Here, we present a series of site-blocking antisense oligonucleotides (sbASOs) designed to outcompete repressive miRNA binding at the MECP2 3' UTR. This strategy aims to increase MeCP2 levels in patients with missense or late-truncating mutations, where the hypomorphic nature of the protein can be offset by enhanced abundance. Our results demonstrate that sbASOs can elevate MeCP2 levels in a dose-dependent manner in SH-SY5Y and patient fibroblast cell lines, plateauing at levels projected to be safe. Confirming in vivo functionality, sbASO administration in wild-type mice led to significant Mecp2 upregulation and the emergence of phenotypes associated with Mecp2 overexpression. In a T158M neural stem cell model of RTT, sbASO treatment significantly increased MeCP2 expression and levels of the downstream effector protein brain-derived neurotrophic factor (BDNF). These findings highlight the potential of sbASO-based therapies for MeCP2-related disorders and advocate for their continued development.

The NEDD4-binding protein N4BP1 degrades mRNA substrates through the coding sequence independent of nonsense-mediated decay

3'UTRs are recognized for their role in regulating mRNA turnover while the turnover of a specific group of mRNAs mediated by coding sequences (CDSs) remains poorly understood. N4BP1 is a critical inflammatory regulator in vivo with a molecular mechanism that is not yet clearly defined. Our study reveals that N4BP1 efficiently degrades its mRNA targets via CDS rather than the 3'-UTR. This CDS-dependent mRNA turnover mechanism appears to be a general feature of N4BP1, as evidenced by testing multiple mRNA substrates, such as Fos-C, Fos-B, Jun-B, and C-X-C motif chemokine ligand 1. Detailed mapping of the motif identified a crucial 33-nt (289-322) sequence near the 5'-end of Fos-C-CDS, where the presence of polyC is necessary for N4BP1-mediated degradation. Functional studies involving domain deletion and point mutations showed that both the K homology and N4BP1, YacP-like nuclease domains are essential for N4BP1 to restrict mRNA substrates. The function of N4BP1 in mRNA turnover is not dependent on nonsense-mediated decay as it efficiently restricts mRNA substrates even in cells deficient in UPF1, UPF3A, and UPF3B. Additionally, the function of N4BP1 is not reliant on LUC7L3 despite its known association with this protein. Our findings suggest that N4BP1 acts as an endoribonuclease to degrade mRNA substrates primarily through CDSs containing a C-rich motif.

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.

miR-124 coordinates metabolic regulators acting at early stages of human neurogenesis

Metabolic dysregulation of neurons is associated with diverse human brain disorders. Metabolic reprogramming occurs during neuronal differentiation, but it is not fully understood which molecules regulate metabolic changes at the early stages of neurogenesis. In this study, we report that miR-124 is a driver of metabolic change at the initiating stage of human neurogenesis. Proteome analysis has shown the oxidative phosphorylation pathway to be the most significantly altered among the differentially expressed proteins (DEPs) in the immature neurons after the knockdown of miR-124. In agreement with these proteomics results, miR-124-depleted neurons display mitochondrial dysfunctions, such as decreased mitochondrial membrane potential and cellular respiration. Moreover, morphological analyses of mitochondria in early differentiated neurons after miR-124 knockdown result in smaller and less mature shapes. Lastly, we show the potential of identified DEPs as novel metabolic regulators in early neuronal development by validating the effects of GSTK1 on cellular respiration. GSTK1, which is upregulated most significantly in miR-124 knockdown neurons, reduces the oxygen consumption rate of neural cells. Collectively, our data highlight the roles of miR-124 in coordinating metabolic maturation at the early stages of neurogenesis and provide insights into potential metabolic regulators associated with human brain disorders characterized by metabolic dysfunctions.

A Comprehensive Analysis of CSN1S2 I and II Transcripts Reveals Significant Genetic Diversity and Allele-Specific Exon Skipping in Ragusana and Amiatina Donkeys

The αs2-casein is a phosphoprotein secreted in the milk of most mammals, and it is the most hydrophilic of all caseins. Contrary to genes found in ruminants, in donkeys two different encoding genes for donkey αs2-casein (CSN1S2 I and CSN1S2 II) have been identified. However, unlike in ruminants, the variability at these loci has not been characterized in detail in donkeys until now. In this study, we analyze the transcript profile of the donkey CSN1S2 I and CSN1S2 II genes, and we identify and describe the variability of these loci in the Ragusana and Amiatina breeds reared in Italy. The analysis of the CSN1S2 I Reverse Transcriptase-Polymerase Chain Reaction (RT-PCR) products and subsequent sequencing showed, in addition to correctly spliced mRNA, seven other minor mRNAs resulting from differential splicing events involving, in various combinations, entire exons (4, 5, 6, and 11), parts of exons (5' or 3' end of exon 17), or the recognition of intronic sequences as an exon (exon 12'). Similarly, the transcription analysis of the CSN1S2 II gene revealed a remarkable variability in splicing events, mainly concerning the alternative insertion of an extra exon 7 (named 7'); the first 33 bp of exon 13; or the alternative skipping of exons 9, 10, 11, 12, and 15, and their combinations. At the mRNA level for CSN1S2 I, seven SNPs were observed, five of which led to amino acid changes: p.T73>A, p.I109>V, p.I130>V, p.I146>T, and p.D217>Y. Similarly, nine SNPs were observed at the CSN1S2 II locus, seven of which are non-synonymous: p.L63>F, p.H70>Q, p.D90>N, p.129A>T, p.H131>Y, p.E144>G, and p.F157>S. In addition, the DNA sequencing of exon 17 and flanking introns of the CSN1S2 I gene revealed a G>A transition at the splice acceptor site of CSN1S2 I exon 17 (FM946022.1:c.375-1G>A), resulting in an allele-specific skipping of the first 15 nucleotides of this exon, which encode the peptide 176NKINQ180, and the recognition of an in-frame cryptic splicing acceptor site: arAACAAAATCAACCAG. A genotyping method based on restriction fragment length polymorphism (XbaI PCR-RFLP) was set up for this SNP. In the total population studied (105 Ragusana and 14 Amiatina donkeys), the A allele had a frequency of 0.2437 with no evidence of deviation from the Hardy-Weinberg equilibrium. This study adds new knowledge regarding the genetic variability of αs2-caseins in donkeys and may contribute significantly to the genetic improvement of milk production for this species.

Dynamic landscape of m6A modifications and related post-transcriptional events in muscle-invasive bladder cancer

BACKGROUND

Muscle-invasive bladder carcinoma (MIBC) is a serious and more advanced stage of bladder carcinoma. N6-Methyladenosine (m6A) is a dynamic and reversible modifications that primarily affects RNA stability and alternative splicing. The dysregulation of m6A in MIBC can be potential target for clinical interventions, but there have been limited studies on m6A modifications in MIBC and their associations with post-transcriptional regulatory processes.

METHODS

Paired tumor and adjacent-normal tissues were obtained from three patients with MIBC following radical cystectomy. The additional paired tissues for validation were obtained from patients underwent transurethral resection. Utilizing Nanopore direct-RNA sequencing, we characterized the m6A RNA methylation landscape in MIBC, with a focus on identifying post-transcriptional events potentially affected by changes in m6A sites. This included an examination of differential transcript usage, polyadenylation signal sites, and variations in poly(A) tail length, providing insights into the broader impact of m6A alterations on RNA processing in MIBC.

RESULTS

The prognostic-related m6A genes and m6A-risk model constructed by machine learning enables the stratification of high and low-risk patients with precision. A novel m6A modification site in the 3' untranslated region (3'UTR) of IGLL5 gene were identified, characterized by a lower m6A methylation ratio, elongated poly(A) tails, and a notable bias in transcript usage. Furthermore, we discovered two particular transcripts, VWA1-203 and CEBPB-201. VWA1-203 displayed diminished m6A methylation levels, a truncated 3'UTR, and an elongated poly(A) tail, whereas CEBPB-201 showed opposite trends, highlighting the complex interplay between m6A modifications and RNA processing. Source code was provided on GitHub ( https://github.com/lelelililele/Nanopore-m6A-analysis ).

CONCLUSIONS

The state-of-the-art Nanopore direct-RNA sequencing and machine learning techniques enables comprehensive identification of m6A modification and provided insights into the potential post-transcriptional regulation mechanisms on the development and progression in MIBC.

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.

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