Patterns of cis-element enrichment reveal potential regulatory modules involved in the transcriptional regulation of anoxia response of japonica rice

Leaf rolling in bread wheat (Triticum aestivum L.) is controlled by the upregulation of a pair of closely linked/duplicate zinc finger homeodomain class transcription factors during moisture stress conditions

Zinc finger-homeodomain (ZF-HDs) class IV transcriptional factors (TFs) is a plant-specific transcription factor and play a key role in stress responses, plant growth, development, and hormonal signaling. In this study, two new leaf rolling TFs genes, namely TaZHD1 and TaZHD10, were identified in wheat using comparative genomic analysis of the target region that carried a major QTL for leaf rolling identified through multi-environment phenotyping and high throughput genotyping of a RIL population. Structural and functional annotation of the candidate ZHD genes with its closest rice orthologs reflects the species-specific evolution and, undoubtedly, validates the notions of remote-distance homology concept. Meanwhile, the morphological analysis resulted in contrasting difference for leaf rolling in extreme RILs between parental lines HD2012 and NI5439 at booting and heading stages. Transcriptome-wide expression profiling revealed that TaZHD10 transcripts showed significantly higher expression levels than TaZHD1 in all leaf tissues upon drought stress. The relative expression of these genes was further validated by qRT-PCR analysis, which also showed consistent results across the studied genotypes at the booting and anthesis stage. The contrasting modulation of these genes under drought conditions and the available evidenced for its epigenetic behavior that might involve the regulation of metabolic and gene regulatory networks. Prediction of miRNAs resulted in five Tae-miRs that could be associated with RNAi mediated control of TaZHD1 and TaZHD10 putatively involved in the metabolic pathway controlling rolled leaf phenotype. Gene interaction network analysis indicated that TaZHD1 and TaZHD10 showed pleiotropic effects and might also involve other functions in wheat in addition to leaf rolling. Overall, the results increase our understanding of TaZHD genes and provide valuable information as robust candidate genes for future functional genomics research aiming for the breeding of wheat varieties tolerant to leaf rolling.

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.

Dynamic genome-wide association analysis and identification of candidate genes involved in anaerobic germination tolerance in rice

BACKGROUND

In Asian rice production, an increasing number of countries now choose the direct seeding mode because of rising costs, labour shortages and water shortages. The ability of rice seeds to undergo anaerobic germination (AG) plays an important role in the success of direct seeding.

RESULTS

In this study, we used 2,123,725 single nucleotide polymorphism (SNP) markers based on resequencing to conduct a dynamic genome-wide association study (GWAS) of coleoptile length (CL) and coleoptile diameter (CD) in 209 natural rice populations. A total of 26 SNP loci were detected in these two phenotypes, of which 5 overlapped with previously reported loci (S1_ 39674301, S6_ 20797781, S7_ 18722403, S8_ 9946213, S11_ 19165397), and two sites were detected repeatedly at different time points (S3_ 24689629 and S5_ 27918754). We suggest that these 7 loci (-log₁₀ (P) value > 7.3271) are the key sites that affect AG tolerance. To screen the candidate genes more effectively, we sequenced the transcriptome of the flooding-tolerant variety R151 in six key stages, including anaerobic (AN) and the oxygen conversion point (AN-A), and obtained high-quality differential expression profiles. Four reliable candidate genes were identified: Os01g0911700 (OsVP1), Os05g0560900 (OsGA2ox8), Os05g0562200 (OsDi19-1) and Os06g0548200. Then qRT-PCR and LC-MS/ MS targeting metabolite detection technology were used to further verify that the up-regulated expression of these four candidate genes was closely related to AG.

CONCLUSION

The four novel candidate genes were associated with gibberellin (GA) and abscisic acid (ABA) regulation and cell wall metabolism under oxygen-deficiency conditions and promoted coleoptile elongation while avoiding adverse effects, allowing the coleoptile to obtain oxygen, escape the low-oxygen environment and germinate rapidly. The results of this study improve our understanding of the genetic basis of AG in rice seeds, which is conducive to the selection of flooding-tolerant varieties suitable for direct seeding.

Molecular Mechanisms Supporting Rice Germination and Coleoptile Elongation under Low Oxygen

Rice germinates under submergence by exploiting the starch available in the endosperm and translocating sugars from source to sink organs. The availability of fermentable sugar under water allows germination with the protrusion of the coleoptile, which elongates rapidly and functions as a snorkel toward the air above. Depending on the variety, rice can produce a short or a long coleoptile. Longer length entails the involvement of a functional transport of auxin along the coleoptile. This paper is an overview of rice coleoptiles and the studies undertaken to understand its functioning and role under submergence.

Differentially Expressed MicroRNAs Link Cellular Physiology to Phenotypic Changes in Rice Under Stress Conditions

Plant microRNAs (miRNAs) and their target genes have important functional roles in nutrition deficiency and stress response. However, the underlying mechanisms relating relative expression of miRNAs and target mRNAs to morphological adjustments are not well defined. By combining miRNA expression profiles, corresponding target genes and transcription factors that bind to computationally identified over-represented cis-regulatory elements (CREs) common in miRNAs and target gene promoters, we implement a strategy that identifies a set of differentially expressed regulatory interactions which, in turn, relate underlying cellular mechanisms to some of the phenotypic changes observed. Integration of experimentally reported individual interactions with identified regulatory interactions explains how (i) during mineral deficiency osa-miR167 inhibits shoot growth but activates adventitious root growth by influencing free auxin content; (ii) during sulfur deficiency osa-miR394 is involved in adventitious root growth inhibition, sulfur and iron homeostasis, and auxin-mediated regulation of sulfur homeostasis; (iii) osa-miR399 contributes to cross-talk between cytokinin and phosphorus deficiency signaling; and (iv) a feed-forward loop involving the osa-miR166, trihelix and HD-ZIP III transcription factors may regulate leaf senescence during drought. This strategy not only identifies various regulatory interactions connecting phenotypic changes with cellular or molecular events triggered by stress, but also provides a framework to deepen our understanding of stress cellular physiology.

MaRAP2-4, a waterlogging-responsive ERF from Mentha, regulates bidirectional sugar transporter AtSWEET10 to modulate stress response in Arabidopsis

As waterlogging and successive events severely influence growth and development of economically important plants, we attempted to characterize the role of a waterlogging-responsive group I (A-6) ethylene response factor (MaRAP2-4) from Mentha arvensis. Waterlogging, ethylene and methyl jasmonate rapidly induced the expression of MaRAP2-4. MaRAP2-4 interacted with multiple cis-elements like dehydration response elements (DRE1/2), anoxia/jasmonic acid response element (JARE) and GCC box showing its involvement in multiple responses. MaRAP2-4 localizes in the nucleus and acts as a transcriptional activator. Truncation and internal deletion identified a 20 amino acids potential transactivation domain (PLPSSVDAKLEAICQSLAIN) in MaRAP2-4. MaRAP2-4 transgenic Arabidopsis showed enhanced waterlogging and subsequent oxidative stress tolerance. Microarray analysis revealed that within up-regulated genes 483, 212 and 132 promoters carry either single or multiple copies of DRE, JARE and GCC cis-element/s, respectively. Within these promoters, a large section belongs to carbohydrate metabolism/transport, including many SWEET transporters. Further analysis showed MaRAP2-4 specifically targets two positions in AtSWEEET10 promoter carrying DRE and/or GCC box that might regulate carbohydrate availability and waterlogging tolerance. These results demonstrate that MaRAP2-4 is a positive regulator of waterlogging tolerance, and as energy-consuming processes such as carbohydrate biosynthesis are reduced under waterlogging-induced hypoxia, sugar transport through SWEETs may be the primary option to make sugar available to the required tissue.

Rice Seed Germination Underwater: Morpho-Physiological Responses and the Bases of Differential Expression of Alcoholic Fermentation Enzymes

