Skip to main content

Identification of BBX gene family and its function in the regulation of microtuber formation in yam


BBX proteins play important roles in all of the major light-regulated developmental processes. However, no systematic analysis of BBX gene family regarding the regulation of photoperiodic microtuber formation has been previously performed in yam. In this study, a systematic analysis on the BBX gene family was conducted in three yam species, with the results, indicating that this gene plays a role in regulating photoperiodic microtuber formation. These analyses included identification the BBX gene family in three yam species, their evolutionary relationships, conserved domains, motifs, gene structure, cis-acting elements, and expressional patterns. Based on these analyses, DoBBX2/DoCOL5 and DoBBX8/DoCOL8 showing the most opposite pattern of expression during microtuber formation were selected as candidate genes for further investigation. Gene expression analysis showed DoBBX2/DoCOL5 and DoBBX8/DoCOL8 were highest expressed in leaves and exhibited photoperiod responsive expression patterns. Besides, the overexpression of DoBBX2/DoCOL5 and DoBBX8/DoCOL8 in potato accelerated tuber formation under short-day (SD) conditions, whereas only the overexpression of DoBBX8/DoCOL8 enhanced the accelerating effect of dark conditions on tuber induction. Tuber number was increased in DoBBX8/DoCOL8 overexpressing plants under dark, as well as in DoBBX2/DoCOL5 overexpressing plants under SD. Overall, the data generated in this study may form the basis of future functional characterizations of BBX genes in yam, especially regarding their regulation of microtuber formation via the photoperiodic response pathway.

Peer Review reports


Yam (Dioscorea spp.) is an important tuberous crop belonging to the Dioscoreaceae family, grown in various regions globally [1, 2]. D. opposita ‘Tiegun’ is a popular Chinese yam cultivar that has been commercially grown in Jiaozuo city of Henan province, China for a long time due to its medicinal and nutritional benefits [3]. However, vegetative propagation of yam through its tubers year after year leads to the spread of viral infections and declining quality [4,5,6]. Researchers have successfully developed microtuber culture techniques, providing virus-free tubers suitable for transportation and storage [7].

Tuberization in yam [8] and other tuberous crops [9, 10] may involve hormomes (e.g., gibberellins (GAs) and abscisic acid (ABA)) and genes in GA and ABA synthesis, catabolism, and signalling [11,12,13]. Tuberization in model plant species, potato (Solanum tuberosum) [14] may also be regulated by photoperiodic genes such as circadian-regulated gene CONSTANS (CO) and the CO/FT module [15, 16]. In potato, StCO inhibited the transcription of the tuberigen StSP6A in leaves and repressed tuberization in a photoperiod-dependent manner. However, it is unknown whether tuberization in yam involves photoperiod-dependent genes, such as the BBX family—ZFP transcription factors.

The BBX family, also known as CONSTANS-LIKE (COL) proteins, are important in photomorphogenesis, seed germination, shade avoidance, photoperiodic regulation of flowering and tuberization [14, 17]. In Arabidopsis, the BBX gene family has 32 members divided into five structural groups based on their domains [18]. AtBBXs are mainly involved in the regulation of flowering. It has been reported that some BBX genes function to regulate flowering and are also involved in the regulation of tuberization. For example, AtBBX1/AtCO promotes flowering under LD conditions, but inhibites tuber development under SD conditions [19, 20]. AtBBX6/COL5 accelerates flowering [21], while over-expression of COL5 from lotus (Nelumbo nucifera) in potato increased tuber weight, but without changing the number of tubers under SD conditions [22]. Besides, AtBBX21 promoting seedling photomorphogenesis, also increased tuber yield in potato [23].

We obtained genomes and transcriptomes of the white Guinea yam (D. rotundata), the greater yam (D. alata), and D. opposita ‘Tiegun’, and analyzed the BBX gene family in the species with a focus on physical and chemical properties, conserved domains, gene structure, chromosome distribution, collinearity, and cis-acting elements. We identified two candidate DoBBX genes in ‘Tiegun’, and investigated their expression patterns in different tissues, under diurnal cycle and continuous light, and at five tuberization stages. The function of DoBBX2/DoCOL5 and DoBBX8/DoCOL8 in regulation of tuberization in potato was identified through heterologous overexpression. This study provides insights into the BBX gene family and the function of DoBBXs in microtuber formation.

Materials and methods

Genome and transcriptome resources

Chromosome-level reference genomes of the white Guinea yam (D. rotundata) [24], and the greater yam (D. alata) [25], A. thaliana version 10 (TAIR 10) [26] and Oryza sativa version 7.0 [27] were downloaded from,,, and, respectively. Our previous transcriptome data were deposited in NCBI Sequence Read Archive (SRA, with accession number SRP061414 [8].

Identification of the BBX gene family in yam

To identify the yam BBX genes, the B-box domain (Pfam00643/cl00034) was downloaded from the PFAM database ( [28] and was used as the query sequence to find the predicted proteins in yam using the HMMER 3.0 program with the threshold of E-value < 10–20 [29]. The predicted proteins containing the B-box domain were used to construct a yam-specific HMM file via hmmbuild from HMMER 3.0 [29]. The yam-specific BBX HMM was used as query against the predicted proteins of yam. The peptide sequences with the threshold of E-value < 10–10 and containing the B-box domain identified by the PFAM database [28] and NCBI-CDD tools ( [30] were selected as candidate proteins. The ClustalW program in MEGA 7.0 and chromosomal locations of candidate genes learned from GFF3 files were used to remove the repetitive sequences and the redundant alternatively spliced sequences [31]. The IDs, characteristics and sequences of BBX proteins from three yam species were listed in Supplementary Tables 13, and were named according to their relationships with homologous genes in A. thaliana.

Multiple sequence alignment and phylogenetic analyses of BBX gene family in yam

All the amino acid residues of DoBBXs, DrBBXs, DaBBXs, AtBBXs and OsBBXs were aligned by MEGA 7.0 [31], and the phylogenetic tree was constructed via the neighbor-joining (NJ) method with 1000 bootstrap replicates.

Chromosomal distribution and gene duplications

The lengths of the chromosomes and the physical locations of DrBBXs and DaBBXs were obtained from the genome annotation information (gff3) of D. rotundata and D. alata, respectively. MapGene2Chromosome v2 software ( was used to map the distribution of DrBBX and DaBBX genes [32]. MCScanX software was used to detect the genome replication gene pairs in D. alata and between different species [33]. Nucleotide sequences with alignment ratios and similarity ratios greater than 75% and with distances between genes on the same chromosome of less than 100 kb were selected as tandem duplications. Moreover, the genes located in the duplicated regions and nucleotide sequences with alignment ratios greater than 75% were selected as resulting from segmental duplications [34]. Tbtools was used to analyze collinearity of the BBX genes and to visualize the duplicated gene pairs [35, 36]. The KaKs_Calculator1.2 software was used to calculate the nonsynonymous (Ka), synonymous (Ks) substitution rates, and Ka/Ks values of the different gene duplication pairs [37]. The Ks values were used to estimate the approximate date of every duplicated event occurred in yam, using the formula: T = Ks/2λ × 10–6 Mya (λ = 6.5 × 10−9) (Supplementary Table 4).

Protein domains, motifs and gene structure analyses of DrBBXs and DaBBXs

The sequence logos of conserved domains of B-box1, B-box2 and CCT were generated by WebLogo (Web-based sequence logo generating application; [38, 39]. EvolView [40, 41] ( was used for phylogenetic tree visualization and conserved domain annotations. The conserved motifs of DrBBXs and DaBBXs proteins were identified by the MEME v4.9.0 ( with the optimum motifs set at ≥ 10 and ≤ 50 amino acids and the maximum number of motifs set at 10 [42]. Gene Structure Display Server ( was used to analyze the intron–exon distribution from the CDS and genomic sequence files of DrBBXs and DaBBXs [43].

Cis-acting elements in the promoter regions of DrBBXs and DaBBXs

We considered the sequence 2000 bp upstream of the initiation codons as the proximal promoter region sequences. Promoters were predicted in the PlantCARE database ( [44] and were categorized for functional groups and visualized using EvolView [40, 41].

Excavation of the differentially expressed BBX genes during microtuber formation in D. opposita ‘Tiegun’

The transcriptome sequencing datasets during tuberization of a traditional Chinese medicinal plant D. opposita ‘Tiegun’ were retrieved form the NCBI Sequence Read Archive (SRA) repository ( The transcriptome data was processed following Kim et al. (2018) [45]. In brief, the SRA files were converted into Fastqa format using the SRAToolkit software [46]. The raw reads were then subjected to quality control using Trimmomatic and FastQC softwares [47, 48]. Subsequently, high-quality reads were aligned to the reference genome from D. alata using the HISAT2 (version: 2.0.5). The mapped reads of were assembled and quantified using the StringTie program [49]. The FPKM was obtained for further analysis. The expression patterns of DoBBX family members during tuberization were examined using the FPKM and the expression heatmap was generated using the pheatmap R package [50].

Plant materials and growth condition

D. opposita ‘Tiegun’ were cultured on the MS medium containing 3 g·L−1 agar and 30 g·L−1sucrose in Engineering Technology Research Center of Nursing and Utilization of Genuine Chinese Crude Drugs in Henan Province (Henan Normal University in Xinxiang, China). The growth conditions were kept at 23 ± 2 °C with 16 h light/8 h dark photoperiod under 38 μm·sec−1·m−2 light intensity, with medium change every four weeks. For microbuer induction, plants from four-week-old D. opposita ‘Tiegun’ were grown on MS medium containing and 60 g·L−1sucrose with shaking (120 r·m−1) at 23 ± 2 °C under dark.

A Chinese potato variety ‘E-potato-3’ (S. tuberosum ‘E3’) kindly provided by Huazhong Agricultural University was used for genetic transformation. The potato seedlings were cultured on the potato medium (P0: MS + 30 g·L−1 sucrose + 3 g·L−1 agar) under the following: the growth conditions 23 ± 2 °C with 16 h light/8 h dark photoperiod (long day condition, LD) under 38 μm·sec−1·m−2 light intensity, with medium change every four weeks. For tuberization, buds at the 2th ~ 4th nodes from the shoot apex from two-week-old ‘E3’ were grown on the P0 medium for two weeks, and then transferred to the potato tuber induction medium (PTI: MS + 6 g·L−1 sucrose + 3 g·L−1 agar). The culture was maintained for 60 days at 18 ~ 20℃ in 38 μm·sec−1·m−2 on a 8 h/16 h day/night cycle (short day condition, SD) for tuberization.

Isolation and qRT-qPCR analyses of DoBBX2 and DoBBX8

All fresh samples were immediately frozen in liquid nitrogen after collection and stored at -80℃ for the expression analysis of DoBBX2 and DoBBX8.

The total RNA was extracted from the leaves at the 3th ~ 4th nodes from the four-week-old D. opposita ‘Tiegun’ using the RNA purification kit (TaKaRa, China). DoBBX2 and DoBBX8 were amplified using primers designed in Primer Premier 5.0 (Premier, Canada). The full CDS sequences of DoBBX2 and DoBBX8 were cloned into PBI121 with the homologous recombination technology (Vazyme, Nanjing, China).

The tissue-specific expression patterns of DoBBX2 and DoBBX8 were examined from root (R), stem (S), leaf (L), main bud (MB) and accessory bud (AB) of four-week-old D. opposita ‘Tiegun’ as well as the MBs 0, 7, 14, 28, 35, and 42 d after microtuber induction. The 3rd leaf from the shoot apex of the 4-week old D. opposita ‘Tiegun’ plants under 16 h light/8 h dark was collected in the dark period every 3 h for the daily expression patterns. Then, the plants were transferred to continuous light for circadian expression analysis. The 3rd leaf was collected every 3 h for 60 h after the initial adaptational period (8 h) for continuous light. The experiments were repeated three times, each time with 10 planets.

cDNAs were synthesized using a HiScriptII 1st Strand cDNA Synthesis Kit (Vazyme, Nanjing, China) and transcript levels were quantified with qPCR run in the Roche Real-Time PCR Detection System with the AceQ qPCR SYBR Green Master Mix (Vazyme, Nanjing, China). All reactions were done in triplicate. The primers for gene isolation and qRT-PCR analysis were listed in Supplementary Table 5.

The SMART ( was employed to determine conserved domains of DoBBX2 and DoBBX8 proteins. The NJ tree was constructed using MEGA 7.0 software with 1,000 bootstrap replicates among full-length amino acid residues of BBX2/COL5 and BBX8/COL8 homologous genes in D. opposita ‘Tiegun’ (Do), D. alata (Da), D. rotundata (Dr), Phoenix dactylifera (Pf), Elaeis guineensis (Eg), Dendrobium officinale (Dof), Phalaenopsis equestris (Peq), O. sativa (Os), A. thaliana (At), S. tuberosum (St), N. nucifera (Nn), Vitis vinifera (Vv), Malus domestica (Md) and Populus trichocarpa (Ptr) (Supplementary Table 6).

Transformation in potato and tuberization

Slices of tuber from ‘E3’ potato were used for Agrobacterium tumefaciens-mediated transformation. The vector containing 35S::DoBBX2 and 35S::DoBBX8 were transformed into A. tumefaciens strain LBA4404 using electroporation, individually. The overexpression (OE) lines were selected on kanamycin (100 mg·L−1) MS medium, verified by semi-RT-PCR and qRT-PCR (specific primers were listed in Supplementary Table 5), and propagated by growing single nodes on P0 medium under LD. Two-week-old seedlings of OE lines and ‘E3’ were transferred to SD or dark for the observation of tuberization, and tuberization rate-related parameters were documented including tuberization time, number of tubers per plants, and tuber weight per tuber. The experiments were repeated three times and each time with 10 plants for each treatment, i.e., SD, or dark.

Statistical analyses

Statistical analysis of collected data was carried out by Excel and SPSS in this chapter. The Fisher LSD test (P < 0.05) were performed to test whether there were significant differences at P = 0.05 in MS Excel 2016 (Microsoft, and SPSS V22 (ANOVA) (IBM SPSS Statistics version 22, SPSS lnc, Chicago, IL, 2014). Boxplot statistics were computed with “paired comparison plot” of Origin 2023 (Origin software; Microcal Software, Northampton, MA).


Identification and characteristics of BBX genes in yam

Identification and phylogenetic analysis of BBX genes in yam

Sequence identity analyses using HMM identified 20, 20, and 19 BBX genes in D. rotundata, D. alata, and D. opposita, respectively (Supplementary Tables 13).

The phylogenetic analysis suggested that there are five groups of the BBX families in D. rotundata, D. alata, and D. opposita (Fig. 1, Supplementary Tables 13). Group IV had the most members (7 Dr, 7 Da and 7 Do), while Group V had the fewest members (1 Dr, 2 Da and 2 Do) (Fig. 1). In addition, the orthologous genes identified among three Dioscorea species are highly conserved.

Fig. 1
figure 1

Phylogenetic relationships of BBXs. The full-length amino acid sequences of BBX genes in A. thaliana, O. sativa, D. rotundata, D. alata, and the translated amino acid sequences of BBX unigenes in D. opposita ‘Tiegun’ transcriptome were used for NJ phylogeny reconstruction with 1000 bootstrap replicates. The green check, purple circle, red triangles, blue stars and green circles represent AtBBXs, OsBBXs, DrBBXs DaBBXs and DoBBXs, respectively. Groups I, II, III, IV, and V are shaded in red, blue, yellow, prey, and purple, respectively

Chromosomal location and gene duplication events of BBX in yam

In D. rotunda, 17 BBX genes are located on 11 of the 20 chromosomes, whereas the other three DrBBXs are found on scaffolds. All of the 20 BBX genes are all located on 12 of the 20 chromosomes in D. alata (Fig. 2a and Supplementary Figure 1). Intraspecific collinearity analysis in D. alata suggests that there are both segmental and tandem duplications in this species, producing four and one BBX gene pairs, respectively (Fig. 2a).

Fig. 2
figure 2

The expansion and evolution analysis of BBX genes in the greater yam. The grey lines and gene names indicate duplicated blocks and red lines and names indicate duplicated gene pairs of BBXs. For each species, chromosome is indicated by the first letter of species name and chromosome number. a Intraspecific colinearity analysis in DaBBXs showing physical positions of 20 DaBBXs gene mapped onto the greater yam chromosome. b Collinearity between D. rotunda and D. alata, D. alata and rice (O. sativa), and D. alata and A. thaliana

Homologous gene analysis between D. alata and D. rotundata showed that 23 Da-Dr orthologous gene pairs were located in collinear regions (Fig. 2b, Supplementary Table 4). While 23 homologous genes pairs are colinear, genes are slightly different in their locations on the chromosomes between D. alata and D. rotundata.

D. alata has more colinear BBX genes with rice (26 pairs) than with A.thaliana (7 pairs). The DaBBX genes showed one-to-one (6), two-to-one (17), and three-to-one (3) corresponding relationships with BBX genes in rice, while showed one-to-one (1), two-to-one (2), and one-to-two (4) corresponding relationships with BBX genes in A. thaliana. The results indicated that loss and expansion of BBX genes were discovered in D. alata when compared with O. sativa and A. thaliana.

Subsequently, the Ka, Ks, and Ka/Ks of the 34 homologous gene pairs in D. alata and D. rotundata were calculated to evaluate their molecular evolutionary rates (Supplementary Table 4). Previous studies have shown that yam possessed a genome-wide duplication event (the value of Ks is approximately 1.21) [25, 51]. The Ks values of DaBBX5-DrBBX6, DaBBX6-DrBBX5, DaBBX10-DrBBX12, DaBBX11-DrBBX11, DaBBX5-DaBBX6, DaBBX10-DaBBX11, DrBBX5-DrBBX6, and DrBBX11-DrBBX12 were 1.6581, 2.1292, 2.3224, 2.5170, 1.8292, 2.0216, 1.994, and 2.7183, respectively. This indicated that these gene pairs were derived from genome-wide duplication events shared by D. alata and D. rotundata. The Ks values for the other paralogs, including DaBBX15-DaBBX16 and DrBBX18-DrBBX19, were 0.0190 and 0.0062, respectively, suggesting that they were derived from the ancient duplication events. The Ka/Ks values of all DrBBX/DaBBX gene pairs were less than 1.0 (Supplementary Table 4), indicating that they were under strong purifying selection during their evolution and a conserved evolutionary pattern was shared among BBX gene family in yam.

Conservative domain, motif and gene structure analyses of DrBBXs and DaBBXs

In the predicted protein domains, there is conserved B-box1 domain in all of DrBBX and DaBBX proteins, in the BBX gene groups I, II, III, IV and V, however, the conserved domains included the B-box1 + B-box2 + CCT, B-box1 + B-box2 + CCT, B-box1 + CCT, B-box1 + B-box2, and B-box1 domain combinations, respectively. Besides, all members in group I and 9 individual members in Group IV contained a valine-proline (VP) motif in their C-terminal region (Fig. 3a, Supplementary Figure 2). Protein sequence alignment and logo analysis show that B-box1, B-box2, CCT domains and VP motifs are highly conserved (Supplementary Figures 2 and 3).

Fig. 3
figure 3

Protein domain diagram, motif phases and gene structure of DrBBX and DaBBX genes mapped on phylogenetic trees. a The conserved domain distribution of DrBBXs and DaBBXs. The groups I, II, III, IV, and V are shaded in red, blue, yellow, grey, and purple, respectively. The red triangles and blue stars represent DaBBXs, and DaBBXs, respectively, and branch lengths are shown close to the branch nodes. b The conserved motifs on DrBBX and DaBBX proteins. c The exon and intron phases of DrBBX and DaBBX genes. The legend is shown in the upper right corner

DrBBXs and DaBBXs from same groups are similar in motif number, exon/intron number, length and arrangement (Fig. 3b,c). Motif 1 + Motif 10 are present in all of the DrBBX and DaBBX groups, Motif 2 is present in group I, II and III, Motif 3 + Motif 4 are present in group I, II and IV, whereas Motif 8 was specific in group I and Motif 6 + Motif 9 was specific in group II. Motif 5 was widely distributed in all of the DrBBX and DaBBX members (except DaBBX20), while Motif 7 was only distributed in DaBBX1, DaBBX7, DaBBX8, DaBBX14, DrBBX1, DrBBX8, DrBBX9, and DrBBX15 (Fig. 3b,c).

Cis-elements in the promoter regions of DrBBXs and DaBBXs

We identified 28 cis-elements. Besides the conventional cis-acting elements (TATA-box, CAAT-box) in the promoter, the other 26 cis-acting elements include 8 light responsive, 8 hormone-responsive, 8 stress responsive, and 2 growth and development groups (Fig. 4). The most frequent elements were Box 4 (the total number was 176) and G-Box (95) (involved in light responsiveness), ABRE (ABA responsive element, 98), and ARE (essential for the anaerobic induction, 92).

Fig. 4
figure 4

Cis-elements of BBX genes in D. rotundata and D. alata mapped on the phylogenetic tree of DrBBXs and DaBBXs, as in Fig. 3. The middle panel is the heat map showing the levels of cis-acting element enrichment. The chart on the right is a proportional map showing the abundance of cis-elements: the green rectangle, light responsive cis-elements; the blue rectangle, hormone responsive cis-elements; the orchid rectangle, stress responsive cis-elements, and the orange rectangle, growth and development related cis-elements

Excavation and functional analysis of BBXs in microtuber formation in D. opposita ‘Tiegun’

Excavation of differentially expressed BBXs during the microtuber formation

Based on the data obtained from RNA-seq during the microtuber formation in D. opposita ‘Tiegun’, the accumulation of DoBBXs during microtuber formation were depicted by heat map (Fig. 5). On the basis of the expression patterns of the BBX genes during the microtuber formation, the DoBBXs were clustered into three main classes: class A, class B and class C (Fig. 5). Six DoBBX members in class B were nearly absent during microtuber formation. Whereas 13 DoBBX genes showed differentially expression patterns during microtuber formation, which included Class A (DoBBX8, DoBBX14 and DoBBX4) were upregulated, and Class C (DoBBX2, DoBBX3, DoBBX7, DoBBX10, DoBBX12, DoBBX15, DoBBX16, and DoBBX18) were down regulated. Among them, DoBBX2 transcripts were the highest at EXP stage, and a significant reduction (~ 10.922-fold) was found from EXP to MTV stages. Besides, we observed the highest increase of DoBBX8 expression (~ 3.880-fold) from EXP to MTV stages. Further study would be focused on DoBBX2 and DoBBX8.

Fig. 5
figure 5

Expression heatmap of DoBBXs genes during the microtuber formation in D. opposita ‘Tiegun’ of different developmental stages. EXP: initial explants (0 d), BUD: axillary bud proliferation after three weeks (21 d), and MTV: microtuber visible after four weeks (28 d)

Spatial expression patterns of DoBBX2 and DoBBX8

DoBBX2 and DoBBX8 were isolated from D. opposita ‘Tiegun’ and submitted to NCBI’s GenBank (DoBBX2/DoCOL5, NCBI Accession No: CP897490; DoBBX8/DoCOL8, NCBI Accession No: CP897491). The full-length sequence of DoBBX2 cDNA was 1332 base pair (bp), ORF was 1041 bp length and encoded a deduced protein of 347 amino acids (aa) residues, with two B-box domains at the N-terminus and one CCT domain at the C-terminus. The full-length cDNA sequence of DoBBX8 was 1811 bp with a 1281 bp of ORF, and the deduced amino acids were 427 aa with one B-box domains at the N-terminus and one CCT domain at the C-terminus (Supplementary Figure 4).

Similar expression patterns were observed for the two DoBBX genes in 5 tissues with the highest expression in leaf (L), moderate in root (R), stem (S) and accessory buds (AB), and the lowest in main buds (MB) (Supplementary Figure 5, Fig. 6a). Specifically, the expression levels of DoBBX2 and DoBBX8 in L was ~ 3.90 times and ~ 7.61 times that in MB, respectively. The transcript levels of DoBBX2 were significantly higher than that of DoBBX8. This finding implied DoBBX2 and DoBBX8 may participate in the regulation of microtuber formation by receiving signals in leaf and then travel to the tissues forming microtubers.

Fig. 6
figure 6

Expression patterns of DoBBX2 and DoBBX8 gene. a DoBBX2 and DoBBX8 mRNA levels relative to DoActin in root (R), stem (S), leaf (L), main bud (MB) and accessory bud (AB). b The expression patterns of DoBBX2 and DoBBX8 gene during the microtuber formation under different light conditions: Black box: dark. White box: light. Grey box: 8-h light adaptations. Red box: continuous light phase. c The diurnal and circadian regulation of DoBBX2 and DoBBX8 gene expression under a 14-h photoperiod and then continuous light. Error bars = SEM

In order to verify the involvement of DoBBX2 and DoBBX8 genes in microtuber formation, the expression levels of these two DoBBXs were detected at five stages during microtuber formation (Fig. 6b). The DoBBX2 expression gradually decreased. However, the expression level of the DoBBX8 increased from EXP stage (0 d) to MTV stage (28 d), then decreased from early to late expansive stages (35d to 42 d). It showed that DoBBX2 and DoBBX8 were involved in microtuber formation through different patterns in D. opposita ‘Tiegun’.

Given that numerous plant BBX proteins are implicated in multiple light-regulated growth and developmental processes, we then sought to study the daily oscillation and circadian expression pattern of DoBBX2 and DoBBX8 (Fig. 6c). The changes in the expression profiles of DoBBX2 and DoBBX8 were observed in the day/night cycle oscillate following about 24-h rhythm. The transcript level of DoBBX2 increased during the day and reached a maximum in the middle of the light phase, whereas that of DoBBX8 had two main peaks throughout the 24-h diurnal cycle, the first with maximal expression in the light phase and the second with maximal expression in the end of light phase (the point of 15 h under light) or in the beginning of dark phase (before -6 h under dark, not shown). In the 8 h light adaptations period, the transcript level of DoBBX2 decreased, while that of DoBBX8 decreased first and then increased. During the continuous 44 h light period, the expression profiles of DoBBX2 were in synchrony with those observed during the diurnal cycle, but the maximum values were decreased. The expression of DoBBX8 was altered by prolonging the total duration at the phase of diurnal cyclic. It suggested that the diurnal expression of DoBBX2 was not directly regulated by light, but by the circadian clock; while diurnal oscillation of DoBBX8 expression required a dark period.

Overexpression of DoBBX2 and DoBBX8 regulates tuberization through photoperiodic pathway in potato

We produced more than 10 transgenic potato lines with 35S::DoBBX2 and 35S:: DoBBX8 constructs as detected using semi-PCR and qRT-PCR (Supplementary Figure 6).

The tuberization in transgenic and control potato lines were investigated under both SD and dark conditions (Fig. 7). Almost all of the transgenic lines (except for DoBBX8-2) formed tuber earlier than 'E3' under SD (Fig. 7a,c,f). The OE lines of DoBBX2-1, DoBBX2-2, DoBBX2-3, DoBBX8-1 and DoBBX8-3 formed tubers after 21, 34, 34, 34, 34 days of being transferred to SD, respectively, while a minimum of 36 days was required for ‘E3’. Under SD for 60 days, we observed a higher percentage of tuberization in the OE lines of DoBBX2 and DoBBX8 than in ‘E3’ (Fig. 7e). In addition, tuber yield was greater in the OE lines of DoBBX2-1, DoBBX2-2 and DoBBX2-3 than in the ‘E3’, with higher average of tubers produced per plant but without changing weight per tuber (Supplementary Figure 7, and Fig. 7g,h). However, the OE of DoBBX8 produced fewer tubers per plant (Fig. 7g,h). It indicated that both DoBBX2 and DoBBX8 advanced tuberization under SD, while number of tubers per plant was higher in OE DoBBX2 lines compared to ‘E3’, but it was slightly reduced in OE DoBBX8 lines.

Fig. 7
figure 7

Tuberization in the control (‘E3’) and transgenic potato plants expressing DoBBX2 and DoBBX8 transgenes (OE DoBBX2 and OE DoBBX8) under the SD and dark conditions. a-d The seedlings and potato tubers of ‘E3’, DoBBX2 and DoBBX8 overexpression lines (DoBBX2-1, DoBBX2-2, DoBBX2-3, DoBBX8-1, DoBBX8-2 and DoBBX8-3) after 60 days of growth under the SD and Dark conditions, bar = 2 cm. e, g, h Boxplots of tuberization rate, the number of tubers per plant and the mass per tuber after 60 days of growth under the SD and Dark conditions. f The boxplots of tuberization time under the SD and Dark. The red line within the boxplot marks the median (50th percentile) while the black box within the box marks the mean

Under dark, the tuberization time of ‘E3’ and OE DoBBX8 lines was significantly earlier than that under SD (Supplementary Figure 7, Fig. 7b,d,f). Under dark for 60 days, the tuberization rate, mass per tuber and yield were all reduced in ‘E3’ and OE lines (Fig. 7e,g,h). The number of tubers per plant in OE DoBBX8-2 and DoBBX8-3 lines had no significant change, while it demonstrated an obvious reduction in the ‘E3’ and OE DoBBX2 lines than that under SD (Fig. 7g). This indicates that although both over-expression of DoBBX2 and DoBBX8 can accelerate tuber formation under dark, only the overexpression of DoBBX8 enhanced the promoting effect under dark on tuber induction.


The BBXs are zinc-finger transcription factors and play vital roles in plant growth, development and response to biotic and abiotic stresses [52,53,54]. BBX gene family has been identified from several plant species, such as Arabidopsis [18], rice [55], orchid [56], bamboo [57], tomato [58], potato [59], pepper [60], grapes [61], cottons [62], pear [63], and apple [64]. However, little information is known about the BBX gene family in yam, which is one of the important tuberous crops. In this study, we performed systematic genome-wide identification and analyses of the BBX gene family in three yam species. Based on the transcriptome data during microtuber formation in D. opposita ‘Tiegun’, we selected two candidate genes for further investigation. Furthermore, we focus on the expression patterns and the potential function of the two candidate genes in potato tuberization.

In the present species, segmental duplication counts for 20% (4 Dr-Dr paralogs) of DrBBX and 25% (5 Da-Da paralogs) of DaBBX (Supplementary Table 4 and Fig. 2). The percentages are smaller than those in rice (60%, 18 OsBBX) [55], Phyllostachys edulis (88.89%, 24 PeBBX) [57], Solanum lycopersicum (40%, 12 SlBBXs) [58], but greater than those in D. officinale (10.53%, 2 DofBBXs) and P. equestris (12.50%, 2 PeqBBXs) [56]. It may be the reason for the difference in BBXs number among these species.

We found five groups of the BBX genes and they varied in conserved domains, motifs and gene structures. Group IV with the largest number of the BBX genes in D. rotundata, D. alata, Arabidopsis and rice, it may have undergone a much greater gene expansion [17]. However, the reasons are unclear. The differences between groups might be linked with a wide functional diversity in BBX gene family [65].

The different cis-regulatory elements in the promoter regions may also be important for functional diversity [66]. The G-box and its variants in the promoter region of BBX gene are the binding sites of central regulators like HY5 and PIFs of photomorphogenesis [67]. The homologous genes of DoBBX2 and DoBBX8 in D. rotundata and D. alata possess the G-box as well as the ABA/GA-responsive cis-elements. 85% of DaBBXs (17) and 90% DrBBXs (18), including homologous genes of DoBBX2 and DoBBX8, have ABRE motifs, whereas 85% of DaBBXs (17) and 85% of DrBBXs (17) contain CGTCA-motif (Fig. 4). The BBX1/CO gene was initially identified in Arabidopsis as an important regulator of flowering in the photoperiodic pathway [19, 68, 69] and also play a role in regulating the tuber formation in potato [16, 53, 55]. Therefore, DoBBX2 and DoBBX8 may have the potential functions of photoperiodic and hormone-regulated tuberization.

Our spatio-temporal expression analyses indicate that DoBBX2 and DoBBX8 proteins may function as components of circadian clock signals during microtuber formation. This is consistent with studies in Arabidopsis and rice [59, 70]. However, specific functions regulated by the circadian clock in DoBBX2 and DoBBX8 have not yet been clarified.

The length of the dark period is critical for tuberization [15]. D. opposita ‘Tiegun’ MBs were previously induced to form microtuber under dark [8]. Tubers of potato plants overexpressing DoBBX8 (the OE DoBBX8 lines) form earlier under dark than under the SD (Fig. 7). Surprisingly, the early tuber formation in the OE DoBBX8 lines occurs in both SD and dark conditions. The overexpression of DoBBX8 also increases the number of tubers per potato plant under dark conditions. However, DoBBX2 is down regulated during the microtuber formation, accelerating tuberization in potato under SD, but delaying tuberization under dark. Besides, DoBBX2 increases the tuber yield per plant both under SD and dark (Supplementary Figure 7). Our results are similar to those of the overexpression of NnCOL5 in potato, which affects the expression levels of NnCOL8, and the related genes in CO-FT and GA signal pathways [22]. Therefore, one of the future studies may focus on elucidating the regulatory relationships between DoBBX2, DoBBX8, hormone-regulated and photoperiod related genes.


In this study, we conducted a systematic genome-wide analysis of the BBX gene family in three yam species and identified 20 DrBBX, 20 DaBBX and 19 DoBBX genes in D. rotundata, D. alata and D. opposita ‘Tiegun’ transcriptome, respectively. The BBX genes form five major groups that are characterized by duplications, conserved domains, motifs, gene structure and cis-elements. Our RNA-seq data suggest that DoBBX2 and DoBBX8 are functional BBX genes during D. opposita ‘Tiegun’ microtuber formation under SD and dark condition. The overexpression of DoBBX2 and DoBBX8 in potato accelerates tuberization under SD and increase the number of tubers per plant under both SD and dark. Overall, this study forms the basis for future functional characterizations of yam BBX gene family, especially regarding regulation of tuberization via the photoperiodic pathway.

Availability of data and materials

The data included in this article and the additional files are available.


  1. Asiedu R, Sartie A. Crops that feed the world 1. Yams Food Security. 2010;2:305–15.

    Article  Google Scholar 

  2. Epping J, Laibach N. An underutilized orphan tuber crop—Chinese yam : a review. Planta. 2020;252:58.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Li M, Li J, Liu W, et al. A protocol for in vitro production of microtubers in Chinese yam (Dioscorea opposita). Biosci Biotechnol Biochem. 2014;78:1005–9.

    Article  CAS  PubMed  Google Scholar 

  4. Yahyaoui E, Casamento D, D’Onghia AM, et al. Fig virus-free production and survival rate improvement using meristem tip culture technique associated with the encapsulation technology. Acta Hort. 2018;1224:81–90.

    Article  Google Scholar 

  5. Bettoni JC, Mathew L, Pathirana R, et al. Eradication of potato virus S, potato virus A, and potato virus M from infected in vitro-grown potato shoots using in vitro therapies. Front Plant Sci. 2022;13:878733.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Yahyaoui E, Marinoni DT, Botta R, et al. Is it possible to produce certified hazelnut plant material in sicily? Identification and recovery of nebrodi genetic resources, in vitro establishment, and innovative sanitation technique from apple mosaic virus. Front Plant Sci. 2021;12:78142.

    Article  Google Scholar 

  7. Li M, Li J, Wang Y, et al. A simple method for microtuber production in Dioscorea opposita using single nodal segments. Pak J Bot. 2015;47:665–8.

    Article  CAS  Google Scholar 

  8. Li J, Zhao X, Dong Y, et al. Transcriptome analysis reveals key pathways and hormone activities involved in early microtuber formation of Dioscorea opposita. Biomed Res Int. 2020;2020:8057929.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Chen P, Yang R, Bartels D, et al. Roles of abscisic acid and gibberellins in stem/root tuber development. Int J Mol Sci. 2022;23:1–15.

    Article  CAS  Google Scholar 

  10. Utsumi Y, Tanaka M, Utsumi C, et al. Integrative omics approaches revealed a crosstalk among phytohormones during tuberous root development in cassava. Plant Mol Biol. 2022;109:249–69.

    Article  CAS  PubMed  Google Scholar 

  11. Zhou Y, Li Y, Gong M, et al. Regulatory mechanism of GA3 on tuber growth by DELLA-dependent pathway in yam (Dioscorea opposita). Plant Mol Biol. 2021;106:433–48.

    Article  CAS  PubMed  Google Scholar 

  12. Wu ZG, Jiang W, Tao ZM, et al. Morphological and stage-specific transcriptome analyses reveal distinct regulatory programs underlying yam (Dioscorea alata L.) bulbil growth. J Exp Bot. 2020;71:1899–914.

    Article  CAS  PubMed  Google Scholar 

  13. Sharma S, Deswal R. Dioscorea alata tuber proteome analysis uncovers differentially regulated growth-associated pathways of tuber development. Plant Cell Physiol. 2021;62:191–204.

    Article  CAS  PubMed  Google Scholar 

  14. González-Schain ND, Díaz-Mendoza M, Zurczak M, et al. Potato CONSTANS is involved in photoperiodic tuberization in a graft-transmissible manner. Plant J. 2012;70:678–90.

    Article  CAS  PubMed  Google Scholar 

  15. Hannapel DJ, Sharma P, Lin T, et al. The multiple signals that control tuber formation. Plant Physiol. 2017;174:845–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Plantenga FDM, Heuvelink E, Rienstra JA, et al. Coincidence of potato CONSTANS (StCOL1) expression and light cannot explain night-break repression of tuberization. Physiol Plant. 2019;167:250–63.

    Article  CAS  PubMed  Google Scholar 

  17. Gangappa SN, Botto JF. The BBX family of plant transcription factors. Trends Plant Sci. 2014;19:460–70.

    Article  CAS  PubMed  Google Scholar 

  18. Khanna R, Kronmiller B, Maszle DR, et al. The Arabidopsis B-box zinc finger family. Plant Cell. 2009;21:3416–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Robson F, Costa MMR, Hepworth SR, et al. Functional importance of conserved domains in the flowering-time gene CONSTANS demonstrated by analysis of mutant alleles and transgenic plants. Plant J. 2001;28:619–31.

    Article  CAS  PubMed  Google Scholar 

  20. Martínez-García JF, Virgós-Soler A, Prat S. Control of photoperiod-regulated tuberization in potato by the Arabidopsis flowering-time gene CONSTANS. Proc Natl Acad Sci USA. 2002;99:15211–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Hassidim M, Harir Y, Yakir E, Kron I, Green RM. Over-expression of CONSTANS-LIKE 5 can induce flowering in short-day grown Arabidopsis. Planta. 2009;230:481–91.

    Article  CAS  PubMed  Google Scholar 

  22. Cao D, Lin Z, Huang L, Damaris RN, et al. A CONSTANS-LIKE gene of Nelumbo nucifera could promote potato tuberization. Planta. 2021;253:1–11.

    Article  CAS  Google Scholar 

  23. Crocco CD, Ocampo GG, Ploschuk EL, et al. Heterologous expression of AtBBX21 enhances the rate of photosynthesis and alleviates photoinhibition in Solanum tuberosum. Plant Physiol. 2018;177:369–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Sugihara Y, Darkwa K, Yaegashi H, et al. Genome analyses reveal the hybrid origin of the staple crop white Guinea yam ( Dioscorea rotundata ). Proc Natl Acad Sci USA. 2020;117:31987–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Bredeson JV, Lyons JB, Oniyinde IO, et al. Chromosome evolution and the genetic basis of agronomically important traits in greater yam. Nat Commun. 2022;13:1–16.

    Article  CAS  Google Scholar 

  26. Philippe L, Berardini TZ, Donghui L, et al. The Arabidopsis information resource (TAIR): improved gene annotation and new tools. Nucleic Acids Res. 2012;40:1202–10.

    Article  CAS  Google Scholar 

  27. Ouyang S, Wei Z, John H, et al. The TIGR rice genome annotation resource: improvements and new features. Nucleic Acids Res. 2007;suppl_1:883–7.

    Article  Google Scholar 

  28. Finn RD, Tate J, Mistry J, et al. The Pfam protein families database. Nucleic Acids Res. 2008;32:D138-141.

    Article  CAS  Google Scholar 

  29. Eddy SR. Multiple alignment using hidden Markov models. Proc Int Conf Intell Syst Mol Biol. 1995;3:114–20 PMID: 7584426.

    CAS  PubMed  Google Scholar 

  30. Aron MB, Derbyshire MK, Gonzales NR, et al. CDD: NCBI’s conserved domain database. Nucleic Acids Res. 2015;43:D222–6.

    Article  CAS  Google Scholar 

  31. Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Jiangtao C, Yingzhen K, Qian W, et al. MapGene2Chrom, a tool to draw gene physical map based on Perl and SVG languages. Hereditas. 2015;37:91–7.

    Article  PubMed  Google Scholar 

  33. Wang Y, Tang H, Debarry JD, et al. MCScanX : a toolkit for detection and evolutionary analysis of gene synteny and collinearity. 2012;40:e49.

    Article  CAS  Google Scholar 

  34. Cannon SB, Mitra A, Baumgarten A, et al. The roles of segmental and tandem gene duplication in the evolution of large gene families in Arabidopsis thaliana. BMC Plant Biol. 2004;4:10.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Chen C, Chen H, Zhang Y, et al. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020;13:1194–202.

    Article  CAS  PubMed  Google Scholar 

  36. Krzywinski M, Schein J, Birol I, et al. Circos: an information aesthetic for comparative genomics. Genome Res. 2009;19:1639–45.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Zhang Z, Li J, Zhao XQ, et al. KaKs_calculator: calculating Ka and Ks through model selection and model averaging. Genom Proteom Bioinform. 2006;4:259–63.

    Article  CAS  Google Scholar 

  38. Crooks G, Hon G, Chandonia J, et al. WebLogo: a sequence logo generator. Genome Res. 2004;14:1188–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Schneider TD, Stephens RM. Sequence logos: a new way to display consensus sequences. Nucleic Acids Res. 1990;18:6097–100.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. He Z, Zhang H, Gao S, et al. Evolview v2: an online visualization and management tool for customized and annotated phylogenetic trees. Nucleic Acids Res. 2016;44:W236–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Zhang H, Gao S, Lercher MJ, et al. EvolView, an online tool for visualizing, annotating and managing phylogenetic trees. Nucleic Acids Res. 2012;40:569–72.

    Article  CAS  Google Scholar 

  42. Bailey TL, Elkan C. Fitting a mixture model by expectation maximization to discover motifs in biopolymers. Proc Int Conf Intell Syst Mol Biol. 1994;2:28–36 PMID: 7584402.

    CAS  PubMed  Google Scholar 

  43. Hu B, Jin J, Guo AY, et al. GSDS 2.0: an upgraded gene feature visualization server. Bioinformatics. 2015;31:1296–7.

    Article  PubMed  Google Scholar 

  44. Lescot M, Déhais P, Thijs G, et al. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002;30:325–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Taemook K, David SH, Lothar H, et al. Octopus-toolkit: a workflow to automate mining of public epigenomic and transcriptomic next-generation sequencing data. Nucleic Acids Res. 2018;46:e53.

    Article  CAS  Google Scholar 

  46. Goldberg DH, Victor JD, Gardner EP, et al. Spike train analysis toolkit: enabling wider application of information-theoretic techniques to neurophysiology. Neuroinformatics. 2009;7:165–78.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Brown J, Pirrung M, McCue LA. FQC Dashboard: integrates FastQC results into a web-based, interactive, and extensible FASTQ quality control tool. Bioinformatics. 2017;33:3137–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Pertea M, Kim D, Pertea GM, et al. Transcript-level expression analysis of RNA-seq experiments with HISAT. StringTie and Ballgown Nat Protoc. 2016;11:1650–67.

    Article  CAS  PubMed  Google Scholar 

  50. Kolde R, Kolde MR. Package “pheatmap.” R Packag. 2018. p. 1.

  51. Li Y, Tan C, Li Z, et al. The genome of Dioscorea zingiberensis sheds light on the biosynthesis, origin and evolution of the medicinally important diosgenin saponins. Hortic Res. 2022;9:uhac165.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Gómez Ocampo G, Ploschuk EL, Mantese A, et al. BBX21 reduces ABA sensitivity, mesophyll conductance and chloroplast electron transport capacity to increase photosynthesis and water use efficient in potato plants cultivated under moderated drought. Plant J. 2021;108:1131–44.

    Article  CAS  PubMed  Google Scholar 

  53. Lyu G, Li D, Li S. Bioinformatics analysis of BBX family genes and its response to UV-B in Arabidopsis thaliana. Plant Signal Behav. 2020;15:1782647.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Steinbach Y. The Arabidopsis thaliana CONSTANS-LIKE 4 (COL4) – a modulator of flowering time. Front Plant Science. 2019;10:651.

    Article  Google Scholar 

  55. Huang J, Zhao X, Weng X, et al. The rice B-box zinc finger gene family: genomic identification, characterization, expression profiling and diurnal analysis. PLoS ONE. 2012;7:15–7.

    Article  CAS  Google Scholar 

  56. Cao Y, Meng D, Han Y, et al. Comparative analysis of B-BOX genes and their expression pattern analysis under various treatments in Dendrobium officinale. BMC Plant Biol. 2019;19:1–16.

    Article  CAS  Google Scholar 

  57. Ma R, Chen J, Huang B, et al. The BBX gene family in Moso bamboo (Phyllostachys edulis): identification, characterization and expression profiles. BMC Genomics. 2021;22:533.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Chu Z, Wang X, Li Y, et al. Genomic organization, phylogenetic and expression analysis of the B-Box gene family in tomato. Front Plant Science. 2016;7:1552.

    Article  Google Scholar 

  59. Talar U, Kieøbowicz-Matuk A, Czarnecka J, et al. Genome-wide survey of B-box proteins in potato (Solanum tuberosum)-Identification, characterization and expression patterns during diurnal cycle, etiolation and deetiolation. PLoS ONE. 2017;12:1–19.

    Article  CAS  Google Scholar 

  60. Ma J, Dai J, Liu X, et al. Genome-wide and expression analysis of B-box gene family in pepper. BMC Genomics. 2021;22:1–18.

    Article  CAS  Google Scholar 

  61. Wang L, Xue J, Dai W, et al. Genome-wide identification, phylogenetic analysis, and expression profiling of CONSTANS-like (COL) genes in Vitis vinifera. J Plant Growth Regul. 2019;38:631–43.

    Article  CAS  Google Scholar 

  62. Feng Z, Li M, Li Y, et al. Comprehensive identification and expression analysis of B-Box genes in cotton. BMC Genomics. 2021;22:439.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Cao Y, Han Y, Meng D, et al. B-BOX genes: genome-wide identification, evolution and their contribution to pollen growth in pear (Pyrus bretschneideri Rehd.). BMC Plant Biol. 2017;17:156.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Shalmani A, Fan S, Jia P, et al. Genome identification of B-BOX gene family members in seven Rosaceae species and their expression analysis in response to flower induction in Malus domestica. Molecules. 2018;23:1763.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Crocco CD, Botto JF. BBX proteins in green plants: Insights into their evolution, structure, feature and functional diversification. Gene. 2013;531:44–52.

    Article  CAS  PubMed  Google Scholar 

  66. Haberer G, Hindemitt T, Meyers BC, et al. Transcriptional similarities, dissimilarities, and conservation of cis-elements in duplicated genes of Arabidopsis. Plant Physiol. 2004;136:3009–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. Lee J, He K, Stolc V, et al. Analysis of transcription factor HY5 genomic binding sites revealed its hierarchical role in light regulation of development. Plant Cell. 2007;19:731–49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Putterill J, Robson F, Lee K, et al. The CONSTANS gene of Arabidopsis promotes flowering and encodes a protein showing similarities to zinc finger transcription factors. Cell. 1995;80:847–57.

    Article  CAS  PubMed  Google Scholar 

  69. Griffiths S, Dunford RP, Coupland G, et al. The evolution of CONSTANS-like gene families in barley, rice, and Arabidopsis. Plant Physiol. 2003;131:1855–67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. An H, Roussot C, Suárez-López P, et al. CONSTANS acts in the phloem to regulate a systemic signal that induces photoperiodic flowering of Arabidopsis. Development. 2004;131:3615–26.

    Article  CAS  PubMed  Google Scholar 

Download references


We are very grateful to Jianhua Li (Hope College in USA) for reading the earlier version of the paper and providing helpful comments. We also very grateful to Conghua Xie (Huazhong Agricultural University) for providing the ‘E3’ potato and for his generous suggestions on the genetic transformation experiments.


This research was supported by China Agriculture Research System (CARS-21), The Zhongyuan high level talents special support plan-Science and Technology Innovation Leading Talents (No.224200510011), Special National Key Research and Development Plan (2021YFD1600100), National Natural Science Foundation of China (No. 81274019), Henan Normal University and Wen County People's Government Jointly Build National Yam Germplasm Resource Garden Project (No.5201049160163), Colleges and Universities in Henan Province Key Scientific Research Project Funding Scheme (No. 21A180013), and the Postdoctor Initiative Foundation of Henan Normal University (No. 5101049470215).

Author information

Authors and Affiliations



M.L. contributed to supervision, project administration, writing review and validation, and X.Z. contributed to the experimental design and supervision, writing review and editing. Q.Y. contributed to writing review and editing. J.Z., J.G. and T.W. gave substantial suggestions to the paper. Y.H., S.Z. and J.W. cultivated sterile tissue potato seedlings. S.L. and H.S. carried out mainly experiments. Y.C. finished the draft the manuscript and comprehensively analyzed data from all experimental results. All authors have read and approved the final manuscript.

Corresponding authors

Correspondence to Mingjun Li or Xiting Zhao.

Ethics declarations

Ethics approval and consent to participate

Not applicable. This study was conducted in accordance with the People’s Republic of China and international authorities’ relevant guidelines and legislation, including the official website of the Committee on Publication Ethics ( and the European Association of Science Editors (EASE) and other institutions’ publishing ethics standards.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1: Supplementary Table 1.

The information of BBX family genes in D. rotundata. Supplementary Table 2. The information of BBX family genes in D. alata. Supplementary Table 3. The information of BBX family genes in D. opposita ‘Tiegun’. Supplementary Table 4. Homologous BBX gene pairs between D. alata, D. rotundata, O. sativa and A. thaliana. The homologous gene pairs were identified by the results of BLAST and collinearity analysis. When both Ka and Ks were equal to 0, Ka/Ks was considered equal to 1. When only Ks was equal to 0, it was marked as Ka >> Ks. Segmental means that the gene might arise from segmental duplication. Tandem means that the homologous sequences in close genomic proximity. Supplementary Table 5. Primers used for gene cloning, semiqRT-PCR and qRT-PCR detection of DoBBX2 and DoBBX8. The double underscores represent the pBI121 plasmid based linker sequences. The restriction endonuclease recognition sites for Sac I and Xba I were included in the respective oligonucleotide primers, which are markered with single underlines. Supplementary Table 6. Accession numbers of the protein sequences used in phylogenetic analysis of DoBBX2 and DoBBX8. Do: D. opposita 'Tiegun', Da: the greater yam (D. alata), Dr: the white Guinea yam (D. rotundata), Pd: P. dactylifera, Eg: E. guineensis, Dof: D. officinale, Peq: P. equestris, Os: O. sativa, At: A. thaliana, St: S. tuberosum, Nn: N. nucifera, Vv: V. vinifera, Md: M. domestica, Ptr: P. trichocarpa.

Additional file 2: Supplementary Figure 1.

Chromosomal distribution of DrBBXs and DaBBXs. Chromosomal mapping was based on the physical position in the white D. rotundata (Dr) and D. alata (Da) chromosomes. The chromosome numbers are presented above each vertical bar. The scale on the left is in base pairs (Mb). Supplementary Figure 2. Multi-sequence alignment of BBX proteins in D. rotundata and D. alata. B-box domains are marked by green boxes, CCT conserved domain is markered by an orchid box, VP motifs are marked by blue boxes. Supplementary Figure 3. Conserved domains analyses of DrBBX and DaBBX proteins were aligned by a WebLogo program using default parameters. The B-box1, B-box2, and CCT conserved domain logos were obtained by aligning all the 40 BBXs sequences, 29 group I, II, IV sequences, and 23 group I, II, III sequences form the white Guinea yam and the greater yam, respectively. The group I’s VP motif logo was obtained by aligning the sequences of four group I DrBBXs and four group I DaBBXs. The group IV’s VP motif logo was obtained from DrBBX13, DrBBX14, DrBBX18, DrBBX19, DaBBX12, DaBBX13, DaBBX15, DaBBX15, and DaBBX16. Supplementary Figure 4. Cloning, domain diagram and phylogenetic analysis of DoBBX2 and DoBBX8. (a) The complete cDNA sequences and amino acid sequences of DoBBX2 and DoBBX8. (b) DoBBX2 and DoBBX8 protein domains. (c) Phylogenetic analysis of DoBBX2 and DoBBX8. Do: D. opposita ‘Tiegun’, Da: the greater yam (D. alata), Dr: the white Guinea yam (D. rotundata), Pd: Phoenix dactylifera, Eg: Elaeis guineensis, Dof: Dendrobium officinale, Peq: Phalaenopsis equestris, Os: O. sativa, At: A. thaliana, St: Solanum tuberosum, Nn: Nelumbo nucifera, Vv: Vitis vinifera, Md: Malus domestica, Ptr: Populus trichocarpa. The red diamond represents DoBBX2/DoCOL5 protein in ‘Tiegun’, and the blue triangle represents DoBBX8/DoCOL8 protein in ‘Tiegun’. Supplementary Figure 5. Morphological structure of 45 d-old D. opposita ‘Tiegun’ plant. (Bar=1 cm). Supplementary Figure 6. The transcript levels of DoBBX2 and DoBBX8 in 15-days-old control pants and overexpression plants. (a) The quantitative qRT-PCR results. (b) The semi quantitative PCR results. Left panel shows the transcript levels of DoBBX2 gene in 15-days-old control pants (‘E3’) and DoBBX2 overexpression plants (DoBBX2-3, DoBBX2-1 and DoBBX2-2). The grouping of gels/blots cropped from different parts of the same gel. Right panel shows the transcript levels of DoBBX8 gene in 15-days-old control pants (‘E3’) and DoBBX8 overexpression plants (DoBBX8-2, DoBBX8-3 and DoBBX8-1). The grouping of gels/blots cropped from different parts of the same gel. The potato ELONGATION FACTOR 1 α (EF1-α) gene was used as reference gene for primary data normalization. Error bars indicate standard deviations (SD) from three biological replicates. Supplementary Figure 7. Tuber yield per plant in controls and the lines overexpressing DoBBX2 and DoBBX8. Data are from three independent biological experiments (n = 3 independent biological experiments × 9 growth chambers × 10 individual plants). The red line within the boxplot marks the median (50th percentile) while the black box within the box marks the mean. Significant differences between two means are indicated by Fisher LSD test (P<0.05).

Additional file 3: Supplementary Info File 1.

The semi quantitative PCR result of OE DoBBX2 and DoBBX8 lines.

Additional file 4: Supplementary Info File 2.

The raw data of potato tuberization index.

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 The Creative Commons Public Domain Dedication waiver ( 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

Chang, Y., Sun, H., Liu, S. et al. Identification of BBX gene family and its function in the regulation of microtuber formation in yam. BMC Genomics 24, 354 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: