Skip to main content

Transcriptome profiling reveals the genes involved in tuberous root expansion in Pueraria (Pueraria montana var. thomsonii)

Abstract

Background

Pueraria is a dry root commonly used in Traditional Chinese Medicine or as food and fodder, and tuberous root expansion is an important agronomic characteristic that influences its yield. However, no specific genes regulating tuberous root expansion in Pueraria have been identified. Therefore, we aimed to explore the expansion mechanism of Pueraria at six developmental stages (P1–P6), by profiling the tuberous roots of an annual local variety “Gange No.1” harvested at 105, 135, 165, 195, 225, and 255 days after transplanting.

Results

Observations of the tuberous root phenotype and cell microstructural morphology revealed that the P3 stage was a critical boundary point in the expansion process, which was preceded by a thickening diameter and yield gain rapidly of the tuberous roots, and followed by longitudinal elongation at both ends. A total of 17,441 differentially expressed genes (DEGs) were identified by comparing the P1 stage (unexpanded) against the P2–P6 stages (expanded) using transcriptome sequencing; 386 differential genes were shared across the six developmental stages. KEGG pathway enrichment analysis showed that the DEGs shared by P1 and P2–P6 stages were mainly involved in pathways related to the “cell wall and cell cycle”, “plant hormone signal transduction”, “sucrose and starch metabolism”, and “transcription factor (TF)”. The finding is consistent with the physiological data collected on changes in sugar, starch, and hormone contents. In addition, TFs including bHLHs, AP2s, ERFs, MYBs, WRKYs, and bZIPs were involved in cell differentiation, division, and expansion, which may relate to tuberous root expansion. The combination of KEGG and trend analyses revealed six essential candidate genes involved in tuberous root expansion; of them, CDC48, ARF, and EXP genes were significantly upregulated during tuberous root expansion while INV, EXT, and XTH genes were significantly downregulated.

Conclusion

Our findings provide new insights into the complex mechanisms of tuberous root expansion in Pueraria and candidate target genes, which can aid in increasing Pueraria yield.

Peer Review reports

Introduction

Pueraria is the dry root of leguminous Pueraria montana var. lobata (Willdenow.) Maesen & S. M. Almeida ex Sanjappa et Predeep or Pueraria montana var. thomsonii (Bentham) M. R. Almeida; the former is well known as Gegen (Chinese name) and the latter known as Fengen (Chinese name) in Traditional Chinese Medicine (TCM) [1]. Pueraria, one of the earliest medicinal herbs used in China [2], is native to Southeast Asia and has been employed as a traditional Chinese medicine, food source, and fodder for thousands of years. Pueraria has been used to treat fever, diarrhea, emesis, cardiac dysfunction, liver injury, weight loss, and toxicosis in humans [3]. In the Pharmacopoeia of the People’s Republic of China (PPRC), there are two primary varieties of Pueraria: Pueraria montana var. lobata (Willdenow.) (Pueraria montana var. lobatae Radix, PLR), and Pueraria montana var. thomsonii (Pueraria montana var. thomsonii Radix, PTR). Historically, the herbal industry and many traditional Chinese medicine dispensers/practitioners believed that these species had similar effects. Therefore, they have been treating these herbs as functionally the same, making no distinction in species used for manufacturing and clinical practice. More recently, PTR has been increasingly used in cooking, emerging as a specialty vegetable that is gradually becoming recognized by the market [4]. Thus, we suggest that the food value of PTR may be greater than its medicinal use. For use as an authentic herb, Pueraria has been thoroughly studied for its chemical composition, pharmacological mechanisms, and clinical applications; however, little research has been conducted to explore its food value, especially at the molecular level.

Since the late 1950s, molecular pharmacognosy of Pueraria, including molecular identification, transcriptome sequencing, cloning, and synthesis of functional genes, has been gradually reported. Research has mainly focused on isoflavone biosynthesis [5,6,7], such as puerarin. Recently, using PacBio-RA-II and Ilumina platforms for simultaneous sequencing analysis of Pueraria [8], 3’- methyltransferase, isoflavone-specific 4’- O methyltransferase, and other genes were more highly expressed in PTR than in PLR. Further, using high-throughput sequencing technology, Nithiwat et al. identified 21 genes that may be involved in isoflavone biosynthesis from the combined tissues of young leaves, mature leaves, tuber cortex, and peel tuber of Thailand Pueraria (Pueraria candollei var. Mirifica) [9]. Additionally, Wang et al. performed molecular cloning and functional characterization of PIUGT43 (a novel glucosyltransferase from PLR). Biochemical analysis showed that PIUGT43 promotes C-glycosylation, which converts soybean sapogenins to puerarin, and that PIUGT43 is active against isoflavonoid soybean sapogenins and genistein. Therefore, the lack of activity against other potential receptors, including flavonoids, confirmed the role of PIUGT43 in puerarin biosynthesis [10].

In addition to isoflavones, the physicochemical properties of Pueraria starch, such as viscosity, thermal instability, retrogradation, gelatinization, and digestibility, comprise current topics of interest in research of Pueraria’s food value [11, 12]. Tuberous root development in Pueraria is a complex process regulated by multiple genes, whose expression is determined by light, temperature, hormones, carbon metabolites, cell expansion and division, and TFs [13, 14]. The ADP-glucose pyrophosphorylase (AGPase) gene, for example, a key enzyme in the starch synthesis of Pueraria, affects the composition of starch by regulating changes in its activity, root length, root diameter, weight per plant, and yield [15]. Recent research has focused only on the molecular mechanism of tuberous root expansion during Pueraria growth and development. For example, Liu et al. screened five specific sequences for controlling tuberous root enlargement in PTR by cDNA-AFLP, which are generally involved in signal transduction, stress resistance response, electron transport, and photosynthetic metabolism [16]. To date, few genes related to tuberous root formation have been identified, and no specific genes regulating tuberous root expansion in Pueraria have been identified. Therefore, more research is needed to reveal the molecular mechanism of Pueraria tuberous root expansion.

In production, tuberous root expansion is a determinant of the yield of PTR cultivated as a medicinal and edible vegetable. However, tuberous root expansion is a complex biological process, and its mechanism is still unclear. It was found that some specific genes and proteins associated with starch and phytohormone synthesis, as well as various TFs, are involved in tuberous root formation and expansion [17,18,19], but there are many genes should be found at the transcriptional level. In this study, we thoroughly investigated the transcriptomic mechanism of tuberous root expansion in PTR. Since we did not find any reports of Pueraria tuberous root expansion, we first observed tuberous root changes during the developmental process of “Gange No. 1”, an annual local PTR variety in Jiangxi province. After then, we referred to the sampling time of other metamorphic organ expansion crops, such as sweet potato [20], potato [21], and yam [22], and divided the tuberous root expansion process of PTR into six stages. Next, we compared the difference in morphological indices, sugar, starch and flavonoid accumulation, endogenous hormone content, and transcriptomes of six different developmental phases. We aimed to explore the gene regulatory networks associated with Pueraria tuberous root expansion and to mine the key genes that mediate tuberous root expansion. The findings could directly guide PTR production, improve the technical level of PTR cultivation, and help to address practical problems in the production process of PTR, such as increasing production or precocious cultivation. Additionally, this study laid the foundation for further research on the puffed mechanism at the molecular gene level, potentially guiding the breeding of PTR to maximize its food and medicinal values.

Results

The phenotype of tuberous root expansion

Samples were collected from field-grown cultivar “Gange No. 1” during its unexpanded and expanded stages (Fig. 1A). From Fig. 1, we did not find any expansion changes in the tuberous root at the P1 stage and defined it as an unexpanded stage (Fig. 1B). Tuberous root expansion initiated from the P2 stage, and went through five expansion periods in total: initial-expansion (the P2 stage) (Fig. 1C), mid-expansion (the P3 stage) (Fig. 1D), mid-late expansion (the P4 stage) (Fig. 1E), late-expansion (the P5 stage) (Fig. 1F), and end-expansion (the P6 stage) (Fig. 1G). We found significant phenotypic differences between the P1 stage (unexpanded stage) and the P2–P6 stages (expanded stages) in multiple tuberous root traits. First, the transverse diameter of the tuberous roots rapidly increased from P1 to P3. Until the P3 stage, the transverse diameter increased by 7.39 cm compared to that of the P1 stage (p < 0.05). Meanwhile, tuberous root weight also increased rapidly from P1 to P3, and compared to that of the P1 stage, root weight increased by 708.2 g in P3 (p < 0.05) (Table 1).

Fig. 1
figure 1

The phenotype of different developmental stages of PTR tuberous root. (A) Field cultivation of “Gange No. 1” in May; (B) P1 represents the unexpanded stage at 105 d after transplanting; (C) P2 represents the initial-expanded stage at 135 d after transplanting; (D) P3 represents the mid-expanded stage at 165 d after transplanting; (E) P4 represents the mid-late expanded stage at 195 d after transplanting; (F) P5 represents the late-expanded stage at 225 d after transplanting; (G) P6 represents the end-expanded stage at 255 d after transplanting

