Cleavage and polyadenylation: Ending the message expands gene regulation

Identification of rare and pathogenic TAL2 gene mutations in B-lineage acute lymphoblastic leukemia (B-ALL) using mutational screening and comprehensive bioinformatics analysis

BACKGROUND

Recent genomic research has identified several genetic factors contributing to B-cell acute lymphoblastic leukemia (B-ALL). However, the exact cause of the disease is still not fully understood. It is known that mutations in the TAL2 gene play important roles in the development of acute lymphoblastic leukemia. This study aimed to analyze the molecular and computational profile of the TAL2 mutations in a group of Iranian B-ALL patients for the first time.

METHODS AND RESULTS

In this study, 188 patients were enrolled, and the TAL2 gene was sequenced to identify gene variations. The study included structural/functional analysis, homology modeling, molecular docking, and molecular dynamics (MD) simulations to assess the potential impact of the missense mutations on the protein's structure. Three nucleotide variations in the exon, three variations in the 3'UTR, and one deletion variant in the 3'UTR were detected in patients. Through in-silico analysis, it was found that the p. Asp35Glu missense mutation is located in the bHLH domain of the TAL2 protein. Also, the structural and functional analyses predicted that this mutation is a pathogenic or likely pathogenic variant in B-ALL patients. Moreover, a multiple nucleotide deletion (g.659_668del) was found in the 3'UTR in most patients. This deletion occurs at the site of poly-A tail attachment and appears to have significant implications.

CONCLUSIONS

These findings offer new insights into the impact of genetic variants in the TAL2 gene on the development of B-ALL and their potential role as tumor biomarkers for the B-ALL. Further research is needed to explore the relationship between specific TAL2 mutations and the clinical presentation of B-ALL.

Cross-disease transcriptomic analysis reveals DOK3 and PAPOLA as therapeutic targets for neuroinflammatory and tumorigenic processes

Objective

Subarachnoid hemorrhage (SAH) and tumorigenesis share numerous biological complexities; nevertheless, the specific gene expression profiles and underlying mechanisms remain poorly understood. This study aims to identify differentially expressed genes (DEGs) that could serve as biomarkers for diagnosis and prognosis.

Methods

Gene expression datasets (GSE122063, GSE13353, GSE161870) were analyzed using machine learning algorithms and logistic regression to identify DEGs associated with both SAH and tumorigenesis. Lasso regression and receiver operating characteristic (ROC) curve analysis were employed to evaluate the classification accuracy of these genes. Validation of critical DEGs was performed through pan-cancer analysis and experimental studies, focusing on the role of DOK3 in modulating inflammation and oxidative stress in U251MG glioblastoma and BV2 microglia cells.

Results

Fifteen common DEGs were identified, with DOK3 and PAPOLA highlighted as crucial genes implicated in SAH and neurodegenerative processes. Experimental validation demonstrated that DOK3 overexpression significantly reduced pro-inflammatory cytokine levels and oxidative stress markers while enhancing antioxidant enzyme activity. Additionally, DOK3 influenced tumorigenic processes such as apoptosis, cell cycle regulation, and proliferation, effectively mitigating LPS-induced cytotoxicity and inflammation in BV2 microglial cells.

Conclusions

DOK3 and PAPOLA play critical roles in both SAH and related neurodegeneration, presenting themselves as potential prognostic biomarkers and therapeutic targets. Notably, DOK3 exhibits potential as an antitumor agent with anti-inflammatory and antioxidative properties, offering therapeutic benefits for both cancer and neuroinflammatory conditions.

Recurrent disruption of tumour suppressor genes in cancer by somatic mutations in cleavage and polyadenylation signals

The expression of eukaryotic genes relies on the precise 3'-terminal cleavage and polyadenylation of newly synthesized pre-mRNA transcripts. Defects in these processes have been associated with various diseases, including cancer. While cancer-focused sequencing studies have identified numerous driver mutations in protein-coding sequences, noncoding drivers - particularly those affecting the cis-elements required for pre-mRNA cleavage and polyadenylation - have received less attention. Here, we systematically analysed somatic mutations affecting 3'UTR polyadenylation signals in human cancers using the Pan-Cancer Analysis of Whole Genomes (PCAWG) dataset. We found a striking enrichment of cancer-specific somatic mutations that disrupt strong and evolutionarily conserved cleavage and polyadenylation signals within tumour suppressor genes. Further bioinformatics and experimental analyses conducted as a part of our study suggest that these mutations have a profound capacity to downregulate the expression of tumour suppressor genes. Thus, this work uncovers a novel class of noncoding somatic mutations with significant potential to drive cancer progression.

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.

Analysis of Haemonchus embryos at single cell resolution identifies two eukaryotic elongation factors as intervention target candidates

Advances in single cell technologies are allowing investigations of a wide range of biological processes and pathways in animals, such as the multicellular model organism Caenorhabditis elegans - a free-living nematode. However, there has been limited application of such technology to related parasitic nematodes which cause major diseases of humans and animals worldwide. With no vaccines against the vast majority of parasitic nematodes and treatment failures due to drug resistance or inefficacy, new intervention targets are urgently needed, preferably informed by a deep understanding of these nematodes' cellular and molecular biology - which is presently lacking for most worms. Here, we created the first single cell atlas for an early developmental stage of Haemonchus contortus - a highly pathogenic, C. elegans-related parasitic nematode. We obtained and curated RNA sequence (snRNA-seq) data from single nuclei from embryonating eggs of H. contortus (150,000 droplets), and selected high-quality transcriptomic data for > 14,000 single nuclei for analysis, and identified 19 distinct clusters of cells. Guided by comparative analyses with C. elegans, we were able to reproducibly assign seven cell clusters to body wall muscle, hypodermis, neuronal, intestinal or seam cells, and identified eight genes that were transcribed in all cell clusters/types, three of which were inferred to be essential in H. contortus. Two of these genes (i.e. Hc-eef-1A and Hc-eef1G), coding for eukaryotic elongation factors (called Hc-eEF1A and Hc-eEF1G), were also demonstrated to be transcribed and expressed in all key developmental stages of H. contortus. Together with these findings, sequence- and structure-based comparative analyses indicated the potential of Hc-eEF1A and/or Hc-eEF1G as intervention targets within the protein biosynthesis machinery of H. contortus. Future work will focus on single cell studies of all key developmental stages and tissues of H. contortus, and on evaluating the suitability of the two elongation factor proteins as drug targets in H. contortus and related nematodes, with a view to finding new nematocidal drug candidates.

LsRTDv1, a reference transcript dataset for accurate transcript-specific expression analysis in lettuce

Accurate quantification of gene and transcript-specific expression, with the underlying knowledge of precise transcript isoforms, is crucial to understanding many biological processes. Analysis of RNA sequencing data has benefited from the development of alignment-free algorithms which enhance the precision and speed of expression analysis. However, such algorithms require a reference transcriptome. Here we generate a reference transcript dataset (LsRTDv1) for lettuce (cv. Saladin), combining long- and short-read sequencing with publicly available transcriptome annotations, and filtering to keep only transcripts with high-confidence splice junctions and transcriptional start and end sites. LsRTDv1 identifies novel genes (mostly long non-coding RNAs) and increases the number of transcript isoforms per gene in the lettuce genome from 1.4 to 2.7. We show that LsRTDv1 significantly increases the mapping rate of RNA-seq data from a lettuce time-series experiment (mock- and Botrytis cinerea-inoculated) and enables detection of genes that are differentially alternatively spliced in response to infection as well as transcript-specific expression changes. LsRTDv1 is a valuable resource for investigation of transcriptional and alternative splicing regulation in lettuce.

Cleavage and polyadenylation factors are potential regulators of adipogenesis

OBJECTIVE

Alternative polyadenylation (APA) is a co-transcriptional process that leads to isoform diversity in the 3' ends of mRNAs. APA is known to occur during differentiation, and its dysregulation is observed in diseases like cancer and autoimmune disorders. It has been previously reported that differentiation of 3T3-L1 cells to adipocytes leads to an overall lengthening of mRNAs, but the proteins involved in this regulation have not been identified. The expression levels of subunits of the cleavage and polyadenylation (C/P) complex can regulate the choice of poly(A) site, which in turn can affect different cellular activities. In this paper, we studied the change in levels of C/P proteins during 3T3-L1 differentiation.

RESULTS

We observed that while the RNA expression of these proteins is unchanged during differentiation, the protein levels of some subunits do change, including a decrease in levels of CPSF73, the nuclease that cuts at the poly(A) site. However, overexpression of CPSF73 alone does not affect the efficiency and rate of differentiation.

DEAD-box ATPase Dbp2 is the key enzyme in an mRNP assembly checkpoint at the 3'-end of genes and involved in the recycling of cleavage factors

mRNA biogenesis in the eukaryotic nucleus is a highly complex process. The numerous RNA processing steps are tightly coordinated to ensure that only fully processed transcripts are released from chromatin for export from the nucleus. Here, we present the hypothesis that fission yeast Dbp2, a ribonucleoprotein complex (RNP) remodelling ATPase of the DEAD-box family, is the key enzyme in an RNP assembly checkpoint at the 3'-end of genes. We show that Dbp2 interacts with the cleavage and polyadenylation complex (CPAC) and localises to cleavage bodies, which are enriched for 3'-end processing factors and proteins involved in nuclear RNA surveillance. Upon loss of Dbp2, 3'-processed, polyadenylated RNAs accumulate on chromatin and in cleavage bodies, and CPAC components are depleted from the soluble pool. Under these conditions, cells display an increased likelihood to skip polyadenylation sites and a delayed transcription termination, suggesting that levels of free CPAC components are insufficient to maintain normal levels of 3'-end processing. Our data support a model in which Dbp2 is the active component of an mRNP remodelling checkpoint that licenses RNA export and is coupled to CPAC release.

ntEmbd: Deep learning embedding for nucleotide sequences

Enabled by the explosion of data and substantial increase in computational power, deep learning has transformed fields such as computer vision and natural language processing (NLP) and it has become a successful method to be applied to many transcriptomic analysis tasks. A core advantage of deep learning is its inherent capability to incorporate feature computation within the machine learning models. This results in a comprehensive and machine-readable representation of sequences, facilitating the downstream classification and clustering tasks. Compared to machine translation problems in NLP, feature embedding is particularly challenging for transcriptomic studies as the sequences are string of thousands of nucleotides in length, which make the long-term dependencies between features from different parts of the sequence even more difficult to capture. This highlights the need for nucleotide sequence embedding methods that are capable of learning input sequence features implicitly. Here we introduce ntEmbd, a deep learning embedding tool that captures dependencies between different features of the sequences and learns a latent representation for given nucleotide sequences. We further provide two sample use cases, describing how learned RNA features can be used in downstream analysis. The first use case demonstrates ntEmbd's utility in classifying coding and noncoding RNA benchmarked against existing tools, and the second one explores the utility of learned representations in identifying adapter sequences in nanopore RNA-seq reads. The tool as well as the trained models are freely available on GitHub at https://github.com/bcgsc/ntEmbd.

Star-PAP controls oncogene expression through primary miRNA 3'-end formation to regulate cellular proliferation and tumour formation

Star-PAP is a non-canonical poly(A) polymerase that is down regulated in breast cancer. While Star-PAP down regulation impairs target mRNA polyadenylation, paradoxically, we see up regulation of a large number of oncogenes on Star-PAP knockdown. Using two breast cancer cells (MCF7 with high Star-PAP, and MDA-MB-231 with negligible Star-PAP level), we discover that Star-PAP negatively regulates oncogene expression and subsequently cellular proliferation. This regulation is compromised with Star-PAP mutant of 3'-end processing function (serine 6 to alanine, S6A phospho-mutation). Concomitantly, xenograft mice model using MDA-MB-231 cells reveals a reduction in the tumour formation on ectopic Star-PAP expression that is ameliorated by S6A mutation. We find that Star-PAP control of target oncogene expression is independent of Star-PAP-mediated alternative polyadenylation or target mRNA 3'-end formation. We demonstrate that Star-PAP regulates target oncogenes through cellular miRNAs (miR-421, miR-335, miR-424, miR-543, miR-205, miR-34a, and miR-26a) that are down regulated in breast cancer. Analysis of various steps in miRNA biogenesis pathway reveals that Star-PAP regulates 3'-end formation and synthesis of primary miRNA (host) transcripts that is dependent on S6 phosphorylation thus controlling mature miRNA generation. Using mimics and inhibitors of two target miRNAs (miR-421 and miR-424) after Star-PAP depletion in MCF7 or ectopic expression in MDA-MB-231 cells, we demonstrate that Star-PAP controls oncogene expression and cellular proliferation through targeting miRNAs that regulates tumour formation. Our study establishes a novel mechanism of oncogene expression independent of alternative polyadenylation through Star-PAP-mediated miRNA host transcript polyadenylation that regulates breast cancer progression.

Genome-wide kinetic profiling of pre-mRNA 3' end cleavage

Cleavage and polyadenylation is necessary for the formation of mature mRNA molecules. The rate at which this process occurs can determine the temporal availability of mRNA for subsequent function throughout the cell and is likely tightly regulated. Despite advances in high-throughput approaches for global kinetic profiling of RNA maturation, genome-wide 3' end cleavage rates have never been measured. Here, we describe a novel approach to estimate the rates of cleavage, using metabolic labeling of nascent RNA, high-throughput sequencing, and mathematical modeling. Using in silico simulations of nascent RNA-seq data, we show that our approach can accurately and precisely estimate cleavage half-lives for both constitutive and alternative sites. We find that 3' end cleavage is fast on average, with half-lives under a minute, but highly variable across individual sites. Rapid cleavage is promoted by the presence of canonical sequence elements and an increased density of polyadenylation signals near a cleavage site. Finally, we find that cleavage rates are associated with the localization of RNA polymerase II at the end of a gene, and faster cleavage leads to quicker degradation of downstream readthrough RNA. Our findings shed light on the features important for efficient 3' end cleavage and the regulation of transcription termination.

Gene expression investigation of four key regulators of polyadenylation and alternative adenylation in the periphery of late-onset Alzheimer's disease patients

BACKGROUND

Alzheimer's disease (AD) is a genetic and sporadic neurodegenerative disease considered by an archetypal cognitive impairment and a decrease in less common cognitive impairment. Notably, the discovery of goals in this paradigm is still a challenge, and understanding basic mechanisms is an important step toward improving disease management. Polyadenylation (PA) and alternative polyadenylation (APA) are two of the most critical RNA processing stages in 3'UTRs that influence various AD-related genes.

METHODS

In this study, we assessed Cleavage and polyadenylation specificity factors 1 and 6 (CPSF1 and CPSF6), cleavage stimulation factor 1 (CSTF1), and WD Repeat Domain 33 (WDR33) genes expression in the periphery of 50 AD patients and 50 healthy individuals with age and gender-matched by quantitative real-time PCR.

RESULTS

Comparing AD patients with healthy people using expression analysis revealed a substantial increase in CSTF1 (posterior beta = 0.773, adjusted P-value = 0.042). Significant positive correlations were found between CSTF1 and CPSF1 (r = 0.365, P < 0.001), WDR33 (r = 0.506, P < 0.001), and CPSF6 (r = 0.446, P < 0.001) expression levels.

CONCLUSION

Although further research is required to determine their potential contribution to AD, our findings offer a fresh perspective on molecular regulatory pathways associated with AD pathogenic mechanisms associated with PA and APA.

Mechanisms and Functions of the RNA Polymerase II General Transcription Machinery during the Transcription Cycle

Central to the development and survival of all organisms is the regulation of gene expression, which begins with the process of transcription catalyzed by RNA polymerases. During transcription of protein-coding genes, the general transcription factors (GTFs) work alongside RNA polymerase II (Pol II) to assemble the preinitiation complex at the transcription start site, open the promoter DNA, initiate synthesis of the nascent messenger RNA, transition to productive elongation, and ultimately terminate transcription. Through these different stages of transcription, Pol II is dynamically phosphorylated at the C-terminal tail of its largest subunit, serving as a control mechanism for Pol II elongation and a signaling/binding platform for co-transcriptional factors. The large number of core protein factors participating in the fundamental steps of transcription add dense layers of regulation that contribute to the complexity of temporal and spatial control of gene expression within any given cell type. The Pol II transcription system is highly conserved across different levels of eukaryotes; however, most of the information here will focus on the human Pol II system. This review walks through various stages of transcription, from preinitiation complex assembly to termination, highlighting the functions and mechanisms of the core machinery that participates in each stage.

TDP-43 loss induces extensive cryptic polyadenylation in ALS/FTD

