本发明涉及一种与瓦氏黄颡鱼生长性状相关的snp分子标记及应用，属于水产动物遗传及分子标记辅助选择育种技术领域。
背景技术：
瓦氏黄颡鱼(pelteobagrusvachelli)是黄颡鱼属鱼类中个体最大的一个物种(体质量达1800g)，与黄颡鱼(p.fulvidraco)并为该属最具经济价值的两个物种，目前养殖热度很高的水产新品系“黄优1号”黄颡鱼是采用瓦氏黄颡鱼作为父本、黄颡鱼作为母本杂交得来，具生长快、体色诱人、抗病力强等优势。在水产新品种培育中，重中之重的环节便是亲本选育，而目前水产新品系“黄优1号”黄颡鱼中作为父本的瓦氏黄颡鱼，相比母本黄颡鱼，在种质改良方面的研究明显滞后，鲜有报道。在瓦氏黄颡鱼繁养殖规模迅速扩大的同时，存在着诸多育种问题，如经济性状退化：原有的捕捞野生亲本进行繁殖制种方式导致野生资源及种质衰退严重，小型化趋势加剧，或多代近亲遗传方式进行制苗，生产的苗种质量良莠不齐，如生长缓慢、头长/体长增大、商品鱼颜色不纯。当务之急需对瓦氏黄颡鱼进行遗传改良，最有效的方法之一就是开发dna分子标记，找到与生长性状连锁的标记，从而为该鱼的分子辅助育种提供遗传依据，这也是瓦氏黄颡鱼及“黄优1号”黄颡鱼繁养殖业可持续发展的必备前提。
相比传统的分子标记(ssr、aflp或rapd等)，snp作为第三代分子标记技术的代表，具有多态性好、广泛分布于全基因组的特点。如今，snp分子标记技术已广泛应用于水产动物生长性状相关snp标记开发中，济鱼类半滑舌鳎(cynoglossussemilaevis)是功能基因关联snps位点开发的一则成功案例，如生长激素受体(ghr1)中snp1357g/a位点gg、ga、aa三种不同基因型与体重、性腺重、三碘甲素及ghr1mrna表达量显著相关，并且gg型甲基化水平明显高于另两个基因型；垂体腺苷酸环化酶激活多肽(pacap)中snp151c/g位点，gg个体具有更高的甲基化水平并且生长速度更快；igf1基因外显子中snp884c/a位点多态性对生长性状、血液生化指标和基因表达有显著影响；天然抗性相关巨噬细胞蛋白(nramp)中snp3125a/g位点等位基因(g)频率和基因型(gg)频率与半滑舌鳎抗鳗弧菌能力呈现极显著相关；igf2基因中两个snp位点与生长性状、血液生化和igf2基因表达显著相关。
技术实现要素：
针对上述现有技术，本发明提供了一种与瓦氏黄颡鱼生长性状相关的snp分子标记，该分子标记可应用于该鱼生长性状的分子标记辅助选择育种，能为进行瓦氏黄颡鱼生长性状关联研究提供有用信息，有利于该鱼生长性状的遗传改良。本发明基于瓦氏黄颡鱼图谱qtl定位结果获取生长性状相关的snp标记，运用目前国际上主流的基因分型方法kasp技术(即竞争性等位基因特异性pcr)开发可直接用于瓦氏黄颡鱼分子标记辅助育种的试剂盒，该技术准确率高，灵活性超强，无需昂贵的双色标记探针，最低的dna样本需求，无需全基因组扩增，快速高效，可以高通量检测大量样本的少数几个标记。其顺利实施必将产生创新性的成果，并加快我国瓦氏黄颡鱼及水产新品系“黄优1号”黄颡鱼良种选育的进程，促进我国黄颡鱼属鱼类良种产业的健康发展。
本发明是通过以下技术方案实现的：
一种与瓦氏黄颡鱼生长性状相关的snp分子标记，该分子标记位于瓦氏黄颡鱼14连锁群41cm位置，其核苷酸序列如seqidno:1所示，在该基因序列的第26位存在一个snp位点：碱基为g或c。
seqidno:1(双下划线所示的碱基“s”，即为snp位点，碱基为g或c)：
用于检测上述snp分子标记的特异性引物，包括正向引物a、正向引物b、反向引物，其核苷酸序列分别如seqidno:2～4所示；在三条引物上进行不同的荧光标记。
正向引物a：5′-catggccagaaaaatctctaatcctg-3′(seqidno:2)。
正向引物b：5′-catggccagaaaaatctctaatcctc-3′(seqidno:3)。
反向引物：5′-gtctcccttaaatgcaagagtgatttcta-3′(seqidno:4)。
优选的，三种荧光标记分别为fam、vic、com。
用于检测上述snp分子标记的试剂盒，包括上述特异性引物，以及检测分子标记所必需的常规试剂。
本发明的发明人经过实验研究发现，瓦氏黄颡鱼14连锁群41cm位置的基因26bp处的snp位点的基因型是gg、gc或cc，与瓦氏黄颡鱼的体重、体宽、全长等生长性状息息相关：cc基因型个体，其体重、体宽、全长这三个生长性状明显优于(所谓优于，是指体重更重，体宽更宽，全长更长)gc基因型个体，gc基因型个体，其体重、体宽、全长这三个生长性状明显优于gg基因型个体。
鉴于上述实验研究成果，上述snp分子标记，可用于检测或辅助检测瓦氏黄颡鱼的生长性状，包括体重、体宽、全长；可用于瓦氏黄颡鱼的辅助选择育种(在以生长性状为选育指标的瓦氏黄颡鱼遗传育种研究过程中，可以优先选择瓦氏黄颡鱼14连锁群41cm位置的基因的26bp处的snp位点的基因型为cc的个体为育种亲本，这对于瓦氏黄颡鱼优良生长性状新品系的选育具有重要的指导意义)。
一种检测瓦氏黄颡鱼生长性状的方法(检测上述snp分子标记的方法)：提取待测瓦氏黄颡鱼的基因组dna(利用少量尾鳍即可提取)，pcr扩增(可利用上述特异性引物或试剂盒进行pcr扩增)，检测瓦氏黄颡鱼14连锁群41cm位置的基因的26bp处的snp位点的基因型是cc、gg还是gc(将所得扩增产物进行荧光信号扫描，根据荧光信号判断瓦氏黄颡鱼14连锁群41cm位置的基因26bp处的snp位点的基因型是cc、gg还是gc)。
进一步地，所述pcr扩增的反应体系为：dna3.5μl，特异性引物混合物0.1μl，2×kaspmastermix3.5μl；特异性引物混合物包括：正向引物a、正向引物b、反向引物和纯化水，正向引物a浓度为11.5μmol/ml，正向引物b浓度为11.5μmol/ml，反向引物浓度为28.25μmol/ml。
进一步地，所述pcr扩增的反应程序为：95℃预变性15min；94℃变性28s，62～57℃梯度pcr，每循环降低-0.5℃，共10个循环；复性28s，95℃变性18s，56℃复性及扩增68s，共28个循环；10℃保存。
本发明的与瓦氏黄颡鱼生长性状相关的snp分子标记，可用于快速鉴定或辅助鉴定优良生长性状(体重、体宽、全长)，不受瓦氏黄颡鱼的年龄、性别等限制，可用于瓦氏黄颡鱼的早期选育，鉴定结果准确可靠，操作简单，成本低，可显著促进瓦氏黄颡鱼的育种进程，解决亲本选育中的盲目性。仅仅使用少量尾鳍即可进行优良生长性状的检测，为分子辅助育种提供科学工具。
具体实施方式
下面结合实施例对本发明作进一步的说明。然而，本发明的范围并不限于下述实施例。本领域的专业人员能够理解，在不背离本发明的精神和范围的前提下，可以对本发明进行各种变化和修饰。
本发明对试验中所使用到的材料以及试验方法进行一般性和/或具体的描述。虽然为实现本发明目的所使用的许多材料和操作方法是本领域公知的，但是本发明仍然在此作尽可能详细描述。
下述实施例中所涉及的仪器、试剂、材料等，若无特别说明，均为现有技术中已有的常规仪器、试剂、材料等，可通过正规商业途径获得。下述实施例中所涉及的实验方法，检测方法等，若无特别说明，均为现有技术中已有的常规实验方法，检测方法等。
下述实施例中的所有引物均由南京集思慧远生物科技有限公司合成。本发明实验所用瓦氏黄颡鱼材料均来自南京师范大学海洋科学与工程学院鱼类遗传育种实验室。
实施例1
本发明的snp分子标记来源于本实验室瓦氏黄颡鱼构建f1子代遗传图谱，后经过生长性状(体重、体宽、全长)qtl定位得来：
(1)挖掘snps分子标记并构建高密度遗传连锁图谱，采用来自不同地理群体、性状分离的亲本，建立相同条件下养殖>3个月的瓦氏黄颡鱼f1代家系，运用微卫星标记筛选出杂合度较高的家系，随机挑选200条全同胞f1代及父母本作为作图群体，记录该群体生长性状(体重、体宽、全长等)数据，然后剪取尾鳍提取dna，基于ddrad-seq技术挖掘snps分子标记，并构建瓦氏黄颡鱼高密度遗传连锁图谱(joinmap)。
(2)生长性状(体重、体宽、全长)进行qtl定位，基于高密度遗传连锁图谱对200条瓦氏黄颡鱼f1代生长性状进行qtl定位(wincartqtl6.0)，发现部分lod值>4与体重、体宽、全长相关的qtl位点，从中筛选出能够有效区分性状、可直接用于瓦氏黄颡鱼分子辅助选择育种的snp标记。
实施例2检测样本的snp分子标记
自2017年7月下旬出膜后，在同一环境(温度、密度等)下养殖105天的450尾瓦氏黄颡鱼个体；将该450尾鱼置于水池(580l)适应2个星期，进行实验前停止投喂饲料2天，使用游标卡尺测量其体宽、全长，使用电子天平测量其体重；随机挑选193条瓦氏黄颡鱼，剪取尾鳍于95％乙醇-20℃保存，用于基因组dna提取，具体步骤如下：
(1)取少量尾鳍，提取瓦氏黄颡鱼的基因组dna；
(2)采用特异性引物组进行pcr扩增，将所得扩增产物进行荧光信号扫描；
(3)根据所述荧光信号判断瓦氏黄颡鱼14连锁群41cm位置26bp处的snp位点的基因型是cc还是gc还是gg。
所述特异性引物组，由fam正向引物、vic正向引物、com反向引物组成，核苷酸序列分别如下所示(如seqidno:2～4所示)。
fam正向引物：5′-catggccagaaaaatctctaatcctg-3′。
vic正向引物：5′-catggccagaaaaatctctaatcctc-3′。
com反向引物：5′-gtctcccttaaatgcaagagtgatttcta-3′。
所述pcr扩增的反应体系为：dna3.5μl，kasp引物混合物0.1μl，2×kaspmastermix3.5μl；特异性引物混合物包括：fam正向引物、vic正向引物、com反向引物和纯化水；fma正向引物浓度为11.5μmol/ml，vic正向引物浓度为11.5μmol/ml，com反向引物浓度为28.25μmol/ml。
所述pcr扩增的反应程序为：95℃预变性15min；94℃变性28s，62～57℃梯度pcr，每循环降低-0.5℃，共10个循环；复性28s，95℃变性18s，56℃复性及扩增68s，共28个循环；10℃保存。
将利用上述方法检测到的瓦氏黄颡鱼基因型，与前期记录得到的体重、体宽、全长进行kruskal-wallistest，得出如下结果(如表1、表2所示)，瓦氏黄颡鱼体重、体宽、全长的kruskal-wallis检验统计量分别为7.336、9.463和6.109，该结果显示：体重、体宽、全长与瓦氏黄颡鱼14连锁群41cm位置的基因的26bp处的snp位点的基因型是cc、gc还是gg显著相关，根据待测瓦氏黄颡鱼的基因型，并经过进一步平均数据统计得到待测三种基因型的瓦氏黄颡鱼体重、体宽、全长关系为：cc>gc>gg。
表1同一环境下养殖的193尾瓦氏黄颡鱼体重、体宽、全长统计数据
表2瓦氏黄颡鱼基因型结合体重、体宽、全长数据的kruskal-wallistest结果
k*：thekruskal-wallis检验统计量；significancelevels:**:0.05，***:0.01。
给本领域技术人员提供上述实施例，以完全公开和描述如何实施和使用所主张的实施方案，而不是用于限制本文公开的范围。对于本领域技术人员而言显而易见的修饰将在所附权利要求的范围内。本说明书引述的所有出版物、专利和专利申请通过引用并入本文，如同这些出版物、专利和专利申请各自特别地和个别地表明通过引用并入本文。
序列表
<110>菏泽学院
<120>一种与瓦氏黄颡鱼生长性状相关的snp分子标记及应用
<141>2018-10-30
<160>4
<170>siposequencelisting1.0
<210>1
<211>150
<212>dna
<213>pelteobagrusvachelli
<400>1
catggccagaaaaatctctaatcctsagcatctgtctcattctctcttcactgaggaaca60
attagaaatcactcttgcatttaagggagactggtgtgaaaattaaacaggatgtgaaaa120
tctataacccttactatctcacttgattta150
<210>2
<211>26
<212>dna
<213>artificialsequence
<400>2
catggccagaaaaatctctaatcctg26
<210>3
<211>26
<212>dna
<213>artificialsequence
<400>3
catggccagaaaaatctctaatcctc26
<210>4
<211>29
<212>dna
<213>artificialsequence
<400>4
gtctcccttaaatgcaagagtgatttcta29