Next Article in Journal
The Gut–Immune–Brain Axis: An Important Route for Neuropsychiatric Morbidity in Inflammatory Bowel Disease
Previous Article in Journal
The Role of Fatty Acid-Binding Protein 4 in the Characterization of Atrial Fibrillation and the Prediction of Outcomes after Catheter Ablation
Previous Article in Special Issue
The PagKNAT2/6b-PagBOP1/2a Regulatory Module Controls Leaf Morphogenesis in Populus
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Characterization and Expression Analysis of Fatty acid Desaturase Gene Family in Poplar

1
Key Laboratory of Landscape Plant Genetics and Breeding, School of Life Sciences, Nantong University, Nantong 226001, China
2
College of Biology and the Environment, Nanjing Forestry University, Nanjing 210037, China
3
College of Arts and Sciences, Arlington International University, Wilmington, DE 19804, USA
4
Department of Horticulture, Faculty of Agriculture and Natural Resources, University of Hormozgan, Bandar Abbas 7916193145, Iran
5
Department of Animal Science, Faculty of Agriculture, Yasouj University, Yasouj 7591874831, Iran
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2022, 23(19), 11109; https://doi.org/10.3390/ijms231911109
Submission received: 25 August 2022 / Revised: 15 September 2022 / Accepted: 16 September 2022 / Published: 21 September 2022

Abstract

:
Fatty acid desaturases (FADs) modulate carbon–carbon single bonds to form carbon–carbon double bonds in acyl chains, leading to unsaturated fatty acids (UFAs) that have vital roles in plant growth and development and their response to environmental stresses. In this study, a total of 23 Populus trichocarpa FAD (PtFAD) candidates were identified from the poplar genome and clustered into seven clades, including FAB2, FAD2, FAD3/7/8, FAD5, FAD6, DSD, and SLD. The exon–intron compositions and conserved motifs of the PtFADs, clustered into the same clade, were considerably conserved. It was found that segmental duplication events are predominantly attributable to the PtFAD gene family expansion. Several hormone- and stress-responsive elements in the PtFAD promoters implied that the expression of the PtFAD members was complicatedly regulated. A gene expression pattern analysis revealed that some PtFAD mRNA levels were significantly induced by abiotic stress. An interaction proteins and gene ontology (GO) analysis indicated that the PtFADs are closely associated with the UFAs biosynthesis. In addition, the UFA contents in poplars were significantly changed under drought and salt stresses, especially the ratio of linoleic and linolenic acids. The integration of the PtFAD expression patterns and UFA contents showed that the abiotic stress-induced PtFAD3/7/8 members mediating the conversion of linoleic and linolenic acids play vital roles in response to osmotic stress. This study highlights the profiles and functions of the PtFADs and identifies some valuable genes for forest improvements.

1. Introduction

Fatty acids (FAs), the hydrophobic interior of cell components, are ubiquitous in plants and participate in many physiological processes [1]. The form of FA in a cell is mainly bound and rarely exists in a free state. Their biosynthesis and regulation play an important role in the basic metabolic activities of cells [2]. The FAs can be divided into saturated FAs (SFAs) and unsaturated FAs (UFAs), according to the presence or absence of a double bond of FAs, and the UFAs are classified into the monounsaturated FAs (MUFAs), with one double bond, and the polyunsaturated FAs (PUFAs), containing two or more double bonds [3]. Also, based on the position of the unsaturated bond, the UFAs can be divided into a series of ω-3, ω-6, ω-7, and ω-9 UFAs [4]. The UFAs are essential for bio-membranes and improve membrane fluidity [5,6]. FA desaturases (FADs) can catalyze single-bond UFAs into double-bond UFAs in fatty acyl chains [7]. Two categories of FADs have been identified in plants; one is the soluble FAD, broadly distributed on plant cell plastids, and the other is the membrane-bound FAD, mainly localized on the endoplasmic reticulum (ER) plastid membranes [8]. It has been reported that the soluble FAD can be identified as stearoyl-ACP desaturase (FAB2 or SAD), which explicitly converts stearoyl-ACP (C18:0) into palmitoleic acid and oleic acid (C18:1). Then, the Δ12 FADs control the oleic acid to linoleic acid, and the Δ15 FADs or Δ6 FADs further convert the linoleic acid into α-linolenic acid or γ-linolenic acid, respectively [9].
There are two spatially independent pathways for glycerolipid biosynthesis in advanced plants, including the prokaryotic and eukaryotic pathways. The prokaryotic pathway is localized in the chloroplast, while the eukaryotic pathway refers to the glyceride biosynthesis in the ER [10]. The glycolytic pathway provides acetyl-coenzyme A (CoA) as the FA precursor [11]. Acetyl-CoA carboxylase (ACC) converts two molecules of the acetyl-CoA into malonyl-CoA, and then the malonyl-CoA can be catalyzed to form the acyl carrier protein (ACP). The C18:0-ACP as the substrate is catalyzed to generate phosphatidylglycerol (PG), and the PG can be converted to diacylglycerol (DG) with the catalyzation phosphatase. For example, the FAD5-8 modulate the PG, sulfoquinovosyl diacylglycerol (SQDG), and digalactosyldiacylglycerol (DGDG) to form C18:3 linolenate [12]. In the eukaryotic pathway, the C16:0-ACP and C18:1-ACP can be catalyzed by the acyl-CoA synthase to form C16:0-CoA and C18:1-CoA, and they are transferred by the acyl-CoA binding protein (ACBP) to the ER for the phosphatidylic acid (PA), phosphatidylcholine (PC), and another glyceride biosynthesis [13]. Therefore, the C18 UFAs are synthesized by the prokaryotic and eukaryotic pathways, while the C16 UFAs can be generated by the eukaryotic pathway [10].
The studies on desaturases in large numbers of plant species have revealed their characterizations and functions. For example, the glycerolipid n-3 desaturase has been revealed to be responsible for lipid composition and has a significant function in the biosynthesis of the trienoic FAs (TFAs), including the C16:3 and C18:3 FAs [14]. The Arabidopsis omega-3 (ω-3) FAD has been demonstrated to be limited to linolenic acid (C18:3) production in seeds [15]. A mutation in the FAD7 in Arabidopsis reduced the TFA contents, and a mutation at the fad8 locus of the Arabidopsis showed lower levels of phospha-tidylglycerol and sulfoquinovosyldiacylglycerol under a low temperature [16]. The Arabidopsis FAD7 and FAD8 were identified as functionally related, and the constitutive expression of the FAD8 in the mutation of the FAD7 plants resulted in a genetic complementation of the mutation [17]. Shah et al. [18] proved that up-regulation of the FAD3 transcript levels leads to increasing levels of the C18:3 FAs and correspondingly decreasing levels of the C18:2 FAs. The FAD2, a plant oleate desaturase, is responsible for PUFA biosynthesis, especially linoleate and linolenate [19]. The FAD2 and FAD3 are exposed on the cytosolic side of the ER membranes and catalyze the conversion of the C18:1 PC-bound oleate acid to the C18:2 linoleic acid and C18:3 linolenic acid, respectively [20]. In addition, UFAs have been identified to play an essential role in responding to abiotic stresses, such as low and high temperature, salt, drought, and heavy metals [21,22]. An increased UFA content in the chloroplast membrane enhanced the Arabidopsis tolerance to chilling, freezing, and oxidative stress [23]. The Oryza sativa FAD2 (OsFAD2) participated in the FA desaturation and promoted the tolerance of rice to low-temperature stress [24]. The Brassica napus FAD3 (BnFAD3) and Arabidopsis thaliana FAD8 (AtFAD8) were proved to promote the linolenic acid/linoleic acid ratio and resistance to osmotic stress [25]. The ratio of linolenic acid to linoleic acid was significantly increased when the cowpea was treated with drought stress [26]. Moreover, FADs are crucial for regulating UFAs, which may be associated with limiting the compositions and contents of the UFAs, especially linoleic and linolenic acids. Previous studies have indicated that FAD genes are essential for plant growth and development and the sustenance of plants resilience to different environmental stresses. The FAD3, FAD7, and FAD8 had close associations with TFAs, the precursor of oxylipin, which modulated Arabidopsis pollen development [27]. The post-translational stability of Arabidopsis FAD8 was considered an essential regulatory mechanism in response to temperature, modulating the TFA contents at high temperatures [28]. The overall mRNA levels of Zea mays FAD7/8 (ZmFAD7 and ZmFAD8) were elevated upon exposure to low temperature, and the ZmFAD7/8 expression was accumulated in the roots when the maize was treated with high-salt stress [29]. A low temperature induced the soybean GmFAD3A gene, and overexpression of it in rice promoted the cold resistance and seed germination rate [30,31]. The AtFAD3 or AtFAD8 was ectopically overexpressed in tobacco, increasing the osmotic stress resistance [9]. The tomato FAD3 has been reported to have crucial roles in resistance to salt stress, while down-regulation of the tomato FAD7 expression dramatically improved its tolerance under high-temperature stress [32]. The Arabidopsis FAD2 and FAD6 were identified to have a positive function in response to salinity stress [33,34]. Also, the mutant of ads2 enhanced sensitivity to chilling and freezing treatments [35], and the mutant of sld showed an increased sensitivity to low-temperature stress [36]. In tobacco, the expression levels of FAD3 and 8 were promoted under drought stress, and the FAD7 overexpression significantly improved the drought tolerance [25]. Thus, the FADs, as essential regulatory genes, can respond to various stresses and are associated with alleviating the damage induced by various stresses.
Plant genomes contributed to the genome-wide analysis of the FAD gene family. In Arabidopsis and rice, 24 and 18 FAD genes were identified, including soluble and membrane-bound FADs [8]. Also, 68 putative FAD members were identified in Brassica napus; 23 FAD genes were retrieved from the Cucumis sativus genome; and 19 FAD members were found in Gossypium raimondii [37,38,39]. In addition, 20 full-length FAD members were identified from Medicago truncatula, and 30 putative FADs were found in the maize genome, including 17 membrane-bound and 13 soluble FADs [13,40]. For Triticum aestivum, 68 FAD genes were retrieved in the wheat genome, and 29 FAD genes were found in the Cucumis sativus genome [41,42]. Poplar, a major source of natural material for industrial production, is widely cultivated worldwide. Here, we identified the genome-wide identification of the Populus trichocarpa FAD members and divided them into two subfamilies, or seven clades, based on the phylogenetic tree. Also, we investigated the chromosome distribution, collinearity relationship, selection pressure, gene structure, and cis-acting elements of the PtFAD members. In addition, we evaluated the expression patterns of the PtFAD genes and UFA contents under abiotic stresses. Those results may lay a critical basis for understanding the biological roles of FADs in response to abiotic stresses and may also contribute to accelerating the process of the genetic improvement of poplar or other species.

2. Results

2.1. Genome-Wide Identification of FAD Members in P. trichocarpa

