CN103667452B - 一种利用ssr分子标记引物体系分析柑橘黄龙病菌遗传多样性的方法 - Google Patents
一种利用ssr分子标记引物体系分析柑橘黄龙病菌遗传多样性的方法 Download PDFInfo
- Publication number
- CN103667452B CN103667452B CN201310564327.7A CN201310564327A CN103667452B CN 103667452 B CN103667452 B CN 103667452B CN 201310564327 A CN201310564327 A CN 201310564327A CN 103667452 B CN103667452 B CN 103667452B
- Authority
- CN
- China
- Prior art keywords
- las
- ssr
- seq
- primer
- 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
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/6844—Nucleic acid amplification reactions
- C12Q1/6858—Allele-specific amplification
Landscapes
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Organic Chemistry (AREA)
- Engineering & Computer Science (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Immunology (AREA)
- Microbiology (AREA)
- Molecular Biology (AREA)
- Analytical Chemistry (AREA)
- Physics & Mathematics (AREA)
- Biotechnology (AREA)
- Biochemistry (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Genetics & Genomics (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明公开了一种利用SSR分子标记引物体系分析柑橘黄龙病菌遗传多样性的方法。所述SSR分子标记引物体系包括33对Las-SSR可用引物，其中18对为SSR核心引物。本发明首次提出了基于SSR分子标记并结合PAGE电泳分析柑橘黄龙病菌亚洲种（Ca.L.asiaticus）菌株多样性的方法，该方法简单、方便、分析效率高，为黄龙病分子流行学的研究奠定了分子基础。
Description
技术领域
本发明涉及分子生物学技术领域。更具体地，涉及一种利用SSR分子标记引物体系分析柑橘黄龙病菌遗传多样性的方法。
背景技术
柑橘黄龙病是危及世界和我国多个柑橘产区的最主要的病害。该病害最早在我国广东发现，至今已有近百年历史。目前，黄龙病在亚洲和非洲广泛分布，在大洋洲、南美洲和北美洲的部分柑橘主产区也严重发生，已报道在50多个国家和地区造成了巨大的损失。对黄龙病病原的认识经历了漫长的时间，目前普遍认为黄龙病是由限制于寄主植物韧皮部、不可体外培养的革兰氏阴性细菌（Candidatus Libriacter spp.）引起的。Ca.L. spp包含3个亚种：亚洲种（Ca. L. asiaticus）、非洲种（Ca.L. africanus）、和美洲种（Ca.L. americanus）。其中亚洲种（Ca.L. asiaticus）为世界范围内流行范围最广、致病力最强的种群。我国发现的所有黄龙病相关病原均属于亚洲种（Ca. L. asiaticus）。
目前已成功地完成了2个亚洲种（Ca.L. asiaticus）菌株的全基因组测序。2009年公布了一个从携毒柑橘木虱Psy62中获得的大小为1.23 Mb的亚洲种（Ca.L. asiaticus）全基因组序列（Genebank ID: gb∣CP001677.5）；2013年，Genebank上又公布了一个来自中国广西的柑橘黄龙病菌亚洲种（Ca. L. asiaticus）菌株的基因组（Genebank ID: gb∣CP004005.1），大小为1.27 Mb。在Psy62菌株的全基因组序列公布以后，Liu等分析了中国云南和广东两个省的柑橘黄龙病菌亚洲种（Ca. L. asiaticus）菌株中的前噬菌体特异引物766F/R的扩增结果，首次报道中国菌株的地区变异性（Liu et al., 2011）。Omp基因的限制性酶切多态性也证明了菌株的多态性与地域相关（Hu et al., 2011）。Wang等分析262个分别来自中国和美国的菌株某区域（CLIBASIA_ 05640- CLIBA SIA_05650）的电泳条带类型谱，发现2个国家之间及中国的云南和广东两省之间的菌株具有多态性，指出基因的插入和缺失突变是菌株遗传多样性的重要依据（Wang et al., 2012）。Zhou报道了基因组中的另两个高度变异的前噬菌体基因，并利用这两个基因分析来自不同国家的32个菌株，得出菌株多态性与地域差异相关的结论（Zhou et al., 2011）。
柑橘黄龙病菌株的地区差异性和遗传多样性一方面曾经引起黄龙病起源的争议，另一方面对于该病害在不同地域间的流行和传播也曾一度引起热议。柑橘生产实践者或研究者甚至发现柑橘黄龙病菌亚洲种在曾经流行过美洲种的地区大势流行，造成严重的危害。因此，基于病原分子多样性的病原流行趋势分析具有重要的现实意义。目前常用的遗传多样性分析都是基于常规PCR技术、琼脂糖凝胶电泳及序列测定来实现的。但是常规PCR通常设置较大的（20μL ~25μL）的体系，浪费大量的试剂和实验材料（尤其是各种PCR酶价格比较昂贵，另外需要从国外取得的实验材料更加珍贵）；并且普通琼脂糖凝胶电泳的分辨率比较低；另外，序列测定需要用到高精度的大型测序仪，批量测序价格昂贵。因此，寻找出一种简单、快速、经济、高效的鉴定柑橘黄龙病菌遗传多样性的方法具有非常重要的理论和实践意义。
简单串联重复序列（Simple sequence repeat，SSR）也叫微卫星标记（Microsatellite），是基于某一物种已知基因组的、以一定数目核苷酸（1~8个碱基对）2次以上重复排列的长达几十到几百个核苷酸的序列。SSR标记因重复性好、稳定性强、多态性高且在基因组中普遍存在的优点被广泛用于植物遗传图谱的构建、品种鉴定、遗传多样性分析和分子标记辅助育种等领域。虽然遗传标记在许多模式植物上已经广泛应用，但在植物病原细菌上还是方兴未艾。植物病原细菌的遗传方式完全不同于植物，细菌的遗传并不涉及菌株的杂交，另一方面，细菌核酸的获得需要克服人工培养的障碍，因此想要利用SSR分子标记的方法来分析植物病原细菌的遗传多样性，还需克服一系列的技术问题。
柑橘黄龙病菌目前还不能进行体外培养，无法提取纯病菌的核酸。目前还没有简单有效的遗传多样性分析方法。
发明内容
本发明的目的在于克服现有柑橘黄龙病菌遗传多样性分析的技术不足，提供一套用于柑橘黄龙病菌遗传多样性分析的SSR分子标记引物体系。
本发明另一目的在于提供一种利用上述SSR分子标记引物体系鉴定和分析柑橘黄龙病菌遗传多样性的方法。
本发明上述目的通过以下技术方案实现：
本发明提供了一种柑橘黄龙病菌遗传多样性分析SSR分子标记引物体系，所述SSR分子标记引物体系包括33对Las-SSR可用引物，分别为Las-ssr1~Las-ssr8、Las-ssr10~Las-ssr14、Las-ssr16~Las-ssr20、Las-ssr22~ Las-ssr36，引物序列如SEQ ID NO.1~16、SEQ ID NO.19~28、SEQ ID NO.31~40、SEQ ID NO.43~72所示；
其中有18对为SSR核心引物，分别是Las-ssr1、Las-ssr3、Las-ssr5、Las-ssr7、Las-ssr12~14、Las-ssr18、Las-ssr22、Las-ssr27~29、Las-ssr31~36，引物序列如SEQ ID NO.1~2、SEQ ID NO.5~6、SEQ ID NO.9~10、SEQ ID NO.13~14、SEQ ID NO.23~28、SEQ ID NO.35~36、SEQ ID NO.43~44、SEQ ID NO.53~58、SEQ ID NO.61~72所示。
所述SSR分子标记引物体系构建方法步骤如下：
S1.柑橘黄龙病菌亚洲种（Ca. L. asiaticus）全基因组序列的获得：从Genebank中下载得到2个柑橘黄龙病菌亚洲种（Ca. L. asiaticus）菌株序列，Genebank ID分别为gb∣CP001677.5和gb∣CP004005.1，分别为采自美国佛罗里达州Psy62菌株和中国广西gxpsy菌株的基因组序列；
S2.全基因组序列中SSR位点的鉴定：利用SSR检索软件Tandem Repeats Finder（version 4.07b）获得了36个含有SSR位点信息的序列；具体是利用Tandem Pepeats Finder (version 4.07b)在线软件扫描分析候选的SSR位点信息，将得到的结果进行筛选，筛选标准：（1）若匹配率（Percent Matches）为100%，则重复单元的重复长度为2 bp~10 bp；（2）若Percent Matches介于90%和100%之间，则重复单元的最小重复长度为10 bp；（3）重复数大于等于2；
S3.根据S2得到的含有SSR位点的序列，利用引物设计软件Primer Premier 5.0合成SSR引物，得到36对Las-SSR引物（如表1所示），引物设计参数主要包括：引物长度为15 bp ~25 bp，退火温度为45℃~65℃，GC含量为40%~65%，PCR扩增产物长度为100 bp ~600 bp，其中Las-ssr28位点的序列总长度大于600 bp，其预测的PCR扩增产物长度为1742 bp；
S4.用普通PCR的方法验证S3得到的36对Las-SSR引物的特异性：选取1个黄龙病柑橘叶片样品（阳性）和1个健康柑橘样品（阴性），剪取0.15 g柑橘叶片中脉，磨碎后用OMEGA公司的植物DNA提取试剂盒（产品号D2485-02）提取基因组总DNA。阳性样品中包含柑橘基因组DNA和柑橘黄龙病菌亚洲种（Ca. L. asiaticus）基因组DNA。利用1%（w/v）的琼脂糖凝胶检测提取的DNA的质量。以上述2个DNA样品为模板，分别用所有设计的SSR引物进行筛选；阳性样品能扩增出目的片段而阴性样品不能扩增出该片段的引物为可用引物。筛选得到33对具有高稳定性、高扩增效率和较强特异性的Las-SSR可用引物（如表1中所示的33对可用引物）；
S5.利用S4筛选到的33对Las-SSR可用引物对柑橘黄龙病菌亚洲种（Ca. L. asiaticus）进行遗传多样性分析：以柑橘黄龙病菌亚洲种（Ca. L. asiaticus）阳性DNA样品为模板进行PCR扩增；预测扩增片段小于500bp的PCR产物用聚丙烯酰胺凝胶电泳，预测片段大于1500bp的PCR产物以1% 琼脂糖凝胶进行电泳分析，并测序比对；筛选出多态性非常好的18对SSR核心引物，其中17对SSR核心引物的PCR扩增片段均小于600bp，用聚丙烯酰胺凝胶电泳即可分辨。只有Las-ssr28预测的PCR扩增产物长度为1742 bp，还需使用琼脂糖凝胶电泳来分析。
步骤S3所述36对Las-SSR引物为Las-ssr1~Las-ssr36，序列如SEQ ID NO.1~ 72所示。
步骤S4所述33对Las-SSR可用引物分别为Las-ssr1~Las-ssr8、Las-ssr10~Las-ssr14、Las-ssr16~Las-ssr20、Las-ssr22~ Las-ssr36，引物序列如SEQ ID NO.1~16、SEQ ID NO.19~28、SEQ ID NO.31~40、SEQ ID NO.43~72所示。
步骤S5所述18对SSR核心引物分别为Las-ssr1、Las-ssr3、Las-ssr5、Las-ssr7、Las-ssr12~14、Las-ssr18、Las-ssr22、Las-ssr27~29、Las-ssr31~36，引物序列如SEQ ID NO.1~2、SEQ ID NO.5~6、SEQ ID NO.9~10、SEQ ID NO.13~14、SEQ ID NO.23~28、SEQ ID NO.35~36、SEQ ID NO.43~44、SEQ ID NO.53~58、SEQ ID NO.61~72所示。
本发明还提供了一种上述SSR分子标记引物体系鉴定和分析柑橘黄龙病菌遗传多样性的方法，其具体步骤是：
S1.实验材料的准备：选取典型斑驳症状的柑橘叶片，提取柑橘样品叶中脉总DNA，采用国际上通用的柑橘黄龙病菌特异引物OI1和OI2c扩增16S rDNA片段，检测为阳性的DNA样品用于SSR遗传多样性分析；
S2.以S1获得的阳性的DNA为模板，利用本发明的SSR分子标记引物体系进行PCR扩增；
S3.对S2所述PCR扩增得到的PCR产物进行电泳、显色，显色结束后，拍照获取谱带信息，比较同一SSR标记在不同菌株中的多态性；
S4.结果统计及遗传多样性分析：读取S3获得的谱带信息，以条带位置为依据，最小片段的带型为“1”，依次递增；收集所有条带类型信息，使用分析软件分别计算SSR标记的多态性位点频率（PPB）、Nei's基因多样性指数（h）和Shannon's信息指数（I），反映菌株的遗传多样性；采用相关软件方法进行分析，获得UPGMA聚类图。
其中，步骤S2所述PCR扩增的体系为15μL：包括9.9μL灭菌ddH2O，1.5μL 含Mg2+的10×Taq buffer、1.5μL的2.5 mM的dNTPs，0.8μL的10μM正反向引物混合液，20 ng样品DNA模板及0.2μL的2.5 M/μL Taq DNA 聚合酶；
扩增程序为：94℃预变性5 min；94℃变性30s、56℃退火40s、72℃延伸45s，循环35次；72℃延伸10min；
步骤S3所述电泳优选为聚丙烯酰胺凝胶电泳（PAGE电泳），所述显色为银染显色；
步骤S4所述分析软件优选为POPGENE1.32软件，所述相关软件方法优选为PowerMark V3.0软件中的N-J聚类方法。
本发明具有以下有益效果：
本发明提供了柑橘黄龙病菌亚洲种（Ca. L. asiaticus）的全套SSR分子标记引物体系并筛选出核心引物，结合PAGE电泳的方法为该病原的遗传多样性研究开辟了新的途径。与现有技术一次只能分析一个遗传位点多样性，本发明一次性涉及全基因组上分布的33个位点，效率大大提高，且足以对菌株遗传多样性进行批量区分，方法简单，步骤少，易于操作，也克服了批量测序价格昂贵的问题。另外，较小的PCR体系（10μL~15μL）就能实现检测目的，节省成本和实验材料，所用的材料均为常规试剂或装置，耗资少；数据的统计分析借助于已开发的软件，可一次性处理所有数据。同时，首次结合PAGE电泳的方法克服了琼脂糖凝胶电泳分辨率低的问题。
另外，本发明采用提取感染黄龙病的叶片总DNA的方法，一并提取得到黄龙病菌DNA，设计的黄龙病菌特异性引物能避免柑橘基因组DNA的干扰，克服了黄龙病菌无法获得纯培养等技术缺陷，具有很好的推广应用价值。
附图说明
图1 为所有36对Las-SSR可用引物的常规PCR筛选结果；其中“+”表示携带柑橘黄龙病菌亚洲种（Ca. L. asiaticus）的砂糖橘叶中脉DNA，“-”表示健康砂糖橘叶中脉DNA。
图2为9对Las-SSR核心引物扩增来自不同地区的32个柑橘黄龙病菌亚洲种（Ca. L. asiaticus）菌株的PAGE电泳多态图。
图3为9对Las-SSR核心引物扩增来自不同地区的32个柑橘黄龙病菌亚洲种（Ca. L. asiaticus）菌株的PAGE电泳多态图。
图4为我国6省的32个柑橘黄龙病菌亚洲种（Ca. L. asiaticus）菌株基于SSR分析的条带多样性的聚类图。
具体实施方式
以下结合附图和具体实施例来进一步说明本发明，但实施例并不对本发明做任何形式的限定。除非特别说明，本发明采用的试剂、设备及方法为本技术领域常规试剂、设备及方法。
实施例1 本发明Las-SSR分子标记引物体系的筛选
1、利用Genebank数据库中的Genome数据库，查找得到gb∣CP001677.5和gb∣CP004005.1两个柑橘黄龙病菌亚洲种（Ca. L. asiaticus）全基因组序列，核苷酸全长分别为1.23 Mb和1.27 Mb。
本实施例所用分析软件为SSR信息查找软件Tandem Repeats Finder（version 4.07b）在线分析系统（http://tandem.bu.edu/trf/trf.html）和引物合成软件Primer Premier 5.0。
实验方法为常规PCR、琼脂糖凝胶电泳、聚丙烯酰胺凝胶电泳（PAGE）。
2、利用SSR检索软件Tandem Repeats Finder（version 4.07b）查找2个基因组中的所有2-500个碱基重复基序的所有SSR位点信息，经筛选后得到36个SSR位点，筛选标准：（1）若匹配率（Percent Matches）为100%，则重复单元的重复长度为4 bp-10 bp；（2）若Percent Matches介于90%和100%之间，则重复单元的最小重复长度为10 bp；（3）重复数大于等于2。
用Primer Premier 5.0软件对这36个位点设计相应的SSR引物，设计得到36对Las-SSR引物（如表1所示）。引物设计参数主要包括：引物长度为15 bp ~25 bp，退火温度为45℃~65℃，GC含量为40%~65%，PCR扩增产物长度为100 bp ~600 bp，其中Las-ssr28位点的序列总长度大于600 bp，其预测的PCR扩增产物长度为1742 bp。
表1本发明中设计合成的柑橘黄龙病菌Las-SSR引物序列
提取1个感染黄龙病菌的砂糖橘叶中脉样品（携带柑橘黄龙病菌亚洲种Ca. L. asiaticus，阳性）和1个健康砂糖橘叶中脉样品（阴性）的基因组DNA，以这2份DNA样品为模板，用上述36对引物进行PCR扩增。
扩增体系15μL：9.9μL灭菌ddH2O，10×Taq buffer（含Mg2+，TianGen）、dNTPs（2.5 mM，TianGen）各1.5μL，0.8μL正反向引物（10μM）混合液，约20 ng样品DNA模板及0.2μL Taq DNA 聚合酶（2.5 M/μL）；PCR扩增程序为：（1）94℃预变性5 min；（2）94℃变性30s、56℃退火40s、72℃延伸45s，该步骤（2）重复35个循环；（3）72℃延伸10min。
PCR产物用1%（w/v）的琼脂糖凝胶进行电泳检测扩增的效果。其中33对引物具有高扩增稳定性、高扩增效率和较强特异性（如附图1所示）。
实施例2 Las-SSR分子标记核心引物的获得及柑橘黄龙病菌遗传多态性分析
选取来自我国福建、广东、广西、江西、云南和浙江6个主要柑橘产区的感染柑橘黄龙病菌的柑橘叶片样品，对每份样品进行拍照记录。利用OMEGA公司E.Z.N.A.TM植物基因组DNA试剂盒（HP plant DNA kit (200), D2485-02, OMEGA bio-tek）提取提取基因组DNA，剪取叶中脉0.15 g置于2.0 mL离心管中，加入500 mL提取液CPL和2颗磁珠，置于MP FastPrep-24磨样机上打碎。用1%（w/v）的琼脂糖凝胶电泳检测基因组DNA的质量。
利用柑橘黄龙病菌亚洲种（Ca. L. asiaticus）16s rDNA的特异引物OI1/OI2c（OI1: 5’- GCGCGTATCCAATACGAGCGGCA-3’；OI2c: 5’- GCCTCGCGACTTCGCAACCCAT-3’）对提取的样品进行分析，筛选出32个阳性样品，即为32个柑橘黄龙病菌亚洲种（Ca. L. asiaticus）菌株（如表2所示）。
表2 本发明中用于引物筛选及遗传多样性分析的柑橘黄龙病菌亚洲种（Ca. L. asiaticus）菌株信息
用实施例一中筛选出来的33对可用引物分别对这32个菌株进行扩增。扩增体系15μL：9.9μL灭菌ddH2O，10×Taq buffer（含Mg2+，TianGen）、dNTPs（2.5 mM，TianGen）各1.5μL，0.8μL正反向引物（10μM）混合液，约20 ng样品DNA模板及0.2μL Taq DNA 聚合酶（2.5 M/μL）；PCR扩增程序为：（1）94℃预变性5 min；（2）94℃变性30s、56℃退火40s、72℃延伸45s，该步骤（2）重复35个循环；（3）72℃延伸10min。
上述PCR产物中加入3×loading buffer缓冲液后，取1μL用6% PAGE胶在110 V电压条件下垂直电泳400 min。
所述6% PAGE胶配制方法：称取9.6 g尿素，溶于35 mL ddH2O，再分别加入16 mL 5×TBE缓冲液（TBE缓冲液配方为54gTris、27.5g硼酸、205g无水NaAc及20mL 1 M Na2EDTA (Ph8.0)，加水定容到1000mL）、16 mL 40%的丙烯酰胺（丙烯酰胺 380.0g，甲叉双丙烯酰胺 20.0 g，加ddH2O定容至1000 mL）、1.0 mL 10%的过硫酰胺（w/v）和65 μL的四甲基乙二胺溶液，摇匀即可灌胶。
垂直电泳跑完之后，取下PAGE胶，经漂洗、银染、显色和漂洗后拍照，收集条带信息。
所述银染的方法为：跑完电泳之后，先用超纯水将胶漂洗2次，每次60 s；之后加入100 mL 1%（w/v）的AgNO3溶液，置于摇床轻摇10min进行银染。
所述显色的方法为：银染后倒掉AgNO3溶液，用超纯水清洗1.5 min；加入显影液（12 g 氢氧化钠，0.15 g 四硼酸钠溶于600 mL ddH2O中，完全溶解后再加3.2 mL甲醛，用ddH2O定容至800 mL，摇匀即可）置于摇床轻摇显色至DNA条带清晰，倒掉显影液，加入超纯水清洗。
显色结束后，将胶置于白瓷盘中，用普通数码相机拍照获取谱带信息；比较同一SSR标记在不同菌株中的多态性。具有多态性且扩增稳定性和特异性强的SSR标记被确认为核心SSR标记。
对32个柑橘黄龙病菌亚洲种（Ca. L. asiaticus）菌株进行PCR扩增结果显示，其中18对引物可稳定扩增出2种或2种以上类型的条带，多态性较好，定为核心引物。所有33对可用引物中，分别有14对、3对和3对引物能扩增出1个、2个和3个条带；分别有8对、3对、2对和1对引物可扩增出1-2个、1-3个、1-4个和1-5个条带。可见18对核心引物的丰度足以满足种间不同样品的遗传多样性分析（核心引物的扩增结果如附图2和图3所示）。
实施例3 中国6省的柑橘黄龙病菌亚洲种（Ca. L. asiaticus）菌株多样性初步分析
1、分析方法
方法同实施例2。
2、结果分析
收集来自中国6省32个菌株的所有32个（不含Las-ssr28）可用SSR引物的条带信息，用不同的阿拉伯数字表示不同的条带类型。
用POPGENE1.32软件分别获得SSR标记的多态性位点频率（PPB）、Nei's基因多样性指数（h）和Shannon's信息指数（I）（如表3和表4所示），反映32个菌株的遗传多样性。
表3 其中32对可用SSR标记的在32个菌株中的多态性位点频率(PPB)
表4 其中32对可用SSR标记的多态性位点Nei's基因多样性指数(h)和Shannon's信息指数(I)
Sample Size Mean = 30. h* Mean = 0.2366 I* Mean = 0.4253 h*St.Dev = 0.2577 I*St.Dev = 0.5002
h=Nei's-1973 gene diversity
I=Shannon's Information index [Lewontin (1972)]
结果表明：Shannon’s信息指数为0.1391-1.8979，平均为0.6805；Nei's基因多样性指数为0.0605-0.8145，平均为0.3785。
同时通过Popgene 1.32计算不同省份菌株之间的遗传距离及遗传相似性，结果表明，不同省（市）之间的菌株遗传相似性为0.4839-0.9677，平均为0.7579；同省菌株之间遗传相似性分别为：江西（0.8137-0.9655），广西（0.6985-0.9677），福建（0.6237-0.9667），云南（0.6433-0.9032），广东（0.8065-0.9525），但其中浙江菌株数低于3株，数据未列出。结果可见，来自福建、广东、广西、江西、云南和浙江各省内部的柑橘黄龙病菌亚洲种（Ca. L. asiaticus）菌株间具有较为丰富的多样性。
以PowerMark V3.0软件采用类平均数聚类方法对32个菌株进行分析，获得N-J聚类图（如附图4所示）。可见这32个菌株可分为5组：I群包含2个省份的4株菌株：福建2株和江西2株，该2省为邻省；II群包含3个省份的6株菌株，其中福建的3株和广西的2株均采自柚子品种；III群包含3个省份的11株菌株，其中包括来自为邻省的广西和云南的分别3株和7株；IV群只包含了广东省的6株菌株；V群包含了2个邻近省份的5株菌株，其中江西3株，浙江2株。以上结果初步表明柑橘黄龙病菌亚洲种（Ca. L. asiaticus）群体与地理分布和品种相关。
SEQUENCE LISTING
<110> 华南农业大学
<120> 一种利用SSR分子标记引物体系分析柑橘黄龙病菌遗传多样性的方法
<130>
<160> 72
<170> PatentIn version 3.3
<210> 1
<211> 22
<212> DNA
<213> Las-ssr1上游引物
<400> 1
tcgcttgaaa gttattattg gg 22
<210> 2
<211> 22
<212> DNA
<213> Las-ssr1下游引物
<400> 2
gcataagaat agaagaaccc ta 22
<210> 3
<211> 17
<212> DNA
<213> Las-ssr2上游引物
<400> 3
tttccaagtc atttccg 17
<210> 4
<211> 17
<212> DNA
<213> Las-ssr2下游引物
<400> 4
caccgaccaa tacaacg 17
<210> 5
<211> 17
<212> DNA
<213> Las-ssr3上游引物
<400> 5
attcatggga agcgtaa 17
<210> 6
<211> 17
<212> DNA
<213> Las-ssr3下游引物
<400> 6
atgctcagaa cactgca 17
<210> 7
<211> 16
<212> DNA
<213> Las-ssr4上游引物
<400> 7
ttttgggaag agttgc 16
<210> 8
<211> 19
<212> DNA
<213> Las-ssr4下游引物
<400> 8
tctttcgtta tttctggta 19
<210> 9
<211> 17
<212> DNA
<213> Las-ssr5上游引物
<400> 9
acagggacgg tgaaaca 17
<210> 10
<211> 17
<212> DNA
<213> Las-ssr5下游引物
<400> 10
tgccttgacg caatcta 17
<210> 11
<211> 16
<212> DNA
<213> Las-ssr6上游引物
<400> 11
gatatgggta tgctcg 16
<210> 12
<211> 17
<212> DNA
<213> Las-ssr6下游引物
<400> 12
tagtgtttct tgtgctg 17
<210> 13
<211> 16
<212> DNA
<213> Las-ssr7上游引物
<400> 13
cgcctataaa tccctt 16
<210> 14
<211> 15
<212> DNA
<213> Las-ssr7下游引物
<400> 14
ctgttgcacc cgact 15
<210> 15
<211> 19
<212> DNA
<213> Las-ssr8上游引物
<400> 15
aagggattct tactatgtg 19
<210> 16
<211> 15
<212> DNA
<213> Las-ssr8下游引物
<400> 16
tcgtgctgct ttgtt 15
<210> 17
<211> 16
<212> DNA
<213> Las-ssr9上游引物
<400> 17
acgatttaca agttca 16
<210> 18
<211> 15
<212> DNA
<213> Las-ssr9下游引物
<400> 18
attattggga gggaa 15
<210> 19
<211> 16
<212> DNA
<213> Las-ssr10上游引物
<400> 19
tgtacgttcg gagaag 16
<210> 20
<211> 16
<212> DNA
<213> Las-ssr10下游引物
<400> 20
ggtctgtatc tggtgc 16
<210> 21
<211> 17
<212> DNA
<213> Las-ssr11上游引物
<400> 21
gtttcacagt ccaaaat 17
<210> 22
<211> 15
<212> DNA
<213> Las-ssr11下游引物
<400> 22
aacgcacaga taagc 15
<210> 23
<211> 16
<212> DNA
<213> Las-ssr12上游引物
<400> 23
gttgcttcgt ttatcc 16
<210> 24
<211> 16
<212> DNA
<213> Las-ssr12下游引物
<400> 24
gcattctgtg cctctt 16
<210> 25
<211> 17
<212> DNA
<213> Las-ssr13上游引物
<400> 25
gtaggagtcc ccgaaat 17
<210> 26
<211> 17
<212> DNA
<213> Las-ssr13下游引物
<400> 26
gcctgtacga ggtttga 17
<210> 27
<211> 17
<212> DNA
<213> Las-ssr14上游引物
<400> 27
ttttggagga actaacg 17
<210> 28
<211> 18
<212> DNA
<213> Las-ssr14下游引物
<400> 28
tcacatgaat atgaaccc 18
<210> 29
<211> 17
<212> DNA
<213> Las-ssr15上游引物
<400> 29
tgggagattc gtgatgt 17
<210> 30
<211> 17
<212> DNA
<213> Las-ssr15下游引物
<400> 30
atagccgttg gaggtag 17
<210> 31
<211> 19
<212> DNA
<213> Las-ssr16上游引物
<400> 31
gtactatgtt tagcctgta 19
<210> 32
<211> 18
<212> DNA
<213> Las-ssr16下游引物
<400> 32
cttctaatca atagaccc 18
<210> 33
<211> 15
<212> DNA
<213> Las-ssr17上游引物
<400> 33
aggggtgttt ctgtc 15
<210> 34
<211> 16
<212> DNA
<213> Las-ssr17下游引物
<400> 34
ataaaatcaa tgtggc 16
<210> 35
<211> 15
<212> DNA
<213> Las-ssr18上游引物
<400> 35
cagatgccga ccttc 15
<210> 36
<211> 18
<212> DNA
<213> Las-ssr18下游引物
<400> 36
gtgagctgta gtccctat 18
<210> 37
<211> 16
<212> DNA
<213> Las-ssr19上游引物
<400> 37
ccataggttg ggaatc 16
<210> 38
<211> 17
<212> DNA
<213> Las-ssr19下游引物
<400> 38
tgccgaataa gtcgttt 17
<210> 39
<211> 15
<212> DNA
<213> Las-ssr20上游引物
<400> 39
ccgtcctgat gtggt 15
<210> 40
<211> 18
<212> DNA
<213> Las-ssr20下游引物
<400> 40
ttgatcttaa tcggtctc 18
<210> 41
<211> 20
<212> DNA
<213> Las-ssr21上游引物
<400> 41
aataagtatt aaaaggtcgg 20
<210> 42
<211> 18
<212> DNA
<213> Las-ssr21下游引物
<400> 42
gctaaggagg gtacagag 18
<210> 43
<211> 16
<212> DNA
<213> Las-ssr22上游引物
<400> 43
gtaaatgatc caccac 16
<210> 44
<211> 17
<212> DNA
<213> Las-ssr22下游引物
<400> 44
caaaaccaag atagact 17
<210> 45
<211> 18
<212> DNA
<213> Las-ssr23上游引物
<400> 45
atggcagacg gtgaaatt 18
<210> 46
<211> 18
<212> DNA
<213> Las-ssr23下游引物
<400> 46
tatcctcaag gcaaacaa 18
<210> 47
<211> 18
<212> DNA
<213> Las-ssr24上游引物
<400> 47
tgcgtcagtt tctcgttt 18
<210> 48
<211> 18
<212> DNA
<213> Las-ssr24下游引物
<400> 48
cgttccaggt gttcattc 18
<210> 49
<211> 15
<212> DNA
<213> Las-ssr25上游引物
<400> 49
tcagttcgtc cttgc 15
<210> 50
<211> 17
<212> DNA
<213> Las-ssr25下游引物
<400> 50
tcatacctat tcaccgt 17
<210> 51
<211> 15
<212> DNA
<213> Las-ssr26上游引物
<400> 51
ggcaagccat acaaa 15
<210> 52
<211> 16
<212> DNA
<213> Las-ssr26下游引物
<400> 52
aagaaagcgt aagcaa 16
<210> 53
<211> 15
<212> DNA
<213> Las-ssr27上游引物
<400> 53
tgccaccgaa tagaa 15
<210> 54
<211> 20
<212> DNA
<213> Las-ssr27下游引物
<400> 54
tacatcgtgt aatagatgaa 20
<210> 55
<211> 19
<212> DNA
<213> Las-ssr28上游引物
<400> 55
catccttgtc gtcatagaa 19
<210> 56
<211> 19
<212> DNA
<213> Las-ssr28下游引物
<400> 56
ggaagtaaat caggtagcc 19
<210> 57
<211> 16
<212> DNA
<213> Las-ssr29上游引物
<400> 57
tccaagctgt tatcca 16
<210> 58
<211> 16
<212> DNA
<213> Las-ssr29下游引物
<400> 58
gcgaaatcag tatgcc 16
<210> 59
<211> 18
<212> DNA
<213> Las-ssr30上游引物
<400> 59
tttcctttgg tttcttct 18
<210> 60
<211> 18
<212> DNA
<213> Las-ssr30下游引物
<400> 60
caccaccaca aacattag 18
<210> 61
<211> 18
<212> DNA
<213> Las-ssr31上游引物
<400> 61
atgagtcggg acggaagg 18
<210> 62
<211> 15
<212> DNA
<213> Las-ssr31下游引物
<400> 62
cgggctaacg cacaa 15
<210> 63
<211> 17
<212> DNA
<213> Las-ssr32上游引物
<400> 63
tgataacaag cacctct 17
<210> 64
<211> 17
<212> DNA
<213> Las-ssr32下游引物
<400> 64
ctaatctgaa cccaaac 17
<210> 65
<211> 17
<212> DNA
<213> Las-ssr33上游引物
<400> 65
tttgtcctgt attcgta 17
<210> 66
<211> 16
<212> DNA
<213> Las-ssr33下游引物
<400> 66
actcccctct tatcaa 16
<210> 67
<211> 18
<212> DNA
<213> Las-ssr34上游引物
<400> 67
acataaaccc aaactagc 18
<210> 68
<211> 18
<212> DNA
<213> Las-ssr34下游引物
<400> 68
ggcattgttc ttgaccgt 18
<210> 69
<211> 17
<212> DNA
<213> Las-ssr35上游引物
<400> 69
gagacggaca ctcaaca 17
<210> 70
<211> 16
<212> DNA
<213> Las-ssr35下游引物
<400> 70
gcaagggatt tatgct 16
<210> 71
<211> 18
<212> DNA
<213> Las-ssr36上游引物
<400> 71
agatattccc aaagcacg 18
<210> 72
<211> 18
<212> DNA
<213> Las-ssr36下游引物
<400> 72
atttggtttg tccgaagg 18
Claims (5)
1.一种用于柑橘黄龙病菌遗传多样性分析的SSR分子标记引物体系，其特征在于，包括33对Las-SSR可用引物，分别为Las-ssr1~Las-ssr8、Las-ssr10~Las-ssr14、Las-ssr16~Las-ssr20、Las-ssr22~ Las-ssr36，引物序列如SEQ ID NO.1~16、SEQ ID NO.19~28、SEQ ID NO.31~40、SEQ ID NO.43~72所示。
2.根据权利要求1所述SSR分子标记引物体系，其特征在于，所述33对Las-SSR可用引物中有18对为SSR核心引物，分别是Las-ssr1、Las-ssr3、Las-ssr5、Las-ssr7、Las-ssr12~14、Las-ssr18、Las-ssr22、Las-ssr27~29、Las-ssr31~36，引物序列如SEQ ID NO.1~2、SEQ ID NO.5~6、SEQ ID NO.9~10、SEQ ID NO.13~14、SEQ ID NO.23~28、SEQ ID NO.35~36、SEQ ID NO.43~44、SEQ ID NO.53~58、SEQ ID NO.61~72所示。
3.一种权利要求1或2所述SSR分子标记引物体系的构建方法，其特征在于，包括如下步骤：
S1. 柑橘黄龙病菌亚洲种全基因组序列的获得：从Genebank中下载得到2个柑橘黄龙病菌亚洲种菌株序列，Genebank ID: gb∣CP001677.5和gb∣CP004005.1；
S2.全基因组序列中SSR位点的鉴定：利用SSR检索软件获得了36个含有SSR位点信息的序列；
S3.根据S2得到的36个含有SSR位点的序列设计得到36对Las-SSR引物；
S4.用普通PCR的方法验证S3得到引物的特异性，筛选得到33对Las-SSR可用引物；
S5.利用S4筛选到的33对Las-SSR可用引物对柑橘黄龙病菌亚洲种进行遗传多样性分析，筛选得到18对SSR核心引物；
其中，步骤S3所述36对Las-SSR引物为Las-ssr1~ Las-ssr36，序列如SEQ ID NO.1~72所示；
步骤S4所述33对Las-SSR可用引物分别为Las-ssr1~Las-ssr8、Las-ssr10~Las-ssr14、Las-ssr16~Las-ssr20、Las-ssr22~ Las-ssr36，引物序列如SEQ ID NO.1~16、SEQ ID NO.19~28、SEQ ID NO.31~40、SEQ ID NO.43~72所示；
步骤S5所述18对SSR核心引物分别为Las-ssr1、Las-ssr3、Las-ssr5、Las-ssr7、Las-ssr12~14、Las-ssr18、Las-ssr22、Las-ssr27~29、Las-ssr31~36，引物序列如SEQ ID NO.1~2、SEQ ID NO.5~6、SEQ ID NO.9~10、SEQ ID NO.13~14、SEQ ID NO.23~28、SEQ ID NO.35~36、SEQ ID NO.43~44、SEQ ID NO.53~58、SEQ ID NO.61~72所示。
4.一种利用权利要求1或2所述SSR分子标记引物体系分析柑橘黄龙病菌遗传多样性的方法，其特征在于，包括如下步骤：
S1.实验材料的准备：选取典型斑驳症状的柑橘叶片，提取柑橘样品叶中脉总DNA，采用国际上通用的柑橘黄龙病菌特异引物OI1和OI2c扩增16S rDNA片段，检测为阳性的DNA样品进行SSR遗传多样性分析；
S2.以S1获得的阳性的DNA为模板，利用权利要求1或2所述的SSR分子标记引物体系进行PCR扩增；
S3.对S2所述PCR扩增得到的PCR产物进行电泳、显色，显色结束后，拍照获取谱带信息，比较同一SSR标记在不同菌株中的多态性；
S4.结果统计及遗传多样性分析：读取S3获得的谱带信息，以条带位置为依据，最小片段的带型为“1”，依次递增；收集所有条带类型信息，使用分析软件分别计算SSR标记的多态性位点频率、Nei's基因多样性指数和Shannon's信息指数，反映菌株的遗传多样性；采用相关软件方法进行分析，获得UPGMA聚类图；
其中，所述的柑橘黄龙病菌为柑橘黄龙病菌亚洲种；
步骤S1所述引物OI1的核苷酸序列为：GCGCGTATCCAATACGAGCGGCA；
OI2c的核苷酸序列为 GCCTCGCGACTTCGCAACCCAT。
5.根据权利要求4所述方法，其特征在于，步骤S3所述电泳为聚丙烯酰胺凝胶电泳。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310564327.7A CN103667452B (zh) | 2013-11-14 | 2013-11-14 | 一种利用ssr分子标记引物体系分析柑橘黄龙病菌遗传多样性的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310564327.7A CN103667452B (zh) | 2013-11-14 | 2013-11-14 | 一种利用ssr分子标记引物体系分析柑橘黄龙病菌遗传多样性的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103667452A CN103667452A (zh) | 2014-03-26 |
| CN103667452B true CN103667452B (zh) | 2014-12-17 |
Family
ID=50306220
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310564327.7A Active CN103667452B (zh) | 2013-11-14 | 2013-11-14 | 一种利用ssr分子标记引物体系分析柑橘黄龙病菌遗传多样性的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103667452B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106282365B (zh) * | 2016-09-09 | 2019-05-31 | 广州大学 | 一种快速对光倒刺鲃不同群体进行分类的试剂盒 |
| CN116083610B (zh) * | 2022-12-15 | 2023-11-07 | 华南农业大学 | 一种基于噬菌体基因为靶标的柑橘黄龙病菌检测引物及其应用 |
-
2013
- 2013-11-14 CN CN201310564327.7A patent/CN103667452B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN103667452A (zh) | 2014-03-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Hassel et al. | Restricted variation in plant barcoding markers limits identification in closely related bryophyte species | |
| CN101851675B (zh) | 检测检疫性实蝇的探针组及检测方法 | |
| CN106701914A (zh) | 一种基于dna条形码的细菌核酸测序鉴定方法 | |
| CN103320521B (zh) | 一种真核浮游植物多样性的快速高通量检测方法 | |
| CN102978198A (zh) | 用于评价动物肠道微生物群落多样性的微生物基因组dna间接提取方法 | |
| CN103981256B (zh) | 一种沙门氏菌crispr分型方法 | |
| CN101845495A (zh) | 检测检疫性实蝇的基因芯片 | |
| CN101712953A (zh) | 用于评价动物肠道微生物群落多样性的dna提取方法 | |
| CN110951911B (zh) | 基于转录组的椴树属est-ssr引物及其筛选方法和应用 | |
| CN103667452B (zh) | 一种利用ssr分子标记引物体系分析柑橘黄龙病菌遗传多样性的方法 | |
| Scali et al. | Vineyards genetic monitoring and Vernaccia di San Gimignano wine molecular fingerprinting | |
| CN103320425A (zh) | 一种用于大规模、高通量基因分型的贝类dna快速提取方法 | |
| CN113186315B (zh) | 一种用于检测水稻细菌性条斑病菌的引物对及检测方法 | |
| CN110970091A (zh) | 标签质控的方法及装置 | |
| CN105567800A (zh) | 一种真菌种属pcr鉴定方法 | |
| Séjalon-Delmas et al. | Molecular tools for the identification of Tuber melanosporum in agroindustry | |
| Kang et al. | High‐molecular‐weight DNA extraction for long‐read sequencing of plant genomes: An optimization of standard methods | |
| CN111554349B (zh) | 一种基于高通量测序的物种鉴定系统和方法 | |
| CN107937395A (zh) | 一种远海梭子蟹多态性微卫星分子标记及鉴定方法与应用 | |
| CN107862177A (zh) | 一种区分鲤群体的单核苷酸多态性分子标记集的构建方法 | |
| CN102534007B (zh) | 温室白粉虱特异性scar标记、特异性引物及快速分子鉴定方法 | |
| CN105039517A (zh) | 一种鉴别5种啮齿动物螺杆菌的pcr-hrm引物和方法 | |
| CN105200051A (zh) | 一种皱纹盘鲍中国和日本群体鉴别用snp标记 | |
| CN107760798A (zh) | 葡萄果实无核主效qtl位点的ssr分子标记 | |
| CN101215606A (zh) | 分析浮游生物群落dna多态性的方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20211215 Address after: Room 503, building 3, No. 133, development avenue, Tongxiang Economic Development Zone, Tongxiang City, Jiaxing City, Zhejiang Province Patentee after: ZHEJIANG MAIZHI NETWORK TECHNOLOGY CO.,LTD. Address before: 510642 No. five, 483 mountain road, Guangzhou, Guangdong, Tianhe District Patentee before: SOUTH CHINA AGRICULTURAL University |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20221219 Address after: 314500 No.1 Gongnong Road, South Gate, Chongfu Town, Tongxiang City, Jiaxing City, Zhejiang Province Patentee after: Jiaxing Zhuoshi Biotechnology Co.,Ltd. Address before: Room 503, building 3, No. 133, development avenue, Tongxiang Economic Development Zone, Tongxiang City, Jiaxing City, Zhejiang Province Patentee before: ZHEJIANG MAIZHI NETWORK TECHNOLOGY CO.,LTD. |