Skip to main content

Morphological, anatomical, and transcriptomics analysis reveals the regulatory mechanisms of cassava plant height development

Abstract

Background

Cassava is one of three major potato crops and the sixth most important food crop globally. Improving yield remains a primary aim in cassava breeding. Notably, plant height significantly impacts the yield and quality of crops; however, the mechanisms underlying cassava plant height development are yet to be elucidated.

Results

In this study, we investigated the mechanisms responsible for cassava plant height development using phenotypic, anatomical, and transcriptomic analyses. Phenotypic and anatomical analysis revealed that compared to the high-stem cassava cultivar, the dwarf-stem cassava cultivar exhibited a significant reduction in plant height and a notable increase in internode tissue xylem area. Meanwhile, physiological analysis demonstrated that the lignin content of dwarf cassava was significantly higher than that of high cassava. Notably, transcriptome analysis of internode tissues identified several differentially expressed genes involved in cell wall synthesis and expansion, plant hormone signal transduction, phenylpropanoid biosynthesis, and flavonoid biosynthesis between the two cassava cultivars.

Conclusions

Our findings suggest that internode tissue cell division, secondary wall lignification, and hormone-related gene expression play important roles in cassava plant height development. Ultimately, this study provides new insights into the mechanisms of plant height morphogenesis in cassava and identifies candidate regulatory genes associated with plant height that can serve as valuable genetic resources for future crop dwarfing breeding.

Peer Review reports

Introduction

Plant height is a crucial trait influencing crop yield and quality. Notably, dwarf crop breeding is a key strategy used to enhance productivity. A reduction in plant height facilitates mechanized harvesting, increases yield and photosynthetic efficiency, enhances fertilizer tolerance, supports dense planting, improves lodging resistance, and enhances nutrient transport to reproductive organs. Moreover, dwarfing reduces the need for chemical growth regulators, thereby protecting the ecological environment. Therefore, understanding the regulatory mechanisms and genetic basis of plant height traits remains a major topic in genetic breeding research. This focus has been particularly extensive in rice, with the application of rice and wheat dwarf varieties significantly increasing crop yield in the 1960s and 1970s, leading to the Green Revolution [1].

Alongside environment factors, gene regulation plays a leading role in crop plant height development. For example, overexpression of the brassinosteroid (BR) degradation–related gene AtBAT1 in transgenic creeping bentgrass was found to induce dwarfing, shorter internodes, smaller leaf angles, and other BR deficiency–related phenotypes [2]. Meanwhile, overexpression of CYP94C2b in rice increases internode length and plant height [3]. Similarly, the GmLHY gene regulates soybean plant height development through regulation of the gibberellin (GA) pathway [4]. In wheat, overexpression of Rht1 and Rht2 genes not only reduces plant height, but also increases tiller numbers and grains per spike, thereby significantly improving yield [1, 5].

Numerous studies have highlighted associations between crop dwarfing and deletions or mutations of genes involved in GA, BR, auxin (IAA), and strigolactone (SL) biosynthesis or signal transduction pathways [6]. For example, mutations in the Arabidopsis BRI1 (Brassinosteroid-Insensitive 1) gene, which encodes a receptor-like kinase, results in the inability to detect BR signals, ultimately leading to a severe dwarf phenotype [7]. In rice, mutations in the sd1 gene, which encodes GA20ox, a key enzyme in the GA synthesis pathway, lead to a reduction in GA content and the shortening of internode length, resulting in a dwarf phenotype [8,9,10]. Similarly, deletion of rice IAA co-receptor genes OsTIR1 and OsAFB2-5 reduce plant height to varying degrees, accompanied by changes in yield, tillering, root traits [11]. Meanwhile, abnormal expression of the key IAA delivery genes ABCB1/PGP1 in the maize mutant brachytic2 (br2) and sorghum mutant dwarf3 (d3) resulted in a dwarfing phenotype [12]. SL has also been found to mediate rice plant height by regulating D53 degradation via the D14-D3 complex, thereby reversing the repression of downstream target genes [13].

Other factors influencing plant dwarfing include transcription factors and cell wall–associated genes. For example, CesA and CSL genes are responsible for cell wall biosynthesis, while expansin (EXP), and XTH genes influence cell wall loosening, thereby affecting cell growth during stem elongation [14]. Notably, plants with CesA or CSL mutations exhibit dwarf phenotypes [15,16,17,18,19]. Meanwhile inhibition of EXP gene expression results in reduced plant height, early flowering, and leaf curl [20,21,22,23]. Transcription factors can respond to both external environmental factors and endogenous hormone signals, regulating plant growth by influencing downstream target genes [24,25,26]. For example, overexpression of NAC transcription factors has been found to reduce abscisic acid (ABA) and BR content in apples, resulting in shorter internodes [27]. However, despite significant advances in understanding the genetic regulation of plant height and identifying several associated genes, the functions and regulatory mechanisms of many plant height–associated genes remains unclear due to the complex nature of this trait. Therefore, further research is needed to explore the roles and regulatory mechanisms of various genes influencing plant height.

Cassava (Manihot esculenta Crantz), a plant belonging to the Manihot genus in the Euphorbiaceae family, ranks among the world’s three top tuber crops alongside potato and sweet potato. With an annual yield exceeding 200 million tons, cassava serves as a staple food for over 700 million people worldwide, making it the sixth most important food crop globally. Beyond its role in food security, cassava serves as an important source of industrial raw material, feed, biomass energy, and food processing [28, 29]. Consequently, cassava has a strong reputation, being termed the “king of starch” and “energy crop”. Notably, cassava exhibits resilience to barren, dry, and acidic soil and can be planted on slope land. Moreover, this crop offers various advantages, including the primary use of asexual reproduction, high photosynthetic efficiency, extensive management, low cost, and high yield potential, making it a promising crop for industrial development. However, research on cassava plant height development remains sparse. Therefore, this study employed two cassava varieties with significant differences in height to investigate the mechanisms underlying cassava plant height development by analyzing morphological structure, anatomical structure, lignin content, cellulose content, and gene expression.

Results

Phenotypic characteristics, cytology, and physiology of XX048 and NZ199 plants

The plant heights of the two cassava varieties XX048 and NZ199 were 420.27 ± 5.2 cm and 219 ± 4.57 cm, respectively; this difference was determined to be statistically significant (Fig. 1). August–September is considered the period of rapid growth in cassava plant height. During this period, XX048 plant height increased by nearly 150 cm, while NZ199 plant height increased by approximately 55 cm. This indicated that XX048 plants grew significantly faster than NZ199 plants, with the difference in height between the two variants being the largest during this period.

To identify the anatomical differences between XX048 and NZ199 plants, we performed paraffin section staining of the topmost internode tissue. The woody area of XX048 was determined to be 414.87 ± 14.28 mm2, while that of NZ199 was 911.10 ± 82.62 mm2; notably, this difference in woody area was statistically significant. Meanwhile, the average length and width of XX048 xylem cells were 27.04 ± 4.37 mm and 22.11 ± 4.36 mm, respectively, while the average length and width of NZ199 xylem cells were 35.25 ± 3.81 mm and 29.12 ± 5.15 mm, respectively (Fig. 2). However, there was no significant difference observed between the length and width of xylem cells in XX048 and NZ199 plants. Overall, these results indicate that the xylem of NZ199 is more developed and has a higher degree of lignification.

Fig. 1
figure 1

Comparison of the phenotypic characteristics of XX048 and NZ199 cassava plants. A: Plant heights of XX048 and NZ199 variants. Scale bars = 0.1 m. B: Plant heights of XX048 and NZ199 across six developmental stages; orange denotes NZ199, while blue represents XX048. Error bars indicate the SD. *** denotes a significant difference of p < 0.001

Fig. 2
figure 2

Anatomical characteristics of stem sections obtained from XX048 and NZ199 cassava plants. A-D: Transverse sections showing the a: epidermis, b: cortex, c: xylem, and d: microtubule column. E-G Analysis of stem sections in each variant. Error bars indicate the SD. Different letters denote significant differences between XX048 and NZ199 plants (P < 0.05)

Next, we measured lignin and cellulose content in the topmost internode tissue of XX048 and NZ199 cassava plants (Fig. 3). The lignin content of NZ199 and XX048 plants is 162.37 ± 1.81 mg/g and 132.01 ± 3.83 mg/g, respectively. Meanwhile, the cellulose content of NZ199 and XX048 is 73.47 ± 2.69 mg/g and 70.91 ± 1.45 mg/g, respectively. Notably, the difference in lignin content between the two variants was determined to be statistically significant, while the difference in cellulose content was not significant. These results indicate differences in the cell wall lignification development in the two cassava variants, with the main stem of NZ199 cassava entering secondary wall lignification development earlier. Notably, these findings are consistent with the cytological results.

Fig. 3
figure 3

Differences in lignin and cellulose content between XX048 and NZ199 cassava plants. Error bars indicate the SD. Different letters associated with the same trait denote significant differences between XX048 and NZ199 variants (P < 0.05)

Transcriptomic analysis

To explore the key genes linked to plant height development in cassava, six cDNA libraries were constructed. After removing splices and low-quality reads from Illumina sequencing data, 121,791,451 clean reads were obtained (Supplementary Table 1); the percentages of Q20 and Q30 base rates were > 97.82% and > 93.65%, respectively. The clean reads of each sample were then compared with the cassava reference genome, and the corresponding comparison efficiency ranged from 94.56 to 94.88%. The sample correlation heat map indicated that the R2 value among three biological repeat samples was > 0.89, with that of most of the samples being > 0.9, indicating that this experiment was highly reproducible and the data were reliable (Fig. 4).

Fig. 4
figure 4

Correlation analysis between sample replicates of XX048 and NZ199 cassava variants

To investigate the expression pattern of genes associated with cassava plant height development, we compared the transcriptome maps of cassava with tall and short stems. Differentially expressed genes (DEGs) between cassava with tall and short stems were identified according to a log2 fold change ≥ 2 and a p – value < 0.01. A total of 2368 DEGs were identified between XX048 and NZ199 variants, of which 1112 genes were up-regulated and 1256 genes were down-regulated (Fig. 5).

Fig. 5
figure 5

Volcano map showing differentially expressed gene expression between XX048 and NZ199 cassava plants

To understand the biological processes associated with the identified DEGs, we classified DEGs into GO functions. These DEGs were categorized into three main GO functional groups (Fig. 6): biological process (BP), cellular component (CC), and molecular function (MF). The three most abundant sub-categories in the BP group were “metabolic process” (753 DEGs), “cellular process” (590 DEGs), and “single-organism process” (458 DEGs). In the CC group, the four most abundant sub-categories were “membrane” (521 DEGs), “membrane part” (480 DEGs), “cell part” (377 DEGs), and “cell” (377 DEGs). Meanwhile, the majority of DEGs in the MF group were associated with “binding” (829 DEGs), “catalytic activity” (795 DEGs), and “transporter activity” (105 DEGs).

Fig. 6
figure 6

Gene Ontology (GO) classification of DEGs between XX048 and NZ199 cassava variants

Finally, KEGG enrichment analysis revealed that the five most abundant pathways associated with these DEGs include “plant-pathogen interaction” (117 DEGs), “plant hormone signal transduction” (47 DEGs), “starch and sucrose metabolism” (46 DEGs), “phenylpropanoid biosynthesis” (43 DEGs), and “carbon metabolism” (41 DEGs). Notably, the “flavonoid biosynthesis”, “brassinosteroid biosynthesis”, and “zeatin biosynthesis” pathways were also enriched, with the flavonoid biosynthesis pathway being the most significantly enriched among these biosynthetic pathways (Fig. 7).

Fig. 7
figure 7

KEGG enrichment analysis of DEGs between XX048 and NZ199 cassava plant variants. The x axis denotes the enrichment score, which represents the ratio of DEGs identified in a specific pathway to all the genes contained in this pathway. The y axis denotes each KEGG pathway. Bubble size indicates the number of genes annotated to that specific KEGG Pathway

DEGs associated with plant hormone biosynthesis and signal transduction pathways

Plant hormone biosynthesis and signal transduction include IAA, GA, cytokinin, BR, ABA, ethylene (ETH), jasmonic acid (JA), and salicylic acid (SA) pathways. Notably, various DEGs involved in IAA, cytokinin, GA, and BR biosynthesis and signal transduction were found to be significantly down-regulated in NZ199 dwarf cassava (Fig. 8, Supplementary Table 2). For example, in IAA biosynthesis and transduction, auxin response factor (ARF), auxin-responsive GH3 gene family (GH3), and small auxin-up RNA (SAUR) were typically down-regulated in NZ199 plants. Meanwhile, in the cytokinin pathway, two-component response regulators in ARR-A and ARR-B families were predominantly down-regulated in NZ199 cassava. In the GA pathway of NZ199 cassava, one DELLA protein was up-regulated and another was down-regulated, while phytochrome-interacting factor 3 (TF) was down-regulated and the GA receptor GID1 was up-regulated. Several BR biosynthesis and signal transduction genes were also down-regulated in NZ199 cassava, including brassinosteroid insensitive 1-associated receptor kinase 1 (BAK1), protein brassinosteroid insensitive 1 (BRI1), and cyclin D3 (CYCD3).

Conversely, many important DEGs involved in ETH, SA, JA, and ABA biosynthesis and signal transduction were significantly up-regulated in dwarf cassava. For example, in the ABA pathway, serine/threonine-protein kinase SRK2 (SnRK2) and ABA-responsive element binding factor (ABF) were typically up-regulated in NZ199 dwarf cassava. Additionally, in the JA pathway, the transcription factor MYC2 was predominantly up-regulated in NZ199 plants. In the ETH and SA pathways, the serine/threonine-protein kinase CTR1, mitogen-activated protein kinase kinase 4/5 (MKK4/5), regulatory protein NPR1, and pathogenesis-related protein 1 (PR1) were up-regulated.

Fig. 8
figure 8

Heatmaps showing the DEGs involved in hormone signal transduction during cassava plant height growth in XX048 and NZ199 cassava variants

DEGs associated with flavonoid biosynthesis and phenylpropanoid biosynthesis

Overall, flavonoid biosynthesis was the most significantly enriched pathway in KEGG enrichment analysis. In this pathway, most associated genes were significantly up-regulated in the XX048 vs. NZ199 groups. Up-regulated DEGs included 4 chalcone synthase genes (CHS, EC: 2.3.1.74), 5 shikimate O-hydroxycinnamoyltransferase (HCT, EC: 2.3.1.133), 1 phlorizin synthase genes (PGT1, EC: 2.4.1.357), 2 chalcone isomerase (EC: 5.5.1.6), 1 naringenin 3-dioxygenase (F3H, EC: 1.14.11.9), 1 bifunctional dihydroflavonol 4-reductase/flavanone 4-reductase (DFR, EC: 1.1.1.219, EC: 1.1.1.234), 1 flavonoid 3’,5’-hydroxylase (CYP75A, EC: 1.14.14.81), 1 flavonol synthase (FLS, EC: 1.14.20.6), 1 anthocyanidin synthase (ANS, EC: 1.14.20.4), 1 anthocyanidin reductase (ANR, EC: 1.3.1.77), and 1 leucoanthocyanidin reductase (LAR, EC: 1.17.1.3). Meanwhile, only 1 chalcone synthase gene (CHS, EC: 2.3.1.74), 1 caffeoyl-CoA O-methyltransferase gene (EC: 2.1.1.104), and 5 phlorizin synthase genes (PGT1, EC: 2.4.1.357) were significantly down-regulated in the XX048 vs. NZ199 groups (Fig. 9, Supplementary Table 3).

Fig. 9
figure 9

Mechanisms of flavonoid biosynthesis in cassava. HCT, shikimate O-hydroxycinnamoyl transferase; CHS, chalcone synthase; PGT1, phlorizin synthase; CYP75A, flavonoid 3’,5’-hydroxylase; F3H, naringenin 3-dioxygenase; DFR, bifunctional dihydroflavonol 4-reductase/flavanone 4-reductase; FLS, flavonol synthase; LAR, leucoanthocyanidin reductase; ANS, anthocyanidin synthase; ANR, anthocyanidin reductase. Red and green indicate up- and downregulated genes, respectively; blue indicates some upregulated and some downregulated genes

Phenylpropanoid biosynthesis is another key pathway identified in KEGG enrichment analysis. The products of this pathway primarily include lignin monomers, including p-coumaroyl, coniferyl alcohol, and sinapyl alcohol. Notably, our study identified 37 DEGs involved in lignin monomer synthesis (Fig. 10, Supplementary Table 4), including 4 phenylalanine ammonia-lyase (PAL), 6 4-coumarate-CoA ligase (4CL), 5 cinnamoyl-CoA reductase (CCR), 4 cinnamyl-alcohol dehydrogenase (CAD), 2 peroxidases (POD), 5 HCT, 5 caffeic acid 3-O-methyltransferase/acetylserotonin O-methyltransferase (COMT), 1 caffeoyl-CoA O-methyltransferase (CCoAOMT), 1 ferulate-5-hydroxylase (F5H), and 4 laccase genes. Of these 37 DEGs, 28 were significantly up-regulated in the XX048 vs. NZ199 groups, while only 9 genes were significantly down-regulated.

Fig. 10
figure 10

Heat maps depicting the differential expression patterns of genes involved in lignin biosynthesis in the XX048 and NZ199 cassava. PAL, phenylalanine ammonia-lyase; 4CL, 4-coumarate–CoA ligase; CCR, cinnamoyl-CoA reductase; CAD, cinnamyl-alcohol dehydrogenase; POD, peroxidase; HCT, shikimate O-hydroxycinnamoyltransferase; COMT, caffeic acid 3-O-methyltransferase/acetylserotonin O-methyltransferase; CCoAOMT, caffeoyl-CoA O-methyltransferase

DEGs associated with cell wall biosynthesis and expansion

The synthesis of new cell wall components is essential to provide mechanical strength to newly elongated cell walls and internodes. In this study, we analyzed the differential expression of genes involved in cell wall synthesis or modification since they may be linked to cassava stem elongation and expansion; such genes include cellulose synthase (CES), EXP, xylosidase (XYL), pectin lyase (PL), pectin esterase (PM), pectin acetyl esterase (PA), and polygalacturonase (PG).

EXPs, nonhydrolytic cell wall-loosening proteins, enable cell expansion while facilitating differentiation and growth. In the present study, two members of the EXP gene family were found to be differentially expressed in tall and dwarf cassava, with significant up-regulation in NZ199 dwarf cassava (Fig. 11, Supplementary Table 5). Xyloglucan endoglycosyltransferase/hydrolase (XTH) also participates in cell expansion by promoting the relaxation and rearrangement of cell wall fibers in growing tissues. Notably, four XTHs were differentially expressed between tall and dwarf cassava, of which 3 were down-regulated and 1 was up-regulated in the tall vs. dwarf groups (Fig. 11, Supplementary Table 5). The degradation of pectin, a complex polymer, has been linked to the initiation of secondary wall lignification and cell wall disintegration. A total of 11 pectin-related genes were significantly differentially expressed between tall and dwarf cassava; this included 5 pectin esterase genes (PM), 3 polygalacturonase genes (PG), 2 pectin lyase genes (PL), and 1 pectin acetyl esterase gene (PA). Additionally, 1 cellulose synthase gene and 2 XYL genes were differentially expressed between the tall and dwarf cassava variants (Fig. 11, Supplementary Table 5).

Fig. 11
figure 11

Heatmap showing the differential expression patterns of genes related to cell wall biosynthesis and expansion in the XX048 and NZ199 cassava variants. XTH, xyloglucan endotransglucosylase/hydrolase; CES, cellulose synthase; XYL, xylosidase; PL, pectate lyase; PM, pectinesterase; PA, pectin acetylesterase; PG, polygalacturonase; EXP, expansin

Transcription factor analysis

In this study, we identified 92 differentially expressed transcription factors between XX048 and NZ199 cassava, across 34 transcription factor families. Among these transcription factors, 56 were down-regulated and 36 were up-regulated (Supplementary Table 6). Notably, these transcription factors were primarily comprised of the bHLH (9), MYB (9), C2H2 (7), NAC (5), MYB-related (5), bZIP (5), AP2 / ERF-ERF (5), and WRKY (4) transcription factor families.

qRT-PCR validation

To validate the accuracy of RNA-Seq results, 11 DEGs were randomly selected for qRT-PCR verification. Overall, the expression patterns of these 11 DEGs were similar to those identified by RNA-Seq, with an R2 value of 0.9804 (Fig. 12), indicating the reliability of our RNA-Seq results.

Fig. 12
figure 12

qRT-PCR validation of the transcriptome data. a-k qRT-PCR validation profles of eleven randomly selected genes. The data was normalized by using Actin as an internal reference. The expression level of XX048 was used as reference state, which was set to 1, and fold change values were shown here. l Pearson’s correlation analysis of the fold changes between qRT-PCR (x axis) and RNA-seq (y axis) data was conducted (R2)

Discussion

Plant height is an important agronomic trait affecting crop yield and quality. Therefore, understanding the molecular mechanism controlling crop plant height is of great significance in the agricultural sector. However, although cassava is considered one of three major potato crops and one of six major food crops globally, the genetic regulation mechanism of cassava plant height is unclear. Consequently, the regulation and genetic improvement of cassava plant height traits remain limited. Therefore, this study aimed to identify candidate genes and gene expression networks that are responsible for regulating cassava plant height development.

Cell wall biosynthesis and expansion play a major role in plant height

