Comparative Transcriptome Analysis Reveals the Mechanisms Underlying Differences in Salt Tolerance Between indica and japonica Rice at Seedling Stage

QTL-Seq identified a genomic region on chromosome 1 for soil-salinity tolerance in F2 progeny of Thai salt-tolerant rice donor line "Jao Khao"

Introduction

Owing to advances in high-throughput genome sequencing, QTL-Seq mapping of salt tolerance traits is a major platform for identifying soil-salinity tolerance QTLs to accelerate marker-assisted selection for salt-tolerant rice varieties. We performed QTL-BSA-Seq in the seedling stage of rice from a genetic cross of the extreme salt-sensitive variety, IR29, and "Jao Khao" (JK), a Thai salt-tolerant variety.

Methods

A total of 462 F2 progeny grown in soil and treated with 160 mM NaCl were used as the QTL mapping population. Two high- and low-bulk sets, based on cell membrane stability (CMS) and tiller number at the recovery stage (TN), were equally sampled. The genomes of each pool were sequenced, and statistical significance of QTL was calculated using QTLseq and G prime (G') analysis, which is based on calculating the allele frequency differences or Δ(SNP index).

Results

Both methods detected the overlapping interval region, wherein CMS-bulk was mapped at two loci in the 38.41-38.85 Mb region with 336 SNPs on chromosome 1 (qCMS1) and the 26.13-26.80 Mb region with 1,011 SNPs on chromosome 3 (qCMS3); the Δ(SNP index) peaks were -0.2709 and 0.3127, respectively. TN-bulk was mapped at only one locus in the overlapping 38.26-38.95 Mb region on chromosome 1 with 575 SNPs (qTN1) and a Δ(SNP index) peak of -0.3544. These identified QTLs in two different genetic backgrounds of segregating populations derived from JK were validated. The results confirmed the colocalization of the qCMS1 and qTN1 traits on chromosome 1. Based on the CMS trait, qCMS1/qTN1 stably expressed 6%-18% of the phenotypic variance in the two validation populations, while qCMS1/qTN1 accounted for 16%-20% of the phenotypic variance in one validation population based on the TN trait.

Conclusion

The findings confirm that the CMS and TN traits are tightly linked to the long arm of chromosome 1 rather than to chromosome 3. The validated qCMS-TN1 QTL can be used for gene/QTL pyramiding in marker-assisted selection to expedite breeding for salt resistance in rice at the seedling stage.

Alternative splicing as a driver of natural variation in abscisic acid response

Abscisic acid (ABA) is a crucial player in plant responses to the environment. It accumulates under stress, activating downstream signaling to implement molecular responses that restore homeostasis. Natural variance in ABA sensitivity remains barely understood, and the ABA pathway has been mainly studied at the transcriptional level, despite evidence that posttranscriptional regulation, namely, via alternative splicing, contributes to plant stress tolerance. Here, we identified the Arabidopsis accession Kn-0 as less sensitive to ABA than the reference Col-0, as shown by reduced effects of the hormone on seedling establishment, root branching, and stomatal closure, as well as by decreased induction of ABA marker genes. An in-depth comparative transcriptome analysis of the ABA response in the two variants revealed lower expression changes and fewer genes affected for the least ABA-sensitive ecotype. Notably, Kn-0 exhibited reduced levels of the ABA-signaling SnRK2 protein kinases and lower basal expression of ABA-reactivation genes, consistent with our finding that Kn-0 contains less endogenous ABA than Col-0. ABA also markedly affected alternative splicing, primarily intron retention, with Kn-0 being less responsive regarding both the number and magnitude of alternative splicing events, particularly exon skipping. We find that alternative splicing introduces a more ecotype-specific layer of ABA regulation and identify ABA-responsive splicing changes in key ABA pathway regulators that provide a functional and mechanistic link to the differential sensitivity of the two ecotypes. Our results offer new insight into the natural variation of ABA responses and corroborate a key role for alternative splicing in implementing ABA-mediated stress responses.

Transcriptomic analysis of salt-tolerant and sensitive high-yield japonica rice (Oryza sativa L.) reveals complicated salt-tolerant mechanisms

Developing and cultivating rice varieties is a potent strategy for reclaiming salinity-affected soils for rice production. Nevertheless, the molecular mechanisms conferring salt tolerance, especially in conventional high-yield japonica rice varieties, remain obscure. In this study, Zhendao 23309 (ZD23309) exhibited significantly less grain yield reduction under a salt stress gradient than the control variety Wuyunjing 30 (WYJ30). High positive correlations between grain yield and dry matter accumulation at the jointing, heading and maturity stages indicated that early salt tolerance performance is a crucial hallmark for yield formation. After a mild salt stress (85 mM NaCl) of young seedlings, RNA sequencing (RNA-seq) of shoot and root separately identified a total of 1952 and 3647 differentially expressed genes (DEGs) in ZD23309, and 2114 and 2711 DEGs in WYJ30, respectively. Gene ontology (GO) analysis revealed numerous DEGs in ZD23309 that play pivotal roles in strengthening salt tolerance, encompassing the response to stimulus (GO:0050896) in shoots and nucleoside binding (GO:0001882) in roots. Additionally, distinct expression patterns were observed in a fraction of genes in the two rice varieties under salt stress, corroborating the efficacy of previously reported salt tolerance genes. Our research not only offers fresh insights into the differences in salt stress tolerance among conventional high-yield rice varieties but also unveils the intricate nature of salt tolerance mechanisms. These findings lay a solid groundwork for deciphering the mechanisms underlying salt tolerance.

The comparative transcriptome analysis of two green super rice genotypes with varying tolerance to salt stress

BACKGROUND

Salinity is one of the main abiotic factors that restrict plant growth, physiology, and crop productivity is salt stress. About 33% of the total irrigated land suffers from severe salinity because of intensive underground water extraction and irrigation with brackish water. Thus, it is important to understand the genetic mechanism and identify the novel genes involved in salt tolerance for the development of climate-resilient rice cultivars.

METHODS AND RESULTS

In this study, two rice genotypes with varying tolerance to salt stress were used to investigate the differential expressed genes and molecular pathways to adapt under saline soil by comparative RNA sequencing at 42 days of the seedling stage. Salt-susceptible (S3) and -tolerant (S13) genotypes revealed 3982 and 3463 differentially expressed genes in S3 and S13 genotypes. The up-regulated genes in both genotypes were substantially enriched in different metabolic processes and binding activities. Biosynthesis of secondary metabolites, phenylpropanoid biosynthesis, and plant signal transduction mechanisms were highly enriched. Salt-susceptible and -tolerant genotypes shared the same salt adaptability mechanism with no significant quantitative differences at the transcriptome level. Moreover, bHLH, ERF, NAC, WRKY, and MYB transcription factors were substantially up-regulated under salt stress. 391 out of 1806 identified novel genes involved in signal transduction mechanisms. Expression profiling of six novel genes further validated the findings from RNA-seq data.

CONCLUSION

These findings suggest that the differentially expressed genes and molecular mechanisms involved in salt stress adaptation are conserved in both salt-susceptible and salt-tolerant rice genotypes. Further molecular characterization of novel genes will help to understand the genetic mechanism underlying salt tolerance in rice.

Comparative transcriptome analysis of gene responses of salt-tolerant and salt-sensitive rice cultivars to salt stress

Salt stress is one unfavorable factor of global climate change that adversely affects rice plant growth and yield. To identify novel salt-tolerant genes and new varieties of salt-tolerant rice, a better understanding of the molecular regulation mechanism of salt tolerance in rice is needed. In this study we used transcriptome analyses to examine changes in gene expression of salt-tolerant and salt-sensitive rice plants. The salt-tolerant cultivar HH11 and salt-sensitive cultivar IR29 were treated with 200 mM NaCl solution for 0 h, 6 h, 24 h and 48 h at the three leaf stage. Physiological parameters and transcriptome were measured and analyzed after each treatment. Activity of SOD and POD, as well as the MDA and protein content of the two rice cultivars generally increased with increasing time of exposure to NaCl. Meanwhile, the APX activity first increased, then decreased in both cultivars, with maximum values seen at 6 h for IR29 and at 24 h for HH11. The GR and GPX activity of HH11 were stronger than that of IR29 in response to salt stress. The H2O2 content first increased at 0-6 h, then decreased at 6-24 h, and then increased again at 24-48 h under salt stress. Compared with IR29, SOD, POD and APX activity of HH11 was more sluggish in response to salt stress, reaching the maximum at 24 h or 48 h. The MDA, H2O2 and proline content of HH11 was lower than that of IR29 under salt stress. Relative to untreated HH11 plants (0 h) and those exposed to salt for 6 h, 24 h, and 48 h (H0-H6, H0-H24 and H0-H48), 7462, 6363 and 6636, differentially expressed genes (DEGs), respectively, were identified. For IR29, the respective total DEGs were 7566, 6075 and 6136. GO and KEGG enrichment analysis showed that metabolic pathways related to antioxidative responses and osmotic balance played vital roles in salt stress tolerance. Sucrose and starch metabolism, in addition to flavonoid biosynthesis and glutathione metabolism, showed positive responses to salt stress. Expression of two SPS genes (LOC_Os01g69030 and LOC_Os08g20660) and two GST genes (LOC_Os06g12290 and LOC_Os10g38740) was up-regulated in both HH11 and IR29, whereas expression of LOC_Os09g12660, a glucose-1-phosphate adenylyltransferase gene, and two SS genes (LOC_Os04g17650 and LOC_Os04g24430) was up-regulated differential expression in HH11. The results showed that HH11 had more favorable adjustment in antioxidant and osmotic activity than IR29 upon exposure to salt stress, and highlighted candidate genes that could play roles in the function and regulation mechanism of salt tolerance in rice.

Integrated Transcriptomic and Metabolomic Analyses Uncover the Differential Mechanism in Saline-Alkaline Tolerance between Indica and Japonica Rice at the Seedling Stage

Saline-alkaline stress is one of the major damages that severely affects rice (Oryza sativa L.) growth and grain yield; however, the mechanism of the tolerance remains largely unknown in rice. Herein, we comparatively investigated the transcriptome and metabolome of two contrasting rice subspecies genotypes, Luohui 9 (abbreviation for Chao2R under study, O. sativa ssp. indica, saline-alkaline-sensitive) and RPY geng (O. sativa ssp. japonica, saline-alkaline-tolerant), to identify the main pathways and important factors related to saline-alkaline tolerance. Transcriptome analysis showed that 68 genes involved in fatty acid, amino acid (such as phenylalanine and tryptophan), phenylpropanoid biosynthesis, energy metabolism (such as Glycolysis and TCA cycle), as well as signal transduction (such as hormone and MAPK signaling) were identified to be specifically upregulated in RPY geng under saline-alkaline conditions, implying that a series of cascade changes from these genes promotes saline-alkaline stress tolerance. The transcriptome changes observed in RPY geng were in high accordance with the specifically accumulation of metabolites, consisting mainly of 14 phenolic acids, 8 alkaloids, and 19 lipids based on the combination analysis of transcriptome and metabolome. Moreover, some genes involved in signal transduction as hub genes, such as PR5, FLS2, BRI1, and NAC, may participate in the saline-alkaline stress response of RPY geng by modulating key genes involved in fatty acid, phenylpropanoid biosynthesis, amino acid metabolism, and glycolysis metabolic pathways based on the gene co-expression network analysis. The present research results not only provide important insights for understanding the mechanism underlying of rice saline-alkaline tolerance at the transcriptome and metabolome levels but also provide key candidate target genes for further enhancing rice saline-alkaline stress tolerance.

Role of non-coding RNAs against salinity stress in Oryza species: Strategies and challenges in analyzing miRNAs, tRFs and circRNAs

Salinity is an imbalanced concentration of mineral salts in the soil or water that causes yield loss in salt-sensitive crops. Rice plant is vulnerable to soil salinity stress at seedling and reproductive stages. Different non-coding RNAs (ncRNAs) post-transcriptionally regulate different sets of genes during different developmental stages under varying salinity tolerance levels. While microRNAs (miRNAs) are well known small endogenous ncRNAs, tRNA-derived RNA fragments (tRFs) are an emerging class of small ncRNAs derived from tRNA genes with a demonstrated regulatory role, like miRNAs, in humans but unexplored in plants. Circular RNA (circRNA), another ncRNA produced by back-splicing events, acts as target mimics by preventing miRNAs from binding with their target mRNAs, thereby reducing the miRNA's action upon its target. Same may hold true between circRNAs and tRFs. Hence, the work done on these ncRNAs was reviewed and no reports were found for circRNAs and tRFs under salinity stress in rice, either at seedling or reproductive stages. Even the reports on miRNAs are restricted to seedling stage only, in spite of severe effects on rice crop production due to salt stress during reproductive stage. Moreover, this review sheds light on strategies to predict and analyze these ncRNAs in an effective manner.

Genetic Dissection of Alkalinity Tolerance at the Seedling Stage in Rice (Oryza sativa) Using a High-Resolution Linkage Map

Although both salinity and alkalinity result from accumulation of soluble salts in soil, high pH and ionic imbalance make alkaline stress more harmful to plants. This study aimed to provide molecular insights into the alkalinity tolerance using a recombinant inbred line (RIL) population developed from a cross between Cocodrie and Dular with contrasting response to alkalinity stress. Forty-six additive QTLs for nine morpho-physiological traits were mapped on to a linkage map of 4679 SNPs under alkalinity stress at the seedling stage and seven major-effect QTLs were for alkalinity tolerance scoring, Na⁺ and K⁺ concentrations and Na⁺:K⁺ ratio. The candidate genes were identified based on the comparison of the impacts of variants of genes present in five QTL intervals using the whole genome sequences of both parents. Differential expression of no apical meristem protein, cysteine protease precursor, retrotransposon protein, OsWAK28, MYB transcription factor, protein kinase, ubiquitin-carboxyl protein, and NAD binding protein genes in parents indicated their role in response to alkali stress. Our study suggests that the genetic basis of tolerance to alkalinity stress is most likely different from that of salinity stress. Introgression and validation of the QTLs and genes can be useful for improving alkalinity tolerance in rice at the seedling stage and advancing understanding of the molecular genetic basis of alkalinity stress adaptation.

Combined transcriptome and metabolome reveal glutathione metabolism plays a critical role in resistance to salinity in rice landraces HD961

Rice (Oryza sativa) is one of the most important food crops around the world, which is sensitive to salt stress, especially in the seedling and booting stage. HD961 is a salt-tolerant rice landrace that grows along coastal beaches and has disease and insect pest resistance, salt tolerance, and vigorous growth characteristics. We performed a combined transcriptome and metabolome analysis to clarify salinity resistance mechanisms in cultivar HD961, which has adapted to salinity soil at the early seedling stage. The results showed that the growth and antioxidant capacity of HD961 were stronger than 9311 under salt stress (SS). Transcriptomic analysis showed that a total of 6,145, 3,309, 1,819, and 1,296 differentially expressed genes (DEGs) were identified in the groups of TH60 (control group vs. 60 mM group of HD961 for transcriptome), TH120 (control group vs. 120 mM group of HD961 for transcriptome), T60 (control group vs. 60 mM group of 9311 for transcriptome), and T120 (control group vs. 120 mM group of 9311 for transcriptome), respectively. Starch and sucrose metabolism and phenylpropanoid biosynthesis were shared in the four treatment groups based on a KEGG enrichment analysis of DEGs. In addition, alpha-linolenic acid metabolism, plant hormone signal transduction, plant-pathogen interaction, and fatty acid elongation were specific and significantly different in HD961. A total of 92, 158, 151, and 179 significantly regulated metabolites (SRMs) responded to SS in MH60 (control group vs. 60 mM group of HD961 for metabolome), MH120 (control group vs. 120 mM group of HD961 for metabolome), M60 (control group vs. 60 mM group of 9311 for metabolome), and M120 (control group vs. 120 mM group of 9311 for metabolome), respectively. The KEGG analysis showed that eight common metabolic pathways were identified in the four treatment groups, of which biosynthesis of amino acids was the most significant. Three specific metabolic pathways were identified in the HD961, including glutathione metabolism, ascorbate and aldarate metabolism, and pantothenate and CoA biosynthesis. Integrative analysis between the transcriptome and metabolome showed that glutathione metabolism was specific and significantly affected under SS in HD961. A total of seven SRMs and 48 DEGs and four SRMs and 15 DEGs were identified in the glutathione metabolism pathway in HD961 and 9311, respectively. The Pearson correlation analysis showed a significant correlation between reduced glutathione and 16 genes (12 upregulated and four downregulated genes), suggesting these genes could be candidates as salt-tolerance regulation genes. Collectively, our data show that glutathione metabolism plays a critical role in response to SS in rice. Moreover, the stronger regulative ability of related common genes and metabolites might contribute to salt resistance in HD961.

De novo assembly of two chromosome-level rice genomes and bin-based QTL mapping reveal genetic diversity of grain weight trait in rice

Following the "green revolution," indica and japonica hybrid breeding has been recognized as a new breakthrough in further improving rice yields. However, heterosis-related grain weight QTLs and the basis of yield advantage among subspecies has not been well elucidated. We herein de novo assembled the chromosome level genomes of an indica/xian rice (Luohui 9) and a japonica/geng rice (RPY geng) and found that gene number differences and structural variations between these two genomes contribute to the differences in agronomic traits and also provide two different favorable allele pools to produce better derived recombinant inbred lines (RILs). In addition, we generated a high-generation (> F₁₅) population of 272 RILs from the cross between Luohui 9 and RPY geng and two testcross hybrid populations derived from the crosses of RILs and two cytoplasmic male sterile lines (YTA, indica and Z7A, japonica). Based on three derived populations, we totally identified eight 1,000-grain weight (KGW) QTLs and eight KGW heterosis loci. Of QTLs, qKGW-6.1 and qKGW-8.1 were accepted as novel KGW QTLs that have not been reported previously. Interestingly, allele genotyping results revealed that heading date related gene (Ghd8) in qKGW-8.1 and qLH-KGW-8.1, can affect grain weight in RILs and rice core accessions and may also play an important role in grain weight heterosis. Our results provided two high-quality genomes and novel gene editing targets for grain weight for future rice yield improvement project.

Genome-Wide Identification and Evolution Analysis of the Gibberellin Oxidase Gene Family in Six Gramineae Crops

The plant hormones gibberellins (GAs) regulate plant growth and development and are closely related to the yield of cash crops. The GA oxidases (GAoxs), including the GA2ox, GA3ox, and GA20ox subfamilies, play pivotal roles in GAs' biosynthesis and metabolism, but their classification and evolutionary pattern in Gramineae crops remain unclear. We thus conducted a comparative genomic study of GAox genes in six Gramineae representative crops, namely, Setaria italica (Si), Zea mays (Zm), Sorghum bicolor (Sb), Hordeum vulgare (Hv), Brachypodium distachyon (Bd), and Oryza sativa (Os). A total of 105 GAox genes were identified in these six crop genomes, belonging to the C19-GA2ox, C20-GA2ox, GA3ox, and GA20ox subfamilies. Based on orthogroup (OG) analysis, GAox genes were divided into nine OGs and the number of GAox genes in each of the OGs was similar among all tested crops, which indicated that GAox genes may have completed their family differentiations before the species differentiations of the tested species. The motif composition of GAox proteins showed that motifs 1, 2, 4, and 5, forming the 2OG-FeII_Oxy domain, were conserved in all identified GAox protein sequences, while motifs 11, 14, and 15 existed specifically in the GA20ox, C19-GA2ox, and C20-GA2ox protein sequences. Subsequently, the results of gene duplication events suggested that GAox genes mainly expanded in the form of WGD/SD and underwent purification selection and that maize had more GAox genes than other species due to its recent duplication events. The cis-acting elements analysis indicated that GAox genes may respond to growth and development, stress, hormones, and light signals. Moreover, the expression profiles of rice and maize showed that GAox genes were predominantly expressed in the panicles of the above two plants and the expression of several GAox genes was significantly induced by salt or cold stresses. In conclusion, our results provided further insight into GAox genes' evolutionary differences among six representative Gramineae and highlighted GAox genes that may play a role in abiotic stress.

Pan-transcriptome assembly combined with multiple association analysis provides new insights into the regulatory network of specialized metabolites in the tea plant Camellia sinensis

Specialized metabolites not only play important roles in biotic and abiotic stress adaptation of tea plants (Camellia sinensis (L.) O. Kuntze) but also contribute to the unique flavor of tea, the most important nonalcoholic beverage. However, the molecular networks and major genes that regulate specialized metabolites in tea plants are not well understood. Here, we constructed a population-level pan-transcriptome of the tea plant leaf using second-leaf transcriptome data from 134 accessions to investigate global expression differences in the population, expression presence or absence variations (ePAVs), and differentially expressed genes (DEGs) between pure Camellia sinensis var. assamica (CSA) and pure Camellia sinensis var. sinensis (CSS) accessions. Next, we used a genome-wide association study, a quantitative trait transcript study, and a transcriptome-wide association study to integrate genotypes, accumulation levels of specialized metabolites, and expression levels of pan-transcriptome genes to identify candidate regulatory genes for flavor-related metabolites and to construct a regulatory network for specialized metabolites in tea plants. The pan-transcriptome contains 30 482 expressed genes, 4940 and 5506 of which were newly annotated from a de novo transcriptome assembly without a reference and a genome reference-based assembly, respectively. DEGs and ePAVs indicated that CSA and CSS were clearly differentiated at the population transcriptome level, and they were closely related to abiotic tolerance and secondary metabolite synthesis phenotypes of CSA and CSS based on gene annotations. The regulatory network contained 212 specialized metabolites, 3843 candidate genes, and 3407 eQTLs, highlighting many pleiotropic candidate genes, candidate gene-rich eQTLs, and potential regulators of specialized metabolites. These included important transcription factors in the AP2/ERF-ERF, MYB, WD40, and bHLH families. CsTGY14G0001296, an ortholog of AtANS, appeared to be directly related to variation in proanthocyanins in the tea plant population, and the CsTGY11G0002074 gene encoding F3'5'H was found to contribute to the biased distribution of catechins between pure CSAs and pure CSSs. Together, these results provide a new understanding of the metabolite diversity in tea plants and offer new insights for more effective breeding of better-flavored tea varieties.