Genome-wide identification and comparative analysis of alternative splicing across four legume species

Alternative splicing shapes the transcriptome complexity in blackgram [Vigna mungo (L.) Hepper]

Vigna mungo, a highly consumed crop in the pan-Asian countries, is vulnerable to several biotic and abiotic stresses. Understanding the post-transcriptional gene regulatory cascades, especially alternative splicing (AS), may underpin large-scale genetic improvements to develop stress-resilient varieties. Herein, a transcriptome based approach was undertaken to decipher the genome-wide AS landscape and splicing dynamics in order to establish the intricacies of their functional interactions in various tissues and stresses. RNA sequencing followed by high-throughput computational analyses identified 54,526 AS events involving 15,506 AS genes that generated 57,405 transcripts isoforms. Enrichment analysis revealed their involvement in diverse regulatory functions and demonstrated that transcription factors are splicing-intensive, splice variants of which are expressed differentially across tissues and environmental cues. Increased expression of a splicing regulator NHP2L1/SNU13 was found to co-occur with lower intron retention events. The host transcriptome is significantly impacted by differential isoform expression of 1172 and 765 AS genes that resulted in 1227 (46.8% up and 53.2% downregulated) and 831 (47.5% up and 52.5% downregulated) transcript isoforms under viral pathogenesis and Fe²⁺ stressed condition, respectively. However, genes experiencing AS operate differently from the differentially expressed genes, suggesting AS is a unique and independent mode of regulatory mechanism. Therefore, it can be inferred that AS mediates a crucial regulatory role across tissues and stressful situations and the results would provide an invaluable resource for future endeavours in V. mungo genomics.

Deciphering the role of alternative splicing as modulators of defense response in the MYMIV- Vigna mungo pathosystem

Alternative splicing (AS) is a crucial regulatory mechanism that impacts transcriptome and proteome complexity under stressful situations. Although its role in abiotic stresses is somewhat understood, our understanding of the mechanistic regulation of pre-mRNA splicing in plant-pathogen interaction is meagre. To comprehend this unexplored immune reprogramming mechanism, transcriptome profiles of Mungbean Yellow Mosaic India Virus (MYMIV)-resistant and susceptible Vigna mungo genotypes were analysed for AS genes that may underlie the resistance mechanism. Results revealed a repertoire of AS-isoforms accumulated during pathogenic infestation, with intron retention being the most common AS mechanism. Identification of 688 differential alternatively spliced (DAS) genes in the resistant host elucidates its robust antiviral response, whereas 322 DAS genes were identified in the susceptible host. Enrichment analyses confirmed DAS transcripts pertaining to stress, signalling, and immune system pathways have undergone maximal perturbations. Additionally, a strong regulation of the splicing factors has been observed both at transcriptional and post-transcriptional levels. qPCR validation of candidate DAS transcripts with induced expression upon MYMIV-infection demonstrated a competent immune response in the resistant background. The AS-impacted genes resulted either in partial/complete loss of functional domains or altered sensitivity to miRNA-mediated gene silencing. A complex regulatory module, miR7517-ATAF2, has been identified in an aberrantly spliced ATAF2 isoform that exposes an intronic miR7517 binding site, thereby suppressing the negative regulator to enhance defense reaction. The present study establishes AS as a non-canonical immune reprogramming mechanism that operates in parallel, thereby offering an alternative strategy for developing yellow mosaic-resistant V. mungo cultivars. This article is protected by copyright. All rights reserved.

OsGRF4AA compromises heat tolerance of developing pollen grains in rice

Extreme high temperature at the meiosis stage causes a severe decrease in spikelet fertility and grain yield in rice. The rice variety grain size on chromosome 2 (GS2) contains sequence variations of OsGRF4 (Oryza sativa growth-regulating factor 4; OsGRF4^AA ), escaping the microRNA miR396-mediated degradation of this gene at the mRNA level. Accumulation of OsGRF4 enhances nitrogen usage and metabolism, and increases grain size and grain yield. In this study, we found that pollen viability and seed-setting rate under heat stress (HS) decreased more seriously in GS2 than in its comparator, Zhonghua 11 (ZH11). Transcriptomic analysis revealed that, following HS, genes related to carbohydrate metabolic processes were expressed and regulated differentially in the anthers of GS2 and ZH11. Moreover, the expression of genes involved in chloroplast development and photosynthesis, lipid metabolism, and key transcription factors, including eight male sterile genes, were inhibited by HS to a greater extent in GS2 than in ZH11. Interestingly, pre-mRNAs of OsGRF4, and a group of essential genes involved in development and fertilization, were differentially spliced in the anthers of GS2 and ZH11. Taken together, our results suggest that variation in OsGRF4 affects proper transcriptional and splicing regulation of genes under HS, and that this can be mediated by, and also feed back to, carbohydrate and nitrogen metabolism, resulting in a reduction in the heat tolerance of rice anthers.

Genome-wide profiling of alternative splicing genes in hybrid poplar (P.alba×P.glandulosa cv.84K) leaves

Alternative splicing (AS) is a post-transcriptional process common in plants and essential for regulation of environmental fitness of plants. In the present study, we focus on the AS events in poplar leaves to understand their effects on plant growth and development. The hybrid poplar (P.alba×P.glandulosa cv.84K) leaves were collected for RNA extraction. The extracted RNA was sequenced using on an Illumina HiSeq™ 2000 platform. Using the Populus trichocarpa genome as the reference, a total of 3810 AS genes were identified (9225 AS events), which accounted for 13.51% of all the expressed genes. Intron retention was the most common AS event, accounting for 43.86% of all the AS events, followed by alternative 3′ splice sites (23.75%), alternative 5′ splice sites (23.71%), and exon skipping (8.68%). Chromosomes 10 had the most condensed AS events (33.67 events/Mb) and chromosome 19 had the least (12.42 events/Mb). Association analysis showed that AS in the poplar leaves was positively correlated with intron length, exon number, exon length, and gene expression level, and was negatively correlated with GC content. AS genes in the poplar leaves were associated mainly with inositol phosphate metabolism and phosphatidylinositol signaling system pathways that would be significant on wooden plant production.

To keep or not to keep: mRNA stability and translatability in root nodule symbiosis

Post-transcriptional control of gene expression allows plants to rapidly adapt to changes in their environment. Under low nitrogen conditions, legume plants engage into a symbiosis with soil bacteria that results in the formation of root nodules, where bacteria are allocated and fix atmospheric nitrogen for the plant's benefit. Recent studies highlighted the importance of small RNA-mediated mechanisms in the control of bacterial infection, nodule organogenesis, and the long-distance signaling that balances plant growth and nodulation. Examples of such mechanisms are shoot-to-root mobile microRNAs and small RNA fragments derived from degradation of bacterial transfer RNAs that repress complementary mRNAs in the host plant. Mechanisms of selective mRNA translation also contribute to rapidly modulate the expression of nodulation genes in a cell-specific manner during symbiosis. Here, the most recent advances made on the regulation of mRNA stability and translatability, and the emerging roles of long non-coding RNAs in symbiosis are summarized.

Analysis and comprehensive comparison of PacBio and nanopore-based RNA sequencing of the Arabidopsis transcriptome

BACKGROUND

The number of studies using third-generation sequencing utilising Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT) is rapidly increasing in many different research areas. Among them, plant full-length single-molecule transcriptome studies have mostly used PacBio sequencing, whereas ONT is rarely used. Therefore, in this study, we examined ONT RNA sequencing methods in plants. We performed a detailed evaluation of reads from PacBio, Nanopore direct cDNA (ONT Dc), and Nanopore PCR cDNA (ONT Pc) sequencing including characteristics of raw data and identification of transcripts. In addition, matched Illumina data were generated for comparison.

RESULTS

ONT Pc showed overall better raw data quality, whereas PacBio generated longer read lengths. In the transcriptome analysis, PacBio and ONT Pc performed similarly in transcript identification, simple sequence repeat analysis, and long non-coding RNA prediction. PacBio was superior in identifying alternative splicing events, whereas ONT Pc could estimate transcript expression levels.

CONCLUSIONS

This paper made a comprehensive comparison of PacBio and nanopore-based RNA sequencing of the Arabidopsis transcriptome, the results indicate that ONT Pc is more cost-effective for generating extremely long reads and can characterise the transcriptome as well as quantify transcript expression. Therefore, ONT Pc is a new cost-effective and worthwhile method for full-length single-molecule transcriptome analysis in plants.

Alternative Splicing Enhances the Transcriptome Complexity of Liriodendron chinense

Alternative splicing (AS) plays pivotal roles in regulating plant growth and development, flowering, biological rhythms, signal transduction, and stress responses. However, no studies on AS have been performed in Liriodendron chinense, a deciduous tree species that has high economic and ecological value. In this study, we used multiple tools and algorithms to analyze transcriptome data derived from seven tissues via hybrid sequencing. Although only 17.56% (8,503/48,408) of genes in L. chinense were alternatively spliced, these AS genes occurred in 37,844 AS events. Among these events, intron retention was the most frequent AS event, producing 1,656 PTC-containing and 3,310 non-PTC-containing transcripts. Moreover, 183 long noncoding RNAs (lncRNAs) also underwent AS events. Furthermore, weighted gene coexpression network analysis (WGCNA) revealed that there were great differences in the activities of transcription and post-transcriptional regulation between pistils and leaves, and AS had an impact on many physiological and biochemical processes in L. chinense, such as photosynthesis, sphingolipid metabolism, fatty acid biosynthesis and metabolism. Moreover, our analysis showed that the features of genes may affect AS, as AS genes and non-AS genes had differences in the exon/intron length, transcript length, and number of exons/introns. In addition, the structure of AS genes may impact the frequencies and types of AS because AS genes with more exons or introns tended to exhibit more AS events, and shorter introns tended to be retained, whereas shorter exons tended to be skipped. Furthermore, eight AS genes were verified, and the results were consistent with our analysis. Overall, this study reveals that AS and gene interaction are mutual-on one hand, AS can affect gene expression and translation, while on the other hand, the structural characteristics of the gene can also affect AS. This work is the first to comprehensively report on AS in L. chinense, and it can provide a reference for further research on AS in L. chinense.

Transcriptomic and Co-Expression Network Profiling of Shoot Apical Meristem Reveal Contrasting Response to Nitrogen Rate between Indica and Japonica Rice Subspecies

Reducing nitrogen (N) input is a key measure to achieve a sustainable rice production in China, especially in Jiangsu Province. Tiller is the basis for achieving panicle number that plays as a major factor in the yield determination. In actual production, excessive N is often applied in order to produce enough tillers in the early stages. Understanding how N regulates tillering in rice plants is critical to generate an integrative management to reduce N use and reaching tiller number target. Aiming at this objective, we utilized RNA sequencing and weighted gene co-expression network analysis (WGCNA) to compare the transcriptomes surrounding the shoot apical meristem of indica (Yangdao6, YD6) and japonica (Nipponbare, NPB) rice subspecies. Our results showed that N rate influenced tiller number in a different pattern between the two varieties, with NPB being more sensitive to N enrichment, and YD6 being more tolerant to high N rate. Tiller number was positively related to N content in leaf, culm and root tissue, but negatively related to the soluble carbohydrate content, regardless of variety. Transcriptomic comparisons revealed that for YD6 when N rate enrichment from low (LN) to medium (MN), it caused 115 DEGs (LN vs. MN), from MN to high level (HN) triggered 162 DEGs (MN vs. HN), but direct comparison of low with high N rate showed a 511 DEGs (LN vs. HN). These numbers of DEG in NPB were 87 (LN vs. MN), 40 (MN vs. HN), and 148 (LN vs. HN). These differences indicate that continual N enrichment led to a bumpy change at the transcription level. For the reported sixty-five genes which affect tillering, thirty-six showed decent expression in SAM at tiller starting phase, among them only nineteen being significantly influenced by N level, and two genes showed significant interaction between N rate and variety. Gene ontology analysis revealed that the majority of the common DEGs are involved in general stress responses, stimulus responses, and hormonal signaling process. WGCNA network identified twenty-two co-expressing gene modules and ten candidate hubgenes for each module. Several genes associated with tillering and N rate fall on the related modules. These indicate that there are more genes participating in tillering regulation in response to N enrichment.

Genome-wide analysis of alternative splicing divergences between Brassica hexaploid and its parents

Compared with its parents, Brassica hexaploid underwent significant AS changes, which may provide diversified gene expression regulation patterns and could enhance its adaptability during evolution Polyploidization is considered a significant evolution force that promotes species formation. Alternative splicing (AS) plays a crucial role in multiple biological processes during plant growth and development. To explore the effects of allopolyploidization on the AS patterns of genes, a genome-wide AS analysis was performed by RNA-seq in Brassica hexaploid and its parents. In total, we found 7913 (27540 AS events), 14447 (70179 AS events), and 13205 (60804 AS events) AS genes in Brassica rapa, Brassica carinata, and Brassica hexaploid, respectively. A total of 920 new AS genes were discovered in Brassica hexaploid. There were 56 differently spliced genes between Brassica hexaploid and its parents. In addition, most of the alternative 5' splice sites were located 4 bp upstream of the dominant 5' splice sites, and most of the alternative 3' splice sites were located 3 bp downstream of the dominant 3' splice sites in Brassica hexapliod, which was similar to B. carinata. Furthermore, we cloned and sequenced all amplicons from the RT-PCR products of GRP7/8, namely, Bol045859, Bol016025 and Bol02880. The three genes were found to produce AS transcripts in a new way. The AS patterns of genes were diverse between Brassica hexaploid and its parents, including the loss and gain of AS events. Allopolyploidization changed alternative splicing sites of pre-mRNAs in Brassica hexaploid, which brought about alterations in the sequences of transcripts. Our study provided novel insights into the AS patterns of genes in allopolyploid plants, which may provide a reference for the study of polyploidy adaptability.