CN112877464A - 一种基于转录组中相邻微卫星序列的植物鉴别方法 - Google Patents
一种基于转录组中相邻微卫星序列的植物鉴别方法 Download PDFInfo
- Publication number
- CN112877464A CN112877464A CN202110265657.0A CN202110265657A CN112877464A CN 112877464 A CN112877464 A CN 112877464A CN 202110265657 A CN202110265657 A CN 202110265657A CN 112877464 A CN112877464 A CN 112877464A
- Authority
- CN
- China
- Prior art keywords
- sequence
- microsatellite
- sequences
- data
- splicing
- 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.)
- Pending
Links
- 108091092878 Microsatellite Proteins 0.000 title claims abstract description 77
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000012163 sequencing technique Methods 0.000 claims abstract description 15
- 108091081062 Repeated sequence (DNA) Proteins 0.000 claims abstract description 7
- 238000001914 filtration Methods 0.000 claims abstract description 6
- 238000012165 high-throughput sequencing Methods 0.000 claims abstract description 4
- 230000004069 differentiation Effects 0.000 claims description 4
- 230000000877 morphologic effect Effects 0.000 abstract description 5
- 108020004414 DNA Proteins 0.000 abstract description 2
- 241000196324 Embryophyta Species 0.000 description 33
- 108091032973 (ribonucleotides)n+m Proteins 0.000 description 19
- 230000002068 genetic effect Effects 0.000 description 7
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 6
- 239000011324 bead Substances 0.000 description 6
- YGSDEFSMJLZEOE-UHFFFAOYSA-N salicylic acid Chemical compound OC(=O)C1=CC=CC=C1O YGSDEFSMJLZEOE-UHFFFAOYSA-N 0.000 description 6
- 210000001519 tissue Anatomy 0.000 description 5
- 238000011282 treatment Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 239000003153 chemical reaction reagent Substances 0.000 description 4
- 239000002299 complementary DNA Substances 0.000 description 4
- 239000012634 fragment Substances 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- 238000002360 preparation method Methods 0.000 description 4
- GEWDNTWNSAZUDX-WQMVXFAESA-N (-)-methyl jasmonate Chemical compound CC\C=C/C[C@@H]1[C@@H](CC(=O)OC)CCC1=O GEWDNTWNSAZUDX-WQMVXFAESA-N 0.000 description 3
- 239000007788 liquid Substances 0.000 description 3
- GEWDNTWNSAZUDX-UHFFFAOYSA-N methyl 7-epi-jasmonate Natural products CCC=CCC1C(CC(=O)OC)CCC1=O GEWDNTWNSAZUDX-UHFFFAOYSA-N 0.000 description 3
- 229910052757 nitrogen Inorganic materials 0.000 description 3
- FJKROLUGYXJWQN-UHFFFAOYSA-N papa-hydroxy-benzoic acid Natural products OC(=O)C1=CC=C(O)C=C1 FJKROLUGYXJWQN-UHFFFAOYSA-N 0.000 description 3
- 229960004889 salicylic acid Drugs 0.000 description 3
- 238000005507 spraying Methods 0.000 description 3
- 241000045403 Astragalus propinquus Species 0.000 description 2
- 108091028043 Nucleic acid sequence Proteins 0.000 description 2
- 238000012408 PCR amplification Methods 0.000 description 2
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 2
- 238000002123 RNA extraction Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000012153 distilled water Substances 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 108020004999 messenger RNA Proteins 0.000 description 2
- 229910052698 phosphorus Inorganic materials 0.000 description 2
- 239000011574 phosphorus Substances 0.000 description 2
- 108090000623 proteins and genes Proteins 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 229910000162 sodium phosphate Inorganic materials 0.000 description 2
- 239000002689 soil Substances 0.000 description 2
- 241000894007 species Species 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 229920002101 Chitin Polymers 0.000 description 1
- 108700031407 Chloroplast Genes Proteins 0.000 description 1
- 235000006200 Glycyrrhiza glabra Nutrition 0.000 description 1
- 244000303040 Glycyrrhiza glabra Species 0.000 description 1
- 235000000100 Hibiscus rosa sinensis Nutrition 0.000 description 1
- 108010070716 Intercellular Signaling Peptides and Proteins Proteins 0.000 description 1
- 102000005755 Intercellular Signaling Peptides and Proteins Human genes 0.000 description 1
- 108020005196 Mitochondrial DNA Proteins 0.000 description 1
- 108700001094 Plant Genes Proteins 0.000 description 1
- 238000003559 RNA-seq method Methods 0.000 description 1
- 241000220317 Rosa Species 0.000 description 1
- 240000008254 Rosa chinensis Species 0.000 description 1
- 235000000664 Rosa chinensis Nutrition 0.000 description 1
- 235000016785 Rosa della China Nutrition 0.000 description 1
- 244000154511 Rosa hybrid cultivar Species 0.000 description 1
- 235000010295 Rosa x kordesii Nutrition 0.000 description 1
- 238000010306 acid treatment Methods 0.000 description 1
- 208000005652 acute fatty liver of pregnancy Diseases 0.000 description 1
- 235000006533 astragalus Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000008014 freezing Effects 0.000 description 1
- 238000007710 freezing Methods 0.000 description 1
- LPLVUJXQOOQHMX-QWBHMCJMSA-N glycyrrhizinic acid Chemical compound O([C@@H]1[C@@H](O)[C@H](O)[C@H](O[C@@H]1O[C@@H]1C([C@H]2[C@]([C@@H]3[C@@]([C@@]4(CC[C@@]5(C)CC[C@@](C)(C[C@H]5C4=CC3=O)C(O)=O)C)(C)CC2)(C)CC1)(C)C)C(O)=O)[C@@H]1O[C@H](C(O)=O)[C@@H](O)[C@H](O)[C@H]1O LPLVUJXQOOQHMX-QWBHMCJMSA-N 0.000 description 1
- 239000001963 growth medium Substances 0.000 description 1
- 235000011477 liquorice Nutrition 0.000 description 1
- 230000002438 mitochondrial effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003147 molecular marker Substances 0.000 description 1
- 235000015097 nutrients Nutrition 0.000 description 1
- 102000054765 polymorphisms of proteins Human genes 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 238000007894 restriction fragment length polymorphism technique Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000013518 transcription Methods 0.000 description 1
- 230000035897 transcription Effects 0.000 description 1
- 230000017260 vegetative to reproductive phase transition of meristem Effects 0.000 description 1
- 229910052902 vermiculite Inorganic materials 0.000 description 1
- 235000019354 vermiculite Nutrition 0.000 description 1
- 239000010455 vermiculite Substances 0.000 description 1
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
Landscapes
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Engineering & Computer Science (AREA)
- Organic Chemistry (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Health & Medical Sciences (AREA)
- Biotechnology (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Immunology (AREA)
- Mycology (AREA)
- Microbiology (AREA)
- Molecular Biology (AREA)
- Botany (AREA)
- Biophysics (AREA)
- Physics & Mathematics (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
本发明公开了一种基于转录组中相邻微卫星序列的植物鉴别方法,包括以下步骤:提取待鉴定植物组织样本的RNA,使用高通量测序的方法对转录组进行测序,将测得的原始序列Raw data经过去掉接头序列,过滤掉不合格的序列后得到有效数据Clean data序列,由有效数据Clean data序列根据已知基因组进行拼接或根据重复序列拼接Unigene,使用MISA软件对拼接除重后得到的序列进行微卫星序列进行分析,得到微卫星序列统计数据,包括微卫星序列的类型和微卫星序列的具体序列;本方法可以灵敏地检测不同植物样本,该方法相较于基于形态的鉴别方法,具有准确可靠,分辨率高的特点,相较于传统基于DNA的分子鉴别方法,具有分辨率高的特点。
Description
技术领域
本发明涉及生物技术领域中区分植物样本,更具体的说是涉及一种基于转录组中相邻微卫星序列的植物鉴别方法。
背景技术
植物样本的精细鉴别在生产实践中有着重要的意义,尤其在经济作物和中药材领域,许多优良的品系遗传差异很小,需要高分辨力的鉴别方法。目前植物的鉴别一般采用形态鉴别法和基于DNA序列的分子鉴别方法。但是这两种方法在进行遗传差异非常小的植物品系,或同种植物由于不同化学处理、环境因素、生长周期等因素的不同,表达不同特定基因的样本时,经常会遇到分辨力不足的情况。
形态学鉴别方法应用范围广,但有些遗传差异并不会表现为形态学上的差异,这就无法运用形态学方式进行鉴别。而基于DNA序列的鉴别方法,如基于线粒体或叶绿体基因保守序列的分析方法,一般可以对属以上的分类单元进行分类。在对于种以下的植物分类单元进行鉴别时,处于成本考虑,一般采用分子标记的方法,如RAPD、SSR、ISSR、RFLP和AFLP等。这样的分子标记方法对于亲缘关系相对较近的植物具有较好的分辨率。但对于亲缘关系非常接近的品系,也会遇到分辨力不足的情况,因为许多存在于调控序列中的基因并不一定与分子标记相关联,从而无法用分子标记的方法鉴别。
然而,在实际的生产实践中,需要对植物组织样本进行更为细致的区分,如亲缘关系十分相近的植物,或无明显性状差异的优良品系;但上述方式就显得无能为力了。
植物中RNA的表达具有多态性,可以提供相较于DNA更多的遗传差异信息。由于目前对于RNA测序的转录组测序服务成本快速下降,该技术也将有更广泛的应用前景。微卫星序列是存在于基因组中的短串重复序列,这些序列可以转录,从而存在于植物的转录组数据中,由于RNA的转录受内部调控序列影响和外部环境影响,所以具有较高的多态性,微卫星序列按照重复单元中碱基的数目和类型一般可以分为单碱基重复、二碱基重复、三碱基重复、四碱基重复、五碱基重复、六碱基重复和相邻微卫星序列。其中相邻微卫星序列是指在两个微卫星序列相邻,之间一般有一段短的碱基序列,将两个微卫星序列隔开,因为中间短的微卫星序列有种属特异性,所以相较于其它微卫星序列类型具有较高的区分度。针对植物转录组中相邻微卫星序列相似度分析,可以揭示植物的亲缘关系,从而作为辅助植物鉴别的一种方法。
发明内容
本发明的目的是提供一种基于转录组中相邻微卫星序列的植物鉴别方法,该鉴别方法,相较于其它微卫星序列类型具有较高的区分度。针对植物转录组中相邻微卫星序列相似度分析,可以作为辅助植物鉴别的一种方法。
为解决上述技术问题,本发明通过一种基于转录组中相邻微卫星序列的植物鉴别方法进行实现,该方法包括以下步骤:
提取待鉴定植物组织样本的RNA,使用高通量测序的方法对转录组进行测序,将测得的原始序列Rawdata经过去掉接头序列,过滤掉不合格的序列后得到有效数据Cleandata序列,由有效数据Clean data序列根据已知基因组进行拼接或根据重复序列拼接Unigene,使用MISA软件对拼接除重后得到的序列进行微卫星序列进行分析,得到微卫星序列统计数据,包括微卫星序列的类型和微卫星序列的具体序列。在不同样本间进行相似度比对时,在微卫星序列的类型中选择有相邻的微卫星序列(c型),统计序列相同的微卫星序列,以相同微卫星序列个数除以除重后总的微卫星序列数得出相似度,利用相似度构建聚类图。从而实现不同植物样本的区分。
本发明设计了基于转录组中相邻微卫星序列的植物鉴别方法。本发明所使用的方法基于植物基因的转录本,具有多态性高、分辨力强和用量少的优点。对于亲缘关系十分接近的植物或不同外界环境因素处理下的同种植物样本的区分鉴别有着重要的意义。
本发明的其他特征和优点将在随后的说明书中阐述,或者,部分特征和优点可以从说明书推知或毫无疑义地确定,或者通过实施本发明的上述技术即可得知。
为使本发明的上述目的、特征和优点能更明显易懂,下文特举较佳实施方式,并配合所附附图,作详细说明如下。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明的最短距离法构建的RC、RJ和RS样本聚类图。
图2为本发明的最短距离法构建的RC、RJ、CKg、HH、CKh和LP样本聚类图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
下述实施例中所使用的实验方法如无特殊说明,均为常规方法。
下述实施例中所用的材料、试剂等,如无特殊说明,均可从商业途径得到。
目前对于植物组织样本进行更为细致的区分,如亲缘关系十分相近的植物,或无明显性状差异的优良品系进行区分鉴别时,是十分困难的,基于这一原因,本发明提供了一种基于转录组中相邻微卫星序列的植物鉴别方法;该方法主要是通过如下步骤加以实施:提取待鉴定植物组织样本的RNA,使用高通量测序的方法对转录组进行测序,将测得的原始序列Rawdata经过去掉接头序列,过滤掉不合格的序列后得到有效数据Clean data序列,由有效数据Clean data序列根据已知基因组进行拼接或根据重复序列拼接Unigene,使用MISA软件对拼接除重后得到的序列进行微卫星序列进行分析,得到微卫星序列统计数据,包括微卫星序列的类型和微卫星序列的具体序列。在不同样本间进行相似度比对时,在微卫星序列的类型中选择有相邻的微卫星序列(c型),统计序列相同的微卫星序列,以相同微卫星序列个数除以除重后总的微卫星序列数得出相似度,利用相似度构建聚类图。从而实现不同植物样本的区分。
为了更加方便本领域技术人员的理解,下面申请人以下面两个实施例加以说明。
实施例1:鉴别正常生长、水杨酸处理和茉莉酸甲酯处理的微型月季。
一、植物材料的准备
选取处于开花期的盆栽微型月季(Rosa hybrid minima),水杨酸处理组(RS)使用100 mL 1 mmol/L的水杨酸进行叶面喷施,茉莉酸甲酯处理组(RJ)使用100 mL 1 mmol/L的茉莉酸甲酯进行叶面喷施,对照组(RC)使用100 mL的蒸馏水,于24 h后取花瓣,取样后液氮速冻,用于总RNA的提取。
二、RNA提取
采用Invitrogen Trizol试剂提取微型月季花瓣的总RNA,提取方法按照试剂说明书操作。RNA的质量和完整性使用Bioanalyzer 2100和RNA 1000 Nano LabChip kit检测。总RNA使用连接有poly-T的磁珠进行两轮纯化。
三、RNA转录组文库的构建和测序
将RNA片段进一步分割后按Illumina的mRNA Seq sample preparation kit说明书反转录来创建cDNA文库。双链产物使用AMPure XP beads进行纯化,将双链产物的粘性末端修复为平末端并加接头,使用AMPureXP beads进行片段选择,将含有U的cDNA第二链酶解后进行PCR扩增以获得最终测序文库。测序文库在质检合格后进行测序,采用IlluminaHiseq平台。
四、序列的组装
测序产生的原始数据(raw data)使用cutadapt去除reads中的接头序列。使用fqtrim过滤不合格的序列,去除poly-A/T和截短后长度小于100 bp或N的含量在5%以上的序列,得到有效数据clean data。使用clean data根据重叠区域单独组装Unigene。最后进行归一化得到Unigene。
五、微卫星序列的统计
使用微卫星序列分析软件MISA分别对不同样本的Unigene,得到不同样本的微卫星序列统计数据。采用两两对比的方式,在微卫星序列的类型中选择有相邻微卫星序列类型(c型),统计其中序列相同的微卫星序列,除去重复序列后,各样本相邻微卫星序列数分别为RC510条,RJ394条,RS419条;RC与RJ相同序列为126条,RC与RS相同序列为127条,RJ与RS相同序列为129条。
以相同微卫星序列个数除以除去重复后总的微卫星序列数得出相似度,
公式为相似度=相同微卫星序列个数/(样本1微卫星序列个数+样本2微卫星序列个数-相同微卫星个数)
样本间相似度如表1所示:
表1 RC、RJ和RS样本间两两相似度
使用R语言的ggplot2包进行绘图,聚类计算函数采用hclust,计算方法采用最短距离法(single),得到聚类图如图1所示。
实施例2
一、植物材料的准备
盆栽蒙古黄芪(Astragalus membranaceus var. Mongolicus),培养基质采用蛭石和灭菌后的河沙,4 kg/盆。LP处理组通过添加NaH2PO4使磷含量达到15 mg/kg土。CKh处理组通过添加NaH2PO4使磷含量达到150 mg/kg土,各处理在培养过程中每周浇水1次,施加200mL无磷霍格兰营养液1次,生长6个月后采集幼苗,洗净后液氮速冻,用于总RNA的提取。
两年生甘草幼苗,CKg处理叶片喷施蒸馏水,HH处理叶面喷施10mg/L高分子量几丁质,24小时后采集叶片,洗净后液氮速冻,用于总RNA的提取。
二、RNA提取
采用Invitrogen Trizol试剂提取植物组织的总RNA,提取方法按照试剂说明书操作。RNA的质量和完整性使用Bioanalyzer 2100和RNA 1000 Nano LabChip kit检测。总RNA使用连接有poly-T的磁珠进行两轮纯化。
三、RNA转录组文库的构建和测序
将RNA片段进一步分割后按Illumina的mRNA Seq sample preparation kit说明书反转录来创建cDNA文库。双链产物使用AMPure XP beads进行纯化,将双链产物的粘性末端修复为平末端并加接头,使用AMPureXP beads进行片段选择,将含有U的cDNA第二链酶解后进行PCR扩增以获得最终测序文库。测序文库在质检合格后进行测序,采用IlluminaHiseq平台。
四、序列的组装
测序产生的原始数据(raw data)使用cutadapt去除reads中的接头序列。使用fqtrim过滤不合格的序列,去除poly-A/T和截短后长度小于100 bp或N的含量在5%以上的序列,得到有效数据clean data。使用clean data根据重叠区域单独组装,最后进行归一化得到Unigene。
五、微卫星序列的统计
使用微卫星序列分析软件MISA分别对不同样本的Unigene,得到不同样本的微卫星序列统计数据。在分析中同时也使用实施例1中的数据RC和RJ。
采用两两对比的方式,在微卫星序列的类型中选择有相邻微卫星序列类型(c型),统计其中序列相同的微卫星序列,除去重复序列后,各样本相邻微卫星序列数分别为RC510条,RJ394条,CKg,H1345条,HH1328条,CKh167条,LP152条;
CKg与CKh相同序列为1条,CKg与RC相同序列为3条,CKh与RC相同序列为3条,CKg与HH相同序列为181条,CKh与LP相同序列为17条,RC与RJ相同序列为126条,CKg与LP相同序列为1条,CKg与RJ相同序列为0条,CKh与HH相同序列为0条,CKh与RJ相同序列为1条,RC与HH相同序列为2条,RC与LP相同序列为0条,HH与LP相同序列为1条,HH与RJ相同序列为3条,LP与RJ相同序列为0条。
以相同微卫星序列个数除以除去重复后总的微卫星序列数得出相似度,
公式为相同微卫星序列个数/(样本1微卫星序列个数+样本2微卫星序列个数-相同微卫星个数)
样本间相似度如表2所示:
表2 CKg、CKh、RC、HH、LP和RJ样本间两两相似度
使用R语言的ggplot2包进行绘图,聚类计算函数采用hclust,计算方法采用最短距离法(single),得到聚类图如图2所示。
最后应说明的是:以上所述实施例,仅为本发明的具体实施方式,用以说明本发明的技术方案,而非对其限制,本发明的保护范围并不局限于此,尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化,或者对其中部分技术特征进行等同替换;而这些修改、变化或者替换,并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应所述以权利要求的保护范围为准。
Claims (1)
1.一种基于转录组中相邻微卫星序列的植物鉴别方法,其特征在于,包括以下步骤:
提取待鉴定植物组织样本的RNA,使用高通量测序的方法对转录组进行测序,将测得的原始序列Raw data经过去掉接头序列,过滤掉不合格的序列后得到有效数据Clean data序列,由有效数据Clean data序列根据重复序列进行拼接,由有效数据Clean data序列根据已知基因组进行拼接或根据重复序列拼接Unigene,使用MISA软件对拼接除重后得到的序列进行微卫星序列进行分析,得到微卫星序列统计数据,包括微卫星序列的类型和微卫星序列的具体序列;对不同样本中进行相似度比对,采用两两对比的方式,在微卫星序列的类型中选择相邻微卫星序列类型,统计其中序列相同的微卫星序列,以相同微卫星序列个数除以除重后总的微卫星序列数得出相似度,利用相似度构建聚类图;从而实现不同植物样本的区分。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110265657.0A CN112877464A (zh) | 2021-03-10 | 2021-03-10 | 一种基于转录组中相邻微卫星序列的植物鉴别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110265657.0A CN112877464A (zh) | 2021-03-10 | 2021-03-10 | 一种基于转录组中相邻微卫星序列的植物鉴别方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112877464A true CN112877464A (zh) | 2021-06-01 |
Family
ID=76041333
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110265657.0A Pending CN112877464A (zh) | 2021-03-10 | 2021-03-10 | 一种基于转录组中相邻微卫星序列的植物鉴别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112877464A (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2005200014A1 (en) * | 1999-12-24 | 2005-02-03 | Agriculture Victoria Services Pty Ltd | Molecular markers in clovers |
CN101684481A (zh) * | 2009-03-30 | 2010-03-31 | 电子科技大学 | 丹参est-ssr分子标记的制备方法、特异引物及其应用 |
US20150082486A1 (en) * | 2011-12-28 | 2015-03-19 | Evogene Ltd. | Isolated polynucleotides and polypeptides, and methods of using same for increasing yield of plants |
CN105349651A (zh) * | 2015-11-18 | 2016-02-24 | 广东省中药研究所 | Est-ssr标记鉴定中药溪黄草品种的方法和引物 |
CN107345256A (zh) * | 2017-08-22 | 2017-11-14 | 山西省农业科学院农作物品种资源研究所 | 一种基于转录组测序开发山黧豆est‑ssr引物组及方法和应用 |
CN109680091A (zh) * | 2019-01-10 | 2019-04-26 | 内蒙古科技大学包头医学院 | 一种基于高通量测序检测丛枝菌根真菌的引物及检测方法 |
CN112397149A (zh) * | 2020-11-11 | 2021-02-23 | 天津现代创新中药科技有限公司 | 无参考基因组序列的转录组分析方法及系统 |
-
2021
- 2021-03-10 CN CN202110265657.0A patent/CN112877464A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2005200014A1 (en) * | 1999-12-24 | 2005-02-03 | Agriculture Victoria Services Pty Ltd | Molecular markers in clovers |
CN101684481A (zh) * | 2009-03-30 | 2010-03-31 | 电子科技大学 | 丹参est-ssr分子标记的制备方法、特异引物及其应用 |
US20150082486A1 (en) * | 2011-12-28 | 2015-03-19 | Evogene Ltd. | Isolated polynucleotides and polypeptides, and methods of using same for increasing yield of plants |
CN105349651A (zh) * | 2015-11-18 | 2016-02-24 | 广东省中药研究所 | Est-ssr标记鉴定中药溪黄草品种的方法和引物 |
CN107345256A (zh) * | 2017-08-22 | 2017-11-14 | 山西省农业科学院农作物品种资源研究所 | 一种基于转录组测序开发山黧豆est‑ssr引物组及方法和应用 |
CN109680091A (zh) * | 2019-01-10 | 2019-04-26 | 内蒙古科技大学包头医学院 | 一种基于高通量测序检测丛枝菌根真菌的引物及检测方法 |
CN112397149A (zh) * | 2020-11-11 | 2021-02-23 | 天津现代创新中药科技有限公司 | 无参考基因组序列的转录组分析方法及系统 |
Non-Patent Citations (11)
Title |
---|
DAOJIE WANG等: ""Comparative transcriptome analyses of drought-resistant and - susceptible Brassica napus L. and development of EST-SSR markers by RNA-Seq"", 《J. PLANT BIOL.》 * |
DAOJIE WANG等: ""Comparative transcriptome analyses of drought-resistant and - susceptible Brassica napus L. and development of EST-SSR markers by RNA-Seq"", 《J. PLANT BIOL.》, vol. 58, 21 July 2015 (2015-07-21), pages 259 - 269 * |
MAISSA M. MORSI,等: "Genetic Variation Analysis for Breeding Purposes in Wild and Cultivated Sunflower Germplasm (Helianthus annuus Sp.)", PAKISTAN JOURNAL OF BIOLOGICAL SCIENCES, vol. 23, no. 12, pages 1539 - 1550 * |
刘畅等: ""苗药八爪金龙转录组测序与次生代谢产物合成相关基因的挖掘"", 《中草药》 * |
刘畅等: ""苗药八爪金龙转录组测序与次生代谢产物合成相关基因的挖掘"", 《中草药》, vol. 52, no. 5, 5 March 2021 (2021-03-05), pages 1434 - 1447 * |
杨进;: "SSR标记技术在植物遗传多样性研究上的应用", 中国农学通报, no. 07, pages 90 - 94 * |
王会敏,等: "夏枯草转录组SSR位点信息分析", 中国现代中药, vol. 20, no. 1, pages 30 * |
王建婷,等: "微型月季抗逆性基因表达与微卫星序列连锁分析", 安徽农业科学, vol. 49, no. 21, pages 124 - 127 * |
魏秀清,等: "莲雾转录组SSR信息分析及其分子标记开发", 园艺学报, vol. 45, no. 3 * |
魏秀清等: ""莲雾转录组 SSR 信息分析及其分子标记开发"", 《园艺学报》 * |
魏秀清等: ""莲雾转录组 SSR 信息分析及其分子标记开发"", 《园艺学报》, vol. 45, no. 3, 25 March 2018 (2018-03-25) * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wan et al. | The Welwitschia genome reveals a unique biology underpinning extreme longevity in deserts | |
Dash et al. | Transcriptome landscape at different developmental stages of a drought tolerant cultivar of flax (Linum usitatissimum) | |
Hiltunen et al. | Maintenance of high genome integrity over vegetative growth in the fairy-ring mushroom Marasmius oreades | |
CN105087574B (zh) | 藜麦est-ssr分子标记及其应用 | |
CN112342313B (zh) | 一种鉴别黄芩产地的snp分子标记物及其方法与应用 | |
CN108517368A (zh) | 利用上位性解析毛白杨LncRNA Pto-CRTG及其靶基因Pto-CAD5互作关系的方法及系统 | |
CN112877464A (zh) | 一种基于转录组中相邻微卫星序列的植物鉴别方法 | |
Ritter | Guava biotechnologies, genomic achievements and future needs | |
CN112322778A (zh) | 一种鉴别内蒙古产黄芩的snp分子标记及其方法与应用 | |
Winsett et al. | Global distribution and molecular diversity of Didymium difforme | |
糟谷大河 et al. | Reexamination of Crepidotus crocophyllus (Basidiomycota, Fungi) in Japan, with reference to its phylogenetic placement | |
CN109207617A (zh) | 一种植物植原体基因组的测序方法 | |
CN110511928B (zh) | 一种牡丹皮的转录组ssr分子标记体系及其应用 | |
CN103667452B (zh) | 一种利用ssr分子标记引物体系分析柑橘黄龙病菌遗传多样性的方法 | |
PRAKRONGRAK et al. | Genetic diversity and phylogenetic analysis of mayfly Caenis (Insecta: Ephemeroptera) using Cytochrome C Oxidase I (COI) and 12s rRNA genes from Thailand | |
CN113684291B (zh) | 一种鉴别东北地区黄芩的ssr分子标记引物及方法 | |
CN112342314B (zh) | 一种鉴别内蒙古杭锦旗甘草的snp分子标记物及其方法与应用 | |
CN111028889A (zh) | 一种获得活体营养型植物病原卵菌无污染基因组的方法 | |
CN113684293B (zh) | 一种鉴别山西和陕西地区黄芩的ssr分子标记引物及方法 | |
KR101634346B1 (ko) | 분자표지자를 이용한 리기다소나무와 테다소나무의 식별 방법 | |
KR102081586B1 (ko) | 양파의 항존 유전자 및 이를 분리하는 방법 | |
CN113073134B (zh) | 一种应用甘蔗着丝粒特异重复序列开展染色体准确计数的新方法 | |
CN113699262B (zh) | 一种鉴别内蒙古地区黄芩的ssr分子标记引物及方法 | |
CN113674003B (zh) | 一种利用宏基因组和机器学习进行蜂蜜蜜源地追踪的方法 | |
CN111549162B (zh) | 鉴别杜仲性别的引物、片段及方法 |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210601 |