Allelic variants of OsSUB1A cause differential expression of transcription factor genes in response to submergence in rice

Integrated meta-analysis and transcriptomics pinpoint genomic loci and novel candidate genes associated with submergence tolerance in rice

BACKGROUND

Due to rising costs, water shortages, and labour shortages, farmers across the globe now prefer a direct seeding approach. However, submergence stress remains a major bottleneck limiting the success of this approach in rice cultivation. The merger of accumulated rice genetic resources provides an opportunity to detect key genomic loci and candidate genes that influence the flooding tolerance of rice.

RESULTS

In the present study, a whole-genome meta-analysis was conducted on 120 quantitative trait loci (QTL) obtained from 16 independent QTL studies reported from 2004 to 2023. These QTL were confined to 18 meta-QTL (MQTL), and ten MQTL were successfully validated by independent genome-wide association studies from diverse natural populations. The mean confidence interval (CI) of the identified MQTL was 3.44 times narrower than the mean CI of the initial QTL. Moreover, four core MQTL loci with genetic distance less than 2 cM were obtained. By combining differentially expressed genes (DEG) from two transcriptome datasets with 858 candidate genes identified in the core MQTL regions, we found 38 common differentially expressed candidate genes (DECGs). In silico expression analysis of these DECGs led to the identification of 21 genes with high expression in embryo and coleoptile under submerged conditions. These DECGs encode proteins with known functions involved in submergence tolerance including WRKY, F-box, zinc fingers, glycosyltransferase, protein kinase, cytochrome P450, PP2C, hypoxia-responsive family, and DUF domain. By haplotype analysis, the 21 DECGs demonstrated distinct genetic differentiation and substantial genetic distance mainly between indica and japonica subspecies. Further, the MQTL7.1 was successfully validated using flanked marker S2329 on a set of genotypes with phenotypic variation.

CONCLUSION

This study provides a new perspective on understanding the genetic basis of submergence tolerance in rice. The identified MQTL and novel candidate genes lay the foundation for marker-assisted breeding/engineering of flooding-tolerant cultivars conducive to direct seeding.

Physiological responses and Ethylene-Response AP2/ERF Factor expression in Indica rice seedlings subjected to submergence and osmotic stress

BACKGROUND

The increased frequency of heavy rains in recent years has led to submergence stress in rice paddies, severely affecting rice production. Submergence causes not only hypoxic stress from excess water in the surrounding environment but also osmotic stress in plant cells. We assessed physiological responses and Ethylene-Response AP2/ERF Factor regulation under submergence conditions alone and with ionic or nonionic osmotic stress in submergence-sensitive IR64 and submergence-tolerant IR64-Sub1 Indica rice cultivars.

RESULTS

Our results indicate that both IR64 and IR64-Sub1 exhibited shorter plant heights and root lengths under submergence with nonionic osmotic stress than normal condition and submergence alone. IR64-Sub1 seedlings exhibited a significantly lower plant height under submergence conditions alone and with ionic or nonionic osmotic stress than IR64 cultivars. IR64-Sub1 seedlings also presented lower malondialdehyde (MDA) concentration and higher survival rates than did IR64 seedlings after submergence with ionic or nonionic osmotic stress treatment. Sub1A-1 affects reactive oxygen species (ROS) accumulation and antioxidant enzyme activity in rice. The results also show that hypoxia-inducible ethylene response factors (ERF)-VII group and alcohol dehydrogenase 1 (ADH1) and lactate dehydrogenase 1 (LDH1) genes exhibited different expression levels under nonionic or ionic osmotic stress during submergence on rice.

CONCLUSIONS

Together, these results demonstrate that complex regulatory mechanisms are involved in responses to the aforementioned forms of stress and offer new insights into the effects of submergence and osmotic stress on rice.

Comprehensive transcriptomic analysis of three varieties with different brown planthopper-resistance identifies leaf sheath lncRNAs in rice

