Number Of Alternative Splicing Events Research Articles

The fungicide pyraclostrobin affects gene expression by altering the DNA methylation pattern in Magnaporthe oryzae.

Rice blast disease caused by Magnaporthe oryzae has long been the main cause of rice (Oryza sativa L.) yield reduction worldwide. The quinone external inhibitor pyraclostrobin is widely used as a fungicide to effectively control the spread of pathogenic fungi, including M. oryzae. However, M. oryzae can develop resistance through multiple levels of mutation, such as target protein cytb mutation G143A/S, leading to a decrease in the effectiveness of the biocide after a period of application. Therefore, uncovering the possible mutational mechanisms from multiple perspectives will further provide feasible targets for drug development. In this work, we determined the gene expression changes in M. oryzae in response to pyraclostrobin stress and their relationship with DNA methylation by transcriptome and methylome. The results showed that under pyraclostrobin treatment, endoplasmic reticulum (ER)-associated and ubiquitin-mediated proteolysis were enhanced, suggesting that more aberrant proteins may be generated that need to be cleared. DNA replication and repair processes were inhibited. Glutathione metabolism was enhanced, while lipid metabolism was impaired. The number of alternative splicing events increased. These changes may be related to the elevated methylation levels of cytosine and adenine in gene bodies. Both hypermethylation and hypomethylation of differentially methylated genes (DMGs) mainly occurred in exons and promoters. Some DMGs and differentially expressed genes (DEGs) were annotated to the same pathways by GO and KEGG, including protein processing in the ER, ubiquitin-mediated proteolysis, RNA transport and glutathione metabolism, suggesting that pyraclostrobin may affect gene expression by altering the methylation patterns of cytosine and adenine. Our results revealed that 5mC and 6mA in the gene body are associated with gene expression and contribute to adversity adaptation in M. oryzae. This enriched the understanding for potential mechanism of quinone inhibitor resistance, which will facilitate the development of feasible strategies for maintaining the high efficacy of this kind of fungicide.

Ubiquitin Ligase SCF-FBXO11 Controls a Network of RNA-Binding Proteins and Splicing in MDS

The F-Box protein FBXO11 is a substrate receptor of Skp-Cullin-F-Box (SCF) ubiquitin ligase complexes. Loss-of-function mutations in FBXO11 are associated with diffuse large B-cell lymphoma, and FBXO11 has been shown to regulate epigenetic drivers of erythroid cell maturation. In our prior study, we found that deletion of FBXO11 improved cell survival under cytokine deprivation in an MDS cell line. Furthermore, we had observed a striking decrease in FBXO11 protein levels in AML patient samples. However, how FBXO11 contributes to myeloid malignancies is unclear. Here, using RNA-sequencing, proteomics, and CRISPR-Cas9 screening, we reveal broad ubiquitylation of RNA-binding proteins (RBPs) and alternative splicing as a mechanism by which deregulation of FBXO11 enhances malignant cell function. In addition, we report the presence of mutations in the uncharacterized N-terminal extension of FBXO11 in both lymphoid and myeloid malignancies. We first observed that heterozygous deletion of FBXO11 increases colony-forming ability of MDS-AML cell lines in vitro whilea2-fold increase in expression results in decreased cell growth and apoptosis. In a transplant model with a retrovirally-driven Runx1 mutation, heterozygous deletion of Fbxo11 exacerbated neutropenia and shortened overall survival. However, depletion of FBXO11 in healthy human and murine HSPCs significantly impaired myeloid progenitor colony-formation, indicating a unique requirement for FBXO11 in myelopoiesis. Based on its tumor suppressive activity, we evaluated whole-exome sequencing data from patient samples of both lymphoid and myeloid malignancies for FBXO11 mutations. This analysis uncovered rare mutations in FBXO11 in myeloid malignancies and revealed a substantial cluster of previously unidentified mutations in the Q-rich, N-terminal extension of FBXO11. We found that this uncharacterized, extended protein isoform is predominantly expressed in both normal and malignant hematopoietic cells. In colony-forming assays, expression of the long isoform and several of the N-terminal mutations inhibited CD34+ progenitor function. On a cellular level, lack of this N-terminal extension in the short isoform resulted in aberrant subcellular localization of FBXO11 within the nucleus. To elucidate the mechanism of FBXO11 in myeloid malignancy, we performed immunoprecipitation of FBXO11 followed by proteomics to identify endogenous FBXO11 complexes in MDS-AML cells. We integrated these data with ubiquitin proteomics of WT and FBXO11-KO MDS cells to identify candidate substrates of SCF-FBXO11-mediated ubiquitylation. FBXO11 complexes are significantly enriched for RNA-binding proteins affecting mRNA splicing, a process commonly deregulated in MDS. FBXO11-KO MDS cells displayed broad depletion of ubiquitylation in a network of RBPs. As expected, targeting FBXO11 did not affect the level of transcripts of interacting proteins; however, there were variable changes in total protein levels indicative of post-transcriptional regulation. We then applied a focused CRISPR-Cas9-KO screen of candidate SCF-FBXO11 substrates in MDS-AML cells depleted of FBXO11. We observed significant bi-directional effects in colony-forming ability mediated by several of the FBXO11-regulated RBPs including TRIM28, HNRNPU, NPM1, and SYNCRIP. Though FBXO11 is significantly decreased in MDS compared to healthy controls, we saw no difference in FBXO11 expression between splicing factor WT and mutant samples, suggesting that the function of FBXO11 in regulating mRNA splicing in MDS is independent of these mutations. Compared to healthy controls, FBXO11-low MDS samples exhibited a greater number of alternative splicing events, predominantly in skipped or mutually-exclusive exons. To experimentally test whether FBXO11 indeed affects RNA splicing, we utilized a bichromatic splicing reporter assay in MDS-AML cells which showed that depletion of FBXO11 is sufficient to alter splicing of the reading frame in the cassette exon, read out as a shift from GFP to RFP. Finally, STRING protein network analysis revealed that the alternative splicing events occurring in FBXO11-low MDS patient samples affect key nodes in protein translation and metabolism. Collectively, our data support a model whereby the change in FBXO11-mediated ubiquitylation of RBPs drives alternative exon usage that exacerbates malignant cell function in MDS.

Comprehensive Profiling of Alternative Splicing and Alternative Polyadenylation during Fruit Ripening in Watermelon (Citrullus lanatus).

Fruit ripening is a highly complicated process that is accompanied by the formation of fruit quality. In recent years, a series of studies have demonstrated post-transcriptional control play important roles in fruit ripening and fruit quality formation. Till now, the post-transcriptional mechanisms for watermelon fruit ripening have not been comprehensively studied. In this study, we conducted PacBio single-molecule long-read sequencing to identify genome-wide alternative splicing (AS), alternative polyadenylation (APA) and long non-coding RNAs (lncRNAs) in watermelon fruit. In total, 6,921,295 error-corrected and mapped full-length non-chimeric (FLNC) reads were obtained. Notably, more than 42,285 distinct splicing isoforms were derived from 5,891,183 intron-containing full-length FLNC reads, including a large number of AS events associated with fruit ripening. In addition, we characterized 21,506 polyadenylation sites from 11,611 genes, 8703 of which have APA sites. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that fructose and mannose metabolism, starch and sucrose metabolism and carotenoid biosynthesis were both enriched in genes undergoing AS and APA. These results suggest that post-transcriptional regulation might potentially have a key role in regulation of fruit ripening in watermelon. Taken together, our comprehensive PacBio long-read sequencing results offer a valuable resource for watermelon research, and provide new insights into the molecular mechanisms underlying the complex regulatory networks of watermelon fruit ripening.

In-depth analysis of alternative splicing landscape in multiple myeloma and potential role of dysregulated splicing factors

Splicing changes are common in cancer and are associated with dysregulated splicing factors. Here, we analyzed RNA-seq data from 323 newly diagnosed multiple myeloma (MM) patients and described the alternative splicing (AS) landscape. We observed a large number of splicing pattern changes in MM cells compared to normal plasma cells (NPC). The most common events were alterations of mutually exclusive exons and exon skipping. Most of these events were observed in the absence of overall changes in gene expression and often impacted the coding potential of the alternatively spliced genes. To understand the molecular mechanisms driving frequent aberrant AS, we investigated 115 splicing factors (SFs) and associated them with the AS events in MM. We observed that ~40% of SFs were dysregulated in MM cells compared to NPC and found a significant enrichment of SRSF1, SRSF9, and PCB1 binding motifs around AS events. Importantly, SRSF1 overexpression was linked with shorter survival in two independent MM datasets and was correlated with the number of AS events, impacting tumor cell proliferation. Together with the observation that MM cells are vulnerable to splicing inhibition, our results may lay the foundation for developing new therapeutic strategies for MM. We have developed a web portal that allows custom alternative splicing event queries by using gene symbols and visualizes AS events in MM and subgroups. Our portals can be accessed at http://rconnect.dfci.harvard.edu/mmsplicing/ and https://rconnect.dfci.harvard.edu/mmleafcutter/.

