基于SNP遗传图谱对甘蓝型油菜部分脂肪酸组成性状的QTL定位

叶桑，崔翠，郜欢欢，雷维，王刘艳，王瑞莉，陈柳依，曲存民，唐章林，李加纳，周清元

基于SNP遗传图谱对甘蓝型油菜部分脂肪酸组成性状的QTL定位

叶桑，崔翠，郜欢欢，雷维，王刘艳，王瑞莉，陈柳依，曲存民，唐章林，李加纳，周清元

（西南大学农学与生物科技学院，重庆 400715）

【】菜籽油在烹饪、食品加工及工业生产中广泛应用，因此，根据生产需要改善菜籽油脂肪酸组分是油菜育种的重要目标。通过对2种环境下甘蓝型油菜主要脂肪酸组成进行QTL定位分析，寻找甘蓝型油菜脂肪酸组分的QTL及影响本群体脂肪酸组分的候选基因。以人工合成甘蓝型油菜10D130和甘蓝型油菜常规品种中双11构建高世代重组自交系（RIL）为研究材料，分别于2016—2017年和2017—2018年2个年度在重庆市北碚区2个不同的环境中设置田间试验，收获自交种子，采用气相色谱法3次重复对种子的脂肪酸组分进行分析。利用油菜6K SNP芯片对该RIL群体进行基因分型，DNA样品预处理及芯片处理严格按照Illumina Inc 公司Infinium HD Assay Ultra操作说明进行。取最小阈值LOD 2.0利用JoinMap4.0软件构建高密度遗传连锁图谱。通过QTL IciMapping V4.1完备区间作图法对油菜主要脂肪酸组成进行QTL定位。2种环境中，两亲本各性状间差异及RIL群体各性状在株系间差异均达到显著或极显著水平，且6种脂肪酸含量在2个环境中均表现为连续分布，适合进行QTL检测。构建用于QTL定位的遗传图谱包含1 897个多态性SNP标记，覆盖甘蓝型油菜基因组3 214.19 cM，平均图距1.69 cM。利用此图谱，在2个环境共检测到位于8条染色体上的23个控制脂肪酸组分QTL位点，与硬脂酸、油酸、亚油酸、亚麻酸、廿碳烯酸和芥酸含量相关的QTL位点分别为6、3、4、5、2和3个，其中在A05、A08和C03染色体上发现多种脂肪酸含量的QTL“富集区”。在A05染色体上检测到亚油酸和亚麻酸含量重叠的主效QTL，亚油酸与亚麻酸表现加性效应相同；在A08和C03上都检测到油酸、廿碳烯酸和芥酸含量重叠的主效QTL，油酸与廿碳烯酸及芥酸表现加性效应相反。与拟南芥脂肪酸代谢基因进行同源性比对分析，在17个QTL置信区间内筛选到22个候选基因，主要通过编码脂肪酸去饱和酶、全羧化酶合酶、碳链延长酶和参与酰基辅酶A生物合成等途径调控脂质的生物合成和代谢。利用甘蓝型油菜6K SNP芯片准确定位了2种环境条件脂肪酸组成的QTL位点，筛选到位于A05、A08和C03染色体上多种脂肪酸QTL的“富集区”，并与拟南芥脂肪酸代谢基因比对出该群体油菜脂肪酸代谢基因，可作为改善油菜籽脂肪酸组成的重要区段及候选基因。

甘蓝型油菜；脂肪酸；遗传图谱；QTL；候选基因

0 引言

【研究意义】油菜是中国主要的油料作物，其种植面积和总产量居世界首位[1]，目前，已经成为继水稻、小麦、玉米之后的第四大作物，是中国的优势油料作物[2]。菜籽油是中国主要的食用油之一，其主要脂肪酸组成包括硬脂酸（C18﹕0）、油酸（C18﹕1）、亚油酸（C18﹕2）、亚麻酸（C18﹕3）、廿碳烯酸（C20﹕1）和芥酸（C22﹕1）等。食用油中油酸具有极高的营养和应用价值，高油酸菜籽油不仅有利于人体心血管健康[3]，还因其热稳定性高[4]，利于纯化储存而被用于煎炸行业[5]以及生产、生活用品的原料[6]。亚油酸是一种必需脂肪酸，可防止动脉粥样硬化，但存在氧化稳定性差、储存时间短的缺点。亚麻酸属多烯类不饱和脂肪酸，因存在3个不饱和键而易氧化变质，故亚麻酸含量较高的菜籽油不耐储存[7]。芥酸因其碳链较长，在人体内不易被消化吸收，不适用于食品加工[8]，然而芥酸是重要的工业原料，被广泛应用于众多工业领域[7]。因此，根据生产需要改善油菜脂肪酸的组成成分是油菜育种的重要目标，定位油菜脂肪酸组分的QTL和筛选影响脂肪酸组分的候选基因在油菜育种中都具有重要的指导作用。【前人研究进展】前人研究表明，在油菜主要脂肪酸的形成过程中，首先由棕榈酸在碳链延长酶的作用下合成硬脂酸[9]，硬脂酸在质体基质中再被硬脂酰-ACP脱氢酶催化形成油酸[10]，而油酸包含2条代谢途径：其一，油酸在碳链延长酶的作用下合成廿碳烯酸，廿碳烯酸进一步延长碳链生成芥酸[11]；其二，油酸在△12减饱和酶作用下生成亚油酸[12]，亚油酸在△15减饱和酶作用下进一步减饱和生成亚麻酸[13]。了解主要脂肪酸的代谢途径是改善油菜脂肪酸组成的前提，通过分子标记定位相关性状QTL并筛选出候选基因有利于更准确了解脂肪酸各组分受到影响的遗传因素。迄今为止，利用分子标记技术的QTL定位分析在油菜相关性状包括抗性[14-16]、育性[17-18]和品质[19-21]等方面广泛发展，为油菜分子的育种改良奠定了基础。油菜脂肪酸含量具有较高的遗传力，对其相关分子标记的确定是改善油菜脂肪酸组成的关键。Burns等[22]利用RFLP分子标记进行油菜脂肪酸组成的QTL定位，获得4个与油酸含量相关的QTL，其中效应较大的QTL位于N3、N8、N18染色体，效应较小的位于N11染色体；获得5个与亚油酸含量相关的QTL，效应相对较大QTL位于N8和N14染色体；获得5个与亚麻酸含量相关的QTL，分别位于N6、N7、N18和N11染色体。张洁夫等[23]采用RAPD、SSR和SRAP 3种分子标记对油菜脂肪酸组成进行QTL定位，获得3个与硬脂酸含量相关的QTL，分别位于N1、N8和N16染色体；获得2个与油酸含量相关的主效QTL，分别位于N8和N13染色体；获得3个与亚油酸含量相关的QTL，其中2个位于N8染色体，1个位于N13染色体；获得3个与亚麻酸含量相关的微效QTL；获得4个与廿碳烯酸含量相关的QTL，分别位于N8、N13和N15染色体；在N8和N13染色体获得2个与芥酸含量相关的主效QTL。Zhao等[24]利用SSR标记构建的遗传图谱找到了7个控制油酸含量的QTL，其中，主效QTL位于N18染色体；找到7个控制亚油酸含量的QTL，其中主效QTL位于N9染色体；并在N14染色体找到1个控制亚麻酸含量的QTL。与上述分子标记技术相比，SNP标记具有数量多、分布广和规模化等优点，已成为大家认可的第三代遗传标记，近几年被广泛应用于油菜产量、品质等性状的分子标记辅助育种及遗传定位等诸多领域。例如，Liu等[25]利用油菜SNP芯片技术构建高密度遗传图谱，对油菜籽木质素、纤维素和半纤维素等性状进行QTL定位。Shen等[26]采用SNP分子标记技术定位了17个油菜分枝角度QTL，并在QTL区间获得27个候选基因。Qu等[27]利用520份油菜资源群体对油菜籽脂肪酸组成进行全基因组关联分析，确定了62个与7种脂肪酸组成显著相关的基因组区域，并鉴定了24个参与脂肪酸生物合成的功能候选基因的同源基因。【本研究切入点】前人对油菜脂肪酸含量的定位研究主要是利用SSR、AFLP、RFLP和SRAP等分子标记，存在标记的多位点现象，给不同连锁图谱之间QTL位点的比较分析带来困难。另外，前人在QTL作图方法上大多采用区间作图法和复合区间作图法，致使QTL效应可能会被侧连标记区间之外的标记变量所吸收[28]。【拟解决的关键问题】本研究利用油菜6K SNP芯片构建高密度遗传连锁图谱，通过完备区间作图法对油菜主要脂肪酸组成进行QTL定位。结合2年的表型鉴定，发掘在不同环境条件下稳定表达的主效QTL，并利用甘蓝型油菜基因组序列，与拟南芥脂肪酸相关代谢基因进行同源性分析，筛选可能的候选基因，旨在探讨油菜主要脂肪酸之间的内在联系，为改善油菜脂肪酸组成进行分子标记辅助选择、克隆相关候选基因奠定基础。

1 材料与方法

1.1 材料

重组自交系群体的母本为人工合成甘蓝型油菜新品系10D130，其芥酸含量50%左右，父本为高产优质的天然甘蓝型油菜常规品种中双11号（简称ZS11），其芥酸含量＜1%。两亲本杂交F2通过单粒法连续自交7代，构建重组自交系群体，以其中的186个材料进行SNP标记分析，构建高密度SNP遗传连锁图谱。所有试验材料均由西南大学重庆市油菜工程技术研究中心提供。

1.2 田间试验

2016—2017年和2017—2018年2个年度，在重庆市北碚区歇马镇油菜实验基地种植重组自交系群体及其亲本，2个环境分别记为17Cq和18Cq。田间播种采用随机区组试验设计，每个小区5行，每行10株。田间管理遵循常规生产方式，在油菜初花期逐株套袋自交，种子完全成熟后每家系收取5株正常植株的种子，自然风干后保存。

1.3 脂肪酸组成的测定

根据文献[29]的方法，利用气相色谱仪GC-2010（Shimazu，JAP）测定2个环境下重组自交系群体种子中硬脂酸、油酸、亚油酸、亚麻酸、廿碳烯酸及芥酸的含量，每份材料重复测定3株，取其平均值。

1.4 遗传连锁图谱构建及QTL分析

参考刘列钊等[30]方法，分别从每个株系的5个幼嫩植株叶片取混合样0.15 g，利用DNA提取试剂盒DP321-03（天根，中国北京）提取DNA并稀释至50 ng·μL-1用于SNP标记分析。严格按照Infinium HD Assay Ultra操作说明书（Illumina Inc 公司）进行DNA样品的预处理（等位扩增、片段化及富集）、与芯片杂交、洗涤、安装流动室、单碱基延伸、染色及包埋。芯片准备好后运用Illumina HiScan扫描仪的iScan Control Software软件扫描，然后利用GenomeStudio genotyping software v2011软件分析扫描结果，获取各个样本的SNP基因型数据，并为获得的SNP标记命名，命名方法以S-95505568为例，S表示SNP，95505568代表GenomeStudio genotyping software生成的相应SNP位点索引号。

连锁图谱的构建采用JoinMap4.0软件[31]，选用Kosambi[32]函数将重组值转换为图距单位（centiMorgan，cM），取最小阈值LOD2.0对所有标记进行分组，每个染色体上标记顺序通过两两标记之间最小重组频率计算，构建用于QTL定位的遗传连锁图谱。

采用软件QTL IciMapping V4.1[33]完备区间作图（inclusive composite interval mapping，ICIM）法，进行2个环境下硬脂酸、油酸、亚油酸、亚麻酸、廿碳烯酸及芥酸含量的QTL定位和检测。进行ICIM分析时，选用1 cM的步长（walking speed），作1 000次排列检验和显著水平0.01确定LOD临界值。软件运行结果可同时给出性状QTL的加性效应和表型贡献率。QTL命名时以q加上性状的英文缩写再加上染色体编号及QTL序号表示，字体斜体，如代表亚油酸含量位于C03染色体上的第一个QTL。同一性状在染色体相同的位置检测到重复的QTL，且加性效应方向一致，认为是同一QTL。

1.5 候选基因筛选

为筛选出与脂肪酸组成相关的候选基因，首先将检测到的QTL置信区间在甘蓝型油菜基因组上对应的序列查询到，然后与Botella等[34]、Qu等[27]搜索出的拟南芥脂肪酸代谢相关基因进行BlastN比对，将E值设定为E-20，最后筛选出每个QTL置信区间内匹配E值小于阈值的候选基因。

2 结果

2.1 重组自交系群体脂肪酸含量的分析

两亲本和RIL群体的表型分析见表1。RIL群体的表型在2个环境中表现一致：6个性状的均值都介于两亲本之间，最大值和最小值都超过亲本值，存在明显的超亲分离现象。另外，两亲本各性状间差异及RIL群体各性状在株系间差异都达到显著或极显著水平，结合变异系数值反映出这几个性状在该群体中具有较广泛的遗传变异。由性状频率分布图（图1）可知，硬脂酸、亚油酸和亚麻酸呈现偏态分布的特征，油酸、廿碳烯酸和芥酸含量则表现为明显的多峰分布，结果表明，6种脂肪酸含量在2个环境中均表现为连续分布，适合进行QTL检测。

由表2可知，在2个环境下油酸含量与硬脂酸及亚油酸含量均表现为极显著正相关（＜0.01），与廿碳烯酸及芥酸含量均表现为极显著负相关（＜0.01），与亚麻酸含量相关性不强。芥酸含量与廿碳烯酸含量下2个亲本的脂肪酸含量在2个环境中均表现极显著正相关（＜0.01），相关系数分别为0.56和0.64；与亚油酸含量在2个环境中均表现极显著负相关（＜0.01），相关系数分别为-0.52和-0.54。结果表明，油酸转化为亚油酸的途径受到芥酸途径的掣肘，即油酸和亚油酸的含量受芥酸含量的变化影响较大。而亚麻酸受油酸和芥酸变化的影响较小，在遗传上具有相对独立的稳定性。

表1 2个环境下油菜主要脂肪酸含量的变化

a、b和A、B分别代表在0.05和0.01水平上差异显著，*和**分别代表在0.05和0.01水平上差异显著。C18:0、C18:1、C18:2、C18:3、C22:1和C22:2分别代表硬脂酸、油酸、亚油酸、亚麻酸、廿碳烯酸和芥酸。下同

a,bandA,Bindicate significance at＜0.05 and＜0.01, * and ** indicate significance at＜0.05 and＜0.01. C18:0, C18:1, C18:2, C18:3, C22:1 and C22:2 represent stearic acid, oleic acid, linoleic acid, linolenic acid, eicosenoic acid and erucic acid, respectively. The same as below

P1、P2和P’1、P’2分别标出了17Cq和18Cq环境

P1, P2and P'1, P'2indicated the fatty acid content of two parents under 17Cq and 18Cq environment respectively

图1 RIL群体硬脂酸、油酸、亚油酸、亚麻酸、廿碳烯酸及芥酸含量的频率分布

Fig. 1 Frequency distribution of stearic acid, oleic acid, linoleic acid, linolenic acid, eicosenoic acid and erucic acid content from the RIL population in two different environments

表2 2个环境下主要脂肪酸含量相关分析

右上部代表17Cq环境，左下部代表18Cq环境

The upper right represents the 17Cq environment and the lower left represents the 18Cq environment

2.2 遗传连锁图谱

利用包含5 058个标记的6K油菜芯片对186个RIL材料进行基因型鉴定，从中筛选出1 897个高质量多态性的SNP标记，约占37.5%，基于这些标记构建遗传连锁图谱。获得用于QTL定位的图谱覆盖甘蓝型油菜基因组3 214.19 cM，平均图距为1.69 cM。每条染色体长度在86.51—298.72 cM，平均长度为169.17 cM。染色体标记数目在24—153，平均数目为99.84个。但是从标记分布来看，各染色体分布不均，其中染色体A03、C03和C04上标记分布较多，分别为139、147和153个；而染色体C05、C08和C09上的标记数目较少，分别只有48、26和24个。此外各染色体标记密度也有较大差异，C04密度最大，平均间距仅1.10 cM，而密度最小的C08平均间距达6.13 cM。

