Comparative Transcriptomics and Co-Expression Networks Reveal Tissue- and Genotype-Specific Responses of qDTYs to Reproductive-Stage Drought Stress in Rice (Oryza sativa L.)

Panicle transcriptome of high-yield mutant indica rice reveals physiological mechanisms and novel candidate regulatory genes for yield under reproductive stage drought stress

BACKGROUND

Reproductive stage drought stress (RDS) is a major global threat to rice production. Due to climate change, water scarcity is becoming an increasingly common phenomenon in major rice-growing areas worldwide. Understanding RDS mechanisms will allow candidate gene identification to generate novel rice genotypes tolerant to RDS.

RESULTS

To generate novel rice genotypes that can sustain yield under RDS, we performed gamma-irradiation mediated mutation breeding in the drought stress susceptible mega rice variety, MTU1010. One of the mutant MM11 (MTU1010 derived mutant11) shows consistently increased performance in yield-related traits under field conditions consecutively for four generations. In addition, compared to MTU1010, the yield of MM11 is sustained in prolonged drought imposed during the reproductive stage under field and in pot culture conditions. A comparative emerged panicle transcriptome analysis of the MTU1010 and MM11 suggested metabolic adjustment, enhanced photosynthetic ability, and hormone interplay in regulating yield under drought responses during emerged panicle development. Regulatory network analysis revealed few putative significant transcription factor (TF)-target interactions involved in integrated signalling between panicle development, yield and drought stress.

CONCLUSIONS

A gamma-irradiate rice mutant MM11 was identified by mutation breeding, and it showed higher potential to sustain yield under reproductive stage drought stress in field and pot culture conditions. Further, a comparative panicle transcriptome revealed significant biological processes and molecular regulators involved in emerged panicle development, yield and drought stress integration. The study extends our understanding of the physiological mechanisms and candidate genes involved in sustaining yield under drought stress.

Comparative transcriptomics and co-expression networks reveal cultivar-specific molecular signatures associated with reproductive-stage cold stress in rice

The resistance of Huaidao5 results from the high constitutive expression of tolerance genes, while that of Huaidao9 is due to the cold-induced resistance in flag leaves and panicles. The regulation mechanism of rice seedlings' cold tolerance is relatively clear, and knowledge of its underlying mechanisms at the reproductive stage is limited. We performed differential expression and co-expression network analyses to transcriptomes from panicle and flag leaf tissues of a cold-tolerant cultivar (Huaidao5), and a sensitive cultivar (Huaidao9), under reproductive-stage cold stress. The results revealed that the expression levels of genes in stress-related pathways such as MAPK signaling pathway, diterpenoid biosynthesis, glutathione metabolism, plant-pathogen interaction and plant hormone signal transduction were constitutively highly expressed in Huaidao5, especially in panicles. Moreover, the Hudaidao5's panicle sample-specific (under cold) module contained some genes related to rice yield, such as GW5L, GGC2, SG1 and CTPS1. However, the resistance of Huaidao9 was derived from the induced resistance to cold in flag leaves and panicles. In the flag leaves, the responses included a series of stress response and signal transduction, while in the panicles nitrogen metabolism was severely affected, especially 66 endosperm-specific genes. Through integrating differential expression with co-expression networks, we predicted 161 candidate genes (79 cold-responsive genes common to both cultivars and 82 cold-tolerance genes associated with differences in cold tolerance between cultivars) potentially affecting cold response/tolerance, among which 85 (52.80%) were known to be cold-related genes. Moreover, 52 (65.82%) cold-responsive genes (e.g., TIFY11C, LSK1 and LPA) could be confirmed by previous transcriptome studies and 72 (87.80%) cold-tolerance genes (e.g., APX5, OsFbox17 and OsSTA109) were located within QTLs associated with cold tolerance. This study provides an efficient strategy for further discovery of mechanisms of cold tolerance in rice.

Transcriptome analysis of drought-responsive and drought-tolerant mechanisms in maize leaves under drought stress

Maize is a major crop essential for food and feed, but its production is threatened by various biotic and abiotic stresses. Drought is one of the most common abiotic stresses, causing severe crop yield reduction. Although several studies have been devoted to selecting drought-tolerant maize lines and detecting the drought-responsive mechanism of maize, the transcriptomic differences between drought-tolerant and drought-susceptible maize lines are still largely unknown. In our study, RNA-seq was performed on leaves of the drought-tolerant line W9706 and the drought-susceptible line B73 after drought treatment. We identified 3147 differentially expressed genes (DEGs) between these two lines. The upregulated DEGs in W9706 were enriched in specific processes, including ABA signaling, wax biosynthesis, CHO metabolism, signal transduction and brassinosteroid biosynthesis-related processes, while the downregulated DEGs were enriched in specific processes, such as stomatal movement. Altogether, transcriptomic analysis suggests that the different drought resistances were correlated with the differential expression of genes, while the drought tolerance of W9706 is due to the more rapid response to stimulus, higher water retention capacity and stable cellular environment under water deficit conditions.

