Differential genome-wide profiling of tandem 3' UTRs among human breast cancer and normal cells by high-throughput sequencing

Star-PAP controlled alternative polyadenylation coupled poly(A) tail length regulates protein expression in hypertrophic heart

Alternative polyadenylation (APA)-mediated 3′-untranslated region (UTR) shortening is known to increase protein expression due to the loss of miRNA regulatory sites. Yet, mRNAs with longer 3′-UTR also show enhanced protein expression. Here, we identify a mechanism by which longer transcripts generated by the distal-most APA site leads to increased protein expression compared to the shorter transcripts and the longer transcripts are positioned to regulate heart failure (HF). A Star-PAP target gene, NQO1 has three poly(A) sites (PA-sites) at the terminal exon on the pre-mRNA. Star-PAP selects the distal-most site that results in the expression of the longest isoform. We show that the NQO1 distal-specific mRNA isoform accounts for the majority of cellular NQO1 protein. Star-PAP control of the distal-specific isoform is stimulated by oxidative stress and the toxin dioxin. The longest NQO1 transcript has increased poly(A) tail (PA-tail) length that accounts for the difference in translation potentials of the three NQO1 isoforms. This mechanism is involved in the regulation of cardiac hypertrophy (CH), an antecedent condition to HF where NQO1 downregulation stems from the loss of the distal-specific transcript. The loss of NQO1 during hypertrophy was rescued by ectopic expression of the distal- but not the proximal- or middle-specific NQO1 mRNA isoforms in the presence of Star-PAP expression, and reverses molecular events of hypertrophy in cardiomyocytes.

A Cancer-Specific Ubiquitin Ligase Drives mRNA Alternative Polyadenylation by Ubiquitinating the mRNA 3' End Processing Complex

Alternative polyadenylation (APA) contributes to transcriptome complexity by generating mRNA isoforms with varying 3' UTR lengths. APA leading to 3' UTR shortening (3' US) is a common feature of most cancer cells; however, the molecular mechanisms are not understood. Here, we describe a widespread mechanism promoting 3' US in cancer through ubiquitination of the mRNA 3' end processing complex protein, PCF11, by the cancer-specific MAGE-A11-HUWE1 ubiquitin ligase. MAGE-A11 is normally expressed only in the male germline but is frequently re-activated in cancers. MAGE-A11 is necessary for cancer cell viability and is sufficient to drive tumorigenesis. Screening for targets of MAGE-A11 revealed that it ubiquitinates PCF11, resulting in loss of CFIm25 from the mRNA 3' end processing complex. This leads to APA of many transcripts affecting core oncogenic and tumor suppressors, including cyclin D2 and PTEN. These findings provide insights into the molecular mechanisms driving APA in cancer and suggest therapeutic strategies.

LanceletDB: an integrated genome database for lancelet, comparing domain types and combination in orthologues among lancelet and other species

Lancelet (amphioxus) represents the most basally divergent extant chordate (cephalochordates) that diverged from the other two chordate lineages (urochordates and vertebrates) more than half a billion years ago. As it occupies a key position in evolution, it is considered as one of the best proxies for understanding the chordate ancestral state. Thus, the construction of a database with multiple lancelet genomes and gene annotation data, including protein domains, is urgently needed to investigate the loss and gain of domains in orthologues among species, especially ancient domain types (non-vertebrate-specific domains) and novel domain combination, which is helpful for providing new insight into the chordate ancestral state and vertebrate evolution. Here, we present an integrated genome database for lancelet, LanceletDB, which provides reference haploid genome sequence and annotation data for lancelet (Branchiostoma belcheri), including gene models and annotation, protein domain types, gene expression pattern in embryogenesis, different expression sequence tag sets and alternative polyadenylation (APA) sites profiled by the sequencing APA sites method. Especially, LanceletDB allows comparison of domain types and combination in orthologues among type species so as to decode the ancient domain types and novel domain combination during evolution. We also integrated the released diploid lancelet genome annotation data (Branchiostoma floridae) to expand LanceletDB and extend its usefulness. These data are available through the search and analysis page, basic local alignment search tool page and genome browser to provide an integrated display.

Global Identification of Human Exosome Substrates Using RNA Interference and RNA Sequencing

The RNA exosome is involved in RNA processing and quality control. In humans, it consists of an enzymatically inactive nine-subunit core, with ribonucleolytic activity contributed by one or two additional components. Moreover, several protein cofactors interact with the exosome to enable and specify its recruitment to a wide range of substrates. A common strategy to identify these substrates has been to deplete an exosome subunit or a cofactor and subsequently interrogate which transcripts become stabilized. Here, we describe an experimental pipeline including siRNA-mediated depletion of the RNA exosome or its cofactors in HeLa cells, confirmation of the knockdown efficiencies, and the manual or high-throughput identification of exosome targets.

PlantAPAdb: A Comprehensive Database for Alternative Polyadenylation Sites in Plants

Alternative cleavage and polyadenylation (APA) is increasingly recognized as an important regulatory mechanism in eukaryotic gene expression and is dynamically modulated in a developmental, tissue-specific, or environmentally responsive manner. Given the functional importance of APA and the rapid accumulation of APA sites in plants, a comprehensive and easily accessible APA site database is necessary for improved understanding of APA-mediated gene expression regulation. We present a database called PlantAPAdb that catalogs the most comprehensive APA site data derived from sequences from diverse 3' sequencing protocols and biological samples in plants. Currently, PlantAPAdb contains APA sites in six species, Oryza sativa (japonica and indica), Arabidopsis (Arabidopsis thaliana), Medicago truncatula, Trifolium pratense, Phyllostachys edulis, and Chlamydomonas reinhardtii APA sites in PlantAPAdb are available for bulk download and can be queried in a Google-like manner. PlantAPAdb provides rich information of the whole-genome APA sites, including genomic locations, heterogeneous cleavage sites, expression levels, and sample information. It also provides comprehensive poly(A) signals for APA sites in different genomic regions according to distinct profiles of cis-elements in plants. In addition, PlantAPAdb contains events of 3' untranslated region shortening/lengthening resulting from APA, which helps to understand the mechanisms underlying systematic changes in 3' untranslated region lengths. Additional information about conservation of APA sites in plants is also available, providing insights into the evolutionary polyadenylation configuration across species. As a user-friendly database, PlantAPAdb is a large and extendable resource for elucidating APA mechanisms, APA conservation, and gene expression regulation.

hnRNPC regulates cancer-specific alternative cleavage and polyadenylation profiles

Alternative cleavage and polyadenylation (APA) can occur at more than half of all human genes, greatly enhancing the cellular repertoire of mRNA isoforms. As these isoforms can have altered stability, localisation and coding potential, deregulation of APA can disrupt gene expression and this has been linked to many diseases including cancer progression. How APA generates cancer-specific isoform profiles and what their physiological consequences are, however, is largely unclear. Here we use a subcellular fractionation approach to determine the nuclear and cytoplasmic APA profiles of successive stages of colon cancer using a cell line-based model. Using this approach, we show that during cancer progression specific APA profiles are established. We identify that overexpression of hnRNPC has a critical role in the establishment of APA profiles characteristic for metastatic colon cancer cells, by regulating poly(A) site selection in a subset of genes that have been implicated in cancer progression including MTHFD1L.

Genome-wide alternative polyadenylation dynamics in response to biotic and abiotic stresses in rice

Alternative polyadenylation (APA) is an important way to regulate gene expression at the post-transcriptional level, and is extensively involved in plant stress responses. However, the systematic roles of APA regulation in response to abiotic and biotic stresses in rice at the genome scale remain unknown. To take advantage of available RNA-seq datasets, using a novel tool APAtrap, we identified thousands of genes with significantly differential usage of polyadenylation [poly(A)] sites in response to the abiotic stress (drought, heat shock, and cadmium) and biotic stress [bacterial blight (BB), rice blast, and rice stripe virus (RSV)]. Genes with stress-responsive APA dynamics commonly exhibited higher expression levels when their isoforms with short 3' untranslated region (3' UTR) were more abundant. The stress-responsive APA events were widely involved in crucial stress-responsive genes and pathways: e.g. APA acted as a negative regulator in heat stress tolerance; APA events were involved in DNA repair and cell wall formation under Cd stress; APA regulated chlorophyll metabolism, being associated with the pathogenesis of leaf diseases under RSV and BB challenges. Furthermore, APA events were found to be involved in glutathione metabolism and MAPK signaling pathways, mediating a crosstalk among the abiotic and biotic stress-responsive regulatory networks in rice. Analysis of large-scale datasets revealed that APA may regulate abiotic and biotic stress-responsive processes in rice. Such post-transcriptome diversities contribute to rice adaption to various environmental challenges. Our study would supply useful resource for further molecular assisted breeding of multiple stress-tolerant cultivars for rice.

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.

Transcriptome Analyses of FY Mutants Reveal Its Role in mRNA Alternative Polyadenylation

A crucial step for mRNA polyadenylation is poly(A) signal recognition by trans-acting factors. The mammalian cleavage and polyadenylation specificity factor (CPSF) complex components CPSF30 and WD repeat-containing protein33 (WDR33) recognize the canonical AAUAAA for polyadenylation. In Arabidopsis (Arabidopsis thaliana), the flowering time regulator FY is the homolog of WDR33. However, its role in mRNA polyadenylation is poorly understood. Using poly(A) tag sequencing, we found that >50% of alternative polyadenylation (APA) events are altered in fy single mutants or double mutants with oxt6 (a null mutant of AtCPSF30), but mutation of the FY WD40-repeat has a stronger effect than deletion of the plant-unique Pro-Pro-Leu-Pro-Pro (PPLPP) domain. fy mutations disrupt AAUAAA or AAUAAA-like poly(A) signal recognition. Notably, A-rich signal usage is suppressed in the WD40-repeat mutation but promoted in PPLPP-domain deficiency. However, fy mutations do not aggravate the altered signal usage in oxt6 Furthermore, the WD40-repeat mutation shows a preference for 3' untranslated region shortening, but the PPLPP-domain deficiency shows a preference for lengthening. Interestingly, the WD40-repeat mutant exhibits shortened primary roots and late flowering with alteration of APA of related genes. Importantly, the long transcripts of two APA genes affected in fy are related to abiotic stress responses. These results reveal a conserved and specific role of FY in mRNA polyadenylation.

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.

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.

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.

APAtrap: identification and quantification of alternative polyadenylation sites from RNA-seq data

Motivation

Alternative polyadenylation (APA) has been increasingly recognized as a crucial mechanism that contributes to transcriptome diversity and gene expression regulation. As RNA-seq has become a routine protocol for transcriptome analysis, it is of great interest to leverage such unprecedented collection of RNA-seq data by new computational methods to extract and quantify APA dynamics in these transcriptomes. However, research progress in this area has been relatively limited. Conventional methods rely on either transcript assembly to determine transcript 3' ends or annotated poly(A) sites. Moreover, they can neither identify more than two poly(A) sites in a gene nor detect dynamic APA site usage considering more than two poly(A) sites.

Results

We developed an approach called APAtrap based on the mean squared error model to identify and quantify APA sites from RNA-seq data. APAtrap is capable of identifying novel 3' UTRs and 3' UTR extensions, which contributes to locating potential poly(A) sites in previously overlooked regions and improving genome annotations. APAtrap also aims to tally all potential poly(A) sites and detect genes with differential APA site usages between conditions. Extensive comparisons of APAtrap with two other latest methods, ChangePoint and DaPars, using various RNA-seq datasets from simulation studies, human and Arabidopsis demonstrate the efficacy and flexibility of APAtrap for any organisms with an annotated genome.

Availability and implementation

Freely available for download at https://apatrap.sourceforge.io.

Contact

liqq@xmu.edu.cn or xhuister@xmu.edu.cn.

Supplementary information

Supplementary data are available at Bioinformatics online.

Alternative polyadenylation of single cells delineates cell types and serves as a prognostic marker in early stage breast cancer

Alternative polyadenylation (APA) in 3’ untranslated regions (3’ UTR) plays an important role in regulating transcript abundance, localization, and interaction with microRNAs. Length-variation of 3’UTRs by APA contributes to efficient proliferation of cancer cells. In this study, we investigated APA in single cancer cells and tumor microenvironment cells to understand the physiological implication of APA in different cell types. We analyzed APA patterns and the expression level of genes from the 515 single-cell RNA sequencing (scRNA-seq) dataset from 11 breast cancer patients. Although the overall 3’UTR length of individual genes was distributed equally in tumor and non-tumor cells, we found a differential pattern of polyadenylation in gene sets between tumor and non-tumor cells. In addition, we found a differential pattern of APA across tumor types using scRNA-seq data from 3 glioblastoma patients and 1 renal cell carcinoma patients. In detail, 1,176 gene sets and 53 genes showed the distinct pattern of 3’UTR shortening and over-expression as signatures for five cell types including B lymphocytes, T lymphocytes, myeloid cells, stromal cells, and breast cancer cells. Functional categories of gene sets for cellular proliferation demonstrated concordant regulation of APA and gene expression specific to cell types. The expression of APA genes in breast cancer was significantly correlated with the clinical outcome of earlier stage breast cancer patients. We identified cell type-specific APA in single cells, which allows the identification of cell types based on 3’UTR length variation in combination with gene expression. Specifically, an immune-specific APA signature in breast cancer could be utilized as a prognostic marker of early stage breast cancer.

Nudt21 regulates the alternative polyadenylation of Pak1 and is predictive in the prognosis of glioblastoma patients

Alternative polyadenylation (APA) has emerged as a prevalent feature associated with cancer development and progression. The advantage of APA to tumor progression is to induce oncogenes through 3'-UTR shortening, and to inactivate tumor suppressor genes via the re-routing of microRNA competition. We previously identified the Mammalian Cleavage Factor I-25 (CFIm25) (encoded by Nudt21 gene) as a master APA regulator whose expression levels directly impact the tumorigenicity of glioblastoma (GBM) in vitro and in vivo. Despite its importance, the role of Nudt21 in GBM development is not known and the genes subject to Nudt21 APA regulation that contribute to GBM progression have not been identified. Here, we find that Nudt21 is reduced in low grade glioma (LGG) and all four subtypes of high grade glioma (GBM). Reduced expression of Nudt21 associates with worse survival in TCGA LGG cohorts and two TCGA GBM cohorts. Moreover, although CFIm25 was initially identified as biochemically associated with both CFIm59 and CFIm68, we observed three CFIm distinct subcomplexes exist and CFIm59 protein level is dependent on Nudt21 expression in GBM cells, but CFIm68 is not, and that only CFIm59 predicts prognosis of GBM patients similar to Nudt21. Through the use of Poly(A)-Click-Seq to characterize APA, we define the mRNAs subject to 3'-UTR shortening upon Nudt21 depletion in GBM cells and observed enrichment in genes important in the Ras signaling pathway, including Pak1. Remarkably, we find that Pak1 expression is regulated by Nudt21 through its 3'-UTR APA, and the combination of Pak1 and Nudt21 expression generates an even stronger prognostic indicator of GBM survival versus either value used alone. Collectively, our data uncover Nudt21 and its downstream target Pak1 as a potential "combination biomarker" for predicting prognosis of GBM patients.

Differential alternative polyadenylation contributes to the developmental divergence between two rice subspecies, japonica and indica

Alternative polyadenylation (APA) is a widespread post-transcriptional mechanism that regulates gene expression through mRNA metabolism, playing a pivotal role in modulating phenotypic traits in rice (Oryza sativa L.). However, little is known about the APA-mediated regulation underlying the distinct characteristics between two major rice subspecies, indica and japonica. Using a poly(A)-tag sequencing approach, polyadenylation (poly(A)) site profiles were investigated and compared pairwise from germination to the mature stage between indica and japonica, and extensive differentiation in APA profiles was detected genome-wide. Genes with subspecies-specific poly(A) sites were found to contribute to subspecies characteristics, particularly in disease resistance of indica and cold-stress tolerance of japonica. In most tissues, differential usage of APA sites exhibited an apparent impact on the gene expression profiles between subspecies, and genes with those APA sites were significantly enriched in quantitative trait loci (QTL) related to yield traits, such as spikelet number and 1000-seed weight. In leaves of the booting stage, APA site-switching genes displayed global shortening of 3' untranslated regions with increased expression in indica compared with japonica, and they were overrepresented in the porphyrin and chlorophyll metabolism pathways. This phenomenon may lead to a higher chlorophyll content and photosynthesis in indica than in japonica, being associated with their differential growth rates and yield potentials. We further constructed an online resource for querying and visualizing the poly(A) atlas in these two rice subspecies. Our results suggest that APA may be largely involved in developmental differentiations between two rice subspecies, especially in leaf characteristics and the stress response, broadening our knowledge of the post-transcriptional genetic basis underlying the divergence of rice traits.

Detection of Differentially Expressed Cleavage Site Intervals Within 3' Untranslated Regions Using CSI-UTR Reveals Regulated Interaction Motifs

The length of untranslated regions at the 3' end of transcripts (3'UTRs) is regulated by alternate polyadenylation (APA). 3'UTRs contain regions that harbor binding motifs for regulatory molecules. However, the mechanisms that coordinate the 3'UTR length of specific groups of transcripts are not well-understood. We therefore developed a method, CSI-UTR, that models 3'UTR structure as tandem segments between functional alternative-polyadenylation sites (termed cleavage site intervals-CSIs). This approach facilitated (1) profiling of 3'UTR isoform expression changes and (2) statistical enrichment of putative regulatory motifs. CSI-UTR analysis is UTR-annotation independent and can interrogate legacy data generated from standard RNA-Seq libraries. CSI-UTR identified a set of CSIs in human and rodent transcriptomes. Analysis of RNA-Seq datasets from neural tissue identified differential expression events within 3'UTRs not detected by standard gene-based differential expression analyses. Further, in many instances 3'UTR and CDS from the same gene were regulated differently. This modulation of motifs for RNA-interacting molecules with potential condition-dependent and tissue-specific RNA binding partners near the polyA signal and CSI junction may play a mechanistic role in the specificity of alternative polyadenylation. Source code, CSI BED files and example datasets are available at: https://github.com/UofLBioinformatics/CSI-UTR.

The mRNA Export Receptor NXF1 Coordinates Transcriptional Dynamics, Alternative Polyadenylation, and mRNA Export

Alternative polyadenylation (APA) produces mRNA isoforms with different 3' UTR lengths. Previous studies indicated that 3' end processing and mRNA export are intertwined in gene regulation. Here, we show that mRNA export factors generally facilitate usage of distal cleavage and polyadenylation sites (PASs), leading to long 3' UTR isoform expression. By focusing on the export receptor NXF1, which exhibits the most potent effect on APA in this study, we reveal several gene features that impact NXF1-dependent APA, including 3' UTR size, gene size, and AT content. Surprisingly, NXF1 downregulation results in RNA polymerase II (Pol II) accumulation at the 3' end of genes, correlating with its role in APA regulation. Moreover, NXF1 cooperates with CFI-68 to facilitate nuclear export of long 3' UTR isoform with UGUA motifs. Together, our work reveals important roles of NXF1 in coordinating transcriptional dynamics, 3' end processing, and nuclear export of long 3' UTR transcripts, implicating NXF1 as a nexus of gene regulation.

Modulation of Auxin Signaling and Development by Polyadenylation Machinery

Polyadenylation influences gene expression by affecting mRNA stability, transport, and translatability. Here, we report that Cleavage stimulation Factor 77 (AtCstF77), a component of the pre-mRNA 3'-end polyadenylation machinery, affects polyadenylation site (PAS) selection in transcripts of some auxin signaling genes in Arabidopsis (Arabidopsis thaliana). Disruption of AtCstF77 reduced auxin sensitivity and decreased the expression of the auxin reporter DR5-GFP Null mutations of cstf77 caused severe developmental defects, but were not lethal as previously reported. cstf77-2 genetically interacted with transport inhibitor response 1 auxin signaling f-box 2 auxin receptor double mutants, further supporting that polyadenylation affects auxin signaling. AtCstF77 was ubiquitously expressed in embryos, seedlings, and adult plants. The AtCstF77 protein was localized in the nucleus, which is consistent with its function in pre-mRNA processing. We observed that PASs in transcripts from approximately 2,400 genes were shifted in the cstf77-2 mutant. Moreover, most of the PAS shifts were from proximal to distal sites. Auxin treatment also caused PAS shifts in transcripts from a small number of genes. Several auxin signaling or homeostasis genes had different PASs in their transcripts in the cstf77-2 mutant. The expression levels of AUXIN RESISTANT 2/INDOLE-3-ACETIC ACID 7 were significantly increased in the cstf77-2 mutant, which can partially account for the auxin resistance phenotype of this mutant. Our results demonstrate that AtCstF77 plays pleiotropic and critical roles in Arabidopsis development. Moreover, disruption of AtCstF64, another component of the polyadenylation machinery, led to developmental defects and reduced auxin response, similar to those of the cstf77-2 mutant. We conclude that AtCstF77 affects auxin responses, likely by controlling PAS selection of transcripts of some auxin signaling components.

The structural basis of CstF-77 modulation of cleavage and polyadenylation through stimulation of CstF-64 activity

Cleavage and polyadenylation (C/P) of mRNA is an important cellular process that promotes increased diversity of mRNA isoforms and could change their stability in different cell types. The cleavage stimulation factor (CstF) complex, part of the C/P machinery, binds to U- and GU-rich sequences located downstream from the cleavage site through its RNA-binding subunit, CstF-64. Less is known about the function of the other two subunits of CstF, CstF-77 and CstF-50. Here, we show that the carboxy-terminus of CstF-77 plays a previously unrecognized role in enhancing C/P by altering how the RNA recognition motif (RRM) of CstF-64 binds RNA. In support of this finding, we also show that CstF-64 relies on CstF-77 to be transported to the nucleus; excess CstF-64 localizes to the cytoplasm, possibly via interaction with cytoplasmic RNAs. Reverse genetics and nuclear magnetic resonance studies of recombinant CstF-64 (RRM-Hinge) and CstF-77 (monkeytail-carboxy-terminal domain) indicate that the last 30 amino acids of CstF-77 increases the stability of the RRM, thus altering the affinity of the complex for RNA. These results provide new insights into the mechanism by which CstF regulates the location of the RNA cleavage site during C/P.

Dual Role of a C-Terminally Truncated Isoform of Large Tumor Suppressor Kinase 1 in the Regulation of Hippo Signaling and Tissue Growth

The considerable amount of experimental evidence has defined the Hippo pathway as a tumor suppressive pathway and increased expression and/or activity of its oncogenic effectors is frequently observed in cancer. However, clinical studies have failed to attribute cancer development and progression to mutations in the pathway. In explaining this conundrum, we investigated the expression and functions of a C-terminally truncated isoform of large tumor suppressor kinase 1 (LATS1) called short LATS1 (sLATS1) in human cell lines and Drosophila. Intriguingly, through overexpression of sLATS1, we demonstrated that sLATS1 either activates or suppresses the activity of Yes-associated protein (YAP), one of the effectors of the Hippo pathway, in a cell type-specific manner. The activation is mediated through inhibition of full-length LATS1, whereas suppression of YAP is accomplished through sLATS1-YAP interaction. In HEK293T cells, the former mechanism may affect the cellular response more dominantly, whereas in U2OS cells and developing tissues in Drosophila, the latter mechanism may be solely carried out. Finally, to find the clinical relevance of this molecule, we examined the expression of sLATS1 in breast cancer patients. The transcriptome analysis showed that the ratio of sLATS1 to LATS1 was increased in tumor tissues comparing to their adjacent normal tissues.

Alternative polyadenylation factors link cell cycle to migration

BACKGROUND

In response to a wound, fibroblasts are activated to migrate toward the wound, to proliferate and to contribute to the wound healing process. We hypothesize that changes in pre-mRNA processing occurring as fibroblasts enter the proliferative cell cycle are also important for promoting their migration.

RESULTS

RNA sequencing of fibroblasts induced into quiescence by contact inhibition reveals downregulation of genes involved in mRNA processing, including splicing and cleavage and polyadenylation factors. These genes also show differential exon use, especially increased intron retention in quiescent fibroblasts compared to proliferating fibroblasts. Mapping the 3' ends of transcripts reveals that longer transcripts from distal polyadenylation sites are more prevalent in quiescent fibroblasts and are associated with increased expression and transcript stabilization based on genome-wide transcript decay analysis. Analysis of dermal excisional wounds in mice reveals that proliferating cells adjacent to wounds express higher levels of cleavage and polyadenylation factors than quiescent fibroblasts in unwounded skin. Quiescent fibroblasts contain reduced levels of the cleavage and polyadenylation factor CstF-64. CstF-64 knockdown recapitulates changes in isoform selection and gene expression associated with quiescence, and results in slower migration.

CONCLUSIONS

Our findings support cleavage and polyadenylation factors as a link between cellular proliferation state and migration.

Differential genome-wide profiling of alternative polyadenylation sites in nasopharyngeal carcinoma by high-throughput sequencing

Background

Alternative polyadenylation (APA) is a widespread phenomenon in the posttranscriptional regulation of gene expression that generates mRNAs with alternative 3′-untranslated regions (3’UTRs). APA contributes to the pathogenesis of various diseases, including cancer. However, the potential role of APA in the development of nasopharyngeal carcinoma (NPC) remains largely unknown.

Methods

A strategy of sequencing APA sites (SAPAS) based on second-generation sequencing technology was carried out to explore the global patterns of APA sites and identify genes with tandem 3’UTRs in samples from 6 NPC and 6 normal nasopharyngeal epithelial tissue (NNET). Sequencing results were then validated using quantitative RT-PCR in a larger cohort of 16 NPC and 16 NNET samples.

Results

The sequencing data showed that the use of tandem APA sites was prevalent in NPC, and numerous genes with APA-switching events were discovered. In total, we identified 195 genes with significant differences in the tandem 3’UTR length between NPC and NNET; including 119 genes switching to distal poly (A) sites and 76 genes switching to proximal poly (A) sites. Several gene ontology (GO) terms were enriched in the list of genes with switched APA sites, including regulation of cell migration, macromolecule catabolic process, protein catabolic process, proteolysis, small conjugating protein ligase activity, and ubiquitin-protein ligase activity.

Conclusions

APA site-switching events are prevalent in NPC. APA-mediated regulation of gene expression may play an important role in the development of NPC, and more detailed studies targeting genes with APA-switching events may contribute to the development of novel future therapeutic strategies for NPC.

Electronic supplementary material

The online version of this article (10.1186/s12929-018-0477-6) contains supplementary material, which is available to authorized users.

Crosstalk between alternative polyadenylation and miRNAs in the regulation of protein translational efficiency

3′ UTRs play important roles in the gene regulation network via their influence on mRNA stability, translational efficiency, and subcellular localization. For a given gene, 3′ UTRs of different lengths generated by alternative polyadenylation (APA) may result in functional differences in regulation. The mechanistic details of how length changes of 3′ UTRs alter gene function remain unclear. By combining APA sequencing and polysome profiling, we observed that mRNA isoforms with shorter 3′ UTRs were bound with more polysomes in six cell lines but not in NIH3T3 cells, suggesting that changing 3′ UTRs to shorter isoforms may lead to a higher gene translational efficiency. By interfering with the expression of TNRC6A and analyzing AGO2-PAR-CLIP data, we revealed that the APA effect on translational efficiency was mainly regulated by miRNAs, and this regulation was cell cycle dependent. The discrepancy between NIH3T3 and other cell lines was due to contact inhibition of NIH3T3. Thus, the crosstalk between APA and miRNAs may be needed for the regulation of protein translational efficiency.

Genome-wide profiling reveals cancer-related genes with switched alternative polyadenylation sites in colorectal cancer

Background

Alternative polyadenylation (APA) is an important post-transcriptional regulation in eukaryotic cells. It plays considerable roles in many biological processes and diseases, such as cell differentiation, proliferation and cancer. Colorectal cancer (CRC) is one of the most common malignancies worldwide, which is among the top five in incidence and mortality of all cancers in China. Although there have been some studies on the APA of CRC, the normal and carcinoma samples used for genome-wide profiling were not matched. The purpose of this study was to obtain genes with switched 3′-untranslated region (UTR) that may be associated with intracellular regulation of CRC by analyzing APA patterns of strict control groups from clinical patients.

Materials and methods

CRC and matched normal tissues were acquired from surgical specimens from three CRC patients. Their libraries of 3′-terminal fragments of mRNA with poly(A) tails were constructed by 3T-seq technology and sequenced by Illumina Hiseq X Ten. APA patterns of cancer and matched normal tissues were analyzed by bioinformatics analysis, and a representative gene, GPI, was verified by quantitative reverse transcription PCR.

Results

Overall, we identified 35,076 poly(A) sites in total. Compared to the matched normal tissues, we detected 350, 405 and 375 genes with significantly APA-mediated 3′-UTR alteration in cancer tissues of three patients, respectively. Forty-seven genes with switched 3′-UTR were shared in all three patients. In addition, most of these genes have shortened 3′-UTRs, some of which were associated with cancers, such as GPI.

Conclusion

Our studies found several genes with switched 3′-UTR in CRC patients, which may provide some important clues for more in-depth study of the cellular regulation in CRC from the perspective of post-transcriptional regulation. It may also help in the search for new biomarkers of CRC.

Integrative 3' Untranslated Region-Based Model to Identify Patients with Low Risk of Axillary Lymph Node Metastasis in Operable Triple-Negative Breast Cancer

BACKGROUND

Sentinel lymph node biopsy is the standard surgical staging approach for operable triple-negative breast cancer (TNBC) with clinically negative axillae. In this study, we sought to develop a model to predict TNBC patients with negative nodal involvement, who would benefit from the exemption of the axillary staging surgery.

MATERIALS AND METHODS

We evaluated 3' untranslated region (3'UTR) profiles using microarray data of TNBC from two Gene Expression Omnibus datasets. Samples from GSE31519 were divided into training set (n = 164) and validation set (n = 163), and GSE76275 was used to construct testing set (n = 164). We built a six-member 3'UTR panel (ADD2, COL1A1, APOL2, IL21R, PKP2, and EIF4G3) using an elastic net model to estimate the risk of lymph node metastasis (LNM). Receiver operating characteristic and logistic analyses were used to assess the association between the panel and LNM status.

RESULTS

The six-member 3'UTR-panel showed a high distinguishing power with an area under the curve of 0.712, 0.729, and 0.708 in the training, validation, and testing sets, respectively. After adjustment by tumor size, the 3'UTR panel retained significant predictive power in the training, validation, and testing sets (odds ratio = 4.93, 4.58, and 3.59, respectively; p < .05 for all). A combinatorial analysis of the 3'UTR panel and tumor size yielded an accuracy of 97.2%, 100%, and 100% in training, validation, and testing set, respectively.

CONCLUSION

This study established an integrative 3'UTR-based model as a promising predictor for nodal negativity in operable TNBC. Although a prospective study is needed to validate the model, our results may permit a no axillary surgery option for selected patients.

IMPLICATIONS FOR PRACTICE

Currently, sentinel lymph node biopsy is the standard approach for surgical staging in breast cancer patients with negative axillae. Prediction estimation for lymph node metastasis of breast cancer relies on clinicopathological characteristics, which is unreliable, especially in triple-negative breast cancer (TNBC)-a highly heterogeneous disease. The authors developed and validated an effective prediction model for the lymph node status of patients with TNBC, which integrates 3'UTR markers and tumor size. This is the first 3'UTR-based model that will help identify TNBC patients with low risk of nodal involvement who are most likely to benefit from exemption axillary surgery.

Androgen-dependent alternative mRNA isoform expression in prostate cancer cells

Background: Androgen steroid hormones are key drivers of prostate cancer. Previous work has shown that androgens can drive the expression of alternative mRNA isoforms as well as transcriptional changes in prostate cancer cells. Yet to what extent androgens control alternative mRNA isoforms and how these are expressed and differentially regulated in prostate tumours is unknown. Methods: Here we have used RNA-Seq data to globally identify alternative mRNA isoform expression under androgen control in prostate cancer cells, and profiled the expression of these mRNA isoforms in clinical tissue. Results: Our data indicate androgens primarily switch mRNA isoforms through alternative promoter selection. We detected 73 androgen regulated alternative transcription events, including utilisation of 56 androgen-dependent alternative promoters, 13 androgen-regulated alternative splicing events, and selection of 4 androgen-regulated alternative 3' mRNA ends. 64 of these events are novel to this study, and 26 involve previously unannotated isoforms. We validated androgen dependent regulation of 17 alternative isoforms by quantitative PCR in an independent sample set. Some of the identified mRNA isoforms are in genes already implicated in prostate cancer (including LIG4, FDFT1 and RELAXIN), or in genes important in other cancers (e.g. NUP93 and MAT2A). Importantly, analysis of transcriptome data from 497 tumour samples in the TGCA prostate adenocarcinoma (PRAD) cohort identified 13 mRNA isoforms (including TPD52, TACC2 and NDUFV3) that are differentially regulated in localised prostate cancer relative to normal tissue, and 3 ( OSBPL1A, CLK3 and TSC22D3) which change significantly with Gleason grade and  tumour stage. Conclusions: Our findings dramatically increase the number of known androgen regulated isoforms in prostate cancer, and indicate a highly complex response to androgens in prostate cancer cells that could be clinically important.

Tandem 3' UTR Patterns and Gene Expression Profiles of Marc-145 Cells During PRRSV Infection

Porcine reproductive and respiratory syndrome virus (PRRSV) causes substantial economic losses to the global pig industry. Alternative polyadenylation (APA) is a mechanism that diversifies gene expression, which is important for tumorigenesis, development, and cell differentiation. However, it is unclear whether APA plays a role in the course of PRRSV infection. To address this issue, in this study we carried out a whole-genome transcriptome analysis of PRRSV-infected Marc-145 African green monkey kidney cells and identified 185 APA switching genes and 393 differentially expressed genes (DEGs). Most of these genes were involved in cellular process, metabolism, and biological regulation, and there was some overlap between the two gene sets. DEGs were found to be more directly involved in the antiviral response than APA genes. These findings provide insight into the dynamics of host gene regulation during PRRSV infection and a basis for elucidating the pathogenesis of PRRSV.

Alternative Polyadenylation of Mammalian Transcripts Is Generally Deleterious, Not Adaptive

Alternative polyadenylation (APA) produces from the same gene multiple mature RNAs with varying 3' ends. Although APA is commonly believed to generate beneficial functional diversity and be adaptive, we hypothesize that most genes have one optimal polyadenylation site and that APA is caused largely by deleterious polyadenylation errors. The error hypothesis, but not the adaptive hypothesis, predicts that, as the expression level of a gene increases, its polyadenylation diversity declines, relative use of the major (presumably optimal) polyadenylation site increases, and that of each minor (presumably nonoptimal) site decreases. It further predicts that the number of polyadenylation signals per gene is smaller than the random expectation and that polyadenylation signals for major but not minor sites are under purifying selection. All of these predictions are confirmed in mammals, suggesting that numerous defective RNAs are produced in normal cells, many phenotypic variations at the molecular level are nonadaptive, and cellular life is noisier than is appreciated.

NUDT21 negatively regulates PSMB2 and CXXC5 by alternative polyadenylation and contributes to hepatocellular carcinoma suppression

Alternative polyadenylation (APA) is an important post-transcriptional regulatory mechanism and involved in many diseases, including cancer. CFIm25, a subunit of the cleavage factor I encoded by NUDT21, is required for 3'RNA cleavage and polyadenylation. Although it has been recently reported to be involved in glioblastoma tumor suppression, its roles and the underlying functional mechanism remain unclear in other types of cancer. In this study, we characterized NUDT21 in hepatocellular carcinoma (HCC). Reduced expression of NUDT21 was observed in HCC tissue compared to adjacent non-tumorous compartment. HCC patients with lower NUDT21 expression have shorter overall and disease-free survival times than those with higher NUDT21 expression after surgery. Knockdown of NUDT21 promotes HCC cell proliferation, metastasis, and tumorigenesis, whereas forced expression of NUDT21 exhibits the opposite effects. We then performed global APA site profiling analysis in HCC cells and identified considerable number of genes with shortened 3'UTRs upon the modulation of NUDT21 expression. In particular, we further characterized the NUDT21-regulated genes PSMB2 and CXXC5. We found NUDT21 knockdown increases usage of the proximal polyadenylation site in the PSMB2 and CXXC5 3' UTRs, resulting in marked increase in the expression of PSMB2 and CXXC5. Moreover, knockdown of PSMB2 or CXXC5 suppresses HCC cell proliferation and invasion. Taken together, our study demonstrated that NUDT21 inhibits HCC proliferation, metastasis and tumorigenesis, at least in part, by suppressing PSMB2 and CXXC5, and thereby provided a new insight into understanding the connection of HCC suppression and APA machinery.

3' UTR shortening represses tumor-suppressor genes in trans by disrupting ceRNA crosstalk

Widespread mRNA 3' UTR shortening through alternative polyadenylation ¹ promotes tumor growth in vivo ² . A prevailing hypothesis is that it induces proto-oncogene expression in cis through escaping microRNA-mediated repression. Here we report a surprising enrichment of 3'UTR shortening among transcripts that are predicted to act as competing-endogenous RNAs (ceRNAs) for tumor-suppressor genes. Our model-based analysis of the trans effect of 3' UTR shortening (MAT3UTR) reveals a significant role in altering ceRNA expression. MAT3UTR predicts many trans-targets of 3' UTR shortening, including PTEN, a crucial tumor-suppressor gene ³ involved in ceRNA crosstalk ⁴ with nine 3'UTR-shortening genes, including EPS15 and NFIA. Knockdown of NUDT21, a master 3' UTR-shortening regulator ² , represses tumor-suppressor genes such as PHF6 and LARP1 in trans in a miRNA-dependent manner. Together, the results of our analysis suggest a major role of 3' UTR shortening in repressing tumor-suppressor genes in trans by disrupting ceRNA crosstalk, rather than inducing proto-oncogenes in cis.

Cell Cycle Regulation by Alternative Polyadenylation of CCND1

Global shortening of 3′UTRs by alternative polyadenylation (APA) has been observed in cancer cells. However, the role of APA in cancer remains unknown. CCND1 is a proto-oncogene that regulates progression through the G1-S phase of the cell cycle; moreover, it has been observed to be switching to proximal APA sites in cancer cells. To investigate the biological function of the APA of CCND1, we edited the weak poly(A) signal (PAS) of the proximal APA site to a canonical PAS using the CRISPR/Cas9 method, which can force the cells to use a proximal APA site. Cell cycle profiling and proliferation assays revealed that the proximal APA sites of CCND1 accelerated the cell cycle and promoted cell proliferation, but UTR-APA and CR-APA act via different molecular mechanisms. These results indicate that PAS editing with CRISPR/Cas9 provides a good method by which to study the biological function of APA.

Widespread intronic polyadenylation diversifies immune cell transcriptomes

Alternative cleavage and polyadenylation (ApA) is known to alter untranslated region (3'UTR) length but can also recognize intronic polyadenylation (IpA) signals to generate transcripts that lose part or all of the coding region. We analyzed 46 3'-seq and RNA-seq profiles from normal human tissues, primary immune cells, and multiple myeloma (MM) samples and created an atlas of 4927 high-confidence IpA events represented in these cell types. IpA isoforms are widely expressed in immune cells, differentially used during B-cell development or in different cellular environments, and can generate truncated proteins lacking C-terminal functional domains. This can mimic ectodomain shedding through loss of transmembrane domains or alter the binding specificity of proteins with DNA-binding or protein-protein interaction domains. MM cells display a striking loss of IpA isoforms expressed in plasma cells, associated with shorter progression-free survival and impacting key genes in MM biology and response to lenalidomide.

MicroRNA-143 Targets ATG2B to Inhibit Autophagy and Increase Inflammatory Responses in Crohn's Disease

BACKGROUND

Dysfunctional autophagy is recognized as a contributing factor in many chronic inflammatory diseases, including Crohn's disease (CD). Genetic analyses have found that microRNA (miRNA) levels are altered in the intestinal tissues of CD patients.

METHODS

The Sequencing Alternative Poly-Adenylation Sites (SAPAS) method was used to compare the 3' end of the total mRNA sequence of 3 surgical specimens of CD patients (including inflamed tissues and corresponding noninflamed tissues in each case). The levels of autophagy-related 2B (ATG2B), LC3, and miR-143 were compared between inflamed tissues and noninflamed tissues using immunoblot and quantitative reverse transcription polymerase chain reaction. Luciferase assays were used to verify the interactions between miR-143 and ATG2B. Autophagy was measured by immunoblot analyses of LC3 and transmission electron microscopy. Inflammatory cytokines and IκBα were analyzed to evaluate the effect of miR-143 on inflammatory response.

RESULTS

The tandem repeat 3'-UTR of ATG2B was longer in inflamed tissues than in corresponding noninflamed tissues and contained an miR-143 target site. miR-143 expression was elevated, whereas ATG2B and LC3-II were downregulated in inflamed tissues. The direct interaction between miR-143 and ATG2B was verified by a 3'-UTR dual-luciferase reporter assay. Constitutive expression of miR-143 or depletion of ATG2B in cultured intestinal epithelial cells inhibited autophagy, reduced IκBα levels, and increased inflammatory responses.

CONCLUSIONS

miR-143 may induce bowel inflammation by regulating ATG2B and autophagy, suggesting that miR-143 might play a critical role in the development of CD. Therefore, miR-143 could be a promising novel target for gene therapy in CD patients.

3' UTR lengthening as a novel mechanism in regulating cellular senescence

Cellular senescence has been viewed as a tumor suppression mechanism and also as a contributor to individual aging. Widespread shortening of 3' untranslated regions (3' UTRs) in messenger RNAs (mRNAs) by alternative polyadenylation (APA) has recently been discovered in cancer cells. However, the role of APA in the process of cellular senescence remains elusive. Here, we found that hundreds of genes in senescent cells tended to use distal poly(A) (pA) sites, leading to a global lengthening of 3' UTRs and reduced gene expression. Genes that harbor longer 3' UTRs in senescent cells were enriched in senescence-related pathways. Rras2, a member of the Ras superfamily that participates in multiple signal transduction pathways, preferred longer 3' UTR usage and exhibited decreased expression in senescent cells. Depletion of Rras2 promoted senescence, while rescue of Rras2 reversed senescence-associated phenotypes. Mechanistically, splicing factor TRA2B bound to a core "AGAA" motif located in the alternative 3' UTR of Rras2, thereby reducing the RRAS2 protein level and causing senescence. Both proximal and distal poly(A) signals showed strong sequence conservation, highlighting the vital role of APA regulation during evolution. Our results revealed APA as a novel mechanism in regulating cellular senescence.

CFIm25 inhibits hepatocellular carcinoma metastasis by suppressing the p38 and JNK/c-Jun signaling pathways

Alternative polyadenylation (APA), a post-transcriptional modification, has been implicated in many diseases, but especially in tumor proliferation. CFIm25, the 25 kDa subunit of human cleavage factor Im (CFIm), is a key factor in APA. We show that CFIm25 expression is reduced in human hepatocellular carcinoma (HCC), and its expression correlates with metastasis. Kaplan-Meier analysis indicated that CFIm25 is related to overall survival in HCC. Moreover, CFIm25 expression is negatively related to the metastatic potential of HCC cell lines. CFIm25 knockdown promotes cell invasion and migration in vitro, while overexpression of CFIm25 inhibits cell invasion and migration in vitro and inhibits intrahepatic and lung metastasis in vivo. Additional studies showed that CFIm25 disrupts epithelial-mesenchymal transition by increasing E-cadherin, that it inhibits HCC cell migration and invasion by blocking the p38 and JNK/c-Jun signaling pathways, and that CFIm25 knockdown increases the transcriptional activity of activating protein-1 (AP-1). These findings indicate that therapy directed at increasing CFIm25 expression is a potential HCC treatment.

The 3'UTR signature defines a highly metastatic subgroup of triple-negative breast cancer

Triple-negative breast cancer (TNBC) is a highly heterogeneous disease with an aggressive clinical course. Prognostic models are needed to chart potential patient outcomes. To address this, we used alternative 3′UTR patterns to improve postoperative risk stratification. We collected 327 publicly available microarrays and generated the 3′UTR landscape based on expression ratios of alternative 3′UTR. After initial feature filtering, we built a 17-3′UTR-based classifier using an elastic net model. Time-dependent ROC comparisons and Kaplan–Meier analyses confirmed an outstanding discriminating power of our prognostic model for TNBC patients. In the training cohort, 5-year event-free survival (EFS) was 78.6% (95% CI 71.2–86.0) for the low-risk group, and 16.3% (95% CI 2.3–30.4) for the high-risk group (log-rank p<0.0001; hazard ratio [HR] 8.29, 95% CI 4.78–14.4), In the validation set, 5-year EFS was 75.6% (95% CI 68.0–83.2) for the low-risk group, and 33.2% (95% CI 17.1–49.3) for the high-risk group (log-rank p<0.0001; HR 3.17, 95% CI 1.66–5.42). In conclusion, the 17-3′UTR-based classifier provides a superior prognostic performance for estimating disease recurrence and metastasis in TNBC patients and it may permit personalized management strategies.

MACC1 is post-transcriptionally regulated by miR-218 in colorectal cancer

Metastasis is a multistep molecular network process, which is lethal for more than 90% of the cancer patients. Understanding the regulatory functions of metastasis-inducing molecules is in high demand for improved therapeutic cancer approaches. Thus, we studied the post-transcriptional regulation of the crucial carcinogenic and metastasis-mediating molecule metastasis associated in colon cancer 1 (MACC1). In silico analysis revealed MACC1 as a potential target of miR-218, a tumor suppressor miRNA. Expression of these two molecules inversely correlated in colorectal cancer (CRC) cell lines. In a cohort of CRC patient tissues (n = 59), miR-218 is significantly downregulated and MACC1 is upregulated compared with normal mucosa. Luciferase reporter assays with a construct of the MACC1-3′-UTR harboring either the wild type or the mutated miR-218 seed sequence confirmed the specificity of the targeting. miR-218 inhibited significantly MACC1 protein expression, and consistently, MACC1-mediated migration, invasion and colony formation in CRC cells. Anti-miR-218 enhanced the MACC1-mediated migration, invasion and colony formation. Similar findings were observed in the gastric cancer cell line MKN-45. Further, we performed methylation-specific PCR of the SLIT2 and SLIT3 promoter, where miR-218 is encoded in intronic regions. The SLIT2 and SLIT3 promoters are hypermethylated in CRC cell lines. miR-218 and SLIT2 expressions correlated positively. Methyltransferase inhibitor 5-Azacytidine induced miR-218 expression and inhibited the expression of its target MACC1. We also determined that MACC1 has alternative polyadenylation (APA) sites, which results in different lengths of 3′-UTR variants in a CRC cell line. Taken together, miR-218 is post-transcriptionally inhibiting the MACC1 expression and its metastasis-inducing abilities.

Aberrant Ki-67 expression through 3'UTR alternative polyadenylation in breast cancers

Ki‐67 (MKI67) is a marker of cellular proliferation of cancer. Here, we show that Ki‐67 is post‐transcriptionally regulated through alternative polyadenylation (APA) and microRNAs in breast cancer. We show that shortening of the Ki‐67 3′UTR results in the loss of the binding sites for the suppressive miRNAs and thus renders the transcript with a shortened 3′UTR insusceptible to miRNA‐mediated suppression. This APA‐mediated shortening of the Ki‐67 3′UTR contributes to increased mRNA stability and enhanced translational efficiency. In summary, our results not only highlight the post‐transcriptional regulation of Ki‐67 involving APA and microRNAs but also suggest that Ki‐67 3′UTR disruption could serve as a molecular marker in breast cancer.

Analyses of alternative polyadenylation: from old school biochemistry to high-throughput technologies

Alternations in usage of polyadenylation sites during transcription termination yield transcript isoforms from a gene. Recent findings of transcriptome-wide alternative polyadenylation (APA) as a molecular response to changes in biology position APA not only as a molecular event of early transcriptional termination but also as a cellular regulatory step affecting various biological pathways. With the development of high-throughput profiling technologies at a single nucleotide level and their applications targeted to the 3'-end of mRNAs, dynamics in the landscape of mRNA 3'-end is measureable at a global scale. In this review, methods and technologies that have been adopted to study APA events are discussed. In addition, various bioinformatics algorithms for APA isoform analysis using publicly available RNA-seq datasets are introduced. [BMB Reports 2017; 50(4): 201-207].

Antisense transcription regulates the expression of sense gene via alternative polyadenylation

Natural antisense transcripts (NAT) and alternative polyadenylation (APA) of messenger RNA (mRNA) are important contributors of transcriptome complexity, each playing a critical role in multiple biological processes. However, whether they have crosstalk and function collaboratively is unclear. We discovered that APA enriched in human sense-antisense (S-AS) gene pairs, and finally focused on RNASEH2C-KAT5 S-AS pair for further study. In cis but not in trans over-expression of the antisense KAT5 gene promoted the usage of distal polyA (pA) site in sense gene RNASEH2C, which generated longer 3' untranslated region (3'UTR) and produced less protein, accompanying with slowed cell growth. Mechanistically, elevated Pol II occupancy coupled with SRSF3 could explain the higher usage of distal pA site. Finally, NAT-mediated downregulation of sense gene's protein level in RNASEH2C-KAT5 pair was specific for human rather than mouse, which lacks the distal pA site of RNASEH2C. We provided the first evidence to support that certain gene affected phenotype may not by the protein of its own, but by affecting the expression of its overlapped gene through APA, implying an unexpected view for understanding the link between genotype and phenotype.

Overexpressed HSF1 cancer signature genes cluster in human chromosome 8q

Background

HSF1 (heat shock factor 1) is a transcription factor that is found to facilitate malignant cancer development and proliferation. In cancer cells, HSF1 mediates a set of genes distinct from heat shock that contributes to malignancy. This set of genes is known as the HSF1 Cancer Signature genes or simply HSF1-CanSig genes. HSF1-CanSig genes function and operate differently than typical cancer-causing genes, yet it is involved in fundamental oncogenic processes.

Results

By utilizing expression data from 9241 cancer patients, we identified that human chromosome 8q21-24 is a location hotspot for the most frequently overexpressed HSF1-CanSig genes. Intriguingly, the strength of the HSF1 cancer program correlates with the number of overexpressed HSF1-CanSig genes in 8q, illuminating the essential role of HSF1 in mediating gene expression in different cancers. Chromosome 8q21-24 is found under selective pressure in preserving gene order as it exhibits strong synteny among human, mouse, rat, and bovine, although the biological significance remains unknown. Statistical modeling, hierarchical clustering, and gene ontology-based pathway analyses indicate crosstalk between HSF1-mediated responses and pre-mRNA 3′ processing in cancers.

Conclusions

Our results confirm the unique role of chromosome 8q mediated by the master regulator HSF1 in cancer cases. Additionally, this study highlights the connection between cellular processes triggered by HSF1 and pre-mRNA 3′ processing in cancers.

Electronic supplementary material

The online version of this article (10.1186/s40246-017-0131-5) contains supplementary material, which is available to authorized users.

Next-generation sequencing: recent applications to the analysis of colorectal cancer

Since the establishment of the Sanger sequencing method, scientists around the world focused their efforts to progress in the field to produce the utmost technology. The introduction of next-generation sequencing (NGS) represents a revolutionary step and promises to lead to massive improvements in our understanding on the role of nucleic acids functions. Cancer research began to use this innovative and highly performing method, and interesting results started to appear in colorectal cancer (CRC) analysis. Several studies produced high-quality data in terms of mutation discovery, especially about actionable or less frequently mutated genes, epigenetics, transcriptomics. Analysis of results is unveiling relevant perspectives aiding to evaluate the response to therapies. Novel evidences have been presented also in other directions such as gut microbiota or CRC circulating tumor cells. However, despite its unquestioned potential, NGS poses some issues calling for additional studies. This review intends to offer a view of the state of the art of NGS applications to CRC through examination of the most important technologies and discussion of recent published results.

