Skip to main content

Genome-wide identification of the AAT gene family in quinoa and analysis of its expression pattern under abiotic stresses

Abstract

Background

Plant amino acid transporters play an important role in the absorption of soil amino acids by roots, the transport of amino acids between xylem and phloem, plant growth and development, and response to abiotic stress.

Result

In this study, we identified 147 AAT genes in the quinoa genome sequence and categorized them into 12 subfamilies on the basis of their similarity and phylogenetic relationships with AAT found in Arabidopsis thaliana. Interestingly, these AAT genes are not evenly distributed on the quinoa chromosomes. Instead, most of these genes are centrally located on the outer edges of the chromosome arms. After performing motif analysis and gene structure analysis, we observed the consistent presence of similar motifs and intron-exon distribution patterns among subfamilies. Tissue expression analysis revealed that CqAAT gene was less expressed in fruits and more expressed in roots, stems, leaves and flowers. Meanwhile, expression analysis under four adversities of high temperature, low temperature, waterlogging, and drought and different treatments of nitrogen, phosphorus, and potash fertilizers found that two genes of the CqGAT subfamily, AUR62031750 and AUR62023955 were up-regulated expressed under abiotic stresses.

Conclusions

In summary, there is a significant differentiation in the tissue expression and stress expression of the CqAAT gene, indicating that CqAATs play a role in regulating growth and development under abiotic stress.

Peer Review reports

Introduction

Quinoa (Chenopodium quinoa Willd.) is an annual dicotyledonous plant of the chenopodium genus in the chenopodium subfamily of the Amaranthaceae family, with higher nutritional content of proteins, carotenoids, and vitamin C in its seeds than in most grains such as rice and wheat. Quinoa is considered by the Food and Agriculture Organization of the United Nations to be the only monocot plant that can meet the basic nutritional needs of the human body, and quinoa is officially recommended as the most suitable fully nutritious food for human beings [1, 2]. Quinoa seedlings and young leaves are edible and, like spinach of the same family, are commonly used as a vegetable, providing important vitamins, nutrients and minerals to humans [3, 4].

Amino acids are not only necessary for protein formation, but also the main form of nitrogen absorbed by the plant root system, which plays an important role in the growth and development of plants. In highbrow plants, the absorption and transportation of amino acids mainly depend on amino acid transporter proteins (AAT), which are widely involved in important physiological processes such as amino acid uptake and transportation, root growth, floral organ development, seed germination, and abiotic stress in plants. The AAT gene family has been systematically identified and characterized in a variety of plant species. There are 63 AAT genes in Arabidopsis thaliana [5], 84 in Phaseolus vulgaris [6], 85 in rice [7], 189 in soybean [8], 296 in common wheat [9], 94 in cereal [10], 104 in Tartary buckwheat [11], and 72 in potato [12], indicating that the AAT gene family is widespread in higher plants. The AAT family in plants consists of the amino acid/Auxin permease (AAAP) family and the amino acid-polyamine-choline (APC) family of transporter proteins.The AAAP gene family is one of the largest gene superfamilies of AAT, and it can be divided into eight branches on the basis of sequence similarity and conserved structural domains, namely, the amino acid permease (AAP), the lysine and histidine transporter (LHT), the gamma- amino butyric acid transporter (GAT), Auxin permease transporter (AUX), proline transporter (ProT), aromatic and neutral amino acid transporter (ANT), and amino acid-like transporter (also classified as ATLa and ATLb) [13,14,15]. Members of the AAAP family belonging to different branches have low sequence similarity, but all AAAP proteins have the same conserved structural domain, PF01490. The APC family consists of four subfamilies: the cationic amino acid transporter protein (CAT), amino acid/polyamine transporter protein (ACT), and polyamine H+ homotransporter protein (PHS), and the tyrosine-specific transporter protein (TTP) family, all ACP proteins have the same conserved structural domain, PF00324.

Amino acid transporter proteins have key roles in plant growth and development. AAP6 in Arabidopsis may be responsible for the uptake of amino acids in xylem, and AAP8 may be important in the uptake of amino acids in seeds [15]. LHT1 is responsible for the uptake of soil amino acids by plants and the transport of amino acids by members of the LHT subfamily is also involved in coordinated interactions between pollen development and pollen gynoecium in plants [1617]. Amino acid transporter proteins are not only responsible for the uptake and transportation of amino acids, but also for the transportation of non-amino acid substances. They may function as receivers of amino acid signals to regulate growth and development, modulate growth and yield traits, change plant shape and enhance disease resistance [18,19,20]. If amino acid transporter proteins are mutated or overexpressed in plants, some amino acids may act as signaling molecules to cause a series of complex signaling reactions to alter gene expression and ultimately lead to phenotypic changes [21]. In agricultural production, planted crops are often subjected to yield reductions due to abiotic stresses, such as high temperature, drought, low temperature, and waterlogging. Many AAT genes can also alleviate the damage of various abiotic stresses on plants. For example, the ProT gene can rapidly partition proline under water stress conditions, thus minimizing the damage caused by water deficit in plants [22]. GmProT1 and GmProT2 are stress-sensitive in roots, and transgenic plants can transport more proline after proline spraying [23]. Increased γ-aminobutyric acid significantly induced AtGAT1 expression in mechanical injury, stress, or sensitive environments [24]. OsAAP6, OsAAP11, and OsANT3 were significantly up-regulated under salt and drought stress, and they could regulate the plant’s tolerance to adversity [7].

Systematic analysis of the AAT gene family in quinoa has not been reported. Quinoa is an important grain, feed and industrial raw material crop in China, as well as the largest crop grown in the world, so the study of quinoa AAT family genes has important theoretical and practical significance. In this study, based on the whole genome level of quinoa, we used bioinformatics to identify the members of quinoa AAT family, and analyzed their evolutionary relationship and tissue expression relationship to reveal the expression pattern of quinoa AAT family genes in different tissues and under adverse conditions.

Results

Identification and physicochemical properties analysis of CqAAT gene family members

Based on bidirectional Blast analysis and HMM (PF00324 and PF01490), a total of 147 AAT genes were identified in the quinoa genome. The physicochemical properties of the protein sequence encoded by the quinoa gene were analyzed on the ExPASy website, as shown in Table S1. The CqAAT protein contains an average of 466 amino acids, ranging in length from 141 (AUR62023532) to 1425 (AUR62040229) amino acids. The molecular weight ranges from 15,381 kD (AUR62009229) to 156543.52 kD (AUR62040229), with an average of 51085.53 kD. The maximum and minimum values differ by 141162.52 kD. The isoelectric points (PI) range from 5.19 (AUR62000014) to 10.12 (AUR62001226), with 105 members of the CqAATs gene family having isoelectric points exceeding 7, indicating that the majority of CqAATs belong to alkaline proteins. The instability index ranged from 21.1 (AUR62027056) to 54.62 (AUR62014264), with 129 genes having an instability index less than 40, indicating that most CqAATs belong to stable proteins. The aliphatic index ranges from 81.86 (AUR62043653) to 129.3 (AUR62027047); The hydrophilicity ranges from − 0.312 (AUR62043653) to 0.923 (AUR62018347), indicating that most quinoa AAT family proteins belong to hydrophobic proteins, with only three hydrophilic proteins (AUR62014264, AUR62043653, AUR62040228). Subcellular localization prediction showed that most proteins were localized to the cell membrane, with 22 proteins localized to Golgi apparatus and 4 proteins localized to Chloroplast (Table S1).

Phylogenetic analysis of AAT gene family members in quinoa

To investigate the evolutionary relationship of the quinoa AAT gene family, 147 quinoa AAT protein sequences and 63 Arabidopsis AAT protein sequences were aligned. Through phylogenetic tree analysis, the 147 CqAAP proteins were divided into 12 subfamilies (Fig. 1): the ATLb subfamily contains 26 AAT members, the ANT subfamily contains 2 AAT members, the ATLa subfamily contains 12 AAT members, the TTP subfamily contains 2 AAT members, the ACT subfamily contains 6 AAT members, the PHS subfamily contains 4 AAT members, the CAT subfamily contains 18 AAT members, the ProT subfamily contains 6 AAT members, the AUX subfamily contains 8 AAT members, and the GAT subfamily contains 14 AAT members. There are 23 AAT members in the AAP family and 26 AAT members in the LHT family; In addition, there are significant differences in the number of branches at all levels and the number of genes in each branch among different subfamilies, indicating the existence of some short specific sequences between different subfamilies, which may endow the CqAAT family with diverse functions.