A BLASTP analysis of the P. trichocarpa genome was executed to investigate the FAD-related proteins. Thirty-two putative FAD genes were identified from the P. trichocarpa genome, and Potri.T034100 and Potri.T180100 were eliminated from the FAD candidates due to their localization on the scaffolds (Supplemental Table S1). The Pfam, including (PF00487) and (PF03405), were also used to identify the PtFAD candidates, and the 23 FAD candidates were identified from the P. trichocarpa genome, when the incomplete PtFADs (Potri.011G137600 and Potri.018G107700) were eliminated (Supplemental Table S1). In conclusion, 23 FAD genes were identified and renamed as PtFAD2s, PtFAD3/7s/8s, PtFAD5s, PtFAD6s, P. trichocarpa delta-8-fatty acid desaturase (PtSLD1/2), P. trichocarpa dihydroceramide desaturase (PtDSD1/2), and PtFAB2s, based on their similarity with AtFADs and OsFADs (Supplemental Table S2). The PtFAD genes ranged from 984 bp (PtFAD2.1) to 1359 bp (PtFAD8.2) in length, and the corresponding amino acids ranged from 327 aa (PtFAD2.1) to 452 aa (PtFAD8.2) in length. The molecular weights of these PtFADs ranged from 38.12 kDa (PtFAD2.1) to 51.81 kDa (PtFAD8.2), with the isoelectric points varying from 5.62 (PtDSD1) to 9.60 (PtFAD5.2) (Supplemental Table S3). Subcellular localization predictions of the PtFADs revealed that they are localized on the chloroplast and ER, and the spatial diversity of these proteins is likely associated with the various functions (Supplemental Table S3).

2.2. Comparative Alignment and Phylogenetic Analysis

A multiple alignment analysis was performed to characterize the conserved domains, and the most distinguished amino acids in the conserved domains, and to determine the conserved domains of the PtFADs. The result of the multi-alignment revealed that the membrane-bound PtFADs contained three relatively conservative histidine boxes (HXXXH, HXXHH, and HXXHH) (Supplemental Figure S1). In addition, to determine the evolutionary relationship of the FADs in different species, including A. thaliana (At), O. sativa (Os), P. trichocarpa (Pt), and T. aestivum (Ta) [42], a phylogenetic tree constructed by the neighbor-joining (NJ) method was used to assess the genetic relationship among the PtFADs, AtFADs, OsFADs, and TaFADs. As shown in Figure 1, all the investigated FADs in this study were classified into two subfamilies: the membrane-bound FADs and soluble FADs (FAB). The membrane-bound FADs were divided into six clades (FAD2, FAD3/7/8, FAD5, FAD6, DSD, and SLD), and the FAB2 subfamily was significantly separated from the membrane-bound FAD subfamily. The phylogenetic tree indicated that the FAB2s and FADs might exist in a common ancestor before the divergence of the monocotyledons and dicotyledons.
The FAD2 clade thought of as ω-6 desaturases included six members from P. trichocarpa, two members from A. thaliana, four members from O. sativa, and fifteen members from T. aestivum. The FAD6 was also considered as a kind of ω-6 desaturases, and the FAD6 clade included one member from A. thaliana, one member from O. sativa, and two members from P. trichocarpa and T. aestivum, respectively (Figure 1). The ω-3 desaturases included FAD3, FAD7, and FAD8. Here, the PtFAD3 included one member, the PtFAD7 comprised PtFAD7.1/7.2, and the PtFAD8 contained PtFAD8.1/8.2 (Figure 1). The SLD clade catalyzed the desaturation of the sphingolipid at the ∆8 position, comprising three members from P. trichocarpa, two members from A. thaliana, two members from O. sativa, and three members from T. aestivum (Figure 1). The sphingolipid ∆4 desaturase DSD comprised two members from P. trichocarpa, one member from A. thaliana, one member from O. sativa, and three members from T. aestivum (Figure 1). In addition, the FAD5 (ADS) clade comprised nine members from A. thaliana and only two from P. trichocarpa (Figure 1). Also, no ADS members were identified from rice and T. aestivum, suggesting they may have experienced evolutionary losses in the O. sativa and P. trichocarpa. In addition, most PtFAB2 members were clustered together with the AtFAB2 members, and the PtFAB2, AtFAB2, OsFAB2, and TaFAB2 members were clustered together, indicating that the FAB2s were not significantly divided into monocotyledonous and dicotyledonous groups and may have independent evolutions in these plant species.

2.3. Chromosome Mapping and Duplication Events of P. trichocarpa FAD Genes

The chromosomal location of the 23 poplar FAD genes revealed their distributions on the 12 chromosomes, except for Chr03, Chr07, Chr09, Chr12, Chr15, and Chr17, where they were not present (Supplemental Figure S2). The PtFAD genes were mapped unevenly on the P. trichocarpa chromosomes. The Chr01, Chr10, and Chr11 had three FAD genes, whereas four genes were located on the Chr06. The Ch08 and Chr16 had two FAD genes, and Ch02, Ch04, Ch05, Ch13, Ch14, and Ch18 each had one FAD gene.
The diversity and expansion of the gene families resulted from tandem and segmental duplication events and we investigated the gene duplication within the PtFAD genes. The 16 segmental duplication events were determined among the PtFAD family members: PtFAD2.1/PtFAD2.4, PtFAD2.1/PtFAD2.3, PtFAD2.3/PtFAD2.4, PtFAD3/PtFAD7.1, PtFAD3/PtFAD8.2, PtFAD3/PtFAD8.1, PtFAD8.1/PtFAD7.1, PtFAD8.1/PtFAD7.2, PtFAD8.1/PtFAD8.2, PtFAD7.1/PtFAD7.2, PtFAD7.1/PtFAD8.2, PtFAD7.2/PtFAD8.2, PtFAD6.1/PtFAD6.2, PtSLD1/PtSLD3, PtDSD1/PtDSD2, PtFAB2.2/PtFAB2.3 (Figure 2), and PtFAB2.3/PtFAB2.4 as tandem duplication was identified in the syntenic PtFAD pairs. Therefore, the segmental duplication events occupied prominent roles in increasing the genetic diversity of the poplar FAD family gene. To understand the driving pressure for the PtFAD family gene evolution, the Ka/Ks values of the syntenic gene pairs were calculated. The values of Ka/Ks among the PtFAD gene pairs were less than one (Supplemental Table S4), which indicated that the PtFAD members have mainly experienced negative selection. The slow evolutionary process contributes to keeping the conserved functions of the PtFAD members.
To further investigate the evolutionary process among the dicotyledonous and monocotyledons, the duplicated gene pairs among the PtFAD, AtFAD, OsFAD, and the Salix purpurea FAD (SpFAD) members were identified (Figure 3). A total of 14 duplicated gene pairs between P. trichocarpa and A. thaliana were determined (Supplemental Table S5). Also, the following 14 gene pairs between P. trichocarpa and O. sativa were found (Supplemental Table S5). In addition, 28 collinearity gene pairs between P. trichocarpa and S. purpurea were found, 3 of which belonged to the syntenic FAB gene pairs, and the rest were membrane-bound FAD gene pairs (Supplemental Table S5). Although there is no identification and characterization of the S. purpurea FAD family members, the collinearity gene pairs might lay a foundation for studying the SpFAD family gene. There were more collinearity gene pairs between P. trichocarpa and S. purpurea, which might be associated with both poplar and willow, belonging to the Salicaceae. The various PtFAD genes have syntenic relationships with three to five willow genes, which suggests that these genes may play essential roles in the evolution of the FAD family gene.

2.4. Exon–Intron Distribution, Conserved Motif, and 3D Structure Analysis

The FAD gene structures were analyzed to gain insight into the exon–intron distributions of the FADs. The 47 FAD genes were clustered into seven clades, with highly similar patterns of gene structures identified for the PtFAD and AtFAD genes clustered into the same clade (Figure 4). For the FAD3/7/8 cluster, the FAD3/7/8 members contained eight exons and seven introns (Supplemental Table S6). Also, the PtFAD3/7.2/8.1 has a longer UTR than other FAD3/7/8 members. All the FAD2 members contained one exon except PtFAD2.1, which had two exons (Supplemental Table S6). In the FAD6 clade, the PtFAD6.1 and 6.2 had eight exons and seven introns, while the AtFAD6 possessed ten exons and nine introns (Supplemental Table S6). For the FAD5 clade, all the FAD5 contained five exons and four introns, and the FAD5 members were more dominantly similar in their exon–intron distributions (Supplemental Table S6). Specifically, the AtDSD1 and PtDSD2 had two exons and an untranslated region (UTR), while the PtDSD1 contained three exons and was not composed of a UTR. In general, the SLD members only contained a single exon and no intron. Most of the FABs harbored three exons in the FAB subfamily, except the PtFAB2.1 and AtFAB2.3/2.4, which had two and four exons, respectively (Supplemental Table S6). Similar gene structures of the FADs were identified in the same clade, suggesting the exon–intron distribution supports the phylogenetic classification of the FADs.
The Multiple Em for Motif Elicitation (MEME) server was used to identify the conserved motifs of the membrane-bound FADs and FAB2s. According to the results, the motifs 1–10 were detected in the FAD family members. Also, the evaluation of motifs 1–10 with the Pfam indicated that motifs 1–7 and motif 10 are related to the FAD gene family (Supplemental Figure S3 and Supplemental Table S7). The amino acids of ten motifs ranged from 21 to 50, while the sites ranged from 12 to 47 (Supplemental Table S7). The FAB2s contain motifs 1, 3, 5, 6, 7, and 10, except for the AtFAB2.4, which harbored the additional motif 2 (Figure 5). The DSD clade had only one motif (motif 3), and the FAD5 and SLD contained two motifs, 2 and 3, except the AtSLD1, which contained motif 4. The clades FAD3/7/8 and FAD2 contained motifs 2, 3, 4, 8, and 9, which indicated that these clades have a close phylogenetic relationship (Figure 5). For the FAD6 clade, motifs 2, 3, and 4 were harbored in the PtFAD6s, while the AtFAD6 was composed of motifs 2, 3, 4, and 8. There is a dominant divergence in the conserved motifs between the membrane-bound FAD and FAB2 (Figure 5). For instance, only motif 3 consensually presented in all the FAD members, while motifs 1, 5, 6, 7, and 10 existed only in the FAB2, suggesting the FAD family members may be involved in two directions in the process of evolution; one is the FAB subfamily, and the other is the membrane-bound FAD subfamily.
The FAD members have been documented to be associated with environmental stresses by introducing a double bond to the formation of UFAs. In order to understand the feathers of the PtFAD protein structures, the 3-dimensional (3D) structures of the PtFADs were generated using the SWISS-MODEL. The 3D structures of the PtFADs and AtFADs showed that the PtFADs and AtFADs are composed of alpha helices, random coils, beta-strands, and the PtFADs and AtFADs clustered into the same clade share similar 3D structures, which provided insight into the function analysis of the PtFAD proteins (Supplemental Figure S4).

2.5. Detection of FAD Cis-Acting Elements

To discover the regulation mechanisms of the FAD member expression patterns, the PlantCARE website was used to investigate the 2 kb upstream promoters of the PtFAD and AtFAD genes. Various cis-acting elements were found, and each FAD gene was composed of divergent cis-acting elements (Supplemental Figure S5). Three different cis-elements were identified in the FAD gene promoters. The first is associated with the plant growth and development, including the light-responsiveness, seed-specific regulation, circadian control, differentiation of the palisade mesophyll cells, and cell cycle regulation. The second is the stress response, including the MeJA-responsiveness, abscisic acid-responsiveness (ABRE), wound-responsiveness, auxin-responsiveness, low-temperature responsiveness, defense and stress-responsiveness, SA-responsiveness, GA-responsiveness, and drought-inducibility. These are related to specific biological processes, including the elements involved in the endosperm expression, zein metabolism regulation, flavonoid biosynthesis regulation, and cell cycle regulation. Most of the PtFAD and AtFAD promoters harbored the elements involved in the light and hormone responses. For example, the PtFAD3s were abundant with the cis-acting elements of hormone responsiveness. Also, the cis-acting elements involved in flavonoid biosynthesis regulation were identified in five AtFADs (AtADS1, AtADS7, AtADS5, AtADS4, AtFAD6), and AtFAB2.2 and one PtFAD (PtDSD2) promoter, showing that these genes may play essential roles in flavonoid biosynthesis. The promoters of the AtDSD1, PtFAD2.1, and PtFAD6.2 contained cis-acting elements participating in the seed-specific regulation, suggesting that these genes may have close associations with the molecular regulation of seed development. The wound-responsive element was detected in the AtFAD7, PtDSD2, and PtFAD2.4, indicating that these genes might be mainly involved inwewound response. In summary, various kinds of the cis-acting elements were identified in the PtFAD and AtFAD genes, indicating that the PtFAD and AtFAD genes may participate in responding to some stimulus.

2.6. Interaction Network and Gene Ontology (GO) Annotation of FADs

To better understand the characterizations and profiles of the FAD genes, interaction networks and GO annotation were performed to gain more insight into the putative function of the FADs. The results showed that 70 Arabidopsis proteins interacted with 24 AtFADs (Supplemental Figure S6). For example, as an SA receptor, NPR1 (non-expressor of pathogenesis-related genes 1) plays an important role in regulating plant disease resistance [43]. KASI, beta-ketoacyl-ACP synthase I, significantly participates in C4–C16 FA biosynthesis and lipid metabolism [44]. CYP38 (Cyclophilin 38) plays an essential role in the assembly and stability of the photosynthetic system II [45]. In addition, 23 PtFADs were predicted to interact with 59 poplar proteins (Supplemental Figure S7). For example, FAB1, phosphatidylinositol-3-monophosphate5-kinases, participates in vacuole/lysosome homeostasis and transporting various proteins to the vacuole [46]. ACC, acetyl-CoA carboxylase, is involved in ethylene biosynthesis [47]. SAY1, steryl acetyl hydrolase 1, plays a part in lipid metabolism and significantly influences lipid homeostasis [48]. AOX, allene oxide synthase, catalyzes the UFAs to generate important JA precursors. In conclusion, the interaction network analysis of the AtFADs and PtFADs predicted that these proteins were closely associated with the FA and lipid metabolisms and participated in various stress responses.
A GO annotation analysis was performed to improve understanding of the pathway involved in the FAD genes. The results indicated that the AtFADs, OsFADs, and PtFADs were dominantly enriched in divergent GO terms, mainly “molecular function” and “biological process” (Supplemental Figure S8). In the category of “biological process”, the AtFADs, OsFADs, and PtFADs were mainly enriched in the lipid metabolic process (GO:0006629), obsolete oxidation-reduction process (GO:0055114), fatty acid metabolic process (GO:0006631), and monocarboxylic acid metabolic process (GO:0032787). In the category of “molecular function”, the AtFADs, OsFADs, and PtFADs were significantly involved in the oxidoreductase activity (GO:0016717, GO:0016705, and GO:0016491), and acyl–acyl-carrier-protein desaturase activity (GO:0045300). It is suspected that the FADs participate in diverse biological processes, including the FA and lipid metabolisms, and occupy various molecular functions containing the redox and desaturation activities. In addition, the protein interaction relationship contributes to exploring the putative signaling pathway or physiological process of the FADs. Therefore, a GO annotation analysis was performed to analyze the FAD and interaction-related protein genes. As shown in Figure 6, the PtFAD and interaction-related protein genes were mainly involved in three GO terms: “molecular function”, “cellular component”, and “biological process”. In the category of cellular components, the PtFAD and interaction-related protein genes significantly enriched the integral components of the membrane (GO:0016021 and GO:0031224), which suggested that the PtFADs may also be associated with the assembly stability of the membrane.

2.7. Expression Profiles of Poplar FAD Genes in Diverse Tissues

The FADs introduce a double bond into the FAs acyl chain leading to the UFAs involved in plant growth and development. To investigate the functional role of the PtFAD genes, a real-time quantitative PCR (RT-qPCR) was used to analyze the expression profiles of the diverse tissues, including the mature leaves (ML), young leaves (YL), upper region of the stems (US), lower region of the stems (LS), and roots. The results showed that each PtFAD gene has one specific expression pattern.
In P. trichocarpa, the relative mRNA levels of the PtFAD8.1/8.2 and PtFAB3.5 were lower and absent from the roots. The relative expression of the PtFAD3/7.1/7.2 was comparatively higher in the YL and lower in the ML. The PtFAB2.1/2.2 was highly expressed in the roots and stems, while lower accumulations were identified in the ML. The PtFAB2.3 occupied the highest expression levels in the LS and the lowest in the ML, while the highest expression level of the PtFAB2.5 was found in the YL and the lowest in the UP (Figure 7A). In ‘Nanlin 895’ (P. deltoides × P. euramericana), the transcript patterns of the PtFAD7.1/8.1 showed a high similarity, with the highest expression of abundance in the UP and the lowest in the roots. Also, the high expression levels of the PtFAB2.1/2.5 were observed in the roots, while low expression abundances were detected in the stems. The PtFAD3/8.2 genes had a high abundance of expression in the US, while the lower expression levels were found in the ML and roots, respectively.
Interestingly, the PtFAB2.2/2.3 showed lower expression levels in all the tissues tested, whereas the PtFAB2.4 showed higher expression accumulations in the YL, US, and roots (Figure 7B). In ‘Shanxinyang’ (P. davidiana × P. bolleana Loucne), similar expression patterns were identified in the PtFAD7.2/8.2, with relatively higher expression accumulations in the YL and lower in the roots. Also, the PtFAD3/7/8 and PtFAB2.2/2.3/2.5 genes have similar expression profiles and are highly expressed in the YL. In addition, the PtFAB2.1 displayed a specific expression profile in a more limited root, while the PtFAB2.1 had a higher expression level in all the tissues tested except the ML (Figure 7C). The expression analysis suggested that the PtFAD members play a key role in the growth and development of the poplar tissues and their functional divergences.

2.8. PtFAD Expression Profiles in Response to Multiple Abiotic Stresses

The external environment greatly influences plant growth and development, and the plant can respond to environmental stress by regulating crucial molecular processes. Various stress-related genes help plants to respond to various stresses under adverse environments. Thus, to investigate whether the PtFAD members play a role in abiotic stress and understand the putative regulation mechanisms, a RT-qPCR was applied to analyze their expression patterns under abiotic stress. As shown in Figure 8A, the expression levels of the PtFADs were significantly changed under the abscisic acid (ABA) treatment, and their expression patterns were classified into two groups. The expression levels of the PtFAD8.2 and PtFAB2.2/2.3/2.5 showed long-lasting and continuous improvement under the ABA treatment. The peak of the PtFAD7.1 transcript level was exhibited at 1 h, while the peak of the PtFAB2.1 expression was identified at 6 h. In addition, the highest expression abundances of the PtFAD3/8.1 were exhibited at 24 h after the ABA treatment. These findings suggested that the PtFAD8.2 and PtFAB2.2/2.3/2.5 play crucial functions in the ABA signal transduction and promote the poplar’s response to environmental stresses.
For the H2O2 treatment, the PtFAD3/8.1/8.2 and PtFAB2.2 displayed peak expression at the 48 h time point. The PtFAB2.4 and PtFAB2.5 had similar trends within the whole time of the H2O2 treatment, with higher expression levels within the time course. The PtFAB2.1 and PtFAB2.3 expression levels were promoted and reached the highest at 12 h and 24 h, respectively. The mRNA level of the PtFAD7.2 was upregulated at all time points, while the PtFAD7.1 expression level at 24–48 h was significantly lower than that at 0–12 h (Figure 8B). The response of the PtFADs to the NaCl treatment was more dominant than to the ABA treatment. All the tested PtFAD expression levels were significantly upregulated during the NaCl treatment, especially the PtFAD3/7s/8s and PtFAB2.3/2.5. The PtFAB2.4 reached its expression peak at 24 h, and the PtFAB2.1 and 2.2 were promoted significantly at 6–12 h (Figure 8C). In response to the PEG6000 treatment, the upregulated expression of the PtFAD3/7.1/8s and PtFAB2.1-2.3, and the PtFAB2.5 was observed, and their expression levels after the PEG6000 treatment were higher than that before the PEG6000 treatment. Also, the PtFAB2.4 expression was down-regulated at 0–24 h, and its expression level was lower than that of the control (Figure 8D). These results fully reflect that some PtFADs can respond to abiotic stress in poplar, which provides insight into the adaptability promotion of poplar in different environments.

2.9. Cluster Network Analysis of PtFADs and PtLACSs

Under the ABA treatment, the cluster dendrogram showed that the PtFAD and PtLACS expression patterns were divided into two categories (Figure 9A). The PtFAD3/7.2/8s, PtFAB2.2/2.3/2.5, and PtLACS6/7/9s clustered together in the same category, and the PtFAD2.3/7.1, PtFAB2.1/2.4, and PtLACS1-4/8.1 clustered together in the other category. Also, a cluster network analysis of the PtFAD and PtLACS expression levels under the H2O2 treatment classified them into two categories (Figure 9B). One is the small cluster dendrogram, including the PtFAD2.3/7.1 and PtLACS3/4, and the other is the big cluster dendrogram, including most of the PtFAD and PtLACS members. For example, the PtFAB2.4/2.5 and PtLACS7/9.1 clustered together, suggesting they may be involved in the LCFA biosynthesis and desaturation. The cluster dendrogram was also divided into four categories (Figure 9C). The PtFAD2.3/3/8.1 and PtLACS9.1, PtFAB2.3-2.5, PtFAD7.2, and PtLACS7/8.1, and the PtFAB2.1/2.2, FAD7.1/8.2, and PtLACS2/6/9.2 clustered together, respectively. The result of the cluster network on the PtFADs and PtLACSs expression under the NaCl treatment indicate that they play essential roles in response to NaCl stress by regulating the FA biosynthesis and desaturation. For the PEG6000 treatment, the dendrogram was clustered into two categories, one of which had the PtLACS1/2/9.1 and PtFAD2.3/2.4; the other contained the rest of the PtLACS and PtFAD members (Figure 9D). For example, the PtLACS9.1 clustered with the PtFAD2.2/7.1/8.1, indicating a close relationship in response to the PEG6000 treatment. The result of the cluster network showed that expression patterns of some PtFADs under the abiotic stress are consistent with those of some PtLACSs, indicating that the PtFADs and PtLACSs involved in the LCFA biosynthesis and desaturation are related to the response to abiotic stress.

2.10. UFAs Analysis before and after Osmotic Stress

The nonadecanoate (C19:0) was considered as an internal standard, and the quantitative investigation of the single compounds was quantified relative to an internal standard through the automatic integration of the peak areas [49]. Also, the Chromeleon 7.0 was performed to analyze the characteristic peaks of the individual components of the SFA and UFA. In addition, the online tool NIST 17 (https://www.sisweb.com/software/ms/nist.htm (accessed on 8 April 2022) database) was used to sort the data of the SFA and UFA in poplars. The contents of the SFA and UFA were dominantly changed by the drought and salt stresses. The contents of the SFA, including myristate, pentadecanoate, palmitate, margarate, octadecanoate, icosanoate, behenate, and lignocerate, were significantly down-regulated when the poplars were treated by the drought and salt stresses. Also, under the drought and salt treatments, the contents of the UFAs, including the cis-9-tetradecenoate, cis-10-pentadecenoate, palmitoleate, cis-10-heptadecenoate, linoleic acid, cis-8,11,14-eicosatrienoate, cis-11-eicosenoate, and nervate, were significantly decreased.
In contrast, the contents of the UFAs, including oleate, linolenic acid, and erucate, were dominantly promoted by the drought and salt stresses (Figure 10). The previous studies proved that the oleic acid and ratio of linolenic acid/linoleic acid are closely associated with the osmotic tolerance and are considered as osmotic-regulating substances against osmotic stress [26,39]. The contents of the oleic acid in poplars cultivated under normal conditions were higher than that in poplars cultivated under the drought and salt treatments (Figure 10). Also, the poplars grown under normal conditions produced total linoleic acid levels 55- or 27-fold higher than those treated with the drought or salt treatment (Figure 10). In addition, the linolenic acid contents were significantly increased when the poplars were treated with the drought and salt treatments (Figure 10), and the value of linolenic acid/linoleic acid in poplars treated by the drought or salt treatment was higher than that under the normal condition. The above results suggested that the UFAs participate in osmotic stress, especially the oleic acid, linoleic acid, and linolenic acid. It is speculated that the accumulations of oleic acid, linoleic acid, and linolenic acid improve the membrane lipid fluidity, which can help the poplar respond to abiotic stress.

3. Discussion

The UFA compounds and contents are considered the key factor for plant growth and development, and the FADs play an essential role in regulating the UFA biosynthesis [50]. The previous studies on the FADs were constrained to specific plants, and not much information was available on the FADs families in poplars [8,13,42]. Given the vital function of the FADs, systematic studies on poplar FADs were performed. In the present study, a total of 23 PtFADs were identified from the P. trichocarpa genome based on the Pfams and alignment with the AtFADs. There is no significant difference in the number of FAD genes among the poplar, Arabidopsis, and rice. However, the number of PtFADs was less than in cotton [51], mustard [52], and wheat [42]. A similar number of FAD genes in poplar, compared to Arabidopsis and rice, may result from their diploid genome. There are more FADs in cotton, mustard, and wheat than in poplar, Arabidopsis, and rice, possibly because of their ploidy. Therefore, the number of the FAD genes did not correspond to the genome size, which indicate that they experienced divergent duplication events during evolution. Gene expansions play an essential role in generating a family gene, and segmental and tandem duplications are commonly associated with the driving forces in family genes [53]. Consistently, the previous studies showed that 22 FAD gene pairs are found in rape, all of which are segmental duplication events [54]. Also, 26 and 126 TaFAD gene pairs belonging to tandem and segmental duplication were identified from the wheat genome, respectively [42]. In addition, 25 segmental duplicated gene pairs and three tandem duplicated gene pairs were in the BjFAD gene family [52]. Here, 17 collinearity gene pairs were determined among the PtFADs, containing 16 gene pairs of segmental duplication event and one gene pair of tandem duplication event, which indicate that the segmental duplication event occupies a prominent role in the PtFAD family gene expansion.
The phylogenetic analysis indicated that the FAD family was significantly divided into two subfamilies, including the soluble and membrane-bound FADs, similar to the previous studies [55]. Also, the membrane-bound FADs were further separated into the AD3/7/8, FAD2, FAD6, FAD5 (ADS), DSD, and SLD clades [42,55]. The rice did not contain the FAD5 clade, and the PtFAD5 clade only contained the PtFAD5.1/5.2, which was smaller than the number of the FAD5 in Arabidopsis. The soluble and membrane-bound FADs were predicted to be localized in diverse positions; the membrane-bound PtFAD members were localized in the ER, and the PtFAB members were chloroplast-localized, like other plants such as sunflower, sesame, and canola [55]. The exon/intron structure variation is related to gene evolution [56]. The PtFAD members clustered into the same clade and shared similar gene distributions and motif compositions, similar to other plant species [42,55]. The oil palm FAD promoters were identified to have some cis-acting elements involved in stress- and light-responsiveness, and their expression was regulated by low temperature and darkness [57]. Also, SA- and JA-responsive elements were identified in the BnFAD2-C5 promoters. Inconsistently, their expression was upregulated by the SA and JA inducible [38]. Both the SeFAD2 and BnFAD promoters contained the ABRE, and their transcript was induced by the ABA [58,59]. Here, some hormones- and stress-responsive elements were identified in the PtFAD promoters, indicating that the PtFADs expression may be affected by hormones and stresses. Also, many transcription factor binding sites, such as the MYB binding sites, were found in the PtFAD promoters, suggesting that these transcription factors regulate the PtFAD expression.
The divergence of the gene expression plays an important role in the family gene, and the PtFAD expression patterns provide an opportunity for exploring their profiles and functions. The mRNA levels of the FAD2/2-1 from Arachis hypogaea were accumulated in the seeds, and the expression levels of the FAD2-2/6 and SLD1 were elevated in the leaves [41]. The BjFAD2s were expressed in all the tissues tested, and the CaFAD2s were mainly accumulated in the flowers and seeds [60]. Some LuFAD2s, LuFAD3s, and LuFAB2s from linseed were highly expressed throughout all stages of the seed development, and the CsFAD genes were constitutively expressed in the cotyledons and leaves [37,61]. Here, the PtFADs were expressed in all the P. trichocarpa tissue tested, and the PtFAD3.1/3.3/3.4 expression was significantly accumulated in the YL. Also, the higher transcript levels of the PtFAB2.1/2.2 were found in the roots and stems. The previous studies revealed that various abiotic stresses affect poplar development and production, and the FADs are involved in the plant’s adaptation to various abiotic stresses [62]. Several BnFAD genes, such as the BnADS4.4/4.8/4.9, BnSLD8, and BnFAD8.1, potentially regulate the membrane adaptation to cadmium stress. Also, the BnADS4.9 expression level was remarkably improved, and the BnFAD7.4 and BnADS4.8 mRNA levels were decreased after the salt treatment [8]. In addition, the transcript levels of the maize FAD2.1/2.2 and SLD1–3 were upregulated under cold stress, while the ZmFAB2.1/2.3/2.12 expression was significantly down-regulated [13]. In this study, the expression patterns of the PtFADs were investigated under various abiotic stresses, and the results revealed that the PtFAD genes occupy different stress resistance for divergent stresses. For example, the transcript levels of the PtFAD8.2 and PtFAB2.2/2.3/2.5 were continuously induced under the ABA treatment; the PtFAD3/7s/8s and PtFAB2.3/2.5 were significantly upregulated within the course of the NaCl treatment; the dominant accumulations of the PtFAB2.4/2.5 expression were identified after the H2O2 treatment; and the PtFAD3/7s/8s and PtFAB2.1-2.3/2.5 expressions were significantly improved when the poplars were treated by the PEG stress. Thus, based on the above evidence, we speculated that the PtFAD expression levels have the dominant changes in response to various stresses and play an essential role in the plant’s stress resistance.
With the development of industrial production with wood as the raw material, the requirement for fast-growth and high-quality poplars has come to the fore again in recent years. Abiotic stress, such as drought and salt, has an adverse effect on tree growth and forest productivity [63]. Thus, the improvement of poplar quality contributes to expanding the planting area of poplars. The FAD catalyzes the formation of the UFAs, including the linolenic and linoleic acids, which are associated with the fluidity of the membrane lipids in response to the defense reactions of plants against environmental stresses. The previous studies have documented that the FADs play key roles in plant FA biosynthesis, development, and environmental stress tolerances [62]. The BnFAD3 were the dominant marker genes for a high oleic acid content, the AtFAD7/8 were associated with the chloroplast trienoic FA contents, and the overexpressing LuFAD3/7 promoted the α-linolenic acid content in transgenic rice and tobacco [64,65]. In addition, the overexpressing CrFAB2 in green alga significantly improved the contents of oleic and linoleic acids, and the overexpression of the cytosolic FAD3 and the plastic FAD8 remarkably upregulated the ratio of linolenic and linoleic acids, which changed the tolerance to various abiotic stresses [25,66]. In the olive fruit mesocarp, the transcript level of olive FAD6 was decreased; the oleic/linoleic ratio was promoted at maturation group 1; and the FA content was decreased at the end of the ripening period [67]. In general, the FAD3/7/8 clade members can convert (Δ9,12) linoleic acid to form (Δ9,12,15) linolenic acid [42]. FAB2, the first desaturase in the desaturation process of the FAs, is responsible for desaturating stearoyl-ACP to form oleoyl-ACP and determines the ratio of SFAs and UFAs [68]. In this study, the NaCl and PEG6000 treatments were considered as the osmotic stress, and the expression levels of the PtFAD3/7/8 and PtFAB2 members significantly increased after the NaCl and PEG6000 treatment. These results indicate that the PtFAD3/7/8 and PtFAB2 members could respond to osmotic stress in poplars and play an essential role in the adaptability promotion of poplars in an adverse environment. Also, the GC–MS showed that the value of linolenic acid/linoleic acid was significantly promoted. Following the integration of the results of the PtFAD3/7/8 and PtFAB2 expression patterns and the SFA and UFA contents, before and after osmotic stress, we proposed that the PtFAD3/7/8 and PtFAB2 have a close association with the osmotic stress response involving regulating the oleic acid content and the ratio of linolenic acid to linoleic acid, and in improving the membrane fluidity in response to osmotic stress. Therefore, based on the gene expression and UFA contents results, we speculated that the PtFAD3/7/8 and PtFAB2 could guide the UFA biosynthesis and modulate the SFA and UFA contents, especially the linolenic and linoleic acids, in response to osmotic stress. Considering the PtFAD3/7/8 and PtFAB2 expression profiles and the SFA and UFA contents, the putative mechanism of the PtFAD and PtFAB involved in the osmotic stress responses is proposed. When the poplar is under osmotic stress, the transcript levels of the PtFAD3/7/8 and PtFAB2 are promoted, and the activities of the PtFAD3/7/8 and PtFAB2 are increased. The increased accumulation of the linolenic acid/linoleic acid value accelerates the membrane fluidity, thus improving the resistance of the poplar to the osmotic stress.

4. Materials and Methods

4.1. Comparative Analysis of FADs from Arabidopsis, Rice, and Poplar

The genome databases of A. thaliana, O. sativa, P. trichocarpa, and S. purpurea were downloaded from Phytozome (http://www.phytozome.net/ (accessed on 8 April 2022)). The previous studies showed that 24 and 18 FAD members are identified from A. thaliana and O. sativa, respectively [8,42], and the AtFADs and OsFADs were obtained from the Arabidopsis Information Resource (TAIR, https://www.arabidopsis.org/ (accessed on 8 April 2022)) and the Rice Gene Annotation Project (http://rice.plantbiology.msu.edu (accessed on 8 April 2022)), respectively. For a further analysis of the poplar soluble- and membrane-bound FADs, the AtFADs and OsFADs, as the queries, were used to search the poplar proteome using the BLASTP program with the default parameters; the identifiers >40% and an Evalue ≤ 1 × 10−10 were used as the criterion. Also, the HMM profile of the FA desaturase (PF00487) and the FA desaturase 2 (PF03405) was used to verify the PtFAD candidates. The candidates were submitted to the NCBI Conserved Domains Database (https://www.ncbi.nlm.nih.gov/Structure/cdd/cdd (accessed on 8 April 2022)) to determine the conserved domain of the PtFAD members. The characterizations and basic information of the PtFADs were identified using the online ExPASy-ProtParam tool (https://web.expasy.org/protparam/ (accessed on 8 April 2022)), and the subcellular localization of the PtFADs was predicted using Cell-PLoc (http://www.csbio.sjtu.edu.cn/bioinf/Cell-PLoc-2/ (accessed on 8 April 2022)).

4.2. Multiple Alignment and Construction of a Phylogenetic Tree

Using the ClustalX software, the AtFADs, OsFADs, TaFADs, and PtFADs were conducted to align with the default parameters. The accession numbers of the PtFADs were renamed, based on their similarity with the AtFADs and OsFADs. Using a NJ method with 1000 bootstrap replicates, the MEGA7.0 was performed to construct the phylogenetic tree on the FADs from the poplar, Arabidopsis, and rice.

4.3. Comprehensive Analysis of Chromosome Locations and Gene Duplications

The general feature format version 3 (gff3) annotation file in the poplar genome database was used to identify the chromosome locations of the PtFAD genes. The chromosomal mapping image of the PtFAD genes was visualized by the TBtools software, based on their starting positions on the poplar chromosomes [69]. The previous studies have defined the criteria for gene duplication, including the length of the matched cover >80% and the identity of the matched regions >80% [8,70]. The tandem and segment duplications, considered the distinguished gene duplications, were performed to analyze the PtFAD members. To further assess the evolutionary pressure of the PtFAD members, the synonymous (Ks) and nonsynonymous (Ka) of the PtFAD gene pairs were calculated, using the TBtools with a simple Ka/Ks calculator. In addition, to investigate the evolutionary relationship among the various species, the MCScanX and BLASTP were performed to determine the gene pairs with a collinearity relationship of the FAD members among the poplar, Arabidopsis, rice, and willow.

4.4. Investigation of Gene Structures, Motif Distributions, and Protein Structures

The poplar coding sequence and genome file were applied to investigate the splicing phase of the AtFADs and PtFADs. The Gene Structure Display Server and TBtools software were used to visualize the structure diagrams of the AtFADs and PtFADs. Also, the MEME was considered for the conserved motifs of the FAD members, and the Pfam and SMART databases were performed to evaluate the function of the conserved motifs. The integrative structures of the AtFADs and PtFADs were constructed by the SWISS-MODEL based on the iterative template-based fragment. The models from the SWISS-MODEL were further refined and visualized by the Chimera software.

4.5. Promoter Analysis and Prediction of Protein Interaction and GO Annotation

The 2-kb upstream promoter sequences of the FADs were extracted, using the coding sequence and genome files from the poplar and Arabidopsis genome databases. The promoter sequences of the FADs were submitted to PlantCARE to identify the cis-regulatory elements. Also, the STRING database was applied to investigate the interaction networks of the AtFAD and PtFAD proteins. To further understand the putative involvement pathways of the PtFAD genes, a GO analysis was performed using the TBtools with a GO-basic file.

4.6. Plant Cultivation and Various Stress Treatments

All the poplar cultivars used in this study were taken from the Key Laboratory of Landscape Plant Genetics and Breeding at Nantong University. The poplar seedlings were cultivated in a greenhouse at 23 °C and 74% humidity. The poplar leaves with wounds were placed in a medium for regeneration, and then the poplar seedlings produced were transported to the shooting medium for a strong seedling. At the generation of the third or fourth leaves, the poplar shootings were selected and transferred to the rooting medium. The P. trichocarpa plants were cultivated in a woody plant medium (WPM) (pH 5.8), supplemented with 0.1 mg/L indole butyric acid (IBA). The ‘Nanlin 895’ (P. deltoides × P. euramericana) plants were propagated in 1/2 Murashige and Skoog (MS) (pH 5.8), and the ‘Shanxinyang’ (P. davidiana × P. bolleana Loucne) plants were cultivated in MS medium (pH 5.8) containing 0.3 mg/L IBA. To identify the transcript profiles of the PtFAD genes in different tissues, the mature leaves (ML), young leaves (YL), upper region of stems (US), lower region of stems (LS), and roots of 2-month-old poplars grown in the greenhouse were collected. In addition, to examine the expression patterns of the PtFAD genes under various stresses, the tissue culture seedlings of the poplar were treated at 0.2 M ABA, 0.2 M NaCl, 2 mM H2O2, and 10% PEG6000, and all the leaf samples were treated by various treatments and collected at 0, 1, 6, 12, 24, and 48 h. All the harvested samples were immediately frozen in liquid nitrogen and stored at −80 °C. Three experimental repetitions were conducted per sample.

4.7. Analysis of PtFAD Gene Expression Patterns

The specific primers of the PtFAD genes were designed to identify the expression patterns using a RT-qPCR (Supplemental Table S8). The total RNA was extracted using the hexadecyl trimethyl ammonium bromide (CTAB) method, as previously described, [71] for all poplar samples. The concentration and quality of the RNA were investigated using a NanoDrop One/OneC spectrophotometer (Thermo Scientific, Waltham, MA, USA). For the reverse-transcription, the PrimeScript™ RT Master Mix (TaKaRa, Kusatsu, Japan) was applied to reverse-transcribe 1 μg of RNA. The gene expression levels were quantified using the ABI 7500 Fast Real-Time PCR System (Applied Biosystem, Waltham, MA, USA) and the UltraSYBR Green I Mixture (CWBIO, Beijing, China). The RT-qPCR amplification reactions were performed in 20 mL volume with the following procedure: pre-denaturation at 95 °C for 5 min, followed by 40 cycles of 95 °C for 10 s, 60 °C for 30 s, and 72 °C for 30 s. The melting curve conditions were 60 °C for 60 s and 95 °C for 15 s. Three experimental repetitions were conducted per selected sample. The 2−ΔΔCt values were estimated to indicate the relative expression levels of the PtFAD genes. The TBtools was used to generate heatmaps according to the 2−ΔΔCt values.

4.8. Linkage Clustering Analysis of PtFAD and PtLACS Genes

Long-chain acyl-CoA synthetases (LACSs) catalyze FAs to form fatty acyl-CoA thioesters and play a vital role in FA metabolism [72]. Based on the expression patterns of the PtFADs and PtLACSs under the ABA, NaCl, H2O2, and PEG6000 treatments, the linkage clustering analysis was performed to explore the putative function of the PtFADs and PtLACSs on the lipid and FA metabolisms. The weighted gene cluster network analysis (WGCNA) package in the R programming language was applied to identify the gene cluster networks [73].

4.9. Analysis of SFA and UFA Contents

For the drought-stress of soil-grown poplars, the irrigation for poplars cultivated in the soil for three months was suspended, and the leaves were harvested after the drought treatment for two weeks. For the salt stress, the 3-month-old poplars were irrigated with a 200 mM NaCl solution for two weeks, and then the poplar leaves were collected. The experiments were performed under long-day (16-h light/8-h dark) conditions at 23 °C. Three experimental repetitions were applied, each time consisting of six lines. In addition, the collected poplar samples were freeze-dried for 48 h, and the samples were immersed in 2 mL of 10% acetylchloromethanol and 1 mL of n-hexane, and a total of 3 mL of reaction solvent was incubated for 2 h at 95 °C. Then, a total of 6 mL of 6% potassium carbonate solution was added to the above mixtures, and the n-hexane phase was obtained after centrifugation for 15 min at 4 °C, then, the n-hexane was removed from the n-hexane phase by vacuum concentration. The collected solvents containing the n-hexane and nonadecanoic acid were injected into a gas chromatography–mass spectrometry (GC–MS) (Thermo ISQ7000, Thermo Fisher, Waltham, MA, USA) or GC–flame ionization detector (Thermo Trace1300, Thermo Fisher, Waltham, MA, USA) with a chromatographic column DB-5 (60 m × 0.25 mm × 0.25 µm). In addition, GC–MS was applied to investigate the UFA contents. The GC–MS was performed as follows: (1) The column temperature was kept at 140 °C for 5 min, increased to 180 °C at the rate of 10 °C/min, and then rose to 210 °C at the rate of 4 °C/min. (2) The temperature was maintained at 310 °C for 30 min at the rate of 10 °C/min, and the capillary GC with a flame ionization detector was used to identify the individual components of the fatty acids in poplars. All data were presented as the mean ± standard deviation. A one-way analysis of variance (ANOVA), followed by the Tukey test, was used to determine the significant difference by the SPSS software (SPSS Inc., Chicago, IL, USA).

5. Conclusions

In this study, 23 PtFAD genes were identified from the P. trichocarpa genome. The phylogenetic tree indicated that the PtFADs could be divided into two subfamilies, including the soluble and membrane-bound FADs. The exon–intron compositions and conserved motifs of the PtFADs were similar within the same clade. A total of 17 syntenic gene pairs of the PtFAD members were found in the poplar genome, which might have a close association with the PtFAD duplication during evolution. Protein interaction analysis revealed that the PtFADs play a crucial role in the UFAs biosynthesis, and the GO analysis showed that the PtFADs are enriched in the molecular function and biological processes. The PtFAD3/7/8 mediated the regulation of the linolenic and linoleic acid levels and were associated with poplar’s resistance to drought and salt stresses. The study on the characterizations and profiles of the PtFAD genes provides insight into the potential functions of the PtFADs and sheds light on more candidate genes for genetic breeding.

Supplementary Materials

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

Author Contributions

J.Z.: funding acquisition; H.W.: writing—review and editing, methodology, software; A.M.: writing—review and editing, visualization; S.X., Y.Z., S.A.-D., M.G.Z., G.L., S.Z., C.Y., Y.C. and F.Z.: investigation, methodology, software, and formal analysis. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Jiangsu Forestry Science and Technology Innovation and Promotion Project (LYKJ [2021]11), Nantong University Scientific Research Start-up Project for Introducing Talents (135421609106), Science Foundation of Nantong (JC2020158), and the Priority Academic Program Development of Jiangsu Higher.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We gratefully thank our team for performing the experiment.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Olzmann, J.A.; Carvalho, P. Dynamics, and functions of lipid droplets. Nat. Rev. Mol. Cell Biol. 2019, 20, 137–155. [Google Scholar] [CrossRef] [PubMed]
  2. Spaggiari, M.; Ricci, A.; Calani, L.; Bresciani, L.; Neviani, E.; Dall’Asta, C.; Lazzi, C.; Galaverna, G. Solid state lactic acid fermentation: A strategy to improve wheat bran functionality. LWT 2020, 118, 108668. [Google Scholar] [CrossRef]
  3. Kim, H.; Rodriguez-Navas, C.; Kollipara, R.K.; Kapur, P.; Pedrosa, I.; Brugarolas, J.; Kittler, R.; Ye, J. Unsaturated fatty acids stimulate tumor growth through stabilization of β-catenin. Cell Rep. 2015, 13, 495–503. [Google Scholar] [CrossRef] [PubMed]
  4. D’Angeli, S.; Altamura, M.M. Unsaturated lipids change in olive tree drupe and seed during fruit development and in response to cold-stress and acclimation. Int. J. Mol. Sci. 2016, 17, 1889. [Google Scholar] [CrossRef]
  5. Gao, J.; Wallis, J.G.; Browse, J. Mutations in the prokaryotic pathway rescue the fatty acid biosynthesis1 mutant in the cold. Plant Physiol. 2015, 169, 442–452. [Google Scholar] [CrossRef]
  6. Li, N.; Xu, C.; Li-Beisson, Y.; Philippar, K. Fatty acid and lipid transport in plant cells. Trends Plant Sci. 2016, 21, 145–158. [Google Scholar] [CrossRef]
  7. Suri, K.; Singh, B.; Kaur, A.; Yadav, M.P.; Singh, N. Influence of microwave roasting on chemical composition, oxidative stability and fatty acid composition of flaxseed (Linum usitatissimum L.) oil. Food Chem. 2020, 326, 126974. [Google Scholar] [CrossRef]
  8. Xu, L.; Zeng, W.; Li, J.; Liu, H.; Yan, G.; Si, P.; Yang, C.; Shi, Y.; He, Q.; Zhou, W. Characteristics of membrane-bound fatty acid desaturase (FAD) genes in Brassica napus L. and their expressions under different cadmium and salinity stresses. Environ. Exp. Bot. 2019, 162, 144–156. [Google Scholar] [CrossRef]
  9. Zhang, Y.; Maximova, S.N.; Guiltinan, M.J. Characterization of a stearoyl-acyl carrier protein desaturase gene family from chocolate tree, Theobroma cacao L. Front. Plant Sci. 2015, 6, 239. [Google Scholar] [CrossRef]
  10. He, M.; Qin, C.X.; Wang, X.; Ding, N.Z. Plant unsaturated fatty acids: Biosynthesis and regulation. Front. Plant Sci. 2020, 11, 390. [Google Scholar] [CrossRef] [Green Version]
  11. Yao, L.; Shen, H.; Wang, N.; Tatlay, J.; Li, L.; Tan, T.W.; Lee, Y.K. Elevated acetyl-CoA by amino acid recycling fuels microalgal neutral lipid accumulation in exponential growth phase for biofuel production. Plant Biotechnol. J. 2017, 15, 497–509. [Google Scholar] [CrossRef] [PubMed]
  12. Browse, J.; Warwick, N.; Somerville, C.R.; Slack, C.R. Fluxes through the prokaryotic and eukaryotic pathways of lipid synthesis in the ‘16: 3′plant Arabidopsis thaliana. Biochem. J. 1986, 235, 25–31. [Google Scholar] [CrossRef] [PubMed]
  13. Zhao, X.; Wei, J.; He, L.; Zhang, Y.; Zhao, Y.; Xu, X.; Wei, Y.; Ge, S.; Ding, D.; Liu, M.; et al. Identification of fatty acid desaturases in maize and their differential responses to low and high temperature. Genes 2019, 10, 445. [Google Scholar] [CrossRef] [PubMed]
  14. Yadav Browse, J.; McCourt, P.; Somerville, C. A mutant of Arabidopsis deficient in C18:3 and C16:3 leaf lipids. Plant Physiol. 1986, 81, 859–864. [Google Scholar] [CrossRef] [PubMed]
  15. Yadav, N.S.; Wierzbicki, A.; Aegerter, M.; Caster, C.S.; Pérez-Grau, L.; Kinney, A.J.; Hitz, W.D.; Booth, J.R., Jr.; Schweiger, B.; Stecca, K.L.; et al. Cloning of higher plant [omega]-3 fatty acid desaturases. Plant Physiol. 1993, 103, 467–476. [Google Scholar] [CrossRef]
  16. McConn, M.; Hugly, S.; Browse, J.; Somerville, C. A mutation at the fad8 locus of Arabidopsis identifies a second chloroplast [omega]-3 desaturase. Plant Physiol. 1994, 106, 1609–1614. [Google Scholar] [CrossRef]
  17. Gibson, S.; Arondel, V.; Iba, K.; Somerville, C. Cloning of a temperature-regulated gene encoding a chloroplast [omega]-3 desaturase from Arabidopsis thaliana. Plant Physiol. 1994, 106, 1615–1621. [Google Scholar] [CrossRef]
  18. Shah, S.; Xin, Z.; Browse, J. Overexpression of the FAD3 desaturase gene in a mutant of Arabidopsis. Plant Physiol. 1997, 114, 1533–1539. [Google Scholar] [CrossRef]
  19. Okuley, J.; Lightner, J.; Feldmann, K.; Yadav, N.; Lark, E.; Browse, J. Arabidopsis FAD2 gene encodes the enzyme that is essential for polyunsaturated lipid synthesis. Plant Cell 1994, 6, 147–158. [Google Scholar]
  20. Dyer, J.M.; Mullen, R.T. Immunocytological localization of two plant fatty acid desaturases in the endoplasmic reticulum. FEBS Lett. 2001, 494, 44–47. [Google Scholar] [CrossRef]
  21. Wu, Q.S.; He, J.D.; Srivastava, A.K.; Zou, Y.N.; Kuča, K. Mycorrhizas enhance drought tolerance of citrus by altering root fatty acid compositions and their saturation levels. Tree Physiol. 2019, 39, 1149–1158. [Google Scholar] [CrossRef] [PubMed]
  22. Song, C.; Yang, Y.; Yang, T.; Ba, L.; Zhang, H.; Han, Y.; Lu, W. MaMYB4 recruits histone deacetylase MaHDA2 and modulates the expression of ω-3 fatty acid desaturase genes during cold stress response in banana fruit. Plant Cell Physiol. 2019, 60, 2410–2422. [Google Scholar] [CrossRef] [PubMed]
  23. Xue, M.; Guo, T.; Ren, M.; Wang, Z.; Tang, K.; Zhang, W.; Wang, M. Constitutive expression of chloroplast glycerol-3-phosphate acyltransferase from Ammopiptanthus mongolicus enhances unsaturation of chloroplast lipids and tolerance to chilling, freezing and oxidative stress in transgenic Arabidopsis. Plant Physiol. Biochem. 2019, 143, 375–387. [Google Scholar] [CrossRef] [PubMed]
  24. Shi, J.; Cao, Y.; Fan, X.; Li, M.; Wang, Y.; Ming, F. A rice microsomal delta-12 fatty acid desaturase can enhance resistance to cold stress in yeast and Oryza sativa. Mol. Breed. 2012, 29, 743–757. [Google Scholar] [CrossRef]
  25. Zhang, M.; Barg, R.; Yin, M.; Gueta-Dahan, Y.; Leikin-Frenkel, A.; Salts, Y.; Shabtai, S.; Ben-Hayyim, G. Modulated fatty acid desaturation via overexpression of two distinct ω-3 desaturases differentially alters tolerance to various abiotic stresses in transgenic tobacco cells and plants. Plant J. 2005, 44, 361–371. [Google Scholar] [CrossRef]
  26. Torres-Franklin, M.L.; Repellin, A.; Huynh, V.B.; d’Arcy-Lameta, A.; Zuily-Fodil, Y.; Pham-Thi, A.T. Omega-3 fatty acid desaturase (FAD3, FAD7, FAD8) gene expression and linolenic acid content in cowpea leaves submitted to drought and after rehydration. Environ. Exp. Bot. 2009, 65, 162–169. [Google Scholar] [CrossRef]
  27. McConn, M.; Browse, J. The critical requirement for linolenic acid is pollen development, not photosynthesis, in an Arabidopsis mutant. Plant Cell 1996, 8, 403–416. [Google Scholar] [CrossRef]
  28. Matsuda, O.; Sakamoto, H.; Hashimoto, T.; Iba, K. A temperature-sensitive mechanism that regulates post-translational stability of a plastidial ω-3 fatty acid desaturase (FAD8) in Arabidopsis leaf tissues. J. Biol. Chem. 2005, 280, 3597–3604. [Google Scholar] [CrossRef]
  29. Berberich, T.; Harada, M.; Sugawara, K.; Kodama, H.; Iba, K.; Kusano, T. Two maize genes encoding ω-3 fatty acid desaturase and their differential expression to temperature. Plant Mol. Biol. 1998, 36, 297–306. [Google Scholar] [CrossRef]
  30. Roman, A.; Andreu, V.; Hernandez, M.L.; Lagunas, B.; Picorel, R.; Martínez-Rivas, J.M.; Alfonso, M. Contribution of the different omega-3 fatty acid desaturase genes to the cold response in soybean. J. Exp. Bot. 2012, 63, 4973–4982. [Google Scholar] [CrossRef]
  31. Wang, X.; Yu, C.; Liu, Y.; Yang, L.; Li, Y.; Yao, W.; Cai, Y.; Yan, X.; Li, S.; Cai, Y.; et al. GmFAD3A, a ω-3 fatty acid desaturase gene, enhances cold tolerance and seed germination rate under low temperature in rice. Int. J. Mol. Sci. 2019, 20, 3796. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Wang, H.S.; Yu, C.; Tang, X.F.; Zhu, Z.J.; Ma, N.N.; Meng, Q.W. A tomato endoplasmic reticulum (ER)-type omega-3 fatty acid desaturase (LeFAD3) functions in early seedling tolerance to salinity stress. Plant Cell Rep. 2014, 33, 131–142. [Google Scholar] [CrossRef] [PubMed]
  33. Zhang, J.T.; Zhu, J.Q.; Zhu, Q.; Liu, H.; Gao, X.S.; Zhang, H.X. Fatty acid desaturase-6 (Fad6) is required for salt tolerance in Arabidopsis thaliana. Biochem. Biophys. Res. Commun. 2009, 390, 469–474. [Google Scholar] [CrossRef]
  34. Zhang, J.; Liu, H.; Sun, J.; Li, B.; Zhu, Q.; Chen, S.; Zhang, H. Arabidopsis fatty acid desaturase FAD2 is required for salt tolerance during seed germination and early seedling growth. PLoS ONE 2012, 7, e30355. [Google Scholar] [CrossRef]
  35. Chen, M.; Thelen, J.J. Acyl-Lipid Desaturase2 is required for chilling and freezing tolerance in Arabidopsis. Plant Cell 2013, 25, 1430–1444. [Google Scholar] [CrossRef]
  36. Chen, M.; Markham, J.E.; Cahoon, E.B. Sphingolipid Δ8 unsaturation is important for glucosylceramide biosynthesis and low-temperature performance in Arabidopsis. Plant J. 2012, 69, 769–781. [Google Scholar] [CrossRef] [PubMed]
  37. Dong, C.J.; Cao, N.; Zhang, Z.G.; Shang, Q.M. Characterization of the fatty acid desaturase genes in cucumber: Structure, phylogeny, and expression patterns. PLoS ONE 2016, 11, e0149917. [Google Scholar] [CrossRef] [PubMed]
  38. Liu, W.; Li, W.; He, Q.; Daud, M.K.; Chen, J.; Zhu, S. Characterization of 19 genes encoding membrane-bound fatty acid desaturases and their expression profiles in Gossypium raimondii under low temperature. PLoS ONE 2015, 10, e0123281. [Google Scholar] [CrossRef]
  39. Xue, Y.; Chen, B.; Wang, R.; Win, A.N.; Li, J.; Chai, Y. Genome-wide survey and characterization of fatty acid desaturase gene family in Brassica napus and its parental species. Appl. Biochem. Biotechnol. 2018, 184, 582–598. [Google Scholar] [CrossRef]
  40. Zhang, Z.; Wei, X.; Liu, W.; Min, X.; Jin, X.; Ndayambaza, B.; Wang, Y. Genome-wide identification and expression analysis of the fatty acid desaturase genes in Medicago truncatula. Biochem. Biophys. Res. Commun. 2018, 499, 361–367. [Google Scholar] [CrossRef]
  41. Chi, X.; Yang, Q.; Lu, Y.; Wang, J.; Zhang, Q.; Pan, L.; Chen, M.; He, Y.; Yu, S. Genome-wide analysis of fatty acid desaturases in soybean (Glycine max). Plant Mol. Biol. Rep. 2011, 29, 769–783. [Google Scholar] [CrossRef] [Green Version]
  42. Hajiahmadi, Z.; Abedi, A.; Wei, H.; Sun, W.; Ruan, H.; Zhuge, Q.; Movahedi, A. Identification, evolution, expression, and docking studies of fatty acid desaturase genes in wheat (Triticum aestivum L.). BMC Genom. 2020, 21, 778. [Google Scholar] [CrossRef] [PubMed]
  43. Zavaliev, R.; Mohan, R.; Chen, T.; Dong, X. Formation of NPR1 condensates promotes cell survival during the plant immune response. Cell 2020, 182, 1093–1108. [Google Scholar] [CrossRef]
  44. Liu, M.; Long, H.; Li, W.; Shi, M.; Cao, H.; Zhang, L.; Tan, X. Boosting C16 fatty acid biosynthesis of Escherichia coli, yeast and tobacco by Tung tree (Vernicia fordii Hemsl.) beta-hydroxyacyl-acyl carrier protein dehydratase gene. Ind. Crop. Prod. 2019, 127, 46–54. [Google Scholar] [CrossRef]
  45. Fu, A.; He, Z.; Cho, H.S.; Lima, A.; Buchanan, B.B.; Luan, S. A chloroplast cyclophilin functions in the assembly and maintenance of photosystem II in Arabidopsis thaliana. Proc. Natl. Acad. Sci. USA 2007, 104, 15947–15952. [Google Scholar] [CrossRef] [PubMed]
  46. Serrazina, S.; Dias, F.V.; Malhó, R. Characterization of FAB 1 phosphatidylinositol kinases in Arabidopsis pollen tube growth and fertilization. New Phytol. 2014, 203, 784–793. [Google Scholar] [CrossRef] [PubMed]
  47. Chen, H.; Sun, J.; Li, S.; Cui, Q.; Zhang, H.; Xin, F.; Wang, H.; Lin, T.; Gao, D.; Wang, S.; et al. An ACC oxidase gene essential for cucumber carpel development. Mol. Plant 2016, 9, 1315–1327. [Google Scholar] [CrossRef]
  48. Peng, Y.J.; Zhang, H.; Feng, M.G.; Ying, S.H. SterylAcetyl hydrolase 1 (BbSay1) links lipid homeostasis to conidiogenesis and virulence in the entomopathogenic fungus Beauveria bassiana. J. Fungi 2022, 8, 292. [Google Scholar] [CrossRef]
  49. Lee, S.B.; Go, Y.S.; Bae, H.J.; Park, J.H.; Cho, S.H.; Cho, H.J.; Lee, D.S.; Park, O.K.; Hwang, I.; Suh, M.C. Disruption of glycosylphosphatidylinositol-anchored lipid transfer protein gene altered cuticular lipid composition, increased plastoglobules, and enhanced susceptibility to infection by the fungal pathogen Alternaria brassicicola. Plant Physiol. 2009, 150, 42–54. [Google Scholar] [CrossRef]
  50. Luttgeharm, K.D.; Kimberlin, A.N.; Cahoon, E.B. Plant sphingolipid metabolism and function. Subcell. Biochem. 2016, 86, 249–286. [Google Scholar]
  51. Feng, J.; Dong, Y.; Liu, W.; He, Q.; Daud, M.K.; Chen, J.; Zhu, S. Genome-wide identification of membrane-bound fatty acid desaturase genes in Gossypium hirsutum and their expressions during abiotic stress. Sci. Rep. 2017, 7, 45711. [Google Scholar] [CrossRef] [PubMed]
  52. Xue, Y.; Chai, C.; Chen, B.; Shi, X.; Wang, B.; Mei, F.; Jiang, M.; Liao, X.; Yang, X.; Yuan, C.; et al. Whole-genome mining and in silico analysis of FAD gene family in Brassica juncea. J. Plant Biochem. Biot. 2020, 29, 149–154. [Google Scholar] [CrossRef]
  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] [PubMed]
  54. Liu, K.; Zhao, S.; Wang, S.; Wang, H.; Zhang, Z. Identification and analysis of the FAD gene family in walnuts (Juglans regia L.) based on transcriptome data. BMC Genom. 2020, 21, 299. [Google Scholar] [CrossRef] [PubMed]
  55. Cheng, C.; Liu, F.; Sun, X.; Wang, B.; Liu, J.; Ni, X.; Hu, C.; Deng, G.; Tong, Z.; Zhang, Y.; et al. Genome-wide identification of FAD gene family and their contributions to the temperature stresses and mutualistic and parasitic fungi colonization responses in banana. Int. J. Biol. Macromol. 2022, 204, 61–676. [Google Scholar] [CrossRef]
  56. Cao, J.; Shi, F. Evolution of the RALF gene family in plants: Gene duplication and selection patterns. Evol. Bioinform. 2012, 8, 271–292. [Google Scholar] [CrossRef]
  57. Chen, L.; Wang, L.; Wang, H.; Sun, R.; You, L.; Zheng, Y.; Yuan, Y.; Li, D. Identification and characterization of a plastidial ω-3 fatty acid desaturase EgFAD8 from oil palm (Elaeis guineensis Jacq.) and its promoter response to light and low temperature. PLoS ONE 2018, 13, e0196693. [Google Scholar] [CrossRef]
  58. Kim, M.J.; Kim, H.; Shin, J.S.; Chung, C.H.; Ohlrogge, J.B.; Suh, M.C. Seed-specific expression of sesame microsomal oleic acid desaturase is controlled by combinatorial properties between negative cis-regulatory elements in the SeFAD2 promoter and enhancers in the 5′-UTR intron. Mol. Genet. Genom. 2006, 276, 351–368. [Google Scholar] [CrossRef]
  59. Xiao, G.; Zhang, Z.Q.; Yin, C.F.; Liu, R.Y.; Wu, X.M.; Tan, T.L.; Chen, S.Y.; Lu, C.M.; Guan, C.Y. Characterization of the promoter and 5′-UTR intron of oleic acid desaturase (FAD2) gene in Brassica napus. Gene 2014, 545, 45–55. [Google Scholar] [CrossRef]
  60. Kim, H.; Go, Y.S.; Kim, A.Y.; Lee, S.; Kim, K.N.; Lee, G.J.; Kim, G.J.; Suh, M.C. Isolation and functional analysis of three microsomal delta-12 fatty acid desaturase genes from Camelina sativa (L.) cv. CAME. J. Plant Biochem. 2014, 41, 146–158. [Google Scholar] [CrossRef]
  61. Rajwade, A.V.; Kadoo, N.Y.; Borikar, S.P.; Harsulkar, A.M.; Ghorpade, P.B.; Gupta, V.S. Differential transcriptional activity of SAD, FAD2 and FAD3 desaturase genes in developing seeds of linseed contributes to varietal variation in α-linolenic acid content. Phytochemistry 2014, 98, 41–53. [Google Scholar] [CrossRef] [PubMed]
  62. Dar, A.A.; Choudhury, A.R.; Kancharla, P.K.; Arumugam, N. The FAD2 gene in plants: Occurrence, regulation, and role. Front. Plant Sci. 2017, 8, 789. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. 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]
  64. Yaeno, T.; Matsuda, O.; Iba, K. Role of chloroplast trienoic fatty acids in plant disease defense responses. Plant J. 2004, 40, 931–941. [Google Scholar] [CrossRef] [PubMed]
  65. Yang, J.; Xing, G.; Niu, L.; He, H.; Guo, D.; Du, Q.; Qian, X.; Yao, Y.; Li, H.; Zhong, X.; et al. Improved oil quality in transgenic soybean seeds by RNAi-mediated knockdown of GmFAD2-1B. Transg. Res. 2018, 27, 155–166. [Google Scholar] [CrossRef] [PubMed]
  66. Hwangbo, K.; Ahn, J.W.; Lim, J.M.; Park, Y.I.; Liu, J.R.; Jeong, W.J. Overexpression of stearoyl-ACP desaturase enhances accumulations of oleic acid in the green alga Chlamydomonas reinhardtii. Plant Biotechnol. Rep. 2014, 8, 135–142. [Google Scholar] [CrossRef]
  67. Moretti, S.; Francini, A.; Hernández, M.L.; Martínez-Rivas, J.M.; Sebastiani, L. Effect of saline irrigation on physiological traits, fatty acid composition and desaturase genes expression in olive fruit mesocarp. Plant Physiol. Bioch. 2019, 141, 423–430. [Google Scholar] [CrossRef]
  68. Dyer, J.M.; Stymne, S.; Green, A.G.; Carlsson, A.S. Highvalue oils from plants. Plant J. 2008, 54, 640–655. [Google Scholar] [CrossRef]
  69. Chen, C.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.; Xia, R. TBtools: An integrative toolkit developed for interactive analyses of big biological data. Mol. Plant 2020, 13, 1194–1202. [Google Scholar] [CrossRef]
  70. Kong, X.; Lv, W.; Jiang, S.; Zhang, D.; Cai, G.; Pan, J.; Li, D. Genome-wide identification and expression analysis of calcium-dependent protein kinase in maize. BMC Genom. 2013, 14, 433. [Google Scholar] [CrossRef]
  71. Meng, D.; He, M.; Bai, Y.; Xu, H.; Dandekar, A.M.; Fei, Z.; Cheng, L. Decreased sorbitol synthesis leads to abnormal stamen development and reduced pollen tube growth via an MYB transcription factor, MdMYB39L, in apple (Malus domestica). New Phytol. 2018, 217, 641–656. [Google Scholar] [CrossRef] [PubMed]
  72. Zhao, H.; Kosma, D.; Lu, S. Functional role of long chain acyl-CoA synthetases (LACS) in plant development and stress responses. Front. Plant Sci. 2021, 12, 276. [Google Scholar]
  73. Langfelder, P.; Horvath, S. WGCNA: An R package for weighted correlation network analysis. BMC Bioinform. 2008, 9, 559. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Phylogenetic relationships of the FAD proteins from Arabidopsis thaliana (At), Oryza sativa (Os), Populus trichocarpa (Pt), and Triticum aestivum (Ta). The PtFAD, AtFAD, and OsFAD proteins were aligned using the ClustalW, and the phylogenetic tree was constructed using the MEGA7.0. The FAD members were divided into seven clades: FAD3/7/8, FAD2, FAD6, DSD, FAD5, SLD, and FAB2, and different clades were indicated by different colors.
