Name: 以乌苏图病毒为例，验证全基因组纳米孔测序 (视频) | JoVE Journal
Uploaded: 2020-03-11T14:30:00.000+00:00
Duration: 5 min 45 s
Description: 了解基于纳米孔测序平台的乌苏图病毒全基因组扩增子测序验证方案。该方法使用MinION和MinIT进行测序，通过Porechop、cutadapt、Minimap2和Samtools进行数据处理，评估错误率并生成一致性序列，为实时疫情支持提供可靠的病毒基因组数据。

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

该项工作获得了欧盟“地平线2020”研究与创新计划的资助，资助协议号为643476（COMPARE）。

注意：将使用的软件工具列表：usearch v11.0.667；muscle v3.8.1551；porechop 0.2.4；cutadapt 2.5；minimap2 2.16-r922；samtools 1.9；trimmomatic 0.39；bbmap 38.33；spades v3.13.1；kma-1.2.8
1. 引物设计
<ol>
	<li>首先从公共或私人数据库中下载或获取一组相关的参考全基因组序列。例如，从NCBI数据库中获取所有全长的USUV基因组（taxid64286）22。USUV编码的基因组长度约为11,000个核苷酸，因此仅获取序列长度在8,000至12,000个核苷酸之间的序列。使用以下搜索条目进行检索： 
	- taxid64286[Organism:noexp] AND 8000[SLEN]:12000[SLEN]。
	<ol>
		<li>点击 Send to | Complete Record | File；选择格式为FASTA并生成文件。</li>
	</ol>
	</li>
	<li>为减少参考序列的数量，需从数据集中去除重复序列或核苷酸序列相似性超过99%的序列。可使用usearch工具中的cluster_fast选项完成此操作23。在命令行中输入： 
	- usearch -cluster_fast All_USUV.fasta -id 0.99 -centroids All_USUV_dedup.fasta</li>
	<li>为设计引物，需对序列进行比对，此步骤使用MUSCLE软件完成24。在命令行中输入： 
	- muscle -in All_USUV_dedup.fasta -out All_USUV_dedup_aligned.fasta -log log_muscle.txt 
	注意：必须手动检查比对结果以发现可能的错误。如有需要，可手动修正错误，并根据大多数全基因组序列的长度对两端进行裁剪。</li>
	<li>使用Primal工具初步筛选可用于全长扩增子测序的引物19。将比对好的序列上传至Primal网站（<a href="http://primal.zibraproject.org/" target="_blank">http://primal.zibraproject.org/</a>），选择理想的扩增子长度以及不同扩增子之间的重叠长度。访问primal.zibraproject.org，填写方案名称，上传比对后的FASTA文件，设置扩增子长度、重叠大小，然后生成引物方案。</li>
	<li>对所有可获得的完整USUV序列（非缩减或去重后的序列集）进行比对。在命令行中输入： 
	- muscle -in All_USUV.fasta -out All_USUV_aligned.fasta -log log_muscle.txt 
	注意：将生成的引物比对到完整的序列比对结果中（不要使用去重后的比对结果），手动修正错误，并最多包含5个简并引物位点。</li>
</ol>
2. 多重 PCR
<ol>
	<li>使用设计的引物进行多重PCR，并进行纳米孔测序和Illumina测序。针对USUV的多重PCR按照先前所述方法进行19,21。</li>
	<li>使用flip-flop版本3.0.6.6+9999d81进行碱基识别。</li>
</ol>
3. 利用纳米孔数据生成共识序列的数据分析
<ol>
	<li>可在单次纳米孔测序运行中对多个样本进行多重检测。完成测序后，对纳米孔数据进行解复用，可使用 Porechop25 实现。为防止污染并提高准确性，需使用 require_two_barcodes 参数。在命令行中输入： 
	- porechop -i Run_USUV.fastq -o Run_USUV_demultiplex --require_two_barcodes</li>
	<li>解复用后，使用 cutadapt26 去除引物序列（正向和反向序列均在文件 Primers_Usutu.fasta 中标明）。此外，去除长度短于 75 个核苷酸的序列。必须去除引物序列，因其可能在共识序列中引入人为偏差。在命令行中输入： 
	- cutadapt -b file:Primers_USUV.fasta -o BC01_trimmed.fastq BC01.fastq -m 75</li>
	<li>可使用 minimap227 将解复用后的测序读段比对至一组不同的参考毒株，并使用 samtools28 生成共识序列。以下示例展示了一个样本（BC01）基于参考序列的比对及共识序列生成流程。在命令行中输入： 
	- minimap2 -ax map-ont Random_Refs_USUV.fasta BC01_trimmed.fastq &#62; BC01.bam 
	- samtools sort BC01.bam &#62; BC01_sorted.bam 
	- bcftools mpileup -Ou -f Random_Refs_USUV.fasta BC01_sorted.bam | bcftools call -mv -Oz -o BC01.vcf.gz 
	- bcftools index BC01.vcf.gz 
	- cat Random_Refs_USUV.fasta | bcftools consensus BC01.vcf.gz &#62; BC01_consensus.fasta</li>
	<li>在进行基于参考序列的比对时，必须使用亲缘关系密切的参考序列。因此，应对生成的共识序列进行 BlastN 搜索，以确定最接近的参考毒株。随后，使用该最近的参考毒株作为参考序列，重复进行基于参考的比对（步骤 3.3 和 3.4）。在命令行中输入： 
	- minimap2 -ax map-ont Ref_USUV_BC01.fasta BC01_trimmed.fastq &#62; BC01_ref.bam 
	- samtools sort BC01_ref.bam &#62; BC01_sorted_ref.bam 
	- bcftools mpileup -Ou -f Ref_USUV_BC01.fasta BC01_sorted_ref.bam | bcftools call -mv -Oz -o BC01_ref.vcf.gz 
	- bcftools index BC01_ref.vcf.gz 
	- cat Ref_USUV_BC01.fasta | bcftools consensus BC01_ref.vcf.gz &#62; BC01_ref_consensus.fasta</li>
