Alternative Splicing Plays a Critical Role in Maintaining Mineral Nutrient Homeostasis in Rice (Oryza sativa)

Integrated analyses of miRNAome and transcriptome reveal zinc deficiency responses in rice seedlings

BACKGROUND

Zinc (Zn) deficiency is one of the most widespread soil constraints affecting rice productivity, but the molecular mechanisms underlying the regulation of Zn deficiency response is still limited. Here, we aim to understand the molecular mechanisms of Zn deficiency response by integrating the analyses of the global miRNA and mRNA expression profiles under Zn deficiency and resupply in rice seedlings by integrating Illumina's high-throughput small RNA sequencing and transcriptome sequencing.

RESULTS

The transcriptome sequencing identified 360 genes that were differentially expressed in the shoots and roots of Zn-deficient rice seedlings, and 97 of them were recovered after Zn resupply. A total of 68 miRNAs were identified to be differentially expressed under Zn deficiency and/or Zn resupply. The integrated analyses of miRNAome and transcriptome data showed that 12 differentially expressed genes are the potential target genes of 10 Zn-responsive miRNAs such as miR171g-5p, miR397b-5p, miR398a-5p and miR528-5p. Some miRNA genes and differentially expressed genes were selected for validation by quantitative RT-PCR, and their expressions were similar to that of the sequencing results.

CONCLUSION

These results provide insights into miRNA-mediated regulatory pathways in Zn deficiency response, and provide candidate genes for genetic improvement of Zn deficiency tolerance in rice.

Phosphoproteomic Profiling Reveals the Importance of CK2, MAPKs and CDPKs in Response to Phosphate Starvation in Rice

Phosphorus is one of the most important macronutrients required for plant growth and development. The importance of phosphorylation modification in regulating phosphate (Pi) homeostasis in plants is emerging. We performed phosphoproteomic profiling to characterize proteins whose degree of phosphorylation is altered in response to Pi starvation in rice root. A subset of 554 proteins, including 546 down-phosphorylated and eight up-phosphorylated proteins, exhibited differential phosphorylation in response to Pi starvation. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis with the differentially phosphorylated proteins indicated that RNA processing, transport, splicing and translation and carbon metabolism played critical roles in response to Pi starvation in rice. Levels of phosphorylation of four mitogen-activated protein kinases (MAPKs), including OsMAPK6, five calcium-dependent protein kinases (CDPKs) and OsCK2β3 decreased in response to Pi starvation. The decreased phosphorylation level of OsMAPK6 was confirmed by Western blotting. Mutation of OsMAPK6 led to Pi accumulation under Pi-sufficient conditions. Motif analysis indicated that the putative MAPK, casein kinase 2 (CK2) and CDPK substrates represented about 54.4%, 21.5% and 4.7%, respectively, of the proteins exhibiting differential phosphorylation. Based on the motif analysis, 191, 151 and 46 candidate substrates for MAPK, CK2 and CDPK were identified. These results indicate that modification of phosphorylation profiles provides complementary information on Pi-starvation-induced processes, with CK2, MAPK and CDPK protein kinase families playing key roles in these processes in rice.

Multiplex CRISPR Mutagenesis of the Serine/Arginine-Rich (SR) Gene Family in Rice

Plant growth responds to various environmental and developmental cues via signaling cascades that influence gene expression at the level of transcription and pre-mRNA splicing. Alternative splicing of pre-mRNA increases the coding potential of the genome from multiexon genes and regulates gene expression through multiple mechanisms. Serine/arginine-rich (SR) proteins, a conserved family of splicing factors, are the key players of alternative splicing and regulate pre-mRNA splicing under stress conditions. The rice (Oryza sativa) genome encodes 22 SR proteins categorized into six subfamilies. Three of the subfamilies are plant-specific with no mammalian orthologues, and the functions of these SR proteins are not well known. The clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) system is a genome engineering tool that cleaves the target DNA at specific locations directed by a guide RNA (gRNA). Recent advances in CRISPR/Cas9-mediated plant genome engineering make it possible to generate single and multiple functional knockout mutants in diverse plant species. In this study, we targeted each rice SR locus and produced single knockouts. To overcome the functional redundancy within each subfamily of SR genes, we utilized a polycistronic tRNA-gRNA multiplex targeting system and targeted all loci of each subfamily. Sanger sequencing results indicated that most of the targeted loci had knockout mutations. This study provides useful resource materials for understanding the molecular role of SR proteins in plant development and biotic and abiotic stress responses.

Serine/Arginine-rich protein family of splicing regulators: New approaches to study splice isoform functions

Serine/arginine-rich (SR) proteins are conserved RNA-binding proteins that play major roles in RNA metabolism. They function as molecular adaptors, facilitate spliceosome assembly and modulate constitutive and alternative splicing of pre-mRNAs. Pre-mRNAs encoding SR proteins and many other proteins involved in stress responses are extensively alternatively spliced in response to diverse stresses. Hence, it is proposed that stress-induced changes in splice isoforms contribute to the adaptation of plants to stress responses. However, functions of most SR genes and their splice isoforms in stress responses are not known. Lack of easy and robust tools hindered the progress in this area. Emerging technologies such as CRISPR/Cas9 will facilitate studies of SR function by enabling the generation of single and multiple knock-out mutants of SR subfamily members. Moreover, CRISPR/Cas13 allows targeted manipulation of splice isoforms from SR and other genes in a constitutive or tissue-specific manner to evaluate functions of individual splice variants. Identification of the in vivo targets of SR proteins and their splice variants using the recently developed TRIBE (Targets of RNA-binding proteins Identified By Editing) and other methods will help unravel their mode of action and splicing regulatory elements under various conditions. These new approaches are expected to provide significant new insights into the roles of SRs and splice isoforms in plants adaptation to diverse stresses.

OsbZIP81, A Homologue of Arabidopsis VIP1, May Positively Regulate JA Levels by Directly Targetting the Genes in JA Signaling and Metabolism Pathway in Rice

Rice (Oryza sativa L.) is one of the most important food crops in the world. In plants, jasmonic acid (JA) plays essential roles in response to biotic and abiotic stresses. As one of the largest transcription factors (TFs), basic region/leucine zipper motif (bZIP) TFs play pivotal roles through the whole life of plant growth. However, the relationship between JA and bZIP TFs were rarely reported, especially in rice. In this study, we found two rice homologues of Arabidopsis VIP1 (VirE2-interacting protein 1), OsbZIP81, and OsbZIP84. OsbZIP81 has at least two alternative transcripts, OsbZIP81.1 and OsbZIP81.2. OsbZIP81.1 and OsbZIP84 are typical bZIP TFs, while OsbZIP81.2 is not. OsbZIP81.1 can directly bind OsPIOX and activate its expression. In OsbZIP81.1 overexpression transgenic rice plant, JA (Jasmonic Acid) and SA (Salicylic acid) were up-regulated, while ABA (Abscisic acid) was down-regulated. Moreover, Agrobacterium, Methyl Jasmonic Acid (MeJA), and PEG6000 can largely induce OsbZIP81. Based on ChIP-Seq and Random DNA Binding Selection Assay (RDSA), we identified a novel cis-element OVRE (Oryza VIP1 response element). Combining ChIP-Seq and RNA-Seq, we obtained 1332 targeted genes that were categorized in biotic and abiotic responses, including α-linolenic acid metabolism and fatty acid degradation. Together, these results suggest that OsbZIP81 may positively regulate JA levels by directly targeting the genes in JA signaling and metabolism pathway in rice.

Brassica napus Infected with Leptosphaeria maculans

Alternative splicing (AS) is a post-transcriptional regulatory process that enhances transcriptome diversity, thereby affecting plant growth, development, and stress responses. To identify the new transcripts and changes in the isoform-level AS landscape of rapeseed (Brassica napus) infected with the fungal pathogen Leptosphaeria maculans, we compared eight RNA-seq libraries prepared from mock-inoculated and inoculated B. napus cotyledons and stems. The AS events that occurred in stems were almost the same as those in cotyledons, with intron retention representing the most common AS pattern. We identified 1892 differentially spliced genes between inoculated and uninoculated plants. We performed a weighted gene co-expression network analysis (WGCNA) to identify eight co-expression modules and their Hub genes, which are the genes most connected with other genes within each module. There are nine Hub genes, encoding nine transcription factors, which represent key regulators of each module, including members of the NAC, WRKY, TRAF, AP2/ERF-ERF, C2H2, C2C2-GATA, HMG, bHLH, and C2C2-CO-like families. Finally, 52 and 117 alternatively spliced genes in cotyledons and stems were also differentially expressed between mock-infected and infected materials, such as HMG and C2C2-Dof; which have dual regulatory mechanisms in response to L. maculans. The splicing of the candidate genes identified in this study could be exploited to improve resistance to L. maculans.

Comprehensive identification of the full-length transcripts and alternative splicing related to the secondary metabolism pathways in the tea plant (Camellia sinensis)

Flavonoids, theanine and caffeine are the main secondary metabolites of the tea plant (Camellia sinensis), which account for the tea's unique flavor quality and health benefits. The biosynthesis pathways of these metabolites have been extensively studied at the transcriptional level, but the regulatory mechanisms are still unclear. In this study, to explore the transcriptome diversity and complexity of tea plant, PacBio Iso-Seq and RNA-seq analysis were combined to obtain full-length transcripts and to profile the changes in gene expression during the leaf development. A total of 1,388,066 reads of insert (ROI) were generated with an average length of 1,762 bp, and more than 54% (755,716) of the ROIs were full-length non-chimeric (FLNC) reads. The Benchmarking Universal Single-Copy Orthologue (BUSCO) completeness was 92.7%. A total of 93,883 non-redundant transcripts were obtained, and 87,395 (93.1%) were new alternatively spliced isoforms. Meanwhile, 7,650 differential expression transcripts (DETs) were identified. A total of 28,980 alternative splicing (AS) events were predicted, including 1,297 differential AS (DAS) events. The transcript isoforms of the key genes involved in the flavonoid, theanine and caffeine biosynthesis pathways were characterized. Additionally, 5,777 fusion transcripts and 9,052 long non-coding RNAs (lncRNAs) were also predicted. Our results revealed that AS potentially plays a crucial role in the regulation of the secondary metabolism of the tea plant. These findings enhanced our understanding of the complexity of the secondary metabolic regulation of tea plants and provided a basis for the subsequent exploration of the regulatory mechanisms of flavonoid, theanine and caffeine biosynthesis in tea plants.