Skip to main content

Comprehensive analysis of Alfin-like transcription factors associated with drought and salt stresses in wheat (Triticum aestivum L.)

Abstract

Background

Alfin-like proteins are a kind of plant-specific transcription factors, and play vital roles in plant growth, development and stress responses.

Results

In this study, a total of 27 Alfin-like transcription factors were identified in wheat. TaAL genes were unevenly distributed on chromosome. Phylogenetic analysis showed TaAL genes were divided into AL-B and AL-C subfamilies, and TaALs with closer evolutionary relationships generally shared more similar exon-intron structures and conserved motifs. The cis-acting element analysis showed MBS, ABRE and CGTCA-motif were the most common in TaAL promoters. The interacting proteins and downstream target genes of TaAL genes were also investigated in wheat. The transcriptome data and real-time PCR results indicated TaAL genes were differentially expressed under drought and salt stresses, and TaAL1-B was significantly up-regulated in response to drought stress. In addition, association analysis revealed that TaAL1-B-Hap-I allelic variation had significantly higher survival rate compared to TaAL1-B-Hap-II under drought stress.

Conclusions

These results will provide vital information to increase our understanding of the Alfin-like gene family in wheat, and help us in breeding better wheat varieties in the future.

Peer Review reports

Introduction

Wheat (Triticum aestivum L.) is one of the most important grain crops worldwide. As a sessile organism, wheat has to suffer from all kinds of adverse environmental factors, such as high salinity and drought, which seriously affect the yield and quality of wheat [1]. Therefore, plants have gradually developed intricate regulatory mechanisms over a long evolutionary period to avoid or defend against adverse environmental factors [2]. Transcription factors (TFs) are essential regulators in plant growth, development and stress responses through regulating the expression of the downstream target genes [3, 4].

The Alfin-like (AL) family proteins are a kind of plant-specific TFs [5]. AL proteins mainly contain two conserved domains, a conserved Alfin domain at the N-terminus consisting of 140 conserved amino acid residues and a PHD-finger domain with conserved Cys4-His-Cys3 (C4HC3) residues at the C-terminus [6, 7]. The PHD-finger are involved in epigenetics and chromatin-mediated transcriptional regulation [8, 9], which is responsible for the nuclear localization of the AL proteins [10, 11].

AL TFs have been identified to participate in various stress responses, such as salt, drought, and powdery mildew resistance [5,6,7, 12]. The first AL gene, Alfin1, was identified in alfalfa (Medicago sativa L.) [13], which could enhance plant root growth, and enhance salt tolerance via improving the expression of salt-induced MsPRP2 gene [14, 15]. Overexpression of Conringia planisiliqua CpAL2 gene conferred Arabidopsis thaliana seedlings with drought and salt resistance [16]. In Arabidopsis, AtAL5 improved the salt, drought and freezing tolerance by binding to G-rich motifs (GTGGNG or GNGGTG) to regulate the expression of downstream stress-related genes [10]. Overexpression of Atriplex hortensis AhAL1 gene enhanced drought and salt resistance in transgenic Arabidopsis plants [7]. Soybean (Glycine max (L.) Merr.) AL gene GmPHD2 improved salt stress tolerance in transgenic Arabidopsis plants [17]. The heterologous expression of MdAL4 gene improved the drought tolerance of transgenic Arabidopsis [6]. AL was also as a histone reader to bind to active histone mark H3K4me2/3 [18, 19]. AL2-PRC1 complexes promoted seed germination by active or repressive H3K4me3-to-H3K27me3 chromatin state switch, leading to the repression of seed developmental genes in Arabidopsis [20].

Alfin-like family genes have been identified in diverse plant species, such as Arabidopsis [10], rice (Oryza sativa L.) [5], apple (Malus domestica) [6], and Chinese cabbage (Brassica rapa) [21]. However, a comprehensive analysis of Alfin-like gene family in wheat has not been investigated. In present study, a genome-wide identification of Alfin-like genes in wheat was performed, and their physicochemical properties, gene structure, conserved motif, evolutionary relationship, chromosome distribution, gene duplication and regulatory network were also analyzed. Furthermore, the tissue-specific and expression patterns under drought and salt stresses of Alfin-like gene were determined. In addition, we used wheat resequencing data from Wheat Union database [22] and drought-related trait data acquired from published study [23] to perform association analysis with TaAL-B gene. These results will provide a valuable foundation for further functional study of Alfin-like genes in wheat.

Results

Identification and phylogenetic analysis of Alfin-like genes

Through a Hidden Markov Model (HMM) search, a total of 27 members of Alfin-like gene family were identified in wheat, and renamed them according to their collinearity analysis and subgenome distribution (Table S1 and Fig. 1). The phylogenetic tree was constructed with AL protein sequences from moss (Physcomitrella patens, 7 PpALs), dicot (Arabidopsis, 7 AtALs), and monocot (Oryza sativa, 10 OsAls; Triticum aestivum, 27 TaALs) (Fig. 1). Based on the phylogenetic relationship of AL proteins, they were classified into three subfamilies: AL-A, AL-B, and AL-C (Fig. 1). AL-B and AL-C subfamilies included twelve and fifteen TaAL protein members, respectively, whereas no TaAL protein member belonged to AL-A subfamily (Fig. 1 and Table S1). AL-A members only existed in ancient plant species, e.g., P. patens. AL-B and AL-C were present in higher plants such as monocotyledonous and dicotyledonous plants [5]. Interestingly, AL-B members in dicot (AtALs) and monocot (OsALs and TaALs) were clustered in separate branches (Fig. 1), suggesting that AL-B proteins occurred after differentiation of the monocot and dicot. In contrast, AL-C members might occur before differentiation of the monocot and dicot. These results suggested that AL proteins evolved with the evolution of plants from ancient plant to higher plant species, and AL-B and AL-C proteins evolved independently in higher plant species.

Fig. 1
figure 1

The neighbor-joining (NJ) phylogenetic tree of AL proteins. The phylogenetic tree was constructed with AL protein sequences from Arabidopsis (At), Oryza sativa (Os), Physcomitrella patens (Pp), and Triticum aestivum (Ta) with bootstrap values of 1000 replicates. The AL proteins are classified into three subfamilies: Alfin-like A (AL-A), Alfin-like A (AL-B), and Alfin-like A (AL-C). AL protein members of different subfamilies are distinguished by different colors

Characterization, gene structure and conserved motif analysis of TaAL genes

The length of the identified 27 TaAL proteins ranged from 219 to 304 amino acid residues, and the molecular weights ranged from 24.95 to 34.6 kDa. The pI values varied from 4.55 (TaAL10-A) to 6.31 (TaAL9-B), with the calculated grand average of hydrophilic index (GRAVY) ranging from − 0.864 (TaAL3-A) to -0.489 (TaAL4-B), suggesting that these 27 TaAL genes encoded highly hydrophilic proteins. The subcellular localization prediction suggested all TaALs were located in the nucleus (Table S1).

To investigate the structural characteristic of TaAL genes, the exon–intron structures and conserved motifs were analyzed (Fig. 2). The exon number of AL-B genes varied from 2 to 5. AL-C subfamily genes included five exons and four introns, except for TaAL2-A contained three exons and two introns (Fig. 2B). Meanwhile, five conserved motifs were identified among TaAL proteins by using the MEME program (Fig. 2C). All TaALs contained motif 1, 2, 4, and 5. Motif 2, 3 and 4 formed Alfin domain, and motif 1 and 5 made up PHD domain (Fig. 2C and Fig. S1). These results also verified the reliability of the identified TaAL proteins. TaALs with closer evolutionary relationships generally shared more similar exon-intron structures and conserved motifs (Fig. 2).

Fig. 2
figure 2

Phylogenetic classification, exon–intron structures and conserved domains of TaAL genes. (A) Phylogenetic analysis of TaAL gene family. (B) The introns, exons and UTR are represented by black lines, blue boxes and yellow boxes, respectively. (C) Each motif is represented using colored boxes

Chromosomal location, collinearity and Ka/Ks analysis of TaAL genes

The chromosomal location and collinearity of TaAL genes were analyzed according to their genomic information (Fig. 3 and Fig. S2). TaAL genes were unevenly distributed on chromosome 1 A, 1B, 1D, 2 A, 2B, 2D, 3 A, 3B, 3D, 4 A, 4B, 4D, 5B, 5D, 7 A, 7D, and Un (Fig. S2). To investigate the expansion and evolution of Alfin-like family genes in wheat, a collinearity analysis of TaAls was performed using TGT (Triticeae-Gene Tribe) website (Fig. 3). A total of thirty-five paralogous gene pairs of TaAL genes were identified in wheat genome. Expect TaAL7-A, TaAL8-A, TaAL10-D and TaAL11-D were tandem duplicate and block duplicate, the other TaAL genes all undergone block duplicate events (Table S2). Ka/Ks (the nonsynonymous and synonymous substitution ratio) values of four paralogous gene pairs (TaAL8-A/TaAL10-D, TaAL8-A/TaAL11-Un, TaAL10-D/TaAL11-D, TaAL10-A/TaAL11-D) were more than 1, suggesting these genes have undergone positive selection after duplication events (Table S2). And the Ka/Ks ratio of the other thirty-one paralogous gene pairs were less than 1, suggesting these thirty-one paralogous gene pairs undergone purifying selection to maintain the function of TaAL gene family (Table S2).

Fig. 3
figure 3

Chromosomal localizations and syntenic relationships among TaAL genes in wheat. Red lines in the highlight indicate the syntenic TaAL gene pairs

Analysis of cis-acting elements of TaAL promoters

The variable cis-acting elements in gene promoters can bind to various transcription factors to regulate gene expression, thus revealing the differences biological function of genes. To study the function of TaAL genes, 1500 bp promoter sequences of TaALs were analyzed through PlantCARE database to identify the cis-acting elements (Fig. 4). A total of 33 cis-acting elements were identified and divided into four categories, including ten stress responsive-elements, six hormone responsive-elements, twelve light responsive-elements, and five growth and development related-elements. The TaAL promoters mainly contained stress responsive- and hormone responsive-elements, especially MBS (MYB binding site), ABRE (ABA-responsive element) and CGTCA-motif (MeJA-responsive element) were the most common in TaAL promoters (Fig. 4). These results suggested that TaAL genes might play crucial roles in response to various stresses in wheat.

Fig. 4
figure 4

Analysis of cis -acting elements in TaALpromoters. (A) Phylogenetic analysis of TaAL gene family. (B) Cis-acting elements in the promoters of TaAL genes. The different colors and numbers of grids indicated the numbers of different promoter elements. (C) The histograms of different colors represented the sum of the cis-acting elements in each category

Expression profiles of TaALs

The expression profiles of the TaAL genes in different tissues and under drought and salt stresses during the wheat seedling stage were determined by using transcriptome data and real-time PCR (Fig. 5 and Fig. S3). The results showed TaAL1-B, TaAL2-B and TaAL5-B were expressed at a higher level in shoot tissues and at a lower level in root tissues (Fig. S3). The transcriptome data revealed that the expression levels of TaAL1-A/B/D, TaAL2-A/B/D, TaAL5-A, TaAL6-B/D and TaAL9-B/D were significantly up-regulated under drought stress, and TaAL3-A/B/D, TaAL4- B/D, TaAL5-B/D, TaAL6-A were down-regulated (Fig. 5A and Table S3). Consistently, real-time PCR results showed that the expression levels of TaAL1-B and TaAL2-B were up-regulated and reached the highest expression level at 36 h and 24 h after PEG-induced drought stress with approximately 2.05- and 1.50-fold compared with the control, respectively. TaAL5-B was down-regulated under PEG stress (Fig. 5B). After salt stress treatment, the expression levels of TaAL1-A/B/D, TaAL3-A/B/D, TaAL6-A/B/D, TaAL4-A and TaAL9-B were obviously up-regulated, and TaAL2-A/B/D, TaAL4-B, TaAL5-A/B/D, and TaAL9-D were down-regulated (Fig. 5C and Table S3). The real-time PCR results also showed that TaAL2-B and TaAL5-B were down-regulated under salt stress (Fig. 5D). These results suggested that TaAL genes had different expression levels under drought and salt stresses and might be involved in the regulatory pathways of drought and salt stresses.

Fig. 5
figure 5

Expression patterns of TaAL genes in response to drought and salt stresses. (A) RNA-seq analysis of the expression profiles of TaAL genes in response to drought stress for 0 h, 2 h, and 7 h, respectively. Fragments per kilobase of exon per million mapped fragments (FPKM) values were used to measure the expression levels of genes. (B) Real-time PCR analysis of the expression profiles of TaAL genes in wheat seedling leaves at 0 h, 12 h, 24 h, and 36 h after PEG stress treatment.(C) RNA-seq analysis of the expression profiles of TaAL genes in response to salt stress for 6 h, 12 h, 24 h, and 36 h, respectively. Fragments per kilobase of exon per million mapped fragments (FPKM) values were used to measure the expression levels of genes. (D) Real-time PCR analysis of the expression profiles of TaAL genes in wheat seedling leaves at 0 h, 12 h, 24 h, and 36 h after salt stress treatment. The values show the mean ± SD. Significant statistical analysis was carried out by Student’s t -test (*p < 0.05; **p < 0.01)

The interacting proteins and downstream target genes analysis of TaAL genes

To determine the regulatory mechanism of TaAL genes in response to drought and salt stresses, the interacting proteins and downstream target genes of TaALs were analyzed (Figs. 6 and 7). A total of 25 interacting proteins of TaAL were identified by using the STRING [24] and plant.MAP [25] database, including Alfin-like transcription factor, beta-catenin-like protein, histone-lysine N-methyltransferase, multiple organellar RNA editing factor, pentatricopeptide repeat-containing protein, protein RBL, protein TRAUCO, and ubiquitin-40 S ribosomal protein (Fig. 6A and Table S4). Cluster analysis of gene expression showed TaAL1-A/D exhibited the most similar expression patterns with interacting protein TRAUCO (TraesCS5D02G141000) under drought and salt stresses. TaAL2-A/B/D had the most similar expression patterns with TRAUCO (TraesCS5A02G132800) and ubiquitin-40 S ribosomal protein (TraesCS1A02G397400). TaAL5-A/B/D showed the most similar expression patterns with ubiquitin-40 S ribosomal protein (TraesCS1D02G109800, and TraesCS1A02G166000) (Fig. 6B). These interacting proteins that have similar expression patterns with TaALs might perform important roles through interacting with TaAL proteins under drought and salt stresses.

Fig. 6
figure 6

Protein–protein interaction (A) and expression profiles analysis of TaAL. Fragments per kilobase of exon per million mapped fragments (FPKM) values were used to measure the expression levels of genes

As transcription factors, TaALs could regulate downstream target genes in response to drought and salt stresses. Therefore, the downstream target genes of TaALs were obtained through identifying the promoters of co-expression genes including G-rich motifs (Table S5 and S6). The results showed the downstream target genes of TaAL2-A (7), TaAL2-B (28), TaAL2-D (22), TaAL3-A (6), TaAL3-B (1), TaAL3-D (10), TaAL4-A (11), TaAL4-D (5), TaAL5-A (7), TaAL5-B (5), TaAL5-D (7), TaAL6-A (27), TaAL6-B (34), TaAL6-D (31), TaAL10-D (1), TaAL10-Un (5) and TaAL11-D (7) were identified (Fig. 7 and Table S5). The expression patterns of TaALs and their downstream target genes under drought and salt stresses were also analyzed (Fig. S4). Cluster analysis of gene expression showed TaAL2-A, TaAL2-B and TaAL2-D had the most similar expression patterns with chromatin remodeling protein (TraesCS5B02G450300), actin-related protein (TraesCS2B02G597500) and E3 ubiquitin-protein ligase (TraesCS7A02G350000), respectively (Fig. 7 and Fig. S4). Additionally, the most similar expression patterns were observed between TaAL5-A with upstream activation factor subunit (TraesCS5D02G411700), TaAL5-B with bZIP transcription factor (TraesCS7A02G268700) and upstream activation factor subunit (TraesCS5B02G406200), TaAL5-D with bZIP transcription factor (TraesCS2A02G352100). These genes that have similar expression patterns with TaALs might perform important roles through regulating the expression by TaAL under drought and salt stresses.

Fig. 7
figure 7

The downstream target gene of TaAL proteins

Association analysis of TaAL1-B gene haplotypes with drought-related traits in wheat

According to transcriptome data and real-time PCR results, TaAL1-B gene was highly expressed under drought stress. Therefore, we detected the variation sites in 1500 bp promoter and CDS (coding sequence) region of TaAL1-B gene by using 681 wheat resequencing data in Wheat Union database (Table S7). The results showed that 7 SNP (single-nucleotide polymorphism) variants were detected, i.e., C/T (-1081 bp) and A/T (-514 bp) were located in promoter region, G/A (3559 bp), A/T (3580 bp), G/A (3602 bp), A/T (3653 bp), and T/G (3719 bp) were located in intron region (Fig. 8A). Based on these SNP variants, two haplotypes were identified and named TaAL1-B-Hap-I and TaAL1-B-Hap-II (Fig. 8A and Table S8). Due to SNP variants, FAR1 (far-red impaired response 1) TF could bind to TaAL1-B-Hap-I promoter at -1081 bp, however, MYB TF bound to TaAL1-B-Hap-II promoter at -1081 bp. Otherwise, MYB TF could bind to TaAL1-B-Hap-I promoter at -514 bp, but no TF was detected to bind to TaAL1-B-Hap-II promoter at -514 bp (Fig. 8B). To analysis the effect of TaAL1-B allele variation on survival rate under drought stress in wheat, we used published data from 39 of these 681 wheat varieties to conduct the associations analysis of TaAL1-B genes with survival rate under drought stress (Table S9). The results showed TaAL1-B-Hap-I allelic variation had significantly higher survival rate compared to TaAL1-B-Hap-II under drought stress (Fig. 8B and Table S9). Therefore, the transcription factor binding site in the promoter of TaAL1-B gene with Hap-I might be involved in regulating wheat drought tolerance.

Fig. 8
figure 8

Association analysis of TaAL1-B gene haplotypes with drought-related traits in wheat. (A) TaAL1-B gene structure and SNP sites of the two haplotypes. (B) Distribution of transcription factor binding site in two haplotype promoters of TaAL1-B gene. (C) Association analysis of two haplotypes of TaAL1-B with survival rate under drought stress. The values show the mean ± SD. Significant statistical analysis was carried out by Student’s t -test (*p < 0.05; **p < 0.01)

Discussion

Alfin-like transcription factors play crucial roles in plant growth, development and stress responses [26]. In this study, 27 TaAL genes were identified and comprehensively analyzed the roles of Alfin-like transcription factors in wheat (Fig. 1 and Table S1). The phylogenetic analysis showed that AL genes were classified into AL-A, AL-B and AL-C subfamilies, and TaAL genes belonged to AL-B and AL-C members (Fig. 1). AL-A subfamily genes only existed in ancient plant species, and AL-B members occurred after differentiation of the monocot and dicot, which was consistent with previous study [5]. Most TaAL genes undergone block duplicate events, which played vital function for expansion of TaAL gene family (Table S2). This was consistent with results of the AL gene family in Malus domestica, Arabidopsis lyrata, Arabidopsis thaliana, and Thellungiella halophila [6, 27].

Previous studies showed that the AL genes participated in regulating plant growth, development and stress responses [5, 6, 16, 28]. In order to investigate the function of TaALs, we analyzed the cis-acting elements in the promoter of TaALs, the stress responsive-elements, hormone responsive-elements, light responsive-elements, and growth and development related-elements were identified (Fig. 4). In light responsive-element category, G-box had the largest number, suggesting the interaction of G-box with TFs facilitated the involvement of TaALs in light signaling. Similarly, TaALs might participate in ABA and MeJA signaling through the interaction of ABRE and CGTCA-motif with TFs. AhAL1 was up-regulated after ABA treatment and could enhance the ABA-mediated stomatal closure [7]. In addition, AL genes were involved in root hair elongation and seed germination [7, 26], which might be regulated by hormone responsive- and growth and development related-elements. Most TaAL promoter included MBS, ABRE and CGTCA-motif (Fig. 4). MBS and ABRE elements were related to drought stress response [29, 30]. Thus, the expression level of TaAL genes were studied under drought and salt stresses. The transcriptome data revealed up-regulation of TaAL1-A/B/D, TaAL6-B/D and TaAL9-B under drought and salt stresses, suggesting their potential crosstalk between drought and salt stresses. The real-time PCR results showed TaAL1-B was significantly up-regulated in response to drought stress (Fig. 5 and Table S3).

The interacting protein and cluster analysis showed that TaALs were most likely to interact with TRAUCO and ubiquitin-40 S ribosomal protein under drought and salt stresses (Fig. 6). TRAUCO protein was a core component of histone methyltransferase complex [31], and AL was a histone reader to bind to histone mark H3K4me2/3 [18, 19]. In Arabidopsis, AL2-PRC1 complexes promoted seed germination through H3K4me3-to-H3K27me3 chromatin state switch in repression of seed developmental genes [20]. Therefore, the interaction of TaAL with TRAUCO protein might regulate transcription initiation of target genes by changing chromatin state to respond drought and salt stresses. TaAL might participate in protein degradation by interacting with ubiquitin-40 S ribosomal protein. The downstream target genes and cluster analysis showed that TaAL might regulate the expression of chromatin remodeling protein, actin-related protein, E3 ubiquitin-protein ligase, upstream activation factor subunit, and bZIP transcription factor gene to respond drought or salt stresses (Fig. 7 and Fig. S4). These target genes played important roles in response to drought or salt stresses [32,33,34,35].

Association analysis of haplotypes with drought-related traits showed that Hap-I allelic variation of TaAL1-B had significantly higher survival rate compared to Hap-II under drought stress, probably due to the difference of transcription factor binding site in the promoter of TaAL1-B haplotypes lead to the difference of wheat drought tolerance (Fig. 8B and Table S9). FAR1 and MYB bound to TaAL1-B-Hap-I and Hap-II promoters, respectively (Fig. 8). In Arabidopsis, FAR1 could be up-regulated by ABA and drought stresses [36]. The far1 mutants had wider stomata, lose water faster, and were more sensitive to drought than the wild type plants [36]. Overexpression of HvFRF9 significantly enhanced the drought tolerance in Arabidopsis by increasing the absorption and transportation of water and the activity of antioxidant enzymes [37]. FAR1 negatively controlled ROS levels by directly or indirectly regulating genes involved in ROS homeostasis [38]. Therefore, the up-regulation of FAR1 might regulate the expression of TaAL1-B-Hap-I to enhance drought resistance by regulating ROS homeostasis. Furthermore, FAR1 was not only involved in light signaling, but also the positive regulator of ABA signaling, which was a hub for the integration of light and ABA signaling pathway [36, 38]. FAR1 might also enhance drought resistance via ABA signaling pathway [38]. Similarly, the SNPs were identified in the promoter region of OsAL7.1, which were significantly associated with grain yield, drought coefficient and seed width [5]. These results lay a foundation for further functional studies of Alfin-like transcription factors under drought and salt stresses.

Materials and methods

Genome-wide identification of Alfin-like family genes

The candidate Alfin-like protein sequences were obtained by using the PFAM ID of Alfin domain (PF12165) [39] to search against wheat protein database in WheatOmics 1.0 website [40]. Then, the InterPro [39] and SMART (Simple Modular Architecture Research Tool) [41] database were used to confirm the candidate Alfin-like proteins. The physiological and biochemical parameters of the Alfin-like proteins were analyzed by WheatOmics 1.0 [40]. The subcellular localization of the Alfin-like proteins were predicted using Plant-mPLoc [42].

Multiple sequence alignment and phylogenetic tree construction

Multiple sequence alignment of Alfin-like amino acid sequences was performed with ClustalW using the default options in MEGA 11 [43] and visualized by ESPript 3.0 [44]. The phylogenetic tree was constructed by using the neighbor-joining (NJ) method with 1000 bootstrap replicates in MEGA 11 software [43] and visualized by Evolview service [45].

Gene structures, conserved motifs and domains analyses

The exon-intron structures of Alfin-like genes were analyzed based on TGT (Triticeae-Gene Tribe) [46]. The conserved motifs and domains of Alfin-like family proteins were annotated using the MEME program [47] and SMART database [41], and visualized by the TBtools [48].

Chromosome localization, synteny and Ka/Ks analyses

The chromosome localization and paralogous gene pairs of Alfin-like genes were identified by using TGT (Triticeae-Gene Tribe) [46]. The gene duplication events were determined by PLAZA database [49]. TBtools was used to calculate the Ka (non-synonymous rate), Ks (synonymous rate), and Ka/Ks (the nonsynonymous and synonymous substitution ratio) values of the paralogous gene pair with the Nei-Gojobori (NG) method [48].

Analysis of cis-acting elements of TaAL promoters

The 1500 bp promoter sequences of TaAL genes were obtained from the wheat genome database [40]. PlantCARE database (https://bioinformatics.psb.ugent.be/webtools/plantcare/html/) was used to identify and count the cis-acting elements on the TaAL genes promoter.

The interacting proteins and downstream target genes prediction

Protein–protein interaction (PPI) was analyzed using the STRING [24] and plant.MAP [25] database. To predict the downstream target genes of TaALs, WheatCENet database [50] was used to obtain the co-expression genes associated with TaALs (PCC > = 0.9), and the co-expression genes that including G-rich motifs (GNGGTG/GTGGNG) [15] on 1.5 kb promoter sequences were potential downstream target genes of TaALs.

Transcriptome data analysis

The transcriptome data of NCBI SRA (SRX2948720, SRX2948721, SRX2948729, SRX2948726, SRX2948731, SRX2948728, SRX2948715, SRX2948717, and SRX2948727) were used to analyze wheat gene expression profiles under drought stress. NCBI SRA data (SRX1162592, SRX1162594, SRX1162596, and SRX1162598) were obtained to analyze wheat gene expression profiles under salt stress.

RNA extraction and real-time PCR

Real-time PCR was performed to detect the expression patterns of TaAL genes according to previous study [51]. The total RNA was isolated using RNApure Plant Kit (Vazyme), and the first-strand cDNA was synthesized from 1 Âµg of total RNA using Prescript III RT ProMix (CISTRO). The real-time PCR was performed using gene-specific primers (Table S10) with 2× Ultra SYBR Green qPCR Mix (CISTRO), and the TaActin gene was selected as a reference control. The real-time PCR cycling parameters were 95 ℃ for 30 s, followed by 45 cycles at 95 ℃ for 5 s and 60 ℃ for 30 s, with a melting curve analysis. All reactions were performed in three technical and biological replicates. The relative expression levels of target genes were calculated using the 2−△△CT method [52].

Association analysis of TaAL1-B gene haplotypes with drought-related traits in wheat

The 681 wheat resequencing data in Wheat Union database [22] (http://wheat.cau.edu.cn/WheatUnion/) was used to analyze the haplotypes of TaAL1-B gene. PlantRegMap database [53] (https://plantregmap.gao-lab.org/index.php) was used to identify transcription factor binding site in the 1500 bp promoter sequences of TaAL1-B. The phenotypic data of 39 wheat survival rate under drought stress were obtained from previous study [23].

Data availability

No datasets were generated or analysed during the current study.

Abbreviations

AL:

Alfin-like

CDS:

Coding sequence

HMM:

Hidden Markov Model

MBS:

MYB binding site

ABRE:

ABA-responsive element

CGTCA-Motif:

MeJA-responsive element

References

  1. Gupta PK, Balyan HS, Sharma S, Kumar R. Genetics of yield, abiotic stress tolerance and biofortification in wheat (Triticum aestivum L). Theor Appl Genet. 2020;133(5):1569–602.

    Article  PubMed  Google Scholar 

  2. Hirayama T, Shinozaki K. Research on plant abiotic stress responses in the post-genome era: past, present and future. Plant J. 2010;61(6):1041–52.

    Article  CAS  PubMed  Google Scholar 

  3. Lai X, Chahtane H, Martin-Arevalillo R, Zubieta C, Parcy F. Contrasted evolutionary trajectories of plant transcription factors. Curr Opin Plant Biol. 2020;54:101–7.

    Article  CAS  PubMed  Google Scholar 

  4. Das Gupta M, Tsiantis M. Gene networks and the evolution of plant morphology. Curr Opin Plant Biol. 2018;45(Pt A):82–7.

    Article  CAS  PubMed  Google Scholar 

  5. Yang Y, Ma X, Xia H, Wang L, Chen S, Xu K, Yang F, Zou Y, Wang Y, Zhu J, et al. Natural variation of Alfin-like family affects seed size and drought tolerance in rice. Plant J. 2022;112(5):1176–93.

    Article  CAS  PubMed  Google Scholar 

  6. Wang P, Lu S, Li W, Ma Z, Mao J, Chen B. Genome-wide characterization of Alfin-Like (AL) genes in apple and functional identification of MdAL4 in response to drought stress. Plant Cell Rep. 2023;42(2):395–408.

    CAS  PubMed  Google Scholar 

  7. Tao JJ, Wei W, Pan WJ, Lu L, Li QT, Ma JB, Zhang WK, Ma B, Chen SY, Zhang JS. An Alfin-like gene from Atriplex hortensis enhances salt and drought tolerance and abscisic acid response in transgenic Arabidopsis. Sci Rep. 2018;8(1):2707.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Mellor J. It takes a PHD to read the histone code. Cell. 2006;126(1):22–4.

    Article  CAS  PubMed  Google Scholar 

  9. Wysocka J, Swigut T, Xiao H, Milne TA, Kwon SY, Landry J, Kauer M, Tackett AJ, Chait BT, Badenhorst P, et al. A PHD finger of NURF couples histone H3 lysine 4 trimethylation with chromatin remodelling. Nature. 2006;442(7098):86–90.

    Article  CAS  PubMed  Google Scholar 

  10. Wei W, Zhang YQ, Tao JJ, Chen HW, Li QT, Zhang WK, Ma B, Lin Q, Zhang JS, Chen SY. The Alfin-like homeodomain finger protein AL5 suppresses multiple negative factors to confer abiotic stress tolerance in Arabidopsis. Plant J. 2015;81(6):871–83.

    Article  CAS  PubMed  Google Scholar 

  11. Lee WY, Lee D, Chung WI, Kwon CS. Arabidopsis ING and Alfin1-like protein families localize to the nucleus and bind to H3K4me3/2 via plant homeodomain fingers. Plant J. 2009;58(3):511–24.

    Article  CAS  PubMed  Google Scholar 

  12. Yan C, Yang N, Li R, Wang X, Xu Y, Zhang C, Wang X, Wang Y. Alfin-like transcription factor VqAL4 regulates a stilbene synthase to enhance powdery mildew resistance in grapevine. Mol Plant Pathol. 2023;24(2):123–41.

    Article  CAS  PubMed  Google Scholar 

  13. Krochko JE, Pramanik SK, Bewley JD. Contrasting storage protein synthesis and Messenger RNA Accumulation during development of zygotic and somatic embryos of Alfalfa (Medicago sativa L). Plant Physiol. 1992;99(1):46–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Winicov I. Alfin1 transcription factor overexpression enhances plant root growth under normal and saline conditions and improves salt tolerance in alfalfa. Planta. 2000;210(3):416–22.

    Article  CAS  PubMed  Google Scholar 

  15. Winicov II, Bastola DR. Transgenic overexpression of the transcription factor alfin1 enhances expression of the endogenous MsPRP2 gene in alfalfa and improves salinity tolerance of the plants. Plant Physiol. 1999;120(2):473–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Zhu Y, Chen Q, Liu X, Qu Y. The Conringia planisiliqua Alfin-like2 gene enhances drought and salt tolerance in Arabidopsis thaliana. Theoretical Experimental Plant Physiol. 2021;33(4):427–41.

    Article  CAS  Google Scholar 

  17. Wei W, Huang J, Hao YJ, Zou HF, Wang HW, Zhao JY, Liu XY, Zhang WK, Ma B, Zhang JS, et al. Soybean GmPHD-type transcription regulators improve stress tolerance in transgenic Arabidopsis plants. PLoS ONE. 2009;4(9):e7209.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Peng L, Wang L, Zhang Y, Dong A, Shen WH, Huang Y. Structural analysis of the Arabidopsis AL2-PAL and PRC1 complex provides mechanistic insight into active-to-repressive chromatin state switch. J Mol Biol. 2018;430(21):4245–59.

    Article  CAS  PubMed  Google Scholar 

  19. Sureshkumar S, Bandaranayake C, Lv J, Dent CI, Bhagat PK, Mukherjee S, Sarwade R, Atri C, York HM, Tamizhselvan P, et al. SUMO protease FUG1, histone reader AL3 and chromodomain protein LHP1 are integral to repeat expansion-induced gene silencing in Arabidopsis thaliana. Nat Plants. 2024;10(5):749–59.

    Article  CAS  PubMed  Google Scholar 

  20. Molitor AM, Bu Z, Yu Y, Shen WH. Arabidopsis AL PHD-PRC1 complexes promote seed germination through H3K4me3-to-H3K27me3 chromatin state switch in repression of seed developmental genes. PLoS Genet. 2014;10(1):e1004091.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Kayum MA, Park JI, Ahmed NU, Jung HJ, Saha G, Kang JG, Nou IS. Characterization and stress-induced expression analysis of Alfin-like transcription factors in Brassica rapa. Mol Genet Genomics. 2015;290(4):1299–311.

    Article  PubMed  Google Scholar 

  22. Wang W, Wang Z, Li X, Ni Z, Hu Z, Xin M, Peng H, Yao Y, Sun Q, Guo W. SnpHub: an easy-to-set-up web server framework for exploring large-scale genomic variation data in the post-genomic era with applications in wheat. Gigascience 2020, 9(6).

  23. Mei F, Chen B, Du L, Li S, Zhu D, Chen N, Zhang Y, Li F, Wang Z, Cheng X, et al. A gain-of-function allele of a DREB transcription factor gene ameliorates drought tolerance in wheat. Plant Cell. 2022;34(11):4472–94.

    Article  PubMed  PubMed Central  Google Scholar 

  24. von Mering C, Huynen M, Jaeggi D, Schmidt S, Bork P, Snel B. STRING: a database of predicted functional associations between proteins. Nucleic Acids Res. 2003;31(1):258–61.

    Article  Google Scholar 

  25. McWhite CD, Papoulas O, Drew K, Cox RM, June V, Dong OX, Kwon T, Wan C, Salmi ML, Roux SJ, et al. A pan-plant protein Complex Map reveals Deep Conservation and Novel assemblies. Cell. 2020;181(2):460–74. e414.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Jin R, Yang H, Muhammad T, Li X, Tuerdiyusufu D, Wang B, Wang J. Involvement of Alfin-Like Transcription Factors in Plant Development and stress response. Genes (Basel) 2024, 15(2).

  27. Song Y, Gao J, Yang F, Kua CS, Liu J, Cannon CH. Molecular evolutionary analysis of the Alfin-like protein family in Arabidopsis lyrata, Arabidopsis thaliana, and Thellungiella Halophila. PLoS ONE. 2013;8(7):e66838.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Chandrika NNP, Sundaravelpandian K, Yu SM, Schmidt W. ALFIN-LIKE 6 is involved in root hair elongation during phosphate deficiency in Arabidopsis. New Phytol. 2013;198(3):709–20.

    Article  CAS  PubMed  Google Scholar 

  29. Yu H, Liu B, Yang Q, Yang Q, Li W, Fu F. Maize ZmLAZ1-3 gene negatively regulates drought tolerance in transgenic Arabidopsis. BMC Plant Biol. 2024;24(1):246.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Li S, Lin YJ, Wang P, Zhang B, Li M, Chen S, Shi R, Tunlaya-Anukit S, Liu X, Wang Z, et al. The AREB1 transcription factor influences histone acetylation to regulate Drought responses and tolerance in Populus trichocarpa. Plant Cell. 2019;31(3):663–86.

    Article  CAS  PubMed  Google Scholar 

  31. Aquea F, Johnston AJ, Canon P, Grossniklaus U, Arce-Johnson P. TRAUCO, a trithorax-group gene homologue, is required for early embryogenesis in Arabidopsis thaliana. J Exp Bot. 2010;61(4):1215–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Ali A, Yun DJ. Chromatin remodeling complex HDA9-PWR-ABI4 epigenetically regulates drought stress response in plants. Plant Signal Behav. 2020;15(10):1803568.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Do BH, Hiep NT, Lao TD, Nguyen NH. Loss-of-function mutation of ACTIN-RELATED PROTEIN 6 (ARP6) impairs Root Growth in response to salinity stress. Mol Biotechnol. 2023;65(9):1414–20.

    Article  CAS  PubMed  Google Scholar 

  34. Li S, Zhang Y, Liu Y, Zhang P, Wang X, Chen B, Ding L, Nie Y, Li F, Ma Z, et al. The E3 ligase TaGW2 mediates transcription factor TaARR12 degradation to promote drought resistance in wheat. Plant Cell. 2024;36(3):605–25.

    Article  PubMed  Google Scholar 

  35. Yang Y, Yu TF, Ma J, Chen J, Zhou YB, Chen M, Ma YZ, Wei WL, Xu ZS. The soybean bZIP transcription factor gene GmbZIP2 confers Drought and Salt resistances in transgenic plants. Int J Mol Sci 2020, 21(2).

  36. Tang W, Ji Q, Huang Y, Jiang Z, Bao M, Wang H, Lin R. FAR-RED ELONGATED HYPOCOTYL3 and FAR-RED IMPAIRED RESPONSE1 transcription factors integrate light and abscisic acid signaling in Arabidopsis. Plant Physiol. 2013;163(2):857–66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. He X, He Y, Dong Y, Gao Y, Sun X, Chen W, Xu X, Su C, Lv Y, Ren B, et al. Genome-wide analysis of FRF gene family and functional identification of HvFRF9 under drought stress in barley. Front Plant Sci. 2024;15:1347842.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Wang H, Wang H. Multifaceted roles of FHY3 and FAR1 in light signaling and beyond. Trends Plant Sci. 2015;20(7):453–61.

    Article  CAS  PubMed  Google Scholar 

  39. Mistry J, Chuguransky S, Williams L, Qureshi M, Salazar GA, Sonnhammer ELL, Tosatto SCE, Paladin L, Raj S, Richardson LJ, et al. Pfam: the protein families database in 2021. Nucleic Acids Res. 2021;49(D1):D412–9.

    Article  CAS  PubMed  Google Scholar 

  40. Ma S, Wang M, Wu J, Guo W, Chen Y, Li G, Wang Y, Shi W, Xia G, Fu D, et al. WheatOmics: a platform combining multiple omics data to accelerate functional genomics studies in wheat. Mol Plant. 2021;14(12):1965–8.

    Article  CAS  PubMed  Google Scholar 

  41. Letunic I, Khedkar S, Bork P. SMART: recent updates, new developments and status in 2020. Nucleic Acids Res. 2021;49(D1):D458–60.

    Article  CAS  PubMed  Google Scholar 

  42. Chou KC, Shen HB. Plant-mPLoc: a top-down strategy to augment the power for predicting plant protein subcellular localization. PLoS ONE. 2010;5(6):e11335.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Tamura K, Stecher G, Kumar S. MEGA11: Molecular Evolutionary Genetics Analysis Version 11. Mol Biol Evol. 2021;38(7):3022–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Gouet P, Robert X, Courcelle E. ESPript/ENDscript: extracting and rendering sequence and 3D information from atomic structures of proteins. Nucleic Acids Res. 2003;31(13):3320–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Subramanian B, Gao S, Lercher MJ, Hu S, Chen WH. Evolview v3: a webserver for visualization, annotation, and management of phylogenetic trees. Nucleic Acids Res. 2019;47(W1):W270–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Chen Y, Song W, Xie X, Wang Z, Guan P, Peng H, Jiao Y, Ni Z, Sun Q, Guo W. A collinearity-incorporating homology inference strategy for connecting emerging assemblies in the Triticeae Tribe as a pilot practice in the plant pangenomic era. Mol Plant. 2020;13(12):1694–708.

    Article  CAS  PubMed  Google Scholar 

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

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

  49. Van Bel M, Diels T, Vancaester E, Kreft L, Botzki A, Van de Peer Y, Coppens F, Vandepoele K. PLAZA 4.0: an integrative resource for functional, evolutionary and comparative plant genomics. Nucleic Acids Res. 2018;46(D1):D1190–6.

    Article  PubMed  Google Scholar 

  50. Li Z, Hu Y, Ma X, Da L, She J, Liu Y, Yi X, Cao Y, Xu W, Jiao Y et al. WheatCENet: a database for comparative co-expression networks analysis of allohexaploid wheat and its progenitors. Genomics Proteom Bioinf 2022.

  51. Liu H, Yang W, Zhao X, Kang G, Li N, Xu H. Genome-wide analysis and functional characterization of CHYR gene family associated with abiotic stress tolerance in bread wheat (Triticum aestivum L). BMC Plant Biol. 2022;22(1):204.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2–∆∆CT method. Methods. 2001;25(4):402–8.

    Article  CAS  PubMed  Google Scholar 

  53. Tian F, Yang DC, Meng YQ, Jin J, Gao G. PlantRegMap: charting functional regulatory maps in plants. Nucleic Acids Res 2020, 48(D1):D1104-D1113.

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by Supported by the Start-up Funding for Scientific Research of Ludong University (20240065), Innovation Project for graduate students of Ludong University (IPGS2024-099), PhD Research Startup Foundation of Henan University of Science and Technology (13480094), Natural Science Foundation of Henan Province (242300421320), Natural Science Talent Project of Weinan Normal University (2022RC25), and National Natural Science Foundation of China (32272159). The funding bodies played no role in not only the design of the study but also the collection and analysis of data and manuscript writing.

Author information

Authors and Affiliations

Authors

Contributions

Designed the experiment: H.L. Performed the experiments and analyzed data: H.L., W.L. and Z.W. Contributed reagents/materials/analysis tools: H.L., N.L., Y.X. and Y.Z. Drafted the manuscript.: H.L. Revised the manuscript: H.L., N.L. Y.X. and Y.Z. All authors read and approved the manuscript.

Corresponding authors

Correspondence to Na Li, Yongfeng Xie or Yanhong Zhao.

Ethics declarations

Ethics approval and consent to participate

The experiments did not involve endangered or protected species. The data collection of plants was carried out with permission of related institution, and complied with national or 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.

Electronic supplementary material

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Liu, H., Liu, W., Wang, Z. et al. Comprehensive analysis of Alfin-like transcription factors associated with drought and salt stresses in wheat (Triticum aestivum L.). BMC Genomics 25, 701 (2024). https://doi.org/10.1186/s12864-024-10557-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-024-10557-y

Keywords