</ol>
4. Illumina 数据分析
<ol>
	<li>这些序列在测序后会自动进行解复用。可使用 Trimmomatic 对读段进行质量控制。29对于成对末端Illumina测序数据，使用常用的截断值——中位PHRED得分为33，并设置最小读长为75，以获得准确且高质量的测序读段。在命令行中输入： 
	- trimmomatic PE -phred33 9_S9_L001_R1_001.fastq.gz 9_S9_L001_R2_001.fastq.gz 9_1P.fastq 9_1U.fastq 9_2P.fastq 9_2U.fastq LEADING:3 TRAILING:3 SLIDINGWINDOW:3:15 MINLEN:75</li>
	<li>使用 cutadapt 去除引物序列（引物序列信息见文件 Primers_Usutu.fasta 中的正向和反向序列），以避免引入人为偏差26此外，使用以下命令删除长度短于75个核苷酸的序列。在命令行中输入： 
	- cutadapt -b o 9_1P_trimmed.fastq -p 9_2P_trimmed.fastq 9_1P.fastq 9_2P.fastq -m 75</li>
	<li>在从头组装之前，可对序列读段进行标准化处理，以实现基因组范围内的均匀覆盖度。这一步骤至关重要，因为诸如 SPAdes 等从头组装软件在组装序列读段时会考虑读段覆盖度。使用 BBMap 软件包中的 BBNorm 将读段标准化至覆盖度为 50。30在命令行中输入： 
	- bbmap/bbnorm.sh target=50 in=9_1P_trimmed.fastq in2=9_2P_trimmed.fastq out=Sample9_FW_norm.fastq out2=Sample9_RE_norm.fastq</li>
	<li>使用 SPAdes 对标准化后的测序读长进行从头组装31使用所有不同的kmer值（21、33、55、77、99和127）进行组装时采用默认设置。在命令行中输入： 
	- spades.py -k 21,33,55,77,99,127 -o Sample9 -1 Sample9.qc.f.fq -2 Sample9.qc.r.fq</li>
	<li>使用 minimap2 以及 Geneious、Bioedit 或 Ugene 等软件将质控后的测序读段比对至获得的共识序列，以校正比对结果。必须检查重叠群（contig）的起始端和末端。
	<ol>
		<li>使用 minimap2 将质控后的测序读段与获得的共识序列进行比对。</li>
		<li>在 Geneious/Bioedit/UGene 中导入比对结果。</li>
		<li>手动检查、校正和审编基因组的起始和终止区域。</li>
	</ol>
	</li>
</ol>
5. 使用Illumina数据作为金标准，确定为补偿纳米孔测序错误特征所需的测序读长覆盖度
<ol>
	<li>选择比对到一个扩增子的序列读段，此处为扩增子26。随后使用 minimap2 将纳米孔测序读段比对至该扩增子。利用 Samtools 仅筛选比对到扩增子26的读段，并将 bam 文件转换为 fastq 格式。在命令行中输入： 
	- minimap2 -ax map-ont -m 150 Amplicon26.fasta BC01_trimmed.fastq &#62; BC01.bam 
	- samtools view -b -F 4 BC01.bam &#62; BC01_mapped.bam 
	- samtools bam2fq BC01_mapped.bam | seqtk seq - -&#62; BC01_mapped.fastq</li>
	<li>对例如200条序列读段的子集进行一千次随机抽样。例如，若将数值改为10，则将进行一千次每次随机抽取10条序列读段的子集。该脚本作为 补充文件1 提供。在命令行中输入： 
	- python Random_selection.py</li>
	<li>将所有随机选取的序列读段比对至扩增子26。使用 KMA32 进行序列读段比对，并立即生成一致性序列。采用纳米孔测序优化设置，即使用 -bcNano 标志。在命令行中输入： 
	- kma index -i Amplicon26.fasta 
	- for file in random_sample*; do 
	- sampleID=${file%.fastq} 
	- kma -i ${sampleID}.fastq -o ${sampleID} -t_db Amplicon26.fasta -mem_mode -mp 5 -mrs 0.0 -bcNano 
	- done</li>
	<li>使用以下命令行指令检查生成的一致性序列： 
	- cat *.fsa &#62; All_genomes.fsa 
	- minimap2 -ax map-ont Amplicon26.fasta All_genomes.fsa &#62; All_genomes.bam 
	- samtools sort All_genomes.bam &#62; All_genomes_sorted.bam 
	- samtools stats All_genomes_sorted.bam &#62; stats.txt
	<ol>
		<li>错误率显示在 stats.txt 文件中标题为 error rate #mismatches / bases mapped 的部分。使用以下命令在屏幕上显示： 
		- grep ^SN stats.txt | cut -f 2-</li>
		<li>插入缺失（indel）的数量显示在标题为 #Indels per cycle 的部分。使用以下命令在屏幕上显示： 
		- grep ^IC stats.txt | cut -f 2-</li>
	</ol>
	</li>
