Differential nucleosome occupancy modulates alternative splicing in Arabidopsis thaliana

At-RS31 orchestrates hierarchical cross-regulation of splicing factors and integrates alternative splicing with TOR-ABA pathways

Alternative splicing is essential for plants, enabling a single gene to produce multiple transcript variants to boost functional diversity and fine-tune responses to environmental and developmental cues. At-RS31, a plant-specific splicing factor in the Serine/Arginine (SR)-rich protein family, responds to light and the Target of Rapamycin (TOR) signaling pathway, yet its downstream targets and regulatory impact remain unknown.To identify At-RS31 targets, we applied individual-nucleotide resolution crosslinking and immunoprecipitation (iCLIP) and RNAcompete assays. Transcriptomic analyses of At-RS31 mutant and overexpressing plants further revealed its effects on alternative splicing.iCLIP identified 4,034 At-RS31 binding sites across 1,421 genes, enriched in CU-rich and CAGA RNA motifs. Comparative iCLIP and RNAcompete data indicate that the RS domain of At-RS31 may influence its binding specificity in planta, underscoring the value of combining in vivo and in vitro approaches. Transcriptomic analysis showed that At-RS31 modulates diverse splicing events, particularly intron retention and exitron splicing, and influences other splicing modulators, acting as a hierarchical regulator.By regulating stress-response genes and genes in both TOR and abscisic acid (ABA) signaling pathways, At-RS31 may help integrate these signals, balancing plant growth with environmental adaptability through alternative splicing.

Rare occurrence of cryptic 5' splice sites by downstream 3' splice site/exon boundary mutations in a heavy-ion-induced egy1-4 allele of Arabidopsis thaliana

Pre-mRNA splicing is a fundamental process in eukaryotic gene expression, and the mechanism of intron definition, involving the recognition of the canonical GU (5'-splice site) and AG (3'-splice site) dinucleotides by splicing factors, has been postulated for most cases of splicing initiation in plants. Splice site mutations have played crucial roles in unraveling the mechanism of pre-mRNA splicing in planta. Typically, splice site mutations abolish splicing events or activate one or more cryptic splice sites surrounding the mutated region. In this report, we investigated the splicing pattern of the EGY1 gene in an Ar-ion-induced egy1-4 allele of Arabidopsis thaliana. egy1-4 has an AG-to-AC mutation in the 3'-end of intron 3, along with 4-bp substitutions and a 5-bp deletion in adjacent exon 4. RT-PCR, cDNA cloning, and amplicon sequencing analyses of EGY1 revealed that while most wild-type EGY1 mRNAs had a single splicing pattern, egy1-4 mRNAs had multiple splicing defects. Almost half of EGY1 transcripts showed 'intron retention' at intron 3, while the other half exhibited activation of 3' cryptic splice sites either upstream or downstream of the original 3'-splice site. Unexpectedly, around 8% of EGY1 transcripts in egy1-4 exhibited activation of cryptic 5'-splice sites positioned upstream of the authentic 5'-splice site of intron 3. Whole genome resequencing of egy1-4 indicated that it has no other known impactful mutations. These results may provide a rare, but real case of activation of cryptic 5'-splice sites by downstream 3'-splice site/exon mutations in planta.

NucMap 2.0: An Updated Database of Genome-wide Nucleosome Positioning Maps Across Species

Nucleosome dynamics plays important roles in many biological processes, such as DNA replication and gene expression. NucMap (https://ngdc.cncb.ac.cn/nucmap) is the first database of genome-wide nucleosome positioning maps across species. Here, we present an updated version, NucMap 2.0, by incorporating more species and MNase-seq samples. In addition, we integrate other related omics data for each MNase-seq sample to provide a comprehensive view of nucleosome positioning, such as gene expression, transcription factor binding sites, histone modifications and DNA methylation. In particular, NucMap 2.0 integrates and pre-analyzes RNA-seq data and ChIP-seq data of human-related samples, which facilitates the interpretation of nucleosome positioning in humans. All processed data are integrated into an in-built genome browser, and users can make comprehensive side-by-side analyses. In addition, more online analytical functions are developed, which allows researchers to identify differential nucleosome regions and explore potential gene regulatory regions. All resources are open access with a user-friendly web interface.

Phytopathogens Reprogram Host Alternative mRNA Splicing

Alternative splicing (AS) is an evolutionarily conserved cellular process in eukaryotes in which multiple messenger RNA (mRNA) transcripts are produced from a single gene. The concept that AS adds to transcriptome complexity and proteome diversity introduces a new perspective for understanding how phytopathogen-induced alterations in host AS cause diseases. Recently, it has been recognized that AS represents an integral component of the plant immune system during parasitic, commensalistic, and symbiotic interactions. Here, I provide an overview of recent progress detailing the reprogramming of plant AS by phytopathogens and the functional implications on disease phenotypes. Additionally, I discuss the vital function of AS of immune receptors in regulating plant immunity and how phytopathogens use effector proteins to target key components of the splicing machinery and exploit alternatively spliced variants of immune regulators to negate defense responses. Finally, the functional association between AS and nonsense-mediated mRNA decay in the context of plant-pathogen interface is recapitulated.

Deep learning with a small dataset predicts chromatin remodelling contribution to winter dormancy of apple axillary buds

Epigenetic changes serve as a cellular memory for cumulative cold recognition in both herbaceous and tree species, including bud dormancy. However, most studies have discussed predicted chromatin structure with respect to histone marks. In the present study, we investigated the structural dynamics of bona fide chromatin to determine how plants recognize prolonged chilling during the initial stage of bud dormancy. The vegetative axillary buds of the 'Fuji' apple, which shows typical low temperature-dependent, but not photoperiod, dormancy induction, were used for the chromatin structure and transcriptional change analyses. The results were integrated using a deep-learning model and interpreted using statistical models, including Bayesian estimation. Although our model was constructed using a small dataset of two time points, chromatin remodelling due to random changes was excluded. The involvement of most nucleosome structural changes in transcriptional changes and the pivotal contribution of cold-driven circadian rhythm-dependent pathways regulated by the mobility of cis-regulatory elements were predicted. These findings may help to develop potential genetic targets for breeding species with less bud dormancy to overcome the effects of short winters during global warming. Our artificial intelligence concept can improve epigenetic analysis using a small dataset, especially in non-model plants with immature genome databases.

Cold-induced deposition of bivalent H3K4me3-H3K27me3 modification and nucleosome depletion in Arabidopsis

Epigenetic regulation of gene expression plays a crucial role in plant development and environmental adaptation. The H3K4me3 and H3K27me3 have not only been discovered in the regulation of gene expression in multiple biological processes but also in responses to abiotic stresses in plants. However, evidence for the presence of both H3K4me3 and H3K27me3 on the same nucleosome is sporadic. Cold-induced deposition of bivalent H3K4me3-H3K27me3 modifications and nucleosome depletion over a considerable number of active genes is documented in potato tubers and provides clues on an additional role of the bivalent modifications. Limited by the available information of genes encoding PcG/TrxG proteins as well as their corresponding mutants in potatoes, the molecular mechanism underlying the cold-induced deposition of the bivalent mark remains elusive. In this study, we found a similar deposition of the bivalent H3K4me3-H3K27me3 mark over 2129 active genes in cold-treated Arabidopsis Col-0 seedlings. The expression levels of the bivalent mark-associated genes tend to be independent of bivalent modification levels. However, these genes were associated with greater chromatin accessibility, presumably to provide a distinct chromatin environment for gene expression. In mutants clf28 and lhp1, failure to deposit H3K27me3 in active genes upon cold treatment implies that the CLF is potentially involved in cold-induced deposition of H3K27me3, with assistance from LHP1. Failure to deposit H3K4me3 during cold treatment in atx1-2 suggests a regulatory role of ATX1 in the deposition of H3K4me3. In addition, we observed a cold-induced global reduction in nucleosome occupancy, which is potentially mediated by LHP1 in an H3K27me3-dependent manner.

Full-length transcriptome sequencing provides insights into alternative splicing under cold stress in peanut

Introduction

Peanut (Arachis hypogaea L.), also called groundnut is an important oil and cash crop grown widely in the world. The annual global production of groundnuts has increased to approximately 50 million tons, which provides a rich source of vegetable oils and proteins for humans. Low temperature (non-freezing) is one of the major factors restricting peanut growth, yield, and geographic distribution. Since the complexity of cold-resistance trait, the molecular mechanism of cold tolerance and related gene networks were largely unknown in peanut.

Methods

In this study, comparative transcriptomic analysis of two peanut cultivars (SLH vs. ZH12) with differential cold tolerance under low temperature (10°C) was performed using Oxford Nanopore Technology (ONT) platform.

Results and discussion

As a result, we identified 8,949 novel gene loci and 95,291 new/novel isoforms compared with the reference database. More differentially expressed genes (DEGs) were discovered in cold-sensitive cultivar (ZH12) than cold-tolerant cultivar (SLH), while more alternative splicing events were found in SLH compared to ZH12. Gene Ontology (GO) analyses of the common DEGs showed that the "response to stress", "chloroplast part", and "transcription factor activity" were the most enriched GO terms, indicating that photosynthesis process and transcription factors play crucial roles in cold stress response in peanut. We also detected a total of 708 differential alternative splicing genes (DASGs) under cold stress compared to normal condition. Intron retention (IR) and exon skipping (ES) were the most prevalent alternative splicing (AS) events. In total, 4,993 transcription factors and 292 splicing factors were detected, many of them had differential expression levels and/or underwent AS events in response to cold stress. Overexpression of two candidate genes (encoding trehalose-6-phosphatephosphatases, AhTPPs) in yeast improves cold tolerance. This study not only provides valuable resources for the study of cold resistance in peanut but also lay a foundation for genetic modification of cold regulators to enhance stress tolerance in crops.

Unveiling unique alternative splicing responses to low temperature in Zoysia japonica through ZjRTD1.0, a high-quality reference transcript dataset

Inadequate reference databases in RNA-seq analysis can hinder data utilization and interpretation. In this study, we have successfully constructed a high-quality reference transcript dataset, ZjRTD1.0, for Zoysia japonica, a widely-used turfgrass with exceptional tolerance to various abiotic stress, including low temperatures and salinity. This dataset comprises 113,089 transcripts from 57,143 genes. BUSCO analysis demonstrates exceptional completeness (92.4%) in ZjRTD1.0, with reduced proportions of fragmented (3.3%) and missing (4.3%) orthologs compared to prior datasets. ZjRTD1.0 enables more precise analyses, including transcript quantification and alternative splicing assessments using public datasets, which identified a substantial number of differentially expressed transcripts (DETs) and differential alternative splicing (DAS) events, leading to several novel findings on Z. japonica's responses to abiotic stresses. First, spliceosome gene expression influenced alternative splicing significantly under abiotic stress, with a greater impact observed during low-temperature stress. Then, a significant positive correlation was found between the number of differentially expressed genes (DEGs) encoding protein kinases and the frequency of DAS events, suggesting the role of protein phosphorylation in regulating alternative splicing. Additionally, our results suggest possible involvement of serine/arginine-rich (SR) proteins and heterogeneous nuclear ribonucleoproteins (hnRNPs) in generating inclusion/exclusion isoforms under low-temperature stress. Furthermore, our investigation revealed a significantly enhanced overlap between DEGs and differentially alternatively spliced genes (DASGs) in response to low-temperature stress, suggesting a unique co-regulatory mechanism governing transcription and splicing in the context of low-temperature response. In conclusion, we have proven that ZjRTD1.0 will serve as a reliable and useful resource for future transcriptomic analyses in Z. japonica.

Exitrons: offering new roles to retained introns-the novel regulators of protein diversity and utility

Exitrons are exonic introns. This subclass of intron retention alternative splicing does not contain a Pre-Terminating stop Codon. Therefore, when retained, they are always a part of a protein. Intron retention is a frequent phenomenon predominantly found in plants, which results in either the degradation of the transcripts or can serve as a stable intermediate to be processed upon induction by specific signals or the cell status. Interestingly, exitrons have coding ability and may confer additional attributes to the proteins that retain them. Therefore, exitron-containing and exitron-spliced isoforms will be a driving force for creating protein diversity in the proteome of an organism. This review establishes a basic understanding of exitron, discussing its genesis, key features, identification methods and functions. We also try to depict its other potential roles. The present review also aims to provide a fundamental background to those who found such exitronic sequences in their gene(s) and to speculate the future course of studies.

R-loops act as regulatory switches modulating transcription of COLD-responsive genes in rice

COLD is a major naturally occurring stress that usually causes complex symptoms and severe yield loss in crops. R-loops function in various cellular processes, including development and stress responses, in plants. However, how R-loops function in COLD responses is largely unknown in COLD susceptible crops like rice (Oryza sativa L.). We conducted DRIP-Seq along with other omics data (RNA-Seq, DNase-Seq and ChIP-Seq) in rice with or without COLD treatment. COLD treatment caused R-loop reprogramming across the genome. COLD-biased R-loops had higher GC content and novel motifs for the binding of distinct transcription factors (TFs). Moreover, R-loops can directly/indirectly modulate the transcription of a subset of COLD-responsive genes, which can be mediated by R-loop overlapping TF-centered or cis-regulatory element-related regulatory networks and lncRNAs, accounting for c. 60% of COLD-induced expression of differential genes in rice, which is different from the findings in Arabidopsis. We validated two R-loop loci with contrasting (negative/positive) roles in the regulation of two individual COLD-responsive gene expression, as potential targets for enhanced COLD resistance. Our study provides detailed evidence showing functions of R-loop reprogramming during COLD responses and provides some potential R-loop loci for genetic and epigenetic manipulation toward breeding of rice varieties with enhanced COLD tolerance.

DNA methylome analysis provides insights into gene regulatory mechanism for better performance of rice under fluctuating environmental conditions: epigenomics of adaptive plasticity

MAIN CONCLUSION

Roots play an important role in adaptive plasticity of rice under dry/direct-sown conditions. However, hypomethylation of genes in leaves (resulting in up-regulated expression) complements the adaptive plasticity of Nagina-22 under DSR conditions. Rice is generally cultivated by transplanting which requires plenty of water for irrigation. Such a practice makes rice cultivation a challenging task under global climate change and reducing water availability. However, dry-seeded/direct-sown rice (DSR) has emerged as a resource-saving alternative to transplanted rice (TPR). Though some of the well-adapted local cultivars are used for DSR, only limited success has been achieved in developing DSR varieties mainly because of a limited knowledge of adaptability of rice under fluctuating environmental conditions. Based on better morpho-physiological and agronomic performance of Nagina-22 (N-22) under DSR conditions, N-22 and IR-64 were grown by transplanting and direct-sowing and used for whole genome methylome analysis to unravel the epigenetic basis of adaptive plasticity of rice. Comparative methylome and transcriptome analyses indicated a large number (4078) of genes regulated through DNA methylation/demethylation in N-22 under DSR conditions. Gene × environment interactions play important roles in adaptive plasticity of rice under direct-sown conditions. While genes for pectinesterase, LRK10, C2H2 zinc-finger protein, splicing factor, transposable elements, and some of the unannotated proteins were hypermethylated, the genes for regulation of transcription, protein phosphorylation, etc. were hypomethylated in CG context in the root of N-22, which played important roles in providing adaptive plasticity to N-22 under DSR conditions. Hypomethylation leading to up-regulation of gene expression in the leaf complements the adaptive plasticity of N-22 under DSR conditions. Moreover, differential post-translational modification of proteins and chromatin assembly/disassembly through DNA methylation in CHG context modulate adaptive plasticity of N-22. These findings would help developing DSR cultivars for increased water-productivity and ecological efficiency.

Is winter coming? Impact of the changing climate on plant responses to cold temperature

Climate change is causing alterations in annual temperature regimes worldwide. Important aspects of this include the reduction of winter chilling temperatures as well as the occurrence of unpredicted frosts, both significantly affecting plant growth and yields. Recent studies advanced the knowledge of the mechanisms underlying cold responses and tolerance in the model plant Arabidopsis thaliana. However, how these cold-responsive pathways will readjust to ongoing seasonal temperature variation caused by global warming remains an open question. In this review, we highlight the plant developmental programmes that depend on cold temperature. We focus on the molecular mechanisms that plants have evolved to adjust their development and stress responses upon exposure to cold. Covering both genetic and epigenetic aspects, we present the latest insights into how alternative splicing, noncoding RNAs and the formation of biomolecular condensates play key roles in the regulation of cold responses. We conclude by commenting on attractive targets to accelerate the breeding of increased cold tolerance, bringing up biotechnological tools that might assist in overcoming current limitations. Our aim is to guide the reflection on the current agricultural challenges imposed by a changing climate and to provide useful information for improving plant resilience to unpredictable cold regimes.

The Global Changes of N6-methyldeoxyadenosine in Response to Low Temperature in Arabidopsis thaliana and Rice

N6-methyldeoxyadenosine (6mA) is a recently discovered DNA modification involved in regulating plant adaptation to abiotic stresses. However, the mechanisms and changes of 6mA under cold stress in plants are not yet fully understood. Here, we conducted a genome-wide analysis of 6mA and observed that 6mA peaks were predominantly present within the gene body regions under both normal and cold conditions. In addition, the global level of 6mA increased both in Arabidopsis and rice after the cold treatment. The genes that exhibited an up-methylation showed enrichment in various biological processes, whereas there was no significant enrichment observed among the down-methylated genes. The association analysis revealed a positive correlation between the 6mA level and the gene expression level. Joint analysis of the 6mA methylome and transcriptome of Arabidopsis and rice unraveled that fluctuations in 6mA levels caused by cold exposure were not correlated to changes in transcript levels. Furthermore, we discovered that orthologous genes modified by 6mA showed high expression levels; however, only a minor amount of differentially 6mA-methylated orthologous genes were shared between Arabidopsis and rice under low-temperature conditions. In conclusion, our study provides information on the role of 6mA in response to cold stress and reveals its potential for regulating the expression of stress-related genes.

PlaASDB: a comprehensive database of plant alternative splicing events in response to stress

BACKGROUND

Alternative splicing (AS) is a co-transcriptional regulatory mechanism of plants in response to environmental stress. However, the role of AS in biotic and abiotic stress responses remains largely unknown. To speed up our understanding of plant AS patterns under different stress responses, development of informative and comprehensive plant AS databases is highly demanded.

DESCRIPTION

In this study, we first collected 3,255 RNA-seq data under biotic and abiotic stresses from two important model plants (Arabidopsis and rice). Then, we conducted AS event detection and gene expression analysis, and established a user-friendly plant AS database termed PlaASDB. By using representative samples from this highly integrated database resource, we compared AS patterns between Arabidopsis and rice under abiotic and biotic stresses, and further investigated the corresponding difference between AS and gene expression. Specifically, we found that differentially spliced genes (DSGs) and differentially expressed genes (DEG) share very limited overlapping under all kinds of stresses, suggesting that gene expression regulation and AS seemed to play independent roles in response to stresses. Compared with gene expression, Arabidopsis and rice were more inclined to have conserved AS patterns under stress conditions.

CONCLUSION

PlaASDB is a comprehensive plant-specific AS database that mainly integrates the AS and gene expression data of Arabidopsis and rice in stress response. Through large-scale comparative analyses, the global landscape of AS events in Arabidopsis and rice was observed. We believe that PlaASDB could help researchers understand the regulatory mechanisms of AS in plants under stresses more conveniently. PlaASDB is freely accessible at http://zzdlab.com/PlaASDB/ASDB/index.html .

Alternative Splicing in the Regulatory Circuit of Plant Temperature Response

As sessile organisms, plants have evolved complex mechanisms to rapidly respond to ever-changing ambient temperatures. Temperature response in plants is modulated by a multilayer regulatory network, including transcriptional and post-transcriptional regulations. Alternative splicing (AS) is an essential post-transcriptional regulatory mechanism. Extensive studies have confirmed its key role in plant temperature response, from adjustment to diurnal and seasonal temperature changes to response to extreme temperatures, which has been well documented by previous reviews. As a key node in the temperature response regulatory network, AS can be modulated by various upstream regulations, such as chromatin modification, transcription rate, RNA binding proteins, RNA structure and RNA modifications. Meanwhile, a number of downstream mechanisms are affected by AS, such as nonsense-mediated mRNA decay (NMD) pathway, translation efficiency and production of different protein variants. In this review, we focus on the links between splicing regulation and other mechanisms in plant temperature response. Recent advances regarding how AS is regulated and the following consequences in gene functional modulation in plant temperature response will be discussed. Substantial evidence suggests that a multilayer regulatory network integrating AS in plant temperature response has been unveiled.

Alternative splicing as a key player in the fine-tuning of the immunity response in Arabidopsis

Plants, like animals, are constantly exposed to abiotic and biotic stresses, which often inhibit plant growth and development, and cause tissue damage, disease, and even plant death. Efficient and timely response to stress requires appropriate co- and posttranscriptional reprogramming of gene expression. Alternative pre-mRNA splicing provides an important layer of this regulation by controlling the level of factors involved in stress response and generating additional protein isoforms with specific features. Recent high-throughput studies have revealed that several defence genes undergo alternative splicing that is often affected by pathogen infection. Despite extensive work, the exact mechanisms underlying these relationships are still unclear, but the contribution of alternative protein isoforms to the defence response and the role of regulatory factors, including components of the splicing machinery, have been established. Modulation of gene expression in response to stress includes alternative splicing, chromatin remodelling, histone modifications, and nucleosome occupancy. How these processes affect plant immunity is mostly unknown, but these facets open new regulatory possibilities. Here we provide an overview of the current state of knowledge and recent findings regarding the growing importance of alternative splicing in plant response to biotic stress.

The Arabidopsis cyclophilin CYP18-1 facilitates PRP18 dephosphorylation and the splicing of introns retained under heat stress

In plants, heat stress induces changes in alternative splicing, including intron retention; these events can rapidly alter proteins or downregulate protein activity, producing nonfunctional isoforms or inducing nonsense-mediated decay of messenger RNA (mRNA). Nuclear cyclophilins (CYPs) are accessory proteins in the spliceosome complexes of multicellular eukaryotes. However, whether plant CYPs are involved in pre-mRNA splicing remain unknown. Here, we found that Arabidopsis thaliana CYP18-1 is necessary for the efficient removal of introns that are retained in response to heat stress during germination. CYP18-1 interacts with Step II splicing factors (PRP18a, PRP22, and SWELLMAP1) and associates with the U2 and U5 small nuclear RNAs in response to heat stress. CYP18-1 binds to phospho-PRP18a, and increasing concentrations of CYP18-1 are associated with increasing dephosphorylation of PRP18a. Furthermore, interaction and protoplast transfection assays revealed that CYP18-1 and the PP2A-type phosphatase PP2A B'η co-regulate PRP18a dephosphorylation. RNA-seq and RT-qPCR analysis confirmed that CYP18-1 is essential for splicing introns that are retained under heat stress. Overall, we reveal the mechanism of action by which CYP18-1 activates the dephosphorylation of PRP18 and show that CYP18-1 is crucial for the efficient splicing of retained introns and rapid responses to heat stress in plants.

Stochastic Variation in DNA Methylation Modulates Nucleosome Occupancy and Alternative Splicing in Arabidopsis thaliana

Plants use complex gene regulatory mechanisms to overcome diverse environmental challenges. For instance, cold stress induces rapid and massive transcriptome changes via alternative splicing (AS) to confer cold tolerance in plants. In mammals, mounting evidence suggests chromatin structure can regulate co-transcriptional AS. Recent evidence also supports co-transcriptional regulation of AS in plants, but how dynamic changes in DNA methylation and the chromatin structure influence the AS process upon cold stress remains poorly understood. In this study, we used the DNA methylation inhibitor 5-Aza-2′-Deoxycytidine (5-aza-dC) to investigate the role of stochastic variations in DNA methylation and nucleosome occupancy in modulating cold-induced AS, in Arabidopsis thaliana (Arabidopsis). Our results demonstrate that 5-aza-dC derived stochastic hypomethylation modulates nucleosome occupancy and AS profiles of genes implicated in RNA metabolism, plant hormone signal transduction, and of cold-related genes in response to cold stress. We also demonstrate that cold-induced remodelling of DNA methylation regulates genes involved in amino acid metabolism. Collectively, we demonstrate that sudden changes in DNA methylation via drug treatment can influence nucleosome occupancy levels and modulate AS in a temperature-dependent manner to regulate plant metabolism and physiological stress adaptation.

Pup1 QTL Regulates Gene Expression Through Epigenetic Modification of DNA Under Phosphate Starvation Stress in Rice

Cytosine methylation, epigenetic DNA modification, is well known to regulate gene expression. Among the epigenetic modifications, 5-methylcytosine (5-mC) has been one of the extensively studied epigenetic changes responsible for regulating gene expression in animals and plants. Though a dramatic change in 5-mC content is observed at the genome level, the variation in gene expression is generally less than that it is expected. Only less is understood about the significance of 5-mC in gene regulation under P-starvation stress in plants. Using whole-genome bisulfite sequencing of a pair of rice [Pusa-44 and its near-isogenic line (NIL)-23 harboring Pup1 QTL] genotypes, we could decipher the role of Pup1 on DNA (de)methylation-mediated regulation of gene expression under P-starvation stress. We observed 13-15% of total cytosines to be methylated in the rice genome, which increased significantly under the stress. The number of differentially methylated regions (DMRs) for hypomethylation (6,068) was higher than those (5,279) for hypermethylated DMRs under the stress, particularly in root of NIL-23. Hypomethylation in CHH context caused upregulated expression of 489 genes in shoot and 382 genes in root of NIL-23 under the stress, wherein 387 genes in shoot and 240 genes in root were upregulated exclusively in NIL-23. Many of the genes for DNA methylation, a few for DNA demethylation, and RNA-directed DNA methylation were upregulated in root of NIL-23 under the stress. Methylation or demethylation of DNA in genic regions differentially affected gene expression. Correlation analysis for the distribution of DMRs and gene expression indicated the regulation of gene mainly through (de)methylation of promoter. Many of the P-responsive genes were hypomethylated or upregulated in roots of NIL-23 under the stress. Hypermethylation of gene body in CG, CHG, and CHH contexts caused up- or downregulated expression of transcription factors (TFs), P transporters, phosphoesterases, retrotransposon proteins, and other proteins. Our integrated transcriptome and methylome analyses revealed an important role of the Pup1 QTL in epigenetic regulation of the genes for transporters, TFs, phosphatases, carbohydrate metabolism, hormone-signaling, and chromatin architecture or epigenetic modifications in P-starvation tolerance. This provides insights into the molecular function of Pup1 in modulating gene expression through DNA (de)methylation, which might be useful in improving P-use efficiency or productivity of rice in P-deficient soil.

Genomically Hardwired Regulation of Gene Activity Orchestrates Cellular Iron Homeostasis in Arabidopsis

Iron (Fe) is an essential micronutrient which plays pivotal roles as electron donor and catalyst across organisms. In plants, variable, often insufficient Fe supply necessitates mechanisms that constantly attune Fe uptake rates and recalibrate cellular Fe homoeostasis. Here, we show that short-term (0.5, 6, and 12 h) exposure of Arabidopsis thaliana plants to Fe deficiency triggered massive changes in gene activity governed by transcription and alternative splicing (AS), regulatory layers that were to a large extent mutually exclusive. Such preclusion was not observed for genes that are directly involved in the acquisition of Fe, which appears to be concordantly regulated by both expression and AS. Generally, genes with lower splice site strengths and higher intron numbers were more likely to be regulated by AS, no dependence on gene architecture was observed for transcriptionally controlled genes. Conspicuously, specific processes were associated with particular genomic features and biased towards either regulatory mode, suggesting that genomic hardwiring is functionally biased. Early changes in splicing patterns were, in many cases, congruent with later changes in transcript or protein abundance, thus contributing to the pronounced transcriptome-proteome discordance observed in plants.

The pattern of DNA methylation alteration, and its association with the changes of gene expression and alternative splicing during phosphate starvation in tomato

DNA methylation is changed and associates with gene expression alterations in plant response to phosphate starvation (Pi-), a common stress that impacts plant growth and productivity. However, in the horticultural model species Solanum lycopersicum (tomato), the dynamics of DNA methylation and its relationship with changes in gene transcription and alternative splicing (AS) under Pi- are unknown. Here, we performed integrative methylome and transcriptome analyses of tomato seedlings under Pi-deficient and -sufficient conditions. We found Pi- caused a slight increase in the overall methylation level, with millions of differentially methylated cytosines (DmCs) and a few hundred differentially methylated regions (DMRs). We also identified thousands of differentially expressed (DE) and differential AS (DAS) genes induced by Pi-, and found that DmCs were more abundant in non-expressed genes than in DE or DAS genes. Moreover, DNA methylation alterations weakly correlated with transcription changes but not with DAS events, and hyper-CHH-DMRs overlapping with transposable elements (TEs) were enriched in a subset of Pi starvation response (PSR) genes. We propose that changes in DNA methylation may be associated with the differential expression of some PSR genes, but that most of these changes probably control the expression of nearby TEs, rather than directly affecting the transcription or AS of PSR genes. Besides, the pattern of methylation changes upon Pi- may largely be shaped by TE distributions. Together, our study provides comprehensive insights into the association of DNA methylation with gene transcription and AS under Pi- in tomato and may contribute to unveiling novel roles of epigenetic mechanisms in plant stress response.

Reading between the Lines: RNA-seq Data Mining Reveals the Alternative Message of the Rice Leaf Transcriptome in Response to Heat Stress

Rice (Oryza sativa L.) is a major food crop but heat stress affects its yield and grain quality. To identify mechanistic solutions to improve rice yield under rising temperatures, molecular responses of thermotolerance must be understood. Transcriptional and post-transcriptional controls are involved in a wide range of plant environmental responses. Alternative splicing (AS), in particular, is a widespread mechanism impacting the stress defence in plants but it has been completely overlooked in rice genome-wide heat stress studies. In this context, we carried out a robust data mining of publicly available RNA-seq datasets to investigate the extension of heat-induced AS in rice leaves. For this, datasets of interest were subjected to filtering and quality control, followed by accurate transcript-specific quantifications. Powerful differential gene expression (DE) and differential AS (DAS) identified 17,143 and 2162 heat response genes, respectively, many of which are novel. Detailed analysis of DAS genes coding for key regulators of gene expression suggests that AS helps shape transcriptome and proteome diversity in response to heat. The knowledge resulting from this study confirmed a widespread transcriptional and post-transcriptional response to heat stress in plants, and it provided novel candidates for rapidly advancing rice breeding in response to climate change.

Stress-Induced Changes in Alternative Splicing Landscape in Rice: Functional Significance of Splice Isoforms in Stress Tolerance

Improvements in yield and quality of rice are crucial for global food security. However, global rice production is substantially hindered by various biotic and abiotic stresses. Making further improvements in rice yield is a major challenge to the rice research community, which can be accomplished through developing abiotic stress-resilient rice varieties and engineering durable agrochemical-independent pathogen resistance in high-yielding elite rice varieties. This, in turn, needs increased understanding of the mechanisms by which stresses affect rice growth and development. Alternative splicing (AS), a post-transcriptional gene regulatory mechanism, allows rapid changes in the transcriptome and can generate novel regulatory mechanisms to confer plasticity to plant growth and development. Mounting evidence indicates that AS has a prominent role in regulating rice growth and development under stress conditions. Several regulatory and structural genes and splicing factors of rice undergo different types of stress-induced AS events, and the functional significance of some of them in stress tolerance has been defined. Both rice and its pathogens use this complex regulatory mechanism to devise strategies against each other. This review covers the current understanding and evidence for the involvement of AS in biotic and abiotic stress-responsive genes, and its relevance to rice growth and development. Furthermore, we discuss implications of AS for the virulence of different rice pathogens and highlight the areas of further research and potential future avenues to develop climate-smart and disease-resistant rice varieties.