Molecular and genetic dissection of recursive splicing

Intronic ratchet points (RPs) are abundant within long introns in the Drosophila genome and consist of juxtaposed splice acceptor and splice donor (SD) sites. Although they appear to encompass zero-nucleotide exons, we recently clarified that intronic recursive splicing (RS) requires a cryptic exon at the RP (an RS-exon), which is subsequently always skipped and thus absent from mRNA. In addition, Drosophila encodes a smaller set of expressed exons bearing features of RS. Here, we investigate mechanisms that regulate the choice between RP and RS-exon SDs. First, analysis of Drosophila RP SD mutants demonstrates that SD competition suppresses inclusion of cryptic exons in endogenous contexts. Second, characterization of RS-exon reporters implicates exonic sequences as influencing choice of RS-exon usage. Using RS-exon swap and mutagenesis assays, we show exonic sequences can determine RS-exon inclusion. Finally, we provide evidence that splicing can suppress utilization of RP SDs to enable RS-exon expression. Overall, multiple factors can influence splicing of Drosophila RS-exons, which usually result in their complete suppression as zero-nucleotide RPs, but occasionally yield translated RS-exons.

Developmental dynamics of voltage-gated sodium channel isoform expression in the human and mouse brain

Background Genetic variants in the voltage-gated sodium channels SCN1A, SCN2A, SCN3A, and SCN8A are leading causes of epilepsy, developmental delay, and autism spectrum disorder. The mRNA splicing patterns of all four genes vary across development in the rodent brain, including mutually exclusive copies of the fifth protein-coding exon detected in the neonate (5N) and adult (5A). A second pair of mutually exclusive exons is reported in SCN8A only (18N and 18A). We aimed to quantify the expression of individual exons in the developing human brain. Methods RNA-seq data from 783 human brain samples across development were analyzed to estimate exon-level expression. Developmental changes in exon utilization were validated by assessing intron splicing. Exon expression was also estimated in RNA-seq data from 58 developing mouse neocortical samples. Results In the mature human neocortex, exon 5A is consistently expressed at least 4-fold higher than exon 5N in all four genes. For SCN2A, SCN3A, and SCN8A, a brain-wide synchronized 5N to 5A transition occurs between 24 post-conceptual weeks (2nd trimester) and 6 years of age. In mice, the equivalent 5N to 5A transition begins at or before embryonic day 15.5. In SCN8A, over 90% of transcripts in the mature human cortex include exon 18A. Early in fetal development, most transcripts include 18N or skip both 18N and 18A, with a transition to 18A inclusion occurring from 13 post-conceptual weeks to 6 months of age. No other protein-coding exons showed comparably dynamic developmental trajectories. Conclusions Exon usage in SCN1A, SCN2A, SCN3A, and SCN8A changes dramatically during human brain development. These splice isoforms, which alter the biophysical properties of the encoded channels, may account for some of the observed phenotypic differences across development and between specific variants. Manipulation of the proportion of splicing isoforms at appropriate stages of development may act as a therapeutic strategy for specific mutations or even epilepsy in general.

Developmental dynamics of voltage-gated sodium channel isoform expression in the human and mouse neocortex

ObjectiveGenetic variants in the voltage-gated sodium channels SCN1A, SCN2A, SCN3A, and SCN8A are leading causes of epilepsy, developmental delay, and autism spectrum disorder. The mRNA splicing patterns of all four genes vary across development in the rodent brain, including mutually exclusive copies of the fifth protein-coding exon detected in the neonate (5N) and adult (5A). A second pair of mutually exclusive exons is reported in SCN8A only (18N and 18A). We aimed to quantify the expression of individual exons in the developing human neocortex.MethodsRNA-seq data from 176 human dorsolateral prefrontal cortex samples across development were analyzed to estimate exon-level expression. Developmental changes in exon utilization were validated by assessing intron splicing. Exon expression was also estimated in RNA-seq data from 58 developing mouse neocortical samples.ResultsIn the mature human neocortex, exon 5A is consistently expressed at least 4-fold higher than exon 5N in all four genes. For SCN2A, SCN3A, and SCN8A a synchronized 5N/5A transition occurs between 24 post-conceptual weeks (2nd trimester) and six years of age. In mice, the equivalent 5N/5A transition begins at or before embryonic day 15.5. In SCN8A, over 90% of transcripts in the mature human cortex include exon 18A. Early in fetal development, most transcripts include 18N or skip both 18N and 18A, with a transition to 18A inclusion occurring from 13 post-conceptual weeks to 6 months of age. No other protein-coding exons showed comparably dynamic developmental trajectories.SignificanceSplice isoforms, which alter the biophysical properties of the encoded channels, may account for some of the observed phenotypic differences across development and between specific variants. Manipulation of the proportion of splicing isoforms at appropriate stages of development may act as a therapeutic strategy for specific mutations or even epilepsy in general.

RNA-seq analysis identifies cytoskeletal structural genes and pathways for meat quality in beef

RNA sequencing (RNA-seq) has allowed for transcriptional profiling of biological systems through the identification of differentially expressed (DE) genes and pathways. A total of 80 steers with extreme phenotypes were selected from the University of Florida multibreed Angus-Brahman herd. The average slaughter age was 12.91±8.69 months. Tenderness, juiciness and connective tissue assessed by sensory panel, along with marbling, Warner-Bratzler Shear Force (WBSF) and cooking loss, were measured in longissimus dorsi muscle. Total RNA was extracted from muscle and one RNA-seq library per sample was constructed, multiplexed, and sequenced based on protocols by Illumina HiSeq-3000 platform to generate 2×101 bp paired-end reads. The overall read mapping rate using the Btau_4.6.1 reference genome was 63%. A total of 8,799 genes were analyzed using two different methodologies, an expression association and a DE analysis. A gene and exon expression association analysis was carried out using a meat quality index on all 80 samples as a continuous response variable. The expression of 208 genes and 3,280 exons from 1,565 genes was associated with the meat quality index (p-value ≤ 0.05). A gene and isoform DE evaluation was performed analyzing two groups with extreme WBSF, tenderness and marbling. A total of 676 (adjusted p-value≤0.05), 70 (adjusted p-value≤0.1) and 198 (adjusted p-value≤0.1) genes were DE for WBSF, tenderness and marbling, respectively. A total of 106 isoforms from 98 genes for WBSF, 13 isoforms from 13 genes for tenderness and 43 isoforms from 42 genes for marbling (FDR≤0.1) were DE. Cytoskeletal and transmembrane anchoring genes and pathways were identified in the expression association, DE and the gene enrichment analyses; these proteins can have a direct effect on meat quality. Cytoskeletal proteins and transmembrane anchoring molecules can influence meat quality by allowing cytoskeletal interaction with myocyte and organelle membranes, contributing to cytoskeletal structure and architecture maintenance postmortem.

DIGGER: exploring the functional role of alternative splicing in protein interactions

Alternative splicing plays a major role in regulating the functional repertoire of the proteome. However, isoform-specific effects to protein-protein interactions (PPIs) are usually overlooked, making it impossible to judge the functional role of individual exons on a systems biology level. We overcome this barrier by integrating protein-protein interactions, domain-domain interactions and residue-level interactions information to lift exon expression analysis to a network level. Our user-friendly database DIGGER is available at https://exbio.wzw.tum.de/digger and allows users to seamlessly switch between isoform and exon-centric views of the interactome and to extract sub-networks of relevant isoforms, making it an essential resource for studying mechanistic consequences of alternative splicing.

Coordinated analysis of exon and intron data reveals novel differential gene expression changes