Comprehensive analysis of pre-mRNA alternative splicing regulated by m6A methylation in pig oxidative and glycolytic skeletal muscles

BackgroundDifferent types of skeletal myofibers exhibit distinct physiological and metabolic properties that are associated with meat quality traits in livestock. Alternative splicing (AS) of pre-mRNA can generate multiple transcripts from an individual gene by differential selection of splice sites. N6-methyladenosine (m6A) is the most abundant modification in mRNAs, but its regulation for AS in different muscles remains unknown. ResultsWe characterized AS events and m6A methylation pattern in pig oxidative and glycolytic muscles. A tota1 of 1294 differential AS events were identified, and differentially spliced genes were significantly enriched in processes related to different phenotypes between oxidative and glycolytic muscles. We constructed the regulatory network between splicing factors and corresponding differential AS events and identified NOVA1 and KHDRBS2 as key splicing factors. AS event was enriched in m6A-modified genes, and the methylation level was positively correlated with the number of AS events in genes. The dynamic change in m6A enrichment was associated with 115 differentially skipping exon (SE-DAS) events within 92 genes involving in various processes, including muscle contraction and myofibril assembly. We obtained 23.4% SE-DAS events (27/115) regulated by METTL3-meditaed m6A and experimentally validated the aberrant splicing of ZNF280D, PHE4DIP, and NEB. The inhibition of m6A methyltransferase METTL3 could induce the conversion of oxidative fiber to glycolytic fiber in PSCs.ConclusionOur study suggested that m6A modification could contribute to significant difference in phenotypes between oxidative and glycolytic muscles by mediating the regulation of AS. These findings would provide novel insights into mechanisms underlying muscle fiber conversion.

The Isolation and Full-Length Transcriptome Sequencing of a Novel Nidovirus and Response of Its Infection in Japanese Flounder (Paralichthys olivaceus).

A novel nidovirus, CSBV Bces-Po19, was isolated from the marine fish, Japanese flounder (Paralichthys olivaceus). The viral genome was 26,597 nucleotides long and shared 98.62% nucleotide identity with CSBV WHQSR4345. PacBio Sequel and Illumina sequencing were used to perform full-length transcriptome sequencing on CSBV Bces-Po19-sensitive (S) and -resistant (R) Japanese flounder. The results of negative staining revealed bacilliform and spherical virions. There were in total 1444 different genes between CSBV Bces-Po19 S and R groups, with 935 being up-regulated and 513 being down-regulated. Metabolism-, immune-, and RNA-related pathways were significantly enriched. Furthermore, CSBV Bces-Po19 infection induced alternative splicing (AS) events in Japanese flounder; the S group had a higher numbers of AS events (12,352) than the R group (11,452). The number of long non-coding RNA (lncRNA) in the S group, on the other hand, was significantly lower than in the R group. In addition to providing valuable information that sheds more light on CSBV Bces-Po19 infection, these research findings provide further clues for CSBV Bces-Po19 prevention and treatment.

The influence of 4-thiouridine labeling on pre-mRNA splicing outcomes.

Metabolic labeling is a widely used tool to investigate different aspects of pre-mRNA splicing and RNA turnover. The labeling technology takes advantage of native cellular machineries where a nucleotide analog is readily taken up and incorporated into nascent RNA. One such analog is 4-thiouridine (4sU). Previous studies demonstrated that the uptake of 4sU at elevated concentrations (>50μM) and extended exposure led to inhibition of rRNA synthesis and processing, presumably induced by changes in RNA secondary structure. Thus, it is possible that 4sU incorporation may also interfere with splicing efficiency. To test this hypothesis, we carried out splicing analyses of pre-mRNA substrates with varying levels of 4sU incorporation (0-100%). We demonstrate that increased incorporation of 4sU into pre-mRNAs decreased splicing efficiency. The overall impact of 4sU labeling on pre-mRNA splicing efficiency negatively correlates with the strength of splice site signals such as the 3' and the 5' splice sites. Introns with weaker splice sites are more affected by the presence of 4sU. We also show that transcription by T7 polymerase and pre-mRNA degradation kinetics were impacted at the highest levels of 4sU incorporation. Increased incorporation of 4sU caused elevated levels of abortive transcripts, and fully labeled pre-mRNA is more stable than its uridine-only counterpart. Cell culture experiments show that a small number of alternative splicing events were modestly, but statistically significantly influenced by metabolic labeling with 4sU at concentrations considered to be tolerable (40 μM). We conclude that at high 4sU incorporation rates small, but noticeable changes in pre-mRNA splicing can be detected when splice sites deviate from consensus. Given these potential 4sU artifacts, we suggest that appropriate controls for metabolic labeling experiments need to be included in future labeling experiments.

