Alternative Splicing as a Regulator of Early Plant Development

Evidence Supporting a Role of Alternative Splicing Participates in Melon (Cucumis melo L.) Fruit Ripening

One key post-transcriptional modification mechanism that dynamically controls a number of physiological processes in plants is alternative splicing (AS). However, the functional impacts of AS on fruit ripening remain unclear. In this research, we used RNA-seq data from climacteric (VED, Harukei 3) and non-climacteric (PI, PS) melon cultivars to explore alternative splicing (AS) in immature and mature fruit. The results revealed dramatic changes in differential AS genes (DAG) between the young and mature fruit stages, particularly in genes involved in fruit development/ripening, carotenoid and capsaicinoid biosynthesis, and starch and sucrose metabolism. Serine/arginine-rich (SR) family proteins are known as important splicing factors in AS events. From the melon genome, a total of 17 SR members were discovered in this study. These genes could be classified into eight distinct subfamilies based on gene structure and conserved motifs. Promoter analysis detected various cis-acting regulatory elements involved in hormone pathways and fruit development. Interestingly, these SR genes exhibited specific expression patterns in reproductive organs such as flowers and ovaries. Additionally, concurrent with the increase in AS levels in ripening fruit, the transcripts of these SR genes were activated during fruit maturation in both climacteric and non-climacteric melon varieties. We also found that most SR genes were under selection during domestication. These results represent a novel finding of increased AS levels and SR gene expression during fruit ripening, indicating that alternative splicing may play a role in fruit maturation.

High Overexpression of SiAAP9 Leads to Growth Inhibition and Protein Ectopic Localization in Transgenic Arabidopsis

Amino acid permeases (AAPs) transporters are crucial for the long-distance transport of amino acids in plants, from source to sink. While Arabidopsis and rice have been extensively studied, research on foxtail millet is limited. This study identified two transcripts of SiAAP9, both of which were induced by NO3- and showed similar expression patterns. The overexpression of SiAAP9L and SiAAP9S in Arabidopsis inhibited plant growth and seed size, although SiAAP9 was found to transport more amino acids into seeds. Furthermore, SiAAP9-OX transgenic Arabidopsis showed increased tolerance to high concentrations of glutamate (Glu) and histidine (His). The high overexpression level of SiAAP9 suggested its protein was not only located on the plasma membrane but potentially on other organelles, as well. Interestingly, sequence deletion reduced SiAAP9's sensitivity to Brefeldin A (BFA), and SiAAP9 had ectopic localization on the endoplasmic reticulum (ER). Protoplast amino acid uptake experiments indicated that SiAAP9 enhanced Glu transport into foxtail millet cells. Overall, the two transcripts of SiAAP9 have similar functions, but SiAAP9L shows a higher colocalization with BFA compartments compared to SiAAP9S. Our research identifies a potential candidate gene for enhancing the nutritional quality of foxtail millet through breeding.

Comprehensive deciphering the alternative splicing patterns involved in leaf morphogenesis of Liriodendron chinense

Alternative splicing (AS), a pivotal post-transcriptional regulatory mechanism, profoundly amplifies diversity and complexity of transcriptome and proteome. Liriodendron chinense (Hemsl.) Sarg., an excellent ornamental tree species renowned for its distinctive leaf shape, which resembles the mandarin jacket. Despite the documented potential genes related to leaf development of L. chinense, the underlying post-transcriptional regulatory mechanisms remain veiled. Here, we conducted a comprehensive analysis of the transcriptome to clarify the genome-wide landscape of the AS pattern and the spectrum of spliced isoforms during leaf developmental stages in L. chinense. Our investigation unveiled 50,259 AS events, involving 10,685 genes (32.9%), with intron retention as the most prevalent events. Notably, the initial stage of leaf development witnessed the detection of 804 differentially AS events affiliated with 548 genes. Although both differentially alternative splicing genes (DASGs) and differentially expressed genes (DEGs) were enriched into morphogenetic related pathways during the transition from fishhook (P2) to lobed (P7) leaves, there was only a modest degree of overlap between DASGs and DEGs. Furthermore, we conducted a comprehensively AS analysis on homologous genes involved in leaf morphogenesis, and most of which are subject to post-transcriptional regulation of AS. Among them, the AINTEGUMENTA-LIKE transcript factor LcAIL5 was characterization in detailed, which experiences skipping exon (SE), and two transcripts displayed disparate expression patterns across multiple stages. Overall, these findings yield a comprehensive understanding of leaf development regulation via AS, offering a novel perspective for further deciphering the mechanism of plant leaf morphogenesis.

Profiling of the gene expression and alternative splicing landscapes of Eucalyptus grandis

Eucalyptus is a widely planted hardwood tree species due to its fast growth, superior wood properties and adaptability. However, the post-transcriptional regulatory mechanisms controlling tissue development and stress responses in Eucalyptus remain poorly understood. In this study, we performed a comprehensive analysis of the gene expression profile and the alternative splicing (AS) landscape of E. grandis using strand-specific RNA-Seq, which encompassed 201 libraries including different organs, developmental stages, and environmental stresses. We identified 10 416 genes (33.49%) that underwent AS, and numerous differentially expressed and/or differential AS genes involved in critical biological processes, such as primary-to-secondary growth transition of stems, adventitious root formation, aging and responses to phosphorus- or boron-deficiency. Co-expression analysis of AS events and gene expression patterns highlighted the potential upstream regulatory role of AS events in multiple processes. Additionally, we highlighted the lignin biosynthetic pathway to showcase the potential regulatory functions of AS events in the KNAT3 and IRL3 genes within this pathway. Our high-quality expression atlas and AS landscape serve as valuable resources for unravelling the genetic control of woody plant development, long-term adaptation, and understanding transcriptional diversity in Eucalyptus. Researchers can conveniently access these resources through the interactive ePlant browser (https://bar.utoronto.ca/eplant_eucalyptus).

Gene expression and alternative splicing contribute to adaptive divergence of ecotypes

Regulation of gene expression is a critical link between genotype and phenotype explaining substantial heritable variation within species. However, we are only beginning to understand the ways that specific gene regulatory mechanisms contribute to adaptive divergence of populations. In plants, the post-transcriptional regulatory mechanism of alternative splicing (AS) plays an important role in both development and abiotic stress response, making it a compelling potential target of natural selection. AS allows organisms to generate multiple different transcripts/proteins from a single gene and thus may provide a source of evolutionary novelty. Here, we examine whether variation in alternative splicing and gene expression levels might contribute to adaptation and incipient speciation of dune-adapted prairie sunflowers in Great Sand Dunes National Park, Colorado, USA. We conducted a common garden experiment to assess transcriptomic variation among ecotypes and analyzed differential expression, differential splicing, and gene coexpression. We show that individual genes are strongly differentiated for both transcript level and alternative isoform proportions, even when grown in a common environment, and that gene coexpression networks are disrupted between ecotypes. Furthermore, we examined how genome-wide patterns of sequence divergence correspond to divergence in transcript levels and isoform proportions and find evidence for both cis and trans-regulation. Together, our results emphasize that alternative splicing has been an underappreciated mechanism providing source material for natural selection at short evolutionary time scales.

Transcriptomic analysis of alternative splicing events for different stages of growth and development in Sistan Yaghooti grape clusters

Sistan Yaghooti grape variety, despite characteristics such as early ripening, is vulnerable to cluster rot due to small berries and dense clusters. In this regard, AS may serve as a regulatory mechanism during developmental processes and in response to environmental signals. RNA-Seq analysis was performed to measure gene expression and the extent of AS events in the cluster growth and development stages of Sistan Yaghooti grape. The number of AS events increased during stages, suggesting that it contributes to the grapevine's adaptability to various stresses. In addition, DEG and DAS genes showed little overlap in cluster growth stages. Functional analysis of 19,194 DAS -gene sets showed that VIT_06s0004g06670 gene is involved in the activation of calcium channels (Ca2+) through the activation of 5 PLC biosynthetic pathways. Among the 27,229 DEG -sets, VIT_07s0005g05320 gene showed higher expression. Interestingly, this gene is involved in the synthesis of an EF -hand domain-containing protein capable of binding to Ca2+ by activating 4 biochemical pathways. These genes increase cytosolic Ca2+ concentration, enhancing plant stress tolerance and resistance to cracking. These results show that AS can respond independently to different types of stress. Among the other DAS genes, the GA2ox gene (VvGA2ox) showed an increase in AS events during cluster development. This gene is critical for initiating the degradation process of GA and plays a crucial role in different stages of seed development. Therefore, it is very likely that this gene is one of the main factors responsible for the density and seedlessness of Sistan Yaghooti grape.

Nonsense-mediated mRNA decay pathway in plants under stress: general gene regulatory mechanism and advances

MAIN CONCLUSION

Nonsense-mediated mRNA decay in eukaryotes is vital to cellular homeostasis. Further knowledge of its putative role in plant RNA metabolism under stress is pivotal to developing fitness-optimizing strategies. Nonsense-mediated mRNA decay (NMD), part of the mRNA surveillance pathway, is an evolutionarily conserved form of gene regulation in all living organisms. Degradation of mRNA-bearing premature termination codons and regulation of physiological RNA levels highlight NMD's role in shaping the cellular transcriptome. Initially regarded as purely a tool for cellular RNA quality control, NMD is now considered to mediate various aspects of plant developmental processes and responses to environmental changes. Here we offer a basic understanding of NMD in eukaryotes by explaining the concept of premature termination codon recognition and NMD complex formation. We also provide a detailed overview of the NMD mechanism and its role in gene regulation. The potential role of effectors, including ABCE1, in ribosome recycling during the translation process is also explained. Recent reports of alternatively spliced variants of corresponding genes targeted by NMD in Arabidopsis thaliana are provided in tabular format. Detailed figures are also provided to clarify the NMD concept in plants. In particular, accumulating evidence shows that NMD can serve as a novel alternative strategy for genetic manipulation and can help design RNA-based therapies to combat stress in plants. A key point of emphasis is its function as a gene regulatory mechanism as well as its dynamic regulation by environmental and developmental factors. Overall, a detailed molecular understanding of the NMD mechanism can lead to further diverse applications, such as improving cellular homeostasis in living organisms.

A SNP variation in the Sucrose synthase (SoSUS) gene associated with sugar-related traits in sugarcane

Background

Sugarcane (Saccharum spp.) is an economically significant crop for both the sugar and biofuel industries. Breeding sugarcane cultivars with high-performance agronomic traits is the most effective approach for meeting the rising demand for sugar and biofuels. Molecular markers associated with relevant agronomic traits could drastically reduce the time and resources required to develop new sugarcane varieties. Previous sugarcane candidate gene association analyses have found single nucleotide polymorphism (SNP) markers associated with sugar-related traits. This study aims to validate these associated SNP markers of six genes, including Lesion simulating disease 1 (LSD), Calreticulin (CALR), Sucrose synthase 1 (SUS1), DEAD-box ATP-dependent RNA helicase (RH), KANADI1 (KAN1), and Sodium/hydrogen exchanger 7 (NHX7), in a diverse population in 2-year and two-location evaluations.

Methods

After genotyping of seven targeted SNP markers was performed by PCR Allelic Competitive Extension (PACE) SNP genotyping, the association with sugar-related traits and important cane yield component traits was determined on a set of 159 sugarcane genotypes. The marker-trait relationships were validated and identified by both t-test analysis and an association analysis based on the general linear model.

Results

The mSoSUS1_SNPCh10.T/C and mSoKAN1_SNPCh7.T/C markers that were designed from the SUS1 and KAN1 genes, respectively, showed significant associations with different amounts of sugar-related traits and yield components. The mSoSUS1_SNPCh10.T/C marker was found to have more significant association with sugar-related traits, including pol, CCS, brix, fiber and sugar yield, with p values of 6.08 × 10-6 to 4.35 × 10-2, as well as some cane yield component traits with p values of 1.61 × 10-4 to 3.35 × 10-2. The significant association is consistent across four environments.

Conclusion

Sucrose synthase (SUS) is considered a crucial enzyme involved in sucrose metabolism. This marker is a high potential functional marker that may be used in sugarcane breeding programs to select superior sugarcane with good fiber and high sugar contents.

Identification of cell-type specific alternative transcripts in the multicellular alga Volvox carteri

BACKGROUND

Cell type specialization is a hallmark of complex multicellular organisms and is usually established through implementation of cell-type-specific gene expression programs. The multicellular green alga Volvox carteri has just two cell types, germ and soma, that have previously been shown to have very different transcriptome compositions which match their specialized roles. Here we interrogated another potential mechanism for differentiation in V. carteri, cell type specific alternative transcript isoforms (CTSAI).

METHODS

We used pre-existing predictions of alternative transcripts and de novo transcript assembly with HISAT2 and Ballgown software to compile a list of loci with two or more transcript isoforms, identified a small subset that were candidates for CTSAI, and manually curated this subset of genes to remove false positives. We experimentally verified three candidates using semi-quantitative RT-PCR to assess relative isoform abundance in each cell type.

RESULTS

Of the 1978 loci with two or more predicted transcript isoforms 67 of these also showed cell type isoform expression biases. After curation 15 strong candidates for CTSAI were identified, three of which were experimentally verified, and their predicted gene product functions were evaluated in light of potential cell type specific roles. A comparison of genes with predicted alternative splicing from Chlamydomonas reinhardtii, a unicellular relative of V. carteri, identified little overlap between ortholog pairs with alternative splicing in both species. Finally, we interrogated cell type expression patterns of 126 V. carteri predicted RNA binding protein (RBP) encoding genes and found 40 that showed either somatic or germ cell expression bias. These RBPs are potential mediators of CTSAI in V. carteri and suggest possible pre-adaptation for cell type specific RNA processing and a potential path for generating CTSAI in the early ancestors of metazoans and plants.

CONCLUSIONS

We predicted numerous instances of alternative transcript isoforms in Volvox, only a small subset of which showed cell type specific isoform expression bias. However, the validated examples of CTSAI supported existing hypotheses about cell type specialization in V. carteri, and also suggested new hypotheses about mechanisms of functional specialization for their gene products. Our data imply that CTSAI operates as a minor but important component of V. carteri cellular differentiation and could be used as a model for how alternative isoforms emerge and co-evolve with cell type specialization.

Advances in alternative splicing identification: deep learning and pantranscriptome

In plants, alternative splicing is a crucial mechanism for regulating gene expression at the post-transcriptional level, which leads to diverse proteins by generating multiple mature mRNA isoforms and diversify the gene regulation. Due to the complexity and variability of this process, accurate identification of splicing events is a vital step in studying alternative splicing. This article presents the application of alternative splicing algorithms with or without reference genomes in plants, as well as the integration of advanced deep learning techniques for improved detection accuracy. In addition, we also discuss alternative splicing studies in the pan-genomic background and the usefulness of integrated strategies for fully profiling alternative splicing.

Eukaryotic splicing machinery in the plant-virus battleground

Plant virual infections are mainly caused by plant-virus parasitism which affects ecological communities. Some viruses are highly pathogen specific that can infect only specific plants, while some can cause widespread harm, such as tobacco mosaic virus (TMV) and cucumber mosaic virus (CMV). After a virus infects the host, undergoes a series of harmful effects, including the destruction of host cell membrane receptors, changes in cell membrane components, cell fusion, and the production of neoantigens on the cell surface. Therefore, competition between the host and the virus arises. The virus starts gaining control of critical cellular functions of the host cells and ultimately affects the fate of the targeted host plants. Among these critical cellular processes, alternative splicing (AS) is an essential posttranscriptional regulation process in RNA maturation, which amplify host protein diversity and manipulates transcript abundance in response to plant pathogens. AS is widespread in nearly all human genes and critical in regulating animal-virus interactions. In particular, an animal virus can hijack the host splicing machinery to re-organize its compartments for propagation. Changes in AS are known to cause human disease, and various AS events have been reported to regulate tissue specificity, development, tumour proliferation, and multi-functionality. However, the mechanisms underlying plant-virus interactions are poorly understood. Here, we summarize the current understanding of how viruses interact with their plant hosts compared with humans, analyze currently used and putative candidate agrochemicals to treat plant-viral infections, and finally discussed the potential research hotspots in the future. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Processing > Splicing Regulation/Alternative Splicing.

Fluid protein fold space and its implications

Fold-switching proteins, which remodel their secondary and tertiary structures in response to cellular stimuli, suggest a new view of protein fold space. For decades, experimental evidence has indicated that protein fold space is discrete: dissimilar folds are encoded by dissimilar amino acid sequences. Challenging this assumption, fold-switching proteins interconnect discrete groups of dissimilar protein folds, making protein fold space fluid. Three recent observations support the concept of fluid fold space: (1) some amino acid sequences interconvert between folds with distinct secondary structures, (2) some naturally occurring sequences have switched folds by stepwise mutation, and (3) fold switching is evolutionarily selected and likely confers advantage. These observations indicate that minor amino acid sequence modifications can transform protein structure and function. Consequently, proteomic structural and functional diversity may be expanded by alternative splicing, small nucleotide polymorphisms, post-translational modifications, and modified translation rates.

Toward a systems view on RNA-binding proteins and associated RNAs in plants: Guilt by association

RNA-binding proteins (RBPs) have a broad impact on most biochemical, physiological, and developmental processes in a plant's life. RBPs engage in an on-off relationship with their RNA partners, accompanying virtually every stage in RNA processing and function. While the function of a plethora of RBPs in plant development and stress responses has been described, we are lacking a systems-level understanding of components in RNA-based regulation. Novel techniques have substantially enlarged the compendium of proteins with experimental evidence for binding to RNAs in the cell, the RNA-binding proteome. Furthermore, ribonomics methods have been adapted for use in plants to profile the in vivo binding repertoire of RBPs genome-wide. Here, we discuss how recent technological achievements have provided novel insights into the mode of action of plant RBPs at a genome-wide scale. Furthermore, we touch upon two emerging topics, the connection of RBPs to phase separation in the cell and to extracellular RNAs. Finally, we define open questions to be addressed to move toward an integrated understanding of RBP function.

Catsnap: a user-friendly algorithm for determining the conservation of protein variants reveals extensive parallelisms in the evolution of alternative splicing

Understanding the evolutionary conservation of complex eukaryotic transcriptomes significantly illuminates the physiological relevance of alternative splicing (AS). Examining the evolutionary depth of a given AS event with ordinary homology searches is generally challenging and time-consuming. Here, we present Catsnap, an algorithmic pipeline for assessing the conservation of putative protein isoforms generated by AS. It employs a machine learning approach following a database search with the provided pair of protein sequences. We used the Catsnap algorithm for analyzing the conservation of emerging experimentally characterized alternative proteins from plants and animals. Indeed, most of them are conserved among other species. Catsnap can detect the conserved functional protein isoforms regardless of the AS type by which they are generated. Notably, we found that while the primary amino acid sequence is maintained, the type of AS determining the inclusion or exclusion of protein regions varies throughout plant phylogenetic lineages in these proteins. We also document that this phenomenon is less seen among animals. In sum, our algorithm highlights the presence of unexpectedly frequent hotspots where protein isoforms recurrently arise to carry physiologically relevant functions. The user web interface is available at https://catsnap.cesnet.cz/.

The physiology of alternative splicing

Alternative splicing is a substantial contributor to the high complexity of transcriptomes of multicellular eukaryotes. In this Review, we discuss the accumulated evidence that most of this complexity is reflected at the protein level and fundamentally shapes the physiology and pathology of organisms. This notion is supported not only by genome-wide analyses but, mainly, by detailed studies showing that global and gene-specific modulations of alternative splicing regulate highly diverse processes such as tissue-specific and species-specific cell differentiation, thermal regulation, neuron self-avoidance, infrared sensing, the Warburg effect, maintenance of telomere length, cancer and autism spectrum disorders (ASD). We also discuss how mastering the control of alternative splicing paved the way to clinically approved therapies for hereditary diseases.

Full-length transcriptome combined with RNA sequence analysis of Fraxinus chinensis

BACKGROUND

The dry root or stem bark of Fraxinus chinensis is a famous herb Qin Pi which is known for its anti-inflammatory, analgesic, anti-tumor, liver protective and diuretic pharmacological effects, the fundamental chemical components are coumarin, phenylethanol glycosides and flavonoids. However, it is difficult to clarify the secondary metabolite synthesis pathway and key genes involved in the pathway because of lack genome information of Fraxinus chinensis.

OBJECTIVE

To generate a complete transcriptome of Fraxinus chinensis and to clarify the differentially expressed genes (DEGs) in leaves and stem barks.

METHODS

In this study, full-length transcriptome analysis and RNA-Seq were combined to characterize Fraxinus chinensis transcriptome.

RESULTS

A total of 69,145 transcripts were acquired and regarded as reference transcriptome, 67,441 transcripts (97.47%) were annotated to NCBI non-redundant protein (Nr), SwissProt, the Kyoto Encyclopedia of Genes and Genomes (KEGG) and eukaryotic orthologous groups (KOG) databases. A total of 18,917 isoforms were annotated to KEGG database and classified to 138 biological pathways. In total, 10,822 simple sequence repeat (SSRs) and 11,319 resistance (R) gene were classified to 18 types, and 3947 transcription factors (TFs) were identified in full-length transcriptome analysis. Additionally, 15,095 DEGs were detected by RNA-seq in leaves and barks, including 4696 significantly up-regulated and 10,399 significantly down-regulated genes. And 254 transcripts were annotated into phenylpropane metabolism pathway containing 86 DEGs and ten of these enzyme genes were verified by qRT-PCR.

CONCLUSION

It laid the foundation for further exploration of the biosynthetic pathway of phenylpropanoids and related key enzyme genes.

Mutation of Arabidopsis SME1 and Sm core assembly improves oxidative stress resilience

Alternative splicing is a key posttranscriptional gene regulatory process, acting in diverse adaptive and basal plant processes. Splicing of precursor-messenger RNA (pre-mRNA) is catalyzed by a dynamic ribonucleoprotein complex, designated the spliceosome. In a suppressor screen, we identified a nonsense mutation in the Smith (Sm) antigen protein SME1 to alleviate photorespiratory H₂O₂-dependent cell death in catalase deficient plants. Similar attenuation of cell death was observed upon chemical inhibition of the spliceosome, suggesting pre-mRNA splicing inhibition to be responsible for the observed cell death alleviation. Furthermore, the sme1-2 mutants showed increased tolerance to the reactive oxygen species inducing herbicide methyl viologen. Both an mRNA-seq and shotgun proteomic analysis in sme1-2 mutants displayed a constitutive molecular stress response, together with extensive alterations in pre-mRNA splicing of transcripts encoding metabolic enzymes and RNA binding proteins, even under unstressed conditions. Using SME1 as a bait to identify protein interactors, we provide experimental evidence for almost 50 homologs of the mammalian spliceosome-associated protein to reside in the Arabidopsis thaliana spliceosome complexes and propose roles in pre-mRNA splicing for four uncharacterized plant proteins. Furthermore, as for sme1-2, a mutant in the Sm core assembly protein ICLN resulted in a decreased sensitivity to methyl viologen. Taken together, these data show that both a perturbed Sm core composition and assembly results in the activation of a defense response and in enhanced resilience to oxidative stress.

A splice site mutation in the FvePHP gene is associated with leaf development and flowering time in woodland strawberry

Leaves and flowers are crucial for the growth and development of higher plants. In this study we identified a mutant with narrow leaflets and early flowering (nlef) in an ethyl methanesulfonate-mutagenized population of woodland strawberry (Fragaria vesca) and aimed to identify the candidate gene. Genetic analysis revealed that a single recessive gene, nlef, controlled the mutant phenotype. We found that FvH4_1g25470, which encodes a putative DNA polymerase α with a polymerase and histidinol phosphatase domain (PHP), might be the candidate gene, using bulked segregant analysis with whole-genome sequencing, molecular markers, and cloning analyses. A splice donor site mutation (C to T) at the 5' end of the second intron led to an erroneous splice event that reduced the expression level of the full-length transcript of FvePHP in mutant plants. FvePHP was localized in the nucleus and was highly expressed in leaves. Silencing of FvePHP using the virus-induced gene silencing method resulted in partial developmental defects in strawberry leaves. Overexpression of the FvePHP gene can largely restore the mutant phenotype. The expression levels of FveSEP1, FveSEP3, FveAP1, FveFUL, and FveFT were higher in the mutants than those in 'Yellow Wonder' plants, probably contributing to the early flowering phenotype in mutant plants. Our results indicate that mutation in FvePHP is associated with multiple developmental pathways. These results aid in understanding the role of DNA polymerase in strawberry development.

A divide-and-conquer approach for genomic prediction in rubber tree using machine learning

Rubber tree (Hevea brasiliensis) is the main feedstock for commercial rubber; however, its long vegetative cycle has hindered the development of more productive varieties via breeding programs. With the availability of H. brasiliensis genomic data, several linkage maps with associated quantitative trait loci have been constructed and suggested as a tool for marker-assisted selection. Nonetheless, novel genomic strategies are still needed, and genomic selection (GS) may facilitate rubber tree breeding programs aimed at reducing the required cycles for performance assessment. Even though such a methodology has already been shown to be a promising tool for rubber tree breeding, increased model predictive capabilities and practical application are still needed. Here, we developed a novel machine learning-based approach for predicting rubber tree stem circumference based on molecular markers. Through a divide-and-conquer strategy, we propose a neural network prediction system with two stages: (1) subpopulation prediction and (2) phenotype estimation. This approach yielded higher accuracies than traditional statistical models in a single-environment scenario. By delivering large accuracy improvements, our methodology represents a powerful tool for use in Hevea GS strategies. Therefore, the incorporation of machine learning techniques into rubber tree GS represents an opportunity to build more robust models and optimize Hevea breeding programs.

Genome-scale analysis of Arabidopsis splicing-related protein kinase families reveals roles in abiotic stress adaptation

Nearly 60 - 80 % of intron-containing plant genes undergo alternative splicing in response to either stress or plant developmental cues. RNA splicing is performed by a large ribonucleoprotein complex called the spliceosome in conjunction with associated subunits such as serine arginine (SR) proteins, all of which undergo extensive phosphorylation. In plants, there are three main protein kinase families suggested to phosphorylate core spliceosome subunits and related splicing factors based on orthology to human splicing-related kinases: the SERINE/ARGININE PROTEIN KINASES (SRPK), ARABIDOPSIS FUS3 COMPLEMENT (AFC), and Pre-mRNA PROCESSING FACTOR 4 (PRP4K) protein kinases. To better define the conservation and role(s) of these kinases in plants, we performed a genome-scale analysis of the three families across photosynthetic eukaryotes, followed by extensive transcriptomic and bioinformatic analysis of all Arabidopsis thaliana SRPK, AFC, and PRP4K protein kinases to elucidate their biological functions. Unexpectedly, this revealed the existence of SRPK and AFC phylogenetic groups with distinct promoter elements and patterns of transcriptional response to abiotic stress, while PRP4Ks possess no phylogenetic sub-divisions, suggestive of functional redundancy. We also reveal splicing-related kinase families are both diel and photoperiod regulated, implicating different orthologs as discrete time-of-day RNA splicing regulators. This foundational work establishes a number of new hypotheses regarding how reversible spliceosome phosphorylation contributes to both diel plant cell regulation and abiotic stress adaptation in plants.

A crosstalk of circadian clock and alternative splicing under abiotic stresses in the plants

The circadian clock is an internal time-keeping mechanism that synchronizes the physiological adaptation of an organism to its surroundings based on day and night transition in a period of 24 h, suggesting the circadian clock provides fitness by adjusting environmental constrains. The circadian clock is driven by positive and negative elements that regulate transcriptionally and post-transcriptionally. Alternative splicing (AS) is a crucial transcriptional regulator capable of generating large numbers of mRNA transcripts from limited numbers of genes, leading to proteome diversity, which is involved in circadian to deal with abiotic stresses. Over the past decade, AS and circadian control have been suggested to coordinately regulate plant performance under fluctuating environmental conditions. However, only a few reports have reported the regulatory mechanism of this complex crosstalk. Based on the emerging evidence, this review elaborates on the existing links between circadian and AS in response to abiotic stresses, suggesting an uncovered regulatory network among circadian, AS, and abiotic stresses. Therefore, the rhythmically expressed splicing factors and core clock oscillators fill the role of temporal regulators participating in improving plant growth, development, and increasing plant tolerance against abiotic stresses.

An Arabidopsis Retention and Splicing complex regulates root and embryo development through pre-mRNA splicing

Pre-mRNA splicing is an important step in the posttranscriptional processing of transcripts and a key regulator of development. The heterotrimeric retention and splicing (RES) complex plays vital roles in the growth and development of yeast, zebrafish, and humans by mediating pre-mRNA splicing of multiple genes. However, whether the RES complex is conserved in plants and what specific functions it has remain unknown. In this study, we identified Arabidopsis (Arabidopsis thaliana) BUD13 (AtBUD13), GROWTH, DEVELOPMENT AND SPLICING 1 (GDS1), and DAWDLE (DDL) as the counterparts of the yeast RES complex subunits Bud site selection protein 13 (Bud13), U2 snRNP component Snu17 (Snu17), and Pre-mRNA leakage protein 1, respectively. Moreover, we showed that RES is an ancient complex evolutionarily conserved in eukaryotes. GDS1 directly interacts with both AtBUD13 and DDL in nuclear speckles. The BUD13 domain of AtBUD13 and the RNA recognition motif domain of GDS1 are necessary and sufficient for AtBUD13-GDS1 interaction. Mutants of AtBUD13, GDS1, and DDL failed to properly splice multiple genes involved in cell proliferation and showed defects in early embryogenesis and root development. In addition, we found that GDS1 and DDL interact, respectively, with the U2 small nuclear ribonucleoproteins auxiliary factor AtU2AF65B and the NineTeen Complex-related splicing factor SKIP, which are essential for early steps of spliceosome assembly and recognition of splice sites. Altogether, our work reveals that the Arabidopsis RES complex is important for root and early embryo development by modulating pre-mRNA splicing.

Genome-wide identification, characterization of Serotonin N-acetyltransferase and deciphering its importance under development, biotic and abiotic stress in soybean

Serotonin N-acetyltransferase (SNAT) is the penultimate enzyme involved in plant melatonin biosynthesis. Identifying its expression under development and stress will reveal the regulatory role in the soybean. To identify and characterize SNAT, we employed genome-wide analysis, gene structure, cis-acting elements, expression, and enzyme activity. We identified seven putative genes by genome-wide analysis and found chloroplast signal peptides in three GmSNATs. To elucidate GmSNATs role, expression datasets of more than a hundred samples related to circadian rhythm, developmental stages, and stress conditions were analysed. Notably, the expression of GmSNAT1 did not show significant expression during biotic and abiotic stress. The GmSNAT1 sequence showed 67.8 and 72.2 % similarities with OsSNAT and AtSNAT, respectively. The Km and Vmax of the purified recombinant GmSNAT1 were 657 μM and 3780 pmol/min/mg, respectively. To further understand the GmSNAT1 role, we supplemented different concentrations of serotonin and melatonin to in-vitro cultures and seed priming. These studies revealed that the GmSNAT1 expression was significantly up-regulated at higher concentrations of serotonin and down-regulated at higher melatonin concentrations. We speculate that a high concentration of melatonin during abiotic, biotic stress, and in-vitro cultures are responsible for regulating GmSNAT1 expression, which may regulate them at the enzyme level during stress in soybean.

Identification and analysis of HD-Zip genes involved in the leaf development of Liriodendron chinense using multidimensional analysis

Homeodomain-leucine zipper (HD-Zip) proteins are plant-specific transcription factors that play important roles in different biological processes, especially leaf development. However, no studies to date have identified the HD-Zip genes in Liriodendron chinense nor characterized their functions. We identified the HD-Zip genes in L. chinense by analysing the phylogeny, chromosome location, structure, conserved motif, cis-regulatory elements, synteny, post-transcriptional regulation and expression patterns of these genes during leaf development. A total of 36 LcHD-Zip genes were identified and divided into four subfamilies (HD-Zip I to IV). Synteny analysis revealed that segmental duplication was the main force driving the expansion of LcHD-Zip genes. These 36 LcHD-Zip genes exhibited 11 different expression patterns. Pattern 1, 2, 3, 4, 6, 7, 8 and 9 genes may play important roles in leaf development, such as leaf initiation, leaf polarity establishment, leaf shape development, phytohormone-mediated leaf growth and leaf epidermal structure formation. Four HD-Zip III genes were targeted by microRNAs (miRNAs), and the miR165/166a-HD-Zip regulatory module formed regulated leaf initiation and leaf polarity establishment. Overall, LcHD-Zip genes play key roles in leaf development of L. chinense. This work provides a foundation for the functional verification of HD-Zip genes identified in this study.

Alternative Splicing and Its Roles in Plant Metabolism

Plant metabolism, including primary metabolism such as tricarboxylic acid cycle, glycolysis, shikimate and amino acid pathways as well as specialized metabolism such as biosynthesis of phenolics, alkaloids and saponins, contributes to plant survival, growth, development and interactions with the environment. To this end, these metabolic processes are tightly and finely regulated transcriptionally, post-transcriptionally, translationally and post-translationally in response to different growth and developmental stages as well as the constantly changing environment. In this review, we summarize and describe the current knowledge of the regulation of plant metabolism by alternative splicing, a post-transcriptional regulatory mechanism that generates multiple protein isoforms from a single gene by using alternative splice sites during splicing. Numerous genes in plant metabolism have been shown to be alternatively spliced under different developmental stages and stress conditions. In particular, alternative splicing serves as a regulatory mechanism to fine-tune plant metabolism by altering biochemical activities, interaction and subcellular localization of proteins encoded by splice isoforms of various genes.

The Function of DNA Demethylase Gene ROS1a Null Mutant on Seed Development in Rice (Oryza Sativa) Using the CRISPR/CAS9 System

The endosperm is the main nutrient source in cereals for humans, as it is a highly specialized storage organ for starch, lipids, and proteins, and plays an essential role in seed growth and development. Active DNA demethylation regulates plant developmental processes and is ensured by cytosine methylation (5-meC) DNA glycosylase enzymes. To find out the role of OsROS1a in seed development, the null mutant of OsROS1a was generated using the CRISPR/Cas9 system. The null mutant of OsROS1a was stable and heritable, which affects the major agronomic traits, particularly in rice seeds. The null mutant of OsROS1a showed longer and narrower grains, and seeds were deformed containing an underdeveloped and less-starch-producing endosperm with slightly irregularly shaped embryos. In contrast to the transparent grains of the wild type, the grains of the null mutant of OsROS1a were slightly opaque and rounded starch granules, with uneven shapes, sizes, and surfaces. A total of 723 differential expression genes (DEGs) were detected in the null mutant of OsROS1a by RNA-Seq, of which 290 were downregulated and 433 were upregulated. The gene ontology (GO) terms with the top 20 enrichment factors were visualized for cellular components, biological processes, and molecular functions. The key genes that are enriched for these GO terms include starch synthesis genes (OsSSIIa and OsSSIIIa) and cellulose synthesis genes (CESA2, CESA3, CESA6, and CESA8). Genes encoding polysaccharides and glutelin were found to be downregulated in the mutant endosperm. The glutelins were further verified by SDS-PAGE, suggesting that glutelin genes could be involved in the null mutant of OsROS1a seed phenotype and OsROS1a could have the key role in the regulation of glutelins. Furthermore, 378 differentially alternative splicing (AS) genes were identified in the null mutant of OsROS1a, suggesting that the OsROS1a gene has an impact on AS events. Our findings indicated that the function on rice endosperm development in the null mutant of OsROS1a could be influenced through regulating gene expression and AS, which could provide the base to properly understand the molecular mechanism related to the OsROS1a gene in the regulation of rice seed development.

The Rice Serine/Arginine Splicing Factor RS33 Regulates Pre-mRNA Splicing during Abiotic Stress Responses

Abiotic stresses profoundly affect plant growth and development and limit crop productivity. Pre-mRNA splicing is a major form of gene regulation that helps plants cope with various stresses. Serine/arginine (SR)-rich splicing factors play a key role in pre-mRNA splicing to regulate different biological processes under stress conditions. Alternative splicing (AS) of SR transcripts and other transcripts of stress-responsive genes generates multiple splice isoforms that contribute to protein diversity, modulate gene expression, and affect plant stress tolerance. Here, we investigated the function of the plant-specific SR protein RS33 in regulating pre-mRNA splicing and abiotic stress responses in rice. The loss-of-function mutant rs33 showed increased sensitivity to salt and low-temperature stresses. Genome-wide analyses of gene expression and splicing in wild-type and rs33 seedlings subjected to these stresses identified multiple splice isoforms of stress-responsive genes whose AS are regulated by RS33. The number of RS33-regulated genes was much higher under low-temperature stress than under salt stress. Our results suggest that the plant-specific splicing factor RS33 plays a crucial role during plant responses to abiotic stresses.

Ubiquitination and Ubiquitin-Like Modifications as Mediators of Alternative Pre-mRNA Splicing in Arabidopsis thaliana

Alternative splicing (AS) is a common post-transcriptional regulatory process in eukaryotes. AS has an irreplaceable role during plant development and in response to environmental stress as it evokes differential expression of downstream genes or splicing factors (e.g., serine/arginine-rich proteins). Numerous studies have reported that loss of AS capacity leads to defects in plant growth and development, and induction of stress-sensitive phenotypes. A role for post-translational modification (PTM) of AS components has emerged in recent years. These modifications are capable of regulating the activity, stability, localization, interaction, and folding of spliceosomal proteins in human cells and yeast, indicating that PTMs represent another layer of AS regulation. In this review, we summarize the recent reports concerning ubiquitin and ubiquitin-like modification of spliceosome components and analyze the relationship between spliceosome and the ubiquitin/26S proteasome pathway in plants. Based on the totality of the evidence presented, we further speculate on the roles of protein ubiquitination mediated AS in plant development and environmental response.

Single-molecule real-time sequencing of the full-length transcriptome of purple garlic (Allium sativum L. cv. Leduzipi) and identification of serine O-acetyltransferase family proteins involved in cysteine biosynthesis

BACKGROUND

Garlic (Allium sativum L.), whose bioactive components are mainly organosulfur compounds (OSCs), is a herbaceous perennial widely consumed as a green vegetable and a condiment. Yet, the metabolic enzymes involved in the biosynthesis of OSCs are not identified in garlic.

RESULTS

Here, a full-length transcriptome of purple garlic was generated via PacBio and Illumina sequencing, to characterize the garlic transcriptome and identify key proteins mediating the biosynthesis of OSCs. Overall, 22.56 Gb of clean data were generated, resulting in 454 698 circular consensus sequence (CCS) reads, of which 83.4% (379 206) were identified as being full-length non-chimeric reads - their further transcript clustering facilitated identification of 36 571 high-quality consensus reads. Once corrected, their genome-wide mapping revealed that 6140 reads were novel isoforms of known genes, and 2186 reads were novel isoforms from novel genes. We detected 1677 alternative splicing events, finding 2902 genes possessing either two or more poly(A) sites. Given the importance of serine O-acetyltransferase (SERAT) in cysteine biosynthesis, we investigated the five SERAT homologs in garlic. Phylogenetic analysis revealed a three-tier classification of SERAT proteins, each featuring a serine acetyltransferase domain (N-terminal) and one or two hexapeptide transferase motifs. Template-based modeling showed that garlic SERATs shared a common homo-trimeric structure with homologs from bacteria and other plants. The residues responsible for substrate recognition and catalysis were highly conserved, implying a similar reaction mechanism. In profiling the five SERAT genes' transcript levels, their expression pattern varied significantly among different tissues.

CONCLUSION

This study's findings deepen our knowledge of SERAT proteins, and provide timely genetic resources that could advance future exploration into garlic's genetic improvement and breeding. © 2021 Society of Chemical Industry.

A Comprehensive Identification and Function Analysis of Serine/Arginine-Rich (SR) Proteins in Cotton (Gossypium spp.)

As one of the most important factors in alternative splicing (AS) events, serine/arginine-rich (SR) proteins not only participate in the growth and development of plants but also play pivotal roles in abiotic stresses. However, the research about SR proteins in cotton is still lacking. In this study, we performed an extensive comparative analysis of SR proteins and determined their phylogeny in the plant lineage. A total of 169 SR family members were identified from four Gossypium species, and these genes could be divided into eight distinct subfamilies. The domain, motif distribution and gene structure of cotton SR proteins are conserved within each subfamily. The expansion of SR genes is mainly contributed by WGD and allopolyploidization events in cotton. The selection pressure analysis showed that all the paralogous gene pairs were under purifying selection pressure. Many cis-elements responding to abiotic stress and phytohormones were identified in the upstream sequences of the GhSR genes. Expression profiling suggested that some GhSR genes may involve in the pathways of plant resistance to abiotic stresses. The WGCNA analysis showed that GhSCL-8 co-expressed with many abiotic responding related genes in a salt-responding network. The Y2H assays showed that GhSCL-8 could interact with GhSRs in other subfamilies. The subcellular location analysis showed that GhSCL-8 is expressed in the nucleus. The further VIGS assays showed that the silencing of GhSCL-8 could decrease salt tolerance in cotton. These results expand our knowledge of the evolution of the SR gene family in plants, and they will also contribute to the elucidation of the biological functions of SR genes in the future.

Analysis of Alternative Splicing During the Combinatorial Response to Simultaneous Copper and Iron Deficiency in Arabidopsis Reveals Differential Events in Genes Involved in Amino Acid Metabolism

Copper (Cu) and iron (Fe) constitute fundamental nutrients for plant biology but are often limited due to low bioavailability. Unlike responses to single Cu or Fe deprivation, the consequences of simultaneous Cu and Fe deficiency have not yet been fully deciphered. Previously, it was demonstrated that Cu and Fe deficiency applied in combination imposes transcriptome, proteome, and metabolome changes different from those triggered under each deficiency individually. Here, we evaluated the effect of alternative splicing (AS) on the transcriptome of rosette leaves under single and simultaneous Cu and Fe deficiency. Differentially spliced genes (DSGs) and differentially expressed genes (DEGs) coincided in number (2,600 approx.) although the overlapping fraction was minimal (15%). Functional annotation of changes exclusively detected under simultaneous Cu and Fe deficiency revealed that DEGs participated in general stress responses and translation, while DSGs were involved in metabolic reactions, especially amino acid biosynthesis. Interestingly, transcripts encoding central features for tryptophan (Trp) and asparagine (Asn) synthesis - two significantly altered metabolites under simultaneous Cu and Fe deficiency - underwent exclusive intron retention events under the double deficiency. However, transcript and protein amounts for these enzymes did not correlate with Trp and Asn concentration. In consequence, we propose that AS might act as a regulatory mechanism to modify the stability and/or functionality of the enzymes and therefore fine-tune amino acid production during the combinatorial response to simultaneous Cu and Fe deficiency.

Alternative Splicing Dynamics During the Lifecycle of Salvia miltiorrhiza Root Revealed the Fine Tuning in Root Development and Ingredients Biosynthesis

Alternative splicing (AS) is an essential post-transcriptional process that enhances the coding and regulatory potential of the genome, thereby strongly influencing multiple plant physiology processes, such as metabolic biosynthesis. To explore how AS affects the root development and synthesis of tanshinones and phenolic acid pathways in Salvia miltiorrhiza roots, we investigated the dynamic landscape of AS events in S. miltiorrhiza roots during an annual life history. Temporal profiling represented a distinct temporal variation of AS during the entire development stages, showing the most abundant AS events at the early seedling stage (ES stage) and troughs in 45 days after germination (DAG) and 120 DAG. Gene ontology (GO) analysis indicated that physiological and molecular events, such as lateral root formation, gravity response, RNA splicing regulation, and mitogen-activated protein kinase (MAPK) cascade, were greatly affected by AS at the ES stage. AS events were identified in the tanshinones and phenolic acids pathways as well, especially for the genes for the branch points of the pathways as SmRAS and SmKSL1. Fifteen Ser/Arg-rich (SR) proteins and eight phosphokinases (PKs) were identified with high transcription levels at the ES stage, showing their regulatory roles for the high frequency of AS in this stage. Simultaneously, a co-expression network that includes 521 highly expressed AS genes, SRs, and PKs, provides deeper insight into the mechanism for the variable programming of AS.

Transcriptomic Complexity in Strawberry Fruit Development and Maturation Revealed by Nanopore Sequencing

The use of alternative transcription start or termination sites (aTSS or aTTS) as well as alternative splicing (AS) produce diverse transcript isoforms, playing indispensable roles in the plant development and environmental adaptations. Despite the advances in the finding of the genome-wide alternatively spliced genes in strawberry, it remains unexplored how AS responds to the developmental cues and what relevance do these outcomes have to the gene function. In this study, we have systematically investigated the transcriptome complexity using long-read Oxford Nanopore Technologies along the four successive developmental stages. The full-length cDNA sequencing results unraveled thousands of previously unexplored transcript isoforms raised from aTSS, aTTS, and AS. The relative contributions of these three processes to the complexity of strawberry fruit transcripts were compared. The aTSS and aTTS were more abundant than the AS. Differentially expressed transcripts unraveled the key transitional role of the white fruit stage. Isoform switches of transcripts from 757 genes were observed. They were associated with protein-coding potential change and domain gain or loss as the main consequences. Those genes with switched isoforms take part in the key processes of maturation in the late stages. A case study using yeast two hybrid analysis supported the functional divergence of the two isoforms of the B-box protein 22. Our results provided a new comprehensive overview of the dynamic transcriptomic landscape during strawberry fruit development and maturation.

How alternative splicing changes the properties of plant proteins

Most plant primary transcripts undergo alternative splicing (AS), and its impact on protein diversity is a subject of intensive investigation. Several studies have uncovered various mechanisms of how particular protein splice isoforms operate. However, the common principles behind the AS effects on protein function in plants have rarely been surveyed. Here, on the selected examples, we highlight diverse tissue expression patterns, subcellular localization, enzymatic activities, abilities to bind other molecules and other relevant features. We describe how the protein isoforms mutually interact to underline their intriguing roles in altering the functionality of protein complexes. Moreover, we also discuss the known cases when these interactions have been placed inside the autoregulatory loops. This review is particularly intended for plant cell and developmental biologists who would like to gain inspiration on how the splice variants encoded by their genes of interest may coordinately work.

Mutually opposing activity of PIN7 splicing isoforms is required for auxin-mediated tropic responses in Arabidopsis thaliana

Advanced transcriptome sequencing has revealed that the majority of eukaryotic genes undergo alternative splicing (AS). Nonetheless, little effort has been dedicated to investigating the functional relevance of particular splicing events, even those in the key developmental and hormonal regulators. Combining approaches of genetics, biochemistry and advanced confocal microscopy, we describe the impact of alternative splicing on the PIN7 gene in the model plant Arabidopsis thaliana. PIN7 encodes a polarly localized transporter for the phytohormone auxin and produces two evolutionarily conserved transcripts, PIN7a and PIN7b. PIN7a and PIN7b, differing in a four amino acid stretch, exhibit almost identical expression patterns and subcellular localization. We reveal that they are closely associated and mutually influence each other's mobility within the plasma membrane. Phenotypic complementation tests indicate that the functional contribution of PIN7b per se is minor, but it markedly reduces the prominent PIN7a activity, which is required for correct seedling apical hook formation and auxin-mediated tropic responses. Our results establish alternative splicing of the PIN family as a conserved, functionally relevant mechanism, revealing an additional regulatory level of auxin-mediated plant development.

Structural property, molecular regulation, and functional diversity of glutamine synthetase in higher plants: a data-mining bioinformatics approach

Glutamine synthetase (GS; E.C.6.3.1.2) is a key enzyme in higher plants with two isozymes, cytosolic GS1 and plastidic GS2, and involves in the assimilation and recycling of NH₄⁺ ions and maintenance of complex traits such as crop nitrogen-use efficiency and yield. Our present understanding of crop nitrogen-use efficiency and its correlation with the functional role of the GS family genes is inadequate, which delays harnessing the benefit of this key enzyme in crop improvement. In this report, we performed a comprehensive investigation on the phylogenetic relationship, structural properties, complex multilevel gene regulation, and expression patterns of the GS genes to enrich present understanding about the enzyme. Our Gene Ontology and protein-protein interactions analysis revealed the functional aspects of GS isozymes in stress mitigation, aging, nucleotide biosynthesis/transport, DNA repair and response to metals. The insight gained here contributes to the future research strategies in developing climate-smart crops for global sustainability.

A dynamic intron retention program regulates the expression of several hundred genes during pollen meiosis

Intron retention is a stage-specific mechanism of functional attenuation of a subset of co-regulated, functionally related genes during early stages of pollen development. To improve our understanding of the gene regulatory mechanisms that drive developmental processes, we performed a genome-wide study of alternative splicing and isoform switching during five key stages of pollen development in field mustard, Brassica rapa. Surprisingly, for several hundred genes (12.3% of the genes analysed), isoform switching results in stage-specific expression of intron-retaining transcripts at the meiotic stage of pollen development. In such cases, we report temporally regulated switching between expression of a canonical, translatable isoform and an intron-retaining transcript that is predicted to produce a truncated and presumably inactive protein. The results suggest a new pervasive mechanism underlying modulation of protein levels in a plant developmental program. The effect is not based on gene expression induction but on the type of transcript produced. We conclude that intron retention is a stage-specific mechanism of functional attenuation of a subset of co-regulated, functionally related genes during meiosis, especially genes related to ribosome biogenesis, mRNA transport and nuclear envelope architecture. We also propose that stage-specific expression of a non-functional isoform of Brassica rapa BrSDG8, a non-redundant member of histone methyltransferase gene family, linked to alternative splicing regulation, may contribute to the intron retention observed.

Medicago ABI3 Splicing Isoforms Regulate the Expression of Different Gene Clusters to Orchestrate Seed Maturation

Seed maturation comprises important developmental processes, such as seed filling and the acquisition of seed germination capacity, desiccation tolerance, longevity, and dormancy. The molecular regulation of these processes is tightly controlled by the LAFL transcription factors, among which ABSCISIC ACID INSENSITIVE 3 (ABI3) was shown to be involved in most of these seed maturation processes. Here, we studied the ABI3 gene from Medicago truncatula, a model legume plant for seed studies. With the transcriptomes of two loss-of-function Medicago abi3 mutants, we were able to show that many gene classes were impacted by the abi3 mutation at different stages of early, middle, and late seed maturation. We also discovered three MtABI3 expression isoforms, which present contrasting expression patterns during seed development. Moreover, by ectopically expressing these isoforms in Medicago hairy roots generated from the abi3 mutant line background, we showed that each isoform regulated specific gene clusters, suggesting divergent molecular functions. Furthermore, we complemented the Arabidopsis abi3 mutant with each of the three MtABI3 isoforms and concluded that all isoforms were capable of restoring seed viability and desiccation tolerance phenotypes even if not all isoforms complemented the seed color phenotype. Taken together, our results allow a better understanding of the ABI3 network in Medicago during seed development, as well as the discovery of commonly regulated genes from the three MtABI3 isoforms, which can give us new insights into how desiccation tolerance and seed viability are regulated.

Genome-Wide Identification and Characterization of Small Peptides in Maize

Small peptides (sPeptides), <100 amino acids (aa) long, are encoded by small open reading frames (sORFs) often found in the 5' and 3' untranslated regions (or other parts) of mRNAs, in long non-coding RNAs, or transcripts from introns and intergenic regions; various sPeptides play important roles in multiple biological processes. In this study, we conducted a comprehensive study of maize (Zea mays) sPeptides using mRNA sequencing, ribosome profiling (Ribo-seq), and mass spectrometry (MS) on six tissues (each with at least two replicates). To identify maize sORFs and sPeptides from these data, we set up a robust bioinformatics pipeline and performed a genome-wide scan. This scan uncovered 9,388 sORFs encoding peptides of 2-100 aa. These sORFs showed distinct genomic features, such as different Kozak region sequences, higher specificity of translation, and high translational efficiency, compared with the canonical protein-coding genes. Furthermore, the MS data verified 2,695 sPeptides. These sPeptides perfectly discriminated all the tissues and were highly associated with their parental genes. Interestingly, the parental genes of sPeptides were significantly enriched in multiple functional gene ontology terms related to abiotic stress and development, suggesting the potential roles of sPeptides in the regulation of their parental genes. Overall, this study lays out the guidelines for genome-wide scans of sORFs and sPeptides in plants by integrating Ribo-seq and MS data and provides a more comprehensive resource of functional sPeptides in maize and gives a new perspective on the complex biological systems of plants.

Daily temperature cycles promote alternative splicing of RNAs encoding SR45a, a splicing regulator in maize

Elevated temperatures enhance alternative RNA splicing in maize (Zea mays) with the potential to expand the repertoire of plant responses to heat stress. Alternative RNA splicing generates multiple RNA isoforms for many maize genes, and here we observed changes in the pattern of RNA isoforms with temperature changes. Increases in maximum daily temperature elevated the frequency of the major modes of alternative splices (AS), in particular retained introns and skipped exons. The genes most frequently targeted by increased AS with temperature encode factors involved in RNA processing and plant development. Genes encoding regulators of alternative RNA splicing were themselves among the principal AS targets in maize. Under controlled environmental conditions, daily changes in temperature comparable to field conditions altered the abundance of different RNA isoforms, including the RNAs encoding the splicing regulator SR45a, a member of the SR45 gene family. We established an "in protoplast" RNA splicing assay to show that during the afternoon on simulated hot summer days, SR45a RNA isoforms were produced with the potential to encode proteins efficient in splicing model substrates. With the RNA splicing assay, we also defined the exonic splicing enhancers that the splicing-efficient SR45a forms utilize to aid in the splicing of model substrates. Hence, with rising temperatures on hot summer days, SR45a RNA isoforms in maize are produced with the capability to encode proteins with greater RNA splicing potential.

Quantifying splice-site usage: a simple yet powerful approach to analyze splicing

RNA splicing, and variations in this process referred to as alternative splicing, are critical aspects of gene regulation in eukaryotes. From environmental responses in plants to being a primary link between genetic variation and disease in humans, splicing differences confer extensive phenotypic changes across diverse organisms (1–3). Regulation of splicing occurs through differential selection of splice sites in a splicing reaction, which results in variation in the abundance of isoforms and/or splicing events. However, genomic determinants that influence splice-site selection remain largely unknown. While traditional approaches for analyzing splicing rely on quantifying variant transcripts (i.e. isoforms) or splicing events (i.e. intron retention, exon skipping etc.) (4), recent approaches focus on analyzing complex/mutually exclusive splicing patterns (5–8). However, none of these approaches explicitly measure individual splice-site usage, which can provide valuable information about splice-site choice and its regulation. Here, we present a simple approach to quantify the empirical usage of individual splice sites reflecting their strength, which determines their selection in a splicing reaction. Splice-site strength/usage, as a quantitative phenotype, allows us to directly link genetic variation with usage of individual splice-sites. We demonstrate the power of this approach in defining the genomic determinants of splice-site choice through GWAS. Our pilot analysis with more than a thousand splice sites hints that sequence divergence in cis rather than trans is associated with variations in splicing among accessions of Arabidopsis thaliana. This approach allows deciphering principles of splicing and has broad implications from agriculture to medicine.

Breeding for Climate Change Resilience: A Case Study of Loblolly Pine (Pinus taeda L.) in North America

Earth's atmosphere is warming and the effects of climate change are becoming evident. A key observation is that both the average levels and the variability of temperature and precipitation are changing. Information and data from new technologies are developing in parallel to provide multidisciplinary opportunities to address and overcome the consequences of these changes in forest ecosystems. Changes in temperature and water availability impose multidimensional environmental constraints that trigger changes from the molecular to the forest stand level. These can represent a threat for the normal development of the tree from early seedling recruitment to adulthood both through direct mortality, and by increasing susceptibility to pathogens, insect attack, and fire damage. This review summarizes the strengths and shortcomings of previous work in the areas of genetic variation related to cold and drought stress in forest species with particular emphasis on loblolly pine (Pinus taeda L.), the most-planted tree species in North America. We describe and discuss the implementation of management and breeding strategies to increase resilience and adaptation, and discuss how new technologies in the areas of engineering and genomics are shaping the future of phenotype-genotype studies. Lessons learned from the study of species important in intensively-managed forest ecosystems may also prove to be of value in helping less-intensively managed forest ecosystems adapt to climate change, thereby increasing the sustainability and resilience of forestlands for the future.

Wheat (Triticum aestivum L.) TaHMW1D Transcript Variants Are Highly Expressed in Response to Heat Stress and in Grains Located in Distal Part of the Spike

High-temperature stress during the grain filling stage has a deleterious effect on grain yield and end-use quality. Plants undergo various transcriptional events of protein complexity as defensive responses to various stressors. The "Keumgang" wheat cultivar was subjected to high-temperature stress for 6 and 10 days beginning 9 days after anthesis, then two-dimensional gel electrophoresis (2DE) and peptide analyses were performed. Spots showing decreased contents in stressed plants were shown to have strong similarities with a high-molecular glutenin gene, TraesCS1D02G317301 (TaHMW1D). QRT-PCR results confirmed that TaHMW1D was expressed in its full form and in the form of four different transcript variants. These events always occurred between repetitive regions at specific deletion sites (5'-CAA (Glutamine) GG/TG (Glycine) or (Valine)-3', 5'-GGG (Glycine) CAA (Glutamine) -3') in an exonic region. Heat stress led to a significant increase in the expression of the transcript variants. This was most evident in the distal parts of the spike. Considering the importance of high-molecular weight glutenin subunits of seed storage proteins, stressed plants might choose shorter polypeptides while retaining glutenin function, thus maintaining the expression of glutenin motifs and conserved sites.

Genome-wide association study identified candidate genes for seed size and seed composition improvement in M. truncatula

Grain legumes are highly valuable plant species, as they produce seeds with high protein content. Increasing seed protein production and improving seed nutritional quality represent an agronomical challenge in order to promote plant protein consumption of a growing population. In this study, we used the genetic diversity, naturally present in Medicago truncatula, a model plant for legumes, to identify genes/loci regulating seed traits. Indeed, using sequencing data of 162 accessions from the Medicago HAPMAP collection, we performed genome-wide association study for 32 seed traits related to seed size and seed composition such as seed protein content/concentration, sulfur content/concentration. Using different GWAS and postGWAS methods, we identified 79 quantitative trait nucleotides (QTNs) as regulating seed size, 41 QTNs for seed composition related to nitrogen (i.e. storage protein) and sulfur (i.e. sulfur-containing amino acid) concentrations/contents. Furthermore, a strong positive correlation between seed size and protein content was revealed within the selected Medicago HAPMAP collection. In addition, several QTNs showed highly significant associations in different seed phenotypes for further functional validation studies, including one near an RNA-Binding Domain protein, which represents a valuable candidate as central regulator determining both seed size and composition. Finally, our findings in M. truncatula represent valuable resources to be exploitable in many legume crop species such as pea, common bean, and soybean due to its high synteny, which enable rapid transfer of these results into breeding programs and eventually help the improvement of legume grain production.

CRISPR-Mediated Engineering across the Central Dogma in Plant Biology for Basic Research and Crop Improvement

The central dogma (CD) of molecular biology is the transfer of genetic information from DNA to RNA to protein. Major CD processes governing genetic flow include the cell cycle, DNA replication, chromosome packaging, epigenetic changes, transcription, posttranscriptional alterations, translation, and posttranslational modifications. The CD processes are tightly regulated in plants to maintain genetic integrity throughout the life cycle and to pass genetic materials to next generation. Engineering of various CD processes involved in gene regulation will accelerate crop improvement to feed the growing world population. CRISPR technology enables programmable editing of CD processes to alter DNA, RNA, or protein, which would have been impossible in the past. Here, an overview of recent advancements in CRISPR tool development and CRISPR-based CD modulations that expedite basic and applied plant research is provided. Furthermore, CRISPR applications in major thriving areas of research, such as gene discovery (allele mining and cryptic gene activation), introgression (de novo domestication and haploid induction), and application of desired traits beneficial to farmers or consumers (biotic/abiotic stress-resilient crops, plant cell factories, and delayed senescence), are described. Finally, the global regulatory policies, challenges, and prospects for CRISPR-mediated crop improvement are discussed.

Identification and Regulation of Tomato Serine/Arginine-Rich Proteins Under High Temperatures

Alternative splicing is an important mechanism for the regulation of gene expression in eukaryotes during development, cell differentiation or stress response. Alterations in the splicing profiles of genes under high temperatures that cause heat stress (HS) can impact the maintenance of cellular homeostasis and thermotolerance. Consequently, information on factors involved in HS-sensitive alternative splicing is required to formulate the principles of HS response. Serine/arginine-rich (SR) proteins have a central role in alternative splicing. We aimed for the identification and characterization of SR-coding genes in tomato (Solanum lycopersicum), a plant extensively used in HS studies. We identified 17 canonical SR and two SR-like genes. Several SR-coding genes show differential expression and altered splicing profiles in different organs as well as in response to HS. The transcriptional induction of five SR and one SR-like genes is partially dependent on the master regulator of HS response, HS transcription factor HsfA1a. Cis-elements in the promoters of these SR genes were predicted, which can be putatively recognized by HS-induced transcription factors. Further, transiently expressed SRs show reduced or steady-state protein levels in response to HS. Thus, the levels of SRs under HS are regulated by changes in transcription, alternative splicing and protein stability. We propose that the accumulation or reduction of SRs under HS can impact temperature-sensitive alternative splicing.

Processing of coding and non-coding RNAs in plant development and environmental responses

Precursor RNAs undergo extensive processing to become mature RNAs. RNA transcripts are subjected to 5′ capping, 3′-end processing, splicing, and modification; they also form dynamic secondary structures during co-transcriptional and post-transcriptional processing. Like coding RNAs, non-coding RNAs (ncRNAs) undergo extensive processing. For example, secondary small interfering RNA (siRNA) transcripts undergo RNA processing, followed by further cleavage to become mature siRNAs. Transcriptome studies have revealed roles for co-transcriptional and post-transcriptional RNA processing in the regulation of gene expression and the coordination of plant development and plant–environment interactions. In this review, we present the latest progress on RNA processing in gene expression and discuss phased siRNAs (phasiRNAs), a kind of germ cell-specific secondary small RNA (sRNA), focusing on their functions in plant development and environmental responses.

Label-Free Quantitative Phosphoproteomics Reveals Signaling Dynamics Involved in Embryogenic Competence Acquisition in Sugarcane

In this study, a label-free quantitative phosphoproteomic analysis was performed to identify and quantify signaling events related to the acquisition of embryogenic competence in sugarcane. Embryogenic and nonembryogenic calli were compared at the multiplication phase, resulting in the identification of 163 phosphoproteins unique to embryogenic calli, 9 unique to nonembryogenic calli, and 51 upregulated and 40 downregulated in embryogenic calli compared to nonembryogenic calli. Data are available via ProteomeXchange with identifier PXD018054. Motif-x analysis revealed the enrichment of [xxxpSPxxx], [RxxpSxxx], and [xxxpSDxxx] motifs, which are predicted phosphorylation sites for several kinases related to stress responses. The embryogenic-related phosphoproteins (those unique and upregulated in embryogenic calli) identified in the present study are related to abscisic acid-induced signaling and abiotic stress response; they include OSK3, ABF1, LEAs, and RD29Bs. On the other hand, the nonembryogenic-related phosphoproteins EDR1 and PP2Ac-2 are negative regulators of abscisic acid signaling, suggesting a relationship between phosphoproteins involved in the abscisic acid and stress responses in the acquisition of embryogenic competence. Moreover, embryogenic-related phosphoproteins associated with epigenetic modifications, such as HDA6, HDA19, and TOPLESS, and with RNA metabolism, including AGO1, DEAH5, SCL30, UB2C, and SR45, were identified to play potential roles in embryogenic competence. These results reveal novel phosphorylation sites for several proteins and identify potential candidate biomarkers for the acquisition of embryogenic competence in sugarcane.

F-Box Gene D5RF Is Regulated by Agrobacterium Virulence Protein VirD5 and Essential for Agrobacterium-Mediated Plant Transformation

We previously reported that the Agrobacterium virulence protein VirD5 possesses transcriptional activation activity, binds to a specific DNA element D5RE, and is required for Agrobacterium-mediated stable transformation, but not for transient transformation. However, direct evidence for a role of VirD5 in plant transcriptional regulation has been lacking. In this study, we found that the Arabidopsis gene D5RF (coding for VirD5 response F-box protein, At3G49480) is regulated by VirD5. D5RF has two alternative transcripts of 930 bp and 1594 bp that encode F-box proteins of 309 and 449 amino acids, designated as D5RF.1 and D5RF.2, respectively. D5RF.2 has a N-terminal extension of 140 amino acids compared to D5RF.1, and both of them are located in the plant cell nucleus. The promoter of the D5RF.1 contains two D5RE elements and can be activated by VirD5. The expression of D5RF is downregulated when the host plant is infected with virD5 deleted Agrobacterium. Similar to VirD5, D5RF also affects the stable but not transient transformation efficiency of Agrobacterium. Some pathogen-responsive genes are downregulated in the d5rf mutant. In conclusion, this study further confirmed Agrobacterium VirD5 as the plant transcription activator and identified Arabidopsis thalianaD5RF.1 as the first target gene of VirD5 in regulation.

Dynamic mRNP Remodeling in Response to Internal and External Stimuli

Signal transduction and the regulation of gene expression are fundamental processes in every cell. RNA-binding proteins (RBPs) play a key role in the post-transcriptional modulation of gene expression in response to both internal and external stimuli. However, how signaling pathways regulate the assembly of RBPs with mRNAs remains largely unknown. Here, we summarize observations showing that the formation and composition of messenger ribonucleoprotein particles (mRNPs) is dynamically remodeled in space and time by specific signaling cascades and the resulting post-translational modifications. The integration of signaling events with gene expression is key to the rapid adaptation of cells to environmental changes and stress. Only a combined approach analyzing the signal transduction pathways and the changes in post-transcriptional gene expression they cause will unravel the mechanisms coordinating these important cellular processes.