Framework for reanalysis of publicly available Affymetrix® GeneChip® data sets based on functional regions of interest

BACKGROUND

Since the introduction of microarrays in 1995, researchers world-wide have used both commercial and custom-designed microarrays for understanding differential expression of transcribed genes. Public databases such as ArrayExpress and the Gene Expression Omnibus (GEO) have made millions of samples readily available. One main drawback to microarray data analysis involves the selection of probes to represent a specific transcript of interest, particularly in light of the fact that transcript-specific knowledge (notably alternative splicing) is dynamic in nature.

RESULTS

We therefore developed a framework for reannotating and reassigning probe groups for Affymetrix® GeneChip® technology based on functional regions of interest. This framework addresses three issues of Affymetrix® GeneChip® data analyses: removing nonspecific probes, updating probe target mapping based on the latest genome knowledge and grouping probes into gene, transcript and region-based (UTR, individual exon, CDS) probe sets. Updated gene and transcript probe sets provide more specific analysis results based on current genomic and transcriptomic knowledge. The framework selects unique probes, aligns them to gene annotations and generates a custom Chip Description File (CDF). The analysis reveals only 87% of the Affymetrix® GeneChip® HG-U133 Plus 2 probes uniquely align to the current hg38 human assembly without mismatches. We also tested new mappings on the publicly available data series using rat and human data from GSE48611 and GSE72551 obtained from GEO, and illustrate that functional grouping allows for the subtle detection of regions of interest likely to have phenotypical consequences.

CONCLUSION

Through reanalysis of the publicly available data series GSE48611 and GSE72551, we profiled the contribution of UTR and CDS regions to the gene expression levels globally. The comparison between region and gene based results indicated that the detected expressed genes by gene-based and region-based CDFs show high consistency and regions based results allows us to detection of changes in transcript formation.

Transcription elongation rate has a tissue-specific impact on alternative cleavage and polyadenylation in Drosophila melanogaster

Alternative polyadenylation (APA) is a mechanism that generates multiple mRNA isoforms with different 3'UTRs and/or coding sequences from a single gene. Here, using 3' region extraction and deep sequencing (3'READS), we have systematically mapped cleavage and polyadenylation sites (PASs) in Drosophila melanogaster, expanding the total repertoire of PASs previously identified for the species, especially those located in A-rich genomic sequences. Cis-element analysis revealed distinct sequence motifs around fly PASs when compared to mammalian ones, including the greater enrichment of upstream UAUA elements and the less prominent presence of downstream UGUG elements. We found that over 75% of mRNA genes in Drosophila melanogaster undergo APA. The head tissue tends to use distal PASs when compared to the body, leading to preferential expression of APA isoforms with long 3'UTRs as well as with distal terminal exons. The distance between the APA sites and intron location of PAS are important parameters for APA difference between body and head, suggesting distinct PAS selection contexts. APA analysis of the RpII215^C4 mutant strain, which harbors a mutant RNA polymerase II (RNAPII) with a slower elongation rate, revealed that a 50% decrease in transcriptional elongation rate leads to a mild trend of more usage of proximal, weaker PASs, both in 3'UTRs and in introns, consistent with the "first come, first served" model of APA regulation. However, this trend was not observed in the head, suggesting a different regulatory context in neuronal cells. Together, our data expand the PAS collection for Drosophila melanogaster and reveal a tissue-specific effect of APA regulation by RNAPII elongation rate.

Isoform Evolution in Primates through Independent Combination of Alternative RNA Processing Events

Recent RNA-seq technology revealed thousands of splicing events that are under rapid evolution in primates, whereas the reliability of these events, as well as their combination on the isoform level, have not been adequately addressed due to its limited sequencing length. Here, we performed comparative transcriptome analyses in human and rhesus macaque cerebellum using single molecule long-read sequencing (Iso-seq) and matched RNA-seq. Besides 359 million RNA-seq reads, 4,165,527 Iso-seq reads were generated with a mean length of 14,875 bp, covering 11,466 human genes, and 10,159 macaque genes. With Iso-seq data, we substantially expanded the repertoire of alternative RNA processing events in primates, and found that intron retention and alternative polyadenylation are surprisingly more prevalent in primates than previously estimated. We then investigated the combinatorial mode of these alternative events at the whole-transcript level, and found that the combination of these events is largely independent along the transcript, leading to thousands of novel isoforms missed by current annotations. Notably, these novel isoforms are selectively constrained in general, and 1,119 isoforms have even higher expression than the previously annotated major isoforms in human, indicating that the complexity of the human transcriptome is still significantly underestimated. Comparative transcriptome analysis further revealed 502 genes encoding selectively constrained, lineage-specific isoforms in human but not in rhesus macaque, linking them to some lineage-specific functions. Overall, we propose that the independent combination of alternative RNA processing events has contributed to complex isoform evolution in primates, which provides a new foundation for the study of phenotypic difference among primates.

Genome-Wide Profiling Reveals That Herbal Medicine Jinfukang-Induced Polyadenylation Alteration Is Involved in Anti-Lung Cancer Activity

Alternative polyadenylation (APA) plays an important role in regulation of genes expression and is involved in many biological processes. As eukaryotic cells receive a variety of external signals, genes produce diverse transcriptional isoforms and exhibit different translation efficiency. The traditional Chinese medicine (TCM) Jinfukang (JFK) has been effectively used for lung cancer treatment. In this study, we investigated whether JFK exerts its antitumor effect by modulating APA patterns in lung cancer cells. We performed a genome-wide APA site profiling analysis in JFK treated lung cancer cells A549 with 3T-seq approach that we reported previously. Comparing with those in untreated A549, in JFK treated A549 we observed APA-mediated 3′ UTRs alterations in 310 genes including 77 genes with shortened 3′ UTRs. In particular, we identified TMEM123, a gene involved in oncotic cell death, which produced transcripts with shortened 3′ UTR and thus was upregulated upon JFK treatment. Taken together, our studies suggest that APA might be one of the antitumor mechanisms of JFK and provide a new insight for the understanding of TCM against cancer.

Defining the 5΄ and 3΄ landscape of the Drosophila transcriptome with Exo-seq and RNaseH-seq

Cells regulate biological responses in part through changes in transcription start sites (TSS) or cleavage and polyadenylation sites (PAS). To fully understand gene regulatory networks, it is therefore critical to accurately annotate cell type-specific TSS and PAS. Here we present a simple and straightforward approach for genome-wide annotation of 5΄- and 3΄-RNA ends. Our approach reliably discerns bona fide PAS from false PAS that arise due to internal poly(A) tracts, a common problem with current PAS annotation methods. We applied our methodology to study the impact of temperature on the Drosophila melanogaster head transcriptome. We found hundreds of previously unidentified TSS and PAS which revealed two interesting phenomena: first, genes with multiple PASs tend to harbor a motif near the most proximal PAS, which likely represents a new cleavage and polyadenylation signal. Second, motif analysis of promoters of genes affected by temperature suggested that boundary element association factor of 32 kDa (BEAF-32) and DREF mediates a transcriptional program at warm temperatures, a result we validated in a fly line where beaf-32 is downregulated. These results demonstrate the utility of a high-throughput platform for complete experimental and computational analysis of mRNA-ends to improve gene annotation.