</ol>

<ol><li>Faria, N. R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Establishment+and+cryptic+transmission+of+Zika+virus+in+Brazil+and+the+Americas%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">Establishment and cryptic transmission of Zika virus in Brazil and the Americas.</a> Nature. 546 (7658), 406-410 (2017).</li>
<li>Bonaldo, M. C., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Genome+analysis+of+yellow+fever+virus+of+the+ongoing+outbreak+in+Brazil+reveals+polymorphisms%5BTitle%5D)+AND+%22Mem%C3%B3rias+do+Instituto+Oswaldo+Cruz%22%5BJournal%5D)">Genome analysis of yellow fever virus of the ongoing outbreak in Brazil reveals polymorphisms.</a> Memórias do Instituto Oswaldo Cruz. 112 (6), 447-451 (2017).</li>
<li>Faria, N. R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Genomic+and+epidemiological+monitoring+of+yellow+fever+virus+transmission+potential%5BTitle%5D)+AND+%22bioRxiv%22%5BJournal%5D)">Genomic and epidemiological monitoring of yellow fever virus transmission potential.</a> bioRxiv. , 299842(2018).</li>
<li>Magi, A., Giusti, B., Tattini, L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Characterization+of+MinION+nanopore+data+for+resequencing+analyses%5BTitle%5D)+AND+%22Briefings+in+Bioinformatics%22%5BJournal%5D)">Characterization of MinION nanopore data for resequencing analyses.</a> Briefings in Bioinformatics. 18 (6), bbw077(2016).</li>
<li>Rang, F. J., Kloosterman, W. P., de Ridder, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((From+squiggle+to+basepair%3A+computational+approaches+for+improving+nanopore+sequencing+read+accuracy%5BTitle%5D)+AND+%22Genome+Biology%22%5BJournal%5D)">From squiggle to basepair: computational approaches for improving nanopore sequencing read accuracy.</a> Genome Biology. 19 (1), 90(2018).</li>
<li>Nanopore Store, R10 flow cells. , https://store.nanoporetech.com/flowcells/spoton-flow-cell-mk-i-r10.html (2019).</li>
<li>Wick, R. R., Judd, L. M., Holt, K. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Performance+of+neural+network+basecalling+tools+for+Oxford+Nanopore+sequencing%5BTitle%5D)+AND+%22Genome+Biology%22%5BJournal%5D)">Performance of neural network basecalling tools for Oxford Nanopore sequencing.</a> Genome Biology. 20 (1), 129(2019).</li>
<li>GitHub - nanoporetech/flappie: Flip-flop basecaller for Oxford Nanopore reads. , https://github.com/nanoporetech/flappie (2019).</li>
<li>Lühken, R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Distribution+of+Usutu+Virus+in+Germany+and+Its+Effect+on+Breeding+Bird+Populations%5BTitle%5D)+AND+%22Emerging+Infectious+Diseases%22%5BJournal%5D)">Distribution of Usutu Virus in Germany and Its Effect on Breeding Bird Populations.</a> Emerging Infectious Diseases. 23 (12), 1994-2001 (2017).</li>
<li>Cadar, D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Widespread+activity+of+multiple+lineages+of+Usutu+virus%2C+Western+Europe%2C+2016%5BTitle%5D)+AND+%22Eurosurveillance%22%5BJournal%5D)">Widespread activity of multiple lineages of Usutu virus, Western Europe, 2016.</a> Eurosurveillance. 22 (4), (2017).</li>
<li>Becker, N., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Epizootic+emergence+of+Usutu+virus+in+wild+and+captive+birds+in+Germany%5BTitle%5D)+AND+%22PLoS+ONE%22%5BJournal%5D)">Epizootic emergence of Usutu virus in wild and captive birds in Germany.</a> PLoS ONE. 7 (2), (2012).</li>
<li>Diagne, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Usutu+Virus+Isolated+from+Rodents+in+Senegal%5BTitle%5D)+AND+%22Viruses%22%5BJournal%5D)">Usutu Virus Isolated from Rodents in Senegal.</a> Viruses. 11 (2), 181(2019).</li>
<li>Bakonyi, T., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Usutu+virus+infections+among+blood+donors%2C+Austria%2C+July+and+August+2017+%E2%80%93+Raising+awareness+for+diagnostic+challenges%5BTitle%5D)+AND+%22Eurosurveillance%22%5BJournal%5D)">Usutu virus infections among blood donors, Austria, July and August 2017 – Raising awareness for diagnostic challenges.</a> Eurosurveillance. 22 (41), (2017).</li>
<li>Zaaijer, H. L., Slot, E., Molier, M., Reusken, C. B. E. M., Koppelman, M. H. G. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Usutu+virus+infection+in+Dutch+blood+donors%5BTitle%5D)+AND+%22Transfusion%22%5BJournal%5D)">Usutu virus infection in Dutch blood donors.</a> Transfusion. , trf.15444(2019).</li>
<li>Cadar, D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Blood+donor+screening+for+West+Nile+virus+%28WNV%29+revealed+acute+Usutu+virus+%28USUV%29+infection%2C+Germany%2C+September+2016%5BTitle%5D)+AND+%22Eurosurveillance%22%5BJournal%5D)">Blood donor screening for West Nile virus (WNV) revealed acute Usutu virus (USUV) infection, Germany, September 2016.</a> Eurosurveillance. 22 (14), 30501(2017).</li>
<li>Pierro, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Detection+of+specific+antibodies+against+West+Nile+and+Usutu+viruses+in+healthy+blood+donors+in+northern+Italy%2C+2010%E2%80%932011%5BTitle%5D)+AND+%22Clinical+Microbiology+and+Infection%22%5BJournal%5D)">Detection of specific antibodies against West Nile and Usutu viruses in healthy blood donors in northern Italy, 2010–2011.</a> Clinical Microbiology and Infection. 19 (10), E451-E453 (2013).</li>
<li>Pecorari, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((First+human+case+of+Usutu+virus+neuroinvasive+infection%2C+Italy%2C+August-September+2009%5BTitle%5D)+AND+%22Euro+surveillance%3A+bulletin+europ%C3%A9en+sur+les+maladies+transmissibles+%3D+European+Communicable+Disease+Bulletin%22%5BJournal%5D)">First human case of Usutu virus neuroinvasive infection, Italy, August-September 2009.</a> Euro surveillance: bulletin européen sur les maladies transmissibles = European Communicable Disease Bulletin. 14 (50), (2009).</li>
<li>Simonin, Y., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Human+Usutu+Virus+Infection+with+Atypical+Neurologic+Presentation%2C+Montpellier%2C+France%2C+2016%5BTitle%5D)+AND+%22Emerging+Infectious+Diseases%22%5BJournal%5D)">Human Usutu Virus Infection with Atypical Neurologic Presentation, Montpellier, France, 2016.</a> Emerging Infectious Diseases. 24 (5), 875-878 (2018).</li>
<li>Quick, J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Multiplex+PCR+method+for+MinION+and+Illumina+sequencing+of+Zika+and+other+virus+genomes+directly+from+clinical+samples%5BTitle%5D)+AND+%22Nature+Protocols%22%5BJournal%5D)">Multiplex PCR method for MinION and Illumina sequencing of Zika and other virus genomes directly from clinical samples.</a> Nature Protocols. 12 (6), 1261-1276 (2017).</li>
<li>Quick, J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Real-time%2C+portable+genome+sequencing+for+Ebola+surveillance%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">Real-time, portable genome sequencing for Ebola surveillance.</a> Nature. 530 (7589), 228-232 (2016).</li>
<li>Oude Munnink, B. B., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Towards+high+quality+real-time+whole+genome+sequencing+during+outbreaks+using+Usutu+virus+as+example%5BTitle%5D)+AND+%22Infection%2C+Genetics+and+Evolution%22%5BJournal%5D)">Towards high quality real-time whole genome sequencing during outbreaks using Usutu virus as example.</a> Infection, Genetics and Evolution. 73, 49-54 (2019).</li>
<li>Benson, D. A., Karsch-Mizrachi, I., Lipman, D. J., Ostell, J., Sayers, E. W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((GenBank%5BTitle%5D)+AND+%22Nucleic+Acids+Research%22%5BJournal%5D)">GenBank.</a> Nucleic Acids Research. 38 (Database issue), D46-D51 (2010).</li>
<li>Edgar, R. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Search+and+clustering+orders+of+magnitude+faster+than+BLAST%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">Search and clustering orders of magnitude faster than BLAST.</a> Bioinformatics. 26 (19), 2460-2461 (2010).</li>
<li>Edgar, R. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((MUSCLE%3A+multiple+sequence+alignment+with+high+accuracy+and+high+throughput%5BTitle%5D)+AND+%22Nucleic+Acids+Research%22%5BJournal%5D)">MUSCLE: multiple sequence alignment with high accuracy and high throughput.</a> Nucleic Acids Research. 32 (5), 1792-1797 (2004).</li>
<li>GitHub - rrwick/Porechop: adapter trimmer for Oxford Nanopore reads. , https://github.com/rrwick/porechop (2018).</li>
<li>Martin, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Cutadapt+removes+adapter+sequences+from+high-throughput+sequencing+reads%5BTitle%5D)+AND+%22EMBnet.journal%22%5BJournal%5D)">Cutadapt removes adapter sequences from high-throughput sequencing reads.</a> EMBnet.journal. 17 (1), 10(2011).</li>
<li>Li, H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Minimap2%3A+pairwise+alignment+for+nucleotide+sequences%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">Minimap2: pairwise alignment for nucleotide sequences.</a> Bioinformatics. 34 (18), 3094-3100 (2018).</li>
<li>Li, H., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+Sequence+Alignment%2FMap+format+and+SAMtools%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">The Sequence Alignment/Map format and SAMtools.</a> Bioinformatics. 25 (16), 2078-2079 (2009).</li>
<li>Bolger, A. M., Lohse, M., Usadel, B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Trimmomatic%3A+a+flexible+trimmer+for+Illumina+sequence+data%5BTitle%5D)+AND+%22Bioinformatics+%28Oxford%2C+England%29%22%5BJournal%5D)">Trimmomatic: a flexible trimmer for Illumina sequence data.</a> Bioinformatics (Oxford, England). 30 (15), 2114-2120 (2014).</li>
<li>BBMap download | SourceForge.net. , https://sourceforge.net/projects/bbmap/ (2019).</li>
<li>Bankevich, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((SPAdes%3A+a+new+genome+assembly+algorithm+and+its+applications+to+single-cell+sequencing%5BTitle%5D)+AND+%22Journal+of+Computational+Biology%22%5BJournal%5D)">SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing.</a> Journal of Computational Biology. 19 (5), 455-477 (2012).</li>
<li>Clausen, P. T. L. C., Aarestrup, F. M., Lund, O. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Rapid+and+precise+alignment+of+raw+reads+against+redundant+databases+with+KMA%5BTitle%5D)+AND+%22BMC+Bioinformatics%22%5BJournal%5D)">Rapid and precise alignment of raw reads against redundant databases with KMA.</a> BMC Bioinformatics. 19 (1), 307(2018).</li>
<li>Brinkmann, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Proficiency+Testing+of+Virus+Diagnostics+Based+on+Bioinformatics+Analysis+of+Simulated+In+Silico+High-Throughput+Sequencing+Data+Sets%5BTitle%5D)+AND+%22Journal+of+Clinical+Microbiology%22%5BJournal%5D)">Proficiency Testing of Virus Diagnostics Based on Bioinformatics Analysis of Simulated In Silico High-Throughput Sequencing Data Sets.</a> Journal of Clinical Microbiology. 57 (8), (2019).</li>
</ol>

