Skip to main content

Genome-wide identification and expression analysis of the xyloglucan endotransglucosylase/hydrolase gene family in poplar

Abstract

Background

Xyloglucan endotransglucosylase/hydrolase (XTH) family plays an important role in cell wall reconstruction and stress resistance in plants. However, the detailed characteristics of XTH family genes and their expression pattern under salt stress have not been reported in poplar.

Results

In this study, a total of 43 PtrXTH genes were identified from Populus simonii × Populus nigra, and most of them contain two conserved structures (Glyco_hydro_16 and XET_C domain). The promoters of the PtrXTH genes contain mutiple cis-acting elements related to growth and development and stress responses. Collinearity analysis revealed that the XTH genes from poplar has an evolutionary relationship with other six species, including Eucalyptus robusta, Solanum lycopersicum, Glycine max, Arabidopsis, Zea mays and Oryza sativa. Based on RNA-Seq analysis, the PtrXTH genes have different expression patterns in the roots, stems and leaves, and many of them are highly expressed in the roots. In addition, there are11 differentially expressed PtrXTH genes in the roots, 9 in the stems, and 7 in the leaves under salt stress. In addition, the accuracy of RNA-Seq results was verified by RT-qPCR.

Conclusion

All the results indicated that XTH family genes may play an important role in tissue specificity and salt stress response. This study will lay a theoretical foundation for further study on molecular function of XTH genes in poplar.

Peer Review reports

Background

As the external supporting structure of cells, cell wall largely determines the shape and size of cells in the process of plant growth and development. The main components of plant cell wall are cellulose, lignin, hemicellulose and pectin [1]. Xyloglucan is the most important hemicellulose in the primary cell wall of dicotyledons [2]. Xyloglucan endotransglucosylase/hydrolase (XTH) is widely present in plant cells, which can catalyze the cleavage and polymerization of xyloglucan molecules, and achieve cell remodeling by regulating the elasticity and ductility of cells [3]. It is a key enzyme in the process of cell wall reconstruction in plants.

XTH family genes contain a typical catalyze enzymatic reaction motif (HDEIDFEFLG), in which the first glutamic acid residue (E) is affinity site and the second is proton donor [4, 5]. XTHs can be used to form a covalent glycosylase intermediate, which can be decomposed by water to carry out hydrolysis reaction, or enter the sugar substrate to produce transglycosylation [3]. According to the conserved motifs, the XTH family can be divided into three groups (I, II, and III), and the group III is divided into IIIA and IIIB group [6]. It has been found that the XTHs in groups I, II and IIIB have significant xyloglucan endohydrolase (XET) activity, while the XTHs in IIIA show xyloglucan endoglucosidase (XEH) activity [7]. Moreover, a small part of outlier group (early diverting) was found in detailed clade [7, 8]. In summary, XTH family is divided into I / II, IIIA, IIIB and the early diverting groups [9].

The expression characteristics of many XTH family genes are different in different tissues during different developmental stages. For example, poplar PtxtXET16–34 was specifically expressed in the developing wood, and overexpression of the gene could increase xyloglucan content in the xylem of the primary wall and promote the growth of vessel element [10]. In Arabidopsis, AtXTH17, AtXTH18, AtXTH19 and AtXTH20 were specifically expressed in roots and play an important role in root elongation and root hair initiation [11, 12]. AtXTH21 participates in the growth of primary roots by changing the cellulose deposition and the extensibility in the cell wall [13]. AtXTH31 can encode genes of root extension region expression and participates in root elongation and growth [14]. In addition, some XTH genes are active in the development of leaf [15], fruit [16] and flower [17]. Many studies have also indicated that the expression of XTHs was regulated by plant hormones. For example, the expression of rice OsXTH8 [18] and Arabidopsis AtXTH21 [13] was up-regulated by gibberellin treatment. The expression of banana MA-XET1 was induced by ethylene and participates in the ripening and softening of peel and pulp [19]. Under the treatment of brassinosteroid (BR), the expression of AtXTH22 and AtXTH24 genes increased significantly in Arabidopsis, thereby promoting the elongation of cell walls [20].

In particular, XTH genes play an important role in abiotic stress in plants. The content of xyloglucan was reduced in AtXTH31 mutant of Arabidopsis, which leads to a decrease in the content of absorbed Al3+, thereby improving the aluminum stress tolerance [14]. Similarly, Arabidopsis xth17 and xth15 mutants had higher aluminum tolerance compared with wild-type plants [21, 22]. CaXTH3 gene in pepper can be induced by drought, high salt and cold, and overexpression of the gene leads to severe fold of Arabidopsis leaves, however, which improves the drought and salt tolerance of transgenic plants [23]. Under low temperature, the expression of persimmon DkXTH6 gene was inhibited, while DkXTH7 showed a high level of transcription [24]. Overexpression of DkXTH1 enhanced the tolerance of transgenic Arabidopsis to salt, ABA and drought stress, and affected the growth of roots and leaves [25].

The first XTH gene was discovered in cowpea [26], then the XTH family genes were established in other species. Among them, a total of 33, 29, 35, and 25 XTH genes were identified in Arabidopsis [11], rice [27], sorghum [28], and tomato [16], respectively. In poplar, Geisler-Lee [29] found 41 XTH genes, and Ye [30] found 42 genes. However, they both did not systematically analyze their expression pattern under salt stress in poplar. In the present research, we found 43 XTHs family genes from poplar. And we determined the evolutionary relationship of poplar XTH family by the analysis of gene structure and homologous collinearity. Meanwhile, the gene expression pattern of XTH family genes was profiled among different tissues under salt stress. This study provides a theoretical foundation for the research on salt-tolerance improvement of poplar plants.

Result

Phylogenetic tree analysis of XTH family

In this study, 43 XTH members were found in poplar genome and named as PtrXTH1 - PtrXTH43 according to their position in the chromosome. The phylogenetic tree of XTH family genes from poplar and Arabidopsis was constructed by MEGA6. As shown in the Fig. 1, the XTH family could be divided into 4 groups (GroupI/II, Group IIIA, Group IIIB and Early diverging group), which was similar to that in Brassica rapa [9]. And the GroupI/IIcontained the most members with 33 members. Group IIIA had 6 members, Group 3B had 3 members, and group 4 had only one member (Fig. 2a).

Fig. 1
figure 1

Phylogenetic analysis of the XTH proteins from poplar and Arabidopsis. The phylogenetic tree was constructed by MEGA6 software, with neighbor-joining method. The evolutionary tree was divided to 4 groups, each color represents one group

Fig. 2
figure 2

Characterization of the genes structure, conserved protein domain and motif of the poplar XTH genes. (a) Phylogenetic tree of poplar XTH genes family, (b) Gene structures, (c) Conserved protein domain, (d) Protein conserved motifs

The physicochemical properties of XTH proteins were also analyzed (Supplementary Table 1). PthXTH2 encoded the largest amino acids (351) with the molecular weight (MW) of 40.77 kDa, while PthXTH11 had the smallest amino acids (117) with the MW of 13.44 kDa. The theoretical isoelectric point (PI) ranged from 4.53 to 9.67, while the instability coefficient ranged from 27.44 to 55.83. In addition, only one protein had grand average of hydropathicity (GRAVY) greater than 0, which was regarded as hydrophobic protein, and the others were less than 0, which were expressed as hydrophilic protein.

For subcellular localization, there were 8 genes predicted to be located in the extracellular, and the other 35 proteins in the plasma membrane, and 80% of XTH proteins had signal peptides. Then we randomly selected PtrXTH2, PtrXTH14 and PtrXTH37 genes to verify the predicted results (Fig. 3). The signal of positive control (35S-GFP) was detected in the Nicotiana benthamiana leaves nucleus and plasma membrane, while 35S-XTH2-GFP, 35S-XTH14-GFP and 35S-XTH37-GFP were detected in the plasma membrane. The results verified the accuracy of the subcellular localization prediction of XTH proteins.

Fig. 3
figure 3

Subcellular localization of PtrXTH2, PtrXTH14 and PtrXTH37 proteins. The control (35S-GFP) and fusion vector (35S-XTH2-GFP, 35S-XTH14-GFP, 35S-XTH37-GFP) were separately transiently expressed in N. benthamiana leaves by agrobacterium-mediated method. Scale bar = 20 um

Gene structure analysis of poplar XTH family genes

In general, the genes that are clustered in one group have a similar structure. In this study, the genes in GroupI/IIhave 2 to 3 introns except PtrXTH11 (1 intron) and PtrXTH1 (4 introns), and the genes in Group IIIA have 3 introns except PtrXTH5 (1 intron). In addition, there were 3 introns in the genes in Group IIIB and Early diverging group, respectively (Fig. 2b). We used PFAM data to analyze the conserved domain of the 43 XTHs. Only three proteins including PtrXTH11, PtrXTH43 and PtrXTH5 had one conserved domain (Glyco_hydro_16), while the other 40 genes had both Glyco_hydro_16 and XET_C domain (Fig. 2c).

In addition, we found 10 related motifs through MEME online software (Fig. 2d and Supplementary Table 2). The longest one contained 50 amino acids (motif 3 and motif 4), and the shortest contained 7 amino acids (motif 7) (Fig. 2d). Glyco_hydro_16 contained the motifs 1, 4, 6, 7 and 8, but not all protein had the 5 motifs. For instance, XET_C contained the motif 5 and motif 10. The results were similar to those in the previous study [9].

Cis-acting elements in PtrXTHs

In this study, the 2000 bp upsteam sequences of the 43 PtrXTHs were analysed by PlantCARE. As shown in Fig. 4 and Supplementary Table 3, we found 16 types of cis-elements in the upstream promoter sequences of PtrXTHs. Among them, there were many elements related to abiotic stresses such as drought-inducibility, defense and stress responsiveness, low-temperature responsiveness, and wound-response. The elements including meristem expression, analytical induction, circadian control, seed-specific regulation, cell cycle regulation and light responsive were related to plant growth and development. The elements such as auxin responsiveness, MeJA-responsiveness, gibberellin-responsive, abscisic acid responsiveness and zeatin metabolism regulation were related to hormones response. In addition, there were multiple elements in response to a certain stress (Supplementary Table 3). For example, I-box, AE-box, and G-Box were involved in light responsive. P-box, TATC-box and GARE-motif were involved in gibberellin-responsive. The results showed that many PtrXTHs were involved in plant growth and stress response.

Fig. 4
figure 4

Cis-elements analysis of the poplar XTH gene promoters. The left side represents the phylogenetic tree of the poplar XTH genes, and the rights with different colored ellipses represent different elements

Chromosomal distribution and synteny analysis of PtrXTHs

We mapped the poplar XTH family genes to poplar chromosomes through genome annotation. As shown in Fig. 5a, the 43 PtrXTHs were not centrally allocated to a certain chromosomes, while they were randomly allocated to as many chromosomes as possible in a scattered manner. For example, there was only one XTH gene allocated on chromosomes 7, 8, 10, 16, 19, and two on the chromosomes 1, 3, 4, 5, 9 and 13, and six on the chromosome 6, which had the most members. Interestingly, there was no PtrXTHs on the chromosomes 12, 15, and 17.

Fig. 5
figure 5

Gene location and repetitive events of the poplar XTH genes. (a) Chromosome location information on the of the poplar XTH genes. The blue arcs represent different pairs of tandem repeat genes. (b) Fragment duplication analysis of the poplar XTH genes. The line and heat map in the outer circle represent gene density on the chromosome, the red lines represents the fragment duplicate gene pair, and the gray lines represent the synteny blocks of the XTH genes in poplar genome

We analyzed the repetitive events of the PtrXTHs through MCscan. As shown in Fig. 5a and Supplementary Table 4, there were 6 pairs of tandem repetitive events among 43 PtrXTHs, of which PtrXTH18 shared two pairs of repetitive events. And 6 pairs of fragment duplication events were also identified in PtrXTHs (Fig. 5b), which were distributed on the different chromosomes.

In order to further study the evolutionary relationship of XTH family, we constructed systematic maps of the XTH genes between poplar and six other species, including four dicotyledonous (E. robusta, S. lycopersicum, G. max and Arabidopsis) and two monocotyledons (Z. mays and O. sativa). As shown in Fig. 6 and Supplementary Table 5, there were 8 homologous pairs in Arabidopsis, 14 pairs in G. max, 12 pairs in E. robusta, 11 pairs in S. lycopersicum, 1 pair in Z. mays, and 0 pairs in O. sativa. Among them, a total of 16 PtrXTHs showed high collinearity relationship with other 6 species (7 in Arabidopsis, 8 in G. max, 9 in E. robusta, 10 in S. lycopersicum, and 1 in Z. mays). In addition, the PtrXTHs shared homologous pairs in the different species (Fig. S1). For example, PtrXTH10 was homologous with the gene from E. robusta, S. lycopersicum and Z. mays, and PtrXTH7, PtrXTH15 and PtrXTH37 were shared in G. max, Arabidopsis and S. lycopersicum.

Fig. 6
figure 6

Collinearity analysis of the XTH genes from poplar and other six species. The blue lines represent XTH syntenic gene pairs between poplar and other species, and the gray lines represent orthologous genesof poplar and other six species

Expression patterns of PtrXTH genes in different tissues

In order to explore the expression patterns of P. simonii × P. nigra XTH genes, the expression level of PtrXTHs in the different tissues was profiled by RNA-Seq. As shown in Fig. 7 and Supplementary Table 6, the XTH genes were divided into 4 groups according to their expression patterns. Among them, the genes in the group 2, group 3 and group 4 were mainly expressed in the leaves, roots and stems, respectively. However, the expression of the genes in the group 1 was mixed. Similarly, some duplicate gene pairs had similar expression patterns. For example, PtrXTH1 and PtrXTH9 were expressed in the leaves, and PtrXTH17 and PtrXTH18 were expressed in the roots. However, some tandem gene pairs showed different expression trends. For instance, the PtrXTH21 was expressed in the leaves, while the PtrXTH24 gene was mainly expressed in the stems. The PtrXTH26 gene was mainly expressed in the roots, while the PtrXTH27 gene was mainly expressed in the leaves.

Fig. 7
figure 7

Expression pattern of the poplar XTH genes in different tissues. Red represents high expression, green represents low expression. The left represents four gene clusters, R, L, and S represents roots, leaves and stems without salt stress, respectively

PtrXTHs expression analysis under salt stress

We also analyzed gene expression patterns of PtrXTHs in different tissues under salt stress by RNA-Seq. As shown in Fig. 8 and Supplementary Table 7, there were 7 differential expression genes (DEGs) in the leaves (4 up- and 3 down-regulated), 9 in the stems (4 up- and 5 down-regulated), and 11 in the roots (5 up- and 6 down-regulated). In addition, we analyzed the shared genes among different tissues. As shown in Fig. 8d, there were 3 shared genes (0 up- and 3 down-regulated) in the roots and stems,3 shared genes (2 up- and 1 down-regulated) in the stems and leaves, and 2 shared genes (0 up- and 2 down-regulated) in the roots and leaves. Furthermore, only one shared genes (0 up- and 1 down-regulated) were shared in the three tissues.

Fig. 8
figure 8

Expression pattern of PtrXTH genes under salt stress. (a-c) Heatmaps of DEGs in leaves, stems and roots, respectivley. R, L, and S represents the roots, leaves and stems without salt stress, respectively. NR, NL, and NS represents the roots, leaves and stems under salt stress, respectively. (d) Venn diagrams of DEGs among the three tissues

To verify the accuracy of RNA-Seq data, the DEGs in the roots, stems and leaves were selected for RT-qPCR. As showed in Fig. 9, the relative expression levels of DEGs verified by RT-qPCR were in consistent with the mRNA abundance profiled by RNA-Seq analysis generally.

Fig. 9
figure 9

Gene expression levels in the three tissues by RNA-Seq and RT-qPCR. Salt and control represent with and without salt stress, respectively

Discussion

XTH family genes play an indispensable role in plants growth and development. In this study, we found 43 XTH genes from P. trichocarpa, which were named as PtrXTH1 to PtrXTH43 according to their location on the chromosomes. The 43 PtrXTHs were divided into 4 subgroups compared with Arabidopsis, which was similar with rice and soybean. Most XTH genes contain two main conserved domains (Glyco_hydro_16 and XET_C domain). However, we found that PtrXTH5 and PtrXTH11 were lack of the XET_C domain. We speculate that there was a loss of XET_C domain during the evolution of the poplar XTH genes. In addition, we predicted that most PtrXTH proteins are located in the plasma membrane, while a small are located in the extracellular, which similar with the tobacco [31]. In our study, the homeopathic transformation experiment of tobacco found that the PtrXTH2, PtrXTH14 and PtrXTH37 proteins are located in the plasma membrane, and possibly in the cell wall. In addition, we found that Arabidopsis AtXTH30 (a homologous protein of PtrXTH2) protein and soybeans GmXTH23 (a homologous protein of PtrXTH14 and PtrXTH37) protein are located in the plasma membrane [32, 33], from which we speculate that they are localized in the plasma membrane.

Plants have experienced many whole-genome replication events in the process of evolution, including tandem repeat, fragment repeat and conversion events [34, 35]. In this study, we analyzed the collinearity events of XTH family genes in poplar. Interestingly, we found that the number of tandem repeat gene pairs was same to that of fragment repeat gene pairs. We speculate that the two evolutionary mechanisms were co-regulated in the process of evolution. In addition, we found that 7 out of 12 repetitive gene pairs had different expression patterns. For example, the fragment repetitive gene pair PtrXTH10 / PtrXTH25 had same structure and conserved motif, however, they displayed obviously different expression patterns. PtrXTH10 was mainly expressed in the roots, while PtrXTH25 was mainly expressed in the leaves. Therefore, we speculate that PtrXTH genes participate in the replication process, which causes gene mutations, leading to the change of gene function and expression patterns. This phenomenon is also occurred in the ERF family [36] and NAC family [37]. In addition, the collinear relationship between PtrXTHs and other six species were analyzed. We found that PtrXTHs had more collinear gene pairs with dicotyledons plants than that with monocotyledons plants.

XTH family genes are involved in the development of plant tissues and organs, and they have expression specificity of tissue and cell. In tobacco, NtXTH4, NtXTH5, NtXTH12, and NtXTH19 had high expression levels in 19 tissues [31], and a similar situation occurred in barley [38]. In our study, we found that the PtrXTHs in groupI such as PtrXTH16 and PtrXTH31 were not expressed in the leaves, roots or stems, which may function in other organizations. The other three groups of PtrXTHs were expressed in the different tissues. However, there were more PtrXTHs expressed in the roots and stems than that in the leaves, which indirectly proves that PtrXTHs may play an important role in the development of root and stem. The previous studies have proved it. For example, Arabidopsis XTH19 and XTH23 (a homologous gene of PtrXTH19) were regulated by BES1 to participate in the development of lateral roots through brassinosteroid signaling pathway, contributing to salt adaptation of lateral roots [39]. The deficiency of XTH9 (a homologous gene of PtrXTH40) in Arabidopsis regulated the secondary wall thickening by triggering integrity signal of cell wall and stimulating the production of xylem cells [40].

During the process of growth and development, plants are subject to multiple abiotic stresses such as high temperature, salinity and drought [41]. XTH genes play an important role in abiotic stress. In our study, we analyzed the expression pattern of poplar XTH family genes under salt stress by RNA-Seq. We found 11 differential expression genes in the roots, 7 in the leaves and 9 in the stems. In addition, RT-qPCR was used to verify the accuracy of RNA-Seq analysis, which indicated the results were consistent in general. Phylogenetic analysis also showed that those genes play an important role in salt stress. For example, Arabidopsis gene XTH30, a homologous gene of PtrXTH2 and PtrXTH8, was up-regulated under salt stress, and the lack of the gene slowed down the decrease of crystalline cellulose content and microtubule depolymerization under salt stress, which had a negative impact on salt tolerance [32]. Overexpression of PeXTH from Populus euphratica (a homologous gene of PtrXTH28) could increase water holding capacity and reduce salt concentration in fleshy tissues and mesophyll cells, which improved salt tolerance of transgenic tobacco [42]. All evidences indicate that PtrXTH genes play an important role in regulating stress response in plants.

Conclusion

In this study, a total of 43 XTH genes were identified from poplar, which were divided into 4 groups. These genes were randomly distributed on the 16 chromosomes and 2 scaffolds in poplar. And as many as 12 pairs of duplication events among poplar XTH family genes were found, including 6 pairs of fragment duplication and 6 pairs of tandem duplication. In addition, we analyzed the homologous evolutionary relationship of XTH genes between poplar and other six species, which indicates the PtrXTHs had stronger homologous relationship with dicotyledonous plants, compared to monocotyledonous plants. Furthermore, we profiled the expression pattern of PtrXTHs in different tissues under salt stress through RNA-Seq and RT-qPCR. There were more PtrXTHs expressed in the roots and stems than that in the leaves, which indicates PtrXTHs may play an important role in the development of root and stem. All the results in this study provide a theoretical basis for function identification of poplar XTH genes.

Materials and methods

Identification of poplar XTH genes

All amino acid sequences of poplar XTH proteins were obtained from the genome of P. trichocarpa v3.0 in Phytozome12 (https://phytozome.jgi.doe.gov/pz/portal.html). The hidden markov models of two typical XTH family protein structures (PF00722 and PF06955) were downloaded from Pfam database (http://pfam.xfam.org/). The HMMER3.0 [43] and SMART database (http://smart.embl-heidelberg.de/) were used to search and filter all potential XTH proteins in poplar. The ExPASy website (http://web.expasy.org/protparam/) was used to calculate physical and chemical parameters of the XTH proteins. SignalP v4.1 server (http://www.cbs.dtu.dk/services/SignalP/) was used to predict the signal peptide cleavage site, and ProtComp 9.0 (http://linux1.softberry.com) was used to predict the subcellular localization of the XTH proteins.

Subcellular localization of PtrXTH2, PtrXTH14 and PtrXTH37 proteins

The full-length coding sequences of the PtrXTH2, PtrXTH14 and PtrXTH37 genes without stop codon were fused with pBI-121-GFP vector driven by 35S promoter, respectively. The fusion vector was transferred into GV3101 agrobacterium by liquid nitrogen freeze-thaw method. One-month-old N. benthamiana seedlings were used for transient transformation. The agrobacterium containing the fusion vector (35S-XTH2-GFP, 35S-XTH14-GFP, 35S-XTH37-GFP and 35S-GFP) was injected into the tobacco leaves with 1 ml syringes. Incubate in the dark for 36 h and observe tobacco cells under LMS800 laser confocal microscope.

Phylogenetic analyses of PtrXTH proteins

The protein sequences of XTH family from poplar and Arabidopsis were downloaded from the Phytozome12 (https://phytozome.jgi.doe.gov/pz/portal.html) and TAIR online websites (https://www.arabidopsis.org/), respectively. ClustalW [44] was used for multiple sequence alignment of the proteins, and MEGA6 (http://www.megasoftware.net/mega6/) was used to construct phylogenetic trees with neighbor-joining (NJ) method (bootstrap analysis for 1000 repetitions). The evolutionary tree was visualized through EvolView online software [45].

Gene structures and conserved motif analyses of PtrXTHs

The poplar genomic sequences were downloaded from the Phytozome12 database (https://phytozome.jgi.doe.gov/pz/portal.html). We used GSDS (http://gsds.cbi.pku.edu.cn/) online software to analyze the gene structure of the XTH family genes. MEME (http://meme-suite.org/tools/meme) online software was used to analyze the conserved motifs with default parameters, and the motifs were visualized with TBtools software [46].

Cis-acting element analysis of PtrXTHs

The upstream 2000 bp sequences of the XTH genes were obtained from the Phytozome12 database The cis-elements in the all sequences were predicted through PlantCRAE (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/), and the cis-elements were screened and visualized by TBtools software [46].

Chromosomal locations, synteny and duplications analyses of PtrXTHs

The genomic data of the poplar was downloaded from phytozome12, and each gene was mapped to the corresponding chromosome position through the position information of the poplar XTH genes. We used the MCscan [47] to calculate the repetitive events of the PtrXTHs family genes among different species. Dual Synteny Plotter was drew to calculate the repetitive events of PtrXTHs between poplar and the other six species (E. robusta, S. lycopersicum, G. max, Arabidopsis, Z. mays and O. sativa) by TBtools software [46].

Plant material and stress treatments

The experimental material used in this study was di-haploid P. simonii × P. nigra and N. benthamiana, which were cultured in the experimental forest of Northeast Forestry University, Harbin, China. The poplar seedlings were grown in 1/2 MS medium, and one-month-old seedlings were stressed with 0 and 150 mM NaCl for 0 and 12 h, respectively. The roots, stems and leaves at different time points were collected for RNA extraction. Meanwhile, the samples were sent to GENEWIZ Company for RNA-Seq with Illiumina platform [48]. The expression patterns of the PtrXTH genes in different tissues, and the DESeq2 [49] was used to screen the differentially expressed genes (DEGs) with absolute value of log2 (fold change) ≥1 and the adjusted p-value ≤0.05. The heatmap was drawn by TBtools [46]. In addition, we used RT-qPCR to verify the accuracy of DEGs screened by RNA-Seq. And the detailed information of RT-qPCR were referred to our previous studies [50]. All the samples were prepared with three biological replicates. All the primers were list in Supplementary Table 8.

Availability of data and materials

All data generated or analyzed during this study are included in this published article and its supplementary information files. The raw sequencing data used during the study have been deposited in NCBI SRA with the accession number SRP267437 (https://trace.ncbi.nlm.nih.gov/Traces/sra/?study=SRP267437).

Abbreviations

XTH:

Xyloglucan endotransglucosylase/hydrolase

XET:

Xyloglucan endohydrolase

XEH:

Xyloglucan endoglucosidase

MW:

Molecular weight

DEGs:

Differentially expressed genes

References

  1. Cosgrove DJ. Growth of the plant cell wall. Nat Rev Mol Cell Biol. 2005;6(11):850–61. https://doi.org/10.1038/nrm1746.

    CAS  Article  PubMed  Google Scholar 

  2. Cosgrove DJ, Jarvis MC. Comparative structure and biomechanics of plant primary and secondary cell walls. Front Plant Sci. 2012;3:204. https://doi.org/10.3389/fpls.2012.00204.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Eklof JM, Brumer H. The XTH gene family: an update on enzyme structure, function, and phylogeny in xyloglucan remodeling. Plant Physiol. 2010;153(2):456–66. https://doi.org/10.1104/pp.110.156844.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  4. Nishitani K. The role of endoxyloglucan transferase in the organization of plant cell walls. Int Rev Cytol. 1997;173:157–206. https://doi.org/10.1016/S0074-7696(08)62477-8.

    CAS  Article  PubMed  Google Scholar 

  5. Van Sandt VST, Guisez Y, Verbelen J-P, Vissenberg K. Analysis of a xyloglucan endotransglycosylase/hydrolase (XTH) from the lycopodiophyte Selaginella kraussiana suggests that XTH sequence characteristics and function are highly conserved during the evolution of vascular plants. J Exp Bot. 2006;57(12):2909–22. https://doi.org/10.1093/jxb/erl064.

    Article  PubMed  Google Scholar 

  6. Campbell P, Braam J. Xyloglucan endotransglycosylases: diversity of genes, enzymes and potential wall-modifying functions. Trends Plant Sci. 1999;4(9):361–6. https://doi.org/10.1016/S1360-1385(99)01468-5.

    CAS  Article  PubMed  Google Scholar 

  7. Baumann MJ, Eklof JM, Michel G, Kallas AM, Teeri TT, Czjzek M, et al. Structural evidence for the evolution of xyloglucanase activity from xyloglucan endo-transglycosylases: biological implications for cell wall metabolism. Plant Cell. 2007;19(6):1947–63. https://doi.org/10.1105/tpc.107.051391.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  8. Michailidis G, Argiriou A, Darzentas N, Tsaftaris A. Analysis of xyloglucan endotransglycosylase/hydrolase (XTH) genes from allotetraploid (Gossypium hirsutum) cotton and its diploid progenitors expressed during fiber elongation. J Plant Physiol. 2009;166(4):403–16. https://doi.org/10.1016/j.jplph.2008.06.013.

    CAS  Article  PubMed  Google Scholar 

  9. Wu D, Liu A, Qu X, Liang J, Song M. Genome-wide identification, and phylogenetic and expression profiling analyses, of XTH gene families in Brassica rapa L. and Brassica oleracea L. BMC Genomics. 2020;21(1):21(1). https://doi.org/10.1186/s12864-020-07153-1.

    CAS  Article  PubMed Central  Google Scholar 

  10. Nishikubo N, Takahashi J, Roos AA, Derba-Maceluch M, Piens K, Brumer H, et al. Xyloglucan endo-Transglycosylase-mediated xyloglucan rearrangements in developing wood of hybrid Aspen. Plant Physiol. 2011;155(1):399–413. https://doi.org/10.1104/pp.110.166934.

    CAS  Article  PubMed  Google Scholar 

  11. Yokoyama R, Nishitani K. A comprehensive expression analysis of all members of a gene family encoding cell-wall enzymes allowed us to predict cis-regulatory regions involved in cell-wall construction in specific organs of Arabidopsis. Plant Cell Physiol. 2001;42(10):1025–33. https://doi.org/10.1093/pcp/pce154.

    CAS  Article  PubMed  Google Scholar 

  12. Osato Y, Yokoyama R, Nishitani K. A principal role for AtXTH18 in Arabidopsis thaliana root growth: a functional analysis using RNAi plants. J Plant Res. 2006;119(2):153–62. https://doi.org/10.1007/s10265-006-0262-6.

    CAS  Article  PubMed  Google Scholar 

  13. Liu Y-B, Lu S-M, Zhang J-F, Liu S, Lu Y-T. A xyloglucan endotransglucosylase/hydrolase involves in growth of primary root and alters the deposition of cellulose in Arabidopsis. Planta. 2007;226(6):1547–60. https://doi.org/10.1007/s00425-007-0591-2.

    CAS  Article  PubMed  Google Scholar 

  14. Zhu XF, Shi YZ, Lei GJ, Fry SC, Zhang BC, Zhou YH, et al. XTH31, encoding an in vitro XEH/XET-active enzyme, regulates aluminum sensitivity by modulating in vivo XET action, Cell Wall xyloglucan content, and aluminum binding capacity in Arabidopsis. Plant Cell. 2012;24(11):4731–47. https://doi.org/10.1105/tpc.112.106039.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  15. Matsui A, Yokoyama R, Seki M, Ito T, Shinozaki K, Takahashi T, et al. AtXTH27 plays an essential role in cell wall modification during the development of tracheary elements. Plant J. 2005;42(4):525–34. https://doi.org/10.1111/j.1365-313X.2005.02395.x.

    CAS  Article  PubMed  Google Scholar 

  16. Miedes E, Lorences EP. Xyloglucan endotransglucosylase/hydrolases (XTHs) during tomato fruit growth and ripening. J Plant Physiol. 2009;166(5):489–98. https://doi.org/10.1016/j.jplph.2008.07.003.

    CAS  Article  PubMed  Google Scholar 

  17. Harada T, Torii Y, Morita S, Onodera R, Hara Y, Yokoyama R, et al. Cloning, characterization, and expression of xyloglucan endotransglucosylase/hydrolase and expansin genes associated with petal growth and development during carnation flower opening. J Exp Bot. 2011;62(2):815–23. https://doi.org/10.1093/jxb/erq319.

    CAS  Article  PubMed  Google Scholar 

  18. Jan A, Yang GX, Nakamura H, Ichikawa H, Kitano H, Matsuoka M, et al. Characterization of a xyloglucan endotransglucosylase gene that is up-regulated by gibberellin in rice. Plant Physiol. 2004;136(3):3670–81. https://doi.org/10.1104/pp.104.052274.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  19. Lu WJ, Nakano R, Kubo Y, Inaba A, Jiang YM. Cloning and expression analysis of an XET cDNA in the peel and pulp of banana fruit ripening and softening. Acta Bot Sin. 2004;46(3):355–62.

    CAS  Google Scholar 

  20. He JX, Fujioka S, Li TC, Kang SG, Seto H, Takatsuto S, et al. Sterols regulate development and gene expression in Arabidopsis. Plant Physiol. 2003;131(3):1258–69. https://doi.org/10.1104/pp.014605.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  21. Zhu XF, Lei GJ, Wang ZW, Shi YZ, Braam J, Li GX, et al. Coordination between Apoplastic and Symplastic detoxification confers plant aluminum resistance. Plant Physiol. 2013;162(4):1947–55. https://doi.org/10.1104/pp.113.219147.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  22. Zhu XF, Wan JX, Sun Y, Shi YZ, Braam J, Li GX, et al. Xyloglucan Endotransglucosylase-Hydrolase17 interacts with xyloglucan Endotransglucosylase-Hydrolase31 to confer xyloglucan Endotransglucosylase action and affect aluminum sensitivity in Arabidopsis. Plant Physiol. 2014;165(4):1566–74. https://doi.org/10.1104/pp.114.243790.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  23. Choi JY, Seo YS, Kim SJ, Kim WT, Shin JS. Constitutive expression of CaXTH3, a hot pepper xyloglucan endotransglucosylase/hydrolase, enhanced tolerance to salt and drought stresses without phenotypic defects in tomato plants (Solanum lycopersicum cv. Dotaerang). Plant Cell Rep. 2011;30(5):867–77. https://doi.org/10.1007/s00299-010-0989-3.

    CAS  Article  PubMed  Google Scholar 

  24. Han Y, Ban QY, Hou YL, Meng K, Suo JT, Rao JP. Isolation and characterization of two persimmon xyloglucan Endotransglycosylase/hydrolase (XTH) genes that have divergent functions in Cell Wall modification and fruit postharvest softening. Front Plant Sci. 2016;7:12. https://doi.org/10.3389/fpls.2016.00624.

    Article  Google Scholar 

  25. Han Y, Han S, Ban Q, He Y, Jin M, Rao J. Overexpression of persimmon DkXTH1 enhanced tolerance to abiotic stress and delayed fruit softening in transgenic plants. Plant Cell Rep. 2017;36(4):583–96. https://doi.org/10.1007/s00299-017-2105-4.

    CAS  Article  PubMed  Google Scholar 

  26. Nishitani K, Tominaga R. Endo-xyloglucan transferase, a novel class of glycosyltransferase that catalyzes transfer of a segment of xyloglucan molecule to another xyloglucan molecule. J Biol Chem. 1992;267(29):21058–64. https://doi.org/10.1016/S0021-9258(19)36797-3.

    CAS  Article  PubMed  Google Scholar 

  27. Yokoyama R, Rose JKC, Nishitani K. A surprising diversity and abundance of xyloglucan endotransglucosylase/hydrolases in rice. Classification and expression analysis. Plant Physiol. 2004;134(3):1088–99. https://doi.org/10.1104/pp.103.035261.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  28. Rai KM, Thu SW, Balasubramanian VK, Cobos CJ, Disasa T, Mendu V. Identification, Characterization, and Expression Analysis of Cell Wall Related Genes in Sorghum bicolor (L.) Moench, a Food, Fodder, and Biofuel Crop. Front Plant Sci. 2016;7:7. https://doi.org/10.3389/fpls.2016.01287.

    Article  Google Scholar 

  29. Geisler-Lee J, Geisler M, Coutinho PM, Segerman B, Nishikubo N, Takahashi J, et al. Poplar carbohydrate-active enzymes. Gene identification and expression analyses. Plant Physiol. 2006;140(3):946–62. https://doi.org/10.1104/pp.105.072652.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  30. Ye X, Yuan S, Guo H, Chen F, Tuskan GA, Cheng Z-M. Evolution and divergence in the coding and promoter regions of the Populus gene family encoding xyloglucan endotransglycosylase/hydrolases. Tree Genet Genomes. 2012;8(1):177–94. https://doi.org/10.1007/s11295-011-0431-1.

    Article  Google Scholar 

  31. Wang M, Xu Z, Ding A, Kong Y. Genome-Wide Identification and Expression Profiling Analysis of the Xyloglucan Endotransglucosylase/Hydrolase Gene Family in Tobacco (Nicotiana tabacum L.). Genes. 2018;9(6):9(6). https://doi.org/10.3390/genes9060273.

    CAS  Article  PubMed Central  Google Scholar 

  32. Yan J, Huang Y, He H, Han T, Di P, Sechet J, et al. Xyloglucan endotransglucosylase-hydrolase30 negatively affects salt tolerance in Arabidopsis. J Exp Bot. 2019;70(19):5495–506. https://doi.org/10.1093/jxb/erz311.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  33. Chen L, Zhang Y, Li Y, Li W. Cloning of GmXTH23 gene and Identification of its drought resistance. Acta Agric Univ Jiangxiensis. 2020;42(5):898–905.

    Google Scholar 

  34. Tuskan GA, DiFazio S, Jansson S, Bohlmann J, Grigoriev I, Hellsten U, et al. The genome of black cottonwood, Populus trichocarpa (Torr. & gray). Science. 2006;313(5793):1596–604. https://doi.org/10.1126/science.1128691.

    CAS  Article  PubMed  Google Scholar 

  35. Flagel LE, Wendel JF. Gene duplication and evolutionary novelty in plants. New Phytol. 2009;183(3):557–64. https://doi.org/10.1111/j.1469-8137.2009.02923.x.

    Article  PubMed  Google Scholar 

  36. Liu M, Sun W, Ma Z, Zheng T, Huang L, Wu Q, et al. Genome-wide investigation of the AP2/ERF gene family in tartary buckwheat (Fagopyum Tataricum). BMC Plant Biol. 2019;19(1):342. https://doi.org/10.1186/s12870-019-1951-3.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  37. Liu M, Ma Z, Sun W, Huang L, Wu Q, Tang Z, et al. Genome-wide analysis of the NAC transcription factor family in Tartary buckwheat (Fagopyrum tataricum). BMC Genomics. 2019;20(1):113. https://doi.org/10.1186/s12864-019-5500-0.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Fu M-M, Liu C, Wu F. Genome-Wide Identification, Characterization and Expression Analysis of Xyloglucan Endotransglucosylase/Hydrolase Genes Family in Barley (Hordeum vulgare). Molecules. 2019;24(10):1935. https://doi.org/10.3390/molecules24101935.

    CAS  Article  PubMed Central  Google Scholar 

  39. Xu PP, Fang S, Chen HY, Cai WM. The brassinosteroid-responsive xyloglucan endotransglucosylase/hydrolase 19 (XTH19) and XTH23 genes are involved in lateral root development under salt stress in Arabidopsis. Plant J. 2020;104(1):59–75. https://doi.org/10.1111/tpj.14905.

    CAS  Article  PubMed  Google Scholar 

  40. Kushwah S, Banasiak A, Nishikubo N, Derba-Maceluch M, Majda M, Endo S, et al. Arabidopsis XTH4 and XTH9 contribute to wood cell expansion and secondary wall formation (1)( OPEN ). Plant Physiol. 2020;182(4):1946–65. https://doi.org/10.1104/pp.19.01529.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  41. Ahuja I, de Vos RCH, Bones AM, Hall RD. Plant molecular stress responses face climate change. Trends Plant Sci. 2010;15(12):664–74. https://doi.org/10.1016/j.tplants.2010.08.002.

    CAS  Article  PubMed  Google Scholar 

  42. Han YS, Wang W, Sun J, Ding MQ, Zhao R, Deng SR, et al. Populus euphratica XTH overexpression enhances salinity tolerance by the development of leaf succulence in transgenic tobacco plants. J Exp Bot. 2013;64(14):4225–38. https://doi.org/10.1093/jxb/ert229.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  43. Finn RD, Clements J, Eddy SR. HMMER web server: interactive sequence similarity searching. Nucleic Acids Res. 2011;39(suppl):W29–37. https://doi.org/10.1093/nar/gkr367.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  44. Thompson JD, Gibson TJ, Higgins DG. Multiple sequence alignment using ClustalW and ClustalX. Curr Protoc Bioinformatics. 2002;Chapter 2:Unit 2.3-Unit 2.3.

    Google Scholar 

  45. Zhang H, Gao S, Lercher MJ, Hu S, Chen W-H. EvolView, an online tool for visualizing, annotating and managing phylogenetic trees. Nucleic Acids Res. 2012;40(W1):W569–72. https://doi.org/10.1093/nar/gks576.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  46. Chen C, Chen H, Zhang Y, Thomas HR, Frank MH, He Y, et al. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020;13(8):1194–202. https://doi.org/10.1016/j.molp.2020.06.009.

    CAS  Article  PubMed  Google Scholar 

  47. Wang Y, Tang H, DeBarry JD, Tan X, Li J, Wang X, et al. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012;40(7):e49. https://doi.org/10.1093/nar/gkr1293.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  48. Yao W, Wang S, Zhou B, Jiang T. Transgenic poplar overexpressing the endogenous transcription factor ERF76 gene improves salinity tolerance. Tree Physiol. 2016;36(7):896–908. https://doi.org/10.1093/treephys/tpw004.

    CAS  Article  PubMed  Google Scholar 

  49. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550. https://doi.org/10.1186/s13059-014-0550-8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  50. Zhang X, Cheng Z, Zhao K, Yao W, Sun X, Jiang T, et al. Functional characterization of poplar NAC13 gene in salt tolerance. Plant Sci. 2019;281:1–8. https://doi.org/10.1016/j.plantsci.2019.01.003.

    CAS  Article  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by the Fundamental Research Funds for the Central Universities (2572018CL03), State Key Laboratory of Tree Genetics and Breeding (Northeast Forestry University) (201205) and Applied technology research and Development Program of Heilongjiang Province (GA20B401).

Author information

Authors and Affiliations

Authors

Contributions

T J and B Z designed research. Z C conducted experiments, data analysis and wrote the manuscript. X Z, K Z, Y G and Q G performed in data analysis. W Y revised the manuscript. All authors read and approved the manuscript.

Corresponding authors

Correspondence to Boru Zhou or Tingbo Jiang.

Ethics declarations

Ethics approval and consent to participate

The di-haploid P. simonii × P. nigra and N. benthamiana were provided by the State Key Laboratory of Tree Genetics and Breeding, Northeast Forestry University, Harbin, China. The pBI-121-GFP vector was a long-term preservation material in our laboratory. The study complies with relevant institutional, national, and international guidelines and legislation.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no conflict of interest.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1: Supplementary Fig. S1.

Upset plot diagram of the poplar XTH genes throughout diverse species. The yellow color represents the number of genes that have collinearity between poplar and other species, the black circles connected by line segments represent genes that are shared by different species, and the black column represents the number of shared genes.

Additional file 2: Supplementary Table 1.

Physicochemical characterization analysis of poplar XTH protein.

Additional file 3: Supplementary Table 2.

Ten conserved motifs sequences found with MEME.

Additional file 4: Supplementary Table 3.

The list of cis-regulatoty elements of PtrXTH gene promoter.

Additional file 5: Supplementary Table 4.

Repetitive events of the PtrXTHs.

Additional file 6: Supplementary Table 5.

The homologous relationships between the poplar XTH genes and the other species.

Additional file 7: Supplementary Table 6.

Expression data of PtrXTH genes in three different tissues under salt and without salt stress.

Additional file 8: Supplementary Table 7.

The fold-changes of differentially expressed PtrXTH genes under salt stress.

Additional file 9: Supplementary Table 8.

The primers sequences used in this study.

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

Verify currency and authenticity via CrossMark

Cite this article

Cheng, Z., Zhang, X., Yao, W. et al. Genome-wide identification and expression analysis of the xyloglucan endotransglucosylase/hydrolase gene family in poplar. BMC Genomics 22, 804 (2021). https://doi.org/10.1186/s12864-021-08134-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-021-08134-8

Keywords

  • Poplar
  • XTH
  • Expression patterns
  • Salt stress