Genome-wide analysis of alternative splicing in Caenorhabditis elegans

Alternative splicing patterns of hnrnp genes in gill tissues of rainbow trout (Oncorhynchus mykiss) during salinity changes

Alternative splicing (AS) plays an important role in various physiological processes in eukaryotes, such as the stress response. However, patterns of AS events remain largely unexplored during salinity acclimation in fishes. In this study, we conducted AS analysis using RNA-seq datasets to explore splicing patterns in the gill tissues of rainbow trout exposed to altered salinity environments, ranging from 0 ‰ (T0) to 30 ‰ (T30). The results revealed 1441, 351, 483, 1051 and 1049 differentially alternatively spliced (DAS) events in 5 pairwise comparisons, including T6 vs. T0, T12 vs. T0, T18 vs. T0, T24 vs. T0, and T30 vs. T0, respectively. These DAS events were derived from 1290, 328, 444, 963 and 948 genes. Enrichment analysis indicated that these DAS genes were related to RNA splicing and processing. Among these, 14 DAS genes were identified as members of the large heterogeneous nuclear RNP (hnRNP) gene family. Alternative 3' splice site (A3SS), exon skipping (SE) and intron retention (RI) events resulted in the fragmentation or even loss of the functional RNA recognition motif (RRM) domains in hnrnpa0, hnrnp1a, hnrnp1b and hnrnpc genes. The incomplete RRM domains would hinder the interactions between hnRNP genes and pre-mRNAs. It would in turn influence the splicing patterns and mRNA stability of downstream target genes in response to salinity changes. The study provides insights into salinity acclimation in gill tissues of rainbow trout and serves as a significant reference on the osmoregulation mechanisms at post-transcription regulation levels in fish.

CCAR-1 works together with the U2AF large subunit UAF-1 to regulate alternative splicing

The Cell Division Cycle and Apoptosis Regulator (CCAR) protein family members have recently emerged as regulators of alternative splicing and transcription, as well as having other key physiological functions. For example, mammalian CCAR2/DBC1 forms a complex with the zinc factor protein ZNF326 to integrate alternative splicing with RNA polymerase II transcriptional elongation in AT-rich regions of the DNA. Additionally, Caenorhabditis elegans CCAR-1, a homolog to mammalian CCAR2, facilitates the alternative splicing of the perlecan unc-52 gene. However, much about the CCAR family's role in alternative splicing is unknown. Here, we have examined the role of CCAR-1 in genome-wide alternative splicing in Caenorhabditis elegans and have identified new alternative splicing targets of CCAR-1 using RNA sequencing. Also, we found that CCAR-1 interacts with the spliceosome factors UAF-1 and UAF-2 using mass spectrometry, and that knockdown of ccar-1 affects alternative splicing patterns, motility, and proteostasis of UAF-1 mutant worms. Collectively, we demonstrate the role of CCAR-1 in regulating global alternative splicing in C. elegans and in conjunction with UAF-1.

Analysis of alternative splicing in chicken macrophages infected with avian pathogenic E. coli (APEC)

Colibacillosis is a complex disease that caused by avian pathogenic Escherichia coli (APEC), resulting in huge economic loss to the global poultry industry and threatening to human health. Alternative splicing (AS) is a universal post-transcriptional regulatory mechanism, which can simultaneously produce many proteins from a single gene to involve in various diseases and individual development. Herein, we characterized genome-wide AS events in wild type macrophages (WT) and APEC infected macrophages (APEC) by high-throughput RNA sequencing technology. A total of 751 differentially expressed (DE) AS genes were identified in the comparison of APEC vs. WT, including 587 of SE, 114 of MXE, 25 of RI, 17 of A3 and 8 of A5 event. Functional analysis showed that these identified DE AS genes were involved in 'Endocytosis', 'p53 signaling pathway', 'MAPK signaling pathway', 'NOD-like receptor signaling pathway', 'Ubiquitin mediated proteolysis' and 'Focal adhesion' immune related pathways. In summary, we comprehensively investigate AS events during APEC infection. This study has expanded our understanding of the process of APEC infection and provided new insights for further treatment options for APEC infection.

Alternative splicing in Atlantic salmon head kidney and SHK-1 cell line during the Piscirickettsia salmonis infection: A comparative transcriptome survey

Piscirickettsia salmonis, an intracellular bacterium in salmon aquaculture, is a big challenge because it is responsible for 54.2% of Atlantic salmon mortalities. In recent years, the high relevance of Alternative Splicing (AS) as a molecular mechanism associated with infectious conditions and host-pathogen interaction processes, especially in host immune activation, has been observed. Several studies have highlighted the role of AS in the host's immune response during viral, bacterial, and endoparasite infection. In the present study, we evaluated AS transcriptome profiles during P. salmonis infection in the two most used study models, SHK-1 cell line and salmon head kidney tissue. First, the SHK-1 cell line was exposed to P. salmonis infection at 0-, 7-, and 14-days post-infection (dpi). Following, total RNA was extracted for Illumina sequencing. On the other hand, RNA-Seq datasets of Atlantic salmon head kidney infected with the same P. salmonis strayingwase used. For both study models, the highest number of differentially alternative splicing (DAS) events was observed at 7 dpi, 16,830 DAS events derived from 9213 DAS genes in SHK-1 cells, and 13,820 DAS events from 7684 DAS genes in salmon HK. Alternative first exon (AF) was the most abundant AS type in the three infection times analyzed, representing 31% in SHK-1 cells and 228.6 in salmon HK; meanwhile, mutually exclusive exon (MX) was the least abundant. Notably, functional annotation of DAS genes in SHK-1 cells infected with P. salmonis showed a high presence of genes related to nucleotide metabolism. In contrast, the salmon head kidney exhibited many GO terms associated with immune response. Our findings reported the role of AS during P. salmonis infection in Atlantic salmon. These studies would contribute to a better understanding of the molecular bases that support the pathogen-host interaction, evidencing the contribution of AS regulating the transcriptional host response.

The Responses of Alternative Splicing during Heat Stress in the Pacific White Shrimp Litopenaeus vannamei

Heat tolerance is increasingly becoming a crucial trait for aquaculture species in the face of rapidly changing climate conditions. Alternative splicing (AS) is a vital mechanism within cells that modulates gene abundance and functional diversity, enabling organisms to effectively respond to diverse stressful conditions, including thermal stress. However, it is still uncertain whether AS contributes to heat tolerance in shrimp. In this study, we conducted an extensive transcriptome analysis on the Pacific white shrimp, Litopenaeus vannamei, revealing a total of 1267, 987, and 130 differential AS events (DAS) in the gill, hepatopancreas, and muscle, respectively, following exposure to heat stress. Among all of the DAS events, exon skipping (ES) was the predominant form of splicing modification observed. Interestingly, a minor portion of DAS genes exhibited overlap across the three tissues, implying that heat stress exerts unique effects on various tissue types. Moreover, the functional enrichment analysis demonstrated that commonly identified DAS genes were primarily associated with the "spliceosome" pathway, indicating that the AS of splicing-related genes played a crucial role in the response to heat stress. Our findings also revealed that heat stress tended to induce longer mRNA isoforms through differential alternative 3' splice site (A3SS) events. Notably, A3SS events exhibited the highest proportion of maintained open reading frames (ORFs) under heat stress. Interestingly, we observed a limited overlap between the genes exhibiting DAS and those showing differential gene expression (DEG), indicating that AS may function as a distinct regulatory mechanism independent of transcriptional regulation in response to heat stress. This is the first comprehensive study on AS in crustacea species under heat stress, which broadens our understanding of the regulatory mechanisms governing the crustaceans' response to environmental stress, providing valuable insights for the aquaculture breeding of shrimp and other aquatic animals.

