一种与大黄鱼抗副溶血弧菌关联的SNP位点、其筛选方法和应用与流程

本发明属于鱼类分子标记筛选技术领域，涉及一种与大黄鱼抗副溶血弧菌关联的单核苷酸多态性(snp)分子标记，特别涉及一种位于大黄鱼ctrp9基因外显子上与抗副溶血弧菌关联的snp位点、其筛选方法和应用。

背景技术：

大黄鱼(larimichthyscrocea)属硬骨鱼纲、鲈形目、石首鱼科、黄鱼属，是我国近海主要养殖经济鱼类，近些年，由于其养殖环境恶化、养殖规模扩大、养殖密度过大等原因，导致其病害频发，尤其是副溶血弧菌(vibrioparahaemolyticus)感染引起的弧菌病，对大黄鱼生长危害极大。弧菌病发病时间短、死亡快、传染性极强，严重威胁着渔业经济发展和人们的健康。因此，培育大黄鱼抗病品系是预防弧菌的重要途径之一。

单核苷酸多态性(singlenucleotidepolymorphism，snp)是指在基因组水平上由单个核苷酸的变异所引起的dna序列多态性，是一种理想的遗传标记。继限制性片段长度多态性(rflp)、微卫星标记(ssr)后，snp作为新一代分子标记，在遗传的多样性分析、关联分析、品种鉴定、高密度遗传连锁图谱的构建以及辅助育种等方面有广泛应用。目前尚未见到有关从免疫相关基因中筛选到与大黄鱼抗副溶血弧菌相关的snp位点，以及应用的报道。

技术实现要素：

本发明针对副溶血弧菌感染大黄鱼患弧菌病，提供一种位于大黄鱼ctrp9基因上与抗副溶血弧菌关联的snp位点，确定大黄鱼免疫相关基因的snp位点与抗副溶血弧菌关联，为大黄鱼抗副溶血弧菌的选育提供分子标记。

本发明中，与大黄鱼抗副溶血弧菌相关的snp位点，位于大黄鱼如seqidno:1所示的ctrp9基因cdna序列上，全长1368bp，编码455个氨基酸，其1147位碱基处存在一个a或g的等位基因突变，且该突变引起第383位氨基酸由异亮氨酸(ile)突变为缬氨酸(val)。

本发明还提供了上述与大黄鱼抗副溶血弧菌相关的snp位点的筛选方法，包括以下步骤：

(ⅰ)用副溶血弧菌感染大黄鱼幼鱼后，提取其rna并反转录成cdna，其核苷酸序列如seqidno:1所示；

(ⅱ)设计引物对如seqidno:2所述的正向引物和如seqidno:3所示的反向引物，以cdna为模板进行pcr扩增反应；

(ⅲ)经琼脂糖凝胶检测，对目的条带清晰无杂带的pcr产物进行测序，筛选得到上述snp位点。

步骤(ⅱ)中，所述pcr扩增反应体系25ul为：2×mix12.5ul，正向引物(10um)和反向引物(10um)各1ul，cdna模板1ul，无菌水9.5ul；和/或

所述pcr扩增反应的程序依次为：95℃预变性5min，94℃变性30s、64℃退火30s、72℃延伸30s(45个循环)，72℃终延伸10min。

本发明还提供了上述与大黄鱼抗副溶血弧菌相关的snp位点在抗副溶血弧菌大黄鱼遗传育种上的用途。

与现有技术相比，本发明的有益效果是：

本发明对大黄鱼免疫相关基因非同义突变snps进行筛选与验证，snps标记不仅可以评估大黄鱼遗传多样性，而且为大黄鱼的抗病选育等工作提供研发基础。同时，本发明检测到大黄鱼ctrp9基因与副溶血弧菌抗性相关的snp位点，为大黄鱼的抗菌选育工作提供新思路，有利于推进大黄鱼的遗传育种进程，提高大黄鱼养殖经济效益。

附图说明

图1是大黄鱼ctrp9基因的snp位点峰型图，自上而下依次为副溶血弧菌易感组和抗副溶血弧菌组。

具体实施方式

下面结合附图和实施例对本发明进一步说明。

以下实例中大黄鱼ctrp9基因上与抗副溶血弧菌关联的snp位点的筛选方法，通过用副溶血弧菌感染大黄鱼幼鱼，提取rna并反转录成cdna，采用引物设计、pcr扩增测序、snp位点筛选获得snp位点，利用spss18.0软件对感染前后抗感易感群体的基因型频率进行差异显著性分析，采用bioedit软件进行分型获得与抗副溶血弧菌相关的snp位点，具体步骤如下：

