Sequenciamento de captura de DNA de ligação de metilo em tecidos de pacientes

Resumo

Aqui apresentamos um protocolo para investigar a metilação do DNA do genoma ampla em grande escala estudos de triagem de pacientes clínicos utilizando o Metil-Binding tecnologia de sequenciamento de DNA Captura (MBDCap-seq ou MBD-seq) eo gasoduto análise bioinformática subsequentes.

Resumo

A metilação é uma das modificações epigenética essenciais para o DNA, o qual é responsável pela regulação precisa dos genes necessários para o desenvolvimento e diferenciação estável de diferentes tipos de tecidos. A desregulação deste processo é frequentemente a característica de diversas doenças como o cancro. Aqui, nós descrevemos uma das técnicas de sequenciação recentes, Metil-Binding sequenciação de ADN de captura (MBDCap-SEQ), usada para quantificar a metilação em vários tecidos normais e de doença para grandes grupos de pacientes. Descreve-se um protocolo detalhado desta abordagem de enriquecimento por afinidade ao longo de uma conduta bioinformática para alcançar quantificação ideal. Esta técnica tem sido utilizada para sequenciar centenas de pacientes em vários tipos de câncer, como parte do projeto methylome 1.000 (Sistema Methylome Câncer).

Introdução

Regulação epigenética de genes por meio de metilação do ADN é um dos mecanismos essenciais necessários para determinar o destino da célula por diferenciação estável de diferentes tipos de tecidos do corpo ^1. A desregulação deste processo tem sido conhecido por causar diversas doenças, incluindo o cancro ^2.

Este processo envolve, principalmente, a adição de grupos metilo no resíduo de citosina nos dinucleótidos CpG de ADN ^3. Existem algumas técnicas diferentes, actualmente utilizados para investigar este mecanismo, tendo cada um as suas próprias vantagens, conforme descrito em muitos estudos ^2-8. Aqui vamos discutir uma destas técnicas chamados metil-Binding Captura de sequenciamento de DNA (MBDCap-seq), onde usamos uma técnica de enriquecimento de afinidade para identificar regiões metiladas do DNA. Esta técnica baseia-se na capacidade de ligação de metilo da proteína MBD2 para enriquecer os fragmentos de DNA genómico contendo locais CpG metilados. Nós utilizamos um kit de enriquecimento de DNA metilado comercialpara o isolamento destas regiões metiladas. O nosso laboratório fez a triagem centenas de amostras de pacientes, utilizando esta técnica e aqui nós fornecemos um protocolo optimizado abrangente, que pode ser usado para investigar grandes grupos de pacientes.

Como é evidente, com qualquer tecnologia de sequenciamento de próxima geração, MBDCap-seq também requer uma abordagem bioinformática específicos, a fim de quantificar com precisão os níveis de metilação através das amostras. Tem havido muitos estudos recentes em um esforço para otimizar o processo de normalização e análise dos dados de sequenciamento ^{9, 10.} Neste protocolo, demonstramos um destes métodos de aplicação de uma abordagem única de recuperação de leitura - LONUT - seguido pela normalização linear de cada amostra, a fim de permitir comparações imparciais toda grande número de amostras de doentes.

Protocolo

Todos os tecidos são obtidos após a aprovação da comissão Institutional Review Board e quando todos os participantes consentiram em ambas as análises moleculares e estudos de acompanhamento. Os protocolos são aprovados pela Comissão de Estudos Humanos da Universidade do Texas Health Science Center em San Antonio.

1. Captação de ligação ao ADN de metilo (MBDCap)

1. Coleta de amostras e isolamento de ADN
   1. Recolha tumor granel ou amostras de tecidos normais a partir de amostras de tecido embebidos em parafina de pacientes.
   2. Usar um DNA Mini Kit comercial para isolar o ADN genómico a partir das amostras de tecido embebidos em parafina de acordo com o protocolo do fabricante.
   3. Use um kit de enriquecimento de DNA metilado com o processo aqui descrito de acordo com o protocolo do fabricante.
2. Lavagem talão inicial
   Nota: Os grânulos são utilizados para acoplar com a proteína MBD-Biotina, que detecta a metilação de ADN no genoma. Grânulos são usualmente suspensoem uma solução de estoque. Use estas etapas para lavar este líquido.
   1. Ressuspender o estoque de contas de estreptavidina (Veja a Tabela de Materiais) cuidado pipetando para cima e para baixo para obter uma suspensão homogénea. Não misture as contas em vórtice.
   2. Para um ug de ADN de entrada, ADD10 uL de esferas para um tubo com 90 ul de 1x Bind Buffer / lavagem. Gire durante 1 min. Não misturar em vórtice.
   3. Colocar o tubo (s) em uma cremalheira magnético durante 1 minuto.
   4. Remover o líquido com uma pipeta e eliminar o líquido.
   5. Adicionar 250 ul de 1x tampão de Bind / lavagem para os grânulos e girar durante 1 min.
   6. Repita os passos 1.2.3 - 1.2.5 duas vezes.
3. O acoplamento da proteína MBD-Biotina para os grânulos
   1. Para 1,2 ug de ADN de entrada, adicionar 7 ul (3,5 ug) de proteína MBD-Biotina a cada tubo.
   2. Adicionar 93 ul de 1x Bind Buffer / lavagem à proteína para fazer um volume final de 100 uL.
   3. Misturar a proteína de grânulos-mistura em um misturador rotativo à temperatura ambiente durante 1 h.
4. DNA fragmentado (entrada)
   1. Sonicar ADN pela adição de 1,2 ug de ADN genómico total em 96 uL de tampão TE e 24? L de 5x tampão de Bind / lavagem para fazer uma solução de 120 uL. Use 30 s ON poder e 30 s OFF energia durante sonicação, 20 - 25 vezes.
   2. Corra 1 uL de solução de ADN sonicada usando um sistema Bioanalyzer com um kit de ADN comercial de alta sensibilidade para verificar o tamanho dos fragmentos de estar na gama de 150-350 pb de acordo com o protocolo do fabricante.
5. Lave as MBD-esferas
   1. Colocar o tubo contendo MBD-esferas em uma cremalheira magnético durante 1 minuto.
   2. Retirar e descartar o líquido com uma pipeta sem tocar as esferas.
   3. Ressuspender as esferas com 250 ul de 1x tampão de Bind / lavagem.
   4. Misturar os grânulos num misturador rotativo à temperatura ambiente durante 5 min.
   5. Repita os passos 1.5.1 - 1.5.4 duas vezes.
6. Reação de captura de DNA fragmentado (1 ug de DNA de entrada)
   1. Transferir a mistura de ADN / tampão ao tubo contendo os MBD-grânulos.
   2. Misturar os grânulos com MBD-ADN em um misturador rotativo durante 1 h à TA. Como alternativa, misture O / N a 4 ^o C.
7. A remoção de DNA não-capturado a partir das pérolas
   1. Depois de misturar o ADN e MBD-grânulos, colocar o tubo na cremalheira magnética durante 1 min a concentrar-se todas as esferas na parede interior do tubo.
   2. Remover o líquido sobrenadante com uma pipeta e salve-o em um tubo de microcentrífuga livre de DNase limpo como DNA não metilado. Armazenar esta amostra em gelo.
   3. Adicionar 200 ul de 1x tampão de Bind / lavagem para os grânulos para lavar as pérolas.
   4. Misturar os grânulos num misturador rotativo durante 3 min.
   5. Colocar o tubo na cremalheira magnético durante 1 minuto.
   6. Retirar o líquido com uma pipeta e salve-o em um tubo de microcentrífuga.
   7. para capture reações de ≤1 mg de DNA de entrada, repita os passos 1.7.3 - 1.7.6 para mais 2 frações de lavagem no total.
8. Eluição em dose única
   1. Misture solução salina diluída com tampão de alto teor salino (1: 1) para obter tampão de eluição (NaCl 1,000 mM).
   2. Ressuspender as pérolas em 200 ul de tampão de eluição (NaCl 1,000 mM).
   3. Incubar os grânulos num misturador rotativo durante 3 min.
   4. Colocar o tubo na cremalheira magnético durante 1 minuto.
   5. Com o tubo no lugar na prateleira magnética, remover o líquido com uma pipeta sem tocar as contas com a ponta da pipeta, e guardá-lo em um tubo de microcentrífuga limpo livre de DNase de 1,5 mL. Armazenar esta amostra em gelo.
   6. Repita os passos 1.8.1 - 1.8.4 uma vez, a segunda coleta no mesmo tubo (volume total será de 400 mL). Armazenar a amostra colectiva no gelo.
9. Precipitação com etanol (limpeza do ADN)
   1. Para cada noncaptured, lavagem e eluição fração from os passos anteriores, adicionar 1 ul de glicogénio (20 mg / mL, incluídos no kit), Volume 10/1 amostra de acetato de sódio 3M, pH 5,2 (por exemplo, 40 ul por 400 ul de amostra) e 2 volumes de amostra de 100% etanol (por exemplo, 800 mL por 400 mL de amostra). O volume total é de 1.241 mL.
   2. Misturar bem e incubar à temperatura de -80 ^° C durante pelo menos 2 h.
   3. Centrifuga-se o tubo durante 15 minutos a 11363 xg, a 4 ^o C.
   4. Cuidadosamente descartar o sobrenadante sem perturbar o sedimento.
   5. Adicionar 500 uL de etanol frio a 70%.
   6. Centrifugue o tubo por 5 min a 11.363 xg a 4 ^o C.
   7. Cuidadosamente descartar o sobrenadante sem perturbar o sedimento.
   8. Repita os passos 1.9.6 - 1.9.7 uma vez e remover qualquer sobrenadante residual.
   9. Ar seco o sedimento para ~ 5 min. Não secar completamente o sedimento.
   10. Ressuspender o sedimento de DNA em 37,5 uL de água sem ADNase ou outravolume apropriado de tampão.
   11. Coloque o DNA em gelo ou armazenar o DNA a -20 ^° C ou inferior até à sua utilização.
   12. Verifique se a quantidade total de DNA é superior a 20 ng, (por exemplo, usar o sistema de quantificação fluorimétrica, veja Materiais Tabela).

2. Sequencing

1. Use os fragmentos de DNA recuperados de procedimento MBDCap para o sequenciamento. Para cada amostra a ser sequenciada, uma quantidade total de 20 - 40 ng de ADN é necessária para a preparação da biblioteca.
2. Para a preparação da biblioteca de ADN-SEQ usar um sistema de construção da biblioteca totalmente automatizado (Ver Tabela Materiais) e seguir o protocolo padrão fornecido pelo fabricante. Seleccione fragmentos de DNA de tamanho 200-400 bp para a preparação da biblioteca. O sistema de preparação biblioteca de automação permite padronizar o procedimento de preparação da biblioteca e minimizar a variação entre as amostras devido à preparação manual.
3. privilegiada uso do adaptadorrs e diluir a concentração de 100x. Deste, use 10 mL para cada amostra.
4. Seguindo o passo 2.2, levar a cabo a reacção, e quantificar a biblioteca de ADN-SEQ usando um sistema de quantificação fluorométrico (Ver Tabela Materiais).
   1. Configurar o procedimento de amplificação PCR. A escolha de PCR Master Mix é melhorar a eficiência da PCR da região CG enriquecida do genoma. Num tubo de PCR, adicionar biblioteca 15? L DNA-SEQ, 25 pL de PCR mestre de mistura (Ver Tabela Materiais), 2 uL do iniciador de PCR mix (ver Materiais Tabela), e 8 mL _de H ₂ O.
5. Seguir a recomendação PCR do fabricante da PCR Master Mix, alterando vezes ciclismo e temperaturas para diferentes matérias-primas: 98 ^o C durante 45 s, ciclos X de 98 ^o C durante 15 s, 65 ^° C durante 30 s, 72 ^° C durante 30 s, e em seguida 72 ^° C durante 1 min. Segure a 4 ^o C.
   1. perform ciclos suficientes para acabar com 2-10 nM de DNA biblioteca (8 - 10 ciclos).
6. Limpar reacção de PCR com 1: 1 de PCR grânulos de purificação de acordo com o protocolo do fabricante (Ver Tabela Materiais).
7. Eluir com 20 - 30 ul de tampão de eluição.
8. Código de barras de cada amostra e piscina 4 amostras juntos em uma pista de uma célula de fluxo. Use o protocolo de leitura padrão de 50 bp única para sequenciação (Ver Tabela de Materiais).

3. Análise Bioinformatics

Nota: Além disso processar os arquivos fastq brutos obtidos a partir do sequenciamento para realizar o controle de qualidade e mapear as sequências de DNA curtas (lê) para o genoma.

1. Use Bowtie alinhador curto ler, ferramenta de mapeamento do genoma em cada arquivo de amostra fastq para mapear a lê para o genoma de referência. Muitas ferramentas de mapeamento de leitura estão disponíveis para realizar esta etapa e pode ser usado com base na preferência individual.
   1. Permitir até 2 inadequações no whole lê enquanto mapeamento para o genoma de referência. Relatar até 20 melhores alinhamentos para cada leitura. Por exemplo, usar gravata borboleta -v 2 --best -k 20 --chunkmbs 200 .
2. Use LONUT ⁹ para recuperar a multiplicar lê mapeado para melhorar a detecção das regiões enriquecidas. Para cada leitura, no máximo, um alinhamento seria recuperado quando se está perto o suficiente para qualquer um dos picos acima chamada. Por exemplo, use: perl lonut.pl -g hg19 -w 250 -p 99 . LONUT irá realizar os seguintes passos.
   1. alinhamentos dividido em mapeados exclusivamente lê e se multiplicam mapeados lê. Pico de chamadas na mapeados exclusivamente lê usando chamador pico CORREIA ^11. A opção -w 250 e -p 99 no comando exemplo são para cinto.
   2. Combine exclusivamente mapeado lê com o recuperado leituras obtidas a partir de LONUT, e o combinado lê, em formato CAMA, vai ser utilizado em análises posteriores. Calcule o número de Combined lê para a normalização futuro.
3. Bin as leituras. Extrato lê na região de interesse usando script perl: perl methyPipeline.pl -up -dn . Para cada arquivo de cama listados na sampleInfoFile, o script perl executa as tarefas abaixo.
   Nota: Ver o arquivo de codificação complementar para o script perl.
   1. Bin as leituras em tamanho de posição fixa, por exemplo, de 100 pb, por 24 cromossomos.
   2. Para cada gene especificado em refGeneFile, extrair os recipientes em torno do local de início da transcrição (SST), até ao comprimento que se estende especificado por -up opção e -dn. Por exemplo, extrair os dados da região que é a montante ea jusante 4 KB de TSS. No bin tamanho de 100 pb, este extraíram os dados têm 81 caixas, ea bin centro corresponde ao TSS.
   3. Para o gene da cadeia anti-sentido, inverter a região extraído da esquerda para a direita, para que caixas de upstream são sempre placed na esquerda com o script perl.
   4. Além disso dividir TSS ± ​​região de 4 kb em três sub-regiões: esquerda, a montante 4 kb a montante de 2 kb; Oriente, a montante ea jusante 2 kb; Direita, a jusante 2 kb a 4 kb.
   5. Para cada amostra, dividir a contagem de leituras em cada bin pelo número total de Combinado mapeados lê calculada em 3.3.3. A normalização vai eliminar a amostra específica ler diferenças e permite comparar as leituras entre amostras diferentes.
4. Executar o script de festa (como indicado no ficheiro de script) para executar o teste estatístico na normalizada leituras obtidas a partir da região esquerda, para detectar a metilação diferencial entre o grupo normal e o caso nas regiões descritas como na região do flanco.
   Nota: Consulte o arquivo de codificação complementar para o script bash. Com base nas regiões que são diferencialmente metilados, existem sete combinações:
   Toda a região: metilação diferencial ao longo ± região 4 kb do TSS
   middlmetilação diferencial tanto em região média e esquerda: eLeft
   metilação diferencial tanto em meio e direita região: MiddleRight
   metilação diferencial em ambas as regiões esquerda e direita: flanco
   Esquerda: metilação diferencial na região deixou apenas
   metilação diferencial na região direito só: Direito
   metilação diferencial na região média apenas: Médio
5. Use o mesmo script bash para visualizar a metilação diferencial em uma trama tornado. Traçar os hiper e hipo genes metilados em um painel separado, e classificar na ordem da combinação região, conforme descrito em 3.4, respectivamente.

Resultados

Usámos MBDCap-SEQ de estudar as alterações de metilação do DNA em um grande número de pacientes de diversos tipos de cancro incluindo cancro da mama ^{12, 13} do endométrio, próstata ^14, e cancros do fígado, entre outros. Aqui demonstramos algumas informações a partir do estudo do cancro da mama publicou recentemente ^12. Neste exemplo, utilizou-se toda a abordagem de sequenciação do genoma para identificar ilhas CpG que são diferencialmente metilados no tumor em relação à perpendicul...

Discussão

A técnica MBDCap-seq é uma abordagem de enriquecimento de afinidade ^3, considerada como uma alternativa de custo eficaz quando se investiga coortes com um grande número de pacientes ^15. O gasoduto apresentado aqui descreve uma abordagem abrangente, desde a colheita de amostra para análise de dados e interpretação. Um dos passos mais importantes é a criação de um procedimento de amplificação por PCR para melhorar a eficiência da PCR das regiões GC enriquecido no genoma como este é onde ocorre a metilação...

Materiais

Name	Company	Catalog Number	Comments
Methylminer DNA enrichment Kit	Invitrogen	ME10025
Dynabeads M-280 Streptavidin	Invitrogen	112-05D
Bioruptor Plus Sonication Device	diagenode	B01020001
3 M sodium acetate, pH 5.2	Sigma	S7899	100 mL
SPRIworks Fragment Library System I	Beckman Coulter	A50100	Fully automated library construction system
Adapter Primers	Bioo Scientific	514104	PCR primer mix
Qubit	Invitrogen	Q32854	Fluorometric Quantitation System
PCR master mix	KAPA scientific	KK2621	PCR master mix
AMPure XP	Beckman Coulter	A63881	PCR Purification beads
EB Buffer	Qiagen	19086
HiSeq 2000 Sequencing System	Illumina