A identificação e caracterização dos plasmídeos é de grande interesse para saúde pública, microbiologia clínica e ecologia microbiana, pois facilitam a manutenção e disseminação de genes adaptativos em ambientes ecológicos 1,5 e servem como plataformas para a evoluçãoacelerada 4. A introdução das tecnologias NGS permitiu o sequenciamento de alto rendimento de genomas de isolados microbianos e de populações, embora a identificação pós-montagem de plasmídeos completos ou de contigs derivados de plasmídeos geralmente exija sequências de leitura longa, pois dados de sequenciamento de leitura curta não conseguem resolver regiões repetidas maiores que o tamanho das leituras. O ONT tornou-se uma abordagem muito popular de sequenciamento de leitura longa, com aplicações que vão além do sequenciamento genômico completo e da metagenômica, incluindo verificação de sequências plasmídicas, detecção e caracterização de modificações epigenéticas, sequenciamento de RNA e avaliação da instabilidade intrínseca de matrizes genéticas em tandem (revisado em12,49).
Aqui, os pesquisadores apresentam um fluxo de trabalho para sequenciamento ONT de plasmídeos, projetado para situações em que a sequência cromossômica já é conhecida ou não tem interesse. A razão é que a extração de plasmídeos é trabalhosa, custosa e demorada (já que a população hospedeira precisa ser cultivada e processada), e sequenciar todo o genoma usando tecnologia de nanoporos é frequentemente suficiente para obter conjuntos de plasmídeos completos enquanto fornece seu contexto genômico. Exemplos em que o protocolo é útil incluem a verificação de uma sequência plasmídica no contexto da engenharia genéticaplasmídica 50,51, a melhoria da precisão das sequências plasmídicas já obtidas pelo WGS24, ou a caracterização de plasmídeos capturados de amostras clínicas ou ambientais por conjugação.
A extração de DNA de amostras clínicas e ambientais pode precisar de alguma personalização. Aqui, os autores relatam que algumas das cepas clínicas que usaram produziram uma matriz de biofilme que entupia os filtros, exigindo múltiplas etapas de pelletagem e, às vezes, a adição de dois filtros. Além disso, para sequenciar plasmídeos conjugativos, é necessário um volume maior de cultura em relação ao sequenciamento de plasmídeos multicópia, pois plasmídeos de médio e grande porte naturalmente tendem a ter um número de cópiasde plasmídeo baixo 2.
O protocolo permitiu o sequenciamento de plasmídeos entre 4 kb e 173 kb de tamanho. Essa faixa está totalmente dentro das especificações do fabricante (entre 2 e 200 kb). Recentemente, um estudo analisou 23.000 sequências de plasmídeos depositadas no NCBI e encontrou diferenças significativas de tamanho dependendo da fonte, com plasmídeos de origem humana tendo o menor tamanho mediano (76 kb), e plasmídeos de solo e plantas os maiores (215 e 427 kb); Plasmídeos de todas as outras fontes tinham tamanho mediano entre 79 e 147 kB, então esse protocolo deve cobrir a maioria dos plasmídeos, sendo possível exceção dos megaplasmídeos52 e plasmídeos de origem soloou vegetal 3.
Os autores incluem uma etapa de enriquecimento de DNA plasmídico para melhorar a precisão das assembleias de DNA plasmídico. A presença de grandes quantidades de DNA cromossômico inunda os poros celulares na célula de fluxo sequenciante com leituras não alvo, reduzindo substancialmente a profundidade do DNA plasmídico. Uma etapa de enriquecimento de DNA plasmídico já foi implementada para tentar a montagem completa do DNA plasmídico em amostras clínicas para a detecção dosARGs 24. Isso mostrou melhorias substanciais na eficácia de completar a montagem de plasmídeos com sequência apenas de nanoporos (78% das montagens completas em relação ao benchmark). Outros exemplos vieramem seguida 50, 53, 54.
Para melhorar a precisão da leitura, dois desses estudos usaram uma técnica desenvolvida pela nanopore chamada "chamada de base pareada em leituras pseudopareadas"50,53. Essa abordagem envolve alinhar os sinais elétricos brutos das cadeias sensoriais e antissenso da mesma molécula de DNA, resultando em uma melhoria na precisão de cerca de uma ordem de magnitude. A chamada de bases pareada em leituras pseudopareadas é compatível com multiplexação, desde que a estratégia permita a identificação das cadeias direta e reversa da mesma molécula. Assim, não pode ser usado em plasmídeos purificados de amostras de complemento plasmídico desconhecido ou contendo uma mistura de plasmídeos. Note que o número de plasmídeos em um genoma normalmente varia entre 0 e 7 plasmídeos (particularmente em isolados clínicos)2,55.
Os pesquisadores notaram um alto nível de contaminação cromossômica nas amostras (Tabela 10). O nível de contaminação de 10 amostras analisadas variou entre 28 e 72% com base no mapeamento de leitura. No entanto, quando os pesquisadores normalizaram essa sequência ao tamanho do cromossomo em relação ao dos plasmídeos presentes, descobriu-se que as sequências plasmídicas eram em média 144 vezes mais abundantes do que as cromossômicas. Traduzir esse número para o nível de enriquecimento de DNA plasmídico exigiria conhecer o número de cópias plasmídicas dos plasmídeos em cada amostra.
Células de fluxo MinION normalmente produzem uma saída de sequenciamento de ~30 GB, possibilitando teoricamente o sequenciamento de bem mais de 96 plasmídeos por execução. Na prática, a multiplexação também é limitada pelo limite de 96 códigos de barras, tornando 96 plasmídeos o limite superior prático por célula de fluxo. No entanto, a contaminação por DNA cromossômico e o desequilíbrio natural do código de barras reduzem a profundidade efetiva de leitura por amostra. Para garantir profundidade de leitura suficiente para cada amostra, uma abordagem mais conservadora é multiplexar aproximadamente 24 plasmídeos por célula de fluxo. A profundidade recomendada para a montagem precisa da sequência plasmídica é entre 50 e70x32,56. Nos exemplos ilustrativos apresentados neste artigo, os autores analisaram 21 amostras multiplexadas em uma célula de fluxo e obtiveram uma profundidade média (109x), que está bem alinhada com o ótimo recomendado para amostras por célula de fluxo.
Note que, com profundidade de leitura excessiva, partições de leitura aprovadas/reprovadas são menos informativas para avaliar a qualidade geral da leitura. Em vez disso, métricas derivadas da distribuição de comprimento fornecem uma representação mais clara da qualidade do conjunto de dados. Uma dessas métricas é o N50, definido como o comprimento de leitura no qual leituras desse comprimento ou mais representam 50% do total de bases. Esse valor é altamente sensível a etapas de controle de qualidade. Por exemplo, filtrar leituras de baixa qualidade ou curtas com ferramentas como Filtlong remove uma fração substancial da cauda mais curta da distribuição, inflando assim o N50. Esse efeito é evidente na Tabela 9.
Os autores usaram o kit Rapid Barcode do ONT, que utiliza um transpósomo para fragmentar o DNA do plasmídeo, e portanto não requer conhecimento prévio do padrão de restrição do plasmídeo. Um estudo anterior constatou que, para plasmídeos pequenos (<20 kb), a fragmentação enzimática é necessária porque pequenos plasmídeos permanecem circularizados durante a maioria das extrações de DNA e, portanto, não contêm extremidades livres para a ligaçãoadaptadora 54. Esse estudo não usou códigos de barras, que geram extremidades livres adicionais por amplificação por PCR, então, neste caso, isso provavelmente não tem importância. O OnRamp, um dos procedimentos avançados para validação rotineira de plasmídeos, evita o uso de códigos de barras ao utilizar leituras de plasmídeos de comprimento total para montagem, simplificando a preparação daamostra 54. Neste caso, os autores preferiram códigos de barras por três razões. A primeira é construir um conjunto independentemente do comprimento do plasmídeo. A segunda é tolerar a contaminação por DNA cromossômico. A terceira é que isso abre a possibilidade de multiplexar várias amostras, como mencionado acima.
Um estudo anterior observou que a agregação de diferentes montadores melhorou o resultado final ao cancelar vieses incorporados nas plataformas de montagemindividuais 24. Essa meta-análise é um recurso incorporado ao programa de montagem que os autores utilizam, o Autocycler34. Em relação à sua versão anterior, chamada Trycycler57, o Autocycler é um pouco mais automatizado; ele utiliza vários subconjuntos das leituras e (como mencionado antes) realiza assemblies com múltiplos montadores para criar uma montagem de consenso mais precisa. No entanto, a maioria dos montadores, mesmo os melhores, como Flye28, Canu29 e Raven30, tem dificuldade com a presença de plasmídeos com homologia de sequências extensa, fundindo erroneamente sequências de plasmídeos ou gerando múltiplas cópias de um único plasmídeo, e tendem a perder pequenosplasmídeos 58,59 (veja também a Figura 4C). Quando isso acontece, cortar leituras de qualidade inferior e remover leituras curtas pode levar à resolução de contigs quiméricos.
Para estabelecer a precisão geral dos protocolos de sequenciamento, os autores compararam os montagens apenas ONT com montagens polidas combinando leituras longas e curtas, que (sem os dados do PacBio) são o padrão ouro atual. Os autores também tinham dados de gel mostrando o número de plasmídeos presentes e seus tamanhos estimados (Tabela 8, terceira coluna). Os pesquisadores montaram com sucesso todos os 13 plasmídeos encontrados nas 10 amostras e obtiveram tamanhos consistentes com as sequências montadas. Isso é um grande sucesso, já que montadores de sequências LR frequentemente falham em montar plasmídeos com sequências estruturalmente complexas. Em um nível mais detalhado, comparações par a par entre os dois conjuntos de sequências indicaram um nível variável de precisão (Tabela 8). Cinco dos conjuntos plasmídicos tinham áreas sem cobertura nos montagens ONT, presumivelmente resultado da sequência adicionada artificialmente durante a montagem de sequências apenas com ONT. Dos outros sete, três apresentaram polimorfismos de nucleotídeo único, com diferenças entre dois e pares variando entre 0,0023 e 0,05%, e quatro apresentaram concordância perfeita. O nível de precisão não parece correlacionar-se com o nível de enriquecimento do DNA plasmídico ou profundidade de sequenciamento, mas a presença de mais de um plasmídeo em uma amostra dada sim, possivelmente porque ter mais de um plasmídeo tende a amplificar o nível de complexidade da sequência (Tabela 8). Em suma, o fluxo de trabalho é projetado para otimizar o rendimento e a precisão das sequências plasmídicas. Embora seja capaz de montar sequências de plasmídeos, produz um nível variável de precisão no nível dos nucleotídeos. O protocolo pode ser usado para sequenciar plasmídeos de uma ampla gama de tamanhos, tolera pelo menos 72% de contaminação por DNA cromossômico e pode acomodar (dependendo do tamanho do plasmídeo e do estado da célula de fluxo) pelo menos 24 plasmídeos em uma única célula de fluxo.