Single Nematode Transcriptomic Analysis, Using Long-Read Technology, Reveals Two Novel Virulence Gene Candidates in the Soybean Cyst Nematode, Heterodera glycines

The soybean cyst nematode (Heterodera glycines, SCN), is the most damaging disease of soybean in North America. While management of this pest using resistant soybean is generally still effective, prolonged exposure to cultivars derived from the same source of resistance (PI 88788) has led to the emergence of virulence. Currently, the underlying mechanisms responsible for resistance breakdown remain unknown. In this study, we combined a single nematode transcriptomic profiling approach with long-read sequencing to reannotate the SCN genome. This resulted in the annotation of 1932 novel transcripts and 281 novel gene features. Using a transcript-level quantification approach, we identified eight novel effector candidates overexpressed in PI 88788 virulent nematodes in the late infection stage. Among these were the novel gene Hg-CPZ-1 and a pioneer effector transcript generated through the alternative splicing of the non-effector gene Hetgly21698. While our results demonstrate that alternative splicing in effectors does occur, we found limited evidence of direct involvement in the breakdown of resistance. However, our analysis highlighted a distinct pattern of effector upregulation in response to PI 88788 resistance indicative of a possible adaptation process by SCN to host resistance.

On the origin and nature of nongenetic information in eumetazoans

Nongenetic information implies all the forms of biological information not related to genes and DNA in general. Despite the deep scientific relevance of the concept, we currently lack reliable knowledge about its carriers and origins; hence, we still do not understand its true nature. Given that genes are the targets of nongenetic information, it appears that a parsimonious approach to find the ultimate source of that information is to trace back the sequential steps of the causal chain upstream of the target genes up to the ultimate link as the source of the nongenetic information. From this perspective, I examine seven nongenetically determined phenomena: placement of locus-specific epigenetic marks on DNA and histones, changes in snRNA expression patterns, neural induction of gene expression, site-specific alternative gene splicing, predator-induced morphological changes, and cultural inheritance. Based on the available evidence, I propose a general model of the common neural origin of all these forms of nongenetic information in eumetazoans.

Identification of sex-specific splicing via comparative transcriptome analysis of gonads from sea cucumber Apostichopus japonicus

Alternative splicing (AS) is an essential post-transcriptional regulation mechanism for sex differentiation and gonadal development, which has rarely been reported in marine invertebrates. Sea cucumber (Apostichopus japonicus) is an economically important marine benthic echinoderm with a potential XX/XY sex determination mechanism, whose molecular mechanism in the gonadal differentiation has not been clearly understood. In this study, we analyzed available RNA-seq datasets of male and female gonads to explore if AS mechanism exerts an essential function in sex differentiation and gonadal development of A. japonicus. In our results, a total of 20,338 AS events from 7219 alternatively spliced genes, and 189 sexually differential alternative splicing (DAS) events from 156 genes were identified in gonadal transcriptome of sea cucumber. Gene Ontology analysis indicated that these DAS genes were significantly enriched in spermatogenesis-related GO terms. Maximal Clique Centrality (MCC) was then applied for protein-protein interaction (PPI) analysis to search for protein interactions and hub DAS gene. Among all DAS genes, we identified 10 DAS genes closely related to spermatogenesis and (or) sperm motility and a hub gene dnah1. Thus, this study revealed that alternative isoforms were generated from certain genes in female and male gonads through alternative splicing, which may provide direct evidence that alternative splicing mechanisms participate in female and male gonads. These results suggested a novel perspective for explaining the molecular mechanisms underlying gonadal differentiation between male and female sea cucumbers.

Re-evaluating the impact of alternative RNA splicing on proteomic diversity

Alternative splicing (AS) constitutes a mechanism by which protein-coding genes and long non-coding RNA (lncRNA) genes produce more than a single mature transcript. From plants to humans, AS is a powerful process that increases transcriptome complexity. Importantly, splice variants produced from AS can potentially encode for distinct protein isoforms which can lose or gain specific domains and, hence, differ in their functional properties. Advances in proteomics have shown that the proteome is indeed diverse due to the presence of numerous protein isoforms. For the past decades, with the help of advanced high-throughput technologies, numerous alternatively spliced transcripts have been identified. However, the low detection rate of protein isoforms in proteomic studies raised debatable questions on whether AS contributes to proteomic diversity and on how many AS events are really functional. We propose here to assess and discuss the impact of AS on proteomic complexity in the light of the technological progress, updated genome annotation, and current scientific knowledge.

Transcriptome Analysis Reveals the Involvement of Alternative Splicing in the Nitrogen Starvation Response of Chlamydomonas reinhardtii

Alternative splicing (AS) is a regulatory mechanism of post-transcriptional regulation that plays an important role in plant response to abiotic stresses. However, corresponding research involving the mechanism of AS in the nitrogen starvation response of C. reinhardtii is rare. This study performed a comprehensive and systematic analysis of AS events in C. reinhardtii at nine time points (0 h, 10 m, 30 m, 1 h, 6 h, 8 h, 24 h, and 48 h) under nitrogen starvation. It used STAR and rMATS tools to identify and quantify the probability of the AS event happening through the transcriptome high-throughput sequencing data. A total of 5806 AS events in 3500 genes were identified, and the retained intron and skipped exon were considered the main AS types. The genes related to the AS event in nitrogen starvation were mainly involved in spliceosome and transporter and enriched in the citrate cycle and fatty acid degradation pathways. These results suggested that AS may play an important role in the nitrogen starvation response in C. reinhardtii, and provided insights into post-transcriptional regulation under nitrogen starvation.

Eomes function is conserved between zebrafish and mouse and controls left-right organiser progenitor gene expression via interlocking feedforward loops

The T-box family transcription factor Eomesodermin (Eomes) is present in all vertebrates, with many key roles in the developing mammalian embryo and immune system. Homozygous Eomes mutant mouse embryos exhibit early lethality due to defects in both the embryonic mesendoderm and the extraembryonic trophoblast cell lineage. In contrast, zebrafish lacking the predominant Eomes homologue A (Eomesa) do not suffer complete lethality and can be maintained. This suggests fundamental differences in either the molecular function of Eomes orthologues or the molecular configuration of processes in which they participate. To explore these hypotheses we initially analysed the expression of distinct Eomes isoforms in various mouse cell types. Next we compared the functional capabilities of these murine isoforms to zebrafish Eomesa. These experiments provided no evidence for functional divergence. Next we examined the functions of zebrafish Eomesa and other T-box family members expressed in early development, as well as its paralogue Eomesb. Though Eomes is a member of the Tbr1 subfamily we found evidence for functional redundancy with the Tbx6 subfamily member Tbx16, known to be absent from eutherians. However, Tbx16 does not appear to synergise with Eomesa cofactors Mixl1 and Gata5. Finally, we analysed the ability of Eomesa and other T-box factors to induce zebrafish left-right organiser progenitors (known as dorsal forerunner cells) known to be positively regulated by vgll4l, a gene we had previously shown to be repressed by Eomesa. Here we demonstrate that Eomesa indirectly upregulates vgll4l expression via interlocking feedforward loops, suggesting a role in establishment of left-right asymmetry. Conversely, other T-box factors could not similarly induce left-right organiser progenitors. Overall these findings demonstrate conservation of Eomes molecular function and participation in similar processes, but differential requirements across evolution due to additional co-expressed T-box factors in teleosts, albeit with markedly different molecular capabilities. Our analyses also provide insights into the role of Eomesa in left-right organiser formation in zebrafish.

Genome annotation of Caenorhabditis briggsae by TEC-RED identifies new exons, paralogs, and conserved and novel operons

The nematode Caenorhabditis briggsae is routinely used in comparative and evolutionary studies involving its well-known cousin Caenorhabditis elegans. The C. briggsae genome sequence has accelerated research by facilitating the generation of new resources, tools, and functional studies of genes. While substantial progress has been made in predicting genes and start sites, experimental evidence is still lacking in many cases. Here, we report an improved annotation of the C. briggsae genome using the trans-spliced exon coupled RNA end determination technique. In addition to identifying the 5' ends of expressed genes, we have discovered operons and paralogs. In summary, our analysis yielded 10,243 unique 5' end sequence tags with matches in the C. briggsae genome. Of these, 6,395 were found to represent 4,252 unique genes along with 362 paralogs and 52 previously unknown exons. These genes included 14 that are exclusively trans-spliced in C. briggsae when compared with C. elegans orthologs. A major contribution of this study is the identification of 492 high confidence operons, of which two-thirds are fully supported by tags. In addition, 2 SL1-type operons were discovered. Interestingly, comparisons with C. elegans showed that only 40% of operons are conserved. Of the remaining operons, 73 are novel, including 12 that entirely lack orthologs in C. elegans. Further analysis revealed that 4 of the 12 novel operons are conserved in Caenorhabditis nigoni. Altogether, the work described here has significantly advanced our understanding of the C. briggsae system and serves as a rich resource to aid biological studies involving this species.

Intron retention coupled with nonsense-mediated decay is involved in cellulase biosynthesis in cellulolytic fungi

BACKGROUND

Knowledge on regulatory networks associated with cellulase biosynthesis is prerequisite for exploitation of such regulatory systems in enhancing cellulase production with low cost. The biological functions of intron retention (IR) and nonsense-mediated mRNA decay (NMD) in filamentous fungi is lack of study, let alone their roles in cellulase biosynthesis.

RESULTS

We found that major cellulase genes (cel7a, cel7b, and cel3a) exhibited concomitant decrease in IR rates and increase in their gene expression in T. reesei under cellulase-producing condition (cellulose and lactose) that was accompanied with a more active NMD pathway, as compared to cellulase non-producing condition (glucose). In the presence of the NMD pathway inhibitor that successfully repressed the NMD pathway, the mRNA levels of cellulase genes were sharply down-regulated, but the rates of IR in these genes were significantly up-regulated. Consistently, the cellulase activities were severely inhibited. In addition, the NMD pathway inhibitor caused the downregulated mRNA levels of two important genes of the target of rapamycin (TOR) pathway, trfkbp12 and trTOR1. The absence of gene trfkbp12 made the cellulase production in T. reesei more sensitive to the NMD pathway inhibitor.

CONCLUSIONS

All these findings suggest that the IR of cellulase genes regulates their own gene expression by coupling with the NMD pathway, which might involve the TOR pathway. Our results provide better understanding on intron retention, the NMD pathway, and cellulase production mechanism in filamentous fungi.

Two-Color Fluorescent Reporters for Analysis of Alternative Splicing

Alternative splicing is a key layer of gene regulation that is frequently modulated in a spatiotemporal manner. As such, it is a major goal to understand the mechanisms controlling alternative splicing in specific cellular contexts. Reporters that recapitulate alternative splicing patterns of endogenous transcripts have served as excellent tools for dissecting regulatory mechanisms of splicing. In this chapter, we describe a two-color fluorescent reporter system that enables the visualization of alternative splicing patterns by microscopy at single-cell resolution in live animals. We present this reporter system in the context of the model nematode C. elegans.

Comparative transcriptome profiles of Schistosoma japonicum larval stages: Implications for parasite biology and host invasion

Schistosoma japonicum is prevalent in Asia with a wide mammalian host range, which leads to highly harmful zoonotic parasitic diseases. Most previous transcriptomic studies have been performed on this parasite, but mainly focus on stages inside the mammalian host. Moreover, few larval transcriptomic data are available in public databases. Here we mapped the detailed transcriptome profiles of four S. japonicum larval stages including eggs, miracidia, sporocysts and cercariae, providing a comprehensive development picture outside of the mammalian host. By analyzing the stage-specific/enriched genes, we identified functional genes associated with the biological characteristic at each stage: e.g. we observed enrichment of genes necessary for DNA replication only in sporocysts, while those involved in proteolysis were upregulated in sporocysts and/or cercariae. This data indicated that miracidia might use leishmanolysin and neprilysin to penetrate the snail, while elastase (SjCE2b) and leishmanolysin might contribute to the cercariae invasion. The expression profile of stem cell markers revealed potential germinal cell conversion during larval development. Additionally, our analysis indicated that tandem duplications had driven the expansion of the papain family in S. japonicum. Notably, all the duplicated cathepsin B-like proteases were highly expressed in cercariae. Utilizing our 3rd version of S. japonicum genome, we further characterized the alternative splicing profiles throughout these four stages. Taken together, the present study provides compressive gene expression profiles of S. japonicum larval stages and identifies a set of genes that might be involved in intermediate and definitive host invasion.

Alternative splicing level related to intron size and organism complexity

Background

Alternative splicing is the process of selecting different combinations of splice sites to produce variably spliced mRNAs. However, the relationships between alternative splicing prevalence and level (ASP/L) and variations of intron size and organism complexity (OC) remain vague. Here, we developed a robust protocol to analyze the relationships between ASP/L and variations of intron size and OC. Approximately 8 Tb raw RNA-Seq data from 37 eumetazoan species were divided into three sets of species based on variations in intron size and OC.

Results

We found a strong positive correlation between ASP/L and OC, but no correlation between ASP/L and intron size across species. Surprisingly, ASP/L displayed a positive correlation with mean intron size of genes within individual genomes. Moreover, our results revealed that four ASP/L-related pathways contributed to the differences in ASP/L that were associated with OC. In particular, the spliceosome pathway displayed distinct genomic features, such as the highest gene expression level, conservation level, and fraction of disordered regions. Interestingly, lower or no obvious correlations were observed among these genomic features.

Conclusions

The positive correlation between ASP/L and OC ubiquitously exists in eukaryotes, and this correlation is not affected by the mean intron size of these species. ASP/L-related splicing factors may play an important role in the evolution of OC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08172-2.

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

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

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

BACKGROUND

The 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.

METHODS

Whole 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.

RESULTS

At 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.

CONCLUSIONS

The 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.

m6 A-mediated alternative splicing coupled with nonsense-mediated mRNA decay regulates SAM synthetase homeostasis

Alternative splicing of pre-mRNAs can regulate gene expression levels by coupling with nonsense-mediated mRNA decay (NMD). In order to elucidate a repertoire of mRNAs regulated by alternative splicing coupled with NMD (AS-NMD) in an organism, we performed long-read RNA sequencing of poly(A)⁺ RNAs from an NMD-deficient mutant strain of Caenorhabditis elegans, and obtained full-length sequences for mRNA isoforms from 259 high-confidence AS-NMD genes. Among them are the S-adenosyl-L-methionine (SAM) synthetase (sams) genes sams-3 and sams-4. SAM synthetase activity autoregulates sams gene expression through AS-NMD in a negative feedback loop. We furthermore find that METT-10, the orthologue of human U6 snRNA methyltransferase METTL16, is required for the splicing regulation in␣vivo, and specifically methylates the invariant AG dinucleotide at the distal 3' splice site (3'SS) in␣vitro. Direct RNA sequencing coupled with machine learning confirms m⁶ A modification of endogenous sams mRNAs. Overall, these results indicate that homeostasis of SAM synthetase in C. elegans is maintained by alternative splicing regulation through m⁶ A modification at the 3'SS of the sams genes.

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

Simple Summary

Studies have shown that alternative splicing (AS) has been utilized in a wide variety of life processes. However, there are very few studies on AS during germ cell development. In this study, we preliminarily investigated the variation of variable shear events during the formation of chicken germ cells through the RNA-seq data analysis of embryonic stem cells (ESCs), gonad PGCs (gPGCs), and spermatogonia stem cells (SSCs), and the critical AS mode for several crucial stage-specific genes, which were identified during germ cell development. The results of this study lay a theoretical foundation for further analysis of the regulation mechanism of key genes involved in germ cell formation.

Abstract

Alternative splicing (AS) is a ubiquitous, co-transcriptional, and post-transcriptional regulation mechanism during certain developmental processes, such as germ cell differentiation. A thorough understanding of germ cell differentiation will help us to open new avenues for avian reproduction, stem cell biology, and advances in medicines for human consumption. Here, based on single-cell RNA-seq, we characterized genome-wide AS events in manifold chicken male germ cells: embryonic stem cells (ESCs), gonad primordial germ cells (gPGCs), and spermatogonia stem cells (SSCs). A total of 38,494 AS events from 15,338 genes were detected in ESCs, with a total of 48,955 events from 14,783 genes and 49,900 events from 15,089 genes observed in gPGCs and SSCs, respectively. Moreover, this distribution of AS events suggests the diverse splicing feature of ESCs, gPGCs, and SSCs. Finally, several crucial stage-specific genes, such as NANOG, POU5F3, LIN28B, BMP4, STRA8, and LHX9, were identified in AS events that were transmitted in ESCs, gPGCs, and SSCs. The gene expression results of the RNA-seq data were validated by qRT-PCR. In summary, we provided a comprehensive atlas of the genome-wide scale of the AS event landscape in male chicken germ-line cells and presented its distribution for the first time. This research may someday improve treatment options for men suffering from male infertility.

Comparative Genome-Wide Alternative Splicing Analysis of Longissimus Dorsi Muscles Between Japanese Black (Wagyu) and Chinese Red Steppes Cattle

Alternative splicing is a ubiquitous regulatory mechanism in gene expression that allows a single gene generating multiple messenger RNAs (mRNAs). Significant differences in fat deposition ability and meat quality traits have been reported between Japanese black cattle (Wagyu) and Chinese Red Steppes, which presented a unique model for analyzing the effects of transcriptional level on marbling fat in livestock. In previous studies, the differentially expressed genes (DGEs) in longissimus dorsi muscle (LDM) samples between Wagyu and other breeds of beef cattle have been reported. In this study, we further investigated the differences in alternative splicing in LDM between Wagyu and Chinese Red Steppes cattle. We identified several alternative splicing types including cassette exon, mutually exclusive exons, alternative 5′ splice site, alternative 3′ splice site, alternative start exon, and intron retention. In total, 115 differentially expressed alternatively spliced genes were obtained, of which 17 genes were enriched in the metabolic pathway. Among the 17 genes, 5 genes, including MCAT, CPT1B, HADHB, SIRT2, and DGAT1, appeared to be the novel spliced candidates that affect the lipid metabolism in cattle. Additionally, another 17 genes were enriched in the Gene Ontology (GO) terms related to muscle development, such as NR4A1, UQCC2, YBX3/CSDA, ITGA7, etc. Overall, altered splicing and expression levels of these novel candidates between Japanese black cattle and Chinese Red Steppes revealed by RNA-seq suggest their potential involvement in the muscle development and fat deposition of beef cattle.

Nanopore long-read RNAseq reveals regulatory mechanisms of thermally variable reef environments promoting heat tolerance of scleractinian coral Pocillopora damicornis

Some scleractinian corals exhibit high thermal adaptability to climate changes, although the mechanism of their adaptation is unclear. This study investigated the adaptability of scleractinian coral Pocillopora damicornis to thermally variable reef environments by applying a nanopore-based RNA sequencing method to characterize different transcription responses that promote heat tolerance of P. damicornis. We identified 1414 novel genes and optimized 6256 mis-annotated loci. Based on full-length transcriptome data, we identified complex alternative polyadenylation and alternative splicing events, which can improve our understanding of the genome annotation and gene structures of P. damicornis. Furthermore, we constructed differentially expressed lncRNA-mRNA co-expression networks, which may play a crucial role in the P. damicornis thermal adaptive response. KEGG function enrichment analysis revealed that P. damicornis from the high-temperature pool had a lower metabolic rate than that from the low-temperature pool. We hypothesize that metabolic readjustment, in the form of a lower metabolic rate, positively correlated with increased heat tolerance in P. damicornis in thermally variable reef environments. Our study provides novel insights into lncRNAs that promote thermally tolerance of scleractinian corals in the thermally variable reef environment, suggesting potential mechanisms for their adaptation to global warming in the future.

Loss of ESRP1 blocks mouse oocyte development and leads to female infertility

Alternative splicing (AS) contributes to gene diversification, but the AS program during germline development remains largely undefined. Here, we interrupted pre-mRNA splicing events controlled by epithelial splicing regulatory protein 1 (ESRP1) and found that it induced female infertility in mice. Esrp1 deletion perturbed spindle organization, chromosome alignment and metaphase-to-anaphase transformation in oocytes. The first polar body extrusion was blocked during oocyte meiosis owing to abnormal activation of spindle assembly checkpoint and insufficiency of anaphase-promoting complex/cyclosome in Esrp1-knockout oocytes. Esrp1-knockout hampered follicular development and ovulation; eventually, premature ovarian failure occurred in six-month-old Esrp1-knockout mouse. Using single-cell RNA-seq analysis, 528 aberrant AS events of maternal mRNA transcripts were revealed and were preferentially associated with microtubule cytoskeletal organization. Notably, we found that loss of ESRP1 disturbed a comprehensive set of gene-splicing sites - including those within Trb53bp1, Rac1, Bora, Kif2c, Kif23, Ndel1, Kif3a, Cenpa and Lsm14b - that potentially caused abnormal spindle organization. Collectively, our findings provide the first report elucidating the ESRP1-mediated AS program of maternal mRNA transcripts, which may contribute to oocyte meiosis and female fertility in mice.

Alternative splicing in endothelial cells: novel therapeutic opportunities in cancer angiogenesis

Alternative splicing (AS) is a pervasive molecular process generating multiple protein isoforms, from a single gene. It plays fundamental roles during development, differentiation and maintenance of tissue homeostasis, while aberrant AS is considered a hallmark of multiple diseases, including cancer. Cancer-restricted AS isoforms represent either predictive biomarkers for diagnosis/prognosis or targets for anti-cancer therapies. Here, we discuss the contribution of AS regulation in cancer angiogenesis, a complex process supporting disease development and progression. We consider AS programs acting in a specific and non-redundant manner to influence morphological and functional changes involved in cancer angiogenesis. In particular, we describe relevant AS variants or splicing regulators controlling either secreted or membrane-bound angiogenic factors, which may represent attractive targets for therapeutic interventions in human cancer.

Global regulatory features of alternative splicing across tissues and within the nervous system of C. elegans

Alternative splicing plays a major role in shaping tissue-specific transcriptomes. Among the broad tissue types present in metazoans, the central nervous system contains some of the highest levels of alternative splicing. Although many documented examples of splicing differences between broad tissue types exist, there remains much to be understood about the splicing factors and the cis sequence elements controlling tissue and neuron subtype-specific splicing patterns. By using translating ribosome affinity purification coupled with deep-sequencing (TRAP-seq) in Caenorhabditis elegans, we have obtained high coverage profiles of ribosome-associated mRNA for three broad tissue classes (nervous system, muscle, and intestine) and two neuronal subtypes (dopaminergic and serotonergic neurons). We have identified hundreds of splice junctions that exhibit distinct splicing patterns between tissue types or within the nervous system. Alternative splicing events differentially regulated between tissues are more often frame-preserving, are more highly conserved across Caenorhabditis species, and are enriched in specific cis regulatory motifs, when compared with other types of exons. By using this information, we have identified a likely mechanism of splicing repression by the RNA-binding protein UNC-75/CELF via interactions with cis elements that overlap a 5′ splice site. Alternatively spliced exons also overlap more frequently with intrinsically disordered peptide regions than constitutive exons. Moreover, regulated exons are often shorter than constitutive exons but are flanked by longer intron sequences. Among these tissue-regulated exons are several highly conserved microexons <27 nt in length. Collectively, our results indicate a rich layer of tissue-specific gene regulation at the level of alternative splicing in C. elegans that parallels the evolutionary forces and constraints observed across metazoa.