Fig. 1
figure 1

Phylogenetic relationship constructed based on amino acid sequences of AAT protein in quinoa and Arabidopsis thaliana. Different subgroups are highlighted in different colors

Gene structure and conserved domain analysis of AAT gene family members in quinoa

To further understand the structural characteristics of CqAATs proteins, 20 conserved motifs were identified from CqAATs using MEME online software (Table 1). Further analyze the gene structure of CqAATs using the GSDS2.0 gene structure display system, and finally draw the gene structure and conserved motif diagram based on the phylogenetic tree relationship (Fig. 2). Based on genome annotation, this study analyzed the intron exon structure of CqAAT family genes. Among the members of the CqAAT family, 10 gene sequences do not contain introns, while the number of introns in the sequences of other family members ranges from 1 to 20. AUR62040229 contains the most exons and introns and 67 genes having no untranslated regions. Similar intron exon structures were observed in most CqAAT genes within the same subfamily, such as CqACT, CqATLb, CqATLa, CqLHT, CqAUX, CqGAT and CqCAT. Conservative motif analysis shows that the number of motifs in the CqAATs gene family members ranges from 1 to 35. Nine genes only contain Motif 11, while AUR62040229, AUR62009231, and AUR62035636 have a higher number of motifs due to the presence of repetitive fragments. It is speculated that during the evolution of genes, certain sequence regions may remain highly conserved due to their importance to gene function. These conserved sequence regions (motifs) are preserved in gene family members and may appear multiple times in the genome through gene replication or evolution. The types and quantities of conserved motifs among members of the same subfamily are very similar, but there are also differences in the patterns of conserved motifs among members of the same subfamily (Fig. 2). Motif 12 only exists in the AAP subfamily; Motif 7 only exists in the ACT, ATLa, and ATLb subfamilies, while Motif 11 only exists in the ACT, PHS, and CAT subfamilies; Motif 10 only exists in the GAT subfamily; Motif 18 only exists in the AUX subfamily; Motif 19 only exists in the AAP and LHT subfamilies; Motif 14 only exists in the ATLa and ATLb subfamilies.

Table 1 Motif sequences of AAT in Quinoa

A B C.

Fig. 2
figure 2

Intraspecific clustering, conserved motifs, and gene structure analysis of the AAT gene family in quinoa. (A) Phylogenetic tree of the AAT gene family in quinoa; (B) 20 conserved AAT protein motifs, with different colored squares representing regions of conserved motifs and black lines representing regions without conserved sequences; (C) The exon/intron structure composition of the AAT gene, with green squares representing upstream/downstream regions, yellow squares representing exons, and black lines representing introns

Chromosome localization and co linear analysis of AAT gene family members in Quinoa

Chromosome localization analysis revealed that 147 CqAAT genes were unevenly distributed on 18 chromosomes (Fig. 3A). The number of AAT genes distributed on different chromosomes ranged from 2 to 13, with chromosome 1 having the most genes and containing 13 AAT genes, followed by chromosomes 4, 8, and 15; Chromosomes 2, 11, and 13 contain the least number of AAT genes, with 2 AAT genes and 12 genes not located on chromosomes. In order to further investigate the evolutionary relationship between quinoa AAT genes and other species, we conducted collinearity analysis on quinoa with Arabidopsis and rice. The results showed that 147 quinoa AAT genes identified 34 pairs of collinear bases with Arabidopsis, involving all chromosomes of quinoa except for chromosome 9 and 18. Only 15 pairs of collinear gene pairs were identified between rice and Poaceae, suggesting a large-scale gene loss event in the AAT family after differentiation between Poaceae and Chenopodiaceae. (Fig. 3B). In order to elucidate the evolutionary relationship between the AAT gene family in quinoa, we constructed a collinear map of the AAT genes in quinoa and found that 44 pairs of AAT genes had collinear relationships among the 18 chromosomes of quinoa (Fig. 3C).

Fig. 3
figure 3

(A) The distribution of 147 CqAAT genes identified on chromosomes. Red on chromosomes represents high density, while blue represents low density. (B) Inter species collinearity analysis of quinoa, rice, and Arabidopsis. The blue lines indicate collinearity between gene pairs. (C) The collinearity map of different AAT genes in quinoa. The red lines indicate collinearity between gene pairs. The first circle from the inner layer to the outer layer represents gene density, and the second circle represents chromosome coordinates

Analysis of the cis acting components of the promoter

The cis acting regulatory elements are transcription factors that regulate gene expression. To further explore the potential biological functions of the AAT gene family genes in quinoa, PlantCARE was used to predict the promoter elements of 147 AAT gene family genes. Classify cis acting elements into five categories, namely phytohormone responsiveness (including ABRE, CGTCA-motif, CARE-motif, etc.), tissue-specific expression (including ARE, RY element), stress responsiveness (including MBS, LTR, GC motif, etc.), light responsiveness (including G-box, GT1 motif, MRE, etc.), and cell cycle regulation (CAT box, MSA like, etc.) (Fig. 4). We found that there are a large number of cis acting elements related to endogenous hormones and light response, and all genes contain cis acting elements related to endogenous hormone response and light response. The number of cis acting elements related to methyl jasmonate and salicylic acid is the highest, and the types of cis acting elements related to auxin gibberellin are the highest. Most genes contain at least one cis acting element related to stress response, with the highest number of cis acting elements involved in low-temperature response. These results indicate that CqAATs are not only involved in regulating plant growth, but also in various physiological processes related to hormone signaling responses and abiotic stress responses.

Fig. 4
figure 4

Analysis of the cis-acting elements of the CqAATs genes. The horizontal axis represents different biological processes with different colors and is divided into different cis acting elements; The vertical axis represents the gene name; The numbers in the figure represent the number of cis acting components

Expression patterns of CqAAT gene family in different tissues

Display the tissue-specific expression pattern of AAT genes in quinoa using a heatmap (Fig. 5), including the transcription levels of AAT in roots, stems, leaves, flowers, and fruits. There are significant differences in the expression levels of AAT genes in various tissues of quinoa. Genes such as AUR62012183, AUR62022819, AUR62042156, AUR62006877 and AUR62000530 have high expression levels in roots, AUR62012183, AUR62000530 and AUR62006877 have high expression levels in stems, AUR62009457, AUR62023139 and AUR62003741 has high expression levels in leaves, and AUR62011899 is not expressed in roots, stems, or leaves. The above indicates that CqAATs genes are expressed in a tissue-specific manner in quinoa. The expression patterns of the AAT gene in the flowers and fruits of two different colored quinoa are basically the same, indicating that the AAT gene is highly conserved in the evolutionary process. Through heatmap, it was found that AUR62023955 had highest expression levels in white quinoa flowers, and AUR62009457 had highest expression levels in yellow quinoa flowers; The genes with the highest expression levels in the fruits of white sweet quinoa and yellow bitter quinoa are AUR62038452 and AUR62043362. AUR62001226, AUR62025868, and AUR62009230 are not expressed in two types of quinoa fruits. We speculate that the AAT gene is also expressed in a tissue-specific manner in the flowers and fruits of quinoa (Table S2).

Fig. 5
figure 5

Expression profiles of CqAAT gene in different tissues. (A) Schematic diagram of different tissues of quinoa. (B) Schematic diagram of the expression of 9 genes in yellow quinoa tissue. Different colors indicate the expression level of each tissue, with red indicating high expression level and blue indicating low expression level. (C) Heatmap of gene expression of CqAATs in different tissues. The horizontal axis represents different tissues; The vertical axis represents the gene name. Red indicates high expression level, while blue indicates low expression level

Expression pattern of CqAATs gene under drought stress

In order to investigate whether CqAATs gene expression is related to drought stress, based on transcriptome data from previous drought stress studies conducted by our research group, a gene expression heatmap was constructed. The results showed that compared to the control group, 13 genes were significantly down regulated and 9 genes were significantly up-regulated in the drought treatment group; Compared with the control group only the expression levels of AUR62027047 and AUR62018121 were significantly increased, while the expression levels of AUR62040086 and AUR62038450 were significantly decreased in rehydration treatment. We found that the expression levels of AUR62035636, AUR62018121, and AUR62017170 increased after drought treatment, while the expression levels decreased after rehydration, which was similar to the control group; The expression levels of genes such as AUR62027047, AUR62035199, AUR62023955, AUR62041352, and AUR62001226 decreased after drought treatment, while their expression levels increased after rehydration, similar to the control group; It is speculated that these genes may be key genes in response to drought stress in quinoa, and they belong to the AAP, ATLb, GAT, AUX, and ProT subgroups (Fig. 6A, Table S3).

Fig. 6
figure 6

Expression profile of AAT gene in quinoa under water stress. Red in the heatmap indicates high expression level, while blue indicates low expression level. (A) Expression profile of AAT gene in quinoa under drought stress. DC group represents 5 days drought control. DR group represents 5 days of drought. RC group represents 1 days of rewatering control. RW group represents 1 days of rewatering. (B) Expression profile of AAT gene in quinoa under waterlogging stress. C1 group represents 7 days waterlogging control; TR group represents 7 days of waterlogging. C2 group represents 7 days recovery control. R group represents 7 days of recovery

Expression pattern of CqAAT gene under waterlogging stress

Based on transcriptome sequencing data, this study analyzed the expression characteristics of the AAT gene family in quinoa in response to waterlogging stress. The results showed that under 7 days of continuous waterlogging stress, three genes, AUR62030901, AUR62027047, and AUR62040086, showed a significant downregulation trend, while 13 genes including AUR62009981, AUR62023955, AUR62031750, and AUR62025861 showed significant upregulation. Compared with the control group, 7 genes were significantly down regulated and 13 genes were significantly up regulated after rehydration treatment. Among them, 9 genes (including AUR62023955, etc.) continued to be up regulated during both waterlogging and rehydration periods, but the difference in expression was reduced compared to the stress stage, indicating that these CqAATs may be involved in the rapid response mechanism of plants to waterlogging stress. From the heatmap, it can be seen that there are still significant differences between the experimental group and the original control group after 7 days of rehydration, indicating that short-term rehydration treatment is not sufficient to fully restore the metabolic network of the plants to steady state. Systematic evolutionary analysis shows that AUR62031750 and AUR62023955, which are consistently highly expressed, belong to the GAT subfamily (Fig. 6B, Table S4).

Expression pattern of CqAAT gene under high temperature stress

Based on transcriptome sequencing data, this study analyzed the response patterns of members of the quinoa AAT gene family to high temperature stress. Differential expression analysis showed that a total of 35 significantly differentially expressed genes were identified in the high emperature sensitive cultivars, of which 22 genes showed upregulated expression and 13 genes were downregulated; 32 DEGs were detected in the high temperature resistant cultivars, including 17 upregulated genes and 15 downregulated genes. It is worth noting that 29 genes in both cultivars showed the same expression pattern under high temperature stress, including 17 upregulated and 12 downregulated genes; However, the sensitive cultivars has five upregulated genes (AUR62002054, AUR62003741, AUR62009231, AUR62009457, AUR62022819), suggesting that these genes may play a negative regulatory role in heat tolerance (Fig. 7A, Table S5).

Fig. 7
figure 7

Expression profile of AAT gene in quinoa under temperature stress. Red in the heatmap indicates high expression level, while blue indicates low expression level. (A) Expression profile of AAT gene in quinoa under high temperature stress. CK group represents normal temperature of heat resistant type. TK group represents heat damage of heat resistant type. CN group represents normal temperature of thermal sensitive type. TN group represents normal temperature of thermal sensitive type. (B) Expression profile of AAT gene in quinoa under low temperature stress. AY1 group represents freeze injury of low temperature sensitive type. BY1 group represents cold damage of low temperature sensitive type. CY1 group represents normal temperature of low temperature sensitive type. AY2 group represents freeze injury of low temperature resistant type. BY2 group represents cold damage of low temperature resistant type. CY2 group represents normal temperature of low temperature resistant type

Expression pattern of CqAAT gene under low temperature stress

Based on transcriptome sequencing data, this study analyzed the dynamic regulatory characteristics of AAT gene family members in quinoa under cold stress (5℃) and freeze injury (-2℃). Among the low temperature sensitive cultivars, there were 10 significantly differentially expressed genes under cold damage, of which 8 were upregulated and 2 were downregulated; There are 12 significantly different genes under freeze injury, of which 7 are upregulated and 5 are downregulated. In the low-temperature tolerant cultivars, there are 16 significantly differentially expressed genes under cold damage, of which 4 are upregulated and 12 are downregulated; There are 13 significantly different genes under freeze injury, of which 7 are upregulated and 6 are downregulated. Both cultivars showed significant upregulation of AUR62041196, AUR62043040, and AUR62043406 under cold and freeze conditions, indicating that these genes play a critical role in quinoa’s response to low temperature stress. The expression levels of AUR62017169 in the low temperature tolerant cultivar increased under cold damage conditions, showing significant differences compared to the low temperature sensitive cultivar. However, their expression levels were consistent under freeze injury conditions, indicating that this genes is key genes in regulating cold stress in quinoa. AUR62017169 belongs to the AAP subfamily, AUR62043040 belongs to the ATLb subfamily, and AUR62041196 and AUR62043406 belong to the CAT subfamily. The GAT subfamily genes AUR62031750 and AUR62023955 exhibit typical freeze injury specific activation patterns. Their expression levels did not show significant differences compared to the control group under cold damage stress, but were significantly upregulated under cold damage conditions. It is speculated that these two genes play a key role in quinoa’s response to freeze injury stress (Fig. 7B, Table S6).

Analysis of the expression pattern of CqAAT gene in different fertilizers and different dosages

Expression pattern analysis of CqAAT gene under different nitrogen applications

To determine the role of CqAAPs in regulating quinoa’s resistance to various nutritional stresses, we studied their transcriptional responses under these conditions. Nitrogen is a key element for plant growth, involved in the synthesis of proteins, nucleic acids, and chlorophyll, and is crucial for plant growth, development, and photosynthesis. In red quinoa, four nitrogen responsive genes (AUR62012331, AUR62002054, AUR62038452, AUR62035199) were significantly upregulated under low nitrogen stress; However, high nitrogen stress significantly upregulated AUR62023955 and AUR62035199. It is worth noting that AUR62023955 shows high expression under both low and high nitrogen stress, suggesting that it may be a core regulatory hub for regulating nitrogen balance. In sharp contrast to red quinoa, no significantly differentially expressed genes were detected in white quinoa under nitrogen stress, indicating that red quinoa may be a high nitrogen sensitive strain, while white quinoa is a high nitrogen tolerant strain (Fig. 8A, Table S7).

Fig. 8
figure 8

Expression profile of AAT gene in quinoa under different fertilizers and at different dosages. R represents red quinoa and W represents white quinoa. Red in the heatmap indicates high expression level, while blue indicates low expression level. (A) Expression profile of AAT gene in quinoa under different nitrogen applications. LN represents low nitrogen; CK represents control group; HN represents high nitrogen. (B) Expression profile of AAT gene in quinoa under different phosphorus applications. LP represents low phosphorus; CK represents control group; HP represents high phosphorus. (C) Expression profile of AAT gene in quinoa under different potassium applications. LK represents low potassium; CK represents control group; HK represents high potassium

Expression pattern analysis of CqAAT gene under different phosphorus applications

Phosphorus is a key element for plant growth and development, involved in energy metabolism, nucleic acid synthesis, and signal transduction. This study reveals the expression patterns of the AAT gene family in red and white quinoa under phosphorus stress. Under low phosphorus stress, three genes AUR62002054, AUR62023955 and AUR62035199 were significantly upregulated in red quinoa; There is a gene AUR62043362 significantly downregulated in white quinoa. Under high phosphorus stress, there were no significant changes in the genes of both strains. AUR62002054 belongs to the ANT subfamily, AUR62023955 belongs to the GAT subfamily, AUR62035199 belongs to the AAP subfamily, and AUR62043362 belongs to the CAT subfamily (Fig. 8B, Table S8).

Expression pattern analysis of CqAAT gene under different potassium applications

Potassium is a key mineral for plant growth, essential for maintaining cell osmotic pressure, regulating stomatal opening and closing, promoting photosynthesis, and enhancing plant stress resistance. Based on transcriptome sequencing results, this study found that there were 7 significantly upregulated genes in red quinoa under low potassium conditions, including AUR62012331, AUR62002054, AUR62038452, AUR62023139, AUR62023955, AUR62035199, AUR62003741; AUR62027047 was significantly downregulated in expression; Under high potassium conditions, AUR62002054 and AUR62023955 were significantly upregulated, while AUR62027047 was significantly downregulated. In white quinoa, AUR62023955 was significantly upregulated and AUR62009231 was significantly downregulated under high potassium conditions. There are significantly more genes in response to different potassium fertilizer application rates in red quinoa, indicating that this cultivar is more sensitive to potassium fertilizer (Fig. 8C, Table S9).

Discussion

Quinoa, as a fully nutritious crop, has important economic value, but the emergence of adverse environments often leads to a decrease in quinoa yield and quality. Therefore, it is of great significance to explore and study some key stress resistant genes for the future cultivation of excellent quinoa varieties, among which transcription factors, as one of the most studied transcription regulatory factors, play an important role in regulating plant response to external environmental stress. In recent years, due to the availability of the whole genome sequence of quinoa, research on the quinoa gene family has been reported, such as VOZ [25], GST [26], WRKY [27], HAK [28]、ACS [29]etc., but there have been no reports on the quinoa AAT gene family. Amino acid transporters (AATs) play important roles in plants, responsible for transmembrane transport of amino acids and participating in multiple physiological processes of plant growth and development, including long-distance amino acid transport, absorption of amino acids from soil, and response to pathogens and abiotic stress [30]. So far, AAT has been recognized and thoroughly studied in many plant species, with 63 AAT genes in Arabidopsis [5], 84 in Phaseolus vulgaris [6], 85 in rice [7], 189 in soybean [8], 296 in common wheat [9], 94 in foxtail millet [10], 104 in Tartary buckwheat [11], and 72 in potatoes [12]. This study was based on the whole genome data of quinoa and used bioinformatics methods to screen 147 CqAATs genes, which were unevenly distributed on 18 chromosomes and divided into 12 subgroups, consistent with Arabidopsis. The overall average hydrophilicity index of the AAT gene in quinoa is 0.537, indicating that it belongs to a hydrophobic protein, which is similar to the hydrophobic properties of general transporters. Ten genes lack introns, which may be closely related to their specific regulation or function. The number of AAT genes on chromosome 1 is the highest, which may be related to gene density and chromosome size. Our analysis found that 121 genes are located on the plasma membrane, and so far all AtAAPs genes are located on the plasma membrane. The entry and exit of amino acids in various organelles are mediated by corresponding transporters and during the transport of amino acids, it is often accompanied by the transmembrane transport of hydrogen ions [31]. Quinoa is a heterozygous tetraploid [32]. Previous studies on orthologs and collinearity have found that it is not simply obtained by doubling the A and B subgenomes provided by donors, but involves complex chromosome fusion processes. The collinearity between quinoa chromosomes in this study also proves this.

The expression of amino acid transporters varies significantly in different plant tissues, and these differences are closely related to their functions in plant growth, development, and nutrient allocation. This study downloaded the gene expression data of various tissues of yellow quinoa and white quinoa from the NCBI database, and the results showed that the AAT gene was expressed in a tissue-specific manner in quinoa, indicating that different AAT proteins may have different functions, and the AAT gene expression was highest in the stem, followed by the roots and leaves, with less gene expression in the fruit. AUR62023139 is highly expressed in all tissues, indicating that these genes are crucial for the overall growth and development of plants. Previous studies have found that tyrosine transporters of the AAT gene family in rice are highly expressed in black rice, but not transcribed in non black genotypes, which are related to regulating grain color formation [33]. Therefore, we analyzed the expression patterns of the AAT gene in the flowers and fruits of two different colored quinoa and found that the expression patterns were basically the same in both quinoa, indicating that the AAT gene was highly conserved in the evolutionary process. Through heat maps, we found that the expression level of AUR62023955 in white quinoa flowers was significantly higher than that in yellow quinoa, and the expression level of AUR62038452 in yellow quinoa fruits was significantly higher than that in white quinoa. The tissue specificity of gene expression may be due to the fact that these genes originate from different subgroups, and the AAT genes of different subgroups may mediate different physiological processes. For example, we found that genes with higher root expression belong to the AAP and ProT subgroups, and the first reported amino acid transporter protein in plants is encoded by the AAPs family genes [34]. AAPs are involved in processes such as root absorption of plant amino acids, xylem loading, phloem loading, exchange between xylem and phloem, and seed loading, affecting the growth and development of plant roots, stems, and leaves [19]. For example, AtAAP1 and AtAAP5 are responsible for root amino acid absorption in Arabidopsis [3536]. Arabidopsis thaliana At AUX1 expression can promote lateral root formation in roots, while rice OsAUX2 and OsAUX5 are specifically expressed in young roots and may be involved in root growth and development processes [37]. In addition, proline as a free radical scavenger, can regulate the redox homeostasis of plant cells and participate in the defense response of plants against biological stress [38]. AtProT2 is mainly expressed in the root epidermis and cortex, and plays a role in the absorption of compatible solutes proline and glycine betaine [39]. AtProT1 is expressed in vascular tissues such as roots, stems, leaves, and flowers, mainly in stems, suggesting that it may play an important role in long-distance transport of proline. AtProT2 is highly expressed in root epidermal and cortical cells, and is extensively induced in leaves under plant stress conditions, suggesting that it may play an important role in plant stress resistance [3940]. Our research found that genes with high expression levels in the flower part of quinoa belong to the ATLa subfamily, while in Arabidopsis, the AtLHT subfamily is mainly responsible for the transport and distribution of lysine and histidine in flower organs, which in turn affects flower development [17, 41]. The genes AUR62038452 and AUR62043362 with the highest expression levels in quinoa fruit belong to the CqCAT subfamily. CAT genes such as AtCAT2, AtCAT3, and AtCAT4 in Arabidopsis are highly expressed in leaves and seeds [42].

Abiotic stress can induce the expression of plant related genes, improve plant adaptability to external environments, enhance stress resistance, and ultimately affect plant quality and yield. The loss of function of amino acid transporters can lead to growth inhibition in plants under adverse conditions, such as reduced efficiency of amino acid allocation and utilization under drought or salt stress, which in turn affects the plant’s stress resistance ability [36]. In Arabidopsis, the expression of the AtAAPs family is inhibited under osmotic pressure stress, while the expression levels of AtAAP4 and AtAAP6, as well as McAAT2 in the ice leafed sunflower (Mesembryanthamum crystallized), are downregulated under salt stress [43]. In addition, the expression levels of AtProT2, McAAT1, and HvProT significantly increased under salt stress conditions, and some genes in the rice OsAAT family were significantly down regulated or up-regulated under salt and drought stress, including OsANT3, OsAAP6, and OsAAP11, indicating that amino acid transporters may play a key role in abiotic stress [7, 44, 45]. Studies have shown that CAT and AAP1 are functionally conserved, and they participate in plant defense responses to abiotic stress by synergistically regulating the transmembrane transport of cationic amino acids. The promoter region of cotton CAT contains cis acting elements related to abiotic stress, and the expression of these genes is significantly affected by treatments such as drought, cold, and salt [46]. Overexpression of TaAAP1 in wheat enhances salt tolerance by regulating ethylene production [47]. We found that the expression levels of AUR62031750 and AUR62023955 in the GAT subfamily of quinoa increased under different abiotic stresses and fertilizer treatments. γ-aminobutyric acid (GABA) is a ubiquitous four carbon non reducing amino acid that participates in regulating plant growth, development, and stress resistance as a signaling substance. When plants are subjected to abiotic stress, GABA accumulates in large quantities. GAT located on the cell membrane can transport GABA across the membrane and transport GABA from the extracellular matrix to the cytoplasm. When plants are subjected to abiotic stress, reducing the accumulation of reactive oxygen species (ROS) can alleviate the damage to cells caused by stress and help maintain the homeostasis of GABA in cells. Previous studies have found that under heavy metal aluminum stress, the LcGAT gene is upregulated and can specifically transport GABA. The accumulation of GABA in leaves promotes the closure of plant stomata, thereby improving plant tolerance [4849]. Only one GAT was found in Arabidopsis, which can transfer GABA with high affinity. When subjected to stress or mechanical damage, plants strongly induce AtGAT1 expression to accumulate more ABA [5051]. Xu et al. studied Arabidopsis and demonstrated that the production of guard cell GABA can reduce stomatal opening and transpiration, thereby improving water use efficiency and drought resistance [52]. Overexpression of PeuGAT3 increases the thickness of xylem cell walls in Arabidopsis and poplar, enhances lignin content in xylem tissues, and increases proline accumulation in poplar leaves, all of which may improve the tolerance of desert poplar to salt and drought stress [53].