Table 1 Statistical analysis of root-related traits in different development stages of PTR

The longitudinal diameter of tuberous roots changed rapidly after the P3 stage, increasing by 17.42 cm (p < 0.05) from P4 to P6. Tuberous root weight gain slowed from the P4 stage onwards, with root weight only increasing by 217.86 g at P6 (p < 0.05) (Table 1). These results indicate that the P3 stage might be an essential boundary point at which tuberous roots of PTR grow faster firstly, especially in terms of lateral growth thickening and weight gain. In contrast, the growth was relatively slow in the later stages, mainly in the longitudinal elongation at both ends.

Microstructural changes in tuberous root expansion

Based on the above experimental results, we considered that the changes in tuberous roots at the P1, P3, and P6 stages were representative, and therefore carried out cytological observations on the tuberous roots at these three developmental stages. The phloem (Ph), vascular cambium (VC), xylem (X), vessel (VE), and pith (P) could be distinguished among the three developmental stages, with the most noticeable distinction between xylem and phloem tissues (Fig. 2).

Fig. 2
figure 2

Cytological observation in different development stages. Tuberous roots from P1, P3, and P6. C: cortex; VC: vascular cambium; P: pith; PR: pith ray; Ph: phloem; Phr: phloem ray; X: xylem; SP: secondary phloem; SX: secondary xylem; XPa: xylem parenchyma; SG: starch granules; WF: wood fiber. (A) Cytological observation of P1, P3, and P6, respectively; (B) Number of total cells per mm2; (C) Number of cells with starch granules per mm2; (D) Number of vessels per mm2. Bar, 1000 μm and 200 μm

In the P1 stage, the shape of the vessels was regular and well-developed, and the lumen diameter was significantly larger than that of thin-walled cells in the xylem (Fig. 2A). The average number of cells per unit area (mm2) reached 656 (Fig. 2B), and the number of intracellular amyloplasts was small and sparsely distributed, mostly at the edges of cells (Fig. 2A and C). In the P3 stage, the proportion of the secondary xylem increased along with tuberous root expansion (Fig. 2A). The number of cells per mm2 decreased and the number of vessels was more differentiated than that in the P1 stage (Fig. 2B and D). The thin-walled cells rapidly became larger, the shape from the original sub-circular constantly increasing mutual extrusion into irregular polygonal and the amyloplasts increased within the cells (Fig. 2C). Additionally, from the middle column of the pith, secondary xylem, newly formed xylem, and starch grains showed a distribution pattern of radial-neatly, dense-disorganized sparse. At P6, the number of cells per mm2 still decreased and the vessels were regular and well developed (Fig. 2B and D), with a wide distribution. At this time, starch mainly accumulated in the parenchyma cells of the secondary xylem, though a few starch granules were also distributed in secondary phloem cells (Fig. 2A). Based on histological observations, tuberous root expansion before the P3 stage may be mainly related to the increase in the number and size of xylem thin-walled cells. In contrast, tuberous root expansion after P3 may be mainly related to the stretch in cell volume. Therefore, starch granules accumulate through PTR tuberous root expansion.

Sugar, starch and flavonoid accumulation

With the tuberous root expansion of Pueraria, the dynamic changes in glucose, fructose, sucrose, and starch content showed different degrees of increase (Fig. 3A–D). Among them, fructose and starch content showed steady growth throughout development, while glucose and sucrose contents fluctuated rapidly in the period from P2–P3 and P4–P6, respectively. Soluble sugars increased from P5 to P6, though there was a slight increase in other stages (Fig. 3E); while reducing sugar content showed a slow increasing trend, with a significant increase at the period from P5–P6 (p < 0.05) (Fig. 3F). In addition, total flavonoid content showed a significant increase with the expansion of the PTR tuberous root (Fig. 3G); while puerarin content increased in the P1–P4 stages, but showed fluctuations occur in the P5 stages (p < 0.05) (Fig. 3H).

Fig. 3
figure 3

Dynamic changes of sugar, starch, and flavonoid contents in different developmental stages. (A) Glucose content; (B) Fructose content; (C) Sucrose content; (D) Starch content; (E) Soluble sugar content; (F) Reducing sugar content; (G) Total flavonoid content; (H) Puerarin content

Endogenous hormone content change

IAA, ZR, and GA3 all decreased during tuberous root expansion, with ZR decreasing more gradually than the others. In contrast, IAA and GA3 decreased more during the P1–P3 stages and tended to level off (Fig. 4A–C). ABA content showed a “decrease–increase–decrease” trend in the process of tuberous root expansion, characterized by a significant decrease before the P3 stage, followed by a rapid increase during the P4 and P5 stages, reaching a maximum in the P5 stage, and a final decline in the P6 stage (Fig. 4D).

Fig. 4
figure 4

Dynamic changes of endogenous hormone contents in different developmental stages. (A) IAA content; (B) ZR content; (C) GA3 content; (D) ABA content

Transcriptome data acquisition

We generated 18 cDNA libraries from unexpanded tuberous roots (P1) and expanded tuberous roots at different stages (P2, P3, P4, P5, and P6) in PTR. Based on the Illumina sequencing results, we assessed gene expression levels for each replicate using principal component analysis (PCA) (Fig. 5A). In this study, two main components, PC1 and PC2, were extracted at 57.3% and 17.7%, respectively, with a cumulative contribution rate of 75.0%. In the PCA score plot, the repeated samples were compactly gathered, indicating that the samples were reproducible, and the data were ready for subsequent analysis. In addition, the P4, P5, and P6 samples clustered far from the P1 samples, indicating that tuberous root expansion induced changes in gene expression between stage types. All Unigenes under the five different groups were annotated and compared using the COG, GO, KEGG, KOG, Swiss-Prot, and eggnog databases (Table 2). We identified 17,441 differentially expressed genes (DEGs) in the six comparison groups. Genes with a false discovery rate (FDR) < 0.05 and an absolute value of log2 ratio ≥ 1 were selected as DEGs. An additional movie file (Table S1) shows this observation in greater detail. The results of the DEGs expression analysis showed that compared with that of the unexpanded stage, the number of differential genes in an expanded stage was relatively high at 255 days (P6 stage) of tuberous root expansion, and the proportion of downregulated genes was higher than that of upregulated genes in all comparison groups (Fig. 5B). Finally, only 386 DEGs (Table S2) were common to all six comparison groups (Fig. 5C).

Fig. 5
figure 5

Comparative transcriptome analysis of P1 stage with P2, P3, P4, P5, and P6 stages of PTR. (A) Scores scatter plot of root transcriptome in unswollen root and swollen root as determined by PCA. (B) The number of regulated genes between the P1 and P2–P6 stage; the red color represents the upregulated genes and the green color represents the downregulated genes. (C) Venn diagram of regulated genes among different developmental stages

Table 2 Reference tag database and major characteristics of DEG libraries

GO and KEGG enrichment analyses of DEGs

To further determine the main biological functions of all DEGs shared by the unexpanded and expanded stages, functional annotation was performed by mapping all the common DEGs to GO terms in the GO database. GO enrichment analysis was implemented using a Bonferroni-corrected p ≤ 0.05. Among the annotated 33,348 Unigenes selected to predict functions by GO annotation, we identified 4,002, 2,623, 6,809, 3,261, and 8,876 DEGs in P1 vs. P2, P1 vs. P3, P1 vs. P4, P1 vs. P5, and P1 vs. P6, assigned to at least one GO term classified into 55, 50, 64, 53, and 69 groups, including 23, 20, 25, 21, and 27 biological processes, 16, 14, 18, 16, and 18 molecular functions, and 16, 16, 21, 16, and 24 cellular components, respectively (Table S3).

Among the DEGs between P1 vs. P2, the “oxidation-reduction process” and “plant-type cell wall organization” were the major terms for biological process, the “cytosolic large ribosomal subunit” and “extracellular region” were the major terms for cellular components, and the “heme binding” was the most represented molecular function term. Between P1 vs. P3 group, the “oxidation-reduction process” was the major term for biological process, the “cytosolic large ribosomal subunit” and “cytosolic small ribosomal subunit” were the major terms for cellular component, and the “structural constituent of the ribosome” was the most represented molecular function term. Between P1 vs. P4 group, the “oxidation-reduction process” and “flavonoid biosynthetic process” were the major terms for biological process, “extracellular region” and “cytosolic large ribosomal subunit” were the major terms for cellular component, and “structural constituent of the ribosome” and “heme binding” were the most represented molecular function terms. Between P1 vs. P5 group, the “oxidation-reduction process” was the major term for biological process, “cytosolic large ribosomal subunit” was the major term for cellular component, and “structural constituent of the ribosome” and “heme binding” were the most represented molecular function terms. Between P1 vs. P6 group, the “plant-type cell wall organization” and “oxidation-reduction process” were the major terms for biological process, “cell wall” was the major term for cellular components, and “structural constituent of the ribosome” was the most represented molecular function term. Notably, the enrichment of DEGs in the P3 stage showed a significant increase in the starch biosynthetic and glycogen metabolic/biosynthetic processes. In contrast, at the P4 stage, the enrichment of DEGs was significantly increased in the flavonoid biosynthetic process.

We performed a pathway enrichment analysis using the KEGG database to further determine the metabolic or signal transduction pathways in which common DEGs may participate in tuberous root expansion. A total of 1,130 (P1 vs. P2), 771 (P1 vs. P3), 1,904 (P1 vs. P4), 860 (P1 vs. P5), and 2,513 (P1 vs. P6) DEGs were assigned to 115, 116, 126, 118, and 125 pathways, respectively, by KEGG pathway enrichment (q ≤ 0.05).Additionally, 11, 10, 11, and 14 pathways were identified as significantly enriched in P1 vs. P2, P1 vs. P3, P1 vs. P4, P1 vs. P5, and P1 vs. P6, respectively (Table 3). Among them, the “biosynthesis of amino acids (ko01230)”, “carbon metabolism (ko01200)”, and “starch and sucrose metabolism(ko00500)” were the top three represented path ways among the DEGs between P1 and P2 (Fig. 6A). “Starch and sucrose metabolism (ko00500)”, “carbon metabolism (ko01200)”, and “plant hormone signal transduction (ko04075)” were the top three path ways among DEGs between P1 and P3 (Fig. 6B). The “carbon metabolism (ko01200)”, “biosynthesis of amino acids (ko01230)”, and “starch and sucrose metabolism (ko00500)” were the top three represented pathways among the DEGs between P1 and P4 (Fig. 6C). The “starch and sucrose metabolism (ko00500)”, “amino sugar and nucleotide sugar metabolism (ko00520)”, and “carbon metabolism (ko01200)” were the top three represented pathways among the DEGs between P1 and P5 (Fig. 6D). The “biosynthesis of amino acids (ko01230)”, “carbon metabolism(ko01200)”, and “starch and sucrose metabolism (ko00500)” were the top three represented path ways among the DEGs between P1 and P6 (Fig. 6E).

Table 3 KEGG enrichment analysis of common differential genes
Fig. 6
figure 6

The top 20 KEGG enrichment pathways between the P1 stage and other expansion stages were obtained by the KEGG database (http://www.kegg.jp/kegg/kegg1.html) analysis. (A–E) P1 vs. P2, P1 vs. P3, P1 vs. P4, P1 vs. P5, and P1 vs. P6. They-axis is the pathway, and the x-axis is the percentage of this pathway of the total Rich Factor. The color depth represents the q-value. The darker the color, the smaller the q-value and the higher the enrichment degree. The size of the dots indicated the number of DEGs in this pathway

Combining the above results, we suggest that the genes involved in the regulation of cell wall and cell division, plant hormone signaling, starch, and sucrose metabolism pathways, and transcription factors might play vital roles in the tuberous root expansion of PRT. Furthermore, these DEGs were annotated using NR, GO, COG, and KEGG annotations, and we found many DEGs involved in these four pathways, indicating that these three pathways played an important role in the process of Pueraria tuberous root expansion. Therefore, we analyzed the genes related to these pathways.

Cell wall and cell cycle

The transcriptome data obtained in this study revealed that 54, 22, 57, 45, and 84 genes were related to the cell wall and cell cycle from the DEGs shared by P1 vs. P2, P1 vs. P3, P1 vs. P4, P1 vs. P5, and P1 vs. P6, respectively (Table S4). Among them, LOB domain-containing protein 4—HL_transcript_188442 (LBD4), and WUSCHEL-related homeobox 4—HL_transcript_36056 (WOX4) were significantly up-regulated at the P3 stage. Moreover, the genes involved in cell division, such as cell division protease (FtsH) in the P2 stage, cell division cycle 2 (CDC2) family in the P4, P5, and P6 stages, and cell division cycle 48 (CDC48) in the cyclin-dependent kinase (CDK) in the P4, P5, and P6 stages, were significantly upregulated in the tuberous root expansion stages. The genes involved in cell extension and expansion, including xyloglucan endotransglucosylase/hydrolase (XTH), expansin (EXP), and extension (EXT), were also significantly changed in the tuberous root expansion stage. These results indicate that the formation and development of tuberous roots are dependent on active meristems and cell division.

In total, 43 genes were identified as DEGs shared between the six stages. Among them, 5 genes were upregulated and 38 genes were downregulated during tuberous root development (Fig. 7A). CDC48, XTH, EXT, EXP, and GDSL/SGNH were the major components of this pathway. HL_transcript_10995 (CDC48) was significantly upregulated, and its expression level increased successively in P2–P6 compared to that of the P1 stage of tuberous root expansion. HL_transcript_18370 and HL_transcript_30528 (XTH), HL_transcript_72441 (EXP), HL_transcript_16215 and HL_transcript_54537 (EXT), and HL_transcript_30674, HL_transcript_54454, and HL_transcript_54561 (GDSL/SGNH) were significantly downregulated, and the expression levels decreased successively during this process.

Fig. 7
figure 7

Heatmap of highly expressed DEGs. Every row shows a related gene of tuberous root expansion. Red, yellow, and blue indicate high, middle, and slow levels of mRNA expression, respectively. (A) Expression of cell wall and cell cycle; (B) Expression of hormone signal; (C) Expression of starch and sucrose metabolism; (D) Expression of TFs

Based on the importance of the P3 stage in the expansion process, we mined 10 genes that were related to the cell wall and cell cycle from the DEGs shared by P3 vs. P1, P3 vs. P2, P3 vs. P4, P3 vs. P5, and P3 vs. P6 (Table 4). In the table, Log2FC is the Fold change value, which indicates the ploidy relationship of gene expression. Finally, we screened a total of 2 DEGs related EXT2-like (HL_transcript_39504) and pectin acetyl esterase 12-like (HL_transcript_40999), after raising the Log2FC value screening threshold from an absolute value greater than 0.5 to 1.

Table 4 DEGs enrichment analysis between P3 and other expansion stages

Plant hormone signal

We screened for key DEGs in auxin, cytokinin, gibberellin, and abscisic acid signaling during the unexpanded stage (P1) and expanded stages (P2–P6) of Pueraria tuberous roots (Table S5), resulting in 495 DEGs with a significantly different expression, including 274 upregulated genes and 221 downregulated genes. The auxin signal transduction pathway was the most active, followed by the abscisic acid signal transduction pathway. Most genes related to auxin signaling (AUX/IAA and ARF) and cytokinin signaling (CRE1) were significantly upregulated. In contrast, genes associated with auxin signaling (GH3), cytokinin signaling (ARR-A and ARR-B), and abscisic acid signaling (PP2C) were significantly downregulated. In this study, the plant hormone signal transduction pathway was one of the most enriched KEGG pathways in the P1 vs. P3 group, indicating that hormone signaling plays a vital role in the P3 stage during tuberous root expansion. In this group, 23 auxin-related genes (AUX/IAA, ARF), 5 cytokinin-related genes (CRE1), 7 gibberellin-related genes (DELLA, PYR/PYL), and 8 abscisic acid-related genes (PP2C) were upregulated in the tuberous root expansion stage, implying that they may be related to cell expansion during the secondary growth of cambium.

Furthermore, 20 genes were identified as DEGs shared between the six stages. Among them, 12 genes were upregulated and 8 genes were downregulated during tuberous root development (Fig. 7B). The genes ARF, PYL8, ARR, and SPY were the major components of this pathway. HL_transcript_72762 (PYL8), HL_transcript_60732 (SPY), HL_transcript_5667, HL_transcript_45339 (ARF), HL_transcript_63852 (ARF), and HL_transcript_72623 (ARF) were significantly upregulated, and their expression levels increased successively in P2–P6 compared to those in the P1 stage of tuberous root expansion. In contrast, HL_transcript_60151 (ARR) was significantly downregulated and its expression levels decreased successively during this process. We only screened a DEG–related ABA catabolic key enzyme (HL_transcript_13040) from the DEGs shared in P3 vs. P1, P3 vs. P2, P3 vs. P4, P3 vs. P5, and P3 vs. P6 (Table 4), but the absolute value of Log2FC value does not exceed 1.

Starch and sucrose metabolism

In this study, changes in the expression of sucrose synthase (SuSy) and invertase (INV) were active in the tuberous roots of Pueraria in P1 vs. P4 and P1 vs. P6. In P1 vs. P4, 25 SuSy genes were significantly upregulated during tuberous root expansion, whereas 8 INV genes were significantly downregulated. In P1 vs. P6, 6 SuSy genes were significantly upregulated during tuberous root expansion, whereas 18 INV genes were significantly downregulated. Furthermore,​ all sucrose phosphate synthase (SPS) genes were upregulated during tuberous root expansion before the P3 stage.

In addition, Pueraria tuberous root expansion before the P3 stage was dominated by the upregulation of granule-bound starch synthase (GBSS), starch branching enzyme (SBE), soluble starch synthase (SSS)-related, and isoamylase (ISA) genes, and later by the upregulation of GBSS, SBE, and ISA genes. Among them, 9 starch-related genes (1 GBSS, 2 SSS, 4 SBE, and 2 ISA) were significantly upregulated in P1 vs. P2, 23 starch-related genes (9 GBSS, 4 SSS, 6 SBE, and 4 ISA) were significantly upregulated in P1 vs. P3, 14 starch-related genes (9 GBSS, 3 SBE, and 2 ISA) were significantly upregulated in P1 vs. P4, 10 starch-related genes (3 GBSS and 7 SBE) were significantly upregulated in P1 vs. P5, and 11 starch-related genes (4 GBSS, 10 SBE, and 1 ISA) were significantly upregulated in P1 vs. P6 (Table S6). These results indicate that many functional genes are involved in the expansion of PRT tuberous roots.

In addition, 11 genes were identified as DEGs shared by the six stages. Among them, two genes were upregulated and nine genes were downregulated during tuberous root development (Fig. 7C). HL_transcript_6305 (INV), the major representative of this pathway, was significantly downregulated, and its expression levels decreased successively in P2–P6 compared to those in the P1 stage of tuberous root expansion. We mined 4 genes related to starch and sucrose metabolism from the DEGs shared by P3 vs. P1, P3 vs. P2, P3 vs. P4, P3 vs. P5, and P3 vs. P6 (Table 4). Finally, we screened a total of 2 DEGs related to β-glucosidase40 (HL_transcript_21640 and HL_transcript_40885), after raising the Log2FC value screening threshold from an absolute value greater than 0.5 to 1.

Transcription factors

In this study, 42 TF genes were identified as DEGs shared by six developmental stages. Among them, 15 TFs were up-regulated, and 20 TFs were down-regulated during the tuberous root expansion (Table S7). bHLH, AP2, ERF, MYB, WRKY, and bZIP TFs were the major represented TF families, which play key regulatory functions in cell differentiation, division and expansion (Fig. 7D). Interestingly, the changes in transcript levels of seven of the eight bHLHs showed an “increasing then decreasing” trend, and the expression levels increased in the P2 stages then successively decreased in the P3 ~ P6 stages in PTR. We mined 10 genes related to TFs from the DEGs shared by P3 vs. P1, P3 vs. P2, P3 vs. P4, P3 vs. P5, and P3 vs. P6 (Table 4). Finally, we screened a total of 2 DEGs related to ERF12 (HL_transcript_47402) and MYB86-like (HL_transcript_49439), after raising the Log2FC value screening threshold from an absolute value greater than 0.5 to 1.

Identification of DEGs involved in tuberous root expansion

To verify the accuracy of the RNA-seq sequencing results, we selected six DEGs (HL_transcript_5667, HL_transcript_6305, HL_transcript_10995, HL_transcript_16215, HL_transcript_18370, and HL_transcript_72441) involved in cell wall and cell cycle, plant hormone signal, as well as starch and sucrose metabolism for qRT-PCR analysis (Fig. 8). The results showed that the relative expression change trends of qRT-PCR were like that of RNA-seq, demonstrating that transcriptome comparison data were reliable.

Fig. 8
figure 8

qRT-PCR validation profiles of 6 randomly selected genes. The data were normalized by using GAPDH as an internal reference. (A) The trend of ARF (auxin response factor, HL_transcript_5667) expression in qRT-PCR and transcriptome. (B) The trend of INV (invertase, HL_transcript_6305) expression in qRT-PCR and transcriptome. (C) The trend of CDC48 (cell division cycle protein 48 homolog, HL_transcript_10995) expression in qRT-PCR and transcriptome. (D) The trend of EXT (extension, HL_transcript_16215) expression in qRT-PCR and transcriptome. (E) The trend of XTH (xyloglucan endotransglucosylase/hydrolase, HL_transcript_18370) expression in qRT-PCR and transcriptome. (F) The trend of EXP (expansion, HL_transcript_72441) expression in qRT-PCR and transcriptome. Black bars represent the data of transcriptome and yellow lines represent the data of qRT-PCR. FPKM values were used to represent the relative expression of genes in the transcriptome. Values are means ± SD (n = 3)

Discussion

The formation and development of tuberous roots directly affect the yield of PTR. Therefore, studying the mechanism of tuberous root expansion is of great importance in improving the yield. The PTR variety used in this study, “Gange No. 1”, is often cultivated as an annual crop, and its entire reproductive period is approximately 260 days. According to the investigation of the traits of Pueraria roots in different reproductive periods, we found that the whole expansion process of “Gange No. 1” could roughly be divided into six stages: the unexpanded, initial-expanded, mid-expanded, mid-late expanded, late-expanded and end-expanded stages. Seedlings were transplanted in a field in the middle of March; vines (stems and leaves) and normal roots grew vigorously from April–June (Fig. 1A). The tuberous roots began to expand in early July (Fig. 1B). Histologically, expansion is caused by the activity of two parts of the meristematic tissue: the primary formative layer (botanically known mainly as the formative vascular layer) and the secondary formative layer (also known as the additional formative layer or paraphyly layer). In this study, we considered that the expansion of the PTR tuberous root was driven by the secondary growth of the xylem (Fig. 2). To play a role in starch accumulation, a large number of thin-walled cells in the secondary xylem divide and differentiate, producing a large number of thin-walled tissues containing starch grains, leading to an increase in tuberous root diameter, which is similar to starch accumulation processes in sweet potato [20] and cassava [23].

Tuberous root expansion relies on an increase in cell volume, in addition to an increase in the number of cells by the division activity of the secondary formative layer [24]. The rate of cell division and the number of cells is usually related to the cell division cycle. For example, genes involved in cell division such as cell division proteases (FtsH), cell division control proteins (CDC48), and cell cycle proteins (CYCT1) were significantly upregulated during radish fleshy root expansion [25]. Additionally, CDK genes, which control the cell division cycle, play an essential role in the expansion of tuber mustard [26]. In our study, the P3 stage (165 days) was the crucial dividing point in the expansion process. Before the P3 stage, round vascular cambium forms and the Pueraria tuberous root begins to thicken. After then, we screened a key gene in the transcriptome data, EXT12-like, which was expressed at a higher level during the P3 stage compared to other expansion stages, and we suggest that it may be an important regulator of tuberous root thickening in PTR. Tuberous root hypertrophy is dominated mainly by increased cell numbers from the formative layer division activity and expanded cell volume, along with a rapid increase in tuberous root weight (Table 1), which is similar to potatoes [21], sweet potatoes [20], and cassava [27, 28].

Subsequently, the weight gain ratio of the tuberous root increased slowly. Tuberous root expansion is dominated by cell volume longitudinal elongation, accompanied by starch and other substance accumulation, resulting in the continuous enlargement of tuberous roots. We found that many vital genes (Table S4), such as LBD4 and WOX4, were significantly upregulated in the P3 stage. The genes involved in cell division and expansion pathways, such as FtsH, CDC, CDK, CDC2, and CDC48 were significantly altered during tuberous root expansion of PTR, especially in the later stages with significant upregulation. EXP1 is a cell wall relaxation protein that is an essential regulator of cell wall elongation during plant cell growth [29, 30], and is closely related to cell size. In sweet potatoes, EXP1 expression levels are higher in the fibrous root stage and lower in the early root storage stage. Antisense EXP1 plants formed storage roots earlier than wild-type plants. Lignification of the fibrous root mid-column was significantly reduced, providing them a higher potential to develop storage roots, suggesting that EXP1 plays a suppressive role in the formation of storage roots in sweet potatoes [31]. In this study, we found a similar result that most EXPs had a significant downregulation trend during tuberous root expansion, indicating that downregulation of the expression of these genes could promote tuberous root expansion of PTR.

Changes in the phytohormone content are also involved in the expansion of tuberous roots as a starchy storage organ [13, 32]. These changes are required to alter cell division patterns and stimulate cell enlargement and metabolism to facilitate the storage of starch and proteins. In this study, the IAA, ZR, and GA3 contents of PTR tuberous roots decreased to different degrees during the expansion process, with IAA and GA3 contents decreasing more and ZR decreasing less. All hormone levels decreased to the lowest value at the end-expansion stage (the P6 stage).

Among them, IAA decreased significantly after the beginning of the expansion and remained at a low level at the mid-expansion and maturity stages, indicating that the low concentration of IAA favored tuberous root expansion, which is consistent with the results of IAA content changes in tuber formation in jicama [33]. Meanwhile, the ZR content of PTR was higher during the initiation stage. However, it significantly decreased with tuberous root expansion, indicating that high concentrations of cytokinin facilitate the initiation of primordium differentiation, which is consistent with observations in onion bulbs [34]. In addition, GA3 content continuously decreased during the entire process of PTR tuberous root expansion. A previous study has shown that a decrease in GA content in potato tubers induces tuber formation [35] while planting under high GA concentrations prevents tuber formation [36]. Similar results were also obtained during the development of sweet potato storage roots [37]. These reports supported our speculation that GA might also be an inhibitor of PTR tuberous root expansion. In transcriptome data, the plant hormone signal transduction pathway was one of the most enriched KEGG pathways in the P1 group compared to the P3 group. Auxin plays a vital role in cambium cell proliferation and cell expansion [38], maintaining the meristem state of cambium cells, and increasing the number of xylem elements [39]. Studies conducted on sweet potato, radish, Rehmannia glutinosa (Di Huang), and Callerya speciosa (Niu Dali) reported that the expression of auxin-related genes is significantly upregulated during tuberous root expansion stages [40,41,42]. In our study, 23 auxin-related genes (10 AUX/IAA and13 ARF) were upregulated in the P1 vs. P3 group (Table S6), indicating that they might be related to cell expansion during the secondary growth of the tuberous root. Cytokinin is involved in the proliferation and development of cambium cells. The expression of auxin-related genes reached the highest level in the rapid growth stage of the tuberous root, which was related to the development and formation of the tuberous root/tuber [39, 43,44,45]. In this study, five cytokinin-related genes (HL_transcript_22447, HL_transcript_59990, HL_transcript_25588, HL_transcript_1520, and HL_transcript_43612) were significantly upregulated during tuberous root expansion, suggesting that cytokinin might promote tuberous root expansion by enhancing cambium development.

In addition, ABA, the classical antagonistic hormone of GA, showed the lowest value at the P3 stage during PRT expansion. From the transcriptome data, we found that most of the PP2C and ABF genes were significantly downregulated in the ABA-responsive pathway during tuber formation. In contrast, PYR/PYL genes were significantly upregulated, suggesting that ABA accumulation during PRT expansion may have a negative feedback regulation on PYR/PYL. Additionally, PP2C expression was enhanced by ABF genes, which is consistent with the literature on ginger rootstock expansion [46]. Although it is known that these hormone-related genes play vital roles during the tuberous root expansion stage, further in-depth studies of their systematic network are needed to assess their role in tuberous root expansion.

The fleshy rootstock becomes an organ for storing carbohydrates (sugars, starch, etc.), accompanied by material and energy changes during development. As a result, the ability to synthesize, accumulate, and metabolize carbohydrates and proteins is an essential criterion for determining the expansion ability of metamorphic rootstocks [47]. Starch is the material basis for the expansion of PTR tuberous roots and it usually accounts for 15–46% of fresh matter weight [48]. In our study, starch grains accumulated significantly with the growth of PTR tuberous roots (Fig. 2), and sugar, sucrose, and starch contents gradually increased, along with total flavonoid and puerarin accumulation (Fig. 3). Combined with the transcriptome data, we found that the pathways “carbon metabolism” and “sucrose and starch metabolism” were present throughout the tuberous root expansion of PRT. These results indicate that carbohydrate metabolism play a vital role in the tuberous root expansion of PTR, which is consistent with our phenotypic data. During the formation and expansion of PTR storage organs, the expression of most SPS, SS, GBSS, SSS, SBE, and ISA genes was significantly upregulated during the process of PRT tuberous root expansion. Most INV genes were significantly downregulated during this process. In contrast, the a-amylase and b-amylase genes first decreased and then increased (Table S6), which is similar to patterns found in most tuber species, such as potato [49], sweet potato [50], radish [51], and cassava [52]. Sucrose, the main product of photosynthesis, supports plant growth and development. Sucrose synthesis in plants is mainly regulated by three enzyme genes: SPS, SuSy, and INV. SPS, the primary source of sucrose synthesis activity [53], was upregulated during PRT tuberous root expansion. As a result, we suggest that SPS might play an essential role in PRT tuberous root expansion, consistent with previous studies in radishes wherein upregulation of SPS played a significant role in the thickening stage of radish taproot [54]. In addition, high transcript levels of β-glucosidase40 in the P3 stage could be related to the increased fiber during expansion, and their regulatory mechanism should be studied in the next step.

The sucrose synthase acts as a bidirectional regulator that synthesizes sucrose and breaks it down into glucose and fructose [55]. Combined with the previous sugar content assay data, SuSy genes upregulated in the pre- and mid-expansion stages of PTR might exercise mainly a catabolic function, consistent with the increase in glucose observed during the P2–P3 stages. In contrast, sucrose synthase primarily exerts a synthetic function during the later phase, consistent with the rapid increase in sucrose during the P4–P6 stages. In addition, INV positively regulates fructose and glucose contents in turnip fleshy roots, thereby accelerating the accumulation of biomass in the roots [56, 57]. Our results most closely resembled those of previous studies on sweet potatoes, which found that INV rapidly decreases to an undetectable level during storage root development [58]. Based on these results, we speculate that INV might not be the primary regulator of glucose and fructose content in PTR tuberous root expansion, and might be related to other genetic characteristics of PTR. Further studies are needed for clarification.

TFs play an important role in the regulation of plant growth and development and secondary metabolism. In this study, we identified 15, 27 TFs that were significantly up-and down-regulated in the tuberous root expansion process. bHLH, AP2, ERF, MYB, WRKY, and bZIP were the main types of TFs, which is different from the results in sweet potato [20] and potato [21], probably because of species differences. Among these TFs, ERF12 and MYB86 are two transcription factors that are prominent in the P3 stage, indicating they may be involved in the regulation of the tuberous root expansion in PTR, but further functional identification studies were needed to confirm the functions of these potential genes.

PTR tuberous root expansion is a complex regulatory process affected by many factors. In this study, through transcriptome analysis combined with phenotype and physiological test results, a hypothetical model of the PTR tuberous root expansion regulatory network was proposed (Fig. 9). The cells in the vascular cambium divide and expand continuously to produce secondary xylem and phloem, resulting in the expansion of tuberous roots. Cell proliferation and size are regulated by several signal transduction pathways (hormone signaling) and metabolic processes (cell wall, sucrose, and starch metabolism). Several genes, including auxin-related genes (Aux/IAA, ARF, CRE1, and PYR/PYL), were highly expressed to promote cell differentiation, division, expansion, sucrose, and starch accumulation in the secondary structure. Additionally, FtsH, CDC, CDK, XTH, EXP, and EXT were involved in the extension and expansion of cell division. Finally, SuSy, SPS, SSS, GBSS, and SBE were involved in the hydrolysis of sucrose and starch synthesis.

Fig. 9
figure 9

A hypothetical model of regulatory network related to tuberous root expansion in PTR

Conclusions

Integrated phenotypic, physiological, and transcriptomic analyses performed here revealed that the P3 stage was a vital boundary point in the expansion period of tuberous roots of PTR. There were 17,441 differential genes shared between the unexpanded (P1) and expanded (P2–P6) stages of PRT. KEGG pathway enrichment analysis revealed that cell wall and cell cycle, plant hormone signal transduction, and starch and sucrose metabolism were the top pathways participating in this process. The qRT-PCR analysis suggested that ARF, INV, CDC48, EXT, XTH, and EXP might play a key role in tuberous root expansion in PTR. Based on the above results, we propose a hypothetical model of a genetic regulatory network associated with tuberous roots in PTR (Fig. 9). The tuberous root expansion of PTR was mainly attributed to cell differentiation, division, expansion, and starch accumulation, all of which are regulated and promoted by specific signal transduction pathways and metabolic processes. These findings not only provide novel insights into the molecular regulation mechanism of tuberous root expansion but also support the theoretical basis for genetic improvement of PTR, which may help enhance PTR yield for food and medicinal purposes.

Methods

Plant materials

The local variety “Gange No. 1” belongs to an annual PTR, and the entire reproductive period is approximately 260 days. The entire planting process of “Gange No. 1” was performed at the Jiangxi Agricultural University Experimental Station in Nanchang (28º 09’ N, 116º 35’ E), Jiangxi Province, China. Pueraria seedlings were transplanted on March 18, 2021. Tuberous roots were collected at 10 am each time across 6 developmental stages at 105 (June 29), 135 (July 29), 165 (August 29), 195 (September 29), 225 (October 29), and 255 days (November 29) post-transplant, named P1, P2, P3, P4, P5, and P6 stages, respectively. For each sample time, three independent biological replicates of tuberous roots were collected. Root samples were frozen immediately in liquid nitrogen and stored at − 80℃ until further processing.

Phenotype assessment

We adopted a bamboo shelf with a nylon mesh cultivation mode. We used single-stem pruning for each PTR plant during cultivation and retained one tuberous root after outcropping. In addition, we used high monopoly cultivation, with the following measurements: width of the ridge bottom (1.1 m), the width of the ridge (60 cm), the height of the ridge (50 cm), and the spacing of the plants (30 cm apart) with planting 1400 seedlings per acre.

For each stage, root diameter, length, and weight were measured and recorded using a ruler and an electronic scale. We measured the root weight after cleaning and drying and recorded the root diameter and length data by measuring the transverse diameters and the most significant longitudinal lengths of the roots. Finally, we performed statistical analyses using SPSS software.

Histological analysis of tuberous roots

Since the P3 stage is an important boundary point during PTR expansion, we collected tuberous roots at three typical expansion stages, P1, P3, and P6, as materials for cytological characterization. Next, the sections were fixed, embedded, cut into slices, and stained according to conventional paraffin sectioning methods. Slices were observed and photographed using a microscope (Nikon ECLIPSE 80i, Tokyo, Japan). Voucher specimens of three PTR developmental stages (P1, P3, and P6) were deposited in the Herbarium of Agronomy College, Jiangxi Agricultural University.

Determination of carbohydrate content

The samples were ground and placed in 50-mL centrifuge tubes to test the main carbohydrate content. The glucose content was determined using the glucose oxidase method [59]. The resorcinol method was used to determine the fructose and sucrose contents [60]. The total starch content was determined using the anthrone sulfate method [61]. The total soluble sugars were determined using the anthrone sulfate method [62]. The reducing sugars were determined colorimetrically using 3, 5-dinitro salicylic acid (DNS) [63].

Determination of hormone contents

We modified the methods used for the extraction and purification of endogenous hormones (IAA, ZR, GA3, and ABA) from those described by Yang et al. [64] and Wang et al. [65]. Culture tissue samples (approximately 0.5 g) were ground using an ice-cooled mortar in 5 mL 80% (v/v) methanol extraction medium containing 1 mM butylated hydroxytoluene and 100 mg polyvinylpyrrolidone (PVP) per gram of fresh material. The extract was incubated at 4 °C for 4 h and centrifuged at 1000 × g for 15 min at the same temperature. The supernatant was passed through a Chromosep C18 column (C18 Sep-Park Cartridge, Waters Corp., Milford, MA, USA) that was prewashed with 1 mL of 80% methanol. Hormone fractions were eluted with 5 mL of 100% methanol, 5 mL of ether, and 5 mL of 100% methanol, then dried under N2, and dissolved in 2 mL of PBS (pH 7.5) containing 0.1% Tween 20 and 0.1% gelatin before using ELISA kits from China Agricultural University (Beijing, China) to estimate hormone levels.

RNA extraction and transcriptome sequencing

Total RNA was extracted according to the instructions in the TRIzol reagent kit (Thermo Fisher Scientific, Waltham, MA, USA). The RNA concentration was measured using a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific). RNA integrity was assessed using an RNA Nano 6000 assay kit from the Agilent Bioanalyzer 2100 system (Agilent Technologies, Santa Clara, CA, USA). The RNA samples were stored at − 80ºC until use. First, mRNA was purified from the total RNA using polyToligo-attached magnetic beads and then broken into short fragments. Random hexamers and RNA fragments were used to prime cDNA synthesis. After purification and connection with adapters, a cDNA library was constructed using PCR amplification. The length of the insert sequence was verified using an Agilent Bioanalyzer 2100 system and the library was quantified using an ABI StepOnePlus Real-Time PCR System. Finally, the qualified cDNA library was sequenced using the Illumina HiSeq 2000 system.