2.3 主要脂肪酸含量的QTL分析

2.3.1 硬脂酸含量的QTL 在2个环境中，硬脂酸含量共定位到6个QTL，解释表型变异5.08%—25.49%，其中主效QTL位点位于A08和C03上。C03上的主效位点在2个环境中能被重复检测到，分别解释表型变异的11.59%和9.12%；在A08上的2个QTL位点分别解释硬脂酸含量表型变异的14.00%和25.49%。且A08和C03上的QTL位点在2个环境中的加性效应均小于0（表3和图2），即加性效应均来源于亲本ZS11。

2.3.2 油酸含量的QTL 在2个环境中共检测到3个与油酸含量相关的 QTL，且3个位点在2个环境中均能被重复检测到，其中位于A08和C03染色体上的是主效位点，可解释油酸含量表型变异的28.28%—43.68%。在2个环境中，位于A05上的QTL位点加性效应分别为3.81和3.91，而A08和C03染色体上主效位点的效应值在-11.42—-8.66（表3和图2），即前者的加性效应来源于亲本10D130，后两者的加性效应则来自亲本ZS11。

2.3.3 亚油酸含量的QTL 在2个环境中，亚油酸含量共检测到4个QTL，除了C03染色体上的以外，其余3个QTL均能在2个环境中被重复检测到。其中主效QTL位于A05和C03上，可解释亚油酸含量表型变异的12.56%—25.25%；A08染色体上重复检测到的QTL位点在2个环境中的遗传贡献率分别为8.91%和10.49%。同时，在2个环境中重复检测到的3个QTL位点的效应值在-2.37—-1.18（表3和图2），即加性效应均来自亲本ZS11。

2.3.4 芥酸含量的QTL 在2个环境中共定位到3个与芥酸含量相关的QTL，其中位于A08和C03上的主效QTL位点能够在2个环境中被重复检测到，可解释表型变异16.45%—67.26%。芥酸含量的3个QTL位点的加性效应均为正值，加性效应均来源于亲本10D130，其中，的效应值最大，为17.53（表3和图2），属主基因座位。其他主要脂肪酸的QTL定位结果见表3和图2。

2.4 主要脂肪酸代谢候选基因筛选

将23个QTL置信区间序列与拟南芥脂肪酸代谢基因进行比对，分别在A01、A02、A05、A08、A09、A10和C03染色体上的17个QTL区间内检测到22个候选基因，匹配E值介于0—4E-17（表4）。位于A01染色体上的基因和分别与拟南芥基因和同源，参与辅酶A的生物合成。A02染色体上的基因、和分别与拟南芥基因、和同源，分别编码长链酰基辅酶A合成酶、β-酮酰辅酶A还原酶和全羧化酶合酶。A05染色体上的基因、、、和分别与拟南芥基因、、、和同源，第一个编码脂肪酸脱氢酶，后4个参与脂质代谢过程。在A09染色体上的基因与拟南芥基因同源，编码脂肪酸去饱和酶。A10染色体上的基因和分别与拟南芥基因和同源，前者是丙酮酸激酶家族蛋白成员，后者是脂质转运蛋白家族成员。A08和C03染色体上存在多个同源基因，其中和与拟南芥基因同源，和与拟南芥基因同源，两者同属3-酮酰辅酶A合成酶家族；和与拟南芥基因同源，编码一种脂肪酰辅酶A还原酶。

3 讨论

3.1 QTL定位结果的比较分析

油菜主要脂肪酸QTL定位研究前人已有报道，本研究的部分结果与已有研究结果一致性较高，但也发现了一些新的QTL区域。例如，本研究通过完备区间作图法发现了控制芥酸含量的3个QTL，其中位于A08和C03染色体上的2个主效QTL在2个环境下被重复检测到，第3个微效QTL也位于A08染色体，与前人将控制芥酸含量的QTL主要集中在A08和C03上的结果一致[23,30,35-36]。油酸的遗传比较复杂，本研究在2个环境中重复检测到位于A08和C03染色体上油酸含量的2个主效QTL，和芥酸的主效QTL重叠，加性效应相反，这与刘列钊等[30]的研究结果相同。不同的是，本研究亲本并非高油酸材料，却在2个环境中重复定位到位于A05染色体上油酸含量的1个微效QTL，该位点同时还是亚油酸的主效QTL位点，这在前人研究中未发现。本研究在该位点上筛选到候选基因，同时检测到该群体（RIL）中的4个高油酸材料的第722位碱基A突变成T导致酪氨酸变成苯丙氨酸，从而影响了蛋白质的空间结构，导致油酸含量提高，这与张宏军等[37]的研究结果一致。而亲本ZS11的正常，结合A05染色体上QTL位点的加性效应来源于亲本10D130，猜测的变异来自亲本10D130。亚麻酸含量是典型的数量性状，易受环境影响[30]。在17Cq环境中，亚麻酸含量检测到1个位于C07染色体上的主效QTL，解释表型变异的11.66%，而在18Cq环境中，检测到一个主效QTL位于A05染色体，解释亚麻酸含量表型变异的23.66%，这与前人主要将亚麻酸含量主效QTL定位在A04和C04染色体不同[12,24,38]，可能的原因是受环境因素影响较大。本研究采用完备区间作图法对甘蓝型油菜RIL群体的一系列油菜籽脂肪酸含量进行QTL定位，一方面进一步验证了已有研究的可靠性，另一方面也说明这些能够被多次重复定位的主效区段和新发现的基因组区段均是控制油菜籽主要脂肪酸含量的重要区域，对剖析油菜籽脂肪酸含量的遗传具有重要意义。

表3 油菜主要脂肪酸含量在2个环境中的QTL位点

1)：标记在遗传连锁图谱上的连锁位置

1): linkage location of markers on the genetic linkage map

表4 在甘蓝型油菜脂肪酸含量QTL置信区间比对拟南芥相关基因获得的候选基因

3.2 主要脂肪酸含量QTL之间的相关性

研究主要对甘蓝型油菜主要脂肪酸，包括硬脂酸、油酸、亚油酸、亚麻酸、廿碳烯酸及芥酸含量进行QTL定位分析。结果表明，控制这些脂肪酸含量的QTL之间存在明显的相关性，主要表现为A05、A08和C03染色体上出现QTL的“富集区”。在A05染色体的相同位置同时存在硬脂酸、油酸、亚油酸和亚麻酸4种脂肪酸的QTL，其中对于亚油酸和亚麻酸是主效QTL，且加性效应方向一致，可解释二者在2个环境中所表现的0.45和0.49的极显著正相关性。在A08和C03染色体上相同的置信区间内同时检测出油酸、廿碳烯酸和芥酸的主效QTL位点，即控制油酸含量的主效QTL同时也是控制廿碳烯酸和芥酸含量的主效QTL，油酸的2个主效位点加性效应值都为负，而廿碳烯酸和芥酸在此位点上的加性效应则表现为正值，充分解释了油酸和廿碳烯酸含量之间、油酸和芥酸含量之间在2个环境中都表现的极显著负相关性以及廿碳烯酸和芥酸含量在2个环境中的极显著正相关性。

根据表型相关性分析和QTL相关性分析认为，油酸的2条代谢途径并不均衡，其“加碳”途径比“脱氢”途径更具优势。一方面表现在A08和C03上芥酸的2个主效QTL对于油酸也是主效位点，可解释油酸表型变异的28.28%—43.68%，而A05上亚油酸的主效QTL对于油酸是微效位点，只解释油酸表型变异的3.8%—5.7%；另一方面，A08染色体上芥酸的主效QTL同时也是亚油酸的QTL，可解释亚油酸表型变异的8.91%，而亚油酸在A05上的主效QTL并不能解释芥酸含量的表型变异，这就表明，“加碳”途径制约着亚油酸的含量，而“脱氢”途径则对芥酸的含量影响不大。在育种工作中，对亚油酸含量的提高，需要阻断油酸向芥酸转化的途径；而高油酸的获得，则需要对2条代谢途径同时进行抑制。所以针对这些区段的深入研究和挖掘为改善油菜脂肪酸组成奠定基础。

3.3 脂肪酸代谢候选基因

本研究在油菜籽脂肪酸组分QTL置信区间内发现22个脂肪酸代谢相关基因。这些候选基因主要集中在A05、A08和C03染色体上，其中位于A05染色体上的基因（）编码脂肪酸脱氢酶，是调控多不饱和脂肪酸合成的关键基因[39]。其中坐落于A08和C03染色体上的同源基因（和）是参与长链脂肪酸合成的关键基因，已被很多学者研究报道[40-42]，这表明我们的QTL定位是可靠的。在A08和C03染色体上筛选的另外2个同源基因（和）和（和），通过各自编码一种醇形成的脂肪酰辅酶A还原酶参与饱和脂肪酸的生物合成[43-44]。位于A01染色体上的基因（）和（）参与辅酶A的生物合成，有研究[45]表明，辅酶A可与醋酸盐结合为脂肪酸从头合成的前体物质乙酰辅酶A。在A02染色体上的候选基因与拟南芥基因同源，已被证实通过催化脂肪酸、ATP和辅酶A形成酰基辅酶A，参与拟南芥籽油的生物合成[46]。此外，本研究还在QTL置信区间内发现许多未知功能基因，可能存在油菜特有的脂肪酸代谢功能基因，有待进一步挖掘。

4 结论

在2个环境中共检测到位于A01、A02、A05、A08、A09、A10、C03和C07染色体上的23个脂肪酸组分QTL位点，与硬脂酸、油酸、亚油酸、亚麻酸、廿碳烯酸和芥酸含量相关的QTL位点分别为6、3、4、5、2和3个，其中在A05、A08和C03染色体上发现多种脂肪酸含量的QTL“富集区”，可作为改善油菜籽脂肪酸组成的重要区段。通过与拟南芥脂肪酸代谢基因的比对分析，在QTL置信区间内筛选到22个脂肪酸代谢相关候选基因。

[1] 官春云. 中国油菜产业发展方向. 粮食科技与经济, 2011, 36(2): 5-6.

GUAN C Y. Development direction of rapeseed industry in China., 2011, 36(2): 5-6.(in Chinese)

[2] 沈金雄, 傅廷栋, 涂金星, 马朝芝. 中国油菜生产及遗传改良潜力与油菜生物柴油发展前景. 华中农业大学学报, 2007, 26(6): 894-899.

SHEN J X, FU T D, TU J X, MA C Z. Potential in production and genetic improvement of rapeseed and prospect for rape oil-based biodiesel in China., 2007, 26(6): 894-899. (in Chinese)

[3] CHANG N W, HUANG P C. Effects of the ratio of polyunsaturated and monounsaturated fatty acid to saturated fatty acid on rat plasma and liver lipid concentrations., 1998, 33: 481-487.

[4] PETUKHOV I, MALCOBNSON L J, PRZYBYLSKI R, ARMSTRONG L. Frying performance of genetically modified canola oils., 1999, 76: 627-632.

[5] MILLER M R, BRIDLE A R, NICHOLS P D, CARTER C G. Increased elongase and desaturase gene expression with stearidonic acid enriched diet does not enhance long-chain (n-3) content of seawater Atlantic salmon (L.)., 2008, 138(11): 2179-2185.

[6] PIAZZA G J, FOGLIA T A. Rapeseed oil for oleochemical uses., 2001, 103: 405-454.

[7] 傅寿仲, 张洁夫, 戚存扣, 浦惠明, 高建芹, 陈新军. 工业专用型高芥酸油菜新品种选育. 作物学报, 2004, 30(5): 409-412.

FU S Z, ZHANG J F, QI C K, PU H M, GAO J Q, CHEN X J. Breeding of high erucic acid rapeseed (L.) for industrial use., 2004, 30(5): 409-412. (in Chinese)

[8] MARTINEZ-RIVAS J M, SPERLING P, LUHS W, HEINZ E. Spatial and temporal regulation of three different microsomal oleate desaturase genes (FAD2) from normal-type and high-oleic varieties of sunflower (L.)., 2001, 8(2): 159-168.

[9] PIDKOWITCH M S, NGUYEN H T, HEILMANN I, ISCHEBECK T, SHANKLIN. Modulating seed beta-ketoacyl-acyl carrier protein synthase II level converts the composition of a temperate seed oil to that of a palm-like tropical oil.AcademySciences, 2007, 104(11): 4742-4747.

[10] SAHA S, ENUGUTTI B, RAJAKUMARI S, RAJASEKHARAN R. Cytosolic triacylglycerol biosynthetic pathway in oilseeds: molecular cloning and expression of peanut cytosolic diacylglycerol acyltransferase., 2006, 141: 1533-1543.

[11] BAUD S, LEPINIEC L. Physiological and developmental regulation of seed oil production., 2010, 49(3): 235-249.

[12] HU X Y, SULLIVAN-GILBERT M, GUPTA M, THOMPSON S A. Mapping of the loci controlling oleic and linoleic acid contents and development of fad2 and fad3 allele-specific markers in canola (L.)., 2006, 113: 497-507.

[13] PHAM A T, SHANNON G J, BILYEU K D. Combinations of mutant FAD2 and FAD3 genes to produce high oleic acid and low linolenic acid soybean oil., 2012, 125(3): 503-515.

[14] YONG H Y, WANG C, BANCROFT I, LI F, WU X, KITASHIBA H, NISHIO T. Identification of a gene controlling variation in the salt tolerance of rapeseed (L.)., 2015, 242(1): 313-26.

[15] FIETCHER R S, HERRMANN D, MULLEN J L, LI Q F, SCHRIDER D R, PRICE N, LIN J J, GROGAN K, KERN A, MCKAY J K. Identification of polymorphisms associated with drought adaptation QTL inby resequencing., 2016, 6(4): 793-803.

[16] BEHLA R, HIRANI A H, ZELMER C D, YU F Q, DILANTHA- FERNANDO W G, MCVETTY P, LI G Y. Identification of common QTL for resistance toin three doubled haploid populations of(L.)., 2017: 213-260.

[17] HUANG J X, XIONG H X, PAN B, NI X Y, ZHANG X Y, ZHAO J Y. Mapping QTL of flowering time and their genetic relationships with seed weight in., 2016, 49(16): 3073-3083.

[18] YE J, YANG Y H, CHEN B, SHI J Q, LUO M Z, ZHAN J P, WANG X F, LIU G H, WANG H Z. An integrated analysis of QTL mapping and RNA sequencing provides further insights and promising candidates for pod number variation in rapeseed (L.)., 2017, 18(1): 18-71.

[19] WEN J, XU J F, LONG Y, WU J G, XU H M, MENG J L, SHI C H. QTL mapping based on the embryo and maternal genetic systems for non-essential amino acids in rapeseed (L.) meal., 2016, 96(2): 465-473

[20] HUANG J X, CHEN F, ZHANG H Z, NI X Y, WANG Y L, LIU H, YAO X T, XU H M, WANG H, MENG J L, ZHAO J Y. Dissection of additive, epistatic and QTL × environment effects involved in oil content variations in rapeseed., 2017, 136(4): 728-737.

[21] RABOANATAHIRY N, CHAO H B, DALIN H, PU S, YAN W, YU L J, WANG B S, LI M T. QTL alignment for seed yield and yield related traits in., 2018, 9: 1127.

[22] BURNS M J, BARNES S R, BOWMAN J G, CLARKE M H, WERNER C P, KEARSEY M J. QTL analysis of an intervarietal set of substitution lines in: (i)Seed oil content and fatty acid composition., 2003, 90(1): 39-48.

