CN111028889B - 一种获得活体营养型植物病原卵菌无污染基因组的方法 - Google Patents
一种获得活体营养型植物病原卵菌无污染基因组的方法 Download PDFInfo
- Publication number
- CN111028889B CN111028889B CN201911223140.4A CN201911223140A CN111028889B CN 111028889 B CN111028889 B CN 111028889B CN 201911223140 A CN201911223140 A CN 201911223140A CN 111028889 B CN111028889 B CN 111028889B
- Authority
- CN
- China
- Prior art keywords
- species
- genome
- pollution
- result
- reads
- 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
- 238000000034 method Methods 0.000 title claims abstract description 50
- 241000233654 Oomycetes Species 0.000 title claims abstract description 45
- 230000001717 pathogenic effect Effects 0.000 title claims abstract description 44
- 238000001727 in vivo Methods 0.000 title claims abstract description 8
- 235000016709 nutrition Nutrition 0.000 title claims abstract description 8
- 241000894007 species Species 0.000 claims abstract description 108
- 108090000623 proteins and genes Proteins 0.000 claims abstract description 31
- 238000009826 distribution Methods 0.000 claims abstract description 9
- 238000012165 high-throughput sequencing Methods 0.000 claims abstract description 5
- 238000012163 sequencing technique Methods 0.000 claims description 38
- 238000011109 contamination Methods 0.000 claims description 16
- 125000003275 alpha amino acid group Chemical group 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 claims description 6
- 102000004169 proteins and genes Human genes 0.000 claims description 5
- 241000894006 Bacteria Species 0.000 claims description 3
- 238000010586 diagram Methods 0.000 claims description 3
- 238000007619 statistical method Methods 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 claims description 2
- 238000012268 genome sequencing Methods 0.000 claims description 2
- 230000003950 pathogenic mechanism Effects 0.000 abstract description 4
- 238000011160 research Methods 0.000 abstract description 4
- 241000196324 Embryophyta Species 0.000 description 41
- 241001466451 Stramenopiles Species 0.000 description 7
- 241001281803 Plasmopara viticola Species 0.000 description 6
- 230000000717 retained effect Effects 0.000 description 4
- 239000003153 chemical reaction reagent Substances 0.000 description 3
- 230000014759 maintenance of location Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 241000233866 Fungi Species 0.000 description 2
- 241000233626 Plasmopara Species 0.000 description 2
- 244000052616 bacterial pathogen Species 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 239000000356 contaminant Substances 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 244000005700 microbiome Species 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 235000015097 nutrients Nutrition 0.000 description 2
- 241000219194 Arabidopsis Species 0.000 description 1
- 241000123650 Botrytis cinerea Species 0.000 description 1
- 241000252254 Catostomidae Species 0.000 description 1
- 241000195493 Cryptophyta Species 0.000 description 1
- 241000206602 Eukaryota Species 0.000 description 1
- 241001465754 Metazoa Species 0.000 description 1
- 241001223281 Peronospora Species 0.000 description 1
- 241001281805 Pseudoperonospora cubensis Species 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000009412 basement excavation Methods 0.000 description 1
- 238000005202 decontamination Methods 0.000 description 1
- 230000003588 decontaminative effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000001963 growth medium Substances 0.000 description 1
- 208000015181 infectious disease Diseases 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000003032 phytopathogenic effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 241000215338 unidentified plant Species 0.000 description 1
- 238000012070 whole genome sequencing analysis Methods 0.000 description 1
Classifications
-
- 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
- G16B30/00—ICT specially adapted for sequence analysis involving nucleotides or amino acids
- G16B30/10—Sequence alignment; Homology search
-
- 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
- G16B30/00—ICT specially adapted for sequence analysis involving nucleotides or amino acids
- G16B30/20—Sequence assembly
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biophysics (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Biotechnology (AREA)
- Evolutionary Biology (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Theoretical Computer Science (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明公开了一种获得活体营养型植物病原卵菌无污染基因组的方法:先鉴定污染源,从NCBI下载候选污染物种对应的全基因组序列,与高通量测序得到的待组装的原始reads比对,去除对比到污染源的reads,得到初步组装结果,再将其分为多个窗口,统计每个窗口的基因组区域的深度平均值和GC含量,通过GC分布将植物病原卵菌和污染源物种基因组的窗口区域区分开,进而获得污染源物种基因组对应的窗口区域,去除map到污染源物种基因组窗口区域的reads序列,保留剩余的reads,组装后得到最终的组装结果。从根本上保证获得的活体营养型植物病原卵菌基因组序列信息的完整性和准确性,为后续致病相关基因、病原卵菌致病机制等的进一步研究,提供了有力保证。
Description
技术领域
本发明涉及生物信息学技术领域,具体涉及一种获得活体营养型植物病原卵菌无污染基因组的方法。
背景技术
卵菌是一类真核生物,包括许多植物病原菌,所致病害给许多农作物和花卉植物造成毁灭性危害。卵菌具有独特分类地位的群体,由于表现出丝状等特性传统上被划分到真菌界中,随着学科的发展,卵菌纲早已从真菌中划分到藻界或茸鞭生物界。
活体营养型植物病原卵菌的基因组中缺失了合成某些生命活动所必须的营养物质的相关基因,因此其相关的营养物质必须通过特殊的吸器从植物细胞获得,所以只能与植物寄生在一起,导致该类型的致病菌无法离开寄主植物存活繁殖,无法对活体营养型植物病原卵菌在离体的无菌培养基上进行培养,这类病原菌有葡萄霜霉菌、拟南芥霜霉菌、黄瓜霜霉菌等。因此在活体营养型植物病原卵菌的菌株培养和收集过程中,必须通过采集寄主植物的组织,如植物叶片对其进行侵染培养,在整个操作过程中,不可避免的会接触外界空气、环境,植物的组织也不可能生长在无菌环境。以上特点,决定了对活体营养型植物病原卵菌的收集和DNA提取过程中,无法避免寄主植物和环境微生物的DNA污染掺入,特别是环境微生物。通常情况下,在最终的高通量测序数据中,有相当一部分测序数据(最多可达到40%~50%左右)为外源污染数据,因此需要有效分析鉴定出外源污染数据,组装活体营养型植物病原卵菌无污染基因组。
目前公开发表的活体营养型植物病原卵菌无污染基因组组装方法有两种:
一、在2018年发表在Genome Biol.Evol.的“A High-QualityGrapevineDownyMildewGenome Assembly Reveals Rapidly Evolving and Lineage-SpecificPutative Host Adaptation Genes”文章,提及了一种相对完整的去污染和组装方法:
1.用所有的原始测序数据进行组装,将所有的contig组装结果与NCBI的NT库比对,并同时将所有的原始reads map到组装结果;
2.根据基因组与NT库的比对结果,鉴定到属于污染物种的contig,并将map到污染物种contig的reads扔掉,保留剩余的reads;
3.再组装,并将1、2步再重复进行,直至没有更多的污染物种contig和reads被发现;
4.将经过上述操作得到的组装结果,结合reads的map深度信息相结合,丢弃掉覆盖深度低于20X,且没有blast比对到植物病原卵菌或同源物种上的contig。
该文章确实提出了一种相对完善的污染鉴定、去除及组装的方法,但作者在文章中同时也指出了存在的问题:该组装方法一定程度上将污染情况降低到一个相对的低水平,但可能同时将一些本来想留下的活体营养型植物病原卵菌的测序数据也扔掉了。这可能是在进行第4步时导致的。由于测序的偏好性,一些活体营养型植物病原卵菌的基因组区域测序深度较低;另外,一些没有比对到植物病原卵菌或同源物种上的contig及其原始reads,可能是一些潜在的尚未被发现的未知植物病原卵菌基因组区域,因此在已知的公共数据库中没有记录;这些基因组区域及其对应的测序数据,也被该方法错误的丢弃了。
另外,该方法除了作者提到的问题之外,仍存在另外两个问题:1)仍无法真正的将所有潜在的污染源序列去除,有部分污染源序列被遗漏;2)存在本该留下的植物病原卵菌序列被错误的去除的可能性。出现以上问题的原因是由于该方法在鉴定污染时,仅仅是以组装结果中,鉴定到的属于污染物种的contig作为参考去除map到这些contig上的reads,而不是以该污染物种对应的完整的基因组序列作为参考。而实际上,由于用于组装的原始数据中,只含有部分的污染物种序列,大部分的序列仍为测序的植物病原卵菌序列,且同时含有其他污染物种的序列;由于污染物种原始数据的测序深度不足导致最终组装结果中污染物种的contig很有可能仅仅是污染物种基因组的一部分,很多污染物种的基因组区域,虽然有测序数据覆盖,但由于测序深度低而没有组装出来;另外,因为原始数据的物种来源混杂,最终的contig组装结果往往不是来自某个单一物种,而是中间还掺杂着其他物种的信息,比如待测序的植物病原卵菌。所以,以这些contig为参考,最终去除及保留的reads,存在两个问题:一是可能去除不够彻底,有些reads属于污染物种,但由于没有被组装到contig中,导致无法通过reads map剔除出来;二是由于contig中可能掺杂其他物种,特别是待测序的植物病原卵菌,导致本该留下的植物病原卵菌reads序列也被错误的剔除出去了。
二、2017年发表在《Scientific Reports》的“A multi-omics study ofthegrapevine-downy mildew(Plasmopara viticola)pathosystem unveils acomplexprotein coding and noncoding-based arms race during infection”文章,也提及了一种污染处理及组装的方法:
1.将所有的数据进行组装,将组装完成后的基因组用AUGUSTUS软件进行基因预测;
2.将预测出的蛋白序列与NCBI的NR库进行blastp比对,获得物种信息;
3.如果整条scaffold中预测出的基因的所属物种均为Stramenopiles类的,被认为是Plasmopara类,即我们需要的组装结果;
4、如果整条scaffold中预测出的基因的所属物种均不是Stramenopiles类的,即不是我们需要的组装结果,这些scaffold将被扔掉;
5.如果整条scaffold中预测出的基因的所属物种既有Stramenopiles类,也有非Stramenopiles类的,对这些scaffold进行截断,根据基因预测的位置信息和测序深度的覆盖信息,将预测为Stramenopiles类或者没有物种信息的,且测序深度覆盖较高的区域,截断并作为候选scaffold保留;
6.对于所有预测出的基因与NCBI的NR库均没有比对结果,或者没有任何基因预测结果的scaffold,将这些scaffold完整的与NCBI的NT库比对,获取其物种信息;将所属物种信息为Stramenopiles类或者没有任何比对结果的scaffold认为是Plasmopara类,将其留下;将剩余的比对结果为非Stramenopiles类的scaffold扔掉;
7.将以上步骤保留下来的scaffold汇总并作为reference,将原始的测序readsmap到这个reference,将map上的序列保留,并用这些reads重新进行组装,得到最终的组装结果。
虽然该文章提出的方法也能在一定程度上解决污染源数据干扰组装的问题,从而获得一个更纯净的基因组组装结果;但从以上方法的步骤中不难看出该方法存在的问题:
由于第一个版本的组装结果,是将包含污染源物种的测序数据一起混合组装的,所以在最终得到的组装结果中,有相当一部分scaffold组装结果是混合了两个或多个物种信息的,而有些区域,是没有基因注释的。但是该方法基于先是基于基因注释信息进行污染鉴定,会直接导致步骤3、4、5存在鉴定错误的情况。
比如步骤3中判定的scaffold中,可能存在非Stramenopiles类且没有基因注释的基因组区域被错误的认定为Stramenopiles类,这些区域被错误的保留了。步骤4中判定的scaffold中,可能存在Stramenopiles类且没有基因注释的基因组区域被错误的认定为非Stramenopiles类,这些区域被错误的剔除了。步骤5中判定的scaffold中最终被保留下来的区域以及被剔除的区域中,同时存在以上步骤3和4中被错误保留或错误剔除的可能。步骤6中根据与NCBI的NT库比对的结果决定scaffold的去留,本质上与基于预测的基因集与NCBI的NR库比对的结果决定scaffold的去留一致,同样存在被错误保留或错误剔除的可能。
同时,很多情况下,在组装结果中,往往是活体营养型植物病原卵菌的基因组序列与污染物种的基因组序列穿插糅合在一起,这导致步骤5中界定两者的边界无法精确;很难找到一个十分准确的碱基位置,在该位置的一边是测序物种序列,另一边紧接着就是污染物种的序列。另外一个不容忽视的问题是有些活体营养型植物病原卵菌的基因组区域,由于测序深度低等原因,未能组装出来,这些区域对应的原始reads信息,无法通过该方法被识别并保留。
总之,通过以上两种方法最终被保留的reads中,都可能仍然存在污染物种的信息,且都存在有些本该属于活体营养型植物病原卵菌的reads数据被错误的丢弃或遗漏的现象。利用这些最终的reads进行二次组装,自然无法最大化的保证组装结果的纯净度,也无法确保基因组的完整度。
本申请的发明人在发明专利《一种快速准确鉴定高通量基因组数据污染源的方法》(CN201610117589)中所提出的方法能够快速准确的鉴定分析出测序数据中的污染源类型及污染程度的情况,但该方法仅仅局限在污染源的鉴定层面,但如何对污染数据进行有效的精准快速去除,并进一步获得准确的基因组组装结果,该方法并没有涉及。
综上所述,提出一套可靠的分析方法,对污染数据进行高效鉴定,并对污染数据进行精准的删除过滤,并最终得到纯净的活体营养型植物病原卵菌基因组组装结果,直接关系到后续基因组分析和致病基因预测的准确性,对后续科学研究具有重要意义。
发明内容
本发明的目的是为解决以上提出的问题,在发明专利《一种快速准确鉴定高通量基因组数据污染源的方法》(CN201610117589)研究的基础上,提供一种获得活体营养型植物病原卵菌无污染基因组的方法。
本发明的方法是通过以下步骤实现的:
(1)使用发明专利《一种快速准确鉴定高通量基因组数据污染源的方法》(CN201610117589)公开的方法进行污染源的鉴定工作,基于获得的污染源对应的物种信息,直接从NCBI下载候选污染物种对应的全基因组序列,并将这些基因组序列合并到一个文件中;
其中,所述污染源的鉴定工作步骤为:
(a)组装denovo测序的原始基因组测序数据,得到组装结果;
(b)将组装结果与NCBI的NT数据库进行blast比对,得到有同源性的序列,作为原始比对数据库;
(c)从原始比对数据库中,提取序列对应的物种信息并排序,将序列对应的物种从多到少进行排序,判断是否存在外源污染;
(d)将组装结果进行基因结构的注释并获得对应的基因所翻译蛋白的氨基酸序列,将氨基酸序列与NCBI的NR数据库进行blast比对;
(e)得到步骤(d)的比对结果后,再依照步骤(c)的方法,提取物种信息和排序,将氨基酸序列对应的物种从多到少进行排序,判断是否存在外源污染;
(f)结合步骤(c)和步骤(e)的结果,根据两步分析统计中均存在的结果,最终确定污染情况及具体污染源;
其中,步骤(a)中采用基于k-mer的开源基因组组装软件进行组装;
其中,步骤(b)中进行blast比对时使用的具体命令为:blastall-p blastn-dBlastDB-i in_file.fasta-m8–b 20-e 1e-5>blast_output;
其中,步骤(c)中,提取物种信息和排序的功能是通过基于Java程序的算法来实现的,具体包含以下两个步骤:
①java PollutedIdentificationgetAnnotationFromNCBIDBByBlastOutblastOutResultncbiDB>result1,本步骤从NCBI的数据库中获取原始比对数据库中每条序列对应的注释信息,注释信息中就包含对应的物种信息;②java PollutedIdentificationgetPollutedSpeciesInfoByNCBIDBBlastOut result1ThreadHold>result2,本步骤设定e value≤1e-5、identity≥80%,将第①步的结果中的物种信息进行统计,凡是大于设定阈值的物种都会被统计,从而得到最终的物种统计结果,最终统计结果的格式为:每一个物种信息占一行,每行两列,第一列为物种信息,第二列为对应的比对到此物种的query序列总数,统计结果以文本文档的格式用记事本或Notepad++软件打开,然后将内容拷贝到excel表格中,按第二列统计数字按从大到小排序,排序后,就可以清晰直观的观察到最终的组装序列对应的物种从多到少的排序,进而根据物种的信息,判断是否有污染,以及主要污染源及其污染情况;
其中,在完成步骤(c)之后,根据排序结果判断,如排在最前面的不是测序物种及近源物种,而是细菌或其他明显为污染的物种(比如测的是植物,而排在前面的是细菌、动物等),且e value≤1e-5、identity≥80%,就说明污染较严重;或排在前面的是近源物种,但其他非近源物种的e value≤1e-5、identity≥80%,且每个占比或总量≥1%,也说明有污染;
其中,在步骤(d)中,使用AUGUSTUS软件将组装结果进行基因结构的注释并获得对应的基因所翻译蛋白的氨基酸序列;进行blast比对时使用的具体命令为:blastall-pblastp-dBlastDB-i in_file.fasta-m8–b 20-e 1e-5>blast_output2;
(2)以步骤(1)合并的污染源全基因组序列文件作为reference,将高通量测序得到的待组装的原始reads map到reference,根据比对结果,将map到污染源reference的reads去除,没有map上的reads保留;
(3)将保留的reads再组装,基于组装结果,重复(1)、(2)步骤,直至没有更多的污染数据被去除,并利用这些最终保留的reads进行组装,获得组装结果;
(4)将获得的组装结果作为reference,将用于组装的reads map比对到该reference,保留原始的比对结果;
(5)将步骤(4)使用的reference分成多个个窗口,根据步骤(4)中的比对结果统计每个窗口的基因组区域的深度平均值和GC含量;
(6)以步骤(5)统计的每个窗口的GC含量为横坐标,以对应的深度平均值为纵坐标,在二维坐标系中标注出所有窗口GC含量和深度平均值的统计结果,最终绘制成一个GC分布图,由于活体营养型植物病原卵菌的GC含量和测序深度与可能的污染源物种存在差异,因此属于植物病原卵菌的统计点会集中分布在一起,与可能的污染源物种的统计点分布在坐标系中会产生明显的分离,因此可通过GC分布将植物病原卵菌和污染源物种基因组的窗口区域区分开,进而获得污染源物种基因组对应的窗口区域;
(7)获得污染源物种基因组的窗口区域,并根据第4步的map结果,去除map到污染源物种基因组窗口区域的reads序列,保留剩余的reads;
(8)将剩余的reads再组装,得到最终的组装结果。
其中,在步骤(5)中,将reference以每500bp作为一个窗口进行统计,依次统计窗口内每一个位点的read map深度,计算深度平均值;同时统计窗口内500bp的基因组区域的GC含量;统计结束后,将统计窗口后移,统计下一个500bp的窗口的GC含量和深度;直至将全基因组的GC含量和深度统计完毕。
其中,深度平均值的计算方法为:500个位点的深度值累加后,再除以500;GC含量的统计方法为:窗口内的G或C的碱基总数,除以500,再乘以100%。
有益效果:
本发明的方法步骤(1)~(3)通过将组装结果与NT和NR库同时比对进行污染源鉴定,下载已鉴定的污染源物种的完整基因组序列作为污染源reference进行污染源reads的去除,取代前人方法中以组装出的代表污染源部分基因组的contig组装结果作为reference作为污染源reads鉴定的方法,从根本上保证最全面的去除属于污染源物种的reads,且彻底杜绝出现将活体营养型植物病原卵菌reads被错误去除的可能性;
本发明的方法步骤(4)~(8)通过结合GC depth分布分析的方法,将一些尚未被测序公布,在公共数据库中没有任何参考信息的污染源物种基因组区域,进一步进行鉴定隔离,并将对应的原始reads分离去除,将最终去除后的reads进行组装,获得最终的组装结果。
本发明的方法彻底解决了活体营养型病原卵菌基因组组装过程中,污染源鉴定效率低、速度慢、鉴定不准确、污染源数据去除不完全不彻底、将植物病原卵菌reads数据错误去除的现象,并且彻底解决了仅仅依赖已知数据库,无法鉴定和去除未知污染源数据的短板。依赖更完整更纯净的活体营养型植物病原卵菌原始reads信息,从根本上保证获得的活体营养型植物病原卵菌基因组序列信息的完整性和准确性,为后续挖掘致病相关基因并对其进行致病机理研究,揭示活体营养型植物病原卵菌的致病机制及其与寄主植物的互作机制,提供有力保证。
本发明的方法鉴定和去除效率高,时间周期短,操作简便,为活体营养型病原卵菌基因组污染鉴定和基因组组装提供了一套切实可操作的方法流程。
具体实施方式
以下的实施例便于更好地理解本发明,但并不限定本发明。下述实施例中的实验方法,如无特殊说明,均为常规方法。下述实施例中所用的试验材料和试剂,如无特殊说明,均为常规生化试剂,购自试剂公司。
实施例:葡萄霜霉菌无污染基因组的组装
(1)对一例葡萄霜霉菌进行PacBio三代全基因组测序,共得到1,836,120条原始reads,共计19.2G的数据量;使用发明专利《一种快速准确鉴定高通量基因组数据污染源的方法》(CN201610117589)公开的方法进行污染源的鉴定工作,鉴定到包括寄主植物在内的246种污染源物种,直接从NCBI下载候选污染物种对应的全基因组序列,并将这些基因组序列合并到一个文件中;
(2)以步骤(1)合并的污染源全基因组序列文件作为reference,将高通量测序得到的待组装的原始reads map到reference,根据比对结果,将map到污染源reference的reads去除,没有map上的reads保留;
(3)将保留的reads再组装,基于组装结果,重复(1)、(2)步骤,直至没有更多的污染数据被去除,共去除污染reads有373,076条,共计达4.6G的数据量,剩余reads有1,463,044条,共计14.7G,去除的污染数据量占原始数据量的24%,利用最终保留的reads进行组装,得到1598条contig,共计98.2M的基因组,contig N50达到404.7K,GC含量为45.21%;
(4)将获得的组装结果作为reference,将用于组装的reads map比对到该reference,保留原始的比对结果;
(5)将步骤(4)使用的reference以每500bp作为一个窗口进行统计,依次统计窗口内每一个位点的read map深度,计算深度平均值;同时统计窗口内500bp的基因组区域的GC含量;统计结束后,将统计窗口后移,统计下一个500bp的窗口的GC含量和深度,直至将全基因组的GC含量和深度统计完毕;
(6)以步骤(5)统计的每个窗口的GC含量为横坐标,以对应的深度平均值为纵坐标,在二维坐标系中标注出所有窗口GC含量和深度平均值的统计结果,最终绘制成一个GC分布图,由于葡萄霜霉菌的GC含量和测序深度与可能的污染源物种存在差异,因此属于葡萄霜霉菌的统计点会集中分布在一起,与可能的污染源物种的统计点分布在坐标系中会产生明显的分离,因此可通过GC分布将葡萄霜霉菌和污染源物种基因组的窗口区域区分开,进而获得污染源物种基因组对应的窗口区域;
其中,深度平均值的计算方法为:500个位点的深度值累加后,再除以500;GC含量的统计方法为:窗口内的G或C的碱基总数,除以500,再乘以100%;
通过该方法分析发现所有的窗口GC分布均集中在一个范围内,并未发生分离,因此未发现新的可能污染物种;说明通过前期的去污染方法,已较好的去除了污染源,最终获得了98.2M的葡萄霜霉菌基因组。
虽然,上文中已经用一般性说明及具体实施方案对本发明作了详尽的描述,但在本发明基础上,可以对之作一些修改或改进,这对本领域技术人员而言是显而易见的。因此,在不偏离本发明精神的基础上所做的这些修改或改进,均属于本发明要求保护的范围。
Claims (4)
1.一种获得活体营养型植物病原卵菌无污染基因组的方法,其特征在于,包括以下步骤:
(1)进行污染源的鉴定工作,基于获得的污染源对应的物种信息,直接从NCBI下载候选污染物种对应的全基因组序列,并将这些基因组序列合并到一个文件中;
其中,所述污染源的鉴定工作步骤为:
(a)组装denovo测序的原始基因组测序数据,得到组装结果;
(b)将组装结果与NCBI的NT数据库进行blast比对,得到有同源性的序列,作为原始比对数据库;
(c)从原始比对数据库中,提取序列对应的物种信息并排序,将序列对应的物种从多到少进行排序,判断是否存在外源污染;
(d)将组装结果进行基因结构的注释并获得对应的基因所翻译蛋白的氨基酸序列,将氨基酸序列与NCBI的NR数据库进行blast比对;
(e)得到步骤(d)的比对结果后,再依照步骤(c)的方法,提取物种信息和排序,将氨基酸序列对应的物种从多到少进行排序,判断是否存在外源污染;
(f)结合步骤(c)和步骤(e)的结果,根据两步分析统计中均存在的结果,最终确定污染情况及具体污染源;
其中,步骤(a)中采用基于k-mer的开源基因组组装软件进行组装;
其中,步骤(b)中进行blast比对时使用的具体命令为:blastall-p blastn-dBlastDB-i in_file.fasta-m8–b 20-e 1e-5>blast_output;
其中,步骤(c)中,提取物种信息和排序的功能是通过基于Java程序的算法来实现的,具体包含以下两个步骤:
①java PollutedIdentificationgetAnnotationFromNCBIDBByBlastOutblastOutResultncbiDB>result1,本步骤从NCBI的数据库中获取原始比对数据库中每条序列对应的注释信息,注释信息中就包含对应的物种信息;②java PollutedIdentificationgetPollutedSpeciesInfoByNCBIDBBlastOut result1ThreadHold>result2,本步骤设定e value≤1e-5、identity≥80%,将第①步的结果中的物种信息进行统计,凡是大于设定阈值的物种都会被统计,从而得到最终的物种统计结果,最终统计结果的格式为:每一个物种信息占一行,每行两列,第一列为物种信息,第二列为对应的比对到此物种的query序列总数,统计结果以文本文档的格式用记事本或Notepad++软件打开,然后将内容拷贝到excel表格中,按第二列统计数字按从大到小排序,排序后,就清晰直观的观察到最终的组装序列对应的物种从多到少的排序,进而根据物种的信息,判断是否有污染,以及主要污染源及其污染情况;
其中,在完成步骤(c)之后,根据排序结果判断,如排在最前面的不是测序物种及近源物种,而是细菌或其他明显为污染的物种,且e value≤1e-5、identity≥80%,就说明污染较严重;或排在前面的是近源物种,但其他非近源物种的e value≤1e-5、identity≥80%,且每个占比或总量≥1%,也说明有污染;
其中,在步骤(d)中,使用AUGUSTUS软件将组装结果进行基因结构的注释并获得对应的基因所翻译蛋白的氨基酸序列;进行blast比对时使用的具体命令为:blastall-p blastp-dBlastDB-i in_file.fasta-m8–b 20-e 1e-5>blast_output2;
(2)以步骤(1)合并的污染源全基因组序列文件作为reference,将高通量测序得到的待组装的原始reads map到reference,根据比对结果,将map到污染源reference的reads去除,没有map上的reads保留;
(3)将保留的reads再组装,基于组装结果,重复(1)、(2)步骤,直至没有更多的污染数据被去除,并利用这些最终保留的reads进行组装,获得组装结果;
(4)将获得的组装结果作为reference,将用于组装的reads map比对到该reference,保留原始的比对结果;
(5)将步骤(4)使用的reference分成多个窗口,根据步骤(4)中的比对结果统计每个窗口的基因组区域的深度平均值和GC含量;
(6)以步骤(5)统计的每个窗口的GC含量为横坐标,以对应的深度平均值为纵坐标,在二维坐标系中标注出所有窗口GC含量和深度平均值的统计结果,最终绘制成一个GC分布图,由于活体营养型植物病原卵菌基因组的GC含量和测序深度与可能的污染源物种存在差异,因此属于植物病原卵菌的统计点会集中分布在一起,与可能的污染源物种的统计点分布在坐标系中会产生明显的分离,因此通过GC分布将植物病原卵菌和污染源物种基因组的窗口区域区分开,进而获得污染源物种基因组对应的窗口区域;
(7)获得污染源物种基因组的窗口区域,并根据第4步的map结果,去除map到污染源物种基因组窗口区域的reads序列,保留剩余的reads;
(8)将剩余的reads再组装,得到最终的组装结果。
2.按照权利要求1所述的方法,其特征在于:在步骤(5)中,将reference以每500bp作为一个窗口进行统计,依次统计窗口内每一个位点的read map深度,计算深度平均值;同时统计窗口内500bp的基因组区域的GC含量;统计结束后,将统计窗口后移,统计下一个500bp的窗口的GC含量和深度;直至将全基因组的GC含量和深度统计完毕。
3.按照权利要求1所述的方法,其特征在于,在步骤(5)中深度平均值的计算方法为:500个位点的深度值累加后,再除以500。
4.按照权利要求1所述的方法,其特征在于,在步骤(5)中GC含量的统计方法为:窗口内的G或C的碱基总数,除以500,再乘以100%。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911223140.4A CN111028889B (zh) | 2019-12-03 | 2019-12-03 | 一种获得活体营养型植物病原卵菌无污染基因组的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911223140.4A CN111028889B (zh) | 2019-12-03 | 2019-12-03 | 一种获得活体营养型植物病原卵菌无污染基因组的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111028889A CN111028889A (zh) | 2020-04-17 |
CN111028889B true CN111028889B (zh) | 2021-04-20 |
Family
ID=70207800
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911223140.4A Active CN111028889B (zh) | 2019-12-03 | 2019-12-03 | 一种获得活体营养型植物病原卵菌无污染基因组的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111028889B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114627966A (zh) * | 2020-12-11 | 2022-06-14 | 厦门大学 | 一种评估用于检测病原体的试剂或试剂组合的有效性的方法及装置 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104164479B (zh) * | 2014-04-04 | 2017-09-19 | 深圳华大基因科技服务有限公司 | 杂合基因组处理方法 |
CN105653899B (zh) * | 2014-09-30 | 2018-02-09 | 深圳华大基因研究院 | 同时确定多种样本的线粒体基因组序列信息的方法和系统 |
CN104630206A (zh) * | 2015-02-05 | 2015-05-20 | 北京诺禾致源生物信息科技有限公司 | 转录组文库的构建方法 |
CN105740650B (zh) * | 2016-03-02 | 2019-04-05 | 广西作物遗传改良生物技术重点开放实验室 | 一种快速准确鉴定高通量基因组数据污染源的方法 |
CN108504750B (zh) * | 2018-04-23 | 2021-10-15 | 深圳华大法医科技有限公司 | 确定菌群snp位点集合的方法、系统及其应用 |
CN109023537B (zh) * | 2018-09-04 | 2021-10-08 | 上海交通大学 | 一种微量dna样品高通量测序文库的构建技术 |
CN109559780A (zh) * | 2018-09-27 | 2019-04-02 | 华中科技大学鄂州工业技术研究院 | 一种高通量测序的rna数据处理方法 |
CN109912699B (zh) * | 2019-05-05 | 2020-09-29 | 南京林业大学 | 樟疫霉效应子蛋白Avh87及其编码基因与应用 |
-
2019
- 2019-12-03 CN CN201911223140.4A patent/CN111028889B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN111028889A (zh) | 2020-04-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111951895B (zh) | 基于宏基因组学的病原分析方法、分析装置、设备及存储介质 | |
CN109273053B (zh) | 一种高通量测序的微生物数据处理方法 | |
Schäferhoff et al. | Towards resolving Lamiales relationships: insights from rapidly evolving chloroplast sequences | |
CN105740650B (zh) | 一种快速准确鉴定高通量基因组数据污染源的方法 | |
CN109346130B (zh) | 一种直接从全基因组重测序数据中得到微单体型及其分型的方法 | |
US20200294628A1 (en) | Creation or use of anchor-based data structures for sample-derived characteristic determination | |
CN115719616B (zh) | 一种病原物种特异性序列的筛选方法及系统 | |
CN113470743A (zh) | 一种基于bd单细胞转录组和蛋白组测序数据的差异基因分析方法 | |
CN111028889B (zh) | 一种获得活体营养型植物病原卵菌无污染基因组的方法 | |
CN103911442A (zh) | 一种基于ssr标记和毛细管电泳技术的水稻指纹图谱构建方法 | |
CN115662516A (zh) | 一种基于二代测序技术的高通量预测噬菌体宿主的分析方法 | |
CN110444253B (zh) | 一种适用于混池基因定位的方法及系统 | |
CN110970091A (zh) | 标签质控的方法及装置 | |
CN111916151A (zh) | 一种苜蓿黄萎病菌的溯源检测方法及应用 | |
CN117089644A (zh) | 用于芦竹品种鉴定的mnp标记位点、引物组合物和试剂盒及其应用 | |
CN110592253A (zh) | 鉴定云南茶树品种的dna组合条形码及其鉴定方法 | |
CN115612751A (zh) | 西瓜全基因组液相芯片及应用 | |
CN114822697A (zh) | 一种利用宏基因组分析溯源土壤耐药基因污染的方法 | |
CN108009401B (zh) | 一种筛选指纹图谱遗传标记的方法 | |
CN101565744A (zh) | 一种三疣梭子蟹多元高通量遗传标记系统及遗传分析方法 | |
Ali et al. | Plant biotechnology research with single-cell transcriptome: recent advancements and prospects | |
CN211578386U (zh) | 一种宏基因组分析装置 | |
CN114277164B (zh) | 一种肺炎链球菌的mnp标记组合、引物对组合、试剂盒及其应用 | |
CN116646010B (zh) | 人源性病毒检测方法及装置、设备、存储介质 | |
Veres et al. | No evidence for phylogenetic structure or environmental filtering of springtail microbiomes |
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 |