The water-, energy-, and labor-intensive system of transplanted puddled rice (Oryza sativa) is steadily being replaced by direct seeding due to the progressive scarcity of these resources. However, the alternate dry direct seeding leads to competition with weeds and poor establishment when soils are flooded. Direct seeded rice capable of anaerobic germination (germination in flooded soil, AG) is ideal, which under rainfed ecosystems would also overcome waterlogging during germination. AG tolerance is associated with faster germination and faster elongation of coleoptiles, with the activities of alcoholic fermentation enzymes replacing aerobic respiration as a source of energy. To better understand the variability in the morpho-physiological responses and in the nature of the alcoholic fermentation enzymes during AG, 21 rice genotypes were studied. The genotypes Khao Hlan On (KHO) and IR42 were used as the tolerant and susceptible checks, respectively. KHO exhibited faster germination, with 82.5% of the coleoptiles emerging out of 10 cm of water within 8 days, whereas IR42 exhibited 20% germination and limited coleoptile growth. Among the test genotypes, four performed well, including two that are drought tolerant. Increased content and activity of the alcoholic fermentation enzymes, alcohol dehydrogenase (ADH1) and acetaldehyde dehydrogenase (ALDH2a and ALDH2b), was noted in KHO under anaerobic than under aerobic conditions and also in comparison with IR42 under AG. Gene transcripts for these enzymes were also more in KHO undergoing AG. However, no major differences were observed between KHO and IR42 in the critical cis-acting regulatory elements, such as the auxin, light, and sugar response elements, in the promoters of ADH1, ALDH2a, and ALDH2b genes. Post-transcriptional and post-translational regulatory mechanisms were implicated for the increased transcript and protein content/activity of the enzymes in KHO by observing four different transcripts of ALDH2a and a unique non-glycosylated form of ADH1 under AG. IR42 lacked the non-glycosylated ADH1 and contained only a truncated form of ALDH2a, which lacked the active site. Additionally, KHO exhibited increased activity and more isoforms for reactive oxygen species detoxifying enzymes under AG compared to IR42. These results highlight the need for a deeper functional understanding of the critical enzymes involved in AG.

Common Stress Transcriptome Analysis Reveals Functional and Genomic Architecture Differences Between Early and Delayed Response Genes

To identify the similarities among responses to diverse environmental stresses, we analyzed the transcriptome response of rice roots to three rhizotoxic perturbations (chromium, ferulic acid and mercury) and identified common early-transient, early-constant and delayed gene inductions. Common early response genes were mostly associated with signal transduction and hormones, and delayed response genes with lipid metabolism. Network component analysis revealed complicated interactions among common genes, the most highly connected signaling hubs being PP2C68, MPK5, LRR-RLK and NPR1. Gene architecture studies revealed different conserved promoter motifs and a different ratio of CpG island distribution between early and delayed genes. In addition, early-transient genes had more exons and a shorter first exon. IMEter was used to calculate the transcription regulation effects of introns, with greater effects for the first introns of early-transient than delayed genes. The higher Ka/Ks (non-synonymous/synonymous mutation) ratio of early-constant genes than early-transient, delayed and the genome median demonstrates the rapid evolution of early-constant genes. Our results suggest that finely tuned transcriptional control in response to environmental stress in rice depends on genomic architecture and signal intensity and duration.

Identification of candidate network hubs involved in metabolic adjustments of rice under drought stress by integrating transcriptome data and genome-scale metabolic network

In this study, we have integrated a rice genome-scale metabolic network and the transcriptome of a drought-tolerant rice line, DK151, to identify the major transcriptional regulators involved in metabolic adjustments necessary for adaptation to drought. This was achieved by examining the differential expressions of transcription factors and metabolic genes in leaf, root and young panicle of rice plants subjected to drought stress during tillering, booting and panicle elongation stages. Critical transcription factors such as AP2/ERF, bZIP, MYB and NAC that control the important nodes in the gene regulatory pathway were identified through correlative analysis of the patterns of spatio-temporal expression and cis-element enrichment. We showed that many of the candidate transcription factors involved in metabolic adjustments were previously linked to phenotypic variation for drought tolerance. This approach represents the first attempt to integrate models of transcriptional regulation and metabolic pathways for the identification of candidate regulatory genes for targeted selection in rice breeding.

Deciphering Cis-Regulatory Element Mediated Combinatorial Regulation in Rice under Blast Infected Condition

Combinations of cis-regulatory elements (CREs) present at the promoters facilitate the binding of several transcription factors (TFs), thereby altering the consequent gene expressions. Due to the eminent complexity of the regulatory mechanism, the combinatorics of CRE-mediated transcriptional regulation has been elusive. In this work, we have developed a new methodology that quantifies the co-occurrence tendencies of CREs present in a set of promoter sequences; these co-occurrence scores are filtered in three consecutive steps to test their statistical significance; and the significantly co-occurring CRE pairs are presented as networks. These networks of co-occurring CREs are further transformed to derive higher order of regulatory combinatorics. We have further applied this methodology on the differentially up-regulated gene-sets of rice tissues under fungal (Magnaporthe) infected conditions to demonstrate how it helps to understand the CRE-mediated combinatorial gene regulation. Our analysis includes a wide spectrum of biologically important results. The CRE pairs having a strong tendency to co-occur often exhibit very similar joint distribution patterns at the promoters of rice. We couple the network approach with experimental results of plant gene regulation and defense mechanisms and find evidences of auto and cross regulation among TF families, cross-talk among multiple hormone signaling pathways, similarities and dissimilarities in regulatory combinatorics between different tissues, etc. Our analyses have pointed a highly distributed nature of the combinatorial gene regulation facilitating an efficient alteration in response to fungal attack. All together, our proposed methodology could be an important approach in understanding the combinatorial gene regulation. It can be further applied to unravel the tissue and/or condition specific combinatorial gene regulation in other eukaryotic systems with the availability of annotated genomic sequences and suitable experimental data.

Upstream regulatory architecture of rice genes: summarizing the baseline towards genus-wide comparative analysis of regulatory networks and allele mining

Dissecting the upstream regulatory architecture of rice genes and their cognate regulator proteins is at the core of network biology and its applications to comparative functional genomics. With the rapidly advancing comparative genomics resources in the genus Oryza, a reference genome annotation that defines the various cis-elements and trans-acting factors that interface each gene locus with various intrinsic and extrinsic signals for growth, development, reproduction and adaptation must be established to facilitate the understanding of phenotypic variation in the context of regulatory networks. Such information is also important to establish the foundation for mining non-coding sequence variation that defines novel alleles and epialleles across the enormous phenotypic diversity represented in rice germplasm. This review presents a synthesis of the state of knowledge and consensus trends regarding the various cis-acting and trans-acting components that define spatio-temporal regulation of rice genes based on representative examples from both foundational studies in other model and non-model plants, and more recent studies in rice. The goal is to summarize the baseline for systematic upstream sequence annotation of the rapidly advancing genome sequence resources in Oryza in preparation for genus-wide functional genomics. Perspectives on the potential applications of such information for gene discovery, network engineering and genomics-enabled rice breeding are also discussed.

Microarray analysis of laser-microdissected tissues indicates the biosynthesis of suberin in the outer part of roots during formation of a barrier to radial oxygen loss in rice (Oryza sativa)

Internal aeration is crucial for root growth in waterlogged soil. A barrier to radial oxygen loss (ROL) can enhance long-distance oxygen transport via the aerenchyma to the root tip; a higher oxygen concentration at the apex enables root growth into anoxic soil. The ROL barrier is formed within the outer part of roots (OPR). Suberin and/or lignin deposited in cell walls are thought to contribute to the barrier, but it is unclear which compound is the main constituent. This study describes gene expression profiles during ROL barrier formation in rice roots to determine the relative responses of suberin and/or lignin biosyntheses for the barrier. OPR tissues were isolated by laser microdissection and their transcripts were analysed by microarray. A total of 128 genes were significantly up- or downregulated in the OPR during the barrier formation. Genes associated with suberin biosynthesis were strongly upregulated, whereas genes associated with lignin biosynthesis were not. By an ab initio analysis of the promoters of the upregulated genes, the putative cis-elements that could be associated with transcription factors, WRKY, AP2/ERF, NAC, bZIP, MYB, CBT/DREB, and MADS, were elucidated. They were particularly associated with the expression of transcription factor genes containing WRKY, AP2, and MYB domains. A semiquantitative reverse-transcription PCR analysis of genes associated with suberin biosynthesis (WRKY, CYP, and GPAT) confirmed that they were highly expressed during ROL barrier formation. Overall, these results suggest that suberin is a major constituent of the ROL barrier in roots of rice.

Metabolic and transcriptional regulatory mechanisms underlying the anoxic adaptation of rice coleoptile

