CN110846435B - 龙须菜鲁龙1号良种的特异性snp标记及其应用 - Google Patents

龙须菜鲁龙1号良种的特异性snp标记及其应用 Download PDF

Info

Publication number
CN110846435B
CN110846435B CN201911291611.5A CN201911291611A CN110846435B CN 110846435 B CN110846435 B CN 110846435B CN 201911291611 A CN201911291611 A CN 201911291611A CN 110846435 B CN110846435 B CN 110846435B
Authority
CN
China
Prior art keywords
sequencing
lulong
genome
centrifuging
dna
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201911291611.5A
Other languages
English (en)
Other versions
CN110846435A (zh
Inventor
隋正红
胡依依
肖保衡
李萍萍
刘源
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Ocean University of China
Original Assignee
Ocean University of China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Ocean University of China filed Critical Ocean University of China
Priority to CN201911291611.5A priority Critical patent/CN110846435B/zh
Publication of CN110846435A publication Critical patent/CN110846435A/zh
Application granted granted Critical
Publication of CN110846435B publication Critical patent/CN110846435B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6876—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
    • C12Q1/6888—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for detection or identification of organisms
    • C12Q1/6895—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for detection or identification of organisms for plants, fungi or algae
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6806—Preparing nucleic acids for analysis, e.g. for polymerase chain reaction [PCR] assay
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2600/00—Oligonucleotides characterized by their use
    • C12Q2600/13—Plant traits
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2600/00—Oligonucleotides characterized by their use
    • C12Q2600/156—Polymorphic or mutational markers
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A40/00—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
    • Y02A40/80—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in fisheries management

Landscapes

  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Organic Chemistry (AREA)
  • Analytical Chemistry (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Engineering & Computer Science (AREA)
  • Zoology (AREA)
  • Wood Science & Technology (AREA)
  • Health & Medical Sciences (AREA)
  • Biotechnology (AREA)
  • Molecular Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Immunology (AREA)
  • Biophysics (AREA)
  • Physics & Mathematics (AREA)
  • Genetics & Genomics (AREA)
  • Biochemistry (AREA)
  • Microbiology (AREA)
  • General Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Botany (AREA)
  • Mycology (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

鲁龙1号龙须菜的特异性SNP标记，具体包括如下步骤：1、基因组DNA提取与质量检测：采用天根植物基因组DNA提取试剂盒进行良种981、2007和鲁龙1号的基因组DNA的提取。2、基因组重测序文库构建与测序。3、测序数据质控：对原始测序数据进行过滤，从而得到高质量的测序数据以保证后续分析的顺利进行。4、比对参考基因组。5、SNP检测与分型。6、特异性SNP标记的筛选。7、特异性SNP标记的引物设计及Sanger测序检验。本发明的有益效果：通过开发鲁龙1号良种的特异性SNP标记，从分子水平上将其与良种981和良种2007区分开来，解决龙须菜在栽培过程中产生的种苗混杂现象，避免对龙须菜栽培产业造成重大损失。

Description

龙须菜鲁龙1号良种的特异性SNP标记及其应用
技术领域
本发明涉及海洋藻类的优良品系鉴定，特别涉及鉴定龙须菜鲁龙1号良种的特异性SNP标记及其应用。
背景技术
龙须菜是一种重要的产琼胶经济海藻。迄今为止，我国成功选育的龙须菜国家级水产品种已有3个，分别为：1998年审定通过的龙须菜981，2014年审定通过的龙须菜2007和2015年审定通过的龙须菜鲁龙1号。其中，鲁龙1号是中国海洋大学海洋生物遗传学与育种教育部重点实验室隋正红教授团队与莆田市水产推广站历时六年培育而成，藻体细长，颜色鲜红，分枝多，生长优势十分明显，并且藻体富有弹性，抗风浪能力强，同时生长速度比其他品种更快，产量更高，经济效益显著。龙须菜新品种的成功选育极大推动了其大规模栽培业的发展。目前龙须菜已经在山东、福建、广东、浙江和辽宁等地广泛栽培，是继海带和紫菜之后的新兴海藻栽培物种。龙须菜栽培成为我国第二大海藻栽培产业。但3个良种在外观形态上难以区分，在大规模栽培过程中容易产生种苗混杂现象。之前通过测定生长性状和生理生化指标来加以区别，但试验周期长，操作繁琐。
分子标记是继形态标记、细胞学标记和生化标记后产生并发展起来的一种新的遗传标记，所反映出的遗传变异种类十分丰富，不易受环境影响，且DNA序列的稳定性使得分子标记具有良好的重复性。分子标记技术可用来进行种质鉴定。单核苷酸多态性(SingleNucleotide Polymorphism,SNP)指由于单个核苷酸的变异所引起的DNA序列多态性。作为最新一代的分子标记，SNP比过去的分子标记更加多态，数目也更多，利于分析。伴随高通量测序技术的发展，利用全基因组测序筛查SNP已成为最直接并且可获得最完整基因组水平单碱基遗传变异信息的方法。
专利号为ZL201310243064.X，专利名称为“鉴定龙须菜981良种的特异性探针及其应用”的发明专利，所提供的技术方案虽然通过AFLP分子标记技术开发出两条981良种的特异性核苷酸序列，但只能用以区分良种981和良种2007，无法区分出良种鲁龙1号。
发明内容
基于上述背景技术，本发明通过在良种981、2007和鲁龙1号中开发最新一代SNP分子标记，筛选获得鲁龙1号良种的特异性SNP标记，从分子水平上将其与良种981和良种2007区分开来，解决龙须菜在栽培过程中产生的种苗混杂现象，避免对龙须菜栽培产业造成重大损失。
为了实现上述目的，本发明采取如下技术方案：
鲁龙1号龙须菜的特异性SNP标记，具体包括如下步骤：
1、基因组DNA提取与质量检测
采用天根植物基因组DNA提取试剂盒进行良种981、2007和鲁龙1号的基因组DNA的提取，具体步骤如下：
(1)缓冲液GP1预热：提取工作开始之前，先用65℃水浴对缓冲液GP1进行预热，并在其中加入0.1％的β-巯基乙醇。而研钵、研杵和药匙则需要利用液氮预冷。
(2)藻体研磨：藻体经无菌处理后，称取约100mg，放入预冷的研砵中，加入液氮充分研磨，使藻体完全破碎成粉末即可。研磨的过程中一定要注意及时添加液氮，并且缓慢加入，防止粉末被快速气化的液氮带出，造成浪费。
(3)藻体裂解：用药匙将藻体粉末转移到(1)中预热的缓冲液GP1中，并迅速颠倒混匀，水浴30～40min进行裂解，水浴过程中每隔5min要进行一次颠倒混匀，保证裂解充分。
(4)离心：裂解结束后，加入700ml氯仿，迅速颠倒混匀，13800g离心5min。
(5)上层透明水相转移：将离心后的上层透明水相转移到新的离心管内，注意吸取时不要混入蛋白层，在新离心管中加入700μl缓冲液GP2，颠倒混匀。
(6)离心：将新离心管内混匀后的液体全部转移到吸附柱CB3中，13800g离心1min，将收集管中的液体弃去。
(7)离心：完成上述操作后，向吸附柱CB3中加入500ml缓冲液GD，13800g离心1min，将收集管中的液体弃去。
(8)离心：向吸附柱CB3中加入600μl漂洗液PW，13800g离心1min，将收集管中的液体弃去。该步骤重复两次。
(9)离心：将吸附柱CB3于13800g离心2min，将收集管中的液体弃去。
(10)吸附柱CB3晾干：将吸附柱CB3置于超净工作台内晾干至无残留酒精气味。需要严格注意吸附柱的管口和底部不能触碰到其他未经灭菌处理的地方，防止污染。
(11)离心制得DNA溶液：将晾干的吸附柱CB3转移到新的离心管中，然后向吸附膜的中间部位垂直悬空滴加洗脱缓冲液TE或ddH2O，体积在50～100ml范围内。室温下静置2～5min，13800g离心2min，离心管中所得液体即为DNA溶液。
(12)去除RNA：加入RNA酶，37℃消化40min，以去除DNA溶液中的RNA。
(13)DNA检测：0.8％琼脂糖凝胶电泳检测所提取的DNA，电泳条件为稳压100V，电泳时间为为40min。使用核酸蛋白质检测仪检测DNA的A 260/A280和A260/A230，并且使用Qubit 2.0测定DNA浓度。
2、基因组重测序文库构建与测序
所提取的基因组DNA经质量检测合格后用于基因组重测序文库构建，具体流程如下：
(1)使用仪器Covaris M220对基因组DNA进行随机打断，电泳回收300～500bp片段。
(2)利用TruSeqTM DNA Sample Prep Kit制备文库：首先在DNA片段两端加接头A和B，然后筛选去除接头自连片段，随后利用琼脂糖凝胶电泳进行片段筛选，保留一端是A接头、一端是B接头的片段，最后利用氢氧化钠变性，产生单链DNA片段。
(3)利用TruSeq PE Cluster Kit制备cluster：DNA片段的一端与引物碱基互补，固定在芯片上，另一端随机与附近的另外一个引物互补，也被固定住，形成“桥”，然后进行PCR扩增，产生DNA簇，最后DNA扩增子线性化成为单链。
(4)利用Illumina HiSeq 4000测序平台进行双末端(PE150)测序。
3、测序数据质控
上机测序后产生大量原始测序数据，其中会包含测序接头序列、低质量读段、N率较高序列及长度过短序列，为保证后续的生物信息分析的准确性，首先对原始测序数据进行过滤，从而得到高质量的测序数据以保证后续分析的顺利进行，具体步骤如下：
(1)去除reads中的接头序列，去除由于接头自连等原因导致没有插入片段的reads。
(2)将序列末端质量较低(质量值小于20)的碱基修剪掉，如剩余序列中仍然有质量值小于10的碱基，则将整条序列剔除，否则保留。
(3)去除含N比率超过10％的reads。
(4)将去除接头及质量修剪后长度小于50bp的序列舍弃。
4、比对参考基因组
以龙须菜已发表的genome survey组装结果序列作为参考基因组，利用BWA(http://bio-bwa.sourceforge.net/)比对软件将高质量reads比对回参考基因。接着利用SAMtools(http://samtools.sourceforge.net/)软件去除PCR-duplication产生的测序片段，再利用本地脚本计算出测序片段的深度与覆盖度，从而获得数据的比对情况。
5、SNP检测与分型
基于测序所得高质量数据与参考基因组序列比对的结果，利用GATK软件来检测SNP位点，结合每个碱基的测序质量，计算每个位点可能的碱基类型，再通过贝叶斯统计原理推断其中具有最大可能性的类型，判断该位点碱基是否具有多态性，然后过滤质量值小于20的位点，最终得到高可信度的SNP数据集。
6、特异性SNP标记的筛选
将鲁龙1号，981和2007这3个样品的SNP数据集去除任一样品中基因型缺失的位点，然后筛选出鲁龙1号与其余2个良种的基因型均不相同的位点，即鲁龙1号的特异性SNP标记，共2个。
7、特异性SNP标记的引物设计及Sanger测序检验
将上述2个SNP位点比对回参考基因组，由于2个SNP之间相差8bp，因此取包含2个SNP在内的前后各200bp共计400bp序列，设计出1对引物如下：
引物正向序列为：AAGGTAGCAGGAGCGTTCGTAAGC；
引物反向序列为：GCGATGGCGTTATTGAACCAACCC；
该引物扩增产物长度为325bp。
以鲁龙1号、981、2007的基因组DNA为模板，分别用该引物进行PCR扩增。
PCR反应条件为：94℃预变性5min；35个循环内94℃变性1min，60℃退火1min，72℃延伸1min；72℃延伸5min；4℃，保存。
PCR扩增后用2％琼脂糖凝胶电泳检测扩增结果，目的条带明显且无非特异性扩增条带，随后进行Sanger测序。
本发明的有益效果为：通过开发鲁龙1号良种的特异性SNP标记，从分子水平上将其与良种981和良种2007区分开来，解决龙须菜在栽培过程中产生的种苗混杂现象，避免对龙须菜栽培产业造成重大损失。
附图说明
图1为鲁龙1号良种的特异性SNP标记图；
图2为PCR反应体系图，以鲁龙1号、981、2007的基因组DNA为模板，分别用该引物进行PCR扩增；
图3示出鲁龙1号特异性SNP标记zcs1的Sanger测序检验结果，图中A、B、C分别代表鲁龙1号、981和2007；A1/A2、B1/B2、C1/C2分别代表每个良种的2个平行样；虚线框内是zcs1在不同样品中的Sanger测序峰，根据峰图读取出zcs1的碱基型并用黑色字体标识；
图4示出鲁龙1号特异性SNP标记zcs2的Sanger测序检验结果，图中A、B、C分别代表鲁龙1号、981和2007；A1/A2、B1/B2、C1/C2分别代表每个良种的2个平行样；虚线框内是zcs2在不同样品中的Sanger测序峰，根据峰图读取出zcs2的碱基型并用黑色字体标识。
具体实施方式
为了便于理解，下面结合附图通过实施例对本发明作进一步详细说明。
鲁龙1号龙须菜的特异性SNP标记，具体包括如下步骤：
1.基因组DNA提取与质量检测
采用天根植物基因组DNA提取试剂盒进行良种981、2007和鲁龙1号的基因组DNA的提取，具体步骤如下：
(1)提取工作开始之前，先用65℃水浴对缓冲液GP1进行预热，并在其中加入0.1％的β-巯基乙醇。而研钵、研杵和药匙则需要利用液氮预冷。
(2)藻体经无菌处理后，称取约100mg，放入预冷的研砵中，加入液氮充分研磨，使藻体完全破碎成粉末即可。研磨的过程中一定要注意及时添加液氮，并且缓慢加入，防止粉末被快速气化的液氮带出，造成浪费。
(3)用药匙将藻体粉末转移到(1)中预热的缓冲液GP1中，并迅速颠倒混匀，水浴30～40min，水浴过程中每隔5min要进行一次颠倒混匀，保证裂解充分。
(4)裂解结束后，加入700ml氯仿，迅速颠倒混匀，13800g离心5min。
(5)将离心后的上层透明水相转移到新的离心管内，注意吸取时不要混入蛋白层，在新离心管中加入700μl缓冲液GP2，颠倒混匀。
(6)将新离心管内混匀后的液体全部转移到吸附柱CB3中，13800g离心1min，将收集管中的液体弃去。
(7)完成上述操作后，向吸附柱CB3中加入500ml缓冲液GD，13800g离心1min，将收集管中的液体弃去。
(8)向吸附柱CB3中加入600μl漂洗液PW，13800g离心1min，将收集管中的液体弃去。该步骤重复两次。
(9)将吸附柱CB3于13800g离心2min，将收集管中的液体弃去。
(10)将吸附柱CB3置于超净工作台内晾干至无残留酒精气味。需要严格注意吸附柱的管口和底部不能触碰到其他未经灭菌处理的地方，防止污染。
(11)将晾干的吸附柱CB3转移到新的离心管中，然后向吸附膜的中间部位垂直悬空滴加洗脱缓冲液TE或ddH2O，体积在50～100ml范围内。室温下静置2～5min，13800g离心2min，离心管中所得液体即为DNA溶液。
(12)加入RNA酶，37℃消化40min，以去除DNA溶液中的RNA。
(13)0.8％琼脂糖凝胶电泳检测所提取的DNA，电泳条件为稳压100V，电泳时间为为40min。使用核酸蛋白质检测仪检测DNA的A 260/A280和A260/A230，并且使用Qubit 2.0测定DNA浓度。
2.基因组重测序文库构建与测序
所提取的基因组DNA经质量检测合格后用于基因组重测序文库构建，具体流程如下：
(1)使用仪器Covaris M220对基因组DNA进行随机打断，电泳回收回收300～500bp片段；
(2)利用TruSeqTM DNA Sample Prep Kit制备文库：首先在DNA片段两端加接头A和B，然后筛选去除接头自连片段，随后利用琼脂糖凝胶电泳进行片段筛选，保留一端是A接头、一端是B接头的片段，最后利用氢氧化钠变性，产生单链DNA片段；
(3)利用TruSeq PE Cluster Kit制备cluster：DNA片段的一端与引物碱基互补，固定在芯片上，另一端随机与附近的另外一个引物互补，也被固定住，形成“桥”，然后进行PCR扩增，产生DNA簇，最后DNA扩增子线性化成为单链；
(4)利用Illumina HiSeq 4000测序平台进行双末端(PE150)测序。
3.测序数据质控
上机测序后产生大量原始测序数据，其中会包含测序接头序列、低质量读段、N率较高序列及长度过短序列，为保证后续的生物信息分析的准确性，首先对原始测序数据进行过滤，从而得到高质量的测序数据以保证后续分析的顺利进行，具体步骤如下：
(1)去除reads中的接头序列，去除由于接头自连等原因导致没有插入片段的reads；
(2)将序列末端质量较低(质量值小于20)的碱基修剪掉，如剩余序列中仍然有质量值小于10的碱基，则将整条序列剔除，否则保留；
(3)去除含N比率超过10％的reads；
(4)将去除接头及质量修剪后长度小于50bp的序列舍弃。
质量过滤中使用软件包括：SeqPrep(https://github.com/jstjohn/SeqPrep)和Sickle(https://github.com/najoshi/sickle)。
4.比对参考基因组
以龙须菜已发表的genome survey组装结果序列作为参考基因组，利用BWA(http://bio-bwa.sourceforge.net/)比对软件将高质量reads比对回参考基因组。接着利用SAMtools(http://samtools.sourceforge.net/)软件去除PCR-duplication产生的测序片段，再利用本地脚本计算出测序片段的深度与覆盖度，从而获得数据的比对情况。
5.SNP检测与分型
基于测序所得高质量数据与参考基因组序列比对的结果，利用GATK(http://www.broadinstitute.org/gatk/)软件来检测SNP位点，结合每个碱基的测序质量，计算每个位点可能的碱基类型，再通过贝叶斯统计原理推断其中具有最大可能性的类型，判断该位点碱基是否具有多态性，如果是，即得到与参考序列相同和不同的碱基类型。过滤质量值小于20的位点，最终得到高可信度的SNP数据集。
6.特异性SNP标记的筛选
如图1所示，将鲁龙1号，981和2007这3个样品的SNP数据集去除任一样品中基因型缺失的位点，然后筛选出鲁龙1号与其余2个良种的基因型均不相同的位点，即鲁龙1号的特异性SNP标记，共2个。
7.特异性SNP标记的引物设计及Sanger测序检验
将上述2个SNP位点比对回参考基因组，由于2个SNP之间相差8bp，因此取包含2个SNP在内的前后各200bp共计400bp序列，设计出1对引物如下：
引物正向序列为：AAGGTAGCAGGAGCGTTCGTAAGC；
引物反向序列为：GCGATGGCGTTATTGAACCAACCC；
该引物扩增产物长度为325bp。
以鲁龙1号、981、2007的基因组DNA为模板，分别用该引物进行PCR扩增，PCR反应体系如图2。
PCR反应条件为：94℃预变性5min；35个循环内94℃变性1min，60℃退火1min，72℃延伸1min；72℃延伸5min；4℃，保存。
PCR扩增后用2％琼脂糖凝胶电泳检测扩增结果，目的条带明显且无非特异性扩增条带，随后进行Sanger测序。特异性标记zcs1的测序结果如图3所示，鲁龙1号的测序峰图在该SNP位点处显示为两个峰C和T，而981和2007均只有一个单峰C，而且两次平行样的结果一致，验证了重测序SNP分型结果，证明zcs1是鲁龙1号的特异性标记。
如图4，鲁龙1号的测序峰图在该SNP位点处显示为一个单峰G，981和2007在该位点均显示两个峰G和T。两次平行实验的结果相一致，验证了重测序SNP分型的结果，说明zcs2是鲁龙1号的特异性标记。可将zcs1和zcs2标记结合起来共同鉴定鲁龙1号。

Claims (6)

1.筛选鲁龙1号龙须菜的特异性SNP标记的方法，其特征在于：具体包括如下步骤：
一.基因组DNA提取与质量检测：采用天根植物基因组DNA提取试剂盒进行良种981、2007和鲁龙1号的提取；
二.基因组重测序文库构建与测序：所提取的基因组DNA经质量检测合格后用于基因组重测序文库构建与测序；
三.测序数据质控：对原始测序数据进行过滤，从而得到高质量的测序数据；
四.比对参考基因组；
五.SNP检测与分型；
六.特异性SNP标记的筛选；
七.特异性SNP标记的引物设计及Sanger测序检验；
步骤六中，将鲁龙1号，981和2007这3个样品的SNP数据集去除任一样品中基因型缺失的位点，然后筛选出鲁龙1号与其余2个良种的基因型均不相同的位点，即鲁龙1号的特异性SNP标记，共2个；
步骤七中，特异性SNP标记的引物设计及Sanger测序检验方法为：将上述2个SNP位点比对回参考基因组，由于2个SNP之间相差8bp，因此取包含2个SNP在内的前后各200bp共计400bp序列，设计出1对引物如下：
引物正向序列为：AAGGTAGCAGGAGCGTTCGTAAGC；
引物反向序列为：GCGATGGCGTTATTGAACCAACCC；
该引物扩增产物长度为325bp；
以鲁龙1号、981、2007的基因组DNA为模板，分别用该引物进行PCR扩增；PCR扩增后用2％琼脂糖凝胶电泳检测扩增结果，目的条带明显且无非特异性扩增条带，随后进行Sanger测序。
2.根据权利要求1所述的筛选鲁龙1号龙须菜的特异性SNP标记的方法，其特征在于：步骤一中，基因组DNA的提取步骤如下：
(1)缓冲液GP1预热：提取工作开始之前，先对缓冲液GP1进行预热，并在其中加入0.1％的β-巯基乙醇；
(2)藻体研磨：藻体经无菌处理后，称取100mg，放入预冷的研砵中，加入液氮充分研磨；
(3)藻体裂解：用药匙将藻体粉末转移到(1)中预热的缓冲液GP1中，并迅速颠倒混匀，水浴30～40min进行充分裂解；
(4)离心：裂解结束后，加入700ml氯仿，迅速颠倒混匀，13800g离心5min；
(5)上层透明水相转移：将离心后的上层透明水相转移到新的离心管内，在新离心管中加入700μl缓冲液GP2，颠倒混匀；
(6)离心：将新离心管内混匀后的液体全部转移到吸附柱CB3中，13800g离心1min，将收集管中的液体弃去；
(7)离心：向吸附柱CB3中加入500ml缓冲液GD，13800g离心1min，将收集管中的液体弃去；
(8)离心：向吸附柱CB3中加入600μl漂洗液PW，13800g离心1min，将收集管中的液体弃去，该步骤重复两次；
(9)离心：将吸附柱CB3于13800g离心2min，将收集管中的液体弃去；
(10)吸附柱CB3晾干：将吸附柱CB3置于超净工作台内晾干；
(11)离心制得DNA溶液：将晾干的吸附柱CB3转移到新的离心管中，然后滴加洗脱缓冲液TE或ddH2O，离心所得液体即为DNA溶液；
(12)去除RNA：加入RNA酶，37℃消化40min，以去除DNA溶液中的RNA；
(13)DNA检测：0.8％琼脂糖凝胶电泳检测所提取的DNA，电泳条件为稳压100V，电泳时间为40min。
3.根据权利要求2所述的筛选鲁龙1号龙须菜的特异性SNP标记的方法，其特征在于：步骤二中，基因组重测序文库构建与测序流程如下：
(1)使用仪器Covaris M220对基因组DNA进行随机打断，电泳回收300～500bp片段；
(2)利用TruSeqTMDNA Sample Prep Kit制备文库；
(3)利用TruSeq PE Cluster Kit制备cluster；
(4)利用Illumina HiSeq 4000测序平台进行双末端(PE150)测序。
4.根据权利要求3所述的筛选鲁龙1号龙须菜的特异性SNP标记的方法，其特征在于：步骤三中，对原始测序数据进行过滤步骤如下：
(1)去除reads中的接头序列，去除由于接头自连原因导致没有插入片段的reads；
(2)将序列末端质量值小于20的碱基修剪掉，如剩余序列中仍然有质量值小于10的碱基，则将整条序列剔除，否则保留；
(3)去除含N比率超过10％的reads；
(4)将去除接头及质量修剪后长度小于50bp的序列舍弃。
5.根据权利要求4所述的筛选鲁龙1号龙须菜的特异性SNP标记的方法，其特征在于：步骤四中，基因组比对参考方法为：以龙须菜已发表的genome survey组装结果序列作为参考基因组，利用BWA比对软件将高质量reads比对回参考基因，接着利用SAMtools软件去除PCR-duplication产生的测序片段，再利用本地脚本计算出测序片段的深度与覆盖度，从而获得数据的比对情况。
6.根据权利要求5所述的筛选鲁龙1号龙须菜的特异性SNP标记的方法，其特征在于：步骤五中，SNP检测与分型方法为：基于测序所得高质量数据与参考基因组序列比对的结果，利用GATK软件来检测SNP位点，结合每个碱基的测序质量，计算每个位点可能的碱基类型，再通过贝叶斯统计原理推断其中具有最大可能性的类型，判断该位点碱基是否具有多态性，过滤质量值小于20的位点，最终得到高可信度的SNP数据集。
CN201911291611.5A 2019-12-16 2019-12-16 龙须菜鲁龙1号良种的特异性snp标记及其应用 Active CN110846435B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911291611.5A CN110846435B (zh) 2019-12-16 2019-12-16 龙须菜鲁龙1号良种的特异性snp标记及其应用

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911291611.5A CN110846435B (zh) 2019-12-16 2019-12-16 龙须菜鲁龙1号良种的特异性snp标记及其应用

Publications (2)

Publication Number Publication Date
CN110846435A CN110846435A (zh) 2020-02-28
CN110846435B true CN110846435B (zh) 2022-09-06

Family

ID=69609235

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911291611.5A Active CN110846435B (zh) 2019-12-16 2019-12-16 龙须菜鲁龙1号良种的特异性snp标记及其应用

Country Status (1)

Country Link
CN (1) CN110846435B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112266974B (zh) * 2020-11-25 2023-09-15 江苏海洋大学 一种用于鉴定龙须菜性别的引物及鉴定方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103397083A (zh) * 2013-06-08 2013-11-20 中国海洋大学 鉴定龙须菜981良种的特异性探针及其应用
CN107354206A (zh) * 2017-07-10 2017-11-17 中国烟草总公司郑州烟草研究院 用于鉴定烤烟南江3号的引物组合及试剂盒、应用及检测方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103397083A (zh) * 2013-06-08 2013-11-20 中国海洋大学 鉴定龙须菜981良种的特异性探针及其应用
CN107354206A (zh) * 2017-07-10 2017-11-17 中国烟草总公司郑州烟草研究院 用于鉴定烤烟南江3号的引物组合及试剂盒、应用及检测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Heterozygous Single Nucleotide Polymorphic Loci in Haploid Gametophytes of Gracilariopsis lemaneiformis (Rhodophyta);Yiyi Hu et al;《Frontiers in Genetics》;20191206;1-9 *
龙须菜"鲁龙1号";周伟等;《中国水产》;20160305(第03期);80-81 *

Also Published As

Publication number Publication date
CN110846435A (zh) 2020-02-28

Similar Documents

Publication Publication Date Title
CN114672587B (zh) 一种与番木瓜果实果糖含量相关的snp分子标记及其扩增引物、检测试剂盒和应用
CN107077537A (zh) 用短读测序数据检测重复扩增
CN106939342B (zh) 一种与谷子米色连锁的snp标记、引物及应用
CN105385755A (zh) 一种利用多重pcr技术进行snp-单体型分析的方法
CN108315439A (zh) 一种与瓦氏黄颡鱼生长相关的snp分子标记及其应用
CN114672588B (zh) 一种SNP分子标记Cpa03g017250：8及其扩增引物、检测试剂盒和应用
CN114672586A (zh) 一种与番木瓜果实宽度性状相关的snp分子标记及其扩增引物、检测试剂盒和应用
CN108531582A (zh) 一种检测人类胚胎α-地中海贫血基因突变的引物组合及方法
CN105969879B (zh) 一种高通量检测AhFAD2A基因突变位点分型的引物组及检测方法
CN106755561A (zh) 一种与大豆根干重相关的qtl、snp分子标记及应用
CN105525012B (zh) 一种花生杂交种的分子鉴定方法
CN110846429A (zh) 一种玉米全基因组InDel芯片及其应用
CN105925680A (zh) 一种四倍体马铃薯高通量测序开发标记的方法及其应用
CN107988385B (zh) 一种检测肉牛PLAG1基因Indel标记的方法及其专用试剂盒
CN113637787B (zh) 一种与油茶单果质量相关的dna片段及其应用
CN110846435B (zh) 龙须菜鲁龙1号良种的特异性snp标记及其应用
CN108642201B (zh) 与谷子株高性状相关的snp标记及其检测引物和应用
CN108660238A (zh) 基于gbs技术的燕麦抗旱性相关snp分子标记及其应用
CN106191253B (zh) 基于gbs技术的北京鸭简化基因测序方法
CN110484629B (zh) 一种与三疣梭子蟹生长性状相关的微卫星标记、其引物及应用
CN105200051B (zh) 一种皱纹盘鲍中国和日本群体鉴别用snp标记
CN117051151A (zh) 一种辣椒5k液相芯片及其应用
CN108707685B (zh) 与谷子分蘖数性状相关的snp标记及其检测引物和应用
CN116064893A (zh) 指状蔷薇珊瑚微卫星标记位点及其引物组和应用
CN113337578A (zh) 一种基于转录组数据高效筛选水产动物阳性snp的方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant