Interspecies Evolution and Networks Investigation of the Auxin Response Protein (AUX/IAA) Family Reveals the Adaptation Mechanisms of Halophytes Crops in Nitrogen Starvation Agroecological Environments

: The maintenance of adaptability to the exposure to agroecological extreme environments is generally a feature after the long-term domestication of crops. Auxin inﬂuences plant growth in all environments. At present, the research on the auxin response factors ( ARFs ) has been very in-depth. However, there is still a large gap in the research on the origin, evolution, and regulatory networks of the Auxin-responsive protein (AUX/IAA) family. Here, we identiﬁed 495 AUX/IAAs from 19 representative species covering aquatic algae to angiosperms and found that they originated from early bryophytes and mainly expanded by polyploidy in angiosperms. In the domesticated crop quinoa, the evolutionary model of the IAA family is relatively independent and forms a robust regulatory network with auxin signals and energy metabolism pathways. In the nitrogen-deﬁcient environment, CqIAAs ( Chenopodium quinoa AUX/IAAs), auxin signals, and TCA pathway genes induced expression in young roots to promote root elongation and could regulate the balance of carbon and nitrogen metabolism to maintain the adaptation of early seedlings in poor environments. Furthermore, a rapidly evolving CqIAA ( AUR62011942 ) not only has different expression levels in two quinoa seeds but also has a signiﬁcant stress response when seedlings face nitrogen deﬁciency stress, which may be a key factor in the adaptive regulation of the barren environment. Our research provides valuable clues for understanding the origin, evolution, and functional innovation of auxin signaling and also provides a reference for future agricultural breeding in the context of global environmental changes. Multiple collinearity scan toolkit (MCScanX) was used to detect duplication types to identify GRASs associated with ancestral polyploidization events (WGD or segmental). The AUX/IAAs related to polyploidy were inferred according to the collinearity persistence of the region around the selected genes. AUX/IAAs were divided into ﬁve groups, including dispersed, proximal, tandem, singleton, and WGD or segmental.


Introduction
In the process of plant evolution from water to land, it experienced complex and violent environmental changes and formed the survival mechanism regulated by various complex environmental signals through singular evolutionary strategies [1]. How land plants adapt and respond to harsh stress without water-surrounded environments is an intriguing and fundamental biological question. Moreover, improving plant resistance is essential for agricultural productivity and environmental sustainability. To survive in a variety of environmental stresses, plants regulate multiple metabolic pathways and have evolved a variety of adaptive characteristics governed by complex signals and regulatory systems, as well as downstream networks controlled by transcription factors (TFs) [2,3]. The auxin signaling pathway is involved in the most basic biological process of plant growth and development, which determines the fate of plants from birth to death [4,5]. In mainly from 400 to 510 nm and from 610 to 720 nm. The hydroponic conditions were set as two groups: the control group was treated with Hoagland complete nutrient solution, and the nitrogen-deficient Hoagland formula was used in the nitrogen-deficient treatment group. We replaced the culture medium every three days. Two weeks after the growth of seedlings, the roots and leaves of quinoa were collected and stored in a refrigerator at −80 • C for subsequent physiological and qPCR experiments. The Arabidopsis genome data were obtained from the TAIR database (https://www. arabidopsis.org/ (accessed on 24 January 2021)), the wheat genome data from the Ensemble database (http://ensemblgenomes.org/, accessed on 24 January 2021), and the genome data of other species including quinoa were downloaded from the Phytozome plant genome database (https://phytozome.jgi.doe.gov/, accessed on 24 January 2021). The Uniport database (https://www.ebi.ac.uk/uniprot/, accessed on 24 January 2021) was used to query all the plant AUX/IAA proteins and obtain their protein sequence files. The AUX/IAA protein sequence was compared with the multispecies proteome data by BLAST method [16]. To cope with the potential issue of low sequence similarity of homologous sequences from distally related species, the hidden Markov model (HMM) profile of the AUX/IAA domain (PF02309) was downloaded from the Pfam database (https://pfam.xfam.org/, accessed on 24 January 2021). The HMM query in the HM-MER v3.1 program was used to retrieve all AUX/IAAs from the 19 genomes with the threshold of an e-value ≤ 10 −5 . All results were submitted to the NCBI-BLAST tool (https://blast.ncbi.nlm.nih.gov/Blast.cgi, accessed on 24 January 2021) for online retrieval to remove entries that did not belong to the family. The protein sequences of lower plants and higher plants' IAA were submitted to the orthovenn2 online tool [17] for orthologous identification. The maximum likelihood (ML) tree of multiple IAA proteins was conducted using the best-fitting in MEGA7, with Arabidopsis AUX/IAA sequences as outgroups. The optimal amino acid substitution model was deduced by the model selection (ML) program in MEGA: retain "JTT with Freqs (JTT+F)" for "Model/Method" and retain "Gamma distributed with sites (G+I)" for "Rates among sites". The default settings for other parameters were retained. The support values were estimated by 500 nonparametric bootstrap iterations. The IAA amino acid sequences of all species from different orthologous were searched in Pfam and CDD databases to find other known domains/motifs in addition to the AUX/IAA domain. In addition, MEME v4.9.0 [18] was used to identify the conserved motifs with the following parameters: width of a motif was 10-50 amino acids, and the number of motifs was 10. In addition, we identify key genes (YUC, TAA, TIR, ARF, PDC, IDH, OGDC, etc.) for quinoa auxin signal and tricarboxylic acid cycle pathway via BLASTp and online database annotation from quinoa genome annotation, combined with sequences of Arabidopsis related genes in the TAIR database. Details of these genes are recorded in Table S12.

Chromosome Distribution, Gene Structure, and Syntenic Analysis in Quinoa and Other Plants
Gff files and sequencing files were used to obtain chromosome localization and gene structure information of CqIAAs and to visualize data through TBtools [19]. To identify homologous blocks in different species, the syntenic analysis maps between the CqIAAs and AUX/IAAs from other plants were constructed with Dual Synteny Plotter software (https://github.com/wyp1125/MCScanx, accessed on 24 January 2021) [20].

Identification of Collinearity and Specific Duplication Events
According to previous studies, gene collinearity analysis was performed using default parameters in MCScanX [21,22]. We used the Duplicate Gene Classifier program in MC-ScanX to speculate on different types of duplicated genes. We found potential anchoring points (E-value < 10 −5 ; top 5 matches) between each possible pair of chromosomes in multiple genomes using BlastP. The loose E-value threshold accommodates a high differentiated Agriculture 2021, 11, 780 4 of 22 evolution of duplicated genes resulting from polyploidization millions of years ago. Each species protein sequences are searched for their own or other species genome sequences. The homologous blocks in each genome and those between the genomes of different species were determined with Colinearscan (maximal gap ≤ 50 genes; p-value < 0.05) [23]. With reference to previous studies, the key parameters for the maximum gap between adjacent genes on chromosomes showing collinearity with genes on the corresponding chromosome sequence was set to 50 intervening genes [24].
2.6. Trend Analysis of Gene Expression STEM v1.3.7 [27] were used to analyze the transcriptome data of different quinoa tissues based on the FPKM values under default parameters. a total of 49 expression pattern matrices were obtained, of which 14 were statistically significant trend modules. Detailed expression data and Profiles information are recorded in Table S8.

Correlation Network Construction of CqIAAs and Other Metabolic Pathways
Referring to previously reported methods, the interaction networks of CqIAAs with Gene of auxin biosynthesis, auxin signaling and TCA cycle pathways genes, auxin signaling, photosynthesis, and flavonoids under all different stresses were constructed by using the R package imsbInfer in GitHub. The associations were calculated using the Pearson correlation measure, and the red and blue edges represent positive and negative correlations.

Calculation of the Synonymous and Nonsynonymous
Using KaKs Calculator 2.0 to calculate the synonymous (Ks) and nonsynonymous (Ka) substitutions of each gene pair and their ratios (Ka/Ks) [28].

Gene Expression Analysis in Two Kinds of Quinoa
Publicly available C. quinoa RNA-Seq datasets were obtained from previous reports and the SRA (Sequence Read Archive, Shenzhen, China, https://www.ncbi.nlm.nih.gov/ sra, accessed on 24 January 2021). These data were analyzed for the expression patterns of the identified AUX/IAAs and related metabolic pathway genes. Materials from fruit (SRS2464890), flower (SRS2464889), leaf (SRS4026093, SRS4026095), stem (SRS2464891), root [29], and materials with different fruit sizes [30] were used to examine the expression patterns and regulatory network. The expression level was normalized, and log2transformed was used for producing the heat map using the heatmap in TBtools. Using R packet limma (https://github.com/wolski/imsbInfer, accessed on 24 January 2021) to analyze the differential expression of RNA-seq data, Foldchange > 2 and p-value < 0.05 were selected as the screening thresholds of significantly different genes. The heat maps and volcanic maps were also drawn by TBtools (v1.09854).

Detection of Physiological Indexes after Stress Treatment
The samples were divided into four groups: large-fruit quinoa control (LF + N), small fruit quinoa control (SF + N), large-fruit quinoa nitrogen deficiency treatment (LF − N), and small-fruit quinoa nitrogen deficiency treatment (SF − N). Each group had 15 plants, and 5 plants with the closest growth were randomly selected as biological repeats, and the average of 5 technical repeats was selected as the final result and recorded in Table S14.
The determination of the physiological parameters mainly refers to the method of Zhang et al. [31]; the specific determination method is as follows: Fresh leaves of 0.5 g were submerged into 90% acetone in an avoiding light environment for 24 h (4 • C). The absorbance of the solution was measured at 646.6 and 663.6 nm with a spectrophotometer (AOE Instruments, Shanghai, China), and then the calculations of chlorophyll content were made with reference to previous reports [31]. The amount of xylogen in the sample was determined by measuring the absorbance of reactants at 280 nm via perchloric acid oxidation method using a xylogen content assay kit AKSU010M (BOXBIO, Beijing, China); using a total pectin content assay kit AKSU069M (BOXBIO, Beijing, China); through the carbazole ratio method, the sample reacted with the carbazole reagent under concentrated sulfuric acid, and the absorbance of the generation was determined at 530 nm to determine the content of the pectin. The fresh leaf or root tissues of different treatments were homogenized with potassium phosphate (PBS) buffer (pH = 7.4) under 0 • C precooling conditions. The homogenate was centrifuged at 20,000× g for 20 min at 4 • C, and the supernatant was collected for testing. Determination of the physiological parameters of plant stress was tested via the kit produced by Nanjing Jiancheng Institute of Bioengineering, taking malondialdehyde (MDA) content as assessment index of lipid peroxidation, and MDA content was determined by the thiobarbituric acid (TBA) method though recording the absorbance at 600, 532, and 450 nm; superoxide dismutase (SOD) activity was determined by nitroblue tetrazolium; catalase (CAT) activity was determined by UV spectrometry method; peroxidase (POD) activity was determined by guaiacol method; hydrogen peroxide (H 2 O 2) content was quantitated by the hydroxylamine method; and the content of ascorbic acid (Vc) was determined by iron reduction/antioxidant and ascorbic acid (FRASC) method [31].

Detection of Gene Expression after Stress Treatment
We determined the expression level of CqIAA and related pathway genes (YUC, TAA, TIR, ARF, IAA, PDC, IDH, OGDC, and MDH) before and after nitrogen deficiency stress in two quinoa young roots by qRT-PCR. Primer 3 (http://frodo.wi.mit.edu/, accessed on 24 January 2021) was used to design primers (Table S16). The elongation factor 1 alpha (EF1α) gene was used as an internal reference gene [32]. In our study, five seedlings with similar growth were selected from each treatment as a biological replicate. Total RNA of two quinoa (LF and SF) roots before and after nitrogen starvation was extracted using RNA Easy Fast DP452 kit (Tiangen, Beijing, China). The cDNA was synthesized with 1 mg RNA samples using a PrimeScript RT reagent kit (TaKaRa, Dalian, China) with gDNA Eraser (TaKaRa) and SYBR Premix Ex Taq II (TaKaRa). The mRNA expression level in four samples (LF + N, SF + N, LF − N, and SF − N) was detected via qPCR experiments and obtained acquired gene expression information by the 2 −∆∆Ct method [33]. During the experiment, each biological replicate was set up with three technical repeats. The mean value of the technical repeats was taken as the final result, and the three biological repeats were statistically tested.

Statistics
All data were analyzed using the Origin Pro 2020 statistics program, and the means were compared by the least significant difference test at the 0.05 level of significance. Error analysis is based on data from 3 biological replicates and 3 technical replicates.

AUX/IAA Family Originated in the Early Stage of Plant Landification and with Conservative Evolutionary Pattern in Higher Plants
Integrating hormone signaling from different plants during landification is complicated because of intense environmental changes. This longstanding problem is even more extreme in terrestrial plants because of their often complex and polyploid genomes as well as past whole-genome duplication events [34]. To further enhance our understanding of the evolution of hormone signals during the process of plant landification, a total of 495 auxin-reactive proteins (AUX/IAA) were identified from 19 plants of agricultural and evolutionary importance. AUX/IAA proteins were not found in two aquatic plants (Volvox carteri and Chlamydomonas reinhardtii) ( Figure 1A,B, Table S1), and this suggests that AUX/IAA protein family may have originated after plants have colonized the land. We also predicted the physical properties of the identified AUX/IAA proteins, and the molecular weights and isoelectric points of different AUX/IAA proteins were different ( Figure 1A, Table S1). Further, we searched for direct homologous proteins from AUX/IAA proteins of lower and higher plants, respectively. According to the conserved ancestral copies shared among different plants, the AUX/IAA proteins in lower plants were divided into nine orthogroups ( Figure 1B). The lower plant orthogroup 01 (LG01) covered all plant AUX/IAA proteins, while in the other eight LGs with some plant AUX/IAA proteins were lost. The integration of Arabidopsis AUX/IAAs (AtIAAs) into a homologous search of AUX/IAAs in lower plants is intended to better understand the transition evolution of IAA proteins in plants. Notably, 15 (60%) Arabidopsis AUX/IAA proteins are in LG01, suggesting that AtIAAs share ancestral genes with early lower plants. The molecular evolution tree based on the sequence alignment of AUX/IAA proteins shows that the LG01 has a more complex topology across multiple lower plants. The structural analysis of protein motifs in orthogroups reveals molecular evolution patterns of low plant IAA family differences. IAA proteins in most lower plants have highly conserved motifs 1, 2, 3, and 5, while we did not find motifs 3 and 5 in some members (85038 and 36236) of the Selaginella moellendorffii. An interesting phenomenon is that complete motifs 1-10 only exist in all Physcomitrella patens family members. However, these unique motifs are once again lost in Arabidopsis IAA proteins ( Figure 1C). Moreover, the GO enrichment results for all identified proteins in the LG01 shared by lower plants showed that the proteins in the LG01 were involved in the biological process of the response to auxin, which further verified the reliability of the family identification results ( Figure 1D,E). The above conclusions reveal the early origin of conserved structures of IAA proteins. The evolution of new motif structures in Physcomitrella patens suggests the structural diversity and functional explosion of IAA proteins in bryophytes. However, as plants evolved, these nonshared motifs were lost in higher plants represented by Arabidopsis ( Figure 1B,C).
Furthermore, we used the same method to systematically classify the IAA proteins of 12 higher plants ( Figure 1F). There are more IAA family members in domesticated crops such as G. max, O. sativa, Triticum aestivum, Zea mays, and C. quinoa; notably, up to 16 G. max IAA proteins were classified into higher plant orthogroups 01 (LG01), and fewer AUX/IAAs in early differentiated species such as coffee and grapes. ( Figure 1F). Only HG04 conforms to the characteristics that cover all higher plants' AUX/IAAs, in sharp contrast to most other HGs. IAA proteins in most species contain conserved motifs 1-10 except for all Fagopyrum tataricum IAA proteins in the HG01 that lost motif 9, one maize IAA protein which lost motif 6, and one quinoa IAA protein which lost motif 5 ( Figure 1G). The IAA protein structure in other HGs is similar to that in HG01, except for the minority protein with nonconservative motifs' loss, the IAA proteins in most species have similar structural composition, suggesting that the evolutionary pattern of IAA proteins in higher plants is conserved ( Figure 1G-K).

Structural and Evolutionary Analysis of the CqIAA Family of Halophytes Crops Quinoa Reveals the Family's Generally Conserved Evolutionary Model and Rapidly Evolving Genes Selected by a Strong Environment
Quinoa has experienced long domestication and complex hybridization events in history and is a good material to study the adaptation mechanism of the agroecological environment after crop domestication [35]. To investigate the molecular evolution of the IAA protein family within the quinoa genome, we constructed phylogenetic trees in 37 CqIAAs and found that CqIAAs evolved into complex topology ( Figure 2A). The promoter region

Structural and Evolutionary Analysis of the CqIAA Family of Halophytes Crops Quinoa Reveals the Family's Generally Conserved Evolutionary Model and Rapidly Evolving Genes Selected by a Strong Environment
Quinoa has experienced long domestication and complex hybridization events in history and is a good material to study the adaptation mechanism of the agroecological environment after crop domestication [35]. To investigate the molecular evolution of the IAA protein family within the quinoa genome, we constructed phylogenetic trees in 37 CqIAAs and found that CqIAAs evolved into complex topology ( Figure 2A). The promoter region sequences of 2000 bp upstream of each IAA were further extracted for predicting potential hormone and environmental response elements ( Figure 2B, Table S3). Among the IAA families, light response elements are the most, followed by stress response elements, hormone response elements, and other environmental response elements ( Figure 2B). Further, we performed gene structure and protein motif analysis for all CqIAAs ( Figure 2C,D). The number of introns in the CqIAAs ranges from 1 to 5, and the CDS region of each gene contains at least one family-specific conserved structure annotated as AUX/IAA. CqIAA proteins with close phylogenetic relationships have similar motif composition, although some protein motifs are developed or lost in some CqIAAs. For example, motif 9 (purple box) exists only in AUR62008800 and AUR62025439, and motif 7 (red box) is detected only in AUR62002810 and AUR62014118 ( Figure 2D). Moreover, motif 2 (yellow box) is present in almost all CqIAA proteins and is their conserved motif. Interestingly, motif 2 in AUR62043556, AUR62016513, AUR62005398, and AUR62004953 are missing ( Figure 2D). The gain and loss of the motif may indicate the evolution of the function, which needs further discussion.

Evolutionary Analysis of IAA Proteins in Typical Polyploid Agricultural Crops Emphasizes the Contribution of Polyploidy to AUX/IAA Family Expansion
The expansion of the gene family, often driven by gene duplication of different ways, has repeatedly played irreplaceable roles in providing the raw material for the adaptive evolution of plants in changing complex environments [36]. Syntenic analysis between chromosomes helps us to understand the duplication of gene clusters within segments caused by ancient polyploidy events. We selected a heterotetraploid crop quinoa, a hexa- Chromosome mapping through Gff file revealed that 37 CqIAAs were unevenly distributed on 12 of 18 chromosomes, and 3 pairs of tandem duplication gene pairs were detected ( Figure 2E) which may suggest that replication events play a key role in the family's expansion. The synonymous (Ks) and non-synonymous (Ka) substitutions of all CqIAA sequences were further calculated. The results showed that the vast majority of CqIAAs were subjected to environmental purification selection, but the Ka/Ks of two pairs of genes (AUR62011942-AUR62016420 and AUR62011942-AUR62023669) was greater than 1 ( Figure 2F, Table S4). Interestingly, most gene pairs with Ka/Ks value close to 1 are AUR62011942 involved in (Table S4), which suggests that the gene may have important functions related to growth and development regulation and evolve rapidly under environmental selection. We also calculated the synonymous substitution (Ks) rate between CqIAAs, in preparation for a subsequent estimate of the expansion time of this family ( Figure 2F, Table S4).

Evolutionary Analysis of IAA Proteins in Typical Polyploid Agricultural Crops Emphasizes the Contribution of Polyploidy to AUX/IAA Family Expansion
The expansion of the gene family, often driven by gene duplication of different ways, has repeatedly played irreplaceable roles in providing the raw material for the adaptive evolution of plants in changing complex environments [36]. Syntenic analysis between chromosomes helps us to understand the duplication of gene clusters within segments caused by ancient polyploidy events. We selected a heterotetraploid crop quinoa, a hexaploid grain T. aestivum, and hybrid polyploid Brassica napus as the materials for genomic collinearity analysis (Figure 2A-C, Table S5). The results showed that 15, 94, and 141 syntenic AUX/IAA gene pairs were found in the three plant genomes, respectively (Figure 2A-C, Table S5). Interestingly, the syntenic pairs in the three subgenomes of hexaploid T. aestivum had a 2:1 ratio, which is consistent with the previously reported characteristics of the wheat MADS family [37], confirming that the AUX/IAA family is preserved with the fusion of different subgenomes of wheat ( Figure 2B). Although quinoa and B. napus are also allopolyploid with different subgenomes, the synteny and evolutionary retention patterns of AUX/IAA family are significantly different than that of T. aestivum.
In addition to synteny gene pairs within species, synteny homologs between species are widespread and often functionally related [38]. We analyzed the syntenic AUX/IAA pairs of quinoa with other higher plants, especially with gramineous plants ( Figure 2D, Table S6). Visualization showed the total syntenic relationships between quinoa and O. sativa (two pairs), T. aestivum (only one pair), S. bicolor (two pairs), Z. mays (two pairs), Arabidopsis (five pairs), and F. tataricum (seven pairs). Moreover, this suggests that species with closer phylogeny have more collinearity IAA pairs, although the syntenic relationship of AUX/IAAs are not strong among multiple species, which suggests that segment duplication is not the main cause of IAA family's expansion in the plant kingdom ( Figure 2D, Table S6).
Early studies have confirmed that another important driver of plant genome evolution is whole-genome duplication (WGD) [39]. We speculated that the expansion of the IAA family after plant landification may be due to WGD events. The Ks values of homologous pairs in syntenic regions of the genome, as well as the mean Ks values of individual syntenic blocks, indicated that WGD events had occurred in the evolutionary history of representative species ( Figure 2E). The Ks density plot of representative species showed that the quinoa genome experienced a late WGD event, but slightly earlier than the WGD event shared by angiosperms, this may be closely related to a hybridization during quinoa domestication, which may be closely related to a hybridization during quinoa domestication [15]. Interestingly, the Ks peak of quinoa genome is close to the Ks peak of CqIAAs ( Figure 2F), which may suggest that the IAA family expands with the doubling of the quinoa genome. To further investigate whether polyploidy events promote the expansion of the IAA family in representative species, we used MCScanX to identify repeat types of AUX/IAAs in eight plants, including quinoa ( Figure 2F, Table S7). As a result, most of the AUX/IAAs in C. quinoa, B. napus, Arabidopsis, O. sativa, and Sorghum bicolor were associated with polyploidy ( Figure 2F, Table S7). The combined results show that the expansion of AUX/IAA family in higher plants represented by quinoa is mainly due to polyploidization events, although AUX/IAAs in some species (T. aestivum and Z. mays) have other expansion patterns (dispersed and tandem duplication).

RNA-Seq Data Analysis of Quinoa Root Provides a Reference for the Network of Hormone Signals and Environmental Regulation Formed by AUX/IAAs
The family of polyploidy reserved transcription factors often contributes to the adaptive evolution of plants to complex environments [40], and the perception and response of these environmental factors are more coupled through complex root signaling networks [40,41]. We collected RNA-seq data sets from different tissues of quinoa from the NABI GEO database for expression trend analysis, and the normalized expression data were used to perform the clustering ( Figure 3A, Table S8). The results showed that 14 statistically significant model profiles (colored profiles) were identified in 49 profiles generated ( Figure 3A). Of these 14 profiles, three special model profiles (profiles 41, 42, and 44) were identified with significantly higher expression levels in roots than in other tissues ( Figure 3B-D). KEGG pathway enrichment analysis of genes in profiles 41, 42, and 44 via KOBAS online database ( Figure 4F-H, Table S8) was conducted. Phenylpropanoid biosynthesis (ko00940), ubiquinone and other terpenoid-quinone biosynthesis (ko00130), plant hormone signal transduction (ko04075), and other entries were significantly enriched in profile 41 ( Figure 4F, Table S8); the enrichment results of amino sugar and nucleotide sugar metabolism (ko00520), plant hormone signal transformation (ko04075), ubiquitinmediated proteolysis (ko04120), and biosynthesis of secondary metabolites (ko01110) in profile 42 have attracted our attention ( Figure 4G); ubiquitin mediated proteolysis (ko04120), plant hormone signal transformation (ko04075 , Table S8), amino sugar, and negative sugar metabolism (ko00520) in the profile 42 also have a significant rich factor ( Figure 4H). IAA is closely related to ARFs' protein-protein interaction and ubiquitination-mediated protein degradation pathway in the auxin signaling pathway, and auxin signal is also an important component of plant hormone signal transduction [5,42,43], which further leads to our strong attention to root auxin signaling process.
Moreover, a coexpression network analysis of root-specific expression genes and 37 IAA genes showed that 437 edges with significant correlation (p < 0.01) were generated when a total of 1019 nodes were input ( Figure 4I, Table S9). When using the expression data of the other four organizations for correlation analysis, there are only 347 significant correlation edges (p < 0.01) generated after the same number of nodes input ( Figure 4I, Table S9). These findings confirm that the IAA family has formed a closer regulatory network in the root system, which is in contrast to the absence of such a correlation in expression data for four other organizations. Gene ontology (GO) analysis helps us to explore the potential functions of genes that are significantly related to IAAs, and we show the top 15 terms in three parts including biological processes (BP), cellular components (CC), and molecular functions (MF). The significantly enriched terms including protein modification process (GO:0036211) (BP terms) cell wall (GO:0005618), plant-type cell wall (GO:0009505) (CC terms) protein binding (GO:0005515), and hydrolase activity (GO:0016787) (MF terms), which have caught our attention because they again remind the potential role of auxin signaling in root growth and development ( Figure 4J, Table S10). Multiple collinearity scan toolkit (MCScanX) was used to detect duplication types to identify GRASs associated with ancestral polyploidization events (WGD or segmental). The AUX/IAAs related to polyploidy were inferred according to the collinearity persistence of the region around the selected genes. AUX/IAAs were divided into five groups, including dispersed, proximal, tandem, singleton, and WGD or segmental. Multiple collinearity scan toolkit (MCScanX) was used to detect duplication types to identify GRASs associated with ancestral polyploidization events (WGD or segmental). The AUX/IAAs related to polyploidy were inferred according to the collinearity persistence of the region around the selected genes. AUX/IAAs were divided into five groups, including dispersed, proximal, tandem, singleton, and WGD or segmental.   Table S8; (E). Heat maps showing in five tissues further confirmed that the transcription level of these genes in roots was significantly higher than that in other tissues; (F-H). KEGG pathway enrichment for three highly expressed gene profiles (profile 41, 42, and 44) in roots; Detailed data are recorded in the Table S8; (I). Correlation network between genes and CqIAAs formation in 3 profiles. The yellow diamonds represents CqIAAs, and blue circle represents the genes in 3 profiles. The red line represents positive correlation; the blue line represents negative correlation. correlation network is obtained by p < 0.01 as threshold screening. Networks on the left and right are constructed using the RNA-seq data root and other organizations, respectively; detailed data are recorded in the Table S9; (J). GO enrichment of CqIAAs' target genes in correlation networks. red, green, and blue column distributions represent three terms. of biological processes, cellular components, and molecular functions. The top 15 entries in each term were chosen for presentation. Detailed data are recorded in the Table S10.

Among Quinoa with Different Fruit Phenotypes, AUX/IAA family Formed a Robust Network of Correlations with Hormone Signaling and Energy Metabolism Pathways
Transcriptome changes of seed continue to affect the development of radicles and the interaction with the environment after germination [44]. For mining genes with potential regulatory functions for root-environment interaction, we first analyzed the RNA-seq data of quinoa seeds to understand the reserve resources available at the seed stage. We selected two quinoa seeds with significant differences (p < 0.001) in size and color, named large fruit (LF) and small-fruit (SF) quinoa according to their size ( Figure 5A, Table S11), and we hierarchically clustered 37 CqIAAs expressed in two quinoa fruits based on the expression level ( Figure 5B, Table S11). Eight IAA genes (AUR62005881, AUR62031867, AUR62004956, AUR62032874, AUR62025439, AUR62011942, AUR62008800, and AUR62001002) have high transcription levels in both large and SF quinoa, while six genes (AUR62035414, AUR62011167, AUR62016513, AUR62021457, AUR62023669, and AUR62002523) were not expressed in either fruit. Notably, AUR62011942 had significantly different expression patterns (Foldchange > 2, p < 0.05) in the two quinoa fruits ( Figure 5B,  Table S8; (E). Heat maps showing in five tissues further confirmed that the transcription level of these genes in roots was significantly higher than that in other tissues; (F-H). KEGG pathway enrichment for three highly expressed gene profiles (profile 41, 42, and 44) in roots; Detailed data are recorded in the Table S8; (I). Correlation network between genes and CqIAAs formation in 3 profiles. The yellow diamonds represents CqIAAs, and blue circle represents the genes in 3 profiles. The red line represents positive correlation; the blue line represents negative correlation. correlation network is obtained by p < 0.01 as threshold screening. Networks on the left and right are constructed using the RNA-seq data root and other organizations, respectively; detailed data are recorded in the Table S9; (J). GO enrichment of CqIAAs' target genes in correlation networks. red, green, and blue column distributions represent three terms. of biological processes, cellular components, and molecular functions. The top 15 entries in each term were chosen for presentation. Detailed data are recorded in the Table S10.

Among Quinoa with Different Fruit Phenotypes, AUX/IAA Family Formed a Robust Network of Correlations with Hormone Signaling and Energy Metabolism Pathways
Transcriptome changes of seed continue to affect the development of radicles and the interaction with the environment after germination [44]. For mining genes with potential regulatory functions for root-environment interaction, we first analyzed the RNA-seq data of quinoa seeds to understand the reserve resources available at the seed stage. We selected two quinoa seeds with significant differences (p < 0.001) in size and color, named large fruit (LF) and small-fruit (SF) quinoa according to their size ( Figure 5A, Table S11), and we hierarchically clustered 37 CqIAAs expressed in two quinoa fruits based on the expression level ( Figure 5B, Table S11). Eight IAA genes (AUR62005881, AUR62031867, AUR62004956, AUR62032874, AUR62025439, AUR62011942, AUR62008800, and AUR62001002) have high transcription levels in both large and SF quinoa, while six genes (AUR62035414, AUR62011167, AUR62016513, AUR62021457, AUR62023669, and AUR62002523) were not expressed in either fruit. Notably, AUR62011942 had significantly different expression patterns (Foldchange > 2, p < 0.05) in the two quinoa fruits ( Figure 5B, Table S11), which may be a potential regulator of the differential phenotype between the two quinoa species.  Table S11); (C). Expression profiles of TCA cycle pathway genes in two quinoa seeds (expression data are recorded in Table S12); (D). Expression profiles of auxin biosynthesis and signal transduction pathway genes in two quinoa seeds (expression data are recorded in Table S12); (E). Correlation network between CqIAAs and different pathway Figure 5. Morphological display, differential expression of CqIAAs, differential expression of auxin signaling and energy metabolism, differential expression of auxin signaling and energy metabolism pathway and sexual network with CqIAAs.  Table S11); (C). Expression profiles of TCA cycle pathway genes in two quinoa seeds (expression data are recorded in Table S12); (D). Expression profiles of auxin biosynthesis and signal transduction pathway genes in two quinoa seeds (expression data are recorded in Table S12); (E). Correlation network between CqIAAs and different pathway genes in two quinoa seeds (correlation data are recorded in Table S13). Yellow diamonds represent CqIAAs; different color circles represent auxin signaling and energy metabolism pathways. The red line represents positive correlation; the blue line represents negative correlation. Correlation network is obtained by p < 0.05 as threshold screening. Networks on the left and right are constructed using RNA-seq data the small-fruit and LF quinoa seeds, respectively.
Energy and hormone metabolisms are the most vigorous biological processes in plant seed germination and seedling growth and development [45,46]. As a result, we further retrieved the gene expression data of the tricarboxylic acid cycle, auxin biosynthesis, and auxin signaling pathway from the transcriptome data sets of two quinoa fruits ( Figure 5C,D, Table S12) to mine target genes that may be closely related to the IAA family. In the TCA cycle, pyruvate dehydrogenase (PDC) gene AUR62037785 upregulated in LF quinoa while AUR62042164 downregulated. Citrate synthase (CSY) is the rate-limiting enzyme in this pathway, but there is no significant difference in its expression in the two quinoa fruits. Furthermore, the expression level of isocitrate dehydrogenase (IDH) gene AUR62027128, succinate dehydrogenase (SDH) gene AUR62043864, and malic dehydrogenase (MDH) gene AUR62032001 in LF quinoa were all higher than SF quinoa. Notably, the α ketoglutarate dehydrogenase (OGDC) gene AUR62028926 coupled carbon metabolism with nitrogen metabolism [47], and it had a higher expression level in SF quinoa than in LF quinoa, suggesting that SF quinoa may have potential nitrogen reserves in the early stages of germination ( Figure 5C, Table S12). Unexpectedly, the expression levels of two key enzyme genes (TAA and YUC) of the auxin biosynthesis pathway in both large and SF quinoa were low ( Figure 5D, Table S12). A recent report confirmed that LRR-RLK was able to respond to soil nitrogen signals and, coupled with auxin signals [48], had higher expression levels in both fruits. Moreover, the GH3, SAUR from the auxin signaling pathway were almost nonexpressed in both fruits, and the proteasome genes expression in the ubiquitination pathway is active, which is in contrast to that of the coding genes of the SCF complex which appear to be silent ( Figure 5D, Table S12). Moreover, using transcriptome data from two quinoa seeds, we constructed a correlation network of genes CqIAAs, auxin biosynthesis, auxin signal transduction, and TCA cycle pathway genes ( Figure 5E, Table S13). Interestingly, CqIAAs form a robust regulatory network with metabolic pathway genes. Hormone signaling and energy metabolism regulation network are more robust in SF quinoa, which suggests that it may have greater adaptability to a complex agroecological environment during germination ( Figure 5E, Table S13).

Metabolic Regulatory Networks Formed by IAA May Help Induce Stress Responses of Root Systems of Quinoa Seedlings in Nitrogen Starvation Agroecological Environment
TCA cycle is the core pathway of cell energy metabolism, which also provides essential precursors for amino acid biosynthesis and maintains metabolic balance with carbon [49]. In order to verify whether the two quinoa seeds with different energy and hormone signal networks formed different environmental responsiveness after germination, we induced the germination of the two quinoa seeds under the same conditions. Quinoa seedlings were treated with water culture under nitrogen-deficient environment when the cotyledons were fully expanded. The results showed that after two weeks of a nitrogen-deficient culture, the roots of both quinoa seedlings were elongated ( Figure 6A, Table S14). Specifically, the root elongation of SF quinoa before and after stress is significantly different, which is in sharp contrast to LF quinoa; there was no significant difference in the increase in root biomass between the two species ( Figure 6B, Table S14). The determination of cell wall components in the root system of two kinds of quinoa showed that the decrease in pectin in SF quinoa root was more obvious than that in the LF quinoa root system; the increase in xylogen in the root of LF quinoa resveratrol was more obvious than that of the SF quinoa ( Figure 6C, Table S14). We also determined the chlorophyll content in the leaves of quinoa. The results showed that the chlorophyll A, chlorophyll B, and total chlorophyll content in the leaves of the two kinds of quinoa decreased after stress, and the decrease in LF quinoa was more significant ( Figure 6D, Table S14). Through the principal component analysis (PCA) of the compound content of two kinds of quinoa under two treatments, it was found that the above indexes of the two kinds of quinoa were close before the stress treatment, while the stress treatment had an obvious influence on the above indexes ( Figure 6E). H 2 O 2 , SOD, GSH, POD, Vc, and CAT levels are often used to reflect the extent of damage to plant cells, and the determination of these physiological indexes can help us to compare the response ability of the two quinoa to nitrogen deficiency stress ( Figure 6F, Table S14). The results showed that the content of H 2 O 2 was upregulated in both root systems and more significantly in LF quinoa after stress treatment; the content change of POD and CAT is similar to that of H 2 O 2 . Notably, both SOD and GSH levels in the stress-treated LF quinoa roots were significantly upregulated and downregulated, but this change was not significant in the SF quinoa roots. Moreover, the changing trend of Vc content in root system of two quinoa was the opposite in the root system of SF quinoa, Vc content was significantly upregulated after stress, which may help to alleviate oxidative damage caused by stress. PCA analysis of stress injury indicators showed that there were some differences in metabolites content in the roots of the two quinoa before stress, but the difference was more obvious after stress ( Figure 6G).
The auxin and TCA cycle pathway genes in quinoa seeds may promote the root elongation of seedlings in a nitrogen starvation environment [49]. In order to explore the causes of differential phenotypes of two quinoa (LF and SF quinoa) under nitrogen starvation stress, we first mined the differentially expressed genes of auxin and TCA cycle pathways by RNA-seq data ( Figure 6H, Table S15). We identified 9 (upregulation) and 10 (downregulation) differentially expressed genes in the IAA protein family, auxin biosynthesis and signal transduction pathways, and TCA cycle pathways ( Figure 6H, Table S15). On these bases, we further detected the expression levels of these pathway genes by the qPCR method ( Figure 6I, Table S15). In total, 12 genes were found to have different expressions before and after nitrogen starvation treatments of the two quinoa ( Figure 6I, Table S15). The riboflavin mono-oxygenase (YUC) genes in the auxin synthesis pathway were more significantly upregulated in LF quinoa, while the TRYPTOPHAN AMINOTRANSFERASEs (TAA) genes were significantly downregulated in SF quinoa; the expression of these three genes was not significantly different before and after treatment with LF quinoa. In the auxin signaling pathway, the TIR1 gene was upregulated in LF quinoa after stress treatment, and TIR2 was upregulated in both quinoas. Interestingly, the encoding gene of ubiquitin-proteasome (UbP) was downregulated in LF quinoa, which is in strong contrast to its upregulation in quinoa. The expression level of PDC, IDH, and OGDC in both quinoas was upregulated in the tricarboxylic acid cycle pathway; MDH was significantly upregulated in LF quinoa, in contrast to its significant downregulation in SF quinoa ( Figure 6I, Table S15). The expression changes of OGDC (AUR62028926) and IAA (AUR62011942) in SF quinoa root after nitrogen starvation treatment are most significant (p < 0.001), respectively acting as key hubs in auxin signal and carbon-nitrogen metabolism balance, which may potentially promote root elongation of quinoa in a nitrogen-deficient environment ( Figure 6I, Table S15). The above conclusions help us to understand the response mechanism of two quinoa seedlings in response to environmental nutritional deficiency and provide candidate genes for breeding varieties with tolerance to nitrogenstarvation soil. lation in SF quinoa ( Figure 6I, Table S15). The expression changes of OGDC (AUR62028926) and IAA (AUR62011942) in SF quinoa root after nitrogen starvation treatment are most significant (p < 0.001), respectively acting as key hubs in auxin signal and carbon-nitrogen metabolism balance, which may potentially promote root elongation of quinoa in a nitrogen-deficient environment ( Figure 6I, TableS15). The above conclusions help us to understand the response mechanism of two quinoa seedlings in response to environmental nutritional deficiency and provide candidate genes for breeding varieties with tolerance to nitrogen-starvation soil. Figure 6. Physiological changes and gene expression differences of two quinoa seedlings in nitrogen-deficient environment. (A). Large-fruit and small-fruit quinoa seedlings under nitrogen deficiency stress. The first row was the big quinoa seedlings; the second row was the small quinoa seedlings; the left side and the right side were the control group and the experimental group, respectively; (B-D). Changes of physiological indexes of two quinoa seedlings under different treatments; ** statistical significance for p < 0.05; *** statistical significance for p < 0.001; (E). Principal component (PCA) analysis of physiological index data of two materials under two treatments. The points of red, green, blue, and purple represent the control group of LF quinoa, the control group of SF quinoa, the experimental group of LF quinoa, and the experimental group of SF quinoa. Five biological repeats in each group (F). Response to stress responses of two quinoa varieties under different treatments (G). Principal component (PCA) analysis of physiological index data of two materials under two treatments. (H). Auxin biosynthesis, auxin signal transduction, and differential expression genes of TCA cycle pathway in two quinoa seeds. The volcanic map is generated at a Foldchange > 2, p value < 0.05 threshold setting. Red and blue dots indicate genes that are significantly upregulated and downregulated in LF quinoa seeds, respectively; (I). Relative expression of genes CqIAAs, auxin pathway, and TCA pathway in different treatments of two quinoa. The genes used for qRT-PCR detection met the following conditions: the expression of quinoa seeds was significantly different and derived from the auxin pathway or TCA cycle pathway.

AUX/IAA Family May Originated after Plant Terrestrialization, Influenced by WGD Events, Formed Robust Networks with Auxin Signaling and Energy Metabolic Genes in Quinoa
The biosynthesis origin and signal transduction of auxin, as core components in the regulation of growth and development of plants, have long been a focus of both plant genome evolution and physiological research [50,51]. A large number of reports have been made on the transcriptional regulation mechanism of auxin signaling pathway in recent years, but they mainly focus on auxin response factor (ARF). Clear evolutionary biological evidence showed that the ARF family originated from early aquatic algae; in particular, the early origin and diversification of ARF in charophytes was recently delineated [52]. However, the auxin-reactive protein (AUX/IAA) family, which acts as an ARF interaction factor, also plays an irreplaceable role in the regulation of the auxin pathway but has little to report in the origin, evolution, and function of the family in depth. In this study, we compared almost all currently published aquatic algal genomes with the conservative structure of the AUX/IAA family, but without finding targets with significant structural consistency. This coincides with the findings of a recent report that genes encoding auxin signal transduction elements AUX/IAA, SAUR, and GH3 are missing in the genomes of Mesostigma viride and Chlorokybus atmophyticus [53]. Typical AUX/IAA family member was found in representative mosses, although they were small numerous and less structurally conserved. Identification of orthologous genes revealed that AUX/IAA families in most higher plants represented by Arabidopsis shared their original ancestors with bryophytes ( Figure 1A-C). The evolution pattern of AUX/IAA is more conservative in higher plants, and the sequences of orthologous proteins in most plants are highly similar. In domesticated crops represented by quinoa, except for a few genes (such as AUR62011942) with a relatively fast evolution rate, most AUX/IAAs have conserved evolution patterns by strong purification selection ( Figure 1G-K and Figure 2). After plant landification, AUX/IAA family expands and may be strongly influenced by several whole-genome duplication (WGD) events shared by angiosperms ( Figure 3E,F). A lack of a strong synteny relationship between quinoa CqIAAs and AUX/IAAs in other species may suggest that the family's evolutionary pattern in quinoa is relatively independent ( Figure 3D).
Most angiosperms have experienced major extinction events in history, such as the third great extinction in the Cretaceous. A large proportion of these species are thought to have escaped extinction destiny through the genome duplication strategy [54]. Due to the more reservations of stress-responsive genes in their different subgenomes, polyploid plants have been speculated to present a stronger ability to survive in the extreme ecological environment [55]. Polyploidized genes often form complex transcriptional regulatory networks to maintain plant homeostasis in response to extreme environments [56]. This study found that the IAA family of polyploid crops represented by quinoa expanded mainly through polyploidization events, and the CqIAAs retained in the quinoa tetraploid genome formed a robust regulatory network with various complex metabolism pathway genes in root. However, there is no strong correlation between CqIAAs and the regulatory network of metabolic pathways in the other four tissues (stem, leaf, flower, and fruit) ( Figure 4I). Comparative transcriptome analysis showed that the robustness of the regulatory network involved by CqIAAs is different in the seed of two different quinoa ( Figure 5), which may cause heterogeneity in energy metabolism and nutrient flow in the germination stage. The above conclusions suggest that the transcription factors family related to growth and development from the quinoa genome does form a closer metabolic regulatory network after polyploidization, which promotes its long-term domestication and adaptive evolution in extreme agroecological environments. Phylogenetic analysis of cross-species from aquatic lower plants to terrestrial angiosperms and multipathway regulatory network mining helps us to have a deeper understanding of the evolutionary origin and functional innovation of the AUX/IAA family.

Root Elongation and Stress Response of Quinoa Seedlings in Nitrogen Starvation Are Helpful for Alleviating Nutritional Deficiency, and Small-Fruit Quinoa Has Better Potential Adaptability in Poor Agricultural Environmental
Nitrogen has diverse and irreplaceable necessities in various stages of plant growth and development. By being the building blocks for participants in basic biological processes such as nucleic acids, proteins, and chlorophyll, they provide important components for cell metabolic activity and structure [57]. In addition, nitrogen as a raw material for tryptophan metabolism plays a role in plant hormone-mediated growth and development regulation [58]. Nitrogen homeostasis in plant cells strongly depends on mitochondrialmediated carbon-nitrogen transformation equilibrium [49]. More and more reports suggest that mitochondrial function and auxin are closely interconnected: mitochondrial metabolic activity regulates auxin signaling, which promotes organogenesis and morphogenesis in plants [59]; in Arabidopsis thaliana, the retrograde signal of mitochondria can be connected to the polar distribution of auxin, which changes the phenotype of plant roots [60,61]; in rice, mutations in mitochondrial structural genes affect auxin biosynthesis and signaling pathways, thus altering the morphological structure of plant vascular bundles [62]. The effect of mitochondria on the auxin pathway is in fact coupled through TCA cycles. Amino acids are important intermediates substrates that transition from carbon metabolism to nitrogen metabolism in the TCA cycle and can enter the tryptophan-dependent IAA biosynthesis pathway [63,64]. In addition, studies have shown that plants can interact with microorganisms to mediate salt stress response via ACC-deaminase and auxin activity [65]. The synthesis of phenol and flavonoid secondary metabolites and the change of endogenous hormones such as abscisic acid and auxin also promoted plant growth and development and stress resistance [66]. This further suggests that metabolism balance plays an important role in crop resistance to environmental stress.
The transcriptional levels of such genes in the TCA cycle and auxin pathway were significantly different in the seeds of the two quinoa materials selected in this study, suggesting that the carbon-nitrogen metabolism of two quinoa species during germination may enter into different metabolic streams. A generally accepted view is that the seed germination stage is most easily disturbed by environmental stress, and the stress responseability of seedlings directly determines the survival rate of plants [67,68]. Our results of the experiment are as expected that under nitrogen deficiency stress, the root of SF quinoa seedlings had obvious elongation, and chlorophyll, cell wall composition, and physiological indexes related to oxidative damage showed better characterization than LF quinoa seedlings ( Figure 6A-G). qPCR detection confirmed that the expression of TCA and auxin-related genes were strongly induced in roots under nitrogen-deficient stress ( Figure 6I). The above conclusions further confirm our inference that the AUX/IAA family in SF quinoa forms a robust regulatory network with TCA cycle and auxin pathway genes, thus enabling the carbon-nitrogen metabolism balance to be more actively deployed in the nitrogen deficiency environment, thus maintaining the physiological homeostasis in the barren environment. Noticeably, the rapidly evolving gene AUX62011942 found by evolution analysis not only has significantly different expression levels in LF quinoa seeds and SF quinoa but also shows a strong response when young roots of SF quinoa face nitrogen stress. Moreover, CqIAA was significantly upregulated under nitrogen stress in SF root, while CqARF was downregulated, in sharp contrast to the consistent downregulation of the two genes during the stress of LF root ( Figure 6I). Phytohormone plays an important role in plants' resistance to environmental stress [69]. Herein, we present a possible mechanism model: before germination, SF quinoa seeds form a hormone and energy metabolism regulation network centered on IAA, and after germination in a nutrientdeficient environment, the complex regulation network affects auxin signal, thus inducing root elongation to seek more nitrogen sources. By coupling the energy metabolism pathway, the intracellular carbon-nitrogen metabolism was rearranged, which together maintained the homeostasis of seedlings in the nitrogen deficiency environment. The above inference needs to be supported by further experimental evidence. To sum up, it is a more efficient genetic and agricultural breeding strategy to excavate the dominant genes preserved by environmental selection at the evolutionary scale, thus providing valuable insights for understanding the evolutionary strategies of crops adapting to the complex agroecological environment in the course of domestication.

Conclusions
Our study examines the AUX/IAA family in 19 plants with agricultural and evolutionary importance, from aquatic algae to angiosperms. The results show that the AUX/IAA family originated after the territization of plants and expanded in angiosperms, which has a conservative evolutionary pattern and finds general conservation of function across species. We also emphasized AUX/IAAs' expansion in angiosperms is due to polyploidization events. We tapped the complex regulatory networks from quinoa transcriptome datasets, which found that AUX/IAA family members show interactions with auxin signaling and energy metabolism genes. More importantly, the stress response of two different quinoa materials under nitrogen starvation was different. CqIAAs, hormone signaling, and energy metabolism pathway genes in quinoa seedlings showed significant responses under nitrogen starvation stress, which provide guidance for understanding the adaptive mechanism of crop adaptation to poor agroecological environments.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/agriculture11080780/s1, Table S1. Basic information on 495 IAA proteins identified from 17 plant genomes, Table S2. Orthogroup identification of multispecies IAA family proteins, Table S3. Prediction of cis-acting elements in promoter regions of 37 quinoa CqIAAs, Table S4. The synonymous (Ks) and nonsynonymous (Ka) substitutions rates of all CqIAAs, Table S5. Synteny IAA gene pairs in the genome of Chenopodium quinoa, Triticum aestivum and Brassica napus, Table S6. Synteny IAA gene pairs between Chenopodium quinoa and other species, Table S7. Identification of repeat types of IAAs in different species, Table S8. Organizational expression matrix, trend analysis of gene expression of 49 profiles and KEGG pathway enrichment analysis of three Gene profiles specifically expressed in roots, Table S9. Coexpression networks formed by CqIAAs in roots and other tissues, Table S10. GO enrichment analysis of target genes associated with IAA in the coexpression network, Table S11. Transcription level (FPKM) of 37 CQIAA in two quinoa species and morphological data of seeds in two quinoa species, Table S12. Transcription level (FPKM) of auxin biosynthesis, auxin signal transduction, and tricarboxylic acid cycle pathway genes, Table S13. Coexpression networks in LF and SF quinoa seeds formed by CqIAA and hormone and energy metabolic pathway genes, Table S14. Physiological indexes of two kinds of quinoa under nitrogen stress, Table S15. Gene expression of roots of two quinoa seedlings under nitrogen deficiency stress, Table S16. Primer sequence for qRT-PCR experiment.
Author Contributions: H.Y. designed and participated in the whole study and wrote the manuscript; S.W. analyzed and visualized plant genome data; Y.Y. carried out plant culture, physiological experiments, and data analysis; H.X. was involved in the comparative genome analysis of multiple species; G.W. was involved in metabolic pathways and regulatory network analysis; J.W., Y.H. and L.J. summarized and analyzed data generated from the experiments; H.C. provided primary guidance and platform support throughout the study and reviewed and approved the manuscript. All authors have read and agreed to the published version of the manuscript.