Global Survey of Alternative Splicing in Rice by Direct RNA Sequencing During Reproductive Development: Landscape and Genetic Regulation

A MYB transcription factor underlying plant height in sorghum qHT7.1 and maize Brachytic 1 loci

Manipulating plant height is an essential component of crop improvement. Plant height was generally reduced through breeding in wheat, rice, and sorghum to resist lodging and increase grain yield but kept high for bioenergy crops. Here, we positionally cloned a plant height quantitative trait locus (QTL) qHT7.1 as a MYB transcription factor controlling internode elongation, cell proliferation, and cell morphology in sorghum. A 740 bp transposable element insertion in the intronic region caused a partial mis-splicing event, generating a novel transcript that included an additional exon and a premature stop codon, leading to short plant height. The dominant allele had an overall higher expression than the recessive allele across development and internode position, while both alleles' expressions peaked at 46 days after planting and progressively decreased from the top to lower internodes. The orthologue of qHT7.1 was identified to underlie the brachytic1 (br1) locus in maize. A large insertion in exon 3 and a 160 bp insertion at the promoter region were identified in the br1 mutant, while an 18 bp promoter insertion was found to be associated with reduced plant height in a natural recessive allele. CRISPR/Cas9-induced gene knockout of br1 in two maize inbred lines showed significant plant height reduction. These findings revealed functional connections across natural, mutant, and edited alleles of this MYB transcription factor in sorghum and maize. This enriched our understanding of plant height regulation and enhanced our toolbox for fine-tuning plant height for crop improvement.

Direct RNA sequencing in plants: Practical applications and future perspectives

The transcriptome serves as a bridge that links genomic variation to phenotypic diversity. A vast number of studies using next-generation RNA sequencing (RNA-seq) over the last 2 decades have emphasized the essential roles of the plant transcriptome in response to developmental and environmental conditions, providing numerous insights into the dynamic changes, evolutionary traces, and elaborate regulation of the plant transcriptome. With substantial improvement in accuracy and throughput, direct RNA sequencing (DRS) has emerged as a new and powerful sequencing platform for precise detection of native and full-length transcripts, overcoming many limitations such as read length and PCR bias that are inherent to short-read RNA-seq. Here, we review recent advances in dissecting the complexity and diversity of plant transcriptomes using DRS as the main technological approach, covering many aspects of RNA metabolism, including novel isoforms, poly(A) tails, and RNA modification, and we propose a comprehensive workflow for processing of plant DRS data. Many challenges to the application of DRS in plants, such as the need for machine learning tools tailored to plant transcriptomes, remain to be overcome, and together we outline future biological questions that can be addressed by DRS, such as allele-specific RNA modification. This technology provides convenient support on which the connection of distinct RNA features is tightly built, sustainably refining our understanding of the biological functions of the plant transcriptome.

Differences in alternative splicing and their potential underlying factors between animals and plants

BACKGROUND

Alternative splicing (AS), a posttranscriptional process, contributes to the complexity of transcripts from a limited number of genes in a genome, and AS is considered a great source of genetic and phenotypic diversity in eukaryotes. In animals, AS is tightly regulated during the processes of cell growth and differentiation, and its dysregulation is involved in many diseases, including cancers. Likewise, in plants, AS occurs in all stages of plant growth and development, and it seems to play important roles in the rapid reprogramming of genes in response to environmental stressors. To date, the prevalence and functional roles of AS have been extensively reviewed in animals and plants. However, AS differences between animals and plants, especially their underlying molecular mechanisms and impact factors, are anecdotal and rarely reviewed.

AIM OF REVIEW

This review aims to broaden our understanding of AS roles in a variety of biological processes and provide insights into the underlying mechanisms and impact factors likely leading to AS differences between animals and plants.

KEY SCIENTIFIC CONCEPTS OF REVIEW

We briefly summarize the roles of AS regulation in physiological and biochemical activities in animals and plants. Then, we underline the differences in the process of AS between plants and animals and especially analyze the potential impact factors, such as gene exon/intron architecture, 5'/3' untranslated regions (UTRs), spliceosome components, chromatin dynamics and transcription speeds, splicing factors [serine/arginine-rich (SR) proteins and heterogeneous nuclear ribonucleoproteins (hnRNPs)], noncoding RNAs, and environmental stimuli, which might lead to the differences. Moreover, we compare the nonsense-mediated mRNA decay (NMD)-mediated turnover of the transcripts with a premature termination codon (PTC) in animals and plants. Finally, we summarize the current AS knowledge published in animals versus plants and discuss the potential development of disease therapies and superior crops in the future.

Concordant Gene Expression and Alternative Splicing Regulation under Abiotic Stresses in Arabidopsis

The current investigation endeavors to identify differentially expressed alternatively spliced (DAS) genes that exhibit concordant expression with splicing factors (SFs) under diverse multifactorial abiotic stress combinations in Arabidopsis seedlings. SFs serve as the post-transcriptional mechanism governing the spatiotemporal dynamics of gene expression. The different stresses encompass variations in salt concentration, heat, intensive light, and their combinations. Clusters demonstrating consistent expression profiles were surveyed to pinpoint DAS/SF gene pairs exhibiting concordant expression. Through rigorous selection criteria, which incorporate alignment with documented gene functionalities and expression patterns observed in this study, four members of the serine/arginine-rich (SR) gene family were delineated as SFs concordantly expressed with six DAS genes. These regulated SF genes encompass cactin, SR1-like, SR30, and SC35-like. The identified concordantly expressed DAS genes encode diverse proteins such as the 26.5 kDa heat shock protein, chaperone protein DnaJ, potassium channel GORK, calcium-binding EF hand family protein, DEAD-box RNA helicase, and 1-aminocyclopropane-1-carboxylate synthase 6. Among the concordantly expressed DAS/SF gene pairs, SR30/DEAD-box RNA helicase, and SC35-like/1-aminocyclopropane-1-carboxylate synthase 6 emerge as promising candidates, necessitating further examinations to ascertain whether these SFs orchestrate splicing of the respective DAS genes. This study contributes to a deeper comprehension of the varied responses of the splicing machinery to abiotic stresses. Leveraging these DAS/SF associations shows promise for elucidating avenues for augmenting breeding programs aimed at fortifying cultivated plants against heat and intensive light stresses.

Strong culm: a crucial trait for developing next-generation climate-resilient rice lines

Lodging, a phenomenon characterized by the bending or breaking of rice plants, poses substantial constraints on productivity, particularly during the harvesting phase in regions susceptible to strong winds. The rice strong culm trait is influenced by the intricate interplay of genetic, physiological, epigenetic, and environmental factors. Stem architecture, encompassing morphological and anatomical attributes, alongside the composition of both structural and non-structural carbohydrates, emerges as a critical determinant of lodging resistance. The adaptive response of the rice culm to various biotic and abiotic environmental factors further modulates the propensity for lodging. Advancements in next-generation sequencing technologies have expedited the genetic dissection of lodging resistance, enabling the identification of pertinent genes, quantitative trait loci, and novel alleles. Concurrently, contemporary breeding strategies, ranging from biparental approaches to more sophisticated methods such as multi-parent-based breeding, gene pyramiding, genomic selection, genome-wide association studies, and haplotype-based breeding, offer perspectives on the genetic underpinnings of culm strength. This review comprehensively delves into physiological attributes, culm histology, epigenetic determinants, and gene expression profiles associated with lodging resistance, with a specialized focus on leveraging next-generation sequencing for candidate gene discovery.

Alternative Splicing Variation: Accessing and Exploiting in Crop Improvement Programs

Alternative splicing (AS) is a gene regulatory mechanism modulating gene expression in multiple ways. AS is prevalent in all eukaryotes including plants. AS generates two or more mRNAs from the precursor mRNA (pre-mRNA) to regulate transcriptome complexity and proteome diversity. Advances in next-generation sequencing, omics technology, bioinformatics tools, and computational methods provide new opportunities to quantify and visualize AS-based quantitative trait variation associated with plant growth, development, reproduction, and stress tolerance. Domestication, polyploidization, and environmental perturbation may evolve novel splicing variants associated with agronomically beneficial traits. To date, pre-mRNAs from many genes are spliced into multiple transcripts that cause phenotypic variation for complex traits, both in model plant Arabidopsis and field crops. Cataloguing and exploiting such variation may provide new paths to enhance climate resilience, resource-use efficiency, productivity, and nutritional quality of staple food crops. This review provides insights into AS variation alongside a gene expression analysis to select for novel phenotypic diversity for use in breeding programs. AS contributes to heterosis, enhances plant symbiosis (mycorrhiza and rhizobium), and provides a mechanistic link between the core clock genes and diverse environmental clues.

