Summary
Aqui, apresentamos um protocolo para permitir que a quantificação exata de metilação de DNA (mtDNA) mitocondrial. Neste protocolo, descrevemos uma digestão enzimática de DNA com BamHI juntamente com um pipeline de análise de bioinformatic que pode ser usado para evitar a superestimação dos níveis de metilação do DNA mitocondrial causada pela estrutura secundária do DNA mitocondrial.
Abstract
Quantificação de metilação do DNA pode ser conseguida usando sequenciamento de bissulfito, que aproveita a propriedade de bissulfito de sódio para converter unmethylated citosina em uracil, num contexto single-stranded DNA. Sequenciamento de bissulfito pode ser direcionado (usando PCR) ou executadas no genoma inteiro e fornece quantificação absoluta de metilação de citosina na única base-resolução. Dada a natureza distinta de DNA nuclear e mitocondrial, nomeadamente a estrutura secundária, devem ser feitas adaptações de métodos de sequenciamento de bissulfito para investigar a metilação de citosina no DNA mitocondrial. Estrutura secundária e terciária de DNA mitocondrial de fato pode levar a bissulfito artefatos levando a falso-positivos devido a acesso de pobres incompleta desnaturação de bissulfito de ADN single-stranded de sequenciamento. Aqui, descrevemos um protocolo usando uma digestão enzimática de DNA com BamHI juntamente com pipeline de análise de bioinformatic para permitir a quantificação exata da citosina níveis de metilação no DNA mitocondrial. Além disso, podemos fornecer diretrizes para projetar os primers de sequenciamento de bissulfito específicos para DNA mitocondrial, para evitar a segmentação segmentos NUclear mitocondrial indesejáveis (NUMTs) inserido no genoma nuclear.
Introduction
O genoma mitocondrial é uma estrutura circular, double-stranded de aproximadamente 16,5 kg base (kb), constitui um pesado e um fio de luz. O genoma mitocondrial está presente em várias cópias dentro de cada célula, herdado maternalmente e codifica componentes essenciais da cadeia respiratória complexos1. Semelhantes aos genomas bacterianas e ao contrário do genoma nuclear, o genoma mitocondrial é organizado em numerosas estruturas secundárias e terciárias, como em estruturas enrolado e supercoiled2, que pode tornar o acesso difícil durante o sequenciamento experiências3.
No núcleo, metilação do DNA é uma marca epigenética extensivamente estudada que desempenha um papel em vários processos, nomeadamente na regulação da expressão gênica. Em genomas de mamíferos, a metilação do DNA ocorre principalmente na posição 5 do anel pirimidina de deoxycytidines, principalmente de dinucleotídeos CG (ou CpG). Metilação de citosina é encontrada em 70% de todos os CpG no genoma de células somáticas e contas para ~ 1% do totalque 4bases de DNA. Methylations de DNA também tem sido descrita em contextos não-CpG, como CpA, CpT e CpC e existem em várias quantidades no DNA nuclear, com valores até 25% de todas as Citosinas metiladas em células-tronco embrionárias5,6,7.
Enquanto o methylation do cytosine do genoma nuclear é amplamente aceito, a existência de metilação de DNA (mtDNA) mitocondrial é ainda controversa. O primeiro estudo investigar methylation DNA mitocondrial foi realizado em pilhas cultivadas onde a metilação do DNA mitocondrial foi prontamente detectada, embora em níveis mais baixos em comparação com nuclear DNA8. Em células humanas e murino, metilação de DNA mitocondrial também foi detectada em níveis baixos (2-5%). Utilizando ensaios contando com 5 metilcitosina captura, tais como imunoprecipitação de DNA metilada (MeDIP) seguida por PCR quantitativo, metilação de DNA mitocondrial também foi detectada em vários rato e humanos e células linhas9,10, 11,12. Usando anticorpos contra 5-metilcitosina em um ensaio de ELISA ou espectrometria de massa, níveis substanciais de metilação do DNA foram detectados purificada fracções mitocondrial13,14,15, 16. no entanto, a maioria dos ensaios nos estudos acima mencionados utilizado técnicas que não foram projetadas para fornecer a quantificação absoluta de metilação do DNA na única base-resolução.
Análise de metilação de DNA quantitativo e resolutiva pode ser alcançado por uma técnica chamada "bissulfito de sequenciamento", que aproveita a propriedade de bissulfito de sódio para converter unmethylated citosina em uracil em single-stranded DNA contexto17 . Usando o sequenciamento de bissulfito, uma constelação de estudos detectou a presença de metilação de citosina em vários níveis. DNA mitocondrial de metilação na região D-loop, a 12S ou região 16S prontamente foi detectado em humano18,19,20,21,22,mouse de23 e24 tecidos e células, no entanto, com uma variabilidade de intrigante, de 1-20% do totais citosinas através de estudos.
Em comparação com estes numerosos estudos, poucos estudos, incluindo do nosso grupo, têm disputado a presença de DNA mitocondrial metilação3,25,26,27 ou questionado a relevância biológica de níveis muito baixos de mtDNA níveis (abaixo de 2%)28. Recentemente, nós relatamos a observação de um artefato de bissulfito-sequenciamento potencial em toda mitocondrial bissulfito sequenciamento3. Nós fornecemos evidências de que a estrutura secundária do DNA mitocondrial pode levar a resultados falsos positivos sequenciamento de bissulfito, desse modo, superestimando os níveis de metilação. Nós fornecemos aqui um protocolo para evitar um artefato de bissulfito-conversão de DNA mitocondrial. Este protocolo utiliza uma simples digestão enzimática do DNA para interromper estruturas secundárias do DNA mitocondrial e permitir o acesso total ao bissulfito, seguindo um protocolo de sequenciamento de bissulfito. Além disso, nós fornecemos um pipeline de bioinformatic acompanhamento para a análise do sequenciamento de bissulfito.
Subscription Required. Please recommend JoVE to your librarian.
Protocol
1. tratamento enzima de restrição
- Linearizar mtDNA tratando o DNA humano total com a enzima de restrição BamHI, que corta na posição 14258 no DNA mitocondrial humano.
Nota: Para mouse DNA, use a enzima de restrição BglII sob condições idênticas como abaixo. - Para cada amostra onde a metilação do DNA mitocondrial é para ser avaliado, preparar um tubo de reacção de 0,2 mL e adicionar o seguinte mistura: 3 μg DNA genômico (quantificado pelo fluorometry), 15 μL tampão 3, 3 μL BamHI e água acima de 150 μL
- Coloque os tubos em um termociclador por 4 h a 37 ° C.
2. bissulfito conversão
- Converta o DNA BamHI-Tratado usando bissulfito de sódio, conforme descrito em outro lugar,29.
- Use 200-500 ng DNA BamHI-Tratado para cada reação de conversão em um volume total de 20 μL. A recuperação do DNA é 50-70%.
- Eluir o DNA de bissulfito-convertido em tampão de eluição 10 μL.
- Quantificar o ADN single-stranded por fluorometry.
- Opcional: Armazenar DNA bissulfito-convertido a-20 ° C antes de prosseguir para o restante do protocolo. Para armazenamento a longo prazo, armazenar DNA bissulfito-convertido a-80 ° C
3. projeto de sequenciamento de Primers de bissulfito
- Desenho de primers para sequenciamento de bissulfito, usando as ferramentas on-line Methprimer30 e31,de BiSearch32.
Nota: Methprimer e BiSearch permitem para procurar sequências de vinculação da primeira demão sobre sequências genomic onde citosinas são convertidas em timinas, similarmente ao pós-bissulfito tratamento. - Para a amplificação da ideal e imparcial, evite dinucleotídeos CpG em primers e projetar o tamanho do amplicon entre bp 100-300.
- Que sejam projetados primers específicos para o DNA mitocondrial, usando a função do BiSearch Busca de Primer. Inserir as primeiras demão já projetado bissulfito convertido, marque a caixa de bissulfito e incluem um genoma de referência e parâmetros PCR.
Nota: Será exibida uma lista de possíveis produtos PCR. - Copie o top sequências de 1 a 5 da primeira demão e colá-los em um formulário de pedidos para a síntese do oligonucleotide.
Nota: Uma lista de humanos e sequências de cartilha de DNA mitocondrial de rato validado no laboratório é exibido na tabela 1. - Iniciadores de teste usando bissulfito convertido PCR após eletroforese em gel de agarose, conforme descrito na etapa 4. Testar e otimizar todos os pares da primeira demão separadamente e verifique o tamanho correto do amplicon aparece sobre o gel de agarose.
Nota: Se multiplexação primários são necessários, adicionar múltiplos primers para uma única reação de PCR e visualizar os produtos PCR na electroforese do gel do agarose ou, um método mais resolutiva como eletroforese em gel de poliacrilamida, para distinguir produtos de semelhante tamanho.
4. bissulfito convertido do PCR e extração do Gel
- Para amplificar as regiões de interesse, realize PCR usando uma começo quente Taq polimerase. Brevemente, mistura 100 ng convertido DNA, 500 cartilhas de sentido e anti-sentido µM, 1 µ l do dNTP mix (10 µ m de cada dNTP) e polimerase em 7,5 U/reação em um volume total de 50 µ l. executar o PCR com as seguintes condições: 5 min a 95 ° C (60 s 94 ° C; 60 s 55 ° C *; 60 s 72 ° C)
35 ciclos; 10 min 72 ° C. Também usar primers separadamente ou multiplexado. - Opcional: Quando as primeiras demão de multiplexação, uso 200 ng convertido DNA e as seguintes condições de ciclismo: 5 min 95 ° C (60 s 94 ° C; 90 s 55 ° C *; 90 s 72 ° C) 35 ciclos; 10 min 72 ° C.
Nota: A temperatura pode variar dependendo da temperatura de fusão de primers. - Produtos PCR sujeitos ao gel eletroforese em agarose 2% em 100 volts.
- Suave luz UV, extrair produtos PCR com um bisturi e adicionar para microtubos. Não exponha os produtos de PCR de luz UV em excesso para evitar danos ao DNA. Evitar o tempo de exposição acima de 30 s.
- Purifica os produtos PCR usando gel de purificação conforme descrito anteriormente3.
- Quantificar o DNA purificado da etapa 4.5 usando fluorometry. Prossiga para o passo 5 ou amostras de armazenamento-20 ° c.
35 ciclos; 10 min 72 ° C. Também usar primers separadamente ou multiplexado.5. bissulfito sequenciamento biblioteca preparação
- Execute a preparação da biblioteca de produtos PCR bissulfito-convertido.
Opcional: multiplex produtos do PCR nesta fase. - Executar reparação final usando até 100 ng do produto do PCR. Adicione 3 µ l final prep enzima e amortecedor da reação reparação final 6,5 µ l (10 x) para 55,5 µ l de produto PCR em um tubo de 0,2 mL. Coloque e incubar o tubo em um termociclador por 30 min a 20 ° C, seguido de 30 minutos a 65 ° C.
- Ligate adaptadores de sequenciamento, misturando o DNA final-reparado com 15 mix µ l Blunt/TA ligase, 2,5 µ l adaptador e potenciador de ligadura 1 µ l. Se usar < 100 produto PCR do ng como entrada, diluir adaptador 10 vezes.
- Coloque a mistura (etapa 5.3) em um termociclador e incubar 15 min a 20 ° C. Pausar o termociclador e adicione 3 µ l usuário enzima ao tubo. Misture e coloque o tubo de volta no termociclador e incubar 15 min a 37 ° C.
- Tamanho selecione o adaptador-ligado DNA usando contas de fase sólida reversível imobilização (SPRI) com razões diferentes das esferas para DNA, dependendo do tamanho do produto do PCR.
- Adicione 13,5 µ l H2O para adaptador-ligado DNA de passo 5.4 e 55 µ l resuspended SPRI missanga. Incube a temperatura ambiente por 5 min e coloque o tubo num suporte magnético. Quando a solução é clara, transferir o sobrenadante para um tubo novo e descartar o tubo contendo os grânulos.
Nota: O sobrenadante contém o DNA de adaptador-ligados e grandes fragmentos indesejados são vinculados aos talões descartados. - Adicione 25 µ l resuspended SPRI grânulos para o sobrenadante da etapa 5.6, misturar e incubar durante 5 min à temperatura ambiente. Colocar o tubo no stand magnético e descartar o sobrenadante, quando a solução é clara.
Nota: O sobrenadante Descartado contém DNA indesejado, Considerando que o DNA de adaptador-ligado está vinculado a grânulos. - Enquanto o tubo é o suporte magnético, adicione 200 µ l 80% de etanol (preparado) para lavar os grânulos. Incubar 30 s e remover e descartar o sobrenadante. Repita esta etapa para um total de 2 lavagens.
- Ar secos grânulos por 5 min, enquanto o tubo é em suporte magnético com uma tampa aberta.
- Retire o tubo do ímã, Eluir o DNA em tampão de eluição 23 µ l e misture. Incube a temperatura ambiente por 2 min.
- Coloque o tubo num suporte magnético e quando a solução é clara, 2 µ l de transferência para um prato PCR de 96 poços para a amplificação do pre-PCR (etapa 5.14).
- Transferi o restante de aproximadamente 21 µ l para um novo tubo de 0,2 mL para amplificação por PCR (etapa 5.17).
Nota: Certifique-se de não transferir qualquer grânulos como grânulos podem inibir reações enzimáticas dos próximos passos. - Prossiga para a etapa 5.14 ou amostras de armazenamento-20 ° c.
- Execute a amplificação do pre-PCR por RT-qPCR para estimar o número de ciclos necessários para amplificar o DNA do adaptador-ligados. Por reação, adicionar o seguinte no PCR de 96 poços da placa contendo 2 µ l DNA (etapa 5.11): cartilha de índice 0,4 µ l, 0,4 µ l universal do PCR da primeira demão, 7,2 µ l H2O, 10 µ l RT-qPCR mestre mix.
- Coloque a placa em uma máquina PCR em tempo real e sujeito a placa para as seguintes condições: 30 s a 98 ° C (10 s 98 ° C; 75 s 65 ° C) 20 ciclos; 5 min 65 ° C.
- Estime o número de ciclos de amplificação necessária para cada amostra subtraindo um valor de Ct para o Ct-valor do pré-amplificação RT-qPCR correspondente ao fim da fase linear da reação de PCR.
- Amplificar o DNA adaptador-ligados da etapa 5.12 misturando: 21 µ l DNA, 25 µ l PCR master mix, 1 µ l índice da primeira demão, 1 µ l do PCR Universal da primeira demão e 2 µ l H2O. Em um termociclador, sujeito cada amostra para as seguintes condições: 30 s a 98 ° C (10 s 98 ° C; 75 s 65 ° C) [Ct calculado da etapa 5.16] ciclos; 5 min 65 ° C.
Nota: Lembre-se de utilizar somente uma primeira demão de índice por amostra para permitir a multiplexação. O índice é inserido durante a PCR. - Purificar os grânulos SPRI DNA amplificados e eluir em tampão de eluição 22 µ l.
- Controlar a qualidade da biblioteca pelo método de gel de eletroforese ou alternativa. Procure por tamanhos de pico de biblioteca correta (tamanho do produto (s) PCR mais adaptador).
- Estimar o tamanho médio de pares de base do produto de biblioteca pela análise do esfregaço: em avançada configurações globais, clique duas vezes em "tabela" sob análise do esfregaço. Definir a região de 100-1000 bp e clicar Okey. Sob a Mesa da região, a região definida aparece e biblioteca de tamanho médio (bp) é calculada.
- Quantificar as bibliotecas por fluorometry. Prossiga para o passo 6 ou bibliotecas de armazenamento-20 ° c.
6. na próxima geração de sequenciamento
- Calcule a concentração de nM de bibliotecas, usando a fórmula:
- Diluir as bibliotecas para o mesmo nM nM de concentração por exemplo 2 (escolher 4 nM, 2 nM, 1 nM, ou 0,5 nM, dependendo da concentração de biblioteca). Todas as bibliotecas para alcançar uma piscina de 2 bibliotecas nM da piscina.
- Desnaturar e diluir bibliotecas seguindo o guia de instrumento de sequenciamento. Em suma: desnaturar o pool nM 2 adicionando 0,2 N NaOH e incubar a temperatura ambiente por 5 min. diluir a piscina desnaturada para uma diluição final de 22:00. Desnaturar e diluir uma biblioteca de controle comercialmente disponível para uma diluição final de 12:05. Para um novo tubo, misture o pool de 22:00 com a biblioteca de controle 12:05 20%.
Nota: Conversão de bissulfito introduz complexidade muito baixa. 20% de controle de biblioteca spike resultará em melhor qualidade de sequenciamento. - Execute um 150 bp emparelhado-final em sequência usando um instrumento de sequenciamento profunda.
7. computacional análise - níveis de metilação de estimativa
- Gerar arquivos de .fastq (aqui chamados de sample1_R1.fastq.gz e sample1_R2.fastq.gz), gerar um índice de bismark do genoma inteiro de interesse (aqui em uma pasta chamada bismarkIndex) e certifique-se de que a sequência genômica de chrM está disponível em fasta formato (aqui em uma pasta chamada chrM).
- Pré-processar leituras para remover adaptadores e viés introduzido por primers: Use a ferramenta Trim_galore33. Remova os dois nucleotídeos da extremidade 5' do lê ambos frente e verso.
trim_galore - clip_R1 2 - clip_R2 2 -o. / --trim1 - emparelhado sample1_R1.fastq.gz sample1_R2.fastq.gz - Alinhar o pre-processed leituras para o genoma inteiro usando Bismark34
Bismark - andorinha..--bam -o. /. / bismarkIndex / o-1./sample1_R1_val_1.fq.gz -2./sample1_R2_val_2.fq.gz
Nota: Se outro alinhador é usado, certifique-se de que lê o que não pode ser alinhado com exclusividade é descartadas. - Extrato lê o mapeamento para o cromossomo M, aqui usando Samtools35
samtools classificar -l 0 - O BAM -o sample1_sorted.bam
sample1_R1_val_1_bismark_bt2_pe.bam
samtools índice sample1_sorted.bam
samtools Ver os chrM de sample1_sorted.bam -u | samtools classificar - n - O BAM -o sample1_chrM.bam - Extrair as informações de metilação
bismark_methylation_extractor -p - CX..--cytosine_report..--abrangente..--genome_folder. /./sample1_chrM.bam chrM - Use o .bedGraph, os COV ou os arquivos de CX_report.txt para uma análise mais aprofundada. Use o arquivo sample1_chrM.CX_report.txt para visualização e teste.
- Execute o recorte ainda mais na etapa de pré-processamento, se várias regiões são interrogadas e um viés de primeira demão pode ser visível nas tramas M-viés geradas durante o mapeamento. Consulte a documentação do Bismark para obter mais informações.
8. computacional análise - teste de diferenças
- Certifique-se de que o arquivo CX_report.txt contém 7 colunas, cromossomo (chrM aqui), posição, vertente, número de leituras metilados, número de leituras unmethylated, C contexto e sequência circundante, para todos os C no cromossomo M.
- Como o CX_report cobre chrM em sua totalidade, verifique o subconjunto para a região (ões) que é amplificado.
- Use testes não-paramétricos para as diferenças entre as amostras, tais como o teste exato de um pescadores do C individual ou um teste de sinal-para toda uma região.
Nota: Quantificação será frequentemente perto de 0% de metilação, que é a razão por que supondo que a normalidade não é apropriado. - Da mesma forma, para visualização, ou Visualizar quantiles ou calcular intervalos de confiança de proporção binomial.
Subscription Required. Please recommend JoVE to your librarian.
Representative Results
Duas etapas neste protocolo são cruciais ao investigar a metilação do DNA mitocondrial. 1) a abertura da estrutura secundária e 2) o projeto de primers de específicas de DNA mitocondriais.
Digerindo o DNA genômico humano com a enzima de restrição BamHI (Figura 1), a estrutura de DNA mitocondrial será cortada no nucleotídeo posição 14.258 e estrutura secundária será aberta.
O comprometimento possível da conversão de bissulfito, com uma estrutura secundária do mtDNA intacta é muito provável que superestima a taxa de unconversion de DNA mitocondrial. Por sequenciamento de bissulfito, a taxa de unconversion do DNA mitocondrial foi investigada, não digerida e digerido o DNA total de células do músculo esquelético humano, em 5 diferentes regiões do genoma mitocondrial: o Loop de deslocamento (D-Loop), o ribossoma fenilalanina e 12S de tRNA (tRNA - F + 12S), Ribossoma 16S (16S), NADH desidrogenase 5 (CD5) e b (CYTB) codificação de mRNA citocroma (Figura 2). Taxa de unconversion variou de 0 - 15,1% de todas as regiões investigadas para o mtDNA não digerido. No entanto, quando o DNA foi digerido antes da conversão de bissulfito, a taxa de unconversion caiu para um máximo de 1% de todas as regiões (Figura 3). Assim, ilustrando a importância da digestão de BamHI antes da conversão de bissulfito para evitar a superestimação dos níveis de metilação do DNA mitocondrial.
Contaminação de DNA nuclear é inevitável quando isolar o genoma mitocondrial e desde que quase todo o genoma mitocondrial foi replicado no genoma nuclear, tem dado origem a inserções nucleares do genoma mitocondrial, denominado NUclear Segmentos mitocondrial (NUMTs)36. Para garantir que os níveis de metilação de DNA analisados correspondem ao DNA mitocondrial e não NUMTs, é de importância ultraperiférica para projetar primers específicos para o genoma mitocondrial. A Figura 4 mostra como verificar a especificidade da primeira demão e tabela 1 exibe humana e sequências de rato mtDNA da primeira demão.
| Nome | Sequência de | Posição | Tamanho do amplicon | Temperatura do recozimento | Strand | ||
| Humana | |||||||
| D-Loop | F: AAATCTATCACCCTATTAAC | ||||||
| R: GTGGAAATTTTTTGTTATGATGT | 6-298 | 292 | 55° C | Antisentido | |||
| D-Loop | F: CATAACAAAAAATTTCCACCAAAC | ||||||
| R: GGGAAAATAATGTGTTAGTT | 279-458 | 179 | 55° C | Antisentido | |||
| 12S + TF | F: TTTATATAACTTACCTCCTC | ||||||
| R: GTGTTTGATGTTTGTTTTTTTTG | 577-765 | 188 | 55° C | Antisentido | |||
| 16S | F: AATAAATTTATAGGTTTTTAAATTATTAAAT | ||||||
| R: TAACTAATAAAATCTTAACATATACTACTC | 2763-2873 | 110 | 53° C | Sentido | |||
| CYTB | F: GGTATTATTTTTTTGTTTGTAATTATAGTA | ||||||
| R: CCTCAAATTCATTAAACTAAATCTATCC | 15091-15243 | 152 | 55° C | Sentido | |||
| Mouse | |||||||
| 12S | F: ACACATACAAACCTCCATAAAC | ||||||
| R: GAGGTATAATTTAGTTAAAT | 111-287 | 176 | 53° C | Antisentido | |||
| 16S | F: AAATTCCAATTCTCCAAACATAC | ||||||
| R: TTTGGATTTTTTTTTTAGGT | 1749-1896 | 147 | 53° C | Antisentido | |||
| CYTB | F: CTACAAAAACACCTAATAACAAAC | ||||||
| R: AGTGTATGGTTAAGAAAAGA | 14129-14307 | 178 | 55° C | Antisentido | |||
| 16S | F: AGAGAAATAGAGTTATTTTATAAATAAG | ||||||
| R: CTAAAAACAAAATTTTAAATCTTAC | 2576-2727 | 151 | 55° C | Sentido | |||
| CD5 | F: TTAAGTTAATTAGGTTTGATAATAGTGA | ||||||
| R: TAAAACCCTATTAAAAATAATATTCCTATA | 12659-12910 | 251. | 55° C | Sentido | |||
| CYTB | F: TGGGTTTTTTTTAGGAGTTTGTTTA | ||||||
| R: CAATAAAAATACTCCAATATTTCAAATTTC | 14242-14504 | 262. | 55° C | Sentido | |||
Tabela 1: sequências de Primer para o ser humano e do rato mtDNA. Observe que recozimento temperaturas do PCR primers variam devido às diferenças da temperatura de fusão de primers.
Figura 1 : Gel de electroforese de DNA BamHI-digerido ou não digerido. DNA genômico de músculo esquelético humano foi tratado ou tratado com BamHI por 4h a 37° C e eletroforese em gel foi executado em um gel de agarose 1,5%. Clique aqui para ver uma versão maior desta figura.
Figura 2 : Mapa do genoma mitocondrial humano. Regiões do genoma mitocondrial investigado pelo sequenciamento de bissulfito e o site de restrição da enzima BamHI são exibidas. O DNA mitocondrial contém 13 mRNAs, 2 RNAs ribossomal e 22 RNAs de transferência. Abreviaturas: PH1: promotor pesado-vertente 1; PH2: Heavy-vertente promover 2; PL: promotor de luz-fio; OH: Origem de replicação do pesado-strand; ÓL: origem de replicação do luz-strand. Clique aqui para ver uma versão maior desta figura.
Figura 3 : BamHI digestão antes da conversão de bissulfito reduz a taxa de unconversion do DNA mitocondrial em células de músculo esquelético humano. Por sequenciamento de bissulfito, não digerida e digerida porcentagem de unconversion DNA mitocondrial foram interrogados em cinco diferentes regiões do genoma mitocondrial em células de músculo esquelético humano (N= 3). A praça cheia representa sites de CpG, e o quadrado aberto representa sites não-CpG. Os resultados são apresentados com um intervalo mín-máx e um teste de sinal foi usado para testar as diferenças unconversion significativas. D-Loop (6-298) P= 1.83E-42; D-Loop (279-458): P= 4.55E-13; tRNA - F + 12S: P= 8.38E-16; 16s: P= 5.91E-06; CD5: P= 1.70E-05; CYTB: P= 2.69E-10. D-Loop (6-298) inclui a origem de replicação e tRNA - F + 12S inclui promotor forte vertente 2. Estes dados foram publicados em outro lugar 3. Clique aqui para ver uma versão maior desta figura.
Figura 4 : O uso de Bisearch para verificar a especificidade da primeira demão, no sentido do genoma mitocondrial. A) sequências da primeira demão bissulfito-convertidos são adicionadas à função BiSearch busca de Primer (http://bisearch.enzim.hu/?m=genompsearch). A bissulfito caixa está marcada, e um genoma de referência está seleccionado. B) exemplo de produtos PCR potenciais gerado no sentido e antisentido correntes. Clique aqui para ver uma versão maior desta figura.
Subscription Required. Please recommend JoVE to your librarian.
Discussion
Aqui, nós fornecemos um protocolo de bissulfito-sequenciamento que é projetado especificamente para interrogar a metilação do DNA mitocondrial. As diferenças com bissulfito-sequenciamento protocolos usados para DNA genômico situa-se na utilização de uma enzima de restrição prévia etapa de digestão e uma análise de bioinformatic acórdão fora falsos positivos decorrentes NUMT sequências.
Nós fornecemos um protocolo para evitar artefatos de bissulfito-sequenciamento ao investigar a metilação do DNA mitocondrial. Artefatos de sequenciamento de bissulfito levando a falsos positivos foram anteriormente relatados37. Remoção insuficiente de proteínas acessórias de DNA por proteinase que k tem sido descritos37. Desnaturação incompleta ou parcial reannealing pode levar a trechos de conversão incompleta, possivelmente diminuindo acesso de bissulfito de single-stranded DNA37. O artefato mais tarde poderia ser em jogo no DNA mitocondrial onde uma estrutura secundária iria prejudicar o acesso de citosinas. O melhor de nosso conhecimento, a adição de uma etapa de digestão enzimática de restrição antes da conversão de bissulfito no contexto de estimar a metilação do DNA mitocondrial tem sido relatada primeiramente em 20163,28. Em nossas mãos também, prévia digestão com uma enzima de restrição single-cortador drasticamente reduzido deteção de citosinas não convertida3,28.
Aqui, nós fornecemos um pipeline de bioinformatic para analisar resultados de sequenciamento de bissulfito e para detectar artefatos de sequenciamento de bissulfito no contexto de sequenciamento de genoma mitocondrial toda bissulfito. Este método origina-se de nossa observação que mtDNA taxa de conversão não correlacionou-se inversamente com profundidade de3,28de sequenciamento. Esta observação sugere que a estrutura secundária e terciária mesmo do mtDNA prejudica a fragmentação do DNA mitocondrial em regiões específicas durante a construção de bibliotecas de sequenciamento e também poderia prejudicar o completo de bissulfito de sódio para DNA mitocondrial. Uma vez que a digestão com uma enzima de um cortador libera estruturas secundárias e terciárias, nossos resultados indicam que a estrutura do DNA mitocondrial é uma fonte de artefato de bissulfito e validar os métodos descritos para a quantificação de metilação de citosina no DNA mitocondrial.
A inevitável contaminação de DNA nuclear em preparação mitocondrial constitui um desafio. Isto leva à contaminação em sequências NUMT em leituras de sequenciamento, que podem influenciar-estimativa de metilação de citosina no DNA mitocondrial, mesmo quando realizando direcionados e sequenciamento de genoma não toda mitocondrial bissulfito. Para evitar tal preconceito, o desenho de primers de sequenciamento de bissulfito exclusivo permite excluir contaminação por metilação em NUMTs. NUMTs abrangem determinadas regiões do ser humano e o genoma mitocondrial de rato com 100% de identidade, e é, portanto, impossível projetar primers exclusivos em certas regiões do DNA mitocondrial. Ao realizar o sequenciamento de genoma mitocondrial toda bissulfito, sequenciamento tempo leituras podem facilitar a discriminação entre o DNA mitocondrial e NUMTs. -Alinhamento de leituras para o genoma inteiro e utilize apenas leituras exclusivamente mapeadas para o genoma mitocondrial podem ainda garantir a detecção de metilação de DNA mitocondrial.
Finalmente, este protocolo pode ser usado para além de genomas mitocondriais, por exemplo, em outras instâncias do sequenciamento de bissulfito, onde a estrutura secundária do DNA representa uma fonte potencial de artefato. A associação entre níveis de metilação e sequenciamento de profundidade pode ser investigada para avaliar a necessidade de uma digestão com uma enzima de restrição antes de sequenciamento de bissulfito complexas de sequências de ADN.
Subscription Required. Please recommend JoVE to your librarian.
Disclosures
Os autores têm sem interesses financeiros concorrentes para declarar.
Acknowledgments
A Novo Nordisk Fundação Centro de pesquisa metabólica básica é um centro de pesquisa independente da Universidade de Copenhagen, parcialmente financiados por uma doação irrestrita da Novo Nordisk Foundation.
Materials
| Name | Company | Catalog Number | Comments |
| BamHI | New England BioLabs | # R0136 | |
| EZ DNA methylation-lightning kit | Zymo Research | # D5030 | |
| Qubit ssDNA assay | Thermo Fisher Scientific | # Q10212 | |
| Qubit assay tubes | Thermo Fisher Scientific | # Q32856 | |
| HotStarTaq plus DNA polymerase kit | Qiagen | # 203603 | |
| QIAquick Gel Extraction Kit | Qiagen | # 28704 | |
| NEBNext Ultra DNA Library Kit for Illumina | New England BioLabs | # E7370S | |
| NEBNext Multiplex Oligoes for Illumina | New England BioLabs | # E7335S | |
| AMPure XP Beads | Beckman Coulter | # A63881 | |
| High Sensitivity DNA chip | Agilent | # 5067-4626 | |
| Qubit high sensitivity dsDNA assay | Thermo Fisher Scientific | # Q33230 | |
| MiSeq reagent kit v2 300 cycles | Illumina | # MS-102-2003 | |
| PhiX control v3 | Illumina | # FC-110-3001 | |
| Sodium hydroxide | Sigma | # S5881 | |
| Thermal cycler C1000 | Biorad | # 1851148 | |
| CFX96 Real-Time PCR detection system | Biorad | #1855195 | |
| Qubit Fluorometer | Thermo Fisher Scientific | # Q33226 | |
| Bioanalyzer 2100 | Agilent | # G2939BA | |
| MiSeq instrument | Illumina | # SY-410-1003 |
References
- Falkenberg, M., Larsson, N. -G., Gustafsson, C. M. DNA replication and transcription in mammalian mitochondria. Annu Rev Biochem. 76, 679-699 (2007).
- Kolesar, J. E., Wang, C. Y., Taguchi, Y. V., Chou, S. -H., Kaufman, B. A. Two-dimensional intact mitochondrial DNA agarose electrophoresis reveals the structural complexity of the mammalian mitochondrial genome. Nucleic Acids Res. 41 (4), e58 (2013).
- Mechta, M., Ingerslev, L. R., Fabre, O., Picard, M., Barres, R. Evidence Suggesting Absence of Mitochondrial DNA Methylation. Front Genet. 8, 166 (2017).
- Ehrlich, M., Gama-Sosa, M. A., et al. Amount and distribution of 5-methylcytosine in human DNA from different types of tissues of cells. Nucleic Acids Res. 10 (8), 2709-2721 (1982).
- Lister, R., Pelizzola, M., et al. Human DNA methylomes at base resolution show widespread epigenomic differences. Nature. 462 (7271), 315-322 (2009).
- Yan, J., Zierath, J. R., Barres, R. Evidence for non-CpG methylation in mammals. Exp Cell Res. 317 (18), 2555-2561 (2011).
- Patil, V., Ward, R. L., Hesson, L. B. The evidence for functional non-CpG methylation in mammalian cells. Epigenetics. 9 (6), 823-828 (2014).
- Nass, M. K. Differential methylation of mitochondrial and nuclear DNA in cultured mouse, hamser and virus trasnforemd hamser cells in vivo and in vitro methylation. J Mol Biol. 80 (1), 155-175 (1973).
- Shock, L. S., Thakkar, P. V., Peterson, E. J., Moran, R. G., Taylor, S. M. DNA methyltransferase 1, cytosine methylation, and cytosine hydroxymethylation in mammalian mitochondria. PNAS. 108 (9), 3630-3635 (2011).
- Bellizzi, D., D'Aquila, P., et al. The Control Region of Mitochondrial DNA Shows an Unusual CpG and Non-CpG Methylation Pattern. DNA Res. 20 (6), 537-547 (2013).
- Jia, Y., Li, R., et al. Maternal Low-Protein Diet Affects Epigenetic Regulation of Hepatic Mitochondrial DNA Transcription in a Sex-Specific Manner in Newborn Piglets Associated with GR Binding to Its Promoter. PLoS ONE. 8 (5), e63855 (2013).
- Jia, Y., Song, H., Gao, G., Cai, D., Yang, X., Zhao, R. Maternal Betaine Supplementation during Gestation Enhances Expression of mtDNA-Encoded Genes through D-Loop DNA Hypomethylation in the Skeletal Muscle of Newborn Piglets. J Agric Food Chem. 63 (46), 10152-10160 (2015).
- Infantino, V., Castegna, A., Iacobazzi, F., Spera, I. Impairment of methyl cycle affects mitochondrial methyl availability and glutathione level in Down's syndrome. Mol Genet Metab. 102 (3), 378-382 (2011).
- Chen, H., Dzitoyeva, S., Manev, H. Effect of valproic acid on mitochondrial epigenetics. Eur J Pharmacol. 690 (1-3), 51-59 (2012).
- Dzitoyeva, S., Chen, H., Manev, H. Effect of aging on 5-hydroxymethylcytosine in brain mitochondria. Neurobiol Aging. 33 (12), 2881-2891 (2012).
- Menga, A., Palmieri, E. M., et al. SLC25A26 overexpression impairs cell function via mtDNA hypermethylation and rewiring of methyl metabolism. FEBS J. 284 (6), 967-984 (2017).
- Frommer, M., McDonald, L. E., et al. A genomic sequencing protocol that yields a positive display of 5-methylcytosine residues in individual DNA strands. PNAS. 89 (5), 1827-1831 (1992).
- Byun, H. -M., Panni, T., et al. Effects of airborne pollutants on mitochondrial DNA Methylation. Part Fibre Toxicol. 10 (1), 1 (2013).
- Janssen, B. G., Byun, H. -M., Gyselaers, W., Lefebvre, W., Baccarelli, A. A., Nawrot, T. S. Placental mitochondrial methylation and exposure to airborne particulate matter in the early life environment: An ENVIRONAGE birth cohort study. Epigenetics. 10 (6), 536-544 (2015).
- Zheng, L. D., Linarelli, L. E., et al. Insulin resistance is associated with epigenetic and genetic regulation of mitochondrial DNA in obese humans. Clin Epigenetics. 7 (1), 1739 (2015).
- Byun, H. -M., Colicino, E., Trevisi, L., Fan, T., Christiani, D. C., Baccarelli, A. A. Effects of Air Pollution and Blood Mitochondrial DNA Methylation on Markers of Heart Rate Variability. J Am Heart Assoc. 5 (4), (2016).
- Bianchessi, V., Vinci, M. C., et al. Mitochondrion. MITOCH. 27 (C), 40-47 (2016).
- Wijst, M. G. P., van Tilburg, A. Y., Ruiters, M. H. J., Rots, M. G. Experimental mitochondria-targeted DNA methylation identifies GpC methylation, not CpG methylation, as potential regulator of mitochondrial gene expression. Sci Rep. 7 (1), 177 (2017).
- Wong, M., Gertz, B., Chestnut, B. A., Martin, L. J. Mitochondrial DNMT3A and DNA methylation in skeletal muscle and CNS of transgenic mouse models of ALS. Front cellr Neurosci. 7, 279 (2013).
- Dawid, I. B. 5-methylcytidylic acid: absence from mitochondrial DNA of frogs and HeLa cells. Science. 184 (4132), 80-81 (1974).
- Gama-Sosa, M. A. Levels and distribution of 5-methylcytosine in Chordate DNA. , 1-201 (1985).
- Hong, E. E., Okitsu, C. Y., Smith, A. D., Hsieh, C. -L. Regionally specific and genome-wide analyses conclusively demonstrate the absence of CpG methylation in human mitochondrial DNA. Mol Cell Biol. 33 (14), 2683-2690 (2013).
- Liu, B., Du, Q., et al. CpG methylation patterns ofhuman mitochondrial DNA. Nature Publishing Group. , 1-10 (2016).
- Donkin, I., Versteyhe, S., et al. Obesity and Bariatric Surgery Drive Epigenetic Variation of Spermatozoa in Humans. Cell Metab. 23 (2), 369-378 (2016).
- Li, L. C., Dahiya, R. MethPrimer: designing primers for methylation PCRs. Bioinformatics. 18 (11), 1427-1431 (2002).
- Tusnády, G. E., Simon, I., Váradi, A., Arányi, T. BiSearch: primer-design and search tool for PCR on bisulfite-treated genomes. Nucleic Acids Res. 33 (1), e9 (2005).
- Arányi, T., Váradi, A., Simon, I., Tusnády, G. E. The BiSearch web server. BMC bioinformatics. 7, 431 (2006).
- Krueger, F. Trim Galore!. , Available from: http://www.bioinformatics.babraham.ac.uk/projects/trim_galore/ (2018).
- Krueger, F., Andrews, S. R. Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics. 27 (11), 1571-1572 (2011).
- Li, H., Handsaker, B., et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 25 (16), 2078-2079 (2009).
- Ramos, A., Barbena, E., et al. Nuclear insertions of mitochondrial origin: Database updating and usefulness in cancer studies. Mitochondrion. 11 (6), 946-953 (2011).
- Warnecke, P. M., Stirzaker, C., Song, J., Grunau, C., Melki, J. R., Clark, S. J. Identification and resolution of artifacts in bisulfite sequencing. Methods (San Diego, Calif). 27 (2), 101-107 (2002).
