Alternative splicing and bioinformatic analysis of human U12-type introns

U12-type introns exist, albeit rarely, in a variety of multicellular organisms. Splicing of U12 intron-containing precursor mRNAs takes place in the U12-type spliceosome that is distinct from the major U2-type spliceosome. Due to incompatibility of these two spliceosomes, alternative splicing involving a U12-type intron may give rise to a relatively complicated impact on gene expression. We studied alternative U12-type intron splicing in an attempt to gain more mechanistic insights. First, we characterized mutually exclusive exon selection of the human JNK2 gene, which involves an unusual intron possessing the U12-type 5′ splice site and the U2-type 3′ splice site. We demonstrated that the long and evolutionary conserved polypyrimidine tract of this hybrid intron provides important signals for inclusion of its downstream alternative exon. In addition, we examined the effects of single nucleotide polymorphisms in the human WDFY1 U12-type intron on pre-mRNA splicing. These results provide mechanistic implications on splice-site selection of U12-type intron splicing. We finally discuss the potential effects of splicing of a U12-type intron with genetic defects or within a set of genes encoding RNA processing factors on global gene expression.     

Estrogen-induced upregulation and 3′-UTR shortening of CDC6

3′-Untranslated region (UTR) shortening of mRNAs via alternative polyadenylation (APA) has important ramifications for gene expression. By using proximal APA sites and switching to shorter 3′-UTRs, proliferating cells avoid miRNA-mediated repression. Such APA and 3′-UTR shortening events may explain the basis of some of the proto-oncogene activation cases observed in cancer cells. In this study, we investigated whether 17 β-estradiol (E2), a potent proliferation signal, induces APA and 3′-UTR shortening to activate proto-oncogenes in estrogen receptor positive (ER+) breast cancers. Our initial probe based screen of independent expression arrays suggested upregulation and 3′-UTR shortening of an essential regulator of DNA replication, CDC6 (cell division cycle 6), upon E2 treatment. We further confirmed the E2- and ER-dependent upregulation and 3′UTR shortening of CDC6, which lead to increased CDC6 protein levels and higher BrdU incorporation. Consequently, miRNA binding predictions and dual luciferase assays suggested that 3′-UTR shortening of CDC6 was a mechanism to avoid 3′-UTR-dependent negative regulations. Hence, we demonstrated CDC6 APA induction by the proliferative effect of E2 in ER+ cells and provided new insights into the complex regulation of APA. E2-induced APA is likely to be an important but previously overlooked mechanism of E2-responsive gene expression.     

Single-cell Iso-Sequencing enables rapid genome annotation for scRNAseq analysis

Single-cell RNA sequencing is a powerful technique that continues to expand across various biological applications. However, incomplete 3′-UTR annotations can impede single-cell analysis resulting in genes that are partially or completely uncounted. Performing single-cell RNA sequencing with incomplete 3′-UTR annotations can hinder the identification of cell identities and gene expression patterns and lead to erroneous biological inferences. We demonstrate that performing single-cell isoform sequencing in tandem with single-cell RNA sequencing can rapidly improve 3′-UTR annotations. Using threespine stickleback fish (Gasterosteus aculeatus), we show that gene models resulting from a minimal embryonic single-cell isoform sequencing dataset retained 26.1% greater single-cell RNA sequencing reads than gene models from Ensembl alone. Furthermore, pooling our single-cell sequencing isoforms with a previously published adult bulk Iso-Seq dataset from stickleback, and merging the annotation with the Ensembl gene models, resulted in a marginal improvement (+0.8%) over the single-cell isoform sequencing only dataset. In addition, isoforms identified by single-cell isoform sequencing included thousands of new splicing variants. The improved gene models obtained using single-cell isoform sequencing led to successful identification of cell types and increased the reads identified of many genes in our single-cell RNA sequencing stickleback dataset. Our work illuminates single-cell isoform sequencing as a cost-effective and efficient mechanism to rapidly annotate genomes for single-cell RNA sequencing.     

Discovery and analysis of evolutionarily conserved intronic splicing regulatory elements

Knowledge of the functional cis-regulatory elements that regulate constitutive and alternative pre-mRNA splicing is fundamental for biology and medicine. Here we undertook a genome-wide comparative genomics approach using available mammalian genomes to identify conserved intronic splicing regulatory elements (ISREs). Our approach yielded 314 ISREs, and insertions of ~70 ISREs between competing splice sites demonstrated that 84% of ISREs altered 5′ and 94% altered 3′ splice site choice in human cells. Consistent with our experiments, comparisons of ISREs to known splicing regulatory elements revealed that 40%–45% of ISREs might have dual roles as exonic splicing silencers. Supporting a role for ISREs in alternative splicing, we found that 30%–50% of ISREs were enriched near alternatively spliced (AS) exons, and included almost all known binding sites of tissue-specific alternative splicing factors. Further, we observed that genes harboring ISRE-proximal exons have biases for tissue expression and molecular functions that are ISRE-specific. Finally, we discovered that for Nova1, neuronal PTB, hnRNP C, and FOX1, the most frequently occurring ISRE proximal to an alternative conserved exon in the splicing factor strongly resembled its own known RNA binding site, suggesting a novel application of ISRE density and the propensity for splicing factors to auto-regulate to associate RNA binding sites to splicing factors. Our results demonstrate that ISREs are crucial building blocks in understanding general and tissue-specific AS regulation and the biological pathways and functions regulated by these AS events.     

Genome-Wide Identification of Alternative Splice Forms Down-Regulated by Nonsense-Mediated mRNA Decay in Drosophila

Alternative mRNA splicing adds a layer of regulation to the expression of thousands of genes in Drosophila melanogaster. Not all alternative splicing results in functional protein; it can also yield mRNA isoforms with premature stop codons that are degraded by the nonsense-mediated mRNA decay (NMD) pathway. This coupling of alternative splicing and NMD provides a mechanism for gene regulation that is highly conserved in mammals. NMD is also active in Drosophila, but its effect on the repertoire of alternative splice forms has been unknown, as has the mechanism by which it recognizes targets. Here, we have employed a custom splicing-sensitive microarray to globally measure the effect of alternative mRNA processing and NMD on Drosophila gene expression. We have developed a new algorithm to infer the expression change of each mRNA isoform of a gene based on the microarray measurements. This method is of general utility for interpreting splicing-sensitive microarrays and high-throughput sequence data. Using this approach, we have identified a high-confidence set of 45 genes where NMD has a differential effect on distinct alternative isoforms, including numerous RNA–binding and ribosomal proteins. Coupled alternative splicing and NMD decrease expression of these genes, which may in turn have a downstream effect on expression of other genes. The NMD–affected genes are enriched for roles in translation and mitosis, perhaps underlying the previously observed role of NMD factors in cell cycle progression. Our results have general implications for understanding the NMD mechanism in fly. Most notably, we found that the NMD–target mRNAs had significantly longer 3′ untranslated regions (UTRs) than the nontarget isoforms of the same genes, supporting a role for 3′ UTR length in the recognition of NMD targets in fly.     

m6A‐mediated alternative splicing coupled with nonsense‐mediated mRNA decay regulates SAM synthetase homeostasis

Alternative splicing of pre‐mRNAs can regulate gene expression levels by coupling with nonsense‐mediated mRNA decay (NMD). In order to elucidate a repertoire of mRNAs regulated by alternative splicing coupled with NMD (AS‐NMD) in an organism, we performed long‐read RNA sequencing of poly(A)⁺ RNAs from an NMD‐deficient mutant strain of Caenorhabditis elegans, and obtained full‐length sequences for mRNA isoforms from 259 high‐confidence AS‐NMD genes. Among them are the S‐adenosyl‐L‐methionine (SAM) synthetase (sams) genes sams‐3 and sams‐4. SAM synthetase activity autoregulates sams gene expression through AS‐NMD in a negative feedback loop. We furthermore find that METT‐10, the orthologue of human U6 snRNA methyltransferase METTL16, is required for the splicing regulation in␣vivo, and specifically methylates the invariant AG dinucleotide at the distal 3′ splice site (3′SS) in␣vitro. Direct RNA sequencing coupled with machine learning confirms m⁶A modification of endogenous sams mRNAs. Overall, these results indicate that homeostasis of SAM synthetase in C. elegans is maintained by alternative splicing regulation through m⁶A modification at the 3′SS of the sams genes.     