Genome-Wide Analysis of Alternative Splicing during Host-Virus Interactions in Chicken.

The chicken is a model animal for the study of evolution, immunity and development. In addition to their use as a model organism, chickens also represent an important agricultural product. Pathogen invasion has already been shown to modulate the expression of hundreds of genes, but the role of alternative splicing in avian virus infection remains unclear. We used RNA-seq data to analyze virus-induced changes in the alternative splicing of Gallus gallus, and found that a large number of alternative splicing events were induced by virus infection both in vivo and in vitro. Virus-responsive alternative splicing events preferentially occurred in genes involved in metabolism and transport. Many of the alternatively spliced transcripts were also expressed from genes with a function relating to splicing or immune response, suggesting a potential impact of virus infection on pre-mRNA splicing and immune gene regulation. Moreover, exon skipping was the most frequent AS event in chickens during virus infection. This is the first report describing a genome-wide analysis of alternative splicing in chicken and contributes to the genomic resources available for studying host–virus interaction in this species. Our analysis fills an important knowledge gap in understanding the extent of genome-wide alternative splicing dynamics occurring during avian virus infection and provides the impetus for the further exploration of AS in chicken defense signaling and homeostasis.

Genome-wide transcriptome analysis of the early developmental stages of Echinococcus granulosus protoscoleces reveals extensive alternative splicing events in the spliceosome pathway

BackgroundThe complex life cycle of Echinococcus granulosus involves numerous environmental conditions within different intermediate and definitive hosts. This requires adaptation at different levels of transcript regulation. Alternative splicing (AS) and the related cellular functions as one of the major fields of post-genomics has been poorly studied in tapeworms. In the present study, we investigated AS events and their potential biological effects in E. granulosus.MethodsWhole transcriptome sequencing data of four groups of protoscoleces were prepared for RNA-seq library construction. Fresh protoscoleces were either used as non-induced controls (NT group) or incubated for 15 min with pepsin (PEP group) and cultivated in a biphasic medium for 12 and 24 h (12 and 24 h groups). The frequency and different types of AS events were identified using rMATS software. Functional annotations and gene ontology of differential AS (DAS) genes were performed using Blast2GO software. AS events were experimentally validated by PCR on the protoscolex cDNAs using specific primers for each gene.ResultsAt least one AS event was found in 38.1% of the genes (3904 out of 10,245) in the protoscoleces during early strobilar development. The genes were associated primarily with cellular and metabolic processes and binding and catalytic activity. KEGG pathway analysis of DAS events revealed a number of genes belonging to different components of the spliceosome complex. These genes tended to belong to common SR proteins, U1-related factors, U2-related factors, complex A-specific factors and other splicing-related proteins.ConclusionsThe high number of AS events in the transcriptome regulatory mechanisms indicates the essential rapid molecular processes required by the parasite for adaptation in different environments.Graphical

The Alternative Splicing Landscape in Multiple Myeloma Is Determined By IGH Translocations and Mutations of RNA Processing Genes