作者无任何利益冲突需要披露。

纳米孔测序技术不断发展，因此需要方法来监测其错误率。本文介绍了一种用于监测纳米孔测序仪错误率的工作流程。该方法在新型流动池发布后，或碱基识别软件有新版本发布时尤为有用。此外，对于希望建立并验证自身测序方案的用户而言，该方法同样具有实用价值。
不同的软件和比对工具可能会产生不同的结果33。在本研究中，我们力求使用常见且可免费获取的软件包，这些软件具有清晰的文档说明。在某些情况下，可能会优先选择商业工具，这类工具通常具有更友好的用户界面，但需要付费使用。未来，如果测序技术或碱基识别软件发生重大更新，该方法可再次应用于同一样本。理想情况下，应在每次碱基识别软件或流通池更新后进行此操作，但鉴于当前技术发展的速度，也可仅在重大更新后进行。
测序错误率的降低使得可 multiplex（多重检测）的样本数量增加。因此，纳米孔测序正逐步接近取代传统的实时PCR用于诊断检测，流感病毒的诊断已实现这一替代。此外，错误率的降低也提高了该技术的可用性，例如用于稀有变异株的鉴定以及高通量、无偏倚的宏基因组测序。
该实验方案中的一个关键步骤是必须具备可靠且高度相近的参考序列。引物的设计基于当前对病毒多样性的认知，可能需要定期更新。在建立基于扩增子的测序方法时，另一个关键点是平衡引物浓度，以获得均匀的扩增子深度。这使得在一次测序运行中可对更多样本进行多重分析，从而显著降低检测成本。