TRIMETHYLGUANOSINE SYNTHASE1 mutations decanalize female germline development in Arabidopsis

Here, we report the characterization of a plant RNA methyltransferase, orthologous to yeast trimethylguanosine synthase1 (Tgs1p) and whose downregulation was associated with apomixis in Paspalum grasses. Using phylogenetic analyses and yeast complementation, we determined that land plant genomes all encode a conserved, specific TGS1 protein. Next, we studied the role of TGS1 in female reproduction using reporter lines and loss-of-function mutants in Arabidopsis thaliana. pAtTGS1:AtTGS1 reporters showed a dynamic expression pattern. They were highly active in the placenta and ovule primordia at emergence but, subsequently, showed weak signals in the nucellus. Although expressed throughout gametophyte development, activity became restricted to the female gamete and was also detected after fertilization during embryogenesis. TGS1 depletion altered the specification of the precursor cells that give rise to the female gametophytic generation and to the sporophyte, resulting in the formation of a functional aposporous-like lineage. Our results indicate that TGS1 participates in the mechanisms restricting cell fate acquisition to a single cell at critical transitions throughout the female reproductive lineage and, thus, expand our current knowledge of the mechanisms governing female reproductive fate in plants.

Nuclear OsFKBP20-1b maintains SR34 stability and promotes the splicing of retained introns upon ABA exposure in rice

Alternative splicing (AS) is a critical means by which plants respond to changes in the environment, but few splicing factors contributing to AS have been reported and functionally characterized in rice (Oryza sativa L.). Here, we explored the function and molecular mechanism of the spliceosome-associated protein OsFKBP20-1b during AS. We determined the AS landscape of wild-type and osfkbp20-1b knockout plants upon abscisic acid (ABA) treatment by transcriptome deep sequencing. To capture the dynamics of translating intron-containing mRNAs, we blocked transcription with cordycepin and performed polysome profiling. We also analyzed whether OsFKBP20-1b and the splicing factors OsSR34 and OsSR45 function together in AS using protoplast transfection assays. We show that OsFKBP20-1b interacts with OsSR34 and regulates its stability, suggesting a role as a chaperone-like protein in the spliceosome. OsFKBP20-1b facilitates the splicing of mRNAs with retained introns after ABA treatment; some of these mRNAs are translatable and encode functional transcriptional regulators of stress-responsive genes. In addition, interacting proteins, OsSR34 and OsSR45, regulate the splicing of the same retained introns as OsFKBP20-1b after ABA treatment. Our findings reveal that spliceosome-associated immunophilin functions in alternative RNA splicing in rice by positively regulating the splicing of retained introns to limit ABA response.

Do not panic: An intron-centric guide to alternative splicing

This review is an attempt to establish concepts of splicing and alternative splicing giving proper relevance to introns, the key actors in this mechanism. It might also work as a guide for those who found their favorite gene undergoes alternative splicing and could benefit from gaining a theoretical framework to understand the possible impacts of this process. This is not a thorough review of all the work in the field, but rather a critical review of some of the most relevant work done to understand the underlying mechanisms of splicing and the key questions that remain unanswered such as: What is the physiological relevance of alternative splicing? What are the functions of the different outcomes? To what extent do different alternative splicing types contribute to the proteome? Intron retention is the most frequent alternative splicing event in plants and, although scientifically neglected, it is also common in animals. This is a heterogeneous type of alternative splicing that includes different sub-types with features that have distinctive consequences in the resulting transcripts. Remarkably, intron retention can be a dead end for a transcript, but it could also be a stable intermediate whose processing is resumed upon a particular signal or change in the cell status. New sequencing technologies combined with the study of intron lariats in different conditions might help to answer key questions and could help us to understand the actual relevance of introns in gene expression regulation.

OsGRF4AA compromises heat tolerance of developing pollen grains in rice

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

MAEL gene contributes to bovine testicular development through the m5C-mediated splicing

Knowledge of RNA molecules regulating testicular development and spermatogenesis in bulls is essential for elite bull selection and an ideal breeding program. Herein, we performed direct RNA sequencing (DRS) to explore the functional characterization of RNA molecules produced in the testicles of 9 healthy Simmental bulls at three testicular development stages (prepuberty, puberty, and postpuberty). We identified 5,043 differentially expressed genes associated with testicular weight. These genes exhibited more alternative splicing at sexual maturity, particularly alternative 3' (A3) and 5' (A5) splice sites usage and exon skipping (SE). The expression of hub genes in testicular developmental stages was also affected by both m6A and m5C RNA modifications. We found m5C-mediated splicing events significantly (p < 0.05) increased MAEL gene expression at the isoform level, likely promoting spermatogenesis. Our findings highlight the complexity of RNA processing and expression as well as the regulation of transcripts involved in testicular development and spermatogenesis.

Integrated omics analysis identified genes and their splice variants involved in fruit development and metabolites production in Capsicum species

To date, several transcriptomic studies during fruit development have been reported; however, no comprehensive integrated study on expression diversity, alternative splicing, and metabolomic profiling was reported in Capsicum. This study analyzed RNA-seq data and untargeted metabolomic profiling from early green (EG), mature green (MG), and breaker (Br) fruit stages from two Capsicum species, i.e., C. annuum (Cann) and C. frutescens (Cfrut) from Northeast India. A total of 117,416 and 96,802 alternatively spliced events (AltSpli-events) were identified from Cann and Cfrut, respectively. Among AltSpli-events, intron retention (IR; 32.2% Cann and 25.75% Cfrut) followed by alternative acceptor (AA; 15.4% Cann and 18.9% Cfrut) were the most abundant in Capsicum. Around 7600 genes expressed in at least one fruit stage of Cann and Cfrut were AltSpli. The study identified spliced variants of genes including transcription factors (TFs) potentially involved in fruit development/ripening (Aux/IAA 16-like, ETR, SGR1, ARF, CaGLK2, ETR, CaAGL1, MADS-RIN, FUL1, SEPALLATA1), carotenoid (PDS, CA1, CCD4, NCED3, xanthoxin dehydrogenase, CaERF82, CabHLH100, CaMYB3R-1, SGR1, CaWRKY28, CaWRKY48, CaWRKY54), and capsaicinoids or flavonoid biosynthesis (CaMYB48, CaWRKY51), which were significantly differentially spliced (DS) between consecutive Capsicum fruit stages. Also, this study observed that differentially expressed isoforms (DEiso) from 38 genes with differentially spliced events (DSE) were significantly enriched in various metabolic pathways such as starch and sucrose metabolism, amino acid metabolism, cysteine cutin suberin and wax biosynthesis, and carotenoid biosynthesis. Furthermore, the metabolomic profiling revealed that metabolites from aforementioned pathways such as carbohydrates (mainly sugars such as D-fructose, D-galactose, maltose, and sucrose), organic acids (carboxylic acids), and peptide groups significantly altered during fruit development. Taken together, our findings could help in alternative splicing-based targeted studies of candidate genes involved in fruit development and ripening in Capsicum crop.

Comparative transcriptomics reveals new insights into melatonin-enhanced drought tolerance in naked oat seedlings

The growth and development of naked oat (Avena nuda L.) seedlings, a grain recognized as nutritious and healthy, is limited by drought. Melatonin plays a positive role in plants under drought stress. However, its function is unclear in naked oats. This study demonstrated that melatonin enhances drought stress tolerance in oat seedlings. Melatonin application alleviated the declining growth parameters of two naked oat varieties, Huazao No.2 (H2) and Jizhangyou No.15 (J15), under drought stress by increasing the chlorophyll content and photosynthetic rate of leaves. Melatonin pretreatment induced differential gene expression in H2 and J15 under drought stress. Subsequently, the differential gene expression responses to melatonin in the two varieties were further analyzed. The key drought response transcription factors and the regulatory effect of melatonin on drought-related transcription factors were assessed, focusing on genes encoding proteins in the ABA signal transduction pathway, including PYL, PP2C, ABF, SNRK2, and IAA. Taken together, this study provides new insights into the effect and underlying mechanism of melatonin in alleviating drought stress in naked oat seedlings.