3' end mRNA processing: molecular mechanisms and implications for health and disease

High-throughput discovery of regulatory effector domains in human RNA-binding proteins

RNA regulation plays an integral role in tuning gene expression and is controlled by thousands of RNA-binding proteins (RBPs). We develop and use a high-throughput recruitment assay (HT-RNA-Recruit) to identify regulatory domains within human RBPs by recruiting over 30,000 protein tiles from 367 RBPs to a reporter mRNA. We discover over 100 unique RNA-regulatory effectors in 86 distinct RBPs, presenting evidence that RBPs contain functionally separable domains that dictate their post-transcriptional control of gene expression, and identify some with unique activity at 5' or 3'UTRs. We identify some domains that downregulate gene expression both when recruited to DNA and RNA, and dissect their mechanisms of regulation. Finally, we build a synthetic RNA regulator that can stably maintain gene expression at desired levels that are predictable by a mathematical model. This work serves as a resource for human RNA-regulatory effectors and expands the synthetic repertoire of RNA-based genetic control tools.

Highlights

HT-RNA-Recruit identifies hundreds of RNA-regulatory effectors in human proteins.Recruitment to 5' and 3' UTRs identifies regulatory domains unique to each position.Some protein domains have both transcriptional and post-transcriptional regulatory activity.We develop a synthetic RNA regulator and a mathematical model to describe its behavior.

Disruption of PABPN1 phase separation by SNRPD2 drives colorectal cancer cell proliferation and migration through promoting alternative polyadenylation of CTNNBIP1

Generally shortened 3' UTR due to alternative polyadenylation (APA) is widely observed in cancer, but its regulation mechanisms for cancer are not well characterized. Here, with profiling of APA in colorectal cancer tissues and poly(A) signal editing, we firstly identified that the shortened 3' UTR of CTNNIBP1 in colorectal cancer promotes cell proliferation and migration. We found that liquid-liquid phase separation (LLPS) of PABPN1 is reduced albeit with higher expression in cancer, and the reduction of LLPS leads to the shortened 3' UTR of CTNNBIP1 and promotes cell proliferation and migration. Notably, the splicing factor SNRPD2 upregulated in colorectal cancer, can interact with glutamic-proline (EP) domain of PABPN1, and then disrupt LLPS of PABPN1, which attenuates the repression effect of PABPN1 on the proximal poly(A) sites. Our results firstly reveal a new regulation mechanism of APA by disruption of LLPS of PABPN1, suggesting that regulation of APA by interfering LLPS of 3' end processing factor may have the potential as a new way for the treatment of cancer.

Identifying human pre-mRNA cleavage and polyadenylation factors by genome-wide CRISPR screens using a dual fluorescence readthrough reporter

Messenger RNA precursors (pre-mRNA) generally undergo 3' end processing by cleavage and polyadenylation (CPA), which is specified by a polyadenylation site (PAS) and adjacent RNA sequences and regulated by a large variety of core and auxiliary CPA factors. To date, most of the human CPA factors have been discovered through biochemical and proteomic studies. However, genetic identification of the human CPA factors has been hampered by the lack of a reliable genome-wide screening method. We describe here a dual fluorescence readthrough reporter system with a PAS inserted between two fluorescent reporters. This system enables measurement of the efficiency of 3' end processing in living cells. Using this system in combination with a human genome-wide CRISPR/Cas9 library, we conducted a screen for CPA factors. The screens identified most components of the known core CPA complexes and other known CPA factors. The screens also identified CCNK/CDK12 as a potential core CPA factor, and RPRD1B as a CPA factor that binds RNA and regulates the release of RNA polymerase II at the 3' ends of genes. Thus, this dual fluorescence reporter coupled with CRISPR/Cas9 screens reliably identifies bona fide CPA factors and provides a platform for investigating the requirements for CPA in various contexts.

Roles for epithelial integrin α3β1 in regulation of the microenvironment during normal and pathological tissue remodeling

Integrin receptors for the extracellular matrix activate intracellular signaling pathways that are critical for tissue development, homeostasis, and regeneration/repair, and their loss or dysregulation contributes to many developmental defects and tissue pathologies. This review will focus on tissue remodeling roles for integrin α3β1, a receptor for laminins found in the basement membranes (BMs) that underlie epithelial cell layers. As a paradigm, we will discuss literature that supports a role for α3β1 in promoting ability of epidermal keratinocytes to modify their tissue microenvironment during skin development, wound healing, or tumorigenesis. Preclinical and clinical studies have shown that this role depends largely on ability of α3β1 to govern the keratinocyte's repertoire of secreted proteins, or the "secretome," including 1) matrix proteins and proteases involved in matrix remodeling and 2) paracrine-acting growth factors/cytokines that stimulate other cells with important tissue remodeling functions (e.g., endothelial cells, fibroblasts, inflammatory cells). Moreover, α3β1 signaling controls gene expression that helps epithelial cells carry out these functions, including genes that encode secreted matrix proteins, proteases, growth factors, or cytokines. We will review what is known about α3β1-dependent gene regulation through both transcription and posttranscriptional mRNA stability. Regarding the latter, we will discuss examples of α3β1-dependent alternative splicing (AS) or alternative polyadenylation (APA) that prevents inclusion of cis-acting mRNA sequences that would otherwise target the transcript for degradation via nonsense-mediated decay or destabilizing AU-rich elements (AREs) in the 3'-untranslated region (3'-UTR). Finally, we will discuss prospects and anticipated challenges of exploiting α3β1 as a clinical target for the treatment of cancer or wound healing.

Causes, functions, and therapeutic possibilities of RNA secondary structure ensembles and alternative states

RNA secondary structure plays essential roles in encoding RNA regulatory fate and function. Most RNAs populate ensembles of alternatively paired states and are continually unfolded and refolded by cellular processes. Measuring these structural ensembles and their contributions to cellular function has traditionally posed major challenges, but new methods and conceptual frameworks are beginning to fill this void. In this review, we provide a mechanism- and function-centric compendium of the roles of RNA secondary structural ensembles and minority states in regulating the RNA life cycle, from transcription to degradation. We further explore how dysregulation of RNA structural ensembles contributes to human disease and discuss the potential of drugging alternative RNA states to therapeutically modulate RNA activity. The emerging paradigm of RNA structural ensembles as central to RNA function provides a foundation for a deeper understanding of RNA biology and new therapeutic possibilities.

EFHC1 gene mutation profile of Turkish JME patients and its association with disease risk

OBJECTIVES

Juvenile myoclonic epilepsy (JME) is a common form of generalized epilepsy with an important genetic component. This cohort study aimed to examine the frequency of EFHC1 gene variants in Turkish JME patients and a healthy control group and evaluate the association between these mutations and disease risk.

METHODS

We screened 72 JME patients with a mean age of 31.8 ± 9.9 (20-65) years and 35 controls with a mean age of 29.1 ± 7.6 (17-50) years from southern Turkey using direct sequencing analyses.

RESULTS

EFCH1 single nucleotide variants were detected in 24 of 72 JME patients and 3 of 35 controls. The most common mutations were R182H in JME patients (p = 0.010) and 3'UTR in the control group (p < 0.001). The R182H mutation is a common variant in JME (95 % CI: 1.232-76.580, p = 0.031) and the 3'UTR mutation may be associated with lower risk of JME in the Turkish population (95 % CI: 13.89-166.67, p < 0.001).

SIGNIFICANCE

Our results indicate that EFHC1 gene variants carry a risk for JME and the 3'UTR variant may have a protective role against JME in the Turkish population. Screening for other genes is needed to further clarify the genetic inheritance of JME in Turkish patients.

Advances in Nonviral mRNA Delivery Materials and Their Application as Vaccines for Melanoma Therapy

Messenger RNA (mRNA) vaccines are promising platforms for cancer immunotherapy because of their potential to encode for a variety of tumor antigens, high tolerability, and capacity to induce strong antitumor immune responses. However, the clinical translation of mRNA cancer vaccines can be hindered by the inefficient delivery of mRNA in vivo. In this review, we provide an overview of mRNA cancer vaccines by discussing their utility in treating melanoma. Specifically, we begin our review by describing the barriers that can impede mRNA delivery to target cells. We then review native mRNA structure and discuss various modification methods shown to enhance mRNA stability and transfection. Next, we outline the advantages and challenges of three nonviral carrier platforms (lipid nanoparticles, polymeric nanoparticles, and lipopolyplexes) frequently used for mRNA delivery. Last, we summarize preclinical and clinical studies that have investigated nonviral mRNA vaccines for the treatment of melanoma. In writing this review, we aim to highlight innovative nonviral strategies designed to address mRNA delivery challenges while emphasizing the exciting potential of mRNA vaccines as next-generation therapies for the treatment of cancers.

The role and impact of alternative polyadenylation and miRNA regulation on the expression of the multidrug resistance-associated protein 1 (MRP-1/ABCC1) in epithelial ovarian cancer

The ATP-binding cassette transporter (ABCC1) is associated with poor survival and chemotherapy drug resistance in high grade serous ovarian cancer (HGSOC). The mechanisms driving ABCC1 expression are poorly understood. Alternative polyadenylation (APA) can give rise to ABCC1 mRNAs which differ only in the length of their 3'untranslated regions (3'UTRs) in a process known as 3'UTR-APA. Like other ABC transporters, shortening of the 3'UTR of ABCC1 through 3'UTR-APA would eliminate microRNA binding sites found within the longer 3'UTRs, hence eliminating miRNA regulation and altering gene expression. We found that the HGSOC cell lines Caov-3 and Ovcar-3 express higher levels of ABCC1 protein than normal cells. APA of ABCC1 occurs in all three cell lines resulting in mRNAs with both short and long 3'UTRs. In Ovcar-3, mRNAs with shorter 3'UTRs dominate resulting in a six-fold increase in protein expression. We were able to show that miR-185-5p and miR-326 both target the ABCC1 3'UTR. Hence, 3'UTR-APA should be considered as an important regulator of ABCC1 expression in HGSOC. Both HGSOC cell lines are cisplatin resistant, and we used erastin to induce ferroptosis, an alternative form of cell death. We showed that we could induce ferroptosis and sensitize the cisplatin resistant cells to cisplatin by using erastin. Knocking down ABCC1 resulted in decreased cell viability, but did not contribute to erastin induced ferroptosis.

Post-transcriptional control of haemostatic genes: mechanisms and emerging therapeutic concepts in thrombo-inflammatory disorders

The haemostatic system is pivotal to maintaining vascular integrity. Multiple components involved in blood coagulation have central functions in inflammation and immunity. A derailed haemostasis is common in prevalent pathologies such as sepsis, cardiovascular disorders, and lately, COVID-19. Physiological mechanisms limit the deleterious consequences of a hyperactivated haemostatic system through adaptive changes in gene expression. While this is mainly regulated at the level of transcription, co- and posttranscriptional mechanisms are increasingly perceived as central hubs governing multiple facets of the haemostatic system. This layer of regulation modulates the biogenesis of haemostatic components, for example in situations of increased turnover and demand. However, they can also be 'hijacked' in disease processes, thereby perpetuating and even causally entertaining associated pathologies. This review summarizes examples and emerging concepts that illustrate the importance of posttranscriptional mechanisms in haemostatic control and crosstalk with the immune system. It also discusses how such regulatory principles can be used to usher in new therapeutic concepts to combat global medical threats such as sepsis or cardiovascular disorders.

The U1 antisense morpholino oligonucleotide (AMO) disrupts U1 snRNP structure to promote intronic PCPA modification of pre-mRNAs

Functional depletion of the U1 small nuclear ribonucleoprotein (snRNP) with a 25 nt U1 AMO (antisense morpholino oligonucleotide) may lead to intronic premature cleavage and polyadenylation of thousands of genes, a phenomenon known as U1 snRNP telescripting; however, the underlying mechanism remains elusive. In this study, we demonstrated that U1 AMO could disrupt U1 snRNP structure both in vitro and in vivo, thereby affecting the U1 snRNP-RNAP polymerase II interaction. By performing chromatin immunoprecipitation sequencing for phosphorylation of Ser2 and Ser5 of the C-terminal domain of RPB1, the largest subunit of RNAP polymerase II, we showed that transcription elongation was disturbed upon U1 AMO treatment, with a particular high phosphorylation of Ser2 signal at intronic cryptic polyadenylation sites (PASs). In addition, we showed that core 3'processing factors CPSF/CstF are involved in the processing of intronic cryptic PAS. Their recruitment accumulated toward cryptic PASs upon U1 AMO treatment, as indicated by chromatin immunoprecipitation sequencing and individual-nucleotide resolution CrossLinking and ImmunoPrecipitation sequencing analysis. Conclusively, our data suggest that disruption of U1 snRNP structure mediated by U1 AMO provides a key for understanding the U1 telescripting mechanism.

The emerging roles of CFIm25 (NUDT21/CPSF5) in human biology and disease

