Designing alternative splicing RNA-seq studies. Beyond generic guidelines

Chitosan induces differential transcript usage of chitosanase 3 encoding gene (csn3) in the biocontrol fungus Pochonia chlamydosporia 123

Background

Pochonia chlamydosporia is an endophytic fungus used for nematode biocontrol that employs its cellular and molecular machinery to degrade the nematode egg-shell. Chitosanases, among other enzymes, are involved in this process. In this study, we improve the genome sequence assembly of P. chlamydosporia 123, by utilizing long Pacific Biosciences (PacBio) sequence reads. Combining this improved genome assembly with previous RNA-seq data revealed alternative isoforms of a chitosanase in the presence of chitosan. This study could open new insights into understanding fungal resistance to chitosan and root-knot nematode (RKN) egg infection processes.

Results

The P. chlamydosporia 123 genome sequence assembly has been updated using long-read PacBio sequencing and now includes 12,810 predicted protein-coding genes. Compared with the previous assembly based on short reads, there are 701 newly annotated genes, and 69 previous genes are now split. Eight of the new genes were differentially expressed in fungus interactions with Meloidogyne javanica eggs or chitosan.

A survey of the RNA-seq data revealed alternative splicing in the csn3 gene that encodes a chitosanase, with four putative splicing variants: csn3_v1, csn3_v2, csn3_v3 and csn3_v4. When P. chlamydosporia is treated with 0.1 mg·mL^− 1 chitosan for 4 days, csn3 is expressed 10-fold compared with untreated controls. Furthermore, the relative abundances of each of the four transcripts are different in chitosan treatment compared with controls. In controls, the abundances of each transcript are nil, 32, 55, and 12% for isoforms csn3_v1, csn3_v2, csn3_v3 and csn3_v4 respectively. Conversely, in chitosan-treated P. chlamydosporia, the abundances are respectively 80, 15%, 2—3%, 2—3%. Since isoform csn3_v1 is expressed with chitosan only, the putatively encoded enzyme is probably induced and likely important for chitosan degradation.

Conclusions

Alternative splicing events have been discovered and described in the chitosanase 3 encoding gene from P. chlamydosporia 123. Gene csn3 takes part in RKN parasitism process and chitosan enhances its expression. The isoform csn3_v1 would be related to the degradation of this polymer in bulk form, while other isoforms may be related to the degradation of chitosan in the nematode egg-shell.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08232-7.

Experimental Design for Time-Series RNA-Seq Analysis of Gene Expression and Alternative Splicing

RNA-sequencing (RNA-seq) is currently the method of choice for analysis of differential gene expression. To fully exploit the wealth of data generated from genome-wide transcriptomic approaches, the initial design of the experiment is of paramount importance. Biological rhythms in nature are pervasive and are driven by endogenous gene networks collectively known as circadian clocks. Measuring circadian gene expression requires time-course experiments which take into account time-of-day factors influencing variability in expression levels. We describe here an approach for characterizing diurnal changes in expression and alternative splicing for plants undergoing cooling. The method uses inexpensive everyday laboratory equipment and utilizes an RNA-seq application (3D RNA-seq) that can handle complex experimental designs and requires little or no prior bioinformatics expertise.

SpliceDetector: a software for detection of alternative splicing events in human and model organisms directly from transcript IDs

In eukaryotes, different combinations of exons lead to multiple transcripts with various functions in protein level, in a process called alternative splicing (AS). Unfolding the complexity of functional genomics through genome-wide profiling of AS and determining the altered ultimate products provide new insights for better understanding of many biological processes, disease progress as well as drug development programs to target harmful splicing variants. The current available tools of alternative splicing work with raw data and include heavy computation. In particular, there is a shortcoming in tools to discover AS events directly from transcripts. Here, we developed a Windows-based user-friendly tool for identifying AS events from transcripts without the need to any advanced computer skill or database download. Meanwhile, due to online working mode, our application employs the updated SpliceGraphs without the need to any resource updating. First, SpliceGraph forms based on the frequency of active splice sites in pre-mRNA. Then, the presented approach compares query transcript exons to SpliceGraph exons. The tool provides the possibility of statistical analysis of AS events as well as AS visualization compared to SpliceGraph. The developed application works for transcript sets in human and model organisms.

Mettl3-mediated m6A regulates spermatogonial differentiation and meiosis initiation

METTL3 catalyzes the formation of N⁶-methyl-adenosine (m⁶A) which has important roles in regulating various biological processes. However, the in vivo function of Mettl3 remains largely unknown in mammals. Here we generated germ cell-specific Mettl3 knockout mice and demonstrated that Mettl3 was essential for male fertility and spermatogenesis. The ablation of Mettl3 in germ cells severely inhibited spermatogonial differentiation and blocked the initiation of meiosis. Transcriptome and m⁶A profiling analysis revealed that genes functioning in spermatogenesis had altered profiles of expression and alternative splicing. Our findings provide novel insights into the function and regulatory mechanisms of Mettl3-mediated m⁶A modification in spermatogenesis and reproduction in mammals.

Nova1 mediates resistance of rat pheochromocytoma cells to hypoxia-induced apoptosis via the Bax/Bcl-2/caspase-3 pathway

Neuro-oncological ventral antigen 1 (Nova1) is a well known brain-specific splicing factor. Several studies have identified Nova1 as a regulatory protein at the top of a hierarchical network. However, the function of Nova1 during hypoxia remains poorly understood. This study aimed to investigate the protective effect of Nova1 against cell hypoxia and to further explore the Bax/Bcl-2/caspase-3 pathway as a potential mechanism. During hypoxia, the survival rate of pheochromocytoma PC12 cells was gradually decreased and the apoptosis rate was gradually increased, peaking at 48 h of hypoxia. At 48 h after transfection of PC12 cells with pCMV-Myc-Nova1, the expression of Nova1 was significantly increased, with wide distribution in the cytoplasm and nucleus. Moreover, the survival rate was significantly increased and the apoptosis rate was significantly decreased. Additionally, the mRNA and protein expression levels of Bax and caspase-3 were significantly increased in the pCMV-Myc group and significantly decreased in the pCMV-Myc-Nova1 group, whereas that of Bcl-2 was significantly decreased in the pCMV-Myc group and significantly increased in the pCMV-Myc-Nova1 group. This study indicated that Nova1 could be linked to resistance to the hypoxia-induced apoptosis of PC12 cells via the Bax/Bcl-2/caspase-3 pathway, and this finding may be of significance for exploring novel mechanisms of hypoxia and the treatment of hypoxia-associated diseases.

Detecting circular RNAs: bioinformatic and experimental challenges

The pervasive expression of circular RNAs (circRNAs) is a recently discovered feature of gene expression in highly diverged eukaryotes. Numerous algorithms that are used to detect genome-wide circRNA expression from RNA sequencing (RNA-seq) data have been developed in the past few years, but there is little overlap in their predictions and no clear gold-standard method to assess the accuracy of these algorithms. We review sources of experimental and bioinformatic biases that complicate the accurate discovery of circRNAs and discuss statistical approaches to address these biases. We conclude with a discussion of the current experimental progress on the topic.