Nano3'RACE: A Method to Analyze Poly(A) Tail Length and Nucleotide Additions at the 3' Extremity of Selected mRNAs Using Nanopore Sequencing

Deadenylation is a major process that regulates gene expression by shaping the length of mRNA poly(A) tails. Deadenylation is controlled by factors in trans that recruit or impede deadenylases, by the incorporation of non-adenosines during poly(A) tail synthesis, and by the posttranscriptional addition of 3' nucleotides to poly(A) tails. Deciphering the regulation of poly(A) tail shortening requires both transcriptome-wide approaches and more targeted methodologies, allowing deep analyses of specific mRNAs. In this chapter, we present Nano3'RACE, a nanopore-based cDNA sequencing method that allows in-depth analysis to precisely measure poly(A) tail length and detect 3' terminal nucleotide addition, such as uridylation, for mRNAs of interest.

Absolute Quantification of mRNA Isoforms in Adult Stem Cells Using Microfluidic Digital PCR

Adult stem cells play key roles in homeostasis and tissue repair. These cells are regulated by a tight control of transcriptional programs. For example, muscle stem cells (MuSCs), located beneath the basal lamina, exist in the quiescent state but can transition to an activated, proliferative state upon injury. The control of MuSC state depends on the expression levels of myogenic transcription factors. Recent studies revealed the presence of different mRNA isoforms, with distinct biological regulation. Quantifying the exact expression levels of the mRNA isoforms encoding these myogenic transcription factors is therefore key to understanding how MuSCs switch between cell states. Previously, quantitative real-time polymerase chain reaction (qRT-PCR) has been used to quantify RNA expression levels. However, qRT-PCR depends on large amounts of RNA input and only measures relative abundance. Here, we present a protocol for the absolute quantification of mRNA isoforms using microfluidic digital PCR (mdPCR). Primary MuSCs isolated from individual skeletal muscles (gastrocnemius and masseter) are lysed, and their RNA is reverse-transcribed into cDNA and copied into double-stranded DNA. Following exonuclease I digestion to remove remaining single-stranded DNA, the samples are loaded onto a mdPCR chip with TaqMan probes targeting the mRNA isoforms of interest, whereupon target molecules are amplified in nanoliter chambers. We demonstrate that mdPCR can give exact molecule counts per cell for mRNA isoforms encoding the myogenic transcription factor Pax3. This protocol enables the absolute quantification of low abundant mRNA isoforms in a fast, precise, and reliable way.

Developmentally regulated alternate 3' end cleavage of nascent transcripts controls dynamic changes in protein expression in an adult stem cell lineage

Alternative polyadenylation (APA) generates transcript isoforms that differ in the position of the 3' cleavage site, resulting in the production of mRNA isoforms with different length 3' UTRs. Although widespread, the role of APA in the biology of cells, tissues, and organisms has been controversial. We identified >500 Drosophila genes that express mRNA isoforms with a long 3' UTR in proliferating spermatogonia but a short 3' UTR in differentiating spermatocytes due to APA. We show that the stage-specific choice of the 3' end cleavage site can be regulated by the arrangement of a canonical polyadenylation signal (PAS) near the distal cleavage site but a variant or no recognizable PAS near the proximal cleavage site. The emergence of transcripts with shorter 3' UTRs in differentiating cells correlated with changes in expression of the encoded proteins, either from off in spermatogonia to on in spermatocytes or vice versa. Polysome gradient fractionation revealed >250 genes where the long 3' UTR versus short 3' UTR mRNA isoforms migrated differently, consistent with dramatic stage-specific changes in translation state. Thus, the developmentally regulated choice of an alternative site at which to make the 3' end cut that terminates nascent transcripts can profoundly affect the suite of proteins expressed as cells advance through sequential steps in a differentiation lineage.

Nanopore sequencing reveals endogenous NMD-targeted isoforms in human cells

BACKGROUND

Nonsense-mediated mRNA decay (NMD) is a eukaryotic, translation-dependent degradation pathway that targets mRNAs with premature termination codons and also regulates the expression of some mRNAs that encode full-length proteins. Although many genes express NMD-sensitive transcripts, identifying them based on short-read sequencing data remains a challenge.

RESULTS

To identify and analyze endogenous targets of NMD, we apply cDNA Nanopore sequencing and short-read sequencing to human cells with varying expression levels of NMD factors. Our approach detects full-length NMD substrates that are highly unstable and increase in levels or even only appear when NMD is inhibited. Among the many new NMD-targeted isoforms that our analysis identifies, most derive from alternative exon usage. The isoform-aware analysis reveals many genes with significant changes in splicing but no significant changes in overall expression levels upon NMD knockdown. NMD-sensitive mRNAs have more exons in the 3΄UTR and, for those mRNAs with a termination codon in the last exon, the length of the 3΄UTR per se does not correlate with NMD sensitivity. Analysis of splicing signals reveals isoforms where NMD has been co-opted in the regulation of gene expression, though the main function of NMD seems to be ridding the transcriptome of isoforms resulting from spurious splicing events.

CONCLUSIONS

Long-read sequencing enables the identification of many novel NMD-sensitive mRNAs and reveals both known and unexpected features concerning their biogenesis and their biological role. Our data provide a highly valuable resource of human NMD transcript targets for future genomic and transcriptomic applications.

Direct long-read RNA sequencing identifies a subset of questionable exitrons likely arising from reverse transcription artifacts

Resistance to CD19-directed immunotherapies in lymphoblastic leukemia has been attributed, among other factors, to several aberrant CD19 pre-mRNA splicing events, including recently reported excision of a cryptic intron embedded within CD19 exon 2. While "exitrons" are known to exist in hundreds of human transcripts, we discovered, using reporter assays and direct long-read RNA sequencing (dRNA-seq), that the CD19 exitron is an artifact of reverse transcription. Extending our analysis to publicly available datasets, we identified dozens of questionable exitrons, dubbed "falsitrons," that appear only in cDNA-seq, but never in dRNA-seq. Our results highlight the importance of dRNA-seq for transcript isoform validation.

Deep sequencing of pre-translational mRNPs reveals hidden flux through evolutionarily conserved alternative splicing nonsense-mediated decay pathways

BACKGROUND

Alternative splicing, which generates multiple mRNA isoforms from single genes, is crucial for the regulation of eukaryotic gene expression. The flux through competing splicing pathways cannot be determined by traditional RNA-Seq, however, because different mRNA isoforms can have widely differing decay rates. Indeed, some mRNA isoforms with extremely short half-lives, such as those subject to translation-dependent nonsense-mediated decay (AS-NMD), may be completely overlooked in even the most extensive RNA-Seq analyses.

RESULTS

RNA immunoprecipitation in tandem (RIPiT) of exon junction complex components allows for purification of post-splicing mRNA-protein particles (mRNPs) not yet subject to translation (pre-translational mRNPs) and, therefore, translation-dependent mRNA decay. Here we compare exon junction complex RIPiT-Seq to whole cell RNA-Seq data from HEK293 cells. Consistent with expectation, the flux through known AS-NMD pathways is substantially higher than that captured by RNA-Seq. Our RIPiT-Seq also definitively demonstrates that the splicing machinery itself has no ability to detect reading frame. We identify thousands of previously unannotated splicing events; while many can be attributed to splicing noise, others are evolutionarily conserved events that produce new AS-NMD isoforms likely involved in maintenance of protein homeostasis. Several of these occur in genes whose overexpression has been linked to poor cancer prognosis.

CONCLUSIONS

Deep sequencing of RNAs in post-splicing, pre-translational mRNPs provides a means to identify and quantify splicing events without the confounding influence of differential mRNA decay. For many known AS-NMD targets, the nonsense-mediated decay-linked alternative splicing pathway predominates. Exon junction complex RIPiT-Seq also revealed numerous conserved but previously unannotated AS-NMD events.

Proteomic and Transcriptomic Analysis Identify Spliceosome as a Significant Component of the Molecular Machinery in the Pituitary Tumors Derived from POU1F1- and NR5A1-Cell Lineages

Background: Pituitary adenomas (PA) are the second most common tumor in the central nervous system and have low counts of mutated genes. Splicing occurs in 95% of the coding RNA. There is scarce information about the spliceosome and mRNA-isoforms in PA, and therefore we carried out proteomic and transcriptomic analysis to identify spliceosome components and mRNA isoforms in PA. Methods: Proteomic profile analysis was carried out by nano-HPLC and mass spectrometry with a quadrupole time-of-flight mass spectrometer. The mRNA isoforms and transcriptomic profiles were carried out by microarray technology. With proteins and mRNA information we carried out Gene Ontology and exon level analysis to identify splicing-related events. Results: Approximately 2000 proteins were identified in pituitary tumors. Spliceosome proteins such as SRSF1, U2AF1 and RBM42 among others were found in PA. These results were validated at mRNA level, which showed up-regulation of spliceosome genes in PA. Spliceosome-related genes segregate and categorize PA tumor subtypes. The PA showed alterations in CDK18 and THY1 mRNA isoforms which could be tumor specific. Conclusions: Spliceosome components are significant constituents of the PA molecular machinery and could be used as molecular markers and therapeutic targets. Splicing-related genes and mRNA-isoforms profiles characterize tumor subtypes.

3'READS + RIP defines differential Staufen1 binding to alternative 3'UTR isoforms and reveals structures and sequence motifs influencing binding and polysome association

Staufen1 (STAU1) is an RNA-binding protein (RBP) that interacts with double-stranded RNA structures and has been implicated in regulating different aspects of mRNA metabolism. Previous studies have indicated that STAU1 interacts extensively with RNA structures in coding regions (CDSs) and 3'-untranslated regions (3'UTRs). In particular, duplex structures formed within 3'UTRs by inverted-repeat Alu elements (IRAlus) interact with STAU1 through its double-stranded RNA-binding domains (dsRBDs). Using 3' region extraction and deep sequencing coupled to ribonucleoprotein immunoprecipitation (3'READS + RIP), together with reanalyzing previous STAU1 binding and RNA structure data, we delineate STAU1 interactions transcriptome-wide, including binding differences between alternative polyadenylation (APA) isoforms. Consistent with previous reports, RNA structures are dominant features for STAU1 binding to CDSs and 3'UTRs. Overall, relative to short 3'UTR counterparts, longer 3'UTR isoforms of genes have stronger STAU1 binding, most likely due to a higher frequency of RNA structures, including specific IRAlus sequences. Nevertheless, a sizable fraction of genes express transcripts showing the opposite trend, attributable to AU-rich sequences in their alternative 3'UTRs that may recruit antagonistic RBPs and/or destabilize RNA structures. Using STAU1-knockout cells, we show that strong STAU1 binding to mRNA 3'UTRs generally enhances polysome association. However, IRAlus generally have little impact on STAU1-mediated polysome association despite having strong interactions with the protein. Taken together, our work reveals complex interactions of STAU1 with its cognate RNA substrates. Our data also shed light on distinct post-transcriptional fates for the widespread APA isoforms in mammalian cells.

Computational Methods for Predicting Functions at the mRNA Isoform Level

Multiple mRNA isoforms of the same gene are produced via alternative splicing, a biological mechanism that regulates protein diversity while maintaining genome size. Alternatively spliced mRNA isoforms of the same gene may sometimes have very similar sequence, but they can have significantly diverse effects on cellular function and regulation. The products of alternative splicing have important and diverse functional roles, such as response to environmental stress, regulation of gene expression, human heritable, and plant diseases. The mRNA isoforms of the same gene can have dramatically different functions. Despite the functional importance of mRNA isoforms, very little has been done to annotate their functions. The recent years have however seen the development of several computational methods aimed at predicting mRNA isoform level biological functions. These methods use a wide array of proteo-genomic data to develop machine learning-based mRNA isoform function prediction tools. In this review, we discuss the computational methods developed for predicting the biological function at the individual mRNA isoform level.

Sierra: discovery of differential transcript usage from polyA-captured single-cell RNA-seq data

High-throughput single-cell RNA-seq (scRNA-seq) is a powerful tool for studying gene expression in single cells. Most current scRNA-seq bioinformatics tools focus on analysing overall expression levels, largely ignoring alternative mRNA isoform expression. We present a computational pipeline, Sierra, that readily detects differential transcript usage from data generated by commonly used polyA-captured scRNA-seq technology. We validate Sierra by comparing cardiac scRNA-seq cell types to bulk RNA-seq of matched populations, finding significant overlap in differential transcripts. Sierra detects differential transcript usage across human peripheral blood mononuclear cells and the Tabula Muris, and 3 'UTR shortening in cardiac fibroblasts. Sierra is available at https://github.com/VCCRI/Sierra .

Isoform-specific NF1 mRNA levels correlate with disease severity in Neurofibromatosis type 1

Background

Neurofibromatosis type 1 (NF1) is characterized by an extreme clinical variability both within and between families that cannot be explained solely by the nature of the pathogenic NF1 gene mutations. A proposed model hypothesizes that variation in the levels of protein isoforms generated via alternative transcript processing acts as modifier and contributes to phenotypic variability.

Results

Here we used real-time quantitative PCR to investigate the levels of two major NF1 mRNA isoforms encoding proteins differing in their ability to control RAS signaling (isoforms I and II) in the peripheral blood leukocytes of 138 clinically well-characterized NF1 patients and 138 aged-matched healthy controls. As expected, expression analysis showed that NF1 isoforms I and II levels were significantly lower in patients than controls. Notably, these differences were more evident when patients were stratified according to the severity of phenotype. Moreover, a correlation was identified when comparing the levels of isoform I mRNA and the severity of NF1 features, with statistically significant lower levels associated with a severe phenotype (i.e., occurrence of learning disability/intellectual disability, optic gliomas and/or other neoplasias, and/or cerebrovascular disease) as well as in patients with cognitive impairment.

Conclusions

The present findings provide preliminary evidence for a role of circuits controlling NF1 transcript processing in modulating NF1 expressivity, and document an association between the levels of neurofibromin isoform I mRNA and the severity of phenotype and cognitive impairment in NF1.

Probing In Vivo Structure of Individual mRNA 3' Isoforms Using Dimethyl Sulfate

The DMS region extraction and deep sequencing (DREADS) procedure was designed to probe RNA structure in vivo and to link this structural information to specific 3' isoforms. Growing cells are treated with the alkylating agent dimethyl sulfate (DMS), which enters easily into cells and modifies RNA molecules at solvent-exposed A and C residues. RNA is isolated, and sequencing libraries are constructed in a manner that preserves the identities of individual mRNA isoforms arising from alternative cleavage/polyadenylation sites. During the cDNA synthesis step of library construction, the progress of reverse transcriptase (RT) is blocked when it encounters a DMS modification on the RNA, leading to disproportionate cDNA termination adjacent to DMS-modified positions. After paired-end deep sequencing, the downstream end of each sequenced fragment is mapped to a specific cleavage/poly(A) site representing an individual mRNA 3' isoform. The upstream mapped end of the sequenced fragment defines where the RT reaction stopped. Over the population of all sequenced fragments derived from a particular isoform, A and C positions that are overrepresented next to the upstream endpoints in the DMS sample (relative to a parallel untreated control) are inferred to have been DMS modified, and hence solvent exposed. This method thus allows in vivo structural information obtained using DMS to be linked to individual mRNA 3' isoforms. © 2019 by John Wiley & Sons, Inc.

Alternative polyadenylation produces multiple 3' untranslated regions of odorant receptor mRNAs in mouse olfactory sensory neurons

BACKGROUND

Odorant receptor genes constitute the largest gene family in mammalian genomes and this family has been extensively studied in several species, but to date far less attention has been paid to the characterization of their mRNA 3' untranslated regions (3'UTRs). Given the increasing importance of UTRs in the understanding of RNA metabolism, and the growing interest in alternative polyadenylation especially in the nervous system, we aimed at identifying the alternative isoforms of odorant receptor mRNAs generated through 3'UTR variation.

RESULTS

We implemented a dedicated pipeline using IsoSCM instead of Cufflinks to analyze RNA-Seq data from whole olfactory mucosa of adult mice and obtained an extensive description of the 3'UTR isoforms of odorant receptor mRNAs. To validate our bioinformatics approach, we exhaustively analyzed the 3'UTR isoforms produced from 2 pilot genes, using molecular approaches including northern blot and RNA ligation mediated polyadenylation test. Comparison between datasets further validated the pipeline and confirmed the alternative polyadenylation patterns of odorant receptors. Qualitative and quantitative analyses of the annotated 3' regions demonstrate that 1) Odorant receptor 3'UTRs are longer than previously described in the literature; 2) More than 77% of odorant receptor mRNAs are subject to alternative polyadenylation, hence generating at least 2 detectable 3'UTR isoforms; 3) Splicing events in 3'UTRs are restricted to a limited subset of odorant receptor genes; and 4) Comparison between male and female data shows no sex-specific differences in odorant receptor 3'UTR isoforms.

CONCLUSIONS

We demonstrated for the first time that odorant receptor genes are extensively subject to alternative polyadenylation. This ground-breaking change to the landscape of 3'UTR isoforms of Olfr mRNAs opens new avenues for investigating their respective functions, especially during the differentiation of olfactory sensory neurons.

The Human TET2 Gene Contains Three Distinct Promoter Regions With Differing Tissue and Developmental Specificities

Tet methylcytosine dioxygenase 2 (TET2) is a tumor suppressor gene that is inactivated in a wide range of hematological cancers. TET2 enzymatic activity converts 5-methylcytosine (5-mC) into 5-hydroxymethylcytosine (5-hmC), an essential step in DNA demethylation. Human TET2 is highly expressed in pluripotent cells and down-regulated in differentiated cells: however, transcriptional regulation of the human TET2 gene has not been investigated in detail. Here we define three promoters within a 2.5 kb region located ∼ 87 kb upstream of the first TET2 coding exon. The three promoters, designated as Pro1, Pro2, and Pro3, generate three alternative first exons, and their presence in TET2 mRNAs varies with cell type and developmental stage. In general, all three TET2 transcripts are more highly expressed in human tissues rich in hematopoietic stem cells, such as spleen and bone marrow, compared to other tissues, such as brain and kidney. Transcripts from Pro2 are expressed by a broad range of tissues and at a significantly higher level than Pro1 or Pro3 transcripts. Pro3 transcripts were highly expressed by embryoid bodies generated from the H9 ES cell line, and the major Pro3 transcript is an alternatively spliced mRNA isoform that produces a truncated TET2 protein lacking the catalytic domain. Our study demonstrates distinct tissue-specific mechanisms of TET2 transcriptional regulation during early pluripotent states and in differentiated cell types.

Experiments with Snails Add to Our Knowledge about the Role of aPKC Subfamily Kinases in Learning

Protein kinase Mζ is considered important for memory formation and maintenance in different species, including invertebrates. PKMζ participates in multiple molecular pathways in neurons, regulating translation initiation rate, AMPA receptors turnover, synaptic scaffolding assembly, and other processes. Here, for the first time, we established the sequence of mRNA encoding PKMζ homolog in land snail Helix lucorum. We annotated important features of this mRNA: domains, putative capping sites, translation starts, and splicing sites. We discovered that this mRNA has at least two isoforms, and one of them lacks sequence encoding C1 domain. C1 deletion may be unique for snail because it has not been previously found in other species. We performed behavioral experiments with snails, measured expression levels of identified isoforms, and confirmed that their expression correlates with one type of learning.

Single-Molecule Real-Time (SMRT) Full-Length RNA-Sequencing Reveals Novel and Distinct mRNA Isoforms in Human Bone Marrow Cell Subpopulations

Hematopoietic cells are continuously replenished from progenitor cells that reside in the bone marrow. To evaluate molecular changes during this process, we analyzed the transcriptomes of freshly harvested human bone marrow progenitor (lineage-negative) and differentiated (lineage-positive) cells by single-molecule real-time (SMRT) full-length RNA-sequencing. This analysis revealed a ~5-fold higher number of transcript isoforms than previously detected and showed a distinct composition of individual transcript isoforms characteristic for bone marrow subpopulations. A detailed analysis of messenger RNA (mRNA) isoforms transcribed from the ANXA1 and EEF1A1 loci confirmed their distinct composition. The expression of proteins predicted from the transcriptome analysis was evaluated by mass spectrometry and validated previously unknown protein isoforms predicted e.g., for EEF1A1. These protein isoforms distinguished the lineage negative cell population from the lineage positive cell population. Finally, transcript isoforms expressed from paralogous gene loci (e.g., CFD, GATA2, HLA-A, B, and C) also distinguished cell subpopulations but were only detectable by full-length RNA sequencing. Thus, qualitatively distinct transcript isoforms from individual genomic loci separate bone marrow cell subpopulations indicating complex transcriptional regulation and protein isoform generation during hematopoiesis.

PD-1 in human NK cells: evidence of cytoplasmic mRNA and protein expression

Under physiological conditions, PD-1/PD-L1 interactions regulate unwanted over-reactions of immune cells and contribute to maintain peripheral tolerance. However, in tumor microenvironment, this interaction may greatly compromise the immune-mediated anti-tumor activity. PD-1⁺ NK cells have been detected in high percentage in peripheral blood and ascitic fluid of ovarian carcinoma patients. To acquire information on PD-1 expression and physiology in human NK cells, we analyzed whether PD-1 mRNA and protein are present in resting, surface PD-1⁻, NK cells from healthy donors. Both different splicing isoforms of PD-1 mRNA and a cytoplasmic pool of PD-1 protein were detected. Similar results were obtained also from both in vitro-activated and tumor-associated NK cells. PD-1 mRNA and protein were higher in CD56^dim than in CD56^bright NK cells. Confocal microscopy analyses revealed that PD-1 protein is present in virtually all NK cells analyzed. The present findings are compatible with a rapid surface expression of PD-1 in NK cells in response to appropriate, still undefined, stimuli.

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.

Investigation of Experimental Factors That Underlie BRCA1/2 mRNA Isoform Expression Variation: Recommendations for Utilizing Targeted RNA Sequencing to Evaluate Potential Spliceogenic Variants

PCR-based RNA splicing assays are commonly used in diagnostic and research settings to assess the potential effects of variants of uncertain clinical significance in BRCA1 and BRCA2. The Evidence-based Network for the Interpretation of Germline Mutant Alleles (ENIGMA) consortium completed a multicentre investigation to evaluate differences in assay design and the integrity of published data, raising a number of methodological questions associated with cell culture conditions and PCR-based protocols. We utilized targeted RNA-seq to re-assess BRCA1 and BRCA2 mRNA isoform expression patterns in lymphoblastoid cell lines (LCLs) previously used in the multicentre ENIGMA study. Capture of the targeted cDNA sequences was carried out using 34 BRCA1 and 28 BRCA2 oligonucleotides from the Illumina Truseq Targeted RNA Expression platform. Our results show that targeted RNA-seq analysis of LCLs overcomes many of the methodology limitations associated with PCR-based assays leading us to make the following observations and recommendations: (1) technical replicates (n > 2) of variant carriers to capture methodology induced variability associated with RNA-seq assays, (2) LCLs can undergo multiple freeze/thaw cycles and can be cultured up to 2 weeks without noticeably influencing isoform expression levels, (3) nonsense-mediated decay inhibitors are essential prior to splicing assays for comprehensive mRNA isoform detection, (4) quantitative assessment of exon:exon junction levels across BRCA1 and BRCA2 can help distinguish between normal and aberrant isoform expression patterns. Experimentally derived recommendations from this study will facilitate the application of targeted RNA-seq platforms for the quantitation of BRCA1 and BRCA2 mRNA aberrations associated with sequence variants of uncertain clinical significance.

Precursor microRNA-122 inhibits synthesis of Insig1 isoform mRNA by modulating polyadenylation site usage

The insulin-induced gene 1 protein (Insig1) inhibits the cholesterol biosynthesis pathway by retaining transcription factor SREBP in the endoplasmic reticulum, and by causing the degradation of HMGCR, the rate-limiting enzyme in cholesterol biosynthesis. Liver-specific microRNA miR-122, on the other hand, enhances cholesterol biosynthesis by an unknown mechanism. We have found that Insig1 mRNAs are generated by alternative cleavage and polyadenylation, resulting in specific isoform mRNA species. During high cholesterol abundance, the short 1.4-kb Insig1 mRNA was found to be preferentially translated to yield Insig1 protein. Precursor molecules of miR-122 down-regulated the translation of the 1.4-kb Insig1 isoform mRNA by interfering with the usage of the promoter-proximal cleavage-polyadenylation site that gives rise to the 1.4-kb Insig1 mRNA. These findings argue that precursor miR-122 molecules modulate polyadenylation site usage in Insig1 mRNAs, resulting in down-regulation of Insig1 protein abundance. Thus, precursor microRNAs may have hitherto undetected novel functions in nuclear gene expression.

From transcriptional complexity to cellular phenotypes: Lessons from yeast

Pervasive transcription has been reported in many eukaryotic organisms, revealing a highly interleaved transcriptome organization that involves thousands of coding and non-coding RNAs. However, to date, the biological impact of transcriptome complexity is still poorly understood. Here I will review how subtle variations of the transcriptome can lead to divergent cellular phenotypes by fine-tuning both its coding potential and regulation. I will discuss strategies that can be used to link molecular variations with divergent biological outcomes. Finally, I will explore the implication of transcriptional complexity for our understanding of gene expression in the context of cell-to-cell phenotypic variability. Copyright © 2017 John Wiley & Sons, Ltd.

The Evolutionary Relationship between Alternative Splicing and Gene Duplication

The protein diversity that exists today has resulted from various evolutionary processes. It is well known that gene duplication (GD) along with the accumulation of mutations are responsible, among other factors, for an increase in the number of different proteins. The gene structure in eukaryotes requires the removal of non-coding sequences, introns, to produce mature mRNAs. This process, known as cis-splicing, referred to here as splicing, is regulated by several factors which can lead to numerous splicing arrangements, commonly designated as alternative splicing (AS). AS, producing several transcripts isoforms form a single gene, also increases the protein diversity. However, the evolution and manner for increasing protein variation differs between AS and GD. An important question is how are patterns of AS affected after a GD event. Here, we review the current knowledge of AS and GD, focusing on their evolutionary relationship. These two processes are now considered the main contributors to the increasing protein diversity and therefore their relationship is a relevant, yet understudied, area of evolutionary study.

Complex tissue-specific patterns and distribution of multiple RAGE splice variants in different mammals

The receptor for advanced glycosylation end products (RAGE) is a multiligand receptor involved in diverse cell signaling pathways. Previous studies show that this gene expresses several splice variants in human, mouse, and dog. Alternative splicing (AS) plays an important role in expanding transcriptomic and proteomic diversity, and it has been related to disease. AS is also one of the main evolutionary mechanisms in mammalian genomes. However, limited information is available regarding the AS of RAGE in a wide context of mammalian tissues. In this study, we examined in detail the different RAGE mRNAs generated by AS from six mammals, including two primates (human and monkey), two artiodactyla (cow and pig), and two rodentia (mouse and rat) in 6–18 different tissues including fetal, adult, and tumor. By nested reverse transcription-polymerase chain reaction (RT-PCR) we identified a high number of splice variants including noncoding transcripts and predicted coding ones with different potential protein modifications affecting mainly the transmembrane and ligand-binding domains that could influence their biological function. However, analysis of RNA-seq data enabled detecting only the most abundant splice variants. More than 80% of the detected RT-PCR variants (87 of 101 transcripts) are novel (different exon/intron structure to the previously described ones), and interestingly, 20–60% of the total transcripts (depending on the species) are noncoding ones that present tissue specificity. Our results suggest that RAGE undergoes extensive AS in mammals, with different expression patterns among adult, fetal, and tumor tissues. Moreover, most splice variants seem to be species specific, especially the noncoding variants, with only two (canonical human Tv1-RAGE, and human N-truncated or Tv10-RAGE) conserved among the six different species. This could indicate a special evolution pattern of this gene at mRNA level.

