Genomewide analysis of mRNA processing in yeast using splicing-specific microarrays

A reliable and quick method for screening alternative splicing variants for low-abundance genes

Alternative splicing (AS) is a universal phenomenon in eukaryotes, and it is still challenging to identify AS events. Several methods have been developed to identify AS events, such as expressed sequence tags (EST), microarrays and RNA-seq. However, EST has limitations in identifying low-abundance genes, while microarray and RNA-seq are high-throughput technologies, and PCR-based technology is needed for validation. To overcome the limitations of EST and shortcomings of high-throughput technologies, we established a method to identify AS events, especially for low-abundance genes, by reverse transcription (RT) PCR with gene-specific primers (GSPs) followed by nested PCR. This process includes two major steps: 1) the use of GSPs to amplify as long as the specific gene segment and 2) multiple rounds of nested PCR to screen the AS and confirm the unknown splicing variants. With this method, we successfully identified three new splicing variants, namely, GenBank Accession No. HM623886 for the bdnf gene (GenBank GeneID: 12064), GenBank Accession No. JF417977 for the trkc gene (GenBank GeneID: 18213) and GenBank Accession No. HM623888 for the glb-18 gene (GenBank GeneID: 172485). In addition to its reliability and simplicity, the method is also cost-effective and labor-intensive. In conclusion, we developed an RT-nested PCR method using gene-specific primers to efficiently identify known and novel AS variants. This approach overcomes the limitations of existing methods for detecting rare transcripts. By enabling the discovery of new isoforms, especially for low-abundance genes, this technique can aid research into aberrant splicing in disease. Future studies can apply this method to uncover AS variants involved in cancer, neurodegeneration, and other splicing-related disorders.

The Schizosaccharomyces pombe DEAD-box protein Mss116 is required for mitoribosome assembly and mitochondrial translation

DEAD-box helicases are important players in mitochondrial gene expression, which is necessary for mitochondrial respiration. In this study, we characterized Schizosaccharomyces pombe Mss116 (spMss116), a member of the family of DEAD-box RNA helicases. Deletion of spmss116 in a mitochondrial intron-containing background significantly reduced the levels of mitochondrial DNA (mtDNA)-encoded cox1 and cob1 mRNAs and impaired mitochondrial translation, leading to a severe respiratory defect and a loss of cell viability during stationary phase. Deletion of mitochondrial introns restored the levels of cox1 and cob1 mRNAs to wide-type (WT) levels but could not restore mitochondrial translation and respiration in Δspmss116 cells. Furthermore, deletion of spmss116 in both mitochondrial intron-containing and intronless backgrounds impaired mitoribosome assembly and destabilization of mitoribosomal proteins. Our findings suggest that defective mitochondrial translation caused by deletion of spmss116 is most likely due to impaired mitoribosome assembly.

Pan-cancer analysis of PLAU indicates its potential prognostic value and correlation with neutrophil infiltration in BLCA

BACKGROUND

PLAU is known as a selected serine protease converting plasminogen to plasmin. The role of PLAU in the development of pan-cancer, especially bladder urothelial carcinoma (BLCA) remains unclear.

METHOD

A variety of online tools and cancer databases, including TCGA, GETx, HPA database, GSCALite, UALCAN, ESTIMATE, CIBERSORT, ssGSEA algorithms and SangerBox website, were applied to investigate the associations between PLAU expression and prognosis, genetic alterations, pathway activation, and tumor immunity in pan-cancer. Through cBioPortal and STITCH platforms, the oncogenic role of PLAU and related targeting medicines in BLCA were also explored. We verified the expression of PLAU in pan-cancer cells and its function in bladder cancer cell lines using wet-lab experiments.

RESULTS

PLAU expression levels were significantly higher in most cancer tissues. PLAU had a certain accuracy in the diagnosis of various types of cancers (90 % AUC > 0.700). In BLCA, PLAU has abundant methylated sites and showed statistical differences in clinical features. PLAU was involved in tumor immune infiltration, and especially positively correlated with neutrophil infiltration. High-expressed PLAU indicated poorer prognosis in the BLCA patients receiving Atezolizumab. A high mRNA and protein expression levels of PLAU were observed in pan-cancer cell lines, especially BLCA cells. Knockdown of PLAU inhibited the invasive, proliferative, and aggressive phenotypes of bladder cancer cells. Immunohistochemical staining validated PLAU's higher expression in BLCA tissues than in adjacent non-cancerous tissues. And overexpression of PLAU was associated with more advanced TNM stage, and high infiltrating depth.

CONCLUSION

Our study revealed that PLAU can serve as a potential therapeutic target and prognostic marker for various malignancies, especially BLCA.

Truncating the spliceosomal 'rope protein' Prp45 results in Htz1 dependent phenotypes

Spliceosome assembly contributes an important but incompletely understood aspect of splicing regulation. Prp45 is a yeast splicing factor which runs as an extended fold through the spliceosome, and which may be important for bringing its components together. We performed a whole genome analysis of the genetic interaction network of the truncated allele of PRP45 (prp45(1-169)) using synthetic genetic array technology and found chromatin remodellers and modifiers as an enriched category. In agreement with related studies, H2A.Z-encoding HTZ1, and the components of SWR1, INO80, and SAGA complexes represented prominent interactors, with htz1 conferring the strongest growth defect. Because the truncation of Prp45 disproportionately affected low copy number transcripts of intron-containing genes, we prepared strains carrying intronless versions of SRB2, VPS75, or HRB1, the most affected cases with transcription-related function. Intron removal from SRB2, but not from the other genes, partly repaired some but not all the growth phenotypes identified in the genetic screen. The interaction of prp45(1-169) and htz1Δ was detectable even in cells with SRB2 intron deleted (srb2Δi). The less truncated variant, prp45(1-330), had a synthetic growth defect with htz1Δ at 16°C, which also persisted in the srb2Δi background. Moreover, htz1Δ enhanced prp45(1-330) dependent pre-mRNA hyper-accumulation of both high and low efficiency splicers, genes ECM33 and COF1, respectively. We conclude that while the expression defects of low expression intron-containing genes contribute to the genetic interactome of prp45(1-169), the genetic interactions between prp45 and htz1 alleles demonstrate the sensitivity of spliceosome assembly, delayed in prp45(1-169), to the chromatin environment.

Genomic and transcriptomic advances in amyotrophic lateral sclerosis

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disorder and the most common motor neuron disease. ALS shows substantial clinical and molecular heterogeneity. In vitro and in vivo models coupled with multiomic techniques have provided important contributions to unraveling the pathomechanisms underlying ALS. To date, despite promising results and accumulating knowledge, an effective treatment is still lacking. Here, we provide an overview of the literature on the use of genomics, epigenomics, transcriptomics and microRNAs to deeply investigate the molecular mechanisms developing and sustaining ALS. We report the most relevant genes implicated in ALS pathogenesis, discussing the use of different high-throughput sequencing techniques and the role of epigenomic modifications. Furthermore, we present transcriptomic studies discussing the most recent advances, from microarrays to bulk and single-cell RNA sequencing. Finally, we discuss the use of microRNAs as potential biomarkers and promising tools for molecular intervention. The integration of data from multiple omic approaches may provide new insights into pathogenic pathways in ALS by shedding light on diagnostic and prognostic biomarkers, helping to stratify patients into clinically relevant subgroups, revealing novel therapeutic targets and supporting the development of new effective therapies.

Upcoming progress of transcriptomics studies on plants: An overview

