Next Article in Journal
Alpha Lipoic Acid as a Protective Mediator for Regulating the Defensive Responses of Wheat Plants against Sodic Alkaline Stress: Physiological, Biochemical and Molecular Aspects
Next Article in Special Issue
Metal(loid)s in Common Medicinal Plants in a Uranium Mining-Impacted Area in Northwestern New Mexico, USA
Previous Article in Journal
Antioxidant and Antiproliferation Activities of Lemon Verbena (Aloysia citrodora): An In Vitro and In Vivo Study
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 11
Issue 6
10.3390/plants11060786
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Identification and Expression Profile Reveal Potential Roles of Peanut ZIP Family Genes in Zinc/Iron-Deficiency Tolerance

by
Zhen Zhang
,
Nannan Chen
,
Zheng Zhang
and
Gangrong Shi
*
College of Life Sciences, Huaibei Normal University, Huaibei 235000, China
*
Author to whom correspondence should be addressed.
Plants 2022, 11(6), 786; https://doi.org/10.3390/plants11060786
Submission received: 9 February 2022 / Revised: 11 March 2022 / Accepted: 14 March 2022 / Published: 16 March 2022
(This article belongs to the Special Issue Microelements in Plant and Soil)
Browse Figures
Versions Notes

Abstract

:
Zinc/iron-regulated transporter-like protein (ZIP) family genes play crucial roles in metal uptake and transport in plants. However, little is known about their functions in peanut. Here, genome-wide analysis identified 30 peanut AhZIP genes that were divided into four classes. Most AhZIPs experienced whole-genome or segmental duplication. AhZIP proteins harbored 3–8 transmembrane domains and a typical ZIP domain, showing considerable homology with BbZIP from Bordetella bronchiseptica. Clustered AhZIPs generally share similar gene/protein structures; however, unique features were found in AhIRT1.2, AhZIP1.2, AhZIP3.5 and AhZIP7.8. RNA-seq data revealed that AhZIP2.1/2.2, AhZIP4.1/4.2 and AhZIP11.1/11.2 were highly and preferentially expressed in roots, nodule and reproductive tissues. RT-qPCR analysis indicated that transcriptional responses of AhZIPs to Fe/Zn deficiency are cultivar dependent. The expressions of AhIRT1.1, AhIRT1.2 and AhZIP6.1 were closely related to Fe uptake and translocation. AhIRT1.1 and AhZIP7.2 expression were significantly correlated with Zn accumulation. The expression of AhIRT1.1, AhIRT1.2, AhZIP3.6, AhZIP6.1 and AhZIP11.1 was associated with Mn uptake and translocation. The results confirmed that AhZIP genes play crucial roles in the uptake and transport of Fe, Zn and Mn in peanut, providing clues to further functionally characterize AhZIP genes in the future.

1. Introduction

Iron (Fe) and zinc (Zn) are essential metal micronutrients for virtually all organisms. Both elements act as catalytic and structural cofactors for a large number of enzymes and play vital roles in many metabolic processes [1]. Thus, the shortage of Fe or Zn inhibits plant growth and development, posing major abiotic stresses in crop production. Most Fe and Zn in the human body are acquired from plant-based diets. However, plants are not a good source of these micronutrients because staple crops have low concentrations of Zn and Fe in edible tissues [2]. In fact, plants often suffer from Zn and Fe deficiencies due to the scarcity of Zn or low availability of Fe in the soil [3,4]. Therefore, billions of people worldwide suffer from deficiencies of these two elements, leading to nutritional disorders [5]. It is of great importance to enhance the content of Fe and Zn in edible parts of crops by improving their capacity for Fe and Zn acquirement.
The zinc/iron-regulated transporter-like protein (ZIP) family plays crucial roles in the uptake and transport of essential or nonessential divalent metals in plants, including Fe and Zn. AtIRT1, the first ZIP protein identified in Arabidopsis thaliana, is responsible for taking up Zn, Mn, Co, Ni and Cd from the rhizosphere to root cells [6,7,8,9]. A recent study indicates that AtIRT1 is also specifically expressed in phloem companion cells and has a role in Fe translocation in aboveground organs [10]. AtIRT2 localizes the vesicle membrane and compartmentalizes Fe into vesicles to prevent its toxicity in the cytosol [11]. AtIRT3 is a plasma membrane localized transporter involved in the uptake of Zn and Fe in Arabidopsis [12]. AtZIP1 is a vacuolar transporter that is responsible for remobilizing Mn and Zn from vacuole to cytoplasm in root cells [13]. AtZIP2 is localized to the plasma membrane and may mediate Mn/Zn uptake into parenchyma cells in the xylem, contributing to xylem loading and transport of Mn/Zn to the shoot [13]. AtZIP3 plays a role in the uptake of Zn and Fe from the soil to the plant roots, while AtZIP4 transports Zn intracellularly or between plant tissues [14,15].
In rice (Oryza sativa), several ZIP members have been functionally characterized. OsIRT1 and OsIRT2 are Fe transporters that directly absorb the Fe2+, Zn2+, Cu2+, and Cd2+ into root cells [16,17,18]. OsZIP1 resides in the endoplasmic reticulum and plasma membrane and functions as a metal exporter in rice under Zn, Cu or Cd excess conditions [19]. OsZIP4 is a Zn transporter that may be involved in the translocation of Zn within plants [20,21]. OsZIP5 and OsZIP8 serve as the plasma membrane-localized transporter involved in Zn uptake and distribution within rice [22,23]. OsZIP7 is a plasma membrane Zn-specific transporter that plays a key role in xylem loading in roots and inter-vascular transfer in nodes to deliver Zn/Cd to developing tissues and grain in rice [24,25,26].
In contrast to the significant progress of functional characterization for individual genes, the systematic genome-wide study of ZIP gene family is limited. Since the first conducted by Guerinot [15], the whole genome identification of ZIP family has been carried out only in a few plants, including maize (Zea mays) [27], potato (Solanum tuberosum) [28], poplar (Populus trichocarpa) [29,30] and trifoliate orange (Poncirus trifoliata) [31]. Based on sequence alignments, the eukaryotic ZIP family is split into four subfamilies (ZIPI, ZIPII, gufA and LIV-1) [32]. Most of them contain 309–476 amino acid residues with eight putative transmembrane domains. The majority of ZIP proteins share a similar membrane topology where the N- and C-terminal ends are extracytoplasmic [32]. A histidine-rich domain (HRD) is contained in the long variable region of the cytoplasmic loop between TM3 and TM4, which is considered the metal binding domain playing roles in metal transport. Amphipathic TM4 and TM5 form cavities through which metal ions pass [15,32].
Peanut (Arachis hypogaea L., 2n = 4x = 40) is the fourth major oil crop widely grown throughout the world. It provides 20% of edible oil and 11% of food protein for global people annually. Peanut is a rich source of micronutrients, including Zn, and thus makes it more important for Zn biofortification [33]. The concentration of Zn in peanut seeds ranged from 11 to 77 mg kg−1, with an average of 45 mg kg−1 [34]. Unfortunately, peanut productivity is always affected by Fe and/or Zn deficiencies in soil because a large proportion of the crop grows in calcareous soils [33,35]. To overcome the deficiencies of Fe and Zn and to enhance their concentrations in seeds, it is necessary to fully understand the mechanism of the uptake, distribution and translocation of Fe and Zn in peanut plants.
Although several transporter genes such as AhIRT1 [36] and AhNramp1 [37] have been identified in peanut, the molecular mechanisms underlying the metal homeostasis remain unknown. Recently, the whole genome sequences of the cultivated peanut (A. hypogaea cv. Tifrunner) as well as the two wild ancestral species, A. duranensis and A. ipaënsis, have been released [38,39]. These studies make it possible for identifying gene families at the whole genome level. Here, 30 genes of the ZIP family were identified from cultivated peanut, and their structure, function and evolution were characterized. Moreover, the expression of ZIP genes in response to Fe/Zn deficiencies was evaluated. Our findings would provide clues to further characterize the functions of ZIP proteins in the uptake and translocation of Fe and Zn in peanut plants.

2. Results

2.1. Identification and Phylogenetic Analysis of the AhZIP Family in Peanut

A total of 30 putative AhZIP genes were identified in peanut, including four AhIRT1, two AhZIP1, two AhZIP2, six AhZIP3, two AhZIP4, four AhZIP6, eight AhZIP7 and two AhZIP11 (Table 1). The length of AhZIP genes varied from 937 bp (AhZIP3.6) to 8962 bp (AhIRT1.1), with CDS lengths from 456 bp (AhIRT1.2) to 1227 bp (AhZIP4.1). The amino acid number of AhZIP proteins ranged from 151 (AhIRT1.2) to 408 bp (AhZIP4.1), and the molecular weight varied from 16.55 kDa (AhIRT1.2) to 43.82 kDa (AhZIP4.1). The instability, GRAVY and aliphatic index of AhZIP proteins ranged from 26.98 (AhZIP2.1) to 47.10 (AhIRT1.2), from 0.167 (AhZIP1.2) to 0.767 (AhZIP11.1) and from 86.03 (AhZIP1.2) to 117.14 (AhZIP3.6), respectively. The isoelectric point (pI) of all AhZIP proteins less than 7 ranged from 5.28 (AhIRT1.2) to 6.86 (AhIRT1.3) (Table 1). TMD numbers of AhZIP proteins ranged from 3 to 8, and most of them were predicted to be plasma membrane localized except AhZIP1.2, AhZIP3.2, AhZIP3.5 and AhZIP7.8, which were predicted to localize to the endomembrane system (Table 1).
To reveal the phylogenetic relationship among AhZIP genes, 69 ZIP protein sequences from peanut, Arabidopsis, rice and trifoliate orange were used to construct phylogeny with the NJ method. As reported in previous studies [29,30,31], ZIP members were divided into four classes: I, II, III and IV (Figure 1). Class I, the largest class, contained 20 AhZIP members. Class II consisted of AhZIP4.1 and AhZIP4.2. Class III was composed of AhZIP2.1/2.2 and AhZIP11.1/11.2. Class IV contained four orthologs of AhZIP6 (Figure 1).

2.2. Gene Structure, Duplication and Ka/Ks of the AhZIP Family

Exon–intron organizations revealed that AhZIP genes belonging to the same phylogenetic groups showed similar exon–intron organizations (Figure 2). Most AhZIP genes contained two exons with one intron. However, eight genes have distinct exon–intron structures, including AhZIP1.1 (3 exons and 2 introns), AhZIP3.1 (3 exons and 3 introns), AhZIP3.2 (2 exons and 2 introns), AhZIP4.1 (4 exons and 3 introns), AhZIP4.2 (2 exons and 4 introns), AhZIP11.1/11.2 (3 exons and 2 introns) and AhZIP7.8 (5 exons and 4 introns) (Figure 2).
The 30 AhZIP genes were located unevenly in 13 chromosomes. A total of 13 and 17 AhZIP genes were identified from subgenomes A (Chr. 01–10) and B (Chr. 11–20), respectively (Figure 3). Chr. 01, 04, 11, 14, 17 and 18 contained three AhZIP genes, Chr. 05, 07, 08 and 15 had two genes in each chromosome, and Chr. 06 and 13 contained only one gene each, while no AhZIP gene was identified in Chr. 02, 03, 09, 10, 12, 19 and 20. Most of the AhZIP genes experienced gene duplication events except AhIRT1.1/1.3, AhZIP7.2/7.8 and AhZIP3.5/3.6, resulting in 16 gene pairs (Figure 3). Among the duplicated genes, 12 collinear blocks resulted from whole-genome duplications (WGDs), and AhZIP6.1/6.3 and AhZIP6.2/6.4 resulted from segmental duplication. No tandem duplication was detected in the AhZIP genes. The Ka/Ks ratios (ratios of the number of nonsynonymous substitutions per nonsynonymous site/the number of synonymous substitutions per synonymous site) of all gene duplication pairs were less than 1 (Table 2), indicating that the AhZIP genes evolved under purifying selection [40].

2.3. Conserved Motifs, Domain Architectures and Models of AhZIP Proteins

Ten conserved motifs were identified in the sequences of AhZIP proteins; among them, motifs 1, 2, 3, 5 and 9 were annotated as zinc transporters, according to the Pfam tools (Figure 4a and Table S1). All AhZIP proteins shared motifs 1 and 2 except AhZIP1.2 (without motif 1) and AhZIP7.8 (without motif 2). The distribution pattern of conserved motifs varied among phylogenetic clades, whereas it was generally similar within the same phylogenetic clades. Five proteins, including AhIRT1.2, AhZIP1.1, AhZIP1.2, AhZIP7.8 and AhZIP3.5, were found to have different motif profiles from their orthologs (Figure 4a), indicating that these proteins might have distinct functions.
All AhZIP proteins contained only one domain named ZIP, which is the typical domain of the family (Figure 4b). All AhZIP proteins were well modelled with the template, 6pgi.1.A (Figure S1 and Table S2), which is the A chain of the BbZIP protein from Bordetella bronchiseptica. BbZIP has been revealed to have a binuclear metal center, where two metal ions were trapped halfway through the membrane and the two metal-binding sites play asymmetric roles within the transport pathway [41]. Sequence identity ranged from 13.27% to 20.16%, the value of GMQE ranged from 0.12 to 0.54 and QMEANDisCo global score ranged from 0.32 to 0.54 (Table S2). These data are suggestive of the high quality of the 3D model predictions on AhZIP proteins.
Multiple sequence alignment showed considerable homology between the BbZIP and AhZIP proteins throughout the TMD, particularly in TM2, TM4 and TM5 (Figure 5). The TMD structure of AhZIP proteins showed great diversity. Among the 30 AhZIP proteins, only 14 have the 3 + 5 TMD structure (Figure 5). AhIRT1.2, AhZIP1.2 and AhZIP3.5 only contained the first three TMDs (TM1–TM3), AhZIP7.8 contained the last four TMDs (TM5–TM8), and the remaining homologs of AhZIP3 and AhZIP7 contained the first five TMDs (TM1–TM5) (Figure 5). A long chain variable region was found between TM3 and TM4 in most AhZIPs, except for AhIRT1.2 and AhZIP7.8, and most of them contain various HRDs such as HXHXH, HHH, HHHHH, HXHXHXH and HHXHXHXH (Figure 5). Additionally, glycine (G) residues are always found near or inside HRDs (Figure 5).

2.4. Expression Profiles of AhZIP Genes in Different Tissues of Peanut

The RNA-seq data showed that all AhZIP genes were tissue-specifically expressed in peanut plants (Table S3). Generally, the gene expression profiles were similar within the same phylogenetic classes. The 30 AhZIP genes could be classified into three distinct groups according to the gene expression patterns (Figure 6). Group 1 consists of two orthologs of ZIP11 that showed the highest gene expression levels in almost all tissues tested. Group 2 contained AhZIP4.1/4.2, AhZIP6.3/6.4 and AhZIP2.1/2.2, representing intermediate levels of gene expression. Group 3 is composed of the remaining 22 genes; these genes were not express in most peanut tissues or showed low expression levels. AhZIP2.1/2.2 is specifically and highly expressed in the root and nodule. By contrast, AhZIP4.1/4.2 was mainly expressed in roots, nodules and reproductive tissues (i.e., peg tip to fruit, seed and pericarp); AhZIP11.1/11.2 was expressed in all tissues but was relatively higher in reproductive tissues (i.e., peg tip to fruit, fruit and pericarp) and vegetative shoot tips (Figure 6).

2.5. Gene Expression of AhZIPs in Response to Fe- and Zn-Deficiency

The expression of AhZIP genes differed between the two cultivars in response to Fe- and Zn-deficiency (Figure 7). Under the control condition, Fenghua 1 showed higher expressions of AhZIP3.6 and AhZIP11.1 than Silihong, while Silihong showed higher expressions of AhIRT1.1, AhIRT1.2, AhZIP1.1, AhZIP1.2, AhZIP6.1 and AhZIP7.8 than Fenghua 1. Under Fe-deficiency condition, Fenghua 1 showed higher expressions of AhZIP1.1 and AhZIP4.1 than Silihong, while Silihong showed higher expressions of AhIRT1.1, AhZIP3.6 and AhZIP7.8 than Fenghua 1 (Figure 7). Under Zn-deficiency condition, Fenghua 1 showed higher expressions of AhZIP2.1, AhZIP3.5 and AhZIP11.1 than Silihong, while Silihong showed higher expressions of AhIRT1.1, AhIRT1.2 and AhZIP3.6 than Fenghua 1 (Figure 7).
All AhZIP genes tested transcriptionally responded to Fe- and Zn-deficiency in peanut roots dependent on cultivars (Figure 7). Fe-deficiency enhanced the expressions of AhIRT1.1, AhIRT1.2 and AhZIP7.2 but reduced those of AhZIP3.5 and AhZIP4.1 for both cultivars, while other genes showed cultivar differences in response to Fe-deficiency. Fe-deficiency inhibited the expressions of AhZIP1.1, AhZIP1.2 and AhZIP6.1 in Silihong, while they were not affected or increased in Fenghua 1. The expressions of AhZIP3.6, AhZIP7.8 and AhZIP11.1 in Fenghua 1 were inhibited by Fe-deficiency, while in Silihong, they were unaffected.
Zn-deficiency upregulated the expression of AhIRT1.2 for both cultivars, while other genes showed cultivar differences in response to Zn-deficiency (Figure 7). The influence of Zn-deficiency on gene expression was more pronounced in Fenghua 1 than in Silihong. The expressions of AhIRT1.1, AhZIP1.1, AhZIP1.2, AhZIP2.1, AhZIP3.5, AhZIP7.8 and AhZIP11.1 were induced by Zn-deficiency in Fenghua 1, while that in Silihong were unchanged or downregulated. In contrast, AhZIP3.6 were significantly induced by Zn-deficiency in Silihong, whereas in Fenghua 1, it was downregulated.

2.6. Metal Accumulation and Translocation in Response to Fe- or Zn-Deficiency

The two peanut cultivars differed from each other in the accumulation of Fe and Mn (Figure 8). The concentration of Fe in roots and shoots as well as the total amount of Fe in plants were significantly higher in Fenghua 1 than in Silihong. Compared with Silihong, Fenghua 1 showed higher root Mn concentrations but lower shoot Mn concentrations, which resulted in a lower root-to-shoot Mn translocation (Figure 8). Fe-deficiency significantly reduced Fe concentrations in roots for both cultivars; however, shoot Fe concentrations were unchanged. This contributed to an increase in Fe translocation in Fe-deficient plants. The total amount of Fe in plants was decreased by Fe-deficiency in Fenghua 1 but was unaffected in Silihong. Zn-deficiency increased shoot Fe concentrations for both cultivars, while root Fe concentrations, total Fe in plants and the percentage of Fe in shoots were not affected.
The concentrations of Zn in roots and shoots as well as the total amount of Zn in plants were significantly enhanced by Fe-deficiency for both cultivars, while the percentage of Zn in shoots was decreased. Zn-deficiency did not change Zn accumulation, but reduced Zn translocation from roots to shoots in Fenghua 1 (Figure 8). Fe-deficiency significantly increased Mn concentrations in roots and shoots as well as the total amount of Mn in plants for both cultivars. The percentage of Mn in shoots declined by Fe-deficiency in Fenghua 1, while it was unaffected in Silihong. Zn-deficiency reduced Mn translocation for both cultivars, while root Mn concentrations in Silihong were increased (Figure 8).

2.7. Relationship of Gene Expression of AhZIPs and Metal Accumulation

To identify the AhZIP genes involved in metal uptake and translocation in peanut roots, a stepwise linear regression analysis was performed. As showed in Table 3, the expression of AhIRT1.1 is significantly correlated with the total amount of Fe in plants, concentrations of Zn and Mn in shoots, and the percentage of Fe and Mn in shoots. The expression of AhIRT1.2 is significantly correlated with Fe concentrations in roots, the total amount of Mn in plants, and the percentage of Mn in shoots. The expression of AhZIP3.6 is significantly correlated with Mn concentrations in shoots and the percentage of Mn in shoots. The expression of AhZIP6.1 is significantly correlated with Fe concentrations in roots and shoots. The expression of AhZIP7.2 is significantly correlated with Zn concentrations in roots and the total amount of Zn in plants. The expression of AhZIP11.1 is significantly correlated with Mn concentrations in shoots.

3. Discussion

In this study, 30 ZIP members were identified in peanut, which is the largest compared with the reported plant species. For example, 12 ZIP genes were identified in most plant species such as potato, trifoliate orange and maize, while 15, 16 and 21 ZIP members were reported in Arabidopsis, rice and P. trichocarpa, respectively [27,28,29,30,31]. Peanut is an allotetraploid species that essentially contains two complete sets of subgenome (A and B) from two diploid ancestral species: A. duranensis (AA) and A. ipaensis (BB) [39]. Here, we showed that all AhZIPs are multicopy genes, and most of them resulted from WGDs (Figure 3). AhZIP genes are unevenly distributed in the two subgenomes of peanut. The subgenome A contained 13 AhZIP genes, while 17 genes were distributed in the subgenome B (Figure 3). This phenomenon indicates that gene loss or gain occurred during the evolutionary process.
In agreement with previous studies [27,30,31], most of the AhZIP proteins were predicted to localize to the plasma membrane, and AhZIP1.2, AhZIP3.2, AhZIP3.5 and AhZIP7.8 localize to the endomembrane system (Table 1). AtIRT1, AtIRT3 and AtZIP2 from Arabidopsis have been confirmed to be localized to the plasma membrane [9,12,13], while AtIRT2 and AtZIP1 localized to the vesicle or vacuolar membranes [11,13]. In rice, OsIRT1, OsIRT2, OsZIP1, OsZIP4, OsZIP5, OsZIP7 and OsZIP8 have been proven to localize to the plasma membrane [16,17,19,20,21,22,23,24,25,26].
AhZIPs showed a wide variation in TMDs, ranging from 3 to 8 TMDs (Table 1). More than half of AhZIPs have 3–5 TMDs, which is not consistent with the results of Guerinot [15], who proposed that ZIP proteins typically contained 8 TMDs. The variation in TMD number in ZIP proteins has been reported in several plant species such as maize (6–13 TMDs), potato (6–9 TMDs), trifoliate orange (6–9 TMDs) and P. trichocarpa (3–13 TMDs) [27,28,30,31]. The ZIP family was generally predicted to have a 3 + 5 TMD architecture [15]. However, the TMD structure is more diverse in AhZIP proteins. Among 30 AhZIP proteins, only 14 have the 3 + 5 TMD structure (Figure 5). The unusual arrangement of TMDs suggests that the peanut AhZIP family has a distinctive evolutionary process and significant divergence of physiological functions.
The HRD within the long cytoplasmic loop between TM3 and TM4 has been postulated to serve as a potential Zn2+ binding site in many ZIPs [15,43]. The histidine residues have been proven to coordinate Zn2+ in the intracellular loop between TM3 and TM4 of human hZIP4 [43]. In the current study, all AhZIP proteins except AhIRT1.2 and AhZIP7.8 have the large intracellular loop and most of them contain HRDs including HXHXH, HHH, HHHHH, HXHXHXH and HHXHXHXH (Figure 5). Additionally, a glycine (G) residue is always found near or inside HRDs in the intracellular loop between TM3 and TM4 of AhZIP proteins. Moreover, almost all TMDs contained a conserved glycine (Figure 5). These findings indicate that glycine residues might be essential for the structure and function of AhZIP proteins. In human hZIP4, six conserved G residues (G330, G374, G512, G526, G535 and G630) were identified to be located at or near the TM–TM interface, suggesting a role in mediating TM packing [41].
All AhZIP proteins contained the typical domain of the family, zip, and are perfectly modeled on the template, 6pgi.1.A (the A chain of BbZIP protein from B. bronchiseptica). BbZIP as a Zn2+ and Cd2+ transporter has a binuclear metal center, where two metal ions are trapped halfway through the membrane, and the two metal-binding sites play asymmetric roles within the transport pathway [41]. A pairwise sequence alignment showed considerable homology between BbZIP and AhZIPs throughout the TMD (Figure 5). BbZIP has been proven to have six Zn2+ binding sites and one Cd2+ binding site, which were formed by 13 amino acid residues including D89 from TM2; M99 and D144 from TM3; H177, N178 and E181 from TM4; Q207, D208 and E211 from TM5; E240 from TM6; H275 and E276 from TM7; and H286 from TM8 [41]. Almost all of them were mapped to conserved amino acid residues of AhZIP proteins, particularly the two metal-binding motifs “177HNLPEG182” from TM4 and “207QDVPEG212” from TM5, which form the binuclear metal center of BbZIP: M1 (H177, E181, Q207 and E211) and M2 (N178, E181, D208, E211, and E239 from TM6) (Figure 5). The similar structures suggest that AhZIPs might have equivalent physiological functions as BbZIP.
The 30 AhZIP members were divided into four classes: I, II, III and IV (Figure 1). The classification concurred with those reported in previous studies [29,30,31]. Class I, as the largest class, contained 20 AhZIP members, including eight orthologs of AhZIP7, six orthologs of AhZIP3, four orthologs of AhIRT1 and two orthologs of AhZIP1. Most of the AhZIPs belonging to class I shared six motifs (4, 8, 2, 3, 5 and 1), except for AhZIP7.8, AhIRT1.2, AhZIP3.5 and AhZIP1.2, which have less motifs (Figure 4a). Motifs 1, 2, 3 and 5 were annotated as zinc transporters (Table S1). The motif composition indicates that the class I proteins might be responsible for metal transport in peanut.
It is noteworthy that several genes belonging to class I significantly differed from other orthologous genes in the gene/protein structures. AhZIP7.8, encoding 159 aa with four TMDs (TM5-TM8) and one motif (motif 1), has a distinctive gene structure (five exons and four introns). AhIRT1.2 contained three motifs and three TMDs, which are greatly less than the other three orthologs. AhZIP1.2 has three TMDs and four motifs, which are less than half of AhZIP1.1. AhZIP1.2 contained two exons and one intron, while AhZIP1.1 contained three exons and two introns. AhZIP3.5 contained three TMDs (TM1–TM3) and five motifs, while other orthologs contained five TMDs (TM1–TM5) and six motifs. Moreover, AhZIP7.8, AhZIP3.5 and AhZIP1.2 were predicted to be located in the endomembrane system. The distinctive features indicate that these genes might differ from their orthologous genes in physiological functions.
Class II consisted of AhZIP4.1 and AhZIP4.2. The two members are significantly similar in protein structure and physiochemical traits, implying similar functions in peanut plants. However, they differed from each other in gene structure. AhZIP4.1 contained four exons and three introns, while AhZIP4.2 contained two exons and four introns. The difference in exon/intron organization indicates a significant gene divergence during the evolutionary process. Class III includes AhZIP2.1/2.2 and AhZIP11.1/11.2. The four members shared a similar motif composition (motifs 1, 2, 6, 7 and 10). However, the gene structure is different between AhZIP2 (two exons and one intron) and AhZIP11 (three exons and two introns). Class IV contained four orthologous genes of AhZIP6, which were derived from both WGD and segmental duplication. The four members are greatly similar in physiochemical features, TMDs, subcellular location, and gene/protein structures, indicating similar physiological functions. The Ks values of AhZIP6.1/6.3 (1.3068) and AhZIP6.2/6.4 (1.3548) were considerably higher than that of AhZIP6.1/6.2 (0.0477) and AhZIP6.3/6.4 (0.0485) (Table 2), indicating that segmental duplication may occur earlier than WGD.
Gene duplication is a major source of novel genes, contributing to the evolution of new functions [44]. However, before the functional divergence, duplicated genes are usually functionally redundant, which may induce gene loss [45]. To avoid gene loss during evolution processes, the expression of duplicated genes is reduced compared with the ancestral gene [45]. In the present study, 22 AhZIP genes, including all members of class I as well as two orthologs of AhZIP6 (AhZIP6.1/6.2), showed low expression in the 22 peanut tissues under normal conditions, and most of them have more than four orthologous genes (Figure 6). Our results are in accordance with those of Qian et al. [45], suggesting that the low expression of these genes might be beneficial for long-term maintenance of duplicate genes and their functional redundancy. Despite this, the RT-qPCR results indicated that these low expressed genes in peanut roots can be induced by Fe- and Zn-deficiency in a cultivar-dependent manner (Figure 7). It was observed that AhZIP2.1/2.2, AhZIP4.1/4.2 and AhZIP11.1/11.2 were preferentially expressed in roots, nodule and reproductive tissues (i.e., peg tip to fruit, seed and pericarp) (Figure 6). These genes might be involved in metal uptake and translocation by the root and might be responsible for the development of pods or seeds by regulating the metal transport.
The expression of AhZIP genes showed wide differences in response to Fe- and Zn-deficiency depending on cultivars (Figure 7). Fe-deficiency enhanced the expressions of AhIRT1.1 and AhIRT1.2 for both cultivars. Moreover, the expression of AhIRT1.2 was also upregulated by Zn-deficiency for both cultivars. Phylogenetic analysis showed that AhIRT1.1 and AhIRT1.2 have a close relationship with AtIRT1 and AtIRT2 (Figure 1). AtIRT1 is a Fe transporter responsible for taking up Zn, Mn, Co, Ni and Cd from the rhizosphere to root cells [6,7,8,9]. AtIRT2 is involved in the compartmentalization of Fe into vesicles to avoid its toxicity in the cytosol [11]. Fe-deficiency significantly increased the concentration of Zn and Mn in roots and shoots as well as the total amount of Zn and Mn in plants for both cultivars (Figure 8). The expression of AhIRT1.1 significantly correlated with the total amount of Fe in plants, concentrations of Zn and Mn in shoots, and the percentage of Fe and Mn in shoots (Table 3). The expression of AhIRT1.2 significantly correlated with Fe concentrations in roots, the total amount of Mn in plants and the percentage of Mn in shoots (Table 3). These findings suggested that AhIRT1.1 and AhIRT1.2 might be responsible for the uptake of Fe, Zn and Mn in peanut plants. PtIRT1 from trifoliate orange, which was closely related to AhIRT1.1 and AhIRT1.2, has been proven to be responsible for Fe, Zn and Mn uptake [31].
Fe-deficiency induced the expression of AhZIP7.2 for both cultivars. The stepwise linear regression analysis revealed that the expression of AhZIP7.2 significantly correlated with Zn concentrations in roots and the total amount of Zn in plants (Table 3). The results indicate that AhZIP7.2 might be involved in Fe/Zn uptake in peanut plants. The expression of AhZIP3.6 was positively correlated with Mn concentrations in shoots and the percentage of Mn in shoots (Table 3), indicating a possible role in Mn translocation from roots to shoots. Phylogenetic analysis showed that AhZIP1/3/7 has a closer relationship with OsZIP3/4/5/8/9 (Figure 1). OsZIP3 is most likely to be localized at the plasma membrane and responsible for xylem unloading of Zn in the nodes of rice [46]. OsZIP4 is a Zn transporter involved in the translocation of Zn within plants [20,21]. OsZIP5 and OsZIP8 serve as the plasma membrane-localized transporter involved in Zn uptake and distribution within rice [22,23]. OsZIP5 and OsZIP9 are tandem duplicates that act synergistically in Zn/Cd uptake [47].
The expression of AhZIP6.1 was repressed by Zn-deficiency, while its response to Fe-deficiency was dependent on cultivars (Figure 7). Fe-deficiency upregulated the expression of AhZIP6.1 in Fenghua 1 but downregulated that in Silihong. The expression of AhZIP6.1 is significantly correlated with Fe concentrations in roots and shoots (Table 3). These data indicate that AhZIP6.1 might be responsible for the transport of Fe in peanut plants. Phylogenetic analysis showed that AhZIP6.1 has a close relationship to OsZIP6 (Figure 1), which is suggested to competitively take up Fe2+ and Co2+ [48]. TcZNT6 of Thlaspi caerulescens is a metal transporter that is responsible for the transport of Zn, Cd, and Fe or Mn [49]. ZIP6 from Arabidopsis halleri mainly expressed in vascular tissues and encodes a Zn/Cd transporter [50].
The expression of AhZIP11.1 in Fenghua 1 was inhibited by iron deficiency but induced by zinc deficiency, whereas its expression in Silihong was not affected (Figure 7). Stepwise linear regression analysis revealed that the expression of AhZIP11.1 was significantly correlated with Mn concentrations in shoots (Table 3). It seems that AhZIP11.1 is involved in Mn transport in peanut. Similarly, PtZIP11 from trifoliate orange has been suggested to be involved in Mn transport but not Zn or Fe [31]. However, AtZIP11 from Arabidopsis [13] and NtZIP11 from Nicotiana tabacum [51] were shown to mediate Zn uptake but not Fe or Mn.
In agreement with previous studies [52], we found that Silihong showed a higher capacity for tolerance of Fe-deficiency than Fenghua 1 (Figure 8). This might be attributed to lower Fe requirement under normal conditions and higher Fe translocation in case of iron deficiency. In addition, the higher expression of AhIRT1.1 in the root of Silihong under iron deficiency might enhance its capacity for Fe acquirement, therefore contributing to tolerance of Fe-deficiency. Moreover, Silihong has a higher capacity for the translocation of Mn from roots to shoots than Fenghua 1 (Figure 8). Stepwise linear regression analysis revealed that the expressions of AhIRT1.1, AhIRT1.2, and AhZIP3.6 were significantly related to the translocation of Mn in peanut plants (Table 3). Silihong showed higher expressions of AhIRT1.1 and AhZIP3.6 than Fenghua 1. Thus, a higher capacity for Mn translocation in Silihong might result from the higher expression of AhIRT1.1 and AhZIP3.6.

4. Materials and Methods

4.1. Identification of ZIP Genes in Peanut

The protein sequences of 15 AtZIPs from Arabidopsis and 12 OsZIPs from rice were retrieved from Phytozome database (https://phytozome-next.jgi.doe.gov/, accessed on 12 August 2021). The sequences obtained were used as queries for BLASTP against the genome of cultivated peanut (cv. Tifrunner) on Phytozome. All candidates were further examined with the hmmscan tool (https://www.ebi.ac.uk/Tools/hmmer/search/hmmscan, accessed on 15 August 2021). The sequences that contain the ZIP domain (Pfam: PF02535) were recognized as AhZIP proteins. Thereafter, physiochemical parameters including molecular weight, amino acid number, GRAVY and pI were analyzed using the ProtParam tool (https://web.expasy.org/protparam/, accessed on 27 August 2021) [53]. The TMDs of AhZIPs were identified using TOPCONS (http://topcons.net/, accessed on 1 September 2021) [54]. Subcellular targeting sites of AhZIP proteins were predicted using Plant-mPLoc (http://www.csbio.sjtu.edu.cn/bioinf/plant-multi/, accessed on 31 August 2021) [55].

4.2. Phylogenetic and Structural Analysis of AhZIP Proteins

Sequences of ZIP proteins of Arabidopsis, rice and trifoliate orange retrieved from Phytozome database were aligned using ClustalW integrated in MEGA-X software (v. 10.2.6). Based on the sequence alignment, a phylogenetic tree was constructed by the neighbor-joining (NJ) method with the No. of no-difference models. A bootstrap test of 1000 replicates was used for estimating the reliability of the phylogenetic tree. The tree was drawn and modified using iTOL (https://itol.embl.de/itol.cgi, accessed on 5 September 2021) [56].
The conserved motifs and domains of AhZIP proteins were analyzed using the MEME (https://meme-suite.org/meme/tools/meme, accessed on 31 January 2022) [57] and Pfam tool (http://pfam.xfam.org/search#tabview=tab1, accessed on 2 February 2022) [58], respectively, and were visualized using TBtools software [59]. Their homology-modelled 3D structures were predicted using the SwissModel (https://swissmodel.expasy.org/, accessed on 5 February 2022) [60].

4.3. Structure, Duplication and Ka/Ks of AhZIP Genes

The exon/intron structures of AhZIP genes were detected with genomic and coding sequences using GSDS v. 2.0 (http://gsds.gao-lab.org/, accessed on 11 September 2021) [61]. Gene collinearity and Ka/Ks were analyzed by the One Step MCScanX and Simple Ka/Ks calculator (NJ) integrated in the TBtools software, respectively [59]. Diagrams of the exon/intron organization and gene duplication event were drawn using the TBtools software [59].

4.4. Gene Expression Analysis Based on RNA-Seq Data

The expression profiles of the AhZIP genes from Tifrunner were identified using RNA-seq data obtained from the PeanutBase database (https://www.peanutbase.org/, accessed on 6 February 2022) [62]. The read counts were transformed to TPM (Transcripts Per Kilobase of exon model per Million mapped reads), and the heatmap diagram was constructed with lg(TPM + 1) using TBtools [59].

4.5. Plant Growth, Metal Determination and RT-qPCR Analysis

Two peanut cultivars differing in Fe-deficiency tolerance, Fenghua 1 (Fe-deficiency sensitive cultivar) and Silihong (Fe-deficiency tolerant cultivar), were used for determining the relationships between the expression of AhZIP genes and metal accumulation in peanut plants [52]. After surface sterilized with 5% sodium hypochlorite (1 min), the seeds were presoaked in distilled water for 24 h, and then, they were sown in sand for germination. Three-day-old uniform seedlings were transferred to polyethylene pots and cultured as previously reported [63]. The seven-day-old seedlings were treated with Fe- (without Fe) or Zn-deficiency (without Zn) in hydroponic cultures, with the normal nutrition solution containing 50 μM Fe-EDTA and 3.8 μM ZnSO4 as the control [64]. The experiment was arranged in a randomized complete design with four replications (pots) for each treatment. Each replication includes three seedlings. During the growing period, pots were randomly arranged and moved daily for minimizing position effects. After 14 days of treatment, plants were harvested and fresh tissues of roots were sampled for RT-qPCR analysis.
The harvested plants were separated into roots and shoots and oven-dried to constant weight. After weighing, plant tissues were ground into powder and were digested with a mixture of HNO3 and HClO4 (3:1, v/v) using the method described by Su et al. [63]. The levels of Fe, Zn and Mn were determined by flame atomic absorbance spectrometry (WFX-110, Beijing Rayleigh Analytical Instrument Company, Beijing, China).
The expression levels of 12 AhZIP genes representing unique gene and protein structures were detected using RT-qPCR with the method described previously [52] and with Ah60S as the endogenous reference gene. The primers are listed in Table S4. Four biological replicates were performed for each treatment, and three technical replicates were performed for each sample. The relative gene expression was calculated using the 2−ΔΔCT method [65].

4.6. Statistical Analysis

Data were subjected to one-way ANOVA, and significant variations among means were determined by the Duncan’s multiple-range test at 0.05 probability. Stepwise linear regression analysis was performed on the expression of AhZIP genes and metal accumulation. All statistical analyses were conducted using IBM SPSS Statistics v. 22 (IBM, New York, NY, USA).

5. Conclusions

In conclusion, a total of 30 ZIP genes were identified in peanut, which were divided into four classes. All AhZIP proteins contained the typical zip domain and are perfectly modeled on the 6pgi.1.A template, suggesting a role of metal transport in peanut. Unlike previous reports, AhZIP proteins showed a wide variation in TMDs (3–8 TMDs) and only 14 AhZIPs have the 3 + 5 TMD structure. Clustered AhZIPs generally share similar gene/protein structures; however, unique features were found in AhIRT1.2, AhZIP1.2, AhZIP3.5 and AhZIP7.8. Most AhZIP genes showed reduced expression under normal conditions, while AhZIP2.1/2.2, AhZIP4.1/4.2 and AhZIP11.1/11.2 are highly and preferentially expressed in roots, nodule and reproductive tissues, suggesting an essential role in pod and seed development. Transcriptional responses of AhZIPs to Fe/Zn deficiency in peanut roots are dependent on cultivar, which might be, at least partially, responsible for the different metal accumulation between Fenghua 1 and Silihong. The findings provide essential information to further functionally characterize AhZIP genes in the uptake and translocation of metal ions in peanut plants, which is great of importance for screening or breeding cultivars for Fe/Zn biofortification.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/plants11060786/s1, Figure S1: Predicted 3D structure of peanut AhZIP proteins using the SwissModel. Models were visualized in rainbow color from N to C termini; Table S1: Analysis of the ten conserved motifs of AhZIP proteins in peanut; Table S2: The best templates of peanut AhZIP proteins selected from the SwissModel template library for building 3D structure models; Table S3: Expression profiles (TPM) of AhZIP genes in different tissues of peanut; Table S4: Primers of peanut AhZIP genes selected for RT-qPCR analysis.

Author Contributions

Conceptualization, Z.Z. (Zhen Zhang) and G.S.; methodology, G.S.; validation, Z.Z. (Zhen Zhang), Z.Z. (Zheng Zhang) and G.S.; formal analysis, Z.Z. (Zhen Zhang), N.C. and G.S.; investigation, Z.Z. (Zhen Zhang), N.C., Z.Z. (Zheng Zhang) and G.S.; resources, G.S.; data curation, Z.Z. (Zhen Zhang), N.C. and G.S.; writing—original draft preparation, Z.Z. (Zhen Zhang), Z.Z. (Zheng Zhang) and G.S.; writing—review and editing, Z.Z. (Zhen Zhang) and G.S.; supervision, G.S.; project administration, G.S.; funding acquisition, Z.Z. (Zheng Zhang) and G.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by grants from the Natural Science Foundation of Anhui Province (2108085MC83), the Natural Science Foundation for Colleges and Universities of Anhui Province (KJ2020ZD83 and KJ2019A0587) and the Innovation Team of Scientific Research Platform of Anhui Province (KJ2015TD001).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Marschner, H. Mineral Nutrition of Higher Plants, 2nd ed.; Academic Press: Boston, MA, USA, 2012. [Google Scholar]
  2. Gómez-Galera, S.; Rojas, E.; Sudhakar, D.; Zhu, C.; Pelacho, A.M.; Capell, T.; Christou, P. Critical evaluation of strategies for mineral fortification of staple food crops. Transgenic Res. 2010, 19, 165–180. [Google Scholar] [CrossRef]
  3. Guerinot, M.L.; Yi, Y. Iron: Nutritious, noxious, and not readily available. Plant Physiol. 1994, 104, 815–820. [Google Scholar] [CrossRef] [Green Version]
  4. Alloway, B.J. Soil factors associated with zinc deficiency in crops and humans. Environ. Geochem. Health 2009, 31, 537–548. [Google Scholar] [CrossRef] [PubMed]
  5. Ma, G.; Jin, Y.; Li, Y.; Zhai, F.; Kok, F.J.; Jacobsen, E.; Yang, X. Iron and zinc deficiencies in China: What is a feasible and cost-effective strategy? Public Health Nutr. 2008, 11, 632–638. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Korshunova, Y.O.; Eide, D.; Clark, W.G.; Guerinot, M.L.; Pakrasi, H.B. The IRT1 protein from Arabidopsis thaliana is a metal transporter with a broad substrate range. Plant Mol. Biol. 1999, 40, 37–44. [Google Scholar] [CrossRef] [PubMed]
  7. Eide, D.; Broderius, M.; Fett, J.; Guerinot, M.L. A novel iron-regulated metal transporter from plants identified by functional expression in yeast. Proc. Natl. Acad. Sci. USA 1996, 93, 5624–5628. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Nishida, S.; Tsuzuki, C.; Kato, A.; Aisu, A.; Yoshida, J.; Mizuno, T. AtIRT1, the primary iron uptake transporter in the root, mediates excess nickel accumulation in arabidopsis thaliana. Plant Cell Physiol. 2011, 52, 1433–1442. [Google Scholar] [CrossRef] [Green Version]
  9. Vert, G.; Grotz, N.; Dédaldéchamp, F.; Gaymard, F.; Guerinot, M.L.; Briat, J.F.; Curie, C. IRT1, an Arabidopsis transporter essential for iron uptake from the soil and for plant growth. Plant Cell. 2002, 14, 1223–1233. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Krausko, M.; Labajová, M.; Peterková, D.; Jásik, J. Specific expression of AtIRT1 in phloem companion cells suggests its role in iron translocation in aboveground plant organs. Plant Signal. Behav. 2021, 16, 1925020. [Google Scholar] [CrossRef] [PubMed]
  11. Vert, G.; Barberon, M.; Zelazny, E.; Séguéla, M.; Briat, J.-F.; Curie, C. Arabidopsis IRT2 cooperates with the high-affinity iron uptake system to maintain iron homeostasis in root epidermal cells. Planta 2009, 229, 1171–1179. [Google Scholar] [CrossRef] [PubMed]
  12. Lin, Y.-F.; Liang, H.-M.; Yang, S.-Y.; Boch, A.; Clemens, S.; Chen, C.-C.; Wu, J.-F.; Huang, J.-L.; Yeh, K.-C. Arabidopsis IRT3 is a zinc-regulated and plasma membrane localized zinc/iron transporter. New Phytol. 2009, 182, 392–404. [Google Scholar] [CrossRef] [PubMed]
  13. Milner, M.J.; Seamon, J.; Craft, E.; Kochian, L.V. Transport properties of members of the ZIP family in plants and their role in Zn and Mn homeostasis. J. Exp. Bot. 2013, 64, 369–381. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Grotz, N.; Fox, T.; Connolly, E.; Park, W.; Guerinot, M.L.; Eide, D. Identification of a family of zinc transporter genes from Arabidopsis that respond to zinc deficiency. Proc. Natl. Acad. Sci. USA 1998, 95, 7220–7224. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Guerinot, M.L. The ZIP family of metal transporters. Biochim. Biophys. Acta (BBA)-Biomembr. 2000, 1465, 190–198. [Google Scholar] [CrossRef] [Green Version]
  16. Ishimaru, Y.; Suzuki, M.; Tsukamoto, T.; Suzuki, K.; Nakazono, M.; Kobayashi, T.; Wada, Y.; Watanabe, S.; Matsuhashi, S.; Takahashi, M.; et al. Rice plants take up iron as an Fe3+-phytosiderophore and as Fe2+. Plant J. 2006, 45, 335–346. [Google Scholar] [CrossRef]
  17. Nakanishi, H.; Ogawa, I.; Ishimaru, Y.; Mori, S.; Nishizawa, N.K. Iron deficiency enhances cadmium uptake and translocation mediated by the Fe2+ transporters OsIRT1 and OsIRT2 in rice. Soil Sci. Plant Nutr. 2006, 52, 464–469. [Google Scholar] [CrossRef]
  18. Lee, S.; An, G. Over-expression of OsIRT1 leads to increased iron and zinc accumulations in rice. Plant Cell Environ. 2009, 32, 408–416. [Google Scholar] [CrossRef] [PubMed]
  19. Liu, X.S.; Feng, S.J.; Zhang, B.Q.; Wang, M.Q.; Cao, H.W.; Rono, J.K.; Chen, X.; Yang, Z.M. OsZIP1 functions as a metal efflux transporter limiting excess zinc, copper and cadmium accumulation in rice. BMC Plant Biol. 2019, 19, 283. [Google Scholar] [CrossRef] [PubMed]
  20. Ishimaru, Y.; Suzuki, M.; Kobayashi, T.; Takahashi, M.; Nakanishi, H.; Mori, S.; Nishizawa, N.K. OsZIP4, a novel zinc-regulated zinc transporter in rice. J. Exp. Bot. 2005, 56, 3207–3214. [Google Scholar] [CrossRef] [PubMed]
  21. Ishimaru, Y.; Masuda, H.; Suzuki, M.; Bashir, K.; Takahashi, M.; Nakanishi, H.; Mori, S.; Nishizawa, N.K. Overexpression of the OsZIP4 zinc transporter confers disarrangement of zinc distribution in rice plants. J. Exp. Bot. 2007, 58, 2909–2915. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Lee, S.; Jeong, H.J.; Kim, S.A.; Lee, J.; Guerinot, M.L.; An, G. OsZIP5 is a plasma membrane zinc transporter in rice. Plant Mol. Biol. 2010, 73, 507–517. [Google Scholar] [CrossRef] [PubMed]
  23. Lee, S.; Kim, S.A.; Lee, J.; Guerinot, M.L.; An, G. Zinc deficiency-inducible OsZIP8 encodes a plasma membrane-localized zinc transporter in rice. Mol. Cells 2010, 29, 551–558. [Google Scholar] [CrossRef] [PubMed]
  24. Gindri, R.G.; Navarro, B.B.; da Cruz Dias, P.V.; Tarouco, C.P.; Nicoloso, F.T.; Brunetto, G.; Berghetti, L.P.; da Silva, L.O.S.; Fett, J.P.; Menguer, P.K.; et al. Physiological responses of rice (Oryza sativa L.) oszip7 loss-of-function plants exposed to varying Zn concentrations. Physiol. Mol. Biol. Plants 2020, 26, 1349–1359. [Google Scholar] [CrossRef] [PubMed]
  25. Ricachenevsky, F.K.; Punshon, T.; Lee, S.; Oliveira, B.H.N.; Trenz, T.S.; Maraschin, F.D.S.; Hindt, M.N.; Danku, J.; Salt, D.E.; Fett, J.P.; et al. Elemental profiling of rice FOX lines leads to characterization of a new Zn plasma membrane trans-porter, OsZIP7. Front. Plant Sci. 2018, 9, 865. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Tan, L.; Zhu, Y.; Fan, T.; Peng, C.; Wang, J.; Sun, L.; Chen, C. OsZIP7 functions in xylem loading in roots and inter-vascular transfer in nodes to deliver Zn/Cd to grain in rice. Biochem. Biophys. Res. Commun. 2019, 512, 112–118. [Google Scholar] [CrossRef] [PubMed]
  27. Mondal, T.K.; Ganie, S.A.; Rana, M.K.; Sharma, T.R. Genome-wide analysis of zinc transporter genes of maize (Zea mays). Plant Mol. Biol. Rep. 2013, 32, 605–616. [Google Scholar] [CrossRef]
  28. Li, X.B.; Suo, H.C.; Liu, J.T.; Wang, L.; Li, C.C.; Liu, W. Genome-wide identification and expression analysis of the potato ZIP gene family under Zn-deficiency. Biol. Plant. 2020, 64, 845–855. [Google Scholar] [CrossRef]
  29. Migeon, A.; Blaudez, D.; Wilkins, O.; Montanini, B.; Campbell, M.M.; Richaud, P.; Thomine, S.; Chalot, M. Genome-wide analysis of plant metal transporters, with an emphasis on poplar. Cell. Mol. Life Sci. 2010, 67, 3763–3784. [Google Scholar] [CrossRef]
  30. Zhang, H.; Zhao, S.; Li, D.; Xu, X.; Li, C. Genome-wide analysis of the zrt, irt-like protein (zip) family and their responses to metal stress in Populus trichocarpa. Plant Mol. Biol. Rep. 2017, 35, 534–549. [Google Scholar] [CrossRef]
  31. Fu, X.Z.; Zhou, X.; Xing, F.; Ling, L.L.; Chun, C.P.; Cao, L.; Aarts, M.G.M.; Peng, L.Z. Genome-wide identification, cloning and functional analysis of the zinc/iron-regulated transporter-like protein (ZIP) gene family in trifoliate orange (Poncirus trifoliata L. Raf.). Front. Plant Sci. 2017, 8, 588. [Google Scholar] [CrossRef] [Green Version]
  32. Gaither, L.A.; Eide, D.J. Eukaryotic zinc transporters and their regulation. BioMetals 2001, 14, 251–270. [Google Scholar] [CrossRef] [PubMed]
  33. Arunachalam, P.; Kannan, P.; Prabukumar, G.; Govindaraj, M. Zinc deficiency in Indian soils with special focus to enrich zinc in peanut. Afr. J. Agric. Res. 2013, 8, 6681–6688. [Google Scholar]
  34. Lal, C.; Singh, A.L. Screening for high zinc density groundnut genotypes in India. In Proceedings of the Zinc Crops 2007 Conference for Improving Crop Production and Human Health, Istanbul, Turkey, 24–26 May 2007. [Google Scholar]
  35. Su, Y.; Zhang, Z.; Su, G.; Liu, J.; Liu, C.; Shi, G. Genotypic differences in spectral and photosynthetic response of peanut to iron deficiency. J. Plant Nutr. 2015, 38, 145–160. [Google Scholar] [CrossRef]
  36. Ding, H.; Duan, L.; Jing, L.; Yan, H.; Zhao, M.; Zhang, F.; Li, W.X. Cloning and functional analysis of the peanut iron trans-porter AhIRT1 during iron deficiency stress and intercropping with maize. J. Plant Physiol. 2010, 167, 996–1002. [Google Scholar] [CrossRef] [PubMed]
  37. Xiong, H.; Kobayashi, T.; Kakei, Y.; Senoura, T.; Nakazono, M.; Takahashi, H.; Nakanishi, H.; Shen, H.; Duan, P.; Guo, X.; et al. AhNRAMP1 iron transporter is involved in iron acquisition in peanut. J. Exp. Bot. 2012, 63, 4437–4446. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Bertioli, D.J.; Cannon, S.B.; Froenicke, L.; Huang, G.; Farmer, A.D.; Cannon, E.K.S.; Liu, X.; Gao, D.; Clevenger, J.; Dash, S.; et al. The genome sequences of Arachis duranensis and Arachis ipaensis, the diploid ancestors of cultivated peanut. Nat. Genet. 2016, 48, 438–446. [Google Scholar] [CrossRef] [PubMed]
  39. Bertioli, D.J.; Jenkins, J.; Clevenger, J.; Dudchenko, O.; Gao, D.; Seijo, G.; Leal-Bertioli, S.C.M.; Ren, L.; Farmer, A.D.; Pandey, M.K.; et al. The genome sequence of segmental allotetraploid peanut Arachis hypogaea. Nat. Genet. 2019, 51, 877–884. [Google Scholar] [CrossRef] [Green Version]
  40. Hurst, L.D. The Ka/Ks ratio: Diagnosing the form of sequence evolution. Trends Genet. 2002, 18, 486–487. [Google Scholar] [CrossRef]
  41. Zhang, T.; Liu, J.; Fellner, M.; Zhang, C.; Sui, D.; Hu, J. Crystal structures of a ZIP zinc transporter reveal a binuclear metal center in the transport pathway. Sci. Adv. 2017, 3, e1700344. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Pattee, H.E.; Johns, E.B.; Singleton, J.A.; Sanders, T.H. Composition changes of peanut fruit parts during maturation1. Peanut Sci. 1974, 1, 57–62. [Google Scholar] [CrossRef]
  43. Bafaro, E.M.; Antala, S.; Nguyen, T.-V.; Dzul, S.P.; Doyon, B.; Stemmler, T.L.; Dempski, R.E. The large intracellular loop of hZIP4 is an intrinsically disordered zinc binding domain. Metallomics 2015, 7, 1319–1330. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Panchy, N.; Lehti-Shiu, M.; Shiu, S.-H. Evolution of gene duplication in plants. Plant Physiol. 2016, 171, 2294–2316. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Qian, W.; Liao, B.-Y.; Chang, A.Y.-F.; Zhang, J. Maintenance of duplicate genes and their functional redundancy by reduced expression. Trends Genet. 2010, 26, 425–430. [Google Scholar] [CrossRef] [Green Version]
  46. Sasaki, A.; Yamaji, N.; Mitani-Ueno, N.; Kashino, M.; Ma, J.F. A node-localized transporter OsZIP3 is responsible for the preferential distribution of Zn to developing tissues in rice. Plant J. 2015, 84, 374–384. [Google Scholar] [CrossRef] [Green Version]
  47. Tan, L.; Qu, M.; Zhu, Y.; Peng, C.; Wang, J.; Gao, D.; Chen, C. Zinc transporter5 and zinc transporter9 function synergistically in zinc/cadmium uptake. Plant Physiol. 2020, 183, 1235–1249. [Google Scholar] [CrossRef] [PubMed]
  48. Kavitha, P.G.; Sam, K.; Mathew, M.K. Functional characterization of a transition metal ion transporter, OsZIP6 from rice (Oryza sativa L.). Plant Physiol. Biochem. 2015, 97, 165–174. [Google Scholar]
  49. Wu, J.; Zhao, F.J.; Ghandilyan, A.; Logoteta, B.; Guzman, M.O.; Schat, H.; Wang, X.; Aarts, M.G.M. Identification and func-tional analysis of two ZIP metal transporters of the hyperaccumulator Thlaspi caerulescens. Plant Soil. 2009, 325, 79. [Google Scholar] [CrossRef] [Green Version]
  50. Spielmann, J.; Ahmadi, H.; Scheepers, M.; Weber, M.; Nitsche, S.; Carnol, M.; Bosman, B.; Kroymann, J.; Motte, P.; Clemens, S.; et al. The two copies of the zinc and cadmium ZIP6 transporter of Arabidopsis halleri have distinct effects on cadmium tolerance. Plant Cell Environ. 2020, 43, 2143–2157. [Google Scholar] [CrossRef] [PubMed]
  51. Kozak, K.; Papierniak, A.; Barabasz, A.; Kendziorek, M.; Palusińska, M.; Williams, L.E.; Antosiewicz, D.M. NtZIP11, a new Zn transporter specifically upregulated in tobacco leaves by toxic Zn level. Environ. Exp. Bot. 2018, 157, 69–78. [Google Scholar] [CrossRef] [Green Version]
  52. Cao, Q.; Xv, C.; Jiang, Q.; Wang, L.; Shi, G. Comparative transcriptome analysis reveals key genes responsible for the homeo-stasis of iron and other divalent metals in peanut roots under iron deficiency. Plant Soil. 2019, 445, 513–531. [Google Scholar] [CrossRef]
  53. Duvaud, S.; Gabella, C.; Lisacek, F.; Stockinger, H.; Ioannidis, V.; Durinx, C. Expasy, the Swiss bioinformatics resource portal, as designed by its users. Nucleic Acids Res. 2021, 49, W216–W227. [Google Scholar] [CrossRef] [PubMed]
  54. Tsirigos, K.D.; Peters, C.; Shu, N.; Käll, L.; Elofsson, A. The TOPCONS web server for consensus prediction of membrane protein topology and signal peptides. Nucleic Acids Res. 2015, 43, W401–W407. [Google Scholar] [CrossRef] [PubMed]
  55. Chou, K.C.; Shen, H.B. Plant-mPLoc: A top-down strategy to augment the power for predicting plant protein subcellular lo-calization. PLoS ONE. 2010, 5, e11335. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Letunic, I.; Bork, P. Interactive Tree of Life (iTOL) v5: An online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021, 49, W293–W296. [Google Scholar] [CrossRef]
  57. Bailey, T.L.; Williams, N.; Misleh, C.; Li, W.W. MEME: Discovering and analyzing DNA and protein sequence motifs. Nucleic Acids Res. 2006, 34, W369–W373. [Google Scholar] [CrossRef]
  58. Mistry, J.; Chuguransky, S.; Williams, L.; Qureshi, M.; Salazar, G.A.; Sonnhammer, E.L.L.; Tosatto, S.C.E.; Paladin, L.; Raj, S.; Richardson, L.J.; et al. Pfam: The protein families database in 2021. Nucleic Acids Res. 2021, 49, D412–D419. [Google Scholar] [CrossRef] [PubMed]
  59. Chen, C.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.; Xia, R. TB tools: An integrative toolkit developed for in-teractive analyses of big biological data. Mol. Plant. 2020, 13, 1194–1202. [Google Scholar] [CrossRef]
  60. Waterhouse, A.; Bertoni, M.; Bienert, S.; Studer, G.; Tauriello, G.; Gumienny, R.; Heer, F.T.; De Beer, T.A.P.; Rempfer, C.; Bordoli, L.; et al. SWISS-MODEL: Homology modelling of protein structures and complexes. Nucleic Acids Res. 2018, 46, W296–W303. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Hu, B.; Jin, J.P.; Guo, A.-Y.; Zhang, H.; Luo, J.C.; Gao, G. GSDS 2.0: An upgraded gene feature visualization server. Bioinformatics 2015, 31, 1296–1297. [Google Scholar] [CrossRef] [Green Version]
  62. Clevenger, J.; Chu, Y.; Scheffler, B.; Ozias-Akins, P. A developmental transcriptome map for allotetraploid Arachis hypogaea. Front. Plant Sci. 2016, 7, 1446. [Google Scholar] [CrossRef] [Green Version]
  63. Su, Y.; Liu, J.; Lu, Z.; Wang, X.; Zhang, Z.; Shi, G. Effects of iron deficiency on subcellular distribution and chemical forms of cadmium in peanut roots in relation to its translocation. Environ. Exp. Bot. 2013, 97, 40–48. [Google Scholar] [CrossRef]
  64. Lu, Z.; Zhang, Z.; Su, Y.; Liu, C.; Shi, G. Cultivar variation in morphological response of peanut roots to cadmium stress and its relation to cadmium accumulation. Ecotoxicol. Environ. Saf. 2013, 91, 147–155. [Google Scholar] [CrossRef]
  65. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR. Methods 2002, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
Plants 11 00786 g001 550
Figure 1. Phylogenetic relationships of ZIP proteins in peanut and other plant species. The species involved in the evolutionary tree include peanut (AhZIP), Arabidopsis thaliana (AtZIP), Oryza sativa (OsZIP), and Poncirus trifoliata (PtZIP). The 30 AhZIP proteins of peanut are marked in red.
Figure 1. Phylogenetic relationships of ZIP proteins in peanut and other plant species. The species involved in the evolutionary tree include peanut (AhZIP), Arabidopsis thaliana (AtZIP), Oryza sativa (OsZIP), and Poncirus trifoliata (PtZIP). The 30 AhZIP proteins of peanut are marked in red.
Plants 11 00786 g001
Plants 11 00786 g002 550
Figure 2. Phylogenetic relationships and exon–intron organization of AhZIP genes from peanut. UTR and CDS represent untranslated regions and coding sequences, respectively.
Figure 2. Phylogenetic relationships and exon–intron organization of AhZIP genes from peanut. UTR and CDS represent untranslated regions and coding sequences, respectively.
Plants 11 00786 g002
Plants 11 00786 g003 550
Figure 3. Chromosomal locations and duplications of peanut AhZIP genes obtained from collinearity analysis.
Figure 3. Chromosomal locations and duplications of peanut AhZIP genes obtained from collinearity analysis.
Plants 11 00786 g003
Plants 11 00786 g004 550
Figure 4. Distributions of the conserved motifs (a) and domains (b) in AhZIP proteins from peanut.
Figure 4. Distributions of the conserved motifs (a) and domains (b) in AhZIP proteins from peanut.
Plants 11 00786 g004
Plants 11 00786 g005 550
Figure 5. Multiple sequence alignment of BbZIP and AhZIP proteins. BbZIP protein from Bordetella bronchiseptica and 30 AhZIP proteins were aligned using ClustalW. The red bars indicate that amino acid residues possibly formed Zn2+ and Cd2+ binding sites in the BbZIP protein. The yellow highlighted motifs indicate histidine-rich domains (HRDs) in the long cytoplasmic loop between TM3 and TM4.
Figure 5. Multiple sequence alignment of BbZIP and AhZIP proteins. BbZIP protein from Bordetella bronchiseptica and 30 AhZIP proteins were aligned using ClustalW. The red bars indicate that amino acid residues possibly formed Zn2+ and Cd2+ binding sites in the BbZIP protein. The yellow highlighted motifs indicate histidine-rich domains (HRDs) in the long cytoplasmic loop between TM3 and TM4.
Plants 11 00786 g005
Plants 11 00786 g006 550
Figure 6. Expression profiles of AhZIP genes across the different tissues. Gene expression is expressed in lg(TPM + 1). Pattee 1, 3, 5, 6, 7, 8 and 10 represent different developmental stages of peanut pods according to Pattee et al. [42], who classified peanut pod maturity into 15 categories.
Figure 6. Expression profiles of AhZIP genes across the different tissues. Gene expression is expressed in lg(TPM + 1). Pattee 1, 3, 5, 6, 7, 8 and 10 represent different developmental stages of peanut pods according to Pattee et al. [42], who classified peanut pod maturity into 15 categories.
Plants 11 00786 g006
Plants 11 00786 g007 550
Figure 7. Expression levels of 12 AhZIP genes in the root of two peanut cultivars in response to Fe/Zn deficiency. Data (means ± SE, n = 4) sharing the same letter(s) above the error bars are not significantly different at the 0.05 level based on the Duncan multiple range test.
Figure 7. Expression levels of 12 AhZIP genes in the root of two peanut cultivars in response to Fe/Zn deficiency. Data (means ± SE, n = 4) sharing the same letter(s) above the error bars are not significantly different at the 0.05 level based on the Duncan multiple range test.
Plants 11 00786 g007
Plants 11 00786 g008 550
Figure 8. The accumulation and translocation of Fe, Mn and Zn in two peanut cultivars in response to Fe- or Zn-deficiency. Data (means ± SE, n = 4) shared the same letter(s) above the error bars are not significantly different at the 0.05 level based on the Duncan multiple range test.
Figure 8. The accumulation and translocation of Fe, Mn and Zn in two peanut cultivars in response to Fe- or Zn-deficiency. Data (means ± SE, n = 4) shared the same letter(s) above the error bars are not significantly different at the 0.05 level based on the Duncan multiple range test.
Plants 11 00786 g008
Table
Table 1. Molecular characterization of AhZIP genes identified in peanut.
Table 1. Molecular characterization of AhZIP genes identified in peanut.
Gene NameGene IDGene Length (bp)CDS (bp)MW a (kDa)Aa bInstabilityAliphatic IndexGRAVY cpI dNo. of TMD eLocation
AhIRT1.1arahy.T4CX6H8962109538.9436433.45109.290.556.38/out-outPM f
AhIRT1.2arahy.B3ZT22237345616.5515147.1109.070.6545.283/out-inPM
AhIRT1.3arahy.234RNS4415108939.2536236.8111.990.5736.868/out-outPM
AhIRT1.4arahy.8VMZ7D4094110139.5636637.5108.630.536.348/out-outPM
AhZIP1.1arahy.XJH13Y5732114040.5137941.77104.010.4776.298/out-outPM
AhZIP1.2arahy.ZLZ7ZM160871425.7423744.0186.030.1676.533/out-inEMS g
AhZIP2.1arahy.VX1J702322106238.7735326.98107.480.4956.528/in-outPM
AhZIP2.2arahy.1Q0IUD2309105938.7135227.91103.350.4726.628/in-outPM
AhZIP3.1arahy.CK2LDM476689432.3829739.12107.710.2546.035/out-inPM
AhZIP3.2arahy.E7VKLQ174783130.1227638.03110.620.2886.625/out-inEMS
AhZIP3.3arahy.ZVRF07155478328.2626033.11107.190.4236.085/out-inPM
AhZIP3.4arahy.30BR38121169925.0923233.22109.660.5685.285/out-inPM
AhZIP3.5arahy.WQ3KQR261888832.2629537.4299.490.2126.453/out-inEMS
AhZIP3.6arahy.05ZQZB93774726.9924831.92117.140.5245.915/out-inPM
AhZIP4.1arahy.KIJ6L72325122743.8240842.9395.220.3246.028/out-outPM
AhZIP4.2arahy.2FF2JF2561110138.8336641.9695.250.346.028/out-outPM
AhZIP6.1arahy.58GIJL218599935.6033235.92106.990.6315.778/out-outPM
AhZIP6.2arahy.78T540281299935.5433235.89107.290.6475.98/out-outPM
AhZIP6.3arahy.0DI5A23196103536.2934437.88106.10.6296.298/out-outPM
AhZIP6.4arahy.E6QUMC3026103536.3134437.66106.370.6306.298/out-outPM
AhZIP7.1arahy.BTX8K3158581929.6727245.6103.20.4006.105/out-inPM
AhZIP7.2arahy.092K8V151481929.6127244.17104.30.4136.585/out-inPM
AhZIP7.3arahy.RP74H9139881029.4526945.52109.440.4126.045/out-inPM
AhZIP7.4arahy.FX0GN0115181929.8827244.84103.570.3476.625/out-inPM
AhZIP7.5arahy.1C0EWX151082229.6827344.23106.780.4296.625/out-inPM
AhZIP7.6arahy.QZI7QE120784330.3328046.91104.110.3996.625/out-inPM
AhZIP7.7arahy.0E1GBK139782229.8827345.75104.620.4086.15/out-inPM
AhZIP7.8arahy.NIU36G386948017.6915946.8399.250.4326.574/out-outEMS
AhZIP11.1arahy.HSP4SF3404103836.6234530.58112.870.7675.738/out-outPM
AhZIP11.2arahy.W56MR22492104736.9434832.71113.050.7495.498/out-outPM
a Molecular weight, b amino acid number, c grand average of hydropathicity, d isoelectric points, e transmembrane domain, f plasma membrane, g endomembrane system.
Table
Table 2. Ka/Ks analysis of all gene duplication pairs for AhZIP genes.
Table 2. Ka/Ks analysis of all gene duplication pairs for AhZIP genes.
Gene PairsDuplicate TypeKa aKs bKa/Ks cPositive Selection
AhZIP6.1/6.4Segmental0.19121.37310.1393No
AhZIP6.1/6.3Segmental0.18841.30680.1442No
AhZIP6.2/6.4Segmental0.18871.35480.1393No
AhZIP6.3/6.2Segmental0.18591.29020.1441No
AhIRT1.2/1.4Whole-genome0.04630.09560.4843No
AhZIP1.1/1.2Whole-genome0.03380.06670.5068No
AhZIP2.1/2.2Whole-genome0.01120.03730.3006No
AhZIP3.1/3.2Whole-genome0.03930.04740.8291No
AhZIP3.3/3.4Whole-genome0.02080.08610.2415No
AhZIP4.1/4.2Whole-genome0.00240.04600.0525No
AhZIP6.3/6.4Whole-genome0.00260.04850.0532No
AhZIP6.1/6.2Whole-genome0.00130.04770.0277No
AhZIP7.3/7.4Whole-genome0.01800.08480.2119No
AhZIP7.1/7.7Whole-genome0.02270.08350.2717No
AhZIP7.5/7.6Whole-genome0.01130.06980.1620No
AhZIP11.1/11.2Whole-genome0.00770.02820.2733No
a The number of nonsynonymous substitutions per nonsynonymous site, b the number of synonymous substitutions per synonymous site, c Ka/Ks ratios.
Table
Table 3. Stepwise linear regression analysis (β value) of metal accumulation and the expression of AhZIP genes in the roots of Fenghua 1 and Silihong (n = 24).
Table 3. Stepwise linear regression analysis (β value) of metal accumulation and the expression of AhZIP genes in the roots of Fenghua 1 and Silihong (n = 24).
Gene Expression a[Fe]root b[Fe]shoot cTotal Fe in
Plants		% of Fe in Shoots[Zn]root d[Zn]shoot eTotal Zn in Plants[Mn]shoot fTotal Mn in Plants% of Mn in Shoots
AhIRT1.1−0.07−0.08−0.52 **0.65 **0.210.60 **0.161.19 ***0.050.95 ***
AhIRT1.2−0.39 *0.020.100.160.300.320.280.060.44 *−0.34 *
AhZIP1.10.100.11−0.14−0.07−0.21−0.17−0.20−0.06−0.270.15
AhZIP1.2−0.010.12−0.24−0.050.130.03−0.12−0.05−0.260.10
AhZIP2.10.010.21−0.300.070.180.16−0.12−0.04−0.110.01
AhZIP3.50.060.27−0.140.090.260.280.09−0.12−0.05−0.14
AhZIP3.6−0.12−0.15−0.320.110.040.10−0.260.57 **−0.250.70 ***
AhZIP4.10.000.13−0.35−0.11−0.06−0.02−0.17−0.18−0.220.05
AhZIP6.1−0.60 **−0.42 *0.020.060.02−0.040.100.130.37−0.15
AhZIP7.20.10−0.040.150.090.49 *0.060.45 *0.070.28−0.10
AhZIP7.80.010.13−0.27−0.080.210.13−0.10−0.03−0.18−0.01
AhZIP11.10.110.230.190.26−0.050.43−0.030.56 *0.20−0.14
a Gene expression is calculated as −ΔΔCT, b Fe concentration in roots, c Fe concentration in shoots, d Zn concentration in roots, e Zn concentration in shoots, f Zn concentration in shoots, * p < 0.05, ** p < 0.01, *** p < 0.001.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, Z.; Chen, N.; Zhang, Z.; Shi, G. Genome-Wide Identification and Expression Profile Reveal Potential Roles of Peanut ZIP Family Genes in Zinc/Iron-Deficiency Tolerance. Plants 2022, 11, 786. https://doi.org/10.3390/plants11060786

AMA Style

Zhang Z, Chen N, Zhang Z, Shi G. Genome-Wide Identification and Expression Profile Reveal Potential Roles of Peanut ZIP Family Genes in Zinc/Iron-Deficiency Tolerance. Plants. 2022; 11(6):786. https://doi.org/10.3390/plants11060786

Chicago/Turabian Style

Zhang, Zhen, Nannan Chen, Zheng Zhang, and Gangrong Shi. 2022. "Genome-Wide Identification and Expression Profile Reveal Potential Roles of Peanut ZIP Family Genes in Zinc/Iron-Deficiency Tolerance" Plants 11, no. 6: 786. https://doi.org/10.3390/plants11060786

APA Style

Zhang, Z., Chen, N., Zhang, Z., & Shi, G. (2022). Genome-Wide Identification and Expression Profile Reveal Potential Roles of Peanut ZIP Family Genes in Zinc/Iron-Deficiency Tolerance. Plants, 11(6), 786. https://doi.org/10.3390/plants11060786

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop