Serviços de Análise de Metilação de DNA: Uma Comparação Prática entre WGBS, EM-Seq, RRBS e Plataformas de Microarray

A metilação do DNA—especificamente a adição de um grupo metilo ao carbono 5 de um anel de citosina (5mC) em dinucleotídeos CpG—é a modificação epigenética mais bem caracterizada no genoma humano. Ela serve como o principal mecanismo para o silenciamento transcricional a longo prazo, a inativação do cromossoma X, a impressão genética e a supressão de elementos transponíveis. Na biopharma moderna, a metilação anormal do DNA não é apenas uma característica do câncer; é um alvo altamente acionável para novas terapias e o biomarcador fundamental para biópsias líquidas de deteção precoce.

No entanto, selecionar a plataforma tecnológica correta para mapear esta modificação é uma decisão notoriamente complexa. O panorama analítico passou por uma mudança tectónica nos últimos cinco anos. Os investigadores já não estão restritos às forças destrutivas do bisulfito de sódio; agora devem navegar por um espectro que vai desde microarrays de alto rendimento até conversões enzimáticas, e até mesmo sequenciação nativa de long-read. Este guia serve como um quadro de decisão consultivo, comparando objetivamente os custos, capacidades e aplicações translacionais de WGBS, EM-Seq, RRBS e Arrays da Illumina.

Para uma visão estratégica mais ampla de como a metilação do DNA se integra com o perfilamento da cromatina em estudos multi-ómicos, recomendamos consultar o nosso Serviços de Análise de Epigenética e Cromatina hub, ou antecipando o nosso próximo guia sobre Integração Multi-Epigenómica.

WGBS: O Padrão Dourado Histórico

A Sequenciação de Bisulfito de Todo o Genoma (WGBS) tem sido há muito o padrão definitivo para o mapeamento do metiloma. Funciona com um princípio químico claro: quando o DNA é tratado com bisulfito de sódio, as citosinas não metiladas são rapidamente desaminadas para uracilo (que são posteriormente lidas como timina (T) durante a PCR e a sequenciação). Em contraste, as citosinas metiladas (5mC) estão quimicamente protegidas da desaminação e permanecem a ser lidas como citosina (C).

Ao mapear estas transições de C para T em relação a um genoma de referência, os bioinformáticos podem determinar o estado de metilação preciso de quase todos os ~28 milhões de locais CpG no genoma humano com resolução de um único nucleótido.

Controlo de Qualidade WGBS: O Spike-In de DNA Lambda

A validade total de um experimento de WGBS baseia-se na suposição de que a conversão por bisulfito foi perfeitamente eficiente. Se a reação química estiver incompleta, as citosinas não metiladas permanecerão como 'C', gerando dados de metilação falsos positivos massivos. Para controlar matematicamente isso, cada reação de WGBS na CD Genomics é adicionada com DNA de fago Lambda não metilado. Como o DNA Lambda não contém 5mC, qualquer 'C' que permanecer nas leituras de Lambda após a sequenciação representa explicitamente a taxa de erro de conversão. Um conjunto de dados de WGBS de qualidade de publicação deve demonstrar estritamente uma eficiência de conversão de Lambda superior a 99,5%.

O Lado Sombrio do Bisulfito: Degradação do DNA

A falha fundamental do WGBS é a dureza da química do bisulfito. O processo requer pH extremo e altas temperaturas, o que destrói fisicamente a espinha dorsal do DNA. Não é incomum que o tratamento com bisulfito degrade mais de 95% do material de DNA inicial. Consequentemente, o WGBS padrão requer quantidades de entrada massivas (frequentemente de 100 ng a 1 μg de DNA genómico puro), tornando-o altamente proibitivo para amostras de tecido preciosas pré-translacionais, células-tronco classificadas por fluxo ou DNA livre de células circulantes (cfDNA).

O ponto cego do 5mC vs. 5hmC: Por que o bisulfito não é perfeito

Para além da degradação do DNA, o sequenciamento padrão com bisulfito apresenta uma profunda lacuna biológica: não consegue diferenciar entre 5-metilcitosina (5mC, uma marca repressiva) e 5-hidroximetilcitosina (5hmC, uma marca intermediária altamente ativa no caminho de desmetilação ativa impulsionado pelas enzimas TET). A química do bisulfito protege TANTO a 5mC como a 5hmC, o que significa que uma leitura padrão de WGBS confunde estes dois sinais epigenéticos opostos num único ponto de dados "metilado".

Na neurobiologia (onde 5hmC está altamente enriquecido no cérebro) e na oncologia (onde as mutações TET são comuns), distinguir estas marcas é crítico. Para resolver isto, os investigadores devem implementar oxBS-Seq (Sequenciação de Bisulfito Oxidativa)Este fluxo de trabalho paralelo utiliza perrutenato de potássio para oxidar 5hmC em 5-formilcitosina (5fC)—que é susceptível à desaminação por bisulfito. Ao subtrair matematicamente os dados de oxBS-Seq dos dados padrão de BS-Seq, os bioinformáticos podem gerar um mapa absoluto e diferenciado de verdadeira 5mC versus verdadeira 5hmC. Embora seja poderoso, isso requer a execução de cada amostra duas vezes, dobrando efetivamente o custo de sequenciação.

Figura 1: O Mecanismo Químico do Bisulfito vs. Bisulfito Oxidativo (oxBS-Seq)

EM-Seq: A Revolução Enzimática para Baixa Entrada e cfDNA

Para contornar a degradação catastrófica do DNA causada pelo bisulfito, a indústria tem mudado agressivamente para Methyl-seq Enzimático (EM-Seq)Em vez de produtos químicos agressivos, o EM-Seq baseia-se em uma cascata enzimática suave de dois passos:

1. ProteçãoA enzima TET2 e um potenciador de oxidação são utilizados para oxidar rapidamente 5mC e 5hmC em estados intermédios protegidos.
2. DesaminaçãoA enzima APOBEC é utilizada para desaminar seletivamente as citosinas não metiladas restantes e desprotegidas em uracilo.

A leitura final é funcionalmente idêntica à WGBS (Cs não metilados tornam-se Ts), mas como a reação é inteiramente enzimática e ocorre a pH neutro, a integridade estrutural do DNA é perfeitamente preservada. Isso resulta em tamanhos de inserção de biblioteca significativamente maiores, eficiências de mapeamento mais altas e uma cobertura genómica muito mais uniforme (eliminando o severo viés de GC frequentemente observado na WGBS).

A Ferramenta Definitiva para Biópsia Líquida

Porque o EM-Seq previne a degradação do ADN, tornou-se o de facto padrão para analisar DNA livre de células circulantes (cfDNA) em biópsias líquidas. O cfDNA já está altamente fragmentado (tipicamente cerca de 166 pb) e existe em quantidades mínimas (frequentemente <10 ng por extração de plasma). Submeter o cfDNA a WGBS destrói-o completamente. O EM-Seq permite que os investigadores oncológicos mapeiem todo o metiloma de um tumor de forma não invasiva, abrindo caminho para painéis de deteção precoce de câncer de próxima geração.

Além disso, uma vez que o EM-Seq preserva fragmentos de cromatina intactos, é especialmente adequado para integrações multi-ómicas. Uma investigação recente que executámos com sucesso envolveu a combinação EM-Seq com ATAC-Seq na mesma subpopulação rara de tumores, permitindo ao cliente correlacionar simultaneamente a acessibilidade da cromatina com o estado de metilação em locais específicos de potenciadores, sem dividir a preciosa amostra.

RRBS: Enriquecimento com Economia

Embora o WGBS e o EM-Seq forneçam paisagens genómicas completas, são indiscutivelmente caros, muitas vezes exigindo uma profundidade de 30x em 3 mil milhões de pares de bases. No entanto, a metilação do DNA não está distribuída de forma uniforme. A vasta maioria das alterações de metilação biologicamente relevantes (especialmente no câncer) ocorre em ilhas CpG (CGIs) — aglomerados densos de dinucleotídeos CpG localizados principalmente em promotores de genes.

Sequenciação de Bisulfito com Representação Reduzida (RRBS) é um método de enriquecimento altamente estratégico projetado para sequenciar apenas estas regiões informativas. O DNA genómico é digerido utilizando uma enzima de restrição insensível à metilação (tipicamente MspI, que corta em motivos CCGG). Como o MspI corta com mais frequência em regiões ricas em CpG, a seleção de tamanho dos fragmentos resultantes (geralmente entre 40 e 220 bp) enriquece fortemente as ilhas de CpG e as regiões promotoras.

Ao concentrar o poder de sequenciação exclusivamente nos ~1-2% do genoma que realmente importa para a regulação transcricional, o RRBS captura mais de 85% de todas as ilhas CpG a aproximadamente 5% a 10% do custo do WGBSPara investigadores que realizam triagens na fase de descoberta em dezenas de modelos de ratos knockout ou linhas celulares tratadas com compostos, o RRBS oferece um equilíbrio imbatível entre eficiência de custos e relevância biológica.

Figura 2: Perfil de Cobertura Genómica: WGBS vs. RRBS

Arrays de Metilação: O Peso Pesado da Genómica Populacional

Para coortes de investigação translacional em grande escala e Estudos de Associação Epigenómica (EWAS), até mesmo RRBS pode tornar-se proibitivamente caro e bioinformaticamente complicado ao escalar para milhares de sujeitos. Nestes contextos epidemiológicos, Microarranjos de Metilação (especificamente as plataformas Illumina Infinium) continuam a ser as melhores.

Plataformas Illumina EPIC 850K e 935K

A evolução dos arrays da Illumina — desde os arrays legados de 27K e 450K até o atual 850K (EPIC) e o altamente avançado Illumina EPIC v2 (935K)—providenciou aos investigadores uma ferramenta de triagem de conteúdo fixo incrivelmente robusta. A matriz EPIC v2 investiga mais de 935.000 locais de CpG altamente curados em promotores, potenciadores e elementos regulatórios identificados pelo ENCODE.

O Relógio Epigenético e EWAS

As matrizes operam hibridizando DNA convertido em bisulfito a sondas ligadas a esferas, produzindo um "valor Beta" contínuo (de 0 a 1) que representa a percentagem absoluta de metilação naquele local específico de CpG. Como os locais-alvo são fixos, as matrizes oferecem consistência absoluta em grandes coortes, evitando completamente as quedas de cobertura estocásticas comuns em métodos baseados em NGS.

Esta imensa consistência torna os arrays a plataforma obrigatória para desenvolver "Relógios Epigenéticos" (como o Relógio de Horvath)—algoritmos preditivos que utilizam o estado de metilação de algumas centenas de locais CpG altamente conservados para calcular a idade biológica do tecido de um sujeito. Se estiver a validar um biomarcador translacional em 5.000 amostras de sangue de sujeitos, um "Serviço de Análise de Metilação de DNA (Illumina 935K)" é a única escolha matematicamente e economicamente sólida.

O Gargalo da Bioinformática: Alinhamentos Assimétricos

Um dos desafios computacionalmente mais intensivos na epigenómica é o alinhamento das leituras de sequenciamento convertidas por bisulfito. No sequenciamento de DNA padrão (como RNA-Seq ou ChIP-Seq), o genoma de referência é estático. No entanto, no WGBS ou RRBS, a química do bisulfito quebra a simetria fundamental de emparelhamento de bases de Watson-Crick. Um 'C' não metilado na cadeia direta é convertido em 'T', o que significa que não irá mais complementar perfeitamente o 'G' na cadeia inversa.

Para superar isto, os bioinformáticos devem utilizar alinhadores especializados, sendo o mais famoso o Bismark. O Bismark realiza matematicamente um alinhamento de três letras. Ele converte computacionalmente cada 'C' no genoma de referência humano em 'T' e, simultaneamente, converte cada 'G' na fita de referência reversa em 'A'. Em seguida, tenta alinhar as leituras de sequenciamento brutas contra estes genomas teoricamente convertidos e massivamente expandidos. Este processo de alinhamento paralelo de 4 fitas requer uma enorme RAM computacional e poder de processamento, resultando frequentemente em gargalos na análise de dados que apanham grupos de investigação despreparados de surpresa.

Nuances Técnicas de Microarrays: Probes Tipo I vs. Probes Tipo II

Embora os arrays da Illumina (como o 935K EPIC v2) sejam altamente padronizados, a análise dos dados brutos (geralmente fornecidos como arquivos IDAT) requer uma profunda expertise técnica para evitar viés sistémico. Os arrays utilizam dois designs de sondas químicas fundamentalmente diferentes para capturar os 935.000 locais CpG:

• Sondas de Tipo IEstes utilizam duas contas físicas distintas por local CpG—uma conta projetada explicitamente para ligar a sequência não metilada (convertida) e uma segunda conta projetada para ligar a sequência metilada (não convertida). Elas são lidas em um único canal de cor.
• Sondas de Tipo IIEstes utilizam apenas uma única pérola física por local CpG. A extremidade da sonda termina exatamente uma base antes do CpG alvo, e o sequenciador lê se um 'G' rotulado (metilado) ou um 'A' (não metilado) é incorporado, utilizando dois canais de cor diferentes (Vermelho/Verde).

Porque as sondas do Tipo I e do Tipo II têm cinéticas de ligação e intervalos dinâmicos de fluorescência diferentes, produzem distribuições de valores Beta matematicamente distintas. Se um bioinformático simplesmente calcular a média dos dados brutos, os resultados estarão altamente distorcidos. Algoritmos de normalização sofisticados, como o BMIQ (Beta Mixture Quantile dilation), devem ser aplicados durante o pré-processamento para alinhar perfeitamente a distribuição estatística das sondas do Tipo II com as sondas do Tipo I, garantindo uma consistência matemática absoluta em todo o conjunto de dados de 935K antes de se tirarem quaisquer conclusões sobre biomarcadores.

A Superar os Limites: Metilação de DNA a Nível de Célula Única (scWGBS)

A WGBS tradicional, tal como o RNA-Seq em massa, faz uma média do sinal de metilação em milhões de células. Embora isso seja suficiente para linhagens celulares homogéneas, falha completamente em capturar a diversidade epigenética dentro de tecidos altamente heterogéneos, como o cérebro mamífero ou um microambiente tumoral sólido.

A fronteira final da perfuração com bisulfito é o Sequenciamento de Bisulfito de Genoma Inteiro em Célula Única (scWGBS) e a sua iteração avançada específica para neurónios, snmC-seq (Sequenciamento de Metilcitosina de Núcleo Único). Adaptar o WGBS a uma única célula é cientificamente desafiador porque a degradação por bisulfito tipicamente destrói mais de 90% do DNA. Se começar com o genoma de uma única célula (aproximadamente 6 picogramas de DNA) e perder 90% dele, a biblioteca resultante está efetivamente vazia.

Para resgatar isto, a CD Genomics emprega fluxos de trabalho de célula única altamente otimizados e de baixa perda. Os núcleos são separados via FACS diretamente em placas de 384 poços contendo tampão de lise. Após uma conversão de bisulfito especializada e mais suave, é utilizada uma amplificação linear com primários aleatórios em vez da PCR padrão. Esta estratégia de amplificação única contorna a necessidade de ligação imediata de adaptadores em DNA altamente fragmentado, resgatando os fragmentos de DNA restantes e permitindo que os bioinformáticos construam um mapa de metilação de alta resolução de um único neurônio ou de uma única célula tumoral circulante.

Multi-Ómicas: A Convergência do EM-Seq e ATAC-Seq

À medida que o campo da epigenómica amadurece, os investigadores estão a exigir cada vez mais insights em múltiplas camadas a partir da mesma amostra. Um pedido frequente dos clientes da biopharma é o desejo de correlacionar a metilação do DNA (que tipicamente silencia a cromatina) com a acessibilidade física da cromatina (mapeada através de ATAC-Seq).

Historicamente, isso exigia dividir a biópsia ao meio—analisando metade para WGBS e metade para ATAC-Seq. No entanto, com o advento do EM-Seq, a verdadeira multi-ômica está a tornar-se uma realidade. Como o EM-Seq utiliza enzimas não destrutivas (TET2/APOBEC) em vez de bisulfito agressivo, o DNA permanece completamente intacto. Protocolos de ponta agora permitem que os investigadores primeiro tagmentem a cromatina aberta usando a transposase Tn5 (a etapa ATAC-Seq) e, em seguida, submetam imediatamente essa mesma biblioteca à conversão enzimática do EM-Seq. Isso resulta numa biblioteca de sequenciamento unificada onde os investigadores podem simultaneamente ler onde a cromatina estava fisicamente aberta E se as ilhas CpG circundantes estavam metiladas, tudo a partir de uma única e preciosa amostra de tecido pré-translacional.

Sequenciação de Metilação Direcionada: O Ponto Ideal para Painéis de Investigação Translacional

Embora o WGBS ofereça uma descoberta abrangente e o RRBS uma triagem económica, nenhum deles é, na prática, adequado para procedimentos de investigação de biomarcadores de rotina (Nota: Os nossos serviços são estritamente RUO) ou monitorização de biomarcadores em estágios avançados de tradução. Num ambiente de investigação translacional, os oncologistas não precisam sequenciar os 28 milhões de locais CpG no genoma humano; eles apenas precisam sequenciar um painel altamente específico e pré-definido de 50 a 500 locais CpG que se sabe estarem definitivamente alterados em um tipo específico de câncer (por exemplo, a metilação do SEPT9 no câncer colorretal).

Para estas aplicações, Sequenciação de Metilação Direcionada representa o ponto ideal. Esta tecnologia utiliza isótopos de RNA biotinilados, projetados sob medida, para "capturar" fisicamente e reter apenas as regiões genómicas de interesse antes da sequenciação. Este enriquecimento direcionado permite que os investigadores concentrem todo o seu orçamento de sequenciação numa fração minúscula do genoma.