The ability of rice to germinate under anoxia by extending the coleoptile is a highly unusual characteristic and a key feature underpinning the ability of rice seeds to establish in such a stressful environment. The process has been a focal point for research for many years. However, the molecular mechanisms underlying the anoxic growth of the coleoptile still remain largely unknown. To unravel the key regulatory mechanisms of rice germination under anoxic stress, we combined in silico modelling with gene expression data analysis. Our initial modelling analysis via random flux sampling revealed numerous changes in rice primary metabolism in the absence of oxygen. In particular, several reactions associated with sucrose metabolism and fermentation showed a significant increase in flux levels, whereas reaction fluxes across oxidative phosphorylation, the tricarboxylic acid cycle and the pentose phosphate pathway were down-regulated. The subsequent comparative analysis of the differences in calculated fluxes with previously published gene expression data under air and anoxia identified at least 37 reactions from rice central metabolism that are transcriptionally regulated. Additionally, cis-regulatory content analyses of these transcriptionally controlled enzymes indicate a regulatory role for transcription factors such as MYB, bZIP, ERF and ZnF in transcriptional control of genes that are up-regulated during rice germination and coleoptile elongation under anoxia.

Differential anoxic expression of sugar-regulated genes reveals diverse interactions between sugar and anaerobic signaling systems in rice

The interaction between the dual roles of sugar as a metabolic fuel and a regulatory molecule was unveiled by examining the changes in sugar signaling upon oxygen deprivation, which causes the drastic alteration in the cellular energy status. In our study, the expression of anaerobically induced genes is commonly responsive to sugar, either under the control of hexokinase or non-hexokinase mediated signaling cascades. Only sugar regulation via the hexokinase pathway was susceptible for O2 deficiency or energy deficit conditions evoked by uncoupler. Examination of sugar regulation of those genes under anaerobic conditions revealed the presence of multiple paths underlying anaerobic induction of gene expression in rice, subgrouped into three distinct types. The first of these, which was found in type-1 genes, involved neither sugar regulation nor additional anaerobic induction under anoxia, indicating that anoxic induction is a simple result from the release of sugar repression by O2-deficient conditions. In contrast, type-2 genes also showed no sugar regulation, albeit with enhanced expression under anoxia. Lastly, expression of type-3 genes is highly enhanced with sugar regulation sustained under anoxia. Intriguingly, the inhibition of the mitochondrial ATP synthesis can reproduce expression pattern of a specific set of anaerobically induced genes, implying that rice cells may sense O2 deprivation, partly via perception of the perturbed cellular energy status. Our study of interaction between sugar signaling and anaerobic conditions has revealed that sugar signaling and the cellular energy status are likely to communicate with each other and influence anaerobic induction of gene expression in rice.

Tolerance of anaerobic conditions caused by flooding during germination and early growth in rice (Oryza sativa L.)

Rice is semi-aquatic, adapted to a wide range of hydrologies, from aerobic soils in uplands to anaerobic and flooded fields in waterlogged lowlands, to even deeply submerged soils in flood-prone areas. Considerable diversity is present in native rice landraces selected by farmers over centuries. Our understanding of the adaptive features of these landraces to native ecosystems has improved considerably over the recent past. In some cases, major genes associated with tolerance have been cloned, such as SUB1A that confers tolerance of complete submergence and SNORKEL genes that control plant elongation to escape deepwater. Modern rice varieties are sensitive to flooding during germination and early growth, a problem commonly encountered in rainfed areas, but few landraces capable of germination under these conditions have recently been identified, enabling research into tolerance mechanisms. Major QTLs were also identified, and are being targeted for molecular breeding and for cloning. Nevertheless, limited progress has been made in identifying regulatory processes for traits that are unique to tolerant genotypes, including faster germination and coleoptile elongation, formation of roots and leaves under hypoxia, ability to catabolize starch into simple sugars for subsequent use in glycolysis and fermentative pathways to generate energy. Here we discuss the state of knowledge on the role of the PDC-ALDH-ACS bypass and the ALDH enzyme as the likely candidates effective in tolerant rice genotypes. Potential involvement of factors such as cytoplasmic pH regulation, phytohormones, reactive oxygen species scavenging and other metabolites is also discussed. Further characterization of contrasting genotypes would help in elucidating the genetic and biochemical regulatory and signaling mechanisms associated with tolerance. This could facilitate breeding rice varieties suitable for direct seeding systems and guide efforts for improving waterlogging tolerance in other crops.