Comprehensive Identification and Expression Analysis of the Members of the B-Box Gene Family of Cotton Involved in Flower bud Differentiation and Responses to Multiple Stresses

doi:10.21203/rs.3.rs-196840/v1

Download PDF

Research Article

Comprehensive Identification and Expression Analysis of the Members of the B-Box Gene Family of Cotton Involved in Flower bud Differentiation and Responses to Multiple Stresses

https://doi.org/10.21203/rs.3.rs-196840/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Background: B-BOX (BBX) proteins are zinc-finger transcription factors with one or two BBX domains and sometimes a CCT domain. These proteins play an essential role in regulating plant growth and development, as well as in resisting abiotic stress. So far, the BBX gene family has been widely studied in other crops. However, no one has systematically studied the BBX gene in cotton.

Results: In the present study, 17, 18, 37 and 33 BBX genes were detected in Gossypium arboreum, G. raimondii, G. hirsutum and G. barbadense, respectively, via genome-wide identification. Phylogenetic analysis showed that all BBX genes were divided into 5 main categories. The protein motifs and exon/intron structures indicated that each group of BBX genes was highly conserved. Collinearity analysis revealed that the amplification of BBX gene family in Gossypium spp. was mainly through segmental replication. Nonsynonymous (Ka)/ synonymous (Ks) substitution ratios indicated that the BBX gene family had undergone purification selection throughout the long-term natural selection process. Moreover, transcriptomic data showed that some GhBBX genes were highly expressed in floral organs. Transcriptome data analysis and qRT-PCR verification showed that different GhBBX genes had different biological functions in flower bud differentiation, abiotic stress and stress response.

Conclusions: Our comprehensive analysis of BBX in G. hirsutum provides a basis for further study on the molecular role of GhBBXs in regulating flowering and cotton resistance to abiotic stress.

Epigenetics & Genomics

G. hirsutum

BBX

flower bud differentiation

phytohormone

stress response

Zinc-finger transcription factors are a kind of vital proteins, which play essential roles in plant growth and development, as well as in responsing to environmental stimuli [1, 2]. B-BOX(BBX) protein is one of zinc-finger transcription factors, which have attracted increased amounts of attention in recent years because of its various functions. BBX proteins are characterized by one or two conserved BBX domains at their N-terminus and sometimes a CCT domain at their C-terminus. BBX domains play an important role in transcriptional regulation and protein-protein interactions [3, 4]. The CCT domain is involved in transcriptional regulation and nuclear transport [5–7]. In Arabidopsis, 32 BBX proteins have been identified. According to the existence of BBX domain and CCT domain, these members can be divided into five subgroups. [3]. A growing body of evidence also shows that BBX proteins play a crucial role in flowering [8, 9], abiotic stress responses [10] and hormonal signaling networks[4].

CO/AtBBX1 was the first BBX gene studied in Arabidopsis; this gene controlled flowering time by regulating the expression of downstream Flowering Locus T (FT) gene [11–13]. Flowering is significantly delayed in CO mutant plants, while overexpression of CO could make plants flower early [14, 15]. Other BBX genes subsequently discovered, such as BBX4, BBX6, BBX7, and BBX32, can also regulate flowering time [16–19]. The mutant plants of BBX4 flower earlier than wild plants, indicating that BBX4 can delay the flowering of plants [16]. The mutant plants of BBX7 could also flower earlier, and the plants with overexpression of BBX7 showing the phenotype of delayed flowering [17]. The overexpression of BBX6 could make the plant flower earlier [19]. Overexpression of BBX32 showed the opposite phenotype, which could delay the flowering of the plant [18].

BBX proteins also participate in abiotic stress responses and hormone signalling networks. For example, compared with the wild type, the overexpression of BBX24 in Arabidopsis has higher salt tolerance, and the root length of BBX24 transgenic plants increased significantly under high-salinity conditions [20]. In Chrysanthemum, overexpression CmBBX24 not only prolonged flowering time, but also enhanced cold and drought resistance [13]. BBX genes also play a role in phytohormone signal transduction. AtBBX18 is a positive regulator of the gibberellin (GA) signalling pathway. Molecular and phenotypic studies have shown that BBX18 promotes hypocotyl growth by increasing bioactive gibberellin levels [21]. While BBX20 is a negative regulator of brassinolide signal pathway. It promotes hypocotyl growth by directly binding BZR1 and inhibiting its expression [22].

Cotton is an important cash crop species. Although members of the BBX family have been identified in Arabidopsis, tomato, pear, and apple [3, 23–25], no comprehensive study of BBX genes in cotton has been reported so far. With the release of the cotton genome [26], we can more systematically study the hypothetical functions of BBX genes in cotton. In the present study, we made a comprehensive analysis of the physical and chemical properties, chromosome distribution, collinearity, gene structures, cis-acting elements and expression patterns of the BBX gene family in G. arboreum, G. barbadense, G. hirsutum and G. raimondii. This research provides basic data for further study on the function of BBX genes in cotton.

Identification, chromosomal distribution and subcellular localization of BBX gene family

To identify the BBX genes in the Gossypium spp. genome and obtain their sequences, a global search of the Gossypium spp. Genomes were carried out by using HMM profiling of the BBX domain (PF00643). After ensuring that the identified members contain conserved domains and delete the repeated sequences, in total, of 17, 18, 37 and 33 putative BBX sequences were identified in G. arboreum, G. raimondii, G. hirsutum and G. barbadense, respectively, via genome-wide identification analysis. In G. hirsutum, 1 BBX was located on scaffold fragments. The BBXs were named according to their location on the chromosomes (Figure 1), and the BBXs located on the scaffold fragments in G. hirsutum is finally named. Table S2 contains detailed location information. The lengths of putative GaBBX protein sequences ranged from 163 aa (GaBBX3) to 374 aa (GaBBX13); GrBBXs, 197 aa (GrBBX16) to 374 aa (GrBBX10); GhBBXs, 166 aa (GhBBX29) to 374 aa (GhBBX32) and GbBBXs, 136 aa (GbBBX3) to 505 aa (GbBBX27). The predicted MW and pI of each BBX are shown in Table S2. The results of subcellular localization showed that most of the BBXs were located in the nucleus, while the others were located in the extracellular, indicating that the nucleus was the main region of biological functions of BBXs.

Based on the genomic location information of 105 BBX genes, we visualized the chromosome distribution of GaBBXs, GrBBXs, GhBBXs and GbBBXs (Figure 1). In G. arboreum, 17 GaBBXs were unevenly distributed on 10 chromosomes. A12 and A13 contained 3 GaBBXs, whereas the other 8 chromosomes, A01, A02, A04, A05, A06, A08, A09 and A11, contained 1 or 2 GaBBXs. In G. raimondii, 18 GrBBXs were located on 9 chromosomes. D02 and D08 contained the most GrBBXs (3), while the other 6 chromosomes contained only 1 or 2 GrBBXs. In G. hirsutum, 37 GhBBXs were unevenly mapped to 21 chromosomes, while, GhBBX37 was located on unassembled scaffolds. At13 contained 4 GhBBXs. At01, Dt01, Dt12 and Dt13 contained 3 GhBBXs, and the other 16 chromosomes contained only 1 or 2 GhBBXs. In G. barbadense, 33 GrBBXs were unevenly mapped to 20 chromosomes. At12, At13, Dt12 and Dt13 contained 3 GrBBXs. The other 16 chromosomes contained only 1 or 2 GrBBXs.

Phylogenetic analysis of the BBX gene family

To investigate the phylogenetic relationships of BBXs, 137 BBX protein sequences (G. arboreum (17), G. raimondii (18), G. hirsutum (37), G. barbadense (33) and A. thaliana (32)) were used to construct a phylogenetic tree based on the NJ method. Members of the BBX family were classified into 5 major groups, I-Ⅴ (Figure 2), and each subgroup was named according to the taxonomic results of previous studies in Arabidopsis [4]. It is worth noting that although AtBBX26 and AtBBX27 belong to group Ⅴ of Arabidopsis according to their structural classification, they are phylogenetically closer to AtBBX12 and AtBBX13, which are in group Ⅱ. As shown in Figure 2, group Ⅱ was the smallest subgroup, containing 7 BBXs. By contrast, group Ⅳ had the most massive numbers of BBX genes, including 69 BBXs. There were 29 BBXs in group I. No cotton species were divided into group Ⅲ or group Ⅴ. In G. hirsutum, BBX gene has 2, 9 and 26 members in group Ⅱ, Ⅰ and Ⅳ, respectively.

Replication events of BBX gene family

From the perspective of the cotton evolution, tetraploid cotton is the result of genome doubling of two diploid cotton hybrids. In terms of the number of genes, the sum of BBX genes in G. arboreum (17) and G. raimondii (19) is about equal to the number of those in G. hirsutum (37) or G. barbadense (33). The results further confirm this view. To explore the replication events of the BBX gene family, MCScanX was used to analyze the collinearity between the A_t and D_t subgenomes of G. hirsutum and their corresponding ancestral A and D diploid genomes (Figure 3). The data showed that most homologous gene pairs of the BBX gene family were amplified by segmental replication, which meant segmental replication played a key role in the evolution of the BBX gene family. However, the genomic evolution of allotetraploid cotton is extremely complex. In the process of evolution, the genome has experienced not only segmental duplication events but also many tandem duplication events. The duplicate types of BBXs in G. hirsutum are shown in detail in Table S3. In G. arboretum and G. raimondii, 1 and 2 tandem duplication events (GaBBX15/GaBBX16 as well as GrBBX1/GrBBX2 and GrBBX17/GrBBX18) were identified, respectively. In G. hirsutum, 4 tandem duplications (GhBBX1/GhBBX2, GhBBX15/GhBBX16, GhBBX18/GhBBX19, GhBBX34/GhBBX35) were discovered.

Ka/Ks ratios were calculated to evaluate the selection pressure of these homologous gene pairs. Among the 95 homologous gene pairs, 90 homologous gene pairs had Ka/Ks values <1, which indicates that most of the homologous gene pairs had undergone purifying selection in the process of evolution, and these genes pairs may play a similar function. Only a few homologous gene pairs have experienced positive selection, which may lead to new biological functions of these genes.

Analysis of gene structure and conservative motif

The results of the phylogenetic analysis showed that 37 GhBBXs could be divided into 5 groups (A-E), which contained 9, 2, 12, 5 and 9 members, respectively (Figure 5I). To better understand the structural characteristics of GhBBXs, the exon/intron structure was analyzed by GSDS (Figure 5III). The GhBBX genes contained 3 to 7 exons, but most of them contained less than 5 exons. Moreover, the conserved motif was further analyzed by MEME program. The GhBBXs in the same group shows similar motif composition, which further validated the classification results (Figure 5II). Except for group A, the order of motif 1 and motif 2 in GhBBX of other groups was the same. Motif 3 existed only in group A and group B, but motif 4 existed in all groups except group A and group B. Motif 5 exists only in group C, while motif 6 only exists in group E. Figure 5 shows that the distribution of conserved motif and exon/intron structure are different among different groups, but they are highly conservative on the same branches. The results showed apparent conservation, which laying a foundation for functional conservatism and providing guidance for follow-up functional research.

Analysis of cis-acting elements in GhBBX promoter regions

To better understand the regulation of GhBBXs gene transcription and expression, the promoter region of GhBBX (genomic DNA sequence 2 kb upstream of the transcription start site) were used to search the PlantCARE database. A variety of cis-elements were found in the GhBBX promoter region. Among the cis-acting elements, the cis-acting elements related to phytohormone and stress response are the focus of our attention. We found abscisic acid (ABA) response element, gibberellin (GA) response element, auxin (IAA) response element, salicylic acid (SA) response element and methyl jasmonate (MeJA) response element in 21, 19, 11, 17 and 17 GhBBX promoters, respectively. In some GhBBX promoters, there are cis-acting element related to multiple phytohormone, while in other GhBBX promoters, there are only cis-acting element related to a single phytohormone response. In terms of stress-related response elements, these cis-acting elements are mainly related to low temperature, drought, anaerobic and other defenses. In the midst of these elements, the anaerobic cis-acting element is the most frequent stress response element, which appears in the promoters of 32 GhBBX genes, followed by the cis-acting element in response to low temperature. It exists in the promoters of 20 GhBBX genes. Thus, it can be seen that GhBBX may respond to stress response and abiotic stress of cotton. In addition, a large number of light response elements were found in the promoter region of GhBBXs, including Box-4, G-box, GT1-motif, TCT-motif and MRE.

Expression patterns of GhBBXs in different tissues

In order to study the expression pattern of GhBBXs in different tissues, we analyzed the transcriptomic data of root, stem, leaf, anther, filament, pistil and petal in TM-1. The results showed that different members of the cotton BBX family showed different expression patterns. According to the expression characteristics and based on hierarchical clustering analysis, 37 GhBBXs were divided into 3 categories (I-III) (Figure. 5). 6 GhBBXs (GhBBX5, 8, 9, 2326, and 28) belonging to group II were highly expressed in nearly all tissues. 10 GhBBXs (GhBBX2, 4, 10, 16, 19, 21, 22, 29, 32 and 35) belonging to group I were poorly expressed in all tissues. The remaining members (GhBBX1, 3, 6, 7, 11, 12, 13, 14, 15, 17, 18, 20, 24, 25, 27, 30, 31, 33, 34, 36 and 37) belonging to group III exhibited slightly higher expression in vegetative organs, while others showed slightly higher expression in floral organs. These differences in expression patterns might be related to the various functions of GhBBXs.

Expression characterization of GhBBXs in cotton flower bud differentiation

Flower bud differentiation is an important sign that a plant is undergoing a transition from vegetative growth to reproductive growth [27]. To explore whether GhBBX gene, which is highly expressed in flower organs, is involved in the process of flower bud differentiation, we analysed the relative expression of these genes in the shoot apex of the early-maturing cotton cultivar CCRI50 and late-maturing cotton cultivar GX11 from one-leaf stage to five-leaf stage. The graphical representation of the expression profiles of 6 genes at 5 different times is shown in Figure 6. Among the 6 BBX genes, the expression levels of GhBBX5, GhBBX8, GhBBX9, and GhBBX26 in the early-maturing material were higher than those in the late-maturing material, and significant or extremely significant differences were detected at 4, 1, 3, and 5 periods, respectively. The expression levels of the other 2 genes, GhBBX23 and GhBBX28, in the early-maturing material were higher than those in the late-maturing materials at the most stages, but at three-leaf stage, the expression levels in the early-maturing material were lower than those in the late-maturing material.

Expression pattern of GhBBXs under multiple stress

The analysis of cis-acting elements in the promoter region showed that GhBBXs might be related to the abiotic stress, and some studies have shown that BBX genes play a positive role in abiotic stress resistance [23, 28]. Therefore, based on the transcriptomic data of TM-1, we further analyzed the expression patterns of GhBBX under PEG and NaCl stresses. The results shown that GhBBXs clustered mainly into three categories (I-III) based on a hierarchical clustering analysis according to expression features (Figure 5B). The GhBBX genes that belong to groups II and III responded to PEG and NaCl, especially at 12 h after treatment, showing significantly downregulated expression. 11 GhBBXs belonging to group I showed universally low expression after stress treatment.

Based on the transcriptomic data, GhBBX5, GhBBX8, GhBBX9, GhBBX17, GhBBX23, GhBBX26, GhBBX28 and GhBBX36 which belong to group II and whose expression levels increased after stress treatment, were selected for qRT-PCR verification. As shown in Figure 7, PEG stress affected the activity of GhBBXs. GhBBX5, GhBBX23 and GhBBX28 responded to PEG treatment at 1 h, after which their expression was downregulated continuously until a minimum point was reached at 12 h. Afterwards, their expression increased at 24 h. The expression of GhBBX8, GhBBX9, GhBBX17 and GhBBX26 increased at the beginning, after which it decreased but then increased again. The minimum expression level occurred at 12 h. GhBBX36 was the only gene whose expression abruptly decreased but then increased, followed by a sudden increase after reaching the minimum at 12 h. Under NaCl stress treatment, GhBBX5, GhBBX8, GhBBX9, GhBBX17, GhBBX23, GhBBX26 and GhBBX36 showed the same expression pattern. The expression of these genes quickly decreased but then increased. Afterwards, their expression reached a minimum at 12 h followed by an increase. GhBBX28 was the only gene whose expression began at 9 h, after which it decreased; after reaching a minimum at 12 h, the expression increased. Taken together, the results show that GhBBXs can respond to plant abiotic stresses, and at the same time, it provides potential candidate genes for further study.

Expression pattern of GhBBXs under phytohormone

In order to study the response of GhBBXs to phytohormone, 9 genes containing corresponding hormone-response elements in their promoters were selected for qRT-PCR experiment to further analyse the transcript levels of BBX genes under three different hormone treatments (Figure 8). After treated with ABA (Figure 8A), the expression of most of these genes increased early (mainly from 0.5 h to 2 h), followed by a decrease. The expression of three genes decreased early, after which it increased and finally decreased with time. Under GA treatment (Figure 8B), 4 expression patterns were found. the expression of most of these genes increased early, peaking at 0.5 h or 1 h after GA application, after which their expression decreased over time. The expression of two genes (GhBBX18 and GhBBX20) was downregulated at all time points, and the expression of one gene, GhBBX11, was upregulated at all time points except the 9 h time point. The expression of GhBBX7 was upregulated early on, after which it decreased and then increased. In response to IAA treatment (Figure 8C), these genes could be divided into 2 major patterns on the basis of their expression characteristics. The expression of most of these genes increased early, followed by a decrease, peaking from 0.5 h to 2 h. Three genes (GhBBX7, GhBBX13 and GhBBX25) presented downregulated expression at all time points, reaching a minimum at 12 h.

Expansion of BBX genes in cotton

Polyploidization is the primary mechanism of species formation and plant adaptation to the external environment [29]. The allotetraploid species G. hirsutum has undergone many gene deplication events, and the number of genes has increased throughout in the whole evolutionary process, which make it a perfect plant species to study the formation of polyploid in the whole evolutionary process [30, 31]. G. hirsutum is an allotetraploid AD genome composed of A genome and D genome. The diploid cotton of A genome is mainly distributed in Asia and Africa, and the cotton of D genome is mainly distributed in America. The AD genome is formed by 1.5 million years of natural hybridization and chromosome doubling of the A genome and the D genome [32]. In this study, in total, 17, 18, 37 and 33 putative BBX sequences were identified in G. arboreum, G. raimondii, G. hirsutum and G. barbadense, respectively, via genome-wide analysis. The sum of BBX genes in G. arboreum (17) and G. raimondii (19) is about equal to the number of those in G. hirsutum (37). Collinearity analysis showed that most homologous gene pairs of the BBX gene family were amplified by segmental replication, which meant segmental replication was the main way for the BBX gene family to evolve.

Conservation and evolutionary of BBX genes

The conservation of biological function was based on the conservation of structural sequences. We analyzed the conservative motifs and exon/intron structure of G. hirsutum, and classified it according to the conservative structure. As shown in figure 4, 9 GhBBX members were classified into class A, and two GhBBX members were classified into class B. All the members in these two categories contain two BBXs plus the CCT domain. The remaining 26 GhBBX members carry both B-BOX1 and B-BOX2, which could be divided into C, D and E class according to their gene structure and motif. The remaining 2 domains, B-BOX1 and BOX1 plus CCT, which appeared in Arabidopsis did not exist in cotton. To sum up, the gene composition of these GhBBX is similar to that of Arabidopsis and there are also differences. We speculate that the loss of the other two gene structures in upland cotton may be caused by evolution.

In the process of evolution, it must be under the pressure of natural selection. We calculated the Ka/Ks ratio between the At and Dt subgenomes of G. hirsutum and their corresponding ancestral A and D diploid genomes to evaluate the selection pressure of homologous genes pairs (Table S3). The results showed that the Ka/Ks values of most BBX homologues gene pairs were less than 1. Therefore, we speculated that most of the BBX genes had been purified selection throughout the long-term evolutionary, resulting in segmental replication, which might cause the gene to maintain its original function unchanged.

Cis-acting elements of BBXs in G. hirsutum

Promoters play an essential role in the regulation of plant gene expression. In this study, we extracted the promoter sequence of 2 kb upstream of the transcriptional initiation site of GhBBXs and statistically analyzed the cis-acting elements. As shown in Table S4, we identified a variety of stress-related cis-acting elements in the cotton BBX gene promoter, and we speculated that GhBBX genes might be related to a variety of external stresses. The result of the qRT-PCR experiment showed that the genes containing these cis-acting elements related to stress enabled plants to respond quickly and regulate the expression of related genes after PEG and salt treatment, so as to avoid further harm to plants. We speculate that they may enhance the resistance to environmental stress by enhancing their resistance to stress in G. hirsutum. At the same time, five hormone response elements were identified in cotton BBX gene promoter, and some genes containing cis-acting elements of ABA, GA and IAA were tested by qRT-PCR experiment. The results showed that the expression of genes containing hormone response elements changed in different degrees after dredging the corresponding hormone treatment. The promoters of some BBX genes contain more than one hormone-responsive cis-acting element. We speculate that these genes might participate in hormone signal transduction and act as transcriptional regulators. Their expression may be regulated by a variety of plant hormones.

Potential roles of GhBBX genes in G. hirsutum

An increasing number of studies have shown that BBX genes are involved in many aspects such as plant growth and development, as well as in resisting abiotic stress. TM-1 expression profiles were used to evaluate the expression patterns of 37 GhBBX genes in different cotton tissues and after external stress treatment. They are classified according to their expression abundance in the transcriptome, which may reflect genes in the same group may have similar functions. The results will provide strong data support for the study of specific genes.

BBX genes are involved in the flowering pathway. In Arabidopsis, BBX6/COL5 promotes flowering by increasing the expression of FT under short-day conditions [19], while BBX7/COL9 negatively regulates flowering by inhibiting the expression of CO and FT [17]. BBX32/EIP6 may inhibit flowering probably in a way independent of CO [33]. According to our findings, GhBBX8 and GhBBX26 are the closest homologues of AtBBX6 in cotton and the expression level in flower organs is significantly higher than that in other organs. Moreover, these genes clustered together with the other four genes, GhBBX5, GhBBX9, GhBBX23 and GhBBX28, in a heat map of tissue expression specificity. Therefore, we speculate that these genes may be related to flowering. Flowering is closely related to flower bud differentiation and early maturity of cotton. We used qRT-PCR experiment to analyse the shoot apex of the early-maturing material CCRI50 from one-leaf stage to five-leaf stage and the late-maturing material GX11 from one-leaf stage to five-leaf stage. These genes presented significant or extremely significant differences in expression between the early-maturing and late-maturing materials at most periods. Taken together, our results suggest that GhBBX genes are involved in the regulation of flowering.

BBX gene is also involved in stress signal transduction. In Arabidopsis, BBX24 is involved in signal transduction in response to salt stress. Compared with wild-type Arabidopsis plants, Arabidopsis plants overexpressing BBX24 are more salt-tolerant [20]. The results of qRT-PCR experiment showed that most BBX family genes began to be expressed at 0.5 h after abiotic stress. This rapid response mechanism can protect plants and reduce the damage caused by abiotic stress. Moreover, the BBX genes have a common feature: regardless of the changes in their expression at previous stages, the minimum expression of each occurred at 12 h. Taken together, these results indicate the GhBBX proteins play an important role in plant resistance to abiotic stress. However, the specific functions and regulatory mechanisms need to be clarified by further investigation.

In this study, members of the BBX gene family were identified in Gossypium arboreum, G. raimondii, G. hirsutum and G. barbadense, respectively. The main purpose of our study was to explore the phylogenetic relationships, collinearity, expansion patterns, gene structures, cis-acting elements and expression patterns of the BBX family. The transcription of GhBBXs in different tissues indicated that some members of GhBBXs might play a role in flower bud differentiation of cotton. Transcriptome and qRT-PCR analysis shown that GhBBX genes were involved in the complex regulatory network of the abiotic stress response. After hormone treatment, it is found that GhBBX gene may play an essential role in regulating plant growth and development by regulating a variety of hormone signal pathways. In general, genome-wide analysis of BBXs has laid a solid foundation for the functional analysis of BBX genes in cotton.

Identification of the BBX gene family