Figure 1. Phylogenetic relationships of the FAD proteins from Arabidopsis thaliana (At), Oryza sativa (Os), Populus trichocarpa (Pt), and Triticum aestivum (Ta). The PtFAD, AtFAD, and OsFAD proteins were aligned using the ClustalW, and the phylogenetic tree was constructed using the MEGA7.0. The FAD members were divided into seven clades: FAD3/7/8, FAD2, FAD6, DSD, FAD5, SLD, and FAB2, and different clades were indicated by different colors.
Ijms 23 11109 g001
Figure 2. Collinearity analysis of PtFAD gene family members. A rectangle represented chromosomes. The chromosomal positions of PtFADs were plotted using the TBtools. Individual PtFAD gene positions were labeled using the PtFAD name, and individual chromosomes were labeled with respective PtFADs in the circle. Each dark and red colored curve indicated the gene duplication event of all collinearity gene pairs and PtFADs across the chromosomes, respectively.
Figure 2. Collinearity analysis of PtFAD gene family members. A rectangle represented chromosomes. The chromosomal positions of PtFADs were plotted using the TBtools. Individual PtFAD gene positions were labeled using the PtFAD name, and individual chromosomes were labeled with respective PtFADs in the circle. Each dark and red colored curve indicated the gene duplication event of all collinearity gene pairs and PtFADs across the chromosomes, respectively.
Ijms 23 11109 g002
Figure 3. Syntenic relationship of FAD genes between P. trichocarpa and three other plant species, including A. thaliana (A), O. sativa (B), and Salix purpurea (C). The collinearity gene pairs and FAD gene pairs were represented in dark and red lines. The chromosomes of P. trichocarpa, A. thaliana, O. sativa, and S. purpurea were remarked by differently colored boxes.
Figure 3. Syntenic relationship of FAD genes between P. trichocarpa and three other plant species, including A. thaliana (A), O. sativa (B), and Salix purpurea (C). The collinearity gene pairs and FAD gene pairs were represented in dark and red lines. The chromosomes of P. trichocarpa, A. thaliana, O. sativa, and S. purpurea were remarked by differently colored boxes.
Ijms 23 11109 g003
Figure 4. Exon–intron distributions of PtFAD and AtFAD members. (A) The unrooted phylogenetic tree was shown, constructed using MEGA7.0 set to the neighbor-joining (NJ) method, and different clades were marked with different color backgrounds. (B) Exons, introns, and untranslated regions (UTR) were indicated by yellow rectangles, grey lines, and green rectangles, respectively. PtFAD and AtFAD members were clustered based on a phylogenetic tree. Base pair: bp.
Figure 4. Exon–intron distributions of PtFAD and AtFAD members. (A) The unrooted phylogenetic tree was shown, constructed using MEGA7.0 set to the neighbor-joining (NJ) method, and different clades were marked with different color backgrounds. (B) Exons, introns, and untranslated regions (UTR) were indicated by yellow rectangles, grey lines, and green rectangles, respectively. PtFAD and AtFAD members were clustered based on a phylogenetic tree. Base pair: bp.
Ijms 23 11109 g004
Figure 5. The conserved motifs of PtFAD and AtFAD members. The distribution of conserved motifs in the PtFAD and AtFAD proteins. The unrooted phylogenetic tree was shown, constructed using MEGA7.0 set to the neighbor-joining (NJ) method, and different clades were marked with different color backgrounds.
Figure 5. The conserved motifs of PtFAD and AtFAD members. The distribution of conserved motifs in the PtFAD and AtFAD proteins. The unrooted phylogenetic tree was shown, constructed using MEGA7.0 set to the neighbor-joining (NJ) method, and different clades were marked with different color backgrounds.
Ijms 23 11109 g005
Figure 6. Gene ontology (GO) analysis of PtFAD and interaction protein genes. GO analysis was performed using TBtools with a go-basic file, and PtFAD and interaction protein genes were enriched in molecular function, cellular component, and biological process.
Figure 6. Gene ontology (GO) analysis of PtFAD and interaction protein genes. GO analysis was performed using TBtools with a go-basic file, and PtFAD and interaction protein genes were enriched in molecular function, cellular component, and biological process.
Ijms 23 11109 g006
Figure 7. The expression pattern of PtFAD genes in various tissues and organs. The RT-qPCR was applied to investigate PtFAD expression patterns in P. trichocarpa (A), ‘Nanlin 895’ (Populus × euramericana cv.) (B), and ‘Shanxinyang’ (P. davidiana × P. bolleana) (C) tissues. The heatmap was generated using TBtools software. The color scale indicated the normalized data, where dark orange represented a high expression level, blue represented a low expression level, and orange represented a medium level. Mature leaves: ML; young leaves: YL; the upper region of stems: US; the lower region of stems: LS.
Figure 7. The expression pattern of PtFAD genes in various tissues and organs. The RT-qPCR was applied to investigate PtFAD expression patterns in P. trichocarpa (A), ‘Nanlin 895’ (Populus × euramericana cv.) (B), and ‘Shanxinyang’ (P. davidiana × P. bolleana) (C) tissues. The heatmap was generated using TBtools software. The color scale indicated the normalized data, where dark orange represented a high expression level, blue represented a low expression level, and orange represented a medium level. Mature leaves: ML; young leaves: YL; the upper region of stems: US; the lower region of stems: LS.
Ijms 23 11109 g007
Figure 8. The expression patterns of PtFAD genes under ABA (A), H2O2 (B), NaCl (C), and PEG6000 (D). The heatmap was generated using TBtools software. The color scale indicated the normalized data, where pink represented a high expression level, light pink represented a low expression level, and light green represented a medium level.
Figure 8. The expression patterns of PtFAD genes under ABA (A), H2O2 (B), NaCl (C), and PEG6000 (D). The heatmap was generated using TBtools software. The color scale indicated the normalized data, where pink represented a high expression level, light pink represented a low expression level, and light green represented a medium level.
Ijms 23 11109 g008
Figure 9. The cluster dendrogram of PtFAD and PtLACS expression patterns under ABA (A), H2O2 (B), NaCl (C), and PEG6000 (D). The data were normalized to PtActin (XM-006370951) expression level. Each PtFAD and PtLACS expression level was calculated based on the corresponding gene mRNA level in poplar without treatments.
Figure 9. The cluster dendrogram of PtFAD and PtLACS expression patterns under ABA (A), H2O2 (B), NaCl (C), and PEG6000 (D). The data were normalized to PtActin (XM-006370951) expression level. Each PtFAD and PtLACS expression level was calculated based on the corresponding gene mRNA level in poplar without treatments.
Ijms 23 11109 g009
Figure 10. The gas chromatography-mass spectrometry (GC-MS) analysis of the SFA and UFA contents of poplar leaves before and after drought or salt stress. (A) The contents of C14-C17 SFAs and UFAs in poplar leaves cultivated under normal conditions and drought or salt stress. (B) The contents of C18-C24 SFAs and UFAs in poplar leaves cultivated under normal conditions and drought or salt stress. Three biological replicates were performed per group. * p < 0.05, ** p < 0.01, *** p < 0.001 compared to control poplar. The kinds of SFAs and UFAs included myristate (C14:0), cis-9-tetradecenoate (C14:1), pentadecanoate (C15:0), cis-10-pentadecenoate (C15:1), palmitate (C16:0), palmitoleate (C16:1), margarate (C17:0), cis-10-heptadecenoate (C17:1), octadecanoate (C18:0), oleate (C18:1), linoleic acid (C18:2), linolenic acid (C18:3), icosanoate (C20:0), eicosenoate (cis-11) (C20:1), eicosadienoate (cis-11, 14) (C20:2), eicosatrienoate (cis-8,11,14) (C20:3), behenate (C22:0), erucate (C22:1), lignocerate (C24:0), and nervate (C24:1).
Figure 10. The gas chromatography-mass spectrometry (GC-MS) analysis of the SFA and UFA contents of poplar leaves before and after drought or salt stress. (A) The contents of C14-C17 SFAs and UFAs in poplar leaves cultivated under normal conditions and drought or salt stress. (B) The contents of C18-C24 SFAs and UFAs in poplar leaves cultivated under normal conditions and drought or salt stress. Three biological replicates were performed per group. * p < 0.05, ** p < 0.01, *** p < 0.001 compared to control poplar. The kinds of SFAs and UFAs included myristate (C14:0), cis-9-tetradecenoate (C14:1), pentadecanoate (C15:0), cis-10-pentadecenoate (C15:1), palmitate (C16:0), palmitoleate (C16:1), margarate (C17:0), cis-10-heptadecenoate (C17:1), octadecanoate (C18:0), oleate (C18:1), linoleic acid (C18:2), linolenic acid (C18:3), icosanoate (C20:0), eicosenoate (cis-11) (C20:1), eicosadienoate (cis-11, 14) (C20:2), eicosatrienoate (cis-8,11,14) (C20:3), behenate (C22:0), erucate (C22:1), lignocerate (C24:0), and nervate (C24:1).
Ijms 23 11109 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wei, H.; Movahedi, A.; Xu, S.; Zhang, Y.; Liu, G.; Aghaei-Dargiri, S.; Ghaderi Zefrehei, M.; Zhu, S.; Yu, C.; Chen, Y.; et al. Genome-Wide Characterization and Expression Analysis of Fatty acid Desaturase Gene Family in Poplar. Int. J. Mol. Sci. 2022, 23, 11109. https://doi.org/10.3390/ijms231911109

AMA Style

Wei H, Movahedi A, Xu S, Zhang Y, Liu G, Aghaei-Dargiri S, Ghaderi Zefrehei M, Zhu S, Yu C, Chen Y, et al. Genome-Wide Characterization and Expression Analysis of Fatty acid Desaturase Gene Family in Poplar. International Journal of Molecular Sciences. 2022; 23(19):11109. https://doi.org/10.3390/ijms231911109

Chicago/Turabian Style

Wei, Hui, Ali Movahedi, Songzhi Xu, Yanyan Zhang, Guoyuan Liu, Soheila Aghaei-Dargiri, Mostafa Ghaderi Zefrehei, Sheng Zhu, Chunmei Yu, Yanhong Chen, and et al. 2022. "Genome-Wide Characterization and Expression Analysis of Fatty acid Desaturase Gene Family in Poplar" International Journal of Molecular Sciences 23, no. 19: 11109. https://doi.org/10.3390/ijms231911109

APA Style

Wei, H., Movahedi, A., Xu, S., Zhang, Y., Liu, G., Aghaei-Dargiri, S., Ghaderi Zefrehei, M., Zhu, S., Yu, C., Chen, Y., Zhong, F., & Zhang, J. (2022). Genome-Wide Characterization and Expression Analysis of Fatty acid Desaturase Gene Family in Poplar. International Journal of Molecular Sciences, 23(19), 11109. https://doi.org/10.3390/ijms231911109

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