AATs can promote the transport of amino acids from source organs to sink organs, thereby improving plant nitrogen absorption and utilization efficiency [54]. When PsAAP1 in pea (Pisum sativum L.) is overexpressed in the phloem and embryo, the allocation of amino acid source pool is improved, significantly enhancing the nitrogen absorption efficiency and utilization efficiency of plants [55]. In Arabidopsis, AtAAP2 encodes an amino acid osmolase that participates in the exchange of amino acids between the xylem and phloem, transporting the amino acids absorbed by the root system to the aboveground parts. When the AtAAP2 gene is knocked out, more nitrogen is allocated to the leaves. At different nitrogen fertilizer levels, the chlorophyll content of leaves is generally higher, Rubisco enzyme activity is stronger, and photosynthetic nitrogen use efficiency (NUE) is also higher [56]. In wheat, the TaATLa1 gene is significantly downregulated under nitrogen deficiency conditions, indicating its possible involvement in nitrogen metabolism and source sink transport of amino acids [57]. The expression of GmAAP6 and GmAAP12 in soybeans is significantly upregulated under low nitrogen conditions, which may affect seed protein content by regulating amino acid transport efficiency [58]. In this study, AUR62023955 and AUR62031750 were significantly upregulated under high nitrogen conditions, suggesting their regulation of quinoa’s response to nitrogen, providing targets for genetic improvement of quinoa nitrogen utilization and protein accumulation.

This study found that AAT in quinoa plays unique roles at different stages and parts of its growth and development, and is involved in quinoa’s response to stress and adversity. Providing a research basis for the regulation mechanism of the AAT gene in quinoa under abiotic stress is of great significance for improving quinoa quality and stress resistance.

Materials and methods

Identification and physicochemical characterization of quinoa AAT gene family members

