Next Article in Journal
Molecular Apomorphies in the Secondary and Tertiary Structures of Length-Variable Regions (LVRs) of 18S rRNA Shed Light on the Systematic Position of the Family Thaumastellidae (Hemiptera: Heteroptera: Pentatomoidea)
Next Article in Special Issue
PeGSTU58, a Glutathione S-Transferase from Populus euphratica, Enhances Salt and Drought Stress Tolerance in Transgenic Arabidopsis
Previous Article in Journal
A Photoactivated Ru (II) Polypyridine Complex Induced Oncotic Necrosis of A549 Cells by Activating Oxidative Phosphorylation and Inhibiting DNA Synthesis as Revealed by Quantitative Proteomics
Previous Article in Special Issue
Comprehensive Analysis of the INDETERMINATE DOMAIN (IDD) Gene Family and Their Response to Abiotic Stress in Zea mays
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 24
Issue 9
10.3390/ijms24097759
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Analysis of the Expansin Gene Family in Populus and Characterization of Expression Changes in Response to Phytohormone (Abscisic Acid) and Abiotic (Low-Temperature) Stresses

by
Zhihui Yin
,
Fangwei Zhou
,
Yingnan Chen
,
Huaitong Wu
* and
Tongming Yin
Key Laboratory for Tree Breeding and Germplasm Improvement, Southern Modern Forestry Collaborative Innovation Center, College of Forestry, Nanjing Forestry University, Nanjing 210037, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(9), 7759; https://doi.org/10.3390/ijms24097759
Submission received: 27 March 2023 / Revised: 19 April 2023 / Accepted: 21 April 2023 / Published: 24 April 2023
(This article belongs to the Special Issue Genomics-Based Tree Breeding to Improve Resilience in the Face of Global Climate Change)
Browse Figures
Versions Notes

Abstract

:
Expansins are a group of cell wall enzyme proteins that help to loosen cell walls by breaking hydrogen bonds between cellulose microfibrils and hemicellulose. Expansins are essential plant proteins that are involved in several key processes, including seed germination, the growth of pollen tubes and root hairs, fruit ripening and abscission processes. Currently, there is a lack of knowledge concerning the role of expansins in woody plants. In this study, we analyzed expansin genes using Populus genome as the study target. Thirty-six members of the expansin gene family were identified in Populus that were divided into four subfamilies (EXPA, EXPB, EXLA and EXLB). We analyzed the molecular structure, chromosome localization, evolutionary relationships and tissue specificity of these genes and investigated expression changes in responses to phytohormone and abiotic stresses of the expansin genes of Populus tremula L. (PtEXs). Molecular structure analysis revealed that each PtEX protein had several conserved motifs and all of the PtEXs genes had multiple exons. Chromosome structure analysis showed that the expansin gene family is distributed on 14 chromosomes. The PtEXs gene family expansion patterns showed segmental duplication. Transcriptome data of Populus revealed that 36 PtEXs genes were differently expressed in different tissues. Cis-element analysis showed that the PtEXs were closely associated with plant development and responses to phytohormone and abiotic stress. Quantitative real-time PCR showed that abscisic acid (ABA) and low-temperature treatment affected the expression of some PtEXs genes, suggesting that these genes are involved in responses to phytohormone and abiotic stress. This study provides a further understanding of the expansin gene family in Populus and forms a basis for future functional research studies.

1. Introduction

Expansins (EXPs) are non-hydrolytic cell wall relaxation proteins that are involved in developmental processes that alter the cell wall in plants. Expansins directly induce the expansion of the cell wall by disrupting non-covalent bonds between cellulose microfibrils and associated matrix polysaccharides [1]. Expansins were first discovered by Cosgrove in 1989 by studying the acid-induced extension of the hypocotyl cell wall of cucumbers [2]. Subsequent studies demonstrated that genes encoding extended proteins expanded rapidly during speciation and formed large gene families in plants that are widely distributed in various plant genomes in the form of gene families [3].
Canonical plant expansins are small proteins of 250–275 amino acids that are thought to have two conserved domains. The N-terminal domain (D1) is a six-stranded double-psi beta-barrel (DPBB) that is characterized by a His-Phe-Asp motif with conserved polar residues. The D1 domain has up to 30% sequence homology with the glycosyl hydrolase 45 family (GH45); however, expansins’ protein do not possess the hydrolase activity of GH45. The second domain (D2; Pollen allerg) contains conserved aromatic amino acids that are suitable for polysaccharide binding. The residues are aligned on the surface of an a-sandwich fold that resembles motifs from family 63 of carbohydrate-binding module domains (CBM63) [4].
Four subfamilies of expansins proteins have been identified based on phylogenetic sequence analysis, specifically, EXPA EXPB EXLA and EXLB, which are recognized in plants [5]. In addition to regulating cell size, different members of the extended protein family are involved in activities including morphogenesis [6,7,8,9,10,11], root hair growth [12,13,14], pollen fertilization [15], fruit softening [16] and role under adversity stress [17,18,19,20]. A genome-wide analysis of expansins previously identified 36 expansin genes in Arabidopsis [21]. Similarly, 56 expansin genes were identified in the rice genome [21] and 75 genes in soybeans [22].
The Populus genus is comprised of around 30 different species, including poplars, aspens and cottonwoods. Populus are widely distributed across the Northern hemisphere and many species and hybrids have been cultivated worldwide [23,24] as they are fast growing and can readily adapt to a wide range of ecological conditions. Populus are a key species for use as roadside trees and in protective forests and so have important ecological, economic and social value [24,25,26].
Plant phytohormones are a group of small, simple organic compounds that play a critical role in plant growth, development, and their response to adversity [27,28,29]. Previous studies have demonstrated that external spraying of salicylic acid enhances rust resistance synthesis in poplars [30]. It was also discovered that gibberellin (GA) treatment reduced the number of adventitious roots in wild-type poplars [31]. Abiotic stresses, which encompass heat, cold, waterlogging, drought, salt, metals, and nutritional deficiencies, are pressures that plants experience that are not caused by the natural environment. Abiotic stresses significantly impact the growth, development and productivity of trees [32,33]. P. cathayana plants exhibited increased antioxidant activity and reduced growth and photosystem II efficiency under salt stress [34]. Low temperatures can cause various types of physiological harm, including a reduction in photosynthetic rate, an increase in reactive oxygen species (ROS), decreased nutrient absorption, altered membrane transport and reduced nutrient absorption [35]. P. cathayana experienced reductions in growth and physiological functions under drought conditions [36].
Previous research studies have explored the structure and function of expansins in several plant species, such as Arabidopsis [37,38,39] and rice [40], yet little is known about the role of expansins in Populus. Recently, the availability of the poplar genome [23,41,42] and the genome database [43] have enabled the analysis of the expansin gene family.
The aim of this study was to investigate specific factors regarding the expansin gene family in poplar, its expression pattern in different tissues, and the gene expression response to phytohormone and abiotic stress. This research provides a theoretical basis for further understanding of the molecular processes and roles of expansin genes in the growth, development and stress resistance of poplar.

2. Results

2.1. Identification of Expansin Genes in Populus and Phylogenetic Relationships

Using the conserved amino acid sequences of the DPBB_1 (Pfam: PF03330) and Pollen_allerg_1 (Pfam: PF01357) domains as queries in HMMsearch, we identified 36 expansin family members in the Populus tremula L. All of the family members were divided into four subfamilies consisting of 27 EXPA members, three EXPB members, two EXLA members, and four EXLB members. These members were named PtEXPA1 to PtEXLB4 based on their chromosomal locations (Table 1). Further analysis showed that the number of amino acid (AA) residues in the proteins ranged from 185 to 596, with an average of 274 residues. The average molecular weight (MW) was 29.88 kD, which ranged from 20.64 to 67.81 kD. The isoelectric point (pI) for the proteins was between 5.26, and 10.22. 6 of the proteins were acidic and 30 were alkaline.
To study the phylogenetic relationships of expansins, a maximum likelihood phylogenic tree was constructed using MEGAX software based on multiple alignments of 146 expansins from Populus tremula L., Oryza sativa L., Arabidopsis thaliana and Carica papaya L. (Figure 1). All of the expansin genes were divided into four subfamilies, specifically, EXPA, EXPB, EXLA and EXLB. The largest subfamily was EXPA and the smallest subfamily was EXLB.

