Elucidation of Novel cis-Regulatory Elements and Promoter Structures Involved in Iron Excess Response Mechanisms in Rice Using a Bioinformatics Approach

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.

Comparative analysis of VMT genes/proteins in selected plant species with emphasis on bread wheat (Triticum aestivum L.)

BACKGROUND

In recent years, the study of molecular basis of uptake, transport and utilization of grain Fe/Zn (GFe/GZn) in wheat has been an active area of research. As a result, it has been shown that a number of transporters are involved in uptake and transport of Fe. In a recent study, knockout of a transporter gene OsVMT (VACUOLAR MUGINEIC ACID TRANSPORTER) in rice was shown to be involved in Fe homoeostasis.

OBJECTIVE

In this study, we analysed VMT genes among six monocots and three dicots with major emphasis on wheat VMT genes (TaVMTs), taking OsVMT gene as a reference.

METHODS AND RESULTS

Using OsVMT gene as a reference, VMT genes were identified and sequence similarities were examined among six monocots and three dicots. Each VMT protein carried one functional domain and 7 to 10 distinct motifs (including 9 novel motifs). The qRT-PCR analysis showed differential expression by all the six TaVMT genes in pairs of contrasting wheat genotypes with high (FAR4 and WB02) and low (K8027 and HD3226) GFe/GZn at two different grain filling stages (14 DAA and 28 DAA). TaVMT1 genes showed up-regulation in high GFe/Zn genotypes relative to low GFe/Zn genotypes, whereas the TaVMT2 genes showed down-regulation or nonsignificant up-regulation in a few cases.

CONCLUSIONS

At 14 DAA, each of the six TaVMT genes exhibited higher expression in wheat genotypes with high GFe and GZn relative to those with low GFe and GZn, suggesting major role of VMT genes in improvement of grain Fe/Zn homoeostasis, thus making TaVMT genes useful for improvement in Fe/Zn in wheat grains.

Identification of oil palm cis-regulatory elements based on DNA free energy and single nucleotide polymorphism density

Transcription control through cis-regulatory elements (CREs) is one of important regulators of gene expression. This study aimed to identify the location of CREs in oil palm (Elaeis guineensis Jacq.) using the combination of DNA free energy and single nucleotide polymorphism (SNP) density approaches. Promoter region sequences were extracted oil palm genome spanning from 1500 nucleotides (nt) upstream to 1000 nt downstream of every annotated transcription start sites (TSS). Free energy profiles of each promoter region were calculated using PromPredict software. Raw reads from the deep sequencing of 59 oil palm origins were used to calculate SNP density of each promoter region. The result showed that the average free energy (AFE) on the upstream region of TSS is about 1.5 kcal/mol higher compared to the downstream region. Using DNA free energy method, 16,281 regions of CREs were predicted. Most of predicted CREs was located between 1 and 500 nt upstream of TSS. Anti-correlation pattern between free energy and SNP density was observed on the predicted regions of CREs. This anti-correlated pattern was also observed on an experimentally determined promoter of the oil palm metallothionein gene, EgMSP1. Considering the increasing use of promoter information on plant biotechnology, an easy and accurate promoter prediction using the combination of free energy and SNP density method could be recommended.

Genome-wide association study unravels genomic regions associated with chlorophyll fluorescence parameters in wheat (Triticum aestivum L.) under different sowing conditions

KEY MESSAGE

Genome-wide association study identified 205 significant marker-trait associations for chlorophyll fluorescence parameters in wheat. Candidate gene mining, in silico expression, and promoter analyses revealed the potential candidate genes associated with the studied parameters. The present study investigated the effect of varied sowing conditions (viz., early, timely, and late) on different chlorophyll fluorescence parameters in diverse wheat germplasm set comprising of 198 lines over two cropping seasons (2020-2021 and 2021-2022). Further, a genome-wide association study was conducted to identify potential genomic regions associated with these parameters. The results revealed significant impacts of sowing conditions on all fluorescence parameters, with the maximum and minimum effects on FI (26.64%) and FV/FM (2.12%), respectively. Among the 205 marker-trait associations (MTAs) identified, 11 high-confidence MTAs were chosen, exhibiting substantial effects on multiple fluorescence parameters, and each explaining more than 10% of the phenotypic variation. Through gene mining of genomic regions encompassing high-confidence MTAs, we identified a total of 626 unique gene models. In silico expression analysis revealed 42 genes with an expression value exceeding 2 TPM. Among them, 10 genes were identified as potential candidate genes with functional relevance to enhanced photosynthetic efficiency. These genes mainly encoded for the following important proteins/products-ankyrin repeat protein, 2Fe-2S ferredoxin-type iron-sulfur-binding domain, NADH-ubiquinone reductase complex-1 MLRQ subunit, oxidoreductase FAD/NAD(P)-binding, photosystem-I PsaF, and protein kinases. Promoter analysis revealed the presence of light-responsive (viz., GT1-motif, TCCC-motif, I-box, GT1-motif, TCT-motif, and SP-1) and stress-responsive (viz., ABRE, AuxRR-core, GARE-motif, and ARE) cis-regulatory elements, which may be involved in the regulation of identified putative candidate genes. Findings from this study could directly help wheat breeders in selecting lines with favorable alleles for chlorophyll fluorescence, while the identified markers will facilitate marker-assisted selection of potential genomic regions for improved photosynthesis.

Magnesium supply alleviates iron toxicity-induced leaf bronzing in rice through exclusion and tissue-tolerance mechanisms

Introduction

Iron (Fe) toxicity is a widespread nutritional disorder in lowland rice causing growth retardation and leaf symptoms referred to as leaf bronzing. It is partly caused by an imbalance of nutrients other than Fe and supply of these is known to mitigate the toxicity. But the physiological and molecular mechanisms involved are unknown.

Methods

We investigated the effect of magnesium (Mg) on Fe toxicity tolerance in a field study in the Central Highlands of Madagascar and in hydroponic experiments with excess Fe (300 mg Fe L-1). An RNA-seq analysis was conducted in a hydroponic experiment to elucidate possible mechanisms underlying Mg effects.

Results and discussion

Addition of Mg consistently decreased leaf bronzing under both field and hydroponic conditions, whereas potassium (K) addition caused minor effects. Plants treated with Mg tended to have smaller shoot Fe concentrations in the field, suggesting enhanced exclusion at the whole-plant level. However, analysis of multiple genotypes showed that Fe toxicity symptoms were also mitigated without a concomitant decrease of Fe concentration, suggesting that increased Mg supply confers tolerance at the tissue level. The hydroponic experiments also suggested that Mg mitigated leaf bronzing without significantly decreasing Fe concentration or oxidative stress as assessed by the content of malondialdehyde, a biomarker for oxidative stress. An RNA-seq analysis revealed that Mg induced more changes in leaves than roots. Subsequent cis-element analysis suggested that NAC transcription factor binding sites were enriched in genes induced by Fe toxicity in leaves. Addition of Mg caused non-significant enrichment of the same binding sites, suggesting that NAC family proteins may mediate the effect of Mg. This study provides clues for mitigating Fe toxicity-induced leaf bronzing in rice.

Predicting transcriptional responses to heat and drought stress from genomic features using a machine learning approach in rice

Plants have evolved various mechanisms to adapt to adverse environmental stresses, such as the modulation of gene expression. Expression of stress-responsive genes is controlled by specific regulators, including transcription factors (TFs), that bind to sequence-specific binding sites, representing key components of cis-regulatory elements and regulatory networks. Our understanding of the underlying regulatory code remains, however, incomplete. Recent studies have shown that, by training machine learning (ML) algorithms on genomic sequence features, it is possible to predict which genes will transcriptionally respond to a specific stress. By identifying the most important features for gene expression prediction, these trained ML models allow, in theory, to further elucidate the regulatory code underlying the transcriptional response to abiotic stress. Here, we trained random forest ML models to predict gene expression in rice (Oryza sativa) in response to heat or drought stress. Apart from thoroughly assessing model performance and robustness across various input training data, the importance of promoter and gene body sequence features to train ML models was evaluated. The use of enriched promoter oligomers, complementing known TF binding sites, allowed us to gain novel insights in DNA motifs contributing to the stress regulatory code. By comparing genomic feature importance scores for drought and heat stress over time, general and stress-specific genomic features contributing to the performance of the learned models and their temporal variation were identified. This study provides a solid foundation to build and interpret ML models accurately predicting transcriptional responses and enables novel insights in biological sequence features that are important for abiotic stress responses.

Temporal regulation of cold transcriptional response in switchgrass

Switchgrass low-land ecotypes have significantly higher biomass but lower cold tolerance compared to up-land ecotypes. Understanding the molecular mechanisms underlying cold response, including the ones at transcriptional level, can contribute to improving tolerance of high-yield switchgrass under chilling and freezing environmental conditions. Here, by analyzing an existing switchgrass transcriptome dataset, the temporal cis-regulatory basis of switchgrass transcriptional response to cold is dissected computationally. We found that the number of cold-responsive genes and enriched Gene Ontology terms increased as duration of cold treatment increased from 30 min to 24 hours, suggesting an amplified response/cascading effect in cold-responsive gene expression. To identify genomic sequences likely important for regulating cold response, machine learning models predictive of cold response were established using k-mer sequences enriched in the genic and flanking regions of cold-responsive genes but not non-responsive genes. These k-mers, referred to as putative cis-regulatory elements (pCREs) are likely regulatory sequences of cold response in switchgrass. There are in total 655 pCREs where 54 are important in all cold treatment time points. Consistent with this, eight of 35 known cold-responsive CREs were similar to top-ranked pCREs in the models and only these eight were important for predicting temporal cold response. More importantly, most of the top-ranked pCREs were novel sequences in cold regulation. Our findings suggest additional sequence elements important for cold-responsive regulation previously not known that warrant further studies.

Consensus genomic regions associated with multiple abiotic stress tolerance in wheat and implications for wheat breeding

In wheat, a meta-analysis was performed using previously identified QTLs associated with drought stress (DS), heat stress (HS), salinity stress (SS), water-logging stress (WS), pre-harvest sprouting (PHS), and aluminium stress (AS) which predicted a total of 134 meta-QTLs (MQTLs) that involved at least 28 consistent and stable MQTLs conferring tolerance to five or all six abiotic stresses under study. Seventy-six MQTLs out of the 132 physically anchored MQTLs were also verified with genome-wide association studies. Around 43% of MQTLs had genetic and physical confidence intervals of less than 1 cM and 5 Mb, respectively. Consequently, 539 genes were identified in some selected MQTLs providing tolerance to 5 or all 6 abiotic stresses. Comparative analysis of genes underlying MQTLs with four RNA-seq based transcriptomic datasets unravelled a total of 189 differentially expressed genes which also included at least 11 most promising candidate genes common among different datasets. The promoter analysis showed that the promoters of these genes include many stress responsiveness cis-regulatory elements, such as ARE, MBS, TC-rich repeats, As-1 element, STRE, LTR, WRE3, and WUN-motif among others. Further, some MQTLs also overlapped with as many as 34 known abiotic stress tolerance genes. In addition, numerous ortho-MQTLs among the wheat, maize, and rice genomes were discovered. These findings could help with fine mapping and gene cloning, as well as marker-assisted breeding for multiple abiotic stress tolerances in wheat.

Consensus genomic regions associated with multiple abiotic stress tolerance in wheat and implications for wheat breeding

In wheat, a meta-analysis was performed using previously identified QTLs associated with drought stress (DS), heat stress (HS), salinity stress (SS), water-logging stress (WS), pre-harvest sprouting (PHS), and aluminium stress (AS) which predicted a total of 134 meta-QTLs (MQTLs) that involved at least 28 consistent and stable MQTLs conferring tolerance to five or all six abiotic stresses under study. Seventy-six MQTLs out of the 132 physically anchored MQTLs were also verified with genome-wide association studies. Around 43% of MQTLs had genetic and physical confidence intervals of less than 1 cM and 5 Mb, respectively. Consequently, 539 genes were identified in some selected MQTLs providing tolerance to 5 or all 6 abiotic stresses. Comparative analysis of genes underlying MQTLs with four RNA-seq based transcriptomic datasets unravelled a total of 189 differentially expressed genes which also included at least 11 most promising candidate genes common among different datasets. The promoter analysis showed that the promoters of these genes include many stress responsiveness cis-regulatory elements, such as ARE, MBS, TC-rich repeats, As-1 element, STRE, LTR, WRE3, and WUN-motif among others. Further, some MQTLs also overlapped with as many as 34 known abiotic stress tolerance genes. In addition, numerous ortho-MQTLs among the wheat, maize, and rice genomes were discovered. These findings could help with fine mapping and gene cloning, as well as marker-assisted breeding for multiple abiotic stress tolerances in wheat.

Heme oxygenase-nitric oxide crosstalk-mediated iron homeostasis in plants under oxidative stress

Plant growth under abiotic stress conditions significantly enhances intracellular generation of reactive oxygen species (ROS). Oxidative status of plant cells is directly affected by the modulation of iron homeostasis. Among mammals and plants, heme oxygenase-1 (HO-1) is a well-known antioxidant enzyme. It catalyzes oxygenation of heme, thereby producing Fe²⁺, CO and biliverdin as byproducts. The antioxidant potential of HO-1 is primarily due to its catalytic reaction byproducts. Biliverdin and bilirubin possess conjugated π-electrons which escalate the ability of these biomolecules to scavenge free radicals. CO also enhances the ROS scavenging ability of plants cells by upregulating catalase and peroxidase activity. Enhanced expression of HO-1 in plants under oxidative stress accompanies sequestration of iron in specialized iron storage proteins localized in plastids and mitochondria, namely ferritin for Fe³⁺ storage and frataxin for storage of Fe-S clusters, respectively. Nitric oxide (NO) crosstalks with HO-1 at multiple levels, more so in plants under oxidative stress, in order to maintain intracellular iron status. Formation of dinitrosyl-iron complexes (DNICs) significantly prevents Fenton reaction during oxidative stress. DNICs also release NO upon dissociation in target cells over long distance in plants. They also function as antioxidants against superoxide anions and lipidic free radicals. A number of NO-modulated transcription factors also facilitate iron homeostasis in plant cells. Plants facing oxidative stress exhibit modulation of lateral root formation by HO-1 through NO and auxin-dependent pathways. The present review provides an in-depth analysis of the structure-function relationship of HO-1 in plants and mammals, correlating them with their adaptive mechanisms of survival under stress.

Transcriptomic Analysis Revealed Reactive Oxygen Species Scavenging Mechanisms Associated With Ferrous Iron Toxicity in Aromatic Keteki Joha Rice

Lowland acidic soils with water-logged regions are often affected by ferrous iron (Fe²⁺) toxicity, a major yield-limiting factor of rice production. Under severe Fe²⁺ toxicity, reactive oxygen species (ROS) are crucial, although molecular mechanisms and associated ROS homeostasis genes are still unknown. In this study, a comparative RNA-Seq based transcriptome analysis was conducted to understand the Fe²⁺ toxicity tolerance mechanism in aromatic Keteki Joha. About 69 Fe homeostasis related genes and their homologs were identified, where most of the genes were downregulated. Under severe Fe²⁺ toxicity, the biosynthesis of amino acids, RNA degradation, and glutathione metabolism were induced, whereas phenylpropanoid biosynthesis, photosynthesis, and fatty acid elongation were inhibited. The mitochondrial iron transporter (OsMIT), vacuolar iron transporter 2 (OsVIT2), ferritin (OsFER), vacuolar mugineic acid transporter (OsVMT), phenolic efflux zero1 (OsPEZ1), root meander curling (OsRMC), and nicotianamine synthase (OsNAS3) were upregulated in different tissues, suggesting the importance of Fe retention and sequestration for detoxification. However, several antioxidants, ROS scavenging genes and abiotic stress-responsive transcription factors indicate ROS homeostasis as one of the most important defense mechanisms under severe Fe²⁺ toxicity. Catalase (CAT), glutathione (GSH), ascorbate peroxidase (APX), monodehydroascorbate reductase (MDHAR), dehydroascorbate reductase (DHAR), and glutathione reductase (GR) were upregulated. Moreover, abiotic stress-responsive transcription factors, no apical meristem (NAC), myeloblastosis (MYB), auxin response factor (ARF), basic helix-loop-helix (bZIP), WRKY, and C2H2-zinc finger protein (C2H2-ZFP) were also upregulated. Accordingly, ROS homeostasis has been proposed as an essential defense mechanism under such conditions. Thus, the current study may enrich the understanding of Fe-homeostasis in rice.

Design and Analysis of Native Photorespiration Gene Motifs of Promoter Untranslated Region Combinations Under Short Term Abiotic Stress Conditions

Quantitative traits are rarely controlled by a single gene, thereby making multi-gene transformation an indispensable component of modern synthetic biology approaches. However, the shortage of unique gene regulatory elements (GREs) for the robust simultaneous expression of multiple nuclear transgenes is a major bottleneck that impedes the engineering of complex pathways in plants. In this study, we compared the transcriptional efficacies of a comprehensive list of well-documented promoter and untranslated region (UTR) sequences side by side. The strength of GREs was examined by a dual-luciferase assay in conjunction with transient expression in tobacco. In addition, we created suites of new GREs with higher transcriptional efficacies by combining the best performing promoter-UTR sequences. We also tested the impact of elevated temperature and high irradiance on the effectiveness of these GREs. While constitutive promoters ensure robust expression of transgenes, they lack spatiotemporal regulations exhibited by native promoters. Here, we present a proof-of-principle study on the characterization of synthetic promoters based on cis-regulatory elements of three key photorespiratory genes. This conserved biochemical process normally increases under elevated temperature, low CO2, and high irradiance stress conditions and results in ∼25% loss in fixed CO2. To select stress-responsive cis-regulatory elements involved in photorespiration, we analyzed promoters of two chloroplast transporters (AtPLGG1 and AtBASS6) and a key plastidial enzyme, AtPGLP using PlantPAN3.0 and AthaMap. Our results suggest that these motifs play a critical role for PLGG1, BASS6, and PGLP in mediating response to elevated temperature and high-intensity light stress. These findings will not only enable the advancement of metabolic and genetic engineering of photorespiration but will also be instrumental in related synthetic biology approaches.