Transcriptome and anatomical studies reveal alterations in leaf thickness under long-term drought stress in tobacco

Drought is one of the foremost environmental factors that limit the growth of plants. Leaf thickness (LT) is an important quantitative trait in plant physiology. The experiment was carried out in a growth room and the plants were divided into two groups such as well-watered and drought-stressed. This work investigated leaf growth in terms of leaf surface growth and expansion rate, leaf stomata traits, LT, anticlinal growth, and leaf cell layers. The results showed that the leaf area and leaf surface expansion rate were decreased by drought stress (DS). Similarly, LT, anticlinal expansion rate, palisade and spongy tissue thickness, and their related expansion rates were also decreased at different days' time points (DTP) of DS. However, a steady increase was observed in the aforementioned parameters after 12 DTP of DS. The stomatal density increased while stomata size decreased at 3 DTP and 12 DTP (low leaf water potential and relative leaf water content at these time points) and vice versa at 24 DTP compared with the well-watered plants indicating adaptations in these traits in response to DS, and thus the leaf water status played a role in the regulation of leaf stomata traits. The cell length decreased in the upper epidermis, palisade and spongy tissues by DS up to 12 DTP led to lower LT while an increase was observed after 12 DTP that resulted in higher LT. The increase in the LT was supported by the upregulation of starch and sucrose metabolism, glycerolipid metabolism, protein processing in endoplasmic reticulum pathways at 18 DTP along with the differentially expressed genes induced that were related to cell wall remodeling (cellulose, expansin, xyloglucans) and cell expansion (auxin response factors and aquaporin). The results explain the response of leaf thickness to drought stress and show alterations in LT and leaf stomatal traits. This study might serve as a valuable source of gene information for functional studies and provide a theoretical basis to understand leaf growth in terms of leaf anatomy and leaf stomatal traits under drought stress.

Genome-wide and molecular characterization of the DNA replication helicase 2 (DNA2) gene family in rice under drought and salt stress

Rice plants experience various biotic (such as insect and pest attack) and abiotic (such as drought, salt, heat, and cold etc.) stresses during the growing season, resulting in DNA damage and the subsequent losses in rice production. DNA Replication Helicase/Nuclease2 (DNA2) is known to be involved in DNA replication and repair. In animals and yeast DNA2 are well characterized because it has the abilities of both helicase and nuclease, it plays a crucial role in DNA replication in the nucleus and mitochondrial genomes. However; they are not fully examined in plants due to less focused on plants damage repair. To fill this research gap, the current study focused on the genome-wide identification and characterization of OsDNA2 genes, along with analyses of their transcriptional expression, duplication, and phylogeny in rice. Overall, 17 OsDNA2 members were reported to be found on eight different chromosomes (2, 3, 4, 6, 7, 9, 10, and 11). Among these chromosomes (Chr), Chr4 contained a maximum of six OsDNA2 genes. Based on phylogenetic analysis, the OsDNA2 gene members were clustered into three different groups. Furthermore, the conserved domains, gene structures, and cis-regulatory elements were systematically investigated. Gene duplication analysis revealed that OsDNA2_2 had an evolutionary relationship with OsDNA2_14, OsDNA2_5 with OsDNA2_6, and OsDNA2_1 with OsDNA2_8. Moreover, results showed that the conserved domain (AAA_11 superfamily) were present in the OsDNA2 genes, which belongs to the DEAD-like helicase superfamily. In addition, to understand the post-transcriptional modification of OsDNA2 genes, miRNAs were predicted, where 653 miRNAs were reported to target 17 OsDNA2 genes. The results indicated that at the maximum, OsDNA2_1 and OsDNA2_4 were targeted by 74 miRNAs each, and OsDNA2_9 was less targeted (20 miRNAs). The three-dimensional (3D) structures of 17 OsDNA2 proteins were also predicted. Expression of OsDNA2 members was also carried out under drought and salt stresses, and conclusively their induction indicated the possible involvement of OsDNA2 in DNA repair under stress when compared with the control. Further studies are recommended to confirm where this study will offer valuable basic data on the functioning of DNA2 genes in rice and other crop plants.

Brassinosteroids induced drought resistance of contrasting drought-responsive genotypes of maize at physiological and transcriptomic levels

The present study investigated the brassinosteroid-induced drought resistance of contrasting drought-responsive maize genotypes at physiological and transcriptomic levels. The brassinosteroid (BR) contents along with different morphology characteristics, viz., plant height (PH), shoot dry weight (SDW), root dry weight (RDW), number of leaves (NL), the specific mass of the fourth leaf, and antioxidant activities, were investigated in two maize lines that differed in their degree of drought tolerance. In response to either control, drought, or brassinosteroid treatments, the KEGG enrichment analysis showed that plant hormonal signal transduction and starch and sucrose metabolism were augmented in both lines. In contrast, the phenylpropanoid biosynthesis was augmented in lines H21L0R1 and 478. Our results demonstrate drought-responsive molecular mechanisms and provide valuable information regarding candidate gene resources for drought improvement in maize crop. The differences observed for BR content among the maize lines were correlated with their degree of drought tolerance, as the highly tolerant genotype showed higher BR content under drought stress.

Comparative physiological, transcriptomic, and WGCNA analyses reveal the key genes and regulatory pathways associated with drought tolerance in Tartary buckwheat

Drought stress is one of the major abiotic stress factors that affect plant growth and crop productivity. Tartary buckwheat is a nutritionally balanced and flavonoid-rich pseudocereal crop and also has strong adaptability to different adverse environments including drought. However, little is known about its drought tolerance mechanism. In this study, we performed comparative physiological and transcriptomic analyses of two contrasting drought-resistant Tartary buckwheat genotypes under nature drought treatment in the reproductive stage. Under drought stress, the drought-tolerant genotype XZSN had significantly higher contents of relative water, proline, and soluble sugar, as well as lower relative electrolyte leakage in the leaves than the drought-susceptible LK3. A total of 5,058 (2,165 upregulated and 2,893 downregulated) and 5,182 (2,358 upregulated and 2,824 downregulated) potential drought-responsive genes were identified in XZSN and LK3 by transcriptome sequencing analysis, respectively. Among the potential drought-responsive genes of XZSN, 1,206 and 1,274 genes were identified to be potential positive and negative contributors for XZSN having higher drought resistance ability than LK3. Furthermore, 851 out of 1,206 positive drought-resistant genes were further identified to be the core drought-resistant genes of XZSN based on WGCNA analysis, and most of them were induced earlier and quicker by drought stress than those in LK3. Functional annotation of the 851 core drought-resistant genes found that a large number of stress-responsive genes were involved in TFs, abscisic acid (ABA) biosynthesis, signal transduction and response, non-ABA signal molecule biosynthesis, water holding, oxygen species scavenging, osmotic adjustment, cell damage prevention, and so on. Transcriptional regulatory network analyses identified the potential regulators of these drought-resistant functional genes and found that the HD-ZIP and MYB TFs might be the key downstream TFs of drought resistance in Tartary buckwheat. Taken together, these results indicated that the XZSN genotype was more drought-tolerant than the LK3 genotype as evidenced by triggering the rapid and dramatic transcriptional reprogramming of drought-resistant genes to reduce water loss, prevent cell damage, and so on. This research expands our current understanding of the drought tolerance mechanisms of Tartary buckwheat and provides important information for its further drought resistance research and variety breeding.

Leaf transcriptome profiling of contrasting sugarcane genotypes for drought tolerance under field conditions

Drought is the most detrimental abiotic stress to sugarcane production. Nevertheless, transcriptomic analyses remain scarce for field-grown plants. Here we performed comparative transcriptional profiling of two contrasting sugarcane genotypes, ‘IACSP97-7065’ (drought-sensitive) and ‘IACSP94-2094’ (drought-tolerant) grown in a drought-prone environment. Physiological parameters and expression profiles were analyzed at 42 (May) and 117 (August) days after the last rainfall. The first sampling was done under mild drought (soil water potential of −60 kPa), while the second one was under severe drought (soil water potential of −75 kPa). Microarray analysis revealed a total of 622 differentially expressed genes in both sugarcane genotypes under mild and severe drought stress, uncovering about 250 exclusive transcripts to ‘IACSP94-2094’ involved in oxidoreductase activity, transcriptional regulation, metabolism of amino acids, and translation. Interestingly, the enhanced antioxidant system of ‘IACSP94-2094’ may protect photosystem II from oxidative damage, which partially ensures stable photochemical activity even after 117 days of water shortage. Moreover, the tolerant genotype shows a more extensive set of responsive transcription factors, promoting the fine-tuning of drought-related molecular pathways. These results help elucidate the intrinsic molecular mechanisms of a drought-tolerant sugarcane genotype to cope with ever-changing environments, including prolonged water deficit, and may be useful for plant breeding programs.

Physiological and Multi-Omics Approaches for Explaining Drought Stress Tolerance and Supporting Sustainable Production of Rice

Drought differs from other natural disasters in several respects, largely because of the complexity of a crop's response to it and also because we have the least understanding of a crop's inductive mechanism for addressing drought tolerance among all abiotic stressors. Overall, the growth and productivity of crops at a global level is now thought to be an issue that is more severe and arises more frequently due to climatic change-induced drought stress. Among the major crops, rice is a frontline staple cereal crop of the developing world and is critical to sustaining populations on a daily basis. Worldwide, studies have reported a reduction in rice productivity over the years as a consequence of drought. Plants are evolutionarily primed to withstand a substantial number of environmental cues by undergoing a wide range of changes at the molecular level, involving gene, protein and metabolite interactions to protect the growing plant. Currently, an in-depth, precise and systemic understanding of fundamental biological and cellular mechanisms activated by crop plants during stress is accomplished by an umbrella of -omics technologies, such as transcriptomics, metabolomics and proteomics. This combination of multi-omics approaches provides a comprehensive understanding of cellular dynamics during drought or other stress conditions in comparison to a single -omics approach. Thus a greater need to utilize information (big-omics data) from various molecular pathways to develop drought-resilient crop varieties for cultivation in ever-changing climatic conditions. This review article is focused on assembling current peer-reviewed published knowledge on the use of multi-omics approaches toward expediting the development of drought-tolerant rice plants for sustainable rice production and realizing global food security.

Variety-specific transcript accumulation during reproductive stage in drought-stressed rice

The divergence of natural stress tolerance mechanisms between species is an intriguing phenomenon. To study it in rice, a comparative transcriptome analysis was carried out in 'heading' stage tissue (flag leaf, panicles and roots) of Nagina 22 (N22; drought-tolerant) and IR64 (drought-sensitive) plants subjected to field drought. Interestingly, N22 showed almost double the number of differentially expressed genes (DEGs) than IR64. Many DEGs colocalized within drought-related QTLs responsible for grain yield and drought tolerance and also associated with drought tolerance and critical drought-related plant traits such as leaf rolling, trehalose content, sucrose and cellulose content. Besides, co-expression analysis of the DEGs revealed several 'hub' genes known to actively regulate drought stress response. Strikingly, 1366 DEGs, including 21 'hub' genes, showed a distinct opposite regulation in the two rice varieties under similar drought conditions. Annotation of these variety-specific DEGs (VS-DEGs) revealed that they are distributed in various biological pathways. Furthermore, 103 VS-DEGs were found to physically interact with over 1300 genes, including 32 that physically interact with other VS-DEGs as well. The promoter region of these genes has sequence variations among the two rice varieties, which might be in part responsible for their unique expression pattern.

Transcriptome profiling of two rice genotypes under mild field drought stress during grain-filling stage

Drought is one of the most critical abiotic stresses that threaten crop production worldwide. This stress affects the rice crop in all stages of rice development; however, the occurrence during reproductive and grain-filling stages has the most impact on grain yield. Although many global transcriptomic studies have been performed during the reproductive stage in rice, very limited information is available for the grain-filling stage. Hence, we intend to investigate how the rice plant responds to drought stress during the grain-filling stage and how the responses change over time under field conditions. Two rice genotypes were selected for RNA-seq analysis: '4610', previously reported as a moderately tolerant breeding line, and Rondo, an elite indica rice cultivar susceptible to drought conditions. Additionally, 10 agronomic traits were evaluated under normal irrigated and drought conditions. Leaf tissues were collected during grain-filling stages at two time points, 14 and 21 days after the drought treatment, from both the drought field and normal irrigated field conditions. Based on agronomic performances, '4610' was less negatively affected than Rondo under mild drought conditions, and expression profiling largely aligned with the phenotypic data. The transcriptomic data indicated that, in general, '4610' had much earlier responses than its counterpart in mitigating the impact of drought stress. Several key genes and gene families related to drought stress or stress-related conditions were found differentially expressed in this study, including transcription factors, drought tolerance genes and reactive oxygen species scavengers. Furthermore, this study identified novel differentially expressed genes (DEGs) without function annotations that may play roles in drought tolerance-related functions. Some of the important DEGs detected in this study can be targeted for future research.