2.2. Gene Structures and Conserved Protein Motifs of PtEXs

Structural analysis can provide valuable information for the classification of gene evolution events and subfamilies. We analyzed the gene structures and conserved motifs of the expansin family members (Figure 2). Nearly all of the genes had UTR sections at both terminals and the structural patterns of the genes in the same subfamily were similar. Most members of the EXPA subfamily had three or four exons except for PtEXPA10 (eight exons), PtEXPA24 and PtEXPA25 (two exons). The EXPB subfamily had four exons and the EXLA subfamily had five exons. Based on the number of exons, the EXLB subfamily can be separated into two groups with four and five exons, respectively.
Members of the PtEXs gene family had 10 different conserved motifs that are summarized in Supplementary Figure S1 and Supplementary Table S1. The fundamental principle was that the motif composition of the peer group was characterized by the same or similar structure, for example, motif 1,4,9 appeared in nearly all members of the EXPA subfamily. In the other three subfamilies, the motifs were similar and motifs 7 and 8 were conserved in the other three subfamilies. Motif 2 was also found in all three subfamilies of EXPA, EXPB and EXLA, but not in EXLB. These results indicated that the conserved motifs may play critical roles in specific functions, or have similar functionality. Whilst the functions of some motifs were not yet clear, the presence of these conserved motifs reflected functional similarities among the PtEXs.

2.3. Chromosome Distributions of the PtEXs

The chromosomal distribution map of the expansin genes was generated based on the genome data, and 36 expansin genes were unevenly distributed on 14 chromosomes (Figure 3). Chromosomes 7, 11, 12, 15, and 18 did not have expansin genes. Chromosome 1 contained the largest number of expansin genes (up to six), while chromosomes 5 and 14 contain only one expansin gene. The number of expansin genes on the other chromosomes ranged from two to four.

2.4. Duplication Events of Expansins

Many gene families in plants occur as a result of tandem or segment duplications. To better understand the evolution of PtEXs genes, we investigated genome duplication events in this family. The synteny relationship between the PtEXs genes was discovered using MCScanX in TBtools software (v1.108 Chengjie Chen, Guangzhou, China). A total of 20 pairs of 23 genes were obtained and the data were visualized using Circos in TBtools software (v1.108 Chengjie Chen, Guangzhou, China) (Figure 4). The results identified 20 pairs of genes with segmental duplication, suggesting that segmental duplication is the main driver of the evolutionary expansion of the PtEXs gene family. The Ka and Ks of the complex base pairs were calculated using TBtools. The calculated Ks values ranged from 0.220531 to 2.402618, suggesting that replication occurred between 132.01 million years ago and 12.11 million years ago. The Ka/Ks replication gene pairs were less than 1, indicating that the PtEXs gene family was subjected to purification selection during the evolutionary process (Table 2).

2.5. Analysis of Cis-Acting Elements

We extracted the 1500 bp upstream of the 36 PtEXs genes from the initiation codon to study the cis-acting elements of expansin genes (Figure 5). A total of nine cis-acting elements were identified, amongst which Methyl jasmonate (MeJA), salicylic acid (SA), abscisic acid (ABA) and gibberellin (GA) are associated with responses to phytohormone stress, drought-inducible and low-temperature response (LTR), which are abiotic stresses, and light-response, related to plant development. Phytohormone stress and abiotic stress cis-elements were abundant. Further analysis was performed on a number of the five main cis-acting elements (Figure 6). Most of the PtEXs genes contained more ABRE acting elements, and PtEXPA3 and PtEXPA13 had the highest number of ABREs, containing six. These data suggest that these genes may play an important role in ABA stress.
Concerning abiotic stresses, we found that some genes with higher ABRE numbers also have cis-acting elements associated with LTR, for example, PtEXPA13 and PtEXPA19. This suggests that the expression of these genes is temperature dependent. We then selected several genes based on the promoter number map for the next step of the analysis.

2.6. Analysis of PtEXs Expression in Populus

To understand the potential functions of PtEXs gene family members, we downloaded the RNA-seq (TPM values) data through the Popgenie (PlantGenIE.org: Home) website (Supplementary Table S2). The RNA-seq data were used to analyze expression levels in different tissues, including buds, petiole mature, the phloem/cambium, roots, twigs, seeds, flowers, leaves and suckers. The data is presented in the form of a heat map shown in Figure 7. The results showed that different expansin members displayed varying levels of expression in various tissues with differences in transcript levels. PtEXPA3/8/13/14/15/16/19/20/21/24/27, PtEXPB3, PtEXLA1 and PtEXLB1 were expressed in almost all of the tissues. PtEXPA5/22/25 was not expressed in the tissues. PtEXPA2/4/7/11 showed similar expression patterns and were only expressed in the buds and seeds, indicating that they may have shared biological processes. PtEXLB3 was only highly expressed in seeds and largely absent in other tissues, suggesting that it might have a role in seed germination. In general, the expression of the majority of expansin genes varied significantly in buds, flowers, leaves, petiole mature, the phloem/cambium, flowers, roots and seeds.

2.7. The Expression of PtEXs in Response to Phytohormone and Abiotic Stresses

Given that some phytohormones and abiotic stress-responsive components were identified in the promoter regions of PtEXs (Figure 6), we selected nine expansin genes based on promoter analysis. These genes were abundant in relevant cis-acting elements. For instance, PtEXPA13/17/18/19 each had at least four ABREs, while PtEXPA27 and PtEXLA2 had an LTR. The expression of PtEXs genes at 0, 3, 6, 12, 24 and 48 h after ABA or low-temperature treatments was assessed. Under ABA stress treatment, there were three types of response patterns for PtEXs gene expression. The first expression pattern involved PtEXs genes being induced to express at a high level at a specific time point, followed by a gradual decrease in expression. This pattern includes PtEXLA2, PtEXPB3 and PtEXPA12, which were significantly induced to express at high levels at 24 h, 12 h and 3 h, respectively. The second expression pattern category consists of genes significantly repressed by ABA, such as the PtEXPA27 gene. The third category is the expression pattern without a significant pattern, possibly because these genes were insensitive to ABA treatment. This category includes PtEXPA3, PtEXPA13, PtEXPA17, PtEXPA18 and PtEXPA19 (Figure 8A). Similarly, there were three response patterns for PtEXs gene expression following low-temperature treatment. The first category consists of genes induced by low temperature, including PtEXPB3, PtEXPA3 and PtEXPA19, all of which were significantly induced to express at high levels after 48 h of treatment. The PtEXPB3 gene was induced to express more than 150-fold, suggesting that it may be an essential candidate for responding to low temperature. The second category includes genes suppressed by low temperature, such as PtEXPA12, PtEXPA17 and PtEXPA27, all of which were significantly repressed after low-temperature treatment. The third category consists of genes with no apparent expression pattern, including PtEXLA2, PtEXPA13 and PtEXPA18 (Figure 8B).

3. Discussion

Plant growth is caused by the proliferation and enlargement of cells and is limited by the cell wall, which restricts the protoplasm. The cell wall restricts the rapid increase of the protoplasm in plant cells, whilst expansin proteins can loosen the cell wall by breaking the hydrogen bonds between cellulose microfibrils and hemicellulose. Expansin proteins play an important role in the growth and development of plants. In this study, we analyzed the Populus tremula L. using bioinformatics tools to gain a better understanding of the expasin gene family.
We identified 36 expansin genes with two conserved domains, DPBB_1 and Pollen_allerg_1. The 36 tremula expansins were grouped into four subfamilies, PtEXPA, PtEXPB, PtEXLA and PtEXLB that were similar to other plants. We discovered an uneven distribution of each gene subfamily among species by investigating and comparing the sizes of the expansin subfamilies in Arabidopsis, Oryza sativa, Chinese jujube and other plants (Table 3). Our results showed that EXPA occupied a high proportion of woody plants, suggesting an important role in cell wall regulation. We also found a higher number of EXPBs in non-woody plants compared to woody plants, suggesting that cell wall regulation by EXPB is more effective in non-woody plants.
One of the main mechanisms driving the evolution of genomes and genetic systems is gene duplication [52]. Previous studies have demonstrated that tandem and segmental duplication are the primary forces behind the growth of gene families [53]. A total of 20 pairs of replication genes in the PtEXs gene family were identified that were all segmental. These data suggest that segmental replication is the main driver of the evolutionary expansion of the PtEXs gene family. Populus has undergone three genome-wide replication events, namely the ancient replication event (100–120 million years ago), the true rose branch replication event and the Populus family replication event (60–65 million years ago) [23]. Recent studies have shown that a whole genome duplication event also occurred in angiosperms 20 million years ago when the temperature and CO2 concentration were low. Most of the retained genes were associated with abiotic stresses, such as salt stress, low-temperature stress and drought stress [54]. Our analysis found that PtEXPA15/PtEXPA19, PtEXLB1/PtEXLB3 and PtEXPA3/PtEXPA17 had 3 pairs of gene duplication times that all coincide with the timing of this Genome-wide duplication event. We hypothesize that these genes play an important role in adversity stress.
Understanding gene expression patterns in tissues is crucial for the mining of functional genes. In previous studies, the expansin genes are expressed in one or more tissues and they are known to be essential for the growth and development of plants. In this study, we analyzed the expression patterns of 36 PtEXs genes in the buds, flowers, leaves, petiole mature, the phloem/cambium, flowers, roots, twigs, seeds and suckers of Populus tremula L. Most PtEXs showed tissue-specific expression. PtEXPA19 and PtEXPA15 were expressed at higher levels in roots and suckers, suggesting that this gene may play a role in the maturation and development of roots. According to earlier studies, AtEXP7 and AtEXP18 can control the initiation of root hairs, in Arabidopsis, RNA interference with the expression of AtEXP7 will result in oppositely shortened root hairs [55], and GmEXP1 ectopic expression of the gene promoted the growth of transgenic tobacco roots [56]. PtEXLA1 and PtEXLA2 were expressed at higher levels in the mature petiole compared to other genes, indicating that the EXLA gene subfamily plays a crucial role in petiole abscission. In Arabidopsis thaliana, the altered expression of expansins can modulate the development of leaves and pedicle abscission [7]. PtEXPA3/8/13/18 were expressed at higher levels in the phloem/cambium, suggesting that these genes may play an important role in wood formation.
Cis-acting elements play a crucial role in transcription and expression in plants [57,58,59]. We evaluated the putative cis-regulatory elements in the 1500-bp putative promoter regions of all of the PtEXs. We found that the promoters of the PtEXs genes mainly contained development-related and adaptation-related elements in response to plant phytohormone and abiotic stress. Previous results have shown that treatment with exogenous auxin after 24 to 48 h increased the expression of one expansin gene (accession no. AF085330) by 50 to 100 fold [60]. Ethylene can also induce the expression of Rp-EXP1 in the leaves of flood-tolerant species [61] and the expression of the LeEXP1 gene is regulated by ethylene in tomato fruits [62].
To further investigate the effect of phytohormone and abiotic stress on the expansin gene family, we selected nine genes, based on the number of cis-acting elements, to perform qRT-PCR experiments on ABA and low-temperature stress. We found that the expressions of some PtEXs were increased or repressed by ABA treatment and low temperatures. These experimental results were analyzed in combination with the number of cis-acting elements. We found that PtEXPA13/18/19 contained five ABREs and were expressed at lower levels following ABA treatment compared to PtEXPA12, which has only some ABRE. Similar data were found in response to low temperature stress. PtEXPB3 with the most significant expression at 48 h did not contain LTR cis-acting elements. Based on these data, we hypothesize that some genes gain new functions during the evolutionary process [63]. Another explanation is that the promoter region may contain cis-elements that enable the gene to respond to low-temperature stress [64].
In conclusion, our data suggest that the Expansin gene family in Populus plays an important role in controlling plant physiology and morphology, and has regulatory roles in response to stress. The genome-wide identification and characterization of the expansin gene family members in Populus serves as a key foundation for further investigating the function of these genes and may be useful in the breeding and genetic advancement of wood plants.

4. Materials and Methods

4.1. Identification of the Expansin Gene Family in Populus

