Name: Validando o sequenciamento de nanoporos de genoma inteiro, usando o vírus Usutu como exemplo (Vídeo) | JoVE Journal
Uploaded: 2020-03-11T14:30:00.000+00:00
Duration: 5 min 45 s
Description: Aprenda a validar o sequenciamento de nanoporos de genoma inteiro usando o vírus Usutu como modelo. Este protocolo detalhado avalia taxas de erro e aprimora a precisão para suporte em surtos de vírus.

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

Este trabalho recebeu financiamento do programa de pesquisa e inovação Horizon 2020 da União Europeia o acordo de subvenção nº 643476 (COMPARE).

NOTA: Lista de ferramentas de software a serem utilizadas: usearch v11.0.667; músculo v3.8.1551; porechop 0.2.4; cutadapt 2.5; minimap2 2.16-r922; samtools 1.9; trimmomatic 0,39; bbmap 38.33; espadas v3.13.1; kma-1.2.8
1. Design de primer
<ol><li>Comece com o download ou recuperação de um conjunto de seqüências de genoma suscitadas de referência relevantes de coleções de dados públicas ou privadas. Por exemplo, recupere todos os genomas USUV de comprimento total (taxid64286) do banco de dados NCBI22. USUV codifica um genoma de cerca de 11.000 nucleotídeos, então só recupere as seqüências com um comprimento de seqüência de 8.000-12.000 nucleotídeos. Faça isso usando a seguinte entrada de pesquisa: - taxid64286[Organism:noexp] E 8000[SLEN]:12000[SLEN].
	<ol>
<li>Clique em Enviar para | Registro completo | Arquivo; usar formato = FASTA e criar o Arquivo.</li>
	</ol>
</li>
	<li>Para reduzir o tamanho do conjunto de seqüências de referência, remova seqüências duplicadas ou seqüências com mais de 99% de identidade de nucleotídeos do conjunto de dados. Faça isso usando a opção cluster fast do usearch23. Na linha de comando digite: - usearch -cluster_fast All_USUV.fasta -id 0,99 -centroids All_USUV_dedup.fasta
</li>
	<li>Para gerar os primers, as sequências precisam ser alinhadas. Isso é feito usando MUSCLE24. Na linha de comando digite: - músculo -em All_USUV_dedup.fasta -out All_USUV_dedup_aligned.fasta -log log_muscle.txt NOTA: É essencial inspecionar manualmente o alinhamento para verificar se há discrepâncias. Estes podem ser corrigidos manualmente se necessário e as extremidades podem ser aparadas de acordo com o comprimento da maioria das seqüências de genoma inteiros.</li>
	<li>Primal é usado para fazer uma seleção de rascunho dos primers que podem ser usados para seqüência de amplicon de comprimento total19. Carregue o alinhamento para o site primitivo<a href="http://primal.zibraproject.org/" target="_blank">(http://primal.zibraproject.org/)</a>e selecione o comprimento de amplicon preferido e sobreponha o comprimento entre os diferentes amplicons. Vá para primal.zibraproject.org, preencha o nome Esquema,carregue o arquivo fasta alinhado, selecione o comprimento do amplicon, sobreponha o tamanho e gere o esquema.</li>
	<li>Alinhe o conjunto completo de seqüências USUV completas disponíveis (não o conjunto reduzido ou deduplicado). Na linha de comando digite: - músculo -em All_USUV.fasta -out All_USUV_aligned.fasta -log log_muscle.txt NOTA: Mapeie os primers gerados contra o alinhamento completo (não use o alinhamento duplicado), corrija manualmente erros e inclua um máximo de 5 posições degenerativas de primer.</li>
</ol>2. Multiplex PCR
<ol><li>Execute o MULTIPLEX PCR usando os primers projetados e sequenciamento de nanoporos e Illumina. O PCR multiplex para USUV foi realizado como descritoanteriormente 19,21.</li>
	<li>Execute a chamada base com a versão 3.0.6.6+9999d81.</li>
</ol>3. Análise de dados para gerar seqüências de consenso a partir de dados de nanoporos
<ol><li>Várias amostras podem ser multiplexadas em uma única seqüência de nanoporos. Depois de executar a sequência de execução, demultiplexe os dados de nanoporos. Use Porechop25 para isso. Para evitar contaminação e aumentar a precisão, utilize a bandeira require_two_barcodes. Na linha de comando digite: - porechop -i Run_USUV.fastq -o Run_USUV_demultiplex --require_two_barcodes
</li>
	<li>Após desmultiplexação, remova seqüências de primer (indicadas no arquivo Primers_Usutu.fasta em ambas as orientações) usando cutadapt26. Além disso, remova sequências com um comprimento menor que 75 nucleotídeos. Os primers devem ser removidos, pois podem introduzir vieses artificiais na seqüência de consenso. Na linha de comando digite: - cutadapt -b arquivo:Primers_USUV.fasta -o BC01_trimmed.fastq BC01.fastq -m 75
</li>
	<li>As leituras de seqüência demultiplexadas podem ser mapeadas contra um painel de diferentes cepas de referência usando minimap227 e uma seqüência de consenso pode ser gerada usando samtools28. Siga o exemplo abaixo que mostra o procedimento de um alinhamento baseado em referência e a geração de seqüência de consenso de uma amostra: BC01. Na linha de comando digite: - minimap2 -ax map-ont Random_Refs_USUV.fasta BC01_trimmed.fastq > BC01.bam - samtools classificar BC01.bam > BC01_sorted.bam - bcftools mpileup -Ou -f Random_Refs_USUV.fasta BC01_sorted.bam | bcftools call -mv -Oz -o BC01.vcf.gz - bcftools índice BC01.vcf.gz - gato Random_Refs_USUV.fasta | bcftools consenso BC01.vcf.gz > BC01_consensus.fasta
</li>
	<li>Para alinhamentos baseados em referência, é essencial que uma seqüência de referência intimamente relacionada seja usada. Portanto, realize uma pesquisa BlastN com a seqüência de consenso gerada para identificar a cepa de referência mais próxima. Depois disso, repita o alinhamento baseado em referência com a linha gemida de referência mais próxima como referência (etapa 3.3 e 3.4). Na linha de comando digite: - minimap2 -ax map-ont Ref_USUV_BC01.fasta BC01_trimmed.fastq > BC01_ref.bam - samtools classificar BC01_ref.bam > 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 - índice bcftools BC01_ref.vcf.gz - gato Ref_USUV_BC01.fasta | consenso bcftools BC01_ref.vcf.gz > BC01_ref_consensus.fasta
</li>
</ol>4. Análise dos dados da Illumina
<ol><li>Essas seqüências são automaticamente demultiplexadas após o sequenciamento. As leituras podem ser controladas pela qualidade usando trimmomatic29. Para seqüências de Illumina de extremidade emparelhada, use a pontuação mediana de phred de corte comumente usada de 33 e um comprimento mínimo de leitura de 75 para obter leituras precisas e de alta qualidade. Na linha de comando digite: - 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 MINLEN:75
</li>
	<li>Remova os primers (indicados no arquivo Primers_Usutu.fasta em ambas as orientações), pois eles podem introduzir vieses artificiais, usando cutadapt26. Além disso, remova seqüências com um comprimento menor que 75 nucleotídeos usando os comandos abaixo. Na linha de comando digite: - cutadapt -b o 9_1P_trimmed.fastq -p 9_2P_trimmed.fastq 9_1P.fastq 9_2P.fastq -m 75
</li>
	<li>Antes da montagem de novo, as leituras da seqüência podem ser normalizadas para uma cobertura uniforme em todo o genoma. Isso é essencial, pois montadoras de novo como SPAdes levam em conta a cobertura de leitura ao montar leituras de seqüência. Normalizar leituras para uma cobertura de leitura de 50 usando BBNorm do pacote BBMap30. Na linha de comando digite: - 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>As leituras normalizadas são montadas de novo usando SPAdes31. As configurações padrão são usadas para o conjunto usando todos os kmers diferentes (21, 33, 55, 77, 99 e 127). Na linha de comando digite: - spades.py -k 21,33,55,77,99,127 -o Amostra9 -1 Amostra9.qc.f.fq -2 Sample9.qc.r.fq
</li>
	<li>Mapeie o QC lê contra a seqüência de consenso obtida usando minimap2 e programas como Geneious, Bioedit ou Ugene para fazer a curadoria do alinhamento. É importante verificar o início e o fim do contig.<ol>
<li>Alinhar as leituras do QC contra o sequenciamento de consenso obtido usando minimap2.</li>
		<li>Importe o alinhamento em Geneious/Bioedit/UGene.</li>
		<li>Inspecione manualmente, corrija e faça a curadoria especialmente o início e o fim do genoma.</li>
	</ol>
</li>
</ol>5. Determinar a cobertura de leitura necessária para compensar o perfil de erro no sequenciamento de nanoporos usando dados de Illumina como padrão-ouro
<ol><li>Selecionar seqüência lê mapeamento para um amplicon, neste caso amplicon 26. Posteriormente, mapeie as leituras de nanoporos contra este amplicon usando minimap2. Use samtools para selecionar apenas o mapeamento de leituras para amplicon 26 e converter o arquivo bam em fastq. Na linha de comando digite: - minimap2 -ax map-ont -m 150 Amplicon26.fasta BC01_trimmed.fastq > BC01.bam - samtools view -b -F 4 BC01.bam > BC01_mapped.bam - samtools bam2fq BC01_mapped.bam | seqtk seq - -> BC01_mapped.fastq
</li>
	<li>Selecionar aleatoriamente subconjuntos de, por exemplo, 200 sequências lê mil vezes. Por exemplo, alterá-lo para 10 resultará na seleção aleatória de mil vezes um subconjunto de 10 leituras de seqüência. O script é fornecido como Arquivo Suplementar 1. Na linha de comando digite: - python Random_selection.py
</li>
	<li>Todas as leituras de seqüência selecionadas aleatoriamente estão alinhadas ao amplicon 26. Use KMA32 para mapear as leituras da seqüência e gerar imediatamente uma seqüência de consenso. Use configurações otimizadas para sequenciamento de nanoporos, indicada pela bandeira -bcNano. Na linha de comando digite: - índice kma -i Amplicon26.fasta - para arquivo em random_sample*; - sampleID=${file%.fastq} - kma -i ${sampleID}.fastq -o ${sampleID} -t_db Amplicon26.fasta -mem_mode -mp 5 -mrs 0.0 -bcNano - feito
</li>
	<li>Inspecione as seqüências de consenso geradas na linha de comando usando: - gato *.fsa > All_genomes.fsa - minimap2 -ax map-ont Amplicon26.fasta All_genomes.fsa > All_genomes.bam - samtools tipo All_genomes.bam > All_genomes_sorted.bam - estatísticas de samtools All_genomes_sorted.bam > stats.txt
	<ol>
<li>A taxa de erro é exibida no stats.txt a taxa de erro de título #mismatches / bases mapeadas. Exiba-o na tela com o seguinte comando: - grep ^SN stats.txt | corte -f 2-
</li>
		<li>A quantidade de indels é exibida a posição #Indels por ciclo. Exiba-o na tela com o seguinte comando: - grep ^IC stats.txt | corte -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>

Os autores não têm nada para revelar.

O sequenciamento de nanoporos está em constante evolução e, portanto, há a necessidade de métodos para monitorar a taxa de erro. Aqui, descrevemos um fluxo de trabalho para monitorar a taxa de erro do sequenciador de nanoporos. Isso pode ser útil após o lançamento de uma nova célula de fluxo, ou se novas versões da chamada base forem lançadas. No entanto, isso também pode ser útil para usuários que desejam configurar e validar seu próprio protocolo de seqüenciamento.
Diferentes softwares e ferramentas de alinhamento podem produzir resultados diferentes33. Neste manuscrito, buscamos usar pacotes de software disponíveis gratuitamente que são comumente usados e que tenham documentação clara. Em alguns casos, a preferência pode ser dada às ferramentas comerciais, que geralmente têm interfaces mais fáceis de usar, mas têm que ser pagas. No futuro, este método pode ser aplicado à mesma amostra no caso de grandes modificações na tecnologia de seqüência ou software de chamada base são introduzidos Preferencialmente isso deve ser feito após cada atualização do chamador base ou flowcell, no entanto, dada a velocidade dos desenvolvimentos atuais, isso também pode ser feito somente após grandes atualizações.
A redução da taxa de erro no sequenciamento permite que um número maior de amostras seja multiplexada. Assim, o sequenciamento de nanoporos está cada vez mais próximo de substituir os PCRs convencionais em tempo real para ensaios diagnósticos, o que já é o caso do diagnóstico do vírus da gripe. Além disso, a redução da taxa de erro aumenta a usabilidade dessa seqüência técnica, por exemplo, para a determinação de variantes menores e para sequenciamento metagenómico imparcial de alto throughput.
Um passo crítico no protocolo é que sequências de referência próximas e confiáveis precisam estar disponíveis. Os primers são baseados no conhecimento atual sobre a diversidade de vírus e podem precisar ser atualizados de vez em quando. Outro ponto crítico ao configurar uma abordagem de seqüenciamento baseada em amplicon é o equilíbrio da concentração da cartilha para obter um equilíbrio uniforme na profundidade de amplicon. Isso permite o multiplexamento de mais amostras em uma sequência de execução e resulta em uma redução significativa de custos.

A evolução rápida das tecnologias de seqüência de terceira geração nos permite avançar para perto do sequenciamento em tempo real durante surtos virais. Essa disponibilidade oportuna de informações genéticas pode ser útil para determinar a origem e evolução dos patógenos virais. Os padrões de ouro nos campos de sequenciamento da próxima geração, no entanto, ainda são os sequenciadores de segunda geração. Essas técnicas dependem de técnicas específicas e demoradas, como amplificação clonal durante uma amplificação de emulsão pcr ou ponte clonal. Os sequenciadores de terceira geração são mais baratos, portáteis e vêm com metodologias simplificadas de preparação de bibliotecas. Especialmente o pequeno tamanho do dispositivo de seqüência e o baixo preço de compra torna-o um candidato interessante para sequenciamento implantável e em campo. Isso poderia, por exemplo, ser visto durante o surto de ebola na Serra Leoa e durante as investigações de surto de arbovirose susceno no Brasil1,2,3. No entanto, a alta taxa de erro relatada4 pode limitar os aplicativos para os quais o sequenciamento de nanoporos pode ser usado.
O sequenciamento de nanoporos está evoluindo rapidamente. Novos produtos estão disponíveis no mercado regularmente. Exemplos disso são, por exemplo, os kits 1D ao quadrado que permitem o sequenciamento de ambos os fios da molécula de DNA, aumentando assim a precisão das chamadas bases5 e o desenvolvimento da célula de fluxo R10 que mede a mudança na corrente em duas instâncias diferentes no poro6. Além disso, ferramentas bioinformática aprimoradas, como melhorias na chamada básica, melhorarão a precisão do basecall7. Um dos mais usados como basecallers (por exemplo, Albacore), foi atualizado pelo menos 12 vezes em um período de 9 meses5. Recentemente, a fabricante também lançou um novo chamado flip-flop, que é implementado no software de nanoporos padrão8. Juntos, todas essas melhorias levarão a sequências mais precisas e diminuirão a taxa de erro do sequenciador nanoporo.
O vírus Usutu (USUV) é um arbovírus transmitido por mosquitos da família Flaviviridae e tem um genoma RNA com um fio positivo de cerca de 11.000 nucleotídeos. O USUV afeta principalmente grandes corujas-cinzentas e blackbirds9,10, embora outras espécies de aves também sejam suscetíveis à infecção pelo USUV11. Recentemente, a USUV também foi identificada em roedores e megeras, embora seu potencial papel na transmissão do vírus permaneça desconhecido12. Em humanos, foram descritas infecções assintomáticas em doadores de sangue13,14,15,16 enquanto infecções usam também foram relatadas como associadas com encefalite ou meningo-encefalite17,18. Nos Países Baixos, o USUV foi detectado pela primeira vez em aves silvestres em 201610 e em doadores de sangue assintomáticos em 201814. Desde a detecção inicial do USUV, surtos foram relatados durante os anos subsequentes e a vigilância, incluindo o sequenciamento de genomas inteiros, está atualmente em andamento para monitorar o surgimento e a disseminação de um arbovírus em uma população anteriormente ingênua.
Semelhante ao que foi descrito para outros vírus, como vírus ebola, zika vírus e vírus da febre amarela3,19,20, desenvolvemos um primer definido para seqüência de comprimento total USUV21. Esta abordagem baseada em reação em cadeia de polimerase (PCR) permite a recuperação de genomas USUV de comprimento total de genomas usuv de tipos de amostras altamente contaminados por hospedeiros, como amostras cerebrais em amostras até um valor de Ct de cerca de 32. Os benefícios de uma abordagem de seqüenciamento baseada em amplicon são uma sensibilidade maior em comparação com o sequenciamento metagenómico e uma especificidade maior. As limitações do uso de uma abordagem baseada em amplicon são que as seqüências devem ser semelhantes para projetar primers que se encaixam em todas as cepas e que os primers são projetados em nosso conhecimento atual sobre a diversidade de vírus.
Dada a evolução constante e melhorias no sequenciamento de terceira geração, há a necessidade de avaliar a taxa de erro do sequenciador regularmente. Aqui, descrevemos um método para avaliar o desempenho do nanoporo diretamente contra o sequenciamento de Illumina usando USUV como exemplo. Este método é aplicado a seqüências geradas com a célula de fluxo R10 mais recente e a chamada base é realizada com a versão mais recente do chamador base flip-flop.

<table><tbody><tr><td>Grânulos Agencourt AMPure XP</td><td>Beckman Coulter</td><td>A63881</td><td></td></tr><tr><td>dNTPs</td><td>Qiagen</td><td>201900</td><td></td></tr><tr><td>FLO-MIN106 R10 flowcell</td><td>Nanopore</td><td>R10 flowcell</td><td></td></tr><tr><td>KAPA Hyperplus kit de preparação libarária</td><td>Roche</td><td>7962436001</td><td></td></tr><tr><td>Library Loading Bead Kit</td><td>Nanopore</td><td>EXP-LLB001</td><td></td></tr><tr><td>Ligation Kit de sequenciamento 1D</td><td>Nanopore</td><td>SQK-LSK109</td><td></td></tr><tr><td>Kit de código de barras nativo 1D 1-12</td><td>Nanopore</td><td>EXP-NBD103</td><td></td></tr><tr><td>Kit de código de barras nativo 1D 13-24</td><td>Nanopore</td><td>EXP-NBD104</td><td></td></tr><tr><td>NEB Blunt/TA Ligase Master Mix</td><td>NEB</td><td>M0367S</td><td></td></tr><tr><td>NEB Next Módulo de ligação rápida</td><td>NEB</td><td>E6056</td><td></td></tr><tr><td>NEB Next Ultra II End Repair / dA-Tailing Módulo</td><td>NEB</td><td>E7546S</td><td></td></tr><tr><td>Protoscript II Transcriptase Reversa</td><td>NEB</td><td>M0368X</td><td></td></tr><tr><td>Q5 Polimerase de Alta Fidelidade</td><td>NEB</td><td>M0491</td><td></td></tr><tr><td>Qubit dsDNA HS Kit de ensaio</td><td>Thermo Fisher</td><td>Q32851</td><td></td></tr><tr><td>Primers aleatórios</td><td>Promega</td><td>C1181</td><td></td></tr><tr><td>RNAsin Ribonuclease Inibidor</td><td>Promega</td><td>N2111</td><td></td></tr></tbody></table>

validando o sequenciamento de nanoporos de genoma inteiro, usando o vírus usutu como exemplo

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

O sequenciamento do genoma inteiro pode ser usado para caracterizar e rastrear surtos virais. Protocolos de sequenciamento de genomas baseados em nanoporos foram descritos para vários vírus diferentes. Essas abordagens utilizam uma abordagem sobreposta baseada em amplicon que pode ser usada para atingir um vírus específico ou grupo de vírus geneticamente relacionados. Além da confirmação da presença do vírus, o sequenciamento pode ser utilizado para estudos de epidemiologia genômica, para rastrear vírus e desvendar origens, reservatórios e modos de transmissão. Para tais aplicações, é fundamental entender possíveis efeitos da taxa de erro associada à plataforma utilizada. A aplicação de rotina em ambientes clínicos e de saúde pública exige que isso seja documentado com todas as mudanças importantes no protocolo. Anteriormente, um protocolo para sequenciamento de vírus de genoma inteiro usutu na plataforma de sequenciamento de nanoporos foi validado (R9.4 flowcell) por comparação direta com o sequenciamento de Illumina. Aqui, descrevemos o método utilizado para determinar a cobertura de leitura necessária, usando como exemplo a comparação entre a célula de fluxo R10 e o sequenciamento de Illumina.

Recentemente, uma nova versão da versão flow cell (R10) foi lançada e ofereceu melhorias ao chamador de base usado para converter o sinal de corrente eletrônica em seqüências de DNA (o chamado flip-flop basecaller). Portanto, sequenciamos o USUV do tecido cerebral de uma coruja usuv-positiva que foi anteriormente seqüenciada em uma célula de fluxo R9.4 e em um instrumento Illumina Miseq21. Aqui, descrevemos o método utilizado para determinar a cobertura de leitura necessária para uma chamada de consenso confiável, em comparação direta com o sequenciamento de Illumina.
Usando a célula de fluxo mais nova em combinação com o flip-flop do chamador base, mostramos que uma cobertura de leitura de 40x resulta em resultados idênticos em comparação com o sequenciamento de Illumina. Uma cobertura de leitura de 30x resulta em uma taxa de erro de 0,0002% que corresponde a um erro em cada 585.000 nucleotídeos seqüenciados, enquanto uma cobertura de leitura de 20x resulta em um erro em cada 63.529 nucleotídeos seqüenciados. Uma cobertura de leitura de 10x resulta em um erro em cada 3.312 nucleotídeos sequenciados, o que significa que mais de três nucleotídeos por genoma USUV completo estão sendo chamados de errados. Com uma cobertura de leitura acima de 30x, não foram observados indels. Uma cobertura de leitura de 20x resultou na detecção de uma posição indel, enquanto uma cobertura de leitura de 10x resultou em indels em 29 posições. Uma visão geral da taxa de erro usando diferentes cortes de cobertura de leitura é mostrada na Tabela 1.
<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always"><tbody>
<tr>
<td>Cobertura</td>
			<td>Erros de iteração 1</td>
			<td>Iteração da taxa de erro 1</td>
			<td>Indels:</td>
			<td>Erros de iteração 2</td>
			<td>Iteração da taxa de erro 2</td>
			<td>Indels:</td>
			<td>Erros de iteração 3</td>
			<td>Iteração da taxa de erro 3</td>
			<td>Indels:</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>Tabela 1: Visão geral da taxa de erro da seqüência de nanoporos. Cada iteração representa mil amostras aleatórias.
Arquivo Suplementar 1: Seleção aleatória. <a href="https://www.jove.com/files/ftp_upload/60906/Random_selection.zip" target="_blank">Clique aqui para ver este arquivo (Clique com o botão direito do mouse para baixar).</a>

Nós validamos anteriormente um protocolo para sequenciamento do vírus Usutu (USUV) baseado em amplicon em uma plataforma de sequenciamento de nanoporos. Aqui, descrevemos os métodos usados com mais detalhes e determinamos a taxa de erro da célula de fluxo Nanopore R10.

Validando o sequenciamento de nanoporos de genoma inteiro, usando o vírus Usutu como exemplo

Whole genome sequencing can be used to characterize and to trace viral outbreaks. Nanopore-based whole genome sequencing protocols have been described for ...

Usamos sequenciamento genético de vírus para tentar entender de onde vêm os vírus e como eles se espalham. E para usar esses dados em tempo real, por exemplo, durante um surto, temos que entender exatamente o quão precisos são os dados que fornecemos. As principais vantagens dessas técnicas são que é relativamente barato, é rápido e permite codificação base em tempo real que pode ser usada para sequenciamento em campo.Este método pode ser usado para a geração oportuna de dados sequenciais que, por exemplo, podem ser usados para determinar a origem de uma infinidade de patógenos. Essa técnica pode ser complexa para usuários iniciantes e é necessário conhecimento básico do Unix. Especialmente para a instalação dos diferentes pacotes de software.A visualização deste método ajudará os usuários a determinar a taxa de erro de seus protocolos específicos de sequenciamento. E para responder a perguntas de pesquisa. A sequência é iniciada em um MinION conectado ao MinIT.Depois de pré-formar uma sequência na plataforma Nanopore, carregue a carga de dados em Porechop. Para evitar contaminação e aumentar a precisão, use a demanda sublinhar dois, ressaltar a bandeira de códigos de barras e digite o comando conforme indicado. Após o des multiplexamento, use cutadapt para remover sequências e sequências de primer com um comprimento menor que 75 nucleotídeos usando o comando como indicado.Use Minimap2 para mapear leituras de sequência de de multiplex contra o painel de cepas de referência distintas usando Minimap2. E usar Samtools para gerar uma sequência de consenso. Para alinhamentos baseados em referência, é essencial pré-formar uma explosão e pesquisar com a sequência de consenso gerada para identificar a cepa de referência mais próxima.Em seguida, repita o alinhamento baseado em referência com a cepa de referência mais próxima como referência. Para determinar a cobertura de leitura necessária para compensar o perfil de erro no sequenciamento de Nanopore, selecione o mapeamento de sequência para um Amplicon e use Minimap2 para mapear as leituras de Nanopore contra este Amplicon. Use Samtools para selecionar apenas o mapeamento de leituras para Amplicon26 e converter o arquivo bam em fastq usando os comandos conforme indicado.Selecione aleatoriamente subconjuntos de, por exemplo, 200 leituras de sequência mil vezes. Todas as leituras de sequência selecionadas aleatoriamente estão alinhadas ao Amplicon26. Use o KMA para mapear as leituras de sequência e para gerar imediatamente uma sequência de consenso usando as configurações otimizadas para sequenciamento de nanopore, conforme indicado pela bandeira bc nano.Para inspecionar as sequências de consenso geradas, use os comandos conforme indicado. Para exibir a taxa de erro no ponto de estatísticas txt sob o título, bases de incompatibilidades de taxa de erro mapeadas, use o comando conforme indicado. Para exibir o número de indels é exibido sob os indels de título por ciclo use o comando conforme indicado.Usando a célula de fluxo mais recente em combinação com o flip-flop de cor base, uma cobertura de leitura de 40 vezes resulta em resultados idênticos em comparação com o sequenciamento de Illumina. Uma cobertura de leitura de 30 vezes resulta em uma taxa de erro de 0,02% que corresponde a um erro em cada 585.000 nucleotídeos sequenciados. Enquanto uma cobertura de leitura de 20 vezes resulta em um erro em cada 63, 529 nucleotídeos sequenciados.Uma leitura de cobertura de 10 vezes resulta em um erro em cada 3.312 nucleotídeos sequenciados. O que significa que mais de 3 nucleotídeos por quatro genomas do vírus Usutu estão sendo chamados de errados. Com uma cobertura de leitura acima de 30 vezes, não são observados indels.Uma leitura de cobertura de 20 vezes resulta nas detecções de uma posição indel, enquanto uma cobertura de leitura de 10 vezes resulta em indels em 29 posições. O mais importante a se lembrar é que, embora nem sempre aparente, as atualizações no software ocorrem com frequência e podem influenciar os resultados. Existem várias opções de softwares para usar.Demonstramos os programas mais usados. Mas diferentes softwares e versões diferentes desses softwares podem resultar em resultados diferentes. Mostramos um exemplo de como verificamos a qualidade dos dados de sequência que geramos.E isso é fundamental para um suporte confiável para surtos.

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

Research

JoVE Journal

Genetics

Genética

O presente artigo descreve os passos necessários para isolar e caracterizar variantes RNA polimerase fidelidade de vírus RNA e como usar dados de freqüência de mutação para confirmar as alterações fidelidade em cultura de tecidos.

Isolamento de variantes de vírus de RNA Fidelity e Caracterização de frequência de mutação de vírus

Genomics vírus individuais (SVG) é um método para isolar e amplificar os genomas de virons individuais. Suspensões virais de uma assembléia mista são classificados usando citometria de fluxo em uma lâmina de microscópio com poços discretos contendo agarose, capturando assim, o virion e reduzindo corte genoma durante o processamento downstream. Amplificação do genoma inteiro é conseguida usando vários amplificação por deslocamento (MDA), resultando em material genómico que é adequado para a...

Isolamento e Análise do Genoma do vírus da única usando &#39;Genomics vírus simples&#39;

Este manuscrito descreve uma técnica para a detecção de mutações de baixa frequência em ctDNA, ER-Seq Este método é diferenciado pela sua utilização exclusiva de correção de erros de bi-direcional, um filtro de fundo especial e eficiente aquisição molecular.

Deteção de mutações raras no CtDNA usando a próxima geração de sequenciamento

Sequenciamento de próxima geração (NGS) é uma poderosa ferramenta para a caracterização do genoma que é limitada pela taxa de erro elevada da plataforma (~0.5–2.0%). Descrevemos os nossos métodos de sequenciamento de correção de erro que permitam obviar a taxa de erro NGS e detectar mutações em fracções de alelo variante tão raras como 0,0001.

Deteção de evento raro usando o DNA de correção de erro e a sequenciação do ARN

Aqui nós apresentamos uma nova abordagem para identificar a planta vírus com genoma de DNA dobro-costa. Utilizamos métodos padrão para extrair o DNA e RNA de folhas contaminadas e realizar o sequenciamento de próxima geração. Ferramentas de Bioinformatic montam sequências em contigs, identificam contigs representando genomas do vírus e atribuir genomas de grupos taxonômicos.

Combinando a análise de DNA em uma extração do Virion bruto com a análise do RNA das folhas infectadas para descobrir novos vírus de genomas

Sequenciamento proviral individual completo (aletas) fornece um método eficiente e de alta produtividade para a amplificação e sequenciamento de single, perto de provírus completos de HIV-1 (intactos e com defeito) e permite a determinação do seu potencial replicação-competência. FLIPS supera as limitações dos ensaios anteriores projetados para sequenciar o reservatório latente do HIV-1.

Amplificação de perto completos provírus HIV-1 para a próxima geração de sequenciamento

O sequenciamento de nanopore é uma nova tecnologia que permite o sequenciamento econômico em locais remotos e configurações pobres em recursos. Aqui, nós apresentamos um protocolo para arranjar em seqüência dos mRNAs do sangue inteiro que é compatível com tais circunstâncias.

Sequenciamento de mRNA de sangue total usando sequenciamento de nanoporos

Fornecemos um protocolo que pode ser geralmente aplicado a aptâmeros selecionados que se ligam apenas a vírus infecciosos e não a vírus que foram tornados não infecciosos por um método de desinfecção ou a quaisquer outros vírus semelhantes. Isso abre a possibilidade de determinar o estado de infectividade em testes portáteis e rápidos.

Seleção In Vitro de Aptamers para Diferenciar Vírus Infecciosos de Não-Infecciosos

Aqui, apresentamos um fluxo de trabalho rápido e econômico para caracterizar genomas do vírus da raiva (RABV) usando a tecnologia de nanoporos. O fluxo de trabalho destina-se a apoiar a vigilância genômica informada em nível local, fornecendo informações sobre linhagens RABV circulantes e sua colocação dentro de filogenias regionais para orientar medidas de controle da raiva.

Sequenciamento do Genoma Completo para Caracterização Rápida do Vírus Rábico Utilizando Tecnologia de Nanoporos

A mutagênese insercional é uma ferramenta essencial na genética direta para identificar elementos genômicos funcionais. Aqui, descrevemos a Tela baseada em inserção para elementos funcionais e transcritos (InSET), um método para detectar locais de integração de lentivírus dentro de uma biblioteca de células de mutagênese de inserção baseada em lentivírus.

Identificação de locais de integração de lentivírus funcionalmente relevantes em uma biblioteca de células de mutagênese insercional

Aprenda a validar o sequenciamento de nanoporos de genoma inteiro usando o vírus Usutu como modelo. Este protocolo detalhado avalia taxas de erro e aprimora a precisão para suporte em surtos de vírus.

Validando o sequenciamento de nanoporos de genoma inteiro, usando o vírus Usutu como exemplo

Explorar mais vídeos