亚麻DMP,基因家族的全基因组鉴定与分析

李雯，赵丽蓉，张建平，刘自刚1，齐燕妮，李闻娟，谢亚萍

（1. 甘肃农业大学农学院，甘肃省干旱生境作物重点实验室，甘肃省作物遗传改良与种质资源创新重点实验室，甘肃 兰州 730070；
2. 甘肃省农业科学院作物研究所，甘肃 兰州 730070）

膜蛋白在细胞增殖与分化、信号转导与识别及物质运输等众多生物过程中发挥重要作用［1-4］。Mori 等［5］发现定位于精子膜的拟南芥（Arabidopsis thaliana）蛋白GEX2 对精子附着并最终完成双受精是必需的。DMPs（domain of unknown function 679 membrane proteins）家族成员是一种绿色植物特有的膜蛋白，通常含有4 个跨膜结构域，具有胞质氨基端和羧基端，参与植物生殖发育及衰老等各种生理过程［6］。目前，已经在拟南芥全基因组水平上鉴定到10 个DMP基因家族成员，不同成员在不同组织器官中具有独特的表达模式，其中AtDMP1，AtDMP2，AtDMP3，AtDMP4和AtDMP7参与不同类型的程序性细胞死亡，包括器官衰老，角果开裂及花器官和角果的脱落，而AtDMP8和AtDMP9在双受精过程中促进配子融合［7］。Zhong 等［8］利用图位克隆鉴定到玉米（Zea mays）DMP基因，并且发现DMP突变植株可诱导母本单倍体的产生。双子叶植物拟南芥AtDMP8和AtDMP9功能缺失也可诱导母本单倍体的产生，而且ZmDMP-like 蛋白在许多双子叶植物中都很保守，序列同源性高达74%，这说明DMP基因突变可以在双子叶植物中触发单倍体诱导［9］。Zhu 等［10］在棉花（Gossypiumspp）中对DMP家族进行了系统分析，发现棉花DMP基因参与植物衰老及生殖过程，其中GhDMP8-A/-D和GbDMP8-A/-D可作为棉花单倍体诱导的候选基因。

亚麻（Linum usitatissimum）属于亚麻科亚麻属，是一种古老的韧皮纤维作物和油料作物［11］。作为人类最早栽培利用的农作物之一，亚麻在全世界被广泛种植，主要的种植国家有中国、加拿大、美国、阿根廷等［12］。按其经济性状与用途可将亚麻分为油用、纤维用和油纤维兼用等3 种类型［13-15］。油用亚麻又称胡麻，富含α-亚麻酸、木酚素、膳食纤维、蛋白质等多种营养物质，其中α-亚麻酸含量达到55%左右，是人类摄取ω-3 脂肪酸的最好来源之一。

随着测序技术的发展，玉米［16］、高粱（Sorghum bicolor）［17］、亚麻［18］、大豆（Glycine max）［19］等多种作物完成了全基因组测序，而基因家族分析也成了挖掘功能基因的有效手段。目前，LBD（lateralorgan boundaries domain）蛋白［20］、天冬氨酸蛋白酶（aspartate protease，AP）［21］、多聚半乳糖醛酸酶（polygalacturonase，PG）［22］等基因家族在很多物种中都被鉴定出来，并进行了系统分析，而DMP基因家族在各种作物中的研究还相对缓慢，亚麻作为重要的油料及经济作物，目前尚无相关报道。本研究旨在鉴定亚麻DMP基因，并对其进行生物信息学及表达模式分析，从而为亚麻育种及功能基因组研究提供基因资源。

1.1 试验材料

本研究选用甘肃省农业科学院育成的胡麻品种“陇亚10 号”“黑亚14 号”“陇亚15 号”（丰产性、抗病性等综合性状优良）［23］作为试验材料。供试材料于2021年3 月种植于甘肃省农业科学院兰州试验基地（36°03′N，103°40′E，海拔1500 m）。在开花期进行标记，分别取开花后10、20、30、40 和50 d 的种子及未成熟花药和成熟花药。在培养皿中萌发陇亚15 号种子，10 d 后选取生长一致的幼苗移至1/2 MS 液体培养基中，3 d 后进行胁迫处理同时更换新的培养基。分别对亚麻幼苗进行盐胁迫（150 mmol·L-1氯化钠）、干旱处理（20%聚乙二醇PEG）、吲哚乙酸处理（15 μmol·L-1IAA）、萘乙酸处理（5 μmol·L-1NAA）、赤霉素处理（5 μmol·L-1GA3）、高温胁迫（45 ℃处理3 h）及低温胁迫（4 ℃处理3 h），分别在处理后0、3、6、9、12、24、48、72 h 及恢复后48 h 取其叶片。所有样品用液氮速冻后均保存于-80 ℃超低温冰箱备用。

1.2 亚麻LuDMP 基因家族的鉴定与克隆

亚麻基因组序列从NCBI 下载（登录号为QMEI02000000），其他基因组注释文件下载地址为https：//doi.org/10.6084/m9. figshare.13614311.v3。利用文献［10］报道的玉米、水稻（Oryza sativa）、葡萄（Vitis vinifera）、大豆及拟南芥DMP基因ID 号，进而通过各自基因组蛋白序列提取相应的DMP 序列，利用Blastp 在亚麻基因组中进行同源比对，并采用Pfam（http：//pfam. sanger. ac. ck/）对结构域进行确认，获得亚麻DMP-like基因。将LuDMP 序列与玉米和拟南芥单倍体诱导相关DMP 蛋白序列进行比对，选择同源性最高的基因作为亚麻单倍体诱导的候选基因并对其进行克隆。以陇亚10 号基因组DNA 和成熟花药的cDNA 为模板分别进行PCR 扩增，PCR 扩增体系为50 μL：Buffer 25 μL；
dNTP 10 μL；
高保真酶KOD 1 μL；
正向引物（10 μmol·L-1）1.5 μL；
反向引物（10 μmol·L-1）1.5 μL；
cDNA 模板（100 ng·μL-1）2 μL；
无RNase 水9 μL。PCR 反应条件：94 ℃2 min；
98 ℃10 s，60 ℃30 s，68 ℃50 s，35 个 循环；
72 ℃7 min。扩增产 物经1% 琼脂糖凝 胶电泳后 酶切回 收，之后 用Zero Blunt® TOPO® PCR Cloning Kit 连接至T 载体并转化DH-5α，将阳性菌液进行测序。测序结果经比对分析后，将正确的基因序列上传至NCBI 数据库。

1.3 亚麻LuDMP 家族序列特征分析

利用ExPASy ProtParam（http：//web.wxoasy.org/protparam/）预测LuDMP 蛋白的分子量和等电点；
并使用Plant-mPLoc（http：//www. csbio. sjtu. edu. cn/bioinf/plant-multi/）和ExPASy ProtScale（https：//web. expasy.org/protscale/）分别预测LuDMP 蛋白的亚细胞定位和亲疏水性，利用TMHMM-2.0（https：// services.healthtech. dtu. dk/service. php？TMHMM-2.0）预测LuDMP 蛋白的跨膜结构域。LuDMP 蛋白结构预测采用SOPMA（https：//npsa-prabi. ibcp. fr/cgi-bin/npsa_automat. pl？ page=/NPSA/npsa_sopma. html）和SWISSMODEL（https：//swissmodel. expasy. org/interactive）。利 用MEME（http：//alternate. meme-suite. org/tools/meme）分析LuDMP 蛋白的保守基序。利用PlantCARE 数据库（http：//bioinformatics. psb. ugent. be/webtools/plantcare/html/）对LuDMP基因启动子的顺式作用元件进行预测。

1.4 亚麻LuDMP 家族成员的染色体分布及复制分析

LuDMP基因染色体定位及可视化采用Map Gene to Chromosome 在线软件。利用MCScanX 分析LuDMP基因的复制事件，之后采用DnaSP 6 计算基因对之间的非同义替换率（nonsynonymous mutation rate，Ka）、同义替换率（synonymous substitution rate，Ks）和Ka/Ks，利用Ka/Ks 分析选择压力。用MYA=Ks/（2×6.1×10-9）×10-6计算基因对之间的分歧时间（million years ago，MYA）［24］。

1.5 DMP 蛋白的系统进化分析

利用MEGA 7.0 软件对亚麻、拟南芥、葡萄、大豆、水稻及玉米6 个物种的DMP 蛋白序列进行比对，比对结果采用邻接法（neighbor-joining；
参数：Bootstrap 1000）构建DMP 蛋白的系统进化树。

1.6 LuDMP 基因的表达分析

用RNA 提取试剂盒（BIOMGA）提取RNA，具体操作步骤按照说明书进行。经NanoDrop-2000 进行浓度和质量检测后，利用TaKaRa 反转录试剂盒将RNA 反转录成cDNA。qRT-PCR 采用TB GreenTMPremix Ex TaqTMII（TaKaRa）荧光定量试剂盒。PCR 反应体系为20 μL：2×mix 10 μL；
上下游引物（10 μmol·L-1）各1 μL；
cDNA 模板2 μL；
无RNase 水6 μL。qRT-PCR 程序为50 ℃2 min；
95 ℃10 min；
95 ℃15 s，60 ℃15 s，72 ℃15 s，40 个循环。每个样品3 次重复，选取GAPDH为内参基因，运用2-ΔΔCt方法［25］计算基因的相对表达量。本研究所用转录组数据从NCBI 下载（登录号为PRJNA505721），利用TBtools 绘制热图，本研究所用引物均见表1（由于LuDMP-1与LuDMP-7序列相似性很高，故二者所用引物相同）。

表1 本研究所用引物Table 1 Primers used in this research

2.1 亚麻LuDMP 基因家族的鉴定与系统进化分析

通过Blastp 比对，并利用Pfam 数据库进行结构域分析，最终在亚麻基因组中鉴定到17 个DMP成员（表2）。根据LuDMP基因在染色体上的位置，将其依次命名为LuDMP-1～LuDMP-17。为深入了解LuDMP 与其他物种DMP 之间的进化关系，本研究将拟南芥、大豆、葡萄、玉米及水稻等5 个物种的DMP 蛋白序列与LuDMP 蛋白序列进行多序列比对并构建系统进化树。结果显示，6 个物种的DMP 蛋白可分为I、II、III、IV 和V 5 个亚族（图1），第III 亚族DMP 数量最多，包含28 个基因，其中LuDMP 11 个，其次为第V、IV、I 亚族，分别含17、14、14 个基因，其中LuDMP 分别为0、2、2 个，第II 亚族DMP 数量最少，包含7 个基因，其中LuDMP 为2 个，每个亚族均含有单子叶与双子叶植物DMP，而且各亚族中的DMP 蛋白具有单子叶或双子叶特异性聚类模式。和其他双子叶植物相比，亚麻DMP基因发生了显著扩张。此外，第IV 亚族中的LuDMP-1和LuDMP-7与单倍体诱导基因AtDMP8、AtDMP9和ZmDMP［8-9］具 有很近的亲缘关系，位于同一个 进化分 支，而 且LuDMP-1 和LuDMP-7 与AtDMP9 及ZmDMP 蛋白的同源性达到67%以上（与AtDMP8 的同源性为57%），这表明它们可能具有相同的功能，可以作为亚麻单倍体诱导的候选基因，而其余LuDMP 成员与ZmDMP 和AtDMP8 的同源性低于46%（表3）。

表3 LuDMP 基因同源性比对Table 3 LuDMP gene homology comparison

图1 DMP 蛋白系统进化树Fig.1 Phylogenetic tree of DMP proteins

表2 亚麻LuDMP 成员基本特征Table 2 Basic characteristics of flax LuDMP members

2.2 亚麻单倍体候选基因LuDMP-1/7 的克隆

对LuDMP-1和LuDMP-7基因进行扩增，只扩增得到LuDMP-1基因的DNA 序列和蛋白质编码区（coding sequence，CDS）序列，大小均为699 bp，与参考序列大小相符（图2），将LuDMP-1基因序列上传至NCBI，GenBank 登录号为OM234689。

图2 LuDMP-1 基因全长扩增Fig.2 Full-length amplification of LuDMP-1M：DL1000 marker；
1：阴性 对 照；
2：DNA 扩 增；
3：cDNA 扩 增。1：Negative control；
2：DNA amplification；
3：cDNA amplification.

2.3 亚麻LuDMP 基因序列分析

亚麻LuDMP 蛋白氨基酸长度为195～240 aa，分子 量 为20946～25830 Da，等 电 点（isoelectric point，PI）为4.76～9.71，所有LuDMP 蛋白的总平均亲水性（grand average of hydropathy，GRAVY）均大于0，即所有LuDMP 蛋白均为疏水性蛋白（表2，图3）。亚细胞定位预测结果表明，除LuDMP-5 和LuDMP-8 定位于叶绿体，其余蛋白均定位于细胞膜，其中LuDMP-1、LuDMP-7 和LuDMP-14 定位于细胞膜和细胞核，LuDMP-11 定位于细胞膜和叶绿体（表2）。此外，LuDMP 蛋白均具有跨膜结构域，数量2～5 个不等（表2）。LuDMP 蛋白二级结构预测结果表明，所有LuDMP 蛋白均含有α-螺旋、β-转角、无规则卷曲以及延伸链（图4）。LuDMP 蛋白种内进化关系如图5A 所示，对LuDMP进行基因结构分析，除LuDMP-9、LuDMP-16和LuDMP-17外，其余LuDMP成员均不含内含子（图5）。对17 个LuDMP 蛋白的保守基序进行分析，共鉴定到8 个motif，其中motif1、motif2、motif3 和motif4 在17 个LuDMP 成员中高度保守，然而一些特定的基序只存在于特定的基因中，如motif5只存在于LuDMP-4、LuDMP-6、LuDMP-13、LuDMP-15、LuDMP-16及LuDMP-17；
motif7 只存在于LuDMP-4、LuDMP-6、LuDMP-13及LuDMP-15；
motif6 存在于除LuDMP-2、LuDMP-3、LuDMP-10及LuDMP-12之外的所有成员中；
而motif8 只存在于LuDMP-1和LuDMP-7中，推测该基序与单倍体诱导相关。

图3 LuDMP 蛋白亲疏水性预测Fig.3 Prediction of LuDMP protein hydrophilicity and hydrophobicity

图4 LuDMP 蛋白二三级结构Fig.4 Secondary and tertiary structure of LuDMP protein

图5 LuDMP 基因系统进化、基因结构与保守基序分析Fig.5 Phylogenetic evolution，gene structure and conservative motifs of LuDMP

2.4 LuDMP 基因染色体定位和复制分析

对17 个LuDMP基因进行染色体定位，结果显示LuDMP-1、LuDMP-7、LuDMP-8、LuDMP-9、LuDMP-10和LuDMP-11分别位于1、7、8、10、12 和14 号染色体，LuDMP-2和LuDMP-3位于2 号染色体，LuDMP-4、LuDMP-5及LuDMP-6位于4 号染色体，LuDMP-12～LuDMP-17均在15 号染色体（图6）。对亚麻LuDMP基因进行复制分析（表4），结果显示有3 对基因发生串联复制事件（LuDMP-4/LuDMP-5、LuDMP-13/LuDMP-14、LuDMP-16/LuDMP-17），8 对 基 因 发 生 片 段 复 制（LuDMP-1/LuDMP-7、LuDMP-2/LuDMP-3、LuDMP-6/LuDMP-9、LuDMP-4/LuDMP-13、 LuDMP-8/LuDMP-9、 LuDMP-8/LuDMP-11、 LuDMP-9/LuDMP-11、 LuDMP-9/LuDMP-15）。通过判断Ka/Ks 进行选择压力分析，由于LuDMP-6/LuDMP-9及LuDMP-9/LuDMP-15的Ks 值为n.a.，导致Ka/Ks 及分歧时间也为n.a.，n.a. 表示这两对基因几乎可发生同义突变的位点都发生了同义突变，表明序列分歧度相当大，进化距离很远。其他发生复制的基因对Ka/Ks 均小于1，表明这些LuDMP基因在进化过程中受到纯化选择。

表4 LuDMP 基因复制分析Table 4 Duplication analysis of LuDMP genes

图6 LuDMP 基因染色体定位Fig.6 Chromosomal locations of LuDMP genes横线代表片段复制，方框代表串联复制。The horizontal line represents segmental duplication，and the box represents tandem duplication. Chr：染色体Chromosome.

2.5 顺式作用元件分析

对亚麻LuDMP基因进行启动子分析，发现LuDMP基因启动子区域除含有TATA-box 等基本作用元件之外，还含有多种与光响应、激素响应、非生物胁迫响应等相关的顺式作用元件（表5）。所有LuDMP基因都含有厌氧诱导必需的作用元件ARE；
每个LuDMP基因至少含有一个或多个光响应元件，如Sp1 等；
除LuDMP-4和LuDMP-11外，其余15 个LuDMP基因均有一个或多个植物激素响应元件，如参与生长素响应的TGA-element，参与赤霉素响应的TATC-box 和GARE-motif，参与水杨酸响应的TCA-element，参与脱落酸响应的顺式作用元件ABRE，以 及 参 与 甲 基 茉 莉 酸（MeJA）响 应 的TGACG-motif 和CGTCA-motif；
除LuDMP-6、LuDMP-7、LuDMP-9、LuDMP-11、LuDMP-13、LuDMP-14和LuDMP-15外，其余LuDMP基因均含有参与玉米醇溶蛋白代谢调控的顺式调控元件O2-site；
参与种子特异性调控的顺式调控元件RY-element 只存在于LuDMP-2中；
除LuDMP-3、LuDMP-4、LuDMP-15、LuDMP-16和LuDMP-17外，其他基因均含有低温响应元件LTR；
参与脱水，低温，盐胁迫的顺式作用元件DRE 只存在于LuDMP-3中；
此外，LuDMP基因启动子区域还存在其他逆境胁迫响应元件，如参与防御和应激反应的TC-rich repeats、参与干旱诱导的MBS。

表5 亚麻LuDMP 启动子顺式元件分析Table 5 Analysis of cis-elements of flax LuDMP promoters

2.6 LuDMP 基因表达模式分析

利用陇亚10 号和黑亚14 号转录组数据分析LuDMP基因在不同品种、不同组织中的表达模式（图7），除LuDMP-9和LuDMP-15外，其他15 个基因均具有表达量数据。由图7 可知，大多数基因在两个亚麻品种中表现为组织特异性表达。其中LuDMP-2、LuDMP-4、LuDMP-5、LuDMP-6、LuDMP-8、LuDMP-11和LuDMP-16在两个品种果实中的表达量高于在茎中的表达量，尤其是LuDMP-4和LuDMP-5在油用品种陇亚10 号中的表达量明显高于纤 维 品 种 黑 亚14 号，而LuDMP-1、LuDMP-3、LuDMP-7、LuDMP-10、LuDMP-12、LuDMP-13、LuDMP-14和LuDMP-17在两个品种茎中的表达量高于果实中的表达量，特别是LuDMP-10和LuDMP-12在黑亚14 号中的表达量高于陇亚10 号。

图7 LuDMP 在不同组织中的表达量热图Fig. 7 Heat map of the expression of LuDMP in different tissues

为进一步验证LuDMP-1/7的功能，本研究分析了该基因在陇亚15 号不同发育时期的种子、未成熟花药及成熟花药中的表达模式。结果显示，LuDMP-1/7基因在各组织中均有表达，但在花药中的表达量明显高于其他组织，尤其在成熟花药中的表达量约为其他组织的6 倍以上（图8）。此外，本研究发现LuDMP-1/7启动子存在激素响应及逆境胁迫相关元件，因此分析了LuDMP-1/7在不同胁迫处理下的表达模式。结果表明（图9），LuDMP-1/7在IAA 处理3 h 后略有降低，6 h 时表达量上升了1.4 倍左右，之后呈先降再升的趋势，处理48 h 后表达量趋于稳定；
NAA 处理后，LuDMP-1/7基因表达量呈先减后增再减的趋势，在24 h 时表达量达到最高，较对照上升1.3 倍左右；
GA3处理后，LuDMP-1/7基因表达量在12 h 时上升了约1.6 倍，之后有所降低，在恢复48 h 时又呈上升的趋势；
用NaCl 处理后，不同时期基因表达均低于对照；
PEG 处理后6 h 表达量达到最高，约为对照的1.3 倍；
在低温（4 ℃）和高温（45 ℃）处理下，LuDMP-1/7基因表达量与对照相比有显著的提高，分别上升了约260倍和600 倍。

图8 LuDMP-1/7 基因在种子及花药中的表达模式Fig.8 Expression patterns of LuDMP-1/7 in seeds and anther不同小写字母表示在0.05 水平具有显著性差异。下同。Different small letters mean significant differences at the 0.05 level.The same below.

图9 胁迫处理下LuDMP-1/7 的表达模式Fig.9 Relative expression patterns of LuDMP-1/7 under stress treatments恢复48 h：48 h after recovery.

DMP 蛋白是一种几乎只存在于绿色植物中的膜蛋白，目前已经在拟南芥［7］和棉花［10］中对其进行了系统分析，此外，在水稻（20 个）、葡萄（6 个）、玉米（15个）、大豆（9 个）等多个物种中分别对其进行了鉴定［10］。本研究在亚麻基因组数据库中筛选获得了17 个DMP基因，不均匀分布于9 条染色体上。将亚麻DMP 蛋白与大豆、拟南芥、葡萄、水稻、玉米等5 个物种的DMP 成员一起构建进化树，发现DMP 蛋白可分为5 个亚族，LuDMP 家族成员在5 个亚族中的数量差异很大，超过65%的LuDMP 成员（11 个）属于第III 亚族。此外，DMP 家族在各亚族中呈现出单子叶或双子叶特异性聚类模式，说明在单子叶植物和双子叶植物分化前DMP基因具有共同的祖先，这与Zhu 等［10］的研究结果一致。串联复制和全基因组复制/片段复制是基因加倍、基因功能特异性和多样化的重要原因之一，也可能导致基因功能冗余［20，26-28］。DMP基因家族的高度扩张主要来源于几次全基因组复制，4 个棉花种中只有3 对基因发生串联复制，其余134 对基因都经历了全基因组复制/片段复制［10］。LuDMP在亚麻进化过程中有3 对基因发生串联复制，8 对基因发生片段复制事件，片段复制是促进LuDMP基因家族扩张的主要原因，这源于亚麻的两次全基因组复制事件［18］。其中第III 亚族的扩张最明显，有9 对基因发生了复制事件，基因占比为73.3%，而亚族I 和IV 各自仅有1 对基因发生了片段复制，说明这两个亚族在进化中更为保守。亚麻中，发生复制的LuDMP基因对Ka/Ks 全部小于1，说明LuDMP基因在亚麻进化过程中经历了纯化选择。对LuDMP 蛋白的亚细胞定位进行预测，发现大多数蛋白定位于细胞膜，与棉花DMP 蛋白预测结果相似［10］，拟南芥中DMP 蛋白定位于液泡膜及内质网［7］，这与DMP 蛋白的功能一致。此外，LuDMP-1、LuDMP-7、LuDMP-11 及LuDMP-14 除了定位于细胞膜，还分别位于细胞核和叶绿体，这说明DMP 蛋白还参与其他生物过程。基因结构分析显示，大多数LuDMP基因不含内含子，与棉花DMP基因家族研究结果相似［10］，说明它们可能具有特殊的功能。保守基序分析结果表明，有4 个motif 存在于所有LuDMP 蛋白，而其他motif 只存在于特定的基因中，这些特有的motif 是导致LuDMP 蛋白功能分化的重要原因。亚麻LuDMP 蛋白具有1 个DMP 结构域及2～5 个跨膜结构域，与棉花中DMP 蛋白的研究结果相同［10］。有研究表明玉米DMP 蛋白第一个跨膜结构域发生单个氨基酸突变可诱导玉米单倍体的产生［8］。拟南芥中也有相似的研究结果，AtDMP8 和AtDMP9 蛋白第一个和第二个跨膜结构域发生突变可诱导拟南芥单倍体的产生［9］。因此，本研究推测与ZmDMP 及AtDMP8 和AtDMP9 亲缘关系相近的亚麻LuDMP-1 和LuDMP-7 蛋白第一个和第二个跨膜结构域也具有相似的功能，可诱导亚麻单倍体的产生。由此可见，DMP 蛋白跨膜结构域在植物育性中起重要作用。

与传统育种方法［29-30］相比，双单倍体育种可在两个世代内获得稳定的二倍体纯系，从而在很大程度上节约时间和成本。利用CRISPR 技术对单倍体诱导相关基因PLA1/MATL和DMP进行编辑可诱导单倍体的产生，这种方法已经成功应用于玉米［31-33］、水稻［34］、小麦（Triticum aestivum）［35］、拟南芥［9］等多个物种。玉米单倍体的产生主要受ZmPLA1/MTL/NLD和ZmDMP的调控［8，31-33］。水稻中OsMATL基因突变可诱导水稻单倍体的产生［34］，AtDMP8和AtDMP9突变可在拟南芥中诱导母本单倍体的产生［9］。在本研究中，亚麻LuDMP-1/7 与AtDMP8、AtDMP9 和ZmDMP 具有很近的亲缘关系，位于同一进化分支，而且序列同源性达到67%以上，这说明它们可能具有相同的功能，参与植物的授粉和受精过程且能够诱导单倍体的产生。此外，利用qRT-PCR 分析LuDMP-1/7基因在亚麻不同组织中的表达模式，发现它们在成熟花药中的表达量最高，进一步说明它们可作为亚麻单倍体诱导的候选基因，但还需进一步对其单倍体诱导功能进行研究。

基因表达通常受到其上游启动子区顺式元件的调控。这些位于基因转录起始位点上游非编码DNA 中的顺式元件调节不同环境下基因的应激表达或组织特异性表达行为。因此，对参与LuDMP基因调控的顺式元件进行分析，有助于了解LuDMP基因的调控机制，并预测其潜在功能。在LuDMP基因启动子区存在许多与生长、逆境胁迫和植物激素响应相关的顺式元件，其中大部分顺式元件参与非生物胁迫和激素信号传导，这与棉花DMP基因顺式作用元件预测结果相似［10］。通过分析基因表达模式可初步预测基因功能。组织特异性分析结果表明，有7 个LuDMP基因在果实中的表达量高于茎中的表达量，特别是LuDMP-4和LuDMP-5在油用品种中的表达量显著高于纤维品种，说明这些基因可能与果实发育或含油率相关；
而其他8 个基因在茎中的表达量较高，尤其是LuDMP-10和LuDMP-12在黑亚14 号中的表达量高于陇亚10 号，推测它们可能参与茎的生长发育（株高、纤维合成等）。此外，由于LuDMP-1/7基因启动子具有激素响应及非生物胁迫相关元件，推测其可能参与亚麻激素响应及逆境胁迫等生理过程。本研究中除NaCl 胁迫外，其他处理下LuDMP-1/7基因表达量均有不同程度的提高，其中基因表达量对低温和高温的响应最明显，分别较对照提高了约260 倍和600 倍，表明该基因参与调控亚麻对激素和逆境胁迫的响应，但具体的功能还需进一步验证。

猜你喜欢亚族单倍体亚麻二穗短柄草CYP72A亚族成员表达分析及亚细胞定位江苏农业科学(2022年17期)2022-09-24甜瓜Trihelix转录因子的鉴定与生物信息学分析江苏农业科学(2022年16期)2022-08-26苦荞蛋白磷酸酶2C家族的鉴定及表达分析四川大学学报（自然科学版）(2021年4期)2021-07-15不同除草剂对玉米单倍体成熟胚的加倍效果山西农业科学(2020年8期)2020-08-13辣椒HD-Zip基因家族鉴定、系统进化及表达分析中国农业科学(2020年5期)2020-03-30亚麻抗白粉病种质资源的鉴定与筛选中国麻业科学(2018年6期)2018-04-09玉米单倍体育性自然恢复研究进展河南农业科学(2018年5期)2018-01-19新申亚麻: 从展会到秀场纺织科学研究(2017年4期)2017-05-17新申亚麻&孙林：亚麻时尚的民俗融合纺织服装流行趋势展望(2016年6期)2016-05-04“逐梦森林”：新申亚麻&孙林时装发布会纺织服装流行趋势展望(2016年2期)2016-05-04

相关热词搜索：基因组 亚麻 基因