刘 聪，张 洋，夏德安，陈雪冰，魏志刚

(1. 东北林业大学，林木遗传育种国家重点实验室，黑龙江 哈尔滨 150000；2. 国家林业与草原局盐碱地研究中心，北京 100091)

植物作为固着生物，需要应对干旱、盐碱、低温、高温以及病虫害侵袭等多种逆境胁迫，而对胁迫信号的感知和转导是植物应对环境胁迫的前提和基础[1]。研究表明，植物体内胁迫信号的转导需要多种信号转导物质参与，如ABA、乙烯、H2O2和NO等激素和小分子化合物[2]以及信号转导蛋白G蛋白[3]、磷脂酶A[4]、磷脂酶C[5]、磷脂酶D（PLD）[6]等激酶，其中，PLD能够水解磷酸二脂键使细胞膜磷脂产生磷脂酸（PA）和可溶性头基，而PA是多种激素信号和胁迫信号转导过程的第二信使，在信号转导、膜迁移、细胞程序性死亡和细胞骨架重排等多种细胞活动中发挥着重要作用[7-10]。研究表明，PLD的典型特征是具有2个间隔250～400个氨基酸的HKD（HxKxxxxD）结构域，虽然2个HKD在蛋白初级结构上相距较远，但它们在三级结构上能形成互作，且这种互作也是其酶活性发挥的关键[11-15]。根据氨基酸序列特征，PLD可分为α、β、γ、δ、ε、ζ、φ 7种类型，其中，PLDα、PLDγ、PLDδ和PLDζ数量较多，PLDβ和PLDε数量较少，而PLDφ发现较晚且较少见[13,16]。根据蛋白的N端结构特征和其酶活性对Ca2+的依赖性，PLD也可分为C2-PLD、PX/PH-PLD和SPPLD 3个亚族，其中，C2-PLD蛋白的N端含有一个高度保守Ca2+依赖型磷脂结合C2结构域，需要结合Ca2+以维持酶活性[12]；PX/PH-PLD的N端特征是一个串连的PX（phox homology）和PH（pleckstrin homology）结构域，活性不依赖于Ca2+[13]；SP-PLD亚家族规模较小，于2007年在水稻（Oryza sativaL.）中被首次发现，其N端含有信号肽，而不是C2、PX或PH结构域[13]。其中，PLDα、PLDβ、PLDγ、PLDδ和PLDε属于C2-PLD亚家族，PLDζ属于PX/PH-PLD亚家族，PLDφ属于SP-PLD亚家族[11]。植物中，首个PLD于1994年从蓖麻（Ricinus CommunisL.）[17]中克隆出来后，随后水稻[18]、拟南芥（Arabidopsis thalianaL.）[19]、玉米（Zea maysL.）以 及 烟 草（Nicotiana tabacumL.）[20]等多个物种中相继被克隆出PLD。随着多个植物全基因组序列的发布[21-27]，人们发现PLD在植物中多以基因家族形式存在，如PLD家族成员拟南芥有12个[12]、水稻有17个[13]、白菜（Brassica rapaLour.）有18个[14]、亚洲棉（Gossypium arboreumL.）有19个[11]、葡萄（Vitis viniferaL.）有11个[28]、大豆（Glycine maxL.）有18个[15]等。

研究发现，PLD在植物响应逆境胁迫反应中有着重要作用，如在拟南芥中，AtPLDα1的沉默会严重减弱植株由ABA诱导的耐旱反应[29]；Atpldα1缺陷型株系中由创伤诱导的 JA合成和相关基因的表达显著降低[30]；Atpldδ缺陷型株系对H2O2胁迫的敏感性显著上升[31]；AtPLDδ在响应真菌（Blumeria graminisDC.）胁迫的防御信号转导过程中有着关键作用，Atpldδ突变体株系中响应微生物的早期应激响应基因的上调表达被显著延迟[32]；C2-PLD突变体株系在响应缺氧胁迫过程中关键信号转导物质PA和胞质钙的含量上升显著受阻[33]；在其他植物中，如在毛白杨（Populus tomentosaCarr.in Rev.Hort.）中过表达AtPLDα1能够显著提高其耐盐性与耐旱性[34]；过表达黄瓜（Cucumis sativusL.）CsPLDα的烟草抗旱性增强[35]；小麦（Triticum aestivumL.）TaPLDα能够响应多种胁迫下表达上调，过表达TaPLDα的拟南芥耐旱能力显著增强[36]。由此可见，PLD家族在植物应对逆境胁迫时具有重要生物学功能。

毛果杨（Populus trichocarpaTorr. & Gray）是首个全基因组测序的木本植物[37]，也是研究木本植物生长发育、材质材性、逆境胁迫响应和其他重要性状的模式植物。然而，目前为止，尚无毛果杨PLD家族鉴定及其在逆境胁迫响应特性的研究报道。本项研究在全基因组水平鉴定出毛果杨PLD家族全部成员的基础上，对其系统进化、同源关系、选择压力在其家族进化中的作用、基因结构、家族各基因的顺式作用元件和蛋白保守基序进行分析，从而全面阐述该基因家族的基本特征。此外，利用实时荧光定量PCR和已报道的毛果杨转录组数据，对其PLD家族基因的组织表达特性以及盐胁迫的响应特性进行分析。本研究为进一步阐明毛果杨PLD家族基因的生物学功能和抗逆基因的挖掘奠定了基础。

1 材料与方法

1.1 毛果杨PLDs家族成员鉴定及基本特征分析

利用报道的拟南芥PLD的氨基酸序列比对Phytozome网站（https://phytozome.jgi.doe.gov）中的毛果杨基因组数据库[12]，同时在毛果杨数据库检索已注释的PLD，将获得的序列结果汇总并剔除重复序列后作为候选基因。为了验证初始结果的可靠性，将候选基因的氨基酸序列上传至HMMER网站 (https://www.ebi.ac.uk/Tools/hmmer)以及NCBI保守结构域数据库（https://www.ncbi.nlm.nih.gov/cdd/）鉴定保守结构域，以含有2个间隔250～400氨基酸的PLD典型HKD结构域（HMM模型登录号为PF00614）为标准进行筛选，并根据保守结构域鉴定蛋白类型，最终鉴定出毛果杨PLD家族全部成员[11,14,38-40]。

从Phytozome网站中获得PtrPLD家族成员的染色体位置、基因序列和开放阅读框长度等信息，并根据基因编码蛋白类型和进化关系进行命名，染色体位置信息经Tbtools软件进行可视化处理[41]。通过ExPasy（https://web.expasy.org）在线预测其等电点与分子质量[42]，通过Plant-mPLoc(http://www.csbio.sjtu.edu.cn/)在线预测其亚细胞定位位置[43]。

1.2 系统进化分析

将鉴定出的毛果杨全部PtrPLDs以及拟南芥AtPLDs编码的氨基酸序列在MEGA X软件的ClustaW程序中进行多重序列比对[44]。采用邻近法（Neighbor-joining）构建系统进化树，其中，bootstrap设置1 000次重复，再经最大似然法（Maximum-likeihood）验证，得到系统发育进化树数据，经Evoview（https://www.evolgenius.info/evolview/）网站可视化处理[12,44-45]。

1.3 同源基因对的Ka和Ks分析

Blast相互比对PtrPLDs的CDS序列，超过300个bp且同源性超过80%为标准确定为同源基因对[46-47]，同源关系经Tbtools软件可视化处理[41]。利用Tbtools软件计算同源基因对之间的同义替换率（Ks）、非同义替换率（Ka）以及Ka/Ks比率，并以此分析PtrPLD家族基因进化过程中的选择压力[41, 46, 48]。

1.4 基因结构及蛋白保守型基序分析

从毛果杨数据库（https://genome.jgi.doe.gov/portal/pages/dynamicOrganismDownload.jsf?organis m=Ptrichocarpa）下载基因组数据后，通过TBtools软件提取各PtrPLD外显子、内含子长度及位置信息并进行可视化处理[41]。

将各PtrPLD氨基酸序列提交到Motif Elicitation网站的MEME程序（http://meme-suite.org/tools/meme）中分析保守基序并进行可视化处理[41,49]。

1.5 启动子区顺式作用元件分析

在 Phytozome网 站（https://phytozome.jgi.doe.gov）中将各PtrPLD起始密码子前2 000 bp的序列作为启动子区域汇总，再上传至Plantcare网站（http://bioinformatics.psb.ugent.be/webtools/plantcare/html）进行顺式作用元件在线分析[46]。

1.6 PtrPLD家族基因的组织表达及其对盐胁迫的响应特性分析

组织表达特异性：在Phytozome网站（https://phytozome.jgi.doe.gov）中下载各PtrPLD在各组织中的表达量数据，再采用qRT-PCR进行验证。野生型毛果杨来自中国科学院上海生命科学研究院，用组织培养的方法扩繁后，将1月大小的组培苗移植到土壤中，在25℃、长日照（光照16 h黑暗8 h）温室中培养3周。分别采集根、茎和叶组织，提取总RNA后反转录获得 cDNA用于qRT-PCR。每组处理重复3次，采用2-ΔΔCT法计算相对表达量并利用Tbtools可视化[41]。

盐胁迫下的响应特性：将1月大小的组培苗移植到土壤中，在25℃、长日照（光照16 h黑暗8 h）温室中培养3周时，将长势一致的幼苗随机分成7组，每组包含5棵毛果杨幼苗。用100 mmol·L-1NaCl处理3、6、12、24、36、72 h，同时用水处理作为对照组。分别采集上述处理组内各植株材料的根、茎和叶组织，提取总RNA后反转录获得cDNA进行qRT-PCR分析。每组处理重复3次，采用2-ΔΔCT法计算相对表达量并利用Tbtools可视化[41]。

RNA提取、反转录及qRT-PCR：利用植物总RNA提取试剂盒（MiniBEST，TaKaRa）提取总RNA，然后采用 PrimeScriptTMRT reagent Kit（Perfect Real Time，TaKaRa）试剂盒反转录 RNA获得 cDNA用于 qRT-PCR。根据荧光定量引物设计原则，利用NCBI Primer-Blast（www.ncbi.nlm.nih.gov/tools/primer）工具设计特异性的PtrPLD家族基因定量引物，以PtrACTIN为内参基因（表1），并通过PCR、电泳验证引物特异性[50]。在赛默飞ABI 7500实时荧光定量PCR仪上进行试验，体系如下:2×TransStart TOP/Tip Green qPCR Supermix 10 μL、上下游混合引物(10 μmol·L-1) 0.4 μL、cDNA 1.5 μL，Passive Reference Dye (50×) 0.4 μL，加ddH2O至20 μL。反应程序：94℃ 30 s; 94℃ 5 s，60℃ 15 s，72℃ 35 s，循环40次; 95℃ 15 s，60℃ 1 min，95℃ 30 s。

表1 毛果杨PtrPLD定量引物序列Table 1 PtrPLD quantitative primer sequence

2 结果与分析

2.1 PtrPLD家族成员及其编码蛋白基本特性

根据已报道的AtPLD家族各基因氨基酸序列比对毛果杨的基因组数据库得到18个候选基因，再经过HKD结构域且二者间隔250～400个氨基酸为标准筛选，剔除Potri.001G112100.1和Potri.T 180000.1（2个基因所编码的蛋白中HKD序列不符合HxKxxxxD特征），共得到16个PtrPLD（表2）。按编码PLD蛋白类型和进化关系命名为PtrPLDα1到PtrPLDξ3（表2）；同时，对PtrPLD家族成员编码蛋白特性进行分析，结果表明，PtrPLD家族基因编码蛋白的氨基酸残基数为 645～1 141、分子量为73.79～129.02 kDa、理论等电点为5.41～6.90、单个HKD结构域跨度在35个氨基酸左右、2个HKD结构域间隔均在300个氨基酸左右（表2）。亚细胞定位预测发现，PtrPLDα1和PtrPLDα2定位在细胞质中，PtrPLDα3、PtrPLDα4、PtrPLDα5和PtrPLDα6定位在内质网和液泡中，PtrPLDβ1定位在叶绿体中，PtrPLDβ2定位在叶绿体和细胞质中，PtrPLDδ1～5和PtrPLDξ1～3均定位在细胞质中（表2）。

表2 毛果杨PLD基因及所编码蛋白的特征Table 2 Characteristics of PLD genes and encoded proteins in Populus trichocarpa

2.2 PtrPLD家族的进化关系

为了解PtrPLD家族成员的进化关系，将16个PtrPLDs与12个AtPLDs编码的蛋白构建进化树。根据拟南芥PLD基因家族进化关系，将PtrPLD家族分为α、β、δ及ζ 4个类型，其中α类型有6个，β类型有2个，δ类型有5个，ζ类型有3个，未鉴定出γ、ε和φ类型的PtrPLD（图1）。结果表明：毛果杨基因组并不存在SP-PLD亚族，PX/PH-PLD亚家族规模较小仅有3个成员，C2-PLD亚家族共有13个成员。

图1 毛果杨与拟南芥PLD基因家族的系统发育进化分析。Fig.1 Phylogenetic analysis of PLD gene family in Populus trichocarpa and Arabidopsis thaliana.

2.3 PtrPLD家族基因扩张

为了分析PtrPLD家族扩张原因，利用Tbtools构建了家族基因染色体定位图以及同源关系图[41]。结果表明：毛果杨16个PtrPLDs不均匀地分布在1、2、3、5、6、7、8、10、13、14和18号染色体上，其中，3号染色体上有4个基因、2号和5号染色体上有2个，其余的染色体上均仅有1个PtrPLD家族基因（图2）。

毛果杨PtrPLD家族共线性分析以及序列Blastn分 析 结 果 表 明：PtrPLDα1和PtrPLDα2、PtrPLDα5和PtrPLDα3、PtrPLDβ1和PtrPLDβ2、PtrPLDδ1和PtrPLDδ4以及PtrPLDδ2、PtrPLDδ3和PtrPLDδ4、PtrPLDξ2和PtrPLDξ3有共线性关系且同源片段长度远大于300 bp同源性超过80%（表3），因此上述基因对是由基因复制事件的进化形成的旁系同源基因[46-47]。其中，3号染色体上出现由PtrPLDα1、PtrPLDα2和PtrPLDα3形成的基因簇（图2），PtrPLDα1和PtrPLDα2同源性极高（表3），因此，PtrPLDα1和PtrPLDα2是由基因串联复制事件演化而来。综上所述，多数PtrPLDs具有基因复制现象，表明基因复制是PtrPLDs家族扩张的主要原因。

表3 同源基因的Ka/Ks比值及同源性Table 3 Ka/Ks values and homologous status of homologous genes

图2 PtrPLDs基因染色体分布Fig.2 Chromosome distribution of PtrPLDs

为了阐明选择压力在PtrPLD家族进化中的作用，利用Tbtools分析了PtrPLD家族同源基因的Ks值、Ka值和Ka/Ks值。结果发现：7对同源基因Ka/Ks均远小于1（表3），表明PtrPLD家族在进化过程中经历了较强的纯化选择。

2.4 PtrPLD家族基因结构以及保守域特征

为了解PtrPLD家族各基因系统进化与其基因结构及其编码蛋白中保守基序之间的关系，对PtrPLD家族基因外显子与内含子的结构以及编码蛋白保守结构域进行了分析。结果表明：PtrPLD家族各基因的内含子和外显子分布方式在保持家族保守性的同时，在进化过程中也产生了明显的分化。除PtrPLDα5和PtrPLDδ3外，其他的家族基因最后一个外显子均紧靠3′端UTR（图3）。同时，多数处于同一进化分支上的基因结构较为相似，如PtrPLDβ1和PtrPLDβ2均仅含有4个外显子；值得注意的是，PtrPLDα1和PtrPLDα2均含有20个外显子，但内含子与外显子在基因中的分布松散，同时，上述2个基因序列长度虽然超过12 000 bp，且远大于家族内其他基因序列长度，然而二者编码氨基酸长度却是家族中最短的；此外，PtrPLDβs基因长度均在5 000 bp以下，但其氨基酸数量最多，长度最长（图3）。

图3 PtrPLD家族基因的外显子-内含子结构分析Fig.3 Exon-intron organization analysis of PtrPLD genes

PtrPLD家族基因编码蛋白具有20个基序，其中Motif 1、Motif 2、Motif 3 和Motif 4为家族共有基序。虽然不同基因编码蛋白所含基序数量与种类存在一定的差异，但相同分支的基因编码蛋白具有相似的基序组合，而不同亚族、类型间的基序有明显的区别（图4）。上述结果表明：PtrPLD家族基因编码蛋白结构在具有保守性的同时也产生了明显分化。

2.5 PtrPLD家族基因启动子元件特征

为了解PtrPLD家族基因可能的生物学功能和调控特性，在利用PlantCARE分析家族成员启动子的同时，结合PtrPLD家族基因系统发育树分析各基因启动子元件组成与基因的进化之间是否具有相关性。结果表明：PtrPLD家族16个基因启动子中共包含2类、12种、154个顺式作用元件（表4）。一大类是非生物胁迫响应元件，如厌氧诱导元件（ARE）、损伤反应元件（WUN-motif）、胁迫响应元件（STRE）、MYB干旱诱导性结合位点（MBS）、DREB/CBF转录因子识别位点（DRE）、抗病和胁迫诱导元件（TC-rich repeats）以及低温响应元件（LTR）；另一大类是植物激素响应元件，如生长素响应元件（TGA-element）、水杨酸诱导元件（TCA-element）、茉莉酸甲酯响应元件（CGTCA-motif 和TGACG-motif）、ABA应答元件（ABRE）。

