Genome-Wide Identiﬁcation and Expression Analysis of GATA Gene Family under Different Nitrogen Levels in Arachis hypogaea L.

: Nitrogen, one of the essential elements, is a key determinant for improving peanut growth and yield. GATA zinc ﬁnger transcription factors have been found to be involved in regulation of nitrogen metabolism. However, a systematic characterization of the GATA gene family and patterns of their expression under different nitrogen levels remains elusive. In this study, a total of 45 GATA genes distributed among 17 chromosomes were identiﬁed in the peanut genome and classiﬁed into three subfamilies I, II and III with 26, 13 and 6 members, respectively, whose physicochemical characteristics, gene structures and conserved motifs were also analyzed. Furthermore, the optimal level of nitrogen fertilizer on the growth of peanut cultivar Yuhua 23 was determined by pod yield and value cost ratio from 2020 to 2022, and the results revealed that 150 kg hm − 2 nitrogen was the best for cultivation of peanut Yuhua 23 because of its highest pod yield and relatively higher VCR of more than four. In addition, expression patterns of peanut GATA genes under different nitrogen levels were detected by real-time quantitative PCR and several GATA genes were signiﬁcantly changed under a nitrogen level of 150 kg hm − 2 . Overall, the above results would be helpful for further understanding biological functions of the GATA gene family in cultivated peanut.


Introduction
GATA transcription factors are an important category of zinc-finger DNA binding proteins in plants, animals, and fungi, which are widely related to multiple biological processes [1]. These proteins can bind to the DNA sequences GATA (A/T) GATA (A/G) and thus be known as GATA transcription factors [2]. GATA transcription factors were originally found in fungi and animals and are often composed of multiple gene family members. Most GATA proteins contain either one or two zinc-finger domains that share the common sequence CX 2 CX 17-18 CX 2 C [2]. GATA transcription factors in animals usually possess two CX 2 CX 17 CX 2 C zinc finger domains, and the C-terminal one could bind to DNA to regulate expression of downstream genes [3]. Most fungal GATA factors have one single CX 2 CX 17 CX 2 C or CX 2 CX 18 CX 2 C domain that is highly similar to zinc-finger domains at the C end of animal GATA factors [4,5]. It was shown that 30 GATA transcription factors were previously found in the Arabidopsis thaliana genome, and these were classified into subfamilies I, II, III and IV [6]. Most plant GATA factors have a single CX 2 CX 18 CX 2 C domain, but some of them have a zinc finger ring with 20 or more amino acid residues in the zinc finger domain [6].
GATA transcription factors have been extensively studied in many animals and fungi. GATA transcription factors in animals play important roles in cell proliferation, differentiation, and development [3]. Fungal GATA transcription factors were involved in regulating light induction, nitrogen metabolism, mating-type switching and lateral body biosynthesis [5]. In Yarrowia lipolytica, nitrogen and lipid metabolism were regulated by GATA transcription factors gzf3 and gzf2 [7]. Moreover, numerous studies have shown that plant GATA transcription factors are related to nitrogen metabolism. The first reported plant GATA protein NTL1 was cloned from tobacco and had been demonstrated as a regulator of nitrogen metabolism [8]. In Oryza sativa, cytokinin-responsive GATA protein Cga1was shown to be involved in plant architecture and chloroplast development under different environmental conditions [9]. Inhibition of the transcription factor GATA4's expression in Arabidopsis thaliana improved tolerance response to nitrogen deficiency [10]. GATA transcription factors in plants also play critical roles in response to different stresses, flowering, and development, and regulating plant hormone signaling. In Brassica napus, a total of 96 GATA family members were identified and displayed different expression patterns under ABA, salinity, drought and cold stresses [11]. Sixty-four putative GATA factors were detected in the soybean genome and exhibited expression diversity under low nitrogen stress [12]. Transcription factor GATA2 is a key player mediating signaling crosstalk between brassinosteroid and light pathways in Arabidopsis thaliana [13,14]. Furthermore, in Arabidopsis thaliana, GATA transcription factors HANABA TARANU, ZIM and Blue Micropylar End 3 have been revealed to regulate shoot apical meristem development, flower development and seed germination, respectively [15][16][17]. In rice, the GATA transcription factor NECK LEAF 1 modulates organogenesis through regulating several pathways during reproductive development [18]. In Triticum aestivum L., over-expression of a GATA-like transcription factor TaZIM-A1 delayed flowering and resulted in a reduction of thousand-kernel weight [19]. More importantly, conserved GATA motifs have been identified in the promoter regions of several nitrogen metabolism-related genes such as nitrite reductase, nitrate reductase and glutamine synthetase [20][21][22].
Cultivated peanuts (Arachis hypogaea L.), one of the most essential oil crops, are widely cultivated in the tropical and subtropical areas of the world and provide edible oils and proteins for humans [23,24]. Nitrogen is one of the necessary elements for crop growth and development [25]. Nitrogen limitation has a negative impact on peanut growth, accumulation of dry matter and grain yield [26,27]. However, until now, members of the GATA gene family in cultivated peanut have not been systematically characterized and their expression patterns under different nitrogen levels have not yet been studied. In this study, a total of 45 GATA family genes were identified in the cultivated peanut genome and a comprehensive analysis of them including physical characteristics, gene structures, conserved motifs, chromosomal distribution, and phylogenetic classification was performed. Furthermore, expression patterns of the identified GATAs in peanut under different nitrogen levels were also investigated. This study will be valuable for further elucidating the functions of the GATA gene family under different nitrogen applications in cultivated peanut.

Identification of GATA Genes in the Arachis hypogaea L. Genome
The annotation files of Arachis hypogaea cv. Tifrunner genome were downloaded from the PeanutBase database (https://peanutbase.org/data/v2/Arachis/hypogaea/ annotations/Tifrunner.gnm1.ann1.CCJH/, accessed on 5 May 2022). To identify all genes of the GATA family in Arachis hypogaea L., the GATA zinc finger domain hidden Markov model (HMM) profile (PF00320) was obtained from the Pfam database in the website of the InterPro database (https://www.ebi.ac.uk/interpro/entry/pfam/, accessed on 5 May 2022) [28,29] and the HMMER program (Version 3.3.2) was used for identifying GATA genes among all the annotated protein sequences in the peanut genome at an E-value smaller than 1e −12 [30]. If one gene had two or more protein transcripts (isoforms), the longest one was used for further analysis. Moreover, the identified GATA genes were further confirmed with the Simple Modular Architecture Research Tool (https://smart.embl.de, accessed on 6 June 2022) [31]. In addition, number of amino acids, molecular weight, instability and aliphatic index, theoretical isoelectric point (pI) and grand average of hydro-pathicity (GRAVY) of all identified GATA proteins were analyzed with the ProtParam tool (https://web.expasy.org/protparam/, accessed on 6 June 2022).

Phylogenetic Classification of GATA Genes in Arachis hypogaea L.
Based on a previous study, the evolutionary relationship of all the known GATA genes from Arabidopsis thaliana was referenced and employed to classify the members of the GATA family in Arachis hypogaea L. [6]. In brief, multiple sequence alignment of all known GATA members from Arabidopsis thaliana and the identified GATA proteins in Arachis hypogaea L. were performed with the MUSCLE program with default settings [32]. Subsequently, a phylogenetic tree was constructed by MEGA 7.0 using the Neighbor-Joining (NJ) method with a bootstrap value of 1000 replicates [33].

Analysis of GATA Gene Structures and Conserved Motifs
Gene structures of the identified GATA members in Arachis hypogaea L. were visualized by using the online Gene Structure Display Server (http://gsds.gao-lab.org/, accessed on 6 July 2022) according to genome annotation [34]. The conserved motifs of the identified GATA members in Arachis hypogaea L. were identified and visualized by online Multiple EM for Motif Elicitation (MEME) (Version 5.4.1) (https://meme-suite.org/meme/tools/meme, accessed on 6 July 2022) with the following parameters: motif counts =10 and motif width between 6 and 50 [35]. DNA binding sites in the upstream 2000 bp regions of the promoter for GATA factors were identified by the online PlantPAN3.0 (http://plantpan.itps.ncku. edu.tw/promoter.php, accessed on 5 January 2023) [36].

Chromosomal Distribution of GATAFamilyGenes in Arachis hypogaea L.
The distribution of identified GATA members in Arachis hypogaea L. genome was determined based on genome annotation data and visualized with TBtoolsv1.108 software [37].

Determination of the Optimal Nitrogen Fertilizer Level for Peanut Growth
The field trials were conducted in the Modern Agricultural Demonstration Garden of Changyuan Branch of Henan Academy of Agricultural Sciences in Changyuan, Henan, China (114 • 38 E, 35 • 08 N) from 2020 to 2022. The soil parameters in 2020 were as follows: organic matter with 12.8 g/kg, available nitrogen with 108.96 mg/kg, available phosphorus with 14.6 mg/kg and available potassium with 135.4 mg/kg. Urea was selected as nitrogen fertilizer. Seeds of peanut cultivars Yuhua 23 were grown under five nitrogen levels (N0: 0kg hm −2 , N1: 75 kg hm −2 , N2: 150 kg hm −2 , N3: 225 kg hm −2 and N4: 300 kg hm −2 ) in 30 plots and each plot was 5.6 m × 2.4 m. Seeds were planted 12 cm apart and the distance between rows was 40 cm. Total pods per plant (TPP) were counted manually. Pods were air-dried outdoors until water content was below 9% and economic pods per plant (EPP), weight of pods per plant (PWP) and hundred-pod weight (HPW) were determined after harvesting. Pod dry weight per harvested plot was then weighted and pod yield was calculated as kg per hm 2 . The value cost ratio (VCR) was computed following the formula: VCR = (Yt − Yc) × P/C, where VCR denotes the value cost ratio of treatment, Yt and Ycare pod yield of treatment and control, respectively, P is price of peanut pod per kg and C is cost of fertilizer per hectare of treatment [38]. Total RNA was extracted from peanut leaves collected at 50 days after sowing by using a Plant RNAprep Pure Kit (TIANGEN, Beijing, China) with three biological replicates of each sample. FastQuant RT kit (TIANGEN, Beijing, China) was used to synthesizethe first-strand cDNA. A SuperReal PreMix Plus (SYBR Green) kit (TIANGEN, Beijing, China) was used for qRT-PCR amplification. Relative expression levels were determined using the 2 −∆∆Ct method and the Ahactin gene was used as the reference gene [39]. Primer sequences used in this study are listed in Table S1.

Genome-Wide Identification and Sequence Characteristics of GATA Gene Family in Arachis hypogaea L.
A total of 45 GATA genes were identified in the Arachis hypogaea L. genome ( Table 1)

Sequence Alignment and Phylogenetic Analysis of Peanut GATA Genes
To elucidate the phylogenetic relationships of GATA gene family members between A. hypogaea L. and A. thaliana, protein sequences of the identified GATAs were used for further building a neighbor-joining phylogenetic tree, as shown in Figure 1. According to the classification of the Arabidopsis GATA gene family, GATA members from A. hypogaea L. and A. thaliana cluster into four phylogenetic subfamilies. Of these, 14 A. thaliana GATA members and 26 A. hypogaea L. GATA proteins cluster to subfamily I, harboring the largest number of GATA in both A. thaliana and A. hypogaea L., 11 A. thaliana GATA members and 13 A. hypogaea L. GATA proteins to subfamily II, and three A. thaliana GATA members and six A. hypogaea L. GATA proteins to subfamily III. Interestingly, there were two A. thaliana GATA members in subfamily IV and it was found that no A. hypogaea L. GATA proteins cluster into subfamily IV. To further analyze the sequence features of 45 identified GATA members in A. hypogaea L., conserved domains of 45 GATA proteins were detected, and the result of protein sequence alignment demonstrated that all of them have only one conserved domain with 18-20 residues among zinc finger loop (C-X 2 -C-X 18-20 -CNAC). All 39 A. hypogaea L. GATA members in subfamilies I and II harbor the C-X 2 -C-X 18 -CNAC conserved domain but the other six members in subfamily III have the C-X 2 -C-X 20 -CNAC rather than C-X 2 -C-X 18 -CNAC domain ( Figure 2). Especially, there is a conserved amino acid motif TPQWRXGPXGXKTL between the second and third cysteine residues in the C-X 2 -C-X 18 -CNAC zinc finger loop of subfamily I and a conserved amino acid motif TX 2 TPLWRXGPXGPKXL between the second and third cysteine residues in the C-X 2 -C-X 18 -CNAC zinc finger loop of subfamily II. In subfamily III, there is also a conserved amino acid motif GX 3 KXTPXMRRGPXGPRXL between the second and third cysteine residues in the C-X 2 -C-X 20 -CNAC zinc finger loop ( Figure 2).

Analysis of Gene Structures and Conserved Motifs of Peanut GATA Family
A detailed illustration of the 45 identified GATA family gene exon-intron structures was made, and the results demonstrate that the numbers of exons of the identified peanut GATA family members vary from 1 (69GBM5.1) to 11 (2L20XB.1 and 86LJVS.1) and GATA genes in each subfamily display similar exon-intron structures (Figure 3). The most GATA members in the subfamily I and II possess two or three exons except for three genes (NH2FYW.1, FRWT0V.1 and V4I8CJ.1) with five exons, one gene (36H91F.1) with four exons and one gene (69GBM5.1) with only one exon. In the subfamily III, all six GATA genes comprise more than seven exons with an average of 8.67 exons per gene.

Analysis of Gene Structures and Conserved Motifs of Peanut GATA Family
A detailed illustration of the 45 identified GATA family gene exon-intron structures was made, and the results demonstrate that the numbers of exons of the identified peanut GATA family members vary from 1 (69GBM5.1) to 11 (2L20XB.1 and 86LJVS.1) and GATA genes in each subfamily display similar exon-intron structures (Figure 3). The most GATA members in the subfamily I and Ⅱ possess two or three exons except for three genes (NH2FYW.1, FRWT0V.1 and V4I8CJ.1) with five exons, one gene (36H91F.1) with four exons and one gene (69GBM5.1) with only one exon. In the subfamily Ⅲ, all six GATA genes comprise more than seven exons with an average of 8.67 exons per gene.

Analysis of Chromosomal Distribution and Gene Duplication of Peanut GATA Genes
As shown in Figure 5

Analysis of Chromosomal Distribution and Gene Duplication of Peanut GATA Genes
As shown in Figure 5   Note: + and − indicate the plus and minus strand of chromosome, respectively.

Optimal Level of Nitrogen Fertilizer on the Growth of Peanut Yuhua 23
The increased application of nitrogen fertilizer significantly improved yield components of peanut, and total pods per plant increased the most, followed by economic pods per plant and hundred-pod weight (Table 3). With increasing level of nitrogen application, the proportion of economic pods to total pods decreased. The increased rates of pod yield in 2020 and 2022 were almost consistent, which were 18.6-19.0%, 33.5-34.6%, 35.9-37.6% and 35.2-37.0% at nitrogen levels of 75 kg hm −2 , 150 kg hm −2 , 225 kg hm −2 and 300 kg hm −2 , respectively. In 2020 and 2022, VCR values were higher than 4 at nitrogen levels of 75 kg hm −2 and 150 kg hm −2 and pod yield increased by 533.6 kg hm −2 and 586.2 kg hm −2 , respectively. Pod yield increased by only about 3% at nitrogen levels of 225 kg hm −2 and 300 kg hm −2 compared with the 150 kg hm −2 and the VCR quickly decreased below 2,

Optimal Level of Nitrogen Fertilizer on the Growth of Peanut Yuhua 23
The increased application of nitrogen fertilizer significantly improved yield components of peanut, and total pods per plant increased the most, followed by economic pods per plant and hundred-pod weight (Table 3). With increasing level of nitrogen application, the proportion of economic pods to total pods decreased. The increased rates of pod yield in 2020 and 2022 were almost consistent, which were 18.6-19.0%, 33.5-34.6%, 35.9-37.6% and 35.2-37.0% at nitrogen levels of 75 kg hm −2 , 150 kg hm −2 , 225 kg hm −2 and 300 kg hm −2 , respectively. In 2020 and 2022, VCR values were higher than 4 at nitrogen levels of 75 kg hm −2 and 150 kg hm −2 and pod yield increased by 533.6 kg hm −2 and 586.2 kg hm −2 , respectively. Pod yield increased by only about 3% at nitrogen levels of 225 kg hm −2 and 300 kg hm −2 compared with the 150 kg hm −2 and the VCR quickly decreased below 2, which indicated that it was not economical to increase nitrogen fertilizer above 150 kg hm −2 . However, VCR was still higher than 4 at 225 kg hm −2 nitrogen application in 2021 but the pod yield was lower than that in 2020 and 2022. So, the 150 kg hm −2 nitrogen application could produce relatively higher peanut yield with VCR more than 4.

Expression Analysis of Peanut GATA Gene Family under Different Nitrogen Levels
To investigate the role of the GATA gene family in improving cultivated peanut growth, expression changes of 42 peanut GATA genes in the early stage of pod development were analyzed under a nitrogen level of 150 kg hm −2 ( Figure 6). In detail, 12 members of subfamily I were significantly up-regulated and the top three of them (Z6CAHI.1, LJYJ4M.1 and MTW1LM.1) were remarkably increased 309.14-, 67.13-and 75.63-fold compared with the control, respectively. Only two members (E1UHJ9.1 and HHX0PC.1) were significantly down-regulated and decreased 3.3-and 3.84-fold, respectively. In the subfamily II, five genes were significantly up-regulated, especially, the member TF5DD7.1, which was 47 times higher than the control and one gene considerably decreased 2.63-fold compared to the control. Four of five detected subfamily III genes were significantly increased with the exception of 0U49HS.1 gene; in particular, the 2L20XB.1 gene had the highest expression and increased 13.82-fold compared to the control. (A-C) indicate expression patterns of peanut subfamilies Ⅰ, Ⅱ and Ⅲ GATA genes, respectively. Asterisks indicate significant differences (Student's t-test). ** and * indicated significant differences with p-values smaller than 0.01 and 0.05, respectively.

Discussion
The GATA gene family has been reported to play important roles in cell proliferation and development, nitrogen and lipid metabolism, and biotic and abiotic stresses [1,3,7]. Although GATA family members have been identified in several plants including Eucalyptus urophylla [40], Cucumis sativus L. [41], Brassica napus [11], Triticum aestivum L. [42] and Glycine max [12], they have not been systematically identified in cultivated peanuts until now. In this study, a total of 45 GATA family members were found in the Arachis hypogaea L. genome and classified into subfamilies I, Ⅱ and Ⅲ with 26, 13 and 6 members, respectively, which is consistent with GATA subfamily I containing the most members, followed by subfamilies Ⅱ and Ⅲ in Salvia miltiorrhiza [43], Eucalyptus urophylla [40], Cucumis sativus L. [41] and Brassica napus [11]. However, there was no GATA gene of subfamily IV found in the Arachis hypogaea L. genome, which is not consistent with the distribution of GATA families in several previously reported dicotyledons, such as Arabidopsis thaliana [6], Brassica napus [11] and Glycine max [12] and is the first report of GATA genes (A-C) indicate expression patterns of peanut subfamilies I, II and III GATA genes, respectively. Asterisks indicate significant differences (Student's t-test). ** and * indicated significant differences with p-values smaller than 0.01 and 0.05, respectively.

Discussion
The GATA gene family has been reported to play important roles in cell proliferation and development, nitrogen and lipid metabolism, and biotic and abiotic stresses [1,3,7]. Although GATA family members have been identified in several plants including Eucalyptus urophylla [40], Cucumis sativus L. [41], Brassica napus [11], Triticum aestivum L. [42] and Glycine max [12], they have not been systematically identified in cultivated peanuts until now. In this study, a total of 45 GATA family members were found in the Arachis hypogaea L. genome and classified into subfamilies I, II and III with 26, 13 and 6 members, respectively, which is consistent with GATA subfamily I containing the most members, followed by subfamilies II and III in Salvia miltiorrhiza [43], Eucalyptus urophylla [40], Cucumis sativus L. [41] and Brassica napus [11]. However, there was no GATA gene of subfamily IV found in the Arachis hypogaea L. genome, which is not consistent with the distribution of GATA families in several previously reported dicotyledons, such as Arabidopsis thaliana [6], Brassica napus [11] and Glycine max [12] and is the first report of GATA genes of subfamily IV not existing in dicotyledons. Interestingly, no GATA members of subfamily IV were identified in monocots moso bamboo and rice [6,44]. Thus, GATA genes of subfamily IV are absent in both dicotyledons and monocots.
Differences in gene structures and conserved motifs among members of the GATA protein family may result in functional divergences. As for gene structure, most (34/39) GATA members of subfamilies I and II in cultivated peanut possess two or three exons, which is similar to subfamilies I and II in cucumber and rapeseed [11,41]. However, all six members of subfamily III have more than seven exons per gene, which was also found in subfamily III in cucumber; this is different from members of the subfamily III in Brassica napus and Salvia miltiorrhiza, some of which have less than five exons or even only one exon [11,41,45]. Furthermore, the conserved domain of all GATA members in subfamilies I and II in cultivated peanut was identified as C-X 2 -C-X 18 -CX 2 C and the conserved domain of six subfamily III proteins was C-X 2 -C-X 20 -C X 2 C, which is consistent with previously identified conserved structures in Ophiorrhiza pumila and Fagopyrum tataricum [46,47]. However, two GATA genes of subfamily II such as BnGATA2.8 and BnGATA2.26 in Brassica napus harbor the N-X 2 -C-X 18 -CX 2 C domain and one subfamily II GATA gene Csa4G286370 in Cucumis sativus L. possesses the C-X 4 -C-X 18 -C-X2-C domain instead of the C-X 2 -C-X 18 -CX 2 C conserved domain [11,41]. In addition, there were three conserved amino acid motifs TPQWRXGPXGXKTL, TX 2 TPLWRXGPXGPKXL and GX 3 KXTPXMRRGPXGPRXL found in the peanut GATA subfamilies I, II and III, respectively. The analysis of gene structures and conserved motifs demonstrates that GATA members of the same subfamily show relatively high conservation in different species and GATA genes between subfamilies have their own special characteristics.
The supply of nitrogen, one of the essential nutrient elements, is a key determinant for peanut growth and yield [48]. Pod yield of peanut cultivar Yuhua 23 increased by 48.8% and 108.6% with nitrogen level of 30 kg ha −1 and 60 kg ha −1 , respectively [49]. In this study, increased nitrogen application from 75 kg hm −2 to 300 kg hm −2 continued to increase peanut pod yield, while it increased by only about 3% under 225 kg hm −2 and 300 kg hm −2 , compared to a 150 kg hm −2 nitrogen level. Importantly, VCR values were higher than four at nitrogen levels of 75 kg hm −2 to 150 kg hm −2 from 2020 and 2022 but VCR quickly decreased below 2 at nitrogen levels of 225 kg hm −2 and 300 kg hm −2 . Although VCR was still higher than 4 at 225 kg hm −2 nitrogen level in 2021, peanut pod yield was not as high as in 2020 and 2022, possibly due to heavy rain during peanut growth stages. A VCR of 2 means 100% return on the cost of purchasing fertilizer [50] and VCR ≥ 4 was regarded as more appropriate to ensure fertilizer cost [51]. So, it is worth noticing that 150 kg hm −2 nitrogen application was the best treatment suitable for peanut Yuhua 23 growth because of its highest pod yield and relatively better VCR of more than four. GATA transcription factors have been shown to participate in nitrogen metabolism in fungi and plants [8,9,12,52,53]. For example, heterologous expression of a putative GATA gene DhGZF3 in Saccharomyces cerevisiae was found to regulate nitrogen metabolic genes [54]. Overexpression of poplar GATA transcription factor PdGNC in Arabidopsis thaliana had pronounced effects on growth rate, chloroplast ultrastructure and photosynthetic capacity under low nitrogen levels [53]. Two soybean GATA members GmGATA44 and GmGATA58 were found to potentially regulate nitrogen metabolism [12]. In our study, most members of the GATA family in cultivated peanut were up-regulated and very few members are downregulated under a 150 kg hm −2 nitrogen level. In particular, the expression level of several GATA genes including Z6CAHI.1, LJYJ4M.1, MTW1LM.1 and TF5DD7.1 was induced more than 40-fold compared to the control and thus it could be speculated that they are involved in nitrogen metabolism regulation because the promoter regions of eight nitrogen metabolism-related genes in the cultivated genome, including two nitrite reductase genes (X7K798 and AXQD4Y), two nitrate reductase genes (7L21K7 and B2KMHD) and four glutamine synthetase genes (0KDC5W, HA9E1Y, PJ2I0S and 9W4G6J), harbored several DNA binding sites for GATA transcription factors (Table S2).

Conclusions
Taken together, a systematic characterization of the GATA gene family and patterns of their expression under different nitrogen levels in Arachis hypogaea L. was performed and the results provide valuable information for further understanding functional differences in GATA transcription factors' response to nitrogen application in cultivated peanut.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/agronomy13010215/s1, Table S1. Primers used in this study; Table S2. DNA binding sites identified for GATA factors of 8 nitrogen metabolism-related genes in the cultivated genome by PlantPAN3.0.