第三代测序技术的快速发展使我们能够在病毒暴发期间接近实时地进行测序。及时获得的遗传信息有助于确定病毒病原体的来源和进化过程。然而，下一代测序领域的金标准仍然是第二代测序仪。这些技术依赖于特定且耗时的步骤，例如乳液PCR中的克隆扩增或克隆桥接扩增。第三代测序仪成本更低、可手持操作，并且具有简化的文库制备方法。特别是测序设备的小型化和较低的购买价格，使其成为可部署、适用于现场测序的理想选择。例如，在塞拉利昂的埃博拉病毒暴发期间以及在巴西正在进行的虫媒病毒暴发调查中已体现出其应用价值1,2,3。然而，据报道，纳米孔测序存在较高的错误率4，这可能会限制其应用范围。
纳米孔测序技术正在快速发展，市场上定期推出新产品。例如，1D²试剂盒能够对DNA分子的两条链进行测序，从而提高碱基识别的准确性5，以及R10流动池的开发，可在纳米孔内两个不同时间点测量电流变化6。此外，生物信息学工具的改进，如碱基识别算法的优化，也将提升碱基识别的准确性7。其中最常用的碱基识别软件之一（如Albacore）在9个月内至少更新了12次5。最近，制造商还发布了一种名为flip-flop的新型碱基识别器，并已集成至默认的纳米孔测序软件中8。这些改进的综合效应将带来更准确的测序结果，并降低纳米孔测序仪的错误率。
乌苏图病毒（USUV）是一种通过蚊子传播的节肢动物传播病毒，属于黄病毒科（Flaviviridae），其基因组为约11,000个核苷酸的单股正链RNA。USUV主要影响灰林鸮和乌鸫9,10，尽管其他鸟类物种也易感USUV感染11。最近，USUV也在啮齿类动物和鼩鼱中被发现，但它们在病毒传播中的潜在作用尚不清楚12。在人类中，已在献血者中发现无症状感染13,14,15,16，同时也有报道指出USUV感染与脑炎或脑膜脑炎相关17,18。在荷兰，USUV于2016年首次在野生鸟类中被检测到10，并于2018年在无症状献血者中发现14。自首次发现USUV以来，此后数年均报告了疫情暴发，目前正在进行包括全基因组测序在内的监测工作，以追踪该节肢动物传播病毒在既往无免疫人群中出现及传播的情况。
与其他病毒（如埃博拉病毒、寨卡病毒和黄热病毒）已报道的方法类似3,19,20，我们开发了一套引物用于对USUV全长基因组进行测序21。这种基于聚合酶链式反应（PCR）的方法能够从宿主污染程度较高的样本（如脑组织样本）中回收全长USUV基因组，适用于Ct值高达约32的样本。基于扩增子的测序方法具有比宏基因组测序更高的灵敏度和特异性。然而，该方法的局限性在于：为覆盖所有毒株，所设计引物需基于序列相似性，且引物设计依赖于当前对病毒多样性的认知水平。
鉴于第三代测序技术的持续发展与改进，有必要定期评估测序仪的错误率。本文以 Usutu 病毒（USUV）为例，介绍一种将纳米孔测序结果直接与 Illumina 测序结果进行比对以评估其性能的方法。该方法应用于使用最新 R10 流通池生成的序列，并采用最新版本的翻转-翻转（flip-flop）碱基识别算法进行碱基识别。