A principal vantagem da Sequenciação de Metilação Direcionada é cobertura ultra-profundaEnquanto uma corrida padrão de WGBS pode gerar uma cobertura de 30x em todo o genoma, um painel de captura direcionada pode facilmente alcançar uma cobertura de 1.000x a 5.000x nas regiões selecionadas. Esta profundidade extrema é absolutamente crítica para biópsias líquidas e monitorização de Doença Residual Mínima (MRD). Ao tentar detectar DNA tumoral circulante (ctDNA) numa amostra de sangue onde 99,9% do cfDNA é derivado de leucócitos saudáveis, uma profundidade de sequenciamento de 30x provavelmente não conseguirá captar o raro fragmento tumoral metilado. Apenas sequenciando o locus 5.000 vezes é que os bioinformáticos podem garantir estatisticamente a deteção de um evento de metilação rara com alta confiança estatística.

Além disso, a captura direcionada integra-se perfeitamente com a química EM-Seq mencionada anteriormente. Ao combinar o enriquecimento de isótopos de RNA personalizados com a conversão enzimática EM-Seq não destrutiva, os desenvolvedores de biopharma podem criar a plataforma de biópsia líquida definitiva: capaz de capturar cfDNA fortemente fragmentado, preservando-o sem degradação química e sequenciando os loci associados ao câncer específicos a profundidades sem precedentes.

A Fronteira de 2026: Sequenciação Direta de Metilação de Terceira Geração

Enquanto os arrays e a sequenciação EM dominam o panorama atual, o futuro definitivo da epigenómica reside na Sequenciação de Longas Leituras de Terceira Geração (TGS), especificamente em plataformas como Oxford Nanopore Technologies (ONT) e PacBio (sequenciação SMRT).

Estas plataformas alcançam algo anteriormente considerado impossível: Deteção direta e nativa de metilação sem qualquer conversão química ou enzimática.

• Sequenciação SMRT da PacBioÀ medida que a polimerase sintetiza a cadeia complementar, apresenta atrasos cinéticos microscópicos (medidos em milissegundos) ao incorporar uma base oposta a uma citosina metilada em comparação com uma citosina não metilada. O software analisa estas cinéticas de pulso para mapear diretamente as modificações 5mC e 5hmC em tempo real.
• Sequenciação por NanoporosÀ medida que uma fita de ADN nativa e não amplificada é puxada através do nanoporo de proteína, o volumoso grupo metilo numa base 5mC provoca uma interrupção de corrente elétrica distintamente diferente da de uma citosina padrão. Algoritmos avançados de aprendizagem automática (como Megalodon ou Remora) decifram esta ondulação elétrica numa chamada de metilação de alta confiança.

Crucialmente, porque estas tecnologias leem fragmentos que têm entre 10kb a 100kb de comprimento, permitem que os investigadores façam a fase de padrões de metilação de forma contínua através de variantes estruturais massivas, regiões centroméricas repetitivas e domínios imprintados específicos de alelos complexos — regiões onde o WGBS de leituras curtas falha completamente. Embora atualmente seja mais caro por gigabase do que o sequenciamento Illumina, a metilação direta por Nanopore está rapidamente a tornar-se a plataforma de escolha para montagens complexas do genoma humano.

O Quadro de Decisão Definitivo

Escolher a plataforma certa é um cálculo estratégico que equilibra orçamento, amostra de entrada e âmbito biológico. A CD Genomics utiliza uma abordagem independente da plataforma, adaptando a tecnologia estritamente à sua questão biológica.

Perguntas Frequentes

Referências:

1. Vaisvila, R., et al. (2021). "A sequenciação metílica enzimática deteta a metilação do DNA com resolução de base única a partir de picogramas de DNA." Pesquisa Genómica, 31(7), 1280-1289.
2. Gu, H., et al. (2011). "Preparação de bibliotecas de sequenciação de bisulfito de representação reduzida para perfilagem de metilação de DNA em escala genómica." Protocolos da Natureza, 6(4), 468-481.
3. Horvath, S. (2013). "Idade de metilação do DNA de tecidos e tipos celulares humanos." Biologia Genómica, 14(10), R115.
4. Simpson, J. T., et al. (2017). "Deteção da metilação de citosina no DNA utilizando sequenciação por nanoporo." Métodos da Natureza, 14(4), 407-410.

Serviços Relacionados

• Sequenciação de Bisulfito de Genoma Inteiro (WGBS)
• Sequenciação de Bisulfito com Representação Reduzida (RRBS)
• Sequenciação de Bisulfito Direcionada
• Microarray de Metilação de DNA
• Análise de Dados Epigenómicos

* Apenas para uso em investigação. Não é para uso em procedimentos de diagnóstico.