Transcriptome sequencing or RNA-Sequencing is a high-resolution, sensitive and high-throughput next-generation sequencing (NGS) approach used to study non-model plants and other organisms. In other words, it is an assembly of RNA transcripts from individual or whole samples of functional and developmental stages. RNA-Seq is a significant technique for identifying gene predictions and mining functional analysis that improves gene ontology understanding mechanisms of biological processes, molecular functions, and cellular components, but there is limited information available on this topic. Transcriptomics research on different types of plants can assist researchers to understand functional genes in better ways and regulatory processes to improve breeding selection and cultivation practices. In recent years, several advancements in RNA-Seq technology have been made for the characterization of the transcriptomes of distinct cell types in biological tissues in an efficient manner. RNA-Seq technologies are briefly introduced and examined in terms of their scientific applications. In a nutshell, it introduces all transcriptome sequencing and analysis techniques, as well as their applications in plant biology research. This review will focus on numerous existing and forthcoming strategies for improving transcriptome sequencing technologies for functional gene mining in various plants using RNA- Seq technology, based on the principles, development, and applications.

DNA-enabled fluorescent-based nanosensors monitoring tumor-related RNA toward advanced cancer diagnosis: A review

As a burgeoning non-invasive indicator for reproducible cancer diagnosis, tumor-related biomarkers have a wide range of applications in early cancer screening, efficacy monitoring, and prognosis predicting. Accurate and efficient biomarker determination, therefore, is of great importance to prevent cancer progression at an early stage, thus reducing the disease burden on the entire population, and facilitating advanced therapies for cancer. During the last few years, various DNA structure-based fluorescent probes have established a versatile platform for biological measurements, due to their inherent biocompatibility, excellent capacity to recognize nucleic and non-nucleic acid targets, obvious accessibility to synthesis as well as chemical modification, and the ease of interfacing with signal amplification protocols. After decades of research, DNA fluorescent probe technology for detecting tumor-related mRNAs has gradually grown to maturity, especially the advent of fluorescent nanoprobes has taken the process to a new level. Here, a systematic introduction to recent trends and advances focusing on various nanomaterials-related DNA fluorescent probes and the physicochemical properties of various involved nanomaterials (such as AuNP, GO, MnO2, SiO2, AuNR, etc.) are also presented in detail. Further, the strengths and weaknesses of existing probes were described and their progress in the detection of tumor-related mRNAs was illustrated. Also, the salient challenges were discussed later, with a few potential solutions.

Pan-cancer analysis of mRNA stability for decoding tumour post-transcriptional programs

Measuring mRNA decay in tumours is a prohibitive challenge, limiting our ability to map the post-transcriptional programs of cancer. Here, using a statistical framework to decouple transcriptional and post-transcriptional effects in RNA-seq data, we uncover the mRNA stability changes that accompany tumour development and progression. Analysis of 7760 samples across 18 cancer types suggests that mRNA stability changes are ~30% as frequent as transcriptional events, highlighting their widespread role in shaping the tumour transcriptome. Dysregulation of programs associated with >80 RNA-binding proteins (RBPs) and microRNAs (miRNAs) drive these changes, including multi-cancer inactivation of RBFOX and miR-29 families. Phenotypic activation or inhibition of RBFOX1 highlights its role in calcium signaling dysregulation, while modulation of miR-29 shows its impact on extracellular matrix organization and stemness genes. Overall, our study underlines the integral role of mRNA stability in shaping the cancer transcriptome, and provides a resource for systematic interrogation of cancer-associated stability pathways.

The role of mRNA stability in shaping the cancer transcriptome is revealed using a statistical analysis of transcriptomic data.

Automation of high-throughput mRNA-seq library preparation: a robust, hands-free and time efficient methodology

Over the last decade, whole transcriptome profiling, also known as RNA-sequencing (RNA-seq), has quickly gained traction as a reliable method for unbiased assessment of gene expression. Integration of RNA-seq expression data into other omics datasets (e.g., proteomics, metabolomics, or epigenetics) solidifies our understanding of cell-specific regulatory patterns, yielding pathways to investigate the key rules of gene regulation. A limitation to efficient, at-scale utilization of RNA-seq is the time-demanding library preparation workflows, which is a 2-day or longer endeavor per cohort/sample size. To tackle this bottleneck, we designed an automated workflow that increases throughput capacity, while minimizing human error to enhance reproducibility. To this end, we converted the manual protocol of the NEBNext Directional Ultra II RNA Library Prep Kit for Illumina on the Beckman Coulter liquid handler, Biomek i7 Hybrid workstation. A total of 84 RNA samples were isolated from two human cell lines and subjected to comparative manual and automated library preparation methods. Qualitative and quantitative results indicated a high degree of similarity between libraries generated manually or through automation. Yet, there was a significant reduction in both hands-on and assay time from a 2-day manual to a 9-hour automated workflow. Using linear regression analysis, we found the Pearson correlation coefficient between libraries generated manually or by automation to be almost identical to a sample being sequenced twice (R²= 0.985 vs 0.983). This demonstrates that high-throughput automated workflows can be of great benefit to genomic laboratories by enhancing efficiency of library preparation, reducing hands-on time and increasing throughput potential.

Dynamic tracking of p21 mRNA in living cells by sticky-flares for the visual evaluation of the tumor treatment effect

Monitoring the expression level of the intracellular tumor suppressor gene p21 mRNA is essential to reveal the progress and prognosis of a tumor. Methods widely reported for the detection of p21 mRNA are the real-time polymerase chain reaction and Northern blot. However, these methods only detect mRNA in vitro and cannot realize the in situ monitoring of the p21 mRNA expression level in living cells. Additionally, the sensor for the real-time tracking and monitoring of the p21 mRNA location without the help of a transfection reagent in living cells is still limited. Herein, a novel sticky-flare was constructed for the dynamic monitoring of the temporal and spatial variations of p21 mRNA in living cells. The nanoprobe consists of AuNP, a recognition sequence modified with Cy5, and a thiol-modified DNA sequence. The thiol oligonucleotide strand could act partially complementary to the Cy5-modified oligonucleotide strand to form a double-stranded DNA linked to AuNP, resulting in the fluorescence quenching of Cy5 due to the energy transfer from Cy5 to the gold sphere. In the presence of p21 mRNA, the Cy5-modified recognition nucleic acid specifically bound to p21 mRNA to form a more stable double chain and escaped from the gold sphere, leading to the recovery of red fluorescence. Our method is better than other methods in its ability to quantify the spatial distribution and expression level of p21 mRNA in living cells and discriminate various tumor cell lines with different p21 mRNA expression levels by the naked eye. Particularly, the sticky-flare probe used in this assay could allow the visual evaluation of the tumor treatment effect and the determination of the tumor progression stage by enabling monitoring of the relative expression level of p21 mRNA in tumor cells after cisplatin treatment. The method reported here is accurate, reliable and needs no auxiliary tools (transfection reagent), and thereby provides a promising route for the prognostic evaluation and drug development of cancer treatment in the future.

DNA-coated gold nanoparticles for tracking hepatocyte growth factor secreted by transplanted mesenchymal stem cells in pulmonary fibrosis therapy

The identification of paracrine factors secreted by transplanted mesenchymal stem cells (MSCs) during the treatment of idiopathic pulmonary fibrosis (IPF) is essential for understanding the role of MSCs in therapy. Herein, we report a facile and efficient strategy for in vivo tracking the secretion of hepatocyte growth factor (HGF) in MSCs during IPF therapy. In our strategy, a novel nanoflare tracer consisting of gold nanoparticles (AuNPs), complementary sequences and dye-labeled recognition sequences is developed. Briefly, the AuNPs are functionalized with oligonucleotide complementary sequences hybridized to the organic dye-labeled recognition sequences, where the organic fluorophores are in close proximity to the AuNPs. In the absence of targets, the dye and AuNPs are separated from each other, inducing the quenching of the fluorescence signal. However, in the presence of targets, the recognition sequences gradually fall off from the AuNPs, causing the fluorescence signal to rise. In brief, in vivo monitoring of the dynamic expression of HGF mRNA in transplanted MSCs during IPF therapy in the current work may provide new insight into the paracrine process of the transplanted MSCs, thereby advancing the MSC-based IPF therapy toward clinical applications.

A broad analysis of splicing regulation in yeast using a large library of synthetic introns

RNA splicing is a key process in eukaryotic gene expression, in which an intron is spliced out of a pre-mRNA molecule to eventually produce a mature mRNA. Most intron-containing genes are constitutively spliced, hence efficient splicing of an intron is crucial for efficient regulation of gene expression. Here we use a large synthetic oligo library of ~20,000 variants to explore how different intronic sequence features affect splicing efficiency and mRNA expression levels in S. cerevisiae. Introns are defined by three functional sites, the 5’ donor site, the branch site, and the 3’ acceptor site. Using a combinatorial design of synthetic introns, we demonstrate how non-consensus splice site sequences in each of these sites affect splicing efficiency. We then show that S. cerevisiae splicing machinery tends to select alternative 3’ splice sites downstream of the original site, and we suggest that this tendency created a selective pressure, leading to the avoidance of cryptic splice site motifs near introns’ 3’ ends. We further use natural intronic sequences from other yeast species, whose splicing machineries have diverged to various extents, to show how intron architectures in the various species have been adapted to the organism’s splicing machinery. We suggest that the observed tendency for cryptic splicing is a result of a loss of a specific splicing factor, U2AF1. Lastly, we show that synthetic sequences containing two introns give rise to alternative RNA isoforms in S. cerevisiae, demonstrating that merely a synthetic fusion of two introns might be suffice to facilitate alternative splicing in yeast. Our study reveals novel mechanisms by which introns are shaped in evolution to allow cells to regulate their transcriptome. In addition, it provides a valuable resource to study the regulation of constitutive and alternative splicing in a model organism.

RNA splicing is a process in which parts of a new pre-mRNA are spliced out of the mRNA molecule to produce eventually a mature mRNA. Those RNA segments that are spliced out are termed introns, and they are found in most genes in eukaryotic organisms. Hence regulation of this process has a major role in the control of gene expression.

The budding yeast S. cerevisiae is a popular model organism for eukaryotic cell biology, but in terms of splicing it differs, as it has only few intron-containing genes. Nevertheless, this species has been used to study basic principles of splicing regulation based on its ~300 introns. Here we used the technology of a large synthetic genetic library to introduce many new intron-containing genes to the yeast genome, to explore splicing regulation at a wider scope than was possible so far.

Reassuringly, our results confirm known regulatory mechanisms, and further expand our understanding of splicing regulation, specifically how the yeast splicing machinery interacts with the end of introns, and how through evolution introns have evolved to avoid unwanted misidentifications of this end. We further demonstrate the potential of the yeast splicing machinery to alternatively splice a two-intron gene, which is common in other eukaryotes but rare in yeast. Our work presents a first-of-its-kind resource for the systematic study of splicing in live cells.

Environmental robustness of the global yeast genetic interaction network

Phenotypes associated with genetic variants can be altered by interactions with other genetic variants (GxG), with the environment (GxE), or both (GxGxE). Yeast genetic interactions have been mapped on a global scale, but the environmental influence on the plasticity of genetic networks has not been examined systematically. To assess environmental rewiring of genetic networks, we examined 14 diverse conditions and scored 30,000 functionally representative yeast gene pairs for dynamic, differential interactions. Different conditions revealed novel differential interactions, which often uncovered functional connections between distantly related gene pairs. However, the majority of observed genetic interactions remained unchanged in different conditions, suggesting that the global yeast genetic interaction network is robust to environmental perturbation and captures the fundamental functional architecture of a eukaryotic cell.

Omics in Systems Biology: Current Progress and Future Outlook

Biological research has undergone tremendous changes over the past three decades. Research used to almost exclusively focus on a single aspect of a single molecule per experiment. Modern technologies have enabled thousands of molecules to be simultaneously analyzed and the way that these molecules influence each other to be discerned. The change is so dramatic that it has given rise to a whole new descriptive suffix (i.e., omics) to describe these fields of study. While genomics was arguably the initial driver of this new trend, it quickly spread to other biological entities resulting in the creation of transcriptomics, proteomics, metabolomics, etc. The development of these "big four omics" created a wave of other omic fields, such as epigenomics, glycomics, lipidomics, microbiomics, and even foodomics; all with the purpose of comprehensively studying all the molecular entities or processes within their respective domain. The large number of omic fields that are invented even led to the term "panomics" as a way to classify them all under one category. Ultimately, all of these omic fields are setting the foundation for developing systems biology; in which the focus will be on determining the complex interactions that occur within biological systems.

Processing and Analysis of RNA-seq Data from Public Resources

Advances in next generation sequencing (NGS) technologies resulted in a broad array of large-scale gene expression studies and an unprecedented volume of whole messenger RNA (mRNA) sequencing data, or the transcriptome (also known as RNA sequencing, or RNA-seq). These include the Genotype Tissue Expression project (GTEx) and The Cancer Genome Atlas (TCGA), among others. Here we cover some of the commonly used datasets, provide an overview on how to begin the analysis pipeline, and how to explore and interpret the data provided by these publicly available resources.

A random forest based biomarker discovery and power analysis framework for diagnostics research

Background

Biomarker identification is one of the major and important goal of functional genomics and translational medicine studies. Large scale –omics data are increasingly being accumulated and can provide vital means for the identification of biomarkers for the early diagnosis of complex disease and/or for advanced patient/diseases stratification. These tasks are clearly interlinked, and it is essential that an unbiased and stable methodology is applied in order to address them. Although, recently, many, primarily machine learning based, biomarker identification approaches have been developed, the exploration of potential associations between biomarker identification and the design of future experiments remains a challenge.

Methods

In this study, using both simulated and published experimentally derived datasets, we assessed the performance of several state-of-the-art Random Forest (RF) based decision approaches, namely the Boruta method, the permutation based feature selection without correction method, the permutation based feature selection with correction method, and the backward elimination based feature selection method. Moreover, we conducted a power analysis to estimate the number of samples required for potential future studies.

Results