表4 PtrPLD家族基因启动子区域的各顺式作用元件分析Table 4 Analysis of cis acting elements in PtrPLD genes promoter region

PtrPLD家族各基因启动子含有元件的种类和数量存在差异，其中，PtrPLDδ启动子中含有的顺式作用元件数量最多，如PtrPLDδ4启动子中具有的元件数量最多（20个），而PtrPLDα6和PtrPLDξ3启动子中仅含有2个元件。上述结果表明，基因在进化的同时，其启动子中元件也会产生相应的分化。

此外，在所有顺式作用元件中ABA应答元件ABRE数量最多，在11个PtrPLD启动子中有48个ABRE，其中，PtrPLDδ1、PtrPLDδ2以及PtrPLDδ4含有的ABRE元件达8个；厌氧诱导元件ARE含量次之，共有19个，分布在PtrPLDα6、PtrPLDδ4、PtrPLDδ2等11个基因启动子中；而DREB/CBF转录因子识别位点DRE数量最少，仅4个，分布在PtrPLDδ1、PtrPLDδ3、PtrPLDδ4和PtrPLDξ3基因启动子中。上述结果表明：PtrPLD家族不同的亚族基因响应植物激素和非生物胁迫的能力存在差异。

2.6 PtrPLD家族基因的组织表达与盐胁迫响应的特性

为了解PtrPLD家族基因的表达模式，对phytozome网站中PtrPLD家族各基因在不同组织的表达量进行分析。结果表明：6个PtrPLDα型蛋白基因中PtrPLDα1、PtrPLDα2和PtrPLDα6在茎中的表达水平较高，PtrPLDα4和PtrPLDα5在根部的表达水平较高，PtrPLDα3在嫩叶和根部的表达水平较高，并未有明显的表达部位偏好；2个PtrPLDβ型蛋白基因在根部的表达水平均较高；5个PtrPLDδ型蛋白基因中只有PtrPLDδ2在根、嫩叶和成熟叶的表达水平较高，其他4个均在根的表达量较高；3个PtrPLDξ型蛋白基因均在茎节的表达水平较高，根部次之，叶中表达量最低（图5A）。同时，利用qRT-PCR对PtrPLD家族各基因在根、茎和叶中的表达量做进一步分析，结果表明：16个PtrPLD家族基因中有9个在根部的表达水平较高，3个在茎部表达水平较高，4个在叶部表达水平较高。qRT-PCR与phytozome网站中的基本相符,大部分的基因在根部的表达水平较高，但不同的是qRT-PCR结果中的PtrPLDα3和PtrPLDα5在叶中表达水平较高；3个PtrPLDξ型蛋白基因中PtrPLDξ1在 叶 部 表 达 水 平 较 高，PtrPLDξ2和PtrPLDξ3均在根部的表达量较高（图5B）；PtrPLDδ2在叶中的表达水平高于根部的。

图5 PtrPLD家族基因组织表达特异性Fig.5 Tissue expression specificity of PtrPLD family genes

为了阐明PtrPLD家族基因的盐胁迫响应特性，利用qRT-PCR对100 mmol·L-1NaCl胁迫处理3、6、12、24、48、72 h的毛果杨幼苗的根、茎、叶中PtrPLD家族各基因的表达进行了研究。结果（图6）表明：NaCl胁迫对16个PtrPLD基因在各组织中表达均产生了影响；在根、茎、叶3种组织中全部PtrPLD家族基因在盐胁迫后12 h内表现出显著的表达上调趋势，且大部分是在6 h或12 h时达到表达高峰，之后快速下调，在48 h到72 h之间又略有上升。较为独特的是，PtrPLDδ4在茎部的表达高峰处于胁迫处理72 h时；PtrPLDα1在根、茎、叶中的48～72 h之间的上调表达趋势相较于其他PtrPLD的均较明显。从响应速度来看，叶和茎部中大多数PtrPLD家族基因响应盐胁迫的快速上升期发生在3 h或6 h，而在根部中的快速上升期发生在6 h或12 h。

图6 不同组织中PtrPLD基因家族在盐胁迫下的表达特性Fig.6 Analysis of salt stress expression of PtrPLD family genes in different tissues

此外，各部位中均有部分PtrPLD家族基因在12 h内的表达高峰前出现下调表达的情况，而不是在胁迫处理后一直上调表达到峰值，如根部中的全部PtrPLDα、PtrPLDδ1、PtrPLDδ3和PtrPLDδ4，茎部的PtrPLDδ2、PtrPLDδ4和PtrPLDδ5，叶部的PtrPLDα2、PtrPLDδ1和PtrPLDδ4。

3 讨论

研究表明，PLD以多基因家族的形式存在于拟南芥[12]、水稻[13]、白菜[14]、亚洲棉[11]等植物中，并在胁迫响应中起到重要作用[10,17,29,34]。然而，目前为止，尚无模式植物毛果杨PLD家族基因基本特征及其在逆境胁迫响应特性的研究报导。本研究通过对毛果杨全基因组分析，共鉴定出16个PtrPLD，其编码蛋白包括3个PX/PH-PLD和13个C2-PLD亚家族成员（图1），这与其它植物中C2-PLD成员数量相对较少的情况相同[12]。PtrPLD家族中共有7对具有旁系同源关系，其中PtrPLDα1和PtrPLDα2是由基因串联复制事件演化而来，其余的6对是由基因组复制事件演化而来，这表明基因复制是PtrPLD基因家族扩张的主要原因（图2和表3）。PtrPLD家族中的7对同源基因的Ka/Ks比值远低于1（表3），表明它们在进化过程中经过了较强的纯化选择，有害的非同义替换在进化过程消失，极少数的无害或有益替换得以保留[48]。上述现象间接证明了PtrPLD家族基因在毛果杨的生命活动中发挥着重要作用。

PtrPLD家族基因编码蛋白含有多个基序，表明该基因家族生物学功能的多样性。同时，所有的PtrPLD家族基因编码蛋白均含有Motif 1～4基序（图4），表明这4个基序是该基因家族的特征基序，且是毛果杨生命活动所必须的。此外，序列比对表明，Motif 1中的HKD基序中的甲硫氨酸高度保守，这种现象也存在于其他植物的PLD蛋白中[11,51-52]，因此，推测Motif 1中的甲硫氨酸是HKD基序功能的关键位点。

PtrPLD家族成员的表达特性分析表明，其中9个家族成员在根部表达水平较高，3个在茎部表达水平较高，4个在叶部表达水平较高（图5）。前期研究认为，PLD基因多涉及环境胁迫响应[10,17,29,34]，而根部是植物感知渗透胁迫的首要部位。据此推测，在毛果杨根部表达水平较高的9个PtrPLD基因可能在根部渗透胁迫信号转导过程中发挥重要作用，而其余7个在茎、叶中表达水平较高的可能在茎部或叶部的胁迫信号转导过程中发挥功能。此外，phytozome网站中数据与qRT-PCR结果中的组织表达特性不同的原因可能是所用材料的树龄不同，同时它们在不同发育阶段中的表达模式存在差异所致。

上游转录因子结合下游基因启动子中顺式作用元件，从而激活逆境响应基因的表达是植物响应环境胁迫的重要环节[53]。因此，对基因启动子区域顺式元件的分析，是鉴定分析基因生物学功能的一项重要内容。本研究发现，PtrPLD家族基因启动子中含有多种与植物激素和逆境胁迫响应相关的两大类响应元件（表4），表明该基因家族广泛地参与毛果杨多种逆境胁迫的响应。PtrPLD家族基因在盐胁迫下的响应特性也证实了上述观点。

4 结论

本研究通过氨基酸序列比对及保守结构域检验，在毛果杨中鉴定了16个PtrPLDs；通过系统发育进化分析将其分为2个亚族，4种类型；利用同源性分析及Ka、Ks分析PtrPLD基因家族说明基因复制是其家族扩张的主要动力，且在进化过程中受到较强地纯化选择；系统发育进化树与基因结构和保守结构域分析表明，同一亚族、同一类型的PtrPLD具有相似的基因结构和蛋白保守结构域；启动子中顺式作用元件分析表明，16个PtrPLDs均能够响应多种激素和胁迫信号；组织表达特异性以及盐胁迫下的表达特性分析表明，同一进化分支上基因仍具有相同表达模式，而不同进化分支的基因存在功能分化。本研究对于PtrPLD家族基因生物学功能的鉴定和抗逆基因资源挖掘具有积极意义。