Abstract
Nitrogen (N) and Water (W) - two resources critical for crop productivity – are becoming increasingly limited in soils globally. To address this issue, we aim to uncover the gene regulatory networks (GRNs) that regulate nitrogen use efficiency (NUE) - as a function of water availability - in Oryza sativa, a staple for 3.5 billion people. In this study, we infer and validate GRNs that correlate with rice NUE phenotypes affected by N-by-W availability in the field. We did this by exploiting RNA-seq and crop phenotype data from 19 rice varieties grown in a 2x2 N-by-W matrix in the field. First, to identify gene-to-NUE field phenotypes, we analyzed these datasets using weighted gene co-expression network analysis (WGCNA). This identified two network modules ("skyblue" & "grey60") highly correlated with NUE grain yield (NUEg). Next, we focused on 90 TFs contained in these two NUEg modules and predicted their genome-wide targets using the N-and/or-W response datasets using a random forest network inference approach (GENIE3). Next, to validate the GENIE3 TF→target gene predictions, we performed Precision/Recall Analysis (AUPR) using nine datasets for three TFs validated in planta. This analysis sets a precision threshold of 0.31, used to "prune" the GENIE3 network for high-confidence TF→target gene edges, comprising 88 TFs and 5,716 N-and/or-W response genes. Next, we ranked these 88 TFs based on their significant influence on NUEg target genes responsive to N and/or W signaling. This resulted in a list of 18 prioritized TFs that regulate 551 NUEg target genes responsive to N and/or W signals. We validated the direct regulated targets of two of these candidate NUEg TFs in a plant cell-based TF assay called TARGET, for which we also had in planta data for comparison. Gene ontology analysis revealed that 6/18 NUEg TFs - OsbZIP23 (LOC_Os02g52780), Oshox22 (LOC_Os04g45810), LOB39 (LOC_Os03g41330), Oshox13 (LOC_Os03g08960), LOC_Os11g38870, and LOC_Os06g14670 - regulate genes annotated for N and/or W signaling. Our results show that OsbZIP23 and Oshox22, known regulators of drought tolerance, also coordinate W-responses with NUEg. This validated network can aid in developing/breeding rice with improved yield on marginal, low N-input, drought-prone soils.
Introduction
Nitrogen (N) and water (W) are essential resources for plant productivity that are becoming increasingly limited in marginal soils world-wide (; ). Moreover, applications of N and W in agriculture are costly resources to society (Williamson, 2011; ; ). Most studies in major crops like rice, examine the effects of N and drought separately (; ; Volante et al., 2017; Zhao et al., 2017). More recently, studies that examine how the interaction between N and W availability affects rice phenotypes and gene regulation have been examined (Swift et al., 2019; ; ; Sevanthi et al., 2021).
Several studies have shown that genes critical to N-uptake, sensing and metabolism have been associated with a drought phenotype. For example, NRT1.1/CHL1/NPF6.3 the a dual-affinity nitrate transporter () is expressed in the guard cells in Arabidopsis. Moreover, nrt1.1/chl1 mutant is more drought tolerant compared to wild-type. The loss of NRT1.1/CHL1 reduced the stomatal opening and transpiration rates which contribute to its drought-tolerant phenotype (). Next, mutants in nitrate reductase in both Arabidopsis (NIA1 and NIA2) and rice (OsNR1.2) exhibit a drought-tolerant phenotype with reduced water loss (; ). Transcription factors (TFs) are also at the center of N-by-W response. NLP7 is a master regulator of nitrogen signaling in Arabidopsis (). The nlp7 mutant shows drought resistant phenotype, similar to nrt1.1/chl1 (). Putting these findings together, it has been hypothesized that NLP7 regulates NRT1.1/CHL1 expression in guard cells and further controls stomatal opening and hence drought tolerance. Another TF in rice, drought and salt tolerance (DST), also bridges between N-assimilation and stomata movement that provides a path to crop improvement under marginal soil (lowN-lowW) ().
On the genome-wide level, our current manuscript explores on the gene regulatory networks (GRN) involved in N-by-W interactions by mining the N-by-W response RNA-seq and phenotype dataset from field grown rice (Swift et al., 2019). In our previous Swift et al 2019 study, we used linear models to discover that N-by-W signaling (N/W, molarity and/or NxW synergistic interactions) significantly correlate with rice field phenotypes, compared to genes that respond only to W-dose or N-moles (Swift et al., 2019). That dataset – which we use in our current analysis includes transcriptomic and phenotypic data for 19 rice varieties that vary in their drought and N-response. These 19 rice varieties were treated in a 2x2 N-by-W matrix of two N-doses (fertilized vs. without N) and W-doses (high vs. low water) in field experiments conducted at the International Rice Research Institute (IRRI) in the Philippines (Swift et al., 2019) (Figure 1). While our Swift et al., 2019 study determined the importance of the N-by-W gene responses (e.g., N/W and NxW) to phenotypic field outcomes in rice, the goal of our present study is to determine the GRNs underlying TF→target gene→phenotype interactions that correlate with NUE phenotypes in the rice N-by-W field study.
Figure 1
To develop sustainable agricultural solutions to feed a growing population, in this study we exploit a systems biology approach to uncover and validate the gene regulatory networks (GRNs) by which rice (Oryza sativa) plants sense and respond to the combination of N- and W- availability to promote crop productivity. To this end, we connected gene-to-NUE phenotype using weighted gene correlation analysis (WGCNA) (). Next, for the target genes that correlate with NUE phenotypes, we identified TF-to-target gene relationships in a gene regulatory network (GRN) using GENIE3 (). We then validated the TF-to-target gene network predictions via precision/recall (AUPR) analysis using validated TF-target gene data obtained in planta using the ConnecTF platform (https://rice.connectf.org). Additionally, we applied the plant cell-based Transient Assay Reporting Genome-wide Effects of Transcription factors (TARGET) system (; ), which we adapted in rice to validate the high-confidence TF-to-gene network for the N-by-W response genes whose expression level correlate with NUE.
Overall, we identified six TFs that regulate genes involved in both N and/or W signaling: OsbZIP23 (LOC_Os02g52780), Oshox22 (LOC_Os04g45810), LOB39 (LOC_Os03g41330), Oshox13 (LOC_Os03g08960), LOC_Os11g38870, LOC_Os06g14670. Two of these TFs are known regulators of drought tolerance - OsbZIP23 and Oshox22 – (Xiang et al., 2008; Zhang et al., 2012; ; Zong et al., 2016). Our present study shows that these TFs involved in drought responses are also responsive to N-by-W interactions. Moreover, we show that these six TFs control N-and/or-W response genes that correlate with NUEg. This information can now be applied to develop/breed rice plants with improved yield, on marginal, low N-input, drought-prone soils and on fields where water and N are limited due to climate change.
Materials and methods
Source of N-by-W response data (transcriptome and phenotype) for 19 rice varieties
Field phenotypic data collection and conditions for 19 rice varieties (Indica and Japonica) can be found in Swift et al., 2019 (Swift et al., 2019). The details of the treatments are in Swift et al., 2019, but as an overview: For the +N treatment, 150 kg/ha dose of (NH4)2SO4 was applied at 23 days after sowing (DAS). The -N treatment had no addition of fertilizer. Plants in the -W condition were covered from rain with a rainout shelter (intermittent watering was applied to ensure growth), while plants in the +W condition received rainfall and normal watering. Water-use-efficiency (WUE) was determined from leaves with carbon isotope discrimination as outlined in Swift et al., 2019 (Swift et al., 2019). The nitrogen usage data was calculated using the Kjeldahl N (KJ N) method which determined the nitrogen content from 1 gram of leaf samples. The total KJ N is determined as in (; ) by converting organic nitrogen forms to NH43+ and then measuring the concentration. To calculate N-uptake, we used the Kjeldahl N percent (KJ N%) and vegetative shoot dry weight (SDW) measurements from Swift et al., 2019 collected from leaf samples. We then used the N uptake measurement to calculate NUEg and NUE biomass (NUEb).
The field transcriptomic data consisted of 19 rice varieties (Indica and Japonica) of varying drought tolerant phenotypes, grown under four N-by-W treatment conditions, with three replicate leaf samples for RNA-seq for a total of 228 RNA-seq samples. Expression counts for 228 RNA-seq samples were normalized with the DESeq2 package (). The TFs and TF families from the N-and/or-W DE gene lists were identified based on the Plant Transcription Factor Database v4.0 categorization (). See data availability in Swift et al., 2019 (Swift et al., 2019) for source phenotypes and transcriptome data.
Potential index (IPO) calculation of NUE under low vs. high N and W conditions
To compare NUEg among the 19 rice varieties, we calculated the potential index (IPO) as similar to Ndiaye et al, 2019 (). For the calculation, each variety's NUEg was compared with the conditional average, using the formula below.
The IPO is the potential index of variety i; Yij is the NUEg of variety i for the condition j where j is HWHN, HWLN, LWHN or LWLN; is the conditional mean of all 19 varieties under condition j. The IPO is a relative value that shows the increase or decrease of a specific variety's NUEg, over the mean values. An IPO > 0 indicates better NUEg, whereas IPO< 0 indicates worse NUEg (Figure 2). The NUEg phenotype data was downloaded from Swift et al, 2019 (Swift et al., 2019).
Figure 2
WGCNA analysis: Gene-to-field phenotype correlation
The normalized counts files for each treatment and genotype were averaged as inputs into WGCNA to match the averaged field phenotypes for each biological replicate. This resulted in 76 transcriptomic and phenotypic values (19 varieties and 4 treatments) as inputs into WGCNA. The transcriptome counts file consists of counts for 22,436 genes in 76 samples. The R package, WGCNA, was used to perform the weighted correlation network analysis using step-by-step network construction and module detection (). We selected a MEDissThres of 0.5 to combine modules correlated with each other. We averaged the absolute value of the NUEg GS, WUE GS, and module membership (MM) scores for the genes in each module to select a cut-off value for highly correlated genes. (Figure 3C and Supplemental Figure 1). Overlapping module gene lists and N-and/or-W DE gene lists were made with Venny 2.1 web tool (). To determine the Z score and p-value of the NUEg and WUE genes that overlap with N-and/or-responsive DE gene lists, we used the Genesect function in Virtual Plant 1.3 () (Figures 3B, D and Supplementary Figure 1).
Figure 3
GENIE3 analysis of GRNs and validation of TF→ target gene predictions by AUPR and "network pruning"
The GENIE3 package in R () was used for network inference analysis. The gene expression data used to make the GENIE3 network consisted of the normalized counts of 228 RNA-seq samples for 10,815 N-and/or-W DE genes from Swift et al., 2019 (Swift et al., 2019) (Figure 1 Step 3). The 90 TFs for GENIE3 were selected from the two WGCNA modules (grey60 and skyblue) that are highly correlated with NUEg and are also N-and/or-W responsive (Figure 3D and Supplementary Data 4). The total unpruned network of 973,260 edges were uploaded to ConnecTF-Rice (rice.connectf.org) for network pruning and AUPR analysis (). This analysis is based on the in planta TF-target gene validation data for OsbZIP23, OsABF1, and OsNAC14 that is housed in the ConnecTF database () (Figure 4 and Supplementary Figure 2). Gene Ontology (GO) biological process analysis was conducted using g:Profiler (https://biit.cs.ut.ee/gprofiler/gost) with settings for only annotated genes and a significance threshold of 0.05 calculated with Benjamini-Hochberg FDR () (Table 1). For this analysis the gene IDs for target genes and genes associated with GO terms were converted between MSU7 and RAPDB gene designations. Cytoscape v3.9.1 was used for network visualization (Paul ) (Figure 5).
Figure 4
Table 1
| Rank. TF Name | Significant overlap of pruned GENIE3 target genes w/1,099 N-and/or-W DE genes in WGCNA modules(grey60&skyblue) #TF2s / #genes (Z score) | Relevant N and/or W GO terms associated with TF-target genes that overlap with N-and/or-W DE genes in WGCNA modules (grey60&skyblue) | TFs with High GS and MM for NUEg &/or WUE in WGCNA | Published TF Function (Reference) |
|---|---|---|---|---|
| 1. OsbZIP23 | 17 TF2s/159 genes (52.2) | "Response to water deprivation" | NUEg & WUE | Drought tolerance (Xiang 2008; ; Zong 2016) |
| 2. Oshox22 | 11 TF2s/93 genes (39.3) | "Response to water deprivation" & "Response to abscisic acid" | NUEg & WUE | Drought tolerance (Zhang et al., 2012) |
| 3. LOB39 | 5 TF2s/53 genes (30.9) | "Nitrate assimilation" | NUEg & WUE | N-responsive gene (; Yang 2017) |
| 4. Oshox13 | 5 TF2s/52 genes (27.6) | "Response to water deprivation" | NUEg & WUE | Unknown/Novel |
| 5. LOC_Os11g38870 | 0 TF2s/37 genes (25.9) | "Nitrate assimilation" | NUEg & WUE | Unknown/Novel |
| 6. LOC_Os06g14670 | 4 TF2s/49 genes (24.2) | "Response to water deprivation" & "Ammonia assimilation cycle" | NUEg & WUE | Unknown/Novel |
| 7. ERF65 | 7 TF2s/53 genes (32.1) | No N and/or W GO terms found | NUEg & WUE | Unknown/Novel |
| 8. OsERF48 | 6 TF2s/57 genes (27.3) | No N and/or W GO terms found | NUEg & WUE | Drought tolerance () |
| 9. OsIRO3 | 2 TF2s/24 genes (16.4) | No N and/or W GO terms found | NUEg & WUE | Iron homeostasis (Wang 2020) |
| 10. LOC_Os03g08470 | 1 TF2/20 gene (15.2) | No N and/or W GO terms found | NUEg & WUE | Unknown |
| 11. OsERF1 | 4 TF2s/25 genes (15.2) | No N and/or W GO terms found | NUEg & WUE | Ethylene response () |
| 12. OsABF1 | 5 TF2s/61 genes (13.6) | No N and/or W GO terms found | NUEg & WUE | Drought tolerance (Zhang 2017) |
| 13. OsIRO2 | 1 TF2/15 genes (13.3) | No N and/or W GO terms found | NUEg & WUE | Iron homeostasis/ N-signaling (; Ueda 2020) |
| 14. OSBZ8 | 1 TF2/19 genes (12.8) | No N and/or W GO terms found | NUEg & WUE | ABA response () |
| 15. RSR1 | 4 TF2s/18 genes (10.2) | No N and/or W GO terms found | NUEg & WUE | Starch biosynthesis () |
| 16. OsSPL9 | 0 TF2s/15 genes (10.0) | No N and/or W GO terms found | NUEg & WUE | Grain yield () |
| 17. EIL4 | 4 TF2s/40 genes (18.4) | No N and/or W GO terms found | NUEg | Unknown/Novel |
| 18. IDEF2 | 5 TF2s/96 genes (10.7) | No N and/or W GO terms found | NUEg | Iron homeostasis () |
| Total 52 TF2s/551 genes | "Response to water deprivation" & "Response to abscisic acid" |
Ranked list of 18 prioritized TFs that correlate with NUEg based on high-confidence edges to N-and/or-W DE genes in WGCNA modules (grey60 and skyblue).
First the TFs were ranked by the Z score for the overlap between the TF→target genes and the 1,099 N and/or W DE genes in WGCNA modules (grey60 & skyblue) associated with NUEg; second, for if the overlapping target genes for each TF were enriched for nitrogen and/or water GO terms; and third, for if the TF was highly associated with both NUEg and WUE.
Figure 5
Plasmid construction for TF-perturbation experiments using TARGET assay in plant cells
The coding sequences of OsABF1 and OsbZIP23 were determined as listed in Phytozome 13 () and were synthesized by GENEWIZ (South Plainfield, NJ) with the GATEWAY cassette for cloning into the p1107 destination plasmid (Supplementary Figure 4). Entry vectors were cloned into the p1107 plasmid using the LR Clonase II reaction according to manufacturer's instructions (Invitrogen). The p1107 plasmid for rice TARGET has a pBeaconRFP_GR () backbone with the following modification. The 35S promoters were replaced with maize Ubiquitin promoter subcloned from pTDM-C (Wu et al., 2016). A biotin ligase recognition peptide (BLRP) was fused at the N-terminal of the GATEWAY cassette, which is followed by the GR protein. All junctions were sequenced and verified for in frame TF-GR fusion proteins. The plasmid map and sequence (.FASTA) are provided in Supplemental Data File 1.
TARGET temporal TF perturbation experiment in rice leaf cells and RNA-sequencing
The rice protocol was adapted from our Arabidopsis TARGET protocol (; ) with some modifications. Rice seeds (Nipponbare) were sterilized by 70% ethanol for 3 mins followed by 50% commercial bleach for 30 min with rotation. The rice seeds were germinated in the dark for 4 days. The germinated rice seeds were transferred to ½ MS plates without sugar for 13 days in the growth chamber, under 16 h light/8 h dark diurnal cycle, at temperatures 27 and 25°C respectively and 70% humidity. On the day of the TARGET experiment, rice shoot tissue was cut into small (1 mm) pieces and stirred with cell-wall digestion solution (1.5% cellulase RS, 0.3% macerozyme R10 (Yakult Honsha), 0.6M mannitol, 10 mM MES (pH 5.7), 1 mM CaCl2, 5 mM b-mercaptoethanol, and 0.1% BSA) in a flask. The flask was vacuumed infiltrated for 20 minutes and shaken at 50 rpm in the dark for 4 hours. Rice shoot protoplasts were filtered through a 40 µm cell strainer (BD Falcon, USA) and spun down for 5 min at 500 g. The rice shoot protoplasts were then washed with 10 mL W5 solution (150 mM NaCl, 1M CaCl2, 1M KCl, 200 mM MES pH 5.7) three times, then resuspended in MMG solution (400 mM D-mannitol, 10 mM MgCl2, 4 mM MES pH 5.7) to 1.0x106 cells/mL. For protoplast transfection with vector, 1.0x105 cells were mixed with 40 µg plasmid DNA and 110 µL 40% PEG solution (40% 4000 PEG (Sigma, 81242), 400 mM D-mannitol, 50 mM CaCl2). The mixture was incubated at room temperature for 10 minutes. After incubation, the protoplasts were washed with W5 solution three times and resuspended in 1 mL WI solution (400 mM D-mannitol, 1M KCl, 200 mM MES pH 5.7). The transfected protoplasts were stored in the dark overnight. The next day, transfected protoplasts were treated with 30 µM cycloheximide (CHX) for 20 minutes (to block translation of secondary TF2 targets genes), before a three-hour 10 mM dexamethasone (DEX) treatment (to induce TF-GR nuclear import). After 3 hours, TF vector and control empty vector transfected protoplasts were FACS sorted for RFP signals into 150 µL TRI regent for RNA extraction (Zymo, R2061) (Supplementary Figure 5). We used Lexogen QuantSeq 3' mRNA-Seq Library Prep Kit FWD for Illumina (Lexogen, 015.2x96) for making RNA-Seq libraries. The libraries were pooled and sequenced on the Illumina NextSeq 500 platform at NYU-CGSB Genomics Core facility.
RNA-seq analysis of TARGET assay for validation of TF-target direct regulated genes
The UMI-incorporated RNA-Seq libraries of TF-transfected and empty vector control were analyzed following Lexogen's guidance (https://www.lexogen.com/quantseq-3mrna-sequencing/). The reads' UMI were extracted from raw fastq files using `extract` command from UMI-tools v1.1.1 (Smith et al., 2017). Then the fastq files were trimmed by fastp 0.21.0 (). The clean fastq files were aligned to MSU7 (Kawahara et al., 2013) genome using STAR 2.7.6a (). The aligned reads with the identical UMI were deduplicated using `dedup` command from UMI-tools v1.1.1 (Smith et al., 2017). The gene counts matrix was generated by featureCounts v2.0.1 () from the deduplicated bam files. The TARGET DE genes for OsABF1 and OsbZIP23 were identified using DESeq2 package () by comparing TF vs Empty Vector with a Benjamin & Hochberg adjusted p-values< 0.05. Differentially expressed (DE) genes identified for OsABF1 and OsbZIP23 are listed in Supplementary Data 9. Overlap between in planta and TARGET data was conducted with Venny 2.1 () and the significance was determined with Genesect in Virtual Plant 1.3 (). The calculations for precisions, recall and F-score for the GENIE3 network was the same as in () (Supplementary Figure 6).
Results
Phenotypic variation in NUEg in 19 rice varieties grown in N-by-W matrix field
The N-by-W response field data set used in our current study consisted of 19 rice varieties treated in a 2x2 matrix of four N-and/or-W treatment conditions (Figure 1) (Swift et al., 2019), comprising: well-watered (HW) with low-or-high N (HWLN, HWHN) (Figure 2A) vs. Low-W (LW) with low-or-high N (LWHN, LWLN) (Figure 2B) (For treatment details see Materials & Methods, and Swift et al., 2019. To refine our focus to NUEg, we examined how each of the 19 rice varieties performed for NUEg in the field (Figure 2). To identify the rice varieties with higher NUEg in the four different N-by-W field conditions, we adapted the Potential Index (IPO) () of NUEg for our N-by-W field dataset (Figure 2). The IPO for NUEg indicated the relative performance of each of the 19 rice varieties, compared to the conditional average (dotted lines). Under the well-watered (HW) condition, none of the rice varieties performed well under both LN and HN conditions (Figure 2A). For example, IR64 showed the highest NUEg values under HWLN, but only average NUEg values under HWHN conditions (Figure 2A). By contrast, Tainung67 showed the highest NUEg values under HWHN, but only average NUEg values under HWLN conditions (Figure 2A). However, under LW treatments, there was one variety, IR108, that performed well under both LWLN and LWHN conditions, with the highest IPO-NUEg compared to the other 18 varieties (Figure 2B). In line with this finding, the IR108 variety has been released under the variety name "Sukha dhan 5" to be used in the drought-prone regions of Nepal (). The IPO analysis reveals that this phenotypic dataset covers a range of rice NUEg values. Therefore, we used this NUEg phenotype data from the 2016 growing season data and the corresponding transcriptome data of Swift et al 2019, for the ensuing network-to-NUE phenotype analysis (Figure 1).
Identification of N-and/or-W responsive DE genes highly correlated with NUEg
To discover the relationships between genes and field phenotypes including NUEg, we used WGCNA () (Figure 1, Step 2, and Figure 3). The WGCNA analysis identified 11 co-expression modules for the 22,436 genes from the rice transcriptome data from the N-by-W field plot (Figure 1, Step 2, Figure 3A, Data in Supplemental Data 1). The genes in each of the WGCNA co-expression modules contribute to a Module Eigengene (ME) value based upon their Module Membership (MM) score. The MM score is the contribution of the individual gene to the ME value of the module (). We used the ME value to determine module correlation with each of the rice phenotypes from the N-by-W field plots (Figure 3A). The ME score for two WGCNA modules, grey60 (3,050 genes) and skyblue (744 genes) was significantly and highly correlated with the NUEg and WUE phenotype data in the N-by-W plot (Figure 3A). The ME value of the grey60 module was negatively correlated with NUEg (-0.71), while the ME value of the skyblue module was positively correlated with NUEg (+0.73) (Figure 3A). However, each WGCNA module contains subsets of genes that can be either positively or negatively correlated with NUEg. In WGCNA, this gene expression-to-phenotype correlation is called Gene Significance (GS), as shown for the plot of MM vs. GS in Figure 3C.
To identify which WGCNA modules had a significant representation of genes responding to N-and/or-W signals, the genes comprising each module were overlapped with the N-and/or-W responsive DE genes from Swift et al 2019 (Swift et al., 2019) (Figure 3B). This analysis uncovered a significant overlap of the N:W- and W- responsive gene lists with the genes in the WGCNA modules - grey60 and skyblue - which are each highly correlated with NUEg and WUE (Figure 3A). This demonstrates that the genes in the WGCNA modules - grey60 and skyblue - not only correlate with the NUEg phenotypes from the N-by-W matrix field plots but are also enriched in genes responsive to N-and/or-W signals (Figure 3B). Additionally, the blue and lightyellow WGCNA modules are enriched in genes that respond to N-moles, but not to the interaction of N and W. While the WGCNA modules - blue and lightyellow - do not correlate significantly with NUEg, each of these modules correlates significantly with chlorophyll concentration (Figure 3A), a trait known to be regulated by N and used to determine N-status and the need for fertilizer in the field ().
Next, we performed two analyses that enabled us to prioritize the N-and/or-W response DE TFs and genes within each of the two WGCNA modules - grey60 and skyblue - that are most highly correlated with the NUEg phenotype (Figures 3C, D). The genes with MM scores closes to -1 or 1 are highly connected to their WGCNA module. In addition, genes with GS scores that have a high absolute value for a specific trait are also more biologically significant (). Therefore, to filter genes in each module that were highly correlated with NUEg, we identified genes with high absolute values for both their MM and GS scores. To do this, we first plotted the absolute values of the MM vs. GS scores for each gene in the WGCNA modules - grey60 and skyblue - which are highly correlated with NUEg (Figure 3C). Next, we calculated the average MM and GS scores for the genes in each of these two modules. This enabled us to set a cut-off and identify genes whose absolute MM and GS values were great than or equal to the average of the genes in each module (Figure 3C, upper right quadrant).
This analysis identified a combined total of 131 TFs and 1,491 genes highly relevant to NUEg in the two WGCNA modules: grey60 (104 TFs & 1,209 genes) and skyblue (27 TFs & 282 genes) (Figure 3C). Next, to identify whether genes highly relevant to NUEg are significantly enriched in N-and/or-W response gene, we performed a Genesect analysis () (Figure 3D). This analysis revealed significant overlaps between the N:W and W responsive gene lists from Swift et al 2019 (Swift et al., 2019), with the genes highly correlated with NUEg (131 TFs and 1,491 genes) from the combined grey60 and skyblue WGCNA modules (Figure 3D). The resulting overlap consisted of 90 TFs and 1,099 genes that are highly associated with NUEg and N-and/or-W responsive (Supplementary Data 2). Next, we determined which of these TFs and genes correlated NUEg were also highly associated with the WUE phenotype. To do this, we conducted the same analysis pipeline as described above for NUEg, in which we determine a new GS cut off value for WUE (Supplementary Figure 1A). This resulted in 79 TFs and 976 genes that are highly correlated with WUE and are N-and/or-W responsive (Supplementary Figure 1B, Supplementary Data 3). We find that 72 (80%) NUEg TFs and 815 (74%) NUEg genes are also highly correlated with WUE, thus suggesting a dual role for these genes/TFs in regulating both N and W responses.
For further analysis, we prioritized 90 TFs from the GENIE3 analysis that are; i) N-and/or-W responsive and ii) highly correlated to NUEg from the combined WGCNA modules - grey60 and skyblue. This analysis resulted in 29 TFs that are N:W-responsive and 61 TFs that are W-responsive (Figures 3C, D).
Validation of TF→target GRN predictions in WGCNA modules associated with NUEg
To predict TF→target gene interactions in GRNs important for NUEg, we used GENIE3, a random forest network inference method (). This analysis will identify potential master TF regulators of the NUEg response amongst the 90 TFs (29 TFs N:W-responsive and 61 TFs W-responsive) (Figure 3D) that are highly correlated with NUEg (e.g., members of WGCNA grey60 and skyblue models) (Figures 3A, C). To identify and rank these 90 TFs from these NUEg modules, we generated a GRN using 90 potential TF-regulators of 10,815 DE (N-and/or-W response genes) from the field N-by-W matrix (Figure 1, Step 1). The output of GENIE3 ranks the TF→target gene predictions in the order of confidence for each of the 90 TFs and the 10,815 DE genes N-and/or-W responsive (Figure 4). In total, the resulting GENIE3 inferred network ranks numerical confidence scores for each TF and target gene, excluding self-regulation of the TF ((90 TF x 10,815 genes) - 90 TFs) = 973,260 TF-target edges (Figure 4 and Supplemental Data 4).
Our next goal was to validate the TF-target gene interactions in our predicted GRN, using TF-target gene data validated in planta. To this end, we used experimentally validated TF-target gene interactions from TF perturbation data in rice, housed in the ConnecTF platform (https://rice.connectf.org) () (Figure 4 and Supplementary Figure 2). The ConnecTF database includes published rice RNA-seq and ChIP-seq data available as of June 2020. To validate the GRN, we uploaded the TF→target gene interactions predicted by the GENIE3 network into ConnecTF and filtered for validated TF-regulation (RNA-seq) and TF-binding (ChIP-seq) data from rice in planta datasets (Figure 4 and Supplementary Figure 2, Supplementary Data 4). We focused our analysis on validated TF-target gene datasets from rice leaf tissue, given that the source RNA-seq data used to make the GENIE3 network was from rice leaves (Supplementary Figure 2).
Our query of the ConnecTF rice TF database identified experimental TF-target gene validation datasets for three TFs in rice leaf tissue from our GENIE3 network (Figure 4 and Supplemental Figure 2). The three validated rice TFs are OsABF1 (Zhang et al., 2017), OsbZIP23 (Zong et al., 2016), and OsNAC14 (Shim et al., 2018). These three validated rice TFs include a total of nine datasets with 10,941 validated target genes from TF-regulation and/or TF-binding data (Figure 4 and Supplementary Figure 2). We then used this in planta data as "gold-standard" data to validate the TF→target gene predictions from our GENIE3 network using Area Under the Precision Recall (AUPR) curve analysis, which is an automated function in the ConnecTF platform (Figure 4). The results show that the AUPR for the TF→target gene predictions (edges) in the rice GENIE3 network were significantly higher than the random TF-target gene edges (P-value<0.001, permutation test) (Figure 4). Given the AUPR curve, we were able to select a precision threshold of 0.31 (e.g., TF→target gene edge score ≥ 0.0581). This cut-off score is equivalent to the TF→target gene predictions being accurate 1/3 of the time and this level of accuracy is comparable to other similar network validation AUPR studies (Varala et al., 2018; ). The GENIE3 network was then pruned for only the high-confidence TF→target gene predictions using this precision cut-off score. This network pruning for precision, resulted in a GRN containing 8,826 high confidence edges connecting 88 TFs and 5,716 target N-and/or-W response DE genes (Figure 4 and Supplementary Data 5).
Prioritization of master TFs that regulate NUEg in response to N-and/or-W signaling
Our next goal was to prioritize candidate N-and/or-W response TFs with a significant influence on NUEg from the pruned GENIE3 network. To this end, we overlapped the pruned high confidence TF→target edges for the 88 TFs in the GENIE3 network with the 1,099 genes from the two WGCNA modules that are highly correlated with NUEg - grey60 & skyblue - N-and/or-W DE genes = 322 N:W response genes + 777 W-response genes) (Supplementary Figure 3). We calculated the significance of the overlapping TF→target genes with the 1,099 NUEg genes. To prioritize the 88 TFs, we ranked them by the Z-score for the overlap (Supplementary Data 6). We found 18 TFs whose high confidence TF→targets gene edges had the highest significant overlap (P-value<0.001, Z score ≥ 10) with the 1,099 genes in the NUEg WGCNA modules – grey60 and skyblue (Table 1). This analysis links 18 TFs→ 551 N-and/or-W response target genes→NUEg. Among the 18 TFs, OsbZIP23 is predicted to regulate the most of the NUEg correlated genes, compared to the other 17 TFs (Table 1). Additionally, we find that 16/18 TFs (all except EIL4 and IDEF4) are also highly corelated with WUE (Table 1 and Supplemental Data 3).
Of these 18 TFs, multiple TFs have published functions in drought tolerance including, OsABF1 (Zhang et al., 2017), OsbZIP23 (Xiang et al., 2008; ; Zong et al., 2016), Oshox22 (Zhang et al., 2012), and OsERF48 (). Of note, OsABF1, OsbZIP23, and Oshox22 are N:W-responsive genes based on the N-and/or-W response DE gene lists from Swift et al 2019 (Supplemental Data 6), suggesting their new function in regulating N:W responses, in addition to drought (Table 1). Published functions for other candidate TFs in the 18 TF list include, N-signaling (LOB39) (; Yang et al., 2017), ABA signaling (OSBZ8) (), ethylene signaling (OsERF1) (), iron homeostasis (IDEF2, OsIRO3, and OsIRO2) (, ; ; Wang et al., 2020), starch biosynthesis (RSR1) (), and grain yield (OsSPL9) () (Table 1). OsIRO2 was also found to regulate NUE in a N-response gene network in rice (Ueda et al., 2020).
Gene ontology for target genes for prioritized TFs
To further determine the mechanism of the prioritized TFs in regulating NUEg, we performed Gene Ontology (GO) analysis on the NUEg target genes from Table 1 regulated by each TF using g:Profiler (Table 1 and Supplemental Data 7) (). For each TF, we focused on the relevant biological process GO terms related to water and nitrogen signaling. We found that the target genes of the TFs, LOB39, LOC_Os11g38870, and LOC_Os06g14670, were enriched for GO terms related to nitrogen including, "nitrate assimilation," and "ammonia assimilation cycle" (Table 1). Further, we found that the target genes of the TFs, OsbZIP23, Oshox22, Oshox13, LOC_Os06g14670, were enriched for GO terms related to drought including, "response to water deprivation," and "response to abscisic acid" (Table 1). LOC_Os06g14670 was the only TF enriched for nitrogen and drought-related GO terms. We did not identify any GO enrichment for the TF-target genes of OsERF48, OsIRO3, LOC_Os03g08470, OSBZ8, RSR1 and IDEF2. However, we did identify some other GO terms of interest for the remaining TFs including, "sulfur compound metabolic process" for EIL4, "cell communication" ERF65, "response to temperature stimulus" for OsABF1, "phosphorus metabolic process" for OsERF1, "iron ion homeostasis" for OsIRO2, and "zinc ion homeostasis" for OsSPL9 (Supplemental Data 8). While these enriched GO terms suggest the relevance of these TFs in other cell processes, we focus on the TFs that regulate the target genes associated with the nitrogen and water related GO terms.
High-confidence GRN of TFs that target nitrogen and drought-related genes
To identify the TFs that regulate both nitrogen and water response from our list of prioritized TFs, we took the subset of the GENIE3 network that includes 18 TFs→ 551 N-and/or-W response target genes associated with NUEg, and identified the target genes from this list of 551 that were part of the GO terms "nitrate assimilation", "ammonia assimilation cycle", "response to water deprivation," and "response to abscisic acid" (Supplemental Data 7). This resulted in a list of 23 target genes regulated by 14 TFs (Supplemental Data 7). We found six TFs that regulated both nitrogen and water related target genes either directly (OsbZIP23, LOB39, LOC_Os11g38870, LOC_Os06g14670, and Oshox13) or indirectly (Oshox22 via regulation of OsbZIP23) (Figure 5). While OsABF1 did not regulate genes related to nitrogen, it is included in the network visualization because it is annotated for the water-related GO terms and is regulated by OsbZIP23 and Oshox22 (Figure 5).
The target genes involved in nitrate and ammonia assimilation that are regulated by the TFs in our high-confidence GRN include validated regulators of NUE, glutamate synthetase 1 (OsGOGAT), and nitrite reductase (OsNiR) (; Yu et al., 2021) (Figure 5). We also find regulation of the putatively expressed nitrate reductase 1 (NIA1) gene, which is necessary for nitrate assimilation (Subudhi et al., 2020). The TFs, OsbZIP23, LOB39 and LOC_Os11g38870 regulate nitrate assimilation genes, while OsbZIP23, Oshox13, and LOC_Os06g14670 regulate the ammonia assimilation gene. OsbZIP23 is the only TF that regulates genes in both nitrate and ammonia assimilation genes (Figure 5).
Furthermore, each TF regulates genes involved in water deprivation and/or ABA signaling (Figure 5). These genes include the TFs OsbZIP46 and OsbZIP72, which are known positive regulators of drought tolerance and function in coordination with OsbZIP23 and OsABF1, two other prioritized TFs in our network (; ; Zhang et al., 2017; Song et al., 2020). We also find regulation of the rice aquaporins, OsPIP1;1, OsPIP1;2, and PIP2A that facilitate water transport (; Xu et al., 2019). In addition, there are genes that regulate multiple components involved in the ABA signaling pathway including, the ABA drought receptors, OsPYL1, OsPYL6 (; ), the clade A type 2C protein phosphatases, OsPP2C51, OsPP2C30 (Zong et al., 2016; ), and the ABA-activated protein kinase, SAPK6 (). Overall, this result demonstrates that a subset of our prioritized candidate TFs regulates both nitrogen and water genes.
Network validation with in vivo TARGET assay
To further validate the nitrogen and drought-related edges in our high-confidence GRN (Figure 5), we performed in vivo Transient Assay Reporting Genome-wide Effects of Transcription factors (TARGET) assays to identify the direct TF-target genes for these TFs. We selected OsbZIP23 and OsABF1 for TARGET assays because we could compare the accuracy of our TARGET results with the available in planta data for these TFs in ConnecTF (). The TARGET TF-perturbation assay in isolated plant cells has been previously used to identify direct TF→regulated target genes in Arabidopsis (; Varala et al., 2018; ). In this paper, we adapt the vectors and the TARGET temporal TF-perturbation assay to rice shoot cells (Supplementary Figure 4).
The TARGET TF-perturbation assay identifies the direct TF→ regulated target gene interactions because; i) there is timed nuclear entry of the TF, and ii) translation of regulated secondary (TF2) transcription factors is blocked by cycloheximide treatment. TF-regulated DE genes are identified by comparison to an empty vector control. The TARGET assay identifies direct TF→target genes as follows: i) the TF is fused to the glucocorticoid receptor (GR) protein that when expressed in the plant cells, ii) the TF-GR fusion is retained in the cytoplasm by HSP90 binding, iii) upon dexamethasone (DEX) treatment, the GR binding is released and the TF is imported into the nucleus where it can regulate expression (; ) (Supplementary Figure 5). iv) Additionally, cycloheximide + DEX treatment inhibits translation of mRNA for TF2s. Therefore only the target genes of the over-expressed TF are identified, when compared to the empty vector control ().
Based on our TARGET assay, OsbZIP23 directly regulates 3,095 target genes, while OsABF1 directly regulates 2,151 target genes in rice shoot protoplasts (Supplementary Figure 6 and Supplemental Data 9). To determine the accuracy of our TARGET results, we took the overlap between the TARGET results and the in planta binding and expression data for each TF from ConnecTF (Zong et al., 2016; Zhang et al., 2017; ). We found a significant overlap between the TARGET and in planta DE genes (Supplemental Figure 6A). This significant overlap suggests that the plant cell-based TF-target data can accurately identify in planta TF-regulated genes. Additionally, we find the TARGET data is as accurate, if not even better, than the in planta data at validating the predicted TF→target genes in the GENIE3 network, with a higher F-score and similar precision and recall values (Supplementary Figure 6B).
Given that the TARGET data was accurate in identifying OsbZIP23 and OsABF1 target genes, we used the TARGET and in planta data to validate the nitrogen and drought-related edges in our high-confidence GRN (Figure 5). We confirm with TARGET that OsbZIP23 directly regulates genes involved in nitrogen and drought responses including, NIA1 involved in nitrate assimilation (Subudhi et al., 2020), ABCG4 involved in abiotic stress responses (), and the rice aquaporin, OsPIP1;2, that improves yield in rice (Xu et al., 2019). Additionally, we confirm with OsbZIP23 TARGET and in planta data that OsbZIP23 regulates drought associated genes OsDhn1 and OsPP2C30 (; ). Furthermore, we confirm the role of OsABF1 in regulating drought signaling, as it regulates the drought-associated gene OsPP2C51 in both TARGET and in planta datasets (Figure 5) (Zong et al., 2016).
Overall, our TARGET results show that the high-confidence edges inferred in our GENIE3 network accurately predict TF→target genes, thus further confirming the role of OsbZIP23 in regulating both NUEg and WUE. In addition, we find a new function for OsbZIP23 in mediating NUEg phenotypes, as previous studies show its role in drought responses (Xiang et al., 2008; ; Zong et al., 2016). Thus, our combined network inference and validation approach reveals new TFs in regulating NUEg (Table 1).
Discussion
In this study, we sought to identify GRNs that control NUEg in response to two key interacting components that regulate rice productivity: N and W. By exploiting transcriptomic and phenotypic data collected from 19 rice varieties grown in a 2x2 N-by-W matrix in the field (Swift et al., 2019), we identified and validated the role GRNs comprised of N-and/or-W response genes for their role in TF→target gene→ NUEg phenotype relationships. The TF to N-by-W response gene information now encoded in this high-confidence GRN correlated to NUEg, can now be applied to develop/breed rice plants with improved yield marginal, low N-input, drought-prone soils – which are increasing in the face of climate change.
High-confidence GRN identifying master regulators of NUEg responsive toN-and/or-W signals
We were able to link the TF→target gene→NUEg phenotype using a combination of four approaches (i) WGCNA-based gene-to-trait co-expression network (), (ii) GENIE3, a random forest machine learning approach to GRN inference for predicting TF-target interactions (), (iii) Experimental validation of GRN predictions and Network "pruning" by AUPR (Varala et al., 2018; ), and (iv) Network validation using TARGET, an approach which uses plant cells to identify direct TF→target gene interactions (; ). Using this pipeline (Figure 1), the WGCNA approach identified two network modules that were highly correlated to NUEg called "grey" and "skyblue". Next, we constructed a GRN for the genes in this module, based on their N-and/or-W response DE genes. Finally, we used experimental data for TF-target genes validated in planta (Zong et al., 2016; Zhang et al., 2017; Shim et al., 2018) as well as ones we generated in rice leaf cells for this study. These validated rice TF datasets were used to conduct precision/recall analysis of our GRN. This enabled us to set a precision cut-off score to prune the network for high confidence TF-target predictions for all TFs in the GRN.
Overall, our GRN analysis and validation identified OsbZIP23 and Oshox22 as top candidate master regulators of NUEg in response to N and W signaling. These two TFs are network hubs, as they regulate the largest number of DE genes (N-and/or-W responsive) that are highly correlated with NUEg in the grey60 and skyblue WGCNA modules (Table 1 and Supplemental Data 6). Further validating their known role in drought, these two TFs have published functions in regulating drought tolerance through the plant hormone abscisic acid (ABA) signaling responses (Xiang et al., 2008; Zhang et al., 2012, 2017; ) (Table 1). Our current study, now links these two well-known drought TFs to regulation by N-and/or-W signaling and NUEg. Our results are also in line with previous studies that show OsbZIP23 activity to be dependent upon phosphorylation by SAPK2 (Zong et al., 2016), an osmotic stress/ABA-activated protein kinase, which promotes nitrate uptake and assimilation under drought stress ().
In addition to the TF hubs (OsbZIP23 and Oshox22), we identify four TFs with novel functions NUEg and WUE gene regulation in our GRN. We identified four TFs (LOB39, Oshox13, LOC_Os11g38870, and LOC_Os06g14670), that regulate genes involved in both N and/or W responses using GO analysis of their predicted TARGET genes in the high-confidence GRN (Table 1 and Figure 5). Unlike OsbZIP23 and Oshox22, the TFs Oshox13, LOC_Os11g38870, and LOC_Os06g14670TFs had until now unknown functions in both nitrogen and drought regulation (Table 1). LOB39 expression is regulated by nitrogen, however it was previously not known to be involved in drought (). OsbZIP23, LOB39 and LOC_Os11g38870 regulate nitrate assimilation genes NIA1 and OsNiR, which is a known to promote nitrogen assimilation and NUE in coordination with OsNLP4 (Figure 5) (Yu et al., 2021). Furthermore, OsbZIP23, Oshox13 and LOC_Os06g14670 regulate the ammonia assimilation gene OsGOGAT1, which improves NUE in low N conditions in coordination with the ammonium transporter OsAMT1;2 (). While it is known that rice prefers ammonia uptake compared to nitrate (Sasakawa and Yamamoto, 1978; ), we find the TFs in this network regulate both pathways, with OsbZIP23 regulating genes involved in both.
We also examined the mechanism of transcriptional regulation between these master TFs in the NUEg GRN by validating TF→target gene interactions using TARGET, a plant cell-based assay that identifies direct TF→TARGET gene interactions (; Varala et al., 2018; ). We find that Oshox22 regulates nitrogen and water responses indirectly via candidate TFs OsbZIP23, and OsABF1 (Figure 5). We then validate the TF→target gene interactions for OsbZIP23 and OsABF1 TFs with the TARGET assay. We confirm that OsbZIP23 regulates both nitrogen and drought response genes, and OsABF1 regulated drought response genes, with TARGET and in planta data.
Overall, these finding supports previous studies that show the regulation of these two essential signals N-and-W are linked (Swift et al., 2019; ; ). Our work presents a path of how ABA/drought induced signaling regulates both N and W responses which ultimately affect crop phenotypes, such as NUEg, the trait of focus in our study.
Validation of GRNs in rice using ConnecTF as a platform to validate and prune for high-confidence networks
In our study, we demonstrate the usefulness of ConnecTF as a platform - now applied to rice - to integrate published TF-binding and TF-expression datasets to identify and validate target genes in GRNs () (Figure 4 and Supplementary Figure 2). While some GRN studies use an arbitrary cut-off value for network pruning as in other network studies (Ueda et al., 2020), we show how TF-perturbation data can be used as a "gold-standard" for GRN validation and "network pruning", using automated AUPR functions in ConnecTF () (Figure 4). We performed Precision/Recall analysis of the GRN for NUEg – using the TF-target gene validation sets for rice housed in the ConnecTF database. This enabled us to empirically select a TF→target precision cut-off value of 0.31 from the AUPR curve. This AUPR cut-off represents that approximately 1/3 of our GENIE3 network predictions are validated (Figure 4). This precision cut-off is comparable to what we find in other network studies in Arabidopsis that use AUPR analysis (Varala et al., 2018; ). Overall, the automated AUPR function in ConnecTF provides an accurate, and facile means to validate GRN predictions in any rice GRN that researchers can load onto the site. Importantly, these cut-off values for TF→target gene validated edges established a cut-off score that can be applied to all TF→target gene edges in the network – including TFs which have not been validated. This enables the generation of a high-confidence network for all TFs in the GRN.
bZIP family TFs as regulators of N and W signaling
In our high-confidence GRN we identify nine bZIP TFs as regulators in our "pruned" network (Supplementary Data 6). Members of the bZIP family of TFs are known to regulate drought stress responses in multiple crops species in addition to rice, including Glycine max, Zea mays and Hordeum vulgare (). Additionally, bZIP family TFs regulate ABA hormone responses, which play a crucial role in regulating the drought response in plants in general (; Zong et al., 2016; Zhang et al., 2017; ). In our high-confidence GRN studies that focus on genes correlated with NUEg, we find that bZIP TFs regulate N-signaling as well as drought responses in rice. In line with our finding, previous studies examining N-responses in rice, identified bZIP transcription factors that regulate NUE (Ueda et al., 2020).
We identified three bZIP family members - OsABF1, OsbZIP23, and OSBZ8 - as top-regulators of N-and/or-W signaling in regulating NUEg (Table 1). Additionally, we find regulation of two other bZIP TFs, OsbZIP72 and OsbZIP46, in our NUEg GRN regulated by Oshox22 and OsbZIP23, respectively (Figure 5). This finding is significant, as OsbZIP23, OsbZIP46, OsbZIP72 are part of the same subgroup-III of bZIP TFs and are known to be coordinated in their regulation of ABA signaling and drought responses (; ; Song et al., 2020). Additionally, ObZIP46 improves drought tolerance in coordination with the ABA-activated protein kinase, SAPK6, which is another target gene in our NUEg GRN (Figure 5) (). Overall, our NUEg GRN results link bZIP TFs in rice as mediating N-and/or-W response genes that control NUEg. We validate the TF→target genes predictions in our high-confidence GRN for NUEg for two bZIP TFs, OsbZIP23 and OsABF1, using the TARGET assay.
Functional validation of TFs in rice: TARGET assay to identify direct TF→target gene interactions in rice cells
The TARGET system allows researchers to identify the validated TF-target gene interactions for any TF of interest using a rapid plant cell based temporal TF perturbation assay (; ). The key to this assay is the inducible TF nuclear localization and its ability to identify direct TF-target genes based on RNA-seq data (). Previously, the TARGET assay has been used to identify direct TF→target gene interactions in Arabidopsis root or shoot cells (; Varala et al., 2018; ). In this study we establish the TARGET system in rice leaf protoplasts (see Methods). We then used the rice TARGET assay to identify the direct regulated target genes of the rice TFs OsbZIP23 and OsABF1 (Supplementary Data 9). Our analysis shows that the TF target genes identified in rice leaf protoplasts using TARGET, are comparable and show a significant overlap with genes identified in planta (Supplementary Figure 6A). Additionally, in this study, we demonstrate that the accuracy of rice TARGET data is comparable to in planta data at validating network predictions (Supplementary Figure 6B). This finding suggests that rice TARGET data can be used to validate GRN predictions in rice, as was shown in Arabidopsis (Varala et al., 2018; ; ; ). In our study, we validated that OsbZIP23 regulates both nitrogen and water-related genes including, NIA1 which is involved in nitrate assimilation (Subudhi et al., 2020), OsDhn1 which is induced by drought (), OsPIP1;2 which is an aquaporin that improves yield (Xu et al., 2019), ABCG4 which is involved in abiotic stress responses (), and OsPP2C30 which a core regulator in the ABA signaling pathway (Zong et al., 2016). Overall, our study supports that the TARGET assay is a fast and reliable approach to identify the direct TF→target genes in rice, bypassing the time-consuming process of developing transgenic rice. Importantly, the rapid rice TARGET TF-perturbation assay, can be used to prioritize rice TFs for more laborious studies in planta.
Our network approach is transferrable to any phenotype in any organisms
The method we applied in this study relies on two inputs: a transcriptome-wide gene expression table and collected phenotypes from the same samples. With the reduced cost of RNA-Seq, especially with the 3′ RNA-sequencing (Weih, 2014; Groen et al., 2020; Weng and Juenger, 2022), it is much more feasible for researchers to obtain transcriptome expression data from many samples. Moreover, the software we used are all open-source and publicly available. This includes WGCNA () for gene-to-phenotype correlation, GENIE3 () for GRN inference and ConnecTF () for network pruning. Putting these together, our network approach is not limited in rice research, but can be applied to any organism for any phenotype or trait.
Conclusions
By using a combination of WGCNA and GENIE3 network methods, we present a gene regulatory network that links TF→target gene→NUEg phenotype to determine the mechanism of N-and/or-W signaling to the regulation of NUEg (Figure 1). We also show how to use TF-validation datasets from rice to validate inferred networks using ConnecTF (https://rice.connectf.org) (). In addition, we apply the cell-based TARGET temporal TF-perturbation system to rice to identify direct TF→target genes interactions and validate inferred gene networks. Overall, we identify a new role for OsbZIP23 and Oshox22 as regulators of the N-and/or-W signaling and regulation of NUEg, in addition to ABA/drought signaling. More broadly, we have identified 18 prioritized TFs and their targets that correlate with NUEg, and results from this network approach can potentially be used to optimize rice varieties to thrive in marginal low-N/arid soils, which are increasing in the face of global climate change.
Funding
This work is supported by NSF-PGRP IOS-1840761 to GC, a Grant from the Zegar Family Foundation (A16-0051) to GC, an NIH Grant RO1-GM121753 to GC, an NIH NIGMS Fellowship F32GM139299 to CS, and JS is an Open Philanthropy Awardee of the Life Sciences Research Foundation.
Acknowledgments
We thank the staff at IRRI for their work on the field studies. We would also like to thank Dr. Manpreet Katari and Will Hinkley for their advice and sharing code for data analysis.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
The data presented in the study are deposited in NCBI repository, BioProject: PRJNA828338.
Author contributions
CS, JH, C-YC, and GC designed the research experiments. CS, JH, C-YC, and H-JS, performed research experiments. JS, AH, MB, VA, and JA contributed data and analysis. CS, JH, C-YC, JS, AH, and GC wrote and edited the paper. All authors contributed to the article and approved the submitted version.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.1006044/full#supplementary-material
References
1
AlvarezJ. M.SchinkeA.-L.BrooksM. D.PasquinoA.LeonelliL.VaralaK.et al. (2020). Transient genome-wide interactions of the master transcription factor NLP7 initiate a rapid nitrogen-response cascade. Nat. Commun.11, 1–13. doi: 10.1038/s41467-020-14979-6
2
AnanthaM. S.PatelD.QuintanaM.SwainP.DwivediJ. L.TorresR. O.et al. (2016). Trait combinations that improve rice yield under drought: Sahbhagi dhan and new drought-tolerant varieties in south Asia. Crop Sci.56, 408–421. doi: 10.2135/cropsci2015.06.0344
3
ArausV.SwiftJ.AlvarezJ. M.HenryA.CoruzziG. M. (2020). A balancing act: how plants integrate nitrogen and water signals. J. Exp. Botany71(15), 4442–445. doi: 10.1093/jxb/eraa054
4
BargmannB. O. R.Marshall-ColonA.EfroniI.RuffelS.BirnbaumK. D.CoruzziG. M.et al. (2013). TARGET: A transient transformation system for genome-wide transcription factor target discovery. Mol. Plant6, 978–980. doi: 10.1093/mp/sst010
5
BremnerJ. M. (1996). Nitrogen-total In SparksD.PageA.HelmkeP.LoeppertR.SoltanpourP.N.TabatabaiM.A.JohnstonC.T.et al Eds. Methods of Soil Analysis. (Sparks, DL). doi: 10.2136/sssabookser5.3.c37
6
BremnerJ.MulvaneyC. (1982). “Nitrogen-total,” in MillerA. L.KeeneyR.H.D. R. Eds. Methods of soil analysis. part 2. chemical and microbiological propertiesMadison, Wisconsin: American Society of Agronomy, Soil Science Society of America, 595–624.
7
BrooksM. D.CirroneJ.PasquinoA. V.AlvarezJ. M.SwiftJ.MittalS.et al. (2019). Network walking charts transcriptional dynamics of nitrogen signaling by integrating validated and predicted genome-wide interactions. Nat. Commun.10, 1–13. doi: 10.1038/s41467-019-09522-1
8
BrooksM. D.JuangC.-L.KatariM. S.AlvarezJ. M.PasquinoA. V.ShihH.-J.et al. (2020). ConnecTF: A platform to build gene networks by integrating transcription factor-target gene interactions. bioRxiv2020, 191627. doi: 10.1101/2020.07.07.191627
9
CastaingsL.CamargoA.PocholleD.GaudonV.TexierY.Boutet-MerceyS.et al. (2009). The nodule inception-like protein 7 modulates nitrate sensing and metabolism in arabidopsis. Plant J.57, 426–435. doi: 10.1111/j.1365-313X.2008.03695.x
10
ChangY.NguyenB. H.XieY.XiaoB.TangN.ZhuW.et al. (2017). Co-Overexpression of the constitutively active form of OsbZIP46 and ABA-activated protein kinase SAPK6 improves drought and temperature stress resistance in rice. Front. Plant Sci.8. doi: 10.3389/fpls.2017.01102
11
ChenS.ZhouY.ChenY.GuJ. (2018). Fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics34, i884–i890. doi: 10.1093/bioinformatics/bty560
12
CirroneJ.BrooksM. D.BonneauR.CoruzziG. M.ShashaD. E. (2020). OutPredict: multiple datasets can improve prediction of expression and inference of causality. Sci. Rep10 (1), 6804. doi: 10.1038/s41598-020-63347-3
13
DeyA.SamantaM. K.GayenS.SenS. K.MaitiM. K. (2016). Enhanced gene expression rather than natural polymorphism in coding sequence of the OsbZIP23 determines drought tolerance and. PLoS One11, 1–26. doi: 10.1371/journal.pone.0150763
14
DobinA.DavisC. A.SchlesingerF.DrenkowJ.ZaleskiC.JhaS.et al. (2013). STAR: Ultrafast universal RNA-seq aligner. Bioinformatics29, 15–21. doi: 10.1093/bioinformatics/bts635
15
D’OdoricoP.ChiarelliD. D.RosaL.BiniA.ZilbermanD.RulliM. C. (2020). The global value of water in agriculture. Proc. Natl. Acad. Sci. United States America117, 21985–21993. doi: 10.1073/pnas.2005835117
16
FageriaN. K.BaligarV. C.JonesC. A. (2010). Growth and mineral nutrition of field crops (3rd ed.). (CRC Press). doi: 10.1201/b10160
17
FuF. F.XueH. W. (2010). Coexpression analysis identifies rice starch regulator1, a rice AP2/EREBP family transcription factor, as a novel rice starch biosynthesis regulator. Plant Physiol.154, 927–938. doi: 10.1104/pp.110.159517
18
GibbsH. K.SalmonJ. M. (2015). Mapping the world’s degraded lands. Appl. Geogr.57, 12–21. doi: 10.1016/j.apgeog.2014.11.024
19
GoodsteinD. M.ShuS.HowsonR.NeupaneR.HayesR. D.FazoJ.et al. (2012). Phytozome: A comparative platform for green plant genomics. Nucleic Acids Res.40, 1178–1186. doi: 10.1093/nar/gkr944
20
GroenSCĆalićIJoly-LopezZPlattsAEChoiJYNatividadMet al. The strength and pattern of natural selection on gene expression in rice. Nature (2020) 578(7796):572–6. doi: 10.1038/s41586-020-1997-2
21
GuoF.-Q.YoungJ.CrawfordN. M. (2003). The nitrate transporter AtNRT1. 1 (CHL1) functions in stomatal opening and contributes to drought susceptibility in arabidopsis. Plant Cell15, 107–117. doi: 10.1105/tpc.006312
22
HachiyaT.SakakibaraH. (2017). Interactions between nitrate and ammonium in their uptake, allocation, assimilation, and signaling in plants. J. Exp. Bot.68, 2501–2512. doi: 10.1093/jxb/erw449
23
HanM.-L.LvQ.-Y.ZhangJ.WangT.ZhangC.-X.TanR.-J.et al. (2022). Decreasing nitrogen assimilation under drought stress by suppressing DST-mediated activation of nitrate reductase 1.2 in rice. Mol. Plant15, 167–178. doi: 10.1016/j.molp.2021.09.005
24
HossainM. A.LeeY.ChoJ.AhnC. H.LeeS. K.JeonJ. S.et al. (2010). The bZIP transcription factor OsABF1 is an ABA responsive element binding factor that enhances abiotic stress signaling in rice. Plant Mol. Biol.72, 557–566. doi: 10.1007/s11103-009-9592-9
25
HsiehP. H.KanC. C.WuH. Y.YangH. C.HsiehM. H. (2018). Early molecular events associated with nitrogen deficiency in rice seedling roots. Sci. Rep.8, 1–23. doi: 10.1038/s41598-018-30632-1
26
HuL.ChenW.YangW.LiX.ZhangC.ZhangX.et al. (2021). OsSPL9 regulates grain number and grain yield in rice. Front. Plant Sci.12. doi: 10.3389/fpls.2021.682018
27
Huynh-ThuV. A.IrrthumA.WehenkelL.GeurtsP. (2010). Inferring regulatory networks from expression data using tree-based methods. PLoS One5, 1–10. doi: 10.1371/journal.pone.0012776
28
HuY.ZhaoL.ChongK.WangT. (2008). Overexpression of OsERF1, a novel rice ERF gene, up-regulates ethylene-responsive genes expression besides affects growth and development in arabidopsis. J. Plant Physiol.165, 1717–1725. doi: 10.1016/j.jplph.2007.12.006
29
JinJ.TianF.YangD. C.MengY. Q.KongL.LuoJ.et al. (2017). PlantTFDB 4.0: Toward a central hub for transcription factors and regulatory interactions in plants. Nucleic Acids Res.45, D1040–D1045. doi: 10.1093/nar/gkw982
30
JoshiR.WaniS. H.SinghB.BohraA.DarZ. A.LoneA. A.et al. (2016). Transcription factors and plants response to drought stress: Current understanding and future directions. Front. Plant Sci.7. doi: 10.3389/fpls.2016.01029
31
JungH.ChungP. J.ParkS. H.RedillasM. C. F. R.KimY. S.SuhJ. W.et al. (2017). Overexpression of OsERF48 causes regulation of OsCML16, a calmodulin-like protein gene that enhances root growth and drought tolerance. Plant Biotechnol. J.15, 1295–1308. doi: 10.1111/pbi.12716
32
KatariM. S.NowickiS. D.AceitunoF. F.NeroD.KelferJ.ThompsonL. P.et al. (2010). VirtualPlant: A software platform to support systems biology research. Plant Physiol.152, 500–515. doi: 10.1104/pp.109.147025
33
KeelerB. L.GourevitchJ. D.PolaskyS.IsbellF.TessumC. W.HillJ. D.et al. (2016). The social costs of nitrogen. Sci. Adv.2(10):e1600219. doi: 10.1126/sciadv.1600219
34
LangfelderP.HorvathS. (2008). WGCNA: An r package for weighted correlation network analysis. BMC Bioinf.9, 559. doi: 10.1186/1471-2105-9-559
35
LeeS. C.KimS. H.KimS. R. (2013). Drought inducible OsDhn1 promoter is activated by OsDREB1A and OsDREB1D. J. Plant Biol.56, 115–121. doi: 10.1007/s12374-012-0377-3
36
LeeS.MarmagneA.ParkJ.FabienC.YimY.KimS.et al. (2020). Concurrent activation of OsAMT1;2 and OsGOGAT1 in rice leads to enhanced nitrogen use efficiency under nitrogen limitation. Plant J.103, 7–20. doi: 10.1111/tpj.14794
37
LiaoY.SmythG. K.ShiW. (2014). FeatureCounts: An efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics30, 923–930. doi: 10.1093/bioinformatics/btt656
38
LiH.HuB.ChuC. (2017). Nitrogen use efficiency in crops: Lessons from arabidopsis and rice. J. Exp. Bot.68, 2477–2488. doi: 10.1093/jxb/erx101
39
LiC.ShenH.WangT.WangX. (2015). ABA regulates subcellular redistribution of OsABI-LIKE2, a negative regulator in ABA signaling, to control root architecture and drought resistance in oryza sativa. Plant Cell Physiol.56, 2396–2408. doi: 10.1093/pcp/pcv154
40
LiuK.-H.HuangC.-Y.TsayY.-F. (1999). CHL1 is a dual-affinity nitrate transporter of arabidopsis involved in multiple phases of nitrate uptake. Plant Cell11, 865–874. doi: 10.1105/tpc.11.5.865
41
LouD.ChenZ.YuD.YangX. (2020). SAPK2 contributes to rice yield by modulating nitrogen metabolic processes under reproductive stage drought stress. Rice13, 35. doi: 10.1186/s12284-020-00395-3
42
LoveM. I.HuberW.AndersS. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol.15, 1–21. doi: 10.1186/s13059-014-0550-8
43
Lozano-JusteJ.LeonJ. (2010). Enhanced abscisic acid-mediated responses in nia1nia2noa1-2 triple mutant impaired in NIA/NR-and AtNOA1-dependent nitric oxide biosynthesis in arabidopsis. Plant Physiol.152, 891–903. doi: 10.1104/pp.109.148023
44
LuG.GaoC.ZhengX.HanB. (2009). Identification of OsbZIP72 as a positive regulator of ABA response and drought tolerance in rice. Planta229, 605–615. doi: 10.1007/s00425-008-0857-3
45
MasudaH.AungM. S.KobayashiT.HamadaT.NishizawaN. K. (2019). Enhancement of iron acquisition in rice by the mugineic acid synthase gene with ferric iron reductase gene and OsIRO2 confers tolerance in submerged and nonsubmerged calcareous soils. Front. Plant Sci.10. doi: 10.3389/fpls.2019.01179
46
MatsudaS.FunabikiA.FurukawaK.KomoriN.KoikeM.TokujiY.et al. (2012). Genome-wide analysis and expression profiling of half-size ABC protein subgroup G in rice in response to abiotic stress and phytohormone treatments. Mol. Genet. Genomics287, 819–835. doi: 10.1007/s00438-012-0719-3
47
NdiayeM.AdamM.GanyoK. K.GuisséA.CisséN.MullerB. (2019). Genotype-environment interaction: Trade-offs between the agronomic performance and stability of dual-purpose sorghum (Sorghum bicolor l. moench) genotypes in Senegal. Agronomy9, 867. doi: 10.3390/agronomy9120867
48
ObertelloM.ShrivastavaS.KatariM. S.CoruzziG. M. (2015). Cross-species network analysis uncovers conserved nitrogen-regulated network modules in rice. Plant Physiol.168, 1830–1843. doi: 10.1104/pp.114.255877
49
OgoY.KobayashiT.ItaiR. N.NakanishiH.KakeiY.TakahashiM.et al. (2008). A novel NAC transcription factor, IDEF2, that recognizes the iron deficiency-responsive element 2 regulates the genes involved in iron homeostasis in plants. J. Biol. Chem.283, 13407–13417. doi: 10.1074/jbc.M708732200
50
OgoY.Nakanishi ItaiR.NakanishiH.KobayashiT.TakahashiM.MoriS.et al. (2007). The rice bHLH protein OsIRO2 is an essential regulator of the genes involved in fe uptake under fe-deficient conditions. Plant J.51, 366–377. doi: 10.1111/j.1365-313X.2007.03149.x
51
OliverosJ. C. (2015) Venny. an interactive tool for comparing lists with venn’s diagrams. Available at: https://bioinfogp.cnb.csic.es/tools/venny/index.html.
52
ParkS. H.JeongJ. S.LeeK. H.KimY. S.Do ChoiY.KimJ. K. (2015). OsbZIP23 and OsbZIP45, members of the rice basic leucine zipper transcription factor family, are involved in drought tolerance. Plant Biotechnol. Rep.9, 89–96. doi: 10.1007/s11816-015-0346-7
53
PlettD. C.RanathungeK.MelinoV. J.KuyaN.UgaY.KronzuckerH. J. (2020). The intersection of nitrogen nutrition and water use in plants: New paths toward improved crop productivity. J. Exp. Bot.71, 4452–4468. doi: 10.1093/jxb/eraa049
54
RaudvereU.KolbergL.KuzminI.ArakT.AdlerP.PetersonH.et al. (2019). G:Profiler: A web server for functional enrichment analysis and conversions of gene lists (2019 Update). Nucleic Acids Res.47, W191–W198. doi: 10.1093/nar/gkz369
55
RoyChoudhuryA.GuptaB.SenguptaD. N. (2008). Trans-acting factor designated OSBZ8 interacts with both typical abscisic acid responsive elements as well as abscisic acid responsive element-like sequences in the vegetative tissues of indica rice cultivars. Plant Cell Rep.27, 779–794. doi: 10.1007/s00299-007-0498-1
56
SakuraiJ.IshikawaF.YamaguchiT.UemuraM.MaeshimaM. (2005). Identification of 33 rice aquaporin genes and analysis of their expression and function. Plant Cell Physiol.46, 1568–1577. doi: 10.1093/pcp/pci172
57
Santosh KumarV. V.YadavS. K.VermaR. K.ShrivastavaS.GhimireO.PushkarS.et al. (2021a). The abscisic acid receptor OsPYL6 confers drought tolerance to indica rice through dehydration avoidance and tolerance mechanisms. J. Exp. Bot.72, 1411–1431. doi: 10.1093/jxb/eraa509
58
Santosh KumarV. V.YadavS. K.VermaR. K.ShrivastavaS.GhimireO.PushkarS.et al. (2021b). The abscisic acid receptor OsPYL6 confers drought tolerance to indica rice through dehydration avoidance and tolerance mechanisms. J. Exp. Bot.72, 1411–1431. doi: 10.1093/jxb/eraa509
59
SasakawaH.YamamotoY. (1978). Comparison of the uptake of nitrate and ammonium by rice seedlings. Plant Physiol.62, 665–669. doi: 10.1104/pp.62.4.665
60
SevanthiA. M.SinhaS. K.RaniM.SainiM. R.KumariS.et al. (2021). Integration of dual stress transcriptomes and major QTLs from a pair of genotypes contrasting for drought and chronic nitrogen starvation identifies key stress responsive genes in rice. Rice14(1):49. doi: 10.1186/s12284-021-00487-8
61
ShannonP.MarkielA.OzierO.BaligaN.WangJ.RamageD.et al. (1971). Cytoscape: A software environment for integrated models. Genome Res.13, 426. doi: 10.1101/gr.1239303.metabolite. 山本隆久.
62
ShimJ. S.OhN.ChungP. J.KimY. S.ChoiY.KimJ. K. (2018). Overexpression of OsNAC14 improves drought tolerance in rice. Front. Plant Sci.9. doi: 10.3389/fpls.2018.00310
63
SmithT.HegerA.SudberyI. (2017). UMI-tools: modeling sequencing errors in unique molecular identifiers to improve quantification accuracy. Genome Res.27, 491–499. doi: 10.1101/gr.209601.116
64
SongS.WangG.WuH.FanX.LiangL.ZhaoH.et al. (2020). OsMFT2 is involved in the regulation of ABA signaling-mediated seed germination through interacting with OsbZIP23/66/72 in rice. Plant J.103, 532–546. doi: 10.1111/tpj.14748
65
SubudhiP. K.GarciaR. S.CoronejoS.TapiaR. (2020). Comparative transcriptomics of rice genotypes with contrasting responses to nitrogen stress reveals genes influencing nitrogen uptake through the regulation of root architecture. Int. J. Mol. Sci.21, 1–23. doi: 10.3390/ijms21165759
66
SwiftJ.AdameM.TranchinaD.HenryA.CoruzziG. M. (2019). Water impacts nutrient dose responses genome-wide to affect crop production. Nat. Commun.10(1):1374. doi: 10.1038/s41467-019-09287-7
67
UedaY.OhtsukiN.KadotaK.TezukaA.NaganoA. J.KadowakiT.et al. (2020). Gene regulatory network and its constituent transcription factors that control nitrogen-deficiency responses in rice. New Phytol.227, 1434–1452. doi: 10.1111/nph.16627
68
VaralaK.Marshall-ColónA.CirroneJ.BrooksM. D.PasquinoA. V.LéranS.et al. (2018). Temporal transcriptional logic of dynamic regulatory networks underlying nitrogen signaling and use in plants. Proc. Natl. Acad. Sci. United States America115, 6494–6499. doi: 10.1073/pnas.1721487115
69
VolanteA.DesiderioF.TondelliA.PerriniR.OrasenG.BiselliC.et al. (2017). Genome-wide analysis of japonica rice performance under limited water and permanent flooding conditions. Front. Plant Sci.8. doi: 10.3389/fpls.2017.01862
70
WangQ.ZengX.SongQ.SunY.FengY.LaiY. (2020). Identification of key genes and modules in response to cadmium stress in different rice varieties and stem nodes by weighted gene co-expression network analysis. Sci. Rep.10, 1–13. doi: 10.1038/s41598-020-66132-4
71
WeihM. (2014). A calculation tool for analyzing nitrogen use efficiency in annual and perennial crops. Agronomy4, 470–477. doi: 10.3390/agronomy4040470
72
WengXJuengerTE. A High-throughput 3'-tag RNA sequencing for large-scale time-series transcriptome studies. Methods Mol Biol (2022) 2398:151–72. doi: 10.1007/978-1-0716-1912-4_13
73
WilliamsonJ. M. (2011). The role of information and prices in the nitrogen fertilizer management decision: New evidence from the agricultural resource management survey. J. Agric. Resource Economics36, 552–572. doi: 10.22004/ag.econ.119180
74
WuT. M.LinK. C.LiauW. S.ChaoY. Y.YangL. H.ChenS. Y.et al. (2016). A set of GFP-based organelle marker lines combined with DsRed-based gateway vectors for subcellular localization study in rice (Oryza sativa l.). Plant Mol. Biol.90, 107–115. doi: 10.1007/s11103-015-0397-8
75
XiangY.TangN.DuH.YeH.XiongL. (2008). Characterization of OsbZIP23 as a key player of the basic leucine zipper transcription factor family for conferring abscisic acid sensitivity and salinity and drought tolerance in rice. Plant Physiol.148, 1938–1952. doi: 10.1104/pp.108.128199
76
XuF.WangK.YuanW.XuW.LiuS.KronzuckerH. J.et al. (2019). Overexpression of rice aquaporin OsPIP1;2 improves yield by enhancing mesophyll CO2 conductance and phloem sucrose transport. J. Exp. Bot.70, 671–681. doi: 10.1093/jxb/ery386
77
YangH. C.KanC. C.HungT. H.HsiehP. H.WangS. Y.HsiehW. Y.et al. (2017). Identification of early ammonium nitrate-responsive genes in rice roots. Sci. Rep.7, 1–16. doi: 10.1038/s41598-017-17173-9
78
YuJ.XuanW.TianY.FanL.SunJ.TangW.et al. (2021). Enhanced OsNLP4-OsNiR cascade confers nitrogen use efficiency by promoting tiller number in rice. Plant Biotechnol. J.19, 167–176. doi: 10.1111/pbi.13450
79
ZhangS.HaiderI.KohlenW.JiangL.BouwmeesterH.MeijerA. H.et al. (2012). Function of the HD-zip I gene Oshox22 in ABA-mediated drought and salt tolerances in rice. Plant Mol. Biol.80, 571–585. doi: 10.1007/s11103-012-9967-1
80
ZhangC.LiC.LiuJ.LvY.YuC.LiH.et al. (2017). The OsABF1 transcription factor improves drought tolerance by activating the transcription of COR413-TM1 in rice. J. Exp. Bot.68, 4695–4707. doi: 10.1093/jxb/erx260
81
ZhaoH.WuD.KongF.LinK.ZhangH.LiG. (2017). The arabidopsis thaliana nuclear factor y transcription factors. Front. Plant Sci.7. doi: 10.3389/fpls.2016.02045
82
ZongW.TangN.YangJ.PengL.MaS.XuY.et al. (2016). Feedback regulation of ABA signaling and biosynthesis by a bZIP transcription factor targets drought-resistance-related genes. Plant Physiol.171, 2810–2825. doi: 10.1104/pp.16.00469
Summary
Keywords
rice, drought, nitrogen, gene regulatory network, network validation, NUE, GENIE3, WGCNA
Citation
Shanks CM, Huang J, Cheng C-Y, Shih H-JS, Brooks MD, Alvarez JM, Araus V, Swift J, Henry A and Coruzzi GM (2022) Validation of a high-confidence regulatory network for gene-to-NUE phenotype in field-grown rice. Front. Plant Sci. 13:1006044. doi: 10.3389/fpls.2022.1006044
Received
28 July 2022
Accepted
01 November 2022
Published
25 November 2022
Volume
13 - 2022
Edited by
Nandula Raghuram, Guru Gobind Singh Indraprastha University, India
Reviewed by
Antonio Lupini, Mediterranea University of Reggio Calabria, Italy; Fabien Chardon, INRA UMR1318 Institut Jean Pierre Bourgin, France
Updates
Copyright
© 2022 Shanks, Huang, Cheng, Shih, Brooks, Alvarez, Araus, Swift, Henry and Coruzzi.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Gloria M. Coruzzi, gloria.coruzzi@nyu.edu
This article was submitted to Plant Physiology, a section of the journal Frontiers in Plant Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.