Alternative Splicing During the Chlamydomonas reinhardtii Cell Cycle

Genome-wide analysis of transcriptome data in Chlamydomonas reinhardtii shows periodic patterns in gene expression levels when cultures are grown under alternating light and dark cycles so that G1 of the cell cycle occurs in the light phase and S/M/G0 occurs during the dark phase. However, alternative splicing, a process that enables a greater protein diversity from a limited set of genes, remains largely unexplored by previous transcriptome based studies in C. reinhardtii. In this study, we used existing longitudinal RNA-seq data obtained during the light-dark cycle to investigate the changes in the alternative splicing pattern and found that 3277 genes (19.75% of 17,746 genes) undergo alternative splicing. These splicing events include Alternative 5′ (Alt 5′), Alternative 3′ (Alt 3′) and Exon skipping (ES) events that are referred as alternative site selection (ASS) events and Intron retention (IR) events. By clustering analysis, we identified a subset of events (26 ASS events and 10 IR events) that show periodic changes in the splicing pattern during the cell cycle. About two-thirds of these 36 genes either introduce a pre-termination codon (PTC) or introduce insertions or deletions into functional domains of the proteins, which implicate splicing in altering gene function. These findings suggest that alternative splicing is also regulated during the Chlamydomonas cell cycle, although not as extensively as changes in gene expression. The longitudinal changes in the alternative splicing pattern during the cell cycle captured by this study provides an important resource to investigate alternative splicing in genes of interest during the cell cycle in Chlamydomonas reinhardtii and other eukaryotes.

Alternative splicing of coq-2 controls the levels of rhodoquinone in animals

Parasitic helminths use two benzoquinones as electron carriers in the electron transport chain. In normoxia, they use ubiquinone (UQ), but in anaerobic conditions inside the host, they require rhodoquinone (RQ) and greatly increase RQ levels. We previously showed the switch from UQ to RQ synthesis is driven by a change of substrates by the polyprenyltransferase COQ-2 (Del Borrello et al., 2019; Roberts Buceta et al., 2019); however, the mechanism of substrate selection is not known. Here, we show helminths synthesize two coq-2 splice forms, coq-2a and coq-2e, and the coq-2e-specific exon is only found in species that synthesize RQ. We show that in Caenorhabditis elegans COQ-2e is required for efficient RQ synthesis and survival in cyanide. Importantly, parasites switch from COQ-2a to COQ-2e as they transit into anaerobic environments. We conclude helminths switch from UQ to RQ synthesis principally via changes in the alternative splicing of coq-2.

mRNA Editing, Processing and Quality Control in Caenorhabditis elegans

While DNA serves as the blueprint of life, the distinct functions of each cell are determined by the dynamic expression of genes from the static genome. The amount and specific sequences of RNAs expressed in a given cell involves a number of regulated processes including RNA synthesis (transcription), processing, splicing, modification, polyadenylation, stability, translation, and degradation. As errors during mRNA production can create gene products that are deleterious to the organism, quality control mechanisms exist to survey and remove errors in mRNA expression and processing. Here, we will provide an overview of mRNA processing and quality control mechanisms that occur in Caenorhabditis elegans, with a focus on those that occur on protein-coding genes after transcription initiation. In addition, we will describe the genetic and technical approaches that have allowed studies in C. elegans to reveal important mechanistic insight into these processes.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Direct full-length RNA sequencing reveals unexpected transcriptome complexity during Caenorhabditis elegans development

Massively parallel sequencing of the polyadenylated RNAs has played a key role in delineating transcriptome complexity, including alternative use of an exon, promoter, 5′ or 3′ splice site or polyadenylation site, and RNA modification. However, reads derived from the current RNA-seq technologies are usually short and deprived of information on modification, compromising their potential in defining transcriptome complexity. Here, we applied a direct RNA sequencing method with ultralong reads using Oxford Nanopore Technologies to study the transcriptome complexity in Caenorhabditis elegans. We generated approximately six million reads using native poly(A)-tailed mRNAs from three developmental stages, with average read lengths ranging from 900 to 1100 nt. Around half of the reads represent full-length transcripts. To utilize the full-length transcripts in defining transcriptome complexity, we devised a method to classify the long reads as the same as existing transcripts or as a novel transcript using sequence mapping tracks rather than existing intron/exon structures, which allowed us to identify roughly 57,000 novel isoforms and recover at least 26,000 out of the 33,500 existing isoforms. The sets of genes with differential expression versus differential isoform usage over development are largely different, implying a fine-tuned regulation at isoform level. We also observed an unexpected increase in putative RNA modification in all bases in the coding region relative to the UTR, suggesting their possible roles in translation. The RNA reads and the method for read classification are expected to deliver new insights into RNA processing and modification and their underlying biology in the future.

Genome-wide evolution of wobble base-pairing nucleotides of branchpoint motifs with increasing organismal complexity

How have the branchpoint motifs evolved in organisms of different complexity? Here we identified and examined the consensus motifs (R₁C₂T₃R₄A₅Y₆, R: A or G, Y: C or T) of 898 fungal genomes. In Ascomycota unicellular yeasts, the G₄/A₄ ratio is mostly (98%) below 0.125 but increases sharply in multicellular species by about 40 times on average, and in the more complex Basidiomycota, it increases further by about 7 times. The global G₄ increase is consistent with A₄ to G₄ transitions in evolution. Of the G₄/A₄-interacting amino acids of the branchpoint binding protein MSL5 (SF1) and the HSH155 (SF3B1), as well as the 5' splice sites (SS) and U2 snRNA genes, the 5' SS G³/A³ co-vary with the G₄ to some extent. However, corresponding increase of the G₄-complementary GCAGTA-U2 gene is rare, suggesting wobble-base pairing between the G₄-containing branchpoint motif and GTAGTA-U2 in most of these species. Interestingly, the G₄/A₄ ratio correlates well with the abundance of alternative splicing in the two phyla, and G₄ enriched significantly at the alternative 3' SS of genes in RNA metabolism, kinases and membrane proteins. Similar wobble nucleotides also enriched at the 3' SS of multicellular fungi with only thousands of protein-coding genes. Thus, branchpoint motifs have evolved U2-complementarity in unicellular Ascomycota yeasts, but have gradually gained more wobble base-pairing nucleotides in fungi of higher complexity, likely to destabilize branchpoint motif-U2 interaction and/or branchpoint A protrusion for alternative splicing. This implies an important role of relaxing the branchpoint signals in the multicellularity and further complexity of fungi.

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.

Combining next-generation sequencing and single-molecule sequencing to explore brown plant hopper responses to contrasting genotypes of japonica rice

BACKGROUND

The brown plant hopper (BPH), Nilaparvata lugens, is one of the major pest of rice (Oryza sativa). Plant defenses against insect herbivores have been extensively studied, but our understanding of insect responses to host plants' resistance mechanisms is still limited. The purpose of this study is to characterize transcripts of BPH and reveal the responses of BPH insects to resistant rice at transcription level by using the advanced molecular techniques, the next-generation sequencing (NGS) and the single-molecule, real-time (SMRT) sequencing.

RESULTS

