O sequenciamento de fragmentos de DNA associado a sítios de restrição (do inglês Restriction site associated DNA sequencing, RADseq) é uma metodologia que descreve o sequenciamento dos marcadores RAD, que são um tipo de marcador molecular baseado em sítios de restrição enzimática ao longo do genoma. Essa metodologia consiste na digestão do DNA genômico por uma ou mais enzimas de restrição e ligação dos fragmentos a adaptadores, o que permite a multiplexação de centenas a milhares de indivíduos, permitindo uma representação reduzida do genoma, dos loci próximos aos sítios de restrição. Dessa forma caracteriza-se por um método de sequenciamento de bibliotecas de representação reduzida (RRLs, do inglês reduced-representation libraries)A técnica consiste basicamente em fragmentar seletivamente o DNA através de uma enzima de restrição, seguida pelo sequenciamento de alto rendimento de uma porção específica dos fragmentos. Os pontos fortes do RAD-Seq incluem sua versatilidade, economia e capacidade de gerar grandes quantidades de dados das mesmas regiões de um genoma em uma grande população de amostras rapidamente e sem a necessidade de um genoma de referência. Andrews et al., 2016 A principal novidade de metodologias como o RADseq em relação ao sequenciamento de DNA de nova geração (NGS) na plataforma Illumina é o método de fragmentação do DNA. No sequenciamento NGS tradicional, o DNA é fragmentado mecanicamente, enquanto que metodologias de representação reduzida do genoma utilizam enzimas de restrição para realizar essa fragmentação. Na metodologia RADseq original, apenas uma enzima de restrição é utilizada, gerando fragmentos com uma extremidade cortada pela enzima e portanto contendo seu sítio de restrição, e a outra extremidade com sequência aleatória, pois é cortada mecanicamente. Adaptadores que contém a sequência do primer de sequenciamento da Illumina seguida por uma sequência de identificação chamada barcode (em português, código de barras) são então ligados às extremidades dos fragmentos. Esses barcodes permitem a identificação de cada amostra, de modo que é possível multiplicar várias amostras em uma mesma reação de sequenciamento. Os fragmentos de DNA gerados são então selecionados por tamanho, reduzindo então a quantidade do genoma a ser sequenciado. Essa preparação do DNA para ser sequenciado é chamada de preparo de biblioteca (ou construção de biblioteca). .

Metodologia Embora cada método baseado em enzimas de restrição tenha suas particularidades no processo de preparo da biblioteca para o sequenciamento, eles compartilham um número de etapas comuns: Extração, quantificação e qualidade do DNA genômico; fragmentação do DNA genômico de todas as amostras com enzimas de restrição e ligação dos adaptadores; amplificação por PCR (RAD-seq e GBS) e seleção de tamanho dos fragmentos (RRL e RAD-seq); sequenciamento e análise das sequências com ou sem suporte de um genoma de referência, e identificação dos SNPs (Davey et al., 2011b; Poland and Rife, 2012). A principal novidade de metodologias como o RADseq em relação ao sequenciamento de DNA de nova geração (NGS) na plataforma Illumina é o método de fragmentação do DNA. No sequenciamento NGS tradicional, o DNA é fragmentado mecanicamente, enquanto que metodologias de representação reduzida do genoma utilizam enzimas de restrição para realizar essa fragmentação. Na metodologia RADseq original, apenas uma enzima de restrição é utilizada, gerando fragmentos com uma extremidade cortada pela enzima e portanto contendo seu sítio de restrição, e a outra extremidade com sequência aleatória, pois é cortada mecanicamente. Adaptadores que contém a sequência do primer de sequenciamento da Illumina seguida por uma sequência de identificação chamada barcode (em português, código de barras) são então ligados às extremidades dos fragmentos. Esses barcodes permitem a identificação de cada amostra, de modo que é possível multiplexar várias amostras em uma mesma reação de sequenciamento. Os fragmentos de DNA gerados são então selecionados por tamanho, reduzindo então a quantidade do genoma a ser sequenciado. Essa preparação do DNA para ser sequenciado é chamada de preparo de biblioteca (ou construção de biblioteca).

Aplicações Métodos de representação reduzida do genoma como o RADseq são muito úteis em estudos na área de genômica evolutiva, ecológica e conservação de espécies. Estudos populacionais necessitam de grande número amostral, sendo necessário, portanto, alternativas mais rápidas e baratas ao sequenciamento do genoma inteiro, mas que sejam representativas. Trabalhos que utilizaram RADseq foram capazes de identificar dezenas de milhares de polimorfismos de nucleotídeo único (SNPs) no genoma de suas espécies de interesse. O mapeamento de SNPs ao longo de genoma permite estimar tamanho efetivo de populações, responder perguntas referentes à adaptação, introgressão, estruturação populacional, e estudos filogenéticos. A metodologia RADseq é voltada especialmente para estudos na área de genômica evolutiva mas também pode ser aplicado no campo da saúde humana, como no estudo da evolução de tumores.O RADseq irá ter uma enorme aplicação em análise genética sobre pesquisas do complexo genômico. Entrando mais no conceito da aplicação dessa técnica, como quando se trata de quadros em que envolvam problemas de saúde humana igual falado sobre a evolução de tumores, uma pesquisa sobre “Diversidade dos Tumores e evolução revelada através do RADseq” mostrou que quando o mesmo foi aplicado tanto em amostras de câncer pancreático humano quanto em melanomas de peixe-zebra e utilizando enzimas de restrição de baixa frequência ou alta frequência, foi identificado com sucesso substituições de nucleotídeos únicos e alterações no número de cópias nos tumores. Com isso, foi mostrada que tais descobertas podem ser desenvolvidas ao efetuar o RADseq em sub linhagens dentro do tumor. Foi também inferido relações filogenéticas entre tumores primários e metástases. Os estudos evolutivos de câncer que se aplicam nas taxas de evolução tumoral e nas relações evolutivas entre linhagens tumorais irão se beneficiar da maleabilidade e eficiência da sequenciação de DNA associada ao local de restrição. Além disso, a aplicação é muito vista em conservação de espécies principalmente se tratando de plantas, igual foi abordada nesse estudo sobre aplicação do RAD-seq para a avaliação da diversidade genética do Panax notoginseng (Araliaceae) domesticado, que se trata de uma planta medicinal chinesa. O RAD-seq serviu para analisar a variação genética e diferenciação desta. Os resultados de diversidade de nucleotídeos e heterozigosidade indicaram que a P. notoginseng tem baixa diversidade genética a nível de espécie e de população e quase nenhuma diferenciação genética foi localizada, e todas as populações eram semelhantes geneticamente devido a um forte fluxo gênico e tempo insuficiente para divisão. Outra aplicação do método visto em plantas foi a “RADseq approaches and applications for forest tree genetics” que detalhou como dados foram gerados a partir de diversas espécies de árvores e que também traz como as publicações sobre o uso do RAD-seq nesse tema aumentou significativamente, sugerindo que ainda tende a crescer cada vez mais e irá aproximar de gerar dados SNP. Também abrangendo áreas onde o método está contribuindo para análise genética de árvores, como mapeamento linkage, estrutura genética e história de populações, genética quantitativa molecular, inferência filogenética e base genética de adaptação e fenótipo.

Seleção de Espécies Métodos de representação reduzida do genoma como o RADseq são muito úteis em estudos na área de genômica evolutiva, ecológica e conservação de espécies. Estudos populacionais necessitam de grande número amostral, sendo necessário, portanto, alternativas mais rápidas e baratas ao sequenciamento do genoma inteiro, mas que sejam representativas. Trabalhos que utilizaram RADseq foram capazes de identificar dezenas de milhares de polimorfismos de nucleotídeo único (SNPs) no genoma de suas espécies de interesse. (informações anteriores). Para uma melhor observação das polimorfias, é necessário um número amostral suficientemente grande para mapeamento genético dos locus ou um mapeamento específico em locus conhecidos sobre o estudo. Como no estudo do esgana-gata de três espinhos, que seu genoma foi sequenciado pois apresenta variação genética quantitativa, então foi realizado mapeamento em locais marcadores potenciais acontecendo como cruz controlada sobre a quantidade, onde as limitações são na diversidade permanente no país de pesquisa. Vale ressaltar que para o uso de RADseq em populações selvagens é diferente, pois apresentam uma distribuição de indivíduos que exclui grandes espécies da consideração para estudo. Também há dificuldade em mapear a estrutura gênica em organismos móveis e em espécies empobrecidas que sofreram alterações demográficas. Nesse caso, a precisão de análises genéticas populacionais depende do número de loci utilizados a partir da análise da variação de disposição de organismos móveis.

Histórico A técnica de marcadores do tipo RADseq (Restriction-site associated DNA) são fragmentos pequenos de DNA adjacentes a sítios de reconhecimento de enzimas de restrição. Inicialmente estes marcadores foram desenvolvidos com o propósito de detectar polimorfismos de sítios únicos (SNPs) por meio de microarranjos (microarrays) (MILLER et al., 2007). Em 2000 iniciaram-se os primeiros conceitos fundamentais sobre Radseq, que envolviam a marcação de locais de restrição no genoma e a subsequente amplificação e sequenciamento dessas regiões. O RADseq foi descrito em 2008 por Nathan Baird e colegas, o mesmo grupo descreveu os marcadores RAD no ano anterior, na Universidade de Oregon, Estados Unidos. Desde então foi amplamente utilizado, principalmente em trabalhos na área de ecologia molecular e genômica evolutiva, e novas abordagens foram criadas, se tornando hoje uma família de métodos. Em 2012, por exemplo, Brant Peterson descreveu uma metodologia variante chamada “double-digest RADseq” (ddRADseq), que inclui uma segunda enzima de restrição, diferentemente do RADseq original que utiliza apenas uma. Dessa forma as duas extremidades de cada fragmento de DNA são conhecidas, o que permite maior redução da biblioteca, e tamanhos de fragmentos mais específicos, reduzindo custo. A partir desta variante e do RADseq original outras novas variantes também surgiram, dando origem a dezenas de metodologias derivadas, como 3RAD, 2bRAD, hyRAD, quaddRAD, entre outras.

Em princípio os marcadores RAD eram utilizados através de SNPs arrays, posteriormente foram adaptadores para o sequenciamento NGS.

Técnicas Laboratoriais

Etapas Laboratoriais envolvidas na técnica RADseq:

Primeiramente o método inclui cinco principais etapas dentre elas estão: Corte do DNA pela restrição de enzima; dessa maneira os fragmentos de cada amostra são ligados a um adaptador P1 que contém uma extremidade adesiva que torna possível a ligação e um MID (identificador molecular que irá identificar e marcar exclusivamente o fragmento); esses fragmentos são agrupados e cortados de forma aleatória para gerar fragmentos de tamanhos entre 300-700pb, ao final todos os fragmentos serão ligados a um adaptador P2 com final divergente em todos os fragmentos com e sem adaptadores P1. Esses fragmentos serão amplificados por PCR com oligonucleotídeos específicos para P1 e P2, isso significa que apenas fragmentos com adaptadores P1 e P2 são amplificados, ou seja, os fragmentos que contém os sítios de restrição.

Trabalhos que utilizaram RADseq

RADSeq: next-generation population genetics <https://academic.oup.com/bfg/article/9/5-6/416/182576>

John W Davey, Mark L Blaxter