We present a number of different RF based stable feature selection methods and compare their performances using simulated, as well as published, experimentally derived, datasets. Across all of the scenarios considered, we found the Boruta method to be the most stable methodology, whilst the Permutation (Raw) approach offered the largest number of relevant features, when allowed to stabilise over a number of iterations. Finally, we developed and made available a web interface (https://joelarkman.shinyapps.io/PowerTools/) to streamline power calculations thereby aiding the design of potential future studies within a translational medicine context.

Conclusions

We developed a RF-based biomarker discovery framework and provide a web interface for our framework, termed PowerTools, that caters the design of appropriate and cost-effective subsequent future omics study.

Schizosaccharomyces pombe Ppr10 is required for mitochondrial translation

The mitochondrial genome encodes key components of the oxidative phosphorylation (OXPHOS) system, whose expression is essential for mitochondrial functions. We have previously shown that deletion of the Schizosaccharomyces pombe ppr10 encoding a pentatricopeptide repeat protein severely reduces the mature levels of intron-containing mitochondrial transcripts cox1 and cob1, and severely impairs mitochondrial translation. In this study, we examined the possibility that the reduced levels of Cox1 and Cob1 proteins in cells were due to lowered levels of cox1 and cob1 mRNAs. We found that deletion of ppr10 did not affect the levels of mature cox1 and cob1 mRNAs in a mitochondrial intronless background. However, synthesis of Cox1 and Cob1 proteins were still severely affected by deletion of ppr10 in a mitochondrial intronless background. Consistent with this, we found that deletion of mitochondrial introns could not rescue the respiratory growth defect of Δppr10 cells. Our results reveal that Ppr10 is not required for the stability of cox1 and cob1 mRNAs, and provide further support for the idea that Ppr10 plays a critical role in mitochondrial translation.

Dbr1 functions in mRNA processing, intron turnover and human diseases

Pre-mRNA processing and mRNA stability play direct roles in controlling protein abundance in a cell. Before the mRNA can be translated into a protein, the introns in the pre-mRNA transcripts need to be removed by splicing, such that exons can be ligated together and can code for a protein. In this process, the function of the RNA lariat debranching enzyme or Dbr1 provides a rate-limiting step in the intron turnover process and possibly regulating the production of translation competent mRNAs. Surprising new roles of Dbr1 are emerging in cellular metabolism which extends beyond intron turnover processes, ranging from splicing regulation to translational control. In this review, we highlight the importance of the Dbr1 enzyme, its structure and how anomalies in its function could relate to various human diseases.

Widespread Transcriptional Readthrough Caused by Nab2 Depletion Leads to Chimeric Transcripts with Retained Introns

Nascent RNA sequencing has revealed that pre-mRNA splicing can occur shortly after introns emerge from RNA polymerase II (RNA Pol II). Differences in co-transcriptional splicing profiles suggest regulation by cis- and/or trans-acting factors. Here, we use single-molecule intron tracking (SMIT) to identify a cohort of regulators by machine learning in budding yeast. Of these, Nab2 displays reduced co-transcriptional splicing when depleted. Unexpectedly, these splicing defects are attributable to aberrant "intrusive" transcriptional readthrough from upstream genes, as revealed by long-read sequencing. Transcripts that originate from the intron-containing gene's own transcription start site (TSS) are efficiently spliced, indicating no direct role of Nab2 in splicing per se. This work highlights the coupling between transcription, splicing, and 3' end formation in the context of gene organization along chromosomes. We conclude that Nab2 is required for proper 3' end processing, which ensures gene-specific control of co-transcriptional RNA processing.

Inhibitory Mechanism of Trichoderma virens ZT05 on Rhizoctonia solani

Trichoderma is a filamentous fungus that is widely distributed in nature. As a biological control agent of agricultural pests, Trichoderma species have been widely studied in recent years. This study aimed to understand the inhibitory mechanism of Trichoderma virens ZT05 on Rhizoctonia solani through the side-by-side culture of T. virens ZT05 and R. solani. To this end, we investigated the effect of volatile and nonvolatile metabolites of T. virens ZT05 on the mycelium growth and enzyme activity of R. solani and analyzed transcriptome data collected from side-by-side culture. T. virens ZT05 has a significant antagonistic effect against R. solani. The mycelium of T. virens ZT05 spirally wraps around and penetrates the mycelium of R. solani and inhibits the growth of R. solani. The volatile and nonvolatile metabolites of T. virens ZT05 have significant inhibitory effects on the growth of R. solani. The nonvolatile metabolites of T. virens ZT05 significantly affect the mycelium proteins of R. solani, including catalase (CAT), superoxide dismutase (SOD), peroxidase (POD), selenium-dependent glutathione peroxidase (GSH-Px), soluble proteins, and malondialdehyde (MDA). Twenty genes associated with hyperparasitism, including extracellular proteases, oligopeptide transporters, G-protein coupled receptors (GPCRs), chitinases, glucanases, and proteases were found to be upregulated during the antagonistic process between T. virens ZT05 and R. solani. Thirty genes related to antibiosis function, including tetracycline resistance proteins, reductases, the heat shock response, the oxidative stress response, ATP-binding cassette (ABC) efflux transporters, and multidrug resistance transporters, were found to be upregulated during the side-by-side culture of T. virens ZT05 and R. solani. T. virens ZT05 has a significant inhibitory effect on R. solani, and its mechanism of action is associated with hyperparasitism and antibiosis.

CORAZON: a web server for data normalization and unsupervised clustering based on expression profiles

Objective

Data normalization and clustering are mandatory steps in gene expression and downstream analyses, respectively. However, user-friendly implementations of these methodologies are available exclusively under expensive licensing agreements, or in stand-alone scripts developed, reflecting on a great obstacle for users with less computational skills.

Results

We developed an online tool called CORAZON (Correlations Analyses Zipper Online), which implements three unsupervised learning methods to cluster gene expression datasets in a friendly environment. It allows the usage of eight gene expression normalization/transformation methodologies and the attribute’s influence. The normalizations requiring the gene length only could be performed to RNA-seq, meanwhile the others can be used with microarray and/or NanoString data. Clustering methodologies performances were evaluated through five models with accuracies between 92 and 100%. We applied our tool to obtain functional insights of non-coding RNAs (ncRNAs) based on Gene Ontology enrichment of clusters in a dataset generated by the ENCODE project. The clusters where the majority of transcripts are coding genes were enriched in Cellular, Metabolic, Transports, and Systems Development categories. Meanwhile, the ncRNAs were enriched in the Detection of Stimulus, Sensory Perception, Immunological System, and Digestion categories. CORAZON source-code is freely available at https://gitlab.com/integrativebioinformatics/corazon and the web-server can be accessed at http://corazon.integrativebioinformatics.me.

Adipogenesis, Osteogenesis, and Chondrogenesis of Human Mesenchymal Stem/Stromal Cells: A Comparative Transcriptome Approach

Adipogenesis, osteogenesis and chondrogenesis of human mesenchymal stem/stromal cells (MSC) are complex and highly regulated processes. Over the years, several studies have focused on understanding the mechanisms involved in the MSC commitment to the osteogenic, adipogenic and/or chondrogenic phenotypes. High-throughput methodologies have been used to investigate the gene expression profile during differentiation. Association of data analysis of mRNAs, microRNAs, circular RNAs and long non-coding RNAs, obtained at different time points over these processes, are important to depict the complexity of differentiation. This review will discuss the results that were highlighted in transcriptome analyses of MSC undergoing adipogenic, osteogenic and chondrogenic differentiation. The focus is to shed light on key molecules, main signaling pathways and biological processes related to different time points of adipogenesis, osteogenesis and chondrogenesis.

Amniotic fluid cell-free transcriptome: a glimpse into fetal development and placental cellular dynamics during normal pregnancy

Background

The amniotic fluid (AF) cell-free transcriptome is modulated by physiologic and pathologic processes during pregnancy. AF gene expression changes with advancing gestation reflect fetal development and organ maturation; yet, defining normal expression and splicing patterns for biomarker discovery in obstetrics requires larger heterogeneous cohorts, evaluation of potential confounding factors, and novel analytical approaches.

Methods

Women with a normal pregnancy who had an AF sample collected during midtrimester (n = 30) or at term gestation (n = 68) were included. Expression profiling at exon level resolution was performed using Human Transcriptome Arrays. Differential expression was based on moderated t-test adjusted p < 0.05 and fold change > 1.25; for differential splicing, a splicing index > 2 and adjusted p < 0.05 were required. Functional profiling was used to interpret differentially expressed or spliced genes. The expression of tissue-specific and cell-type specific signatures defined by single-cell genomics was quantified and correlated with covariates. In-silico validation studies were performed using publicly available datasets.

Results

1) 64,071 genes were detected in AF, with 11% of the coding and 6% of the non-coding genes being differentially expressed between midtrimester and term gestation. Expression changes were highly correlated with those previously reported (R > 0.79, p < 0.001) and featured increased expression of genes specific to the trachea, salivary glands, and lung and decreased expression of genes specific to the cardiac myocytes, uterus, and fetal liver, among others. 2) Single-cell RNA-seq signatures of the cytotrophoblast, Hofbauer cells, erythrocytes, monocytes, T and B cells, among others, showed complex patterns of modulation with gestation (adjusted p < 0.05). 3) In 17% of the genes detected, we found differential splicing with advancing gestation in genes related to brain development processes and immunity pathways, including some that were missed based on differential expression analysis alone.

Conclusions

This represents the largest AF transcriptomics study in normal pregnancy, reporting for the first time that single-cell genomic signatures can be tracked in the AF and display complex patterns of expression during gestation. We also demonstrate a role for alternative splicing in tissue-identity acquisition, organ development, and immune processes. The results herein may have implications for the development of fetal testing to assess placental function and fetal organ maturity.

