Bisulfite amplicon sequencing (BSAS) is a method for quantifying cytosine methylation in targeted genomic regions of interest. This method uses bisulfite conversion paired with PCR amplification of target regions prior to next-generation sequencing to produce absolute quantitation of DNA methylation at a base-specific level.
The role of epigenetic processes in the control of gene expression has been known for a number of years. DNA methylation at cytosine residues is of particular interest for epigenetic studies as it has been demonstrated to be both a long lasting and a dynamic regulator of gene expression. Efforts to examine epigenetic changes in health and disease have been hindered by the lack of high-throughput, quantitatively accurate methods. With the advent and popularization of next-generation sequencing (NGS) technologies, these tools are now being applied to epigenomics in addition to existing genomic and transcriptomic methodologies. For epigenetic investigations of cytosine methylation where regions of interest, such as specific gene promoters or CpG islands, have been identified and there is a need to examine significant numbers of samples with high quantitative accuracy, we have developed a method called Bisulfite Amplicon Sequencing (BSAS). This method combines bisulfite conversion with targeted amplification of regions of interest, transposome-mediated library construction and benchtop NGS. BSAS offers a rapid and efficient method for analysis of up to 10 kb of targeted regions in up to 96 samples at a time that can be performed by most research groups with basic molecular biology skills. The results provide absolute quantitation of cytosine methylation with base specificity. BSAS can be applied to any genomic region from any DNA source. This method is useful for hypothesis testing studies of target regions of interest as well as confirmation of regions identified in genome-wide methylation analyses such as whole genome bisulfite sequencing, reduced representation bisulfite sequencing, and methylated DNA immunoprecipitation sequencing.
Il a été plus d'un demi-siècle, depuis le premier rapport de modifications de l'ADN d'origine naturelle sous forme de méthylation de cytosine 1. La méthylation des cytosines est une base cytosine nucleotide modifié où la 5 ème carbone sur l'anneau de base a un groupe méthyle (5-mC), se produisant le plus fréquemment dans le dinucléotide CpG motif dans les génomes de mammifères. La présence fonctionnelle de 5-mC dans les promoteurs de gènes est généralement associée à la répression transcriptionnelle, tandis que l'absence est associée à une activité de transcription 2.
D'énormes progrès ont été réalisés dans notre compréhension du rôle de la méthylation de l'ADN dans le développement 3, transgénérationnelle propagation de profils épigénétiques 4, la pathogenèse du cancer 5,6, et un certain nombre d'autres domaines de recherche. Beaucoup de ces progrès sont venus au cours des dernières années que de nouvelles méthodes ont été développées pour le profil de méthylation de l'ADN 7. Avec l'avènement etpopularisation de séquençage de prochaine génération (NGS) techniques il ya maintenant un certain nombre de méthodes pour le profil méthylation de l'ADN dans un génome entier ou de grandes régions d'un génome huit. Ces approches ouvrent la possibilité d'études de découverte qui permettent de quantifier les modèles et les différences de méthylation de l'ADN méthylation ADN. En plus de ces approches hypothèse génératrices, il existe un besoin important pour cible / ADN hypothèse de test méthodes méthylation de quantification.
Pyroséquençage, PCR spécifique de la méthylation et séquençage Sanger directe de l'ADN converti au bisulfite ont été les méthodes les plus utilisées pour l'analyse des régions ciblées (ce est à dire, une région promotrice d'un seul gène ou d'un îlot CpG) 9-12. Toutes ces méthodes reposent sur la conversion au bisulfite, une réaction de désamination chimique où la cytosine non méthylées sont converties en de de l'uracile, la cytosine méthylés tandis que de 13,14 restent intacts. L'amplification par PCR de bisulfite-Résultats d'ADN converties en uracile de remplacement de la thymine avec 15 permettant la méthylation différentielle à lire comme une différence de base. Bien que très utile, les limites de ces méthodes comprennent faible précision quantitatives, faible longueur de lecture, faible capacité, et d'échantillon. Pour répondre à ces limitations nous avons développé une méthode que nous avons appelé bisulfite Amplicon séquençage (BSAS) 16. L'objectif de cette méthode est d'être capable d'analyser des régions génomiques cibles d'intérêt dans un grand nombre d'échantillons avec une précision quantitative élevé.
BSAS utilise des éléments de méthodes existantes (conversion au bisulfite et l'amplification de la PCR spécifique à la région) et les combine avec une simple construction de la prochaine génération bibliothèque (transposome médiation) 17 et paillasse NGS 18 (Figure 1). Cette approche fournit un procédé de débit quantitatif plus élevé et à examiner la méthylation de la cytosine dans ne importe quelle région d'intérêt. Nous présentons ici laméthode en détail et décrire des approches pour le dépannage et vérification de la qualité du processus BSAS.
BSAS permet de concentré, la quantification précise de méthylation d'ADN avec des capacités de débit élevées à la fois dans le nombre d'échantillons, le nombre de cibles. En outre, cette génération bibliothèque à l'aide de tagmentation transposome médiation nécessite l'ADN d'entrée moins, est plus rapide et contient moins d'étapes que cisaillement traditionnelle et les techniques de ligature adaptateur ultérieures. Ces améliorations par rapport aux méthodes existantes permettent précise le niveau base de l'analyse de méthylation d'ADN de toute cible d'intérêt. En outre, BSAS fournit une nouvelle méthode de confirmation de régions d'intérêt (régions de méthylation de manière différentielle (DMR), îlots CpG, régions réglementaires) qui ont été identifiés à l'aide de génomes entiers bisulfite 22 ou capturer 23 approches. En utilisant des approches de confirmation orthogonales et plus quantitativement méthodes précises améliore la rigueur analytique des études épigénomiques. La profondeur de lecture élevé atteint avec BSAS permet une quantification précise dans un intervalle de confiance étroit, rectsconsulting dans la quantification précise 16. En outre, la capacité d'exécuter de nombreux échantillons en une seule expérience peut augmenter la taille de l'échantillon pour les confirmations de méthodes du génome entier, ce qui augmentera la puissance statistique; une faiblesse de nombreuses études épigénétiques actuelles. Surtout, BSAS fournit une quantification spécifique de méthylation CpG-base, ce qui permet la génération d'hypothèses vérifiables pour la recherche en épigénétique. Les hypothèses comme l'augmentation de la méthylation à un élément de réponse spécifique provoque une diminution de la liaison d'un facteur de transcription spécifique. BSAS, combinées avec les données d'expression d'ARNm couplés, fournit un procédé d'obtention à la fois la régulation épigénétique et l'expression de l'ARNm au sein d'une seule pièce de tissu ou population de cellules. Mesures paires de régulation et l'expression augmentent également la rigueur scientifique et l'impact des résultats.
Les étapes critiques au sein du protocole BSAS comprennent la conception primaire, l'optimisation de l'amplicon, et la génération de bibliothèque /QC. Biais PCR est introduit si amorces sont pas bien conçus et amplifient l'efficacité à différents 8. L'utilisation de normes de méthylation, comme 100% et 0% cytosine méthylés normes générés par voie enzymatique, peut être utilisé pour déterminer le degré, le cas échéant, de partialité méthylation de quantification, et est un contrôle méthodologique appropriée lors de la quantification méthylation 16. De même, il faut prendre soin lors de l'optimisation des réactions de PCR pour générer un seul amplicon de haute qualité. Si aucun amplicon est présent en utilisant les lignes directrices proposées décrites dans le protocole, des mesures d'optimisation comprennent l'augmentation du nombre de cycles PCR ou le montant d'entrée bsDNA. Se il ya des produits de PCR multiples, y compris des dimères d'amorces, des mesures d'optimisation comprennent quantité décroissante d'entrée bsDNA, diminuant les concentrations d'entrée molaires d'amorces PCR, la hausse des températures de recuit, ou en diminuant le nombre de cycles PCR. L'un ou une combinaison de ces procédures d'optimisation donne des résultats suffisants pour générer une seule higamplicon h-qualité. Alternativement, la refonte amorces est une bonne approche pour jeux d'amorces qui ne produisent pas un seul amplicon. Dans les régions où la refonte ne est pas approprié ou possible, dégénérer jeux d'amorces dont les deux peuvent travailler de la cytosine et la thymine de sites CpG dans amorces sens et de l'adénine et la guanine de sur les sites CpG dans amorces inverses. Cependant, ces jeux d'amorces doivent poursuite de l'optimisation pour se assurer qu'aucun biais PCR est en cours, un processus en dehors du champ d'application de ce protocole. bibliothèques de qualité sont extrêmement important pour le succès. Assurez-vous que QC mesures sont effectuées pour déterminer la taille, la qualité et la quantité des bibliothèques avant séquençage. Les réactions de séquençage sont sujettes à l'échec avec la participation des bibliothèques de faible qualité. Biais dans la méthylation quantification peut être atténué en veillant à ce que les fréquences de méthylation sont présents à la fois avant et arrière lit séquençage. De plus, les scores de qualité de l'alignement des bases de sites CpG devraient être ≥Q30 pour la haute confidence dans la quantification. Alors que d'autres ont montré précédemment peu ou pas de biais dans les réactions de tagmentation de bisulfite ADN converti 19, assurant que les mesures décrites ci-dessus permet de confirmation de faible quantification biais de méthylation.
BSAS mesure actuellement méthylation des sites CpG, cependant, de futures extensions de cette méthode permettra mesures paires de CpG 5-HMC avec l'ajout d'étapes expérimentales distinction entre 5-mC et 5-HMC telles que l'introduction perruthénate de potassium avant la conversion au bisulfite 24. Cela permettra d'accroître l'impact de l'utilité BSAS en permettant pour la méthylation jumelé, deux 5-mC et 5-HMC, et des données d'expression afin de déterminer les rôles régulation de l'expression des gènes de méthylation de l'ADN. Préparation de bibliothèque Transposome médiée par des amplicons PCR ne entraîne une diminution de la profondeur de séquençage aux extrémités des régions amplifiées. Par conséquent, les régions d'un grand intérêt devraient être conçus pour résider dans le milieu de amplicons. Commentjamais, parce que BSAS génère séquençage lire profondeurs de> 1000 X, la précision quantitative de 5-mC mesures près des extrémités des régions amplifiées reste élevé.
The authors have nothing to disclose.
The authors thank Colleen Van Kirk for mouse retina and cerebellum samples, and Peter Gregory for figure generation. The authors also thank Dr. Allison Gillaspy and the Laboratory for Molecular Biology and Cytometry Research core for access to the MiSeq. This work was supported by NIH grants DA029405, EY021716, and AG026607 and the Donald W. Reynolds Foundation.
Name | Company | Cat # | Comments |
Agilent 2100 Bioanlyzer | Agilent | G29393AA | |
Agilent RNA 6000 Nano kit | Agilent | 5067-1511 | |
Agilent High Sensitivity DNA kit | Agilent | 5067-4626 | |
Typhoon 9200 | Amersham/GE Healthcare Life Sciences | ||
Agencourt AMPure XP – PCR Purification | Beckman Coulter | A63880 | 5 mL |
PowerPac Basic Power Supply | Bio Rad | 164-5050 | |
twin.tech PCR plate 96 | Eppendorf | 951020362 | semi-skirted |
Heatsealing Film | Eppendorf | 0030127838 | |
Heatsealing Foil | Eppendorf | 0030127854 | |
Heat Sealer | Eppendorf | 5390000.024 | |
0.1-10ul Tips | Eppendorf | 0030073.002 | |
2-200ul Tips | Eppendorf | 0030073.045 | |
50-1000ul Tips | Eppendorf | 0030073.100 | |
MixMate | Eppendorf | 5353000.014 | |
1.5 mL Microcentrifuge Tubes | Eppendorf | 0030121.023 | |
Centrifuge 5424 | Eppendorf | 5424000.010 | |
2.0 mL Microcentrifuge Tube | Eppendorf | 022363352 | |
Nextera XT DNA Sample Preparation Kit | Illumina | FC-131-1024 | 24 rxn |
Nextera XT Index Kit | Illumina | FC-131-1001 | 24 indexes |
MiSeq Desktop Sequencer | Illumina | ||
MiSeq Reagent Kit v2 | Illumina | MS-102-2002 | 300cycle |
KAPA SYBR FAST Universal qPCR kit | KAPA Biosystems | KK4824 | |
Quant-iT PicoGreen dsDNA Assay kit | Life Technologies | P7589 | |
Quant-iT RiboGreen RNA Assay kit | Life Technologies | R11490 | |
ProFlex 96-well PCR system | Life Technologies | 4484075 | |
Novex 6% TBE DNA gel | Life Technologies | EC6261BOX | 10-well |
Novex TBE Running Buffer | Life Technologies | LC6675 | 5X (1 L) |
Novex Hi-Density TBE Sample Buffer | Life Technologies | LC6678 | 5X (10 mL) |
1 Kb Plus DNA Ladder | Life Technologies | 10787-018 | |
Xcell SureLock Mini-Cell | Life Technologies | EI0001 | |
UltraPure 10mg/mL Ethidium Bromide | Life Technologies | 15585-011 | |
7900HT Fast Real-Time PCR system | Life Technologies | 4329001 | 384well block |
DynaMag-96 Side Skirted | Life Technologies | 12027 | |
SpectroMax M2 Multi-Mode Microplate Reader | Molecular Devices/VWR | 89429-532 | |
AllPrep DNA/RNA Mini Kit | Qiagen | 80204 | |
Stainless Steel Beads | Qiagen | 69989 | 5 mm |
CLC Genomics Workbench | Qiagen | Software for methylation pipeline | |
QIAQuick PCR Purification kit | Qiagen | 28104 | 50 rxn |
TissueLyser II | Retsch/Qiagen | 85300 | |
Nuclease-Free Water | Sigma | W4502 | |
TRIS hydrochloride | Sigma | PHG0002-100G | |
Tween-20 | Sigma | P9416-50ML | |
NanoDrop 2000 | Thermo Scientific | ||
384-well Full-Skirted Plate, Standard | Thermo Scientific | AB-1384 | |
Absolute qPCR Seal | Thermo Scientific | AB-1170 | |
EZ DNA Methylation-Lightning kit | Zymo Research | D5030 | 50 rxn |
ZymoTaq DNA Polymerase | Zymo Research | E2001 | 50 rxn |