BACKGROUND

Long non-coding RNAs (lncRNAs) have been brought great attention for their crucial roles in diverse biological processes. However, systematic identification of lncRNAs associated with specialized rice pest, brown planthopper (BPH), defense in rice remains unexplored.

RESULTS

In this study, a genome-wide high throughput sequencing analysis was performed using leaf sheaths of susceptible rice Taichung Native 1 (TN1) and resistant rice IR36 and R476 with and without BPH feeding. A total of 2283 lncRNAs were identified, of which 649 lncRNAs were differentially expressed. During BPH infestation, 84 (120 in total), 52 (70 in total) and 63 (94 in total) of differentially expressed lncRNAs were found only in TN1, IR36 and R476, respectively. Through analyzing their cis-, trans-, and target mimic-activities, not only the lncRNAs targeting resistance genes (NBS-LRR and RLKs) and transcription factors, but also the lncRNAs acting as the targets of the well-studied stress-related miRNAs (miR2118, miR528, and miR1320) in each variety were identified. Before the BPH feeding, 238 and 312 lncRNAs were found to be differentially expressed in TN1 vs. IR36 and TN1 vs. R476, respectively. Among their putative targets, the plant-pathogen interaction pathway was significantly enriched. It is speculated that the resistant rice was in a priming state by the regulation of lncRNAs. Furthermore, the lncRNAs extensively involved in response to BPH feeding were identified by Weighted Gene Co-expression Network Analysis (WGCNA), and the possible regulation networks of the key lncRNAs were constructed. These lncRNAs regulate different pathways that contribute to the basal defense and specific resistance of rice to the BPH.

CONCLUSION

In summary, we identified the specific lncRNAs targeting the well-studied stress-related miRNAs, resistance genes, and transcription factors in each variety during BPH infestation. Additionally, the possible regulating network of the lncRNAs extensively responding to BPH feeding revealed by WGCNA were constructed. These findings will provide further understanding of the regulatory roles of lncRNAs in BPH defense, and lay a foundation for functional research on the candidate lncRNAs.

A NAC transcription factor OsNAC3 positively regulates ABA response and salt tolerance in rice

BACKGROUND

NAC (NAM, ATAF and CUC) transcription factors (TFs) play vital roles in plant development and abiotic stress tolerance. Salt stress is one of the most limiting factors for rice growth and production. However, the mechanism underlying salt tolerance in rice is still poorly understood.

RESULTS

In this study, we functionally characterized a rice NAC TF OsNAC3 for its involvement in ABA response and salt tolerance. ABA and NaCl treatment induced OsNAC3 expression in roots. Immunostaining showed that OsNAC3 was localized in all root cells. OsNAC3 knockout decreased rice plants' sensitivity to ABA but increased salt stress sensitivity, while OsNAC3 overexpression showed an opposite effect. Loss of OsNAC3 also induced Na⁺ accumulation in the shoots. Furthermore, qRT-PCR and transcriptomic analysis were performed to identify the key OsNAC3 regulated genes related to ABA response and salt tolerance, such as OsHKT1;4, OsHKT1;5, OsLEA3-1, OsPM-1, OsPP2C68, and OsRAB-21.

CONCLUSIONS

This study shows that rice OsNAC3 is an important regulatory factor in ABA signal response and salt tolerance.

Promoter Architecture and Transcriptional Regulation of Genes Upregulated in Germination and Coleoptile Elongation of Diverse Rice Genotypes Tolerant to Submergence

Rice has the natural morphological adaptation to germinate and elongate its coleoptile under submerged flooding conditions. The phenotypic deviation associated with the tolerance to submergence at the germination stage could be due to natural variation. However, the molecular basis of this variation is still largely unknown. A comprehensive understanding of gene regulation of different genotypes that have diverse rates of coleoptile elongation can provide significant insights into improved rice varieties. To do so, publicly available transcriptome data of five rice genotypes, which have different lengths of coleoptile elongation under submergence tolerance, were analyzed. The aim was to identify the correlation between promoter architecture, associated with transcriptional and hormonal regulation, in diverse genotype groups of rice that have different rates of coleoptile elongation. This was achieved by identifying the putative cis-elements present in the promoter sequences of genes upregulated in each group of genotypes (tolerant, highly tolerant, and extremely tolerant genotypes). Promoter analysis identified transcription factors (TFs) that are common and unique to each group of genotypes. The candidate TFs that are common in all genotypes are MYB, bZIP, AP2/ERF, ARF, WRKY, ZnF, MADS-box, NAC, AS2, DOF, E2F, ARR-B, and HSF. However, the highly tolerant genotypes interestingly possess binding sites associated with HY5 (bZIP), GBF3, GBF4 and GBF5 (bZIP), DPBF-3 (bZIP), ABF2, ABI5, bHLH, and BES/BZR, in addition to the common TFs. Besides, the extremely tolerant genotypes possess binding sites associated with bHLH TFs such as BEE2, BIM1, BIM3, BM8 and BAM8, and ABF1, in addition to the TFs identified in the tolerant and highly tolerant genotypes. The transcriptional regulation of these TFs could be linked to phenotypic variation in coleoptile elongation in response to submergence tolerance. Moreover, the results indicate a cross-talk between the key TFs and phytohormones such as gibberellic acid, abscisic acid, ethylene, auxin, jasmonic acid, and brassinosteroids, for an altered transcriptional regulation leading to differences in germination and coleoptile elongation under submergence. The information derived from the current in silico analysis can potentially assist in developing new rice breeding targets for direct seeding.

Flood resilience loci SUBMERGENCE 1 and ANAEROBIC GERMINATION 1 interact in seedlings established underwater

Crops with resilience to multiple climatic stresses are essential for increased yield stability. Here, we evaluate the interaction between two loci associated with flooding survival in rice (Oryza sativa L.). ANAEROBIC GERMINATION 1 (AG1), encoding trehalose 6-phosphate phosphatase 7 (TPP7), promotes mobilization of endosperm reserves to enhance the elongation of a hollow coleoptile in seeds that are seeded directly into shallow paddies. SUBMERGENCE 1 (SUB1), encoding the ethylene-responsive transcription factor SUB1A-1, confers tolerance to complete submergence by dampening carbohydrate catabolism, to enhance recovery upon desubmergence. Interactions between AG1/TPP7 and SUB1/SUB1A-1 were investigated under three flooding scenarios using four near-isogenic lines by surveying growth and survival. Pyramiding of the two loci does not negatively affect anaerobic germination or vegetative-stage submergence tolerance. However, the pyramided AG1 SUB1 genotype displays reduced survival when seeds are planted underwater and maintained under submergence for 16 d. To better understand the roles of TPP7 and SUB1A-1 and their interaction, temporal changes in carbohydrates and shoot transcriptomes were monitored in the four genotypes varying at the two loci at four developmental timeponts, from day 2 after seeding through day 14 of complete submergence. TPP7 enhances early coleoptile elongation, whereas SUB1A-1 promotes precocious photoautotrophy and then restricts underwater elongation. By contrast, pyramiding of the AG1 and SUB1 slows elongation growth, the transition to photoautotrophy, and survival. mRNA-sequencing highlights time-dependent and genotype-specific regulation of mRNAs associated with DNA repair, cell cycle, chromatin modification, plastid biogenesis, carbohydrate catabolism and transport, elongation growth, and other processes. These results suggest that interactions between AG1/TPP7 and SUB1/SUB1A-1 could impact seedling establishment if paddy depth is not effectively managed after direct seeding.

Identification of main-effect quantitative trait loci (QTLs) for low-temperature stress tolerance germination- and early seedling vigor-related traits in rice (Oryza sativa L.)

An attempt was made in the current study to identify the main-effect and co-localized quantitative trait loci (QTLs) for germination and early seedling growth traits under low-temperature stress (LTS) conditions in rice. The plant material used in this study was an early backcross population of 230 introgression lines (ILs) in BCIF7 generation derived from the Weed Tolerant Rice-1 (WTR-1) (as the recipient) and Haoannong (HNG) (as the donor). Genetic analyses of LTS tolerance revealed a total of 27 main-effect quantitative trait loci (M-QTLs) mapped on 12 chromosomes. These QTLs explained more than 10% of phenotypic variance (PV), and average PV of 12.71% while employing 704 high-quality SNP markers. Of these 27 QTLs distributed on 12 chromosomes, 11 were associated with low-temperature germination (LTG), nine with low-temperature germination stress index (LTGS), five with root length stress index (RLSI), and two with biomass stress index (BMSI) QTLs, shoot length stress index (SLSI) and root length stress index (RLSI), seven with seed vigor index (SVI), and single QTL with root length (RL). Among them, five significant major QTLs (qLTG(I) 1 , qLTGS(I) 1-2 , qLTG(I) 5 , qLTGS(I) 5 , and qLTG(I) 7 ) mapped on chromosomes 1, 5, and 7 were associated with LTG and LTGS traits and the PV explained ranged from 16 to 23.3%. The genomic regions of these QTLs were co-localized with two to six QTLs. Most of the QTLs were growth stage-specific and found to harbor QTLs governing multiple traits. Eight chromosomes had more than four QTLs and were clustered together and designated as promising LTS tolerance QTLs (qLTTs), as qLTT 1 , qLTT 2 , qLTT 3 , qLTT 5 , qLTT 6 , qLTT 8 , qLTT 9 , and qLTT 11 . A total of 16 putative candidate genes were identified in the major M-QTLs and co-localized QTL regions distributed on different chromosomes. Overall, these significant genomic regions of M-QTLs are responsible for multiple traits and this suggested that these could serve as the best predictors of LTS tolerance at germination and early seedling growth stages. Furthermore, it is necessary to fine-map these regions and to find functional markers for marker-assisted selection in rice breeding programs for cold tolerance.

Analysis of the distribution of assimilation products and the characteristics of transcriptomes in rice by submergence during the ripening stage

Background

Research on the submergence stress of rice has concentrated on the quiescence strategy to survive in long-term flooding conditions based on Submergence-1A (SUB1A). In the case of the ripening period, it is important that submergence stress can affect the quality as well as the survival of rice. Therefore, it is essential to understand the changes in the distribution of assimilation products in grain and ripening characteristics in submergence stress conditions. However, such studies have been insufficient at the physiological and molecular biological levels.

Results

We confirmed that the distribution rate of assimilation products in grain was decreased by submergence treatment. These results were caused by an increase in the distribution rate of assimilation products to the stem according to escape strategy. To understand this phenomenon at the molecular level, we analyzed the relative expression levels of genes related to sucrose metabolism, and found that the sucrose phosphate synthase gene (OsSPS), which induces the accumulation of sucrose in tissues, was decreased in the seeds and leaves, but not in the stems. Furthermore, the sucrose transporter gene (OsSUT) related to sucrose transport decreased in the seeds and leaves, but increased in stems. We also analyzed the biological metabolic processes related to starch and sucrose synthesis, carbon fixation, and glycolysis using the KEGG mapper with selected differentially expressed genes (DEGs) in seeds, stems, and leaves caused by submergence treatment. We found that the expression of genes for each step related to starch and D-glucose synthesis was down-regulated in the seeds and leaves but up-regulated in the stem.

Conclusion

The results of this study provide basic data for the development of varieties and corresponding technologies adapted to submergence conditions, through understanding the action network of the elements that change in the submergence condition, as well as information regarding useful DEGs.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5320-7) contains supplementary material, which is available to authorized users.