(1)大黄鱼暂养、驯化

出膜40天、1～2cm的大黄鱼苗由象山港水产引种育种有限公司赠送，放到1.0m(内径)×1.2m(高)的养殖桶里暂养3～4天，水体高1m，温度26℃，每天投饵2次(早上7:00、下午6:30)，换水2次，每次换水90％，清理底部食物残渣及粪便。为了防止大黄鱼患白点病，换水时加入5ml甲醛，浓度为0.005‰。待鱼苗适应环境、进食稳定后，将暂养的大黄鱼移入实验缸中再进行实验。

(2)人工感染

实验水体为长80cm、宽45cm的玻璃缸，加入海水70l，每缸放苗500尾，水温控制在26℃，鱼苗稳定后，进行实验；设置对照组和副溶血弧菌(atcc17802)感染实验组，浓度为1.5×10⁷cfu/ml，实验组和对照组均设置3个平行组。各缸每隔24h投喂一次饵料，清理掉食物残渣和粪便，换水1/3，并补充感染菌以维持其恒定的浓度。每天清理、记录各缸中死苗数，rna保护液保存死亡样本用于后续pcr验证。

(3)转录组分析

取实验第七天的对照组、副溶血弧菌感染患病的、副溶血弧菌感染未患病的样本各三条，trizol法提取样本总rna，混合后质检、建库，用truseqpeclusterkit在cbot中进行clustergeneration，然后在illuminahiseqtm2500中进行双向测序。测序得的rawreads去除低质量和错误碱基后得到cleanreads，以模式生物斑马鱼转录组数据为参考序列进行比对，对测序数据的kegg通路、差异表达基因和snp进行分析。从转录组中筛选免疫相关的kegg通路，根据通路中的免疫相关因子从差异表达基因中筛选差异的免疫基因，并筛选免疫基因的有义突变位点。

从抗副溶血弧菌组与对照组中筛选出127个免疫差异基因，其中有53个免疫基因上调，73个下调基因；53个免疫上调基因中共发现29个基因编码区存在146个snp位点，其中非同义突变位点有28个。

(4)pcr扩增及测序

选择对照组、副溶血弧菌感染第四天死亡(易感)、第七天存活(抗感)的大黄鱼，各10尾、共30尾，分别提取整鱼的总rna。以提取的rna为模版用takara公司的primescript^tmii1standcdnasynthesiskit反转录成cdna，-20℃保存备用。

设计引物对序列如下所示：

primerf：agcctaaccttcctgttcctttctatactatcattt

primerr：ggtggtgttgccatcgggacctt。

采用上述引物以cdna为模板进行pcr扩增，选择诺唯赞公司的2×taqplusmastermix(dyeplus)试剂进行pcr反应，反应体系25ul：2×mix12.5ul，primerf(10um)和primerr(10um)各1ul，cdna模板1ul，无菌水9.5ul。扩增反应器为eppendorfpcr仪(艾本德，德国)。pcr程序：95℃预变性5min，94℃变性30s、64℃退火30s、72℃延伸30s(45个循环)，72℃终延伸10min。pcr产物经琼脂糖凝胶检测、目的条带清晰无杂带，剩余产物送上海生工生物工程股份有限公司进行测序。

(5)基因型的判定及数据分析

根据测序结果用bioedit软件观察测序峰型图，进而确定snps位点的基因型；单峰为纯合子基因型，双峰为杂合子基因型。以上数据用spss18.0软件对对照组、第四天死亡样本(易感)、第七天存活样本(抗感)进行单因素方差分析(p<0.05：显著性差异；p<0.01：极显著差异)。结果发现：ctrp9基因型频率在副溶血弧菌抗感易感群体中存在极显著性差异，其他免疫基因如nf-kappa-binhibitorzeta-like、glutathiones-transferaseomega1、galectin-4、tyrosine-proteinkinasestyk1-like、apobec1complementationfactor、interferonregulatoryfactor2-bindingprotein1-like、arginase基因在副溶血弧菌实验组样本中不存在显著性差异。

(6)扩大样本检测ctrp9基因snp位点

为了更准确验证ctrp9基因snp位点在副溶血弧菌感染组样本中的情况，增加对照组样本至20尾、副溶血弧菌易感组至18尾、抗副溶血弧菌组至16尾，对ctrp9基因的两个snp位点进一步验证。

如图1所示的副溶血弧菌易感组18个样本和抗副溶血弧菌组16个样本一代测序峰型图，snp位点ctrp9-1147-a/g在副溶血弧菌易感组中的基因型频率分别为aa22.2％、ag77.8％、gg0，在抗副溶血弧菌组中aa37.5％、ag25％、gg37.5％。ctrp9基因的ile突变为val位点，副溶血弧菌易感组与抗副溶血弧菌组间存在极显著性差异，但是抗副溶血弧菌组与对照组不存在显著性差异，如表1所示。

表1：ctrp9基因在副溶血弧菌实验组snp位点统计

由表1可知，snp位点ctrp9-1147-a/g在副溶血弧菌抗感易感组基因型频率存在极显著性差异x²＝11.879，p<0.01。表明snp位点ctrp9-1147-a/g在副溶血弧菌抗感易感群体中存在极显著性差异，此位点与大黄鱼抗副溶血弧菌相关联。以上数据分析显示snp位点ctrp9-1147-a/g在副溶血弧菌抗感和易感群体中基因型频率分别为gg37.5％和0，存在极显著性差异(p<0.01)。

snp位点ctrp9-31147-gg基因型可以作为大黄鱼抗副溶血弧菌的候选基因，用于指导大黄鱼的分子标记辅助选择育种。与现有大黄鱼抗菌育种技术相比，从基因角度出发，本发明筛选到与大黄鱼抗副溶血弧菌相关的snp位点，并把此位点运用在大黄鱼育种工作中，有望从遗传角度根本性的解决大黄鱼患弧菌病问题。

序列表

<110>上海海洋大学

<120>一种与大黄鱼抗副溶血弧菌关联的snp位点、其筛选方法和应用

<141>2019-11-18

<160>3

<170>siposequencelisting1.0

<210>1

<211>1368

<212>dna

<213>2ambystomalateralexambystomajeffersonianum

<400>1

atgctaaacaacactgccccagggactggctttggtgtggacaacaccattgattggaga60

ggctcgggctcgaactcgggctcggactcgggctcggactcggactcggactcagatgaa120

agtggggggcacttcccctttaatgcacgtgcccgtcatcaccacatgccggggccaccc180

agcaacatgacccggatgatgatgagagaagatgacatgcccattttacctgatctatca240

atctgtgacatgctgttgagctcacctaacccaccacctatcaacgagattccttttttc300

tgcctctgttcccactgtaagggtaccgtgggacccaaaggagaccgcggagaccgaggt360

ctcccaggtcagaagggcgacgttggggaacaggggatgcctggatctatgggtttcact420

gggatgaagggcgagcgaggcttcaaaggagaaaaaggtgacagaggaacggcaggcccc480

caaggtcaacaaggcccccagggagaaaccggcatatgccctgcttcctgcgatagtgtc540

cagggtccaccaggtgaacaaggcccctctgggccagccggttctcggggcctgcctgga600

atcaagggagaaaggggacccatgggttttaagggtgataaaggtgacctgggcatacct660

ggcaaccctggggtggacggccagaagggtgatcaaggtatgcaggggatatgtgagtgc720

acagatggggaagacggcgctaatggggcagcaggagcaaagggggataaaggggataaa780

ggtgacactggcccccggggtgcacagggtactatggggctaaaaggcaacatgggtgcc840

atgggtatgatggggcagcctgggccgtgctccccaaccatcaaatcagcattctccgcg900

gcaatcaaccagtcgtttccagagcctaaccttcctgttcctttctatactatcatttat960

aataaggaaaaccatttcgatccaagaggcatgtacacagcccctgttaatggcacatac1020

gtcttcaccttcaacatggcagttgctgaaaggactcttaaggttggcctgttccaaaac1080

ttctggcctataatcaaagcaacggaagtggtaggaaacgagcagtcatccgtgagccag1140

actatcatcctccacctcaaaaggtttgaccgggtgtggctgcaggtgaagagcagcacc1200

accaacggcatgtacaccgacagtgagagcagcagcacgttctctgggtggctgcagttc1260

cctgacagctgcgatccgcccctgggcagacagatgtttcctccacaagatgaccgagac1320

tttagctgggaaggtcccgatggcaacaccacccctccaccatactag1368

<210>2

<211>36

<212>dna

<213>2ambystomalateralexambystomajeffersonianum

<400>2

agcctaaccttcctgttcctttctatactatcattt36

<210>3

<211>23

<212>dna

<213>2ambystomalateralexambystomajeffersonianum

<400>3

ggtggtgttgccatcgggacctt23