The genomes sequences of G. arboreum, G. barbadense, G. hirsutum and G. raimondii were obtained from the CottonFGD (http://www.cottonfgd.org/) [34]. The BBX protein sequences of Arabidopsis thaliana were obtained from the TAIR website (https://www.arabidopsis.org/) [35]. The hidden Markov model (HMM) configuration file of the BBX conservative domain (PF00643) was downloaded from the Pfam database (https://pfam.xfam.org/) [36]. HMMsearch program was used to search the BBXs in cotton by using the HMM file. The threshold of E-value was 1.0 E^-10[37]. The SMART database (http://smart.embl-heidelberg.de/) was subsequently used to identify whether the candidate protein sequences contain the special domain of BBX family [38]. The ExPASy website (http://web.expasy.org/protparam/) was used to predict the amino acid (aa) sequence, molecular weight (MW) and isoelectric point (pI) of the identified BBX proteins, and the subcellular localization was predicted through the CELLO version 2.5 server (http://cello.life.nctu.edu.tw/) [39, 40].

Phylogenetic analysis and sequence alignment

Multiple sequence alignment of BBX proteins from the four Gossypium species and Arabidopsis was performed by ClustalW. By using the MEGA 7.0 program and neighbour-joining (NJ) method, a rootless phylogenetic tree was constructed. Bootstrap replications (1000) were used to evaluate the reliability of the interior branches [41].

Chromosome location and gene duplication analysis

The physical chromosome locations of all BBX genes were obtained from the genome sequence databases of the four Gossypium species, and MapChart software was used for visualization [42]. The genome replication gene pairs of G. hirsutum, G. raimondii and G. arboreum were detected by MCScanX software [43]. Circos was used to map the segmental duplication events on chromosomes [44]. TBtools was used to calculate the non-synonymous mutation rate (Ka) and synonymous mutation rate (Ks) of BBX genes replication [45].

Conserved sequence and gene structure analysis

The exon-intron structure was analyzed on the GSDS (http://gsds.cbi.pku.edu.cn/) based on inputting GFF gene annotation files [46]. The motifs of BBX protein sequences were confirmed by the MEME (http://meme-suite.org/tools/meme), using the following parameters: up to six motifs, the best width of 6 to 50 [46].

Promoter region cis-acting elements analysis

The promoter sequence 2000 bp upstream of each GhBBX gene coding region was retrieved from CRIBI (http://genomes.cribi.unipd.it/). The predicted cis-acting elements were searched by using the PlantCARE online program (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/) [47].

Gene expression pattern analysis

The expression levels of GhBBXs in different tissues and after abiotic stress responses were obtained from previously reported transcriptome data of G. hirsutum TM-1, which were obtained from the NCBI Sequence Read Archive (SRA:PRJNA490626) (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA490626) [48]. TBtools was used to draw a heat map to show the expression patterns of GhBBX genes [45].

Plant Materials and Treatments

The early-maturing cotton cultivar CCRI50 and the late-maturing cotton cultivar Guoxinmian11(GX11) were grown in the field in Anyang, Henan Province, China. Shoot apex samples were taken at one, two, three, four and five-leaf stages. The samples were collected and frozen immediately in liquid nitrogen, after which they were stored at −80 °C.

The G. hirsutum cultivar Texas Marker-1 (TM-1) was used for stress and hormone treatments of seedlings in a climate-controlled greenhouse (light/dark cycle: 16 h at 28 °C/8 h at 22 °C). The seedlings were grown in hydroponic solution (pH 6.2) at the four-week-old of the three-leaf stage were treated with consisting of sodium chloride (300 mmol/L), polyethylene glycol (PEG) 6000 (30%), 5 µM abscisic acid (ABA), 5 µM gibberellic acid (GA) or 5 µM auxin (IAA), respectively. Control plants were only treated with hydroponic solution. Leaves were harvested at 0 h, 1 h, 3 h, 9 h, 12 h, 24 h after stress treatments, and leaves were harvested at 0 h, 0.5 h, 1 h, 2 h, 3 h, 6 h, 9 h, and 12 h after hormone treatments. The samples were collected and frozen immediately in liquid nitrogen, after which they were stored at −80 °C.

RNA extraction and qRT-PCR analysis

Total RNA was extracted from the collected samples using the RNA-prep Pure Plant Kit (TIANGEN, Beijing, China). Reverse transcription experiment was carried out using a PrimeScript RT Reagent Kit (Takara, Dalian, China). The gene-specific primers for qRT-PCR (Table S1) were designed with Oligo 7.0 software; Cotton histone 3 (GenBank accession No. AF094716) was used as an internal control. qRT-PCR (Promega, Madison, WI, USA) (three replications) was performed on an ABI 7500 real-time PCR system (Applied Biosystems, USA). The scheme was carried out as follows: step 1, 95 °C for 30 s; step 2, 40 cycles of 95 °C for 5 s followed by 60 °C for 34 s; and step 3, melting curve analysis. The relative expression levels of GhBBXs was calculated by 2^-ΔΔCT method, and statistical analyzed by T-tests.

Ethics approval and consent to participate

Not applicable

Consent for publication

Not applicable

Data availability

The data included in this article and the additional files are available. The transcriptome datasets of G. hirsutum TM-1 are under the accession number in PRJNA490626 NCBI.

Competing interests

The authors declare that they have no competing interests.

Funding

This work was supported by National Natural Science Foundation of China (32072112）Central Public-interest Scientific Institution Basal Research Fund (1610162020010307), Agricultural Science and Technology Innovation Program of Chinese Academy of Agricultural Sciences

Author contributions

HTW and ZF conceived and designed the study and prepared the manuscript. ZF, ML, YL, and XY performed the experiments. HW, LM, XF, and JL assisted with the analysis and interpretation of the data. SY participated in the design of the experiments and provided a critical review. All authors have read, edited, and approved the current version of the manuscript.

Riechmann JL, Heard J, Martin G, Reuber L, Jiang C, Keddie J, Adam L, Pineda O, Ratcliffe OJ, Samaha RR et al: Arabidopsis transcription factors: genome-wide comparative analysis among eukaryotes. Science 2000, 290(5499):2105-2110.
Zhao T, Wu T, Zhang J, Wang Z, Pei T, Yang H, Li J, Xu X: Genome-Wide Analyses of the Genetic Screening of C2H2-Type Zinc Finger Transcription Factors and Abiotic and Biotic Stress Responses in Tomato (Solanum lycopersicum) Based on RNA-Seq Data. Front Genet 2020, 11:540.
Khanna R, Kronmiller B, Maszle DR, Coupland G, Holm M, Mizuno T, Wu SH: The Arabidopsis B-box zinc finger family. Plant Cell 2009, 21(11):3416-3420.
Gangappa SN, Botto JF: The BBX family of plant transcription factors. Trends Plant Sci 2014, 19(7):460-470.
Griffiths S, Dunford RP, Coupland G, Laurie DA: The evolution of CONSTANS-like gene families in barley, rice, and Arabidopsis. Plant Physiol 2003, 131(4):1855-1867.
Yan H, Marquardt K, Indorf M, Jutt D, Kircher S, Neuhaus G, Rodriguez-Franco M: Nuclear localization and interaction with COP1 are required for STO/BBX24 function during photomorphogenesis. Plant Physiol 2011, 156(4):1772-1782.
Gendron JM, Pruneda-Paz JL, Doherty CJ, Gross AM, Kang SE, Kay SA: Arabidopsis circadian clock protein, TOC1, is a DNA-binding transcription factor. Proc Natl Acad Sci U S A 2012, 109(8):3167-3172.
Gonzalez-Schain ND, Diaz-Mendoza M, Zurczak M, Suarez-Lopez P: Potato CONSTANS is involved in photoperiodic tuberization in a graft-transmissible manner. Plant J 2012, 70(4):678-690.
Lee YS, Jeong DH, Lee DY, Yi J, Ryu CH, Kim SL, Jeong HJ, Choi SC, Jin P, Yang J et al: OsCOL4 is a constitutive flowering repressor upstream of Ehd1 and downstream of OsphyB. Plant J 2010, 63(1):18-30.
Ding L, Wang S, Song ZT, Jiang Y, Han JJ, Lu SJ, Li L, Liu JX: Two B-Box Domain Proteins, BBX18 and BBX23, Interact with ELF3 and Regulate Thermomorphogenesis in Arabidopsis. Cell Rep 2018, 25(7):1718-1728 e1714.
Putterill J, Robson F, Lee K, Simon R, Coupland G: The CONSTANS gene of Arabidopsis promotes flowering and encodes a protein showing similarities to zinc finger transcription factors. Cell 1995, 80(6):847-857.
Robson F, Costa MM, Hepworth SR, Vizir I, Pineiro M, Reeves PH, Putterill J, Coupland G: Functional importance of conserved domains in the flowering-time gene CONSTANS demonstrated by analysis of mutant alleles and transgenic plants. Plant J 2001, 28(6):619-631.
Yang Y, Ma C, Xu Y, Wei Q, Imtiaz M, Lan H, Gao S, Cheng L, Wang M, Fei Z et al: A Zinc Finger Protein Regulates Flowering Time and Abiotic Stress Tolerance in Chrysanthemum by Modulating Gibberellin Biosynthesis. Plant Cell 2014, 26(5):2038-2054.
Samach A, Onouchi H, Gold SE, Ditta GS, Schwarz-Sommer Z, Yanofsky MF, Coupland G: Distinct roles of CONSTANS target genes in reproductive development of Arabidopsis. Science 2000, 288(5471):1613-1616.
Valverde F, Mouradov A, Soppe W, Ravenscroft D, Samach A, Coupland G: Photoreceptor regulation of CONSTANS protein in photoperiodic flowering. Science 2004, 303(5660):1003-1006.
Datta S, Hettiarachchi GH, Deng XW, Holm M: Arabidopsis CONSTANS-LIKE3 is a positive regulator of red light signaling and root growth. Plant Cell 2006, 18(1):70-84.
Cheng XF, Wang ZY: Overexpression of COL9, a CONSTANS-LIKE gene, delays flowering by reducing expression of CO and FT in Arabidopsis thaliana. Plant J 2005, 43(5):758-768.
Tripathi P, Carvallo M, Hamilton EE, Preuss S, Kay SA: Arabidopsis B-BOX32 interacts with CONSTANS-LIKE3 to regulate flowering. Proc Natl Acad Sci U S A 2017, 114(1):172-177.
Hassidim M, Harir Y, Yakir E, Kron I, Green RMJP: Over-expression of CONSTANS-LIKE 5 can induce flowering in short-day grown Arabidopsis. 2009, 230(3):481-491.
Nagaoka S, Takano T: Salt tolerance-related protein STO binds to a Myb transcription factor homologue and confers salt tolerance in Arabidopsis. J Exp Bot 2003, 54(391):2231-2237.
Wang Q, Zeng J, Deng K, Tu X, Zhao X, Tang D, Liu X: DBB1a, involved in gibberellin homeostasis, functions as a negative regulator of blue light-mediated hypocotyl elongation in Arabidopsis. Planta 2011, 233(1):13-23.
Fan XY, Sun Y, Cao DM, Bai MY, Luo XM, Yang HJ, Wei CQ, Zhu SW, Sun Y, Chong K et al: BZS1, a B-box protein, promotes photomorphogenesis downstream of both brassinosteroid and light signaling pathways. Mol Plant 2012, 5(3):591-600.
Chu Z, Wang X, Li Y, Yu H, Li J, Lu Y, Li H, Ouyang B: Genomic Organization, Phylogenetic and Expression Analysis of the B-BOX Gene Family in Tomato. Front Plant Sci 2016, 7:1552.
Cao Y, Han Y, Meng D, Li D, Jiao C, Jin Q, Lin Y, Cai Y: B-BOX genes: genome-wide identification, evolution and their contribution to pollen growth in pear (Pyrus bretschneideri Rehd.). BMC Plant Biol 2017, 17(1):156.
Liu X, Li R, Dai Y, Chen X, Wang X: Genome-wide identification and expression analysis of the B-box gene family in the Apple (Malus domestica Borkh.) genome. Mol Genet Genomics 2018, 293(2):303-315.
Jaillon O, Aury JM, Noel B, Policriti A, Clepet C, Casagrande A, Choisne N, Aubourg S, Vitulo N, Jubin C et al: The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature 2007, 449(7161):463-467.
Yan B, Hou J, Cui J, He C, Li W, Chen X, Li M, Wang W: The Effects of Endogenous Hormones on the Flowering and Fruiting of Glycyrrhiza uralensis. Plants (Basel) 2019, 8(11).
Wei H, Wang P, Chen J, Li C, Wang Y, Yuan Y, Fang J, Leng X: Genome-wide identification and analysis of B-BOX gene family in grapevine reveal its potential functions in berry development. BMC Plant Biol 2020, 20(1):72.
Schemske RDWJARoE, Systematics: PATHWAYS, MECHANISMS, AND RATES OF POLYPLOID FORMATION IN FLOWERING PLANTS. 1998, 29:467-501.
Paterson AH, Wendel JF, Gundlach H, Guo H, Jenkins J, Jin D, Llewellyn D, Showmaker KC, Shu S, Udall J et al: Repeated polyploidization of Gossypium genomes and the evolution of spinnable cotton fibres. Nature 2012, 492(7429):423-427.
Zhu YX, Li FG: The Gossypium raimondii genome, a huge leap forward in cotton genomics. J Integr Plant Biol 2013, 55(7):570-571.
Li F, Fan G, Lu C, Xiao G, Zou C, Kohel RJ, Ma Z, Shang H, Ma X, Wu J et al: Genome sequence of cultivated Upland cotton (Gossypium hirsutum TM-1) provides insights into genome evolution. Nat Biotechnol 2015, 33(5):524-530.
Park HY, Lee SY, Seok HY, Kim SH, Sung ZR, Moon YH: EMF1 interacts with EIP1, EIP6 or EIP9 involved in the regulation of flowering time in Arabidopsis. Plant Cell Physiol 2011, 52(8):1376-1388.
Zhu T, Liang C, Meng Z, Sun G, Meng Z, Guo S, Zhang R: CottonFGD: an integrated functional genomics database for cotton. BMC Plant Biol 2017, 17(1):101.
Reiser L, Subramaniam S, Li D, Huala E: Using the Arabidopsis Information Resource (TAIR) to Find Information About Arabidopsis Genes. Curr Protoc Bioinformatics 2017, 60:1 11 11-11 11 45.
Finn RD, Coggill P, Eberhardt RY, Eddy SR, Mistry J, Mitchell AL, Potter SC, Punta M, Qureshi M, Sangrador-Vegas A et al: The Pfam protein families database: towards a more sustainable future. Nucleic Acids Res 2016, 44(D1):D279-285.
Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, Madden TL: BLAST+: architecture and applications. BMC Bioinformatics 2009, 10:421.
Letunic I, Doerks T, Bork P: SMART: recent updates, new developments and status in 2015. Nucleic Acids Res 2015, 43(Database issue):D257-260.
Artimo P, Jonnalagedda M, Arnold K, Baratin D, Csardi G, de Castro E, Duvaud S, Flegel V, Fortier A, Gasteiger E et al: ExPASy: SIB bioinformatics resource portal. Nucleic Acids Res 2012, 40(Web Server issue):W597-603.
Yu CS, Lin CJ, Hwang JK: Predicting subcellular localization of proteins for Gram-negative bacteria by support vector machines based on n-peptide compositions. Protein Sci 2004, 13(5):1402-1406.
Kumar S, Stecher G, Tamura K: MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol Biol Evol 2016, 33(7):1870-1874.
Voorrips RE: MapChart: software for the graphical presentation of linkage maps and QTLs. J Hered 2002, 93(1):77-78.
Wang Y, Tang H, Debarry JD, Tan X, Li J, Wang X, Lee TH, Jin H, Marler B, Guo H et al: MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res 2012, 40(7):e49.
Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA: Circos: an information aesthetic for comparative genomics. Genome Res 2009, 19(9):1639-1645.
Chen C, Chen H, Zhang Y, Thomas HR, Frank MH, He Y, Xia R: TBtools: An Integrative Toolkit Developed for Interactive Analyses of Big Biological Data. Mol Plant 2020, 13(8):1194-1202.
Hu B, Jin J, Guo AY, Zhang H, Luo J, Gao G: GSDS 2.0: an upgraded gene feature visualization server. Bioinformatics 2015, 31(8):1296-1297.
Lescot M, Dehais P, Thijs G, Marchal K, Moreau Y, Van de Peer Y, Rouze P, Rombauts S: PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res 2002, 30(1):325-327.
Zhang T, Hu Y, Jiang W, Fang L, Guan X, Chen J, Zhang J, Saski CA, Scheffler BE, Stelly DM et al: Sequencing of allotetraploid cotton (Gossypium hirsutum L. acc. TM-1) provides a resource for fiber improvement. Nat Biotechnol 2015, 33(5):531-537.

No competing interests reported.

TableS1qPCRprimer.xlsx
Additional file 1: Table S1. Primer pairs used in the qRT-PCR analysis.
TableS2basicinformation.xlsx
Additional file 2: Table S2. Detailed physicochemical characteristics and Gene location of BBX proteins of A. thaliana, G. arboreum, G. raimondii, G. hirsutum and G. barbadense. The pI and MW were computed by ExPASy website. Subcellular localizations of BBX proteins were predicted by CELLO version 2.5 server. Gene locus information for BBXs in G. arboreum, G. raimondii, G. hirsutum and G. barbadense. was obtained from gene annotation.
TableS3genepair.xlsx
Additional file 3: Table S3 | Ka/Ks ratios and duplicate types of gene pairs between the A and D subgenomes of allotetraploid cottons and their corresponding ancestral A and D diploid genomes, respectively. 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. WGD or segmental means that the gene might arise from whole-genome duplication or segmental duplication.
TableS4numberofciselements.xlsx
Additional file 4: Table S4 | Statistical results of phytohormone and stress response cis-acting elements in the promoter segments of BBXs. These cis-acting elements were identified using PlantCARE software with the upstream 2000-bp sequences of cotton BBXs.

Download PDF

Editorial decision: Major revision
05 Mar, 2021
Reviews received at journal
02 Mar, 2021
Reviewers agreed at journal
09 Feb, 2021
Reviewers agreed at journal
05 Feb, 2021
Reviewers invited by journal
04 Feb, 2021
Editor assigned by journal
04 Feb, 2021
Editor invited by journal
04 Feb, 2021
Submission checks completed at journal
04 Feb, 2021
First submitted to journal
01 Feb, 2021

You are reading this latest preprint version

Comprehensive Identification and Expression Analysis of the Members of the B-Box Gene Family of Cotton Involved in Flower bud Differentiation and Responses to Multiple Stresses

Status:

Version 1

Abstract

Figures

Background

Results

Identification, chromosomal distribution and subcellular localization of BBX gene family

Phylogenetic analysis of the BBX gene family

Replication events of BBX gene family

Analysis of gene structure and conservative motif

Analysis of cis-acting elements in GhBBX promoter regions

Expression patterns of GhBBXs in different tissues

Expression characterization of GhBBXs in cotton flower bud differentiation

Expression pattern of GhBBXs under multiple stress

Expression pattern of GhBBXs under phytohormone

Discussion

Expansion of BBX genes in cotton

Conservation and evolutionary of BBX genes

Cis-acting elements of BBXs in G. hirsutum

Potential roles of GhBBX genes in G. hirsutum

Conclusion

Materials And Methods

Identification of the BBX gene family

Phylogenetic analysis and sequence alignment

Chromosome location and gene duplication analysis

Conserved sequence and gene structure analysis

Promoter region cis-acting elements analysis

Gene expression pattern analysis

Plant Materials and Treatments

RNA extraction and qRT-PCR analysis

Declarations

Ethics approval and consent to participate

Consent for publication

Data availability

Competing interests

Funding

Author contributions

References

Additional Declarations

Supplementary Files

Status:

Version 1