Alternative Polyadenylation in Human Diseases

PABPN1 loss-of-function causes APA-shift in oculopharyngeal muscular dystrophy

Alternative polyadenylation (APA) at the 3' UTR of transcripts contributes to the cell transcriptome. APA is suppressed by the nuclear RNA-binding protein PABPN1. Aging-associated reduced PABPN1 levels in skeletal muscles lead to muscle wasting. Muscle weakness in oculopharyngeal muscular dystrophy (OPMD) is caused by short alanine expansion in PABPN1 exon1. The expanded PABPN1 forms nuclear aggregates, an OPMD hallmark. Whether the expanded PABPN1 affects APA and how it contributes to muscle pathology is unresolved. To investigate these questions, we developed a procedure including RNA library preparation and a simple pipeline calculating the APA-shift ratio as a readout for PABPN1 activity. Comparing APA-shift results to previously published PAS utilization and APA-shift results, we validated this procedure. The procedure was then applied on the OPMD cell model and on RNA from OPMD muscles. APA-shift was genome-wide in the mouse OPMD model, primarily affecting muscle transcripts. In OPMD individuals, APA-shift was enriched with muscle transcripts. In an OPMD cell model APA-shift was not significant. APA-shift correlated with reduced expression levels of a subset of PABPN1 isoforms, whereas the expression of the expanded PABPN1 did not correlate with APA-shift. PABPN1 activity is not affected by the expression of expanded PABPN1, but rather by reduced PABPN1 expression levels. In muscles, PABPN1 activity initially affects muscle transcripts. We suggest that muscle weakness in OPMD is caused by PABPN1 loss-of-function leading to APA-shift that primarily affects in muscle transcripts.

An AluYa5 Insertion in the 3'UTR of COL4A1 and Cerebral Small Vessel Disease

Importance

Cerebral small vessel diseases (CSVDs) account for one-fifth of stroke cases. Numerous familial cases remain unresolved after routine screening of known CSVD genes.

Objective

To identify novel genes and mechanisms associated with familial CSVD.

Design, Setting, and Participants

This 2-stage study involved linkage analysis and a case-control study; linkage analysis and whole exome and genome sequencing were used to identify candidate gene variants in 2 large families with CSVD (9 patients with CSVD). Then, a case-control analysis was conducted on 246 unrelated probands, including probands from these 2 families and 244 additional probands. All probands (clinical onset

Main Outcomes and Measures

A pathogenic AluYa5 insertion was identified within the COL4A1 3'UTR in the 2 large families with CSVD. Reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR), Western blot, and long-read RNA sequencing were used to investigate outcomes associated with the insertion using patient fibroblasts. Clinical and magnetic resonance imaging features of probands with variants and available relatives were assessed.

Results

Among 246 probands (141 females [57.3%]; median [IQR] age at referral, 56 [49-64] years), 7 patients of French ancestry carried the insertion. This insertion was absent in 467 healthy French individuals in a control group (odds ratio, ∞; 95% CI, 2.78 to ∞; P = 5 × 10-4) and 10 847 individuals from the gnomAD structural variant database (odds ratio, ∞; 95% CI, 64.77 to ∞; P = 2.42 × 10-12). In these 7 patients' families, 19 family members with CSVD carried the insertion. RT-qPCR and Western blot showed an upregulation of COL4A1 mRNA (10.6-fold increase; 95% CI, 1.4-fold to 17.1-fold increase) and protein levels (2.8-fold increase; 95% CI, 2.1-fold to 3.5-fold increase) in patient vs control group fibroblasts. Long-read RNA sequencing data showed that the insertion was associated with perturbation in the use of canonical COL4A1 polyadenylation signals (approximately 87% of isoforms transcribed from the wild type allele vs 5% of isoforms transcribed from the allele with the insertion used the 2 distal canonical polyadenylation signals). The main clinical feature of individuals with CSVD was the recurrence of pontine ischemic lesions starting at an early age (17 of 19 patients [89.5%]).

Conclusions and relevance

This study found a novel mechanism associated with COL4A1 upregulation and a highly penetrant adult-onset CSVD. These findings suggest that quantitative alterations of the cerebrovascular matrisome are associated with CSVD pathogenesis, with diagnostic and therapeutic implications.

Collapse

Key Words Collapse

MESH Headings

• Adult
• Female
• Humans
• Middle Aged
• 3' Untranslated Regions/genetics
• Alleles
• Case-Control Studies
• Cerebral Small Vessel Diseases/genetics
• Collagen Type IV/metabolism
• Protein Isoforms
• Mutagenesis, Insertional

Collapse

Grants Collapse

Affiliation(s)

Chaker Aloui NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France
Lisa Neumann NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France
Françoise Bergametti NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France
Eric Sartori Service de Neurologie, Centre Hospitalier Bretagne Sud, Lorient, France
Marc Herbreteau Service de Neurologie, Centre Hospitalier Bretagne Sud, Lorient, France
Arnaud Maillard Assistance Publique-Hôpitaux de Paris, Service de Génétique Moléculaire Neurovasculaire, Hôpital Saint-Louis, Paris, France
Thibault Coste NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France Assistance Publique-Hôpitaux de Paris, Service de Génétique Moléculaire Neurovasculaire, Hôpital Saint-Louis, Paris, France
Hélène Morel NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France Assistance Publique-Hôpitaux de Paris, Service de Génétique Moléculaire Neurovasculaire, Hôpital Saint-Louis, Paris, France
Dominique Hervé NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France Assistance Publique-Hôpitaux de Paris, Service de Neurologie, Hôpital Lariboisière, Paris, France
Hugues Chabriat NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France Assistance Publique-Hôpitaux de Paris, Service de Neurologie, Hôpital Lariboisière, Paris, France
Serge Timsit Service de Neurologie Vasculaire, Centre Hospitalier Régional Universitaire de Brest, Brest, France
Irina Viakhireva Service de Neurologie Vasculaire, Centre Hospitalier Régional Universitaire de Brest, Brest, France
Yves Denoyer Service de Neurologie, Centre Hospitalier Bretagne Sud, Lorient, France Université de Rennes, Laboratoire Traitement du Signal et de l'Image, Institut National de la Santé Et de la Recherche Médicale Unité Mixte de Recherche 1099, Rennes, France
Rémi Allibert Service de Neurologie, Unité Neurovasculaire, Centre Hospitalier Universitaire de Saint Etienne, Saint Etienne, France
Florence Demurger Service de Neurologie, Unité Neurovasculaire, Centre Hospitalier Bretagne Atlantique, Vannes, France
Cedric Gollion Service de Neurologie, Centre Hospitalier Universitaire de Toulouse, Toulouse, France
Patrick Vermersch Univ. Lille, Institut National de la Santé Et de la Recherche Médicale Unité Mixte de Recherche 1172 LilNCog, Centre Hospitalier Universitaire Lille, Fédérations Hospitalo-Universitaire Precise, Lille, France
Florence Marchelli Assistance Publique-Hôpitaux de Paris, Service de Génétique Moléculaire Neurovasculaire, Hôpital Saint-Louis, Paris, France
Corinne Blugeon GenomiqueENS, Institut de Biologie de l’Ecole Normale Supérieur, Département de biologie, École Normale Supérieure, Centre National de la Recherche Scientifique, Institut National de la Santé Et de la Recherche Médicale, Université Paris Sciences et Lettres, Paris, France
Sophie Lemoine GenomiqueENS, Institut de Biologie de l’Ecole Normale Supérieur, Département de biologie, École Normale Supérieure, Centre National de la Recherche Scientifique, Institut National de la Santé Et de la Recherche Médicale, Université Paris Sciences et Lettres, Paris, France
Claire Tourtier-Bellosta Service de Neurologie, Centre Hospitalier Bretagne Sud, Lorient, France
Alexis Brouazin Service de neurologie, Centre Hospitalier de Cornouaille, Quimper, France
Anne-Louise Leutenegger NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France
Eva Pipiras NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France Assistance Publique-Hôpitaux de Paris, Hôpitaux Jean Verdier et Armand Trousseau, Université Sorbonne Paris Nord, Bobigny, France
Elisabeth Tournier-Lasserve NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France Assistance Publique-Hôpitaux de Paris, Service de Génétique Moléculaire Neurovasculaire, Hôpital Saint-Louis, Paris, France

Collapse

Alternative polyadenylation reprogramming of MORC2 induced by NUDT21 loss promotes KIRC carcinogenesis

Alternative polyadenylation (APA), a posttranscriptional mechanism of gene expression via determination of 3'UTR length, has an emerging role in carcinogenesis. Although abundant APA reprogramming is found in kidney renal clear cell carcinoma (KIRC), which is one of the major malignancies, whether APA functions in KIRC remains unknown. Herein, we found that chromatin modifier MORC2 gained oncogenic potential in KIRC among the genes with APA reprogramming, and moreover, its oncogenic potential was enhanced by 3'UTR shortening through stabilization of MORC2 mRNA. MORC2 was found to function in KIRC by downregulating tumor suppressor DAPK1 via DNA methylation. Mechanistically, MORC2 recruited DNMT3A to facilitate hypermethylation of the DAPK1 promoter, which was strengthened by 3'UTR shortening of MORC2. Furthermore, loss of APA regulator NUDT21, which was induced by DNMT3B-mediated promoter methylation, was identified as responsible for 3'UTR shortening of MORC2 in KIRC. Additionally, NUDT21 was confirmed to act as a tumor suppressor mainly depending on downregulation of MORC2. Finally, we designed an antisense oligonucleotide (ASO) to enhance NUDT21 expression and validated its antitumor effect in vivo and in vitro. This study uncovers the DNMT3B/NUDT21/APA/MORC2/DAPK1 regulatory axis in KIRC, disclosing the role of APA in KIRC and the crosstalk between DNA methylation and APA.

Alternative polyadenylation regulation in cardiac development and cardiovascular disease

Cleavage and polyadenylation of pre-mRNAs is a necessary step for gene expression and function. Majority of human genes exhibit multiple polyadenylation sites, which can be alternatively used to generate different mRNA isoforms from a single gene. Alternative polyadenylation (APA) of pre-mRNAs is important for the proteome and transcriptome landscape. APA is tightly regulated during development and contributes to tissue-specific gene regulation. Mis-regulation of APA is linked to a wide range of pathological conditions. APA-mediated gene regulation in the heart is emerging as a new area of research. Here, we will discuss the impact of APA on gene regulation during heart development and in cardiovascular diseases. First, we will briefly review how APA impacts gene regulation and discuss molecular mechanisms that control APA. Then, we will address APA regulation during heart development and its dysregulation in cardiovascular diseases. Finally, we will discuss pre-mRNA targeting strategies to correct aberrant APA patterns of essential genes for the treatment or prevention of cardiovascular diseases. The RNA field is blooming due to advancements in RNA-based technologies. RNA-based vaccines and therapies are becoming the new line of effective and safe approaches for the treatment and prevention of human diseases. Overall, this review will be influential for understanding gene regulation at the RNA level via APA in the heart and will help design RNA-based tools for the treatment of cardiovascular diseases in the future.

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.

Comprehensive analysis of alternative polyadenylation regulators concerning CD276 and immune infiltration in bladder cancer

Alternative polyadenylation (APA) is emerging as a crucial regulatory mechanism in bladder cancer (BC), while it remains elusive whether APA influences the tumor immune microenvironment (TIME) in BC. We identified two distinct subtypes of BC by APA-related regulatory genes expression profiles. The two subtypes have different pathological grades, prognostic outcomes, tumor immune infiltration characteristics, and pathway enrichment. Subsequently, CPSF3 was identified as a potential immune infiltration-related gene in BC. Highly expressed CPSF3 was positively correlated with unfavorable prognosis and high CD276 expression in BC. Moreover, we verified the expression of CPSF3 in BC tissues and cell lines by qRT-PCR. In conclusion, the study indicates that APA regulatory factors play an important role in immune infiltration of BC, and that CPSF3 was a potentially prognostic marker and immunotherapy target for BC.

Alternative polyadenylation associated with prognosis and therapy in colorectal cancer

Colorectal cancer (CRC) is among the most widely spread cancers globally. Aberrant alternative polyadenylation (APA) plays a role in cancer onset and its progression. Consequently, this study focused on highlighting the role of APA events and signals in the prognosis of patients with CRC. The APA events, RNA sequencing (RNA-seq), somatic mutations, copy number variants (CNVs), and clinical information of the CRC cohort were obtained from The Cancer Genome Atlas (TCGA) database and UCSC (University of California-Santa Cruz) Xena database. The whole set was sorted into two sets: a training set and a test set in a ratio of 7:3. 197 prognosis-related APA events were collected by performing univariate Cox regression signature in patients with CRC. Subsequently, a signature for APA events was established by least absolute shrinkage and selection operator (LASSO) and multivariate Cox analysis. The risk scores were measured for individual patients on the basis of the signature and patients were sorted into two groups; the high-risk group and the low-risk group as per their median risk scores. Kaplan–Meier curves, principal component analysis (PCA), and time-dependent receiver operator characteristic (ROC) curves revealed that the signature was able to predict patient prognosis effectively and further validation was provided in the test set and the whole set. The high-risk and low-risk groups displayed various distributions of mutations and CNVs. Tumor mutation burden (TMB) alone and in combination with the signature predicted the prognosis of CRC patients, but the gene frequencies of TMBs and CNVs did not change in the low- and high-risk groups. Moreover, immunotherapy and chemotherapy treatments showed different responses to PD-1 inhibitors and multiple chemotherapeutic agents in the low and high-risk groups based on the tumor immune dysfunction and exclusion (TIDE) and genomics of drugs sensitivity in cancer (GDSC) databases. This study may help in understanding the potential roles of APA in CRC, and the signature for prognosis-related APA events can work as a potential predictor for survival and treatment in patients with CRC.

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.

Shortening of HO1 3'UTRs by Alternative Polyadenylation Suppresses Adipogenesis in 3T3-L1

Appropriately increasing intramuscular fat content can help improve meat quality, so it is necessary to explore the internal molecular mechanism of preadipocyte differentiation. The role of heme oxygenase 1 (HO1) in cell oxidative stress, energy metabolism, cell proliferation, and differentiation has gradually been revealed. Here, we used 3'RACE to identify the full-length 3' untranslated region (3'UTR) of HO1 and found that a very short 3'UTR variant was produced by alternative polyadenylation (APA). HO1 with a long 3'UTR variant was identified as a direct target of miR155-5P and miR377-3P. Our experimental results verified the inhibitory effect of HO1 on preadipocyte differentiation. In addition, our research confirms that by escaping microRNA inhibitory effects, the HO1 3'UTR short variant produced by APA has a higher level of expression. Thus, the HO1 3'UTR short variant has a stronger inhibitory effect on the preadipocyte differentiation than the HO1 3'UTR long variants in 3T3-L1.

A computational pipeline to infer alternative poly-adenylation from 3' sequencing data

An increasing number of investigations have established alternative polyadenylation (APA) as a key mechanism of gene regulation through altering the length of 3' untranslated region (UTR) and generating distinct mRNA termini. Further, appreciation for the significance of APA in disease contexts propelled the development of several 3' sequencing techniques. While these RNA sequencing technologies have advanced APA analysis, the intrinsic limitation of 3' read coverage and lack of appropriate computational tools constrain precise mapping and quantification of polyadenylation sites. Notably, Poly(A)-ClickSeq (PAC-seq) overcomes limiting factors such as poly(A) enrichment and 3' linker ligation steps using click-chemistry. Here we provide an updated PolyA-miner protocol, a computational approach to analyze PAC-seq or other 3'-Seq datasets. As a key practical constraint, we also provide a detailed account on the impact of sequencing depth on the number of detected polyadenylation sites and APA changes. This protocol is also updated to handle unique molecular identifiers used to address PCR duplication potentially observed in PAC-seq.

Generation of 3'UTR knockout cell lines by CRISPR/Cas9-mediated genome editing

In addition to the protein code, messenger RNAs (mRNAs) also contain untranslated regions (UTRs). 3'UTRs span the region between the translational stop codon and the poly(A) tail. Sequence elements located in 3'UTRs are essential for pre-mRNA processing. 3'UTRs also contain elements that can regulate protein abundance, localization, and function. At least half of all human genes use alternative cleavage and polyadenylation (APA) to further diversify the regulatory potential of protein functions. Traditional gene editing approaches are designed to disrupt the production of functional proteins. Here, we describe a method that allows investigators to manipulate 3'UTR sequences of endogenous genes for both single- 3'UTR and multi-3'UTR genes. As 3'UTRs can regulate individual functions of proteins, techniques to manipulate 3'UTRs at endogenous gene loci will help to disentangle multi-functionality of proteins. Furthermore, the ability to directly examine the impact of gene regulatory elements in 3'UTRs will provide further insights into their functional significance.

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.

Alternative polyadenylation: methods, mechanism, function, and role in cancer

Occurring in over 60% of human genes, alternative polyadenylation (APA) results in numerous transcripts with differing 3'ends, thus greatly expanding the diversity of mRNAs and of proteins derived from a single gene. As a key molecular mechanism, APA is involved in various gene regulation steps including mRNA maturation, mRNA stability, cellular RNA decay, and protein diversification. APA is frequently dysregulated in cancers leading to changes in oncogenes and tumor suppressor gene expressions. Recent studies have revealed various APA regulatory mechanisms that promote the development and progression of a number of human diseases, including cancer. Here, we provide an overview of four types of APA and their impacts on gene regulation. We focus particularly on the interaction of APA with microRNAs, RNA binding proteins and other related factors, the core pre-mRNA 3'end processing complex, and 3'UTR length change. We also describe next-generation sequencing methods and computational tools for use in poly(A) signal detection and APA repositories and databases. Finally, we summarize the current understanding of APA in cancer and provide our vision for future APA related research.

Animal-APAdb: a comprehensive animal alternative polyadenylation database

Alternative polyadenylation (APA) is an important post-transcriptional regulatory mechanism that recognizes different polyadenylation signals on transcripts, resulting in transcripts with different lengths of 3′ untranslated regions and thereby influencing a series of biological processes. Recent studies have highlighted the important roles of APA in human. However, APA profiles in other animals have not been fully recognized, and there is no database that provides comprehensive APA information for other animals except human. Here, by using the RNA sequencing data collected from public databases, we systematically characterized the APA profiles in 9244 samples of 18 species. In total, we identified 342 952 APA events with a median of 17 020 per species using the DaPars2 algorithm, and 315 691 APA events with a median of 17 953 per species using the QAPA algorithm in these 18 species, respectively. In addition, we predicted the polyadenylation sites (PAS) and motifs near PAS of these species. We further developed Animal-APAdb, a user-friendly database (http://gong_lab.hzau.edu.cn/Animal-APAdb/) for data searching, browsing and downloading. With comprehensive information of APA events in different tissues of different species, Animal-APAdb may greatly facilitate the exploration of animal APA patterns and novel mechanisms, gene expression regulation and APA evolution across tissues and species.

Emerging cellular and molecular determinants of idiopathic pulmonary fibrosis

Idiopathic pulmonary fibrosis (IPF), the most common form of idiopathic interstitial pneumonia, is a progressive, irreversible, and typically lethal disease characterized by an abnormal fibrotic response involving vast areas of the lungs. Given the poor knowledge of the mechanisms underpinning IPF onset and progression, a better understanding of the cellular processes and molecular pathways involved is essential for the development of effective therapies, currently lacking. Besides a number of established IPF-associated risk factors, such as cigarette smoking, environmental factors, comorbidities, and viral infections, several other processes have been linked with this devastating disease. Apoptosis, senescence, epithelial-mesenchymal transition, endothelial-mesenchymal transition, and epithelial cell migration have been shown to play a key role in IPF-associated tissue remodeling. Moreover, molecules, such as chemokines, cytokines, growth factors, adenosine, glycosaminoglycans, non-coding RNAs, and cellular processes including oxidative stress, mitochondrial dysfunction, endoplasmic reticulum stress, hypoxia, and alternative polyadenylation have been linked with IPF development. Importantly, strategies targeting these processes have been investigated to modulate abnormal cellular phenotypes and maintain tissue homeostasis in the lung. This review provides an update regarding the emerging cellular and molecular mechanisms involved in the onset and progression of IPF.

PolyA-miner: accurate assessment of differential alternative poly-adenylation from 3'Seq data using vector projections and non-negative matrix factorization

Almost 70% of human genes undergo alternative polyadenylation (APA) and generate mRNA transcripts with varying lengths, typically of the 3′ untranslated regions (UTR). APA plays an important role in development and cellular differentiation, and its dysregulation can cause neuropsychiatric diseases and increase cancer severity. Increasing awareness of APA’s role in human health and disease has propelled the development of several 3′ sequencing (3′Seq) techniques that allow for precise identification of APA sites. However, despite the recent data explosion, there are no robust computational tools that are precisely designed to analyze 3′Seq data. Analytical approaches that have been used to analyze these data predominantly use proximal to distal usage. With about 50% of human genes having more than two APA isoforms, current methods fail to capture the entirety of APA changes and do not account for non-proximal to non-distal changes. Addressing these key challenges, this study demonstrates PolyA-miner, an algorithm to accurately detect and assess differential alternative polyadenylation specifically from 3′Seq data. Genes are abstracted as APA matrices, and differential APA usage is inferred using iterative consensus non-negative matrix factorization (NMF) based clustering. PolyA-miner accounts for all non-proximal to non-distal APA switches using vector projections and reflects precise gene-level 3′UTR changes. It can also effectively identify novel APA sites that are otherwise undetected when using reference-based approaches. Evaluation on multiple datasets—first-generation MicroArray Quality Control (MAQC) brain and Universal Human Reference (UHR) PolyA-seq data, recent glioblastoma cell line NUDT21 knockdown Poly(A)-ClickSeq (PAC-seq) data, and our own mouse hippocampal and human stem cell-derived neuron PAC-seq data—strongly supports the value and protocol-independent applicability of PolyA-miner. Strikingly, in the glioblastoma cell line data, PolyA-miner identified more than twice the number of genes with APA changes than initially reported. With the emerging importance of APA in human development and disease, PolyA-miner can significantly improve data analysis and help decode the underlying APA dynamics.

SNP2APA: a database for evaluating effects of genetic variants on alternative polyadenylation in human cancers

Alternative polyadenylation (APA) is an important post-transcriptional regulation that recognizes different polyadenylation signals (PASs), resulting in transcripts with different 3' untranslated regions, thereby influencing a series of biological processes and functions. Recent studies have revealed that some single nucleotide polymorphisms (SNPs) could contribute to tumorigenesis and development through dysregulating APA. However, the associations between SNPs and APA in human cancers remain largely unknown. Here, using genotype and APA data of 9082 samples from The Cancer Genome Atlas (TCGA) and The Cancer 3'UTR Altas (TC3A), we systematically identified SNPs affecting APA events across 32 cancer types and defined them as APA quantitative trait loci (apaQTLs). As a result, a total of 467 942 cis-apaQTLs and 30 721 trans-apaQTLs were identified. By integrating apaQTLs with survival and genome-wide association studies (GWAS) data, we further identified 2154 apaQTLs associated with patient survival time and 151 342 apaQTLs located in GWAS loci. In addition, we designed an online tool to predict the effects of SNPs on PASs by utilizing PAS motif prediction tool. Finally, we developed SNP2APA, a user-friendly and intuitive database (http://gong_lab.hzau.edu.cn/SNP2APA/) for data browsing, searching, and downloading. SNP2APA will significantly improve our understanding of genetic variants and APA in human cancers.

Unraveling the RNA Binding Properties of the Iron-Sulfur Zinc Finger Protein CPSF30

Cleavage and polyadenylation specificity factor 30 (CPSF30) is a "zinc finger" protein that plays a crucial role in the transition of pre-mRNA to RNA. CPSF30 contains five conserved CCCH domains and a CCHC "zinc knuckle" domain. CPSF30 activity is critical for pre-mRNA processing. A truncated form of the protein, in which only the CCCH domains are present, has been shown to specifically bind AU-rich pre-mRNA targets; however, the RNA binding and recognition properties of full-length CPSF30 are not known. Herein, we report the isolation and biochemical characterization of full-length CPSF30. We report that CPSF30 contains one 2Fe-2S cluster in addition to five zinc ions, as measured by inductively coupled plasma mass spectrometry, ultraviolet-visible spectroscopy, and X-ray absorption spectroscopy. Utilizing fluorescence anisotropy RNA binding assays, we show that full-length CPSF30 has high binding affinity for two types of pre-mRNA targets, AAUAAA and polyU, both of which are conserved sequence motifs present in the majority of pre-mRNAs. Binding to the AAUAAA motif requires that the five CCCH domains of CPSF30 be present, whereas binding to polyU sequences requires the entire, full-length CPSF30. These findings implicate the CCHC "zinc knuckle" present in the full-length protein as being critical for mediating polyU binding. We also report that truncated forms of the protein, containing either just two CCCH domains (ZF2 and ZF3) or the CCHC "zinc knuckle" domain, do not exhibit any RNA binding, indicating that CPSF30/RNA binding requires several ZF (and/or Fe-S cluster) domains working in concert to mediate RNA recognition.

The DEAD-box RNA helicase SHI2 functions in repression of salt-inducible genes and regulation of cold-inducible gene splicing

Gene regulation is central for growth, development, and adaptation to environmental changes in all living organisms. Many genes are induced by environmental cues, and the expression of these inducible genes is often repressed under normal conditions. Here, we show that the SHINY2 (SHI2) gene is important for repressing salt-inducible genes and also plays a role in cold response. The shi2 mutant displayed hypersensitivity to cold, abscisic acid (ABA), and LiCl. Map-based cloning demonstrates that SHI2 encodes a DEAD- (Asp-Glu-Ala-Asp) box RNA helicase with similarity to a yeast splicing factor. Transcriptomic analysis of the shi2 mutant in response to cold revealed that the shi2 mutation decreased the number of cold-responsive genes and the magnitude of their response, and resulted in the mis-splicing of some cold-responsive genes. Under salt stress, however, the shi2 mutation increased the number of salt-responsive genes but had a negligible effect on mRNA splicing. Our results suggest that SHI2 is a component in a ready-for-transcription repressor complex important for gene repression under normal conditions, and for gene activation and transcription under stress conditions. In addition, SHI2 also serves as a splicing factor required for proper splicing of cold-responsive genes and affects 5' capping and polyadenylation site selection.

Stiff matrix instigates type I collagen biogenesis by mammalian cleavage factor I complex-mediated alternative polyadenylation

Alternative polyadenylation (APA) is a widespread and important mechanism in regulation of gene expression. Dysregulation of the 3' UTR cleavage and polyadenylation represents a common characteristic among many disease states, including lung fibrosis. In this study, we investigated the role of mammalian cleavage factor I-mediated (CFIm-mediated) APA in regulating extracellular matrix production in response to mechanical stimuli from stiffened matrix simulating the fibrotic lungs. We found that stiff matrix downregulated expression of CFIm68, CFIm59 and CFIm25 subunits and promoted APA in favor of the proximal poly(A) site usage in the 3' UTRs of type I collagen (COL1A1) and fibronectin (FN1) in primary human lung fibroblasts. Knockdown and overexpression of each individual CFIm subunit demonstrated that CFIm68 and CFIm25 are indispensable attributes of stiff matrix-induced APA and overproduction of COL1A1, whereas CFIm did not appear to mediate stiffness-regulated FN1 APA. Furthermore, expression of the CFIm subunits was associated with matrix stiffness in vivo in a bleomycin-induced mouse model of pulmonary fibrosis. These data suggest that stiff matrix instigates type I collagen biogenesis by selectively targeting mRNA transcripts for 3' UTR shortening. The current study uncovered a potential mechanism for regulation of the CFIm complex by mechanical cues under fibrotic conditions.

Template-switching artifacts resemble alternative polyadenylation

BACKGROUND

Alternative polyadenylation is commonly examined using cDNA sequencing, which is known to be affected by template-switching artifacts. However, the effects of such template-switching artifacts on alternative polyadenylation are generally disregarded, while alternative polyadenylation artifacts are attributed to internal priming.

RESULTS

Here, we analyzed both long-read cDNA sequencing and direct RNA sequencing data of two organisms, generated by different sequencing platforms. We developed a filtering algorithm which takes into consideration that template-switching can be a source of artifactual polyadenylation when filtering out spurious polyadenylation sites. The algorithm outperformed the conventional internal priming filters based on comparison to direct RNA sequencing data. We also showed that the polyadenylation artifacts arise in cDNA sequencing at consecutive stretches of as few as three adenines. There was no substantial difference between the lengths of poly(A) tails at the artifactual and the true transcriptional end sites even though it is expected that internal priming artifacts have shorter poly(A) tails than genuine polyadenylated reads.

CONCLUSIONS

Our findings suggest that template switching plays an important role in the generation of spurious polyadenylation and support the need for more rigorous filtering of artifactual polyadenylation sites in cDNA data, or that alternative polyadenylation should be annotated using native RNA sequencing.

Alternative polyadenylation of mRNA and its role in cancer

Alternative polyadenylation (APA) is a molecular process that generates diversity at the 3′ end of RNA polymerase II transcripts from over 60% of human genes. APA is derived from the existence of multiple polyadenylation signals (PAS) within the same transcript, and results in the differential inclusion of sequence information at the 3′ end. While APA can occur between two PASs allowing for generation of transcripts with distinct coding potential from a single gene, most APA occurs within the untranslated region (3′UTR) and changes the length and content of these non-coding sequences. APA within the 3′UTR can have tremendous impact on its regulatory potential of the mRNA through a variety of mechanisms, and indeed this layer of gene expression regulation has profound impact on processes vital to cell growth and development. Recent studies have particularly highlighted the importance of APA dysregulation in cancer onset and progression. Here, we review the current knowledge of APA and its impacts on mRNA stability, translation, localization and protein localization. We also discuss the implications of APA dysregulation in cancer research and therapy.

The Length of the Expressed 3' UTR Is an Intermediate Molecular Phenotype Linking Genetic Variants to Complex Diseases

In the last decades, genome-wide association studies (GWAS) have uncovered tens of thousands of associations between common genetic variants and complex diseases. However, these statistical associations can rarely be interpreted functionally and mechanistically. As the majority of the disease-associated variants are located far from coding sequences, even the relevant gene is often unclear. A way to gain insight into the relevant mechanisms is to study the genetic determinants of intermediate molecular phenotypes, such as gene expression and transcript structure. We propose a computational strategy to discover genetic variants affecting the relative expression of alternative 3′ untranslated region (UTR) isoforms, generated through alternative polyadenylation, a widespread posttranscriptional regulatory mechanism known to have relevant functional consequences. When applied to a large dataset in which whole genome and RNA sequencing data are available for 373 European individuals, 2,530 genes with alternative polyadenylation quantitative trait loci (apaQTL) were identified. We analyze and discuss possible mechanisms of action of these variants, and we show that they are significantly enriched in GWAS hits, in particular those concerning immune-related and neurological disorders. Our results point to an important role for genetically determined alternative polyadenylation in affecting predisposition to complex diseases, and suggest new ways to extract functional information from GWAS data.

Alternative cleavage and polyadenylation in health and disease

Most human genes have multiple sites at which RNA 3' end cleavage and polyadenylation can occur, enabling the expression of distinct transcript isoforms under different conditions. Novel methods to sequence RNA 3' ends have generated comprehensive catalogues of polyadenylation (poly(A)) sites; their analysis using innovative computational methods has revealed how poly(A) site choice is regulated by core RNA 3' end processing factors, such as cleavage factor I and cleavage and polyadenylation specificity factor, as well as by other RNA-binding proteins, particularly splicing factors. Here, we review the experimental and computational methods that have enabled the global mapping of mRNA and of long non-coding RNA 3' ends, quantification of the resulting isoforms and the discovery of regulators of alternative cleavage and polyadenylation (APA). We highlight the different types of APA-derived isoforms and their functional differences, and illustrate how APA contributes to human diseases, including cancer and haematological, immunological and neurological diseases.

Tissue-specific mechanisms of alternative polyadenylation: Testis, brain, and beyond (2018 update)

Alternative polyadenylation (APA) is how genes choose different sites for 3′ end formation for mRNAs during transcription. APA often occurs in a tissue‐ or developmental stage‐specific manner that can significantly affect gene activity by changing the protein product generated, the stability of the transcript, its localization within the cell, or its translatability. Despite the important regulatory effects that APA has on tissue‐specific gene expression, only a few examples have been characterized mechanistically. In this 2018 update to our 2010 review, we examine mechanisms for the control of APA and update our understanding of the older mechanisms since 2010. We once postulated the existence of tissue‐specific factors in APA. However, while a few tissue‐specific polyadenylation factors are known, the emerging conclusion is that the majority of APA is accomplished by altering levels of core polyadenylation proteins. Examples of those core proteins include CSTF2, CPSF1, and subunits of mammalian cleavage factor I. But despite support for these mechanisms, no one has yet documented any of these proteins changing in either a tissue‐specific or developmental manner. Given the profound effect that APA can have on gene expression and human health, improved understanding of tissue‐specific APA could lead to numerous advances in gene activity control.

This article is categorized under:

RNA Processing > 3′ End Processing

RNA in Disease and Development > RNA in Development