Plant cell walls not only provide mechanical support for cell components and protects against various external stimuli, but is also essential for cell growth. Lignin and cellulose are key components of cell walls; therefore, their content is considered the primary factor affecting plant height. Notably, prior studies have established that lower lignin content is associated with increased plant height [30]. Changes in lignin content have been linked to the expression of genes related to lignin biosynthesis. Genes involved in lignin biosynthesis, such as PAL, HCT, 4CL, CCR, and CAD, play crucial regulatory roles, with lignin synthesis being positively correlated with the expression of these genes [31, 32]. In the present study, various PAL, HCT, 4CL, CCR, and CAD genes were found to be significantly up-regulated in dwarf cassava (Fig. 10), consistent with the observed physiological results indicating higher lignin content in dwarf cassava compared to tall cassava (Fig. 3). Similar results have also been observed in rice, with the stems of dwarf rice mutants displaying higher lignin content and higher lignin-related gene expression levels [33].

POD and laccase are also key components involved in lignin biosynthesis. Specifically, these enzymes catalyze cross-linking between phenolic groups in cell wall structural proteins, pectin, and other cell wall polymers, thereby inducing cell wall hardening, inhibiting cell elongation, and facilitating cell wall lignification. Notably, the present study revealed significant up-regulation in POD and laccase gene expression in dwarf cassava (Fig. 10). Ultimately, these results demonstrated that the cell wall lignification development of the two cassava varieties differed. In particular, the dwarf cassava entered the secondary wall lignification development state relatively earlier, resulting in the termination of cell elongation and the development of a dwarf phenotype.

Expansion is a process involving non-hydrolytic cell wall relaxation of protein in plant cell walls. It is considered to play an important role in regulating cell wall ductility and is a key factor in determining cell expansion [22]. The role of EXPs has been confirmed in various plant species. For example, overexpression of OsEXP4 in transgenic rice can significantly increase plant height [23]. Similarly, overexpression of ClEXPA2 in tobacco also significantly increases internode length [34]. However, although EXPs have long been thought to play an active role in promoting plant height and tissue expansion, this is not always the case. For example, Rochange et al. demonstrated that EXP overexpression in transgenic plants resulted in the development of smaller leaves and shorter internodes due to the reduced sensitivity of the cell wall to EXP [35]. Notably, a similar phenomenon was also observed in the current study. Specifically, two EXP genes were significantly up-regulated in dwarf cassava, indicating the normal transcription of EXPs (Fig. 11). Therefore, we hypothesize that a decrease in cell wall sensitivity to EXP is one of the primary factors affecting plant height dwarfing in cassava.

The XTH gene encodes a cell wall modifying enzyme that is found in higher plant cells; specifically, this enzyme alters cell wall structure by catalyzing the hydrolysis or transfer of xyloglucan molecules. Prior studies have revealed that XTH expression is closely related to plant height development. For example, OsXTH8 expression is positively correlated with plant height in rice [36]. Meanwhile, mutations in AtXTH9 have been found to result in shorter internodes length in Arabidopsis spp. [37]. Similarly, our results demonstrated that the expression level of XTH genes was significantly higher in high-stem cassava than that in dwarf-stem cassava (Fig. 11), indicating that these genes may be involved in cassava plant height development.

Role of plant hormones in plant height development

As a physiological signal regulating plant growth and development, plant hormones play an important role in plant morphogenesis, growth, and metabolism. Notably, previous studies have revealed that dwarfing mutations are closely associated with plant hormones.

IAA, a key growth regulator, plays an important role in all stages of plant growth and development. Thus, abnormal expression of IAA transport key genes can lead to plant dwarfing. In the current study, the IAA response gene CH3 (Manes.12G087200.v7.0) was found to be significantly down-regulated in dwarf cassava, indicating a reduction in the IAA response. ARF transcription factors also mediate IAA responses and modulate stem development by regulating cell division [38]. In the present study, ARF expression was significantly down-regulated in dwarf cassava (Fig. 8), indicating inhibition of cell division. SAUR, the largest family of IAA response genes, also plays a key role in dwarfing; notably the gain-of-function of SAUR19–24 in seedlings has been found to result in increased hypocotyl length and leaf size, while the loss-of-function of lead to reduced hypocotyl length and leaf size [39]. Similarly, in the present study, a large number of SAUR genes were significantly down-regulated in dwarf cassava (Fig. 8); similar results have been observed in Sophora davidii and peanut plants [40, 41]. Overall, our findings indicate that these genes influence plant height during cassava growth and development.

Notably, GAs also significantly influence plant stem elongation. The plant hormone GA, the receptor GID1, and the inhibitor DELLA collectively represent the GA signal transduction pathway When activated by GA, the GID1 receptor binds to DELLA through the E3 ubiquitination pathway, leading to DELLA protein degradation and the obstruction of DELLA-mediated growth inhibition [42]. Thus, a loss-of-function mutation in the GID1 gene has been found to result in a dwarf phenotype [43]. Conversely, in the present study, GID1 was significantly up-regulated in dwarf cassava (Fig. 8). This contrasts with the understanding that loss-of-function mutations in GID1 cause plants to exhibit dwarf phenotypes and GA insensitivity [42]. Nonetheless, our finding aligns with the established GA-triggered negative feedback loop observed in in Lepidium sativum that results in the repression of GID1 gene transcription [44]. Similar results have been observed in studies on bamboo. For example, higher expression levels of PeGID1 have been detected in the internodes of dwarf Moso bamboo than in wild-type bamboo [45]. DELLA is a transcription regulatory factor that responds to GA and acts as an inhibitory factor for GA-induced growth. In a prior study, the DELLA protein gene PhSLR1 in moso bamboo was transferred into Arabidopsis thaliana; notably, these transgenic plants were found to be shorter than their wild-type counterparts [46]. In the present study, our findings aligned with those of previous studies, with one DELLA protein gene being up-regulated in dwarf cassava (Fig. 8). Therefore, cassava plant height dwarfing may be attributed to GA signal transduction inhibition by DELLA protein.

BRs are crucial positive regulators of stem elongation. Thus, the deletion or mutation of genes, transcription factors, and enzymes related to BR synthesis and signal transduction pathways reduce plant cell elongation, resulting in dwarfing. BRI1 is the primary component involved in BR signal detection; therefore, mutations in the BRI1 gene in Arabidopsis and rice spp. result in the development of dwarfing phenotypes [47,48,49]. Additionally, BAK1 also plays a key role in BR signal transduction. Consequently, BAK1 loss-of-function mutations have been shown to reduced BR sensitivity, leading to the development of shorter stems [50]. Meanwhile, inhibition of CYCD3, a D-type plant cyclin gene that promotes plant cell division, has been found to result in significant dwarfing in Arabidopsis spp. [51]. In the present study, the expression levels of 6 BRI1s, 3 BAK1s, and 1 CYCD3 genes were down-regulated in dwarf cassava (Fig. 8), indicating a reduction in BR sensitivity in dwarf cassava; therefore, the inhibition of cell elongation may serve as an important factor in cassava dwarfing.

Role of flavonoids in plant height development

In recent years, many studies have shown the relationship between plant height development is related to flavonoids. Transcriptome and metabolome analysis has identified a series of genes involved in the flavonoid biosynthesis pathway that are differentially expressed in dwarf and tall mutants of Sophora davidii. A total of 8 differential metabolites were identified between these mutants, of which 6 were flavonoids, indicating that flavonoid metabolites play an important role in plant height morphogenesis in S. davidii [40]. Similarly, the dwarfing mechanism of Polish dwarf wheat was identified via transcriptome and proteome analysis; these results demonstrated that genes involved in flavonoid synthesis are up-regulated in dwarf wheat, resulting in increased flavonoid content, which limits IAA transport and leads to dwarfing [52]. Overall, these findings demonstrate that flavonoids regulate plant height development by controlling the polar transport of IAA. Similar to results of previous studies, our study revealed that the flavonoid biosynthesis was the most significantly enriched pathway in cassava; in particular, a series of genes involved in the flavonoid biosynthesis pathway (F3H, CHS, FLS, HCT, and F3′5′H) were found to be significantly up-regulated in dwarf cassava (Fig. 9). The up-regulated expression of these genes increases the flavonoid content of dwarf cassava, thereby inhibiting the polar transport of IAA; this inhibited IAA transport reduces cell ductility and ultimately limits cell expansion, resulting in dwarfing. However, the exact mechanisms involved in this process require further investigation.

Role of transcription factors in plant height development

Some transcription factors, such as GRF, PRE, bZIP, TCP, and WRKY, have been shown to play important roles in plant height development. In the present study, we identified a total of 92 transcription factors differentially expressed between XX048 and NZ199 cassava, of which bHLH and MYB genes were the predominant transcription factor families.

bHLHs regulate various growth and development processes, including seed germination, floral organ development, and stem and lateral root growth. In Arabidopsis spp., the mutation of bHLH transcription factor SPATULA has been found to limit cell proliferation and expansion, resulting in shortened internode length [53, 54]. Meanwhile, maize double mutants (bHLH12/14) exhibit early internode differentiation and abnormal branching, resulting in dwarfism and reduced stem veins [55]. In the present study, most bHLH transcription factors were significantly down-regulated in dwarf cassava; ultimately, this may have resulted in the inhibition of cell proliferation and expansion.

MYB transcription factors are a large class of plant transcription factors involved in physiological and biochemical processes, including root development, cell development, secondary cell wall synthesis, and seed development. In rice, overexpression of the MYB transcription factor OsMPH1 increases plant height by increasing the length of internode cells [56]. MYB transcription factors are also involved in plant hormone signal transduction. For example, MYB33 positively regulates GA signal transduction during plant growth [57]. bZIP is also involved in plant height growth and development. In Arabidopsis spp., the bZIP transcription factor influences plant height growth and development by controlling GA activity [58, 59]. In the present study, most MYB and bZIP transcription factors were significantly down-regulated in dwarf cassava. Ultimately, these genes may be involved in GA signal transduction in cassava, serving as positive regulators of plant height development. However, the regulatory effects of these transcription factors on cassava plant height development requires further verification.

Regulatory networks associated with cassava plant height development

Cassava plant height development is a complex regulatory process influenced by various factors. Based on the results of this study and previous research, we proposed a regulatory network model related to cassava plant height development (Fig. 13). Notably, increased cell wall lignification, reduced expression of plant hormone signal transduction–related genes, and excessive accumulation of flavonoid compounds are key factors involved in cassava plant height dwarfing. Furthermore, down-regulation of genes involved in cell wall synthesis and expansion, alongside downregulation of related transcriptional regulatory factors, may also influence height development in cassava. Specifically, the expression of genes related to the hormone signaling pathway and cell wall synthesis is down-regulated in dwarf cassava, resulting in the inhibition of cell wall biosynthesis and expansion. Additionally, some transcription factors influence cell development by participating in plant hormone signal transduction. The up-regulated expression of certain enzymes involved in lignin synthesis promotes cell wall lignification, thereby preventing cell elongation and growth and leading to dwarfing. Meanwhile, the up-regulated expression of flavonoid biosynthesis–related genes promotes the accumulation of flavonoids; these compounds inhibit IAA transport, which inhibits cell division and expansion and results in plant height dwarfing. In summary, the DEGs identified in this study play an important role in the regulatory network of cassava plant height development. Nonetheless, further functional identification studies are needed to elucidate the function of these candidate genes.

Fig. 13
figure 13

A hypothetical model illustrating the regulatory network related to plant height growth in cassava

Conclusion

This study analyzed transcriptome, physiological, and anatomical data to identify candidate genes related to cassava plant height development, ultimately proposing a genetic regulatory network model for cassava plant height development based on these results. Overall, we revealed the molecular mechanisms underlying the influence of plant hormone signaling, transcription factors, cell wall synthesis and elongation, flavonoid biosynthesis, and cell wall lignification on cassava plant height development. Therefore, our study provides a theoretical basis for understanding the molecular regulatory mechanisms of cassava plant height development; ultimately, these insights can be utilized to improve cassava crop traits via genetic breeding.

Materials and methods

Plant materials, growth conditions, and measurements of agronomic traits

In this study, the industrial cassava variants Xinxuan048 (XX048: high-stem) and Nanzhi199 (NZ199: dwarf-stem) were used. All plants were grown at the test base of the Guangxi South Asia Tropical Agricultural Science Research Institute, Longzhou County, Chongzuo City, China, under similar conditions. This test base, situated at an altitude of 138 m, experiences a southern subtropical monsoon climate, with an average annual temperature ≥ 23.1 °C, annual rainfall ≥ 1266 mm, annual sunshine of approximately 1521 h, an average annual relative humidity of approximately 88%, and an annual frost-free period of 361 days. The test base consisted of loam soil with medium fertility. The cassava were planted in April 2022 using a completely randomized block design, with each variety treated as a single treatment and repeated three times. The planting distance was 1 m between rows and 0.8 m between individual plants. Field management followed standard agricultural practices.

Two cassava lines, XX048 and NZ199, with uniform growth were selected to evaluate plant height during the stem growth period. The topmost internode tissue was collected for cytological observation [40] when the plant height difference between XX048 and NZ199 was the largest (September 2022, during the tuberous root expansion period). The microstructure of XX048 and NZ199 cassava was then compared, and the content of cellulose and lignin, the primary components of the cell wall was measured. Transcriptome sequencing was also conducted for further comparison.

Cytological observation

In September 2022, 5 plants with consistent growth were randomly selected from each planting plot. The topmost internodes tissues from XX048 and NZ199 were collected and named XX048-1, XX048-2, XX048-3, NZ199-1, NZ199-2 and NZ199-3. These samples were immediately fixed in formalin-aceto-alcohol for over 24 h. Before dehydration, plant lignocellulose was softened with an acidic softening solution for several days to a week. The samples were then dehydrated with gradient ethanol (75%, 85%, 90%, 95%, 100%, 100%) sequentially (4 h, 2 h, 2 h, 2 h, 2 h, 2 h, 2 h), followed by infiltration with alcohol benzene, transparent agent I, and transparent agent II for 30 min each, before performing paraffin embedding. Subsequently, 8 μm thick serial sections were cut using a 340E microtome (ThermoFisher Scientific, China) and fixed on slides. Finally, the sections were stained with toluidine blue, dehydrated, made transparent, sealed with neutral gum, and observed with an Olympus BX53F2 photomicroscope (Olympus, Tokyo, Japan).

Determination of lignin and cellulose content

In September 2022, 5 plants with consistent growth were randomly selected from each planting plot. The topmost internodes tissues from XX048 and NZ199 were then collected. The collected samples were dried at 80 °C to a constant weight, crushed, and sieved through 40 mesh. Subsequently. lignin and cellulose contents were determined using a kit from Suzhou Keming Biotechnology Co. Ltd., (China), following the manufacturer’s instructions.

Sample collection, RNA extraction, RNA sequencing, and data analysis

In September 2022, 5 plants with consistent growth were randomly selected from each planting plot. The topmost internodes tissues from XX048 and NZ199 were then collected and stored at − 80 °C for total RNA extraction.

Total RNA was extracted using the RNAprep Pure Plant Kit (Tiangen, Beijing, China), according to the manufacturer’s instructions. RNA concentration and purity were measured using a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Wilmington, DE). Meanwhile, RNA integrity was assessed using the RNA Nano 6000 Assay Kit with the Agilent Bioanalyzer 2100 system (Agilent Technologies, CA, USA).

For RNA sample preparations, 1 µg of RNA per sample was used as the input material. Sequencing libraries were then generated using the Hieff NGS Ultima Dual-mode mRNA Library Prep Kit for Illumina (Yeasen Biotechnology Co., Ltd., Shanghai, China), following the manufacturer’s recommendations; index codes were added to attribute sequences to each sample. Briefly, mRNA was purified from total RNA samples using poly-T oligo-attached magnetic beads. Subsequently, first-strand and second-strand cDNA synthesis were performed. Blunt ends were then formed using exonuclease and polymerase enzymes. After adenylation of 3′ ends, the NEBNext Adaptor with a hairpin loop structure was ligated to prepare for hybridization. The library fragments were purified using the AMPure XP system (Beckman Coulter, Beverly, USA). Then, 3 µl of USER Enzyme (NEB, USA) was used to treat size-selected, adaptor-ligated cDNA at 37℃ for 15 min, followed by 5 min at 95℃. Subsequently, PCR was performed with Phusion High-Fidelity DNA polymerase, universal PCR primers, and Index (X) Primer. PCR products were purified using the AMPure XP system and library quality was assessed with the Agilent Bioanalyzer 2100 system. Clean reads were obtained by removing reads containing adapters, ploy-N sequences, and low-quality reads from raw data. Q20, Q30, GC-content, and sequence duplication levels of the clean reads were calculated. All subsequent analyses were based on the high-quality clean data, including mapping to the cassava genome sequence. Gene expression levels were estimated as fragments per kilobase of transcript per million fragments mapped. DEGs were identified with an adjusted P-value < 0.01 and a log2 fold change > 1 using DESeq2. Finally, gene functions were annotated based on the following databases: NR (NCBI non-redundant protein sequences), Pfam (Protein family), KOG/COG (Clusters of Orthologous Groups of proteins), Swiss-Prot (A manually annotated and reviewed protein sequence database), KO (KEGG Ortholog database), and GO (Gene Ontology).

Fluorescence quantitative verification

Real-time quantitative polymerase chain reaction (qPCR) was performed as described by Cai et al. [60]. The primer sequences used in this study are shown in Supplementary Table S7. Relative expression levels were calculated using the 2−ΔΔCT method with actin as the endogenous reference gene.

Statistical analysis

Plant height data were collected and processed using Excel 2010 software. The statistical analysis of all differences between the means was performed by using ANONA and DUNCAN with the IBM SPSS statistics 20 program. Image-Pro Plus 6.0 software was used for data analysis of the slices, and Excel 2010 was used for mapping. All heatmapping was done using the OECloud tool on https://cloud.oebiotech.com.

Data availability

The raw sequencing data have been submitted to the NCBI SRA database (PRJNA976406).

Abbreviations

BR:

Brassinosteroid

GA:

Gibberellin

IAA:

Auxin

SL:

Strigolactone

SA:

Salicylic acid

JA:

Jasmonic acid

ABA:

Abscisic acid

XTH:

Xyloglucan endotransglucosylase/hydrolases

CesA:

Cellulose synthase

CSL:

Cellulose synthase-like

CYCD3:

Cyclin D3

SnRK2:

Serine/threonine-protein kinase SRK2

ETH:

Ethylene

CHS:

Chalcone synthase

HCT:

Shikimate O-hydroxycinnamoyltransferase

PGT1:

Phlorizin synthase

F3H:

Naringenin 3-dioxygenase

DFR:

Bifunctional dihydroflavonol 4-reductase/flavanone 4-reductase

CYP75A:

Flavonoid 3’,5’-hydroxylase

ANS:

Anthocyanidin synthase

ANR:

Anthocyanidin reductase

LAR:

Leucoanthocyanidin reductase

PAL:

Phenylalanine ammonia-lyase

4CL:

4-coumarate–CoA ligase

CCR:

Cinnamoyl-CoA reductase

CAD:

Cinnamyl-alcohol dehydrogenase

POD:

Peroxidase

HCT:

Shikimate O-hydroxycinnamoyltransferase

COMT:

Caffeic acid 3-O-methyltransferase/acetylserotonin O-methyltransferase

CCoAOMT:

Caffeoyl-CoA O-methyltransferase

CES:

Cellulose synthase

EXP:

Expansin

XYL:

Xylosidase

PL:

Pectin lyase

PM:

Pectin esterase

PA:

Pectin acetyl esterase

PG:

Polygalacturonase

XTH:

Xyloglucan endoglycosyltransferase/hydrolase

References

  1. Peng J, Richards DE, Hartley NM, Murphy GP, Devos KM, Flintham JE, et al. Green revolution’ genes encode mutant gibberellin response modulators. Nature. 1999;400(6741):256–61.

    Article  CAS  PubMed  Google Scholar 

  2. Han YJ, Kim YS, Hwang OJ, Roh J, Ganguly K, Kim SK, et al. Overexpression of Arabidopsis thaliana brassinosteroid-related acyltransferase 1 gene induces brassinosteroid-deficient phenotypes in creeping bentgrass. PLoS ONE. 2017;12(10):e0187378.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Kurotani KI, Hattori T, Takeda S. Overexpression of a CYP94 family gene CYP94C2b increases internode length and plant height in rice. Plant Signal Behav. 2015;10(7):e1046667.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Cheng Q, Dong L, Su T, Li T, Gan Z, Nan H, et al. CRISPR/Cas9-mediated targeted mutagenesis of GmLHY genes alters plant height and internode length in soybean. BMC Plant Biol. 2019;19(1):562.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Hedden P. The genes of the Green Revolution. Trends Genet. 2003;19(1):5–9.

    Article  CAS  PubMed  Google Scholar 

  6. Wang Y, Li J. Molecular basis of plant architecture. Annu Rev Plant Biol. 2008;59:253–79.

    Article  CAS  PubMed  Google Scholar 

  7. Clouse SD, Langford M, McMorris TC. A brassinosteroid-insensitive mutant in Arabidopsis thaliana exhibits multiple defects in growth and development. Plant Physiol. 1996;111(3):671–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Monna L, Kitazawa N, Yoshino R, Suzuki J, Masuda H, Maehara Y, et al. Positional cloning of rice semidwarfing gene, sd-1: rice green revolution gene encodes a mutant enzyme involved in gibberellin synthesis. DNA Res. 2002;9(1):11–7.

    Article  CAS  PubMed  Google Scholar 

  9. Sasaki A, Ashikari M, Ueguchi-Tanaka M, Itoh H, Nishimura A, Swapan D, et al. Green revolution: a mutant gibberellin-synthesis gene in rice. Nature. 2002;416(6882):701–2.

    Article  CAS  PubMed  Google Scholar 

  10. Spielmeyer W, Ellis MH, Chandler PM. Semidwarf (sd-1), green revolution rice, contains a defective gibberellin 20-oxidase gene. Proc Natl Acad Sci USA. 2002;99(13):9043–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Guo F, Huang Y, Qi P, Lian G, Hu X, Han N, et al. Functional analysis of auxin receptor OsTIR1/OsAFB family members in rice grain yield, tillering, plant height, root system, germination, and auxinic herbicide resistance. New Phytol. 2021;229(5):2676–92.

    Article  CAS  PubMed  Google Scholar 

  12. Multani DS, Briggs SP, Chamberlin MA, Blakeslee JJ, Murphy AS, Johal GS. Loss of an MDR transporter in compact stalks of maize br2 and sorghum dw3 mutants. Science. 2003;302(5642):81–4.

    Article  CAS  PubMed  Google Scholar 

  13. Jiang L, Liu X, Xiong G, Liu H, Chen F, Wang L, et al. DWARF 53 acts as a repressor of strigolactone signaling in rice. Nature. 2013;504(7480):401–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Guo F, Ma J, Hou L, Shi S, Sun J, Li G, et al. Transcriptome profiling provides insights into molecular mechanism in Peanut semi-dwarf mutant. BMC Genomics. 2020;21(1):211.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Arioli T, Peng L, Betzner AS, Burn J, Wittke W, Herth W, et al. Molecular analysis of cellulose biosynthesis in Arabidopsis. Science. 1998;279(5351):717–20.

    Article  CAS  PubMed  Google Scholar 

  16. Fagard M, Desnos T, Desprez T, Goubet F, Refregier G, Mouille G, et al. PROCUSTE1 encodes a cellulose synthase required for normal cell elongation specifically in roots and dark-grown hypocotyls of Arabidopsis. Plant Cell. 2000;12(12):2409–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Burton RA, Gibeaut DM, Bacic A, Findlay K, Roberts K, Hamilton A, et al. Virus-induced silencing of a plant cellulose synthase gene. Plant Cell. 2000;12(5):691–706.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Li M, Xiong G, Li R, Cui J, Tang D, Zhang B, et al. Rice cellulose synthase-like D4 is essential for normal cell-wall biosynthesis and plant growth. Plant J. 2009;60(6):1055–69.

    Article  CAS  PubMed  Google Scholar 

  19. Ding Z, Lin Z, Li Q, Wu H, Xiang C, Wang J. DNL1, encodes cellulose synthase-like D4, is a major QTL for plant height and leaf width in rice (Oryza sativa L). Biochem Biophys Res Commun. 2015;457(2):133–40.

    Article  CAS  PubMed  Google Scholar 

  20. Choi D, Lee Y, Cho HT, Kende H. Regulation of expansin gene expression affects growth and development in transgenic rice plants. Plant Cell. 2003;15(6):1386–98.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Cho HT, Cosgrove DJ. Altered expression of expansin modulates leaf growth and pedicel abscission in Arabidopsis thaliana. Proc Natl Acad Sci U S A. 2000;97(17):9783–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Zenoni S, Reale L, Tornielli GB, Lanfaloni L, Porceddu A, Ferrarini A, et al. Downregulation of the Petunia hybrida alpha-expansin gene PhEXP1 reduces the amount of crystalline cellulose in cell walls and leads to phenotypic changes in petal limbs. Plant Cell. 2004;16(2):295–308.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Goh HH, Sloan J, Dorca-Fornell C, Fleming A. Inducible repression of multiple expansin genes leads to growth suppression during leaf development. Plant Physiol. 2012;159(4):1759–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Bai MY, Fan M, Oh E, Wang ZY. A triple helix-loop-helix/basic helix-loop-helix cascade controls cell elongation downstream of multiple hormonal and environmental signaling pathways in Arabidopsis. Plant Cell. 2012;24(12):4917–a4929.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Bai MY, Shang JX, Oh E, Fan M, Bai Y, Zentella R, et al. Brassinosteroid, gibberellin, and phytochrome impinge on a common transcription module in Arabidopsis. Nat Cell Biol. 2012;14(8):810–b817.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Zheng K, Wang Y, Zhang N, Jia Q, Wang X, Hou C, et al. Involvement of PACLOBUTRAZOL RESISTANCE6/KIDARI, an atypical bHLH transcription factor, in auxin responses in Arabidopsis. Front Plant Sci. 2017;8:1813.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Jia D, Gong X, Li M, Li C, Sun T, Ma F. Overexpression of a Novel Apple NAC Transcription Factor Gene, MdNAC1, confers the dwarf phenotype in Transgenic Apple (Malus domestica). Genes (Basel). 2018;9(5):229.

    Article  PubMed  Google Scholar 

  28. Parmar A, Sturm B, Hensel O. Crops that feed the world: production and improvement of cassava for food, feed, and industrial uses. Food Security: Sci Sociol Econ Food Prod Access Food. 2017;9(3):1–21.

    Google Scholar 

  29. Ferraro V, Piccirillo C, Tomlins K, Pintado ME. Cassava (Manihot esculenta Crantz) and Yam (Dioscorea spp.) crops and their derived foodstuffs: Safety, Security and Nutritional Value. Crit Rev Food Sci Nutr. 2016;56(16):2714–27.

    Article  CAS  PubMed  Google Scholar 

  30. Ma N. Mapping of plant height genes and analysis of candidate genes in Brassica napus L. Northwest A&F University.2021.001715.

  31. Jia N, Liu J, Sun Y, Tan P, Cao H, Xie Y, et al. Citrus sinensis MYB transcription factors CsMYB330 and CsMYB308 regulate fruit juice sac lignification through fine-tuning expression of the Cs4CL1 gene. Plant Sci. 2018;277:334–43.

    Article  CAS  PubMed  Google Scholar 

  32. Wang YX, Teng RM, Wang WL, Wang Y, Shen W, Zhuang J. Identification of genes revealed differential expression profiles and lignin accumulation during leaf and stem development in tea plant (Camellia sinensis (L.) O. Kuntze). Protoplasma. 2019;256(2):359–70.

    Article  CAS  PubMed  Google Scholar 

  33. Ke SW, Luan X, Liang JY, Hung YH, Hsieh TF. Rice OsPEX1, an extensin-like protein, affects lignin biosynthesis and plant growth. Plant Mol Biol. 2019;100(1–2):151–61.

    Article  CAS  PubMed  Google Scholar 

  34. Wang G, Gao Y, Wang J, Yang L, Song R, Li X, et al. Overexpression of two cambium-abundant Chinese fir (Cunninghamia lanceolata) α-expansin genes ClEXPA1 and ClEXPA2 affect growth and development in transgenic tobacco and increase the amount of cellulose in stem cell walls. Plant Biotechnol J. 2011;9(4):486–502.

    Article  CAS  PubMed  Google Scholar 

  35. Rochange SF, Wenzel CL, McQueen-Mason SJ. Impaired growth in transgenic plants over-expressing an expansin isoform. Plant Mol Biol. 2001;46(5):581–9.

    Article  CAS  PubMed  Google Scholar 

  36. Jan A, Yang G, Nakamura H, Ichikawa H, Kitano H, Matsuoka M, et al. Characterization of a xyloglucan endotransglucosylase gene that is up-regulated by gibberellin in rice. Plant Physiol. 2004;136(3):3670–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Hyodo H, Yamakawa S, Takeda Y, Tsuduki M, Yokota A, Nishitani K, et al. Active gene expression of a xyloglucan endotransglucosylase/hydrolase gene, XTH9, in inflorescence apices is related to cell elongation in Arabidopsis thaliana. Plant Mol Biol. 2003;52(2):473–82.

    Article  CAS  PubMed  Google Scholar 

  38. Campanoni P, Nick P. Auxin-dependent cell division and cell elongation. 1-Naphthaleneacetic acid and 2,4-dichlorophenoxyacetic acid activate different pathways. Plant Physiol. 2005;137(3):939–48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Spartz AK, Lee SH, Wenger JP, Gonzalez N, Itoh H, Inzé D, et al. The SAUR19 subfamily of SMALL AUXIN UP RNA genes promote cell expansion. Plant J. 2012;70(6):978–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Zhao X, Sun XF, Zhao LL, Huang LJ, Wang PC. Morphological, transcriptomic and metabolomic analyses of Sophora davidii mutants for plant height. BMC Plant Biol. 2022;22(1):144.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Guo F, Hou L, Ma C, Li G, Lin R, Zhao Y, Wang X. Comparative transcriptome analysis of the peanut semi-dwarf mutant 1 reveals regulatory mechanism involved in plant height. Gene. 2021;791:145722.

    Article  CAS  PubMed  Google Scholar 

  42. Sun TP. Gibberellin-GID1-DELLA: a pivotal regulatory module for plant growth and development. Plant Physiol. 2010;154(2):567–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Ueguchi-Tanaka M, Ashikari M, Nakajima M, Itoh H, Katoh E, Kobayashi M, et al. GIBBERELLIN INSENSITIVE DWARF1 encodes a soluble receptor for gibberellin. Nature. 2005;437(7059):693–8.

    Article  CAS  PubMed  Google Scholar 

  44. Voegele A, Linkies A, Müller K, Leubner-Metzger G. Members of the gibberellin receptor gene family GID1 (GIBBERELLIN INSENSITIVE DWARF1) play distinct roles during Lepidium sativum and Arabidopsis thaliana seed germination. J Exp Bot. 2011;62(14):5131–47.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Wang T, Liu L, Wang X, Liang L, Yue J, Li L. Comparative analyses of anatomical structure, phytohormone levels, and Gene expression profiles reveal potential dwarfing mechanisms in Shengyin Bamboo (Phyllostachys edulis f. tubaeformis). Int J Mol Sci. 2018;19(6):1697.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Pan XF. Cloning and function analysis of genesrelated to gibberellin acid signaling in Phyllostachys edulis. Hangzhou: Zhe jiang A&F University; 2016.

    Google Scholar 

  47. Noguchi T, Fujioka S, Choe S, Takatsuto S, Yoshida S, Yuan H, et al. Brassinosteroid-insensitive dwarf mutants of Arabidopsis accumulate brassinosteroids. Plant Physiol. 1999;121(3):743–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Kim EJ, Russinova E. Brassinosteroid signalling. Curr Biol. 2020;30(7):R294–8.

    Article  CAS  PubMed  Google Scholar 

  49. Nakamura A, Fujioka S, Sunohara H, Kamiya N, Hong Z, Inukai Y et al. The role of OsBRI1 and its homologous genes, OsBRL1 and OsBRL3, in rice. Plant Physiol. 2006;140(2):580–590.

  50. Li J, Wen J, Lease KA, Doke JT, Tax FE, Walker JC. BAK1, an Arabidopsis LRR receptor-like protein kinase, interacts with BRI1 and modulates brassinosteroid signaling. Cell. 2002;110(2):213–22.

    Article  CAS  PubMed  Google Scholar 

  51. Jin YS, Seong ES, Qu GZ, Han W, Yoon BS, Wang MH. Cloning and morphological properties of Nicgl;CYCD3;1 gene in genetic tumors from interspecific hybrid of N. Langsdorffii and N. Glauca. J Plant Physiol. 2008;165(3):317–23.

    Article  CAS  PubMed  Google Scholar 

  52. Wang Y, Xiao X, Wang X, Zeng J, Kang H, Fan X, et al. RNA-Seq and iTRAQ reveal the Dwarfing mechanism of dwarf Polish wheat (Triticum polonicum L). Int J Biol Sci. 2016;12(6):653–66.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Makkena S, Lamb RS. The bHLH transcription factor SPATULA is a key regulator of organ size in Arabidopsis thaliana. Plant Signal Behav. 2013;8(5):e24140a.

    Article  Google Scholar 

  54. Makkena S, Lamb RS. The bHLH transcription factor SPATULA regulates root growth by controlling the size of the root meristem. BMC Plant Biol. 2013;13:1b.

    Article  Google Scholar 

  55. Tsuda K, Abraham-Juarez MJ, Maeno A, Dong Z, Aromdee D, Meeley R, et al. KNOTTED1 cofactors, BLH12 and BLH14, regulate Internode Patterning and Vein Anastomosis in Maize. Plant Cell. 2017;29(5):1105–18.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Zhang Y, Yu C, Lin J, Liu J, Liu B, Wang J, et al. OsMPH1 regulates plant height and improves grain yield in rice. PLoS ONE. 2017;12(7):e0180825.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Gocal GF, Sheldon CC, Gubler F, Moritz T, Bagnall DJ, MacMillan CP, et al. GAMYB-like genes, flowering, and gibberellin signaling in Arabidopsis. Plant Physiol Plant Physiol. 2001;127(4):1682–93.

    Article  CAS  PubMed  Google Scholar 

  58. Ali Z, Sarwat SS, Karim I, Faridi R, Jaskani MJ, Khan AA. Functions of plant’s bZIP transcription factors. Pak J Agr Sci. 2016;53(2):303–14.

    Google Scholar 

  59. akoby M, Weisshaar B, Dröge-Laser W, Vicente-Carbajosa J, Tiedemann J, Kroj T, et al. bZIP transcription factors in Arabidopsis. Trends Plant Sci. 2002;7(3):106–11.

    Article  Google Scholar 

  60. Cai Z, Cai Z, Huang J, Wang A, Ntambiyukuri A, Chen B, et al. 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(1):473.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We would like to thank Editage (www.editage.cn) for English language editing.

Funding

This work was supported by the National Key Research and Development Program of China (2023YFD1600601), the Guangxi Natural Science Foundation Program of China (2023GXNSFBA026150, 2024GXNSFBA010382), and the Basic Research Business Special Project of Guangxi Academy of Agricultural Sciences of China (Guinongke 2021YT157).

Author information

Authors and Affiliations

Authors

Contributions

ZQC: Writing Original Draft, Data analysis, Conceptualization; LXR, WLW, WH, HYY and HXC : Data analysis; ZHL, ZLH, RLH and XFZ: Prepared experimental material; CHZ, XL and TYL: manuscript review; HRL and LFH: Conceptualization, supervision, Writing- Review &Editing.  All authors reviewed the manuscrip.All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Longfei He or Hengrui Li.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Method declaration

All methods used in this study were implemented in accordance with relevant guidelines.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

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

Cai, Z., Ruan, L., Wei, W. et al. Morphological, anatomical, and transcriptomics analysis reveals the regulatory mechanisms of cassava plant height development. BMC Genomics 25, 699 (2024). https://doi.org/10.1186/s12864-024-10599-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-024-10599-2

Keywords