CN107201408A - 一种基于转录组测序开发剑麻ssr引物的方法 - Google Patents
一种基于转录组测序开发剑麻ssr引物的方法 Download PDFInfo
- Publication number
- CN107201408A CN107201408A CN201710577729.9A CN201710577729A CN107201408A CN 107201408 A CN107201408 A CN 107201408A CN 201710577729 A CN201710577729 A CN 201710577729A CN 107201408 A CN107201408 A CN 107201408A
- Authority
- CN
- China
- Prior art keywords
- ssr
- sisal hemp
- sequence
- transcript profile
- label primer
- 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.)
- Granted
Links
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/6811—Selection methods for production or design of target specific oligonucleotides or binding molecules
-
- 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
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B25/00—ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
-
- 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
Landscapes
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Organic Chemistry (AREA)
- Engineering & Computer Science (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Physics & Mathematics (AREA)
- Biotechnology (AREA)
- Analytical Chemistry (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Genetics & Genomics (AREA)
- General Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Immunology (AREA)
- Microbiology (AREA)
- Biochemistry (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Theoretical Computer Science (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Medical Informatics (AREA)
- Botany (AREA)
- Mycology (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明公开了一种基于转录组测序开发剑麻SSR标记引物的方法，属于分子生物学技术领域；其制备方法包括Total RNA 的提取与转录组文库的构建、转录组测序及序列质量分析、转录本组装与SSR标记分析、SSR引物设计、引物PCR扩增与扩增产物电泳检测五个步骤，本发明通过转录组测序的方法，获得剑麻某一特定时期的转录组序列，然后通过生物信息学分析软件进行SSR标记分析，开发剑麻SSR标记引物，标记开发效率高，填补了目前剑麻SSR标记稀缺的空白。
Description
技术领域
本发明涉及一种基于转录组测序开发剑麻SSR标记引物的方法，属于分子生物学技术领域。
背景技术
剑麻属龙舌兰属多年生草本植物，原产于墨西哥等热带、亚热带地区，我国保存的种质资源严重缺乏，仅有100余份。现有的种质资源主要是通过麻园选育、杂交和引种获得，遗传背景模糊，命名极不规范，给种质资源的保存和应用以及育种工作带来了极大的不便。因此需要一种快速有效的方法对现有的种质资源进行遗传多样性分析和资源鉴定。
微卫星DNA（Microsatellite DNA）又称简单重复序列（Simple Sequence Repeat,SSR）是由1-6个核苷酸为重复单位串联而成的长达几十个核苷酸的重复序列。其侧翼通常都是保守性高的单拷贝序列，因此可根据侧翼序列设计引物进行PCR扩增，由于重复片段大小或重复次数的差异，显示微卫星位点的多态性。与其它分子标记相比，SSR标记具有多态性高、共显性遗传、技术简单、重复性好，特异性强以及操作便利等优点，广泛应用于遗传多样性分析、种质资源鉴定以及遗传图谱构建等领域。但对于SSR标记的应用前提是，必须首先要从该物种中获取重复序列两侧的序列信息，并设计引物，而后才能被应用。
SSR标记可分为基因组SSR(gSSR)和表达序列标签SSR（EST-SSR），与gSSR标记相比，通过表达序列标签（Expressed sequence tag，EST）序列开发SSR标记更经济，效率更高，而且在不同属内的通用性更好。随着分子生物学和测序技术的快速发展，EST序列数据急剧增加，加上生物信息学的飞速发展使得大批量数据处理成为可能，因此给SSR标记的开发提供了大量的序列信息资源和技术支持。并且EST-SSR标记是基于某一时期的表达标签序列，能直接与功能基因相关，在分子标记辅助育种如重要性状相关标记关联分析、分离和新基因的鉴定等方面均有极高的应用价值。
目前剑麻的全基因组信息未知，SSR标记缺乏，剑麻SSR标记引物未见报道。随着高通量测序技术和生物信息学的飞速发展，对于无参考基因组的转录组分析，技术非常成熟。因此，利用某一发育时期或逆境胁迫下的剑麻为材料，开展剑麻转录组研究成为可能，这不仅可以挖掘剑麻自身的功能基因，同时对解决剑麻种质资源背景模糊，SSR标记缺乏等实际问题具有十分重要的现实意义。
发明内容
本发明的目的是针对现有技术的不足，提供一种基于转录组测序开发剑麻SSR标记引物的方法，通过转录组测序的方法，获得剑麻某一特定时期的转录组序列，然后通过生物信息学分析软件开发剑麻SSR标记，为后续利用SSR标记引物进行剑麻遗传多样性研究、种质资源的鉴定提供可靠的技术手段。
为实现上述目的，本方法采用的技术方案是：
一种基于转录组测序开发剑麻SSR标记引物的方法，包括如下步骤：
（1）Total RNA 的提取与转录组文库的构建
采用Trizol裂解法提取剑麻叶片总RNA，用带有Oligo(dT)的磁珠富集mRNA，反转录并合成双链cDNA，经QiaQuick PCR Purification Kit 纯化后，在cDNA末端添加腺嘌呤核苷进行末端修复并连接测序接头，然后用琼脂糖凝胶电泳回收目标片段，最后对回收片段进行PCR扩增，扩增产物即为转录组文库；
（2）转录组测序及序列质量分析
用Illumina 2500测序平台进行转录组测序，所得原始序列以Fasta Q格式保存；由于原始序列含有低质量的序列，在进行数据分析之前先对原始数据进行质量分析和过滤。
（3）转录本组装与SSR分析
利用Trinity软件将所得有效序列拼接成一个完整的转录组，作为后续分析的参考序列，取每条基因中最长的转录本作为Unigene，采用MISA1.0软件对每个Unigene进行简单序列重复（SSR）分析；
（4）SSR标记引物设计、扩增与检测
采用引物设计软件Primer3进行SSR引物设计，从所设计的SSR引物中随机选取100对引物，首先以热麻1号基因组DNA为模板，采用Touch-down PCR程序进行扩增，扩增产物经3%琼脂糖凝胶电泳，初步筛选有稳定产物的SSR引物，然后以6份剑麻种质DNA为模板，采用Touch-down PCR程序，对有稳定产物的SSR引物进行PCR扩增，扩增产物经3%琼脂糖凝胶电泳，进一步筛选多态性的SSR标记引物，检测SSR标记引物的有效性。
步骤（1）所述的接头序列为SEQ ID No.1- SEQ ID No.2。
步骤（2）所述的低质量序列是指不确定碱基N比例大于10%的序列和低质量碱基含量大于50%的序列；低质量碱基指Q≤5的碱基。
步骤（4）所述的SSR标记筛选参数为：单核苷酸的重复次数大于或等于10次，二核苷酸重复次数大于或等于6次，三核苷酸、四核苷酸、五核苷酸、六核苷酸的重复次数大于或等于5次。
步骤（4）所述的SSR标记引物设计参数为：引物长度18-25bp，退火温度Tm 56-65℃，预期产物长度为100-300bp。
步骤（4）所述的Touch-down PCR扩增程序为：首先94℃ 15s，60℃ 15s，72℃ 30s，16个循环，每个循环退火温度降低0.7℃；然后进入下一个扩增阶段：94℃ 15s，50℃ 15s，72℃ 30s，15个循环，最后72℃延伸60min，扩增产物4℃保存备用。
步骤（4）所述的有多态性的18对SSR标记引物组序列为SEQ ID No.3- SEQ IDNo.38。
步骤（4）所述的6份剑麻种质为热麻1号、H.11648、番麻、普通剑麻、桂幅4号和广西76416。
本发明的有益效果为：
（1）通过转录组测序的方法，获得剑麻某一特定时期的转录组序列，然后通过生物信息学分析软件开发剑麻SSR标记，开发了一种基于转录组测序开发剑麻SSR标记引物的方法，开发效率更高。
（2）本发明为剑麻SSR标记开发提供了一条新的有效途径，填补了目前剑麻SSR标记引物稀缺的空白。
附图说明
图1 剑麻SSR标记引物开发流程图；
图2 琼脂糖凝胶电泳检测部分剑麻SSR标记引物扩增产物；
图3 剑麻SSR标记引物多态性筛选。
具体实施方式
下面通过实例对本发明做进一步详细说明，这些实例仅用来说明本发明，并不限制本发明的范围。未加特殊说明，转录组文库构建与测序均按标准流程进行，所有试剂盒操作均按试剂盒说明书进行，所有的试剂均为生物试剂。
本发明所提供的18对SSR标记引物组均来自热麻1号转录组序列，其引物核酸序列分别为SEQ ID No.3- SEQ ID No.38。所用的植物材料为热麻1号，病原菌为烟草疫霉。
一．RNA提取与转录组文库构建
取保存于本实验室的烟草疫霉接种于马铃薯培养基（PDA）上，28℃培养1周后，接种热麻1号叶片，取不同接种时间的热麻1号叶片进行转录组测序。具体步骤如下：
用灭菌的大头针将叶片正面刺伤，取直径为5mm的菌饼，将菌饼的菌丝生长面贴在伤口的位置，用无菌湿棉花保湿，然后用保鲜膜包裹叶片，置于25-30℃条件下培养，并分别在接种前、接种24小时、36小时、48小时和72小时取叶片，液氮速冻后保存于-80℃备用。
采用Trizol 裂解法提取热麻1号叶片总RNA，具体操作按说明书进行。采用Nanodrop-2000分光光度计和Bioanalyzer2100生物分析仪对RNA质量进行检测。质量合格的RNA样品将用于转录组文库构建。文库构建采用Illumina公司的文库构建试剂盒进行。首先用带有Oligo（dT)的磁珠富集mRNA，并将mRNA片段化并反转录成第一链cDNA，然后合成双链cDNA，经QiaQuick PCR Purification Kit 纯化后，在cDNA末端添加腺嘌呤核苷进行末端修复并连接测序接头，然后用琼脂糖凝胶电泳回收目标片段，最后对回收片段进行PCR扩增，扩增产物即为转录组文库，将用于后续的转录组测序。
二．测序与转录组分析
采用Illumina 2500测序平台进行转录组测序，所得原始序列以Fasta Q格式保存。由于原始序列含有接头污染和低质量的序列，因此，为了防止这些序列对后续分析产生不利影响，在进行序列分析前，先要去掉测序时的接头序列，不确定碱基（N）比例大于10%的序列以及低质量碱基（Q≤5）含量大于50%的序列，所得序列即为有效的转录组序列，并用于后续的序列分析。
利用Trinity（版本为v2012-06-08，参数为默认参数）软件将所得有效序列拼接成一个完整的转录组，作为后续分析的参考序列，取每条基因中最长的转录本作为Unigene，热麻1号总计获得了103,326个转录本和70,110条Unigene序列。转录本和Unigene平均长度分别为726bp和645bp。转录本和Unigene具体数目分布如表1所示。
表1 转录本和Unigene拼接长度频数分布
三．SSR位点查找及引物设计
以热麻1号70,110条Unigene为材料，利用MISA1.0软件对Unigene序列进行SSR位点查找，查找含有1、2、3、4、5和6碱基重复的SSR位点，且重复次数依次不小于10、6、5、5、5和5次。总计查找到了13,175个SSR位点，SSR密度分布出现频率最高的依次为单核苷酸重复、二核苷酸重复和三核苷酸重复，分别为5001个、4339个和3676个。使用Primer3.0软件对SSR候选位点进行引物设计。引物设计参数为，引物长度18-26bp，GC含量40%-60%，退火温度Tm值55-65℃（上下游引物的Tm值相差不能大于5℃）；PCR目标产物在100-300bp；尽量避免产生引物二聚体，发夹结构、错配等。总计设计了11,946对SSR引物。其中长度为20bp的引物最多，为10,270对，占总引物的85.97%。
四．SSR标记引物的有效性验证
从10,270对引物中随机选取100对引物进行PCR扩增，筛选有稳定产物条带的SSR标记引物，然后以6份剑麻种质基因组DNA为模板，对有稳定产物的SSR标记引物进行扩增，筛选有多态性的SSR标记引物，检测SSR标记引物的有效性。引物由生工生物工程（上海）有限公司合成，引物序列见核苷酸序列表(序列分别如SEQ ID No.3-42)。
序列说明：SED ID No.3-38为热麻1号SSR标记引物序列，其中SED ID No.3和4为一对SSR引物，SED ID No.5和6为一对SSR引物，依次类推，18对引物的退火温度和扩增产物大小见表2。
表2 18对剑麻SSR引物的退火温度和扩增产物大小
1.DNA提取
采用天泽公司的柱式植物DNAout试剂盒提取热麻1号、H.11648、番麻、普通剑麻、桂幅4号和广西76416等6份剑麻种质的基因组DNA。取叶片1克，经液氮速冻后快速研磨成粉末，先加入750µL 65℃预热的裂解液，充分混匀后65℃预热5-10分钟，室温13,000rpm离心5min，取上清液500 µL分别进行抽提、漂洗后过柱，最后用100 µL洗脱液洗脱2次，洗脱液即为提取的DNA， 4℃保存备用。
2.PCR反应体系的建立
20µLPCR反应体系中各组份的浓度及使用量：2×NOVA Taq-Plus PCR Forest Mix(江苏愚公生命科技有限公司):2µL；引物(10μmol/L) F:0.2 µL，R:0.2µL；DNA:50ng(1µL)；灭菌ddH2O :8.6µL。NOVA Taq-Plus PCR Forest Mix购自江苏愚公生命科技有限公司。
Touch-down PCR扩增程序：94℃(15s),60℃(15s)(△℃=-0.7，即每增加一个循环，退火温度降低0.7℃)，72℃（30s）(16个循环)，然后进入下一个扩增阶段，94℃（15s）,50℃(15s), 72℃(30s)(15个循环)，最后72℃（60min），扩增产物4℃保存备用。
3.PCR产物的检测
采用3%琼脂糖凝胶（100mL 1XTBE缓冲液中加入3g琼脂糖）电泳检测PCR产物。将20µLPCR产物全部上样，先150V电泳10min，然后120V电泳40min，电泳结果于凝胶成像系统拍照保存。
本发明采用Touch-down PCR程序，从100对热麻1号SSR标记引物中，用琼脂糖凝胶电泳筛选出了66对有稳定产物条带的SSR标记引物（附图2），然后用6份剑麻种质，从66对SSR标记引物中，筛选出了18对有多态性的SSR标记引物（附图3），证实了通过转录组测序开发剑麻SSR标记引物是切实可行的。
序列表
<110>中国热带农业科学院南亚热带作物研究所
<120>一种基于转录组测序开发剑麻SSR引物的方法
<141>2017-07-14
<160>38
<170>SIPOSequenceListing 1.0
<210>1
<211>58
<212>DNA
<213>人工序列
<400>1
aatgatacgg cgaccaccga gatctacact ctttccctac acgacgctct tccgatct 58
<210>2
<211>57
<212>DNA
<213>人工序列
<400>2
gatcggaaga gcacacgtct gaactccagt cacatctcgt atgccgtctt ctgcttg 57
<210>3
<211>20
<212>DNA
<213>人工序列
<400>3
tcgcgtgcac caacaatttc 20
<210>4
<211>20
<212>DNA
<213>人工序列
<400>4
gtagcggatg taggagacgc 20
<210>5
<211>20
<212>DNA
<213>人工序列
<400>5
tgcttcgact cctgcttctg 20
<210>6
<211>20
<212>DNA
<213>人工序列
<400>6
agtggtggcc gtggaaatag 20
<210>7
<211>20
<212>DNA
<213>人工序列
<400>7
gtgtgtgtgt gtgtgttggg 20
<210>8
<211>20
<212>DNA
<213>人工序列
<400>8
ggccgaatcc tttccactca 20
<210>9
<211>20
<212>DNA
<213>人工序列
<400>9
cgctcgtcct cttctttcgt 20
<210>10
<211>23
<212>DNA
<213>人工序列
<400>10
tcccatccaa tagtccccca 20
<210>11
<211>20
<212>DNA
<213>人工序列
<400>11
ggtcatgatg aaggccacca 20
<210>12
<211>20
<212>DNA
<213>人工序列
<400>12
gcgaacctgc attgctgaat 20
<210>13
<211>20
<212>DNA
<213>人工序列
<400>13
ccttaggctc cctgctgttc 20
<210>14
<211>23
<212>DNA
<213>人工序列
<400>14
ccacaagagc cgctaccatc 20
<210>15
<211>20
<212>DNA
<213>人工序列
<400>15
aacaaccaga gcccaaacca 20
<210>16
<211>20
<212>DNA
<213>人工序列
<400>16
ggggaggtgg tttggtgatc 20
<210>17
<211>20
<212>DNA
<213>人工序列
<400>17
ggttagggtt cttggtgggg 20
<210>18
<211>20
<212>DNA
<213>人工序列
<400>18
gcttcctgat cttcttgttg gc 22
<210>19
<211>20
<212>DNA
<213>人工序列
<400>19
aaaatccatg aggcggctga 20
<210>20
<211>20
<212>DNA
<213>人工序列
<400>20
tagtagctag gcccaggcaa 20
<210>21
<211>20
<212>DNA
<213>人工序列
<400>21
acagcacgag aaatgagctc a 21
<210>22
<211>20
<212>DNA
<213>人工序列
<400>22
ccgatccggc gtaattctct 20
<210>23
<211>20
<212>DNA
<213>人工序列
<400>23
gccttctccc acggaatcaa 20
<210>24
<211>20
<212>DNA
<213>人工序列
<400>24
tgtggagtgt gatgggagtg 20
<210>25
<211>20
<212>DNA
<213>人工序列
<400>25
tggagggtga tggatagggg 20
<210>26
<211>20
<212>DNA
<213>人工序列
<400>26
gatgaggcca tcgttttggt 20
<210>27
<211>20
<212>DNA
<213>人工序列
<400>27
agagttgcca gatgtgtgca 20
<210>28
<211>20
<212>DNA
<213>人工序列
<400>28
aggtgggatt cttgcggatg 20
<210>29
<211>20
<212>DNA
<213>人工序列
<400>29
tcaaaagcaa cgaacagcgg 20
<210>30
<211>20
<212>DNA
<213>人工序列
<400>30
cgacttcctc atcgatgcga 20
<210>31
<211>20
<212>DNA
<213>人工序列
<400>31
gcaggccctg tagtttgact 20
<210>32
<211>20
<212>DNA
<213>人工序列
<400>32
ttcgtgccca gtttctcctc 20
<210>33
<211>20
<212>DNA
<213>人工序列
<400>33
atcttcaggt ttccgctgca 20
<210>34
<211>20
<212>DNA
<213>人工序列
<400>34
ccgagagaga gcgagagaga 20
<210>35
<211>20
<212>DNA
<213>人工序列
<400>35
accgcattca tcggtctctc 20
<210>36
<211>20
<212>DNA
<213>人工序列
<400>36
ggtcctcgct ctgatcttgg 20
<210>37
<211>20
<212>DNA
<213>热人工序列
<400>37
attgcttgaa gatggctgct 20
<210>38
<211>20
<212>DNA
<213>人工序列
<400>38
catgcatacc ttcctccccc 20
Claims (8)
1.一种基于转录组测序开发剑麻SSR标记引物的方法，其特征在于：包括如下步骤：
（1）Total RNA 的提取与转录组文库的构建
采用Trizol裂解法提取剑麻叶片总RNA，用带有Oligo的磁珠富集mRNA，反转录并合成双链cDNA，经QiaQuick PCR Purification Kit 纯化后，在cDNA末端添加腺嘌呤核苷进行末端修复并连接测序接头，然后用琼脂糖凝胶电泳回收目标片段，最后对回收片段进行PCR扩增，扩增产物即为转录组文库；
（2）转录组测序及序列质量分析
采用Illumina 2500测序平台进行转录组测序，所得原始序列以Fasta Q格式保存；由于原始序列含有接头污染和低质量的序列，在进行数据分析之前先对原始数据进行质量分析和过滤；
（3）转录本组装与SSR标记分析
利用Trinity软件将所得有效序列拼接成一个完整的转录组，作为后续分析的参考序列，取每条基因中最长的转录本作为Unigene，采用MISA1.0软件对每个Unigene进行简单序列重复（SSR）分析；
SSR标记引物设计、扩增与检测
采用引物设计软件Primer3进行SSR标记引物设计，从所设计的SSR标记引物中随机选取100对引物，首先以热麻1号基因组DNA为模板，采用Touch-down PCR程序进行扩增，扩增产物经3%琼脂糖凝胶电泳，初步筛选有稳定扩增产物的SSR引物，然后以6份剑麻种质DNA为模板，采用Touch-down PCR程序，对有稳定扩增产物的SSR引物进行PCR扩增，扩增产物经3%琼脂糖凝胶电泳，进一步筛选多态性的SSR标记引物，检测SSR标记引物的有效性。
2.根据权利要求1所述的一种基于转录组测序开发剑麻SSR标记引物的方法，其特征在于：步骤（1）所述的接头序列为SEQ ID No.1- SEQ ID No.2。
3.根据权利要求1所述的一种基于转录组测序开发剑麻SSR标记引物的方法，其特征在于：步骤（2）所述的低质量序列为带有测序接头的序列，不确定碱基N比例大于10%的序列和低质量碱基含量大于50%的序列；低质量碱基指Q≤5的碱基。
4.根据权利要求1所述的一种基于转录组测序开发剑麻SSR标记引物的方法，其特征在于：步骤（4）所述的SSR标记筛选参数为：单核苷酸的重复次数大于或等于10次，二核苷酸重复次数大于或等于6次，三核苷酸、四核苷酸、五核苷酸、六核苷酸的重复次数大于或等于5次。
5.根据权利要求1所述的一种基于转录组测序开发剑麻SSR标记引物的方法，其特征在于：步骤（4）所述的SSR标记引物设计参数为：引物长度18-25bp，退火温度Tm 56-65℃，预期产物长度为100-300bp。
6.根据权利要求1所述的一种基于转录组测序开发剑麻SSR标记引物的方法，其特征在于：步骤（4）所述的Touch-down PCR扩增程序为：首先94℃ 15s，60℃ 15s，72℃ 30s，16个循环，每个循环退火温度降低0.7℃；然后进入下一个扩增阶段：94℃ 15s，50℃ 15s，72℃30s，15个循环，最后72℃延伸60min，扩增产物4℃保存备用。
7.根据权利要求1所述的一种基于转录组测序开发剑麻SSR标记引物的方法，其特征在于：步骤（4）所述的有多态性的18对SSR标记引物组序列为SEQ ID No.3- SEQ ID No.38。
8.根据权利要求1所述的一种基于转录组测序开发剑麻SSR标记引物的方法，其特征在于：步骤（4）所述的6份剑麻种质为热麻1号、H.11648、番麻、普通剑麻、桂幅4号和广西76416。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710577729.9A CN107201408B (zh) | 2017-07-15 | 2017-07-15 | 一种基于转录组测序开发剑麻ssr引物的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710577729.9A CN107201408B (zh) | 2017-07-15 | 2017-07-15 | 一种基于转录组测序开发剑麻ssr引物的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107201408A true CN107201408A (zh) | 2017-09-26 |
| CN107201408B CN107201408B (zh) | 2020-10-02 |
Family
ID=59911995
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710577729.9A Active CN107201408B (zh) | 2017-07-15 | 2017-07-15 | 一种基于转录组测序开发剑麻ssr引物的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107201408B (zh) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108034696A (zh) * | 2018-02-02 | 2018-05-15 | 中南大学 | 一种基于转录组测序开发ssr引物的方法 |
| CN108611434A (zh) * | 2018-05-18 | 2018-10-02 | 中国热带农业科学院南亚热带作物研究所 | 剑麻内参基因act/gapdh及其应用 |
| CN108754018A (zh) * | 2018-07-27 | 2018-11-06 | 大连民族大学 | 一种刺五加目的基因ssr分子标记的筛选方法及应用 |
| CN109207560A (zh) * | 2018-09-19 | 2019-01-15 | 广西壮族自治区林业科学研究院 | 一种基于转录组测序开发叶子花ssr引物的方法 |
| CN109280694A (zh) * | 2018-10-23 | 2019-01-29 | 浙江省萧山棉麻研究所 | 基于全长转录组测序开发的姜荷花多态性ssr引物及试剂盒 |
| CN110029185A (zh) * | 2019-04-11 | 2019-07-19 | 广西壮族自治区亚热带作物研究所（广西亚热带农产品加工研究所） | 一种在转录组测序水平上进行剑麻叶纤维发育的研究方法 |
| CN111705160A (zh) * | 2020-07-31 | 2020-09-25 | 广西壮族自治区亚热带作物研究所（广西亚热带农产品加工研究所） | 龙舌兰麻cpSSR标记引物及其应用 |
| CN113215220A (zh) * | 2021-05-31 | 2021-08-06 | 福建农林大学 | 一种基于转录组测序的橄榄ssr分子标记开发方法 |
| CN114898802A (zh) * | 2022-07-14 | 2022-08-12 | 臻和(北京)生物科技有限公司 | 基于血浆游离dna甲基化测序数据的末端序列频率分布特征确定方法、评价方法及装置 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120283110A1 (en) * | 2011-04-21 | 2012-11-08 | Jay Shendure | Methods for retrieval of sequence-verified dna constructs |
| CN103233075A (zh) * | 2013-05-09 | 2013-08-07 | 南京农业大学 | 一种基于转录组测序开发菊属植物ssr引物的方法 |
| CN103642912A (zh) * | 2013-11-29 | 2014-03-19 | 中国农业科学院作物科学研究所 | 基于转录组测序开发绿豆ssr引物的方法 |
| CN106754886A (zh) * | 2017-01-19 | 2017-05-31 | 北京林业大学 | 基于转录测序获得黑果枸杞ssr引物的方法 |
-
2017
- 2017-07-15 CN CN201710577729.9A patent/CN107201408B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120283110A1 (en) * | 2011-04-21 | 2012-11-08 | Jay Shendure | Methods for retrieval of sequence-verified dna constructs |
| CN103233075A (zh) * | 2013-05-09 | 2013-08-07 | 南京农业大学 | 一种基于转录组测序开发菊属植物ssr引物的方法 |
| CN103642912A (zh) * | 2013-11-29 | 2014-03-19 | 中国农业科学院作物科学研究所 | 基于转录组测序开发绿豆ssr引物的方法 |
| CN106754886A (zh) * | 2017-01-19 | 2017-05-31 | 北京林业大学 | 基于转录测序获得黑果枸杞ssr引物的方法 |
Non-Patent Citations (1)
| Title |
|---|
| 张燕梅等: "剑麻 ISSR 反应体系的建立与分析", 《中国麻业科学》 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108034696B (zh) * | 2018-02-02 | 2020-07-28 | 中南大学 | 一种基于转录组测序开发ssr引物的方法 |
| CN108034696A (zh) * | 2018-02-02 | 2018-05-15 | 中南大学 | 一种基于转录组测序开发ssr引物的方法 |
| CN108611434A (zh) * | 2018-05-18 | 2018-10-02 | 中国热带农业科学院南亚热带作物研究所 | 剑麻内参基因act/gapdh及其应用 |
| CN108611434B (zh) * | 2018-05-18 | 2022-02-08 | 中国热带农业科学院南亚热带作物研究所 | 剑麻内参基因act/gapdh及其应用 |
| CN108754018A (zh) * | 2018-07-27 | 2018-11-06 | 大连民族大学 | 一种刺五加目的基因ssr分子标记的筛选方法及应用 |
| CN109207560A (zh) * | 2018-09-19 | 2019-01-15 | 广西壮族自治区林业科学研究院 | 一种基于转录组测序开发叶子花ssr引物的方法 |
| CN109280694B (zh) * | 2018-10-23 | 2021-09-10 | 浙江省萧山棉麻研究所 | 基于全长转录组测序开发的姜荷花多态性ssr引物及试剂盒 |
| CN109280694A (zh) * | 2018-10-23 | 2019-01-29 | 浙江省萧山棉麻研究所 | 基于全长转录组测序开发的姜荷花多态性ssr引物及试剂盒 |
| CN110029185A (zh) * | 2019-04-11 | 2019-07-19 | 广西壮族自治区亚热带作物研究所（广西亚热带农产品加工研究所） | 一种在转录组测序水平上进行剑麻叶纤维发育的研究方法 |
| CN111705160A (zh) * | 2020-07-31 | 2020-09-25 | 广西壮族自治区亚热带作物研究所（广西亚热带农产品加工研究所） | 龙舌兰麻cpSSR标记引物及其应用 |
| CN111705160B (zh) * | 2020-07-31 | 2023-04-07 | 广西壮族自治区亚热带作物研究所（广西亚热带农产品加工研究所） | 龙舌兰麻cpSSR标记引物及其应用 |
| CN113215220A (zh) * | 2021-05-31 | 2021-08-06 | 福建农林大学 | 一种基于转录组测序的橄榄ssr分子标记开发方法 |
| CN114898802A (zh) * | 2022-07-14 | 2022-08-12 | 臻和(北京)生物科技有限公司 | 基于血浆游离dna甲基化测序数据的末端序列频率分布特征确定方法、评价方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107201408B (zh) | 2020-10-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107201408A (zh) | 一种基于转录组测序开发剑麻ssr引物的方法 | |
| CN107345256B (zh) | 一种基于转录组测序开发山黧豆est-ssr引物组及方法和应用 | |
| CN106947827B (zh) | 一种获得鳙性别特异分子标记及其筛选方法和应用 | |
| CN112359093B (zh) | 血液中游离miRNA文库制备和表达定量的方法及试剂盒 | |
| CN103667275A (zh) | 油茶ssr分子标记 | |
| CN106754886A (zh) | 基于转录测序获得黑果枸杞ssr引物的方法 | |
| CN105238781A (zh) | 基于转录组序列开发的李ssr标记引物对及其应用 | |
| CN108103168A (zh) | 一种ffpe样本中dna质量的评估方法 | |
| CN107475449A (zh) | 一种适用于矮缩病毒科和双生病毒科病毒基因组拼接的转录组测序方法 | |
| CN108192893B (zh) | 基于转录组测序开发艾纳香ssr引物的方法 | |
| CN103013986B (zh) | 利用est序列的冗余性开发辣椒ssr标记及其方法 | |
| CN108998553A (zh) | 一种快速筛选多态性微卫星位点目标引物的方法及引物 | |
| CN104195245B (zh) | 基于转录组测序开发的红掌ssr引物对及试剂盒 | |
| CN110951911A (zh) | 基于转录组的椴树属est-ssr引物及其筛选方法和应用 | |
| CN104988148B (zh) | 一种沼泽型水牛ssr引物及其应用 | |
| CN111944922B (zh) | 基于芍药转录组序列开发的est-ssr引物组及其应用 | |
| CN104278021A (zh) | 一种铁皮石斛引物的制备方法 | |
| CN107828858A (zh) | 一种基于转录组测序开发鬼针草植物ssr引物的方法 | |
| CN112662806A (zh) | 一种钻喙兰属ssr分子标记引物组合物及其应用 | |
| KR20100061910A (ko) | 벼 품종을 구분하기 위한 특이 프라이머 및 이의 용도 | |
| CN106399475B (zh) | 一种获取rDNA ITS2序列用来对铁皮枫斗鉴定的方法 | |
| CN113151543B (zh) | 一种利用ssr标记快速鉴定马蹄莲种苗类型的引物组、方法、试剂盒及其应用 | |
| CN111235303B (zh) | 一种鉴别大米草和互花米草的方法 | |
| KR101807623B1 (ko) | 방풍, 식방풍 및 해방풍의 엽록체 게놈 및 핵 리보솜 dna 서열의 완전 해독 기반 종 식별 마커와 프라이머 세트 및 이의 용도 | |
| Coate | Beyond transcript concentrations: quantifying polyploid expression responses per biomass, per genome, and per cell with RNA-Seq |
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 |