CRISPRpas: programmable regulation of alternative polyadenylation by dCas9

Most human protein-coding genes produce alternative polyadenylation (APA) isoforms that differ in 3′ UTR size or, when coupled with splicing, have variable coding sequences. APA is an important layer of gene expression program critical for defining cell identity. Here, by using a catalytically dead Cas9 and coupling its target site with polyadenylation site (PAS), we develop a method, named CRISPRpas, to alter APA isoform abundance. CRISPRpas functions by enhancing proximal PAS usage, whose efficiency is influenced by several factors, including targeting strand of DNA, distance between PAS and target sequence and strength of the PAS. For intronic polyadenylation (IPA), splicing features, such as strengths of 5′ splice site and 3′ splice site, also affect CRISPRpas efficiency. We show modulation of APA of multiple endogenous genes, including IPA of PCF11, a master regulator of APA and gene expression. In sum, CRISPRpas offers a programmable tool for APA regulation that impacts gene expression.     

Genomic Hallmarks of Genes Involved in Chromosomal Translocations in Hematological Cancer

Reciprocal chromosomal translocations (RCTs) leading to the formation of fusion genes are important drivers of hematological cancers. Although the general requirements for breakage and fusion are fairly well understood, quantitative support for a general mechanism of RCT formation is still lacking. The aim of this paper is to analyze available high-throughput datasets with computational and robust statistical methods, in order to identify genomic hallmarks of translocation partner genes (TPGs). Our results show that fusion genes are generally overexpressed due to increased promoter activity of 5′ TPGs and to more stable 3′-UTR regions of 3′ TPGs. Furthermore, expression profiling of 5′ TPGs and of interaction partners of 3′ TPGs indicates that these features can help to explain tissue specificity of hematological translocations. Analysis of protein domains retained in fusion proteins shows that the co-occurrence of specific domain combinations is non-random and that distinct functional classes of fusion proteins tend to be associated with different components of the gene fusion network. This indicates that the configuration of fusion proteins plays an important role in determining which 5′ and 3′ TPGs will combine in specific fusion genes. It is generally accepted that chromosomal proximity in the nucleus can explain the specific pairing of 5′ and 3′ TPGS and the recurrence of hematological translocations. Using recently available data for chromosomal contact probabilities (Hi-C) we show that TPGs are preferentially located in early replicated regions and occupy distinct clusters in the nucleus. However, our data suggest that, in general, nuclear position of TPGs in hematological cancers explains neither TPG pairing nor clinical frequency. Taken together, our results support a model in which genomic features related to regulation of expression and replication timing determine the set of candidate genes more likely to be translocated in hematological tissues, with functional constraints being responsible for specific gene combinations.     

Identification of a Novel Transcript and Regulatory Mechanism for Microsomal Triglyceride Transfer Protein

Microsomal triglyceride transfer protein (MTP) is essential for the assembly of triglyceride-rich apolipoprotein B-containing lipoproteins. Previous studies in our laboratory identified a novel splice variant of MTP in mice that we named MTP-B. MTP-B has a unique first exon (1B) located 2.7 kB upstream of the first exon (1A) for canonical MTP (MTP-A). The two mature isoforms, though nearly identical in sequence and function, have different tissue expression patterns. In this study we report the identification of a second MTP splice variant (MTP-C), which contains both exons 1B and 1A. MTP-C is expressed in all the tissues we tested. In cells transfected with MTP-C, protein expression was less than 15% of that found when the cells were transfected with MTP-A or MTP-B. In silico analysis of the 5’-UTR of MTP-C revealed seven ATGs upstream of the start site for MTP-A, which is the only viable start site in frame with the main coding sequence. One of those ATGs was located in the 5’-UTR for MTP-A. We generated reporter constructs in which the 5’-UTRs of MTP-A or MTP-C were inserted between an SV40 promoter and the coding sequence of the luciferase gene and transfected these constructs into HEK 293 cells. Luciferase activity was significantly reduced by the MTP-C 5’-UTR, but not by the MTP-A 5’-UTR. We conclude that alternative splicing plays a key role in regulating MTP expression by introducing unique 5’-UTRs, which contain elements that alter translation efficiency, enabling the cell to optimize MTP levels and activity.