Secondary structures involving the poly(A) tail and other 3' sequences are major determinants of mRNA isoform stability in yeast

In Saccharomyces cerevisiae, previous measurements of mRNA stabilities have been determined on a per-gene basis. We and others have recently shown that yeast genes give rise to a highly heterogeneous population of mRNAs due to extensive alternative 3’ end formation. Typical genes can have fifty or more distinct mRNA isoforms with 3’ endpoints differing by as little as one and as many as hundreds of nucleotides. In our recent paper [Geisberg et al. Cell (2014) 156: 812-824] we measured half-lives of individual mRNA isoforms in Saccharomyces cerevisiae by using the anchor away method for the rapid removal of Rpb1, the largest subunit of RNA Polymerase II, from the nucleus, followed by direct RNA sequencing of the cellular mRNA population over time. Combining these two methods allowed us to determine half-lives for more than 20,000 individual mRNA isoforms originating from nearly 5000 yeast genes. We discovered that different 3’ mRNA isoforms arising from the same gene can have widely different stabilities, and that such half-life variability across mRNA isoforms from a single gene is highly prevalent in yeast cells. Determining half-lives for many different mRNA isoforms from the same genes allowed us to identify hundreds of RNA sequence elements involved in the stabilization and destabilization of individual isoforms. In many cases, the poly(A) tail is likely to participate in the formation of stability-enhancing secondary structures at mRNA 3’ ends. Our results point to an important role for mRNA structure at 3’ termini in governing transcript stability, likely by reducing the interaction of the mRNA with the degradation apparatus.

Survivin isoforms and clinicopathological characteristics in colorectal adenocarcinomas using real-time qPCR

AIM: To investigate three isoforms of survivin in colorectal adenocarcinomas.

METHODS: We used the LightCycler Technology (Roche), along with a common forward primer and reverse primers specific for the splice variants and two common hybridization probes labeled with fluorescein and LightCycler-Red fluorophore (LC-Red 640). Real time quantitative polymerase chain reaction (PCR) was performed on cDNAs from 52 tumor specimens from colorectal cancer patients and 10 unrelated normal colorectal tissues. In the patients group, carcinoembryonic antigen (CEA) and CA19-9 tumor markers were also measured immunochemically.

RESULTS: Wild type survivin mRNA isoform was expressed in 48% of the 52 tumor samples, survivin-2b in 38% and survivin-ΔΕx3 in 29%, while no expression was found in normal tissues. The mRNA expression of wild type survivin presented a significant correlation with the expression of the ratio of survivin-2b, survivin-ΔΕx3, survivin-2b/wild type survivin and survivin-ΔΕx3/wild type survivin (P < 0.001). The mRNA expression of wild-survivin and survivin-ΔΕx3 was related with tumor size and invasion (P = 0.006 and P < 0.005, respectively). A significant difference was found between survivin-2b and morphologic cancer type. Also, the ratio of survivin-ΔEx3/wild-survivin was significantly associated with prognosis. No association was observed between the three isoforms and grade, metastasis, Dukes stage and gender. The three isoforms were not correlated with CEA and CA19-9.

CONCLUSION: Survivin isoforms may play a role in cell apoptosis and their quantification could provide information about clinical management of patients suffering from colorectal cancer.