<table><tbody><tr><td>Agencourt AMPure XP 磁珠</td><td>Beckman Coulter</td><td>A63881</td><td></td></tr><tr><td>dNTP 混合物</td><td>Qiagen</td><td>201900</td><td></td></tr><tr><td>FLO-MIN106 R10 流式芯片</td><td>Nanopore</td><td>R10 flowcell</td><td></td></tr><tr><td>KAPA Hyperplus 文库制备试剂盒</td><td>Roche</td><td>7962436001</td><td></td></tr><tr><td>文库加载磁珠试剂盒</td><td>Nanopore</td><td>EXP-LLB001</td><td></td></tr><tr><td>1D 连接测序试剂盒</td><td>Nanopore</td><td>SQK-LSK109</td><td></td></tr><tr><td>1D 原始条形码试剂盒 1-12</td><td>Nanopore</td><td>EXP-NBD103</td><td></td></tr><tr><td>1D 原始条形码试剂盒 13-24</td><td>Nanopore</td><td>EXP-NBD104</td><td></td></tr><tr><td>NEB 平末端/TA 连接酶主混合液</td><td>NEB</td><td>M0367S</td><td></td></tr><tr><td>NEB Next 快速连接模块</td><td>NEB</td><td>E6056</td><td></td></tr><tr><td>NEB Next Ultra II 末端修复 / dA 加尾模块</td><td>NEB</td><td>E7546S</td><td></td></tr><tr><td>Protoscript II 逆转录酶</td><td>NEB</td><td>M0368X</td><td></td></tr><tr><td>Q5 高保真聚合酶</td><td>NEB</td><td>M0491</td><td></td></tr><tr><td>Qubit dsDNA HS 检测试剂盒</td><td>Thermo Fisher</td><td>Q32851</td><td></td></tr><tr><td>随机引物</td><td>Promega</td><td>C1181</td><td></td></tr><tr><td>RNAsin 核糖核酸酶抑制剂</td><td>Promega</td><td>N2111</td><td></td></tr></tbody></table>

以乌苏图病毒为例验证全基因组纳米孔测序

Validating Whole Genome Nanopore Sequencing, using Usutu Virus as an Example

全基因组测序可用于病毒暴发的特征分析与溯源。已有多种基于纳米孔测序平台的全基因组测序方案被用于不同病毒的研究。这些方法采用基于重叠扩增子的策略，可靶向特定病毒或一组遗传关系相近的病毒。除了确认病毒的存在外，测序还可用于基因组流行病学研究，追踪病毒传播路径，并揭示其来源、宿主库及传播方式。对于此类应用，理解所用测序平台相关错误率的潜在影响至关重要。在临床和公共卫生领域的常规应用中，每次对方案进行重要修改时，都必须对这些影响进行记录。此前，已通过与Illumina测序直接比较的方法，验证了在纳米孔测序平台（R9.4流动池）上进行乌苏图病毒全基因组测序的方案。本文以R10流动池与Illumina测序的比较为例，描述了确定所需测序读段覆盖深度的方法。

