Summary
Nous présentons ici un protocole pour permettre une quantification précise de la méthylation de l’ADN (ADN mitochondrial) mitochondriale. Dans ce protocole, nous décrivons une digestion enzymatique de l’ADN avec BamHI couplé avec un pipeline d’analyse bio-informatique qui peut être utilisé pour éviter une surestimation des niveaux de méthylation de l’ADNmt causée par la structure secondaire de l’ADNmt.
Abstract
Quantification de la méthylation de l’ADN peut être réalisée en utilisant le séquençage de bisulfite, qui tire parti de la propriété de bisulfite de sodium pour convertir non méthylé cytosine en uracile, dans un contexte d’ADN simple brin. Le séquençage du bisulfite peut être ciblée (par PCR) ou effectuée sur l’ensemble du génome et fournit une quantification absolue de la méthylation de la cytosine à la résolution de base unique. Compte tenu de la nature distincte de l’ADN nucléaire et mitochondrial, notamment dans la structure secondaire, adaptations des méthodes de séquençage de bisulfite pour étudier la méthylation de la cytosine dans l’ADN mitochondrial se référera. Structure secondaire et tertiaire de l’ADNmt peut en effet conduire à bisulfite de séquençage des artefacts conduisant à des faux positifs en raison de la piètre accès incomplet dénaturation de bisulfite à l’ADN simple brin. Nous décrivons ici un protocole utilisant une digestion enzymatique de l’ADN avec BamHI couplé avec pipeline analyse bioinformatique pour permettre une quantification précise de la cytosine niveaux de méthylation de l’ADN mitochondrial. En outre, nous fournissons des lignes directrices pour concevoir les amorces de séquençage de bisulfite spécifiques d’ADN mitochondrial, afin d’éviter de cibler les segments nucléaire MiTochondrial indésirables (NUMTs) inséré dans le génome nucléaire.
Introduction
Le génome mitochondrial est une structure circulaire double brin de base environ 16,5 kilos (Ko), constituant un lourd et un brin de lumière. Le génome mitochondrial est présent en plusieurs exemplaires dans chaque cellule, maternellement hérité et encode les éléments essentiels de la chaîne respiratoire des complexes1. Semblables aux génomes bactériens et à la différence du génome nucléaire, le génome mitochondrial est organisé dans nombreuses structures secondaires et tertiaires, comme structures enroulé et surenroulé2, qui peuvent rendre l’accès difficile pendant le séquencement expériences,3.
Dans le noyau, la méthylation de l’ADN est une marque épigénétique étudiée qui joue un rôle dans de nombreux processus, notamment dans la régulation de l’expression des gènes. Dans les génomes de mammifères, la méthylation de l’ADN se produit principalement sur la position de l’anneau de pyrimidine de deoxycytidines, pour la plupart sur des dinucléotides CG (ou CpG) 5. Méthylation de la cytosine se trouve à 70 % de tous les CpG dans le génome des cellules somatiques et représente environ 1 % du total que des bases de l’ADN4. Méthylations de l’ADN ont également été recensés dans des contextes non-CpG, tels que le CpA, CpT et CpC et existent en différentes quantités dans l’ADN nucléaire, avec des valeurs jusqu'à 25 % de tous les cytosines méthylées dans les cellules souches embryonnaires5,6,7.
Alors que la méthylation de la cytosine du génome nucléaire est largement acceptée, l’existence de méthylation de l’ADN (ADN mitochondrial) mitochondriale est encore controversée. La première étude d’instruction ADNmt méthylation a été réalisée dans des cellules cultivées où la méthylation de l’ADN mitochondrial a été facilement détectée, mais à des niveaux inférieurs par rapport à l’ADN nucléaire8. Dans les cellules humaines et murines, méthylation de l’ADN mitochondrial a été également détectée à faibles doses (2-5 %). À l’aide d’essais en s’appuyant sur 5 méthylcytosine capture comme immunoprécipitation d’ADN méthylée (MeDIP) suivie d’une PCR quantitative, méthylation de l’ADN mitochondrial a aussi détectée dans diverses souris et les humains et les cellules lignes9,10, 11,,12. Utilisant des anticorps contre la 5-méthylcytosine dans une analyse ELISA ou spectrométrie de masse, des niveaux importants de méthylation de l’ADN ont été détectés des fractions mitochondriales purifiée13,14,15, 16. Toutefois, la plupart des essais dans les études susmentionnées a utilisé des techniques qui ne visaient pas à procurer une quantification absolue de méthylation de l’ADN à la résolution de base unique.
Analyse quantitative et résolutoire méthylation d’ADN peut être réalisée par une technique nommée « bisulfite de séquençage », qui tire parti de la propriété de bisulfite de sodium pour convertir non méthylé cytosine en uracile en simple brin ADN contexte17 . Moyen du séquençage de bisulfite, une constellation d’études a détecté la présence de la méthylation de la cytosine à différents niveaux. Méthylation de l’ADN mitochondrial dans la région de la boucle D, de la 12 s ou de la région de 16 s était facilement détectable à l’humain18,19,20,21,22,23 et souris24 tissus et cellules, cependant, avec une variabilité intrigante, de 1 à 20 % du totales cytosines dans toutes les études.
En comparaison avec ces nombreuses études, seules quelques études, y compris de notre groupe, ont contesté la présence de l’ADN mitochondrial méthylation3,25,26,27 ou conteste la pertinence biologique de niveaux très faibles d’ADNmt niveaux (inférieur à 2 %)28. Récemment, nous avons rapporté l’observation d’un artefact de bisulfite de séquençage potentiel au bisulfite ensemble mitochondrial séquençage3. On a démontré que la structure secondaire de l’ADN mitochondrial pourrait conduire à de faux positifs dans le séquençage de bisulfite, surestimant ainsi des niveaux de méthylation. Nous fournissons ici un protocole visant à prévenir un artefact de bisulfite de conversion de l’ADNmt. Ce protocole utilise une simple digestion enzymatique de l’ADN afin de perturber les structures secondaires de l’ADNmt et autorise l’accès complet au bisulfite selon un protocole de séquençage de bisulfite. En outre, nous fournissons un accompagnement pipeline bioinformatiques pour l’analyse du séquençage de bisulfite.
Subscription Required. Please recommend JoVE to your librarian.
Protocol
1. traitement de l’Enzyme de Restriction
- Linéariser ADNmt en traitant l’ADN humain total avec l’enzyme de restriction BamHI, qui coupe à position 14258 dans l’ADN mitochondrial humain.
Remarque : Pour souris ADN, utiliser l’enzyme de restriction BglII dans des conditions identiques comme ci-dessous. - Pour chaque échantillon où la méthylation de l’ADN mitochondrial doit être évaluée, préparer un tube de réaction de 0,2 mL et ajouter le mélange : 3 suivant μg ADN génomique (quantifié par fluorimétrie), 15 μl tampon de 3, 3 μL BamHI et l’eau jusqu'à 150 μL
- Placer les tubes dans un thermocycleur pendant 4 h à 37 ° C.
2. Conversion de bisulfite
- Convertir BamHI-a traité l’ADN à l’aide de bisulfite de sodium comme décrit ailleurs29.
- Utiliser 200-500 ng BamHI-a traité l’ADN de chaque réaction de conversion dans un volume total de 20 μL. La récupération de l’ADN est de 50 à 70 %.
- Éluer l’ADN Bisulfite-converti en tampon d’élution 10 μL.
- Quantifier l’ADN monocaténaire par fluorimétrie.
- En option : Stocker ADN bisulfite-converti à-20 ° C avant de passer le reste du protocole. Pour le stockage à long terme, stocker l’ADN bisulfite-converti à-80 ° C
3. conception du Bisulfite de séquençage des amorces
- Concevoir des amorces pour l’ordonnancement de bisulfite en utilisant les outils en ligne Methprimer30 et BiSearch31,32.
Remarque : Methprimer et BiSearch permettent de rechercher des séquences de liaison d’amorces sur les séquences génomiques où les cytosines sont converties en thymines, de la même façon au bisulfite après traitement. - D’amplification optimale et non biaisée, éviter dinucléotides dans les amorces et la taille de l’amplicon entre 100-300 bp de conception.
- S’assurer que les amorces sont spécifiques de l’ADN mitochondrial, à l’aide de la fonction de BiSearch Recherche de Primer. Insérer les amorces déjà conçu bisulfite converti, cocher la case de bisulfite et incluent un génome de référence et les paramètres de la PCR.
NOTE : Une liste de produits PCR possibles s’affichera. - Copiez le haut 1-5 séquences d’amorces et les coller dans un formulaire de commande pour la synthèse d’oligonucléotide.
NOTE : Une liste de l’homme et des séquences d’amorces d’ADNmt souris qui a été validé en laboratoire est affiché dans le tableau 1. - Des amorces de test à l’aide de bisulfite converti PCR Après électrophorèse sur gel d’agarose, tel que décrit à l’étape 4. Tester et optimiser toutes les paires d’amorces séparément et que la taille de l’amplicon correct apparaît sur le gel d’agarose.
Remarque : Si le multiplexage amorces sont nécessaires, ajoutez plusieurs amorces à une simple réaction de PCR et visualiser les produits PCR sur gel d’agarose ou, une méthode plus résolutoire comme l’électrophorèse sur gel de polyacrylamide, pour distinguer des produits similaires taille.
4. bisulfite converti PCR et Gel Extraction
- Pour amplifier les régions d’intérêt, effectuer la PCR à l’aide d’une polymérase de Taq de démarrage à chaud. En bref, mix 100 ng converti ADN, 500 amorces sens et anti-sens µM, 1 µL dNTP mix (10 µM de chaque dNTP) et polymérase à 7,5 U/réaction dans un volume total de 50 µL. PCR exécuter dans les conditions suivantes : 5 min à 95 ° C (60 s 94 ° C ; 60 s 55 ° C * ; 60 s 72 ° C)
35 cycles ; 10 min 72 ° C. Utiliser les amorces soit séparément ou multiplexées. - Facultatif : Lorsque le multiplexage des amorces, utilisation 200 ng converti ADN ainsi que les conditions suivantes de cyclisme : 5 min 95 ° C (60 s 94 ° C ; 90 s 55 ° C * ; 90 s 72 ° C) 35 cycles ; 10 min 72 ° C.
Remarque : La température peut varier selon la température de fusion d’amorces. - Produits PCR soumis au gel d’électrophorèse sur agarose à 2 % à 100 volts.
- Sous la douce lumière UV, extraire des produits PCR avec un scalpel et ajouter aux microtubes. Ne pas exposer les produits PCR à la lumière UV excessive afin d’éviter des dommages à l’ADN. Éviter les temps d’exposition plus de 30 s.
- Purifier les produits PCR à l’aide de gel de purification comme décrit plus haut3.
- Quantifier l’ADN purifié de l’étape 4,5 à l’aide de fluorométrie. Passez à l’étape 5 ou Banque d’échantillons à-20 ° C.
35 cycles ; 10 min 72 ° C. Utiliser les amorces soit séparément ou multiplexées.5. bisulfite séquençage bibliothèque préparation
- Effectuer la préparation de la bibliothèque des produits PCR du bisulfite-converti.
En option : multiplex produits PCR à ce stade. - Effectuer fin-réparation à l’aide de jusqu'à 100 ng de produit PCR. Ajouter 3 préparation enzymatique µL fin et tampon de réaction réparation fin µL 6.5 (x 10) à 55,5 produit de PCR µL dans un tube de 0,2 mL. Placer et incuber les tubes dans un thermocycleur pendant 30 min à 20 ° C, suivie de 30 min à 65 ° C.
- Ligaturer les adaptateurs de séquençage en mélangeant l’ADN de fin-réparé alliant 15 µL Blunt/TA ligase, 2,5 µL adaptateur et exhausteur de ligature 1 µL. Si vous utilisez < 100 produit de PCR ng comme entrée, diluer adaptateur 10 fois.
- Placez le mélange (étape 5.3) dans un thermocycleur et incuber pendant 15 min à 20 ° C. Mettre en pause le thermocycleur et ajouter 3 µl utilisateur Enzyme au tube. Mélanger et placer le tube dans le thermocycleur et incuber 15 min à 37 ° C.
- Taille sélectionnez l’ADN adaptateur-ligaturé à l’aide de perles Solid Phase réversible immobilisation (SPRI) avec des proportions variables de perles à l’ADN, selon la taille de produit PCR.
- Ajouter 13,5 µL H2O à ADN ligaturé à l’adaptateur d’étape 5.4 et 55 µL resuspendues SPRI Perles. Incuber à température ambiante pendant 5 min et place le tube sur un support magnétique. La solution est claire, transférer le surnageant dans un nouveau tube et jeter le tube contenant les perles.
Remarque : Le surnageant contient l’ADN adaptateur-ligaturé et grands fragments non désirés sont liés aux talons mis au rebut. - Ajouter 25 µL resuspendues SPRI perles au surnageant de point 5.6, mélanger, laisser incuber 5 min à température ambiante. Placer le tube sur le support magnétique et éliminer le surnageant lorsque la solution est claire.
Remarque : Le surnageant mis au rebut contient des séquences d’ADN non désirées, alors que l’ADN adaptateur-ligaturé est lié aux perles. - Bien que le tube sur le support magnétique, ajouter 200 µL 80 % d’éthanol (fraîchement préparé) pour laver les perles. Incuber 30 s et retirer et jeter le surnageant. Répétez cette étape pour un total de 2 lavages.
- Alors que le tube soit sur support magnétique avec un couvercle ouvert d’air secs perles pendant 5 min.
- Retirer le tube de l’aimant, éluer l’ADN dans 23 µL de tampon d’élution et mélanger. Incuber à température ambiante pendant 2 min.
- Placer le tube sur un support magnétique et lorsque la solution est claire, transférer 2 µL dans une plaque à 96 puits PCR pour l’amplification PCR-pré (étape 5.14).
- Transférer le µL environ 21 reste dans un nouveau tube de 0,2 mL pour l’amplification PCR (étape 5.17).
Remarque : Veillez à ne pas transférer les perles Perles peuvent inhiber les réactions enzymatiques des prochaines étapes. - Passez à l’étape 5.14 ou Banque d’échantillons à-20 ° C.
- Effectuer l’amplification pre-PCR par RT-qPCR pour estimer le nombre de cycles nécessaire pour amplifier l’ADN ligaturé à l’adaptateur. Par réaction, ajouter ce qui suit dans la PCR 96 puits contenant l’ADN 2 µL (étape 5.11) sur plaque : 0,4 amorce d’Index µL, 0,4 amorce PCR universelle µL, 7,2 µL H2O, 10 µL RT-qPCR master mix.
- Placer la plaque dans une machine PCR en temps réel et la plaque sous réserve des conditions suivantes : 30 s à 98 ° C (10 s 98 ° C ; 75 s 65 ° C) 20 cycles ; 5 min 65 ° C.
- Estimer le nombre de cycles d’amplification nécessaire pour chaque échantillon en soustrayant un Ct-valeur à la valeur Ct de la RT-qPCR de pré-amplification correspondant à la fin de la phase linéaire de la réaction PCR.
- Amplifier l’ADN ligaturé à l’adaptateur de l’étape 5.12 en mélangeant : 21 µL ADN, 25 µL PCR master mix, amorce de Index 1 µL, amorces de PCR universelle 1 µL et µL 2 H2O. Dans un thermocycleur, soumettre chaque échantillon aux conditions suivantes : 30 s à 98 ° C (10 s 98 ° C ; 75 s 65 ° C) [Ct calculée à l’étape 5.16] cycles ; 5 min 65 ° C.
Remarque : N’oubliez pas d’utiliser uniquement une amorce d’Index par échantillon permettant de multiplexage. L’index est inséré au cours de la PCR. - Purifier les perles SPRI de l’ADN amplifiés et éluer dans 22 µL de tampon d’élution.
- Contrôle qualité de la bibliothèque par la méthode de gel électrophorèse ou alternative. Recherchez les tailles de pointe bibliothèque correcte (taille de produit PCR et adaptateur).
- Estimer la taille moyenne de paires de bases du ou des produits bibliothèque par l’analyse du frottis : dans les paramètres globaux avancés, double-cliquez sur «tableau» selon l’analyse du frottis. Définir la zone de 100-1000 bp, puis cliquez sur OK. Sous la Table de la région, la région définie apparaît et bibliothèque de taille moyenne (bp) est calculée.
- Quantifier les bibliothèques par fluorimétrie. Passez à l’étape 6 ou bibliothèques de magasin à-20 ° C.
6. prochaine génération séquençage
- Calculer la concentration de nM de bibliothèques en utilisant la formule :
- Diluer les bibliothèques à la même nM de concentration par exemple 2 nM (choisissez 4 nM, 2 nM, 1 nM, soit 0,5 nM selon les concentrations de bibliothèque). Mettre en commun toutes les bibliothèques pour réaliser un bassin de 2 bibliothèques nM.
- Dénaturer et diluer les bibliothèques suivant le guide d’instrument de séquençage. En bref : dénaturer la piscine nM 2 en ajoutant 0,2 N NaOH et incuber à température ambiante pendant 5 min. diluer le pool dénaturé à une dilution finale de 22:00. Dénaturer et diluer une bibliothèque de contrôles disponible dans le commerce à une dilution finale de 12:05. Dans un nouveau tube, côtoient la piscine de 22:00 20 % bibliothèque de contrôles 12:05.
Remarque : Conversion de Bisulfite introduit une complexité très faible. 20 % contrôle bibliothèque spike-dans se traduira de meilleure qualité de séquençage. - Exécuter une séquence de bout jumelé bp 150 à l’aide d’un instrument de séquençage en profondeur.
7. calcul analyse - estimation méthylation niveaux
- Générer des fichiers .fastq (ici appelés sample1_R1.fastq.gz et sample1_R2.fastq.gz), création d’un index de bismark du génome entier d’intérêt (ici dans un dossier appelé bismarkIndex) et s’assurer que la séquence génomique de chrM est disponible en fasta format (ici dans un dossier appelé Arca).
- Pré-traiter lectures pour éliminer les adaptateurs et biais introduit par les amorces : utilisez l’outil Trim_galore33. Retirez deux nucléotides de l’extrémité 5' de la se de lit avant et arrière.
trim_galore--clip_R1 2--clip_R2 2 -o. /--jumelé trim1--sample1_R1.fastq.gz sample1_R2.fastq.gz - Aligner les lectures prétraitées à l’ensemble du génome à l’aide de Bismark34
Bismark--queue d’aronde--bam -o. /. / bismarkIndex / -1./sample1_R1_val_1.fq.gz -2./sample1_R2_val_2.fq.gz
Remarque : Si un autre dispositif d’alignement est utilisé, s’assurer que les lectures qui ne peuvent pas être particulièrement bien alignés sont ignorées. - Extrait lit cartographie du chromosome M, ici en utilisant Samtools35
samtools Trier -l 0 - O BAM -o sample1_sorted.bam
sample1_R1_val_1_bismark_bt2_pe.BAM
samtools index sample1_sorted.bam
samtools Découvre -u sample1_sorted.bam chrM | samtools tri - n - O BAM -o sample1_chrM.bam - Extraire les informations de méthylation
bismark_methylation_extractor -p--CX--cytosine_report--complet--genome_folder. /./sample1_chrM.bam chrM - Utiliser le .bedGraph, le .cov ou les fichiers CX_report.txt pour une analyse plus approfondie. Utilisez le fichier sample1_chrM.CX_report.txt pour la visualisation et l’analyse.
- Effectuer d’écrêtage plus loin dans l’étape de prétraitement, si plusieurs régions sont interrogées et un biais d’apprêt peut être visible dans les parcelles de M-partialité générées au cours du mappage. Reportez-vous à la documentation de Bismark pour plus d’informations.
8. computational Analysis - Test des différences
- Vérifiez que le fichier CX_report.txt contient 7 colonnes, chromosome (ici RMca), position, strand, nombre de lectures méthylés, le nombre de lectures non méthylés, C contexte et séquence environnant, pour tous les C sur le chromosome M.
- Comme le CX_report couvre chrM dans son intégralité, vérifier le sous-ensemble à l’ou les régions qui est amplifié.
- Utilisez les tests non paramétriques des différences entre les échantillons, tels que le test exact d’un pêcheurs pour individuel C ou un signe-test pour l’ensemble de la région.
NOTE : Quantification sera souvent proche de 0 % de méthylation, qui est la raison pourquoi l’hypothèse de normalité n’est pas adaptée. - De même, pour la visualisation, soit visualiser les quantiles ou calculer des intervalles de confiance de proportion binomiale.
Subscription Required. Please recommend JoVE to your librarian.
Representative Results
Deux étapes dans le présent protocole sont cruciales lors d’enquêtes sur la méthylation de l’ADN mitochondrial. 1) l’ouverture de la structure secondaire et 2) le design des amorces de spécifique de l’ADN mitochondriales.
Par digestion de l’ADN du génome humain avec l’enzyme de restriction BamHI (Figure 1), la structure de l’ADN mitochondriale sera réduite à la position de nucléotides 14 258 et la structure secondaire s’ouvrira.
Les altérations possibles de la conversion du bisulfite avec une structure secondaire de l’ADNmt intacte sont très probable à surestimer le taux d’unconversion d’ADNmt. Par séquençage de bisulfite, le taux d’unconversion d’ADN mitochondrial a été étudié, en non digérée et digéré l’ADN total des cellules de muscle squelettique humain, à 5 différentes régions du génome mitochondrial : le déplacement de boucle (D-Loop), le ribosome de phénylalanine et 12 s ARNt (ARNt - F + 12 s), 16 s ribosome (16), NADH déshydrogénase 5 (ND5) et gène de codage des ARNm du cytochrome b (CYTB) (Figure 2). Taux de unconversion varie entre 0 - 15,1 % dans toutes les régions étudiées pour ADNmt non digéré. Toutefois, lorsque l’ADN est digéré avant la conversion de bisulfite, le taux d’unconversion est tombé à un maximum de 1 % dans toutes les régions (Figure 3). Ainsi, illustrant l’importance de la digestion BamHI avant la conversion du bisulfite pour éviter une surestimation des niveaux de méthylation de l’ADN mitochondrial.
Contamination de l’ADN nucléaire est inévitable lorsque isoler le génome mitochondrial, et étant donné que presque l’ensemble du génome mitochondrial ont été répliquées dans le génome nucléaire, il a donné lieu à des insertions nucléaires du génome mitochondrial appelé nucléaire Segments mitochondriale (NUMTs)36. Pour s’assurer que les niveaux de méthylation de l’ADN analysés correspondent à l’ADN mitochondrial et pas NUMTs, il est important d’ultrapériphérique pour concevoir des amorces spécifiques pour le génome mitochondrial. La figure 4 montre comment vérifier la spécificité de l’apprêt et les affichages tableau 1 humains et les séquences d’amorces de souris ADNmt.
| Nom | Séquence | Position | Taille de l’amplicon | Température de recuit | Strand | ||
| Humaine | |||||||
| D-Loop | F: AAATCTATCACCCTATTAAC | ||||||
| R : GTGGAAATTTTTTGTTATGATGT | 6-298 | 292 | 55° C | Antisens | |||
| D-Loop | F: CATAACAAAAAATTTCCACCAAAC | ||||||
| R : GGGAAAATAATGTGTTAGTT | 279-458 | 179 | 55° C | Antisens | |||
| 12 S + TF | F: TTTATATAACTTACCTCCTC | ||||||
| R : GTGTTTGATGTTTGTTTTTTTTG | 577-765 | 188 | 55° C | Antisens | |||
| 466 | F: AATAAATTTATAGGTTTTTAAATTATTAAAT | ||||||
| R : TAACTAATAAAATCTTAACATATACTACTC | 2763-2873 | 110 | 53° C | Sens | |||
| CYTB | F: GGTATTATTTTTTTGTTTGTAATTATAGTA | ||||||
| R : CCTCAAATTCATTAAACTAAATCTATCC | 15091-15243 | 152 | 55° C | Sens | |||
| Souris | |||||||
| 426 | F: ACACATACAAACCTCCATAAAC | ||||||
| R : GAGGTATAATTTAGTTAAAT | 111-287 | 176 | 53° C | Antisens | |||
| 466 | F: AAATTCCAATTCTCCAAACATAC | ||||||
| R : TTTGGATTTTTTTTTTAGGT | 1749-1896 | 147 | 53° C | Antisens | |||
| CYTB | F: CTACAAAAACACCTAATAACAAAC | ||||||
| R : AGTGTATGGTTAAGAAAAGA | 14129-14307 | 178 | 55° C | Antisens | |||
| 466 | F: AGAGAAATAGAGTTATTTTATAAATAAG | ||||||
| R : CTAAAAACAAAATTTTAAATCTTAC | 2576-2727 | 151 | 55° C | Sens | |||
| ND5 | F: TTAAGTTAATTAGGTTTGATAATAGTGA | ||||||
| R : TAAAACCCTATTAAAAATAATATTCCTATA | 12659-12910 | 251 | 55° C | Sens | |||
| CYTB | F: TGGGTTTTTTTTAGGAGTTTGTTTA | ||||||
| R : CAATAAAAATACTCCAATATTTCAAATTTC | 14242-14504 | 262 | 55° C | Sens | |||
Tableau 1 : séquences d’amorces pour l’homme et la souris ADNmt. Notez que recuit températures de PCR amorces varient en raison des différences dans la température de fusion d’amorces.
Figure 1 : Électrophorèse de l’ADN digéré BamHI ou non digérée en. L’ADN génomique de muscle squelettique humain était non traitée ou traitée avec BamHI pendant 4h à 37° C et électrophorèse sur gel a été exécuté sur un gel d’agarose de 1,5 %. S’il vous plaît cliquez ici pour visionner une version agrandie de cette figure.
Figure 2 : Carte du génome mitochondrial humain. Régions du génome mitochondrial étudiées par séquençage du bisulfite et du site d’enzyme de restriction BamHI sont affichées. L’ADN mitochondrial contient 13 mRNA, 2 ARN ribosomal et 22 ARN de transfert. Abréviations : PH1: promoteur lourds-volet 1 ; PH2: Heavy-brin promouvoir 2 ; PL: promoteur de lumière-volet ; OH : Origine de réplication de lourds-strand ; OL: origine de réplication de lumière-strand. S’il vous plaît cliquez ici pour visionner une version agrandie de cette figure.
Figure 3 : BamHI digestion avant la conversion du bisulfite réduit le taux d’unconversion d’ADN mitochondrial dans les cellules de muscle squelettique humain. Par séquençage de bisulfite, pourcentage unconversion non digéré et digéré l’ADN mitochondrial ont été interrogés à cinq différentes régions du génome mitochondrial dans les cellules de muscle squelettique humain (N= 3). Le carré plein représente sites CpG et la place ouverte non-CpG sites. Les résultats sont présentés avec un intervalle mini-maxi et un test du signe a été utilisé pour vérifier les différences importantes unconversion. D-Loop (6-298) P= 1.83E-42 ; D-Loop (279-458) : P= 4.55E-13 ; ARNt - F + 12 ans : P= 8.38e-16 ; 16 s : P= 5.91E-06 ; ND5 : P= 1.70E-05 ; CYTB : P= 2.69E-10. D-Loop (6-298) comprend l’origine de réplication et ARNt - F + 12 s comprend promoteur brin lourd 2. Ces données ont été publiées ailleurs 3. S’il vous plaît cliquez ici pour visionner une version agrandie de cette figure.
Figure 4 : L’utilisation de Bisearch pour vérifier la spécificité des amorces vers le génome mitochondrial. A) séquences d’amorces Bisulfite-convertis sont ajoutés à la fonction BiSearch recherche Primer (http://bisearch.enzim.hu/?m=genompsearch). La boîte de bisulfite est cochée, et un génome de référence est sélectionné. B) exemple des produits PCR potentiels générés sur les chaînes de sens et antisens. S’il vous plaît cliquez ici pour visionner une version agrandie de cette figure.
Subscription Required. Please recommend JoVE to your librarian.
Discussion
Ici, nous fournissons un protocole de bisulfite de séquençage qui est spécifiquement conçu pour interroger la méthylation de l’ADN mitochondrial. Les différences avec le bisulfite de séquençage protocoles utilisés pour l’ADN génomique réside dans l’utilisation d’une étape de digestion préalable des enzymes de restriction et une analyse bioinformatique arrêt des faux positifs découlant de séquences NUMT.
Nous fournissons un protocole pour éviter les artefacts de bisulfite de séquençage lors d’enquêtes sur la méthylation de l’ADN mitochondrial. Artefacts de séquençage de bisulfite conduisant à des faux positifs ont été rapportés antérieurement37. Élimination insuffisante des protéines accessoires ADN de protéinase que k a été décrit37. Dénaturation incomplète ou partielle reannealing peut conduire à bouts droits de conversion incomplète, éventuellement en réduisant l’accès de bisulfite d’ADN simple brin37. L’artefact ultérieure pourrait être en jeu dans l’ADN mitochondrial où une structure secondaire porterait atteinte à l’accès à des cytosines. Au meilleur de nos connaissances, l’ajout d’une étape de digestion d’enzymes de restriction avant la conversion du bisulfite dans le contexte de l’estimation de méthylation de l’ADN mitochondrial a été tout d’abord signalé en 20163,28. Dans nos mains aussi, digestion préalable avec une enzyme de restriction single-coupeur a considérablement réduit détection des cytosines non converti3,28.
Ici, nous fournissons un pipeline de bioinformatique pour analyser les résultats du séquençage du bisulfite et de détecter des artefacts de séquençage de bisulfite dans le cadre du séquençage de génome entier mitochondrial bisulfite. Cette méthode provenance de notre observation qu’ADNmt non conversion taux était inversement corrélée avec le séquençage de profondeur3,28. Cette observation suggère que la structure secondaire et même supérieur d’ADNmt altère la fragmentation de l’ADN mitochondrial dans des régions spécifiques lors de la construction de bibliothèques de séquençage et pourrait également altérer le plein accès de bisulfite de sodium de l’ADNmt. Puisque la digestion avec un coupe-une enzyme libère les structures secondaires et tertiaires, nos résultats indiquent que la structure de l’ADN mitochondrial est une source d’artefact du bisulfite et valider les méthodes décrites ci-après pour la quantification de la méthylation de la cytosine dans l’ADN mitochondrial.
La contamination inévitable de l’ADN nucléaire en préparation mitochondriale constitue un défi. Cela conduit à la contamination dans les séquences NUMT dans lectures de séquençage, qui peut biaiser l’estimation de la méthylation de la cytosine dans l’ADN mitochondrial, même lorsque l’exécution ciblée et séquençage de génome pas ensemble mitochondrial bisulfite. Pour éviter ce biais, la conception des amorces de séquençage de bisulfite unique permet d’exclure contamination par méthylation au NUMTs. NUMTs couvrent certaines régions de l’être humain et le génome mitochondrial de la souris avec 100 % d’identité, et, par conséquent, il est impossible de concevoir des amorces uniques à certaines régions de l’ADN mitochondrial. Lorsque vous effectuez le séquençage génome mitochondrial ensemble bisulfite, séquençage longtemps lectures peuvent faciliter la discrimination entre l’ADN mitochondrial et NUMTs. Alignement de lectures pour le génome entier et n’utiliser que les lectures unique mappés sur le génome mitochondrial peuvent en outre assurent la détection de la méthylation de l’ADN mitochondrial.
Enfin, ce protocole peut être utilisé au-delà des génomes mitochondriaux, par exemple, dans d’autres cas de séquençage de bisulfite où la structure secondaire de l’ADN représente une source potentielle d’artefact. L’association entre niveaux de méthylation et séquençage de profondeur peut être étudiée afin d’évaluer la nécessité d’une digestion avec une enzyme de restriction avant séquençage du bisulfite de séquences d’ADN complexes.
Subscription Required. Please recommend JoVE to your librarian.
Disclosures
Les auteurs n’ont aucun intérêt financier concurrentes à déclarer.
Acknowledgments
Novo Nordisk Foundation Centre for Basic Research métabolique est un centre de recherche indépendant à l’Université de Copenhague, partiellement financé par un don sans restriction de la 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.
- 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).