The protein sequence information of poplar was downloaded from the PopGenie database (https://plantgenie.org/ accessed on 5 December 2022) to identify the expansin genes. Utilizing conserved domains of DPBB 1 (PF03330) and Pollen allerg 1 (PF01357) derived from the Pfam database (https://pfam.xfam.org/ accessed on 5 December 2022), a hidden Markov model (HMM) of the expansins was constructed. A search for Hidden Markov Model protein sequences with PF03330 and PF01357 in the Populus protein database was performed using HMMsearch with a threshold of e value < 10−5 [65]. The results of the initial screening were manually compared with the SMART (Simple Modular Architecture Research Tool) and the NCBI-CDD databases (NCBI conservative domain database) to confirm that the identified members contained both structural domains. For each protein sequence, we analyzed the molecular weight (MW) amino acids (AA) and isoelectric point (pI) on the ExPaSy (https://web.expasy.org/compute_pi/ accessed on 5 December 2022).

4.2. Phylogenetic Analysis

To analyze phylogenetic relationships, the expansin protein sequences for Oryza sativa L., Arabidopsis thaliana and Carica papaya L. were obtained from EXPANSIN CENTRAL (http://www.personal.psu.edu/fsl/ExpCentral/ accessed on 12 December 2022). Multiple sequence alignments of the identified Populus tremula L., Oryza sativa L., Arabidopsis thaliana and Carica papaya L. expansins were executed using the MEGA X software [66] and a phylogenetic tree was constructed using the Maximum likelihood method. The bootstrap replicates were set to 1000 and all other parameters were left at their default values.

4.3. Gene Structure, Chromosomal Locations and Cis-Regulatory Elements

Using the obtained genome annotation data, TBtools was used to map the exon gene structure and chromosome locations [67]. The conserved motifs were identified using the MEME tool (http://meme-suite.org/tools/meme accessed on 15 December 2022). Plant CARE (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/ accessed on 15 December 2022) was used to identify cis-regulatory elements in the 1.5 kb upstream sequences of each expansin gene. The cis-acting elements were visualized using TBtools.

4.4. Gene Duplication Analyses

The PtEXs gene family tandem and segmental duplications were investigated using the McscanX in TBtools software (v1.108 Chengjie Chen, Guangzhou, China) [68]. The Advanced Circos function of the TBtools software was used to visualize segmental duplication relationships [67]. TBtools’ simple Ka/Ks calculator was used to calculate Ka/Ks values for collinear pairs. The evolution time (T) was calculated according to the Ks value: T = Ks/2λ, λ = 9.1 × 10−9 [63].

4.5. Transcriptome Data Analysis

The transcriptional data of PtEXs from 9 tissues/organs in Populus were collected from the PopGenie database (PlantGenIE.org: Home) and a correlation heatmap was analyzed using TBtools.

4.6. Plant Materials and Treatments

45-day-old Populus davidiana × P. bolleana tissue-cultured seedlings were grown by the Key Laboratory of Forest Tree Genetic Breeding, Nanjing Forestry University. The laboratory conditions were as follows: temperature 20 ± 5 °C and substantial daily sunshine (16-h light from 07:00 to 23:00 h, 1000–2000 lx), and 70–80% humidity. Subsequently, we used 36 uniformly growing plants, dividing them into three groups of 12 plants each, with a total of three biological replicates. In each biological replicate, six plants were randomly selected for ABA treatment and low-temperature treatment, respectively. At each of the time points 0, 3, 6, 12, 24 and 48 h after the stress treatment, we collected the second to fourth node leaves of one plant, counting from the top to the bottom. For ABA stress treatment, the plants were fully stressed by mashing the medium. A 100 μM [64,69] ABA solution was poured onto the crushed tissue culture medium so that the plant roots were immersed in the solution. For the low-temperature treatment, the plants were placed in a Percival incubator (Percival, CU-22L) with the temperature set to 10 °C [70,71]. The leaves of the plants were harvested at 6 time points (0, 3, 6, 12, 24 and 48 h) [64,69,72] and immediately frozen at −80 °C.

4.7. RNA Extraction and qRT-PCR Analysis

RNA was extracted from the leaves of Populus davidiana × P. bolleana using an RNA extraction kit (TIANGEN, Nanjing, China). RNA was detected by 1% agarose gel electrophoresis. RNA was reverse transcribed to synthesize the first strand using a 1-step kit (TIANGEN, Nanjing, China). The cDNA was diluted at a concentration of 1:10. Primers were designed based on the sequences of Populus davidiana × P. bolleana in the CDS database using Primer 5 software, and the amplification efficiency of each primer pair was analyzed (Supplementary Figure S2 and Supplementary Table S3). The ubiquitin gene (UBQ, gene ID Potri.001G418500) was used as a reference control gene [73]. The qRT-PCR reaction was carried out under the following conditions; 1 cycle at 98 °C for 3 min, then 40 cycles at 95 °C for 15 s, 60 °C for 30 s, and 72 °C for 30 s. The 2−ΔΔCt method was used to evaluate genes expression levels [74]. Each sample was analyzed as 3 biological replicates with 3 technical replicates.

4.8. Statistical Analysis

For statistical analysis we use GraphPad Prism v8.0.2 software. One-way ANOVA was used to compare the differences between means. Statistically significant difference was considered at * p < 0.05 ** p < 0.01. The gene expression during 0-h stress treatment was used as a control for significant analysis.

5. Conclusions

In this study, we used bioinformatics tools to systematically analyze the expansin gene family in Populus. We identified 36 expansin genes that were divided into four subfamilies, EXPA, EXPB, EXLA and EXLB. The gene structure and conserved domain maps were compared according to different families. PtEXs in the same family had similar gene structures and conserved domains. Collinearity, Ka and Ks analysis of the evolution of the PtEXs gene family showed that segmental duplication was the main driving force for the expansion of the genes, and duplicate genes were subject to strong purifying selection in the evolution process. Transcriptome data demonstrated the tissue specific expression of expansins. Combined cis-acting elements and qRT-PCR analysis indicated that some genes may play important roles in stress. Our results will provide a theoretical basis for genetic improvements of the expansin gene family in poplar.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms24097759/s1.

Author Contributions

Conceptualization, T.Y. and H.W.; software, Z.Y. and F.Z.; validation, Z.Y., F.Z. and Y.C.; data curation, Z.Y. and F.Z.; writing—original draft preparation, Z.Y.; writing—review and editing, H.W.; visualization, Z.Y.; supervision, T.Y. and H.W.; project administration, T.Y.; funding acquisition, T.Y., Y.C. and H.W. All authors have read and agreed to the published version of the manuscript.

Funding

This study was financially supported by the National Key Research and Development Program of China (grant number 2021YFD2201100); the Natural Science Foundation of China (grant numbers 32071795); the Youth Science and Technology Innovation Fund of Nanjing Forestry University (grant number CX2018008).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data generated in the current work is provided in the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. McQueen-Mason, S.; Cosgrove, D.J. Disruption of hydrogen bonding between plant cell wall polymers by proteins that induce wall extension. Proc. Natl. Acad. Sci. USA 1994, 91, 6574–6578. [Google Scholar] [CrossRef] [Green Version]
  2. Cosgrove, D.J. Characterization of long-term extension of isolated cell walls from growing cucumber hypocotyls. Planta 1989, 177, 121–130. [Google Scholar] [CrossRef]
  3. Lee, Y.; Choi, D.; Kende, H. Expansins: Ever-expanding numbers and functions. Curr. Opin. Plant Biol. 2001, 4, 527–532. [Google Scholar] [CrossRef] [PubMed]
  4. Georgelis, N.; Yennawar, N.H.; Cosgrove, D.J. Structural basis for entropy-driven cellulose binding by a type-A cellulose-binding module (CBM) and bacterial expansin. Proc. Natl. Acad. Sci. USA 2012, 109, 14830–14835. [Google Scholar] [CrossRef] [Green Version]
  5. Kende, H.; Bradford, K.; Brummell, D.; Cho, H.-T.; Cosgrove, D.; Fleming, A.; Gehring, C.; Lee, Y.; Mcqueen-Mason, S.; Rose, J.; et al. Nomenclature for members of the expansin superfamily of genes and proteins. Plant Mol. Biol. 2004, 55, 311–314. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Keller, E.; Cosgrove, D.J. Expansins in growing tomato leaves. Plant J. 1995, 8, 795–802. [Google Scholar] [CrossRef] [PubMed]
  7. Cho, H.-T.; Cosgrove, D.J. Altered expression of expansin modulates leaf growth and pedicel abscission in Arabidopsis thaliana. Proc. Natl. Acad. Sci. USA 2000, 97, 9783–9788. [Google Scholar] [CrossRef] [Green Version]
  8. Pien, S.; Wyrzykowska, J.; McQueen-Mason, S.; Smart, C.; Fleming, A. Local expression of expansin induces the entire process of leaf development and modifies leaf shape. Proc. Natl. Acad. Sci. USA 2001, 98, 11812–11817. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Cho, H.-T.; Kende, H. Expansins and Internodal Growth of Deepwater Rice. Plant Physiol. 1997, 113, 1145–1151. [Google Scholar] [CrossRef] [Green Version]
  10. Choi, D.; Lee, Y.; Cho, H.-T.; Kende, H. Regulation of expansin gene expression affects growth and development in transgenic rice plants. Plant Cell 2003, 15, 1386–1398. [Google Scholar] [CrossRef] [Green Version]
  11. Choi, D.; Cho, H.-T.; Lee, Y. Expansins: Expanding importance in plant growth and development. Physiol. Plant. 2006, 126, 511–518. [Google Scholar] [CrossRef]
  12. Cho, H.-T.; Cosgrove, D.J. Regulation of root hair initiation and expansin gene expression in arabidopsis. Plant Cell 2002, 14, 3237–3253. [Google Scholar] [CrossRef] [Green Version]
  13. Wu, Y.; Thorne, E.T.; Sharp, R.E.; Cosgrove, D.J. Modification of expansin transcript levels in the maize primary root at low water potentials. Plant Physiol. 2001, 126, 1471–1479. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Guo, W.; Zhao, J.; Li, X.; Qin, L.; Yan, X.; Liao, H. A soybean β-expansin gene GmEXPB2 intrinsically involved in root system architecture responses to abiotic stresses. Plant J. 2011, 66, 541–552. [Google Scholar] [CrossRef]
  15. Pezzotti, M.; Feron, R.; Mariani, C. Pollination modulates expression of the PPAL gene, a pistil-specific β-expansin. Plant Mol. Biol. 2002, 49, 187–197. [Google Scholar] [CrossRef]
  16. Rose, J.K.C.; Lee, H.H.; Bennett, A.B. Expression of a divergent expansin gene is fruit-specific and ripening-regulated. Proc. Natl. Acad. Sci. USA 1997, 94, 5955–5960. [Google Scholar] [CrossRef] [Green Version]
  17. Geilfus, C.-M.; Zörb, C.; Mühling, K.H. Salt stress differentially affects growth-mediating β-expansins in resistant and sensitive maize (Zea mays L.). Plant Physiol. Biochem. 2010, 48, 993–998. [Google Scholar] [CrossRef]
  18. Han, Y.Y.; Li, A.X.; Li, F.; Zhao, M.R.; Wang, W. Characterization of a wheat (Triticum aestivum L.) expansin gene, TaEXPB23, involved in the abiotic stress response and phytohormone regulation. Plant Physiol. Biochem. 2012, 54, 49–58. [Google Scholar] [CrossRef]
  19. Kuluev, B.; Avalbaev, A.; Mikhaylova, E.; Nikonorov, Y.; Berezhneva, Z.; Chemeris, A. Expression profiles and hormonal regulation of tobacco expansin genes and their involvement in abiotic stress response. J. Plant Physiol. 2016, 206, 1–12. [Google Scholar] [CrossRef] [PubMed]
  20. Dotto, M.; Pombo, M.A.; Martínez, G.A.; Civello, P.M. Heat treatments and expansin gene expression in strawberry fruit. Sci. Hortic. 2011, 130, 775–780. [Google Scholar] [CrossRef]
  21. Sampedro, J.; Cosgrove, D.J. The expansin superfamily. Genome Biol. 2005, 6, 242. [Google Scholar] [CrossRef] [Green Version]
  22. Zhu, Y.; Wu, N.; Song, W.; Yin, G.; Qin, Y.; Yan, Y.; Hu, Y. Soybean (Glycine max) expansin gene superfamily origins: Segmental and tandem duplication events followed by divergent selection among subfamilies. BMC Plant Biol. 2014, 14, 93. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Lin, Y.-C.; Wang, J.; Delhomme, N.; Schiffthaler, B.; Sundström, G.; Zuccolo, A.; Nystedt, B.; Hvidsten, T.R.; De La Torre, A.; Cossu, R.M.; et al. Functional and evolutionary genomic inferences in Populus through genome and population sequencing of American and European aspen. Proc. Natl. Acad. Sci. USA 2018, 115, E10970–E10978. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Zhou, F.; Chen, Y.; Wu, H.; Yin, T. Genome-Wide Comparative Analysis of R2R3 MYB Gene Family in Populus and Salix and Identification of Male Flower Bud Development-Related Genes. Front. Plant Sci. 2021, 12, 721558. [Google Scholar] [CrossRef] [PubMed]
  25. Biselli, C.; Vietto, L.; Rosso, L.; Cattivelli, L.; Nervo, G.; Fricano, A. Advanced Breeding for Biotic Stress Resistance in Poplar. Plants 2022, 11, 2032. [Google Scholar] [CrossRef]
  26. Zhang, Y.; Zhou, F.; Wang, H.; Chen, Y.; Yin, T.; Wu, H. Genome-Wide Comparative Analysis of the Fasciclin-like Arabinogalactan Proteins (FLAs) in Salicacea and Identification of Secondary Tissue Development-Related Genes. Int. J. Mol. Sci. 2023, 24, 1481. [Google Scholar] [CrossRef]
  27. Santner, A.; A Calderon-Villalobos, L.I.; Estelle, M. Plant hormones are versatile chemical regulators of plant growth. Nat. Chem. Biol. 2009, 5, 301–307. [Google Scholar] [CrossRef]
  28. Popko, J.; Hänsch, R.; Mendel, R.-R.; Polle, A.; Teichmann, T. The role of abscisic acid and auxin in the response of poplar to abiotic stress. Plant Biol. 2010, 12, 242–258. [Google Scholar] [CrossRef] [PubMed]
  29. Hu, J.; Su, H.; Cao, H.; Wei, H.; Fu, X.; Jiang, X.; Song, Q.; He, X.; Xu, C.; Luo, K. AUXIN RESPONSE FACTOR7 integrates gibberellin and auxin signaling via interactions between DELLA and AUX/IAA proteins to regulate cambial activity in poplar. Plant Cell 2022, 34, 2688–2707. [Google Scholar] [CrossRef] [PubMed]
  30. Ullah, C.; Tsai, C.; Unsicker, S.B.; Xue, L.; Reichelt, M.; Gershenzon, J.; Hammerbacher, A. Salicylic acid activates poplar defense against the biotrophic rust fungus Melampsora larici-populina via increased biosynthesis of catechin and proanthocyanidins. New Phytol. 2019, 221, 960–975. [Google Scholar] [CrossRef] [Green Version]
  31. Mauriat, M.; Petterle, A.; Bellini, C.; Moritz, T. Gibberellins inhibit adventitious rooting in hybrid aspen and Arabidopsis by affecting auxin transport. Plant J. 2014, 78, 372–384. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Yang, Y.; Shao, Y.; Chaffin, T.A.; Lee, J.H.; Poindexter, M.R.; Ahkami, A.H.; Blumwald, E.; Stewart, C.N.J. Performance of abiotic stress-inducible synthetic promoters in genetically engineered hybrid poplar (Populus tremula × Populus alba). Front. Plant Sci. 2022, 13, 1011939. [Google Scholar] [CrossRef]
  33. Harfouche, A.; Meilan, R.; Altman, A. Molecular and physiological responses to abiotic stress in forest trees and their relevance to tree improvement. Tree Physiol. 2014, 34, 1181–1198. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Wu, N.; Li, Z.; Wu, F.; Tang, M. Comparative photochemistry activity and antioxidant responses in male and female Populus cathayana cuttings inoculated with arbuscular mycorrhizal fungi under salt. Sci. Rep. 2016, 6, 37663. [Google Scholar] [CrossRef] [PubMed]
  35. Zhou, B.; Kang, Y.; Leng, J.; Xu, Q. Genome-Wide Analysis of the miRNA–mRNAs Network Involved in Cold Tolerance in Populus simonii × P. nigra. Genes 2019, 10, 430. [Google Scholar] [CrossRef] [Green Version]
  36. Li, Y.; Duan, B.; Chen, J.; Korpelainen, H.; Niinemets, Ü.; Li, C. Males exhibit competitive advantages over females of Populus deltoides under salinity stress. Tree Physiol. 2016, 36, 1573–1584. [Google Scholar] [CrossRef] [Green Version]
  37. Zhang, X.-Q.; Wei, P.-C.; Xiong, Y.-M.; Yang, Y.; Chen, J.; Wang, X.-C. Overexpression of the Arabidopsis α-expansin gene AtEXPA1 accelerates stomatal opening by decreasing the volumetric elastic modulus. Plant Cell Rep. 2011, 30, 27–36. [Google Scholar] [CrossRef]
  38. Lee, H.W.; Kim, J. EXPANSINA17 up-regulated by LBD18/ASL20 promotes lateral root formation during the auxin response. Plant Cell Physiol. 2013, 54, 1600–1611. [Google Scholar] [CrossRef] [Green Version]
  39. Goh, H.-H.; Sloan, J.; Dorca-Fornell, M.D.C.; Fleming, A. Inducible repression of multiple expansin genes leads to growth suppression during leaf development. Plant Physiol. 2012, 159, 1759–1770. [Google Scholar] [CrossRef] [Green Version]
  40. Yu, Z.M.; Kang, B.; He, X.W.; Lv, S.L.; Bai, Y.H.; Ding, W.N.; Chen, M.; Cho, H.-T.; Wu, P. Root hair-specific expansins modulate root hair elongation in rice. Plant J. 2011, 66, 725–734. [Google Scholar] [CrossRef]
  41. Bai, S.; Wu, H.; Zhang, J.; Pan, Z.; Zhao, W.; Li, Z.; Tong, C. Genome Assembly of Salicaceae Populus deltoides (Eastern Cottonwood) I-69 Based on Nanopore Sequencing and Hi-C Technologies. J. Hered. 2021, 112, 303–310. [Google Scholar] [CrossRef] [PubMed]
  42. Wu, H.; Yao, D.; Chen, Y.; Yang, W.; Zhao, W.; Gao, H.; Tong, C. De Novo Genome Assembly of Populus simonii Further Supports That Populus simonii and Populus trichocarpa Belong to Different Sections. G3 Genes Genomes Genet. 2020, 10, 455–466. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Sjödin, A.; Street, N.R.; Sandberg, G.; Gustafsson, P.; Jansson, S. The Populus Genome Integrative Explorer (PopGenIE): A new resource for exploring the Populus genome. New Phytol. 2009, 182, 1013–1025. [Google Scholar] [CrossRef]
  44. Sampedro, J.; Lee, Y.; Carey, R.E.; DePamphilis, C.; Cosgrove, D.J. Use of genomic history to improve phylogeny and understanding of births and deaths in a gene family. Plant J. 2005, 44, 409–419. [Google Scholar] [CrossRef]
  45. Hou, L.; Zhang, Z.; Dou, S.; Zhang, Y.; Pang, X.; Li, Y. Genome-wide identification, characterization, and expression analysis of the expansin gene family in Chinese jujube (Ziziphus jujuba Mill.). Planta 2019, 249, 815–829. [Google Scholar] [CrossRef] [PubMed]
  46. Zhang, S.; Xu, R.; Gao, Z.; Chen, C.; Jiang, Z.; Shu, H. A genome-wide analysis of the expansin genes in Malus × Domestica. Mol. Genet. Genom. 2014, 289, 225–236. [Google Scholar] [CrossRef]
  47. Zhang, W.; Yan, H.; Chen, W.; Liu, J.; Jiang, C.; Jiang, H.; Zhu, S.; Cheng, B. Genome-wide identification and characterization of maize expansin genes expressed in endosperm. Mol. Genet. Genom. 2014, 289, 1061–1074. [Google Scholar] [CrossRef]
  48. Feng, X.; Li, C.; He, F.; Xu, Y.; Li, L.; Wang, X.; Chen, Q.; Li, F. Genome-Wide Identification of Expansin Genes in Wild Soybean (Glycine soja) and Functional Characterization of Expansin B1 (GsEXPB1) in Soybean Hair Root. Int. J. Mol. Sci. 2022, 23, 5407. [Google Scholar] [CrossRef]
  49. Santo, S.D.; Vannozzi, A.; Tornielli, G.B.; Fasoli, M.; Venturini, L.; Pezzotti, M.; Zenoni, S. Genome-wide analysis of the expansin gene superfamily reveals grapevine-specific structural and functional characteristics. PLoS ONE 2013, 8, e62206. [Google Scholar] [CrossRef]
  50. Lu, Y.; Liu, L.; Wang, X.; Han, Z.; Ouyang, B.; Zhang, J.; Li, H. Genome-wide identification and expression analysis of the expansin gene family in tomato. Mol. Genet. Genom. 2016, 291, 597–608. [Google Scholar] [CrossRef]
  51. Ding, A.; Marowa, P.; Kong, Y. Genome-wide identification of the expansin gene family in tobacco (Nicotiana tabacum). Mol. Genet. Genom. 2016, 291, 1891–1907. [Google Scholar] [CrossRef]
  52. Moore, R.C.; Purugganan, M.D. The early stages of duplicate gene evolution. Proc. Natl. Acad. Sci. USA 2003, 100, 15682–15687. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Cannon, S.B.; Mitra, A.; Baumgarten, A.; Young, N.D.; May, G. The roles of segmental and tandem gene duplication in the evolution of large gene families in Arabidopsis thaliana. BMC Plant Biol. 2004, 4, 10. [Google Scholar] [CrossRef] [Green Version]
  54. Wu, S.; Han, B.; Jiao, Y. Genetic Contribution of Paleopolyploidy to Adaptive Evolution in Angiosperms. Mol. Plant 2020, 13, 59–71. [Google Scholar] [CrossRef]
  55. Ma, Z.; Bielenberg, D.G.; Brown, K.M.; Lynch, J.P. Regulation of root hair density by phosphorus availability in Arabidopsis thaliana. Plant Cell Environ. 2001, 24, 459–467. [Google Scholar] [CrossRef] [Green Version]
  56. Lee, D.-K.; Ahn, J.H.; Song, S.-K.; Choi, Y.D.; Lee, J.S. Expression of an expansin gene is correlated with root elongation in soybean. Plant Physiol. 2003, 131, 985–997. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Lü, S.; Gu, H.; Yuan, X.; Wang, X.; Wu, A.-M.; Qu, L.; Liu, J.-Y. The GUS reporter-aided analysis of the promoter activities of a rice metallothionein gene reveals different regulatory regions responsible for tissue-specific and inducible expression in transgenic Arabidopsis. Transgenic Res. 2007, 16, 177–191. [Google Scholar] [CrossRef] [PubMed]
  58. Walther, D.; Brunnemann, R.; Selbig, J. The regulatory code for transcriptional response diversity and its relation to genome structural properties in a. thaliana. PLoS Genet. 2007, 3, e11. [Google Scholar] [CrossRef]
  59. Le, D.T.; Nishiyama, R.; Watanabe, Y.; Vankova, R.; Tanaka, M.; Seki, M.; Ham, L.H.; Yamaguchi-Shinozaki, K.; Shinozaki, K.; Tran, L.-S.P. Identification and expression analysis of cytokinin metabolic genes in soybean under normal and drought conditions in relation to cytokinin levels. PLoS ONE 2012, 7, e42411. [Google Scholar] [CrossRef] [Green Version]
  60. Hutchison, K.W.; Singer, P.B.; McInnis, S.; Diaz-Sala, C.; Greenwood, M.S. Expansins are conserved in conifers and expressed in hypocotyls in response to exogenous auxin1. Plant Physiol. 1999, 120, 827–832. [Google Scholar] [CrossRef] [Green Version]
  61. Vriezen, W.H.; De Graaf, B.; Mariani, C.; Voesenek, L. Submergence induces expansin gene expression in flooding-tolerant Rumex palustris and not in flooding-intolerant R. acetosa. Planta 2000, 210, 956–963. [Google Scholar] [CrossRef]
  62. Xue, Z.; Kou, X.; Luo, Y.; Zhu, B.; Xu, W. Effect of ethylene on polygalacturonase, lipoxygenase and expansin in ripening of tomato fruits. Trans. Tianjin Univ. 2009, 15, 173–177. [Google Scholar] [CrossRef]
  63. Lynch, M.; Conery, J.S. The evolutionary fate and consequences of duplicate genes. Science 2000, 290, 1151–1155. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Liu, R.; Wu, M.; Liu, H.; Gao, Y.; Chen, J.; Yan, H.; Xiang, Y. Genome-wide identification and expression analysis of the NF-Y transcription factor family in Populus. Physiol. Plant. 2021, 171, 309–327. [Google Scholar] [CrossRef] [PubMed]
  65. Finn, R.D.; Bateman, A.; Clements, J.; Coggill, P.; Eberhardt, R.Y.; Eddy, S.R.; Heger, A.; Hetherington, K.; Holm, L.; Mistry, J.; et al. Pfam: The protein families database. Nucleic Acids Res. 2014, 42, D222–D230. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef]
  67. Chen, C.J.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.H.; Xia, R. TBtools: An Integrative Toolkit Developed for Interactive Analyses of Big Biological Data. Mol. Plant 2020, 13, 1194–1202. [Google Scholar] [CrossRef]
  68. Wang, Y.; Tang, H.; DeBarry, J.D.; Tan, X.; Li, J.; Wang, X.; Lee, T.-H.; Jin, H.; Marler, B.; Guo, H.; et al. MCScanX: A toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012, 40, e49. [Google Scholar] [CrossRef] [Green Version]
  69. Zhang, C.; Wang, W.; Wang, D.; Hu, S.; Zhang, Q.; Wang, Z.; Cui, L. Genome-Wide Identification and Characterization of the WRKY Gene Family in Scutellaria baicalensis Georgi under Diverse Abiotic Stress. Int. J. Mol. Sci. 2022, 23, 4225. [Google Scholar] [CrossRef]
  70. Kong, R.S.; Way, D.A.; Henry, H.A.L.; Smith, N.G. Stomatal conductance, not biochemistry, drives low temperature acclimation of photosynthesis in Populus balsamifera, regardless of nitrogen availability. Plant Biol. 2022, 24, 766–779. [Google Scholar] [CrossRef]
  71. Eom, S.H.; Ahn, M.-A.; Kim, E.; Lee, H.J.; Lee, J.H.; Wi, S.H.; Kim, S.K.; Bin Lim, H.; Hyun, T.K. Plant Response to Cold Stress: Cold Stress Changes Antioxidant Metabolism in Heading Type Kimchi Cabbage (Brassica rapa L. ssp. Pekinensis). Antioxidants 2022, 11, 700. [Google Scholar] [CrossRef]
  72. Lv, T.; Liu, Q.; Xiao, H.; Fan, T.; Zhou, Y.; Wang, J.; Tian, C.-E. Genome-wide identification and analysis of the IQM gene family in soybean. Front. Plant Sci. 2022, 13, 1093589. [Google Scholar] [CrossRef]
  73. Gutierrez, L.; Mauriat, M.; Guénin, S.; Pelloux, J.; Lefebvre, J.-F.; Louvet, R.; Rusterucci, C.; Moritz, T.; Guerineau, F.; Bellini, C.; et al. The lack of a systematic validation of reference genes: A serious pitfall undervalued in reverse transcription-polymerase chain reaction (RT-PCR) analysis in plants. Plant Biotechnol. J. 2008, 6, 609–618. [Google Scholar] [CrossRef] [Green Version]
  74. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
Ijms 24 07759 g001 550
Figure 1. A phylogenetic tree of the expansins from Populus tremula L., Oryza sativa L., Arabidopsis thaliana and Carica papaya L. MEGAX was used to construct a Maximum likelihood phylogenetic tree with 1000 bootstrap replications. Circles, squares, stars and triangle represent the expansins of Oryza sativa L., Arabidopsis thaliana, Carica papaya L. and Populus tremula L. Purple, pink, orange and green colors represent the EXPA, EXLA, EXLB and EXPB subfamilies.
Figure 1. A phylogenetic tree of the expansins from Populus tremula L., Oryza sativa L., Arabidopsis thaliana and Carica papaya L. MEGAX was used to construct a Maximum likelihood phylogenetic tree with 1000 bootstrap replications. Circles, squares, stars and triangle represent the expansins of Oryza sativa L., Arabidopsis thaliana, Carica papaya L. and Populus tremula L. Purple, pink, orange and green colors represent the EXPA, EXLA, EXLB and EXPB subfamilies.
Ijms 24 07759 g001
Ijms 24 07759 g002 550
Figure 2. The phylogenetic relationships, exon-intron and motif structures of poplar expansin genes. (A) The phylogenetic tree was constructed based on the PtEXs sequences. According to phylogenetic relationships, 36 PtEXs were clustered into PtEXPA, PtEXPB, PtEXLA and PtEXLB groups. (B) The 10 motifs of expansin proteins are distinguished by different colors. (C) Gene structures of the expansin genes in poplar.
Figure 2. The phylogenetic relationships, exon-intron and motif structures of poplar expansin genes. (A) The phylogenetic tree was constructed based on the PtEXs sequences. According to phylogenetic relationships, 36 PtEXs were clustered into PtEXPA, PtEXPB, PtEXLA and PtEXLB groups. (B) The 10 motifs of expansin proteins are distinguished by different colors. (C) Gene structures of the expansin genes in poplar.
Ijms 24 07759 g002
Ijms 24 07759 g003 550
Figure 3. Chromosomal locations of the identified PtEXs in Populus. The color gradient from blue to red on the chromosomes indicates the gene density (from low to high).
Figure 3. Chromosomal locations of the identified PtEXs in Populus. The color gradient from blue to red on the chromosomes indicates the gene density (from low to high).
Ijms 24 07759 g003
Ijms 24 07759 g004 550
Figure 4. Collinearity mapping of expansin genes in the Populus genome. The red, blue, green and orange lines show homologous gene pairs representing the EXPA, EXPB, EXLA and EXLB subfamilies, respectively. From the outside to the inside, the first circle represents chromosome coordinates and the second circle represents gene density distribution.
Figure 4. Collinearity mapping of expansin genes in the Populus genome. The red, blue, green and orange lines show homologous gene pairs representing the EXPA, EXPB, EXLA and EXLB subfamilies, respectively. From the outside to the inside, the first circle represents chromosome coordinates and the second circle represents gene density distribution.
Ijms 24 07759 g004
Ijms 24 07759 g005 550
Figure 5. A 1500 bp upstream sequence was used to predict the cis-elements. Each of the 9 predicted cis-elements is represented by a different colored box.
Figure 5. A 1500 bp upstream sequence was used to predict the cis-elements. Each of the 9 predicted cis-elements is represented by a different colored box.
Ijms 24 07759 g005
Ijms 24 07759 g006 550
Figure 6. Number of major cis-acting elements of 36 Populus expansin genes. The blank box indicates that the quantity was 0. Different colors indicate different number of cis-acting elements.
Figure 6. Number of major cis-acting elements of 36 Populus expansin genes. The blank box indicates that the quantity was 0. Different colors indicate different number of cis-acting elements.
Ijms 24 07759 g006
Ijms 24 07759 g007 550
Figure 7. The expression profiles of 36 Populus expansin genes in different tissues. The legend represents the logarithmic normalized TPM.
Figure 7. The expression profiles of 36 Populus expansin genes in different tissues. The legend represents the logarithmic normalized TPM.
Ijms 24 07759 g007
Ijms 24 07759 g008 550
Figure 8. The expression profiles of the 9 PtEXs in 45-day-old Populus davidiana × P. bolleana tissue-cultured seedlings. (A) ABA and (B) low-temperature treatments. The expression of the ubiquitin (UBQ) housekeeping gene was used as a control. The x-axis represents the time and the y-axis represents the level of expression. The data are presented as the mean ± SD (n = 3). An asterisk indicates that the expression level after stress was significantly different to the level before the stress (* p < 0.05, ** p < 0.01).
Figure 8. The expression profiles of the 9 PtEXs in 45-day-old Populus davidiana × P. bolleana tissue-cultured seedlings. (A) ABA and (B) low-temperature treatments. The expression of the ubiquitin (UBQ) housekeeping gene was used as a control. The x-axis represents the time and the y-axis represents the level of expression. The data are presented as the mean ± SD (n = 3). An asterisk indicates that the expression level after stress was significantly different to the level before the stress (* p < 0.05, ** p < 0.01).
Ijms 24 07759 g008
Table
Table 1. Summary of information relating to the 36 PtEXs proteins.
Table 1. Summary of information relating to the 36 PtEXs proteins.
Gene NameGene IdpIMw/kDAAs
EXPAPtEXPA1Potra2n1c11.18.9930.19278
EXPAPtEXPA2Potra2n1c960.19.4929.15269
EXPAPtEXPA3Potra2n1c2087.19.4829.02267
EXPAPtEXPA4Potra2n1c3505.19.6925.33231
EXPAPtEXPA5Potra2n2c4734.110.2220.64185
EXPAPtEXPA6Potra2n2c4737.18.9827.76256
EXPAPtEXPA7Potra2n2c6293.19.127.83256
EXPAPtEXPA8Potra2n4c9142.18.8830.02275
EXPAPtEXPA9Potra2n4c9551.18.6625.81241
EXPAPtEXPA10Potra2n4c10239.18.7767.81596
EXPAPtEXPA11Potra2n5c10644.19.2727.79256
EXPAPtEXPA12Potra2n6c14375.19.1427.62255
EXPAPtEXPA13Potra2n6c14588.19.6234.96323
EXPAPtEXPA14Potra2n8c17125.19.52.79258
EXPAPtEXPA15Potra2n8c17409.19.5326.72251
EXPAPtEXPA16Potra2n9c18642.18.0927.80250
EXPAPtEXPA17Potra2n9c19851.19.433.06302
EXPAPtEXPA18Potra2n10c20623.19.3730.40279
EXPAPtEXPA19Potra2n10c20953.19.930.34279
EXPAPtEXPA20Potra2n13c24992.17.5226.70250
EXPAPtEXPA21Potra2n13c25752.19.0125.59242
EXPAPtEXPA22Potra2n16c30169.19.3435.69323
EXPAPtEXPA23Potra2n16c30491.18.3727.38254
EXPAPtEXPA24Potra2n17c30722.18.8929.94274
EXPAPtEXPA25Potra2n17c31117.18.6928.32257
EXPAPtEXPA26Potra2n17c31179.19.4428.55264
EXPAPtEXPA27Potra2n19c33898.16.0126.54250
EXPBPtEXPB1Potra2n13c25121.17.6128.71262
EXPBPtEXPB2Potra2n14c26821.15.6828.91274
EXPBPtEXPB3Potra2n19c33436.18.5128.60262
EXLAPtEXLA1Potra2n4c10039.18.2527.74252
EXLAPtEXLA2Potra2n9c18880.18.9629.88274
EXLBPtEXLB1Potra2n1c1267.15.2630.20275
EXLBPtEXLB2Potra2n1c1318.16.5731.54290
EXLBPtEXLB3Potra2n3c7743.15.7133.60302
EXLBPtEXLB4Potra2n3c7783.16.7127.61250
Table
Table 2. Ks and Ka analysis of duplicated gene pairs.
Table 2. Ks and Ka analysis of duplicated gene pairs.
Gene 1Gene 2KaKsKa/KsDuplication TypeT (MYA) 1
PtEXLB2PtEXLB40.0783120.3403650.230082Segmental Duplication18.7
PtEXPA19PtEXPA150.0287160.3695110.077713Segmental Duplication20.3
PtEXLB1PtEXLB30.1085420.3654740.296989Segmental Duplication20.08
PtEXPA7PtEXPA110.0518960.3259830.1592Segmental Duplication17.91
PtEXLA1PtEXLA20.0750970.4314830.174043Segmental Duplication23.7
PtEXPA3PtEXPA140.0999921.261390.079272Segmental Duplication69.3
PtEXPA14PtEXPA170.0976091.5495210.062993Segmental Duplication85.13
PtEXPA3PtEXPA170.0376180.3670910.102475Segmental Duplication20.16
PtEXPA23PtEXPA270.0943012.3999050.039294Segmental Duplication131.86
PtEXPA23PtEXPA120.0635630.489160.129944Segmental Duplication26.87
PtEXPB1PtEXPB30.0377490.2205310.171174Segmental Duplication12.11
PtEXPA3PtEXPA180.1310031.8171190.072094Segmental Duplication99.84
PtEXPA18PtEXPA140.0812690.4064430.199952Segmental Duplication22.33
PtEXPA18PtEXPA170.1798241.5640930.11497Segmental Duplication85.93
PtEXPA19PtEXPA260.1458771.03020.141601Segmental Duplication56.6
PtEXPA26PtEXPA90.0554520.4482070.123719Segmental Duplication24.62
PtEXPA26PtEXPA150.1081331.4712180.073499Segmental Duplication80.83
PtEXPA3PtEXPA130.1197061.6857950.071008Segmental Duplication92.62
PtEXPA13PtEXPA170.1242161.5911130.078069Segmental Duplication87.42
PtEXPA22PtEXPA140.4042782.4026180.168266Segmental Duplication132.01
1 million years ago (Mya).
Table
Table 3. Summary of each expansin subfamily in 10 plant species.
Table 3. Summary of each expansin subfamily in 10 plant species.
SpeciesEXPAEXPBEXLAEXLBTotalReference
Populus tremula L.27 (75%)3 (8.3%)2 (5.5%)4 (11.1%)36In this study
Arabidopsis25 (71.4%)6 (17.1%)3 (8.6%)1 (2.9%)35[44]
Oryza sativa34 (58.6%)19(32.8%)4 (6.9%)1 (1.7%)58[44]
Chinese jujube19 (63.3%)3 (10.0%)1 (3.3%)7 (23.3%)30[45]
Apple34 (82.9%)1 (2.4%)2 (4.9%)4 (9.8%)41[46]
maize36 (40.9%)48 (54.5%)4 (4.5%)0 (0%)88[47]
Soybean49 (65.3%)9 (12.0%)12(2.7%)25(20.0%)72[48]
Grapevine20 (69.0%)4 (13.8%)1 (3.4%)4 (13.8%)29[49]
Tomato25 (65.8%)8 (21.1%)1 (2.6%)4 (10.5%)38[50]
Tobacco36 (69.2%)6 (11.5%)3 (5.8%)7 (13.5%)52[51]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yin, Z.; Zhou, F.; Chen, Y.; Wu, H.; Yin, T. Genome-Wide Analysis of the Expansin Gene Family in Populus and Characterization of Expression Changes in Response to Phytohormone (Abscisic Acid) and Abiotic (Low-Temperature) Stresses. Int. J. Mol. Sci. 2023, 24, 7759. https://doi.org/10.3390/ijms24097759

AMA Style

Yin Z, Zhou F, Chen Y, Wu H, Yin T. Genome-Wide Analysis of the Expansin Gene Family in Populus and Characterization of Expression Changes in Response to Phytohormone (Abscisic Acid) and Abiotic (Low-Temperature) Stresses. International Journal of Molecular Sciences. 2023; 24(9):7759. https://doi.org/10.3390/ijms24097759

Chicago/Turabian Style

Yin, Zhihui, Fangwei Zhou, Yingnan Chen, Huaitong Wu, and Tongming Yin. 2023. "Genome-Wide Analysis of the Expansin Gene Family in Populus and Characterization of Expression Changes in Response to Phytohormone (Abscisic Acid) and Abiotic (Low-Temperature) Stresses" International Journal of Molecular Sciences 24, no. 9: 7759. https://doi.org/10.3390/ijms24097759

APA Style

Yin, Z., Zhou, F., Chen, Y., Wu, H., & Yin, T. (2023). Genome-Wide Analysis of the Expansin Gene Family in Populus and Characterization of Expression Changes in Response to Phytohormone (Abscisic Acid) and Abiotic (Low-Temperature) Stresses. International Journal of Molecular Sciences, 24(9), 7759. https://doi.org/10.3390/ijms24097759

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop