Summary
Aquí, presentamos un protocolo para permitir la cuantificación exacta de la metilación de la DNA (mtDNA) mitocondrial. En este protocolo, describimos una digestión enzimática del ADN con BamHI juntada con una tubería de análisis bioinformáticos que se puede utilizar para evitar la sobreestimación de los niveles de metilación de mtDNA causada por la estructura secundaria del ADNmt.
Abstract
Cuantificación de la metilación del ADN puede lograrse utilizando secuenciación de bisulfito, que aprovecha la propiedad de bisulfito de sodio para convertir unmethylated citosina en uracilo, en un contexto de ADN monocatenario. Secuenciación de bisulfito puede ser dirigida (mediante PCR) o en todo el genoma y proporciona la cuantificación absoluta de la metilación de la citosina en la sola base-resolución. Dada la naturaleza distinta del ADN nuclear y mitocondrial, en particular en la estructura secundaria, deben realizarse adaptaciones de métodos de secuenciación de bisulfito para la investigación de metilación de citosina en el ADNmt. Estructura secundaria y terciaria del mtDNA puede de hecho conducir a bisulfito de secuenciación artefactos llevando a falsos positivos debido a la deficiente acceso desnaturalización incompleta de bisulfito a ADN monocatenario. Aquí, describimos un protocolo con una digestión enzimática del ADN BamHI juntada con tubería de análisis bioinformáticos para permitir la cuantificación exacta de la citosina niveles de metilación en el ADNmt. Además, proporcionamos las pautas para el diseño de los cebadores de secuencia de bisulfito específicos de mtDNA, para evitar dirigidos a segmentos NUclear mitocondrial indeseables (NUMTs) insertan en el genoma nuclear.
Introduction
El genoma mitocondrial es una estructura circular, doble cadena de aproximadamente 16,5 kilos base (kb), lo que constituye de un pesado y un filamento de la luz. El genoma mitocondrial está presente en múltiples copias en cada célula, maternal heredado y codifica componentes esenciales de la cadena respiratoria complejos1. Similar a los genomas bacterianos y a diferencia del genoma nuclear, el genoma mitocondrial se organiza en numerosas estructuras secundarias y terciarias, como en las estructuras en espiral y superenrollado2, que puede hacer que acceso difícil durante la secuencia experimentos3.
En el núcleo, la metilación del ADN es una marca epigenética ampliamente estudiada que interviene en numerosos procesos, en particular en la regulación de la expresión génica. En genomas mamíferos, la metilación del ADN ocurre principalmente en la posición 5 del anillo de pirimidina de deoxycytidines, sobre todo en los dinucleótidos CG (o CpG). Metilación de la citosina se encuentra en el 70% de los CpG en el genoma de las células somáticas y ~ 1% del total que bases de ADN4. Metilaciones del ADN se ha descrito también en contextos no-CpG, CpC, CpA y CpT y existen en cantidades diferentes en el ADN nuclear, con valores de hasta un 25% de todos los cytosines metilados en células madre embrionarias5,6,7.
Mientras que la metilación de la citosina del genoma nuclear es ampliamente aceptada, la existencia de la metilación de la DNA (mtDNA) mitocondrial es aún controversial. La metilación de mtDNA investigadora primer estudio fue realizada en células cultivadas donde mtDNA metilación fue detectada fácilmente, aunque en niveles más bajos en comparación con el ADN nuclear8. En las células humanas y murinas, metilación del mtDNA también fue detectada en niveles bajos (2-5%). Usando análisis basándose en 5 METILCITOSINA captura como inmunoprecipitación de ADN metilada (MeDIP) seguida de PCR cuantitativa, metilación del mtDNA también fue detectada en varios ratón y humanos y células líneas9,10, 11,12. Usando los anticuerpos contra 5-METILCITOSINA en un ensayo ELISA o espectrometría de masas, niveles considerables de metilación del ADN fueron detectados de fracciones purificadas mitocondriales13,14,15, 16. sin embargo, la mayoría de los ensayos en los estudios antes mencionados utiliza técnicas que no fueron diseñadas para proporcionar la cuantificación absoluta de la metilación del ADN en la sola resolución de base.
Análisis de metilación de DNA cuantitativo y resolutivo puede lograrse mediante una técnica llamada "bisulfito de secuenciación", que se aprovecha de la propiedad de bisulfito de sodio para convertir unmethylated citosina en uracilo en single-stranded DNA contexto17 . Mediante secuenciación de bisulfito, una constelación de estudios ha detectado la presencia de metilación de citosina en los distintos niveles. MtDNA de la metilación de la región D-loop, el 12S o la región 16S fue detectado fácilmente en humanos18,19,20,21,22,23 y ratón24 tejidos y células, sin embargo, con una variabilidad intrigante, de 1-20% del total cytosines en estudios.
En comparación con estos numerosos estudios, sólo unos pocos estudios, como de nuestro grupo, han disputado la presencia de mtDNA metilación3,25,26,27 o cuestionado la importancia biológica de niveles muy bajos de mtDNA (por debajo de 2%)28. Recientemente, informó de la observación de un artefacto de bisulfito de-secuencias posibles en todo mitocondrial bisulfito secuencia3. Proporcionamos evidencia que la estructura secundaria del ADN mitocondrial podría dar lugar a falsos positivos en secuenciación de bisulfito, así sobreestimar los niveles de metilación. Aquí proporcionamos un protocolo para evitar que un artefacto de bisulfito-conversión del ADNmt. Este protocolo utiliza una simple digestión enzimática del ADN para desarticular estructuras secundarias de mtDNA y permitir un acceso completo a siguiendo un protocolo de secuencia de bisulfito bisulfito. Además, ofrecemos un acompañamiento pipeline bioinformática para el análisis de secuenciación de bisulfito.
Subscription Required. Please recommend JoVE to your librarian.
Protocol
1. tratamiento enzimático de la restricción
- Alinear el mtDNA por tratar el DNA humana total con la enzima de restricción BamHI, que corta en la posición 14258 en el ADN mitocondrial humano.
Nota: Para ADN de ratón, utilice la enzima de restricción BglII en condiciones idénticas como abajo. - Para cada muestra donde metilación del mtDNA es evaluarse, preparar un tubo de reacción de 0.2 mL y añadir el siguiente mezcla: 3 μg ADN genómico (cuantificado por fluorometría), 15 μL tampón 3, 3 μL BamHI y agua hasta de 150 μL
- Colocar los tubos en un termociclador durante 4 h a 37 ° C.
2. bisulfito conversión
- Convertir el ADN BamHI-tratado con bisulfito de sodio como se describe en otra parte29.
- Utilizar 200-500 ng de ADN tratados con BamHI para cada reacción de conversión en un volumen total de 20 μL. La recuperación de ADN es 50-70%.
- Eluir la DNA bisulfito-convertido en tampón de elución de 10 μL.
- Cuantificar el ADN por fluorometría.
- Opcional: Almacenar DNA bisulfito-convertido a-20 ° C antes de continuar el resto del protocolo. Para el almacenamiento a largo plazo, almacenar DNA bisulfito-convertido a-80 ° C
3. diseño de bisulfito de cebadores de secuencia
- Diseño de iniciadores para la secuenciación de bisulfito usando las herramientas en línea Methprimer30 y31,de BiSearch32.
Nota: Methprimer y BiSearch permiten para buscar secuencias de primer enlace en secuencias genómicas donde cytosines se convierten en thymines, semejantemente al bisulfito post tratamiento. - Para la amplificación óptima e imparcial, evitar dinucleótidos CpG en las cartillas y diseñar el tamaño de amplicón entre 100-300 bp.
- Que sean diseñados primers específicos para el ADN mitocondrial usando la función de BiSearch Primer búsqueda. Inserte los cebadores ya diseñados bisulfito convertido, marque la casilla de bisulfito e incluyen un genoma de referencia y parámetros de la polimerización en cadena.
Nota: Se mostrará una lista de posibles productos PCR. - Copiar la parte superior de secuencias de la cartilla de 1-5 y pegar en un formulario de pedido para síntesis de oligonucleótidos.
Nota: Una lista de los humanos y secuencias de la cartilla de mtDNA ratón que ha sido validado en el laboratorio se muestra en la tabla 1. - Cartillas de prueba utilizando bisulfito convierten PCR después de electroforesis en gel de agarosa como se describe en el paso 4. Probar y optimizar por separado todos los pares de la cartilla y el tamaño de amplicón correcta aparece en el gel de agarosa.
Nota: Si necesitan cebadores multiplexación, añadir cartillas múltiples a una sola reacción de PCR y visualizar los productos PCR en electroforesis en gel de agarosa o, un método más resolutivo como electroforesis en gel de poliacrilamida, para distinguir productos de similares tamaño.
4. bisulfito convertido PCR y extracción del Gel
- Para amplificar las regiones de interés, realizar PCR utilizando una polimerasa de Taq hot start. Brevemente, mezcla 100 ng convertido ADN, 500 μm cebadores sentido y anti-sentido, mezcla de dNTP 1 μl (10 μm de cada dNTP) y polimerasa en 7,5 U reacción en un volumen total de 50 μl. ejecutar PCR con las siguientes condiciones: 5 min a 95 ° C (60 s 94 ° C, 60 s 55 ° C *; 60 s 72 ° C)
35 ciclos; 10 min 72 ° C. Utilizar cebadores ya sea por separado o multiplexados. - Opcional: Al multiplexación cartillas, uso 200 ng convierte en ADN y las siguientes condiciones bicicleta: 5 min 95 ° C (60 s 94 ° C; 90 s 55 ° C *; 90 s 72 ° C) 35 ciclos; 10 min 72 ° C.
Nota: La temperatura puede variar dependiendo de la temperatura de fusión de los primers. - Productos PCR sujetos a gelificarse electroforesis en agarosa al 2% de 100 voltios.
- Bajo la suave luz de UV, extraer productos de PCR con un bisturí y añadir a microtubos. No exponga los productos PCR a la luz UV excesiva para evitar daños en el ADN. Evitar el tiempo de exposición de más de 30 s.
- Purificar los productos PCR usando gel de purificación como se describió anteriormente3.
- Cuantificar el ADN purificado de paso 4.5 mediante fluorometría. Proceda al paso 5 o del almacén de muestras a-20 ° C.
35 ciclos; 10 min 72 ° C. Utilizar cebadores ya sea por separado o multiplexados.5. secuencia biblioteca elaboración de bisulfito
- Realizar la preparación de la biblioteca de productos PCR convertido a bisulfito.
Opcional: múltiplex productos de PCR en esta etapa. - Realizar reparación final con hasta 100 ng del producto de PCR. Agregar 3 μl final preparación enzima y 6,5 μl finales reparación reacción de tampón (10 x) a 55,5 μl del producto de PCR en un tubo de 0.2 mL. Lugar e incubar los tubos en un termociclador durante 30 min a 20 ° C, seguida por 30 min a 65 ° C.
- Ligan adaptadores de secuencia mezclando el ADN reparado final con 15 μl Blunt/TA ligasa mezclar, 2,5 μl adaptador y potenciador de la ligadura de 1 μl. Si utiliza < 100 producto PCR de ng como entrada, adaptador de diluir 10 veces.
- Colocar la mezcla (paso 5.3) en un termociclador e incubar 15 min a 20 ° C. Hacer una pausa en el termociclador y agregar 3 μl usuario enzima al tubo. Mezcla y vuelva a colocar el tubo en el termociclador e incubar 15 min a 37 ° C.
- Tamaño Seleccione el ADN ligada adaptador usando granos de sólido fase Reversible inmovilización (SPRI) con diferentes proporciones de granos al ADN dependiendo de tamaño del producto PCR.
- Añadir 13,5 μl de H2O a ADN ligada adaptador de paso 5.4 y 55 μl resuspendió SPRI granos. Incubar a temperatura ambiente durante 5 minutos y lugar el tubo sobre un soporte magnético. Cuando la solución esté clara, transferir el sobrenadante a un tubo nuevo y deseche el tubo que contiene los granos.
Nota: El sobrenadante contiene el ADN ligada adaptador y grandes fragmentos no deseados están limitados a los granos descartados. - Añadir 25 μl resuspendió SPRI cuentas al sobrenadante del paso 5.6, mezclar e incube durante 5 min a temperatura ambiente. Coloque el tubo en el soporte magnético y descarte el sobrenadante cuando la solución es clara.
Nota: El sobrenadante desechado contiene ADN no deseado, mientras que el ADN ligada por el adaptador está enlazado a granos. - Mientras el tubo está en el soporte magnético, agregar 200 μL 80% de etanol (recién preparado) para lavar los granos. Incubar 30 s y quite y descarte el sobrenadante. Repita este paso para un total de 2 lavados.
- Aire seco granos durante 5 minutos mientras que el tubo está en soporte magnético con una tapa abierta.
- Retirar el tubo del imán, eluir el ADN en 23 μl de tampón de elución y mezclar. Incubar a temperatura ambiente durante 2 minutos.
- Coloque el tubo sobre un soporte magnético y cuando la solución esté clara, transferir 2 μL en una placa PCR de 96 pocillos para la amplificación de pre-PCR (paso 5.14).
- Transferir el resto aproximadamente 21 μl a un nuevo tubo 0.2 mL para amplificación por PCR (paso 5.17).
Nota: No olvide transfiera cualquier granos como granos pueden inhibir reacciones enzimáticas de los próximos pasos. - Proceder a paso 5.14 o almacén de muestras a-20 ° C.
- Realizar la amplificación por RT-qPCR para estimar el número de ciclos necesarios para amplificar el ADN ligada adaptador pre-PCR. Por reacción, añadir lo siguiente en la PCR de 96 pocillos de la placa que contiene de 2 μl ADN (paso 5.11): 0,4 μL índice la cartilla, cartilla PCR universal de 0.4 μL, 7,2 μl de H2O, 10 μl RT-qPCR master mezclar.
- Colocar la placa en una máquina PCR en tiempo real y someter la placa a las siguientes condiciones: 30 s a 98 ° C (10 s 98 ° C, 75 s 65 ° C) 20 ciclos; 5 min 65 ° C.
- Estimar el número de ciclos de amplificación necesarios para cada muestra al restar un valor de Ct para el Ct-valor de la amplificación previa RT-qPCR correspondiente al final de la fase lineal de la reacción de PCR.
- Amplificar el ADN ligada adaptador de paso 5.12 mezclando: 21 μl ADN, 25 μl PCR master mix, 1 μl índice la cartilla, cartilla de PCR Universal de 1 μl y 2 μl de H2O. En un termociclador, tema cada muestra a las siguientes condiciones: 30 s a 98 ° C (10 s 98 ° C, 75 s 65 ° C) [Ct calculado en el paso 5.16] ciclos; 5 min 65 ° C.
Nota: Recuerde utilizar solamente un primer índice por muestra para permitir la multiplexación. El índice se inserta durante la PCR. - Purificar las cuentas SPRI ADN amplificadas y eluir en 22 μL de tampón de elución.
- Control de calidad de la biblioteca por gel electroforesis o alternativa método. Buscar tamaños de pico biblioteca correcta (tamaño de producto PCR más adaptador).
- Estimar el tamaño promedio de pares de base de los productos de la biblioteca por el análisis de frotis: en avanzada configuración global, haga doble clic en "tabla" bajo análisis de frotis. Definir la región de 100-1000 bp y haga clic en Aceptar. Debajo de la Mesa de la región, la región definida aparece y se calcula el tamaño promedio de la biblioteca (bp).
- Cuantificar las bibliotecas por fluorometría. Proceda al paso 6 o bibliotecas de tienda a-20 ° C.
6. la siguiente secuencia de generación
- Calcular la concentración de nM de bibliotecas mediante la fórmula:
- Diluir las bibliotecas para el mismo nM de concentración por ejemplo de 2 nM (elegir 4 nM, 2 nM, 1 nM y 0,5 nM dependiendo de las concentraciones de la biblioteca). Todas las bibliotecas para conseguir una piscina de 2 bibliotecas nM de la piscina.
- Desnaturalizar y diluir las bibliotecas siguiendo la guía de instrumento de la secuencia. En Resumen: desnaturalizar la piscina nM 2 mediante la adición de 0.2 N NaOH e incubar a temperatura ambiente durante 5 min diluir la piscina desnaturalizada a una dilución final de 22:00. Desnaturalizar y diluir una biblioteca de Control disponible en el mercado a una dilución final de 12:05. A un tubo nuevo, mezcla de la piscina de 22:00 con la biblioteca de Control 20% 12:05.
Nota: Conversión de bisulfito presenta muy baja complejidad. 20% control biblioteca spike en resultará en una mejor calidad de la secuencia. - Ejecute una secuencia de extremo apareado de 150 bp usando un instrumento de secuenciación profunda.
7. Computación Análisis - niveles de metilación de estimación
- Generar archivos de .fastq (aquí llamados sample1_R1.fastq.gz y sample1_R2.fastq.gz), generar un índice de bismark de la totalidad del genoma de interés (aquí en una carpeta llamada bismarkIndex) y asegúrese de que la secuencia genomic de chrM está disponible en fasta formato (aquí en una carpeta llamada chrM).
- Procesamiento previo de Lee para eliminar el sesgo introducido por los iniciadores y adaptadores: utilizar la herramienta Trim_galore33. Quite dos nucleótidos del extremo 5' de Lee hacia adelante y hacia atrás.
trim_galore--clip_R1 2--clip_R2 2 -o. /--trim1 - junto sample1_R1.fastq.gz sample1_R2.fastq.gz - Alinee el preprocesado lee el genoma entero con Bismark34
Bismark, cola de pato - bam -o. /. / bismarkIndex / -1./sample1_R1_val_1.fq.gz -2./sample1_R2_val_2.fq.gz
Nota: Si se utiliza otro alineador, asegúrese de que dice que no puede ser alineado únicamente se descarta. - Extracto de Lee traz al cromosoma M, utilizando aquí Samtools35
samtools tipo -l 0 - O BAM -o sample1_sorted.bam
sample1_R1_val_1_bismark_bt2_pe.BAM
samtools índice sample1_sorted.bam
samtools ve -u sample1_sorted.bam chrM | samtools tipo - n - O BAM -o sample1_chrM.bam - Extraer la información de la metilación
bismark_methylation_extractor -p--CX--cytosine_report--completo--genome_folder. / chrM./sample1_chrM.bam - Use el .bedGraph, la COV o los archivos de CX_report.txt para su posterior análisis. Utilice el archivo sample1_chrM.CX_report.txt para visualización y pruebas.
- Realizar más recortes en el paso de preprocesamiento, si varias regiones son interrogadas y un sesgo de la cartilla puede ser visible en las parcelas M-sesgo generadas durante la asignación. Consulte la documentación de Bismark para obtener más información.
8. computación de análisis - prueba de las diferencias
- Asegúrese de que el archivo CX_report.txt contiene 7 columnas, cromosoma (aquí chrM), posición, filamento, número de lecturas metilados, el número de lecturas unmethylated, C contexto y secuencia circundante, para cada C en el cromosoma M.
- Como la CX_report cubre chrM en su totalidad, compruebe el subconjunto a la región o regiones que se amplifica.
- Utilizar pruebas no paramétricas para las diferencias entre las muestras, como la prueba exacta de pescadores individuales c o una prueba de signo para toda la región.
Nota: Cuantificación a menudo estará cerca de 0% de metilación, que es la razón por qué no es apropiado asumir normalidad. - Del mismo modo, para la visualización, visualizar cuantiles o calcular intervalos de confianza de la proporción binomial.
Subscription Required. Please recommend JoVE to your librarian.
Representative Results
Dos pasos en este protocolo son cruciales al investigar metilación del mtDNA. 1) la apertura de la estructura secundaria y 2) el diseño de iniciadores específicos de ADN mitocondriales.
Por digerir el ADN genómico humano con la enzima de restricción BamHI (figura 1), se reducirá la estructura del ADN mitocondrial en la posición del nucleótido 14.258 y se abrirá la estructura secundaria.
El posible deterioro del valor de conversión de bisulfito con una estructura secundaria de mtDNA intacto es muy probable que sobrestimar la tasa de unconversion de mtDNA. Ordenando del bisulfito, la tasa de unconversion de mtDNA fue investigada en sin digerir y digiere el ADN total de músculo esquelético humano células, en 5 regiones diferentes del genoma mitocondrial: el lazo de desplazamiento (D-Loop), el ribosoma de ARNt de la fenilalanina y 12S (tRNA - F + 12S), ribosoma 16S (16S), NADH deshidrogenasa 5 (ND5) y citocromo b (CYTB) codificación de mRNA del gen (figura 2). Tasa de unconversion entre 0 - 15.1% en todas las regiones investigadas para el mtDNA no digerida. Sin embargo, cuando el ADN fue digerido antes de la conversión de bisulfito, la unconversion tasa se redujo a un máximo de 1% en todas las regiones (figura 3). Así, ilustrando la importancia de la digestión BamHI antes de la conversión de bisulfito para evitar la sobreestimación de los niveles de metilación de mtDNA.
Contaminación de ADN nuclear es inevitable al aislar el genoma mitocondrial y desde casi el entero genoma mitocondrial se ha replicado en el genoma nuclear, ha dado lugar a inserciones nucleares del genoma mitocondrial denominado NUclear Segmentos mitocondrial (NUMTs)36. Para asegurar que los niveles de metilación de ADN analizados corresponden a ADN mitocondrial y no NUMTs, es de más importancia para el diseño de primers específicos del genoma mitocondrial. La figura 4 muestra cómo comprobar la especificidad de la cartilla y tabla 1 pantallas humanas y secuencias de la cartilla de mtDNA de ratón.
| Nombre | Secuencia de | Posición | Tamaño de amplicón | Temperatura de recocido | Strand | ||
| Humano | |||||||
| 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 | |||
| 062 | F: AATAAATTTATAGGTTTTTAAATTATTAAAT | ||||||
| R: TAACTAATAAAATCTTAACATATACTACTC | 2763-2873 | 110 | 53° C | Sentido | |||
| CYTB | F: GGTATTATTTTTTTGTTTGTAATTATAGTA | ||||||
| R: CCTCAAATTCATTAAACTAAATCTATCC | 15091-15243 | 152 | 55° C | Sentido | |||
| Ratón | |||||||
| 12S | F: ACACATACAAACCTCCATAAAC | ||||||
| R: GAGGTATAATTTAGTTAAAT | 111-287 | 176 | 53° C | Antisentido | |||
| 062 | F: AAATTCCAATTCTCCAAACATAC | ||||||
| R: TTTGGATTTTTTTTTTAGGT | 1749-1896 | 147 | 53° C | Antisentido | |||
| CYTB | F: CTACAAAAACACCTAATAACAAAC | ||||||
| R: AGTGTATGGTTAAGAAAAGA | 14129-14307 | 178 | 55° C | Antisentido | |||
| 062 | F: AGAGAAATAGAGTTATTTTATAAATAAG | ||||||
| R: CTAAAAACAAAATTTTAAATCTTAC | 2576-2727 | 151 | 55° C | Sentido | |||
| ND5 | F: TTAAGTTAATTAGGTTTGATAATAGTGA | ||||||
| R: TAAAACCCTATTAAAAATAATATTCCTATA | 12659 12910 | 251 | 55° C | Sentido | |||
| CYTB | F: TGGGTTTTTTTTAGGAGTTTGTTTA | ||||||
| R: CAATAAAAATACTCCAATATTTCAAATTTC | 14242-14504 | 262 | 55° C | Sentido | |||
Tabla 1: secuencias de la cartilla para el ser humano y ratón mtDNA. Cuenta que el recocido temperatura de PCR cartillas varían debido a diferencias en la temperatura de fusión de los primers.
Figura 1 : Gel de electroforesis de ADN BamHI digerido o no digerido. ADN genómico de músculo esquelético humano fue no tratado o tratado con BamHI de 4h a 37° C y electroforesis en gel se ejecutó en un gel de agarosa al 1.5%. Haga clic aquí para ver una versión más grande de esta figura.
Figura 2 : Mapa del genoma mitocondrial humano. Se muestran las regiones del genoma mitocondrial investigados por secuenciación de bisulfito y el sitio de enzimas de restricción BamHI. El ADN mitocondrial contiene 13 mRNAs, 2 ARN ribosomal y 22 RNAs de transferencia. Abreviaturas: PH1: promotor de la pesado-cadena 1; PH2: cadena pesada promover 2; PL: promotor del filamento de la luz; OH: Origen de replicación de la cadena pesada; OL: origen de replicación del filamento de la luz. Haga clic aquí para ver una versión más grande de esta figura.
Figura 3 : Digestión BamHI antes de la conversión de bisulfito reduce la tasa de unconversion de mtDNA en células del músculo esquelético humano. Por secuenciación de bisulfito, digerido y no digerido porcentaje de unconversion mitocondrial de ADN fueron interrogados en cinco regiones diferentes del genoma mitocondrial en células del músculo esquelético humano (N= 3). El cuadrado completo representa sitios de CpG, y la plaza abierta representa sitios de CpG no. Los resultados se presentan con un intervalo de min-max y una prueba de la muestra se utilizó para probar las diferencias significativas unconversion. 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; ND5: P= 1.70E-05; CYTB: P= 2.69E-10. D-Loop (6-298) incluye el origen de replicación y tRNA - F + 12S promotor cadena pesada 2. Estos datos han sido publicados en otros lugares 3. Haga clic aquí para ver una versión más grande de esta figura.
Figura 4 : El uso de Bisearch para verificar la especificidad de la cartilla hacia el genoma mitocondrial. A) secuencias de primer convertido de bisulfito se agregan a la función BiSearch (http://bisearch.enzim.hu/?m=genompsearch) la búsqueda de la cartilla. Está marcada la casilla de bisulfito, y un genoma de referencia seleccionado. B) ejemplo de potenciales productos PCR generados en las cadenas de sentido y antisentido. Haga clic aquí para ver una versión más grande de esta figura.
Subscription Required. Please recommend JoVE to your librarian.
Discussion
Presentamos un protocolo de secuenciación del bisulfito que está específicamente diseñado para interrogar a metilación del mtDNA. Las diferencias con la secuenciación del bisulfito de protocolos utilizados para la extracción de ADN genómico se encuentra en la utilización de un paso de digestión de la enzima de la restricción anterior y un análisis Bioinformático predominante hacia fuera falso positivo derivados NUMT secuencias.
Proporcionamos un protocolo para evitar artefactos de secuenciación del bisulfito al investigar metilación del mtDNA. Artefactos de la secuencia de bisulfito a falsos positivos fueron previamente divulgados37. Eliminación insuficiente de proteínas accesorias de ADN por proteinasa que k ha sido descrito37. Desnaturalización incompleta o parcial reannealing puede llevar a tramos de la conversión incompleta, posiblemente al reducir el acceso de bisulfito a single-stranded DNA37. El artefacto más adelante podría estar en juego en el ADNmt donde una estructura secundaria deterioraría acceso a cytosines. A lo mejor de nuestro conocimiento, la adición de un paso de digestión enzimático de la restricción antes de la conversión de bisulfito en el contexto de estimar la metilación del mtDNA ha divulgado primero en 20163,28. En nuestras manos, previa digestión con una enzima de la restricción de solo corte dramáticamente baja detección de inconversos cytosines3,28.
Aquí, ofrecemos una tubería de bioinformática para analizar los resultados de la secuencia de bisulfito y detectar artefactos de secuenciación de bisulfito en el contexto de la secuenciación del genoma mitocondrial entero bisulfito. Este método origina a partir de nuestra observación que mtDNA no-tasa de conversión se correlacionó inversamente con secuencia profundidad3,28. Esta observación sugiere que la estructura secundaria y hasta terciaria del mtDNA deteriora la fragmentación del mtDNA en regiones específicas durante la construcción de bibliotecas de la secuencia y también podría afectar el acceso completo de bisulfito de sodio para el mtDNA. Puesto que la digestión con una enzima de un cortador lanza estructuras secundarias y terciarias, nuestros resultados indican que la estructura de mtDNA es una fuente de artefacto de bisulfito y validar los métodos descritos para la cuantificación de la metilación de la citosina en el ADNmt.
La inevitable contaminación de ADN nuclear en preparación mitocondrial constituye un desafío. Esto conduce a la contaminación en secuencias NUMT en Lee de la secuencia, que pueden sesgar la estimación de la metilación de la citosina en el ADNmt, incluso cuando realiza y secuenciación del genoma mitocondrial no todo bisulfito. Para evitar estos sesgos, el diseño de los primers de secuencia única bisulfito permite para excluir la contaminación con la metilación en NUMTs. NUMTs abarcan algunas regiones del ser humano y el genoma mitocondrial de ratón con 100% de la identidad, y por lo tanto, es imposible diseñar cartillas únicas en algunas regiones de mtDNA. Cuando se realiza la secuenciación del genoma mitocondrial entero bisulfito, secuencia larga Lee puede facilitar la discriminación entre el ADNmt y NUMTs. Alineación Lee al genoma entero y utilice únicamente Lee únicamente asignado al genoma mitocondrial puede aún garantizar detección de metilación del mtDNA.
Por último, este protocolo puede utilizarse más allá de los genomas mitocondriales, por ejemplo, en otros casos donde la estructura secundaria del ADN representa una fuente potencial de artefacto de la secuencia de bisulfito. La asociación entre los niveles de metilación y profundidad de la secuencia puede investigarse para evaluar la necesidad de una digestión con una enzima de la restricción antes de secuenciación de bisulfito complejo de secuencias de ADN.
Subscription Required. Please recommend JoVE to your librarian.
Disclosures
Los autores tienen intereses financieros que compiten para declarar.
Acknowledgments
Novo Nordisk Fundación Centro para la investigación metabólica básica es un centro independiente de investigación en la Universidad de Copenhague, parcialmente financiado por una donación sin restricciones de la Fundación de Novo Nordisk.
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).