Nuclear depletion and cytoplasmic aggregation of the RNA-binding protein TDP-43 is the hallmark of ALS, occurring in over 97% of cases. A key consequence of TDP-43 nuclear loss is the de-repression of cryptic exons. Whilst TDP-43 regulated cryptic splicing is increasingly well catalogued, cryptic alternative polyadenylation (APA) events, which define the 3' end of last exons, have been largely overlooked, especially when not associated with novel upstream splice junctions. We developed a novel bioinformatic approach to reliably identify distinct APA event types: alternative last exons (ALE), 3'UTR extensions (3'Ext) and intronic polyadenylation (IPA) events. We identified novel neuronal cryptic APA sites induced by TDP-43 loss of function by systematically applying our pipeline to a compendium of publicly available and in house datasets. We find that TDP-43 binding sites and target motifs are enriched at these cryptic events and that TDP-43 can have both repressive and enhancing action on APA. Importantly, all categories of cryptic APA can also be identified in ALS and FTD post mortem brain regions with TDP-43 proteinopathy underlining their potential disease relevance. RNA-seq and Ribo-seq analyses indicate that distinct cryptic APA categories have different downstream effects on transcript and translation. Intriguingly, cryptic 3'Exts occur in multiple transcription factors, such as ELK1, SIX3, and TLX1, and lead to an increase in wild-type protein levels and function. Finally, we show that an increase in RNA stability leading to a higher cytoplasmic localisation underlies these observations. In summary, we demonstrate that TDP-43 nuclear depletion induces a novel category of cryptic RNA processing events and we expand the palette of TDP-43 loss consequences by showing this can also lead to an increase in normal protein translation.

WDR33 alternative polyadenylation is dependent on stochastic poly(a) site usage and splicing efficiencies

Transcripts from the human WDR33 gene, which encodes a central component of the mRNA polyadenylation (PA) machinery, are subject to alternative polyadenylation (APA) within promoter-proximal introns/exons. This APA, which itself involves usage of multiple PA sites, results in the production of two non-canonical protein isoforms, V2 and V3, that are functionally completely unrelated to the full-length protein, with roles in innate immunity. The mechanism and regulation of WDR33 APA are unclear. Here, we report that levels of the PA factor CFIm25 modulate V2 and V3 expression, and that PA site usage of both V2 and V3 varies in distinct immune responses. Using newly developed assays to measure splicing and PA site strength, we show that splicing of V2-associated intron 6 is inefficient, allowing V2 to be produced using weak PA sites. Usage of V3's strong PA sites, on the other hand, is relatively low, reflecting the high efficiency of intron 7 splicing coupled with dependency on usage of an alternative 3' splice site within the intron. Overall, our findings demonstrate that usage of WDR33 alternative PA sites is stochastic, dependent on a complex interplay between splicing and PA, and thus provide new insights into mechanisms underlying APA.

Efficient 3'-pairing renders microRNA targeting less sensitive to mRNA seed accessibility

MicroRNAs (miRNAs) are short RNAs that post-transcriptionally regulate gene expression by binding to specific sites in mRNAs. Site recognition is primarily mediated by the seed region (nucleotides g2-g8 in the miRNA), but pairing beyond the seed (3'-pairing) is important for some miRNA:target interactions. Here, we use SHAPE, luciferase reporter assays and transcriptomics analyses to study the combined effect of 3'-pairing and secondary structures in mRNAs on repression efficiency. Using the interaction between miR-34a and its SIRT1 binding site as a model, we provide structural and functional evidence that 3'-pairing can compensate for low seed-binding site accessibility, enabling repression of sites that would otherwise be ineffective. We show that miRNA 3'-pairing regions can productively base-pair with nucleotides far upstream of the seed-binding site and that both hairpins and unstructured bulges within the target site are tolerated. We use SHAPE to show that sequences that overcome inaccessible seed-binding sites by strong 3'-pairing adopt the predicted structures and corroborate the model using luciferase assays and high-throughput modelling of 8177 3'-UTR targets for six miRNAs. Finally, we demonstrate that PHB2, a target of miR-141, is an inaccessible target rescued by efficient 3'-pairing. We propose that these results could refine predictions of effective target sites.

Versatile RNA: overlooked gems of the transcriptome

The critical role of RNA, its use and targetability concerning different aspects of human health are gaining more attention because our understanding of the versatility of RNA has dramatically evolved over the last decades. We now appreciate that RNA is far more critical than a messenger molecule and possesses many complicated functions. As a multifunctional molecule with its sequence, flexible structures and enzymatic abilities, RNA is genuinely powerful. Mammalian transcriptomes consist of a dynamically regulated plethora of coding and noncoding RNA types. However, some aspects of RNA metabolism remain to be explored. In this Viewpoint, we focus on the transcriptome's unconventional and possibly overlooked aspects to emphasize the importance of RNA in mammalian systems.

Near chromosome-level genome assembly of the microsporidium Hamiltosporidium tvaerminnensis

Microsporidia are intracellular parasitic fungi whose genomes rank among the smallest of all known eukaryotes. A number of outstanding questions remain concerning the evolution of their large-scale variation in genome architecture, responsible for genome size variation of more than an order of magnitude. This genome report presents the first near-chromosomal assembly of a large-genome microsporidium, Hamiltosporidium tvaerminnensis. Combined Oxford Nanopore, Pacific Biosciences (PacBio), and Illumina sequencing led to a genome assembly of 17 contigs, 11 of which represent complete chromosomes. Our assembly is 21.64 Mb in length, has an N50 of 1.44 Mb, and consists of 39.56% interspersed repeats. We introduce a novel approach in microsporidia, PacBio Iso-Seq, as part of a larger annotation pipeline for obtaining high-quality annotations of 3,573 protein-coding genes. Based on direct evidence from the full-length Iso-Seq transcripts, we present evidence for alternative polyadenylation and variation in splicing efficiency, which are potential regulation mechanisms for gene expression in microsporidia. The generated high-quality genome assembly is a necessary resource for comparative genomics that will help elucidate the evolution of genome architecture in response to intracellular parasitism.

Mapping the splicing landscape of the human immune system

Most human genes code for more than one transcript. Different ratios of transcripts of the same gene can be found in different cell types or states, indicating differential use of transcription start sites or differential splicing. Such differential transcript use (DTUs) events provide an additional layer of regulation and protein diversity. With the exceptions of PTPRC and CIITA, there are very few reported cases of DTU events in the immune system. To rigorously map DTUs between different human immune cell types, we leveraged four publicly available RNA sequencing datasets. We identified 282 DTU events between five human healthy immune cell types that appear in at least two datasets. The patterns of the DTU events were mostly cell-type-specific or lineage-specific, in the context of the five cell types tested. DTUs correlated with the expression pattern of potential regulators, namely, splicing factors and transcription factors. Of the several immune related conditions studied, only sepsis affected the splicing of more than a few genes and only in innate immune cells. Taken together, we map the DTUs landscape in human peripheral blood immune cell types, and present hundreds of genes whose transcript use changes between cell types or upon activation.

Extension of mRNA poly(A) tails and 3'UTRs during neuronal differentiation exhibits variable association with post-transcriptional dynamics

Differentiation of neural progenitor cells into mature neuronal phenotypes relies on extensive temporospatial coordination of mRNA expression to support the development of functional brain circuitry. Cleavage and polyadenylation of mRNA has tremendous regulatory capacity through the alteration of mRNA stability and modulation of microRNA (miRNA) function, however the extent of utilization in neuronal development is currently unclear. Here, we employed poly(A) tail sequencing, mRNA sequencing, ribosome profiling and small RNA sequencing to explore the functional relationship between mRNA abundance, translation, poly(A) tail length, alternative polyadenylation (APA) and miRNA expression in an in vitro model of neuronal differentiation. Differential analysis revealed a strong bias towards poly(A) tail and 3'UTR lengthening during differentiation, both of which were positively correlated with changes in mRNA abundance, but not translation. Globally, changes in miRNA expression were predominantly associated with mRNA abundance and translation, however several miRNA-mRNA pairings with potential to regulate poly(A) tail length were identified. Furthermore, 3'UTR lengthening was observed to significantly increase the inclusion of non-conserved miRNA binding sites, potentially enhancing the regulatory capacity of these molecules in mature neuronal cells. Together, our findings suggest poly(A) tail length and APA function as part of a rich post-transcriptional regulatory matrix during neuronal differentiation.

mRNA 3'UTR lengthening by alternative polyadenylation attenuates inflammatory responses and correlates with virulence of Influenza A virus

Changes of mRNA 3'UTRs by alternative polyadenylation (APA) have been associated to numerous pathologies, but the mechanisms and consequences often remain enigmatic. By combining transcriptomics, proteomics and recombinant viruses we show that all tested strains of IAV, including A/PR/8/34(H1N1) (PR8) and A/Cal/07/2009 (H1N1) (Cal09), cause APA. We mapped the effect to the highly conserved glycine residue at position 184 (G184) of the viral non-structural protein 1 (NS1). Unbiased mass spectrometry-based analyses indicate that NS1 causes APA by perturbing the function of CPSF4 and that this function is unrelated to virus-induced transcriptional shutoff. Accordingly, IAV strain PR8, expressing an NS1 variant with weak CPSF binding, does not induce host shutoff but only APA. However, recombinant IAV (PR8) expressing NS1(G184R) lacks binding to CPSF4 and thereby also the ability to cause APA. Functionally, the impaired ability to induce APA leads to an increased inflammatory cytokine production and an attenuated phenotype in a mouse infection model. Investigating diverse viral infection models showed that APA induction is a frequent ability of many pathogens. Collectively, we propose that targeting of the CPSF complex, leading to widespread alternative polyadenylation of host transcripts, constitutes a general immunevasion mechanism employed by a variety of pathogenic viruses.

Towards understandings of serine/arginine-rich splicing factors

Serine/arginine-rich splicing factors (SRSFs) refer to twelve RNA-binding proteins which regulate splice site recognition and spliceosome assembly during precursor messenger RNA splicing. SRSFs also participate in other RNA metabolic events, such as transcription, translation and nonsense-mediated decay, during their shuttling between nucleus and cytoplasm, making them indispensable for genome diversity and cellular activity. Of note, aberrant SRSF expression and/or mutations elicit fallacies in gene splicing, leading to the generation of pathogenic gene and protein isoforms, which highlights the therapeutic potential of targeting SRSF to treat diseases. In this review, we updated current understanding of SRSF structures and functions in RNA metabolism. Next, we analyzed SRSF-induced aberrant gene expression and their pathogenic outcomes in cancers and non-tumor diseases. The development of some well-characterized SRSF inhibitors was discussed in detail. We hope this review will contribute to future studies of SRSF functions and drug development targeting SRSFs.

Comparative transcriptome characterization of esophageal squamous cell carcinoma and adenocarcinoma

Background

Esophageal cancers are primarily categorized as esophageal squamous cell carcinoma (ESCC) and esophageal adenocarcinoma (EAC). While various (epi) genomic alterations associated with tumor development in ESCC and EAC have been documented, a comprehensive comparison of the transcriptomes in these two cancer subtypes remains lacking.

Methods

We collected 551 gene expression profiles from publicly available sources, including normal, ESCC, and EAC tissues or cell lines. Subsequently, we conducted a systematic analysis to compare the transcriptomes of these samples at various levels, including gene expression, promoter activity, alternative splicing (AS), alternative polyadenylation (APA), and gene fusion.

Results

Seven distinct cluster gene expression patterns were identified among the differentially expressed genes in normal, ESCC, and EAC tissues. These patterns were enriched in the PI3K-Akt signaling pathway and the activation of extracellular matrix organization and exhibited repression of epidermal development. Notably, we observed additional genes or unique expression levels enriched in these shared pathways and biological processes related to tumor development and immune activation. In addition to the differentially expressed genes, there was an enrichment of lncRNA co-expression networks and downregulation of promoter activity associated with the repression of epidermal development in both ESCC and EAC. This indicates a common feature between these two cancer subtypes. Furthermore, differential AS and APA patterns in ESCC and EAC appear to partially affect the expression of host genes associated with bacterial or viral infections in these subtypes. No gene fusions were observed between ESCC and EAC, thus highlighting the distinct molecular mechanisms underlying these two cancer subtypes.

Conclusions

We conducted a comprehensive comparison of ESCC and EAC transcriptomes and uncovered shared and distinct transcriptomic signatures at multiple levels. These findings suggest that ESCC and EAC may exhibit common and unique mechanisms involved in tumorigenesis.

Gene co-expression network identifies critical genes, pathways and regulatory motifs mediating the progression of rift valley fever in Bostaurus

Rift Valley Fever (RVF) is a mosquito-borne viral disease caused by the Rift Valley Fever Virus. The disease is a zoonosis that largely affects domestic animals, including sheep, goats, and cattle, resulting in severe morbidity and mortality marked by massive storm abortions. To halt human and livestock deaths due to RVF, the development of efficacious vaccines and therapeutics is a compelling and urgent priority. We sought to identify potential key modules (gene clusters), hub genes, and regulatory motifs involved in the pathogenesis of RVF in Bos taurus that are amenable to inhibition. We analyzed 39 Bos taurus RNA-Seq samples using the weighted gene co-expression network analysis (WGCNA) R package and uncovered significantly enriched modules containing genes with potential pivotal roles in RVF progression. Moreover, regulatory motif analysis conducted using the Multiple Expectation Maximization for Motif Elicitation (MEME) suite identified motifs that probably modulate vital biological processes. Gene ontology terms associated with identified motifs were inferred using the GoMo human database. The gene co-expression network constructed in WGCNA using 5000 genes contained seven (7) modules, out of which four were significantly enriched for terms associated with response to viruses, response to interferon-alpha, innate immune response, and viral defense. Additionally, several biological pathways implicated in developmental processes, anatomical structure development, and multicellular organism development were identified. Regulatory motifs analysis identified short, repeated motifs whose function(s) may be amenable to disruption by novel therapeutics. Predicted functions of identified motifs include tissue development, embryonic organ development, and organ morphogenesis. We have identified several hub genes in enriched co-expressed gene modules and regulatory motifs potentially involved in the pathogenesis of RVF in B. taurus that are likely viable targets for disruption by novel therapeutics.

RNA-dependent interactome allows network-based assignment of RNA-binding protein function

RNA-binding proteins (RBPs) form highly diverse and dynamic ribonucleoprotein complexes, whose functions determine the molecular fate of the bound RNA. In the model organism Sacchromyces cerevisiae, the number of proteins identified as RBPs has greatly increased over the last decade. However, the cellular function of most of these novel RBPs remains largely unexplored. We used mass spectrometry-based quantitative proteomics to systematically identify protein-protein interactions (PPIs) and RNA-dependent interactions (RDIs) to create a novel dataset for 40 RBPs that are associated with the mRNA life cycle. Domain, functional and pathway enrichment analyses revealed an over-representation of RNA functionalities among the enriched interactors. Using our extensive PPI and RDI networks, we revealed putative new members of RNA-associated pathways, and highlighted potential new roles for several RBPs. Our RBP interactome resource is available through an online interactive platform as a community tool to guide further in-depth functional studies and RBP network analysis (https://www.butterlab.org/RINE).

Viral proteogenomic and expression profiling during productive replication of a skin-tropic herpesvirus in the natural host

Efficient transmission of herpesviruses is essential for dissemination in host populations; however, little is known about the viral genes that mediate transmission, mostly due to a lack of natural virus-host model systems. Marek's disease is a devastating herpesviral disease of chickens caused by Marek's disease virus (MDV) and an excellent natural model to study skin-tropic herpesviruses and transmission. Like varicella zoster virus that causes chicken pox in humans, the only site where infectious cell-free MD virions are efficiently produced is in epithelial skin cells, a requirement for host-to-host transmission. Here, we enriched for heavily infected feather follicle epithelial skin cells of live chickens to measure both viral transcription and protein expression using combined short- and long-read RNA sequencing and LC/MS-MS bottom-up proteomics. Enrichment produced a previously unseen breadth and depth of viral peptide sequencing. We confirmed protein translation for 84 viral genes at high confidence (1% FDR) and correlated relative protein abundance with RNA expression levels. Using a proteogenomic approach, we confirmed translation of most well-characterized spliced viral transcripts and identified a novel, abundant isoform of the 14 kDa transcript family via IsoSeq transcripts, short-read intron-spanning sequencing reads, and a high-quality junction-spanning peptide identification. We identified peptides representing alternative start codon usage in several genes and putative novel microORFs at the 5' ends of two core herpesviral genes, pUL47 and ICP4, along with strong evidence of independent transcription and translation of the capsid scaffold protein pUL26.5. Using a natural animal host model system to examine viral gene expression provides a robust, efficient, and meaningful way of validating results gathered from cell culture systems.

A Putative TRAPα Protein of Microsporidia Nosema bombycis Exhibits Non-Canonical Alternative Polyadenylation in Transcripts

Microsporidia are obligate intracellular eukaryotic parasites that have significantly reduced genomes and that have lost most of their introns. In the current study, we characterized a gene in microsporidia Nosema bombycis, annotated as TRAPα (HNbTRAPα). The homologous of TRAPα are a functional component of ER translocon and facilitates the initiation of protein translocation in a substrate-specific manner, which is conserved in animals but absent from most fungi. The coding sequence of HNbTRAPα consists of 2226 nucleotides, longer than the majority of homologs in microsporidia. A 3′ RACE analysis indicated that there were two mRNA isoforms resulting from non-canonical alternative polyadenylation (APA), and the polyadenylate tail was synthesized after the C951 or C1167 nucleotide, respectively. Indirect immunofluorescence analysis showed two different localization characteristics of HNbTRAPα, which are mainly located around the nuclear throughout the proliferation stage and co-localized with the nuclear in mature spores. This study demonstrated that the post-transcriptional regulation mechanism exists in Microsporidia and expands the mRNA isoform repertoire.

Genome-wide screening for genetic variants in polyadenylation signal (PAS) sites in mouse selection lines for fatness and leanness

Alternative polyadenylation (APA) determines mRNA stability, localisation, translation and protein function. Several diseases, including obesity, have been linked to APA. Studies have shown that single nucleotide polymorphisms in polyadenylation signals (PAS-SNPs) can influence APA and affect phenotype and disease susceptibility. However, these studies focussed on associations between single PAS-SNP alleles with very large effects and phenotype. Therefore, we performed a genome-wide screening for PAS-SNPs in the polygenic mouse selection lines for fatness and leanness by whole-genome sequencing. The genetic variants identified in the two lines were overlapped with locations of PAS sites obtained from the PolyASite 2.0 database. Expression data for selected genes were extracted from the microarray expression experiment performed on multiple tissue samples. In total, 682 PAS-SNPs were identified within 583 genes involved in various biological processes, including transport, protein modifications and degradation, cell adhesion and immune response. Moreover, 63 of the 583 orthologous genes in human have been previously associated with human diseases, such as nervous system and physical disorders, and immune, endocrine, and metabolic diseases. In both lines, PAS-SNPs have also been identified in genes broadly involved in APA, such as Polr2c, Eif3e and Ints11. Five PAS-SNPs within 5 genes (Car, Col4a1, Itga7, Lat, Nmnat1) were prioritised as potential functional variants and could contribute to the phenotypic disparity between the two selection lines. The developed PAS-SNPs catalogue presents a key resource for planning functional studies to uncover the role of PAS-SNPs in APA, disease susceptibility and fat deposition.

Comparative genomics and interactomics of polyadenylation factors for the prediction of new parasite targets: Entamoeba histolytica as a working model

Protein-protein interactions (PPI) play a key role in predicting the function of a target protein and drug ability to affect an entire biological system. Prediction of PPI networks greatly contributes to determine a target protein and signal pathways related to its function. Polyadenylation of mRNA 3'-end is essential for gene expression regulation and several polyadenylation factors have been shown as valuable targets for controlling protozoan parasites that affect human health. Here, by using a computational strategy based on sequence-based prediction approaches, phylogenetic analyses, and computational prediction of PPI networks, we compared interactomes of polyadenylation factors in relevant protozoan parasites and the human host, to identify key proteins and define potential targets for pathogen control. Then, we used Entamoeba histolytica as a working model to validate our computational results. RT-qPCR assays confirmed the coordinated modulation of connected proteins in the PPI network and evidenced that silencing of the bottleneck protein EhCFIm25 affects the expression of interacting proteins. In addition, molecular dynamics simulations and docking approaches allowed to characterize the relationships between EhCFIm25 and Ehnopp34, two connected bottleneck proteins. Interestingly, the experimental identification of EhCFIm25 interactome confirmed the close relationships among proteins involved in gene expression regulation and evidenced new links with moonlight proteins in E. histolytica, suggesting a connection between RNA biology and metabolism as described in other organisms. Altogether, our results strengthened the relevance of comparative genomics and interactomics of polyadenylation factors for the prediction of new targets for the control of these human pathogens.

Macrophage differentiation is marked by increased abundance of the mRNA 3' end processing machinery, altered poly(A) site usage, and sensitivity to the level of CstF64

Regulation of mRNA polyadenylation is important for response to external signals and differentiation in several cell types, and results in mRNA isoforms that vary in the amount of coding sequence or 3' UTR regulatory elements. However, its role in differentiation of monocytes to macrophages has not been investigated. Macrophages are key effectors of the innate immune system that help control infection and promote tissue-repair. However, overactivity of macrophages contributes to pathogenesis of many diseases. In this study, we show that macrophage differentiation is characterized by shortening and lengthening of mRNAs in relevant cellular pathways. The cleavage/polyadenylation (C/P) proteins increase during differentiation, suggesting a possible mechanism for the observed changes in poly(A) site usage. This was surprising since higher C/P protein levels correlate with higher proliferation rates in other systems, but monocytes stop dividing after induction of differentiation. Depletion of CstF64, a C/P protein and known regulator of polyadenylation efficiency, delayed macrophage marker expression, cell cycle exit, attachment, and acquisition of structural complexity, and impeded shortening of mRNAs with functions relevant to macrophage biology. Conversely, CstF64 overexpression increased use of promoter-proximal poly(A) sites and caused the appearance of differentiated phenotypes in the absence of induction. Our findings indicate that regulation of polyadenylation plays an important role in macrophage differentiation.

UBE3D Regulates mRNA 3'-End Processing and Maintains Adipogenic Potential in 3T3-L1 Cells

We have previously described the role of an essential Saccharomyces cerevisiae gene, important for cleavage and polyadenylation 1 (IPA1), in the regulation of gene expression through its interaction with Ysh1, the endonuclease subunit of the mRNA 3'-end processing complex. Through a similar mechanism, the mammalian homolog ubiquitin protein ligase E3D (UBE3D) promotes the migratory and invasive potential of breast cancer cells, but its role in the regulation of gene expression during normal cellular differentiation has not previously been described. In this study, we show that CRISPR/Cas9-mediated knockout of Ube3d in 3T3-L1 cells blocks their ability to differentiate into mature adipocytes. Consistent with previous studies in other cell types, Ube3d knockout leads to decreased levels of CPSF73 and global changes in cellular mRNAs indicative of a loss of 3'-end processing capacity. Ube3d knockout cells also display decreased expression of known preadipogenic markers. Overexpression of either UBE3D or CPSF73 rescues the differentiation defect and partially restores protein levels of these markers. These results support a model in which UBE3D is necessary for the maintenance of the adipocyte-committed state via its regulation of the mRNA 3'-end processing machinery.

Slight Variations in the Sequence Downstream of the Polyadenylation Signal Significantly Increase Transgene Expression in HEK293T and CHO Cells

Compared to transcription initiation, much less is known about transcription termination. In particular, large-scale mutagenesis studies have, so far, primarily concentrated on promoter and enhancer, but not terminator sequences. Here, we used a massively parallel reporter assay (MPRA) to systematically analyze the influence of short (8 bp) sequence variants (mutations) located downstream of the polyadenylation signal (PAS) on the steady-state mRNA level of the upstream gene, employing an eGFP reporter and human HEK293T cells as a model system. In total, we evaluated 227,755 mutations located at different overlapping positions within +17..+56 bp downstream of the PAS for their ability to regulate the reporter gene expression. We found that the positions +17..+44 bp downstream of the PAS are more essential for gene upregulation than those located more distal to the PAS, and that the mutation sequences ensuring high levels of eGFP mRNA expression are extremely T-rich. Next, we validated the positive effect of a couple of mutations identified in the MPRA screening on the eGFP and luciferase protein expression. The most promising mutation increased the expression of the reporter proteins 13-fold and sevenfold on average in HEK293T and CHO cells, respectively. Overall, these findings might be useful for further improving the efficiency of production of therapeutic products, e.g., recombinant antibodies.

TENT5 cytoplasmic noncanonical poly(A) polymerases regulate the innate immune response in animals

Innate immunity is the first line of host defense against pathogens. Here, through global transcriptome and proteome analyses, we uncover that newly described cytoplasmic poly(A) polymerase TENT-5 (terminal nucleotidyltransferase 5) enhances the expression of secreted innate immunity effector proteins in Caenorhabditis elegans. Direct RNA sequencing revealed that multiple mRNAs with signal peptide-encoding sequences have shorter poly(A) tails in tent-5-deficient worms. Those mRNAs are translated at the endoplasmic reticulum where a fraction of TENT-5 is present, implying that they represent its direct substrates. Loss of tent-5 makes worms more susceptible to bacterial infection. Notably, the role of TENT-5 in innate immunity is evolutionarily conserved. Its orthologs, TENT5A and TENT5C, are expressed in macrophages and induced during their activation. Analysis of macrophages devoid of TENT5A/C revealed their role in the regulation of secreted proteins involved in defense response. In summary, our study reveals cytoplasmic polyadenylation to be a previously unknown component of the posttranscriptional regulation of innate immunity in animals.

Transcriptomic complexity of the human malaria parasite Plasmodium falciparum revealed by long-read sequencing

The Plasmodium falciparum human malaria parasite genome is incompletely annotated and does not accurately represent the transcriptomic diversity of this species. To address this need, we performed long-read transcriptomic sequencing. 5' capped mRNA was enriched from samples of total and nuclear-fractionated RNA from intra-erythrocytic stages and converted to cDNA library. The cDNA libraries were sequenced on PacBio and Nanopore long-read platforms. 12,495 novel isoforms were annotated from the data. Alternative 5' and 3' ends represent the majority of isoform events among the novel isoforms, with retained introns being the next most common event. The majority of alternative 5' ends correspond to genomic regions with features similar to those of the reference transcript 5' ends. However, a minority of alternative 5' ends showed markedly different features, including locations within protein-coding regions. Alternative 3' ends showed similar features to the reference transcript 3' ends, notably adenine-rich termination signals. Distinguishing features of retained introns could not be observed, except for a tendency towards shorter length and greater GC content compared with spliced introns. Expression of antisense and retained intron isoforms was detected at different intra-erythrocytic stages, suggesting developmental regulation of these isoform events. To gain insights into the possible functions of the novel isoforms, their protein-coding potential was assessed. Variants of P. falciparum proteins and novel proteins encoded by alternative open reading frames suggest that P. falciparum has a greater proteomic repertoire than the current annotation. We provide a catalog of annotated transcripts and encoded alternative proteins to support further studies on gene and protein regulation of this pathogen.

A Survey on Methods for Predicting Polyadenylation Sites from DNA Sequences, Bulk RNA-seq, and Single-cell RNA-seq

Alternative polyadenylation (APA) plays important roles in modulating mRNA stability, translation, and subcellular localization, and contributes extensively to shaping eukaryotic transcriptome complexity and proteome diversity. Identification of poly(A) sites (pAs) on a genome-wide scale is a critical step toward understanding the underlying mechanism of APA-mediated gene regulation. A number of established computational tools have been proposed to predict pAs from diverse genomic data. Here we provided an exhaustive overview of computational approaches for predicting pAs from DNA sequences, bulk RNA sequencing (RNA-seq) data, and single-cell RNA sequencing (scRNA-seq) data. Particularly, we examined several representative tools using bulk RNA-seq and scRNA-seq data from peripheral blood mononuclear cells and put forward operable suggestions on how to assess the reliability of pAs predicted by different tools. We also proposed practical guidelines on choosing appropriate methods applicable to diverse scenarios. Moreover, we discussed in depth the challenges in improving the performance of pA prediction and benchmarking different methods. Additionally, we highlighted outstanding challenges and opportunities using new machine learning and integrative multi-omics techniques, and provided our perspective on how computational methodologies might evolve in the future for non-3' untranslated region, tissue-specific, cross-species, and single-cell pA prediction.

The Role of RNA-Binding Proteins in Hematological Malignancies

Hematological malignancies comprise a plethora of different neoplasms, such as leukemia, lymphoma, and myeloma, plus a myriad of dysplasia, such as myelodysplastic syndromes or anemias. Despite all the advances in patient care and the development of new therapies, some of these malignancies remain incurable, mainly due to resistance and refractoriness to treatment. Therefore, there is an unmet clinical need to identify new biomarkers and potential therapeutic targets that play a role in treatment resistance and contribute to the poor outcomes of these tumors. RNA-binding proteins (RBPs) are a diverse class of proteins that interact with transcripts and noncoding RNAs and are involved in every step of the post-transcriptional processing of transcripts. Dysregulation of RBPs has been associated with the development of hematological malignancies, making them potential valuable biomarkers and potential therapeutic targets. Although a number of dysregulated RBPs have been identified in hematological malignancies, there is a critical need to understand the biology underlying their contribution to pathology, such as the spatiotemporal context and molecular mechanisms involved. In this review, we emphasize the importance of deciphering the regulatory mechanisms of RBPs to pinpoint novel therapeutic targets that could drive or contribute to hematological malignancy biology.

BaRTv2: a highly resolved barley reference transcriptome for accurate transcript-specific RNA-seq quantification

Accurate characterisation of splice junctions (SJs) as well as transcription start and end sites in reference transcriptomes allows precise quantification of transcripts from RNA-seq data, and enables detailed investigations of transcriptional and post-transcriptional regulation. Using novel computational methods and a combination of PacBio Iso-seq and Illumina short-read sequences from 20 diverse tissues and conditions, we generated a comprehensive and highly resolved barley reference transcript dataset from the European 2-row spring barley cultivar Barke (BaRTv2.18). Stringent and thorough filtering was carried out to maintain the quality and accuracy of the SJs and transcript start and end sites. BaRTv2.18 shows increased transcript diversity and completeness compared with an earlier version, BaRTv1.0. The accuracy of transcript level quantification, SJs and transcript start and end sites have been validated extensively using parallel technologies and analysis, including high-resolution reverse transcriptase-polymerase chain reaction and 5'-RACE. BaRTv2.18 contains 39 434 genes and 148 260 transcripts, representing the most comprehensive and resolved reference transcriptome in barley to date. It provides an important and high-quality resource for advanced transcriptomic analyses, including both transcriptional and post-transcriptional regulation, with exceptional resolution and precision.

Avian influenza viruses suppress innate immunity by inducing trans-transcriptional readthrough via SSU72

Innate immunity plays critical antiviral roles. The highly virulent avian influenza viruses (AIVs) H5N1, H7N9, and H5N6 can better escape host innate immune responses than the less virulent seasonal H1N1 virus. Here, we report a mechanism by which transcriptional readthrough (TRT)-mediated suppression of innate immunity occurs post AIV infection. By using cell lines, mouse lungs, and patient PBMCs, we showed that genes on the complementary strand (“trans” genes) influenced by TRT were involved in the disruption of host antiviral responses during AIV infection. The trans-TRT enhanced viral lethality, and TRT abolishment increased cell viability and STAT1/2 expression. The viral NS1 protein directly bound to SSU72, and degradation of SSU72 induced TRT. SSU72 overexpression reduced TRT and alleviated mouse lung injury. Our results suggest that AIVs infection induce TRT by reducing SSU72 expression, thereby impairing host immune responses, a molecular mechanism acting through the NS1-SSU72-trans-TRT-STAT1/2 axis. Thus, restoration of SSU72 expression might be a potential strategy for preventing AIV pandemics.

Colorectal cancer concurrent gene signature based on coherent patterns between genomic and transcriptional alterations

Background

The aim of the study was to enhance colorectal cancer prognostication by integrating single nucleotide polymorphism (SNP) and gene expression (GE) microarrays for genomic and transcriptional alteration detection; genes with concurrent gains and losses were used to develop a prognostic signature.

Methods

The discovery dataset comprised 32 Taiwanese colorectal cancer patients, of which 31 were assayed for GE and copy number variations (CNVs) with Illumina Human HT-12 BeadChip v4.0 and Omni 25 BeadChip v1.1. Concurrent gains and losses were declared if coherent manners were observed between GE and SNP arrays. Concurrent genes were also identified in The Cancer Genome Atlas Project (TCGA) as the secondary discovery dataset (n = 345).

Results

The “universal” concurrent genes, which were the combination of z-transformed correlation coefficients, contained 4022 genes. Candidate genes were evaluated within each of the 10 public domain microarray datasets, and 1655 (2000 probe sets) were prognostic in at least one study. Consensus across all datasets was used to build a risk predictive model, while distinct relapse-free/overall survival patterns between defined risk groups were observed among four out of five training datasets. The predictive accuracy of recurrence, metastasis, or death was between 61 and 86% (cross-validation area under the receiver operating characteristic (ROC) curve: 0.548-0.833) from five independent validation studies.

Conclusion

The colorectal cancer concurrent gene signature is prognostic in terms of recurrence, metastasis, or mortality among 1746 patients. Genes with coherent patterns between genomic and transcriptional contexts are more likely to provide prognostication for colorectal cancer.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-09627-9.

The role of RNA binding proteins in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is one of the leading causes of overall cancer deaths worldwide with limited therapeutic options. Due to the heterogeneity of HCC pathogenesis, the molecular mechanisms underlying HCC development are not fully understood. Emerging evidence indicates that RNA-binding proteins (RBPs) play a vital role throughout hepatocarcinogenesis. Thus, a deeper understanding of how RBPs contribute to HCC progression will provide new tools for early diagnosis and prognosis of this devastating disease. In this review, we summarize the tumor suppressive and oncogenic roles of RBPs and their roles in hepatocarcinogenesis. The diagnostic and therapeutic potential of RBPs in HCC, including their limitations, are also discussed.

Who let the DoGs out? - biogenesis of stress-induced readthrough transcripts

Readthrough transcription caused by inefficient 3'-end cleavage of nascent mRNAs has emerged as a hallmark of the mammalian cellular stress response and results in the production of long noncoding RNAs known as downstream-of-gene (DoG)-containing transcripts. DoGs arise from around 10% of human protein-coding genes and are retained in the nucleus. They are produced minutes after cell exposure to stress and can be detected hours after stress removal. However, their biogenesis and the role(s) that DoGs or their production play in the cellular stress response are incompletely understood. We discuss findings that implicate host and viral proteins in the mechanisms underlying DoG production, as well as the transcriptional landscapes that accompany DoG induction under different stress conditions.

Differential Effects on the Translation of Immune-Related Alternatively Polyadenylated mRNAs in Melanoma and T Cells by eIF4A Inhibition

Targeting the translation initiation complex eIF4F, which binds the 5' cap of mRNAs, is a promising anti-cancer approach. Silvestrol, a small molecule inhibitor of eIF4A, the RNA helicase component of eIF4F, inhibits the translation of the mRNA encoding the signal transducer and activator of transcription 1 (STAT1) transcription factor, which, in turn, reduces the transcription of the gene encoding one of the major immune checkpoint proteins, i.e., programmed death ligand-1 (PD-L1) in melanoma cells. A large proportion of human genes produce multiple mRNAs differing in their 3'-ends through the use of alternative polyadenylation (APA) sites, which, when located in alternative last exons, can generate protein isoforms, as in the STAT1 gene. Here, we provide evidence that the STAT1α, but not STAT1β protein isoform generated by APA, is required for silvestrol-dependent inhibition of PD-L1 expression in interferon-γ-treated melanoma cells. Using polysome profiling in activated T cells we find that, beyond STAT1, eIF4A inhibition downregulates the translation of some important immune-related mRNAs, such as the ones encoding TIM-3, LAG-3, IDO1, CD27 or CD137, but with little effect on the ones for BTLA and ADAR-1 and no effect on the ones encoding CTLA-4, PD-1 and CD40-L. We next apply RT-qPCR and 3'-seq (RNA-seq focused on mRNA 3' ends) on polysomal RNAs to analyze in a high throughput manner the effect of eIF4A inhibition on the translation of APA isoforms. We identify about 150 genes, including TIM-3, LAG-3, AHNAK and SEMA4D, for which silvestrol differentially inhibits the translation of APA isoforms in T cells. It is therefore crucial to consider 3'-end mRNA heterogeneity in the understanding of the anti-tumor activities of eIF4A inhibitors.

Rapid factor depletion highlights intricacies of nucleoplasmic RNA degradation

Turnover of nucleoplasmic transcripts by the mammalian multi-subunit RNA exosome is mediated by two adaptors: the Nuclear EXosome Targeting (NEXT) complex and the Poly(A) tail eXosome Targeting (PAXT) connection. Functional analyses of NEXT and PAXT have largely utilized long-term factor depletion strategies, facilitating the appearance of indirect phenotypes. Here, we rapidly deplete NEXT, PAXT and core exosome components, uncovering the direct consequences of their acute losses. Generally, proteome changes are sparse and largely dominated by co-depletion of other exosome and adaptor subunits, reflecting possible subcomplex compositions. While parallel high-resolution 3′ end sequencing of newly synthesized RNA confirms previously established factor specificities, it concomitantly demonstrates an inflation of long-term depletion datasets by secondary effects. Most strikingly, a general intron degradation phenotype, observed in long-term NEXT depletion samples, is undetectable upon short-term depletion, which instead emphasizes NEXT targeting of snoRNA-hosting introns. Further analysis of these introns uncovers an unusual mode of core exosome-independent RNA decay. Our study highlights the accumulation of RNAs as an indirect result of long-term decay factor depletion, which we speculate is, at least partly, due to the exhaustion of alternative RNA decay pathways.

Diverse functional elements in RNA predicted transcriptome-wide by orthogonal RNA structure probing

RNA molecules can fold into complex structures and interact with trans-acting factors to control their biology. Recent methods have been focused on developing novel tools to measure RNA structure transcriptome-wide, but their utility to study and predict RNA-protein interactions or RNA processing has been limited thus far. Here, we extend these studies with the first transcriptome-wide mapping method for cataloging RNA solvent accessibility, icLASER. By combining solvent accessibility (icLASER) with RNA flexibility (icSHAPE) data, we efficiently predict RNA-protein interactions transcriptome-wide and catalog RNA polyadenylation sites by RNA structure alone. These studies showcase the power of designing novel chemical approaches to studying RNA biology. Further, our study exemplifies merging complementary methods to measure RNA structure inside cells and its utility for predicting transcriptome-wide interactions that are critical for control of and regulation by RNA structure. We envision such approaches can be applied to studying different cell types or cells under varying conditions, using RNA structure and footprinting to characterize cellular interactions and processing involving RNA.

Disassembling the Nature of Capsid: Biochemical, Genetic, and Imaging Approaches to Assess HIV-1 Capsid Functions

The human immunodeficiency virus type 1 (HIV-1) capsid and its disassembly, or capsid uncoating, has remained an active area of study over the past several decades. Our understanding of the HIV-1 capsid as solely a protective shell has since shifted with discoveries linking it to other complex replication events. The interplay of the HIV-1 capsid with reverse transcription, nuclear import, and integration has led to an expansion of knowledge of capsid functionality. Coincident with advances in microscopy, cell, and biochemistry assays, several models of capsid disassembly have been proposed, in which it occurs in either the cytoplasmic, nuclear envelope, or nuclear regions of the cell. Here, we discuss how the understanding of the HIV-1 capsid has evolved and the key methods that made these discoveries possible.

Aberrant Post-Transcriptional Regulation of Protein Expression in the Development of Chronic Obstructive Pulmonary Disease

Chronic obstructive pulmonary disease (COPD) is an incurable and prevalent respiratory disorder that is characterized by chronic inflammation and emphysema. COPD is primarily caused by cigarette smoke (CS). CS alters numerous cellular processes, including the post-transcriptional regulation of mRNAs. The identification of RNA-binding proteins (RBPs), microRNAs (miRNAs), and long non-coding RNAs (lncRNAs) as main factors engaged in the regulation of RNA biology opens the door to understanding their role in coordinating physiological cellular processes. Dysregulation of post-transcriptional regulation by foreign particles in CS may lead to the development of diseases such as COPD. Here we review current knowledge about post-transcriptional events that may be involved in the pathogenesis of COPD.

Alternative polyadenylation: An enigma of transcript length variation in health and disease

Alternative polyadenylation (APA) is a molecular mechanism during a pre-mRNA processing that involves usage of more than one polyadenylation site (PA-site) generating transcripts of varying length from a single gene. The location of a PA-site affects transcript length and coding potential of an mRNA contributing to both mRNA and protein diversification. This variation in the transcript length affects mRNA stability and translation, mRNA subcellular and tissue localization, and protein function. APA is now considered as an important regulatory mechanism in the pathophysiology of human diseases. An important consequence of the changes in the length of 3'-untranslated region (UTR) from disease-induced APA is altered protein expression. Yet, the relationship between 3'-UTR length and protein expression remains a paradox in a majority of diseases. Here, we review occurrence of APA, mechanism of PA-site selection, and consequences of transcript length variation in different diseases. Emerging evidence reveals coordinated involvement of core RNA processing factors including poly(A) polymerases in the PA-site selection in diseases-associated APAs. Targeting such APA regulators will be therapeutically significant in combating drug resistance in cancer and other complex diseases. This article is categorized under: RNA Processing > 3' End Processing RNA in Disease and Development > RNA in Disease Translation > Regulation.

On the Cutting Edge: Regulation and Therapeutic Potential of the mRNA 3' End Nuclease

Cleavage of nascent transcripts is a fundamental process for eukaryotic mRNA maturation and for the production of different mRNA isoforms. In eukaryotes, cleavage of mRNA precursors by the highly conserved endonuclease CPSF73 is critical for mRNA stability, export from the nucleus, and translation. As an essential enzyme in the cell, CPSF73 surprisingly shows promise as a drug target for specific cancers and for protozoan parasites. In this review, we cover our current understanding of CPSF73 in cleavage and polyadenylation, histone pre-mRNA processing, and transcription termination. We discuss the potential of CPSF73 as a target for novel therapeutics and highlight further research into the regulation of CPSF73 that will be critical to understanding its role in cancer and other diseases.

Novel insights into the pervasive role of RNA structure in post-transcriptional regulation of gene expression in plants

RNA folding is an intrinsic property of RNA that serves a key role in every step of post-transcriptional regulation of gene expression, from RNA maturation to translation in plants. Recent developments of genome-wide RNA structure profiling methods have transformed research in this area enabling focus to shift from individual molecules to the study of tens of thousands of RNAs. Here, we provide a comprehensive review of recent advances in the field. We discuss these new insights of RNA structure functionality within the context of post-transcriptional regulation including mRNA maturation, translation, and RNA degradation in plants. Notably, we also provide an overview of how plants exhibit different RNA structures in response to environmental changes.

Emerging roles of alternative cleavage and polyadenylation (APA) in human disease

In the messenger RNA (mRNA) maturation process, the 3'-end of pre-mRNA is cleaved and a poly(A) sequence is added, this is an important determinant of mRNA stability and its cellular functions. More than 60%-70% of human genes have three or more polyadenylation (APA) sites and can be cleaved at different sites, generating mRNA transcripts of varying lengths. This phenomenon is termed as alternative cleavage and polyadenylation (APA) and it plays role in key biological processes like gene regulation, cell proliferation, senescence, and also in various human diseases. Loss of regulatory microRNA binding sites and interactions with RNA-binding proteins leading to APA are largely investigated in human diseases. However, the functions of the core APA machinery and related factors during disease conditions remain largely unknown. In this review, we discuss the roles of polyadenylation machinery in relation to brain disease, cardiac failure, pulmonary fibrosis, cancer, infectious conditions, and other human diseases. Collectively, we believe this review will be a useful avenue for understanding the emerging role of APA in the pathobiology of various human diseases.