The mammalian cleavage factor I subunit CFIm25 (NUDT21) binds to the UGUA sequences of precursor RNAs. Traditionally, CFIm25 is known to facilitate 3' end formation of pre-mRNAs resulting in the formation of polyadenylated transcripts. Recent studies suggest that CFIm25 may be involved in the cyclization and hence generation of circular RNAs (circRNAs) that contain UGUA motifs. These circRNAs act as competing endogenous RNAs (ceRNAs) that disrupt the ceRNA-miRNA-mRNA axis. Other emerging roles of CFIm25 include regulating both alternative splicing and alternative polyadenylation (APA). APA generates different sized transcripts that may code for different proteins, or more commonly transcripts that code for the same protein but differ in the length and sequence content of their 3' UTRs (3' UTR-APA). CFIm25 mediated global changes in 3' UTR-APA affect human physiology including spermatogenesis and the determination of cell fate. Deregulation of CFIm25 and changes in 3' UTR-APA have been implicated in several human diseases including cancer. In many cancers, CFIm25 acts as a tumor suppressor. However, there are some cancers where CFIm25 has the opposite effect. Alterations in CFIm25-driven 3' UTR-APA may also play a role in neural dysfunction and fibrosis. CFIm25 mediated 3' UTR-APA changes can be used to generate specific signatures that can be used as potential biomarkers in development and disease. Due to the emerging role of CFIm25 as a regulator of the aforementioned RNA processing events, modulation of CFIm25 levels may be a novel viable therapeutic approach. This article is categorized under: RNA Processing > 3' End Processing RNA in Disease and Development > RNA in Disease.

Prx5 of Cristaria plicata has antioxidant function and is regulated by Nrf2/ARE signaling pathway

Cristaria plicata is one of the more important freshwater pearl bivalves in China, which is susceptible to pathogen infection, and greatly impacts the ability of breeding pearls. Nrf2/ARE signaling pathway and its downstream target gene Prx5 have endogenous antioxidant functions to protect cells from oxidative damage. The full-length cDNA of Prx5 was cloned from C. Plicata, which was 1420 bp, encoding a total of 189 amino acids and had two conserved cysteine residues (Cys78 and Cys179). The amino acid sequence of CpPrx5 was highly similar to Prx5 of other species. Real-time fluorescence quantitative PCR showed that CpPrx5 was distributed in various tissues of mussels, and the highest expression was in hepatopancreas. The expression of CpPrx5 up-regulated in hepatopancreas and gills after LPS, PGN and Poly:I:C stimulation. The recombinant plasmid DE3-PGEX-4T-1-CpPrx5 was expressed in Escherichia coli BL21 and showed antioxidant activity. With the increase of CpPrx5 protein concentration, the superhelical form of DNA was protected. The expression of CpPrx5 was up-regulated after interference CpKeap1 and down-regulated after interference CpNrf2. Gel block assay showed that CpNrf2 and CpMafK proteins blocked CpPrx5 promoter. Subcellular localization showed that CpPrx5 was located in 293T nucleus and cytoplasm and CpMafK was located in 293T nucleus. GST-Pull down verified that CpMafK and CpPrx5 could bind in vitro. These results indicated that Prx5 had antioxidant function and could protects DNA from oxidative damage, and participated in transcriptional regulation by combining with the transcription factor MafK. In addition, MafK could combine with Nrf2 to regulate the downstream target gene Prx5.

Pathophysiological Role and Diagnostic Potential of R-Loops in Cancer and Beyond

R-loops are DNA-RNA hybrids that play multifunctional roles in gene regulation, including replication, transcription, transcription-replication collision, epigenetics, and preserving the integrity of the genome. The aberrant formation and accumulation of unscheduled R-loops can disrupt gene expression and damage DNA, thereby causing genome instability. Recent links between unscheduled R-loop accumulation and the abundance of proteins that modulate R-loop biogenesis have been associated with numerous human diseases, including various cancers. Although R-loops are not necessarily causative for all disease entities described to date, they can perpetuate and even exacerbate the initially disease-eliciting pathophysiology, making them structures of interest for molecular diagnostics. In this review, we discuss the (patho) physiological role of R-loops in health and disease, their surprising diagnostic potential, and state-of-the-art techniques for their detection.

Deciphering the impact of genetic variation on human polyadenylation using APARENT2

BACKGROUND

3'-end processing by cleavage and polyadenylation is an important and finely tuned regulatory process during mRNA maturation. Numerous genetic variants are known to cause or contribute to human disorders by disrupting the cis-regulatory code of polyadenylation signals. Yet, due to the complexity of this code, variant interpretation remains challenging.

RESULTS

We introduce a residual neural network model, APARENT2, that can infer 3'-cleavage and polyadenylation from DNA sequence more accurately than any previous model. This model generalizes to the case of alternative polyadenylation (APA) for a variable number of polyadenylation signals. We demonstrate APARENT2's performance on several variant datasets, including functional reporter data and human 3' aQTLs from GTEx. We apply neural network interpretation methods to gain insights into disrupted or protective higher-order features of polyadenylation. We fine-tune APARENT2 on human tissue-resolved transcriptomic data to elucidate tissue-specific variant effects. By combining APARENT2 with models of mRNA stability, we extend aQTL effect size predictions to the entire 3' untranslated region. Finally, we perform in silico saturation mutagenesis of all human polyadenylation signals and compare the predicted effects of [Formula: see text] million variants against gnomAD. While loss-of-function variants were generally selected against, we also find specific clinical conditions linked to gain-of-function mutations. For example, we detect an association between gain-of-function mutations in the 3'-end and autism spectrum disorder. To experimentally validate APARENT2's predictions, we assayed clinically relevant variants in multiple cell lines, including microglia-derived cells.

CONCLUSIONS

A sequence-to-function model based on deep residual learning enables accurate functional interpretation of genetic variants in polyadenylation signals and, when coupled with large human variation databases, elucidates the link between functional 3'-end mutations and human health.

The Bidirectional Link Between RNA Cleavage and Polyadenylation and Genome Stability: Recent Insights From a Systematic Screen

The integrity of the genome is governed by multiple processes to ensure optimal survival and to prevent the inheritance of deleterious traits. While significant progress has been made to characterize components involved in the DNA Damage Response (DDR), little is known about the interplay between RNA processing and the maintenance of genome stability. Here, we describe the emerging picture of an intricate bidirectional coupling between RNA processing and genome integrity in an integrative manner. By employing insights from a recent large-scale RNAi screening involving the depletion of more than 170 components that direct (alternative) polyadenylation, we provide evidence of bidirectional crosstalk between co-transcriptional RNA 3′end processing and the DDR in a manner that optimizes genomic integrity. We provide instructive examples illustrating the wiring between the two processes and show how perturbations at one end are either compensated by buffering mechanisms at the other end, or even propel the initial insult and thereby become disease-eliciting as evidenced by various disorders.

DNA Double-Strand Breaks as Pathogenic Lesions in Neurological Disorders

The damage and repair of DNA is a continuous process required to maintain genomic integrity. DNA double-strand breaks (DSBs) are the most lethal type of DNA damage and require timely repair by dedicated machinery. DSB repair is uniquely important to nondividing, post-mitotic cells of the central nervous system (CNS). These long-lived cells must rely on the intact genome for a lifetime while maintaining high metabolic activity. When these mechanisms fail, the loss of certain neuronal populations upset delicate neural networks required for higher cognition and disrupt vital motor functions. Mammalian cells engage with several different strategies to recognize and repair chromosomal DSBs based on the cellular context and cell cycle phase, including homologous recombination (HR)/homology-directed repair (HDR), microhomology-mediated end-joining (MMEJ), and the classic non-homologous end-joining (NHEJ). In addition to these repair pathways, a growing body of evidence has emphasized the importance of DNA damage response (DDR) signaling, and the involvement of heterogeneous nuclear ribonucleoprotein (hnRNP) family proteins in the repair of neuronal DSBs, many of which are linked to age-associated neurological disorders. In this review, we describe contemporary research characterizing the mechanistic roles of these non-canonical proteins in neuronal DSB repair, as well as their contributions to the etiopathogenesis of selected common neurological diseases.

Modulation of mRNA 3'-End Processing and Transcription Termination in Virus-Infected Cells

Eukaryotic mRNA 3´-end processing is a multi-step process beginning with pre-mRNA transcript cleavage followed by poly(A) tail addition. Closely coupled to transcription termination, 3´-end processing is a critical step in the regulation of gene expression, and disruption of 3´-end processing is known to affect mature mRNA levels. Various viral proteins interfere with the 3´-end processing machinery, causing read-through transcription and altered levels of mature transcripts through inhibition of cleavage and polyadenylation. Thus, disruption of 3´-end processing contributes to widespread host shutoff, including suppression of the antiviral response. Additionally, observed features of read-through transcripts such as decreased polyadenylation, nuclear retention, and decreased translation suggest that viruses may utilize these mechanisms to modulate host protein production and dominate cellular machinery. The degree to which the effects of read-through transcript production are harnessed by viruses and host cells remains unclear, but existing research highlights the importance of host 3´-end processing modulation during viral infection.

Implications of Poly(A) Tail Processing in Repeat Expansion Diseases

Repeat expansion diseases are a group of more than 40 disorders that affect mainly the nervous and/or muscular system and include myotonic dystrophies, Huntington’s disease, and fragile X syndrome. The mutation-driven expanded repeat tract occurs in specific genes and is composed of tri- to dodeca-nucleotide-long units. Mutant mRNA is a pathogenic factor or important contributor to the disease and has great potential as a therapeutic target. Although repeat expansion diseases are quite well known, there are limited studies concerning polyadenylation events for implicated transcripts that could have profound effects on transcript stability, localization, and translation efficiency. In this review, we briefly present polyadenylation and alternative polyadenylation (APA) mechanisms and discuss their role in the pathogenesis of selected diseases. We also discuss several methods for poly(A) tail measurement (both transcript-specific and transcriptome-wide analyses) and APA site identification—the further development and use of which may contribute to a better understanding of the correlation between APA events and repeat expansion diseases. Finally, we point out some future perspectives on the research into repeat expansion diseases, as well as APA studies.

Analysis Polyadenylation Signal Usage in Sus scrofa

RNA polyadenylation is an important step in the messenger RNA (mRNA) maturation process, and the first step is recognizing the polyadenylation signal (PAS). The PAS type and distribution is a key determinant of post-transcriptional mRNA modification and gene expression. However, little is known about PAS usage and alternative polyadenylation (APA) regulation in livestock species. Recently, sequencing technology has enabled the generation of a large amount of sequencing data revealing variation in poly(A) signals and APA regulation in Sus scrofa. We identified 62,491 polyadenylation signals in Sus scrofa using expressed sequence tag (EST) sequences combined with RNA-seq analysis. The composition and usage frequency of polyadenylation signal in Sus scrofa is similar with that of human and mouse. The most highly conserved polyadenylation signals are AAUAAA and AUUAAA, used for over 63.35% of genes. In addition, we also analyzed the U/GU-rich downstream sequence (DSE) element, located downstream of the cleavage site. Our results indicate that APA regulation was widely occurred in Sus scrofa, as in other organisms. Our result was useful for the accurate annotation of RNA 3' ends in Sus scrofa and the analysis of polyadenylation signal usage in Sus scrofa would give the new insights into the mechanisms of transcriptional regulation.

Architectural and functional details of CF IA proteins involved in yeast 3'-end pre-mRNA processing and its significance for eukaryotes: A concise review

In eukaryotes, maturation of pre-mRNA relies on its precise 3'-end processing. This processing involves co-transcriptional steps regulated by sequence elements and other proteins. Although, it holds tremendous importance, defect in the processing machinery will result in erroneous pre-mRNA maturation leading to defective translation. Remarkably, more than 20 proteins in humans and yeast share homology and execute this processing. The defects in this processing are associated with various diseases in humans. We shed light on the CF IA subunit of yeast Saccharomyces cerevisiae that contains four proteins (Pcf11, Clp1, Rna14 and Rna15) involved in this processing. Structural details of various domains of CF IA and their roles during 3'-end processing, like cleavage and polyadenylation at 3'-UTR of pre-mRNA and other cellular events are explained. Further, the chronological development and important discoveries associated with 3'-end processing are summarized. Moreover, the mammalian homologues of yeast CF IA proteins, along with their key roles are described. This knowledge would be helpful for better comprehension of the mechanism associated with this marvel; thus opening up vast avenues in this area.

mRNA - A game changer in regenerative medicine, cell-based therapy and reprogramming strategies

After thirty years of intensive research shaping and optimizing the technology, the approval of the first mRNA-based formulation by the EMA and FDA in order to stop the COVID-19 pandemic was a breakthrough in mRNA research. The astonishing success of these vaccines have brought the mRNA platform into the spotlight of the scientific community. The remarkable persistence of the groundwork is mainly attributed to the exceptional benefits of mRNA application, including the biological origin, immediate but transitory mechanism of action, non-integrative properties, safe and relatively simple manufacturing as well as the flexibility to produce any desired protein. Based on these advantages, a practical implementation of in vitro transcribed mRNA has been considered in most areas of medicine. In this review, we discuss the key preconditions for the rise of the mRNA in the medical field, including the unique structural and functional features of the mRNA molecule and its vehicles, which are crucial aspects for a production of potent mRNA-based therapeutics. Further, we focus on the utility of mRNA tools particularly in the scope of regenerative medicine, i.e. cell reprogramming approaches or manipulation strategies for targeted tissue restoration. Finally, we highlight the strong clinical potential but also the remaining hurdles to overcome for the mRNA-based regenerative therapy, which is only a few steps away from becoming a reality.

Regulation of closely juxtaposed proto-oncogene c-fms and HMGXB3 gene expression by mRNA 3' end polymorphism in breast cancer cells

Sense-antisense mRNA pairs generated by convergent transcription is a way of gene regulation. c-fms gene is closely juxtaposed to the HMGXB3 gene in the opposite orientation, in chromosome 5. The intergenic region (IR) between c-fms and HMGXB3 genes is 162 bp. We found that a small portion (∼4.18%) of HMGXB3 mRNA is transcribed further downstream, including the end of the c-fms gene generating antisense mRNA against c-fms mRNA. Similarly, a small portion (∼1.1%) of c-fms mRNA is transcribed further downstream, including the end of the HMGXB3 gene generating antisense mRNA against the HMGXB3 mRNA. Insertion of the strong poly(A) signal sequence in the IR results in decreased c-fms and HMGXB3 antisense mRNAs, resulting in up-regulation of both c-fms and HMGXB3 mRNA expression. miR-324-5p targets HMGXB3 mRNA 3' UTR, and as a result, regulates c-fms mRNA expression. HuR stabilizes c-fms mRNA, and as a result, down-regulates HMGXB3 mRNA expression. UALCAN analysis indicates that the expression pattern between c-fms and HMGXB3 proteins are opposite in vivo in breast cancer tissues. Together, our results indicate that the mRNA encoded by the HMGXB3 gene can influence the expression of adjacent c-fms mRNA, or vice versa.

Lipofection with Synthetic mRNA as a Simple Method for T-Cell Immunomonitoring

The quantification of T-cell immune responses is crucial for the monitoring of natural and treatment-induced immunity, as well as for the validation of new immunotherapeutic approaches. The present study presents a simple method based on lipofection of synthetic mRNA in mononuclear cells as a method to determine in vitro T-cell responses. We compared several commercially available transfection reagents for their potential to transfect mRNA into human peripheral blood mononuclear cells and murine splenocytes. We also investigated the impact of RNA modifications in improving this method. Our results demonstrate that antigen-specific T-cell immunomonitoring can be easily and quickly performed by simple lipofection of antigen-coding mRNA in complex immune cell populations. Thus, our work discloses a convenient solution for the in vitro monitoring of natural or therapy-induced T-cell immune responses.

Modulation of alternative cleavage and polyadenylation events by dCas9-mediated CRISPRpas

The CRISPR/Cas9 technology is revolutionizing genomic engineering. The high efficiency and selectivity of the system have inspired the development of various derived tools for gene regulation at different levels, such as transcriptional activation or inhibition, epigenetic modification, splicing, and base editing. Cleavage and polyadenylation (CPA) is an essential 3' end maturation step for almost all eukaryotic mRNAs. CPA is tightly coupled with transcriptional termination, and its activity impacts gene expression. Over half of all human genes display alternative polyadenylation (APA), where multiple cleavage and polyadenylation sites (PASs) lead to mRNA isoforms with variable termini. APA isoforms often have distinct metabolisms, and their relative abundance can change drastically in different cells. Here, we describe a method based on delivering a catalytically dead Cas9 (dCas9) to genomic regions nears the PAS, which alters APA site usage in 3'UTRs or introns. This method, named CRISPRpas, allows investigators to examine functional significance of APA isoforms of individual genes. We also describe using the bioinformatics program APAlyzer to examine APA events of interest with RNA-seq data.

Dynamics of alternative splicing during somatic cell reprogramming reveals functions for RNA-binding proteins CPSF3, hnRNP UL1, and TIA1

BACKGROUND

Somatic cell reprogramming is the process that allows differentiated cells to revert to a pluripotent state. In contrast to the extensively studied rewiring of epigenetic and transcriptional programs required for reprogramming, the dynamics of post-transcriptional changes and their associated regulatory mechanisms remain poorly understood. Here we study the dynamics of alternative splicing changes occurring during efficient reprogramming of mouse B cells into induced pluripotent stem (iPS) cells and compare them to those occurring during reprogramming of mouse embryonic fibroblasts.

RESULTS

We observe a significant overlap between alternative splicing changes detected in the two reprogramming systems, which are generally uncoupled from changes in transcriptional levels. Correlation between gene expression of potential regulators and specific clusters of alternative splicing changes enables the identification and subsequent validation of CPSF3 and hnRNP UL1 as facilitators, and TIA1 as repressor of mouse embryonic fibroblasts reprogramming. We further find that these RNA-binding proteins control partially overlapping programs of splicing regulation, involving genes relevant for developmental and morphogenetic processes.

CONCLUSIONS

Our results reveal common programs of splicing regulation during reprogramming of different cell types and identify three novel regulators of this process and their targets.

TRENDseq-A highly multiplexed high throughput RNA 3' end sequencing for mapping alternative polyadenylation

Alternative polyadenylation (APA) is a widespread and highly dynamic mechanism of gene regulation. It affects more than 70% of all genes, resulting in transcript isoforms with distinct 3' end termini. APA thereby considerably expands the diversity of the transcriptome 3' end (TREND). This leads to mRNA isoforms with profoundly different physiological effects, by affecting protein output, production of distinct protein isoforms, or modulating protein localization. APA is globally regulated in various conditions, including developmental and adaptive programs. Since perturbations of APA can disrupt biological processes, ultimately resulting in most devastating disorders, querying the APA landscape is crucial to decipher underlying mechanisms, resulting consequences and potential diagnostic and therapeutic implications. Here we provide a detailed step-by-step protocol for TRENDseq, a method for transcriptome-wide high-throughput sequencing of polyadenylated RNA 3' ends in a highly multiplexed fashion. TRENDseq exploits linear amplification of the starting material to improve sensitivity while significantly reducing the amount of input material. It thereby represents a powerful tool to study APA in numerous experimental set-ups and/or limited human samples in a highly multiplexed and reproducible manner.

The Detection and Bioinformatic Analysis of Alternative 3' UTR Isoforms as Potential Cancer Biomarkers

Alternative transcript cleavage and polyadenylation is linked to cancer cell transformation, proliferation and outcome. This has led researchers to develop methods to detect and bioinformatically analyse alternative polyadenylation as potential cancer biomarkers. If incorporated into standard prognostic measures such as gene expression and clinical parameters, these could advance cancer prognostic testing and possibly guide therapy. In this review, we focus on the existing methodologies, both experimental and computational, that have been applied to support the use of alternative polyadenylation as cancer biomarkers.

Systematic refinement of gene annotations by parsing mRNA 3' end sequencing datasets

Alternative cleavage and polyadenylation generates mRNA 3' isoforms in a cell type-specific manner. Due to finite available RNA sequencing data of organisms with vast cell type complexity, currently available gene annotation resources are incomplete, which poses significant challenges to the comprehensive interpretation and quantification of transcriptomes. In this chapter, we introduce 3'GAmES, a stand-alone computational pipeline for the identification and quantification of novel mRNA 3'end isoforms from 3'mRNA sequencing data. 3'GAmES expands available repositories and improves comprehensive gene-tag counting by cost-effective 3' mRNA sequencing, faithfully mirroring whole-transcriptome RNAseq measurements. By employing R and bash shell scripts (assembled in a Singularity container) 3'GAmES systematically augments cell type-specific 3' ends of RNA polymerase II transcripts and increases the sensitivity of quantitative gene expression profiling by 3' mRNA sequencing. Public access: https://github.com/AmeresLab/3-GAmES.git.

Understanding RNA Binding by the Nonclassical Zinc Finger Protein CPSF30, a Key Factor in Polyadenylation during Pre-mRNA Processing

Cleavage and polyadenylation specificity factor 30 (CPSF30) is a zinc finger protein that regulates pre-mRNA processing. CPSF30 contains five CCCH domains and one CCHC domain and recognizes two conserved 3' pre-mRNA sequences: an AU hexamer and a U-rich motif. AU hexamer motifs are common in pre-mRNAs and are typically defined as AAUAAA. Variations within the AAUAAA hexamer occur in certain pre-mRNAs and can affect polyadenylation efficiency or be linked to diseases. The effects of disease-related variations on CPSF30/pre-mRNA binding were determined using a construct of CPSF30 that contains just the five CCCH domains (CPSF30-5F). Bioinformatics was utilized to identify the variability within the AU hexamer sequence in pre-mRNAs. The effects of this sequence variability on CPSF30-5F/RNA binding affinities were measured. Bases at positions 1, 2, 4, and 5 within the AU hexamer were found to be important for RNA binding. Bioinformatics revealed that the three bases flanking the AU hexamer at the 5' and 3' ends are twice as likely to be adenine or uracil as guanine and cytosine. The presence of A and U residues in these flanking regions was determined to promote higher-affinity CPSF30-5F/RNA binding than G and C residues. The addition of the zinc knuckle domain to CPSF30-5F (CPSF30-FL) restored binding to AU hexamer variants. This restoration of binding is connected to the presence of a U-rich sequence within the pre-mRNA to which the zinc knuckle binds. A mechanism of differential RNA binding by CPSF30, modulated by accessibility of the two RNA binding sites, is proposed.

Targeting RNA-binding proteins in acute and chronic leukemia

RNA-binding proteins (RBPs) play a crucial role in cellular physiology by regulating RNA processing, translation, and turnover. In neoplasms, RBP support of cancer-relevant expression of alternatively spliced, modified, and stabilized mRNA transcripts is essential to self-renewal, proliferation, and adaptation to stress. In this review, we assess the impact of key families of RBPs in leukemogenesis, review progress in targeting those proteins with small molecules, and discuss how multilevel composition of posttranscriptional regulation of gene expression could be used for potential therapies in acute and chronic leukemia.

Hypothalamic Norepinephrine Concentration and Heart Mass in Hypertensive ISIAH Rats Are Associated with a Genetic Locus on Chromosome 18

The relationship between activation of the sympathetic nervous system and cardiac hypertrophy has long been known. However, the molecular genetic basis of this association is poorly understood. Given the known role of hypothalamic norepinephrine in the activation of the sympathetic nervous system, the aim of the work was to carry out genetic mapping using Quantitative Trait Loci (QTL) analysis and determine the loci associated both with an increase in the concentration of norepinephrine in the hypothalamus and with an increase in heart mass in Inherited Stress-Induced Arterial Hypertension (ISIAH) rats simulating the stress-sensitive form of arterial hypertension. The work describes a genetic locus on chromosome 18, in which there are genes that control the development of cardiac hypertrophy associated with an increase in the concentration of norepinephrine in the hypothalamus, i.e., genes involved in enhanced sympathetic myocardial stimulation. No association of this locus with the blood pressure was found. Taking into consideration previously obtained results, it was concluded that the contribution to the development of heart hypertrophy in the ISIAH rats is controlled by different genetic loci, one of which is associated with the concentration of norepinephrine in the hypothalamus (on chromosome 18) and the other is associated with high blood pressure (on chromosome 1). Nucleotide substitutions that may be involved in the formation or absence of association with blood pressure in different rat strains are discussed.

TREND-DB-a transcriptome-wide atlas of the dynamic landscape of alternative polyadenylation

Alternative polyadenylation (APA) profoundly expands the transcriptome complexity. Perturbations of APA can disrupt biological processes, ultimately resulting in devastating disorders. A major challenge in identifying mechanisms and consequences of APA (and its perturbations) lies in the complexity of RNA 3′ end processing, involving poorly conserved RNA motifs and multi-component complexes consisting of far more than 50 proteins. This is further complicated in that RNA 3′ end maturation is closely linked to transcription, RNA processing and even epigenetic (histone/DNA/RNA) modifications. Here, we present TREND-DB (http://shiny.imbei.uni-mainz.de:3838/trend-db), a resource cataloging the dynamic landscape of APA after depletion of >170 proteins involved in various facets of transcriptional, co- and post-transcriptional gene regulation, epigenetic modifications and further processes. TREND-DB visualizes the dynamics of transcriptome 3′ end diversification (TREND) in a highly interactive manner; it provides a global APA network map and allows interrogating genes affected by specific APA-regulators and vice versa. It also permits condition-specific functional enrichment analyses of APA-affected genes, which suggest wide biological and clinical relevance across all RNAi conditions. The implementation of the UCSC Genome Browser provides additional customizable layers of gene regulation accounting for individual transcript isoforms (e.g. epigenetics, miRNA-binding sites and RNA-binding proteins). TREND-DB thereby fosters disentangling the role of APA for various biological programs, including potential disease mechanisms, and helps identify their diagnostic and therapeutic potential.

Comparative accuracies of genetic values predicted for economically important milk traits, genome-wide association, and linkage disequilibrium patterns of Canadian Holstein cows

Genomic selection methodologies and genome-wide association studies use powerful statistical procedures that correlate large amounts of high-density SNP genotypes and phenotypic data. Actual 305-d milk (MY), fat (FY), and protein (PY) yield data on 695 cows and 76,355 genotyping-by-sequencing-generated SNP marker genotypes from Canadian Holstein dairy cows were used to characterize linkage disequilibrium (LD) structure of Canadian Holstein cows. Also, the comparison of pedigree-based BLUP, genomic BLUP (GBLUP), and Bayesian (BayesB) statistical methods in the genomic selection methodologies and the comparison of Bayesian ridge regression and BayesB statistical methods in the genome-wide association studies were carried out for MY, FY, and PY. Results from LD analysis revealed that as marker distance decreases, LD increases through chromosomes. However, unexpected high peaks in LD were observed between marker pairs with larger marker distances on all chromosomes. The GBLUP and BayesB models resulted in similar heritability estimates through 10-fold cross-validation for MY and PY; however, the GBLUP model resulted in higher heritability estimates than BayesB model for FY. The predictive ability of GBLUP model was significantly lower than that of BayesB for MY, FY, and PY. Association analyses indicated that 28 high-effect markers and markers on Bos taurus autosome 14 located within 6 genes (DOP1B, TONSL, CPSF1, ADCK5, PARP10, and GRINA) associated significantly with FY.

Aptamer-Mediated Control of Polyadenylation for Gene Expression Regulation in Mammalian Cells

Small aptamer-based regulatory devices can be designed to control a range of RNA-dependent cellular processes and emerged as promising tools for fine-tuning gene expression in synthetic biology. Here, we design a conceptually new riboswitch device that allows for the conditional regulation of polyadenylation. By making use of ligand-induced sequence occlusion, the system efficiently controls the accessibility of the eukaryotic polyadenylation signal. Undesirable 3'-extended read-through products are counteracted by the downstream insertion of a microRNA target site. We demonstrate the modularity of the system with regard to sensor aptamers and polyadenylation signals used and combine the newly designed riboswitch with well-known aptazymes to yield superior composite systems. In addition, we show that the switches can be used to control alternative polyadenylation. The presented genetic switches require very little coding space and can be easily optimized by rational adjustments of the thermodynamic stability. The polyadenylation riboswitch extends the repertoire of RNA-based regulators and opens new possibilities for the generation of complex synthetic circuits.

USP6NL mediated by LINC00689/miR-142-3p promotes the development of triple-negative breast cancer

Background

Triple-negative breast cancer (TNBC), in part because of the high metastasis rate, is one of the most prevalent causes of malignancy-related mortality globally. Ubiquitin specific peptidase 6 N-terminal like (USP6NL) has been unmasked to be implicated in some human cancers. However, the precise biological function of USP6NL in TNBC has not been defined.

Methods

RNA expression was examined by real-time quantitative PCR (RT-qPCR), while USP6NL protein level was tested through western blot. Besides, cell proliferation was assessed by using colony formation assay, whereas cell apoptosis estimated by flow cytometry analysis, JC-1 assay and TUNEL assay. Transwell assays were adopted to detect the migration and invasion of indicated TNBC cells. Immunofluorescence (IF) assay evaluated epithelial-mesenchymal transitions (EMT) progress in TNBC. Further, RNA immunoprecipitation (RIP), RNA pull down and luciferase reporter assays were implemented for measuring the mutual interplay among USP6NL, miR-142-3p and long intergenic non-protein coding RNA 689 (LINC00689).

Results

Elevated USP6NL level was uncovered in TNBC cells. RNA interference-mediated knockdown of USP6NL inhibited TNBC cell growth, motility and EMT. Further, USP6NL was proved as the target of a tumor-inhibitor miR-142-3p, and LINC00689 augmented USP6NL expression by absorbing miR-142-3p. Importantly, miR-142-3p deficiency or USP6NL overexpression fully abolished the inhibitory effect of LINC00689 silence on TNBC cellular behaviors.

Conclusion

All data revealed the important role of USP6NL/LINC00689/miR-142-3p signaling in TNBC. The findings might provide a new and promising therapeutic biomarker for treating patients with TNBC.

Systematic identification of functional SNPs interrupting 3'UTR polyadenylation signals

Alternative polyadenylation (APA) is emerging as a widespread regulatory layer since the majority of human protein-coding genes contain several polyadenylation (p(A)) sites in their 3’UTRs. By generating isoforms with different 3’UTR length, APA potentially affects mRNA stability, translation efficiency, nuclear export, and cellular localization. Polyadenylation sites are regulated by adjacent RNA cis-regulatory elements, the principals among them are the polyadenylation signal (PAS) AAUAAA and its main variant AUUAAA, typically located ~20-nt upstream of the p(A) site. Mutations in PAS and other auxiliary poly(A) cis-elements in the 3’UTR of several genes have been shown to cause human Mendelian diseases, and to date, only a few common SNPs that regulate APA were associated with complex diseases. Here, we systematically searched for SNPs that affect gene expression and human traits by modulation of 3’UTR APA. First, focusing on the variants most likely to exert the strongest effect, we identified 2,305 SNPs that interrupt the canonical PAS or its main variant. Implementing pA-QTL tests using GTEx RNA-seq data, we identified 330 PAS SNPs (called PAS pA-QTLs) that were significantly associated with the usage of their p(A) site. As expected, PAS-interrupting alleles were mostly linked with decreased cleavage at their p(A) site and the consequential 3’UTR lengthening. However, interestingly, in ~10% of the cases, the PAS-interrupting allele was associated with increased usage of an upstream p(A) site and 3’UTR shortening. As an indication of the functional effects of these PAS pA-QTLs on gene expression and complex human traits, we observed for few dozens of them marked colocalization with eQTL and/or GWAS signals. The PAS-interrupting alleles linked with 3’UTR lengthening were also strongly associated with decreased gene expression, indicating that shorter isoforms generated by APA are generally more stable than longer ones. Last, we carried out an extended, genome-wide analysis of 3’UTR variants and detected thousands of additional pA-QTLs having weaker effects compared to the PAS pA-QTLs.

mRNA molecules that encode for proteins end with a long stretch of adenosines, called poly(A) tail. The poly(A) tail contributes to the stability of the mRNA molecules, their translation to proteins and their import from the nucleus to the cytoplasm. The process of adding this tail to the mRNAs is called polyadenylation, and the termination site on the mRNAs at which the poly(A) tail is added is called the poly(A) site. In recent years it became evident that the vast majority of mRNAs of human genes contain several alternative poly(A) sites and their usage generates different mRNA isoforms that differ in their stability and translation efficiency. Therefore, alternative polyadenylation (APA) is emerging as a novel and important, yet underexplored, mechanism that regulate gene expression. The choice between alternative p(A) sites in an mRNA molecule is regulated by regulatory sequences located within a region in the mRNA called the 3’ untranslated region (3’UTR). A major challenge in present human genetics research is to understand how common genetic variants affect individuals’ health. In our study, we systematically identified dozens of genetic variants that affect the choice between alternative p(A) sites and demonstrated that by that, these variants influence the expression level of the target genes. Our results help to illuminate a novel mechanism by which genetic variants that are common in the population affect different traits including our risk for developing diseases.

An Improved Anticancer Drug-Response Prediction Based on an Ensemble Method Integrating Matrix Completion and Ridge Regression

In this study, we proposed an ensemble learning method, simultaneously integrating a low-rank matrix completion model and a ridge regression model to predict anticancer drug response on cancer cell lines. The model was applied to two benchmark datasets, including the Cancer Cell Line Encyclopedia (CCLE) and the Genomics of Drug Sensitivity in Cancer (GDSC). As previous studies suggest, the dual-layer integrated cell line-drug network model was one of the best models by far and outperformed most state-of-the-art models. Thus, we performed a head-to-head comparison between the dual-layer integrated cell line-drug network model and our model by a 10-fold crossvalidation study. For the CCLE dataset, our model has a higher Pearson correlation coefficient between predicted and observed drug responses than that of the dual-layer integrated cell line-drug network model in 18 out of 23 drugs. For the GDSC dataset, our model is better in 26 out of 28 drugs in the phosphatidylinositol 3-kinase (PI3K) pathway and 26 out of 30 drugs in the extracellular signal-regulated kinase (ERK) signaling pathway, respectively. Based on the prediction results, we carried out two types of case studies, which further verified the effectiveness of the proposed model on the drug-response prediction. In addition, our model is more biologically interpretable than the compared method, since it explicitly outputs the genes involved in the prediction, which are enriched in functions, like transcription, Src homology 2/3 (SH2/3) domain, cell cycle, ATP binding, and zinc finger.

DeepPASTA: deep neural network based polyadenylation site analysis

MOTIVATION

Alternative polyadenylation (polyA) sites near the 3' end of a pre-mRNA create multiple mRNA transcripts with different 3' untranslated regions (3' UTRs). The sequence elements of a 3' UTR are essential for many biological activities such as mRNA stability, sub-cellular localization, protein translation, protein binding and translation efficiency. Moreover, numerous studies in the literature have reported the correlation between diseases and the shortening (or lengthening) of 3' UTRs. As alternative polyA sites are common in mammalian genes, several machine learning tools have been published for predicting polyA sites from sequence data. These tools either consider limited sequence features or use relatively old algorithms for polyA site prediction. Moreover, none of the previous tools consider RNA secondary structures as a feature to predict polyA sites.

RESULTS

In this paper, we propose a new deep learning model, called DeepPASTA, for predicting polyA sites from both sequence and RNA secondary structure data. The model is then extended to predict tissue-specific polyA sites. Moreover, the tool can predict the most dominant (i.e. frequently used) polyA site of a gene in a specific tissue and relative dominance when two polyA sites of the same gene are given. Our extensive experiments demonstrate that DeepPASTA signisficantly outperforms the existing tools for polyA site prediction and tissue-specific relative and absolute dominant polyA site prediction.

AVAILABILITY AND IMPLEMENTATION

https://github.com/arefeen/DeepPASTA.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Emerging Roles of RNA 3'-end Cleavage and Polyadenylation in Pathogenesis, Diagnosis and Therapy of Human Disorders

A crucial feature of gene expression involves RNA processing to produce 3′ ends through a process termed 3′ end cleavage and polyadenylation (CPA). This ensures the nascent RNA molecule can exit the nucleus and be translated to ultimately give rise to a protein which can execute a function. Further, alternative polyadenylation (APA) can produce distinct transcript isoforms, profoundly expanding the complexity of the transcriptome. CPA is carried out by multi-component protein complexes interacting with multiple RNA motifs and is tightly coupled to transcription, other steps of RNA processing, and even epigenetic modifications. CPA and APA contribute to the maintenance of a multitude of diverse physiological processes. It is therefore not surprising that disruptions of CPA and APA can lead to devastating disorders. Here, we review potential CPA and APA mechanisms involving both loss and gain of function that can have tremendous impacts on health and disease. Ultimately we highlight the emerging diagnostic and therapeutic potential CPA and APA offer.

Aire Gene Influences the Length of the 3' UTR of mRNAs in Medullary Thymic Epithelial Cells

Aire is a transcriptional controller in medullary thymic epithelial cells (mTECs) modulating a set of peripheral tissue antigens (PTAs) and non-PTA mRNAs as well as miRNAs. Even miRNAs exerting posttranscriptional control of mRNAs in mTECs, the composition of miRNA-mRNA networks may differ. Under reduction in Aire expression, networks exhibited greater miRNA diversity controlling mRNAs. Variations in the number of 3'UTR binding sites of Aire-dependent mRNAs may represent a crucial factor that influence the miRNA interaction. To test this hypothesis, we analyzed through bioinformatics the length of 3'UTRs of a large set of Aire-dependent mRNAs. The data were obtained from existing RNA-seq of mTECs of wild type or Aire-knockout (KO) mice. We used computational algorithms as FASTQC, STAR and HTSEQ for sequence alignment and counting reads, DESEQ2 for the differential expression, 3USS for the alternative 3'UTRs and TAPAS for the alternative polyadenylation sites. We identified 152 differentially expressed mRNAs between these samples comprising those that encode PTAs as well as transcription regulators. In Aire KO mTECs, most of these mRNAs featured an increase in the length of their 3'UTRs originating additional miRNA binding sites and new miRNA controllers. Results from the in silico analysis were statistically significant and the predicted miRNA-mRNA interactions were thermodynamically stable. Even with no in vivo or in vitro experiments, they were adequate to show that lack of Aire in mTECs might favor the downregulation of PTA mRNAs and transcription regulators via miRNA control. This could unbalance the overall transcriptional activity in mTECs and thus the self-representation.

Single-Nucleotide Polymorphisms (SNPs) Both Associated with Hypertension and Contributing to Accelerated-Senescence Traits in OXYS Rats

Aging is a major risk factor of numerous human diseases. Adverse genetic variants may contribute to multiple manifestations of aging and increase the number of comorbid conditions. There is evidence of links between hypertension and age-related diseases, although the genetic relationships are insufficiently studied. Here, we investigated the contribution of hypertension to the development of accelerated-senescence syndrome in OXYS rats. We compared transcriptome sequences of the prefrontal cortex, hippocampus, and retina of OXYS rats with the genotypes of 45 rat strains and substrains (which include models with hypertension) to find single-nucleotide polymorphisms (SNPs) both associated with hypertension and possibly contributing to the development of age-related diseases. A total of 725 polymorphisms were common between OXYS rats and one or more hypertensive rat strains/substrains being analyzed. Multidimensional scaling detected significant similarities between OXYS and ISIAH rat genotypes and significant differences between these strains and the other hypertensive rat strains/substrains. Nonetheless, similar sets of SNPs produce a different phenotype in OXYS and ISIAH rats depending on hypertension severity. We identified 13 SNPs causing nonsynonymous amino-acid substitutions having a deleterious effect on the structure or function of the corresponding proteins and four SNPs leading to functionally significant structural rearrangements of transcripts in OXYS rats. Among them, SNPs in genes Ephx1, Pla2r1, and Ccdc28b were identified as candidates responsible for the concomitant manifestation of hypertension and signs of accelerated aging in OXYS rats.

PARPs and ADP-ribosylation in RNA biology: from RNA expression and processing to protein translation and proteostasis

In this review, Kim et al. discuss the importance of PARP family members and ADPRylation in gene regulation, mRNA processing, and protein abundance.

ADP-ribosylation (ADPRylation) is a posttranslational modification of proteins discovered nearly six decades ago, but many important questions remain regarding its molecular functions and biological roles, as well as the activity of the ADP-ribose (ADPR) transferase enzymes (PARP family members) that catalyze it. Growing evidence indicates that PARP-mediated ADPRylation events are key regulators of the protein biosynthetic pathway, leading from rDNA transcription and ribosome biogenesis to mRNA synthesis, processing, and translation. In this review we describe the role of PARP proteins and ADPRylation in all facets of this pathway. PARP-1 and its enzymatic activity are key regulators of rDNA transcription, which is a critical step in ribosome biogenesis. An emerging role of PARPs in alternative splicing of mRNAs, as well as direct ADPRylation of mRNAs, highlight the role of PARP members in RNA processing. Furthermore, PARP activity, stimulated by cellular stresses, such as viral infections and ER stress, leads to the regulation of mRNA stability and protein synthesis through posttranscriptional mechanisms. Dysregulation of PARP activity in these processes can promote disease states. Collectively, these results highlight the importance of PARP family members and ADPRylation in gene regulation, mRNA processing, and protein abundance. Future studies in these areas will yield new insights into the fundamental mechanisms and a broader utility for PARP-targeted therapeutic agents.

CPSF3-dependent pre-mRNA processing as a druggable node in AML and Ewing's sarcoma

The post-genomic era has seen many advances in our understanding of cancer pathways, yet resistance and tumor heterogeneity necessitate multiple approaches to target even monogenic tumors. Here, we combine phenotypic screening with chemical genetics to identify pre-messenger RNA endonuclease cleavage and polyadenylation specificity factor 3 (CPSF3) as the target of JTE-607, a small molecule with previously unknown target. We show that CPSF3 represents a synthetic lethal node in a subset of acute myeloid leukemia (AML) and Ewing's sarcoma cancer cell lines. Inhibition of CPSF3 by JTE-607 alters expression of known downstream effectors in AML and Ewing's sarcoma lines, upregulates apoptosis and causes tumor-selective stasis in mouse xenografts. Mechanistically, it prevents the release of newly synthesized pre-mRNAs, resulting in read-through transcription and the formation of DNA-RNA hybrid R-loop structures. This study implicates pre-mRNA processing, and specifically CPSF3, as a druggable target providing an avenue to therapeutic intervention in cancer.

Inference of the human polyadenylation code

Motivation

Processing of transcripts at the 3′-end involves cleavage at a polyadenylation site followed by the addition of a poly(A)-tail. By selecting which site is cleaved, the process of alternative polyadenylation enables genes to produce transcript isoforms with different 3′-ends. To facilitate the identification and treatment of disease-causing mutations that affect polyadenylation and to understand the sequence determinants underlying this regulatory process, a computational model that can accurately predict polyadenylation patterns from genomic features is desirable.

Results

Previous works have focused on identifying candidate polyadenylation sites and classifying tissue-specific sites. By training on how multiple sites in genes are competitively selected for polyadenylation from 3′-end sequencing data, we developed a deep learning model that can predict the tissue-specific strength of a polyadenylation site in the 3′ untranslated region of the human genome given only its genomic sequence. We demonstrate the model’s broad utility on multiple tasks, without any application-specific training. The model can be used to predict which polyadenylation site is more likely to be selected in genes with multiple sites. It can be used to scan the 3′ untranslated region to find candidate polyadenylation sites. It can be used to classify the pathogenicity of variants near annotated polyadenylation sites in ClinVar. It can also be used to anticipate the effect of antisense oligonucleotide experiments to redirect polyadenylation. We provide analysis on how different features affect the model’s predictive performance and a method to identify sensitive regions of the genome at the single-based resolution that can affect polyadenylation regulation.

Supplementary information

Supplementary data are available at Bioinformatics online.

The role of RNA G-quadruplexes in human diseases and therapeutic strategies

G-quadruplexes (GQs) are four-stranded secondary structures formed by G-rich nucleic acid sequence(s). DNA GQs are present abundantly in the genome and affect a wide range of processes associated with DNA. Recent studies show that RNA GQs are present in different transcripts, including coding and noncoding areas of mRNA, telomeric RNA as well as in other premature and mature noncoding RNAs. When present at specific locations within the RNAs, GQs play important roles in key biological functions, including the regulation of gene expression and telomere homeostasis. RNA GQs regulate pre-mRNA processing, such as splicing and polyadenylation. Evidently, among other processes, RNA GQs also control mRNA translation, miRNA and piRNA biogenesis, and RNA localization. The regulatory mechanisms controlled by RNA GQs mainly involve binding to RNA binding protein that modulate GQ conformation or serve as an entity in recruiting additional protein regulators to act as a block element to the processing machinery. Here we provide an overview of the ever-increasing number of discoveries revealing the role of RNA GQs in biology and their relevance in human diseases and therapeutics. This article is categorized under: RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems RNA in Disease and Development > RNA in Disease.

Genome-wide profiling reveals alternative polyadenylation of mRNA in human non-small cell lung cancer

BACKGROUND

Lung cancer is the second most common cancer with an extremely poor overall survival rate. Post-transcriptional regulation of gene expression play many important roles in human cancer, and one of the potential mechanisms underlying this is alternative mRNA maturation at its 3' untranslated regions (3'-UTRs).

METHODS

Cancer tissues and paired adjacent normal lung tissues from 26 patients diagnosed with non-small cell lung cancer (NSCLC) were analyzed by in vitro transcription-sequencing alternative polyadenylation sites (IVT-SAPAS). 41,773,101 reads in average were obtained from each paired sample. A potential regulation of Cleavage Stimulation Factor Subunit 2 (CSTF2) on 3'UTR length of genes was tested in H460 cells.

RESULTS

1439 (10.26%) genes showed up-regulated expression and 1364 (9.72%) genes showed down-regulated expression in lung cancer tissue versus normal lung tissue, and shorten 3'UTR in cancer tissue was detected in cancer tissues collected from 96.2% (25/26) patients, indicating lung cancer tend to have shortened 3'UTRs of these identified genes. KEGG analysis showed 1855 genes with shorten 3'UTR were enriched in mTOR signaling, ubiquitin mediated proteolysis and RNA degradation. Knocking down CSTF2 expression in H460 cells results in 3'UTR elongation of genes that was identified to be with shortened length in cancer tissues.

CONCLUSION

Alternative polyadenylation (APA) site-switching of 3'UTRs is prevalent in NSCLC, and CSTF2 may serve as an oncogene regulates the 3'UTR length of cancer related genes in NSCLC.

A Deep Neural Network for Predicting and Engineering Alternative Polyadenylation

Alternative polyadenylation (APA) is a major driver of transcriptome diversity in human cells. Here, we use deep learning to predict APA from DNA sequence alone. We trained our model (APARENT, APA REgression NeT) on isoform expression data from over 3 million APA reporters. APARENT's predictions are highly accurate when tasked with inferring APA in synthetic and human 3'UTRs. Visualizing features learned across all network layers reveals that APARENT recognizes sequence motifs known to recruit APA regulators, discovers previously unknown sequence determinants of 3' end processing, and integrates these features into a comprehensive, interpretable, cis-regulatory code. We apply APARENT to forward engineer functional polyadenylation signals with precisely defined cleavage position and isoform usage and validate predictions experimentally. Finally, we use APARENT to quantify the impact of genetic variants on APA. Our approach detects pathogenic variants in a wide range of disease contexts, expanding our understanding of the genetic origins of disease.