Introduction: Alternative Splicing (AS) plays a key role in regulating numerous cellular processes in both normal and malignant cells. Previous studies have revealed mutations in the spliceosome complex, such as SF3B1, can cause increased AS frequencies in multiple myeloma (MM) patients, and patients with increased levels of AS are associated with a poor prognosis. Other frequently mutated genes involved in RNA processing include DIS3 and FAM46C, thus, systematically investigating other causes of AS abnormalities and pathologies in MM patients is highly necessary.Materials and Methods: RNA-seq data from 598 newly diagnosed MM patients from the MMRF CoMMpass study were utilized to generate AS comparisons. They were previously annotated for cytogenetic, copy number, and mutation data (version IA16). RNA-seq data were aligned to HG38 using STAR and Salmon. SUPPA2 was used for calling AS differences. For each identified AS event, the splicing level was defined by Percentage of Spliced-In (PSI) while the mean difference of splicing levels between two groups was measured by ΔPSI (dPSI) and by the P-value from independent T-tests against PSIs in the two groups. Filtering thresholds were determined to find high-quality differentially spliced events and were filtered for those also present in normal PCs (GSE110486). Geneset enrichment analysis was performed to identify dysregulated pathways caused by differential splicing and differential expression. Survival analysis was performed on clinical annotations of 598 NDMM patients while the Logrank test and Cox regression were used to evaluate the risk of AS and other genomic factors. Kaplan-Meier curves were plotted for various subgroups.Results: We compared 16 major cytogenetic subgroups, including Ig translocations (t(4;14), t(14;16), t(11;14)), hyperdiploidy, mutations in KRAS, NRAS, BRAF, FAM46C, SF3B1, DIS3 and TP53, combined events (t(4;14) plus DIS3 mutation), as well as those with biallelic abnormalities (DIS3, FAM46C, and TP53). Samples with SF3B1 hotspot mutations identified the greatest number of AS events (n=862), and samples with any SF3B1 mutation had approximately half as many. IGH translocations had an equivalent number of AS events to those with SF3B1 mutations, with t(14;16) having the most (n=587) followed by t(11;14) (n=366), and t(4;14) (n=256). We observed an increased number of significant AS events in bi-allelic DIS3 and FAM46C groups (n=404 and 171) compared to their mono-allelic abnormalities (n=114 and 35). As DIS3 mutations are enriched in the t(4;14) subgroup we also examined that interaction and found significantly more AS events (n=481; p<0.01) in the combination compared to either event alone. As expected, KRAS, NRAS and BRAF mutations did not have enrichment for AS events (n=2, 15, 23, respectively).The majority of AS events were unique to each subgroup, exemplifying the AS heterogeneity in these subgroups. Among overlapped events, an alternative first (AF) exon in ACACA was consistently more spliced in t(14;16), t(11;14) and t(4;14) groups (dPSI=0.18, 0.10, 0.12, P=2x10 -5, 2x10 -9, 5x10 -5). ACACA encodes an enzyme that significantly affects MM cell growth and viability, suggesting that similar regulations exist in the three translocation groups. Unique events were also detected including an AF event in MIB2 (E3 Ubiquitin Protein Ligase 2) in the t(11;14) group (dPSI=0.17, P=7x10 -14), and a skipped exon in UBXN4 (related to ER stress) in t(14;16) group (dPSI=0.1, P=3x10 -4). AS heterogeneity also leads to functional heterogeneity in the three groups. Besides commonly downregulated RNA catabolic processes, cell adhesion, migration and mobility related pathways are enriched pathways in t(14;16); cell growth related pathways in t(11;14); and ERK related pathways in t(4;14). High-risk events were identified through survival analysis and included a retained intron in RPS16 in the t(14;16) (Hazard Ratio (HR)=18.81, p=0.004). Similarly, high risk was associated with an AF event in DDX39B in t(11;14) (HR=2.62, p=0.001) and an AF event in COPA in t(4;14) (HR=6.29, p=0.001).Conclusion: AS is defined by multiple genomic events, including primary translocations and mutations in RNA processing genes, DIS3 and FAM46C, and interactions between genomic markers can increase AS. AS events contribute to outcome and some high risk AS events may serve as prognostic marker or potentially novel targets. DisclosuresWalker: Sanofi: Speakers Bureau; Bristol Myers Squibb: Research Funding.

Phosphoproteome analysis reveals the involvement of protein dephosphorylation in ethylene-induced corolla senescence in petunia

BackgroundSenescence represents the last stage of flower development. Phosphorylation is the key posttranslational modification that regulates protein functions, and kinases may be more required than phosphatases during plant growth and development. However, little is known about global phosphorylation changes during flower senescence.ResultsIn this work, we quantitatively investigated the petunia phosphoproteome following ethylene or air treatment. In total, 2170 phosphosites in 1184 protein groups were identified, among which 2059 sites in 1124 proteins were quantified. To our surprise, treatment with ethylene resulted in 697 downregulated and only 117 upregulated phosphosites using a 1.5-fold threshold (FDR < 0.05), which showed that ethylene negatively regulates global phosphorylation levels and that phosphorylation of many proteins was not necessary during flower senescence. Phosphoproteome analysis showed that ethylene regulates ethylene and ABA signalling transduction pathways via phosphorylation levels. One of the major targets of ethylene-induced dephosphorylation is the plant mRNA splicing machinery, and ethylene treatment increases the number of alternative splicing events of precursor RNAs in petunia corollas.ConclusionsProtein dephosphorylation could play an important role in ethylene-induced senescence, and ethylene treatment increased the number of AS precursor RNAs in petunia corollas.

Integrated SMRT Technology with UMI RNA-Seq Reveals the Hub Genes in Stamen Petalody in Camellia oleifera

Male sterility caused by stamen petalody is a key factor for a low fruit set rate and a low yield of Camellia oleifera but can serve as a useful genetic tool because it eliminates the need for artificial emasculation. However, its molecular regulation mechanism still remains unclear. In this study, transcriptome was sequenced and analyzed on two types of bud materials, stamen petalody mutants and normal materials, at six stages of stamen development based on integrated single-molecule real-time (SMRT) technology with unique molecular identifiers (UMI) and RNA-seq technology to identify the hub genes responsible for stamen petalody in C. oleifera. The results show that a large number of alternative splicing events were identified in the transcriptome. A co-expression network analysis of MADSs and all the differentially expressed genes between the mutant stamens and the normal materials showed that four MADS transcription factor genes, CoSEP3.1, CoAGL6, CoSEP3.2, and CoAP3, were predicted to be the hub genes responsible for stamen petalody. Among these four, the expression patterns of CoAGL6 and CoSEP3.2 were consistently high in the mutant samples, but relatively low in the normal samples at six stages, while the patterns of CoSEP3.1 and CoAP3 were initially low in mutants and then were upregulated during development but remained relatively high in the normal materials. Furthermore, the genes with high connectivity to the hub genes showed significantly different expression patterns between the mutant stamens and the normal materials at different stages. qRT-PCR results showed a similar expression pattern of the hub genes in the RNA-seq. These results lay a solid foundation for the directive breeding of C. oleifera varieties and provide references for the genetic breeding of ornamental Camellia varieties.

The Alter Retina: Alternative Splicing of Retinal Genes in Health and Disease.

Alternative splicing of mRNA is an essential mechanism to regulate and increase the diversity of the transcriptome and proteome. Alternative splicing frequently occurs in a tissue- or time-specific manner, contributing to differential gene expression between cell types during development. Neural tissues present extremely complex splicing programs and display the highest number of alternative splicing events. As an extension of the central nervous system, the retina constitutes an excellent system to illustrate the high diversity of neural transcripts. The retina expresses retinal specific splicing factors and produces a large number of alternative transcripts, including exclusive tissue-specific exons, which require an exquisite regulation. In fact, a current challenge in the genetic diagnosis of inherited retinal diseases stems from the lack of information regarding alternative splicing of retinal genes, as a considerable percentage of mutations alter splicing or the relative production of alternative transcripts. Modulation of alternative splicing in the retina is also instrumental in the design of novel therapeutic approaches for retinal dystrophies, since it enables precision medicine for specific mutations.

SMRT- and Illumina-based RNA-seq analyses unveil the ginsinoside biosynthesis and transcriptomic complexity in Panax notoginseng

Panax notoginseng is one of the most widely used traditional Chinese herbs with particularly valued roots. Triterpenoid saponins are mainly specialized secondary metabolites, which medically act as bioactive components. Knowledge of the ginsenoside biosynthesis in P. notoginseng, which is of great importance in the industrial biosynthesis and genetic breeding program, remains largely undetermined. Here we combined single molecular real time (SMRT) and Second-Generation Sequencing (SGS) technologies to generate a widespread transcriptome atlas of P. notoginseng. We mapped 2,383 full-length non-chimeric (FLNC) reads to adjacently annotated genes, corrected 1,925 mis-annotated genes and merged into 927 new genes. We identified 8,111 novel transcript isoforms that have improved the annotation of the current genome assembly, of which we found 2,664 novel lncRNAs. We characterized more alternative splicing (AS) events from SMRT reads (20,015 AS in 6,324 genes) than Illumina reads (18,498 AS in 9,550 genes), which contained a number of AS events associated with the ginsenoside biosynthesis. The comprehensive transcriptome landscape reveals that the ginsenoside biosynthesis predominantly occurs in flowers compared to leaves and roots, substantiated by levels of gene expression, which is supported by tissue-specific abundance of isoforms in flowers compared to roots and rhizomes. Comparative metabolic analyses further show that a total of 17 characteristic ginsenosides increasingly accumulated, and roots contained the most ginsenosides with variable contents, which are extraordinarily abundant in roots of the three-year old plants. We observed that roots were rich in protopanaxatriol- and protopanaxadiol-type saponins, whereas protopanaxadiol-type saponins predominated in aerial parts (leaves, stems and flowers). The obtained results will greatly enhance our understanding about the ginsenoside biosynthetic machinery in the genus Panax.

Genome-wide identification and expression analysis of the SR gene family in longan (Dimocarpus longan Lour.)

Longan (Dimocarpus longan Lour.) is an important commercial fruit tree in southern China. The embryogenesis of longan affects the quality and yield of fruit. A large number of alternative splicing events occurs during somatic embryogenesis (SE), which is regulated by serine/arginine-rich (SR) proteins. However, the functions of SR proteins in longan are poorly understood. In this study, 21 Dlo-SR gene family members belonging to six subfamilies were identified, among which Dlo-RSZ20a, Dlo-SR30, Dlo-SR17, Dlo-SR53 and Dlo-SR32 were localized in the nucleus, Dlo-RSZ20b, Dlo-RSZ20c, Dlo-RSZ20d, Dlo-SC18, Dlo-RS2Z29, Dlo-SCL41, and Dlo-SR33 were localized in chloroplasts, and Dlo-RS43, Dlo-SC33, Dlo-SC37, Dlo-RS2Z33, Dlo-RS2Z16, Dlo-RS2Z24, Dlo-SCL43, Dlo-SR112, and Dlo-SR59 were localized in the nucleus and chloroplasts. The Dlo-SR genes exhibited differential expression patterns in different tissues of longan. The transcript levels of Dlo-RSZ20a, Dlo-SC18, Dlo-RS2Z29, DLo-SR59, Dlo-SR53, and Dlo-SR17 were low in all analyzed tissues, whereas Dlo-RS43, Dlo-RS2Z16, Dlo-RS2Z24, and Dlo-SR30 were highly expressed in all tissues. To clarify their function during SE, the transcript levels of Dlo-SR genes were analyzed at different four stages of SE, comprising non-embryonic callus (NEC), friable-embryogenic callus (EC), incomplete compact pro-embryogenic culture (ICpEC) and globular embryo (GE). Interestingly, the transcript levels of Dlo-RS2Z29 and Dlo-SR112 were increased in embryogenic cells compared with the NEC stage, whereas transcript levels of Dlo-RSZ20a, Dlo-RS43, Dlo-SC37, and Dlo-RS2Z16 were especially increased at the GE stage compared with the other stages. Alternative splicing events of Dlo-SR mRNA precursors (pre-mRNAs) was detected during SE, with totals of 41, 29, 35, and 44 events detected during NEC, EC, ICpEC, and GE respectively. Protein–protein interaction analysis showed that SR proteins were capable of interaction with each other. The results indicate that the alternative splicing of Dlo-SR pre-mRNAs occurs during SE and that Dlo-SR proteins may interact to regulate embryogenesis of longan.

Wide-ranging analysis of survival-related alternative splicing events in invasive breast carcinoma.

Invasive breast carcinoma (BRCA) is a serious disease that threatens the survival time of those affected. Alternative splicing (AS) involved in BRCA pathogenesis may be a potential therapeutic target. However, to the best of our knowledge, a systematic analysis of survival-related alternative splicing events (SREs) has not yet been reported. The aim of the present study was to identify SREs and analyze their potential biological functions as BRCA prognostic biomarkers. An UpSet plot demonstrated AS global characteristics. Cox's proportional hazards regression model quantitatively demonstrated the prognostic relevance of AS events. Functional enrichment analysis investigated the potential pathways through which AS events affect BRCA progression. The receiver operating characteristic curve model determined the clinical significance of AS events represented using percent-spliced-in (PSI) values. The regulatory network of splicing factors (SFs) and AS events laid the foundation for studying the role of SFs in BRCA. The present study identified 1,215 SREs and their distribution characteristics, suggesting that AS events in exon skipping (ES) primarily exerted normal physiological functions, while AS events in alternative terminator sites had the most significant prognostic effect. The present study demonstrated that survival-associated genes are involved primarily in certain biological processes of ribosomal proteins. In the diagnostic model, the alternative acceptor site, alternative donor site, alternative promoter site and ES performed well. ELAVL4 was the key gene associated with prognosis and SREs. In conclusion, a number of AS events affect BRCA initiation, progression and prognosis. The PSI value of AS events has the potential to diagnose BRCA and predict a prognosis; however, this must be confirmed in additional studies.

Asian Zika Virus Isolate Significantly Changes the Transcriptional Profile and Alternative RNA Splicing Events in a Neuroblastoma Cell Line.

The alternative splicing of pre-mRNAs expands a single genetic blueprint to encode multiple, functionally diverse protein isoforms. Viruses have previously been shown to interact with, depend on, and alter host splicing machinery. The consequences, however, incited by viral infection on the global alternative slicing (AS) landscape are under-appreciated. Here, we investigated the transcriptional and alternative splicing profile of neuronal cells infected with a contemporary Puerto Rican Zika virus (ZIKVPR) isolate, an isolate of the prototypical Ugandan ZIKV (ZIKVMR), and dengue virus 2 (DENV2). Our analyses revealed that ZIKVPR induced significantly more differential changes in expressed genes compared to ZIKVMR or DENV2, despite all three viruses showing equivalent infectivity and viral RNA levels. Consistent with the transcriptional profile, ZIKVPR induced a higher number of alternative splicing events compared to ZIKVMR or DENV2, and gene ontology analyses highlighted alternative splicing changes in genes associated with mRNA splicing. In summary, we show that ZIKV affects cellular RNA homeostasis not only at the transcriptional levels but also through the alternative splicing of cellular transcripts. These findings could provide new molecular insights into the neuropathologies associated with this virus.

Endophytic fungus, Fusarium sp. reduces alternative splicing events in rice plants under salinity stress

Alternative splicing (AS) of transcripts is a well-recognized phenomenon in plants. Several studies have shown that the frequency of AS events are significantly enhanced in plants subjected to abiotic stresses, presumably as an adaptation to increase transcript diversity. In this study, we examine the effect of a salt tolerant endophytic fungus, Fusarium sp. on the AS events of IR-64 rice plants, under salinity stress. RNA seq data generated from rice plants, treated and not treated with Fusarium sp. respectively and subjected to 150 mM salinity stress was used to analyze the number and type of AS events and their chromosomal distribution using approriate bioinformatics pipeline. Besides, we also annotated the locus IDs of genes for gene enrichment and KEGG pathway analysis. Our results show an unequivocal decrease in the number of AS events under salinity stress as influenced by the endophyte. The average AS events per gene also decreased from 2.28 to 2.11 upon colonization by the endophyte. This is the first documented evidence of an endophyte-induced alteration in the frequency of alternative splicing event in plants subjected to salinity stress. The exact mechanisms through which the splicing activities are restrained, however need to be further investigated.

Coding and Non-Coding RNA Abnormalities in Bipolar Disorder.

The molecular mechanisms underlying bipolar disorder (BPD) have remained largely unknown. Postmortem brain tissue studies comparing BPD patients with healthy controls have produced a heterogeneous array of potentially implicated protein-coding RNAs. We hypothesized that dysregulation of not only coding, but multiple classes of RNA (coding RNA, long non-coding (lnc) RNA, circular (circ) RNA, and/or alternative splicing) underlie the pathogenesis of BPD. Using non-polyadenylated libraries we performed RNA sequencing in postmortem human medial frontal gyrus tissue from BPD patients and healthy controls. Twenty genes, some of which not previously implicated in BPD, were differentially expressed (DE). PCR validation and replication confirmed the implication of these DE genes. Functional in silico analyses identified enrichment of angiogenesis, vascular system development and histone H3-K4 demethylation. In addition, ten lncRNA transcripts were differentially expressed. Furthermore, an overall increased number of alternative splicing events in BPD was detected, as well as an increase in the number of genes carrying alternative splicing events. Finally, a large reservoir of circRNAs populating brain tissue not affected by BPD is described, while in BPD altered levels of two circular transcripts, cNEBL and cEPHA3, are reported. cEPHA3, hitherto unlinked to BPD, is implicated in developmental processes in the central nervous system. Although we did not perform replication analyses of non-coding RNA findings, our findings hint that RNA dysregulation in BPD is not limited to coding regions, opening avenues for future pharmacological investigations and biomarker research.

Genome-Wide Analysis of Alternative Splicing Provides Insights Into Stress Response of the Pacific White Shrimp Litopenaeus vanname.

Alternative splicing (AS) can enhance transcript diversity dramatically and play an important role in stress adaptation. Limited researches of AS have been reported in the Pacific white shrimp (Litopenaeus vannamei), which is an important aquaculture species in the world. Here, we performed a genome-wide identification of AS events in L. vannamei based on eight transcriptomes. We identified 38,781 AS events in the shrimp genome, and some of them were validated by polymerase chain reaction experiments. These AS events correspond to 9,209 genes, accounting for 36% of protein-coding genes in the shrimp genome. The number of AS events increased after virus or bacteria infection and low salinity stress. Type 1 AS genes (AS was initially activated) were mainly enriched in substance and energy metabolism, such as carbon metabolism and amino metabolism. However, type 2 AS genes (AS events changed) displayed specific enrichment under different stress challenges. Specifically, type 2 AS genes under biotic stresses were mainly enriched in the pathogenic pathway and immune network, and the AS genes under low salinity stress were significantly enriched for betalain biosynthesis. In summary, our study indicates that AS events are complex in shrimp and may be related to stress adaptation. These results will provide valuable resource for functional genomic studies on crustaceans.