Skip to main content

Genome-wide identification and analysis of WD40 proteins reveal that NtTTG1 enhances drought tolerance in tobacco (Nicotiana tabacum)

Abstract

Background

WD40 proteins, which are highly prevalent in eukaryotes, play important roles in plant development and stress responses. However, systematic identification and exploration of WD40 proteins in tobacco have not yet been conducted.

Results

In this study, a total of 399 WD40 regulatory genes were identified in common tobacco (Nicotiana tabacum). Gene structure and motif analysis revealed structural and functional diversity among different clades of tobacco WD40 regulatory genes. The expansion of tobacco WD40 regulatory genes was mainly driven by segmental duplication and purifying selection. A potential regulatory network of NtWD40s suggested that NtWD40s might be regulated by miRNAs and transcription factors in various biological processes. Expression pattern analysis via transcriptome analysis and qRT-PCR revealed that many NtWD40s exhibited tissue-specific expression patterns and might be involved in various biotic and abiotic stresses. Furthermore, we have validated the critical role of NtTTG1, which was located in the nuclei of trichome cells, in enhancing the drought tolerance of tobacco plants.

Conclusions

Our study provides comprehensive information to better understand the evolution of WD40 regulatory genes and their roles in different stress responses in tobacco.

Peer Review reports

Background

WD40 proteins, also known as WD40 domain-containing proteins, act as scaffolds for protein‒protein or protein‒DNA interactions, providing a platform for the recruitment of a variety of molecules to form functional complexes [1,2,3]. WD40 proteins are rarely found in prokaryotes but are highly prevalent in eukaryotes [4]. Each WD40 protein has several WD40 repeats, which contain 40–60 amino acid residue units as their main feature. Commonly, this unit contains glycine-histidine (GH) at the N-terminus and tryptophan-aspartate (Trp-Asp) at the C-terminus [5]. These conserved residues form a strong hydrogen bond network and stabilize the WD40 repeat fold [6, 7]. In addition, a single WD40 repeat contains a four-stranded antiparallel β–sheet [8]. Generally, WD40 protein motifs include 5–8 repeats, most often 7 repeats, which form a relatively stable β-propeller structure for various protein–protein interaction [9]. Therefore, investigating WD40 proteins is helpful for understanding their interactions with other proteins involved in various biological processes.

To date, WD40 proteins have been systematically explored in multiple plant species [10, 11], with 230 WD40s in Arabidopsis thaliana [11], 200 WD40s in rice (Oryza sativa) [12], 191 WD40s in cucumber (Cucumis sativus) [11], 315 WD40s in mango (Mangifera indica L.) [13], 579 WD40s in cotton (Gossypium hirsutum) [14], 743 WD40s in wheat (Triticum aestivum) [15], 231 WD40s in Cucurbita maxim [16], 177 WD40s in Sugar beet [17], and 178 WD40s in potato (Solanum tuberosum) [18]. WD40 proteins play important roles in diverse biological processes, such as cell cycle regulation [19, 20], apoptosis[21], autophagy [22, 23], gene transcription [24], signal transduction [25, 26], histone modification [27, 28], and chromatin assembly [29, 30]. A WD40-repeat protein from the Recretohalophyte Limonium bicolor could enhance trichome formation and salt tolerance in Arabidopsis [31]. Germinating modulator of rice pollen (GORI), encoding the WD40 domain protein, is required for pollen tube germination and elongation in rice [32]. The TRANSPARENT TESTA GLABRA 1 (TTG1) protein from mango enhances root growth and abiotic tolerance in Arabidopsis [13]. Liu et al. have revealed that WD40-REPEAT 5a plays a role in drought stress tolerance by regulating nitric oxide accumulation in Arabidopsis [33]. Tian et al. have reported that the full-length allele TaWD40-4B.1C, which possesses nonsense nucleotide variation, enhancing the drought tolerance and grain yield of wheat plants under drought stress. TaWD40-4B.1C could interact with canonical catalases, promote their oligomerization and activities, and reduce H2O2 levels under drought stress. The knockdown of catalase genes erases the role of TaWD40-4B.1C in drought tolerance. The introgression of TaWD40-4B.1C enhances the drought tolerance of the cultivar harboring TaWD40-4B.1 [34]. The well-known MYB-bHLH-WD40 (MBW) complex plays a vital role in regulating the anthocyanin biosynthesis pathway in plants. For example, in barley (Hordeum vulgare L.), Ant13 encodes a WD40-type regulatory protein that is needed for the transcriptional activation of a set of structural genes encoding enzymes involved in flavonoid biosynthesis at the leaf sheath base and in grains [35]. For this cascading MBW regulatory model in rice, Sun et al. has reported that S1 (bHLH) acts as the master gene by activating the expression of C1 (MYB), after which C1 activates the expression of WA1 (WD40) during anthocyanin biosynthesis [36]. Yue et al. have demonstrated that the MBW complex associated with flavonoid metabolism in strawberry fruits is FaMYB5/FaMYB10-FaEGL3 (bHLH)-FaLWD1/FaLWD1-like (WD40). Moreover, it directly targets the promoters of flavonoid 3′-hydroxylase (F3'H), leucoanthocyanidin reductase (LAR) and H+-ATP (autoinhibited H+-ATPase isoform 10, AHA10), thereby promoting the accumulation of flavonoids [37]. Given the important roles of WD40 proteins in various plant biological processes, additional studies need to be performed.

Cultivated tobacco (Nicotiana tabacum), an allotetraploid plant, is an important model species for plant genetic research. TTG proteins, which are WD40 proteins, are essential regulators of plant trichome and seed development [38,39,40]. In tobacco, NtTTG1 can regulate immune responses by recognizing a biotic elicitor and involving the transcription factor NPR1. In contrast, NtTTG2 is an immune repressor that particularly suppresses the NPR1 functional pathway to confer disease susceptibility. Transient expression of TTG1 from Raphanus sativus results in an increase in anthocyanin accumulation in tobacco leaves [41]. However, WD40 studies in tobacco plants are still scarce. Moreover, genome-wide investigations of WD40s in tobacco have not been reported until now. In this study, a comprehensive analysis of WD40 proteins in tobacco was conducted, including phylogenetic relationship analysis, duplication analysis, functional element analysis, and expression pattern analysis in different tissues and in response to different biotic and abiotic treatments. Furthermore, the significant role of one WD40 protein, NtTTG1, in the drought response was validated. Our study provides comprehensive insights into the biological functions of WD40 proteins in tobacco and unveils important roles of WD40 proteins under different stress treatments.

Results

Genome-wide identification of WD40 regulatory genes in tobacco

A total of 399 WD40 genes were identified based on comprehensive genome-wide analysis and were named NtWD40-001 to NtWD40-399 based on their location in the genome (Table S1). A total of 200 NtWD40 genes were unevenly distributed on 24 chromosomes, while the others were mapped onto scaffolds (Fig. 1). Chromosome 4 contained the highest number of WD40 genes (21), followed by chromosomes 17, 15, and 4 with 19, 15, and 15 genes, respectively. There were only 2 genes on chromosome 8. Subcellular localization prediction revealed that the tobacco WD40 regulatory genes were located mainly in the nucleus (216/399: 54.14%), cytoplasm (81/399: 20.3%) and chloroplast (70/399: 17.54%) (Table S1). Only a few proteins were located in the extracellular space (NtWD40-218 and NtWD40-339), vacuoles (NtWD40-251 and NtWD40-352) and peroxisomes (NtWD40-014 and NtWD40-151). The lengths and physicochemical properties varied greatly among the tobacco WD40 regulatory proteins. The lengths of the tobacco WD40 proteins ranged from 87 (NtWD40-050) to 2,378 (NtWD40-037) amino acids. The molecular weights (MWs) ranged from 9.871 (NtWD40-050) to 262.21 kDa (NtWD40-037), and the theoretical isoelectric points (pIs) ranged from 4.151 (NtWD40-067) to 10.427 (NtWD40-308) (Table S1).

Fig. 1
figure 1

Chromosome distribution for 200 NtWD40 genes. The scale bar on the left indicates the length (Mb) of the tobacco chromosome. The chromosome number is shown on the left side of each chromosome

Phylogenetic, gene structure and conserved motif analysis

To further analyze the evolutionary relationships of the WD40 regulatory genes in tobacco, 399 tobacco WD40 proteins and 230 Arabidopsis WD40 proteins were used to construct a phylogenetic tree. The tobacco WD40 regulatory genes could be divided into 6 major distinct clusters (Clusters A–F), containing 25, 30, 59, 42, 104 and 139 proteins, respectively (Fig. 2; Table S1). Genes with similar structures were distributed in the same cluster. Figure 2 showed that the structures of proteins from between Arabidopsis and tobacco were similar. Specifically, NtWD40-389, TTG member named NtTTG1, was located in the F cluster (Fig. 2).

Fig. 2
figure 2

Phylogenetic classification of NtWD40 proteins. Different colors in outer circle denote clusters AF. The blue color in inner circle represents tobacco and the green color in inner circle represents Arabidopsis

Interestingly, different NtWD40 clusters might vary significantly from each other (Fig. S1). For example, the protein lengths and molecular weights of the proteins in the F cluster were much greater than those of proteins in the other clusters. The WD40 genes in Cluster A trended to be more hydrophobicity than those in the other clusters. Clusters C and D had more nuclear‒localized members than those in the other clusters. To obtain further insights into WD40 evolution, we examined the exon‒intron structures of the NtWD40 genes (Table S1) and found that 27 NtWD40s had only one exon and no introns, 173 NtWD40s contained 2–9 exon, 171 NtWD40s contained 10–20 exons, and 28 NtWD40s contained more than 20 exons. Interestingly, the NtWD40 genes in Cluster B had more exons than those in the other clusters (Fig. S1). The 10 most conserved motifs for the tobacco WD40 regulatory genes were investigated using the MEME program and annotated using InterProScan (Table S2). Eight motifs (1, 2, 3, 4, 5, 6, 9 and 10) were annotated as WD40 repeats (Fig. S2A) and were present in most of the NtWD40s (99.7%, 99.7%, 90.2%, 94.7%, 78.4%, 80.2%, 77.7%, and 74.7%, respectively) (Fig. S2B), suggesting that these motifs had been preserved for a long period of time.

Gene duplication events and collinearity analysis of the WD40 regulatory genes in tobacco

Since gene duplication was an important mechanism for the evolution of novel gene functions, different duplication events of NtWD40 regulatory genes were investigated in this study. A total of 206 duplication events involving 342 paralogs were detected, including 41 segmental duplications, 4 tandem duplications and 161 uncertain duplication events (Fig. 3A). The ratios of nonsynonymous (Ka) to synonymous (Ks) substitution rates were calculated to explore the evolutionary trajectory of duplicated tobacco WD40 regulatory genes (Table S3). Ka/Ks = 1 was regarded as neutral selection, Ka/Ks < 1 indicated purifying selection, and Ka/Ks > 1 referred to positive selection. Among the 206 duplicated events, with the exception of 9 uncertain events, the Ka/Ks ratios of all the analyzed duplication events were less than 1 (Fig. 3B), indicating the predominance of strong purifying selection on these genes throughout their evolutionary history. This suggested that functional differentiation following gene duplication was limited by purifying selection.

Fig. 3
figure 3

Duplication and collinearity analysis of NtWD40 genes. A Tandem and segmental duplicated gene pairs are linked by the red in same chromosome and cyan lines between chromosomes, respectively. B Ka/Ks distribution for different duplication events. C Micro-collinearity relationship between WD40 regulatory genes in tobacco and those in Arabidopsis and tomato

To further explore the evolutionary patterns of WD40 regulatory genes among different plants, two representative species, Arabidopsis and tomato, were selected for comparative collinearity analysis with tobacco. A total of 58 pairs of genes were collinear between tobacco and Arabidopsis, and 125 pairs were observed between tobacco and tomato (Fig. 3C). The number of collinear gene pairs between tobacco and tomato was much greater than that between tobacco and Arabidopsis. The number of WD40 genes with collinearity between the different plants clearly varied widely, suggesting that the expansion of WD40 genes might have occurred during Solanaceae differentiation.

Cis-acting elements and regulatory networks associated with the NtWD40s

To further understand the functions of the NtWD40s, GO and KEGG analyses were performed. The GO results revealed that the NtWD40 regulatory genes might play important roles in the development of different tissues, including shoot apical meristem development; the positive regulation of developmental growth; mucilage extrusion from the seed coat; and various stress responses, including cellular responses to biotic stimuli and cellular responses to external biotic stimuli (Fig. 4A, 4B and 4C). Similar results were also observed in the KEGG pathway analysis, including circadian rhythm-plant nucleocytoplasmic transport, cell cycle, and ubiquitin-mediated proteolysis (Fig. 4D; Table S4).

Fig. 4
figure 4

GO and KEGG enrichment for NtWD40 genes. A Bar plot for top 20 enriched biological process (BP) terms. B Bar plot for top 20 enriched molecular function (MF) terms. C Bar plot for top 20 enriched cellular component (CC) terms. D Bar plot for top 10 enriched KEGG pathways

Considering that tobacco WD40s might participate in the development of different tissues and stress responses, we explored cis-acting elements (CREs) in the promoter regions of the NtWD40 genes to better understand the potential genetic regulatory mechanism involved. These identified CREs can be further divided into four types (Table S5): light-responsive (46.6%), phytohormone-responsive (28.3%), stress-responsive (9.3%) and development related elements (15.8%). The most abundant elements were light-responsive elements, including Box 4, G-Box, and the GT1 motif. ABRE, CGTCA-motif and TGACG-motif responsive elements accounted for the majority of phytohormone-responsive elements, with 657, 448, and 448, respectively. These large numbers of phytohormone-responsive elements indicated that hormone signals may regulate the expression of NtWD40 regulatory genes. Notably, all four CRE types were observed in the majority of WD40 genes (313/399:78.45%), indicating various regulatory roles for the NtWD40 genes (Fig. S3).

Due to the enrichment of numerous CREs in the promoter regions of the NtWD40 genes, we speculated that the corresponding transcription factors (TFs) may directly regulate the NtWD40 genes in tobacco. Hence, we investigated the relationship between TFs and NtWD40s using PlantTFDB. A total of 395 TF members from 39 families might play important roles in the regulation of NtWD40s (Fig. 5; Table S6). Among these TFs, ERF, MADS, Dof, MYB, C2HE, AP2, and bHLH were the most abundant. In addition, potential miRNA-binding sites for the NtWD40s were also explored using PsRNATarget. Overall, 75 miRNAs in 27 miRNA families might have potential regulatory relationships with NtWD40s (Fig. 5; Table S6). Most of the miRNAs had multiple target genes, of which the nta-miR159 and nta-miR162 families could target 5 and 3 NtWD40s, respectively. In contrast, several NtWD40s could be targeted by several miRNAs. For example, NtWD40-003 could be targeted by two miRNAs, nta-miR160 and nta-miR161, and NtWD40-012 could be targeted by three miRNAs, nta-miR408, nta-miR6146, and nta-miR6155. Hence, additional studies were needed to determine the regulatory relationship between NtWD40s and TFs/miRNAs.

Fig. 5
figure 5

Visualization of miRNAs and transcription factor (TF) that could target NtWD40 genes

Tissue expression analysis of NtWD40 regulatory genes

To further study the potential function of the NtWD40 genes, the expression patterns of eight representative tissues were explored. With the exception of 12 NtWD40s that were not expressed or were weakly expressed in the eight tissues, the remaining 387 genes were expressed in at least one tissue (Table S7). Interestingly, numerous NtWD40 regulatory genes, particularly NtWD40-035, NtWD40-088, NtWD40-298, and NtWD40-395, exhibited high expression in the roots, suggesting that these genes might play crucial roles in root development (Fig. 6A). Similarly, several NtWD40 genes, such as NtWD40-291 and NtWD40-375, were highly expressed specifically at flowering stage. In addition, several NtWD40 genes might be expressed in different tissues at different developmental stages. For example, some NtWD40 genes, such as NtWD40-008, NtWD40-161, NtWD40-275, and NtWD40-359, were unexpressed or expressed at low levels in young leaves and senescent flowers, but were expressed at higher levels in other developmental stages, revealing the differing functions of the NtWD40 genes during leaf and flower development.

Fig. 6
figure 6

Expression profile of NtWD40 gene in different tissues. A Heatmap of the expression of NtWD40s in different tissues using public RNA-seq data, scaled by column. B qRT-PCR analysis of nine selected NtWD40s in different tissues at three developmental stages. Data are presented as mean ± SD. C qRT-PCR analysis of nine selected NtWD40s in trichome tissues. Data are presented as mean ± SD

Nine NtWD40 genes (NtWD40-060, 074, 085, 135, 141, 161, 175, 336, and NtTTG1) with specific expression patterns in different tissues were randomly selected for further qRT‒PCR analysis, and the results were similar to those of the RNA‒seq analyses (Fig. 6B). NtWD40-060, NtTTG1, NtWD40-175, and NtWD40-141 were highly expressed at the flowering stage, especially in petals and stamens. In contrast, the expression of NtWD40-074 decreased significantly at the flowering stage. NtWD40-161 was highly expressed specifically in the seeds, NtWD40-336 was highly expressed in leaves at the flowering stage and in the seeds at the reproductive organ stage, and NtWD40-135 was highly expressed in the buds at the flowering stage. To further explore the expression patterns of the NtWD40s in trichomes, trichomes from the leaf, petiole and stem were explored among the different NtWD40s (Fig. 6C). The results showed that except for NtWD40-161, most of the NtWD40 genes examined, especially for NtTTG1, NtWD40-085, NtWD40-175, and NtWD40-060, were highly expressed in trichome tissues compared with expression in leaf and stem, indicating these NtWD40s might be involved in trichome development.

Expression analysis of NtWD40 genes in response to various abiotic and biotic stresses

To further explore the regulatory role of NtWD40 in various stress treatments, the expression patterns of the NtWD40s were studied using publicly available transcriptome data. As shown in Fig. 7A, the NtWD40s exhibited distinct response patterns to the different stress treatments. The expression of most of the NtWD40 genes did not significantly change in the leaves of tobacco plants subjected to cucumber mosaic virus (CMV), salt or cadmium treatment (Table S8). Notably, the response patterns of the NtWD40 regulatory genes differed between leaf and root tissues under similar stresses, such as high temperature, low temperature, and cadmium treatment. For example, the majority of NtWD40 genes were upregulated in leaves under low temperature, while the opposite patterns were observed in roots, in which most of the genes were downregulated. Similarly, the expression of numerous NtWD40 genes significantly decreased in cadmium-treated tobacco root tissue but exhibited no change or were upregulated in the leaves, suggesting a tissue-specific response of NtWD40s under different stress treatments. In addition, different NtWD40s might exhibit opposite response patterns to various stress treatments. For example, the expression of several genes (such as NtWD40-161 and NtWD40-175) was significantly upregulated after salt treatment, while that of other genes, such as NtWD40-049 and NtWD40-074, was significantly downregulated. Under drought treatment, most NtWD40 genes, such as NtWD40-175, NtWD40-336, NtWD40-074, NtWD40-085, and NtWD40-135, were significantly upregulated, while some, such as NtWD40-141, NtWD40-189, and NtWD40-254, were significantly downregulated. Interestingly, NtWD40-212, which was not expressed in any tissues of the wild type, was significantly upregulated in both leaves and roots under both high and low temperatures, indicating its specific role in the temperature response.

Fig. 7
figure 7

Expression profile of NtWD40 gene in various stress treatment. A Heatmap of the expression of NtWD40s under high temperature (HT), low temperature (LT), R. solanacearum (RS), black shank (BS), Cucumber Mosaic Virus (CMV), cadmium, salt, drought, Potato Virus Y (PVY) and topping treatment using public RNA-seq dataset. The expression change is indicated by the ratio of TPM value between the treatment and control (CK). B Fold change of the nine NtWD40s between expression of control and drought treatment. Note: Data are presented as mean ± SD, with one, two, and three asterisks denoting statistical significance at p < 0.05, p < 0.01, and p < 0.001, respectively

To compare the roles of the NtWD40s among the different stress treatments, nine genes with distinct response patterns were selected for further investigation of their expression patterns via qRT‒PCR. Under drought treatment, NtWD40-074, 135, 161, and 336 were significantly upregulated in both leaf and root tissues (Fig. 7B). Similarly, NtWD40-141 and NtWD40-189 exhibited opposite expression patterns between leaf and root tissue, indicating tissue-specific responses to drought stress. In contrast, NtWD40-389 was significantly downregulated in both leaf and root tissues. Similarly, NtWD40-074, 085, 135, 161, and 336 were also upregulated under salt treatment (Fig. S4). Additionally, we observed gradually increasing expression patterns for these four genes in both root and leaf tissues under salt treatment, with the highest expression occurring at 8 h. NtWD40-389 presented opposite expression patterns between the leaf and root tissues, indicating tissue-specific responses to salt stress. In response to cold treatment, the expression of NtWD40-135/189 in leaves and roots significantly increased after 39 h and decreased after 89 h, but the NtWD40-389 (NtTTG1) only shows this phenomenon in leaves. (Fig. S5A). In contrast, the NtWD40-060/085 ratio gradually decreased under cold treatment. Interestingly, compared with the expression of other hormones, the expression of NtWD40-389 (NtTTG1) was significantly induced by IAA treatment, indicating its potential function in plant growth (Fig. S5B). In contrast, the expression of NtWD40-074 significantly decreased under all the hormone treatments except IAA. NtWD40-175 exhibited different response patterns to different hormones. When tobacco plants were infected with powdery mildew, most genes, such as NtWD40-060, 074, 085, 135, 141, 175, and 389, were upregulated, although NtWD40-189 and NtWD40-336 were not (Fig. S5C). In summary, the distinct expression patterns for NtWD40 genes indicated that they might be involved in different abiotic and biotic stress responses.

Silencing of NtTTG1 enhanced drought tolerance of tobacco

We demonstrated that one of the TTG proteins, NtWD40-389 (named NtTTG1), was highly expressed in trichome tissue and exhibited significant downregulation under drought stress and upregulation under IAA treatment. To further explore the influence of NtWD40s on tobacco drought tolerance, we constructed VIGS-silenced plants harboring NtTTG1. Subcellular localization analysis indicated that NtTTG1 was preferentially localized in the nuclei of trichomes (Fig. 8A), which was consistent with the predicted results (Table S1) as well as findings in other species [42]. Compared with those of wild-type plants, the transcript abundance of NtTTG1 was significantly lower in positively silenced plants (Fig. 8B and 8C). After five days of drought treatment, the wild-type plants exhibited more severe wilting, but the middle and upper leaves of the NtTTG1-silenced plants had spread out (Fig. 8B), indicating that silencing NtTTG1 promoted the growth of the tobacco plants and enhanced their drought tolerance.

Fig. 8
figure 8

Subcellular localization of NtTTG1 and drought tolerance assessment in NtTTG1-silenced tobacco plants. A Subcellular localization of NtTTG1. The scale bars represent 20 μM. B Phenotype of NtTTG1-silenced tobacco plants after 5 days of drought treatment. C Relative expression of NtTTG1 detected by qRT-PCR before drought treatment. Note: Data are presented as mean ± SD, with one, two, and three asterisks denoting statistical significance at p < 0.05, p < 0.01, and p < 0.001, respectively

Discussion

WD40 proteins have been reported to play important roles in plant development and various physiological processes [31, 39]. However, studies related to the WD40 genes in the tobacco family are rare. With the progress in producing high-quality tobacco genome data, it has become possible to investigate tobacco WD40 proteins at the genome-wide level. This study was the first comprehensive and systematic report on the characteristics of the NtWD40 gene family in tobacco. A total of 399 NtWD40 proteins were identified and divided into 6 major distinct clusters (Clusters A–F), containing 25, 30, 59, 42, 104 and 139 proteins (Fig. 2). The number of NtWD40 genes in tobacco (399) was greater than that in most other species, such as the dicotyledonous species Arabidopsis (230), cucumber (191), Cucurbita maxima (231), sugar beet (177), and the monocotyledonous species rice (200). However, the number of DEGs was lower than that in cotton (579) and wheat (743), which implied considerable expansion of the WD40 gene family in tobacco, cotton and wheat compared with that in other plants. Common wheat, which has a relatively large genome size (17 GB), indeed contained the highest number of WD40 regulatory genes. However, for plant species whose genome size was smaller than that of tobacco, such as terrestrial cotton (2.5 GB), there was a greater abundance of WD40 regulatory genes. These results suggested that there might be no correlation between genome size and the number of WD40-regulated genes. In plants, polyploidy and gene duplication are often the main triggers for the expansion of gene families. Tobacco is a heterotetraploid plant that resulted from interspecific hybridization between N. tomentosiformis (2n = 24) and N. sylvestris (2n = 24) nearly 200,000 years ago. Similarly, common wheat and terrestrial cotton experienced genome-wide tripling duplication (WGT) and genome-wide duplication events, respectively [43, 44]. Commonly, duplication and homologous recombination are the main drivers of gene family evolution, and changes in family size and gene family distribution are associated with tandem or fragment duplication [11, 29]. Gene family expansion occurs mainly through three modes: fragment duplication of multiple genes, tandem duplication of single genes, and genome-wide duplication [45]. A total of 206 duplication events were found, including 41 segmental duplications, 4 tandem duplications and 161 uncertain duplication events (Table S3). Fragment duplication was considered the main driver of NtWD40 evolution, as 19.9% (41) of the duplication were fragment events (Table S3). Ka and Ks analyses of the duplicated gene pairs suggested that both were subject to purifying selection, which is consistent with findings in studies of wheat and Qingke [15, 46].

Cis-regulatory elements (CREs) are important for regulating gene activity in various biological processes. In the present study, by analyzing the promoter regions of NtWD40s in tobacco, numerous CREs related to the plant stress response, hormone response, and growth were identified. Similar results were also found for the promoter regions of WD40 regulatory genes in other plant species, such as Cucurbita maxima, Qingke, and cotton [14, 16, 46]. Based on these cis-elements, we found that 394 TF members from 39 families might play important roles in the regulation of NtWD40s. Many TF families implicated in stress responses, including WRKY, MYB, NAC, and bHLH, have been identified. The interactions between WD40 proteins and MYB and bHLH transcription factors have been widely investigated in flavonoid biosynthesis [42, 46]. In strawberry fruits, FaMYB5/FaMYB10-FaEGL3 (bHLH)-FaLWD1/FaLWD1-like (WD40) directly targeted the promoters of F3'H, LAR, and AHA10, thus committing to flavonoid accumulation. Similarly, in Raphanus sativus, RsTTG1, a WD40 protein, could interact with the bHLH transcription factor RsTT8 to regulate anthocyanin and proanthocyanidin biosynthesis [41]. Additionally, the TTG1-bHLH-MYB complex might be directly involved in the regulation of root growth and trichome patterning in Arabidopsis leaves [47]. Based on the results of cis-element analysis in our study, a WD40-bHLH-MYB interaction model was constructed. Finally, only two NtWD40 genes (NtWD40-060 and NtWD40-389, named NtTTG1) were identified in the WD40-bHLH-MYB model; these genes could also be regulated by three bHLH family members (Fig. S6). Moreover, several MYB and bHLH families could interact with each other through corresponding cis-responsive elements, in line with the findings of previous studies [48]. However, no interaction was detected between the WD40 and MYB families. Hence, we speculated that the bHLH families might serve as intermediate bridges for the WD40 and MYB families, but additional studies were needed to validate this hypothesis.

In addition, to elucidate the functions and regulatory mechanisms of miRNAs and WD40s, we predicted miRNAs that might target tobacco WD40 genes. Twenty-seven miRNA families, consisting of 75 miRNAs, may have regulatory relationships with NtWD40s (Fig. 5; Table S6). Previous studies have revealed that moderate miR156 transcript level were sufficient to enhance drought resilience in alfalfa by silencing SPL13 and increasing WD40-1 expression, whereas higher miR156 transcript levels resulted in drought susceptibility [49]. Similarly, our results also indicated that NtWD40-379, which was a guanine nucleotide-binding protein subunit beta-like protein (G-protein β), contained binding sites for nta-miR156. Previous studies have proved that plant G-protein β subunits could positively regulate drought tolerance by increasing the detoxification of ROS [50]. Hence, we speculated that NtWD40-379, regulated by nta-miR156, might also play important roles in the drought tolerance of tobacco plants.

The analysis of tissue expression patterns in tobacco showed that many WD40 genes were expressed in distinct tissue expression patterns. A large number of WD40 regulatory genes were highly expressed in flower tissue (Fig. 6A), suggesting their specific roles in flower development. Chen et al. reported that convergent selection of a WD40 protein enhanced grain yield in maize and rice [51]. The WD40 domain protein GORI was an integrative scaffold that was needed for pollen tube growth in rice [32]. Our qRT‒PCR results also indicated that many genes were highly expressed at the flowering stage, suggesting that these genes might be involved in flower growth and seed development. For example, NtWD40-175, a JINGUBANG protein, was a homolog of GORI and was highly expressed at the flowering stage in our study (Fig. 6B).

Moreover, by analyzing plant response patterns under different stress treatments, we found that the expression of some NtWD40s did not significantly change in response to stress (Fig. 7A), suggesting that these genes may respond to plant stress through other pathways. It was worth noting that the same gene type may exhibit varying stress response patterns. For instance, NtWD40-024 exhibited robust responses to the biotic stresses R. solanacearum (RS) and black shank (BS) but displayed insensitivity to CMV or PVY infection (Fig. 7A). However, further investigations were needed to determine the disease response mechanisms of WD40 regulatory genes in tobacco. Specifically, we focused on the response of tobacco NtWD40s to drought and salt stresses. For the randomly selected genes, many genes, such as NtWD40-074 and NtWD40-389, exhibited opposite expression patterns under drought and salt treatments, indicating that different NtWD40s may have different mechanisms for these two stresses. Silencing of NtWD40-389 (NtTTG1) also confirmed its important regulatory role in drought stress (Fig. 8). In summary, the expression patterns of the NtWD40s under growth, development, and stress conditions were diverse, indicating their functional diversity. However, additional studies were needed to explore the detailed mechanisms of candidate NtWD40s in tobacco plants.

Conclusion

The present study provides the first comprehensive and systematic identification and exploration of the WD40 gene family in tobacco. A total of 399 WD40 genes were identified and further classified into 6 clusters. The gene structure, conserved domains, and motifs of the NtWD40 genes were subsequently evaluated. Moreover, segmental duplication was found to be the major driver of the expansion of tobacco WD40 regulatory genes. Many NtWD40 genes were demonstrated to be expressed in a tissue-specific manner and to be involved in various biotic and abiotic stresses. Notably, NtTTG1-silenced tobacco plants exhibited significantly enhanced drought tolerance. These results could be useful for enhancing tobacco resistance to different stresses. Our study provided a valuable basis for further functional verification of WD40s in tobacco.

Methods

Identification and characterization of WD40 regulatory genes in tobacco

To identify all members of WD40 regulatory genes in tobacco, HMMER (v3.3.2) was employed using hidden markov model (HMM) based on the WD40 conserved domain (PF00400). All candidate sequences were confirmed using the Simple Modular Architecture Research Tool (SMART) (http://smart.embl-heidelberg.de/). TBtools was used to calculate the isoelectric points and molecular weights of tobacco WD40 proteins [52]. WoLF PSORT Online software (https://wolfpsort.hgc.jp/)was used to investigate the subcellular localization of tobacco WD40 proteins [53].

Phylogenetic classification, gene structures and conserved motifs analysis

The FastTre software was used to construct the maximum likelihood tree of the tobacco WD40 proteins, with a bootstrap value of 1000, and visualized by ITOL (https://itol.embl.de/) [54, 55]. The detailed protein sequences were listed in Table S9. To reveal the exon–intron organization of the WD40 genes in tobacco, TBtools software was used to determine the gene structures of each WD40 gene [52]. To analyze the conserved motifs in WD40 regulatory genes, MEME software was employed with maximum number of motifs of 10 [56]. Subsequently, InterProScan was used to annotate these identified motifs [57].

The chromosomal distributions, gene duplication and collinearity analysis

All tobacco WD40 regulatory genes were mapped to their corresponding chromosomes or scaffolds. Then, the origin, loss and expansion of WD40 regulatory genes in tobacco were inferred through phylogenetic tree. The duplication events of the tobacco WD40 genes were analyzed by BLAST (v2.13.0, e-value = 1e − 5), based on the similarity of amino acid sequences. Gene pairs with the similarity greater than 75% were considered as duplicated [58]. Gene pairs located on the same chromosome or scaffold are defined as tandem repeats, while those located on different chromosomes are defined as segmental repeats, and the rest are defined as uncertain. TBtools was used to calculate the non-synonymous (Ka) and synonymous (Ks) substitution rates of the duplicated events, and selective pressure was evaluated using the Ka/Ks ratio [52]. Collinearity analysis of WD40 genes was performed by TBtools and MCScanX between different species [59].

Enrichment analysis, promoter analysis and interaction network prediction

To explore the function of the tobacco WD40 genes, all protein sequences were aligned to the emapper database (v5.0.2), with the species range set to Viridiplantae [60]. Subsequently, the R package ClusterProfiler was used for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of tobacco WD40 genes [61,62,63,64]. For cis-acting element analysis, 2,000 bp upstream sequences from the transcription start site were extracted and submitted to the PlantCare database [65]. Mature miRNA sequences were downloaded from the miRBase database, and the psRNATarget database was used to search the regulatory relationship between miRNAs and tobacco WD40 gene [66, 67]. The regulatory relationship between tobacco WD40 genes and transcription factors (TFs) was retrieved from the Plant Transcription Factor Database PlantTFDB (http://planttfdb.cbi.pku.edu.cn) by searching for promoter sequences.

Expression profiles based on transcriptome data

To investigate tissue-specific expression patterns of WD40 regulatory genes, transcriptome data from eight representative tobacco tissues (root, stem, young leaf, mature leaf, senescent leaf, immature flower, mature flower, and senescent flower) were utilized. To explore the response of WD40 regulatory genes to different biotic and abiotic stresses, tobacco transcriptome samples subjected to various biotic stresses such as Potato Virus Y (PVY), black shank (BS), Cucumber mosaic virus (CMV), and Ralstonia solanacearum (RS), as well as abiotic stresses including high/low temperature, drought, salt, cadmium, and topping were obtained from the NCBI SRA database. The accession numbers for all the transcriptome data utilized in this study could be found in Table S10. The Salmon software (v1.10.1) was applied to perform quantification on all clean read [68], the R package ‘tximport’ (https://github.com/mikelove/tximport) was used to calculate the transcripts per kilobase million (TPM) values, and TBtools software was used for visualization of expression profiling [52].

Plant materials and stress treatment

The cultivated tobacco variety K326 was used to explore the expression of WD40 genes in various tissues and different stresses. Seedlings were cultured in plastic pots with a photoperiod of 16 h at 28 °C during the day and 23 °C at night. Root, stem, leaf, bud tissues at seedling stage (SS), vigorous growth stage (VG), and flowering stage (FS), as well as petal, stamen, pistil and seed samples at reproductive organs (RO) were sampled. To further investigate the expression of tobacco WD40 genes in trichome, trichome samples from leaf, petiole and stem were also collected. For salinity stress, 28 d old seedlings were treated with 400 mM NaCl for 3 h, 8 h, 12 h and 24 h, respectively. Same stage seedlings were also treated with plant hormones for 6 h, including ABA (10 μ M), IAA (10 μ M), Salicylic acid (SA) (10 μ M) and jasmonic acid (JA) (50 μ M and 100 μ M) [69, 70]. In addition, drought treatments were carried out by not watering for 7 days. Untreated seedlings were regarded as controls. All samples were independently collected with three independent replicates.

RNA isolation and relative expression analysis

Total RNA was extracted using SuperPure Plantpoly RNA Kit (Gene Answer, Beijing, China) and DNA contamination was eliminated with RNase-free DNase I (Gene Answer). First-strand cDNA synthesis was performed using 1 μg of total RNA as a template by using reverse transcription kit (Trans, Beijing, China), and diluted to a concentration of 50 ng/µL. Real-Time Quantitative Reverse Transcription PCR (qRT-PCR) of tobacco WD40 genes was performed with SYBR Green kit (Trans, Beijing, China) with a reaction system of 20 μL. The PCR procedure used a two-step method consisting of an initial step of 95 °C for 30 s, followed by 40 cycles of 95 °C for 10 s and 60 °C for 30 s. The relative gene expression was calculated by 2−∆∆Ct method and GAPDH was used as the internal reference gene. All qRT-PCR specific primers were listed in Table S11. Three independent biological replicates were conducted for qRT-PCR analysis.

Subcellular localization analysis

To verify the subcellular localization of WD40 genes in tobacco, primers were designed for segment-length coding sequences (CDS) of NtWD40-389 (which named as NtTTG1 for its particularly high sequence similarity with TTG1 gene in Arabidopsis and other plants). The subcellular localization vector NtTTG1-PC1300s-GFP was constructed by homologous recombination (Fig. S7A). Generally, the plasmids (NtTTG1-PC1300s-GFP and PC1300s-GFP) were transformed into Agrobacterium GV3101, and positive colonies were verified and expanded to OD value 1 at 28 °C. These colonies were re-suspended with MgCl2 + As + MES. After 3 h in darkness, they were injected into leaves of Nicotiana benthamiana. After injection, tobacco plants were cultivated in darkness for 1 d, and in normal condition for 1 d. Finally, GFP signals were detected using a confocal microscopy system (Nikon C2-ER, Japan).

Functional assessment of NtTTG1 by virus-induced gene silencing

Wild-type N. benthamiana plants were cultivated in a climate-controlled chamber at a temperature ranging from 20–22 °C and exposed to a light period of 16 h. The virus-induced gene silencing (VIGS) vector NtTTG1-pYY13 was constructed similar with subcellular localization analysis (Fig. S7B). The plasmid (NtTTG1-pYY13) was transformed into Agrobacterium GV3101, and positive colonies were verified and expanded to OD value 1 at 28 °C. In addition, TRV1, TRV2 and PDS were also expanded. After these colonies re-suspended with MgCl2 + As + MES and kept in darkness for 3 h, NtTTG1-PYY13, TRV2 and PDS were mixed with TRV1 and injected into N. benthamiana leaves, respectively. After in darkness for 1 d, the injected plants were cultivated in a well-regulated growth chamber for 14 d. When PDS-silenced plants showed a severe phenotype with complete photo bleaching in all newly emerging leaves, qRT-PCR analysis was performed on uninfected leaves for each group. Subsequently, silenced plants were subjected to drought treatment, and corresponding phenotype were observed.

Availability of data and materials

All data generated or analyzed in this study are included in the materials and methods section of this article.

References

  1. Liu M, Wang Y, Teng F, Mai X, Wang X, Su MY, Stjepanovic G. Structure of the DDB1-AMBRA1 E3 ligase receptor complex linked to cell cycle regulation. Nat Commun. 2023;14(1):7631.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  2. Kim YJ, Jung KH. WD40-domain protein GORI is an integrative scaffold that is required for pollen tube growth in rice. Plant Signal Behav. 2023;18(1):2082678.

    Article  MathSciNet  PubMed  Google Scholar 

  3. Jain BP, Pandey S. WD40 Repeat Proteins: Signalling Scaffold with Diverse Functions. Protein J. 2018;37(5):391–406.

    Article  CAS  PubMed  Google Scholar 

  4. Hu XJ, Li T, Wang Y, Xiong Y, Wu XH, Zhang DL, Ye ZQ, Wu YD. Prokaryotic and Highly-Repetitive WD40 Proteins: A Systematic Study. Sci Rep. 2017;7(1):10585.

    Article  ADS  PubMed  PubMed Central  Google Scholar 

  5. Fong HK, Hurley JB, Hopkins RS, Miake-Lye R, Johnson MS, Doolittle RF, Simon MI. Repetitive segmental structure of the transducin beta subunit: homology with the CDC4 gene and identification of related mRNAs. Proc Natl Acad Sci USA. 1986;83(7):2162–6.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  6. Wu XH, Chen RC, Gao Y, Wu YD. The effect of Asp-His-Ser/Thr-Trp tetrad on the thermostability of WD40-repeat proteins. Biochemistry. 2010;49(47):10237–45.

    Article  CAS  PubMed  Google Scholar 

  7. Xu C, Min J. Structure and function of WD40 domain proteins. Protein Cell. 2011;2(3):202–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Wall MA, Coleman DE, Lee E, Iñiguez-Lluhi JA, Posner BA, Gilman AG, Sprang SR. The structure of the G protein heterotrimer Gi alpha 1 beta 1 gamma 2. Cell. 1995;83(6):1047–58.

    Article  CAS  PubMed  Google Scholar 

  9. Andrade MA, Perez-Iratxeta C, Ponting CP. Protein repeats: structures, functions, and evolution. J Struct Biol. 2001;134(2–3):117–31.

    Article  CAS  PubMed  Google Scholar 

  10. He S, Tong X, Han M, Hu H, Dai F. Genome-Wide Identification and Characterization of WD40 Protein Genes in the Silkworm, Bombyx mori. Int J Mol Sci. 2018;19(2):527.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Li Q, Zhao P, Li J, Zhang C, Wang L, Ren Z. Genome-wide analysis of the WD-repeat protein family in cucumber and Arabidopsis. Mol Gen Genom. 2014;289(1):103–24.

    Article  CAS  Google Scholar 

  12. Ouyang Y, Huang X, Lu Z, Yao J. Genomic survey, expression profile and co-expression network analysis of OsWD40 family in rice. BMC Genomics. 2012;13:100.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Tan L, Salih H, Htet NNW, Azeem F, Zhan R. Genomic analysis of WD40 protein family in the mango reveals a TTG1 protein enhances root growth and abiotic tolerance in Arabidopsis. Sci Rep. 2021;11(1):2266.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Salih H, Gong W, Mkulama M, Du X. Genome-wide characterization, identification, and expression analysis of the WD40 protein family in cotton. Genome. 2018;61(7):539–47.

    Article  CAS  PubMed  Google Scholar 

  15. Kong D, Li M, Dong Z, Ji H, Li X. Identification of TaWD40D, a wheat WD40 repeat-containing protein that is associated with plant tolerance to abiotic stresses. Plant Cell Rep. 2015;34(3):395–410.

    Article  CAS  PubMed  Google Scholar 

  16. Chen C, Yang Y, Pan L, Xia W, Xu L, Hua B, Zhang Z, Miao M. Genome-Wide Identification of WD40 Proteins in Cucurbita maxima Reveals Its Potential Functions in Fruit Development. Genes. 2023;14(1):220.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Wu Z, Zhang T, Li J, Chen S, Grin IR, Zharkov DO, Yu B, Li H. Genome-wide analysis of WD40 protein family and functional characterization of BvWD40-82 in sugar beet. Front Plant Sci. 2023;14:1185440.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Liu Z, Liu Y, Coulter JA, Shen B-c, Li Y, Li C, Cao Z-M, Zhang J: The WD40 Gene Family in Potato (Solanum Tuberosum L.): Genome-Wide Analysis and Identification of Anthocyanin and Drought-Related WD40s. Agronomy 2020.

  19. Yu H. Cdc20: a WD40 activator for a cell cycle degradation machine. Mol Cell. 2007;27(1):3–16.

    Article  CAS  PubMed  Google Scholar 

  20. Stirnimann CU, Petsalaki E, Russell RB, Müller CW. WD40 proteins propel cellular networks. Trends Biochem Sci. 2010;35(10):565–74.

    Article  CAS  PubMed  Google Scholar 

  21. Saeki M, Irie Y, Ni L, Yoshida M, Itsuki Y, Kamisaki Y. Monad, a WD40 repeat protein, promotes apoptosis induced by TNF-alpha. Biochem Biophys Res Commun. 2006;342(2):568–72.

    Article  CAS  PubMed  Google Scholar 

  22. Fimia GM, Stoykova A, Romagnoli A, Giunta L, Di Bartolomeo S, Nardacci R, Corazzari M, Fuoco C, Ucar A, Schwartz P, et al. Ambra1 regulates autophagy and development of the nervous system. Nature. 2007;447(7148):1121–5.

    Article  ADS  CAS  PubMed  Google Scholar 

  23. Xiong Y, Contento AL, Bassham DC. AtATG18a is required for the formation of autophagosomes during nutrient stress and senescence in Arabidopsis thaliana. Plant J. 2005;42(4):535–46.

    Article  CAS  PubMed  Google Scholar 

  24. Pickles LM, Roe SM, Hemingway EJ, Stifani S, Pearl LH. Crystal structure of the C-terminal WD40 repeat domain of the human Groucho/TLE1 transcriptional corepressor. Structure. 2002;10(6):751–61.

    Article  CAS  PubMed  Google Scholar 

  25. Chen RH, Miettinen PJ, Maruoka EM, Choy L, Derynck R. A WD-domain protein that is associated with and phosphorylated by the type II TGF-beta receptor. Nature. 1995;377(6549):548–52.

    Article  ADS  CAS  PubMed  Google Scholar 

  26. Datta PK, Moses HL. STRAP and Smad7 synergize in the inhibition of transforming growth factor beta signaling. Mol Cell Biol. 2000;20(9):3157–67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Suganuma T, Pattenden SG, Workman JL. Diverse functions of WD40 repeat proteins in histone recognition. Genes Dev. 2008;22(10):1265–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Song JJ, Garlick JD, Kingston RE. Structural basis of histone H4 recognition by p55. Genes Dev. 2008;22(10):1313–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Li H, Luan S. The cyclophilin AtCYP71 interacts with CAF-1 and LHP1 and functions in multiple chromatin remodeling processes. Mol Plant. 2011;4(4):748–58.

    Article  CAS  PubMed  Google Scholar 

  30. Kim DK, Redon CE, Aladjem MI, Kim HK, Jang SM. Molecular double clips within RepID WD40 domain control chromatin binding and CRL4-substrate assembly. Biochem Biophys Res Commun. 2021;567:208–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Yuan F, Leng B, Zhang H, Wang X, Han G, Wang B. A WD40-Repeat Protein From the Recretohalophyte Limonium bicolor Enhances Trichome Formation and Salt Tolerance in Arabidopsis. Front Plant Sci. 2019;10:1456.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Kim YJ, Kim MH, Hong WJ, Moon S, Kim EJ, Silva J, Lee J, Lee S, Kim ST, Park SK, et al. GORI, encoding the WD40 domain protein, is required for pollen tube germination and elongation in rice. Plant J. 2021;105(6):1645–64.

    Article  CAS  PubMed  Google Scholar 

  33. Liu WC, Li YH, Yuan HM, Zhang BL, Zhai S, Lu YT. WD40-REPEAT 5a functions in drought stress tolerance by regulating nitric oxide accumulation in Arabidopsis. Plant Cell Environ. 2017;40(4):543–52.

    Article  CAS  PubMed  Google Scholar 

  34. Tian G, Wang S, Wu J, Wang Y, Liu S, Han D, Xia G, Wang M. Allelic variation of TaWD40–4B.1 contributes to drought tolerance by modulating catalase activity in wheat. Nat Commun. 2023;14(1):1200.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  35. Shoeva OY, Mukhanova MA, Zakhrabekova S, Hansson M. Ant13 Encodes Regulatory Factor WD40 Controlling Anthocyanin and Proanthocyanidin Synthesis in Barley (Hordeum vulgare L.). J Agric Food Chem. 2023;71(18):6967–77.

    Article  CAS  PubMed  Google Scholar 

  36. Sun X, Zhang Z, Li J, Zhang H, Peng Y, Li Z. Uncovering Hierarchical Regulation among MYB-bHLHWD40 Proteins and Manipulating Anthocyanin Pigmentation in Rice. Int J Mol Sci. 2022;23(15):8203.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Yue M, Jiang L, Zhang N, Zhang L, Liu Y, Lin Y, Zhang Y, Luo Y, Zhang Y, Wang Y, et al. Regulation of flavonoids in strawberry fruits by FaMYB5/FaMYB10 dominated MYB-bHLH-WD40 ternary complexes. Front Plant Sci. 2023;14:1145670.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Larkin JC, Walker JD, Bolognesi-Winfield AC, Gray JC, Walker AR. Allele-specific interactions between ttg and gl1 during trichome development in Arabidopsis thaliana. Genetics. 1999;151(4):1591–604.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Ge J, Li B, Shen D, Xie J, Long J, Dong H. Tobacco TTG2 regulates vegetative growth and seed production via the predominant role of ARF8 in cooperation with ARF17 and ARF19. BMC Plant Biol. 2016;16(1):126.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Ishida T, Hattori S, Sano R, Inoue K, Shirano Y, Hayashi H, Shibata D, Sato S, Kato T, Tabata S, et al. Arabidopsis TRANSPARENT TESTA GLABRA2 is directly regulated by R2R3 MYB transcription factors and is involved in regulation of GLABRA2 transcription in epidermal differentiation. Plant Cell. 2007;19(8):2531–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Lim SH, Kim DH, Lee JY. RsTTG1, a WD40 Protein, Interacts with the bHLH Transcription Factor RsTT8 to Regulate Anthocyanin and Proanthocyanidin Biosynthesis in Raphanus sativus. Int J Mol Sci. 2022;23(19):11973.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Zhao M: Regulation of Arabidopsis trichome patterning and anthocyanin biosynthesis by the TTG1 -bHLH -MYB complex. Dissertations & Theses Gradworks 2007.

  43. 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–30.

    Article  PubMed  Google Scholar 

  44. International Wheat Genome Sequencing Consortium (IWGSC): Shifting the limits in wheat research and breeding using a fully annotated reference genome. Science (New York, NY) 2018, 361(6403).

  45. Panchy N, Lehti-Shiu M, Shiu SH. Evolution of Gene Duplication in Plants. Plant Physiol. 2016;171(4):2294–316.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Chen L, Cui Y, Yao Y, An L, Bai Y, Li X, Yao X, Wu K. Genome-wide identification of WD40 transcription factors and their regulation of the MYB-bHLH-WD40 (MBW) complex related to anthocyanin synthesis in Qingke (Hordeum vulgare L. var. nudum Hook. f.). BMC Genomics. 2023;24(1):166.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Liang G, He H, Li Y, Ai Q, Yu D. MYB82 functions in regulation of trichome development in Arabidopsis. J Exp Bot. 2014;65(12):3215–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Wang B, Luo Q, Li Y, Du K, Wu Z, Li T, Shen WH, Huang CH, Gan J, Dong A. Structural insights into partner selection for MYB and bHLH transcription factor complexes. Nature plants. 2022;8(9):1108–17.

    Article  ADS  CAS  PubMed  Google Scholar 

  49. Feyissa BA, Arshad M, Gruber MY, Kohalmi SE, Hannoufa A. The interplay between miR156/SPL13 and DFR/WD40-1 regulate drought tolerance in alfalfa. BMC Plant Biol. 2019;19(1):434.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Liu C, Xu Y, Long D, Cao B, Hou J, Xiang Z, Zhao A. Plant G-protein β subunits positively regulate drought tolerance by elevating detoxification of ROS. Biochem Biophys Res Commun. 2017;491(4):897–902.

    Article  CAS  PubMed  Google Scholar 

  51. Chen W, Chen L, Zhang X, Yang N, Guo J, Wang M, Ji S, Zhao X, Yin P, Cai L, et al. Convergent selection of a WD40 protein that enhances grain yield in maize and rice. Science. 2022;375(6587):eabg7985.

    Article  CAS  PubMed  Google Scholar 

  52. 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–202.

    Article  CAS  PubMed  Google Scholar 

  53. Horton P, Park KJ, Obayashi T, Fujita N, Harada H, Adams-Collier CJ, Nakai K. WoLF PSORT: protein localization predictor. Nucleic Acids Res. 2007;35(Web Server issue):W585–587.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Price MN, Dehal PS, Arkin AP. FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol Biol Evol. 2009;26(7):1641–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Letunic I, Bork P. Interactive Tree Of Life (iTOL) v4: recent updates and new developments. Nucleic Acids Res. 2019;47(W1):W256–w259.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Bailey TL, Boden M, Buske FA, Frith M, Grant CE, Clementi L, Ren J, Li WW, Noble WS. MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 2009;37(Web Server issue):W202–208.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Quevillon E, Silventoinen V, Pillai S, Harte N, Mulder N, Apweiler R, Lopez R. InterProScan: protein domains identifier. Nucleic Acids Res. 2005;33(Web Server issue):W116–120.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, Madden TL. BLAST+: architecture and applications. BMC Bioinformatics. 2009;10:421.

    Article  PubMed  PubMed Central  Google Scholar 

  59. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Cantalapiedra CP, Hernández-Plaza A, Letunic I, Bork P, Huerta-Cepas J. eggNOG-mapper v2: Functional Annotation, Orthology Assignments, and Domain Prediction at the Metagenomic Scale. Mol Biol Evol. 2021;38(12):5825–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Wu T, Hu E, Xu S, Chen M, Guo P, Dai Z, Feng T, Zhou L, Tang W, Zhan L, et al. clusterProfiler 4.0: A universal enrichment tool for interpreting omics data. Innovation (Camb). 2021;2(3):100141.

    CAS  PubMed  Google Scholar 

  62. Kanehisa M, Furumichi M, Sato Y, Kawashima M, Ishiguro-Watanabe M. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2023;51(D1):D587–d592.

    Article  CAS  PubMed  Google Scholar 

  63. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28(11):1947–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Lescot M, Déhais P, Thijs G, Marchal K, Moreau Y, Van de Peer Y, Rouzé 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–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Kozomara A, Birgaoanu M, Griffiths-Jones S. miRBase: from microRNA sequences to function. Nucleic Acids Res. 2019;47(D1):D155–d162.

    Article  CAS  PubMed  Google Scholar 

  67. Dai X, Zhuang Z, Zhao PX. psRNATarget: a plant small RNA target analysis server (2017 release). Nucleic Acids Res. 2018;46(W1):W49–w54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Patro R, Duggal G, Love MI, Irizarry RA, Kingsford C. Salmon provides fast and bias-aware quantification of transcript expression. Nat Methods. 2017;14(4):417–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Torun H, Novák O, Mikulík J, Strnad M, Ayaz FA. The Effects of Exogenous Salicylic Acid on Endogenous Phytohormone Status in Hordeum vulgare L. under Salt Stress. Plants (Basel, Switzerland). 2022;11(5):618.

  70. Braybrook SA. Analyzing Cell Wall Elasticity After Hormone Treatment: An Example Using Tobacco BY-2 Cells and Auxin. Methods Mol Biol. 2017;1497:125–33.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We would like to acknowledge all researchers in our laboratory for their help.

Funding

This work was supported by the Zhengzhou Tobacco Research Institute [CNTC: 110202101008(JY-08), 110202201001(JY-01), 110202202038], Natural Science Foundation of HeNan (232300420220); Beijing Life Science Academy [BLSA: 2023000CC0110].

Author information

Authors and Affiliations

Authors

Contributions

#These authors contributed equally. JJJ and PJC conceived and designed the experiments. HS, PL, TJM and JFZ performed bioinformatics data analysis. LJM and ZCQ did qRT-PCR experiments, the subcellular localization and VIGS experiments. LJM, HS, WZ, NL, HL, PJC and JJJ wrote the manuscript and all authors read and approved the final version.

Corresponding author

Correspondence to Jingjing Jin.

Ethics declarations

Ethics approval and consent to participate

The seeds of K326 and N. benthamiana used in this study were preserved by Tobacco Research Institute of Chinese Academy of Agricultural Sciences. These two materials were also globally common varieties, and we ensured that we obtained the necessary permits and permissions when collecting the plant material. The collection and use of these seeds, including experimental research and field studies, were performed completely in accordance with relevant institutional, national, and international guidelines and legislation.

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

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Meng, L., Su, H., Qu, Z. et al. Genome-wide identification and analysis of WD40 proteins reveal that NtTTG1 enhances drought tolerance in tobacco (Nicotiana tabacum). BMC Genomics 25, 133 (2024). https://doi.org/10.1186/s12864-024-10022-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-024-10022-w

Keywords