CN114592073B

CN114592073B - 影响鸭脂肪酸组成相关基因elovl3的snp

Info

Publication number: CN114592073B
Application number: CN202210336713.XA
Authority: CN
Inventors: 侯卓成; 李晓琴
Original assignee: China Agricultural University
Current assignee: China Agricultural University
Priority date: 2022-04-01
Filing date: 2022-04-01
Publication date: 2023-03-14
Anticipated expiration: 2042-04-01
Also published as: ZA202212446B; WO2022233345A9; CN114592073A; WO2022233345A3; WO2022233345A2

Abstract

本申请提供了影响鸭脂肪酸组成相关基因ELOVL3的SNP及其在选育肉鸭、鉴定肉质等方面的应用，本申请中的SNP位于ELOVL3基因上游‑619bp处，基因型为AA、GG或GA。本发明结合GWAS分析及细胞试验鉴别到了位于ELOVL3基因上游区域影响其表达活性的SNP，进一步检测确定AA型个体具有比GG型个体更高的多不饱和脂肪酸碳链延长能力，选择对提高鸭肉脂肪酸含量有利的优势基因型进行留种，可迅速提高鸭肉长链多不饱和脂肪酸含量以及肉品质，将加快优质鸭育种改良的进程，为鸭的产业发展带来极大的推动作用和客观的经济效益。

Description

影响鸭脂肪酸组成相关基因ELOVL3的SNP

技术领域

本申请属于分子生物学领域和家禽选育领域。具体地，本申请提供了影响鸭脂肪酸组成相关基因ELOVL3的SNP及其应用。

背景技术

肉禽是我国畜禽生产的重要组成部分，具有易饲养，生产性能高的特点，是世界第二大消费肉类。随着现代社会对健康和多样化食品的需求，家禽产品因其营养均衡、安全保健和风味独特等优势，市场占有率逐步上升，在生活水平提高、生活方式变迁等因素的推动下，家禽产品消费结构不断升级，消费者对家禽产品质量安全和营养价值更加青睐。

单核苷酸多态性(SNP)由基因组水平上由单个核苷酸变异引起，占所有已知多态性的90％以上，平均每300个碱基对中就有1个。全基因组关联分析(GWAS)基于SNP间的连锁不平衡来鉴别影响表型和基因之间关系，可以有效挖掘与主选性状相关的分子标记，并应用于分子标记辅助选择(MAS)、全基因组选择(GS)来实现目标性状的早期选择。相比传统育种方法，MAS和GS选育使性状相关等位基因快速纯和，不但加快了选育效率和遗传选择进展，也可在配套系创制中避免商品代性状分离，更高效地选育出性状优良的畜禽品种。

不饱和脂肪酸(unsaturated fatty acid)是碳链中含有双键的一类脂肪酸，因含有双键的数量不同，分为单不饱和脂肪酸(monounsaturated fatty acid,MUFA)和多不饱和脂肪酸(ployunsaturated fatty acid,PUFA)。饱和脂肪酸(SFA)含量高不易被身体吸收，研究发现，SFA长期摄入量过多会导致高血脂症，形成动脉粥样硬化，增加心脑血管疾病的发生风险。相反，临床研究证明，多不饱和脂肪酸(PUFA)有利于人体健康，因其能调节人体脂质代谢，减少血液中胆固醇和甘油三酯含量，从而降低血液黏稠度并改善血液微循环，对心脑血管疾病的治疗和预防有重要价值；同时还被证明，其能增强人体免疫防御功能以及排除过量SFA形成的多余脂肪，具有减肥的功效。长链多不饱和脂肪酸(long chainployunsaturated fatty acid,LCPUFA)是指链长为18-22个碳原子的PUFA，LCPUFA广泛存在于生物体内，其中最受关注且功能研究比较清楚的3个成员分别是AA、EPA和DHA。LCPUFA对人体具有重要的生理功能，具有软化血管、健脑益智、改善视力等功效。其中，n-3系列PUFAs是发挥这些功能的主要成分，其对健康的积极影响主要通过调节和抑制类二十烷酸的酸合成途径，进而改变炎症反应和相关分子蛋白表达活性、调节与正常及病理细胞功能有关的各种信号通路中分子及酶活性，且n-3脂肪酸可并入膜磷脂直接影响基因的表达，所有这些通路均高度相关，n-3脂肪酸对健康和疾病的生物有益作用是靠多种协调机制完成的。

然而，对于PUFA，人体自身大多不能合成，需要从外界摄取。由此，多不饱和脂肪酸在食品、健康辅料、化妆品和制药中的潜在市场是十分巨大的。目前，不饱和脂肪酸仅能从高等植物种子和深海鱼油中获取。但是受限于其易氧化的性质和工艺复杂的提炼过程，无法满足市场及广大人民需求。因此，长期以来，人们在不断探求可替代的生物资源。因禽肉具有普遍性，禽类养殖产量庞大，且已有研究表明，禽肉中的PUFAs高于畜肉，可作为人类多不饱和脂肪酸摄入的重要来源。禽类研究发现，鸭肉中PUFAs含量高于鸡肉，而在表达数据中，脂肪酸延伸家族基因中ELOVL3的表达鸭显著高于鸡。且ELOVL3在白色脂肪组织中具有显著表达，该基因具有调节饱和超长链脂肪酸和甘油三酯的内源性合成等作用。目前，国内外对其它家畜脂肪酸组成的研究较多，但有关禽肉的脂肪酸组成及合成机制报道还较少。启动子区域决定转录起始，调控基因表达的强弱，研究基因的启动子转录调控具有非常重要的作用，关于ELOVL3的启动子结构尚未解析，所以通过对鸭ELOVL3基因启动子活性分析，明确该基因启动子调控区，同时预测核心启动子区的候选转录因子结合位点，并结合全基因组关联分析及选择信号分析结果，找到ELOVL3调控区的候选因果突变，通过实验加以验证，可为后续高多不饱和脂肪酸比例禽类新品种的选育工作提供理论支撑和基础数据，对提高其食用价值和推动相关产业发展具有重要意义。

发明内容

为解决上述问题，一方面，本申请提供了影响鸭脂肪酸组成相关基因ELOVL3的SNP，所述SNP位于ELOVL3基因上游-619bp处。

进一步地，所述SNP处核苷酸为AA、GG或GA。

另一方面，本申请提供了上述SNP在判断鸭中ELOVL3表达水平中的应用。

另一方面，本申请提供了上述SNP在判断鸭肉长链多不饱和脂肪酸含量中的应用。

另一方面，本申请提供了使用上述SNP在选育鸭肉中长链多不饱和脂肪酸含量高的鸭中的应用。

进一步地，所述方法包括检测ELOVL3基因上游-619bp处的核苷酸。

进一步地，ELOVL3基因上游-619bp处的核苷酸为AA基因型时个体有更高的中长链多不饱和脂肪碳链延长能力，其肉中有更高的中长链多不饱和脂肪含量。

进一步地，所述应用包括使用ELOVL3基因上游-619bp处的核苷酸为AA基因型的个体选育品种。

进一步地，检测ELOVL3基因上游-619bp处的核苷酸使用的方法选自，测序法、限制性片段长度多态性分析、熔解曲线分析、单链构象异构多态性分析、探针扩增阻滞突变、色谱方法。

另一方面，本发明提供了扩增鸭多不饱和脂肪酸组成相关基因ELOVL3编码区引物及完整cDNA序列的RACE引物对，其核苷酸序列为：

CDS区域：SEQ ID NO.1和2；

cDNA RACE：SEQ ID NO.3和4。

另一方面，本发明提供了扩增鸭多不饱和脂肪酸组成相关基因ELOVL3上游调控区不同长度的截断引物对，其核苷酸序列为：SEQ ID NO.5-12。

另一方面，本发明提供了ELOVL3基因的核心启动子区域，其位于ELOVL3基因上游区域为-765bp～+62bp处。

另一方面，本申请还提供了用于进行上述应用或者包含上述引物的试剂盒。

本申请中检测多态性的方法不限于上面所列的测序法、限制性片段长度多态性分析、熔解曲线分析、单链构象异构多态性分析、探针扩增阻滞突变、色谱方法，其他已知或研究中的方法也可以应用于本发明。

根据所用方法不同，本申请的试剂盒中可以包含各种检测所需的试剂，包括但不限于引物、探针、缓冲液、酶等。