[23] 张洁夫, 戚存扣, 浦惠明, 陈松, 陈锋, 高建芹, 陈新军, 顾慧, 傅寿仲. 甘蓝型油菜主要脂肪酸组成的QTL定位. 作物学报, 2008, 34(1): 54-60.

ZHANG J F, QI C K, PU H M, CHEN S, CHEN F, GAO J Q, CHEN X J, GU H, FU S Z. QTL identification for fatty acid content in rapeseed (L.)., 2008, 34(1): 54-60. (in Chinese)

[24] ZHAO J Y, DIMOV Z, BECKER H C, ECKE W, MOLLERS C. Mapping QTL controlling fatty acid composition in a doubled haploid rapeseed population segregating for oil content., 2008, 21(1): 115-125.

[25] LIU L Z, QU C M, WITTKOP B, YI B, XIAO Y, HE Y J, SNOWDON R J, LI J N. A high-density SNP map for accurate mapping of seed fibre QTL inL.., 2013, 8: e83052.

[26] SHEN Y S, YANG Y, XU E S, GE X H, XIANG Y, LI Z Y. Novel and major QTL for branch angle detected by using DH population from an exotic introgression in rapeseed (L.)., 2018, 131(1): 67-78.

[27] QU C M, JIA L D, FU F Y, ZHAO H Y, LU K, WEI L J, XU X F, LIANG Y, LI S M, WANG R, LI J N. Genome-wide association mapping and Identification of candidate genes for fatty acid composition inL. using SNP markers., 2017, 18(1): 232

[28] 王健康. 数量性状基因的完备区间作图方法. 作物学报, 2009, 35(2): 1-7.

WANG J K. Inclusive composite interval mapping of quantitative trait genes., 2009, 35(2): 1-7. (in Chinese)

[29] RUCKER B, ROBBELEN G. Impact of low linolenic acid content on seed yield of winter oilseed rape (L.)., 1996, 115(4): 226-230.

[30] 刘列钊, 李加纳. 利用甘蓝型油菜高密度SNP遗传图谱定位油酸、亚麻酸及芥酸含量QTL位点. 中国农业科学, 2014, 47(1): 24-32.

LIU L Z, LI J N. QTL Mapping of oleic acid, linolenic acid and erucic acid content inby using the high density SNP genetic map., 2014, 47(1): 24-32. (in Chinese)

[31] VAN OOIJEN J W. JoinMap®4: Software for the calculation of genetic linkage maps in experimental populations. Netherlands:Wagningen University, 2006. hops://www.kyazma.nl/index.php/JoinMap/.

[32] KOSAMBI D D. The estimation of map distances from recombination values., 1944, 12: 172-175.

[33] LI H H, YE G Y, WANG J K. A modified algorithm for the improvement of composite interval mapping., 2007, 175(1): 361-374.

[34] BOTELLA C, SAUTRON E, BOUDIERE L, MICHAUD M, DUBOTS E, YAMARYO-BOTTE Y, ALBRIEUX C, MARECHAL E, BLOCK M A, JOUHET J. ala10, a phospholipid flippase, controls fad2/fad3 desaturation of phosphatidylcholine in the ER and affects chloroplast lipid composition in., 2016, 170(3): 1300-1314.

[35] QIU D, MORGAN J, SHI J, LONG Y, LIU J, LI R, ZHUANG X, WANG Y, TAN X, DIETRICH E, WEIHMANN T, EVERETT C, VANSTRAELEN S, BECKETT P, FRASER F, TRICK M, BARNES S, WILMER J, SCHMIDT R, LI J, LI D, MENG J, BANCROFT I. A comparative linkage map of oilseed rape and its use for QTL analysis of seed oil and erucic acid content., 2006, 114(1): 67-80.

[36] 王欣娜, 阎星颖, 李加纳, 刘列钊. 油菜种子油脂分子标记及QTL定位的研究进展. 中国农学通报, 2012, 28(24): 1-7.

WANG X N, YAN X Y, LI J N, LIU L Z. The research progress of lipid molecular markers and QTL mapping in rapeseed., 2012, 28(24): 1-7. (in Chinese)

[37] 张宏军, 肖钢, 谭太龙, 李栒, 官春云. EMS处理甘蓝型油菜()获得高油酸材料. 中国农业科学, 2008, 41(12): 4016-4022.

ZHANG H J, XIAO G, TAN T L, LI X, GUAN C Y. High oleate material of rapeseed () produced by EMS treatment., 2008, 41(12): 4016-4022. (in Chinese)

[38] 陈伟, 范楚川, 钦洁, 郭振华, 傅廷栋, 周永明. 分子标记辅助选择改良甘蓝型油菜种子油酸和亚麻酸含量. 分子植物育种, 2011, 9(2): 190-197.

CHEN W, FAN C C, QIN J, GUO Z H, FU T D, ZHOU Y M. Genetic improvement of oleic and linolenic acid content through marker- assisted selection inseeds., 2011, 9(2): 190-197. (in Chinese)

[39] OKULEY J, LIGHTNER J, FELDMANN K, YADAV N, LARK E, BROWSE J.gene encodes the enzyme that is essential for polyunsaturated lipid synthesis., 1994, 6(1): 147-58.

[40] LI F, CHEN B, XU K, WU J, SONG W, BANCROFT I, HARPER A L, TRICK M, LIU S, GAO G. Genome-wide association study dissects the genetic architecture of seed weight and seed quality in rapeseed (L.)., 2014, 21(4): 355-367.

[41] WU G, WU Y, XIAO L, LI X, LU C. Zero erucic acid trait of rapeseed (L.) results from a deletion of four base pairs in the fatty acid elongase 1 gene., 2008, 116(4): 491-499.

[42] WANG N, WANG Y J, TIAN F, KING G J, ZHANG C Y, LONG Y, SHI L, MENG J L. A functional genomics resource for: development of an EMS mutagenized population and discovery of FAE1 point mutations by TILLING., 2008, 180(4): 751-765.

[43] TRESCH S, HEILMANN M, CHRISTIANSEN N, LOOSER R, GROSSMANN K. Inhibition of saturated very-long-chain fatty acid biosynthesis by mefluidide and perfluidone, selective inhibitors of 3-ketoacyl-CoA synthases., 2012, 76: 162-171.

[44] ROWLAND O, ZHENG H, HEPWORTH S R, LAM P, JETTER R, KUNST L. FAR3 encodes an alcohol-forming fatty acyl-coenzyme A reductase involved in cuticular wS production in., 2006, 142(3): 866-877.

[45] NIKOLAU B J, OHLROGGE J B, WURTELE E S. Plant biotin- containing carboxylases., 2003, 414(2): 211-222.

[46] ZHAO L F, KATAVIC V, LI F L, HAUGHN G W, KUNST L. Insertional mutant analysis reveals that long-chain acyl-CoA synthetase 1 (LACS1), but not LACS8, functionally overlaps with LACS9 inseed oil biosynthesis., 2010, 64(6): 1048-1058.

QTL Identification for Fatty Acid Content inUsing the High Density SNP Genetic Map

YE Sang, CUI Cui, GAO HuanHuan, LEI Wei, WANG LiuYan, WANG RuiLi, CHEN LiuYi, QU CunMin, TANG ZhangLin, LI JiaNa, ZHOU QingYuan

(College of Agronomy and Biotechnology, Southwest University, Chongqing 400715)

【】Rapeseed oil is widely used in cooking, food processing and industrial production. Therefore, improving the fatty acid composition of rapeseed oil according to the specific production objective is an important goal of rapeseed breeding. In this study, QTL mapping of main fatty acid composition inunder two environments was conducted, which was designed to search for the QTL and related candidate genes of fatty acid components in.【】The high generation recombinant inbred lines (RILs) were used as experiment materials, which were derived from synthetic 10D130 and conventional variety Zhongshuang11, and field experiments were conducted in 2016-2017 and 2017-2018 with two different environments in Beibei District of Chongqing City, respectively. After self-pollinated seeds were harvested each year, fatty acid components of seeds were measured by GC with three technical repeats. Then, the RIL population genotype was analyzed with the rapeseed 6K SNP chip array. The DNA preparation and the chip preparation were processed strictly according to Infinium HD Assay Ultra manual of Illumina Inc. The SNP linkage map was constructed by using JoinMap 4.0 program with minimum LOD 2.0. QTL mapping of main fatty acid composition was conducted by composite interval mapping using software Windows QTL IciMapping V4.1.【】In the two environments, the differences of parents’ traits and RILs population’s traits reached significant or extremely significant levels, and the contents of six fatty acids showed continuous distribution, which were suitable for the detection of QTLs. The reference SNP genetic map contains 1 897 polymorphic SNP markers, covering 3 214.19 cM ofgenome with an average map distance of 1.69 cM. Twenty-three QTLs loci of fatty acid components on 8 chromosomes were detected in two environments. The QTLs loci related to stearic acid, oleic acid, linoleic acid, linolenic acid, eicosapentaenoic acid and erucic acid contents were 6, 3, 4, 5, 2 and 3, respectively. And the “enrichment regions” of multiple fatty acid contents were found on chromosomes A05, A08 and C03. Meantime, the main QTL with overlapping linoleic acid and linolenic acid content was detected on A05 chromosome, showing the same additive effect. And the main effect QTLs, oleic acid with overlapping contents of eicosenic acid and erucic acid were detected on A08 and C03, which was opposite to the additive effect of eicosenic acid and erucic acid. Twenty-two candidate genes of fatty acid metabolic were found underlying confidence intervals of seventeen QTLs by comparing with homologous genes in. These genes regulate lipid biosynthesis and metabolism through encoding fatty acid desaturase, total carboxylase synthase, carbon chain lengthening enzyme and participating in acyl coenzyme A biosynthesis. 【】The fatty acid composition QTL under two environments were mapped accurately with the 6K SNP chip of rapeseed, and the “enrichment regions” of multiple fatty acid QTLs on chromosomes A05, A08 and C03 were screened. Compared withfatty acid metabolic genes, candidate genes for fatty acid metabolism in this population were detected, which could be used for improving fatty acid composition in rapeseed.

; fatty acid; genetic map; QTL; candidate gene

2019-05-28；

2019-07-11

10.3864/j.issn.0578-1752.2019.21.002

国家重点研发计划（2018YFD0100500）、农业部现代农业产业技术体系（CARS-12）、重庆市技术创新与应用发展（cstc2019jscx-msxmX0383）

叶桑，Tel：13068300612；E-mail：837276825@qq.com。崔翠，Tel：13883787860；E-mail：cuigreeny@163.com。叶桑和崔翠为同等贡献作者。通信作者周清元，Tel：13883388890；E-mail：zhouqy2005@163.com

（责任编辑 李莉）