RBM17 Interacts with U2SURP and CHERP to Regulate Expression and Splicing of RNA-Processing Proteins

RNA splicing entails the coordinated interaction of more than 150 proteins in the spliceosome, one of the most complex of the cell’s molecular machines. We previously discovered that the RNA-binding motif protein 17 (RBM17), a component of the spliceosome, is essential for survival and cell maintenance. Here, we find that it interacts with the spliceosomal factors U2SURP and CHERP and that they reciprocally regulate each other’s stability, both in mouse and in human cells. Individual knockdown of each of the three proteins induces overlapping changes in splicing and gene expression of transcripts enriched for RNA-processing factors. Our results elucidate the function of RBM17, U2SURP, and CHERP and link the activity of the spliceosome to the regulation of downstream RNA-binding proteins. These data support the hypothesis that, beyond driving constitutive splicing, spliceosomal factors can regulate alternative splicing of specific targets.

De Maio et al. find that the splicing factor RBM17 establishes a physical and functional relation with U2SURP and CHERP. Knockdown of these U2 snRNP-associated spliceosomal components reveals their synergistic activity toward regulation of a given set of transcripts rather than a more predictable transcriptome-wide inhibition of splicing.

Smu1 and RED are required for activation of spliceosomal B complexes assembled on short introns

Human pre-catalytic spliceosomes contain several proteins that associate transiently just prior to spliceosome activation and are absent in yeast, suggesting that this critical step is more complex in higher eukaryotes. We demonstrate via RNAi coupled with RNA-Seq that two of these human-specific proteins, Smu1 and RED, function both as alternative splicing regulators and as general splicing factors and are required predominantly for efficient splicing of short introns. In vitro splicing assays reveal that Smu1 and RED promote spliceosome activation, and are essential for this step when the distance between the pre-mRNA’s 5′ splice site (SS) and branch site (BS) is sufficiently short. This Smu1-RED requirement can be bypassed when the 5′ and 3′ regions of short introns are physically separated. Our observations suggest that Smu1 and RED relieve physical constraints arising from a short 5′SS-BS distance, thereby enabling spliceosomes to overcome structural challenges associated with the splicing of short introns.

Human spliceosome components Smu1 and RED regulate alternative splicing. Here the authors show that Smu1 and RED are also required for constitutive splicing of short introns.

AtBUD13 affects pre-mRNA splicing and is essential for embryo development in Arabidopsis

Pre-mRNA splicing is an important step for gene expression regulation. Yeast Bud13p (bud-site selection protein 13) regulates the budding pattern and pre-mRNA splicing in yeast cells; however, no Bud13p homologs have been identified in plants. Here, we isolated two mutants that carry T-DNA insertions at the At1g31870 locus and shows early embryo lethality and seed abortion. At1g31870 encodes an Arabidopsis homolog of yeast Bud13p, AtBUD13. Although AtBUD13 homologs are widely distributed in eukaryotic organisms, phylogenetic analysis revealed that their protein domain organization is more complex in multicellular species. AtBUD13 is expressed throughout plant development including embryogenesis and AtBUD13 proteins is localized in the nucleus in Arabidopsis. RNA-seq analysis revealed that AtBUD13 mutation predominantly results in the intron retention, especially for shorter introns (≤100 bases). Within this group of genes, we identified 52 genes involved in embryogenesis, out of which 22 are involved in nucleic acid metabolism. Our results demonstrate that AtBUD13 plays critical roles in early embryo development by effecting pre-mRNA splicing.

Postmortem brain tissue as an underutilized resource to study the molecular pathology of neuropsychiatric disorders across different ethnic populations

In recent years, large scale meta-analysis of genome-wide association studies (GWAS) have reliably identified genetic polymorphisms associated with neuropsychiatric disorders such as schizophrenia (SCZ), bipolar disorder (BPD) and major depressive disorder (MDD). However, the majority of disease-associated single nucleotide polymorphisms (SNPs) appear within functionally ambiguous non-coding genomic regions. Recently, increased emphasis has been placed on identifying the functional relevance of disease-associated variants via correlating risk polymorphisms with gene expression levels in etiologically relevant tissues. For neuropsychiatric disorders, the etiologically relevant tissue is brain, which requires robust postmortem sample sizes from varying genetic backgrounds. While small sample sizes are of decreasing concern, postmortem brain databases are composed almost exclusively of Caucasian samples, which significantly limits study design and result interpretation. In this review, we highlight the importance of gene expression and expression quantitative loci (eQTL) studies in clinically relevant postmortem tissue while addressing the current limitations of existing postmortem brain databases. Finally, we introduce future collaborations to develop postmortem brain databases for neuropsychiatric disorders from Chinese and Asian subpopulations.

RNA-Seq approach for accurate characterization of splicing efficiency of yeast introns

Introns in different genes, or even different introns within the same gene, often have different splice sites and differ in splicing efficiency (SE). One expects mass-transcribed genes to have introns with higher SE than weakly transcribed genes. However, such a simple expectation cannot be tested directly because variable SE for these genes is often not measured. Mechanistically, SE should depend on signal strength at key splice sites (SS) such as 5'SS, 3'SS and branchpoint site (BPS), i.e., SE = F(5'SS, 3'SS, BPS). However, without SE, we again cannot model how these splice sites contribute to SE. Here I present an RNA-Seq approach to quantify SE for each of the 304 introns in yeast (Saccharomyces cerevisiae) genes, including 24 in the 5'UTR, by measuring 1) number of reads mapped to exon-exon junctions (N_EE) as a proxy for the abundance of spliced form, and 2) number of reads mapped to exon-intron junction (N_EI5 and N_EI3 at 5' and 3' ends of intron) as a proxy for the abundance of unspliced form. The total mRNA is N_Total = N_EE + p * N_EI5 + (1-p) * N_EI3, with the simplest p = 0.5 but statistical methods were presented to estimate p from data. An estimated p is needed because N_EI5 is expected to be smaller than N_EI3 due to 1) step 1 splicing occurs before step 2 so EI5 is broken before EI3, 2) enrichment of poly(A) mRNA by oligo-dT, and 3) 5' degradation. SE is defined as the proportion (N_EE/N_Total). Application of the method shows that ribosomal protein messages are efficiently and mostly cotranscriptionally spliced. Yeast genes with long introns are also spliced efficiently. HAC1/YFL031W is poorly spliced partly because its splicing involves a nonspliceosome mechanism and partly because Ire1p, which participate in splicing HAC1, is hardly expressed. Many putative yeast genes have low SE, and some splice sites are incorrectly annotated.

Structures of the Catalytically Activated Yeast Spliceosome Reveal the Mechanism of Branching

Pre-mRNA splicing is executed by the spliceosome. Structural characterization of the catalytically activated complex (B^∗) is pivotal for understanding the branching reaction. In this study, we assembled the B^∗ complexes on two different pre-mRNAs from Saccharomyces cerevisiae and determined the cryo-EM structures of four distinct B^∗ complexes at overall resolutions of 2.9-3.8 Å. The duplex between U2 small nuclear RNA (snRNA) and the branch point sequence (BPS) is discretely away from the 5'-splice site (5'SS) in the three B^∗ complexes that are devoid of the step I splicing factors Yju2 and Cwc25. Recruitment of Yju2 into the active site brings the U2/BPS duplex into the vicinity of 5'SS, with the BPS nucleophile positioned 4 Å away from the catalytic metal M2. This analysis reveals the functional mechanism of Yju2 and Cwc25 in branching. These structures on different pre-mRNAs reveal substrate-specific conformations of the spliceosome in a major functional state.

Spliceosome Profiling Visualizes Operations of a Dynamic RNP at Nucleotide Resolution

Tools to understand how the spliceosome functions in vivo have lagged behind advances in the structural biology of the spliceosome. Here, methods are described to globally profile spliceosome-bound pre-mRNA, intermediates, and spliced mRNA at nucleotide resolution. These tools are applied to three yeast species that span 600 million years of evolution. The sensitivity of the approach enables the detection of canonical and non-canonical events, including interrupted, recursive, and nested splicing. This application of statistical modeling uncovers independent roles for the size and position of the intron and the number of introns per transcript in substrate progression through the two catalytic stages. These include species-specific inputs suggestive of spliceosome-transcriptome coevolution. Further investigations reveal the ATP-dependent discard of numerous endogenous substrates after spliceosome assembly in vivo and connect this discard to intron retention, a form of splicing regulation. Spliceosome profiling is a quantitative, generalizable global technology used to investigate an RNP central to eukaryotic gene expression.

Networks of mRNA Processing and Alternative Splicing Regulation in Health and Disease

mRNA processing events introduce an intricate layer of complexity into gene expression processes, supporting a tremendous level of diversification of the genome's coding and regulatory potential, particularly in vertebrate species. The recent development of massive parallel sequencing methods and their adaptation to the identification and quantification of different RNA species and the dynamics of mRNA metabolism and processing has generated an unprecedented view over the regulatory networks that are established at this level, which contribute to sustain developmental, tissue specific or disease specific gene expression programs. In this chapter, we provide an overview of the recent evolution of transcriptome profiling methods and the surprising insights that have emerged in recent years regarding distinct mRNA processing events - from the 5' end to the 3' end of the molecule.

Spliceosome factors target timeless (tim) mRNA to control clock protein accumulation and circadian behavior in Drosophila

Transcription-translation feedback loops that comprise eukaryotic circadian clocks rely upon temporal delays that separate the phase of active transcription of clock genes, such as Drosophila period (per) and timeless (tim), from negative feedback by the two proteins. However, our understanding of the mechanisms involved is incomplete. Through an RNA interference screen, we found that pre-mRNA processing 4 (PRP4) kinase, a component of the U4/U5.U6 triple small nuclear ribonucleoprotein (tri-snRNP) spliceosome, and other tri-snRNP components regulate cycling of the molecular clock as well as rest:activity rhythms. Unbiased RNA-Sequencing uncovered an alternatively spliced intron in tim whose increased retention upon prp4 downregulation leads to decreased TIM levels. We demonstrate that the splicing of tim is rhythmic with a phase that parallels delayed accumulation of the protein in a 24 hr cycle. We propose that alternative splicing constitutes an important clock mechanism for delaying the daily accumulation of clock proteins, and thereby negative feedback by them.

Editorial note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

The prognostic significance of global aberrant alternative splicing in patients with myelodysplastic syndrome

Aberrant alternative splicing (AS) is a hallmark of cancer development. However, there are limited data regarding its clinical implications in myelodysplastic syndrome (MDS). In this study, we performed an in-depth analysis of global AS in 176 primary MDS patients with 20 normal marrow transplant donors as reference. We found that 26.9% of the expressed genes genome-wide were aberrantly spliced in MDS patients compared with normal donors. These aberrant AS genes were related to pathways involved in cell proliferation, cell adhesion and protein degradation. A higher degree of global aberrant AS was associated with male gender and U2AF1 mutation, and predicted shorter overall survival and time to leukemic change. Moreover, it was an independent unfavorable prognostic factor irrespective of age, revised international prognostic scoring system (IPSS-R) risk, and mutations in SRSF2, ZRSR2, ASXL1, TP53, and EZH2. With LASSO-Cox regression method, we constructed a simple prognosis prediction model composed of 13 aberrant AS genes, and demonstrated that it could well stratify MDS patients into distinct risk groups. To our knowledge, this is the first report demonstrating significant prognostic impacts of aberrant splicing on MDS patients. Further prospective studies in larger cohorts are needed to confirm our observations.

RES complex is associated with intron definition and required for zebrafish early embryogenesis

Pre-mRNA splicing is a critical step of gene expression in eukaryotes. Transcriptome-wide splicing patterns are complex and primarily regulated by a diverse set of recognition elements and associated RNA-binding proteins. The retention and splicing (RES) complex is formed by three different proteins (Bud13p, Pml1p and Snu17p) and is involved in splicing in yeast. However, the importance of the RES complex for vertebrate splicing, the intronic features associated with its activity, and its role in development are unknown. In this study, we have generated loss-of-function mutants for the three components of the RES complex in zebrafish and showed that they are required during early development. The mutants showed a marked neural phenotype with increased cell death in the brain and a decrease in differentiated neurons. Transcriptomic analysis of bud13, snip1 (pml1) and rbmx2 (snu17) mutants revealed a global defect in intron splicing, with strong mis-splicing of a subset of introns. We found these RES-dependent introns were short, rich in GC and flanked by GC depleted exons, all of which are features associated with intron definition. Using these features, we developed and validated a predictive model that classifies RES dependent introns. Altogether, our study uncovers the essential role of the RES complex during vertebrate development and provides new insights into its function during splicing.

Gene Regulatory Network Perturbation by Genetic and Epigenetic Variation

Gene regulatory networks underlie biological function and cellular physiology. Alternative splicing (AS) is a fundamental step in gene regulatory networks and plays a key role in development and disease. In addition to the identification of aberrant AS events, an increasing number of studies are focusing on molecular determinants of AS, including genetic and epigenetic regulators. We review here recent efforts to identify various deregulated AS events as well as their molecular determinants that alter biological functions, and discuss clinical features of AS and their druggable potential.

Poly-adenine-mediated fluorescent spherical nucleic acid probes for live-cell imaging of endogenous tumor-related mRNA

Identification of tumor-related mRNA in living cells hold great promise for early cancer diagnosis and pathological research. Herein, we present poly-adenine (polyA)-mediated fluorescent spherical nucleic acid (FSNA) probes for intracellular mRNA detection with regulable sensitivities by programmably adjusting the loading density of DNA on gold nano-interface. Gold nanoparticles (AuNPs) functionalized with polyA-tailed recognition sequences were hybridized to fluorescent "reporter" strands to fabricate fluorescence-quenched FSNA probes. While exposed to target gene, the "reporter" strands were released from FSNA through strand displacement and fluorescence was recovered. With polyA20 tail as the attaching block, the detection limit of FSNA probes was calculated to be 0.31 nM, which is ~55 fold lower than that of thiolated probes without surface density regulation. Quantitative intracellular mRNA detection and imaging could be achieved with polyA-mediated FSNA probes within 2 hours, indicating their application potential in rapid and sensitive intracellular target imaging.

Transcriptome and alternative splicing analysis of nucleus pulposus cells in response to high oxygen tension: Involvement of high oxygen tension in the pathogenesis of intervertebral disc degeneration

High oxygen tension caused by neovascularization in the microenvironment of intervertebral discs (IVDs) is associated with the pathogenesis of IVD degeneration (IDD). Pre-mRNAs undergo alternative splicing (AS) to produce structurally and functionally diverse mRNA and proteins. However, the precise role of high oxygen tension in IDD and the relationship between AS and high oxygen tension in disc cells remain unknown. To investigate the effect of high oxygen tension on disc cells, Affymetrix Rat Transcriptome Array 1.0 was used to determine differentially expressed genes (DEGs) and alternative splicing genes (ASGs) in rat nucleus pulposus (NP) cells treated with 20% O₂. NP cells at 1% O₂ served as the control. PCR was used for validation. GO and KEGG pathway analysis was performed. Furthermore, the reactive oxygen species (ROS) production, growth, cell cycle and matrix metabolism of NP cells were also investigated. In total, 2499 DEGs and 8451 ASGs were identified. Various GO terms and KEGG pathways were potently associated with IDD, including autophagy, mTOR signaling pathway and angiogenesis. Especially, high oxygen tension increased ROS production in NP cells. It also accelerated the matrix metabolism of NP cells and induced NP cell cycle arrest to retard cell growth. This study, for the first time, analyzes the transcriptome and AS of NP cells in response to high oxygen tension, indicating that high oxygen tension is involved in the establishment and progression of IDD through its wide effects on the viability and function of disc cells.

The fission yeast Schizosaccharomyces pombe Mtf2 is required for mitochondrial cox1 gene expression

Mitochondrial gene expression is essential for adenosine triphosphate synthesis via oxidative phosphorylation, which is the universal energy currency of cells. Here, we report the identification and characterization of a homologue of Saccharomyces cerevisiae Mtf2 (also called Nam1) in Schizosaccharomyces pombe. The Δmtf2 mutant with the intron-containing mitochondrial DNA (mtDNA) exhibited impaired growth on a rich medium containing the non-fermentable carbon source glycerol, suggesting that mtf2 is involved in mitochondrial function. mtf2 deletion in a mitochondrial intron-containing background resulted in a barely detectable level of the cox1 mRNA and a reduction in the level of the cob1 mRNA, and severely impaired cox1 translation. In contrast, mtf2 deletion in a mitochondrial intron-less background did not affect the levels of cox1 and cob1 mRNAs. However, Cox1 synthesis could not be restored to the control level in the Δmtf2 mutant with intron-less mtDNA. Our results suggest that unlike its counterpart in S. cerevisiae which plays a general role in synthesis of mtDNA-encoded proteins, S. pombe Mtf2 primarily functions in cox1 translation and the effect of mtf2 deletion on splicing of introns in mtDNA is likely due to a deficiency in the synthesis of intron-encoded maturases.

Thyroid transcriptome analysis reveals different adaptive responses to cold environmental conditions between two chicken breeds

Selection for cold tolerance in chickens is important for improving production performance and animal welfare. The identification of chicken breeds with higher cold tolerance and production performance will help to target candidates for the selection. The thyroid gland plays important roles in thermal adaptation, and its function is influenced by breed differences and transcriptional plasticity, both of which remain largely unknown in the chicken thyroid transcriptome. In this study, we subjected Bashang Long-tail (BS) and Rhode Island Red (RIR) chickens to either cold or warm environments for 21 weeks and investigated egg production performance, body weight changes, serum thyroid hormone concentrations, and thyroid gland transcriptome profiles. RIR chickens had higher egg production than BS chickens under warm conditions, but BS chickens produced more eggs than RIRs under cold conditions. Furthermore, BS chickens showed stable body weight gain under cold conditions while RIRs did not. These results suggested that BS breed is a preferable candidate for cold-tolerance selection and that the cold adaptability of RIRs should be improved in the future. BS chickens had higher serum thyroid hormone concentrations than RIRs under both environments. RNA-Seq generated 344.3 million paired-end reads from 16 sequencing libraries, and about 90% of the processed reads were concordantly mapped to the chicken reference genome. Differential expression analysis identified 46-1,211 genes in the respective comparisons. With regard to breed differences in the thyroid transcriptome, BS chickens showed higher cell replication and development, and immune response-related activity, while RIR chickens showed higher carbohydrate and protein metabolism activity. The cold environment reduced breed differences in the thyroid transcriptome compared with the warm environment. Transcriptional plasticity analysis revealed different adaptive responses in BS and RIR chickens to cope with the cold, and showed higher responsiveness in BS compared with RIR chickens, suggesting greater adaptability of the thyroid in BS chickens. Moreover, 10,053 differential splicing events were revealed among the groups, with RNA splicing and processing, gene expression, transport, and metabolism being the main affected biological processes, identifying a valuable alternative splicing repertoire for the chicken thyroid. A short isoform of TPO (encoding thyroid peroxidase) containing multiple open reading frames was generated in both breeds by skipping exons 4 and 5 in the cold environment. These findings provide novel clues for future studies of the molecular mechanisms underlying cold adaptation and/or acclimation in chickens.

Regulatory Potential of the RNA Processing Machinery: Implications for Human Disease

Splicing and nuclear export of mRNA are critical steps in the gene expression pathway. While RNA processing factors can perform general, essential functions for intron removal and bulk export of mRNA, emerging evidence highlights that the core RNA splicing and export machineries also display regulatory potential. Here, we discuss recent insights into how this regulatory potential can selectively alter gene expression and regulate important biological processes. We also highlight the participation of RNA processing pathways in the cellular response to DNA damage at multiple levels. These findings have important implications for the contribution of selective mRNA processing and export to the development of human cancers and neurodegenerative disorders.

Advantages of RNA-seq compared to RNA microarrays for transcriptome profiling of anterior cruciate ligament tears

Microarrays and RNA-seq are at the forefront of high throughput transcriptome analyses. Since these methodologies are based on different principles, there are concerns about the concordance of data between the two techniques. The concordance of RNA-seq and microarrays for genome-wide analysis of differential gene expression has not been rigorously assessed in clinically derived ligament tissues. To demonstrate the concordance between RNA-seq and microarrays and to assess potential benefits of RNA-seq over microarrays, we assessed differences in transcript expression in anterior cruciate ligament (ACL) tissues based on time-from-injury. ACL remnants were collected from patients with an ACL tear at the time of ACL reconstruction. RNA prepared from torn ACL remnants was subjected to Agilent microarrays (N = 24) and RNA-seq (N = 8). The correlation of biological replicates in RNA-seq and microarrays data was similar (0.98 vs. 0.97), demonstrating that each platform has high internal reproducibility. Correlations between the RNA-seq data and the individual microarrays were low, but correlations between the RNA-seq values and the geometric mean of the microarrays values were moderate. The cross-platform concordance for differentially expressed transcripts or enriched pathways was linearly correlated (r = 0.64). RNA-Seq was superior in detecting low abundance transcripts and differentiating biologically critical isoforms. Additional independent validation of transcript expression was undertaken using microfluidic PCR for selected genes. PCR data showed 100% concordance (in expression pattern) with RNA-seq and microarrays data. These findings demonstrate that RNA-seq has advantages over microarrays for transcriptome profiling of ligament tissues when available and affordable. Furthermore, these findings are likely transferable to other musculoskeletal tissues where tissue collection is challenging and cells are in low abundance. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:484-497, 2018.

Comparative evaluation of isoform-level gene expression estimation algorithms for RNA-seq and exon-array platforms

Given that the majority of multi-exon genes generate diverse functional products, it is important to evaluate expression at the isoform level. Previous studies have demonstrated strong gene-level correlations between RNA sequencing (RNA-seq) and microarray platforms, but have not studied their concordance at the isoform level. We performed transcript abundance estimation on raw RNA-seq and exon-array expression profiles available for common glioblastoma multiforme samples from The Cancer Genome Atlas using different analysis pipelines, and compared both the isoform- and gene-level expression estimates between programs and platforms. The results showed better concordance between RNA-seq/exon-array and reverse transcription-quantitative polymerase chain reaction (RT-qPCR) platforms for fold change estimates than for raw abundance estimates, suggesting that fold change normalization against a control is an important step for integrating expression data across platforms. Based on RT-qPCR validations, eXpress and Multi-Mapping Bayesian Gene eXpression (MMBGX) programs achieved the best performance for RNA-seq and exon-array platforms, respectively, for deriving the isoform-level fold change values. While eXpress achieved the highest correlation with the RT-qPCR and exon-array (MMBGX) results overall, RSEM was more highly correlated with MMBGX for the subset of transcripts that are highly variable across the samples. eXpress appears to be most successful in discriminating lowly expressed transcripts, but IsoformEx and RSEM correlate more strongly with MMBGX for highly expressed transcripts. The results also reinforce how potentially important isoform-level expression changes can be masked by gene-level estimates, and demonstrate that exon arrays yield comparable results to RNA-seq for evaluating isoform-level expression changes.

Impact of Alternative Splicing on the Human Proteome

Alternative splicing is a critical determinant of genome complexity and, by implication, is assumed to engender proteomic diversity. This notion has not been experimentally tested in a targeted, quantitative manner. Here, we have developed an integrative approach to ask whether perturbations in mRNA splicing patterns alter the composition of the proteome. We integrate RNA sequencing (RNA-seq) (to comprehensively report intron retention, differential transcript usage, and gene expression) with a data-independent acquisition (DIA) method, SWATH-MS (sequential window acquisition of all theoretical spectra-mass spectrometry), to capture an unbiased, quantitative snapshot of the impact of constitutive and alternative splicing events on the proteome. Whereas intron retention is accompanied by decreased protein abundance, alterations in differential transcript usage and gene expression alter protein abundance proportionate to transcript levels. Our findings illustrate how RNA splicing links isoform expression in the human transcriptome with proteomic diversity and provides a foundation for studying perturbations associated with human diseases.

The pre-mRNA retention and splicing complex controls expression of the Mediator subunit Med20

The heterotrimeric pre-mRNA retention and splicing (RES) complex, consisting of Bud13p, Snu17p and Pml1p, promotes splicing and nuclear retention of a subset of intron-containing pre-mRNAs. Yeast cells deleted for individual RES genes show growth defects that are exacerbated at elevated temperatures. Although the growth phenotypes correlate to the splicing defects in the individual mutants, the underlying mechanism is unknown. Here, we show that the temperature sensitive (Ts) growth phenotype of bud13Δ and snu17Δ cells is a consequence of inefficient splicing of MED20 pre-mRNA, which codes for a subunit of the Mediator complex; a co-regulator of RNA polymerase II transcription. The MED20 pre-mRNA splicing defect is less pronounced in pml1Δ cells, explaining why they grow better than the other 2 RES mutants at elevated temperatures. Inactivation of the cytoplasmic nonsense-mediated mRNA decay (NMD) pathway in the RES mutants leads to accumulation of MED20 pre-mRNA, indicating that inefficient nuclear retention contributes to the growth defect. Further, the Ts phenotype of bud13Δ and snu17Δ cells is partially suppressed by the inactivation of NMD, showing that the growth defects are augmented by the presence of a functional NMD pathway. Collectively, our results demonstrate an important role of the RES complex in maintaining the Med20p levels.

Ubiquitination of hnRNPA1 by TRAF6 links chronic innate immune signaling with myelodysplasia

Toll-like receptor (TLR) activation contributes to premalignant hematologic conditions, such as myelodysplastic syndromes (MDS). TRAF6, a TLR-effector with ubiquitin (Ub) ligase activity, is overexpressed in MDS hematopoietic stem/progenitor cells (HSPC). Here we show that TRAF6 overexpression in mouse HSPC resulted in impaired hematopoiesis and bone marrow failure. Through the use of a global Ub screen, we identified hnRNPA1, an RNA-binding protein and auxiliary splicing factor, as a substrate of TRAF6. TRAF6 ubiquitination of hnRNPA1 regulated alternative splicing of Arhgap1, which resulted in Cdc42 activation and accounted for hematopoietic defects in TRAF6-expressing HSPC. These results implicate Ub signaling in coordinating RNA processing by TLR pathways during an immune response and in premalignant hematologic diseases, such as MDS.

Workflow for Genome-Wide Determination of Pre-mRNA Splicing Efficiency from Yeast RNA-seq Data

Pre-mRNA splicing represents an important regulatory layer of eukaryotic gene expression. In the simple budding yeast Saccharomyces cerevisiae, about one-third of all mRNA molecules undergo splicing, and splicing efficiency is tightly regulated, for example, during meiotic differentiation. S. cerevisiae features a streamlined, evolutionarily highly conserved splicing machinery and serves as a favourite model for studies of various aspects of splicing. RNA-seq represents a robust, versatile, and affordable technique for transcriptome interrogation, which can also be used to study splicing efficiency. However, convenient bioinformatics tools for the analysis of splicing efficiency from yeast RNA-seq data are lacking. We present a complete workflow for the calculation of genome-wide splicing efficiency in S. cerevisiae using strand-specific RNA-seq data. Our pipeline takes sequencing reads in the FASTQ format and provides splicing efficiency values for the 5′ and 3′ splice junctions of each intron. The pipeline is based on up-to-date open-source software tools and requires very limited input from the user. We provide all relevant scripts in a ready-to-use form. We demonstrate the functionality of the workflow using RNA-seq datasets from three spliceosome mutants. The workflow should prove useful for studies of yeast splicing mutants or of regulated splicing, for example, under specific growth conditions.

Advances in analyzing RNA diversity in eukaryotic transcriptomes: peering through the Omics lens

Alternative splicing, polyadenylation, and chemical modifications of RNA generate astonishing complexity within eukaryotic transcriptomes. The last decade has brought numerous advances in sequencing technologies that allow biologists to investigate these phenomena with greater depth and accuracy while reducing time and cost. A commensurate development in biochemical techniques for the enrichment and analysis of different RNA variants has accompanied the advancement of global sequencing analysis platforms. Here, we present a detailed overview of the latest biochemical methods, along with bioinformatics pipelines that have aided in identifying different RNA variants. We also highlight the ongoing developments and challenges associated with RNA variant detection and quantification, including sample heterogeneity and isolation, as well as 'Omics' big data handling.

The power of fission: yeast as a tool for understanding complex splicing

Pre-mRNA splicing is an essential component of eukaryotic gene expression. Many metazoans, including humans, regulate alternative splicing patterns to generate expansions of their proteome from a limited number of genes. Importantly, a considerable fraction of human disease causing mutations manifest themselves through altering the sequences that shape the splicing patterns of genes. Thus, understanding the mechanistic bases of this complex pathway will be an essential component of combating these diseases. Dating almost to the initial discovery of splicing, researchers have taken advantage of the genetic tractability of budding yeast to identify the components and decipher the mechanisms of splicing. However, budding yeast lacks the complex splicing machinery and alternative splicing patterns most relevant to humans. More recently, many researchers have turned their efforts to study the fission yeast, Schizosaccharomyces pombe, which has retained many features of complex splicing, including degenerate splice site sequences, the usage of exonic splicing enhancers, and SR proteins. Here, we review recent work using fission yeast genetics to examine pre-mRNA splicing, highlighting its promise for modeling the complex splicing seen in higher eukaryotes.

Transcriptomic comparison of Drosophila snRNP biogenesis mutants reveals mutant-specific changes in pre-mRNA processing: implications for spinal muscular atrophy

Survival motor neuron (SMN) functions in the assembly of spliceosomal small nuclear ribonucleoproteins (snRNPs) that catalyze pre-mRNA splicing. Here, we used disruptions in Smn and two additional snRNP biogenesis genes, Phax and Ars2, to classify RNA processing differences as snRNP-dependent or gene-specific in Drosophila Phax and Smn mutants exhibited comparable reductions in snRNAs, and comparison of their transcriptomes uncovered shared sets of RNA processing changes. In contrast, Ars2 mutants displayed only small decreases in snRNA levels, and RNA processing changes in these mutants were generally distinct from those identified in Phax and Smn animals. Instead, RNA processing changes in Ars2 mutants support the known interaction of Ars2 protein with the cap-binding complex, as splicing changes showed a clear bias toward the first intron. Bypassing disruptions in snRNP biogenesis, direct knockdown of spliceosomal proteins caused similar changes in the splicing of snRNP-dependent events. However, these snRNP-dependent events were largely unaltered in three Smn mutants expressing missense mutations that were originally identified in human spinal muscular atrophy (SMA) patients. Hence, findings here clarify the contributions of Phax, Smn, and Ars2 to snRNP biogenesis in Drosophila, and loss-of-function mutants for these proteins reveal differences that help disentangle cause and effect in SMA model flies.