有益效果：ELOVL3是动物体内控制长链脂肪酸碳链延伸的关键酶。本发明结合GWAS分析及细胞试验鉴别到了位于ELOVL3基因上游区域影响其表达活性的SNP，进一步检测个体表型效应后将其用于MAS中，选择对提高鸭肉脂肪酸含量有利的优势基因型进行留种，可迅速提高鸭肉长链多不饱和脂肪酸含量以及肉品质，将加快优质鸭育种改良的进程，为鸭的产业发展带来极大的推动作用和客观的经济效益。本发明首次得到与多不饱和脂肪酸组成相关基因ELOVL3表达活性相关的一个SNP分子标记，为该基因上游-619位点的A<G突变,表型优势基因型AA的个体，为早期选择提供试验数据依据，从而大大地提高新品系选育效率，为优质鸭的标记辅助选择提供了新的方法。

附图说明

图1A和图1B为北京鸭肌肉脂肪含量的全基因组关联分析(GWAS)曼哈顿图及该性状在7号染色体上关联分析的曼哈顿图。横坐标表示鸭的染色体编号；纵坐标表示SNP位点的-logP值；

图2为鸭ELOVL3基因完整的cDNA结构示意图，包括三个部分：5’UTR(5’非编码序列)；CDs(编码序列)；3’UTR(3’非编码序列)；

图3为北京鸭ELOVL3截断片段电泳结果；

图4为鸭ELOVL3启动子区域截断片段双荧光酶相对活性测定结果，横坐标为相对荧光活性值；纵坐标为启动子片段截断区域；

图5为鸭ELOVL3基因核心启动子区域-619位点(A<G)单点突变序列测序情况；

图6为鸭ELOVL3核心启动子区域-619位点不同碱基情况下的双荧光酶相对活性测定结果。横坐标为相对荧光活性值；纵坐标为核心启动子区域-619位点的不同碱基情况。

具体实施方式

为了使本发明的目的、技术方案更加清楚明白，以下结合实施例进一步说明本发明的内容，但不应理解为对本发明的限制。在不背离本发明精神和实质的情况下，对本发明方法、步骤或条件所作的修改或替换，均属于本发明的范围。若未特别指明，实施例所用技术手段为本领域技术人员所熟知的常规手段。

实施例1鸭多不饱和脂肪酸组成相关基因ELOVL3完整cDNA序列的扩增

根据Ensembl(http://asia.ensembl.org/)数据库中获取鸭7号染色体ELOVL3(ENSAPLG00020013672)的cDNA序列，分别在5’UTR和3’UTR区域设计扩增引物。

引物序列为：5’-GGACAGGCCAGAGGTCACT-3’(正向，SEQ ID NO.1)，5’-GAACAGGCGGCATGGTCACT-3’(反向，SEQ ID NO.2)，以野鸭肝脏RNA反转录获得的cDNA为模板扩增包含该基因完整蛋白编码区域(CDs)的序列。

PCR：20μl反应体系：2x

PCR SuperMix 10μl，上下游引物(10μmol/L))各0.5μl，DNA模板1μl，超纯水补8μl。

PCR扩增程序：94℃预变性5min；94℃变性30s，60℃退火30s，72℃延伸1min，共35个循环；72℃后延长10min，最后4℃保存。

以扩增得到的CDs序列为模板，首先分别设计5’RACE特异性引物GPS1和3’RACE特异性引物GPS2，引物前端根据线性化pUC19载体序列分别添加15bp同源序列GATTACGCCAAGCTT，引物序列为：

5’-GATTACGCCAAGCTTGCCGCCGCCAGGCACCATCTCCTTGTA-3’(GPS1，SEQ ID NO.3)；

5’-GATTACGCCAAGCTTCGCGTCTTCCATGGGATTCAAGCAGTCGG-3’(GPS2，SEQ ID NO.4)。

按照5’/3’RACE试剂盒(

RACE 5’/3’Kit,Takara bio,USA)操作指南分别进行5’和3’RACE反应。

获得目的长度片段后利用纯化试剂盒(NucleoSpin Gel and PCR Clean-Up Kit,Takara bio,USA)对目的片段进行纯化。

通过使用无缝克隆试剂盒(

HD Cloning Kit,Takara bio,USA)同源重组连接至线性化的pUC19载体上，测序获得序列信息，拼接获得ELOVL3基因的完整cDNA结构，共包含1724碱基对，其中，5’UTR区含有143bp；CDs区810bp；3’UTR区771bp，如图2。以上引物合成及测序均由北京擎科生物技术有限公司完成。

实施例2鸭多不饱和脂肪酸组成相关基因ELOVL3的核心启动子区域确认

从Ensembl(http://asia.ensembl.org/)数据库获取鸭7号染色体ELOVL3(ENSAPLG00020013672)上游2500bp序列，利用在线网站Primer3(https://bioinfo.ut.ee/primer3-0.4.0/)进行不同长度片段引物设计，并添加同源重组臂及酶切位点，扩增不同长度的缺失片段，截断引物序列分别为：

上游引物(5’-3’)：

ELOVL3-2095-F1:

ggcctaactggccggtaGGTACCAGTGAAGAGCTCCACGCACT(SEQ ID NO.5)

ELOVL3-1856-F2:

ggcctaactggccggtaGGTACCTGGGAAACAAGTCATGTCCA(SEQ ID NO.6)

ELOVL3-1518-F3:

ggcctaactggccggtaGGTACCCAGGGCAATTCCTCTAGCAT(SEQ ID NO.7)

ELOVL3-1249-F4:

ggcctaactggccggtaGGTACCCCCTGGGACAGCTCACAGTA(SEQ ID NO.8)

ELOVL3-765-F5:ggcctaactggccggtaGGTACCCACATCCTTCCATGCCACAT(SEQ IDNO.9)

ELOVL3-572-F6:ggcctaactggccggtaGGTACCAAAGATGATGCCCAAATTGC(SEQ IDNO.10)

ELOVL3-389-F7:ggcctaactggccggtaGGTACCGATAACTGGGTGCAGGTTCC(SEQ IDNO.11)

下游引物(5’-3’)：

ELOVL3-+62-R:cgaggccagatcttgatatcCTCGAGCTCCGCTCCACCTCATACTC(SEQ IDNO.12)

利用基因组提取试剂盒(TIANamp Genomic DNA Kit,Tiangen bio,CHINA)提取北京鸭血液基因组作为模板，使用PCR mix(2x

PCR SuperMix,Transgen bio,CHINA)进行PCR扩增，

PCR：20μl反应体系：2x

取PCR产物5μl进行凝胶电泳，检测产物片段大小符合目的片段大小，鉴定结果如图3，并将PCR产物送北京擎科生物技术有限公司测序，测序结果确认PCR产物为目标序列则片段扩增成功。

片段扩增成功后利用纯化试剂盒(TaKaRa MiniBEST AgaroseGel DNAExtraction Kit Ver.4.0,Takara bio,USA)纯化目的片段。

利用软件SnapGene 5.3.1分析pGL4.10限制性内切酶位点，选择利用Kpn I和XhoI(NEB,US)对其进行双酶切并纯化(TaKaRa MiniBEST DNAFragment Purification KitVer.4.0,Takara bio,USA)线性化载体。

利用同源重组试剂(

Seamless Cloning and Assembly Kit,Transgen bio,CHINA)将各片段构建至线性化的pGL4.10报告载体上。

鸡永生化细胞系ICP在37℃、5％CO2浓度和95％空气湿度条件下进行培养，培养基为含10％胎牛血清和1％双抗的DMEM/F12，细胞状态良好时进行细胞计数铺于24孔板中，铺板密度为在24h能生长至80％。

铺板24h后，细胞贴壁覆盖度达80％左右时进行转染，每孔表达载体(0.475μg)、pRL-TK(0.025μg)、脂质体(1.0μL)和Opti-MEM共50μl混匀孵育10min，转染到ICP1细胞中，6h后更换新鲜培养基。

转染48h后收集细胞样品，加入100μL 1x的细胞裂解液裂解5min后，离心取20μL上清液，利用多功能酶标仪(Infinite F200,CH)测定并计算萤火虫荧光素酶与海肾荧光素酶的相对活性(Dual Luciferase Reporter Assay Kit,Vazyme bio,CHINA)，试验进行3次系统重复，每次包含3次技术重复。测定结果如图4所示。

上述结果表明，ELOVL3基因启动子核心区域为-765bp～+62bp区域，此结果为后续探究ELOVL3基因表达活性的影响因素提供了重要基础。

实施例3与上述鸭脂肪酸组成相关基因ELOVL3核心启动子活性相关的SNP突变效应验证

获得核心启动子区域候选突变后以原始表达质粒为模板，利用在线网站(https://crm.vazyme.com/cetool/singlepoint.html)引入单碱基突变的扩增引物进行反向扩增，引物序列为：

5’-ACgCAATTCTTGTAAAGAGAGACATGACTATG-3’(正向，SEQ ID NO.13)，5’-CTTTACAAGAATTGcGTAATAATACTTTGCTGAATTCTCCTGC-3’(反向，SEQ ID NO.14)，针对扩增产物进行DpnI消化，去除甲基化模板质粒后进行重组连接反应，最终进行产物转化、涂板及克隆鉴定(Mut

II Fast Mutagenesis Kit V2,Vazyme bio,CHINA)。克隆菌液送北京擎科生物技术有限公司测序，测序结果如图5，确认单点突变构建成功。提取相应表达质粒备用。

转染48h后收集细胞样品，加入100μL 1x的细胞裂解液裂解5min后，离心取20μL上清液，利用多功能酶标仪(Infinite F200,CH)测定并计算萤火虫荧光素酶与海肾荧光素酶的相对活性(Dual Luciferase Reporter Assay Kit,Vazyme bio,CHINA)，试验进行3次系统重复，每次包含3次技术重复。测定结果如图6所示。

上述试验结果表明，当ELOVL3基因上游-619位点为AA基因型时，其核心启动子区域有更高的转录活性，表明AA型个体具有比GG型个体更高的多不饱和脂肪酸碳链延长能力，为其在检测和提升鸭肉长链多不饱和脂肪酸含量的应用提供了重要依据。

SEQUENCE LISTING

<110> 中国农业大学

<120> 影响鸭脂肪酸组成相关基因ELOVL3的SNP

<160> 14

<170> PatentIn version 3.5

<210> 1

<211> 19

<212> DNA

<213> artificial

<400> 1

ggacaggcca gaggtcact 19

<210> 2

<211> 20

<212> DNA

<213> artificial

<400> 2

gaacaggcgg catggtcact 20

<210> 3

<211> 42

<212> DNA

<213> artificial

<400> 3

gattacgcca agcttgccgc cgccaggcac catctccttg ta 42

<210> 4

<211> 44

<212> DNA

<213> artificial

<400> 4

gattacgcca agcttcgcgt cttccatggg attcaagcag tcgg 44

<210> 5

<211> 43

<212> DNA

<213> artificial

<400> 5

ggcctaactg gccggtaggt accagtgaag agctccacgc act 43

<210> 6

<211> 43

<212> DNA

<213> artificial

<400> 6

ggcctaactg gccggtaggt acctgggaaa caagtcatgt cca 43

<210> 7

<211> 43

<212> DNA

<213> artificial

<400> 7

ggcctaactg gccggtaggt acccagggca attcctctag cat 43

<210> 8

<211> 43

<212> DNA

<213> artificial

<400> 8

ggcctaactg gccggtaggt accccctggg acagctcaca gta 43

<210> 9

<211> 43

<212> DNA

<213> artificial

<400> 9

ggcctaactg gccggtaggt acccacatcc ttccatgcca cat 43

<210> 10

<211> 43

<212> DNA

<213> artificial

<400> 10

ggcctaactg gccggtaggt accaaagatg atgcccaaat tgc 43

<210> 11

<211> 43

<212> DNA

<213> artificial

<400> 11

ggcctaactg gccggtaggt accgataact gggtgcaggt tcc 43

<210> 12

<211> 46

<212> DNA

<213> artificial

<400> 12

cgaggccaga tcttgatatc ctcgagctcc gctccacctc atactc 46

<210> 13

<211> 32

<212> DNA

<213> artificial

<400> 13

acgcaattct tgtaaagaga gacatgacta tg 32

<210> 14

<211> 43

<212> DNA

<213> artificial

<400> 14

ctttacaaga attgcgtaat aatactttgc tgaattctcc tgc 43

Claims

1.选育鸭肉中长链多不饱和脂肪酸含量高的鸭的方法，其特征在于，所述方法包括检测SNP处的基因型并选用所述SNP处为AA基因型的个体培育品种；其中所述SNP位于Ensembl中序列ENSAPLG00020013672所示ELOVL3基因上游-619 bp处，所述SNP处基因型为AA、GG或GA。

2.根据权利要求1所述的方法，其中检测所述SNP处的基因型的方法选自测序法、限制性片段长度多态性分析、熔解曲线分析、单链构象异构多态性分析、探针扩增阻滞突变或色谱方法。