Transcriptome analysis

Low-quality reads (more than 20% of bases with quality ≤ 10) and reads with adapters (more than 5% unknown nucleotides) were filtered to generate only clean reads for data collection. De Novo assembly and redundant sequence removal were performed using Trinity and Tgirl, respectively.

The resultant transcripts were searched against the NCBI nonredundant nucleotide (Nt), the NCBI nonredundant protein sequences (Nr), and SwissProt protein databases for functional annotation using the BLAST algorithm with an E-value cutoff of 1e-5. The available categories of these unique sequences were further analyzed using Clusters of Orthologous Groups of Proteins (GO), Gene Ontology (GO), and Kyoto Encyclopedia of Transcripts and Genome (KEGG) databases using the BLAST and Blast2 GO programs. The clean reads were mapped to the reference using Bowtie 2 to estimate the transcript expression profiles. The expression levels were calculated as fragments per kilobase of exon per million fragments (FPKM) using RSEM. Candidate transcripts involved in swollen root formation were selected from previous reports and databases with FPKM values of the transcripts converted to log10 values (FPKM ≥ 5). They were visualized in a heatmap to identify different expression profiles between P1 vs. P2, P1 vs. P3, P1 vs. P4, P1 vs. P5, and P1 vs. P6.

Identification of candidate genes involved in tuberous root expansion

The expression of the six Unigenes identified by RNA-seq was validated using RT-qPCR. Gene-specific primers were designed based on Unigene sequences (Table 5). The 18SrRNA genes were amplified as an endogenous loading control. Finally, the expression of each gene was confirmed using at least three rounds of independent qRT-PCR.

Table 5 The primers of 6 DEGs validated by RT-qPCR analysis

Statistical analysis

The experimental data were analyzed in triplicate. The data were analyzed using SPSS 17.0. The measurement data were expressed as mean ± standard deviation (SD). Duncan’s multiple range test was used for inter-group and intra-group comparison. p < 0.05 denotes statistical significance.

Data Availability

RNA sequencing data have been deposited in the sequence read archive - SRA (https://www.ncbi.nlm.nih.gov/sra/PRJNA923481). Data supporting the findings of this study are available at https://github.com/xiaoxf820/my-document.

Abbreviations

DEGs:

Differentially expressed genes

IAA:

Auxin

CTK:

Cytokinin

ABA:

Abscisic acid

GA:

Gibberellin

LBD:

LOB domain-containing protein

WOX4:

WUSCHEL-related homeobox 4

XTH:

Xyloglucanendotrans-glucosylase/hydrolase

FtsH:

Cell division protease

CDC2:

Cell division cycle2-like protein

CDC20:

Cell division cycle20-like protein

CDC48:

Cell division cycle protein 48 homolog

CDK:

Cyclin-dependent kinase

CDKI:

Cyclin-dependent kinase inhibitor

EXP:

Expansin

EXT:

Extension

AUX1:

Auxin influx carrie

AUX/IAA:

Auxin-responsive protein IAA

GH3:

Auxin responsive GH3 gene family

TIR1:

Transport inhibitor response

SAUR:

Small auxin up RNA

CRE1:

Cytokinin receptor

ARR-B:

Two-component response regulator ARR-B family

ARR-A:

Two-component response regulator ARR-A family

GID2:

F-box protein GID2

PIF3:

Phytochrome-interacting factor 3

PYR/PYL:

Abscisic acid receptor

PP2C:

Protein phosphatase 2C

ABF:

ABA responsive element binding factor

SuSy:

Sucrose synthase

SPS:

Sucrose phosphate synthase

INV:

Invertase

SS:

Starch synthase

GBSS:

Granule-bound starch synthase

SSS:

Soluble starch synthase

SBE:

Starch branching enzyme

AGPase:

ADP-glucose pyrophosphorylase

ISA:

Isoamylase

References

  1. Commission CP. Pharmacopoeia of the People’s Republic of China. The Chemical Industry Publishing House; 2020.

  2. Prasain JK, Peng N, Rajbhandari R, Michaelwyss J. The chinese Pueraria root extract (Pueraria lobata) ameliorates impaired glucose and lipid metabolism in obese mice. Phytomedicine. 2013;20(1):17–23.

    Article  Google Scholar 

  3. Song W, Li Y, Qiao X, Qian Y, Ye M. Chemistry of the chinese herbal medicine Puerariae Radix (Ge-Gen): a review. J Chin Pharm Sci. 2014;23:347–60.

    Article  CAS  Google Scholar 

  4. Wei S, Chen Y, Xu X. Progress on the pharmacological research of puerarin: a review. Chin J Nat Med. 2014;12:407–14.

    CAS  PubMed  Google Scholar 

  5. Wang X, Li S, Li J, Li C, Zhang Y. De novo transcriptome sequencing in Pueraria lobata to identify putative genes involved in isoflavones biosynthesis. Plant Cell Rep. 2015;34(5):733–43.

    Article  CAS  PubMed  Google Scholar 

  6. Mocan A, Carradori S, Locatelli M, Secci D, Cesa S, Mollica A, Riga S, Angeli A, Supuran CT, Celia C, Marizio LD. Bioactive isoflavones from Pueraria lobata root and starch: different extraction techniques and carbonic anhydrase inhibition. Food ChemToxicol. 2018;112:441–7.

    Article  CAS  Google Scholar 

  7. Wang C, Xu N, Cui S. Comparative transcriptome analysis of roots, stems, and leaves of Pueraria lobata (Willd.) Ohwi: identification of genes involved in isoflavonoid biosynthesis. Peer J. 2021;9:e10885.

    Article  PubMed  PubMed Central  Google Scholar 

  8. He M, Yao Y, Li Y, Yang M, Li Y, Wu B, Yu D. Comprehensive transcriptome analysis reveals genes potentially involved in isoflavone biosynthesis in Pueraria thomsonii Benth. PLoS ONE. 2019;14(6):e0217593.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Nithiwat S, Kittitya T, Pinidphon P, Khwanlada K, Julie Z, Siriporn W, Hubert S, Mami Y, Kazuki S, Wanchai D, Sornkanok V, Supaart S. Transcriptome analysis of Pueraria candollei var. Mirifica for gene discovery in the biosyntheses of isoflavones and miroestrol. BMC Plant Biol. 2019;19(1):581.

    Article  Google Scholar 

  10. Wang X, Li C, Zhou C, Li J, Zhang Y. Molecular characterization of the C-glucosylation for puerarin biosynthesis in Pueraria lobata. Plant J. 2017;90(3):535–46.

    Article  CAS  PubMed  Google Scholar 

  11. Reddy CK, Fei L, Xu B. Morphology, crystallinity, pasting, thermal and quality characteristics of starches from adzuki bean (Vigna angularis L.) and edible kudzu (Pueraria thomsonii Benth). Int J Biol Macromol. 2017;105:354–62.

    Article  CAS  PubMed  Google Scholar 

  12. Liu D, Ma L, Zhou Z, Liang Q, Xie Q, Ou K, Liu Y, Su Y. Starch and mineral element accumulation during root tuber expansion period of Pueraria thomsonii Benth. Food Chem. 2021;343:128445.

    Article  CAS  PubMed  Google Scholar 

  13. Zierer W, Rüscher D, Sonnewald U, Sonnewald S. Tuber and tuberous root development. Annu Rev Plant Biol. 2021;72:551–80.

    Article  CAS  PubMed  Google Scholar 

  14. Zhang H, Wang Z, Li X, Gao X, Dai Z, Cui Y, Zhi Y, Liu Q, Zhai H, Gao S, Zhao N, He S. The IbBBX24-IbTOE3-IbPRX17 module enhances abiotic stress tolerance by scavenging reactive oxygen species in sweet potato. New Phytol. 2022;233:1133–52.

    Article  CAS  PubMed  Google Scholar 

  15. Guo L, Zhang J, Xiao D, Huang Y, Zeng D, He L, Wang A. Relationship between photosynthesis characteristics and starch accumulation of Pueraria lobata during different growth stages. J South Agri. 2018;49(4):662–8.

    CAS  Google Scholar 

  16. Liu D, Wu C, Fan S. cDNA-AFLP analysis of gene expression of Kudzu root expansion-related genes analysis. Acta Agricultural Universitatis Jiangxiensis. 2012;34(2):369.

    Google Scholar 

  17. Wang Z, Fang B, Chen J, Zhang X, Luo Z, Huang L, Chen X, Li Y. De novo assembly and characterization of root transcriptome using Illumina paired-end sequencing and development of cSSR markers in sweet potato (Ipomoea batatas). BMC Genomics. 2010;11:726.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Wang Z, Fang B, Chen X, Liao M, Chen J, Zhang X, Huang L, Luo Z, Yao Z, Li Y. Temporal patterns of gene expression associated with tuberous root formation and development in sweet potato (Ipomoea batatas). BMC Plant Biol2015, 15:180.

  19. Dong T, Zhu M, Yu J, Han R, Tang C, Xu T, Liu J, Li Z. RNA-Seq and iTRAQ reveal multiple pathways involved in storage root formation and development in sweet potato (Ipomoea batatas L). BMC Plant Biol. 2019;19(1):136.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Cai ZQ, Cai ZP, Huang JL, Wang AQ, Ntambiyukuri A, Chen BM, Zheng GH, Li HF, Huang YM, Zhan J, Xiao D, He LF. Transcriptomic analysis of tuberous root in two sweet potato varieties reveals the important genes and regulatory pathways in tuberous root development. BMC Genomics. 2022;23:473.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Kloosterman B, De Koeyer D, Griffiths R, Flinn B, Steuernagel B, Scholz U, Sonnewald S, Sonnewald U, Bryan GJ, Prat S. Genes driving potato tuber initiation and growth: identification based on transcriptional changes using the POCI array. Funct Integr Genomics. 2008;8(4):329–40.

    Article  CAS  PubMed  Google Scholar 

  22. Zhou Y, Luo S, Hameed S, Xiao D, Zhan J, Wang A, He L. Integrated mRNA and miRNA transcriptome analysis reveals a regulatory network for tuber expansion in chinese yam (Dioscorea opposita). BMC Genomics. 2020;21(1):117.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Naconsie M, Lertpanyasampatha M, Viboonjun U, Netrphan S, Narangajavana J. Cassava root membrane proteome reveals activities during storage root maturation. J Plant Res. 2015;129:51–65.

    Article  PubMed  Google Scholar 

  24. Polyn S, Willems A, De Veylder L. Cell cycle entry, maintenance, and exit during plant development. Curr Opin Plant Biol. 2015;23:1–7.

    Article  CAS  PubMed  Google Scholar 

  25. Yu R, Wang J, Xu L, Wang Y, Wang R, Zhu X, Sun X, Luo X, Xie Y, EverlyneM. Transcriptome profiling of taproot reveals complex regulatory networks during taproot thickening in radish (Raphanus sativus L). Front Plant Sci. 2016;7:1210.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Shi H, Wang L, Sun L, Dong L, Liu B, Chen L. Cell division and endoreduplication play important roles in stem swelling of tuber mustard (Brassica juncea Coss. Var. Tumida Tsen et Lee). Plant Biol. 2012;14(6):956–63.

    Article  CAS  PubMed  Google Scholar 

  27. Siebers T, Catarino B, Agusti J. Identification and expression analyses of new potential regulators of xylem development and cambium activity in cassava (Manihot esculenta). Planta. 2017;245(3):539–48.

    Article  CAS  PubMed  Google Scholar 

  28. Mehdi R, Lamm CE, Bodampalli Anjanappa R, Mudsam C, Saeed M, Klima J, Kraner ME, Ludewig F, Knoblauch M, Gruissem W, Sonnewald U, Zierer W. Symplasmic phloem unloading and radial post-phloem transport via vascular rays in tuberous roots of Manihot esculenta. J Exp Bot. 2019;70:5559–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Cosgrove DJ, Li LC, Cho HT, Hoffmann-Benning S, Moore RC, Blecker D. The growing world of expansins. Plant Cell Physiol. 2002;43(12):1436–44.

    Article  CAS  PubMed  Google Scholar 

  30. Li Y, Jones L, McQueen-Mason S. Expansins and cell growth. Curr Opin Plant Biol. 2003;6(6):603–10.

    Article  CAS  PubMed  Google Scholar 

  31. Noh SA, Lee HS, Kim YS, Peak KH, Shin JS, Bae JM. Down-regulation of the IbEXP1 gene enhanced storage root development in sweet potato. J Exp Bot. 2013;64(1):129–42.

    Article  CAS  PubMed  Google Scholar 

  32. Bishopp A, Help H, El-Showk S, Weijers D, Scheres B, Friml J, Benková E, Mähönen AP, Helariutta Y. A mutually inhibitory interaction between auxin and cytokinin specifies vascular pattern in roots. Curr Biol. 2011;21:917–26.

    Article  CAS  PubMed  Google Scholar 

  33. Li LL, Shao TY, Yang H, Chen M, Gao X, Long X, Shao H, Liu Z, Rengel Z. The endogenous plant hormones and ratios regulate sugar and dry matter accumulation Jerusalem artichoke in salt-soil. Sci Total Environ. 2017;578:40–6.

    Article  CAS  PubMed  Google Scholar 

  34. Rahman M, Isenberg F. The role of exogenous plant regulation in the dormancy of onion bulbs. J Agri Sci. 1974;82(1):113–6.

    Article  Google Scholar 

  35. Roumeliotis E, Kloosterman B, Oortwijn M, Kohlen W, Bouwmeester HJ, Visser RG, Bachem CW. The effects of auxin and strigolactones on tuber initiation and stolon architecture in potato. J Exp Bot. 2012;63(12):4539–47.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Ševčíková H, Mašková P, Tarkowská D, Mašek T, Lipavská H. Carbohydrates and gibberellins relationship in potato tuberization. J Plant Physiol. 2017;214:53–63.

    Article  PubMed  Google Scholar 

  37. Dong T, Zhu M, Yu J, Han R, Tang C, Xu T, Liu J, Li Z. RNA-Seq and iTRAQ reveal multiple pathways involved in storage root formation and development in sweet potato (Ipomoea batatas L). BMC Plant Biol. 2019;19(1):1–16.

    Article  Google Scholar 

  38. Noh SA, Lee HS, Huh EJ, Huh GH, Paek KH, Shin JS, Bae JM. SRD1 is involved in the auxin-mediated initial thickening growth of storage root by enhancing proliferation of metaxylem and cambium cells in sweet potato (Ipomoea batatas). J Exp Bot. 2010;61(5):1337–49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Nieminen K, Immanen J, Laxell M, Kauppinen L, Tarkowski P, Dolezal K, Tähtiharju S, Elo A, Decourteix M, Ljung A, Bhalerao R, Keinonen K, Albert VA, Helariutta Y. Cytokinin signaling regulates cambial development in poplar. Proc Natl Acad Sci U S A. 2008;105(50):20032–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Li M, Yang Y, Li X, Gu L, Wang F, Feng F, Tian Y, Wang F, Wang X, Lin W, Chen X, Zhang Z. Analysis of integrated multiple ’omics’ datasets reveals the mechanisms of initiation and determination in the formation of tuberous roots in Rehmannia glutinosa. J Exp Bot. 2015;66(19):5837–51.

    Article  CAS  PubMed  Google Scholar 

  41. Yu R, Wang J, Xu L, Wang Y, Wang R, Zhu X, Sun X, Luo X, Xie Y, Everlyne M, Liu L. Transcriptome profling of taproot reveals complex regulatory networks during taproot thickening in radish (Raphanus sativus L). Front Plant Sci. 2016;7:1210.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Yao S, Lan Z, Huang R, Tan Y, Huang D, Gu J, Pan C. Hormonal and transcriptional analyses provides new insights into the molecular mechanisms underlying root thickening and isofavonoid biosynthesis in Callerya speciosa (champ. Ex Benth.) Schot. Sci Rep. 2021;11(1):9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Hejátko J, Ryu H, Kim GT, Dobesová R, Choi S, Choi SM, Soucek P, Horák J, Pekárová B, Palme K, Brzobohaty B, Hwang I. The histidine kinases CYTOKININ-INDEPENDENT 1 and ARABIDOPSIS HISTIDINE KINASE 2 and 3 regulate vascular tissue development in Arabidopsis shoots. Plant Cell. 2009;21(7):2008–21.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Kumar D, Wareing PF. Factors controlling stolon development in the potato plant. New Phytol. 1972;71(4):639–48.

    Article  CAS  Google Scholar 

  45. Yasunori K. Changes in levels of butanol- and water-soluble cytokinins during the life cycle of potato tubers. Plant Cell Physiol. 1982;23(5):843–9.

    Google Scholar 

  46. Lv Y, Li Y, Liu X, Xu K. Identification of ginger (Zingiber Officinale Roscoe) reference genesfor gene expression analysis. Front Genet. 2020;11:586098.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Liu M, Bassetti N, Petrasch S, Zhang N, Bucher J, Shen S, Zhao J, Bonnema G. What makes turnips: anatomy, physiology and transcriptome during early stages of itshypocotyl-tuber development. Hortic Res. 2019;6:38.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Liu D, Ma L, Zhou Z, Liang Q, Xie Q, Ou K, Liu Y, Su Y. Starch and mineral element accumulation during root tuber expansion of Pueraria thomosonii Benth. Food Chem. 2021;343:128445.

    Article  CAS  PubMed  Google Scholar 

  49. Abelenda JA, Navarro C, Prat S. From the model to the crop: genes controlling tuber formation in potato. Curr Opin Biotechnol. 2011;22(2):287–92.

    Article  CAS  PubMed  Google Scholar 

  50. Zhang K, Wu Z, Tang D, Luo K, Lu H, Liu Y, Dong J, Wang X, Lv C, Wang J, Lu K. Comparative transcriptome analysis reveals critical function of sucrose metabolism related-enzymes in starch accumulation in the storage root of sweet potato. Front Plant Sci. 2017;8:914.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Mitsui Y, Shimomura M, Komatsu K, Namiki N, Shibata-Hatta M, Imai M, Katayose Y, Mukai Y, Kanamori H, Kurita K, Kagami T, Wakatsuk A, Ohyanagi H, Ikawa H, Minaka N, Nakagawa K, Shiwa Y, Sasaki T. The radish genome and comprehensive gene expression profile of tuberous root formation and development. Sci Rep. 2015;5:10835.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Wang X, Chang L, Tong Z, Wang D, Yin Q, Wang D, Jin X, Yang Q, Wang L, Yong S, Huang Q, Guo A, Peng M. Proteomics profiling reveals carbohydrate metabolic enzymes and 14-3-3 proteins play important roles for starch accumulation during cassava root tuberization. Sci Rep. 2016;6:19643.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Solís-Guzmán M, Argüello-Astorga G, López-Bucio J, Ruiz-Herrera L, López-Meza J, Sánchez-Calderón L, Carreón-Abud Y, Martínez-Trujillo M. Arabidopsis thaliana sucrose phosphate synthase (sps) genes are expressed differentially in organs and tissues, and their transcription is regulated by osmotic stress. Gene Expr Patterns. 2017;25:92–101.

    Article  PubMed  Google Scholar 

  54. Yu R, Jing W, Liang X, Yan W, W R, Zhu X, Sun X, Luo X, Xie Y, Everlyne M, Liu L. Transcriptome profiling of taproot reveals complex regulatory networks during taprootthickening in radish (Raphanus sativus L). Front Plant Sci. 2016;7:1210.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Duan Y, Yang L, Zhu H, Zhou J, Sun H, Gong H. Structure and expression analysis of sucrose phosphate synthase, sucrose synthase and invertase gene families in Solanum lycopersicum. Int J Mol Sci. 2021;22(9):4698.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Li XJ, Yang JL, Hao B, Lu YC, Qian ZL, Li Y, Ye S, Tang J, Chen M, Long G, Zhao Y, Zhang G, Chen J, Fan W, Yang S. Comparative transcriptome and metabolome analyses provide new insights into the molecular mechanisms underlying taproot thickening in Panax notoginseng. BMC Plant Biol. 2019;19(1):451.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Gupta AK, Singh J, Kaur N. Sink development, sucrose metabolizing enzymes and carbohydrate status in turnip (Brassica rapa L). Acta Physiol Plant. 2001;23:31–6.

    Article  CAS  Google Scholar 

  58. Li XQ, Zhang D. Gene expression activity and pathway selection for sucrose metabolism in developing storage root of sweet potato. Plant Cell Physiol. 2003;44(6):630–6.

    Article  CAS  PubMed  Google Scholar 

  59. Luo X, Xu G, Li Y, Li B, Tu J. Kuzu root (Ge-Gen) regulates on glucose and lipid metabolism to ameliorating insulin resistance on 3T3-L1 adipocytes. Zhongguo Zhong Yao Za Zhi. 2016;41(14):26876–2694.

    Google Scholar 

  60. Begović L, Pospihalj T, Lončarić P, Štolfa Čamagajevac I, Cesar V, Leljak-Levanić D. Distinct accumulation and remobilization of fructans in barley cultivars contrasting for photosynthetic performance and yield. Theor Exp Plant Phys. 2020;32:1–12.

    Google Scholar 

  61. Wu Y, Shen Y. Dormancy in Tilia miqueliana is attributable to permeability barries and mechanical constraints in the endosperm and seed coat. Braz J Bot. 2021;3:725–40.

    Article  Google Scholar 

  62. Yang S, Na C. Effect of two kinds of antistaling agent on the pigments of cutting branch for Cornus alba L. Biotechnol J Int. 2019;22(4):1–7.

    CAS  Google Scholar 

  63. Hsiao H, Chen RLC, Chou C, Cheng T. Hand-held colorimetry sensor platform for determining salivary α-amylase activity and its applications for stress assessment. Sens (Basel). 2019;19(7):1571.

    Article  CAS  Google Scholar 

  64. Yang J, Zhang J, Wang Z, Zhu Q. Hormonal changes in the grains of rice subjected to water stress during grain filling. Plant Physiol. 2001;127:315–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Wang H, Zhou Y, Gilmer S, Whitwill S, Fowke LC. Expression of the plant cyclin-dependent kinase inhibitor ICK1 affects cell division, plant growth and morphology. Plant J. 2000;24:613–23.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

None.

Funding

This research was funded by Jiangxi Modern Agricultural Industrial Technology System Construction Project (JXARS-16) and the National Natural Science Foundation of China (31860560).

Author information

Authors and Affiliations

Authors

Contributions

X.X., W.C., Z.W., and F.S. designed the research. H.Y., Z.H., and G.F. carried out the experiments. X.X., S.S., H.Y., and F.S. analyzed the data. X.X. and S.S. wrote the manuscript. All authors have read and approved the manuscript.

Corresponding authors

Correspondence to Xiao Xufeng, Wu Caijun or Fan Shuying.

Ethics declarations

Ethics approval and consent to participate

The plant samples used in this study were grown and collected at Jiangxi Agricultural University, Nanchang, China. The sample collection complies with relevant institutional, national, and international guidelines and legislation. In addition, the deposition number of voucher specimens was PTRP1, PTRP3, and PTRP6, respectively.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have 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 4.0 International License, which permits use, sharing, adaptation, 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 changes were made. 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/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Xufeng, X., Yuanfeng, H., Ming, Z. et al. Transcriptome profiling reveals the genes involved in tuberous root expansion in Pueraria (Pueraria montana var. thomsonii). BMC Plant Biol 23, 338 (2023). https://doi.org/10.1186/s12870-023-04303-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-023-04303-x

Keywords