The whole genome sequence of quinoa was downloaded from NCBI (https://www.ncbi.nlm.nih.gov/). The Arabidopsis AAT family protein sequences were downloaded from TAIR (https://www.arabidopsis.org/). The quinoa family genes were characterized by 2 methods: first, the quinoa AAT genes were characterized by a bidirectional BLAST method, using the Arabidopsis protein sequences as seed files, and BLAST analysis was performed in TBTools [59]. Second, AAT protein Hidden Markov Model (HMM) seed files (PF00324 and PF01490) were obtained from the Pfam database (http://pfam.xfam.org/), and the NCBI-CDD (https://www.ncbi.nlm.nih.gov/cdd) database was used to analyze the candidate genes’ Protein sequences were analyzed to exclude candidate proteins that did not contain the AA-permease structural domain and the Aa-trans structural domain and to finalize the quinoa AAT family members. Physicochemical properties of the self-identified protein sequences were analyzed using the ExPASy website (https://www.expasy.org/) [60]. Subcellular localization was performed using the Cell-PLoc2.0 website (http://www.csbio.sjtu.edu.cn/bioinf/Cell-PLoc-2/) [61].

Phylogenetic analysis of quinoa AAT protein family members

Multiple sequence comparison of AAT protein sequences from Arabidopsis thaliana and quinoa was performed using MEGA11 software, and phylogenetic trees were constructed using the neighbor-joining (NJ) method with the parameter Bootstrap set to 1000 times [62]. The evolutionary tree was landscaped using the online tool iTOL (https://itol.embl.de/) [63].

Gene structure and conserved motif analysis of quinoa AAT gene family members

TBtools was utilized to identify introns and exons of the CqAAT gene [59]. Conserved motifs of quinoa AAT family protein sequences were analyzed using MEME online software, and the number of output conserved motifs was set to 20 [64]. The results of the above analysis were integrated and visualized by TBtools software.

Chromosomal localization and covariance analysis of quinoa AAT gene family members

Sequences of Arabidopsis and rice AAT family proteins were downloaded from TAIR (https://www.arabidopsis.org/) and Rice Genome Annotation Project (http://rice.uga.edu/index.shtml). Information on the distribution of the quinoa AAT gene on chromosomes was extracted from the gff3 file of the quinoa genome, and gene duplication events within quinoa species were analyzed by the Multiple Covariance Scanning Toolkit (MCScanX) and visualized by the Advanced Circos tool in TBtools. To explore the co-lineage relationship of AAT genes in quinoa and other selected species, co-lineage analysis maps were constructed using MCScanX and visualized using Multiple SyntenyPlot in TBtools [59].

Analysis of cis acting components

A sequence of 2000 bp upstream of the transcription start site of AAT gene family members was extracted from the quinoa genome database as the promoter region of the genes, and the cis acting elements were predicted and analyzed using the PlantCARE online website (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/) [65]. Use Excel to filter and analyze the types and quantities of cis acting elements, and use a heatmap to display the types and quantities of cis acting elements in the AAT gene [59].

Tissue expression analysis of AAT gene family members in quinoa

Download tissue expression data of various organs and tissues during the growth and development of quinoa from the GEO public database on the NCBI website(https://www.ncbi.nlm.nih.gov/sra)(Search numbers SRP278144, SRP226463). This data includes the roots, stems, leaves, flowers, and fruits of yellow quiona and the flowers and fruits of white quinoa. Obtain the corresponding expression level (FPKM) values of each member of the AAT gene family through screening. Then use the Heatmap tool in TBtool software to draw an organizational expression heatmap [59](Table S10).

Expression analysis of AAT gene in quinoa under different abiotic stress conditions

In order to investigate whether CqAAT gene expression is related to adversity stress, a gene expression heatmap was constructed to analyze the expression of AAT genes in quinoa leaves under different stresses based on quinoa transcriptome data from the previous high temperature, low temperature, drought, and waterlogging stress studies of our group [66,67,68] (Tables S11-S14). Perform differential expression analysis between sample groups using DESeq2 [6970] to obtain a set of differentially expressed genes between two biological conditions. After differential analysis, it is necessary to use the Benjamin Hochberg method to perform multiple hypothesis testing correction on the hypothesis testing probability (P value) to obtain the False Discovery Rate (FDR). The screening criteria for differentially expressed genes are | log2Fold Change|>=1 and FDR < 0.05. The material is planted in a greenhouse at the Yunnan Agricultural University base in Xundian County, Kunming(102°41′E, 25°20′N). Samples were taken during the seedling stage (Six leaf stage) under high temperature, low temperature, and drought stress, while samples were taken during the grain filling stage under waterlogging stress. The sampling site is the leaf and the sample information for different stress treatments is shown in Tables 2, 3, 4 and 5.

Table 2 Samples information on waterlogging stress
Table 3 Samples information on drought stress
Table 4 Samples information on high temperature stress
Table 5 Samples information on low temperature stress

Expression analysis of quinoa AAT gene under different fertilization conditions

In order to investigate whether the expression of CqAAT gene responds to different fertilization conditions, a gene expression heatmap was constructed based on transcriptome data of red quinoa (Dianli-1299) and white quinoa (Dianli-71) leaves published in our research group. The expression of AAT gene in quinoa leaves under different nitrogen, phosphorus, and potassium supply conditions was analyzed [71,72,73](Table S15-S17). Perform differential expression analysis between sample groups using DESeq2 [6970] to obtain a set of differentially expressed genes between two biological conditions. After differential analysis, it is necessary to use the Benjamin Hochberg method to perform multiple hypothesis testing correction on the hypothesis testing probability (P value) to obtain the False Discovery Rate (FDR). The screening criteria for differentially expressed genes are | log2Fold Change|>=1 and FDR < 0.05. The material is planted in a greenhouse at the Yunnan Agricultural University base in Xundian County, Kunming(102°41′E, 25°20′N). Sampling at the six leaf stage, with the sampling site being the leaves. The sample information for different fertilizers and dosages is shown in Table 6.

Table 6 Sample information on different fertilization conditions

Conclusions

In this study, we identified 147 AAT genes in the quinoa genome sequence. Multiple analyses, including the construct phylogenetic tree, gene exon–intron structure and motif, indicated that the CqAAT gene family was divided into 12 subgroups. Organizational expression analysis revealed significant differences in the expression of the CqAAT gene, with less expression in fruits and more expression in roots, stems, leaves, and flowers. The expression analysis of two genes AUR62031750 and AUR62023955 in the CqGAT subfamily was found to be upregulated under various abiotic stresses, including high temperature, low temperature, waterlogging, and drought, as well as different nitrogen, phosphorus, and potassium fertilizer application rates. This indicates that they play a role in the growth and development regulation of CqAATs under abiotic stress. Quinoa is rich in amino acids, and in-depth research on the function of CqAATs is needed in the future, which is of great significance for improving the quality and breeding efficiency of quinoa, as well as enhancing its stress resistance and adaptability.

Data availability

The datasets generated and/or analyzed during the current study are available in the [NCBI] (National Center for Biotechnology Information) repository, [PRJNA939366],[PRJNA1039692 ],[PRJNA857812], [PRJNA946047], [PRJNA839110].

References

  1. Bhargava A, Shukla S, Ohri D. Genetic variability and interrelationship among various morphological and quality traits in quinoa (Chenopodium quinoa willd). Field Crops Res. 2007;101(1):104–16. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.fcr.2006.10.001.

    Article  Google Scholar 

  2. Vega-Gálvez A, Miranda M, Vergara J, Uribe E, Puente L, Martínez EA. Nutrition facts and functional potential of quinoa (Chenopodium quinoa willd.), an ancient Andean grain: a review. J Sci Food Agric. 2010;90(15):2541–7. https://doiorg.publicaciones.saludcastillayleon.es/10.1002/jsfa.4158.

    Article  CAS  PubMed  Google Scholar 

  3. Le L, Gong X, An Q, Xiang D, Zou L, Peng L, Wu X, Tan M, Nie Z, Wu Q, Zhao G, Wan Y. Quinoa sprouts as potential vegetable source: Nutrient composition and functional contents of different quinoa sprout varieties. Food Chem. 2021;357:129752. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.foodchem.2021.129752. Advance online publication.

    Article  CAS  PubMed  Google Scholar 

  4. Pathan S, Siddiqui RA. Nutritional Composition and Bioactive Components in Quinoa (Chenopodium quinoa Willd). Greens: Rev Nutrients. 2022;14(3):558. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/nu14030558.

    Article  CAS  Google Scholar 

  5. Rentsch D, Schmidt S, Tegeder M. Transporters for uptake and allocation of organic nitrogen compounds in plants. FEBS Lett. 2007;581(12):2281–9. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.febslet.2007.04.013.

    Article  CAS  PubMed  Google Scholar 

  6. Nanjareddy K, Guerrero-Carrillo MF, Lara M, Arthikala MK. Genome-wide identification and comparative analysis of the Amino Acid Transporter (AAT) gene family and their roles during Phaseolus vulgaris symbioses. Funct Integr Genom. 2024;24(2):47. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s10142-024-01331-0.

    Article  CAS  Google Scholar 

  7. Zhao H, Ma H, Yu L, Wang X, Zhao J. Genome-wide survey and expression analysis of amino acid transporter gene family in rice (Oryza sativa L). PLoS ONE. 2012;7(11):e49210. https://doiorg.publicaciones.saludcastillayleon.es/10.1371/journal.pone.0049210.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Cheng L, Yuan HY, Ren R, Zhao SQ, Han YP, Zhou QY, Ke DX, Wang YX, Wang L. Genome-Wide Identification, Classification, and Expression Analysis of Amino Acid Transporter Gene Family in Glycine Max. Front Plant Sci. 2016;7:515. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fpls.2016.00515.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Tian R, Yang Y, Chen M. Genome-wide survey of the amino acid transporter gene family in wheat (Triticum aestivum L.): Identification, expression analysis and response to abiotic stress. Int J Biol Macromol. 2020;162:1372–87. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.ijbiomac.2020.07.302.

    Article  CAS  PubMed  Google Scholar 

  10. Yang Y, Chai Y, Liu J, Zheng J, Zhao Z, Amo A, Cui C, Lu Q, Chen L, Hu YG. Amino acid transporter (AAT) gene family in foxtail millet (Setaria italica L.): widespread family expansion, functional differentiation, roles in quality formation and response to abiotic stresses. BMC Genomics. 2021;22(1):519. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12864-021-07779-9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Yang Y, Wang X, Zheng J, Men Y, Zhang Y, Liu L, Han Y, Hou S, Sun Z. Amino acid transporter (AAT) gene family in Tartary buckwheat (Fagopyrum tataricum L. Gaertn.): Characterization, expression analysis and functional prediction. Int J Biol Macromol. 2022;217:330–44. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.ijbiomac.2022.07.059.

    Article  CAS  PubMed  Google Scholar 

  12. Ma H, Cao X, Shi S, Li S, Gao J, Ma Y, Zhao Q, Chen Q. Genome-wide survey and expression analysis of the amino acid transporter superfamily in potato (Solanum tuberosum L). Plant Physiol Biochem. 2016;107:164–77. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.plaphy.2016.06.007.

    Article  CAS  PubMed  Google Scholar 

  13. A WNF, Bruno André b A, D. R.A, S. K.A, M. T.& A, K. B., et al. Amino acid transport in plants. Trends Plantence. 1998;3(5):188–95. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/S1360-1385(98)01231-X.

    Article  Google Scholar 

  14. Saier MH Jr, Yen MR, Noto K, Tamang DG, Elkan C. The Transporter Classification Database: recent advances. Nucleic Acids Res. 2009;37:D274–8. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/nar/gkn862.

    Article  CAS  PubMed  Google Scholar 

  15. Okumoto S, Schmidt R, Tegeder M, Fischer WN, Rentsch D, Frommer WB, Koch W. High affinity amino acid transporters specifically expressed in xylem parenchyma and developing seeds of Arabidopsis. J Biol Chem. 2002;277(47):45338–46. https://doiorg.publicaciones.saludcastillayleon.es/10.1074/jbc.M207730200.

    Article  CAS  PubMed  Google Scholar 

  16. Chen L, Bush DR. LHT1, a lysine- and histidine-specific amino acid transporter in arabidopsis. Plant Physiol. 1997;115(3):1127–34. https://doiorg.publicaciones.saludcastillayleon.es/10.1104/pp.115.3.1127.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Foster J, Lee YH, Tegeder M. Distinct expression of members of the lht amino acid transporter family in flowers indicates specific roles in plant reproduction. Sex Plant Reprod. 2008;21(2):143–52. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s00497-008-0074-z.

    Article  CAS  Google Scholar 

  18. Guo N, Zhang S, Gu M, Xu G. Function, transport, and regulation of amino acids: what is missing in rice? Crop J. 2021;9(3):530–42. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.cj.2021.04.002.

    Article  Google Scholar 

  19. Ji Y, Huang W, Wu B, Fang Z, Wang X. The amino acid transporter AAP1 mediates growth and grain yield by regulating neutral amino acid uptake and reallocation in Oryza sativa. J Exp Bot. 2020;71(16):4763–77. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/jxb/eraa256.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Yao X, Nie J, Bai R, Sui X. Amino Acid Transporters in Plants: Identification and Function. Plants (Basel Switzerland). 2020;9(8):972. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/plants9080972.

    Article  CAS  PubMed  Google Scholar 

  21. Forde BG, Lea PJ. Glutamate in plants: metabolism, regulation, and signalling. J Exp Bot. 2007;58(9):2339–58. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/jxb/erm121.

    Article  CAS  PubMed  Google Scholar 

  22. Lehmann S, Funck D, Szabados L, Rentsch D. Proline metabolism and transport in plant development. Amino Acids. 2010;39(4):949–62. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s00726-010-0525-3.

    Article  CAS  PubMed  Google Scholar 

  23. Na GUO, Dong XUE, Wei, ZHANG, et al. Overexpression of gmprot1 and gmprot2 increases tolerance to drought and salt stresses in transgenic arabidopsis. J Integr Agric. 2016;15(8):1727–43. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/S2095-3119(15)61288-6.

    Article  CAS  Google Scholar 

  24. Breitkreuz KE, Shelp BJ, Fischer WN, Schwacke R, Rentsch D. Identification and characterization of GABA, proline and quaternary ammonium compound transporters from Arabidopsis thaliana. FEBS Lett. 1999;450(3):280–4. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/s0014-5793(99)00516-5.

    Article  CAS  PubMed  Google Scholar 

  25. Shi P, Jiang R, Li B, Wang D, Fang D, Yin M, Yin M, Gu M. Genome-Wide Analysis and Expression Profiles of the VOZ Gene Family in Quinoa (Chenopodium quinoa). Genes. 2022;13(10):1695. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/genes13101695.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Tiwari S, Vaish S, Singh N, Basantani M, Bhargava A. (2024). Correction: Genome-wide identification and characterization of glutathione S-transferase gene family in quinoa (Chenopodium quinoa Willd.). 3 Biotech, 14(11), 265. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s13205-024-04076-6

  27. Yue H, Chang X, Zhi Y, Wang L, Xing G, Song W, Nie X. Evolution and Identification of the WRKY Gene Family in Quinoa (Chenopodium quinoa). Genes. 2019;10(2):131. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/genes10020131.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Chen Y, Lin Y, Zhang S, Lin Z, Chen S, Wang Z. Genome-Wide Identification and Characterization of the HAK Gene Family in Quinoa (Chenopodium quinoa Willd.) and Their Expression Profiles under Saline and Alkaline Conditions. Plants. 2023;12(21):3747. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/plants12213747.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Yin L, Zhang X, Gao A, Cao M, Yang D, An K, Guo S, Yin H. Genome-Wide Identification and Expression Analysis of 1-Aminocyclopropane-1-Carboxylate Synthase (ACS) Gene Family in Chenopodium quinoa. Plants. 2023;12(23):4021. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/plants12234021.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Margheritis E, Imperiali FG, Cinquetti R, Vollero A, Terova G, Rimoldi S, Girardello R, Bossi E. Amino acid transporter B(0)AT1 (slc6a19) and ancillary protein: impact on function. Pflug Arch: Eur J Physiol. 2016;468(8):1363–74. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s00424-016-1842-5.

    Article  CAS  Google Scholar 

  31. Tegeder M, Rentsch D. Uptake and partitioning of amino acids and peptides. Mol Plant. 2010;3(6):997–1011. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/mp/ssq047.

    Article  CAS  PubMed  Google Scholar 

  32. Jarvis DE, Ho YS, Lightfoot DJ, Schmöckel SM, Li B, Borm TJ, Ohyanagi H, Mineta K, Michell CT, Saber N, Kharbatia NM, Rupper RR, Sharp AR, Dally N, Boughton BA, Woo YH, Gao G, Schijlen EG, Guo X, Momin AA, Tester M. The genome of Chenopodium quinoa. Nature. 2017;542(7641):307–12. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/nature21370.

    Article  CAS  PubMed  Google Scholar 

  33. Li B, Jia Y, Xu L, Zhang S, Long Z, Wang R, Guo Y, Zhang W, Jiao C, Li C, Xu Y. Transcriptional convergence after repeated duplication of an amino acid transporter gene leads to the independent emergence of the black husk/pericarp trait in barley and rice. Plant Biotechnol J. 2024;22(5):1282–98. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/pbi.14264.

    Article  CAS  PubMed  Google Scholar 

  34. Lee YH, Foster J, Chen J, Voll LM, Weber AP, Tegeder M. AAP1 transports uncharged amino acids into roots of Arabidopsis. Plant journal: cell Mol biology. 2007;50(2):305–19. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.1365-313X.2007.03045.x.

    Article  CAS  Google Scholar 

  35. Svennerstam H, Ganeteg U, Näsholm T. Root uptake of cationic amino acids by Arabidopsis depends on functional expression of amino acid permease 5. New Phytol. 2008;180(3):620–30. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.1469-8137.2008.02589.x.

    Article  CAS  PubMed  Google Scholar 

  36. Zhang L, Tan Q, Lee R, Trethewy A, Lee YH, Tegeder M. Altered xylem-phloem transfer of amino acids affects metabolism and leads to increased seed yield and oil content in Arabidopsis. Plant Cell. 2010;22(11):3603–20. https://doiorg.publicaciones.saludcastillayleon.es/10.1105/tpc.110.073833.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Ugartechea-Chirino Y, Swarup R, Swarup K, Péret B, Whitworth M, Bennett M, Bougourd S. The AUX1 LAX family of auxin influx carriers is required for the establishment of embryonic root cell organization in Arabidopsis thaliana. Ann Botany. 2010;105(2):277–89. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/aob/mcp287.

    Article  CAS  Google Scholar 

  38. Szabados L, Savouré A. Proline: a multifunctional amino acid. Trends Plant Sci. 2010;15(2):89–97. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.tplants.2009.11.009.

    Article  CAS  PubMed  Google Scholar 

  39. Grallath S, Weimar T, Meyer A, Gumy C, Suter-Grotemeyer M, Neuhaus JM, Rentsch D. The AtProT family. Compatible solute transporters with similar substrate specificity but differential expression patterns. Plant Physiol. 2005;137(1):117–26. https://doiorg.publicaciones.saludcastillayleon.es/10.1104/pp.104.055079.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Lehmann S, Gumy C, Blatter E, Boeffel S, Fricke W, Rentsch D. In planta function of compatible solute transporters of the AtProT family. J Exp Bot. 2011;62(2):787–96. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/jxb/erq320.

    Article  CAS  PubMed  Google Scholar 

  41. Lee YH, Tegeder M. Selective expression of a novel high-affinity transport system for acidic and neutral amino acids in the tapetum cells of Arabidopsis flowers. Plant journal: cell Mol biology. 2004;40(1):60–74. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.1365-313X.2004.02186.x.

    Article  CAS  Google Scholar 

  42. Yang H, Krebs M, Stierhof YD, Ludewig U. Characterization of the putative amino acid transporter genes AtCAT2, 3 &4: the tonoplast localized AtCAT2 regulates soluble leaf amino acids. J Plant Physiol. 2014;171(8):594–601. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.jplph.2013.11.012.

    Article  CAS  PubMed  Google Scholar 

  43. Rentsch D, Hirner B, Schmelzer E, Frommer WB. Salt stress-induced proline transporters and salt stress-repressed broad specificity amino acid permeases identified by suppression of a yeast amino acid permease-targeting mutant. Plant Cell. 1996;8(8):1437–46. https://doiorg.publicaciones.saludcastillayleon.es/10.1105/tpc.8.8.1437.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Ueda A, Shi W, Sanmiya K, Shono M, Takabe T. Functional analysis of salt-inducible proline transporter of barley roots. Plant Cell Physiol. 2001;42(11):1282–9. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/pcp/pce166.

    Article  CAS  PubMed  Google Scholar 

  45. Popova OV, Dietz KJ, Golldack D. Salt-dependent expression of a nitrate transporter and two amino acid transporter genes in Mesembryanthemum crystallinum. Plant Mol Biol. 2003;52(3):569–78. https://doiorg.publicaciones.saludcastillayleon.es/10.1023/a:1024802101057.

    Article  CAS  PubMed  Google Scholar 

  46. Chen X, Wu Z, Yin Z, Zhang Y, Rui C, Wang J, Malik WA, Lu X, Wang D, Wang J, Guo L, Wang S, Zhao L, Qaraevna Z, Chen B, Wang C, X., Ye W. Comprehensive genomic characterization of cotton cationic amino acid transporter genes reveals that GhCAT10D regulates salt tolerance. BMC Plant Biol. 2022;22(1):441. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12870-022-03829-w.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Wang K, Zhai M, Cui D, Han R, Wang X, Xu W, Qi G, Zeng X, Zhuang Y, Liu C. Genome-Wide Analysis of the Amino Acid Permeases Gene Family in Wheat and TaAAP1 Enhanced Salt Tolerance by Accumulating Ethylene. Int J Mol Sci. 2023;24(18):13800. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/ijms241813800.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Li L, Dou N, Zhang H, Wu C. The versatile GABA in plants. Plant Signal Behav. 2021;16(3):1862565. https://doiorg.publicaciones.saludcastillayleon.es/10.1080/15592324.2020.1862565.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Hu L, Fan R, Wang P, Hao Z, Yang D, Lu Y, Shi J, Chen J. Identification, Phylogenetic and Expression Analyses of the AAAP Gene Family in Liriodendron chinense Reveal Their Putative Functions in Response to Organ and Multiple Abiotic Stresses. Int J Mol Sci. 2022;23(9):4765. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/ijms23094765.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Meyer A, Eskandari S, Grallath S, Rentsch D. AtGAT1, a high affinity transporter for gamma-aminobutyric acid in Arabidopsis thaliana. J Biol Chem. 2006;281(11):7197–204. https://doiorg.publicaciones.saludcastillayleon.es/10.1074/jbc.M510766200.

    Article  CAS  PubMed  Google Scholar 

  51. Yang Y, Chai Y, Liu J, Zheng J, Zhao Z, Amo A, Cui C, Lu Q, Chen L, Hu YG. Amino acid transporter (AAT) gene family in foxtail millet (Setaria italica L.): widespread family expansion, functional differentiation, roles in quality formation and response to abiotic stresses. BMC Genomics. 2021;22(1):519. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12864-021-07779-9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Xu B, Long Y, Feng X, Zhu X, Sai N, Chirkova L, Betts A, Herrmann J, Edwards EJ, Okamoto M, Hedrich R, Gilliham M. Author Correction: GABA signalling modulates stomatal opening to enhance plant water use efficiency and drought resilience. Nat Commun. 2024;15(1):1737. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41467-024-46158-2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Bai X, Xu J, Shao X, Luo W, Niu Z, Gao C, Wan D. A Novel Gene Coding γ-Aminobutyric Acid Transporter May Improve the Tolerance of Populus euphratica to Adverse Environments. Front Plant Sci. 2019;10:1083. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fpls.2019.01083.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Dong K, Ye Z, Hu F, Shan C, Wen D, Cao J. Improvement of plant quality by amino acid transporters: A comprehensive review. Plant Physiol Biochem. 2024;215:109084. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.plaphy.2024.109084.

    Article  CAS  PubMed  Google Scholar 

  55. Tegeder M, Tan Q, Grennan AK, Patrick JW. Amino acid transporter expression and localisation studies in pea (Pisum sativum). Funct Plant Biol. 2007;34(11):1019–28. https://doiorg.publicaciones.saludcastillayleon.es/10.1071/FP07107.

    Article  CAS  PubMed  Google Scholar 

  56. Wan Y, King R, Mitchell RAC, Hassani-Pak K, Hawkesford MJ. Spatiotemporal expression patterns of wheat amino acid transporters reveal their putative roles in nitrogen transport and responses to abiotic stress. Sci Rep. 2017;7(1):5461. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41598-017-04473-3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Chen H, Liu Y, Zhang J, Chen Y, Dai C, Tian R, Liu T, Chen M, Yang G, Wang Z, Li H, Cao X, Gao X. Amino acid transporter gene TaATLa1 from Triticum aestivum L. improves growth under nitrogen sufficiency and is down regulated under nitrogen deficiency. Planta. 2022;256(4):65. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s00425-022-03978-0.

    Article  CAS  PubMed  Google Scholar 

  58. Zhang Y, Wang L, Song B-H, Zhang D, Zhang H. Genome-Wide Identification, Characterization, and Expression Analysis of the Amino Acid Permease Gene Family in Soybean. Agronomy. 2024;14(1):52. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/agronomy14010052.

    Article  CAS  Google Scholar 

  59. Chen C, Wu Y, Li J, Wang X, Zeng Z, Xu J, Liu Y, Feng J, Chen H, He Y, Xia R. TBtools-II: A one for all, all for one bioinformatics platform for biological big-data mining. Mol Plant. 2023;16(11):1733–42. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.molp.2023.09.010.

    Article  CAS  PubMed  Google Scholar 

  60. Duvaud S, Gabella C, Lisacek F, Stockinger H, Ioannidis V, Durinx C. Expasy, the Swiss Bioinformatics Resource Portal, as designed by its users. Nucleic Acids Res. 2021;49(W1):W216–27. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/nar/gkab225.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Chou KC, Shen HB. Cell-PLoc: a package of Web servers for predicting subcellular localization of proteins in various organisms. Nat Protoc. 2008;3(2):153–62. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/nprot.2007.494.

    Article  CAS  PubMed  Google Scholar 

  62. Tamura K, Stecher G, Kumar S. Mol Biol Evol. 2021;38(7):3022–7. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/molbev/msab120. MEGA11: Molecular Evolutionary Genetics Analysis Version 11.

  63. Letunic I, Bork P. Interactive Tree of Life (iTOL) v6: recent updates to the phylogenetic tree display and annotation tool. Nucleic Acids Res. 2024;52:W78–82. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/nar/gkae268.

    Article  PubMed  PubMed Central  Google Scholar 

  64. Bailey TL, Johnson J, Grant CE, Noble WS. The MEME Suite. Nucleic Acids Res. 2015;43W1:W39–49. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/nar/gkv416.

    Article  CAS  Google Scholar 

  65. Lescot M, Déhais P, Thijs G, Marchal K, Moreau Y, Van de Peer Y, Rouzé P, Rombauts S. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002;30(1):325–7. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/nar/30.1.325.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Huan X, Li L, Liu Y, Kong Z, Liu Y, Wang Q, Liu J, Zhang P, Guo Y, Qin P. Integrating transcriptomics and metabolomics to analyze quinoa (Chenopodium quinoa Willd.) responses to drought stress and rewatering. Front Plant Sci. 2022;13:988861. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fpls.2022.988861.

    Article  PubMed  PubMed Central  Google Scholar 

  67. Xie H, Wang Q, Zhang P, Zhang X, Huang T, Guo Y, Liu J, Li L, Li H, Qin P. Transcriptomic and Metabolomic Analysis of the Response of Quinoa Seedlings to Low Temperatures. Biomolecules. 2022;12(7):977. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/biom12070977.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Xie H, Zhang P, Jiang C, Wang Q, Guo Y, Zhang X, Huang T, Liu J, Li L, Li H, Wang H, Qin P. Combined transcriptomic and metabolomic analyses of high temperature stress response of quinoa seedlings. BMC Plant Biol. 2023;23(1):292. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12870-023-04310-y.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13059-014-0550-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Varet H, Brillet-Guéguen L, Coppée JY, Dillies MA. PLoS ONE. 2016;11(6):e0157022. https://doiorg.publicaciones.saludcastillayleon.es/10.1371/journal.pone.0157022. SARTools: A DESeq2- and EdgeR-Based R Pipeline for Comprehensive Differential Analysis of RNA-Seq Data.

  71. Li H, Wang Q, Huang T, Liu J, Zhang P, Li L, Xie H, Wang H, Liu C, Qin P. Transcriptome and Metabolome Analyses Reveal Mechanisms Underlying the Response of Quinoa Seedlings to Nitrogen Fertilizers. Int J Mol Sci. 2023;24(14):11580. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/ijms241411580.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Wang Q, Guo Y, Huang T, Zhang X, Zhang P, Xie H, Liu J, Li L, Kong Z, Qin P. Transcriptome and Metabolome Analyses Revealed the Response Mechanism of Quinoa Seedlings to Different Phosphorus Stresses. Int J Mol Sci. 2022;23(9):4704. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/ijms23094704.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Huang T, Zhang X, Wang Q, Guo Y, Xie H, Li L, Zhang P, Liu J, Qin P. Metabolome and transcriptome profiles in quinoa seedlings in response to potassium supply. BMC Plant Biol. 2022;22(1):604. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12870-022-03928-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We gratefully acknowledge the fnancial support of the Yunnan Province Academician Workstation(202405AF140012), the Yunnan Province’s “Xing Dian Talent Support Program” (XDYC-CYCX-2022-0031) and Yunnan Provincial Department of Education Science Research Fund Project [2025Y0493].

Funding

This research was funded by the Yunnan Province Academician Workstation(202405AF140012), the Yunnan Province’s “Xing Dian Talent Support Program” (XDYC-CYCX-2022-0031) and Yunnan Provincial Department of Education Science Research Fund Project [2025Y0493].

Author information

Authors and Affiliations

Authors

Contributions

H.L. Writing-Original Draft, Methodology. C.J. Writing-Original Draft, Formal analysis. J.L. Writing-Original Draft, Formal analysis. P.Z. Conceptualization, Writing-Review & Editing. L.L. Formal analysis, Investigation. L.H. Methodology, Visualization. R.L. collected feld samples. X.W. Formal analysis, Investigation. G. J. Investigation. L.Z. Methodology, Visualization. Y. B. Investigation. P. Q. Supervision, Project administration, Funding acquisition.

Corresponding author

Correspondence to Peng Qin.

Ethics declarations

Ethics approval and consent to participate

Experimental research on plants complies with relevant institutional, national, and international guidelines and legislation.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic Supplementary Material

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, H., Jiang, C., Liu, J. et al. Genome-wide identification of the AAT gene family in quinoa and analysis of its expression pattern under abiotic stresses. BMC Genomics 26, 298 (2025). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12864-025-11491-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12864-025-11491-3

Keywords