Genome-Wide Identiﬁcation and Expression Analysis of AMT and NRT Gene Family in Pecan ( Carya illinoinensis ) Seedlings Revealed a Preference for NH 4+ -N

: Nitrogen (N) is a major limiting factor for plant growth and crop production. The use of N fertilizer in forestry production is increasing each year, but the loss is substantial. Mastering the regulatory mechanisms of N uptake and transport is a key way to improve plant nitrogen use efﬁciency (NUE). However, this has rarely been studied in pecans. In this study, 10 AMT and 69 NRT gene family members were identiﬁed and systematically analyzed from the whole pecan genome using a bioinformatics approach, and the expression patterns of AMT and NRT genes and the uptake characteristics of NH 4+ and NO 3 − in pecan were analyzed by aeroponic cultivation at varying NH 4+ /NO 3 − ratios (0/0, 0/100,25/75, 50/50, 75/25,100/0 as CK, T1, T2, T3, T4, and T5). The results showed that gene duplication was the main reason for the ampliﬁcation of the AMT and NRT gene families in pecan, both of which experienced purifying selection. Based on qRT-PCR results, CiAMTs were primarily expressed in roots, and CiNRTs were majorly expressed in leaves, which were consistent with the distribution of pecan NH 4+ and NO 3 − concentrations in the organs. The expression levels of CiAMTs and CiNRTs were mainly signiﬁcantly upregulated under N deﬁciency and T4 treatment. Meanwhile, T4 treatment signiﬁcantly increased the NH 4+ , NO 3 − , and NO 2 − concentrations as well as the Vmax and Km values of NH 4+ and NO 3 − in pecans, and Vmax/Km indicated that pecan seedlings preferred to absorb NH 4+ . In summary, considering the single N source of T5, we suggested that the NH 4+ /NO 3 − ratio of 75:25 was more beneﬁcial to improve the NUE of pecan, thus increasing pecan yield, which provides a theoretical basis for promoting the scale development of pecan and provides a basis for further identiﬁcation of the functions of AMT and NRT genes in the N uptake and transport process of pecan.


Introduction
Ammonium nitrogen (NH 4 + ) and nitrate nitrogen (NO 3 − ) are the two main forms of nitrogen (N) that plants can absorb and use. The absorption, transport, and assimilation of NH 4 + and NO 3 − in plants have been studied extensively. Physiologically, NH 4 + and NO 3 − uptake and translocation systems were identified mainly by root absorption kinetics, which were classified into two types: high-affinity transport systems (HATs) and low-affinity transport pathways (LATs) [1]. Vmax (maximum ion uptake rate) and Km (Mie's constant) are the two main parameters in the kinetic equation for root nutrient uptake, which can quantitatively characterize the plant uptake of nutrient ions. In general, a larger Vmax value indicates that the plant has a great uptake potential for a certain ion, and the number of the ion transport carrier protein on the cell membrane determines the size of Vmax. Km indicates the affinity between the ion absorbed by the root system and the uptake site (transport carrier), and the greater the affinity, the smaller the Km value [2]. Previous studies have found that plants have a preference for the intake and utilization of NH 4 + and NO 3 − ; when both N forms are present, plants will preferentially intake and utilize one of them [3]. At the molecular level, many members of the ammonium transporter (AMT) and nitrate transporter (NRT) gene families have been cloned and characterized, and the uptake and utilization of NH 4 + and NO 3 − by plants is regulated by multiple genes in the N transporter protein family [4,5].
AMTs are carrier proteins that actively transport NH 4 + across biological cell membranes, and they were divided into three main subfamilies: AMT/MEP/Rh (Ammonium transporter/Methylamine permease/Rhesus protein) [6]. The AMT family can usually be divided into two subfamilies, namely AMT1 and AMT2 [7]. Both AMT1 and AMT2 show high affinity for NH 4 + , but members of the AMT1 subfamily play a more important role in high-affinity NH 4 + uptake [8]. AMT2 has a more complex gene structure and protein profile than AMT1 [9]. The AMT1 subfamily of Arabidopsis thaliana has five members, while AMT2 has low homology with the five AMT1 subfamily members and belongs to the MEP subfamily [10]. Expression of AtAMT1;1 in roots is significantly correlated with N uptake, whereas expression of AtAMT1;2 in roots is insensitive to changes in N concentration [11]. AtAMT1;4 mediate the intake of NH 4 + at the pollen plasma membrane [12]. NH 4 + is regulated into the cytoplasm through AtAMT1;1 and AtAMT2;1, which are primarily expressed in leaves [13]. In addition, AtAMT2;1 may play a role in the transport of NH 4 + from the roots to the ground [10]. A study by Couturier et al. on AMT proteins in poplar (Populus L.) found that the expression of PtrAMT1;2 was influenced by intracellular N concentration. Most PtrAMT1s were preferentially expressed in roots, while most PtrAMT2s were majorly expressed in stems, and PtrAMT3;1 was expressed only in senescing leaves [14].
The absorption systems, corresponding genes, and regulatory mechanisms of NO 3 − by plants are different from those of NH 4 + . The uptake of NO 3 by plants is the process of pumping protons out of the cell through the H + -ATPase in the plasma membrane, creating pH and electrical (∆Ψ) gradients across the plasma membrane that allows the NO 3 transporter to take up NO 3 − into the cell [15]. Four families of transporters are known to contribute to nitrate uptake and transport in plants: the nitrate transporter protein 1 (NRT1/PTR/NPF), nitrate transporter 2 (NRT2), chloride channel (CLC), and slow anion-associated channel homolog (SLC/SLAH) family [16]. The NRT1 and NRT2 are responsible for LATS and HATS, respectively [17]. Members of the NRT1 subfamily are responsible for transporting NO 3 − , hormones, glucosinolates and dipeptides [18]. In Arabidopsis, AtNPF6;3/AtNRT1;1/CHL1 was the first NRT1 subfamily member to be identified and cloned, and this protein is an amphiphilic NRT [19]. Except for AtNRT1;1, most NRT1 exhibited low affinity [20]. NPF4;6/NRT1;2 and NPF2;7/NAXT1 were also shown to be involved in root NO 3 − uptake, with NPF4.6 acting on NO 3 − influx [21] and NPF2.7 involved in NO 3 − efflux [22]. Other NRT1s are mainly related to the internal transport of NO 3 − in processes such as xylem and phloem loading and transport to leaves or seeds [23]. NRT2 belongs to the major facilitator superfamily (MFS) [24], NRT3/NAR2 is a high-affinity NRT [25]. NRT2.1, NRT2.2, NRT2.4, and NRT2.5 were known to be associated with root NO 3 − influx, and AtNRT2.1 expression was induced at low NO 3 − levels and repressed at high NO 3 − concentrations [26]. AtNRT2.4 and AtNRT2.5 are associated with root NO 3 − absorbance during severe N deficiency [27]. The NRT3 subfamily plays a significant role in NO 3 − transport by regulating the activity of NRT2, but they are not transport proteins themselves [28].
Pecan (Carya illinoinensis (Wangenh.) K. Koch) belongs to the Juglandaceae family and is one of the world's famous nut tree species, indigenous to the United States and Mexico [29]. Pecans have a good development prospect in China because of the high content of various nutrients and the important economic value of the kernels [30]. However, as an important economic tree species introduced for many years in China, pecan still has the problem of insufficient yield, and low nitrogen use efficiency (NUE) is an important factor leading to the underproduction of pecan. The key way to improve NUE is to master the N absorption and utilization pattern of plants and the molecular regulation mechanism, but there is no systematic and in-depth research on pecan in this subject. In this study, we determined the absorption characteristics of NH 4 + and NO 3 − in different N forms. We identified 10 AMT and 69 NRT gene family members in pecan and classified them into different evolutionary subfamilies. Then, we analyzed their gene structure, replication events and expression patterns to explore their functions in varying N forms, and we provide a theoretical basis for improving the NUE of pecan.

Identification and Sequence Analysis of AMT and NRT Gene Family Members in Pecan
A total of 11 AMT candidates were identified that contained AMT or AMT-like repeats, and 95 NRT candidates were identified that contained NRT or NRT-like repeats. After validation of AMT and NRT structural domains by Pfam and NCBI CD-search, 10 were identified as AMT gene family members, and 69 were identified as NRT gene family members. The AMT and NRT genes were renamed according to the Arabidopsis gene names as well as the NCBI blastp results for subsequent analysis. Table S3 provided details of AMTs and NRTs.
Sequence analysis of pecan AMT and NRT gene family members revealed that the number of exons in AMTs ranged from one to five, and in NRTs from two to nine. The CDS length of AMTs ranged from 471 to 1542 bp, and in NRTs from 1053 to 2826 bp. Most AMTs (9/10) and NRTs (56/69) were stable proteins with a low protein instability index (instability index < 40). GRAVY analysis showed that the hydration of AMT and NRT proteins in pecan was greater than 0, indicating that these proteins are hydrophobic. Subcellular localization predictions showed that most AMTs (9/10) localized to the cell membrane and most NRTs (64/69) localized to the vesicles.

Phylogenetic Analysis of AMTs and NRTs in Different Species
To decipher the evolutionary relationships and functional associations of pecan AMTs and NRTs, phylogenetic trees were constructed using pecan AMT and NRT proteins and other plants, respectively (Figures 1 and 2). According to the phylogenetic trees, all AMT proteins were divided into two distinct evolutionary branches: AMT1 and AMT2 with strong support (Bootstrap = 100%). The AMT2 was further divided into three clusters: AMT2a, AMT2b, and AMT2c. AMT1 was the largest branch and included 34 AMTs (four CiAMTs), AMT2a included two CiAMTs, AMT2b included three CiAMTs, and AMT2c had no pecan AMT gene family members.
All NRT proteins were divided into three main clades: NRT1/PTR, NRT2, and NRT3 subfamily. The NRT1 formed four subclasses, named NRT1a, NRT1b, NRT1c, and NRT1d, and included 62 CiNRTs and 52 AtNRTs. NRT2 included five CiNRTs and seven AtNRTs, while NRT3 included two CiNRTs and two AtNRTs. The evolutionary tree had 48 sister pairs, the majority of which were paralogous proteins, 38 pairs in total (21 pairs in pecan and 17 pairs in Arabidopsis), and 10 pairs of orthologous proteins. Only the NRT3 subfamily had no orthologous proteins, while all other clades contained orthologous and paralogous proteins.
only the MSF_1 structural domain; NRT3 contained only the NAR2 structural domain.
In addition, we analyzed the exons/introns of the pecan AMT and NRT genes to study the structural diversity. The results showed that the AMT1s contained no intron, and the AMT2s contained 2-4 introns ( Figure 3D); the NRT1s had 2-9 introns, the NRT2s and NRT3s contained 1-2 introns ( Figure 4D). In conclusion, the phylogenetic correlation between gene structure and prediction strongly supports a close evolutionary relationship between paired genes within the same subfamily.  . AMT proteins from seven species were divided into two subfamilies (AMT1 and AMT2). The AMT2 was further divided into three clusters (AMT2a, AMT2b, and AMT2c).

Phylogenetic Tree, Conserved Motif, Conserved Domain, and Gene Structural Analyses of CiAMTs and CiNRTs
To better understand the evolutionary relationships and functional associations of pecan AMT and NRT genes, we constructed unrooted phylogenetic trees with pecan AMT and NRT proteins, respectively ( Figures 3A and 4A). According to MEME analysis, we found 10 motifs in most of the pecan AMT and NRT proteins ( Figures 3B and 4B, Table S5). Motifs 1 to 7 were found in all AMT subfamilies, suggesting that these motifs may be characteristic motifs associated with members of the AMT gene family. Only 2-3 motifs in NRT2 were identical to the NRT1 subfamily, while NRT3 had no identical motifs to the NRT1 subfamily. We analyzed the structural domains of pecan AMT and NRT proteins using Pfam search and found that only one conserved Ammonium_trasp domain existed in all pecan AMT proteins, and they were all located at similar positions, while three conserved structural domains existed in pecan NRT proteins ( Figures 3C and 4C, Table S4).    and Arabidopsis (Arabidopsis thaliana). NRT proteins from two species were divided into three subfamilies (NRT1, NRT2 and NRT3). The NRT1 was further divided into four clusters (NRT1a, NRT1b, NRT1c, and NRT1d).  In addition, we analyzed the exons/introns of the pecan AMT and NRT genes to study the structural diversity. The results showed that the AMT1s contained no intron, and the AMT2s contained 2-4 introns ( Figure 3D); the NRT1s had 2-9 introns, the NRT2s and NRT3s contained 1-2 introns ( Figure 4D). In conclusion, the phylogenetic correlation between gene structure and prediction strongly supports a close evolutionary relationship between paired genes within the same subfamily.

Synteny Analysis of CiAMTs and CiNRTs
To determine the replication events of pecan AMT and NRT genes, we performed a synteny analysis of the pecan genome. The results showed that there were five duplicated gene pairs among the ten members of the pecan AMT gene family, two of which originated from tandem duplication and three from segmental duplication ( Figure S1). There are 103 duplicated gene pairs among 69 members of the pecan NRT gene family, of which 9 originated from tandem duplication and 94 from segmental duplication. This suggested that fragment replication events played an important role in the expansion of the AMT and NRT gene families in pecan.
To examine the selection type of duplicate gene pairs in the pecan AMT and NRT gene families, the Ka/Ks ratios were analyzed for duplication events (Table S6). Ka/Ks < 1 means the gene is subjected to purifying selection, Ka/Ks > 1 means the gene underwent positive selection, and Ka/Ks = 1 means neutral evolution. The Ka/Ks values of all duplicate gene pairs were less than 1, indicating that the amplification of the pecan AMT and NRT genes was mainly influenced by purifying selection.

Effect of N Forms on Quantitative qRT-PCR Analysis of AMT and NRT Gene Expression Levels in Pecan
The results of qRT-PCR analysis of CiAMTs showed that the relative expression levels of CiAMTs were significantly affected by different N forms ( Figure 5). The relative expression levels of almost all CiAMTs were higher in roots than in leaves, indicating that CiAMTs mainly worked in roots. In the roots, CiAMT1.1 was significantly upregulated under T3 and T4 (p < 0.05), CiAMT1.3a was significantly upregulated only under T5 (p < 0.05), and CiAMT1.3b, CiAMT1.4 and CiAMT3.1a were all significantly upregulated only under T4 (p < 0.05). CiAMT2.1 was significantly upregulated under T1, T3, and T4 (p < 0.05), CiAMT2.2 was significantly upregulated under T4 and T5 (p < 0.05), and CiAMT3.1b was significantly upregulated under T1, T2, and T3 (p < 0.05   Based on the results of previous studies [16], some members of the Arabidopsis and rice NRT gene families have been shown to be associated with NO 3 − uptake and transport. According to the phylogenetic tree, we selected 16 pecan CiPawNRTs corresponding to them and further analyzed them using qRT-PCR ( Figure 6). The qRT-PCR results showed that CiNRTs showed different expression patterns in different N forms as well as in different pecan organs. Except for CiNPF2.13, CiNPF4.6, CiNPF5.5, CiNPF6.3a, CiNPF6.4, CiNRT2.5, and CiNRT2.7, which showed relative higher expression in leaves than in roots under most N form treatments, most of the other CiNRTs showed higher expression in roots than in leaves. In leaves, CiNPF2. 4  To further investigate the response mechanisms of pecan AMT and NRT genes to N forms, we measured the concentrations of NH 4 + , NO 3 − , and NO 2 − in pecan under different N forms (Figure 7). The results showed that there was no significant difference in the NH 4 + concentrations of pecan in all organs under different N forms. Except for the T5 treatment, NH 4 + concentrations under all other treatments showed greater in roots than in leaves and stem (p < 0.05) and no significant difference between leaves and stems. The variability of NO 3 − concentrations in each organ of pecan varied among treatments, and T5 showed significantly greater than CK and T3 in leaves (p < 0.05), and no significant differences were found between T1, T2, T4, and other treatments. In the stems, no significant differences were found between treatments. In the roots, it showed that T4 was significantly greater than CK (p < 0.05), and T4 was not significantly different from the other treatments. Except for T3 and T5 treatment, the NO 3 − concentrations of pecan under all other treatments showed greater leaves than roots (p < 0.05) and no significant difference between stems, leaves, and roots. There was no significant difference in NO 2 − concentrations in all organs of pecan under different treatments, while the variability of NO 2 − concentrations in different organs under each treatment varied. CK, T1, and T5 showed no significant difference between stems and roots, and both of them were significantly greater than leaves (p < 0.05), T2 and T3 showed no significant difference between leaves and roots, and both of them were significantly greater than stems (p < 0.05), while T4 showed no significant variation in different organs. Figure 6. Relative expression of selected pecan CiNRT genes under varying NH4 + :NO3 − ratios. The expression levels of selected CiNRT in pecan leaves and roots after varying NH4 + :NO3 − ratio treatments were quantified by qRT-PCR, with Actin as the reference gene. Different capital letters indicate significant differences in leaves (p < 0.05), and different lowercase letters indicate significant differences in roots (p < 0.05).

Effect of N Forms on NH4 + , NO3 − and NO2 − Concentration in Pecan
To further investigate the response mechanisms of pecan AMT and NRT genes to N forms, we measured the concentrations of NH4 + , NO3 − , and NO2 − in pecan under different N forms (Figure 7). The results showed that there was no significant difference in the NH4 + concentrations of pecan in all organs under different N forms. Except for the T5 treatment, NH4 + concentrations under all other treatments showed greater in roots than in leaves and stem (p < 0.05) and no significant difference between leaves and stems. The variability of NO3 − concentrations in each organ of pecan varied among treatments, and T5 showed significantly greater than CK and T3 in leaves (p < 0.05), and no significant differences were found between T1, T2, T4, and other treatments. In the stems, no significant differences were found between treatments. In the roots, it showed that T4 was significantly greater than CK (p < 0.05), and T4 was not significantly different from the other treatments. Except for T3 and T5 treatment, the NO3 − concentrations of pecan under all other treatments showed greater leaves than roots (p < 0.05) and no significant difference between stems, leaves, and roots. There was no significant difference in NO2 − concentrations in all organs The expression levels of selected CiNRT in pecan leaves and roots after varying NH 4 + :NO 3 − ratio treatments were quantified by qRT-PCR, with Actin as the reference gene. Different capital letters indicate significant differences in leaves (p < 0.05), and different lowercase letters indicate significant differences in roots (p < 0.05).
At the mean level, pecan NH 4 + concentrations were significantly greater in the T4 than in the other treatments (p < 0.05), with no significant differences between the other treatments. The pecan NO 3 concentrations revealed no significant difference between T1, T4, and T5, and both of them were significantly greater than CK, T2, and T3 (p < 0.05). There was no significant difference between T2 and T1, T3 and T4, T2 was significantly greater than CK (p < 0.05), and no significant difference between CK and T3. The variability of pecan NO 2 and NH 4 + concentrations was consistent, suggesting that T4 was significantly greater than the other treatments (p < 0.05), with no significant differences between the other treatments.

Effect of N Forms on the Uptake Kinetics of NH 4 + and NO 3 − in Pecan
We also determined the kinetic properties of pecan NH 4 + and NO 3 − uptake under different N forms (Figure 8). The results showed that the uptake rates of NH 4 + and NO 3 − under T2, T3, and T4 were still in a significant upward trend at the ion concentration of 2000 µmol·L −1 ; the uptake rates of NH 4 + and NO 3 − under CK, T1, and T5 leveled off at the medium ion concentration of 1000 µmol·L −1 .
ity of pecan NO2and NH4 + concentrations was consistent, suggesting that T4 was significantly greater than the other treatments (p < 0.05), with no significant differences between the other treatments.

Effect of N Forms on the Uptake Kinetics of NH4 + and NO3 − in Pecan
We also determined the kinetic properties of pecan NH4 + and NO3 − uptake under different N forms (Figure 8). The results showed that the uptake rates of NH4 + and NO3 − under T2, T3, and T4 were still in a significant upward trend at the ion concentration of 2000 μmol·L −1 ; the uptake rates of NH4 + and NO3 − under CK, T1, and T5 leveled off at the medium ion concentration of 1000 μmol·L −1 . The root uptake rates were processed data according to the Hofstee transformation equation to obtain the maximum uptake rate (Vm) and the Mee's constant (Km) of pecan for different N forms with significant coefficients of determination R 2 (Table 1). Different N forms showed significant effects on Vmax and Km of NH4 + and NO3 − (p < 0.05), both showing T4 > T3 > T2 > CK > T1 > T5. This indicates that the involvement of a certain proportion of NH4 + accelerated the uptake of NH4 + and NO3 − by pecans, but the affinity decreased, while NH4 + above a certain proportion decreased the uptake rate and increased the affinity. Except for the T1 treatment, the Vmax/Km of NH4 + was greater than that of NO3 − under all treatments, indicating that the rate of NH4 + uptake by pecan was greater than that of NO3 − . The root uptake rates were processed data according to the Hofstee transformation equation to obtain the maximum uptake rate (Vm) and the Mee's constant (Km) of pecan for different N forms with significant coefficients of determination R 2 (Table 1). Different N forms showed significant effects on Vmax and Km of NH 4 + and NO 3 − (p < 0.05), both showing T4 > T3 > T2 > CK > T1 > T5. This indicates that the involvement of a certain proportion of NH 4 + accelerated the uptake of NH 4 + and NO 3 − by pecans, but the affinity decreased, while NH 4 + above a certain proportion decreased the uptake rate and increased the affinity. Except for the T1 treatment, the Vmax/Km of NH 4 + was greater than that of NO 3 − under all treatments, indicating that the rate of NH 4 + uptake by pecan was greater than that of NO 3 − .

Functional Differentiation of AMT and NRT Gene Family in Pecan Genome
The biological functions of AMT and NRT gene families in pecan are poorly understood. Therefore, we report for an earlier time identification of 10 AMT and 69 NRT gene family members from pecan, confirming that direct orthologs of AMT and NRT proteins should be highly conserved evolutionarily throughout the plant.
Studies on Arabidopsis AMT proteins have shown that AtAMTs had a prominent role in NH 4 + assimilation at the cell membranes [10,11,31], and studies of flowering Chinese cabbage (Brassica campestris) also indicated that BcAMT2 was located on the plasma membrane [32], which was consistent with the predicted results of subcellular localization in this study (Table S3), and this is the optimal cellular structure for maintaining stable NH 4 + concentrations in plants. AtNPF5.11, AtNPF5.12, and AtNPF5.16 were identified to function in the process of NO 3 − from the vesicle to the cytoplasm, thereby regulating NO 3 − distribution between roots and shoots [33]. In contrast, the predicted subcellular localization of the pecan NRT genes in this study showed that most of the NRT1s were localized to the vacuoles (Table S3), suggesting that pecan NRT1s may mostly act in the transport of NO 3 − between the vacuoles and the cytoplasm. There were also some NRT1s localized on the cell membrane, while all NRT2s were all localized to the cell membrane, suggesting that this fraction of NRT proteins may mediate assimilation and efflux of NO 3 − in plants as well as inter-subcellular transport, which was consistent with the findings of cassava (Manihot esculenta) [34]. According to the predicted results of subcellular localization, the subcellular localization of all NRT3s differed significantly from NRT1 and NRT2 in pecan, with CiNRT3.1 localized to the nucleus and CiNRT3.2 localized on the cell membrane, cell wall, chloroplast, and vacuoles (Table S3), indicating that the pecan NRT3s had more complex structures and functions.

Effect of N Forms on the Absorption Characteristics of NH 4 + and NO 3 − in Pecan
NH 4 + and NO 3 − are the two main forms of N absorbed and utilized by plants, and previous studies have shown that a nutrient mixture of NH 4 + and NO 3 − can improve crop yield and quality compared to a single source of N [35]. Therefore, understanding the best NH 4 + :NO 3 − ratios provides the possibility to improve the N utilization efficiency of plants. We investigated the effect of varying NH 4 + :NO 3 − ratios on the absorption and transport of N by measuring the concentrations of different N forms in the pecan organs (Figure 7).
This study showed that the concentrations of different N forms in pecan were tissuespecific. The NH 4 + concentrations of pecan were significantly higher in roots than in leaves and stems, while NO 3 − and NO 2 − concentrations were significantly higher in leaves than in stems and roots, which was consistent with the finding that NH 4 + was majorly assimilated in roots and NO 3 − was mostly translocated to leaves for storage or assimilation [36]. T4 significantly increased the total NH 4 + concentrations of pecan, primarily in the roots, suggesting that T4 was more favorable to promote NH 4 + absorbed by pecan roots, but had no effect on NH 4 + transport between organs. All NH 4 + :NO 3 − ratio treatments increased the total NO 3 − concentrations, primarily in the leaves, indicating that the feeding of NH 4 + and NO 3 − mainly promoted the translocation of pecan NO 3 − from roots to leaves, and possibly the acclimation of NO 3 − in the leaves. Compared to T1, T3 reduced the total pecan NO 3 − concentrations, which was in agreement with the results of studies in pepper (Capsicum annuum) [37]. T4 significantly increased the total NO 2 − concentrations of pecan, indicating that T4 promoted the conversion of NO 3 − to NO 2 − . Previous conclusions on the effects of the simultaneous presence of NH 4 + and NO 3 − on each other's uptake were varied, with some suggesting that they have a facilitative or inhibitory effect [38,39], and that plant absorption of NH 4 + and NO 3 − is limited by the maximum uptake threshold [40]. In this study, the uptake rates of NH 4 + and NO 3 − gradually saturated with increasing substrate concentrations under a single N source, while they remained on an increasing trend under mixed N source treatments, with the T4 treatment being the most obvious, indicating that NH 4 + and NO 3 − promoted each other's intake in this study, as evidenced by the magnitude of Vmax values of pecan. However, single N source treatment increased the affinity of NH 4 + and NO 3 − , and the combined application of NH 4 + and NO 3 − decreased their affinity instead, which was generally consistent with the results of Kamminga-Van Wijk and Prins [41]. Vmax/Km is also commonly used to indicate plant preference for NH 4 + and NO 3 − uptake [42], suggesting that pecans may be more biased toward NH 4 + absorption.

Effect of N Forms on the Expression of CiAMTs and CiNRTs
The ingestion and transport of NH 4 + and NO 3 − in plants is majorly mediated by AMT and NRT gene family members. Therefore, we investigated the expression modes of these genes in varying organs and N forms to further understand the effect of N forms on N uptake and transport in pecan.
In the present study, the relative expression levels of almost all CiAMTs were higher in roots than in leaves, indicating that the addition of NH 4 + and NO 3 − majorly stimulated the expression of CiAMTs in roots. Studies have shown that AtAMT1.1, AtAMT1.2, AtAMT1.3, and AtAMT1.5 are the principal transporter proteins that take up high-affinity NH 4 + into Arabidopsis roots, with AtAMT1.1 and AtAMT1.3 responsible for approximately two-thirds of the high-affinity NH 4 + uptake capacity in the root; the CiAMT1s may have the same function [31,43]. The results of this study showed that all CiAMT1s were upregulated in roots of T4, except for CiAMT1.3a, which was upregulated in roots under T5, suggesting that T4 may have improved the absorption of high-affinity NH 4 + by the pecan root. In this study, CiAMT1s were mostly upregulated in leaves under N lack, but numerous studies have shown that expression levels of AMT1s were upregulated in roots under N limited conditions [14,43,44]. This may be due to the variation caused by the longer duration of N deficiency in this study, or it is possible that the expression of CiAMTs in leaves was not only affected by N lack, but may also be involved in some other regulatory mechanisms [45]. Previous studies have shown that peach (Prunus persica) PpeAMT3;4 was primarily expressed in roots [46], and the expression pattern of CiAMT2s was similar to that of CiAMT1s, being expressed mostly in roots and almost significantly upregulated at all timepoints under T4.
Studies on Arabidopsis suggested that AtNPF6.3/AtNRT1.1 was not only an amphiphilic NO 3 − transport protein but might also act as a NO 3 − sensor under low NO 3 − conditions [47]. The expression level of CiNPF6.3a was significantly upregulated in roots under T4, suggesting that it may carry both translocation and signaling functions at this time. AtNPF4.6/AtNRT1.2 was not only a low-affinity NO 3 − transporter protein that plays a role in NO 3 − influx [21], but also an abscisic acid (ABA) transporter protein that positively regulated the ABA response [48]. CiNPF4.6 expression was significantly upregulated in roots under T3 and in leaves under T5, perhaps because CiNPF4.6 worked primarily on NO 3 − influx under T3, while it worked mostly on ABA regulation under T5. MtNPF6.8/MtNRT1.3 of Medicago truncatula was an amphipathic NO 3 − transport protein and was upregulated by the absence of NO 3 − [49]. This was the same as the results of our study, where CiNPF6.4/CiNRT1.3 expression was significantly upregulated in roots under CK and in leaves under T5. AtNPF7.3/AtNRT1.5 mediated root-to-stem transport, and the same conclusion was found for ZxNPF7.3/ZxNRT1.5 in the Zygophyllum xanthoxylum, which also contributed to the uptake of NO 3 − [50]. The expression levels of CiNPF7.3a/CiNRT1.3a were significantly upregulated in roots, leaves under CK, and roots under T4, indicating that both N scarcity and T4 may promote NO 3 − transport from roots to stems. Moreover, AtNPF7.2/AtNRT1.8 was phylogenetically similar to AtNPF7.3 [51], while CiNPF7.2/CiNRT1.8 was also significantly upregulated in roots under T4 and may function similarly to CiNPF7.3. NPF2.13/NRT1.7 and NPF1.2/NRT1.11 were proven to be involved in the transfer and redistribution of NO 3 − from xylem or the NO 3 − containing tissues to the phloem [23], CiNPF2.13 was largely induced by N limitation and T3, while CiNPF1.2a was mainly induced by N limitation, and T4. NPF2.11/NRT1.10 was identified to be involved in thioglucoside transport [52], and CiNPF2.11 expression was significantly upregulated under CK, suggesting that CiNPF2.11 may resist N deficiency stress by regulating the concentrations of thioglucosides. AtNRT2.5 was a plasma-membrane localized high-affinity NO 3 − transporter protein that mediated NO 3 − acquisition, and reactivation under N deficient conditions [27], and CiNRT2.5 in this study exhibited the same expression pattern. The transcript levels of NRT2.7 in Fraxinus mandshurica were both upregulated in leaves due to N limitation [53], whereas CiNRT2.7 in this study was significantly expressed in roots under CK and in leaves under T1, which may be caused by species differences.

Plant Materials and Experimental Design
In the study, aeroponic cultivation trials were carried out in the greenhouse of the campus of Nanjing Forestry University from 18 April to 9 June 2021 and from 4 May to 21 June 2022. The plant materials and experimental design refer to Chen et al. [29]. In the case of the same N supply, the five ammonia-to-nitrate ratios (NH 4 + :NO 3 − ) were 100:0, 75:25, 50:50, 25:75, and 0:100, corresponding to T1, T2, T3, T4, and T5, respectively. The nutrient solution without N was used as the control (CK), and each treatment was repeated 3 times, each with 6 seedlings. Regulation of the NH 4 + :NO 3 − ratios for each treatment was achieved with specific source compounds (Table S1). Samples were taken after 45 days of treatment for further determination.

Identification of AMT and NRT Genes in Pecan
To identify the pecan AMT and NRT genes, we obtained all of the protein sequences of pecan from the Phytozome v13 database (https://phytozome-next.jgi.doe.gov/info/ CillinoinensisPawnee_v1_1 (accessed on 1 November 2021)). The hidden Markov model (HMM) profiles of the AMT domain (PF00909) and NRT domain (PF07690, PF00854, PF16974), downloaded from the Pfam database (http://pfam.xfam.org/ (accessed on 1 February 2022)) [54], were used to identify the pecan AMT and NRT proteins by using HMM search through HMMER3.0 program (www.hmmer.org (accessed on 1 February 2022)) with default parameters [55]. As multiple AMT and NRT proteins were corresponding to a specific gene in several cases, only one protein sequence corresponding to each gene was retained for further detailed analysis. The Pfam database was used to confirm the presence of these conserved domains of the screened genes. The biophysical properties such as amino acid length (AA), molecular weights (MWs), theoretical isoelectric points (pIs), and grand average of hydration (GRAVY) of pecan AMT and NRT proteins were estimated by ExPASy ProtParam server (http://web.expasy.org/protparam (accessed on 1 February 2022)) [56]. Transmembrane helices (TMHs) were determined using the TMHMM tools (http://www.cbs.dtu.dk/services/TMHMM-2.0/ (accessed on 1 February 2022)), prediction of subcellular localization information for pecan AMT and NRT proteins using the Cell-PLoc 2.0 software (http://www.csbio.sjtu.edu.cn/bioinf/Cell-PLoc-2/ (accessed on 1 February 2022)) [57].

Phylogenetic Analysis
The Arabidopsis AMT and NRT protein sequences were downloaded from the TAIR database (https://www.arabidopsis.org/ (accessed on 1 November 2021)) [58]. The poplar, apple, pear, tomato, and rice AMT protein sequences were downloaded from the Phytozome database. Phylogenetic trees were constructed by MEGA 7.0 with the following settings: the neighbor-joining (NJ) method, 1000 bootstrap replicates, the Jones-Taylor-Thornton (JTT) model, and pairwise deletion [59]. The evolutionary tree was visualized using the online tool Evolview (http://evolgenius.info/ (accessed on 1 March 2022)) [60]. Because CiAMT3.1c is a partial sequence, it was not included in the phylogenetic tree.

Analysis of Gene Structure, Conservative Motifs, and Domains
The pecan genome sequence files and the General Feature Format (GFF) annotation files were downloaded from the Phytozome database. We used the online MEME tool (http://MEME-suite.org/ (accessed on 1 April 2022)) for topic prediction, keeping the maximum number of topics at 10 [61]. TBtools visualizes AMT and NRT gene structures, conservative motifs, and domains [62].

Synteny Analysis
For detecting syntenic blocks, the whole genome sequence file and the GFF annotation file of pecan were used to identify all duplication events in the pecan genome using MCScanX software [63]. Then, the AMT and NRT gene families were analyzed for synteny and visualized using TBtools.

Estimation of the Ka/Ks Values
Multiple sequence alignment of full-length coding sequences (CDS) in the AMT and NRT gene families in pecan was performed using MEGA 7 software and further used to calculate nonsynonymous (Ka) and synonymous (Ks), and the Ka/Ks ratios. Ks values were commonly used to determine the time since gene duplication, and the selection pressure of duplication event was determined by the Ka/Ks ratio.

Cis-Regulatory Elements Analysis
To investigate potential cis-regulatory elements in the promoters of the pecan AMT and NRT genes, a 1000 bp region upstream of the AMT and NRT genes was retrieved from the pecan genome sequences. Then, the cis-regulatory elements were predicted using PlantCARE software (//bioinformatics.psb.ugent.be/webtools/plantcare/html/ (accessed on 1 April 2022)) [64] and screened and visualized by TBtools software (Figures S2 and S3).

RNA Collection and qRT-PCR Expression Analysis
According to the manufacturer's protocol, the total RNA was extracted from the leaves and roots of pecan using a Universal Plant Total RNA Extraction Kit (Bioteke, Beijing, China) and stored at −80 • C until further use. The purity and integrity of the isolated total RNA were analyzed by agarose gel electrophoresis and Nanodrop 2000 spectrophotometer (Thermo Scientific, Wilmington, NC, USA). First-strand cDNA was synthesized using a cDNA Synthesis Kit (HiScript ® RIII RT SuperMix for qPCR +gDNA wiper, Vazyme, Nanjing, China). The qRT-PCR was performed on a 7500 Real-Time PCR system (Applied Biosystems TM , Foster City, CA, USA) using a Taq Pro Universal SYBR qPCR Master Mix (Vazyme, Nanjing, China). The specific primers were synthesized by Tsingke Biotechnology Ltd. (Nanjing, China), and the details of the primers were provided in Supplementary  Table S2. The PCR parameters applied here were as follows: 95 • C for 30 s, followed by 40 cycles of 5 s at 95 and 30 s at 60 • C. The Actin gene was used as an internal reference gene [66], and the relative expression levels of pecan AMT and NRT genes were determined using the 2 −∆∆Ct method [67]. Values represent mean calculated from three biological replicates and three technological repeats.

Measurement of NH 4 + , NO 3 − and NO 2 − Concentration in Pecan
The NH 4 + concentrations of pecan roots, stems and leaves were determined according to the Berthelot reaction [68]. NO 3 − concentrations were determined according to the method suggested by Patterson et al. [69]. The NO 2 − concentrations were determined by the method of Ogawa et al. [70].

Kinetic Characterization of NH 4 + and NO 3 − Uptake in Pecan
The kinetic characteristics of NH 4 + and NO 3 − uptake in pecan seedlings were determined by the conventional depletion method, the ion concentrations in the solutions to be measured were determined after preparing different concentrations of NH 4 + and NO 3 − for 24 h incubation of the plants, the net rate of ion uptake per unit fresh root per unit time was calculated, and the kinetic parameters of uptake were mathematically derived. The concentration of NH 4 + and NO 3 − was determined by referring to 2.10.

Data Analysis
Before analysis of variance (ANOVA), data were checked for normality and homogeneity of variances. One-way ANOVA was performed to test the effects of different N forms on NH 4 + , NO 3 − , NO 2 − concentrations, absorption kinetic characteristics and relative gene expression of pecan seedlings. Differences were considered significant at p < 0.05.
The kinetic parameters of NH 4 + and NO 3 − uptake were calculated using the Hofstee transformation of the Michaelis-Menten kinetic equation: V = C × Vmax/(Km + C), C represents the ion concentration, V indicates the net ion uptake rate, Vm indicates the maximum uptake rate, and the Km value represents the root uptake site for ion affinity. Non-linear regression fitting and graphing were performed using SPSS to obtain the Vmax and Km. The α value represents the competitive ability of the plant root system for nutrient uptake, α = Vmax/Km. All statistical analyses were performed with SPSS 23.0 software (Version 23.0, Chicago, IL, USA). All charts were drawn with Excel (Version 2019, Redmond, WA, USA) and SigmaPlot (Version 14.0, Barcelona, Spain).

Conclusions
We identified 10 AMT and 69 NRT genes in the pecan genome, and the analysis showed that the biophysical properties, gene structure, and expression levels of CiAMTs and CiNRTs were strongly associated with pecan NH 4 + and NO 3 − uptake. Combining the effects of different N form treatments on the expression levels of CiAMTs and CiNRTs, the N concentrations of pecan, and the uptake rate of NH 4 + and NO 3 − , we concluded that pecan preferred NH 4 + and that the NH 4 + /NO 3 − ratio of 75:25 was more favorable to improve the N uptake capacity of pecan seedlings. This study provides a basis for further identification of the functions of AMT and NRT genes in N uptake and transport of pecan, and it provides a theoretical basis for the application of an optimal proportion of N fertilizer to improve NUE, thereby increasing pecan yield and promoting pecan industrialization.