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.
Ha pasado más de medio siglo desde el primer informe de origen natural modificaciones del ADN en forma de citosina metilación 1. Metilación de citosina es una base de nucleótido citosina modificada donde el 5º de carbono en el anillo de base tiene un grupo metilo (5-mC), que se producen con más frecuencia en el motivo dinucleótido CpG en genomas de mamíferos. La presencia funcional de 5-mC en los promotores de genes se asocia generalmente con la represión transcripcional, mientras que la ausencia se asocia con la actividad transcripcional 2.
Grandes avances se han hecho en nuestra comprensión del papel de la metilación del ADN en desarrollo 3, propagación TransGeneration de perfiles epigenéticos 4, la patogénesis del cáncer de 5,6, y una serie de otras áreas de investigación. Muchos de estos avances han llegado en los últimos años a medida que se han desarrollado nuevas metodologías para perfilar la metilación del ADN 7. Con el advenimiento ypopularización de secuenciación (NGS) técnicas de última generación en la actualidad hay una serie de métodos para perfilar la metilación del ADN a través de todo un genoma o grandes regiones de un genoma 8. Estos enfoques abren la posibilidad de realizar estudios de descubrimiento que cuantifican los patrones y las diferencias en la metilación del ADN de metilación del ADN. Además de estos métodos de generación de hipótesis, hay una necesidad significativa de ADN pruebas de métodos de cuantificación de metilación específica / hipótesis.
Pirosecuenciación, PCR específica de metilación y la secuenciación directa de ADN Sanger convertida bisulfito han sido los métodos más utilizados para el análisis de regiones específicas (es decir, una región promotora de un gen o de una isla CpG) 9-12. Todos estos métodos se basan en la conversión de bisulfito, una reacción química donde la desaminación de la citosina no metilada se convierten en uracilo de, mientras que la citosina metilados permanecen intactos 13,14. Amplificación por PCR de bisulfite-Los resultados de ADN convertidos en sustitución de de uracilo con timina 15 que permiten la metilación diferencial a leer como una diferencia de base. Aunque muy útil, las limitaciones a estos métodos incluyen baja precisión cuantitativa, la longitud de lectura corta, y bajo rendimiento de la muestra. Para hacer frente a estas limitaciones, hemos desarrollado un método que hemos llamado bisulfito Amplicon Secuenciación (BSAS) 16. El objetivo de este método es ser capaz de analizar diana regiones genómicas de interés en un gran número de muestras con alta precisión cuantitativa.
BSAS utiliza elementos de los métodos existentes (la conversión de bisulfito y la amplificación por PCR específica de la región) y los combina con una construcción simple de próxima generación biblioteca (transposome mediada) 17 y de sobremesa NGS 18 (Figura 1). Este enfoque proporciona un método rendimiento más cuantitativa y superior para examinar la metilación de citosina en cualquier región de interés. Aquí presentamos elmétodo en detalle y describir enfoques para la resolución de problemas y la calidad de comprobar el proceso BSAS.
BSAS permite centrado, precisa cuantificación de metilación del ADN con capacidades de alto rendimiento, tanto en número de muestras y número de objetivos. Además, esta generación biblioteca utilizando tagmentation transposome mediada requiere menos entrada de ADN, es más rápido y contiene menos pasos que esquila tradicional y técnicas de unión adaptador posteriores. Estas mejoras sobre los métodos existentes permiten precisa de nivel de base el análisis de metilación del ADN de cualquier diana de interés. Por otra parte, BSAS proporciona un nuevo método para la confirmación de las regiones de interés (regiones de metilación diferencialmente (DMR), islas CpG, regiones reguladoras) que han sido identificados utilizando todo el genoma de bisulfito 22 o capturar 23 enfoques. Utilizando métodos de confirmación ortogonales y más cuantitativamente métodos precisos mejora el rigor analítico de estudios epigenómicos. La profundidad máxima de lectura alcanzado con BSAS permite la cuantificación precisa en un intervalo de confianza estrecho, reresultante en la cuantificación precisa 16. Además, la capacidad de ejecutar muchas muestras en un solo experimento puede aumentar el tamaño de la muestra para la confirmación de los métodos de todo el genoma, lo que aumentará el poder estadístico; una debilidad de muchos estudios epigenéticos actuales. Es importante destacar que, BSAS proporciona una cuantificación de metilación CpG específicos de base, lo que permite la generación de hipótesis comprobables para la investigación epigenética. Hipótesis como el aumento de la metilación en un elemento de respuesta específica se traducirá en disminución de la unión de un factor de transcripción específico. BSAS, combinado con los datos de expresión de ARNm emparejados, proporciona un método para obtener tanto la regulación epigenética y la expresión del ARNm dentro de una sola pieza de tejido o población de células. Mediciones pareadas de regulación y expresión también aumentan el rigor científico y el impacto de los hallazgos.
Los pasos críticos dentro del protocolo BSAS incluyen primer diseño, optimización de amplificación, y la generación de biblioteca /QC. PCR sesgo se introduce si cebadores no están diseñados adecuadamente y amplifican en diferentes eficiencias 8. El uso de estándares de metilación, como 100% y 0% de citosina metilados normas enzimáticamente generadas, se puede utilizar para determinar el grado, si alguno, de sesgo cuantificación de metilación, y es un control metodológico adecuado cuando la cuantificación de metilación 16. Del mismo modo, se debe tener cuidado cuando la optimización de las reacciones de PCR para generar un único amplicón de alta calidad. Si no está presente amplicón utilizando las directrices sugeridas se detallan en el protocolo, pasos de optimización incluyen un número creciente de ciclo de PCR o la cantidad de entrada bsDNA. Si hay múltiples productos de PCR, incluyendo dímeros de cebadores, pasos de optimización incluyen la disminución de la cantidad de entrada bsDNA, la disminución de las concentraciones de entrada de cebadores PCR molares, aumento de las temperaturas de recocido, o disminuyendo el número de ciclos de PCR. Uno o una combinación de estos procedimientos de optimización da resultados suficientes para generar un único higamplicón h calidad. Alternativamente, el rediseño de cebadores es un buen enfoque para conjuntos de cebadores que no producen un solo amplicón. Para aquellas regiones donde rediseño no es adecuado o posible, degenerar conjuntos de primer incluyendo tanto pueden trabajar de citosina y timina de CpG en los sitios en los cebadores directos y adenina y guanina de de en los sitios CpG en cebadores inversos. Sin embargo, estos conjuntos de cebadores necesitan una mayor optimización para asegurar que no se está produciendo PCR sesgo, un proceso fuera del alcance de este protocolo. Bibliotecas de calidad son extremadamente importantes para el éxito. Asegúrese de que las medidas de control de calidad se llevan a cabo para determinar el tamaño, la calidad y cantidad de las bibliotecas antes de la secuenciación. Las reacciones de secuenciación son propensos a fallas con el aporte de las bibliotecas de baja calidad. Sesgo en la cuantificación de metilación puede ser mitigado por garantizar que las frecuencias de metilación están presentes tanto en avance y retroceso secuenciación lee. Además, las puntuaciones de calidad de las bases de la alineación de los sitios CpG deben ser ≥Q30 para alta confidence en la cuantificación. Mientras que otros han demostrado previamente poco o ningún sesgo en las reacciones tagmentation de ADN convertido bisulfito 19, garantizando los parámetros descritos anteriormente permite la confirmación de la baja cuantificación metilación sesgo.
BSAS actualmente mide la metilación en sitios CpG, sin embargo, las futuras ampliaciones de este método permitirá que las mediciones pareadas de CpG 5-HMC con la adición de etapas experimentales distinguiendo entre 5-mC y 5-HMC como la introducción de perrutenato de potasio antes de la conversión de bisulfito 24. Esto aumentará el impacto de utilidad BSAS permitiendo para la metilación emparejado, tanto los datos de expresión 5-mC y 5-HMC, y con el fin de determinar la expresión de genes funciones de regulación de la metilación del ADN. Preparación biblioteca mediada por transposome de amplicones de PCR da lugar a la disminución de la profundidad de secuenciación en los extremos de las regiones amplificadas. Por lo tanto, las regiones de alto interés deben ser diseñados para residir en el medio de amplicones. Cómonunca, porque BSAS genera secuenciación leer profundidades de> 1000 X, la precisión cuantitativa de 5-mC mediciones cerca de los extremos de las regiones amplificadas sigue siendo alta.
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 |