最近，新型流动槽版本（R10）已发布，该版本改进了用于将电子电流信号转换为DNA序列的碱基识别软件（即所谓的“翻转-翻转”碱基识别器）。因此，我们对先前在R9.4流动槽和Illumina MiSeq测序仪上测序过的一只USUV阳性猫头鹰脑组织中的USUV进行了重新测序21。本文描述了通过与Illumina测序结果直接比较，确定可靠获得一致序列所需测序深度的方法。
通过使用新型流动池并结合翻转-翻转（flip-flop）碱基识别方法，我们证明，40x的测序深度可获得与Illumina测序相当的完全一致结果。30x测序深度的错误率为0.0002%，相当于每测序585,000个核苷酸出现一个错误；而20x测序深度则为每测序63,529个核苷酸出现一个错误。10x测序深度下，每测序3,312个核苷酸出现一个错误，这意味着在完整的USUV基因组中，平均有超过三个核苷酸被错误识别。当测序深度高于30x时，未观察到任何插入缺失（indel）。20x测序深度检测到一个indel位点，而10x测序深度则在29个位点上检测到indel。不同测序深度阈值下的错误率概览见表1。
<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always">
	<tbody>
		<tr>
			<td>覆盖度</td>
			<td>第一次迭代错误数</td>
			<td>第一次迭代错误率</td>
			<td>插入缺失：</td>
			<td>第二次迭代错误数</td>
			<td>第二次迭代错误率</td>
			<td>插入缺失：</td>
			<td>第三次迭代错误数</td>
			<td>第三次迭代错误率</td>
			<td>插入缺失：</td>
		</tr>
		<tr>
			<td>10×</td>
			<td>100</td>
			<td>0.0274%</td>
			<td>4</td>
			<td>116</td>
			<td>0.0297%</td>
			<td>18</td>
			<td>110</td>
			<td>0.0282%</td>
			<td>7</td>
		</tr>
		<tr>
			<td>20×</td>
			<td>4</td>
			<td>0.0010%</td>
			<td>0</td>
			<td>6</td>
			<td>0.0015%</td>
			<td>1</td>
			<td>7</td>
			<td>0.0018%</td>
			<td>0</td>
		</tr>
		<tr>
			<td>30x</td>
			<td>2</td>
			<td>0.0005%</td>
			<td>0</td>
			<td>0</td>
			<td>0.0000%</td>
			<td>0</td>
			<td>0</td>
			<td>0.0000%</td>
			<td>0</td>
		</tr>
		<tr>
			<td>40x</td>
			<td>0</td>
			<td>0.0000%</td>
			<td>0</td>
			<td>0</td>
			<td>0.0000%</td>
			<td>0</td>
			<td>0</td>
			<td>0.0000%</td>
			<td>0</td>
		</tr>
		<tr>
			<td>50×</td>
			<td>0</td>
			<td>0.0000%</td>
			<td>0</td>
			<td>0</td>
			<td>0.0000%</td>
			<td>0</td>
			<td>0</td>
			<td>0.0000%</td>
			<td>0</td>
		</tr>
	</tbody>
</table>
表1：纳米孔测序错误率概述。 每次迭代代表一千个随机样本。
补充文件 1：随机选择。 <a href="https://www.jove.com/files/ftp_upload/60906/Random_selection.zip" target="_blank">请点击此处查看该文件（右键点击可下载）。</a>

我们此前已验证了一种基于扩增子的乌苏图病毒（USUV）全基因组测序方案，该方案在纳米孔测序平台上进行。本文中，我们对所使用的方法进行了更详细的描述，并确定了纳米孔R10流通池的错误率。

全基因组测序可用于表征和追踪病毒暴发。基于纳米孔的全基因组测序方案已被描述用于……

我们利用病毒基因测序来探究病毒的来源及其传播方式。为了能够实时使用这些数据（例如在疫情暴发期间），我们必须准确了解所提供数据的精确程度。这些技术的主要优势在于成本相对较低、速度快，并且能够实现实时碱基编码，可用于现场测序。该方法可用于及时生成序列数据，例如用于确定多种病原体的来源。对于初次使用者而言，该技术可能较为复杂，需要具备基本的 Unix 知识，尤其是在安装不同软件包时。本方法的可视化将帮助用户确定其特定测序方案的错误率，并解答科研问题。测序运行在连接至 MinIT 的 MinION 上启动。在 Nanopore 平台上完成测序后，将数据加载到 Porechop 中。为防止污染并提高准确性，使用 require_two_barcodes 标志，并按照所示命令输入。去多重后，使用 cutadapt 去除引物序列以及长度短于 75 个核苷酸的序列，具体命令如所示。使用 Minimap2 将去多重化的序列读段比对到不同参考株的面板上。然后使用 Samtools 生成一致性序列。对于基于参考序列的比对，必须使用生成的一致性序列进行 blast 搜索，以鉴定最接近的参考株。然后以最接近的参考菌株作为参考序列，重复基于参考的比对。为了确定为补偿Nanopore测序错误谱所需的测序深度，选择比对到一个扩增子的序列，并使用Minimap2将Nanopore测序读段比对至该扩增子。使用Samtools仅筛选比对到Amplicon26的读段，并按照所示命令将bam文件转换为fastq格式。例如，随机选择 200 条序列读长的子集，重复一千次。将所有随机选出的序列读长比对至 Amplicon26。使用 KMA 软件，依据 BC nano 标志所示的 Nanopore 测序优化参数，对序列读长进行比对并立即生成一致性序列。要检查生成的共识序列，请使用所示命令。要显示统计信息文件（stats dot txt）中“错误率：错配碱基/比对”标题下的错误率，请使用所示命令。要显示“每循环插入缺失数”标题下的插入缺失数量，请使用所示命令。使用新型流动池结合碱基颜色翻转技术，40倍测序深度的结果与Illumina测序结果一致。30倍测序深度的错误率为0.02%，相当于每测序585,000个核苷酸出现一个错误。而20倍测序深度则相当于每测序63,529个核苷酸出现一个错误。10 倍的测序深度会导致每 3,312 个测序的核苷酸中出现一个错误，这意味着每四个乌苏图病毒基因组中就有超过 3 个核苷酸被错误识别。当测序深度超过 30 倍时，则不会观察到插入缺失突变（indels）。20倍的测序深度可检测到1个插入缺失位点，而10倍的测序深度则可检测到29个插入缺失位点。需要牢记最重要的一点是，尽管软件更新并不总是显而易见，但其更新频繁，可能对结果产生影响。目前有多种可用的软件选择。我们已经演示了最常用的程序。但不同的软件以及同一软件的不同版本可能会导致不同的结果。我们向您展示了一个示例，说明如何检查我们生成的序列数据的质量。这对于可靠的疫情应对至关重要。

validating-whole-genome-nanopore-sequencing-using-usutu-virus-as-an

Research

科研

JoVE Journal

Genetics

遗传学

以乌苏图病毒为例，验证全基因组纳米孔测序

本文描述了分离和鉴定RNA病毒RNA聚合酶保真性变异体所需的步骤，以及如何利用突变频率数据在组织培养中确认保真性的改变。...

RNA病毒保真性变异体的分离及病毒突变频率的表征

单病毒基因组学（SVG）是一种用于分离和扩增单个病毒颗粒基因组的方法。通过流式细胞术将混合病毒悬液分选到含有琼脂糖的载玻片微孔中，从而捕获单个病毒颗粒，并在后续处理过程中减少基因组剪切。利用多重置换扩增（MDA）实现全基因组扩增，获得适用于测序的基因组材料。...

单病毒基因组学用于单个病毒颗粒的分离与基因组分析

本文描述了一种用于检测循环肿瘤DNA（ctDNA）中低频突变的技术——ER-Seq。该方法的独特之处在于采用了双向错误校正、特殊的背景过滤以及高效的分子富集。...

利用下一代测序技术检测循环肿瘤DNA中的稀有突变

下一代测序（NGS）是一种强大的基因组表征工具，但其应用受到该平台较高错误率（约0.5–2.0%）的限制。我们介绍了错误校正测序的方法，可消除NGS的错误率影响，检测到低至0.0001变异等位基因分数的突变。...

利用纠错型DNA和RNA测序进行稀有事件检测

本文介绍了一种鉴定具有双链DNA基因组植物病毒的新方法。我们采用标准方法从感染的叶片中提取DNA和RNA，并进行高通量测序。生物信息学工具将测序结果组装成重叠群，识别代表病毒基因组的重叠群，并将这些基因组归类到相应的分类群中。...

结合粗提病毒颗粒中的DNA分析与感病叶片中的RNA分析以发现新病毒基因组

全长单个前病毒测序（FLIPS）是一种高效、高通量的方法，可用于扩增和测序单个接近全长的（完整和有缺陷的）HIV-1前病毒，并可评估其潜在的复制能力。FLIPS克服了以往用于测序潜伏HIV-1储存库的检测方法的局限性。...

用于下一代测序的近全长HIV-1前病毒扩增

纳米孔测序是一种新型技术，可在偏远地区和资源匮乏的环境中实现低成本测序。本文介绍了一种适用于此类条件的全血 mRNA...

使用纳米孔测序技术对全血 mRNA 进行测序

我们提供了一种通用的实验方案，可用于筛选仅结合具有感染性的病毒、而不结合因消毒方法而失去感染性的病毒或其他类似病毒的适配体。这为在便携式快速检测中确定病毒的感染状态提供了可能。...

体外筛选适配体以区分感染性与非感染性病毒

本文介绍了一种利用纳米孔技术对狂犬病病毒（RABV）基因组进行快速且经济高效的表征工作流程。该工作流程旨在支持基于基因组信息的本地化监测，提供有关流行RABV谱系及其在区域系统发育树中位置的信息，以指导狂犬病防控措施的制定。...

利用纳米孔技术进行狂犬病病毒快速全基因组测序

插入突变是正向遗传学中用于鉴定功能性基因组元件的重要工具。本文中，我们描述了一种基于插入的功能元件与转录本筛选方法（InSET），该方法可用于检测基于慢病毒的插入突变细胞文库中的慢病毒整合位点。...

插入突变细胞文库中功能相关慢病毒整合位点的鉴定

了解基于纳米孔测序平台的乌苏图病毒全基因组扩增子测序验证方案。该方法使用MinION和MinIT进行测序，通过Porechop、cutadapt、Minimap2和Samtools进行数据处理，评估错误率并生成一致性序列，为实时疫情支持提供可靠的病毒基因组数据。

以乌苏图病毒为例验证全基因组纳米孔测序

探索更多视频