Detection and visualization of differential splicing in RNA-Seq data with JunctionSeq

Local and systemic transcriptome and spliceome reprogramming induced by the root-knot nematode Meloidogyne incognita in tomato

Root-knot nematodes (Meloidogyne spp.) are widely spread root parasites that infect thousands of vascular plant species. These highly polyphagous nematodes engage in sophisticated interactions with host plants that results in the formation of knot-like structures known as galls whose ontogeny remains largely unknown. Here, we determined transcriptome changes and alternative splicing variants induced by Megalaima incognita in galls and neighboring root cells at two distinct infective stages. M. incognita induced substantial transcriptome changes in tomato roots both locally in galls and systemically in neighboring cells. A considerable parallel regulation of gene expression in galls and neighboring cells were detected, indicative of effective intercellular communications exemplified by suppression of basal defense responses particularly during the early stage of infection. The transcriptome analysis also revealed that M. incognita exerts a tight control over the cell cycle process as a whole that results in an increase of ploidy levels in the feeding sites and accelerated mitotic activity of the gall cells. Alternative splicing analysis indicated that M. incognita significantly modulates pre-mRNA splicing as a total of 9064 differentially spliced events from 2898 genes were identified where intron retention and exon skipping events were largely suppressed. Furthermore, a number of differentially spliced events were functionally validated using transgenic hairy root system and found to impact gall formation and nematode egg mass production. Together, our data provide unprecedented insights into the transcriptome and spliceome reprogramming induced by M. incognita in tomato with respect to gall ontogeny and nematode parasitism.

Canine hyper-sociability structural variants associated with altered three-dimensional chromatin state

Strong selection on complex traits can lead to skewed trait means and reduced trait variability in populations. An example of this phenomenon can be evidenced in allele frequency changes and skewed trait distributions driven by persistent human-directed selective pressures in domesticated species. Dog domestication is linked to several genomic variants; however, the functional impacts of these variants may not always be straightforward when found in non-coding regions of the genome. Four polymorphic transposable elements (TE) found within non-coding sites along a 5 Mb region on canine CFA6 have evolved due to directional selection associated with heightened human-directed hyper-sociability in domesticated dogs. We found that the polymorphic TE in intron 17 of the canine GTF2I gene, which was previously reported to be negatively correlated with canid human-directed hyper-sociability, is associated with altered chromatin looping and hence distinct cis-regulatory landscapes. We reported supporting evidence of an E2F1-DNA binding peak concordant with the altered loop and higher expression of GTF2I exon 18, indicative of alternative splicing. Globally, we discovered differences in pathways regulating the extra-cellular matrix with respect to TE copy number. Overall, we reported evidence suggesting an intriguing molecular convergence between the emergence of hypersocial behaviors in dogs and the same genes that, when hemizygous, produce human Williams Beuren Syndrome characterized by cranio-facial defects and heightened social behaviors. Our results additionally emphasize the often-overlooked potential role of chromatin architecture in social evolution.

betAS: intuitive analysis and visualization of differential alternative splicing using beta distributions

Next-generation RNA sequencing allows alternative splicing (AS) quantification with unprecedented resolution, with the relative inclusion of an alternative sequence in transcripts being commonly quantified by the proportion of reads supporting it as percent spliced-in (PSI). However, PSI values do not incorporate information about precision, proportional to the respective AS events' read coverage. Beta distributions are suitable to quantify inclusion levels of alternative sequences, using reads supporting their inclusion and exclusion as surrogates for the two distribution shape parameters. Each such beta distribution has the PSI as its mean value and is narrower when the read coverage is higher, facilitating the interpretability of its precision when plotted. We herein introduce a computational pipeline, based on beta distributions accurately modeling PSI values and their precision, to quantitatively and visually compare AS between groups of samples. Our methodology includes a differential splicing significance metric that compromises the magnitude of intergroup differences, the estimation uncertainty in individual samples, and the intragroup variability, being therefore suitable for multiple-group comparisons. To make our approach accessible and clear to both noncomputational and computational biologists, we developed betAS, an interactive web app and user-friendly R package for visual and intuitive differential splicing analysis from read count data.

Comprehensive benchmark of differential transcript usage analysis for static and dynamic conditions

RNA sequencing offers unique insights into transcriptome diversity, and a plethora of tools have been developed to analyze alternative splicing. One important task is to detect changes in the relative transcript abundance in differential transcript usage (DTU) analysis. The choice of the right analysis tool is non-trivial and depends on experimental factors such as the availability of single- or paired-end and bulk or single-cell data. To help users select the most promising tool for their task, we performed a comprehensive benchmark of DTU detection tools. We cover a wide array of experimental settings, using simulated bulk and single-cell RNA-seq data as well as real transcriptomics datasets, including time-series data. Our results suggest that DEXSeq, edgeR, and LimmaDS are better choices for paired-end data, while DSGseq and DEXSeq can be used for single-end data. In single-cell simulation settings, we showed that satuRn performs better than DTUrtle. In addition, we showed that Spycone is optimal for time series DTU/IS analysis based on the evidence provided using GO terms enrichment analysis.

Spliceosome malfunction causes neurodevelopmental disorders with overlapping features

Pre-mRNA splicing is a highly coordinated process. While its dysregulation has been linked to neurological deficits, our understanding of the underlying molecular and cellular mechanisms remains limited. We implicated pathogenic variants in U2AF2 and PRPF19, encoding spliceosome subunits in neurodevelopmental disorders (NDDs), by identifying 46 unrelated individuals with 23 de novo U2AF2 missense variants (including 7 recurrent variants in 30 individuals) and 6 individuals with de novo PRPF19 variants. Eight U2AF2 variants dysregulated splicing of a model substrate. Neuritogenesis was reduced in human neurons differentiated from human pluripotent stem cells carrying two U2AF2 hyper-recurrent variants. Neural loss of function (LoF) of the Drosophila orthologs U2af50 and Prp19 led to lethality, abnormal mushroom body (MB) patterning, and social deficits, which were differentially rescued by wild-type and mutant U2AF2 or PRPF19. Transcriptome profiling revealed splicing substrates or effectors (including Rbfox1, a third splicing factor), which rescued MB defects in U2af50-deficient flies. Upon reanalysis of negative clinical exomes followed by data sharing, we further identified 6 patients with NDD who carried RBFOX1 missense variants which, by in vitro testing, showed LoF. Our study implicates 3 splicing factors as NDD-causative genes and establishes a genetic network with hierarchy underlying human brain development and function.

Comprehensive transcriptomic analysis shows disturbed calcium homeostasis and deregulation of T lymphocyte apoptosis in inclusion body myositis

OBJECTIVE

Inclusion body myositis (IBM) has an unclear molecular etiology exhibiting both characteristic inflammatory T-cell activity and rimmed-vacuolar degeneration of muscle fibers. Using in-depth gene expression and splicing studies, we aimed at understanding the different components of the molecular pathomechanisms in IBM.

METHODS

We performed RNA-seq on RNA extracted from skeletal muscle biopsies of clinically and histopathologically defined IBM (n = 24), tibial muscular dystrophy (n = 6), and histopathologically normal group (n = 9). In a comprehensive transcriptomics analysis, we analyzed the differential gene expression, differential splicing and exon usage, downstream pathway analysis, and the interplay between coding and non-coding RNAs (micro RNAs and long non-coding RNAs).

RESULTS

We observe dysregulation of genes involved in calcium homeostasis, particularly affecting the T-cell activity and regulation, causing disturbed Ca2+-induced apoptotic pathways of T cells in IBM muscles. Additionally, LCK/p56, which is an essential gene in regulating the fate of T-cell apoptosis, shows increased expression and altered splicing usage in IBM muscles.

INTERPRETATION

Our analysis provides a novel understanding of the molecular mechanisms in IBM by showing a detailed dysregulation of genes involved in calcium homeostasis and its effect on T-cell functioning in IBM muscles. Loss of T-cell regulation is hypothesized to be involved in the consistent observation of no response to immune therapies in IBM patients. Our results show that loss of apoptotic control of cytotoxic T cells could indeed be one component of their abnormal cytolytic activity in IBM muscles.

Transcriptional Activation, Deactivation and Rebound Patterns in Cortex, Hippocampus and Amygdala in Response to Ketamine Infusion in Rats

Ketamine, an N-methyl-D-aspartate (NMDA)-receptor antagonist, is a recently revitalized treatment for pain and depression, yet its actions at the molecular level remain incompletely defined. In this molecular-pharmacological investigation in the rat, we used short- and longer-term infusions of high dose ketamine to stimulate neuronal transcription processes. We hypothesized that a progressively stronger modulation of neuronal gene networks would occur over time in cortical and limbic pathways. A continuous intravenous administration paradigm for ketamine was developed in rat consisting of short (1 h) and long duration (10 h, and 10 h + 24 h recovery) infusions of anesthetic concentrations to activate or inhibit gene transcription in a pharmacokinetically controlled fashion. Transcription was measured by RNA-Seq in three brain regions: frontal cortex, hippocampus, and amygdala. Cellular level gene localization was performed with multiplex fluorescent in situ hybridization. Induction of a shared transcriptional regulatory network occurred within 1 h in all three brain regions consisting of (a) genes involved in stimulus-transcription factor coupling that are induced during altered synaptic activity (immediate early genes, IEGs, such as c-Fos, 9–12 significant genes per brain region, p < 0.01 per gene) and (b) the Nrf2 oxidative stress-antioxidant response pathway downstream from glutamate signaling (Nuclear Factor Erythroid-Derived 2-Like 2) containing 12–25 increasing genes (p < 0.01) per brain region. By 10 h of infusion, the acute results were further reinforced and consisted of more and stronger gene alterations reflecting a sustained and accentuated ketamine modulation of regional excitation and plasticity. At the cellular level, in situ hybridization localized up-regulation of the plasticity-associated gene Bdnf, and the transcription factors Nr4a1 and Fos, in cortical layers III and V. After 24 h recovery, we observed overshoot of transcriptional processes rather than a smooth return to homeostasis suggesting an oscillation of plasticity occurs during the transition to a new phase of neuronal regulation. These data elucidate critical molecular regulatory actions during and downstream of ketamine administration that may contribute to the unique drug actions of this anesthetic agent. These molecular investigations point to pathways linked to therapeutically useful attributes of ketamine.

Probing Isoform Switching Events in Various Cancer Types: Lessons From Pan-Cancer Studies

Alternative splicing is an essential regulatory mechanism for gene expression in mammalian cells contributing to protein, cellular, and species diversity. In cancer, alternative splicing is frequently disturbed, leading to changes in the expression of alternatively spliced protein isoforms. Advances in sequencing technologies and analysis methods led to new insights into the extent and functional impact of disturbed alternative splicing events. In this review, we give a brief overview of the molecular mechanisms driving alternative splicing, highlight the function of alternative splicing in healthy tissues and describe how alternative splicing is disrupted in cancer. We summarize current available computational tools for analyzing differential transcript usage, isoform switching events, and the pathogenic impact of cancer-specific splicing events. Finally, the strategies of three recent pan-cancer studies on isoform switching events are compared. Their methodological similarities and discrepancies are highlighted and lessons learned from the comparison are listed. We hope that our assessment will lead to new and more robust methods for cancer-specific transcript detection and help to produce more accurate functional impact predictions of isoform switching events.

Multiparametric Evaluation of Post-MI Small Animal Models Using Metabolic ([18F]FDG) and Perfusion-Based (SYN1) Heart Viability Tracers

Cardiovascular diseases (CVD), with myocardial infarction (MI) being one of the crucial components, wreak havoc in developed countries. Advanced imaging technologies are required to obtain quick and widely available diagnostic data. This paper describes a multimodal approach to in vivo perfusion imaging using the novel SYN1 tracer based on the fluorine-18 isotope. The NOD-SCID mice were injected intravenously with SYN1 or [¹⁸F] fluorodeoxyglucose ([¹⁸F]-FDG) radiotracers after induction of the MI. In all studies, the positron emission tomography–computed tomography (PET/CT) technique was used. To obtain hemodynamic data, mice were subjected to magnetic resonance imaging (MRI). Finally, the biodistribution of the SYN1 compound was performed using Wistar rat model. SYN1 showed normal accumulation in mouse and rat hearts, and MI hearts correctly indicated impaired cardiac segments when compared to [¹⁸F]-FDG uptake. In vivo PET/CT and MRI studies showed statistical convergence in terms of the size of the necrotic zone and cardiac function. This was further supported with RNAseq molecular analyses to correlate the candidate function genes’ expression, with Serpinb1c, Tnc and Nupr1, with Trem2 and Aldolase B functional correlations showing statistical significance in both SYN1 and [¹⁸F]-FDG. Our manuscript presents a new fluorine-18-based perfusion radiotracer for PET/CT imaging that may have importance in clinical applications. Future research should focus on confirmation of the data elucidated here to prepare SYN1 for first-in-human trials.

A prenylated dsRNA sensor protects against severe COVID-19

Inherited genetic factors can influence the severity of COVID-19, but the molecular explanation underpinning a genetic association is often unclear. Intracellular antiviral defenses can inhibit the replication of viruses and reduce disease severity. To better understand the antiviral defenses relevant to COVID-19, we used interferon-stimulated gene (ISG) expression screening to reveal that 2′-5′-oligoadenylate synthetase 1 (OAS1), through ribonuclease L, potently inhibits severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). We show that a common splice-acceptor single-nucleotide polymorphism (Rs10774671) governs whether patients express prenylated OAS1 isoforms that are membrane-associated and sense-specific regions of SARS-CoV-2 RNAs or if they only express cytosolic, nonprenylated OAS1 that does not efficiently detect SARS-CoV-2. In hospitalized patients, expression of prenylated OAS1 was associated with protection from severe COVID-19, suggesting that this antiviral defense is a major component of a protective antiviral response.

Alternative RNA splicing in stem cells and cancer stem cells: Importance of transcript-based expression analysis

Alternative ribonucleic acid (RNA) splicing can lead to the assembly of different protein isoforms with distinctive functions. The outcome of alternative splicing (AS) can result in a complete loss of function or the acquisition of new functions. There is a gap in knowledge of abnormal RNA splice variants promoting cancer stem cells (CSCs), and their prospective contribution in cancer progression. AS directly regulates the self-renewal features of stem cells (SCs) and stem-like cancer cells. Notably, octamer-binding transcription factor 4A spliced variant of octamer-binding transcription factor 4 contributes to maintaining stemness properties in both SCs and CSCs. The epithelial to mesenchymal transition pathway regulates the AS events in CSCs to maintain stemness. The alternative spliced variants of CSCs markers, including cluster of differentiation 44, aldehyde dehydrogenase, and doublecortin-like kinase, α6β1 integrin, have pivotal roles in increasing self-renewal properties and maintaining the pluripotency of CSCs. Various splicing analysis tools are considered in this study. LeafCutter software can be considered as the best tool for differential splicing analysis and identification of the type of splicing events. Additionally, LeafCutter can be used for efficient mapping splicing quantitative trait loci. Altogether, the accumulating evidence re-enforces the fact that gene and protein expression need to be investigated in parallel with alternative splice variants.

Epigenetic differences in an identical genetic background modulate alternative splicing in A. thaliana

How stable and temperature-dependent variations in DNA methylation and nucleosome occupancy influence alternative splicing (AS) remains poorly understood in plants. To answer this, we generated transcriptome, whole-genome bisulfite, and MNase sequencing data for an epigenetic Recombinant Inbred Line (epiRIL) of A. thaliana at normal and cold temperature. For comparative analysis, the same data sets for the parental ecotype Columbia (Col-0) were also generated, whereas for DNA methylation, previously published high confidence methylation profiles of Col-0 were used. Significant epigenetic differences in an identical genetic background were observed between Col-0 and epiRIL lines under normal and cold temperatures. Our transcriptome data revealed that differential DNA methylation and nucleosome occupancy modulate expression levels of many genes and AS in response to cold. Collectively, DNA methylation and nucleosome levels exhibit characteristic patterns around intron-exon boundaries at normal and cold conditions, and any perturbation in them, in an identical genetic background is sufficient to modulate AS in Arabidopsis.

A Global Landscape of Miniature Inverted-Repeat Transposable Elements in the Carrot Genome

Miniature inverted-repeat transposable elements (MITEs) are the most abundant group of Class II mobile elements in plant genomes. Their presence in genic regions may alter gene structure and expression, providing a new source of functional diversity. Owing to their small size and lack of coding capacity, the identification of MITEs has been demanding. However, the increasing availability of reference genomes and bioinformatic tools provides better means for the genome-wide identification and analysis of MITEs and for the elucidation of their contribution to the evolution of plant genomes. We mined MITEs in the carrot reference genome DH1 using MITE-hunter and developed a curated carrot MITE repository comprising 428 families. Of the 31,025 MITE copies spanning 10.34 Mbp of the carrot genome, 54% were positioned in genic regions. Stowaways and Tourists were frequently present in the vicinity of genes, while Mutator-like MITEs were relatively more enriched in introns. hAT-like MITEs were relatively more frequently associated with transcribed regions, including untranslated regions (UTRs). Some carrot MITE families were shared with other Apiaceae species. We showed that hAT-like MITEs were involved in the formation of new splice variants of insertion-harboring genes. Thus, carrot MITEs contributed to the accretion of new diversity by altering transcripts and possibly affecting the regulation of many genes.

Polee: RNA-Seq analysis using approximate likelihood

The analysis of mRNA transcript abundance with RNA-Seq is a central tool in molecular biology research, but often analyses fail to account for the uncertainty in these estimates, which can be significant, especially when trying to disentangle isoforms or duplicated genes. Preserving uncertainty necessitates a full probabilistic model of the all the sequencing reads, which quickly becomes intractable, as experiments can consist of billions of reads. To overcome these limitations, we propose a new method of approximating the likelihood function of a sparse mixture model, using a technique we call the Pólya tree transformation. We demonstrate that substituting this approximation for the real thing achieves most of the benefits with a fraction of the computational costs, leading to more accurate detection of differential transcript expression and transcript coexpression.

Quantifying splice-site usage: a simple yet powerful approach to analyze splicing

RNA splicing, and variations in this process referred to as alternative splicing, are critical aspects of gene regulation in eukaryotes. From environmental responses in plants to being a primary link between genetic variation and disease in humans, splicing differences confer extensive phenotypic changes across diverse organisms (1–3). Regulation of splicing occurs through differential selection of splice sites in a splicing reaction, which results in variation in the abundance of isoforms and/or splicing events. However, genomic determinants that influence splice-site selection remain largely unknown. While traditional approaches for analyzing splicing rely on quantifying variant transcripts (i.e. isoforms) or splicing events (i.e. intron retention, exon skipping etc.) (4), recent approaches focus on analyzing complex/mutually exclusive splicing patterns (5–8). However, none of these approaches explicitly measure individual splice-site usage, which can provide valuable information about splice-site choice and its regulation. Here, we present a simple approach to quantify the empirical usage of individual splice sites reflecting their strength, which determines their selection in a splicing reaction. Splice-site strength/usage, as a quantitative phenotype, allows us to directly link genetic variation with usage of individual splice-sites. We demonstrate the power of this approach in defining the genomic determinants of splice-site choice through GWAS. Our pilot analysis with more than a thousand splice sites hints that sequence divergence in cis rather than trans is associated with variations in splicing among accessions of Arabidopsis thaliana. This approach allows deciphering principles of splicing and has broad implications from agriculture to medicine.

Human photoreceptor cells from different macular subregions have distinct transcriptional profiles

The human neural retina is a light sensitive tissue with remarkable spatial and cellular organization. Compared with the periphery, the central retina contains more densely packed cone photoreceptor cells with unique morphologies and synaptic wiring. Some regions of the central retina exhibit selective degeneration or preservation in response to retinal disease and the basis for this variation is unknown. In this study, we used both bulk and single-cell RNA sequencing to compare gene expression within concentric regions of the central retina. We identified unique gene expression patterns of foveal cone photoreceptor cells, including many foveal-enriched transcription factors. In addition, we found that the genes RORB1, PPFIA1 and KCNAB2 are differentially spliced in the foveal, parafoveal and macular regions. These results provide a highly detailed spatial characterization of the retinal transcriptome and highlight unique molecular features of different retinal regions.

Gene expression and alternative splicing dynamics are perturbed in female head transcriptomes following heterospecific copulation

BACKGROUND

Despite the growing interest in the female side of copulatory interactions, the roles played by differential expression and alternative splicing mechanisms of pre-RNA on tissues outside of the reproductive tract have remained largely unknown. Here we addressed these questions in the context of con- vs heterospecific matings between Drosophila mojavensis and its sister species, D. arizonae. We analyzed transcriptional responses in female heads using an integrated investigation of genome-wide patterns of gene expression, including differential expression (DE), alternative splicing (AS) and intron retention (IR).

RESULTS

Our results indicated that early transcriptional responses were largely congruent between con- and heterospecific matings but are substantially perturbed over time. Conspecific matings induced functional pathways related to amino acid balance previously associated with the brain's physiology and female postmating behavior. Heterospecific matings often failed to activate regulation of some of these genes and induced expression of additional genes when compared with those of conspecifically-mated females. These mechanisms showed functional specializations with DE genes mostly linked to pathways of proteolysis and nutrient homeostasis, while AS genes were more related to photoreception and muscle assembly pathways. IR seems to play a more general role in DE regulation during the female postmating response.

CONCLUSIONS

We provide evidence showing that AS genes substantially perturbed by heterospecific matings in female heads evolve at slower evolutionary rates than the genome background. However, DE genes evolve at evolutionary rates similar, or even higher, than those of male reproductive genes, which highlights their potential role in sexual selection and the evolution of reproductive barriers.

Comprehensive Identification and Alternative Splicing of Microexons in Drosophila

Interrupted exons in the pre-mRNA transcripts are ligated together through RNA splicing, which plays a critical role in the regulation of gene expression. Exons with a length ≤ 30 nt are defined as microexons that are unique in identification. However, microexons, especially those shorter than 8 nt, have not been well studied in many organisms due to difficulties in mapping short segments from sequencing reads. Here, we analyzed mRNA-seq data from a variety of Drosophila samples with a newly developed bioinformatic tool, ce-TopHat. In addition to the Flybase annotated, 465 new microexons were identified. Differentially alternatively spliced (AS) microexons were investigated between the Drosophila tissues (head, body, and gonad) and genders. Most of the AS microexons were found in the head and two AS microexons were identified in the sex-determination pathway gene fruitless.

Differential RNA splicing as a potentially important driver mechanism in multiple myeloma

Disruption of the normal splicing patterns of RNA is a major factor in the pathogenesis of a number of diseases. Increasingly research has shown the strong influence that splicing patterns can have on cancer progression. Multiple Myeloma is a molecularly heterogeneous disease classified by the presence of key translocations, gene expression profiles and mutations but the splicing patterns in MM remains largely unexplored. We take a multifaceted approach to define the extent and impact of alternative splicing in MM. We look at the spliceosome component, SF3B1, with hotspot mutations (K700E and K666T/Q) shown to result in an increase in alternative splicing in other cancers. We discovered a number of differentially spliced genes in comparison of the SF3B1 mutant and wild type samples that included, MZB1, DYNLL1, TMEM14C and splicing related genes DHX9, CLASRP, and SNRPE. We identified a broader role for abnormal splicing showing clear differences in the extent of novel splice variants in the different translocation groups. We show that a high number of novel splice loci is associated with adverse survival and an ultra-high risk group. The enumeration of patterns of alternative splicing has the potential to refine MM classification and to aid in the risk stratification of patients.

RUNX1/RUNX1T1 mediates alternative splicing and reorganises the transcriptional landscape in leukemia

The fusion oncogene RUNX1/RUNX1T1 encodes an aberrant transcription factor, which plays a key role in the initiation and maintenance of acute myeloid leukemia. Here we show that the RUNX1/RUNX1T1 oncogene is a regulator of alternative RNA splicing in leukemic cells. The comprehensive analysis of RUNX1/RUNX1T1-associated splicing events identifies two principal mechanisms that underlie the differential production of RNA isoforms: (i) RUNX1/RUNX1T1-mediated regulation of alternative transcription start site selection, and (ii) direct or indirect control of the expression of genes encoding splicing factors. The first mechanism leads to the expression of RNA isoforms with alternative structure of the 5'-UTR regions. The second mechanism generates alternative transcripts with new junctions between internal cassettes and constitutive exons. We also show that RUNX1/RUNX1T1-mediated differential splicing affects several functional groups of genes and produces proteins with unique conserved domain structures. In summary, this study reveals alternative splicing as an important component of transcriptome re-organization in leukemia by an aberrant transcriptional regulator.

Hearing impairment due to Mir183/96/182 mutations suggests both loss and gain of function effects

The microRNA miR-96 is important for hearing, as point mutations in humans and mice result in dominant progressive hearing loss. Mir96 is expressed in sensory cells along with Mir182 and Mir183, but the roles of these closely-linked microRNAs are as yet unknown. Here we analyse mice carrying null alleles of Mir182, and of Mir183 and Mir96 together to investigate their roles in hearing. We found that Mir183/96 heterozygous mice had normal hearing and homozygotes were completely deaf with abnormal hair cell stereocilia bundles and reduced numbers of inner hair cell synapses at four weeks old. Mir182 knockout mice developed normal hearing then exhibited progressive hearing loss. Our transcriptional analyses revealed significant changes in a range of other genes, but surprisingly there were fewer genes with altered expression in the organ of Corti of Mir183/96 null mice compared with our previous findings in Mir96 Dmdo mutants, which have a point mutation in the miR-96 seed region. This suggests the more severe phenotype of Mir96 Dmdo mutants compared with Mir183/96 mutants, including progressive hearing loss in Mir96 Dmdo heterozygotes, is likely to be mediated by the gain of novel target genes in addition to the loss of its normal targets. We propose three mechanisms of action of mutant miRNAs; loss of targets that are normally completely repressed, loss of targets whose transcription is normally buffered by the miRNA, and gain of novel targets. Any of these mechanisms could lead to a partial loss of a robust cellular identity and consequent dysfunction.

Tracing key genes associated with the Pinctada margaritifera albino phenotype from juvenile to cultured pearl harvest stages using multiple whole transcriptome sequencing

BACKGROUND

Albino mutations are commonly observed in the animal kingdom, including in bivalves. In the black-lipped pearl oyster Pinctada margaritifera, albino specimens are characterized by total or partial absence of colouration resulting in typical white shell phenotype expression. The relationship of shell colour with resulting cultured pearl colour is of great economic interest in P. margaritifera, on which a pearl industry is based. Hence, the albino phenotype provides a useful way to examine the molecular mechanisms underlying pigmentation.

RESULTS

Whole transcriptome RNA-sequencing analysis comparing albino and black wild-type phenotypes at three stages over the culture cycle of P. margaritifera revealed a total of 1606, 798 and 187 differentially expressed genes in whole juvenile, adult mantle and pearl sac tissue, respectively. These genes were found to be involved in five main molecular pathways, tightly linked to known pigmentation pathways: melanogenesis, calcium signalling pathway, Notch signalling pathway, pigment transport and biomineralization. Additionally, significant phenotype-associated SNPs were selected (N = 159), including two located in the Pif biomineralization gene, which codes for nacre formation. Interestingly, significantly different transcript splicing was detected between juvenile (N = 1366) and adult mantle tissue (N = 313) in, e.g., the tyrosinase Tyr-1 gene, which showed more complex regulation in mantle, and the Notch1 encoding gene, which was upregulated in albino juveniles.

CONCLUSION

This multiple RNA-seq approach provided new knowledge about genes associated with the P. margaritifera albino phenotype, highlighting: 1) new molecular pathways, such as the Notch signalling pathway in pigmentation, 2) associated SNP markers with biomineraliszation gene of interest like Pif for marker-assisted selection and prevention of inbreeding, and 3) alternative gene splicing for melanin biosynthesis implicating tyrosinase.

Transcriptomic analysis reveals differential gene expression, alternative splicing, and novel exons during mouse trophoblast stem cell differentiation

BACKGROUND

Differentiation of mouse trophoblast stem cells (TSCs) to trophoblast giant cells (TGCs) has been widely used as a model system to study placental development and function. While several differentially expressed genes, including regulators of TSC differentiation, have been identified, a comprehensive analysis of the global expression of genes and splice variants in the two cell types has not been reported.

RESULTS

Here, we report ~ 7800 differentially expressed genes in TGCs compared to TSCs which include regulators of the cell cycle, apoptosis, cytoskeleton, cell mobility, embryo implantation, metabolism, and various signaling pathways. We show that several mitotic proteins, including Aurora A kinase, were downregulated in TGCs and that the activity of Aurora A kinase is required for the maintenance of TSCs. We also identify hitherto undiscovered, cell-type specific alternative splicing events in 31 genes in the two cell types. Finally, we also report 19 novel exons in 12 genes which are expressed in both TSCs and TGCs.

CONCLUSIONS

Overall, our results uncover several potential regulators of TSC differentiation and TGC function, thereby providing a valuable resource for developmental and molecular biologists interested in the study of stem cell differentiation and embryonic development.

Alternative splicing is highly variable among Daphnia pulex lineages in response to acute copper exposure

BACKGROUND

Despite being one of the primary mechanisms of gene expression regulation in eukaryotes, alternative splicing is often overlooked in ecotoxicogenomic studies. The process of alternative splicing facilitates the production of multiple mRNA isoforms from a single gene thereby greatly increasing the diversity of the transcriptome and proteome. This process can be important in enabling the organism to cope with stressful conditions. Accurate identification of splice sites using RNA sequencing requires alignment to independent exonic positions within the genome, presenting bioinformatic challenges, particularly when using short read data. Although technological advances allow for the detection of splicing patterns on a genome-wide scale, very little is known about the extent of intraspecies variation in splicing patterns, particularly in response to environmental stressors. In this study, we used RNA-sequencing to study the molecular responses to acute copper exposure in three lineages of Daphnia pulex by focusing on the contribution of alternative splicing in addition to gene expression responses.

RESULTS

By comparing the overall gene expression and splicing patterns among all 15 copper-exposed samples and 6 controls, we identified 588 differentially expressed (DE) genes and 16 differentially spliced (DS) genes. Most of the DS genes (13) were not found to be DE, suggesting unique transcriptional regulation in response to copper that went unnoticed with conventional DE analysis. To understand the influence of genetic background on gene expression and alternative splicing responses to Cu, each of the three lineages was analyzed separately. In contrast to the overall analysis, each lineage had a higher proportion of unique DS genes than DE genes suggesting that genetic background has a larger influence on DS than on DE. Gene Ontology analysis revealed that some pathways involved in stress response were jointly regulated by DS and DE genes while others were regulated by only transcription or only splicing.

CONCLUSIONS

Our findings suggest an important role for alternative splicing in shaping transcriptome diversity in response to metal exposure in Daphnia, highlighting the importance of integrating splicing analyses with gene expression surveys to characterize molecular pathways in evolutionary and environmental studies.

Splicing factor DHX15 affects tp53 and mdm2 expression via alternate splicing and promoter usage

DHX15, a DEAH box containing RNA helicase, is a splicing factor required for the last step of splicing. Recent studies identified a recurrent mutational hotspot, R222G, in DHX15 in ∼ 6% of acute myeloid leukemia (AML) patients that carry the fusion protein RUNX1-RUNX1T1 produced by t (8;21) (q22;q22). Studies using yeast mutants showed that substitution of G for the residue equivalent to R222 leads to loss of its helicase function, suggesting that it is a loss-of-function mutation. To elucidate the role of DHX15 during development, we established the first vertebrate knockout model with CRISPR/Cas9 in zebrafish. Our data showed that dhx15 expression is enriched in the brain, eyes, pectoral fin primordia, liver and intestinal bulb during embryonic development. Dhx15 deficiency leads to pleiotropic morphological phenotypes in homozygous mutant embryos starting at 3 days post fertilization (dpf) that result in lethality by 7 dpf, revealing an essential role during embryonic development. RNA-seq analysis suggested important roles of Dhx15 in chromatin and nucleosome assembly and regulation of the Mdm2-p53 pathway. Interestingly, exons corresponding to the alternate transcriptional start sites for tp53 and mdm2 were preferentially expressed in the mutant embryos, leading to significant upregulation of their alternate isoforms, Δ113p53 (orthologous to Δ133p53 isoform in human) and mdm2-P2 (isoform using distal promoter P2), respectively. We speculate that these alterations in the Mdm2-p53 pathway contribute to the development of AML in patients with t(8;21) and somatically mutated DHX15.

Genome-Wide Characterization of Alternative Splicing Events and Their Responses to Cold Stress in Tilapia

Alternative splicing (AS) is an important post-transcriptional regulatory mechanism for cells to generate transcript variability and proteome diversity. No systematic investigation of AS events among different tissues in response to stressors is available for tilapia currently. In this study, AS among different tissues was identified and the cold stress-related AS events were explored in a Nile tilapia (Oreochromis niloticus) line based on 42 RNA-seq datasets using a bioinformatics pipeline. 14,796 (82.76%; SD = 2,840) of the expression genes showed AS events. The two most abundant AS types were alternative transcription start site (TSS) and terminal site (TTS) in tilapia. Testis, brain and kidney possess the most abundant AS gene number, while the blood, muscle and liver possess the least number in each tissue. Furthermore, 208 differentially alternative splicing (DAS) genes in heart and 483 DAS in brain in response to cold stress. The number of AS types for alternative exon end, exon skipping and retention of single intron increased significantly under cold stress. GO enrichment and pathway overrepresentation analysis indicated that many DAS genes, e.g., genes in circadian clock pathway, may influence expression of downstream genes under cold stress. Our study revealed that AS exists extensively in tilapia and plays an important role in cold adaption.

Systematic evaluation of differential splicing tools for RNA-seq studies

Differential splicing (DS) is a post-transcriptional biological process with critical, wide-ranging effects on a plethora of cellular activities and disease processes. To date, a number of computational approaches have been developed to identify and quantify differentially spliced genes from RNA-seq data, but a comprehensive intercomparison and appraisal of these approaches is currently lacking. In this study, we systematically evaluated 10 DS analysis tools for consistency and reproducibility, precision, recall and false discovery rate, agreement upon reported differentially spliced genes and functional enrichment. The tools were selected to represent the three different methodological categories: exon-based (DEXSeq, edgeR, JunctionSeq, limma), isoform-based (cuffdiff2, DiffSplice) and event-based methods (dSpliceType, MAJIQ, rMATS, SUPPA). Overall, all the exon-based methods and two event-based methods (MAJIQ and rMATS) scored well on the selected measures. Of the 10 tools tested, the exon-based methods performed generally better than the isoform-based and event-based methods. However, overall, the different data analysis tools performed strikingly differently across different data sets or numbers of samples.

Visualization and analysis of RNA-Seq assembly graphs

RNA-Seq is a powerful transcriptome profiling technology enabling transcript discovery and quantification. Whilst most commonly used for gene-level quantification, the data can be used for the analysis of transcript isoforms. However, when the underlying transcript assemblies are complex, current visualization approaches can be limiting, with splicing events a challenge to interpret. Here, we report on the development of a graph-based visualization method as a complementary approach to understanding transcript diversity from short-read RNA-Seq data. Following the mapping of reads to a reference genome, a read-to-read comparison is performed on all reads mapping to a given gene, producing a weighted similarity matrix between reads. This is used to produce an RNA assembly graph, where nodes represent reads and edges similarity scores between them. The resulting graphs are visualized in 3D space to better appreciate their sometimes large and complex topology, with other information being overlaid on to nodes, e.g. transcript models. Here we demonstrate the utility of this approach, including the unusual structure of these graphs and how they can be used to identify issues in assembly, repetitive sequences within transcripts and splice variants. We believe this approach has the potential to significantly improve our understanding of transcript complexity.

Broad regulation of gene isoform expression by Wnt signaling in cancer

Differential gene isoform expression is a ubiquitous mechanism to enhance proteome diversity and maintain cell homeostasis. Mechanisms such as splicing that drive gene isoform variability are highly dynamic and responsive to changes in cell signaling pathways. Wnt/β-catenin signaling has profound effects on cell activity and cell fate and is known to modify several splicing events by altering the expression of individual splicing factors. However, a global assessment of how extensively Wnt signaling regulates splicing and other mechanisms that determine mRNA isoform composition in cancer is lacking. We used deep time-resolved RNA-seq in two independent in vivo Wnt-addicted tumor models during treatment with the potent Wnt inhibitor ETC-159 and examined Wnt regulated splicing events and splicing regulators. We found 1025 genes that underwent Wnt regulated variable exon usage leading to isoform expression changes. This was accompanied by extensive Wnt regulated changes in the expression of splicing regulators. Many of these Wnt regulated events were conserved in multiple human cancers, and many were linked to previously defined cancer-associated splicing quantitative trait loci. This suggests that the Wnt regulated splicing events are components of fundamental oncogenic processes. These findings demonstrate the wide-ranging effects of Wnt signaling on the isoform composition of the cell and provides an extensive resource of expression changes of splicing regulators and gene isoforms regulated by Wnt signaling.

Are Small Nucleolar RNAs "CRISPRable"? A Report on Box C/D Small Nucleolar RNA Editing in Human Cells

CRISPR technologies are nowadays widely used for targeted knockout of numerous protein-coding genes and for the study of various processes and metabolic pathways in human cells. Most attention in the genome editing field is now focused on the cleavage of protein-coding genes or genes encoding long non-coding RNAs (lncRNAs), while the studies on targeted knockout of intron-encoded regulatory RNAs are sparse. Small nucleolar RNAs (snoRNAs) present a class of non-coding RNAs encoded within the introns of various host genes and involved in post-transcriptional maturation of ribosomal RNAs (rRNAs) in eukaryotic cells. Box C/D snoRNAs direct 2'-O-methylation of rRNA nucleotides. These short RNAs have specific elements in their structure, namely, boxes C and D, and a target-recognizing region. Here, we present the study devoted to CRISPR/Cas9-mediated editing of box C/D snoRNA genes in Gas5. We obtained monoclonal cell lines carrying mutations in snoRNA genes and analyzed the levels of the mutant box C/D snoRNA as well as the 2'-O-methylation status of the target rRNA nucleotide in the obtained cells. Mutations in SNORD75 in the obtained monoclonal cell line were shown to result in aberrant splicing of Gas5 with exclusion of exons 3 to 5, which was confirmed by RT-PCR and RNA-Seq. The obtained results suggest that SNORD75 contains an element for binding of some factors regulating maturation of Gas5 pre-lncRNA. We suggest that METTL3/METTL14 is among such factors, and m⁶A-methylation pathways are involved in regulation of Gas5 splicing. Our results shell light on the role of SNORDs in regulating splicing of the host gene.

Unstable TTTTA/TTTCA expansions in MARCH6 are associated with Familial Adult Myoclonic Epilepsy type 3

Familial Adult Myoclonic Epilepsy (FAME) is a genetically heterogeneous disorder characterized by cortical tremor and seizures. Intronic TTTTA/TTTCA repeat expansions in SAMD12 (FAME1) are the main cause of FAME in Asia. Using genome sequencing and repeat-primed PCR, we identify another site of this repeat expansion, in MARCH6 (FAME3) in four European families. Analysis of single DNA molecules with nanopore sequencing and molecular combing show that expansions range from 3.3 to 14 kb on average. However, we observe considerable variability in expansion length and structure, supporting the existence of multiple expansion configurations in blood cells and fibroblasts of the same individual. Moreover, the largest expansions are associated with micro-rearrangements occurring near the expansion in 20% of cells. This study provides further evidence that FAME is caused by intronic TTTTA/TTTCA expansions in distinct genes and reveals that expansions exhibit an unexpectedly high somatic instability that can ultimately result in genomic rearrangements.

PRMT5 methylome profiling uncovers a direct link to splicing regulation in acute myeloid leukemia

Protein arginine methyltransferase 5 (PRMT5) has emerged as a promising cancer drug target, and three PRMT5 inhibitors are currently in clinical trials for multiple malignancies. In this study, we investigated the role of PRMT5 in human acute myeloid leukemia (AML). Using an enzymatic dead version of PRMT5 and a PRMT5-specific inhibitor, we demonstrated the requirement of the catalytic activity of PRMT5 for the survival of AML cells. We then identified PRMT5 substrates using multiplexed quantitative proteomics and investigated their role in the survival of AML cells. We found that the function of the splicing regulator SRSF1 relies on its methylation by PRMT5 and that loss of PRMT5 leads to changes in alternative splicing of multiple essential genes. This explains the requirement of PRMT5 for leukemia cell survival. We show that PRMT5 regulates binding of SRSF1 to mRNAs and proteins and provide potential biomarkers for the treatment response to PRMT5 inhibitors.

PathwaySplice: an R package for unbiased pathway analysis of alternative splicing in RNA-Seq data

Summary

Pathway analysis of alternative splicing would be biased without accounting for the different number of exons or junctions associated with each gene, because genes with higher number of exons or junctions are more likely to be included in the 'significant' gene list in alternative splicing. We present PathwaySplice, an R package that (i) Performs pathway analysis that explicitly adjusts for the number of exons or junctions associated with each gene; (ii) visualizes selection bias due to different number of exons or junctions for each gene and formally tests for presence of bias using logistic regression; (iii) supports gene sets based on the Gene Ontology terms, as well as more broadly defined gene sets (e.g. MSigDB) or user defined gene sets; (iv) identifies the significant genes driving pathway significance and (v) organizes significant pathways with an enrichment map, where pathways with large number of overlapping genes are grouped together in a network graph.

Availability and implementation

https://bioconductor.org/packages/release/bioc/html/PathwaySplice.html.

Supplementary information

Supplementary data are available at Bioinformatics online.

RNA sequencing of transplant-stage idiopathic pulmonary fibrosis lung reveals unique pathway regulation

Idiopathic pulmonary fibrosis (IPF), the scarring of lung parenchyma resulting in the loss of lung function, remains a fatal disease with a significant unmet medical need. Patients with severe IPF often develop acute exacerbations resulting in the rapid deterioration of lung function, requiring transplantation. Understanding the pathophysiological mechanisms contributing to IPF is key to develop novel therapeutic approaches for end-stage disease. We report here RNA-sequencing analyses of lung tissues from a cohort of patients with transplant-stage IPF (n=36), compared with acute lung injury (ALI) (n=11) and nondisease controls (n=19), that reveal a robust gene expression signature unique to end-stage IPF. In addition to extracellular matrix remodelling pathways, we identified pathways associated with T-cell infiltration/activation, tumour development, and cholesterol homeostasis, as well as novel alternatively spliced transcripts that are differentially regulated in the advanced IPF lung versus ALI or nondisease controls. Additionally, we show a subset of genes that are correlated with percent predicted forced vital capacity and could reflect disease severity. Our results establish a robust transcriptomic fingerprint of an advanced IPF lung that is distinct from previously reported microarray signatures of moderate, stable or progressive IPF and identifies hitherto unknown candidate targets and pathways for therapeutic intervention in late-stage IPF as well as biomarkers to characterise disease progression and enable patient stratification.

DIEGO: detection of differential alternative splicing using Aitchison's geometry

Motivation

Alternative splicing is a biological process of fundamental importance in most eukaryotes. It plays a pivotal role in cell differentiation and gene regulation and has been associated with a number of different diseases. The widespread availability of RNA-Sequencing capacities allows an ever closer investigation of differentially expressed isoforms. However, most tools for differential alternative splicing (DAS) analysis do not take split reads, i.e. the most direct evidence for a splice event, into account. Here, we present DIEGO, a compositional data analysis method able to detect DAS between two sets of RNA-Seq samples based on split reads.

Results

The python tool DIEGO works without isoform annotations and is fast enough to analyze large experiments while being robust and accurate. We provide python and perl parsers for common formats.

Availability and implementation

The software is available at: www.bioinf.uni-leipzig.de/Software/DIEGO.

Supplementary information

Supplementary data are available at Bioinformatics online.

Altered Levels of Proteins and Phosphoproteins, in the Absence of Early Causative Transcriptional Changes, Shape the Molecular Pathogenesis in the Brain of Young Presymptomatic Ki91 SCA3/MJD Mouse

Spinocerebellar ataxia type 3 (SCA3/MJD) is a polyQ neurodegenerative disease where the presymptomatic phase of pathogenesis is unknown. Therefore, we investigated the molecular network of transcriptomic and proteomic triggers in young presymptomatic SCA3/MJD brain from Ki91 knock-in mouse. We found that transcriptional dysregulations resulting from mutant ataxin-3 are not occurring in young Ki91 mice, while old Ki91 mice and also postmitotic patient SCA3 neurons demonstrate the late transcriptomic changes. Unlike the lack of early mRNA changes, we have identified numerous early changes of total proteins and phosphoproteins in 2-month-old Ki91 mouse cortex and cerebellum. We discovered the network of processes in presymptomatic SCA3 with three main groups of disturbed processes comprising altered proteins: (I) modulation of protein levels and DNA damage (Pabpc1, Ddb1, Nedd8), (II) formation of neuronal cellular structures (Tubb3, Nefh, p-Tau), and (III) neuronal function affected by processes following perturbed cytoskeletal formation (Mt-Co3, Stx1b, p-Syn1). Phosphoproteins downregulate in the young Ki91 mouse brain and their phosphosites are associated with kinases that interact with ATXN3 such as casein kinase, Camk2, and kinases controlled by another Atxn3 interactor p21 such as Gsk3, Pka, and Cdk kinases. We conclude that the onset of SCA3 pathology occurs without altered transcript level and is characterized by changed levels of proteins responsible for termination of translation, DNA damage, spliceosome, and protein phosphorylation. This disturbs global cellular processes such as cytoskeleton and transport of vesicles and mitochondria along axons causing energy deficit and neurodegeneration also manifesting in an altered level of transcripts at later ages.

Induction of interferon-stimulated genes and cellular stress pathways by morpholinos in zebrafish

The phenotypes caused by morpholino-mediated interference of gene function in zebrafish are often not observed in the corresponding mutant(s). We took advantage of the availability of a relatively large collection of transcriptomic datasets to identify common signatures that characterize morpholino-injected animals (morphants). In addition to the previously reported activation of tp53 expression, we observed increased expression of the interferon-stimulated genes (ISGs), isg15 and isg20, the cell death pathway gene casp8, and other cellular stress response genes including phlda3, mdm2 and gadd45aa. Studies involving segmentation stage embryos were more likely to show upregulation of these genes. We also found that the expression of these genes could be upregulated by increasing doses of an egfl7 morpholino, or even high doses of the standard control morpholino. Thus, these data show that morpholinos can induce the expression of ISGs in zebrafish embryos and further our understanding of morpholino effects.

Aberrant Clonal Hematopoiesis following Lentiviral Vector Transduction of HSPCs in a Rhesus Macaque

Lentiviral vectors (LVs) are used for delivery of genes into hematopoietic stem and progenitor cells (HSPCs) in clinical trials worldwide. LVs, in contrast to retroviral vectors, are not associated with insertion site-associated malignant clonal expansions and, thus, are considered safer. Here, however, we present a case of markedly abnormal dysplastic clonal hematopoiesis affecting the erythroid, myeloid, and megakaryocytic lineages in a rhesus macaque transplanted with HSPCs that were transduced with a LV containing a strong retroviral murine stem cell virus (MSCV) constitutive promoter-enhancer in the LTR. Nine insertions were mapped in the abnormal clone, resulting in overexpression and aberrant splicing of several genes of interest, including the cytokine stem cell factor and the transcription factor PLAG1. This case represents the first clear link between lentiviral insertion-induced clonal expansion and a clinically abnormal transformed phenotype following transduction of normal primate or human HSPCs, which is concerning, and suggests that strong constitutive promoters should not be included in LVs.

The LL-100 panel: 100 cell lines for blood cancer studies

For many years, immortalized cell lines have been used as model systems for cancer research. Cell line panels were established for basic research and drug development, but did not cover the full spectrum of leukemia and lymphoma. Therefore, we now developed a novel panel (LL-100), 100 cell lines covering 22 entities of human leukemia and lymphoma including T-cell, B-cell and myeloid malignancies. Importantly, all cell lines are unequivocally authenticated and assigned to the correct tissue. Cell line samples were proven to be free of mycoplasma and non-inherent virus contamination. Whole exome sequencing and RNA-sequencing of the 100 cell lines were conducted with a uniform methodology to complement existing data on these publicly available cell lines. We show that such comprehensive sequencing data can be used to find lymphoma-subtype-characteristic copy number aberrations, mRNA isoforms, transcription factor activities and expression patterns of NKL homeobox genes. These exemplary studies confirm that the novel LL-100 panel will be useful for understanding the function of oncogenes and tumor suppressor genes and to develop targeted therapies.

Analysis of differential gene expression and alternative splicing is significantly influenced by choice of reference genome

RNA-seq analysis has enabled the evaluation of transcriptional changes in many species including nonmodel organisms. However, in most species only a single reference genome is available and RNA-seq reads from highly divergent varieties are typically aligned to this reference. Here, we quantify the impacts of the choice of mapping genome in rice where three high-quality reference genomes are available. We aligned RNA-seq data from a popular productive rice variety to three different reference genomes and found that the identification of differentially expressed genes differed depending on which reference genome was used for mapping. Furthermore, the ability to detect differentially used transcript isoforms was profoundly affected by the choice of reference genome: Only 30% of the differentially used splicing features were detected when reads were mapped to the more commonly used, but more distantly related reference genome. This demonstrated that gene expression and splicing analysis varies considerably depending on the mapping reference genome, and that analysis of individuals that are distantly related to an available reference genome may be improved by acquisition of new genomic reference material. We observed that these differences in transcriptome analysis are, in part, due to the presence of single nucleotide polymorphisms between the sequenced individual and each respective reference genome, as well as annotation differences between the reference genomes that exist even between syntenic orthologs. We conclude that even between two closely related genomes of similar quality, using the reference genome that is most closely related to the species being sampled significantly improves transcriptome analysis.

The C. elegans embryonic transcriptome with tissue, time, and alternative splicing resolution

We have used RNA-seq in Caenorhabditis elegans to produce transcription profiles for seven specific embryonic cell populations from gastrulation to the onset of terminal differentiation. The expression data for these seven cell populations, covering major cell lineages and tissues in the worm, reveal the complex and dynamic changes in gene expression, both spatially and temporally. Also, within genes, start sites and exon usage can be highly differential, producing transcripts that are specific to developmental periods or cell lineages. We have also found evidence of novel exons and introns, as well as differential usage of SL1 and SL2 splice leaders. By combining this data set with the modERN ChIP-seq resource, we are able to support and predict gene regulatory relationships. The detailed information on differences and similarities between gene expression in cell lineages and tissues should be of great value to the community and provides a framework for the investigation of expression in individual cells.

Medullary Respiratory Circuit Is Reorganized by a Seasonally-Induced Program in Preparation for Hibernation

Deep hibernators go through several cycles of profound drops in body temperature during the winter season, with core temperatures sometimes reaching near freezing. Yet unlike non-hibernating mammals, they can sustain breathing rhythms. The physiological processes that make this possible are still not understood. In this study, we focused on the medullary Ventral Respiratory Column of a facultative hibernator, the Syrian hamster. Using shortened day-lengths, we induced a "winter-adapted" physiological state, which is a prerequisite for hibernation. When recording electrophysiological signals from acute slices in the winter-adapted pre-Bötzinger complex (preBötC), spike trains showed higher spike rates, amplitudes, complexity, as well as higher temperature sensitivity, suggesting an increase in connectivity and/or synaptic strength during the winter season. We further examined action potential waveforms and found that the depolarization integral, as measured by the area under the curve, is selectively enhanced in winter-adapted animals. This suggests that a shift in the ion handling kinetics is also being induced by the winter-adaptation program. RNA sequencing of respiratory pre-motor neurons, followed by gene set enrichment analysis, revealed differential regulation and splicing in structural, synaptic, and ion handling genes. Splice junction analysis suggested that differential exon usage is occurring in a select subset of ion handling subunits (ATP1A3, KCNC3, SCN1B), and synaptic structure genes (SNCB, SNCG, RAB3A). Our findings show that the hamster respiratory center undergoes a seasonally-cued alteration in electrophysiological properties, likely protecting against respiratory failure at low temperatures.

EGFR Controls Hair Shaft Differentiation in a p53-Independent Manner

Epidermal growth factor receptor (EGFR) signaling controls skin development and homeostasis in mice and humans, and its deficiency causes severe skin inflammation, which might affect epidermal stem cell behavior. Here, we describe the inflammation-independent effects of EGFR deficiency during skin morphogenesis and in adult hair follicle stem cells. Expression and alternative splicing analysis of RNA sequencing data from interfollicular epidermis and outer root sheath indicate that EGFR controls genes involved in epidermal differentiation and also in centrosome function, DNA damage, cell cycle, and apoptosis. Genetic experiments employing p53 deletion in EGFR-deficient epidermis reveal that EGFR signaling exhibits p53-dependent functions in proliferative epidermal compartments, as well as p53-independent functions in differentiated hair shaft keratinocytes. Loss of EGFR leads to absence of LEF1 protein specifically in the innermost epithelial hair layers, resulting in disorganization of medulla cells. Thus, our results uncover important spatial and temporal features of cell-autonomous EGFR functions in the epidermis.

Bioinformatics challenges and perspectives when studying the effect of epigenetic modifications on alternative splicing

It is widely known that epigenetic modifications are important in regulating transcription, but several have also been reported in alternative splicing. The regulation of pre-mRNA splicing is important to explain proteomic diversity and the misregulation of splicing has been implicated in many diseases. Here, we give a brief overview of the role of epigenetics in alternative splicing and disease. We then discuss the bioinformatics methods that can be used to model interactions between epigenetic marks and regulators of splicing. These models can be used to identify alternative splicing and epigenetic changes across different phenotypes.

This article is part of a discussion meeting issue ‘Frontiers in epigenetic chemical biology’.

RNA-sequencing in ophthalmology research: considerations for experimental design and analysis

High-throughput, massively parallel sequence analysis has revolutionized the way that researchers design and execute scientific investigations. Vast amounts of sequence data can be generated in short periods of time. Regarding ophthalmology and vision research, extensive interrogation of patient samples for underlying causative DNA mutations has resulted in the discovery of many new genes relevant to eye disease. However, such analysis remains functionally limited. RNA-sequencing accurately snapshots thousands of genes, capturing many subtypes of RNA molecules, and has become the gold standard for transcriptome gene expression quantification. RNA-sequencing has the potential to advance our understanding of eye development and disease; it can reveal new candidates to improve our molecular diagnosis rates and highlight therapeutic targets for intervention. But with a wide range of applications, the design of such experiments can be problematic, no single optimal pipeline exists, and therefore, several considerations must be undertaken for optimal study design. We review the key steps involved in RNA-sequencing experimental design and the downstream bioinformatic pipelines used for differential gene expression. We provide guidance on the application of RNA-sequencing to ophthalmology and sources of open-access eye-related data sets.

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.

Expanding the Boundaries of RNA Sequencing as a Diagnostic Tool for Rare Mendelian Disease

Gene-panel and whole-exome analyses are now standard methodologies for mutation detection in Mendelian disease. However, the diagnostic yield achieved is at best 50%, leaving the genetic basis for disease unsolved in many individuals. New approaches are thus needed to narrow the diagnostic gap. Whole-genome sequencing is one potential strategy, but it currently has variant-interpretation challenges, particularly for non-coding changes. In this study we focus on transcriptome analysis, specifically total RNA sequencing (RNA-seq), by using monogenetic neuromuscular disorders as proof of principle. We examined a cohort of 25 exome and/or panel "negative" cases and provided genetic resolution in 36% (9/25). Causative mutations were identified in coding and non-coding exons, as well as in intronic regions, and the mutational pathomechanisms included transcriptional repression, exon skipping, and intron inclusion. We address a key barrier of transcriptome-based diagnostics: the need for source material with disease-representative expression patterns. We establish that blood-based RNA-seq is not adequate for neuromuscular diagnostics, whereas myotubes generated by transdifferentiation from an individual's fibroblasts accurately reflect the muscle transcriptome and faithfully reveal disease-causing mutations. Our work confirms that RNA-seq can greatly improve diagnostic yield in genetically unresolved cases of Mendelian disease, defines strengths and challenges of the technology, and demonstrates the suitability of cell models for RNA-based diagnostics. Our data set the stage for development of RNA-seq as a powerful clinical diagnostic tool that can be applied to the large population of individuals with undiagnosed, rare diseases and provide a framework for establishing minimally invasive strategies for doing so.

Whole transcriptome sequencing and biomineralization gene architecture associated with cultured pearl quality traits in the pearl oyster, Pinctada margaritifera

Background

Cultured pearls are unique gems produced by living organisms, mainly molluscs of the Pinctada genus, through the biomineralization properties of pearl sac tissue. Improvement of P. margaritifera pearl quality is one of the biggest challenges that Polynesian research has faced to date. To achieve this goal, a better understanding of the complex mechanisms related to nacre and pearl formation is essential and can now be approached through the use of massive parallel sequencing technologies. The aim of this study was to use RNA-seq to compare whole transcriptome expression of pearl sacs that had producing pearls with high and low quality. For this purpose, a comprehensive reference transcriptome of P. margaritifera was built based on multi-tissue sampling (mantle, gonad, whole animal), including different living stages (juvenile, adults) and phenotypes (colour morphotypes, sex).

Results

Strikingly, few genes were found to be up-regulated for high quality pearls (n = 16) compared to the up-regulated genes in low quality pearls (n = 246). Biomineralization genes up-regulated in low quality pearls were specific to prismatic and prism-nacre layers. Alternative splicing was further identified in several key biomineralization genes based on a recent P. margaritifera draft genome.

Conclusion

This study lifts the veil on the multi-level regulation of biomineralization genes associated with pearl quality determination.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5443-5) contains supplementary material, which is available to authorized users.

Correction of gene model annotations improves isoform abundance estimates: the example of ketohexokinase ( Khk)

Next generation sequencing protocols such as RNA-seq have made the genome-wide characterization of the transcriptome a crucial part of many research projects in biology. Analyses of the resulting data provide key information on gene expression and in certain cases on exon or isoform usage. The emergence of transcript quantification software such as Salmon has enabled researchers to efficiently estimate isoform and gene expressions across the genome while tremendously reducing the necessary computational power. Although overall gene expression estimations were shown to be accurate, isoform expression quantifications appear to be a more challenging task. Low expression levels and uneven or insufficient coverage were reported as potential explanations for inconsistent estimates. Here, through the example of the ketohexokinase ( Khk) gene in mouse, we demonstrate that the use of an incorrect gene annotation can also result in erroneous isoform quantification results. Manual correction of the input Khk gene model provided a much more accurate estimation of relative Khk isoform expression when compared to quantitative PCR (qPCR measurements). In particular, removal of an unexpressed retained intron and a proper adjustment of the 5’ and 3’ untranslated regions both had a strong impact on the correction of erroneous estimates. Finally, we observed a better concordance in isoform quantification between datasets and sequencing strategies when relying on the newly generated Khk annotations. These results highlight the importance of accurate gene models and annotations for correct isoform quantification and reassert the need for orthogonal methods of estimation of isoform expression to confirm important findings.

Identification and characterization of long non-coding RNAs in porcine granulosa cells exposed to 2,3,7,8-tetrachlorodibenzo-p-dioxin

Background

Long non-coding RNAs (lncRNAs) may regulate gene expression in numerous biological processes including cellular response to xenobiotics. The exposure of living organisms to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), a persistent environmental contaminant, results in reproductive defects in many species including pigs. The aims of the study were to identify and characterize lncRNAs in porcine granulosa cells as well as to examine the effects of TCDD on the lncRNA expression profile in the cells.

Results

One thousand six hundred sixty-six lncRNAs were identified and characterized in porcine granulosa cells. The identified lncRNAs were found to be shorter than mRNAs. In addition, the number of exons was lower in lncRNAs than in mRNAs and their exons were longer. TCDD affected the expression of 22 lncRNAs (differentially expressed lncRNAs [DELs]; log₂ fold change  ≥ 1, P-adjusted < 0.05) in the examined cells. Potential functions of DELs were indirectly predicted via searching their target cis- and trans-regulated protein-coding genes. The co-expression analysis revealed that DELs may influence the expression of numerous genes, including those involved in cellular response to xenobiotics, dioxin metabolism, endoplasmic reticulum stress and cell proliferation. Aryl hydrocarbon receptor (AhR) and cytochrome P450 1A1 (CYP1A1) were found among the trans-regulated genes.

Conclusions

These findings indicate that the identified lncRNAs may constitute a part of the regulatory mechanism of TCDD action in granulosa cells. To our knowledge, this is the first study describing lncRNAs in porcine granulosa cells as well as TCDD effects on the lncRNA expression profile. These results may trigger new research directions leading to better understanding of molecular processes induced by xenobiotics in the ovary.