Characterization Of Alternative Splicing Research Articles

Characterization of Alternative Splicing in High-Risk Wilms' Tumors.

The significant heterogeneity of Wilms' tumors between different patients is thought to arise from genetic and epigenetic distortions that occur during various stages of fetal kidney development in a way that is poorly understood. To address this, we characterized the heterogeneity of alternative mRNA splicing in Wilms' tumors using a publicly available RNAseq dataset of high-risk Wilms' tumors and normal kidney samples. Through Pareto task inference and cell deconvolution, we found that the tumors and normal kidney samples are organized according to progressive stages of kidney development within a triangle-shaped region in latent space, whose vertices, or "archetypes", resemble the cap mesenchyme, the nephrogenic stroma, and epithelial tubular structures of the fetal kidney. We identified a set of genes that are alternatively spliced between tumors located in different regions of latent space and found that many of these genes are associated with the epithelial-to-mesenchymal transition (EMT) and muscle development. Using motif enrichment analysis, we identified putative splicing regulators, some of which are associated with kidney development. Our findings provide new insights into the etiology of Wilms' tumors and suggest that specific splicing mechanisms in early stages of development may contribute to tumor development in different patients.

IsoSplitter: identification and characterization of alternative splicing sites without a reference genome.

Long-read transcriptome sequencing is designed to sequence full-length RNA molecules and advantageous for identifying alternative splice isoforms; however, in the absence of a reference genome, it is difficult to accurately locate splice sites because of the diversity of patterns of alternative splicing (AS). Based on long-read transcriptome data, we developed a versatile tool, IsoSplitter, to reverse-trace and validate AS gene “split sites” with the following features: (i) IsoSplitter initially invokes a modified SIM4 program to find transcript split sites; (ii) each split site is then quantified, to reveal transcript diversity, and putative isoforms are grouped into gene clusters; (iii) an optional step for aligning short reads is provided, to validate split sites by identifying unique junction reads, and revealing and quantifying tissue-specific alternative splice isoforms. We tested IsoSplitter AS prediction using data sets from multiple model and nonmodel plant species and showed that the IsoSplitter pipeline is efficient to handle different transcriptomes with high accuracy. Furthermore, we evaluated the IsoSplitter pipeline compared with that of the splice junction identification tools, Program to Assemble Spliced Alignments (PASA software needs a reference genome for AS identification) and AStrap, using data from the model plant Arabidopsis thaliana. We found that IsoSplitter determined more than twice as many AS events than AStrap analysis, and 94.13% of the IsoSplitter predicted AS events were also identified by the PASA analysis. Starting from a simple sequence file, IsoSplitter is an assembly-free tool for identification and characterization of AS. IsoSplitter is developed and implemented in Python 3.5 using the Linux platform and is freely available at https://github.com/Hengfu-Yin/IsoSplitter.

Characterization of Alternative Splicing (AS) Events during Chicken (Gallus gallus) Male Germ-Line Stem Cell Differentiation with Single-Cell RNA-seq.

Simple SummaryStudies have shown that alternative splicing (AS) has been utilized in a wide variety of life processes. However, there are very few studies on AS during germ cell development. In this study, we preliminarily investigated the variation of variable shear events during the formation of chicken germ cells through the RNA-seq data analysis of embryonic stem cells (ESCs), gonad PGCs (gPGCs), and spermatogonia stem cells (SSCs), and the critical AS mode for several crucial stage-specific genes, which were identified during germ cell development. The results of this study lay a theoretical foundation for further analysis of the regulation mechanism of key genes involved in germ cell formation.Alternative splicing (AS) is a ubiquitous, co-transcriptional, and post-transcriptional regulation mechanism during certain developmental processes, such as germ cell differentiation. A thorough understanding of germ cell differentiation will help us to open new avenues for avian reproduction, stem cell biology, and advances in medicines for human consumption. Here, based on single-cell RNA-seq, we characterized genome-wide AS events in manifold chicken male germ cells: embryonic stem cells (ESCs), gonad primordial germ cells (gPGCs), and spermatogonia stem cells (SSCs). A total of 38,494 AS events from 15,338 genes were detected in ESCs, with a total of 48,955 events from 14,783 genes and 49,900 events from 15,089 genes observed in gPGCs and SSCs, respectively. Moreover, this distribution of AS events suggests the diverse splicing feature of ESCs, gPGCs, and SSCs. Finally, several crucial stage-specific genes, such as NANOG, POU5F3, LIN28B, BMP4, STRA8, and LHX9, were identified in AS events that were transmitted in ESCs, gPGCs, and SSCs. The gene expression results of the RNA-seq data were validated by qRT-PCR. In summary, we provided a comprehensive atlas of the genome-wide scale of the AS event landscape in male chicken germ-line cells and presented its distribution for the first time. This research may someday improve treatment options for men suffering from male infertility.

Single-Cell Alternative Splicing Analysis with Expedition Reveals Splicing Dynamics during Neuron Differentiation

Alternative splicing (AS) generates isoform diversity for cellular identity and homeostasis in multicellularlife. Although AS variation has been observed amongsingle cells, little is known about the biological or evolutionary significance of such variation. We developed Expedition, a computational framework consisting of outrigger, a de novo splice graph transversal algorithm to detect AS; anchor, a Bayesian approach to assign modalities; and bonvoyage, a visualization tool using non-negative matrix factorization to display modality changes. Applying Expedition to single pluripotent stem cells undergoing neuronal differentiation, we discover that up to 20% of AS exons exhibit bimodality. Bimodal exons are flanked by more conserved intronic sequences harboring distinct cis-regulatory motifs, constitute much of cell-type-specific splicing, are highly dynamic during cellular transitions, preserve reading frame, and reveal intricacy of cell states invisible to conventional gene expression analysis. Systematic AS characterization in single cells redefines our understanding of AS complexity in cell biology.

A Comprehensive Analysis of Alternative Splicing in Paleopolyploid Maize.

Identifying and characterizing alternative splicing (AS) enables our understanding of the biological role of transcript isoform diversity. This study describes the use of publicly available RNA-Seq data to identify and characterize the global diversity of AS isoforms in maize using the inbred lines B73 and Mo17, and a related species, sorghum. Identification and characterization of AS within maize tissues revealed that genes expressed in seed exhibit the largest differential AS relative to other tissues examined. Additionally, differences in AS between the two genotypes B73 and Mo17 are greatest within genes expressed in seed. We demonstrate that changes in the level of alternatively spliced transcripts (intron retention and exon skipping) do not solely reflect differences in total transcript abundance, and we present evidence that intron retention may act to fine-tune gene expression across seed development stages. Furthermore, we have identified temperature sensitive AS in maize and demonstrate that drought-induced changes in AS involve distinct sets of genes in reproductive and vegetative tissues. Examining our identified AS isoforms within B73 × Mo17 recombinant inbred lines (RILs) identified splicing QTL (sQTL). The 43.3% of cis-sQTL regulated junctions are actually identified as alternatively spliced junctions in our analysis, while 10 Mb windows on each side of 48.2% of trans-sQTLs overlap with splicing related genes. Using sorghum as an out-group enabled direct examination of loss or conservation of AS between homeologous genes representing the two subgenomes of maize. We identify several instances where AS isoforms that are conserved between one maize homeolog and its sorghum ortholog are absent from the second maize homeolog, suggesting that these AS isoforms may have been lost after the maize whole genome duplication event. This comprehensive analysis provides new insights into the complexity of AS in maize.

Identification and characterization of alternative splicing in parasitic nematode transcriptomes

BackgroundAlternative splicing (AS) of mRNA is a vital mechanism for enhancing genomic complexity in eukaryotes. Spliced isoforms of the same gene can have diverse molecular and biological functions and are often differentially expressed across various tissues, times, and conditions. Thus, AS has important implications in the study of parasitic nematodes with complex life cycles. Transcriptomic datasets are available from many species, but data must be revisited with splice-aware assembly protocols to facilitate the study of AS in helminthes.MethodsWe sequenced cDNA from the model worm Caenorhabditis elegans using 454/Roche technology for use as an experimental dataset. Reads were assembled with Newbler software, invoking the cDNA option. Several combinations of parameters were tested and assembled transcripts were verified by comparison with previously reported C. elegans genes and transcript isoforms and with Illumina RNAseq data.ResultsThoughtful adjustment of program parameters increased the percentage of assembled transcripts that matched known C. elegans sequences, decreased mis-assembly rates (i.e., cis- and trans-chimeras), and improved the coverage of the geneset. The optimized protocol was used to update de novo transcriptome assemblies from nine parasitic nematode species, including important pathogens of humans and domestic animals. Our assemblies indicated AS rates in the range of 20-30%, typically with 2-3 transcripts per AS locus, depending on the species. Transcript isoforms from the nine species were translated and searched for similarity to known proteins and functional domains. Some 21 InterPro domains, including several involved in nucleotide and chromatin binding, were statistically correlated with AS genetic loci. In most cases, the Roche/454 data explored in this study are the only sequences available from the species in question; however, the recently published genome of the human hookworm Necator americanus provided an additional opportunity to validate our results.ConclusionsOur optimized assembly parameters facilitated the first survey of AS among parasitic nematodes. The nine transcriptome assemblies, their protein translations, and basic annotations are available from Nematode.net as a resource for the research community. These should be useful for studies of specific genes and gene families of interest as well as for curating draft genome assemblies as they become available.

Snapshot liver transcriptome in hepatocellular carcinoma

Lately, advances in high throughput technologies in biomedical research have led to a dramatic increase in the accessibility of molecular insights at different levels of cancer biology such as genome, epigenome, transcriptome, proteome, and others. Among the diverse biological layers, the transcriptome has been most extensively studied especially due to the successful and broad introduction of the microarray technology. The future prospect of broad disposability of deep sequencing technology will furthermore lead to a more sensitive detection of lowly expressed transcripts and to an increase in the number of newly identified transcripts, but also to increase the discovery and characterization of alternative splicing. A common goal of large scale transcriptome profiling methods is the stratification of patients, eventually leading to personalized prognostic predictions and therapeutic strategies. Also, the observation that diverse etiologies of the mostly underlying liver cirrhosis (alcohol, HBV, HCV a.o.) are represented by different transcriptome profiles may be valuable for the identification of essential tumor targets [3] .I n

Gene Structure of the Murine Calcium Channel β3 Subunit, cDNA and Characterization of Alternative Splicing and Transcription Products

The beta3 subunit of high-voltage-gated calcium channels is a peripheral membrane protein that copurifies with neural N-type calcium channels. Murine genomic clones containing the full coding sequence of beta3 were isolated and the exons were mapped and sequenced. The murine calcium channel beta3 subunit is encoded by a unique gene composed of 13 translated exons that spread over approximately 8 kb of genomic sequence. Alternatively spliced transcripts of the beta3 gene were identified and characterized. The primary structure of beta3 is highly conserved between the murine, human, rabbit and rat proteins (98% identity). The intron placement within that primary structure correlates with the previously postulated exon positions in transcripts encoding the members of the calcium channel beta subunit family and confirm a close evolutionary relationship of the beta3, beta1, beta2 and beta4 subunit genes.

The human peripheral renzodiazepine receptor gene: Cloning and characterization of alternative splicing in normal tissues and in a patient with congenital lipoid adrenal hyperplasia

The mitochondrial benzodiazepine receptor (mBzR)appears to be a key factor in the flow of cholesterol into mitochondia to permit the initiation of steroid hormone synthesis. The mBzR consists of three components; the 18-kDa component on the outer mitochondrial membrane appears to contain the benzodiazepine binding site, and is hence often termed the peripheral benzodiazepine receptor (PBR). Using a cloned human PBR cDNA as probe, we have cloned the human PBR gene. The 13-kb gene is divided into four exons, with exon 1 encoding only a short 5′ untranslated segment. The 5′ flanking DNA lacks TATA and CAAT boxes but contains a cluster of SP-1 binding sites, typical of “housekeeping” genes. The encoded PBR mRNA is alternately spliced into two forms: “authentic” PBR mRNA retains all four exons, while a short form termed PBRS lacks exon 2. While PBR-S contains a 102-codon open reading frame with a typical initiator sequence, the reading frame differs from that of PBR, so that the encoded protein is unrelated to PBR. RT-PCR and RNase protection experiments confirm that both PBR and PBR-S are expressed in all tissues examined and that expression PBR-S is about 10 times the level of PBR. Expression of PBR cDNA in pCMV5 vectors transfected into COS-1 cells resulted in increased binding of [3H]PK11195, but expression of PBR-S did not. It has been speculated that patients with congenital lipoid adrenal hyperplasia, who cannot make any steroids, might have a genetic lesion in mBzR. RT-PCR analysis of testicular RNA from such a patient, sequencing of the cDNA, and blotting analysis of genomic DNA all indicate that the gene and mRNA for the PBR component of mBzR are normal in this disease.

Structure of the rat plasma membrane Ca(2+)-ATPase isoform 3 gene and characterization of alternative splicing and transcription products. Skeletal muscle-specific splicing results in a plasma membrane Ca(2+)-ATPase with a novel calmodulin-binding domain.

We have isolated the rat gene encoding isoform 3 of the plasma membrane Ca(2+)-ATPase (PMCA3) and have determined its exon/intron organization. The PMCA3 gene contains 24 exons and spans approximately 70 kilobases. In addition, we have analyzed the splicing and polyadenylation patterns leading to the production of an alternative 4.5-kilobase (PMCA3) skeletal muscle mRNA that differs from the previously characterized 7.5-kilobase brain mRNA (Greeb, J., and Shull, G. E. (1989) J. Biol. Chem. 264, 18569-18576). cDNA cloning, Northern blot hybridization, and polymerase chain reaction analyses of the 4.5-kilobase mRNA demonstrate (i) the inclusion of a novel 68-nucleotide exon (exon 22) that is specific for skeletal muscle and significantly alters the calmodulin-binding domain and (ii) the utilization of an alternative polyadenylation site following exon 23 which eliminates the last coding exon (exon 24) and 3'-untranslated sequence of the 7.5-kilobase mRNA. We have also identified a 42-nucleotide exon (exon 8) that is included in the skeletal muscle PMCA3 mRNAs, but may be either included or excluded in the brain mRNAs. Exon 8 is inserted immediately before the sequence encoding a putative phospholipid binding domain and thus may alter regulatory interactions of the enzyme with acidic phospholipids.