Transcriptome Sequencing Reveals the Mechanism of Auxin Regulation during Root Expansion in Carrot

Carrot is an important vegetable with roots as the edible organ. A complex regulatory network controls root growth, in which auxin is one of the key players. To clarify the molecular mechanism on auxin regulating carrot root expansion, the growth process and the indole-3-acetic acid (IAA) content in the roots were measured in this experiment. It was found that the rapid expansion period of the root was from 34 to 41 days after sowing and the IAA content was the highest during this period. The root growth then slowed down and the IAA levels decreased. Using the transcriptome sequencing database, we analyzed the expression of IAA-metabolism-related genes and found that the expression of most of the IAA synthesis genes, catabolism genes, and genes related to signal transduction was consistent with the changes in IAA content during root expansion. Among them, a total of 31 differentially expressed genes (DEGs) were identified, including 10 IAA synthesis genes, 8 degradation genes, and 13 genes related to signal transduction. Analysis of the correlations between the DEGs and IAA levels showed that the following genes were closely related to root development: three synthesis genes, YUCCA10 (DCAR_012429), TAR2 (DCAR_026162), and AMI1 (DCAR_003244); two degradation genes, LPD1 (DCAR_023341) and AACT1 (DCAR_010070); and five genes related to signal transduction, IAA22 (DCAR_012516), IAA13 (DCAR_012591), IAA27 (DCAR_023070), IAA14 (DCAR_027269), and IAA7 (DCAR_030713). These results provide a reference for future studies on the mechanism of root expansion in carrots.


Introduction
Carrot (Daucus carota var.sativa DC.) is a biennial plant with nutritious fleshy roots that are widely consumed.Rich in carotenoids, vitamins, lignin, and a variety of other nutritional elements, carrots can enhance immunity, prevent cancers and other diseases, and lower cholesterol in humans [1].Since roots are important edible organs in carrots, it is of great significance to study root expansion.
Fleshy root expansion is a complex process and includes morphogenesis and dry matter accumulation.It is influenced and regulated by multiple factors, including environmental factors and endogenous hormones [2].Previous studies have shown that the endogenous hormones in plants mainly affect the growth and development of the roots by controlling the division and differentiation of root cells; hence, they play an important role in root expansion [3].Auxin, the first type of plant hormone discovered, is widely distributed in plants but mostly concentrated in root tip meristem, cambium, and young seeds.This hormone regulates plant growth and development, cell elongation, and division and promotes xylem and phloem differentiation [4].The formation of sweet potato root tubers is caused by auxin-induced activation of the root cambium and xylem cell proliferation, thus promoting the enlargement of the root tubers [5].Moreover, IAA participates in the formation of lateral roots and the change in root architecture of sweet potato under the low K + stress [6].
Root development is closely related to the activity of secondary cambium, which is determined by plant hormones.Studies have shown that indole-3-acetic acid (IAA) exists in a dynamic balance in plants, and its content and distribution in the roots directly determine cell division [7].In Zingiber, the levels of the hormones such as IAA, gibberellic acid (GA), and jasmonic acid (JA) were significantly increased with increasing rhizome diameter [8].However, the role of IAA in root formation remains contentious.Some researchers believe that IAA affects the activity of lignin and regulates its synthesis in carrot, thus affecting root formation [9,10].
In this study, the carrot roots were investigated at five growth stages.By combining the IAA content and transcriptome sequencing (RNA-Seq) results, the mechanism of IAA regulation expansion was assessed from multiple perspectives.The findings lay the theoretical foundation for the study of carrot root expansion mechanisms in the future, with a view to providing valuable reference data for high-yield breeding of other root crops.

Analysis of the Carrot Root Morphological Indices
The root length and diameter of the carrots gradually increased after sowing (Figure 1A,B).At the early stages, it was mainly the root length that increased, and the rate peaked at 34 days after sowing, reaching 53.6%; however, this decreased rapidly at 41 days after sowing and maintained flat thereafter (Figure 1C).The increase in the root diameter began 34 days after sowing and the growth rate increased rapidly from 34 to 41 days, reaching 130.5% before slowing down (Figure 1D).
proliferation, thus promoting the enlargement of the root tubers [5].Moreover, IAA participates in the formation of lateral roots and the change in root architecture of sweet potato under the low K + stress [6].
Root development is closely related to the activity of secondary cambium, which is determined by plant hormones.Studies have shown that indole-3-acetic acid (IAA) exists in a dynamic balance in plants, and its content and distribution in the roots directly determine cell division [7].In Zingiber, the levels of the hormones such as IAA, gibberellic acid (GA), and jasmonic acid (JA) were significantly increased with increasing rhizome diameter [8].However, the role of IAA in root formation remains contentious.Some researchers believe that IAA affects the activity of lignin and regulates its synthesis in carrot, thus affecting root formation [9,10].
In this study, the carrot roots were investigated at five growth stages.By combining the IAA content and transcriptome sequencing (RNA-Seq) results, the mechanism of IAA regulation expansion was assessed from multiple perspectives.The findings lay the theoretical foundation for the study of carrot root expansion mechanisms in the future, with a view to providing valuable reference data for high-yield breeding of other root crops.

Analysis of the Carrot Root Morphological Indices
The root length and diameter of the carrots gradually increased after sowing (Figure 1A,B).At the early stages, it was mainly the root length that increased, and the rate peaked at 34 days after sowing, reaching 53.6%; however, this decreased rapidly at 41 days after sowing and maintained flat thereafter (Figure 1C).The increase in the root diameter began 34 days after sowing and the growth rate increased rapidly from 34 to 41 days, reaching 130.5% before slowing down (Figure 1D).

Comparison of the IAA Content in Carrot Roots at Different Growth Stages
The content of IAA in the carrot roots was measured at different time points (Figure 2A).The IAA content was higher 34 and 41 days after sowing and peaked at 41 days, followed by a gradual decrease.The IAA content differed significantly between 34 and 41 and between 41 and 48 days.Additionally, the correlation analysis between the growth indicators and IAA content showed that IAA content was negatively correlated with root length (−0.70) and root diameter (−0.83), indicating that IAA is associated with root expansion (Figure 2B).

Comparison of the IAA Content in Carrot Roots at Different Growth Stages
The content of IAA in the carrot roots was measured at different time poin 2A).The IAA content was higher 34 and 41 days after sowing and peaked a followed by a gradual decrease.The IAA content differed significantly between and between 41 and 48 days.Additionally, the correlation analysis between t indicators and IAA content showed that IAA content was negatively correlated length (−0.70) and root diameter (−0.83), indicating that IAA is associated with ro sion (Figure 2B).

Quality Assessment of Transcriptome Sequencing Data
A total of 15 cDNA libraries were constructed in order to analyze the gene e in the carrot roots at different growth stages.After filtering out the contami quences and low-quality sequences from the raw data, the high-quality reads w the range of 91.13-94.52%and of a good sequencing quality (Table S1).The effi the high-quality reads obtained from each library against the carrot referenc (http://plants.ensembl.org/Daucus_carota/Info/Index(accessed on 14 Janua ranged from 86.97 to 88.41%, confirming the reliability of the data and their suit subsequent experimental analysis.

Sample Repeatability Test
The repeatability between biological replicates reflects the reliability of the cipal component analysis provides a visual demonstration of whether the biolog licated samples are sufficiently similar and whether the differences between g sufficiently large.Samples of the carrot roots from different growth stages were together and clearly distinguished from other stages, suggesting that the d within the groups were small and the differences between the groups were (Figure S1).Spearman's correlation coefficient used to evaluate the correlation the biological replicates of all samples (Figure S2) further confirmed the robustn sequencing data.

Analysis of Transcriptome Sequencing 2.3.1. Quality Assessment of Transcriptome Sequencing Data
A total of 15 cDNA libraries were constructed in order to analyze the gene expression in the carrot roots at different growth stages.After filtering out the contaminating sequences and low-quality sequences from the raw data, the high-quality reads were all in the range of 91.13-94.52%and of a good sequencing quality (Table S1).The efficiency of the highquality reads obtained from each library against the carrot reference genome (http:// plants.ensembl.org/Daucus_carota/Info/Index(accessed on 14 January 2022)) ranged from 86.97 to 88.41%, confirming the reliability of the data and their suitability for subsequent experimental analysis.

Sample Repeatability Test
The repeatability between biological replicates reflects the reliability of the data.Principal component analysis provides a visual demonstration of whether the biologically replicated samples are sufficiently similar and whether the differences between groups are sufficiently large.Samples of the carrot roots from different growth stages were clustered together and clearly distinguished from other stages, suggesting that the differences within the groups were small and the differences between the groups were significant (Figure S1).Spearman's correlation coefficient used to evaluate the correlation between the biological replicates of all samples (Figure S2) further confirmed the robustness of the sequencing data.

Real-Time Quantitative PCR Validation of the RNA-Seq Data
Six genes were randomly selected for real-time quantitative PCR (RT-qPCR) to verify the reliability of the gene expression results obtained with RNA-Seq.As shown in Figure S3, the fragments per kilobase of exon model per million mapped fragments (FPKM) values of DCAR_012429, DCAR_030523, DCAR_003708, DCAR_026162, DCAR_016234, and DCAR_018164 were consistent with the fluorescence quantitative results, indicating that the transcriptome data were reliable.

Expression Analysis of IAA-Metabolism-Related Genes
Studies have shown that the aromatic amino acid l-tryptophan (Trp) is the central precursor for IAA biosynthesis in plants [11].Trp-dependent IAA biosynthesis involves multiple parallel pathways, leading to the synthesis of IAA as the final product (Figure 3).The indole-3-pyruvic acid (IPyA) pathway is the main pathway for the synthesis of IAA in plants and consists of a two-step reaction in which Trp is first deamidated by tryptophan aminotransferase of arabidopsis 1 (TAA1) into IPyA.It is then decarboxylated into IAA and catalyzed by a flavin-containing monooxygenase of the YUCCA (YUC) family [12].The Trp pathway and the IPyA pathway share the same first step, and the intermediate product, IPyA, is decarboxylated into IAA in the presence of the key enzymes aldehyde dehydrogenase (ALDH, EC: 1.2.1.3)and indole-3-acetaldehyde oxidase (AAO, EC: 1.2.3.7).In the indole-3-acetaldoxime (IAOx) pathway, the conversion of Trp into its derivative IAOx is catalyzed by the cytochrome P450 monooxygenases CYP79B2 and CYP79B3 [13].Indole-3-acetonitrile (IAN) is then finally converted into IAA by the action of the plant nitrilases family (NITs) [14].In addition, Trp generates the intermediate product indole-3acetamide (IAM) via the IAM pathway, which, in turn, catalyzes the conversion of IAM into IAA by AMIDASE1 (AMI1) [15,16].DCAR_018164 were consistent with the fluorescence quantitative results, indicating that the transcriptome data were reliable.

Expression Analysis of IAA Synthesis Genes
The synthesis of IAA in plants is divided into four pathways: the tryptamine pathway, the IPyA pathway, the IAOx pathway, and the IAM pathway (Figure 3).Analysis of the genes encoding the key enzymes in the IAA synthesis pathway identified 49 genes involved in the process (Table 1).Most of the genes encoding TAA1, a key enzyme of the tryptamine and IPyA pathways, were expressed at higher levels at the early growth stages than at later growth stages, which coincided with the elevated IAA content during early growth.Subsequently, the gene expression level decreased, as did the IAA content, resulting in slower carrot root expansion.Within this, the expression of DCAR_026162 varied significantly during the root expansion.In the tryptamine pathway, the expression changes in most of the genes encoding the key enzyme ALDH were consistent with the changes in the IAA content.Among them, DCAR_011215, DCAR_024255, DCAR_031309, and DCAR_009596 were significantly differentially expressed during root expansion.In the IPyA pathway, there were 19 genes encoding YUC [19], of which DCAR_012429 and DCAR_000740 were significantly differentially expressed during root expansion.In the IAM pathway, eight genes encoding AMI1 had higher expression levels at the early growth stages than the later growth stages.Of these, DCAR_003244 and DCAR_003706 were significantly up-regulated during root expansion, which was consistent with the IAA content.These genes play a crucial role in IAA synthesis and are responsible for the changes in IAA content.However, the expression of IAOx-pathway-related genes was not detected, suggesting that the IAOx pathway may not be involved in the process of carrot root expansion.Therefore, we propose that the main IAA synthesis pathways in carrot root expansion are the tryptamine, IPyA, and IAM pathways.

Expression Analysis of Auxin Degradation Genes
IAA catabolism includes aminobenzoate degradation, nicotinamide metabolism, benzoate degradation, and glycolysis.The key enzymes involved in the IAA catabolic pathway are shown in Figure 3. Analysis of expression of the genes encoding these key enzymes revealed that a total of 31 genes were involved in the IAA catabolic process (Table S2).The DEGs and some genes with annotation of the most likely orthologue genes in the degradation and their expressions level are shown in Table 2.Among these, DCAR_013925 and DCAR_009536, which encode CAT, were significantly down-regulated in H34 vs. H41, which is the opposite trend to that detected for the IAA content.The DAO-encoding genes were all down-regulated during root expansion, among which DCAR_016234 and DCAR_016235 were significantly down-regulated.This is consistent with the results of Porco et al., which suggested that AtDAO1 in Arabidopsis regulates the length of the root hairs and controls the auxin levels in vivo [20].In this experiment, four DAOencoding genes (DCAR_016234, DCAR_002760, DCAR_006175, and DCAR_006176) were up-regulated in H41 vs. H48 to maintain lower IAA content, required for carrot root growth in a suitable range.DCAR_016716 and DCAR_022577, DCAR_024929, DCAR_024930, DCAR_002529, DCAR_026968, and DCAR_029358, which encode OADH, another key enzyme in the catabolic pathway of IAA, and DCAR_019674, DCAR_006315, DCAR_010070, and Daucus_carota_newGene_5365, which encode ACAT, were down-regulated during root expansion, which is consistent with the changes in IAA content.This suggests that they may be involved in glycolysis and decompose excess auxin in the carrot root to maintain auxin concentration homeostasis.

Expression Analysis of Genes Related to Auxin Signal Transduction
The response to auxin in plants is mainly controlled by auxin signal transduction.Typically, the auxin signal transduction pathway is composed of the responses of transport inhibitor response 1 (TIR1), auxin/indole-3-acetic acid (Aux/IAA), and auxin response factor (ARF) (Figure 4).Among these, the TIR1 protein is an important component of the SCF TIR1 complex, which is considered an auxin receptor, and its mutant tir1 exhibits defective hypocotyl and lateral root formation [21].Aux/IAA is a class of auxin-responsive proteins.At low auxin concentrations, the SCF TIR1 complex interacts with Aux/IAA to form the Aux/IAA-ARF transcriptional regulatory element.When its concentration is elevated, auxin mediates the Aux/IAA proteins to undergo ubiquitination degradation, thereby releasing ARF, which promotes the transcription of the corresponding genes and regulates plant growth and development [22].
. Mol The response to auxin in plants is mainly controlled by auxin signal transductio Typically, the auxin signal transduction pathway is composed of the responses transport inhibitor response 1 (TIR1), auxin/indole-3-acetic acid (Aux/IAA), and auxin sponse factor (ARF) (Figure 4).Among these, the TIR1 protein is an important compone of the SCF TIR1 complex, which is considered an auxin receptor, and its mutant tir1 exhib defective hypocotyl and lateral root formation [21].Aux/IAA is a class of auxin-responsi proteins.At low auxin concentrations, the SCF TIR1 complex interacts with Aux/IAA to fo the Aux/IAA-ARF transcriptional regulatory element.When its concentration is elevate auxin mediates the Aux/IAA proteins to undergo ubiquitination degradation, thereby leasing ARF, which promotes the transcription of the corresponding genes and regula plant growth and development [22].A total of 78 genes were identified to be involved in the transcriptional regulation the three proteins TIR1, Aux/IAA, and ARF in this experiment (Table S2) and the DE and key genes are shown in Table 3.The expression of DCAR_028134, the gene encodi TIR1, was significantly higher at the early growth stages (34 d and 41 d after sowing) th at the later growth stages (48 d, 55 d, and 62 d after sowing), which was consistent w the trend in IAA content.There were 33 AUX/IAA-encoding genes identified (Figure Table S3), which were more sensitive to the auxin concentration.The expression of genes (Daucus_carota_newGene_6588, DCAR_010763, DCAR_012516, DCAR_0125 DCAR_018241, DCAR_018447, DCAR_002354, DCAR_023070, DCAR_0242 DCAR_026353, DCAR_027352, DCAR_027269, DCAR_030712, and DCAR_030713) w significantly higher at the early growth stages than the later growth stages, consistent w the trend in the change in IAA content.This indicates that it might be the AUX/IAA p teins that responded to the changes in IAA content and regulated the expression of t transcription factors in the roots.ARF is a key transcription factor in the auxin signal tran duction pathway and works together with the AUX/IAA proteins to regulate the expr sion of auxin-related genes [23].DCAR_015554, DCAR_014788, DCAR_0174 DCAR_020464, DCAR_005400, and DCAR_006053, which encode ARF, showed significa differences in expression, and all of them were higher at the early growth stages.The e A total of 78 genes were identified to be involved in the transcriptional regulation of the three proteins TIR1, Aux/IAA, and ARF in this experiment (Table S2) and the DEGs and key genes are shown in Table 3.The expression of DCAR_028134, the gene encoding TIR1, was significantly higher at the early growth stages (34 d and 41 d after sowing) than at the later growth stages (48 d, 55 d, and 62 d after sowing), which was consistent with the trend in IAA content.There were 33 AUX/IAA-encoding genes identified (Figure S4, Table S3), which were more sensitive to the auxin concentration.The expression of 14 genes (Dau-cus_carota_newGene_6588, DCAR_010763, DCAR_012516, DCAR_012591, DCAR_018241, DCAR_018447, DCAR_002354, DCAR_023070, DCAR_024202, DCAR_026353, DCAR_027352, DCAR_027269, DCAR_030712, and DCAR_030713) was significantly higher at the early growth stages than the later growth stages, consistent with the trend in the change in IAA content.This indicates that it might be the AUX/IAA proteins that responded to the changes in IAA content and regulated the expression of the transcription factors in the roots.ARF is a key transcription factor in the auxin signal transduction pathway and works together with the AUX/IAA proteins to regulate the expression of auxin-related genes [23].DCAR_015554, DCAR_014788, DCAR_017496, DCAR_020464, DCAR_005400, and DCAR_006053, which encode ARF, showed significant differences in expression, and all of them were higher at the early growth stages.The expression of the six ARF genes was consistent with the trend in the change in IAA content.In Arabidopsis, AtARF5, AtARF6, AtARF7, AtARF8, and AtARF19 together form a transcriptional activator that acts on flower, seedling, and root development as well as hypocotyl elongation [24], and the double mutants of arf19 and arf7 exhibit auxin resistance [25], suggesting that the ARF transcription factor controls the auxin response during carrot root expansion.

Discovery of Genes Related to IAA Regulation during Root Expansion
Screening of the differentially expressed genes (DEGs) with false discovery rate (FDR) < 0.01 and fold change (FC) ≥ 1.5 found a total of eight IAA synthesis genes, seven catabolism genes, and twelve genes related to auxin signal transduction that were differentially expressed during root expansion (Table 4).In order to obtain the key genes for auxin regulation of the root expansion, all of the above DEGs were correlated with the IAA content and a total of 11 DEGs with correlation coefficients ≥ 0.8 were identified (Figure 5).Their expression was significantly correlated with the changes in IAA content.In addition, by analyzing the expression patterns of the genes in the IAA metabolic pathways, we found three synthetic genes, YUCCA10 (DCAR_012429), TAR2 (DCAR_026162), and AMI1 (DCAR_003244); one catabolic gene AACT1 (DCAR_010070); and five signaling genes IAA22 (DCAR_012516), IAA13 (DCAR_012591), IAA27 (DCAR_023070), IAA14 (DCAR_027269), and IAA7 (DCAR_030713) were differentially expressed in H41 vs. H48, and one catabolic gene LPD1 (DCAR _023341) was differentially expressed in H55 vs. H62 (Table 5), and these genes were consistent with the expression of genes involved in IAA metabolism in the above analysis.In Arabidopsis, AtARF5, AtARF6, AtARF7, AtARF8, and AtARF19 together form a transcriptional activator that acts on flower, seedling, and root development as well as hypocotyl elongation [24], and the double mutants of arf19 and arf7 exhibit auxin resistance [25], suggesting that the ARF transcription factor controls the auxin response during carrot root expansion.

Discovery of Genes Related to IAA Regulation during Root Expansion
Screening of the differentially expressed genes (DEGs) with false discovery rate (FDR) < 0.01 and fold change (FC) ≥ 1.5 found a total of eight IAA synthesis genes, seven catabolism genes, and twelve genes related to auxin signal transduction that were differentially expressed during root expansion (Table 4).In order to obtain the key genes for auxin regulation of the root expansion, all of the above DEGs were correlated with the IAA content and a total of 11 DEGs with correlation coefficients ≥ 0.8 were identified (Figure 5).Their expression was significantly correlated with the changes in IAA content.In addition, by analyzing the expression patterns of the genes in the IAA metabolic pathways, we found three synthetic genes, YUCCA10 (DCAR_012429), TAR2 (DCAR_026162), and AMI1 (DCAR_003244); one catabolic gene AACT1 (DCAR_010070); and five signaling genes IAA22 (DCAR_012516), IAA13 (DCAR_012591), IAA27 (DCAR_023070), IAA14 (DCAR_027269), and IAA7 (DCAR_030713) were differentially expressed in H41 vs. H48, and one catabolic gene LPD1 (DCAR _023341) was differentially expressed in H55 vs. H62 (Table 5), and these genes were consistent with the expression of genes involved in IAA metabolism in the above analysis.DCAR_012429, DCAR_026162, and DCAR_003244 encode the YUCCA, TAA1, and AMI1 enzymes, respectively, which are involved in the tryptamine, IPyA, and IAM pathways of IAA synthesis.Expression of these genes was down-regulated in H41 vs. H48, indicating that the enzyme activity decreased, and the IAA content gradually decreased on 48 d after sowing.These results are consistent with the observed changes in the IAA content, indicating that the down-regulated gene expression led to a decrease in the IAA content and the root expansion rate consequently decreased gradually.Expression levels of DCAR_012429 and DCAR_003244 were up-regulated in H34 vs. H41, indicating that the rate of IAA synthesis increased and the IAA content also increased from 34 d to 41 d after sowing.At this stage, the carrot roots were in a phase of rapid growth.DCAR_010070 and DCAR_023341 encode the ACAT and OADH enzymes, respectively, and they were up-regulated in H41 vs. H48 and in H55 vs. H62.Acetyl-CoA was continuously produced to participate in glucose metabolism and maintain the dynamic IAA balance in the plants.The DEGs involved in auxin signal transduction were all genes encoding the AUX/IAA protein and were down-regulated in H41 vs. H48.The activities of the AUX/IAA-encoding genes were affected at once when the IAA content changed, indicating that the AUX/IAA transcription factor family regulated the synthesis of IAA from 41 d to 48 d after sowing and led to a drastic decrease in the IAA content (Table 5 and Figure 6).Therefore, it is suggested that the above 11 genes are candidate genes that play a crucial role in IAA regulation of carrot root expansion, but their specific mechanisms of action need to be further explored.

Discussion
During the root formation of carrot, changes in endogenous hormone content a main regulatory factors.The process of maintaining continuous root growth invo complex network of hormonal regulation, in which auxin plays an important role.concentrations of auxin in the roots have the effect of localizing and establishing roo cell populations, which can promote cell differentiation, whereas the process of t peated division of daughter cells requires appropriately reduced auxin concentr [26].In the present study, the IAA content at the early growth stages of root develop was much higher than at the later growth stages, which confirms that auxin plays a portant role in root initiation and formation processes.However, Perrin et al. showe carotenoids are differentially accumulated in phloem and xylem carrot root depend the genotype [27].Therefore, whether there are differences in auxin content betwee ferent tissues in carrot root and its specific mechanism will be a valuable research tion.In addition, Khadr et al. showed that increased IAA content may cause an in in the number of xylem cells [28]; however, the specific mechanism needs to be f studied.
To further clarify the role of IAA in the root expansion of carrots, we analyze associated metabolic processes.No genes related to the IAOx pathway were detected gesting that this pathway may not exist in carrot root growth.Studies have show

Discussion
During the root formation of carrot, changes in endogenous hormone content are the main regulatory factors.The process of maintaining continuous root growth involves a complex network of hormonal regulation, in which auxin plays an important role.High concentrations of auxin in the roots have the effect of localizing and establishing root stem cell populations, which can promote cell differentiation, whereas the process of the repeated division of daughter cells requires appropriately reduced auxin concentrations [26].In the present study, the IAA content at the early growth stages of root development was much higher than at the later growth stages, which confirms that auxin plays an important role in root initiation and formation processes.However, Perrin et al. showed that carotenoids are differentially accumulated in phloem and xylem carrot root depending on the genotype [27].Therefore, whether there are differences in auxin content between different tissues in carrot root and its specific mechanism will be a valuable research direction.In addition, Khadr et al. showed that increased IAA content may cause an increase in the number of xylem cells [28]; however, the specific mechanism needs to be further studied.
To further clarify the role of IAA in the root expansion of carrots, we analyzed the associated metabolic processes.No genes related to the IAOx pathway were detected, suggesting that this pathway may not exist in carrot root growth.Studies have shown that IAOx plays a role in the synthesis of IAA in Arabidopsis thaliana and this pathway is specific to Brassica crops [29,30].
In the tryptamine and the IPyA pathways, TAR2 (DCAR_026162) encodes L-tryptophanpyruvate transaminase, which mainly acts on the elongation of the hypocotyls, lateral root formation, gravity growth [31], and conversion of tryptophan into indolepyruvate in plants.Subsequent research has confirmed that TAA and members of the YUCCA gene family are the first pathways jointly involved in IAA synthesis and the most important pathway in the IAA-dependent synthesis of tryptophan, controlling changes in IAA levels [12].In the present study, the expression of the YUCCA gene DCAR_012429 41 days after sowing was significantly higher than at 48 days and regulated the rapid growth of the carrot roots, and it is a rate-limiting enzyme in the IPyA pathway.Therefore, we speculate that YUCCA10 (DCAR_012429) is a candidate gene for IAA metabolism in the process of root expansion.
ALDH-encoding genes ALDH2B4 (DCAR_011215 and DCAR_031309) and ALDH3H1 (DCAR_009596) are specifically expressed in carrot roots, and they mainly respond to abiotic stress.They are considered plant antidotes [32].AMI1, a key enzyme in the IAM pathway, hydrolyzes indole-3-acetamide to produce IAA, and this specific amidase is widespread among plants.Genes encoding this enzyme are actively expressed in tissues with a high IAA content, which, in plants, are mostly vigorously growing tissues [33,34], consistent with our results.DCAR_003708 and DCAR_003244 were highly expressed alongside elevated IAA levels at the early growth and developmental stages in the carrots.In IAA degradation, CAT, DAO, OADH, and ACAT are key enzymes that metabolize excess IAA and maintain the IAA balance in plants, thereby promoting root expansion [35].
IAA is synthesized in plants and specifically binds to the auxin receptor SCF TIR1 in the nucleus, thus starting the auxin signal transduction pathway.Recently, a study showed that TIR1/AFB-Aux/IAA pathway rapidly modulates root cell elongation [36].In our experiment, only the AUX/IAA-encoding genes are differentially expressed, as they are more sensitive to changes in IAA content.When the expression of the gene YUCCA10 (DCAR_012429) was up-regulated in H34 vs. H41, resulting in an increase in IAA content, signal transduction was actively expressed in the plant and the corresponding genes were up-regulated, thereby promoting root expansion.

Plant Materials and Cultivation
The 'Dongfanghong' carrot variety was selected as the experimental material and coated carrot seeds were sown in a plastic greenhouse at the experimental station of the College of Horticulture, Shanxi Agricultural University, China, on 7 August 2020.Ridge cultivation was employed with a width of 60 cm and a height of 20 cm.Two rows were sown on each ridge and the carrot plants were traditionally managed.

Phenotyping and Sampling of the Carrot Roots at Different Growth Stages 4.2.1. Determination of Morphological Indices in the Carrot Roots
To understand the dynamics of the carrots during enlargement, the root length and diameter were measured weekly from 13 to 62 days after sowing.

Determination of the IAA Content at Different Growth Stages
Approximately 2 g of each sample was taken in triplicate and IAA content was determined using high-performance liquid chromatography (HPLC), as described by Zhong et al. [37], using a Hypersil ODS C18 column (250 mm × 4.6 mm, 5 µm; Thermo Scientific, Waltham, MA, USA) and an Agilent 1200 UPLC system (Agilent Technologies Inc., Santa Clara, CA, USA).Isocratic elution was performed using a mobile phase of 20:80 methanol/acetic acid aqueous solution (pH 3.6) at a flow rate of 1 mL•min −1 with an injection volume of 20 µL.The detection wavelength was 254 nm.Approximately 1 g of each sample was taken in triplicate, an RNAprep Pure Plant Kit (DP432, TIANGEN, Beijing, China) was used to extract the total RNA, and the purity of the RNA was tested.Qualified RNA samples were added to magnetic beads with Oligo (dT) to enrich and fragment the mRNAs, which served as templates to synthesize and purify cDNAs.The cDNAs were then end-repaired, poly A tails added, and sequencing connectors attached.Appropriately sized fragments were selected for the construction of a cDNA library using PCR amplification.After quality testing (effective concentration of the library > 2 nM), the Illumina platform was used for sequencing.

Analysis of IAA-Metabolism-Related Genes
From the RNA-Seq results, the expression of genes encoding key enzymes in the IAA synthesis and degradation pathways (tryptophan metabolism, ko00380) and auxin signal transduction pathway (plant hormone signal transduction, ko04075) in different periods was analyzed to explore the molecular basis of the IAA content changes during carrot root expansion.GO, KEGG, Ensemble Plants, Panther, and Dicots PLAZA database were used to annotate IAA metabolism and signal-transduction-related genes (Table S2).

Screening of IAA-Metabolism-Related Genes during Root Expansion
The screening threshold for genes involved in IAA metabolism was set as FDR < 0.01 and FC ≥ 1.5, and the DEseq package was used to identify DEGs [42].Expression levels of the DEGs and IAA levels were then correlation analyzed to obtain candidate genes related to IAA regulation in root expansion.

Correlation Analysis between IAA Content and IAA Metabolism and Signal-Transduction-Related DEGs
In order to obtain the key genes for auxin regulation during root expansion, with p ≤ 0.05 as the significant threshold, the correlation index was calculated by Pearson's correlation method.The Gidio Bioinformatics Cloud Tool Platform (http://www.omicsshare.com/tools (accessed on 14 January 2024)) was used to perform the correlation analysis between IAA content and the FPKM values of DEGs related to IAA metabolism and signal transduction (Table 4).
cDNAs were synthesized using reverse transcription of the qualified RNA samples using a PrimeScriptTM RT reagent kit (Perfect Real Time, RR037A, TaKaRa, Beijing, China) and diluted to 80 ng•µL −1 .An ABI 7500 fluorescence quantitative PCR instrument was used with TB Green ® Premix Ex Taq™ II (Tli RNaseH Plus, RR820A, Takara), and three biological replicates were included.Thermal cycling was performed at 94 • C for 5 min, followed by 40 cycles at 94 • C for 30 s, 56 • C for 30 s, and 72 • C for 30 s. Relative gene expression levels were calculated using the 2 −∆∆CT method [43].

Statistical Analysis
Statistical analyses were performed using Microsoft Excel and SPSS version 28.0.1.1.Differences were analyzed by an ANOVA test, followed by the Duncan's new multiple range test for root length and diameter of carrot and IAA content.Clustering of gene expression was performed using R version 3.6.1 pheatmap package with Ward method on Euclidian distances.

Conclusions
Dynamic analysis of the root growth in carrots showed that the rapid expansion period of roots was from 34 to 41 days after sowing.Determination of the IAA content at five different growth stages during root expansion revealed that the IAA content was the highest on 41 d, but the growth of the roots subsequently slowed down and the IAA levels also decreased.Analysis of the expression of genes encoding the key enzymes in IAA metabolism identified 21 IAA-synthesis-related genes, 10 IAA-degradation-related genes, and 22 genes related to auxin signal transduction that were consistent with the IAA content.Among them, 12 IAA synthesis DEGs, 8 IAA degradation DEGs, and 13 DEGs related to auxin signal transduction were obtained.Analysis of the correlations between the expression of these genes and the IAA levels showed that the following genes were most closely related to root expansion: three IAA synthesis genes, YUCCA10 (DCAR_012429), TAR2 (DCAR_026162), and AMI1 (DCAR_003244); two degradation genes, LPD1 (DCAR_023341) and AACT1 (DCAR_010070); and five genes related to auxin signal transduction, IAA22 (DCAR_012516), IAA13 (DCAR_012591), IAA27 (DCAR_023070), IAA14 (DCAR_027269), and IAA7 (DCAR_030713).These results provide a reference for a further exploration of the mechanism of root expansion in future work.

Figure 1 .
Figure 1.Root phenotype of carrot at different growth stages.(A) Appearance; (B) cross d, 41 d, 48 d, 55 d, and 62 d after sowing; (C) root length (cm) and growth rate; (D) ro (mm) and growth rate.Different letters indicate significant difference.Duncan's test, p <

2. 3 . 3 .
Real-Time Quantitative PCR Validation of the RNA-Seq Data Six genes were randomly selected for real-time quantitative PCR (RT-qPCR the reliability of the gene expression results obtained with RNA-Seq.As shown S3, the fragments per kilobase of exon model per million mapped fragments (F

Figure 3 .
Figure 3. IAA metabolism process and key enzymes.The pink, blue, green, and yellow represent the tryptamine pathway, IPyA pathway, IAOx pathway, and IAM pathway, respectively.Red letters: key enzymes.Dotted lines: valid pathways.

Figure 3 .
Figure 3. IAA metabolism process and key enzymes.The pink, blue, green, and yellow represent the tryptamine pathway, IPyA pathway, IAOx pathway, and IAM pathway, respectively.Red letters: key enzymes.Dotted lines: valid pathways.

Figure 4 .
Figure 4. Auxin signal transduction pattern diagram.Auxin-inducible genes have AuxRE in th promoters, which are bound by dimers of the ARF protein family.Gene expression is prevented recruitment of Aux/IAA transcriptional repressors to these promoters via their interaction with ARFs.Solid arrows: verified pathway; dashed arrows: predicted pathways.DEGs encoding TI ARF, and AUX/IAA are listed.

Figure 4 .
Figure 4. Auxin signal transduction pattern diagram.Auxin-inducible genes have AuxRE in their promoters, which are bound by dimers of the ARF protein family.Gene expression is prevented by recruitment of Aux/IAA transcriptional repressors to these promoters via their interaction with the ARFs.Solid arrows: verified pathway; dashed arrows: predicted pathways.DEGs encoding TIR1, ARF, and AUX/IAA are listed.

Figure 5 .
Figure 5. Correlation heatmap of IAA content and FPKM values of all DEGs related to IAA metabolism and signal transduction.Values are Pearson correlation coefficients; |r| ≥ 0.8 indicates gene expression significantly correlated with IAA content; yellow indicates genes significantly correlated with changes in IAA content.

Figure 5 .
Figure 5. Correlation heatmap of IAA content and FPKM values of all DEGs related to IAA metabolism and signal transduction.Values are Pearson correlation coefficients; |r| ≥ 0.8 indicates gene expression significantly correlated with IAA content; yellow indicates genes significantly correlated with changes in IAA content.

Figure 6 .
Figure 6.Expression analysis of DEGs during auxin metabolism.(A) IAA synthesis and degra pathway; (B) auxin signal transduction pathway.The yellow highlights represent the genes expression levels are closely related to the change in IAA content, namely candidate genes.

Figure 6 .
Figure 6.Expression analysis of DEGs during auxin metabolism.(A) IAA synthesis and degradation pathway; (B) auxin signal transduction pathway.The yellow highlights represent the genes whose expression levels are closely related to the change in IAA content, namely candidate genes.

Table 1 .
Expression of key enzyme genes in IAA synthesis process.

Table 2 .
Expression of key enzyme genes in IAA degradation process.

Table 3 .
Expression of genes related to auxin signal transduction.

Table 4 .
DEGs of tryptophan metabolism in root expansion at different stages.

Table 5 .
DEGs of root expansion at different stages.