The current study obtained 24,891 collapsed isoforms of full-length transcripts, and 20,662 were mapped to known annotated genes, including 17,175 novel transcripts. The current study also identified 915 fusion genes, 1794 novel genes, 2435 long non-coding RNAs (lncRNAs), and 20,356 alternative splicing events. Moreover, analysis of differentially expressed genes (DEGs) revealed that genes involved in metabolic and cell proliferation processes were significantly enriched in up-regulated and down-regulated sets, respectively, in BPH fed on resistant rice relative to BPH fed on susceptible wild type rice. Furthermore, the FoxO signaling pathway was involved and genes related to BPH starvation response (Nlbmm), apoptosis and autophagy (caspase 8, ATG13, BNIP3 and IAP), active oxygen elimination (catalase, MSR, ferritin) and detoxification (GST, CarE) were up-regulated in BPH responses to resistant rice.

CONCLUSIONS

The current study provides the first demonstrations of the full diversity and complexity of the BPH transcriptome, and indicates that BPH responses to rice resistance, might be related to starvation stress responses, nutrient transformation, oxidative decomposition, and detoxification. The current result findings will facilitate further exploration of molecular mechanisms of interaction between BPH insects and host rice.

The regulation properties of RNA secondary structure in alternative splicing

The RNA secondary structure is important for many functional processes in the cell. The secondary and tertiary structures of cellular RNAs are essential for the activity of these molecules in processes such as transcription, splicing, translation, and localization. New high-throughput analytical methods, including next generation sequencing, have allowed for the in-depth characterization of the 'RNA structurome': a new term describing how the RNA structure controls the activity of RNA by itself and how it regulates the expression of genes. In this review, we present many examples of the influence of structural motifs of RNA, long range interactions and global RNA structure on the alternative splicing processes. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

Conserved functions of RNA-binding proteins in muscle

Animals require different types of muscle for survival, for example for circulation, motility, reproduction and digestion. Much emphasis in the muscle field has been placed on understanding how transcriptional regulation generates diverse types of muscle during development. Recent work indicates that alternative splicing and RNA regulation are as critical to muscle development, and altered function of RNA-binding proteins causes muscle disease. Although hundreds of genes predicted to bind RNA are expressed in muscles, many fewer have been functionally characterized. We present a cross-species view summarizing what is known about RNA-binding protein function in muscle, from worms and flies to zebrafish, mice and humans. In particular, we focus on alternative splicing regulated by the CELF, MBNL and RBFOX families of proteins. We discuss the systemic nature of diseases associated with loss of RNA-binding proteins in muscle, focusing on mis-regulation of CELF and MBNL in myotonic dystrophy. These examples illustrate the conservation of RNA-binding protein function and the marked utility of genetic model systems in understanding mechanisms of RNA regulation.

Predominant patterns of splicing evolution on human, chimpanzee and macaque evolutionary lineages

Although splicing is widespread and evolves rapidly among species, the mechanisms driving this evolution, as well as its functional implications, are not yet fully understood. We analyzed the evolution of splicing patterns based on transcriptome data from five tissues of humans, chimpanzees, rhesus macaques and mice. In total, 1526 exons and exon sets from 1236 genes showed significant splicing differences among primates. More than 60% of these differences represent constitutive-to-alternative exon transitions while an additional 25% represent changes in exon inclusion frequency. These two dominant evolutionary patterns have contrasting conservation, regulation and functional features. The sum of these features indicates that, despite their prevalence, constitutive-to-alternative exon transitions do not substantially contribute to long-term functional transcriptome changes. Conversely, changes in exon inclusion frequency appear to be functionally relevant, especially for changes taking place in the brain on the human evolutionary lineage.

Young genes have distinct gene structure, epigenetic profiles, and transcriptional regulation

Species-specific, new, or "orphan" genes account for 10%-30% of eukaryotic genomes. Although initially considered to have limited function, an increasing number of orphan genes have been shown to provide important phenotypic innovation. How new genes acquire regulatory sequences for proper temporal and spatial expression is unknown. Orphan gene regulation may rely in part on origination in open chromatin adjacent to preexisting promoters, although this has not yet been assessed by genome-wide analysis of chromatin states. Here, we combine taxon-rich nematode phylogenies with Iso-Seq, RNA-seq, ChIP-seq, and ATAC-seq to identify the gene structure and epigenetic signature of orphan genes in the satellite model nematode Pristionchus pacificus Consistent with previous findings, we find young genes are shorter, contain fewer exons, and are on average less strongly expressed than older genes. However, the subset of orphan genes that are expressed exhibit distinct chromatin states from similarly expressed conserved genes. Orphan gene transcription is determined by a lack of repressive histone modifications, confirming long-held hypotheses that open chromatin is important for new gene formation. Yet orphan gene start sites more closely resemble enhancers defined by H3K4me1, H3K27ac, and ATAC-seq peaks, in contrast to conserved genes that exhibit traditional promoters defined by H3K4me3 and H3K27ac. Although the majority of orphan genes are located on chromosome arms that contain high recombination rates and repressive histone marks, strongly expressed orphan genes are more randomly distributed. Our results support a model of new gene origination by rare integration into open chromatin near enhancers.

Alternative splicing links histone modifications to stem cell fate decision

BACKGROUND

Understanding the embryonic stem cell (ESC) fate decision between self-renewal and proper differentiation is important for developmental biology and regenerative medicine. Attention has focused on mechanisms involving histone modifications, alternative pre-messenger RNA splicing, and cell-cycle progression. However, their intricate interrelations and joint contributions to ESC fate decision remain unclear.

RESULTS

We analyze the transcriptomes and epigenomes of human ESC and five types of differentiated cells. We identify thousands of alternatively spliced exons and reveal their development and lineage-dependent characterizations. Several histone modifications show dynamic changes in alternatively spliced exons and three are strongly associated with 52.8% of alternative splicing events upon hESC differentiation. The histone modification-associated alternatively spliced genes predominantly function in G2/M phases and ATM/ATR-mediated DNA damage response pathway for cell differentiation, whereas other alternatively spliced genes are enriched in the G1 phase and pathways for self-renewal. These results imply a potential epigenetic mechanism by which some histone modifications contribute to ESC fate decision through the regulation of alternative splicing in specific pathways and cell-cycle genes. Supported by experimental validations and extended datasets from Roadmap/ENCODE projects, we exemplify this mechanism by a cell-cycle-related transcription factor, PBX1, which regulates the pluripotency regulatory network by binding to NANOG. We suggest that the isoform switch from PBX1a to PBX1b links H3K36me3 to hESC fate determination through the PSIP1/SRSF1 adaptor, which results in the exon skipping of PBX1.

CONCLUSION

We reveal the mechanism by which alternative splicing links histone modifications to stem cell fate decision.

The matrices and constraints of GT/AG splice sites of more than 1000 species/lineages

To provide a resource for the splice sites (SS) of different species, we calculated the matrices of nucleotide compositions of about 38 million splice sites from >1000 species/lineages. The matrices are enriched of aGGTAAGT (5'SS) or (Y)₆N(C/t)AG(g/a)t (3'SS) overall; however, they are quite diverse among hundreds of species. The diverse matrices remain prominent even under sequence selection pressures, suggesting the existence of diverse constraints as well as U snRNAs and other spliceosomal factors and/or their interactions with the splice sites. Using an algorithm to measure and compare the splice site constraints across all species, we demonstrate their distinct differences quantitatively. As an example of the resource's application to answering specific questions, we confirm that high constraints of particular positions are significantly associated with transcriptome-wide, increased occurrences of alternative splicing when uncommon nucleotides are present. More interestingly, the abundance of alternative splicing in 16 species correlates with the average constraint index of splice sites in a bell curve. This resource will allow users to assess specific sequences/splice sites against the consensus of every Ensembl-annotated species, and to explore the evolutionary changes or relationship to alternative splicing and transcriptome diversity. Web-search or update features are also included.

Mutually exclusive alternative splicing of pre-mRNAs

Pre-mRNA alternative splicing is an important mechanism used to expand protein diversity in higher eukaryotes, and mutually exclusive splicing is a specific type of alternative splicing in which only one of the exons in a cluster is included in functional transcripts. The most extraordinary example of this is the Drosophila melanogaster Down's syndrome cell adhesion molecule gene (Dscam), which potentially encodes 38,016 different isoforms through mutually exclusive splicing. Mutually exclusive splicing is a unique and challenging model that can be used to elucidate the evolution, regulatory mechanism, and function of alternative splicing. The use of new approaches has not only greatly expanded the mutually exclusive exome, but has also enabled the systematic analyses of single-cell alternative splicing during development. Furthermore, the identification of long-range RNA secondary structures provides a mechanistic framework for the regulation of mutually exclusive splicing (i.e., Dscam splicing). This article reviews recent insights into the identification, underlying mechanism, and roles of mutually exclusive splicing. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems.

In vitro and in vivo defensive effect of probiotic LAB against Pseudomonas aeruginosa using Caenorhabditis elegans model

This study aimed to investigate in vitro and in vivo the probiotic characteristics of lactic acid bacteria (LAB) isolated from Korean traditional fermented foods. Caenorhabditis elegans  (C. elegans) was used for analytical assays of fertility, chemotaxis, life-span, worm-killing and bacterial colonization in the intestinal lumen of the worm. All 35 strains of LAB reduced fertility and slowed development in the worms. The worm-killing assay showed that LAB significantly increased the lifespan (P < 0.05) and reduced the susceptibility to virulent PA14; however, the heat-killed LAB did not. The bacterial colonization assay revealed that LAB proliferated and protected the gut of the worm against infection by Pseudomonas aeruginosa PA14. In addition, specific LAB Pediococcus acidilactici(P. acidilactici DM-9), Pediococcus brevis (L. brevis SDL1411), and Pediococcus pentosaceus (P. pentosaceus SDL1409) strains showed acid resistance (66-91%), resistance to pepsin (64-67%) and viability in simulated intestinal fluid (67-73%) based on in vitro probiotic analyses. Taken together, these results suggest that C. elegans may be a tractable model for screening efficient probiotics.

Functions for fission yeast splicing factors SpSlu7 and SpPrp18 in alternative splice-site choice and stress-specific regulated splicing

Budding yeast spliceosomal factors ScSlu7 and ScPrp18 interact and mediate intron 3'ss choice during second step pre-mRNA splicing. The fission yeast genome with abundant multi-intronic transcripts, degenerate splice signals and SR proteins is an apt unicellular fungal model to deduce roles for core spliceosomal factors in alternative splice-site choice, intron retention and to study the cellular implications of regulated splicing. From our custom microarray data we deduce a stringent reproducible subset of S. pombe alternative events. We examined the role of factors SpSlu7 or SpPrp18 for these splice events and investigated the relationship to growth phase and stress. Wild-type log and stationary phase cells showed ats1+ exon 3 skipped and intron 3 retained transcripts. Interestingly the non-consensus 5'ss in ats1+ intron 3 caused SpSlu7 and SpPrp18 dependent intron retention. We validated the use of an alternative 5'ss in dtd1+ intron 1 and of an upstream alternative 3'ss in DUF3074 intron 1. The dtd1+ intron 1 non-canonical 5'ss yielded an alternative mRNA whose levels increased in stationary phase. Utilization of dtd1+ intron 1 sub-optimal 5' ss required functional SpPrp18 and SpSlu7 while compromise in SpSlu7 function alone hampered the selection of the DUF3074 intron 1 non canonical 3'ss. We analysed the relative abundance of these splice isoforms during mild thermal, oxidative and heavy metal stress and found stress-specific splice patterns for ats1+ and DUF3074 intron 1 some of which were SpSlu7 and SpPrp18 dependent. By studying ats1+ splice isoforms during compromised transcription elongation rates in wild-type, spslu7-2 and spprp18-5 mutant cells we found dynamic and intron context-specific effects in splice-site choice. Our work thus shows the combinatorial effects of splice site strength, core splicing factor functions and transcription elongation kinetics to dictate alternative splice patterns which in turn serve as an additional recourse of gene regulation in fission yeast.

Differential Gene Expression Profiles and Alternative Isoform Regulations in Gill of Nile Tilapia in Response to Acute Hypoxia

Fish often encounters exposures to acute environmental hypoxia either spatially or temporally. Gill organ plays important roles in response to hypoxic stress in fish. Few studies focus on the molecular regulation mechanisms of gills under hypoxic stress. In this study, we investigated the transcriptomic response to 12-h acute hypoxia in gill of a hypoxia tolerant fish, Nile tilapia Oreochromis niloticus through RNA sequencing (RNA-Seq). We sequenced messenger RNA from three control samples and three hypoxia-treated samples. Bioinformatics analysis identified 239 differentially expressed genes (DEG) and 34 genes (DUES) that had significant differential alternative isoform regulation events in at least one exonic region in gill in response to acute hypoxia. The spatiotemporal expression analysis in five tissues (heart, liver, brain, gill, and spleen) sampled at three time points (6, 12, and 24 h) under hypoxia treatment confirmed the significant association of differential exon usages in two DUES genes (TLDC2 and SSX2IPA) with hypoxia conditions. Further functional analysis suggested several energy and immune response-related pathways, e.g., metabolic pathway and antigen processing and presentation, contained the most abundant DEG genes. We found that some GO biological processes for DEG genes were significantly enriched under hypoxic stress, such as glycolysis, metabolic process, generation of precursor metabolites and energy, and cholesterol metabolic process. Our findings suggest abundant differential gene expression changes and alternative isoform regulation events in genes involved in the hypoxia response in gill. Our results provide a basis for exploring the gene regulation mechanism under hypoxic stress in fish.

Deep learning of the splicing (epi)genetic code reveals a novel candidate mechanism linking histone modifications to ESC fate decision

Alternative splicing (AS) is a genetically and epigenetically regulated pre-mRNA processing to increase transcriptome and proteome diversity. Comprehensively decoding these regulatory mechanisms holds promise in getting deeper insights into a variety of biological contexts involving in AS, such as development and diseases. We assembled splicing (epi)genetic code, DeepCode, for human embryonic stem cell (hESC) differentiation by integrating heterogeneous features of genomic sequences, 16 histone modifications with a multi-label deep neural network. With the advantages of epigenetic features, DeepCode significantly improves the performance in predicting the splicing patterns and their changes during hESC differentiation. Meanwhile, DeepCode reveals the superiority of epigenomic features and their dominant roles in decoding AS patterns, highlighting the necessity of including the epigenetic properties when assembling a more comprehensive splicing code. Moreover, DeepCode allows the robust predictions across cell lineages and datasets. Especially, we identified a putative H3K36me3-regulated AS event leading to a nonsense-mediated mRNA decay of BARD1. Reduced BARD1 expression results in the attenuation of ATM/ATR signalling activities and further the hESC differentiation. These results suggest a novel candidate mechanism linking histone modifications to hESC fate decision. In addition, when trained in different contexts, DeepCode can be expanded to a variety of biological and biomedical fields.

The combinatorial control of alternative splicing in C. elegans

Normal development requires the right splice variants to be made in the right tissues at the right time. The core splicing machinery is engaged in all splicing events, but which precise splice variant is made requires the choice between alternative splice sites—for this to occur, a set of splicing factors (SFs) must recognize and bind to short RNA motifs in the pre-mRNA. In C. elegans, there is known to be extensive variation in splicing patterns across development, but little is known about the targets of each SF or how multiple SFs combine to regulate splicing. Here we combine RNA-seq with in vitro binding assays to study how 4 different C. elegans SFs, ASD-1, FOX-1, MEC-8, and EXC-7, regulate splicing. The 4 SFs chosen all have well-characterised biology and well-studied loss-of-function genetic alleles, and all contain RRM domains. Intriguingly, while the SFs we examined have varied roles in C. elegans development, they show an unexpectedly high overlap in their targets. We also find that binding sites for these SFs occur on the same pre-mRNAs more frequently than expected suggesting extensive combinatorial control of splicing. We confirm that regulation of splicing by multiple SFs is often combinatorial and show that this is functionally significant. We also find that SFs appear to combine to affect splicing in two modes—they either bind in close proximity within the same intron or they appear to bind to separate regions of the intron in a conserved order. Finally, we find that the genes whose splicing are regulated by multiple SFs are highly enriched for genes involved in the cytoskeleton and in ion channels that are key for neurotransmission. Together, this shows that specific classes of genes have complex combinatorial regulation of splicing and that this combinatorial regulation is critical for normal development to occur.

Alternative splicing (AS) is a highly regulated process that is crucial for normal development. It requires the core splicing machinery, but the specific choice of splice site during AS is controlled by splicing factors (SFs) such as ELAV or RBFOX proteins that bind to specific sequences in pre-mRNAs to regulate usage of different splice sites. AS varies across the C. elegans life cycle and here we study how diverse SFs combine to regulate AS during C. elegans development. We selected 4 RRM-containing SFs that are all well studied and that have well-characterised loss-of-function genetic alleles. We find that these SFs regulate many of the same targets, and that combinatorial interactions between these SFs affect both individual splicing events and organism-level phenotypes including specific effects on the neuromuscular system. We further show that SFs combine to regulate splicing of an individual pre-mRNA in two distinct modes—either by binding in close proximity or by binding in a defined order on the pre-mRNA. Finally, we find that the genes whose splicing are most likely to be regulated by multiple SFs are genes that are required for the proper function of the neuromuscular system. These genes are also most likely to have changing AS patterns across development, suggesting that their splicing regulation is highly complex and developmentally regulated. Taken together, our data show that the precise splice variant expressed at any point in development is often the outcome of regulation by multiple SFs.

Quantitative RNA-seq meta-analysis of alternative exon usage in C. elegans

Almost 20 years after the completion of the C. elegans genome sequence, gene structure annotation is still an ongoing process with new evidence for gene variants still being regularly uncovered by additional in-depth transcriptome studies. While alternative splice forms can allow a single gene to encode several functional isoforms, the question of how much spurious splicing is tolerated is still heavily debated. Here we gathered a compendium of 1682 publicly available C. elegans RNA-seq data sets to increase the dynamic range of detection of RNA isoforms, and obtained robust measurements of the relative abundance of each splicing event. While most of the splicing reads come from reproducibly detected splicing events, a large fraction of purported junctions is only supported by a very low number of reads. We devised an automated curation method that takes into account the expression level of each gene to discriminate robust splicing events from potential biological noise. We found that rarely used splice sites disproportionately come from highly expressed genes and are significantly less conserved in other nematode genomes than splice sites with a higher usage frequency. Our increased detection power confirmed trans-splicing for at least 84% of C. elegans protein coding genes. The genes for which trans-splicing was not observed are overwhelmingly low expression genes, suggesting that the mechanism is pervasive but not fully captured by organism-wide RNA-seq. We generated annotated gene models including quantitative exon usage information for the entire C. elegans genome. This allows users to visualize at a glance the relative expression of each isoform for their gene of interest.

Acute hypoxia stress induced abundant differential expression genes and alternative splicing events in heart of tilapia

Hypoxia is one of the critical environmental stressors for fish in aquatic environments. Although accumulating evidences indicate that gene expression is regulated by hypoxia stress in fish, how genes undergoing differential gene expression and/or alternative splicing (AS) in response to hypoxia stress in heart are not well understood. Using RNA-seq, we surveyed and detected 289 differential expressed genes (DEG) and 103 genes that undergo differential usage of exons and splice junctions events (DUES) in heart of a hypoxia tolerant fish, Nile tilapia, Oreochromis niloticus following 12h hypoxic treatment. The spatio-temporal expression analysis validated the significant association of differential exon usages in two randomly selected DUES genes (fam162a and ndrg2) in 5 tissues (heart, liver, brain, gill and spleen) sampled at three time points (6h, 12h, and 24h) under acute hypoxia treatment. Functional analysis significantly associated the differential expressed genes with the categories related to energy conservation, protein synthesis and immune response. Different enrichment categories were found between the DEG and DUES dataset. The Isomerase activity, Oxidoreductase activity, Glycolysis and Oxidative stress process were significantly enriched for the DEG gene dataset, but the Structural constituent of ribosome and Structural molecule activity, Ribosomal protein and RNA binding protein were significantly enriched only for the DUES genes. Our comparative transcriptomic analysis reveals abundant stress responsive genes and their differential regulation function in the heart tissues of Nile tilapia under acute hypoxia stress. Our findings will facilitate future investigation on transcriptome complexity and AS regulation during hypoxia stress in fish.

Identification and analysis of alternative splicing events in Phaseolus vulgaris and Glycine max

BACKGROUND

The vast diversification of proteins in eukaryotic cells has been related with multiple transcript isoforms from a single gene that result in alternative splicing (AS) of primary transcripts. Analysis of RNA sequencing data from expressed sequence tags and next generation RNA sequencing has been crucial for AS identification and genome-wide AS studies. For the identification of AS events from the related legume species Phaseolus vulgaris and Glycine max, 157 and 88 publicly available RNA-seq libraries, respectively, were analyzed.

RESULTS

We identified 85,570 AS events from P. vulgaris in 72% of expressed genes and 134,316 AS events in 70% of expressed genes from G. max. These were categorized in seven AS event types with intron retention being the most abundant followed by alternative acceptor and alternative donor, representing ~75% of all AS events in both plants. Conservation of AS events in homologous genes between the two species was analyzed where an overrepresentation of AS affecting 5'UTR regions was observed for certain types of AS events. The conservation of AS events was experimentally validated for 8 selected genes, through RT-PCR analysis. The different types of AS events also varied by relative position in the genes. The results were consistent in both species.

CONCLUSIONS

The identification and analysis of AS events are first steps to understand their biological relevance. The results presented here from two related legume species reveal high conservation, over ~15-20 MY of divergence, and may point to the biological relevance of AS.

Differential alternative splicing coupled to nonsense-mediated decay of mRNA ensures dietary restriction-induced longevity

Alternative splicing (AS) coupled to nonsense-mediated decay (AS-NMD) is a conserved mechanism for post-transcriptional gene regulation. Here we show that, during dietary restriction (DR), AS is enhanced in Caenorhabditis elegans and mice. A splicing mediator hrpu-1 regulates a significant part of these AS events in C. elegans; knocking it down suppresses DR-mediated longevity. Concurrently, due to increased AS, NMD pathway genes are upregulated and knocking down UPF1 homologue smg-2 suppresses DR lifespan. Knockdown of NMD during DR significantly increases the inclusion of PTC-containing introns and the lengths of the 3′UTRs. Finally, we demonstrate that PHA-4/FOXA transcriptionally regulates the AS-NMD genes. Our study suggests that DR uses AS to amplify the proteome, supporting physiological remodelling required for enhanced longevity. This increases the dependence on NMD, but also helps fine-tune the expression of metabolic and splicing mediators. AS-NMD may thus provide an energetically favourable level of dynamic gene expression control during dietary restriction.

Alternative splicing coupled to nonsense-mediated decay (AS-NMD) is a conserved mechanism for post-transcriptional gene regulation. Here, the authors provide evidence that AS-NMD is enhanced during dietary restriction (DR) and is required for DR-mediated longevity assurance in C. elegans.