RNA-Seq expression analysis currently relies primarily upon exon expression data. The recognized role of introns during translation, and the presence of substantial RNA-Seq counts attributable to introns, provide the rationale for the simultaneous consideration of both exon and intron data. We describe here a method for the coordinated analysis of exon and intron data by investigating their relationship within individual genes and across samples, while taking into account changes in both variability and expression level. This coordinated analysis of exon and intron data offers strong evidence for significant differences that distinguish the profiles of the exon-only expression data from the combined exon and intron data. One advantage of our proposed method, called matched change characterization for exons and introns (MEI), is its straightforward applicability to existing archived data using small modifications to standard RNA-Seq pipelines. Using MEI, we demonstrate that when data are examined for changes in variability across control and case conditions, novel differential changes can be detected. Notably, when MEI criteria were employed in the analysis of an archived data set involving polyarthritic subjects, the number of differentially expressed genes was expanded by sevenfold. More importantly, the observed changes in exon and intron variability with statistically significant false discovery rates could be traced to specific immune pathway gene networks. The application of MEI analysis provides a strategy for incorporating the significance of exon and intron variability and further developing the role of using both exons and intron sequencing counts in studies of gene regulatory processes.

P120 Catenin Isoforms Differentially Associate with Breast Cancer Invasion and Metastasis

Tumor metastasis is the endpoint of tumor progression and depends on the ability of tumor cells to locally invade tissue, transit through the bloodstream and ultimately to colonize secondary organs at distant sites. P120 catenin (p120) has been implicated as an important regulator of metastatic dissemination because of its roles in cell–cell junctional stability, cytoskeletal dynamics, growth and survival. However, conflicting roles for p120 in different tumor models and steps of metastasis have been reported, and the understanding of p120 functions is confounded by the differential expression of p120 isoforms, which differ in N-terminal length, tissue localization and, likely, function. Here, we used in silico exon expression analyses, in vitro invasion assays and both RT-PCR and immunofluorescence of human tumors. We show that alternative exon usage favors expression of short isoform p120-3 in 1098 breast tumors and correlates with poor prognosis. P120-3 is upregulated at the invasive front of breast cancer cells migrating as collective groups in vitro. Furthermore, we demonstrate in histological sections of 54 human breast cancer patients that p120-3 expression is maintained throughout the metastatic cascade, whereas p120-1 is differentially expressed and diminished during invasion and in metastases. These data suggest specific regulation and functions of p120-3 in breast cancer invasion and metastasis.

RNA-seq analysis identifies cytoskeletal structural genes and pathways for meat quality in beef

Background RNA sequencing (RNA-seq) has allowed for transcriptional profiling of biological systems through identification of differentially expressed (DE) genes and pathways. Results A total of 80 steers were selected from the multibreed Angus-Brahman herd of the University of Florida. Sensory panel tenderness, juiciness and connective tissue as well as marbling, WBSF and cooking loss were assessed in longissimus dorsi muscle. Nuclear RNA was extracted from muscle and an RNA-seq library for each sample was constructed, multiplexed, and sequenced based on protocols by Illumina HiSeq 3000 PE100 platform to generate 2 × 101 bp paired-end reads. On average, 34.9 million high-quality paired reads were uniquely mapped to the Btau_4.6.1 reference genome and a total of 8,799 genes were analyzed. Including all 80 animals, gene and exon expression analysis was carried out using a meat quality index as a continuous response variable. The expression of 208 genes and 3,280 exons from 1,565 genes was associated with the meat quality index (p-value ≤ 0.05). Out of the 80 samples sequenced, 40 animals with extreme low and high WBSF, tenderness and marbling values were selected for a differential expression (DE) analysis for gene and isoforms. A total of 676 (adjusted p-value ≤ 0.05), 70 (adjusted p-value ≤ 0.1) and 198 (adjusted p-value ≤ 0.1) genes were DE for WBSF, tenderness and marbling, respectively. A total of 106 isoforms from 98 genes for WBSF, 13 isoforms from 13 genes for tenderness and 43 isoforms from 42 genes for marbling (FDR ≤ 0.1) were DE. Conclusion A number of cytoskeletal and transmembrane anchoring related genes and pathways were identified in the expression, DE and gene enrichment analyses, and these proteins can have a direct effect on meat quality. Cytoskeletal proteins and transmembrane anchoring molecules can influence meat quality by allowing cytoskeletal filament interaction with myocyte and organelle membranes, contributing to cytoskeletal structure, microtubule network stability, and cellular architecture maintenance during the postmortem.