Comprehensive mapping of the alternative polyadenylation site usage and its dynamics at single cell resolution

Alternative polyadenylation (APA) plays an important role in post-transcriptional gene regulation such as transcript stability and translation efficiency. However, our knowledge about APA dynamics at single cell level is largely unexplored. Here we developed single cell polyadenylation sequencing (scPolyA-seq), a strand-specific approach for sequencing 3’ end of transcripts, to investigate the landscape of APA at single cell level. By analyzing several cell lines, we found many genes using multiple polyA sites in bulk data are prone to use only one polyA site in each single cell. Interestingly, cell cycle was significantly enriched in genes showing high variation of polyA site usages. We further identified 414 genes showing polyA site usage switch after cell synchronization. Genes showing cell cycle associated polyA site usage switch were grouped into 6 clusters, with cell phase specific functional categories enriched in each cluster. Furthermore, scPolyA-seq could facilitate study of APA in various biological processes.

Investigating the structural changes due to adenosine methylation of the Kaposi's sarcoma-associated herpes virus ORF50 transcript

Kaposi's sarcoma-associated herpes virus (KSHV) is a human oncovirus. KSHV relies on manipulating the host cell N6- methyl adenosine (m6A) RNA modification pathway. Methylation within a RNA stem loop of the open reading frame 50 (ORF50) enhances transcript stability via the recruitment of the m6A reader, SND1. In this contribution we explore the energy landscapes of the unmethylated and methylated RNA stem loops of ORF50 to investigate the effect of methylation on the structure of the stem loop. We observe a significant shift upon methylation between an open and closed configuration of the top of the stem loop. In the unmethylated stem loop the closed configuration is much lower in energy, and, as a result, exhibits higher occupancy.

Modulation of AggR levels reveals features of virulence regulation in enteroaggregative E. coli

Enteroaggregative Escherichia coli (EAEC) strains are one of the diarrheagenic pathotypes. EAEC strains harbor a virulence plasmid (pAA2) that encodes, among other virulence determinants, the aggR gene. The expression of the AggR protein leads to the expression of several virulence determinants in both plasmids and chromosomes. In this work, we describe a novel mechanism that influences AggR expression. Because of the absence of a Rho-independent terminator in the 3′UTR, aggR transcripts extend far beyond the aggR ORF. These transcripts are prone to PNPase-mediated degradation. Structural alterations in the 3′UTR result in increased aggR transcript stability, leading to increased AggR levels. We therefore investigated the effect of increased AggR levels on EAEC virulence. Upon finding the previously described AggR-dependent virulence factors, we detected novel AggR-regulated genes that may play relevant roles in EAEC virulence. Mutants exhibiting high AggR levels because of structural alterations in the aggR 3′UTR show increased mobility and increased pAA2 conjugation frequency. Furthermore, among the genes exhibiting increased fold change values, we could identify those of metabolic pathways that promote increased degradation of arginine, fatty acids and gamma-aminobutyric acid (GABA), respectively. In this paper, we discuss how the AggR-dependent increase in specific metabolic pathways activity may contribute to EAEC virulence.

M6A-Mediated Upregulation of LINC00106 Promotes Stemness and Metastasis Properties of Hepatocellular Carcinoma via Sponging Let7f

Background: Hepatocellular carcinoma (HCC) cells exhibit the stemness property, which makes the patient with HCC prone to tumor recurrence and metastasis. Despite the prominent regulatory role of long non-coding RNAs (lncRNAs) in tumor stemness, the roles and molecular mechanisms of LINC00106 in HCC are poorly understood.Methods: LINC00106, let7f and periostin expression levels in tissue specimens and cell lines were assessed through qRT-PCR and immunohistochemistry (IHC). Various in vivo and in vitro assays, namely sphere/colony formation, proportion of side population cells (SP%), invasion, migration, western blot, and murine xenograft model were employed for assessing the stemness and metastatic properties of HCC cells. Luciferase reporter assays, RNA-seq, RNA pull-down, RNA immunoprecipitation (RIP) were conducted to clarificate the target gene and analyze the underlying mechanisms.Results: LINC00106 was prominently upregulated in tissues and cell lines of HCC. Patients having a high LINC00106 level exhibited a poor outcome. Under in vivo and in vitro conditions, the stemness and metastatic properties of HCC cells were augmented by LINC00106. Additionally, LINC00106 was found to sponge let7f to upregulate periostin, which lead to the activation of periostin-associated PI3K-AKT signaling pathway. Moreover, m6A methylation was found to cause LINC00106 upregulation while maintaining LINC00106 RNA transcript stability.Conclusion: m6A methylation triggers the upregulation of LINC00106, which promotes the stemness and metastasis properties in HCC cells by sponging let7f, thereby resulting in periostin activation. The findings indicate the potential of LINC00106 as a diagnostic marker and therapeutic target for HCC.

Protein Binding to Cis-Motifs in mRNAs Coding Sequence Is Common and Regulates Transcript Stability and the Rate of Translation

Protein binding to the non-coding regions of mRNAs is relatively well characterized and its functionality has been described in many examples. New results obtained by high-throughput methods indicate that binding to the coding sequence (CDS) by RNA-binding proteins is also quite common, but the functions thereof are more obscure. As described in this review, CDS binding has a role in the regulation of mRNA stability, but it has also a more intriguing role in the regulation of translational efficiency. Global approaches, which suggest the significance of CDS binding along with specific examples of CDS-binding RBPs and their modes of action, are outlined here, pointing to the existence of a relatively less-known regulatory network controlling mRNA stability and translation on yet another level.

The Function of microRNAs in Pulmonary Embolism: Review and Research Outlook

Pulmonary embolism (PE) is a common pathologic condition that frequently occurs in patients with deep venous thrombosis. Severe PE may critically suppress cardiopulmonary function, thereby threatening the life of patients. Chronic pulmonary hypertension caused by PE may lead to deterioration of respiratory dysfunction, resulting in complete disability. MicroRNAs (miRNAs) are a group of abundantly expressed non-coding RNAs that exert multiple functions in regulating the transcriptome via post-transcriptional targeting of mRNAs. Specifically, miRNAs bind to target mRNAs in a matching mechanism between the miRNA seed sequence and mRNA 3ʹ UTR, thus modulating the transcript stability or subsequent translation activity by RNA-induced silencing complex. Current studies have reported the function of miRNAs as biomarkers of PE, revealing their mechanism, function, and targetome in venous thrombophilia. This review summarizes the literature on miRNA functions and downstream mechanisms in PE. We conclude that various related miRNAs play important roles in PE and have great potential as treatment targets. For clinical application, we propose that miRNA biomarkers combined with traditional biomarkers or miRNA signatures generated from microchips may serve as a great predictive tool for PE occurrence and prognosis. Further, therapies targeting miRNAs or their upstream/downstream molecules need to be developed more quickly to keep up with the progress of routine treatments, such as anticoagulation, thrombolysis, or surgery.

The snoGloBe interaction predictor enables a broader study of box C/D snoRNA functions and mechanisms

Box C/D small nucleolar RNAs (snoRNAs) are a conserved class of noncoding RNA known to serve as guides for the site-specific 2'-O-ribose methylation of ribosomal RNAs and the U6 small nuclear RNA, through direct base pairing with the target. In recent years however, several examples of box C/D snoRNAs regulating different levels of gene expression including transcript stability and splicing have been reported. These regulatory interactions typically require direct binding of the target but do not always involve the guide region. Supporting these new box C/D snoRNA functions, high-throughput RNA-RNA interaction datasets detect many interactions between box C/D snoRNAs and messenger RNAs. To facilitate the study of box C/D snoRNA functionality, we created snoGloBe, a box C/D snoRNA machine learning target predictor based on a gradient boosting classifier and considering snoRNA and target sequence and position as well as target type. SnoGloBe convincingly outperforms general RNA duplex predictors and PLEXY, the only box C/D snoRNA-specific target predictor available. The study of snoGloBe human transcriptome-wide predictions identifies enrichment in snoRNA interactions in exons and on exon-intron junctions. Some specific snoRNAs are predicted to target groups of functionally-related transcripts on common regulatory elements and the exact position of the predicted targets strongly overlaps binding sites of RNA-binding proteins involved in relevant molecular functions. SnoGloBe was also applied to predicting interactions between human box C/D snoRNAs and the SARS-CoV-2 transcriptome, identifying known and novel interactions. Overall, snoGloBe is a timely new tool that will accelerate our understanding of C/D snoRNA targets and function.

Members of the CUGBP Elav-Like Family of RNA-Binding Proteins are Expressed in Distinct Populations of Primary Sensory Neurons

