We beschrijven hier, het protocol en de uitvoering van Methyl-Seq, een epigenomic-platform, met behulp van een rat-model om epigenetische veranderingen die zijn gekoppeld aan chronische stress blootstelling te identificeren. Resultaten tonen aan dat de rat Methyl-Seq platform is voor het opsporen van methylation van verschillen die uit stress blootstelling bij ratten voortvloeien van.
Zoals genoom van een breder scala dieren beschikbaar komen, is er een toenemende behoefte aan tools die dynamische Epigenetische veranderingen in deze dierlijke modellen kunt vastleggen. De rat is een bepaald model dier waar een epigenetische instrument veel farmacologische en behavioral studies inzichtelijke mechanistische informatie te verschaffen kan aanvullen. Te dien einde aangepast we de SureSelect Target Capture System (hierna: Methyl-Seq) voor de rat, die DNA methylering niveaus over het genoom van de rat beoordelen kan. Het ontwerp van de rat gericht initiatiefnemers, CpG eilanden eiland oevers en GC-rijke regio’s uit alle RefSeq genen.
Ter uitvoering van het platform op een rat-experiment, werden mannelijke Sprague-Dawley ratten blootgesteld aan chronische variabele spanning voor 3 weken, waarna de bloedmonsters werden verzameld voor genomic DNA-extractie. Methyl-Seq bibliotheken werden opgebouwd uit de rat DNA-monsters door scheren, adapter afbinding doel verrijking, Bisulfiet Omzetting en multiplex. Bibliotheken werden sequenced op een volgende-generatie sequencing-platform en de gesequenceerd leest werden geanalyseerd om te achterhalen van DMRs tussen DNA van beklemtoonde en onbeklemtoonde ratten. Top kandidaat DMRs werden onafhankelijk gevalideerd door bisulfiet pyrosequencing te bevestigen de robuustheid van het platform.
Resultaten tonen aan dat de rat Methyl-Seq-platform een nuttig epigenetische hulpmiddel dat vangen methylering veranderingen bij blootstelling is kan te benadrukken.
Voorschotten in high-throughput sequencing hebben geleid tot een schat aan de genomic opeenvolgingen voor zowel model als niet-modelorganismen. De beschikbaarheid van dergelijke sequenties heeft onderzoek in genetica, vergelijkende genomica en transcriptomics vergemakkelijkt. Bijvoorbeeld, beschikbaar genomic opeenvolgingen zijn zeer nuttig voor het uitlijnen van sequencing gegevens van ChIP-Seq experimenten die DNA op basis van zijn associatie met histone modificaties1of bisulfiet sequencing verrijken, welke maatregelen methylation van DNA door opsporen van uracil gevormd uit Bisulfiet Omzetting van unmethylated cytosines2. Echter zijn er vertragingen bij de tenuitvoerlegging van epigenomic-platforms die beschikbaar genomic sequencing gegevens in hun ontwerp als gevolg van een gebrek aan geannoteerde gegevens van soortspecifieke regulerende sequenties die invloed op de genfunctie uitoefenen kunnen opnemen.
In het bijzonder is methylation van DNA een van de meest bestudeerde epigenetische aanpassingen op DNA dat kan hefboomwerking beschikbare genomic gegevens voor het bouwen van een methylomic platform. Een voorbeeld hiervan is een array gebaseerde platform voor de menselijke methylome3, die wijd verbeid in verschillende disciplines van oncologie psychiatrie4,5 gebruikt is. Soortgelijke platforms voor niet-menselijke dierlijke modellen zijn helaas schaars is, aangezien er vrijwel geen gebruikte platforms die hebben geprofiteerd van de genomic volgorde in hun oorspronkelijke ontwerp.
Een gemeenschappelijke methode voor de beoordeling van het methylomic landschap van niet-menselijke dierlijke modellen is verminderde vertegenwoordiging bisulfiet sequencing (RRBS)6. Deze aanpak overwint de kosten van de gehele-genoom bisulfiet sequencing waarmee, terwijl het verstrekken van een uitgebreide methylomic landschap, lagere Lees-diepte dekking als gevolg van de kosten en de beperkte functionele informatie in grote gebieden van de gen-armen van het genoom2 . RRBS houdt de samenvatting van de beperking en grootte-selectie van genomic DNA te verrijken voor zeer GC-rijke reeksen zoals CpG eilanden die vaak worden gevonden in de buurt van gene initiatiefnemers en dacht dat een rol te spelen in gen verordening7. Terwijl de methode RRBS is gebruikt in een aantal belangrijke studies, is haar afhankelijkheid van de enzymen van de beperking niet zonder opvallende uitdagingen en beperkingen. Bijvoorbeeld, is verrijking van GC-rijke sequenties in RRBS volledig afhankelijk van de aanwezigheid van specifieke opeenvolgingen erkend door de restrictie-enzym en de daaropvolgende maat keuze door elektroforese. Dit betekent dat elke genomic gebieden die geen deze beperkingsplaatsen bevatten zijn uitgesloten tijdens de selectie van de grootte. Kruis-soorten vergelijkingen zijn ook uitdagend tenzij de dezelfde beperkingsplaatsen in de dezelfde loci onder de verschillende soorten aanwezig zijn.
Een benadering voor het overwinnen van de beperkingen van RRBS is het gebruik van een methode van de verrijking die van de gepubliceerde genomic sequentie in het ontwerp van het platform profiteert. Het arraygebaseerd menselijke platform gebruikt primer sondes ontworpen tegen specifieke CpGs voor allele-specifieke (CG vs. TG na bisulfiet omzetting) doel gloeien en primer extensie. Haar ontwerp weerspiegelt niet alleen de beschikbare menselijke genomic volgorde, maar experimenteel geverifieerd regelgevende regio’s verworven van meerdere lijnen van onderzoek, zoals ENCODE en ENSEMBL8. Ondanks het ruime gebruik ervan in menselijke methylomic onderzoeken bestaat een soortgelijk platform niet voor model dieren. Daarnaast plaatst de indeling arraygebaseerd aanzienlijke beperking op de oppervlakte die beschikbaar is voor de plaatsing van de sonde. In de afgelopen jaren, hebben inspanningen geleverd om te combineren van de doel-specificiteit door vangen sonde ontwerp en de functie van de hoge-doorvoer van volgende-generatie sequencing geboden. Een dergelijke inspanning heeft geleid tot de sequentiebepaling gebaseerde target verrijking-systeem voor de muis-genoom (muis Methyl-Seq), die werd gebruikt als aanduiding van de hersenen-specifieke of glucocorticoide-geïnduceerde verschillen in methylering9,10. Soortgelijke platforms voor andere model en niet-model dieren zijn nodig om epigenomic onderzoek in deze dieren.
Hier tonen we de uitvoering van dit nieuwe platform te analyseren van de methylomic op de rat. De rat heeft gediend als een belangrijke diermodel in farmacologie, metabolisme, neuroendocrinologie en gedrag. Bijvoorbeeld, is er een toenemende behoefte te begrijpen van de onderliggende mechanismen die aanleiding tot drug toxiciteit, overgewicht, stress-respons of drugsverslaving geven. Een high-throughput platform geschikt voor het vastleggen van de wijzigingen van de methylomic die zijn gekoppeld aan deze voorwaarden zou het vergroten van ons begrip van de mechanismen. Sinds het genoom van de rat ontbreekt nog steeds een aantekening voor regelgevende gebieden, we opgenomen niet-redundante initiatiefnemers, CpG eilanden, eiland oevers11, en eerder aangeduid GC-rijke sequenties in de rat Methyl-Seq platform12.
Om te beoordelen succesvol ontwerp en tenuitvoerlegging van het platform SureSelect Target verrijking (generiek Methyl-Seq genoemd) voor de rat genoom, we een rat model van chronische variabele spanning (CVS)13 te identificeren differentieel methylated gebruikt regio’s tussen onbeklemtoonde en gestresste dieren. Onze platformontwerp, protocol en uitvoering kunnen nuttig zijn voor onderzoekers die willen een uitgebreide en onbevooroordeelde epigenetische onderzoek in te stellen op een organisme waarvan genomic volgorde is al beschikbaar, maar blijft slecht geannoteerde.
In deze studie, we ontworpen en geïmplementeerd de Methyl-Seq-platform voor het genoom van de rat. Door aan te tonen het nut ervan met een rat-model van stress, toonden we aan dat de pijpleiding van experimentele en analytische differentieel gemethyleerd regio’s tussen twee groepen van de vergelijking kan bieden.
Om te garanderen een succesvolle implementatie van het platform, moeten verschillende kritische stappen in acht worden genomen. Eerste, initiële DNA kwaliteit en kwantiteit heeft een aanzienlijke invloed op de kwaliteit en kwantiteit van de definitieve Methyl-Seq-bibliotheek. We gebruikten een Fluorimeter, in plaats van een spectrofotometer, om ervoor te zorgen dat onze DNA meting weerspiegeld de hoeveelheid double-stranded DNA aanwezig. De bioanalyzer werd gebruikt voor het meten van de moleculaire grootte en hoeveelheid van DNA na scheren en na adapter afbinding. Controleren of de moleculaire grootte “shift” tussen deze stappen is cruciaal voor het bevestigen van de aanwezigheid van adapters aan het einde van elk fragment van DNA die adapter-gemedieerde PCR in de opeenvolgende stappen zal ondergaan. De hoeveelheid DNA resterende aan het eind van de adapter afbinding stap is ook belangrijk, sinds ten minste 100 ng van het product van de bibliotheek is die bij deze stap nodig is om voldoende hoeveelheid is beschikbaar na de doelgroep verrijking en Bisulfiet Omzetting stappen. Een definitieve ultragevoelige meting werd uitgevoerd op de geconstrueerde Methyl-Seq-bibliotheek, zodat de bibliotheek correct kan worden verdund voor latere clustering op de volgende-generatie-sequencer. Ten slotte was bisulfiet pyrosequencing werkzaam als een hoogst-kwantitatieve, onafhankelijke methode te beoordelen van de juistheid van de analytische pijpleiding. De uiteindelijke validatie met behulp van de oorspronkelijke monsters en replicatie met behulp van extra dieren zijn cruciale stappen om ervoor te zorgen dat het experiment biologisch belangrijke veranderingen in DNA-methylering kan detecteren.
Wij omvatten ook verschillende richtsnoeren in geval van afwijking van het protocol of als problemen ondervindt. Eerst, is het mogelijk om te verliezen teveel DNA tijdens einde-reparatie, adapter afbinding of magnetische kraal zuivering stappen. Anderzijds beginnen bedragen van DNA zou kleine (< 200 ng) als gevolg van de beperkte weefsel/DNA beschikbaarheid of de implementatie van verschillende verrijking methoden zoals het sorteren van de cel van de fluorescentie geactiveerd. Verhoging van het aantal cyclus tijdens de twee bibliotheek amplificatie stappen mogelijk kunnen compenseren voor het overmatig verlies van DNA of lage DNA bedrag in het gehele protocol bibliotheek bouw beginnen. Echter worden niet meer dan een extra 2-3 cycli aanbevolen, zoals overmatige sjabloon versterking dreigt te leiden tot een toename van het aantal dubbele leest wordt sequenced. Deze dubbele vermeldingen zijn uitgesloten tijdens de uitlijning stap om te voorkomen dat vooringenomenheid in percentage methylering berekeningen. Ten tweede, als de gemiddelde grootte van de DNA niet door meer dan 30 bps toeneemt, controleren om ervoor te zorgen dat de reagentia zijn nieuwe, als de polymerase van DNA van de T4, Klenow, en/of T4 ligase oude. Verkrijgbare vervanging reagentia kunnen worden gebruikt.
Bovendien is het mogelijk dat de voorspelde DMRs pyrosequencing, waar DNA methylering verschillen bestaan niet of zijn aanzienlijk minder dan voorspeld door de analyse niet kunnen valideren. Slechte validatie van kandidaat-regio’s is een probleem te vaak voor veel genoom-brede analyses, zoals wanneer pyrosequencing resultaten differentiële methylering niet bevestigen of de effect grootte is veel kleiner dan die voorspeld door de analyse. BSmooth is een analytische pakket dat “glad” de methylering niveaus in een venster van meerdere CpGs. Voor de huidige experiment betrokken BSmooth een DMR waarvan methylering niveaus zijn gevalideerd door bisulfiet pyrosequencing. Er zal waarschijnlijk wel verschillen vastgesteld tussen methylering niveaus voorspeld door BSmooth en geverifieerd door pyrosequencing. De verschillen voortkomen uit de smoothing functie die de gemiddelde methylering waarden in alle van de CpGs binnen een DMR schat, met inbegrip van opeenvolgende CpGs die door meer dan 50% in methylation van DNA afwijken of CpGs waarvan methylation van de waarden werden uitgesloten wegens sub drempel Lees diepte. R-pakketten zoals MethylKit24 kunnen worden gebruikt voor identificatie van kleinere vensters of zelfs één CpGs waarvan methylering niveaus sterk met die gevalideerd door pyrosequencing correleren en CpGs. Uitvoering van de verschillende pakketten en testen hun voorspelde regio’s of de CpGs van differentiële methylering door pyrosequencing zorgt robuustheid van gegevens. Als alternatief, originele Methyl-Seq bibliotheken kunnen worden resequenced en toegevoegd aan de gelezen bestanden toe Lees diepte. Aangezien de gehaltebepaling methylering semi-kwantitatieve en gedicteerde door het aantal leest [(# of CpGs) / (# van TpGs + CpGs)], verhoging van de Lees diepte voor een bepaalde CpG de nauwkeurigheid van de waarde van percentage methylering verhoogt. In deze studie vonden wij alleen CpGs waarvan methylation van de waarden werden bepaald door ten minste tien leest en een algehele Lees dekking van 19 x voor elke CpG bereikt.
De rat Methyl-Seq platform is niet zonder haar beperkingen. Hoewel het rendabeler dan bisulfiet geheel-genoom sequencing, is het aanzienlijk duurder dan andere methoden. Niettemin was de grootste deel van de kosten voor de aankoop van rijstroken op de sequencer en niet voor het systeem vastleggen. Afhankelijk van de diepte van Lees nodig met kruis-weefsel vergelijkingen vereisen minder als gevolg van grote (25-70%) verschillen12 in methylation van DNA, kan de kosten worden verminderd door multiplexing meer monsters per rijstrook en met behulp van een hogere capaciteit platform. Ook kost de monsterverwerking meer tijd dan andere methoden. Terwijl gelijkaardig aan andere pulldown benaderingen die gebruikmaken van de volgende generatie sequencing, toevoegen de toegevoegde Bisulfiet Omzetting en zuivering stappen aan de werklast. Globaal, is het Methyl-Seq-platform is een kosteneffectief alternatief voor geheel-genoom sequentiebepaling en cijferreeks basenpaar op meer dan 2,3 miljoen CpGs, dat is aanzienlijk meer dan de vehiculumcontrolegroep door microarray gebaseerde platformen. Tot op heden, de verkrijgbare mens en muis zijn Methyl-Seq platformen gebruikt om het document alcohol-afhankelijke veranderingen in de makaak hersenen25,26, neurologische genen in de muis hersenen9, en de blood-brain doelstellingen van glucocorticoïden10. Verder is maakt de mogelijkheid om zich te richten op specifieke regio’s, ongeacht de volgorde erkenning door enzymen van de beperking het een ideaal platform voor cross-soorten vergelijkingen. Voor deze studie ontwierpen we de Methyl-Seq-platform voor de rat, waarvoor vele farmacologische, metabole en gedrags experimenten worden uitgevoerd zonder het voordeel van een genoom-brede methylomic tool. Onze gegevens tonen aan dat het kan worden gebruikt voor het detecteren van DMRs in een rat-model van stress en andere fysiologische parameters zoals algemene plasma CORT niveaus gecorreleerd.
De Methyl-Seq-platform is ideaal voor epigenetische experimenten in dieren met gesequenceerd genomen die wellicht niet genoeg experimenteel bewijsmateriaal documenteren van de regelgevende gebieden. Wanneer dergelijke regio’s ter beschikking worden gesteld, kunnen extra gebieden worden speciaal ontworpen en gekoppeld aan de huidige versie. Verder, het platform is ideaal voor vergelijkende genomica, aangezien de verrijking van het doel niet wordt belemmerd door restrictie-enzym erkenning. Bijvoorbeeld, kan de regio van de promotor van een gen van belang worden vastgelegd ongeacht of het een specifieke beperkingsplaats havens. Ook kunnen eventuele regelgevende gebieden, zoals die welke in de muis of mens, die zijn bewaard in het genoom van belang worden vastgelegd.
The authors have nothing to disclose.
Deze studie werd gefinancierd door de NIH grant MH101392 (RSL) en steun van de volgende awards en stichtingen: een NARSAD Young Investigator Award, Margaret Ann prijs onderzoeker Fonds, het James Wah Mood Disorders geleerde fonds via de Charles T. Bauer Foundation, Baker Foundation en de Stichting van het Project-Match (RSL).
Radioimmuno assay (RIA) | MP Biomedicals | 7120126 | Corticosterone, 125I labeled |
Master Pure DNA Purification Kit | Epicentre/Illumina | MC85200 | |
Thermal-LOK 2-Position Dry Heat Bath | USA Scientific | 2510-1102 | Used with 1.5 mL tubes |
Vortex Genie 2 | Fisher | 12-812 | Vortex Mixer |
Ethyl alcohol, Pure | Sigma-Aldrich | E7023 | 100% Ethanol, molecular grade |
Centrifuge 5424 R | Eppendorf | – | Must be capable of 20000 x g |
Qubit 2.0 | ThermoFisher Scientific | Q32866 | Fluorometer |
Qubit dsDNA BR Assay Kit | ThermoFisher Scientific | Q32850 | |
Qubit dsDNA HS Assay Kit | ThermoFisher Scientific | Q32851 | High sensitivity DNA detection reagents |
Qubit Assay Tubes | ThermoFisher Scientific | Q32856 | |
SureSelectXT Rat Methyl-Seq Reagent Kit | Agilent Technologies | G9651A | Reagents for preparing the Methyl-Seq library |
SureSelect Rat Methyl-Seq Capture Library | Agilent Technologies | 931143 | RNA baits for enrichment of rat targets |
IDTE, pH 8.0 | IDT DNA | 11-05-01-09 | 10 mM TE, 0.1 mM EDTA |
DNA LoBind Tube 1.5 mL | Eppendorf | 22431021 | |
Covaris E-series or S-series | Covaris | – | Isothermal sonicator |
microTUBE AFA Fiber Pre-Slit Snap-Cap 6x16mm (25) | Covaris | 520045 | |
Water, Ultra Pure (Molecular Biology Grade) | Quality Biological | 351-029-721 | |
Veriti 96 Well-Thermal Cycler | Applied Biosystems | 4375786 | |
AMPure XP Beads | Beckman Coulter | A63880 | DNA-Binding magnetic beads |
96S Super Magnet | ALPAQUA | A001322 | Magnetic plate for purification steps |
2200 TapeStation | Agilent Technologies | G2965AA | Electrophresis-based bioanalyzer |
D1000 ScreenTape | Agilent Technologies | 5067-5582 | |
D1000 ScreenTape High Sensitivity | Agilent Technologies | 5067-5584 | |
D1000 Reagents | Agilent Technologies | 5067-5583 | |
D1000 Reagents High Sensitivity | Agilent Technologies | 5067-5585 | |
DNA110 SpeedVac | ThermoFisher Scientific | – | Vacuum Concentrator |
Dynabeads MyOne Streptavidin T1 magnetic beads | Invitrogen | 65601 | Streptavidin magnetic beads |
Labquake Tube Rotator | ThermoFisher Scientific | 415110Q | Nutator Mixer is also acceptable |
EZ DNA Methylation-Gold Kit | Zymo Research | D5006 | Bisulfite conversion kit. Contains Binding, Wash, Desulphonation, and Elution buffers |
Illumina Hi-Seq 2500 | Illumina | – | Next-generation sequencing machine |
PCR and Pyrosequencing Primers | IDT DNA | Variable | |
Taq DNA Polymerase with ThermoPol Buffer – 2,000 units | New England BioLabs | M0267L | |
Deoxynucleotide (dNTP) Solution Set | New England BioLabs | N0446S | |
Pyromark MD96 | QIAGEN | – | Pyrosequencing machine |
Ethyl Alcohol 200 Proof | Pharmco-Aaper | 111000200 | 70% Ethanol solution |
Sodium Hydroxide Pellets | Sigma-Aldrich | 221465 | 0.2 M NaOH denature buffer solution |
Tris (Base) from J.T. Baker | Fisher Scientific | 02-004-508 | 10 mM Tris Acetate Buffer wash buffer solution |
PyroMark Gold Q96 Reagents (50×96) | QIAGEN | 972807 | Reagents required for pyrosequencing |
PyroMark Annealing Buffer | QIAGEN | 979009 | |
PyroMark Binding Buffer (200 ml) | QIAGEN | 979006 | |
Streptavidin Sepharose High Performance Beads | GE Healthcare | 17-5113-01 | Streptavidin-coated sepharose beads |
PyroMark Q96 HS Plate | QIAGEN | 979101 | Pyrosequencing assay plate |
Eppendorf Thermomixer R | Fisher Scientific | 05-400-205 | Plate mixer. 96-well block sold separately (cat. No 05-400-207) |
SureDesign Website | Agilent Technologies | – | Target capture design software (https://earray.chem.agilent.com/suredesign/) |
UCSC Genome Browser | University of California Santa Cruz | – | rat Nov 2004 rn4 assembly |
Agilent Methyl-Seq Protocol | Agilent Technologies | – | https://www.agilent.com/cs/library/usermanuals/public/G7530-90002.pdf |