Primary sensory Dorsal Root Ganglia (DRG) neurons are diverse, with distinct populations that respond to specific stimuli. Previously, we observed that functionally distinct populations of DRG neurons express mRNA transcript variants with different 3` untranslated regions (3`UTRs). 3`UTRs harbor binding sites for interaction with RNA-binding proteins (RBPs) critical for targeting mRNAs to subcellular domains, modulating transcript stability and regulating the rate of translation. In the current study we sought to determine if 3`UTR-binding proteins are restricted to specific DRG neuron populations. Analysis of publicly available single-cell RNA-Sequencing (scRNA-Seq) data generated from adult mice revealed that 17 3`UTR-binding RBPs were enriched in specific populations of DRG neurons. This included 4 members of the CUGBP Elav-Like Family (CELF). CELF2 and CELF4 were enriched in peptidergic, CELF6 in both peptidergic and nonpeptidergic and CELF3 in tyrosine hydroxylase-expressing neurons. CELF4 is a known regulator of neural excitability, likely through modulation of protein synthesis via binding to interaction sites within the 3`UTRs of mRNAs. Immunofluorescence studies showed 60% of CELF4+ neurons are small diameter C fibers and 33% medium diameter myelinated (likely Aδ) fibers. Co-expression analyses using transcriptomic data and quantitative immunofluorescence revealed that CELF4 is enriched in nociceptive neurons that express GFRA3, CGRP and the capsaicin receptor TRPV1. Finally, genes with CELF4 binding motifs expressed in CELF4+ neurons are significantly associated with gene ontology (GO) terms such as RNA-binding and translation. We propose that CELF4 may therefore control a novel regulon that coordinates the translation of mRNAs encoding components of the protein translation apparatus in nociceptors.

5-methylcytosine modification by Plasmodium NSUN2 stabilizes mRNA and mediates the development of gametocytes

5-methylcytosine (m5C) is an important epi-transcriptomic modification involved in mRNA stability and translation efficiency in various biological processes. However, it remains unclear if m5C modification contributes to the dynamic regulation of the transcriptome during the developmental cycles of Plasmodium parasites. Here, we characterize the landscape of m5C mRNA modifications at single nucleotide resolution in the asexual replication stages and gametocyte sexual stages of rodent (P. yoelii) and human (P. falciparum) malaria parasites. While different representations of m5C-modified mRNAs are associated with the different stages, the abundance of the m5C marker is strikingly enhanced in the transcriptomes of gametocytes. Our results show that m5C modifications confer stability to the Plasmodium transcripts and that a Plasmodium ortholog of NSUN2 is a major mRNA m5C methyltransferase in malaria parasites. Upon knock-out of P. yoelii nsun2 (pynsun2), marked reductions of m5C modification were observed in a panel of gametocytogenesis-associated transcripts. These reductions correlated with impaired gametocyte production in rodent and human malaria parasites. Restoration of the nsun2 gene in the knock-out parasites rescued the gametocyte production phenotype as well as m5C modification of the gametocytogenesis-associated transcripts. Together with the mRNA m5C profiles for two species of Plasmodium, our findings demonstrate a major role for NSUN2-mediated m5C modifications in mRNA transcript stability and sexual differentiation in malaria parasites.

Characterization of the human TARDBP gene promoter

The expression of TDP-43, the main component of neuronal intracellular inclusions across a broad spectrum of ALS and FTD disorders, is developmentally regulated and studies in vivo have shown that TDP-43 overexpression can be toxic, even before observation of pathological aggregates. Starting from these observations, the regulation of its expression at transcriptional level might represent a further key element for the pathogenesis of neurodegenerative diseases. Therefore, we have characterized the human TARDBP promoter, in order to study the transcriptional mechanisms of expression. Mapping of cis-acting elements by luciferase assays in different cell outlined that the activity of the promoter seems to be higher in SH-SY5Y, Neuro2A, and HeLa than in HEK293. In addition, we tested effects of two SNPs found in the promoter region of ALS patients and observed no significant effect on transcription levels in all tested cell lines. Lastly, while TDP-43 overexpression did not affect significantly the activity of its promoter (suggesting that TDP-43 does not influence its own transcription), the presence of the 5′UTR sequence and of intron-1 splicing seem to impact positively on TDP-43 expression without affecting transcript stability. In conclusion, we have identified the region spanning nucleotides 451–230 upstream from the transcription start site as the minimal region with a significant transcription activity. These results lay an important foundation for exploring the regulation of the TARDBP gene transcription by exogenous and endogenous stimuli and the implication of transcriptional mechanisms in the pathogenesis of TDP-43 proteinopathies.