$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Como prueba de concepto, ChIP-Seq se completó para seis donantes humanos con tres conjuntos de tipos de células inmunes: células T CD4 ingenuas (CD4), monocitos clásicos (MO) y células asesinas naturales (NK), enriquecidas por clasificación FACS como se describió antes13. El procedimiento subrayado consta de nueve procedimientos distintos como se representa en la Figura 1.

Figura 1: Diagrama de flujo general para el procedimiento. (A) Una caricatura del procedimiento general (generado en BioRender). (B) Diagrama de flujo para todos los pasos principales para el protocolo y el tiempo práctico y total estimado asociado con cada día. La secuenciación podría ocurrir al final del Día 5 o más tarde con múltiples rondas. La línea de tiempo también se puede escalonar a lo largo de la semana, donde los días secuenciales 3-4 se pueden completar varias veces en una semana para generar 48 muestras de ChIP. Haga clic aquí para ver una versión más grande de esta figura.
Después del aislamiento celular por citometría de flujo13,las células clasificadas se centrifugaron y las células se fijaron y almacenaron como se describió anteriormente. Una vez que se recogieron todas las muestras, las muestras se lisaron y prepararon para el cizallamiento cromático en lotes de 12 como se describió anteriormente. Para cada muestra, se completó el número de ciclos para alcanzar la sonicación óptima10. La medición cuantitativa, así como las mediciones del tamaño del fragmento de cromatina cortada mostraron una gran reproducibilidad de nuestro método en los tres conjuntos de células inmunes(Figura 2A). Las diferentes células inmunes humanas se sonicaron en lotes separados y rindieron de manera muy consistente con > el 70% de la muestra entre 100 y 500 pb durante 14 ciclos (16 s ON, 32 s OFF por ciclo). En este punto, las muestras con fragmentos grandes después de la sonicación (< el 70% de la muestra entre 100 - 500 pb) se consideraron fallidas. Estas muestras podrían sonicarse durante 1-2 ciclos adicionales o fueron descartadas y reemplazadas más tarde con células de otro pellet. Nuestro método mostró que ninguna de las muestras requirió más sonicación o fue eliminada, lo que sugiere una robustez absoluta del procedimiento.

Figura 2: Ejemplos de control de calidad de presecuenciación. (A) Los geles de agarosa al 1,2% muestran reproducibilidad de la sonicación. Muestras de sonicación para 6 donantes en tres tipos de células: células T CD4 ingenuas (CD4), monocitos clásicos (MO) y células asesinas naturales (NK). Las muestras fueron sonicadas durante 14 ciclos (16 s ON, 32 s OFF por ciclo). Para cada muestra, se cargaron alrededor de 200 ng de cromatina desenlazada en un gel de agarosa al 1%. Las muestras se consideraron buenas si más del 70 % de los fragmentos están dentro de 100-500 pb. (B) Top - Análisis de curvas de amplificación qPCR para determinar el número óptimo de ciclos para la amplificación (Ct donde hay 1/2 de la intensidad máxima). Las muestras ideales tienen un Ct de aproximadamente 15 y la amplificación se puede completar hasta 2 ciclos más del Ct medido. La flecha es un ejemplo de un mal ejemplo donde el Ct es mayor que 18. Abajo - Se muestra un ejemplo de un conjunto pobre de muestras que tienen un Ct mayor que 18. Estas muestras también mostraron una menor intensidad de fluorescencia. (C) Izquierda - Los rastros de electroforesis del analizador de fragmentos mostraron la distribución de las bibliotecas etiquetadas finales después de la amplificación y la selección de tamaño. Las muestras con más del 85 % de la biblioteca de fragmentos dentro de 200-1.000 pb se consideraron como buenas muestras. La medición de la intensidad máxima de la fluorescencia también se considera un parámetro importante de control de calidad, de hecho, si la señal es baja, es poco probable que la muestra se secuencia bien. Derecha: se muestran ejemplos de muestras positivas en CD4, MO y NK. Haga clic aquí para ver una versión más grande de esta figura.
Después de la cuantificación, las muestras se ejecutaron en un manipulador de líquidos ChIP con anticuerpos H3K27ac, seguido de la etiquetación con la enzima transposasa Tn5. Para determinar el número apropiado de ciclos de amplificación por qPCR, se utilizó el 10% de las muestras etiquetadas. Para la determinación del número de ciclos para la amplificación de las muestras, encontramos el ciclo en el que la intensidad de la muestra es la mitad del máximo promedio para la determinación del ciclo (Figura 2B). Las muestras con valores de Ct de más de 18 no tuvieron un buen desempeño después de la secuenciación y, por lo tanto, su valor de Ct fue indicativo de una muestra de ChIP fallida. Estas muestras generalmente también produjeron una menor cantidad de ADN después de la amplificación. Las muestras (entrada de 100.000 células) con un valor de Ct igual o inferior a 15 fueron ideales y las muestras entre 15 y 18 fueron aceptables pero menos consistentes después de la secuenciación. Para muestras con menos de 100,000 celdas de entrada, los valores de Ct generalmente se encontraron entre 15 y 18, pero no necesitaron más de 18 ciclos para producir suficiente producto para la secuenciación.
Después de la amplificación del ADN etiquetado, las bibliotecas se purificaron y seleccionaron el tamaño para obtener una distribución de tamaño ideal, que oscilaba entre 200 y 1.000 pb, para la secuenciación NextGen. La evaluación de la distribución del tamaño en cada una de las bibliotecas se completó porque se obtuvieron los mejores datos de secuenciación cuando más del 85% de los fragmentos de ADN oscilaban entre 200 y 1.000 pb(Figura 2C). En particular, a medida que se cargaba la misma cantidad de ADN (medida por cuantificación de fluorescencia), se notó que las muestras con menor intensidad de fluorescencia generalmente se secuenciaban mal(Figura 2C).
Post secuenciación, se aplicaron controles de calidad estándar basados en las directrices ENCODE ChIP-Seq5,14,15.

Figura 3: Reproducibilidad de las muestras de células inmunes. (A) H3K27ac rastrea (UCSC Genome Browser, intensidad máxima, función de suavizado de 4, todos con eje Y igualmente escalado) para 6 donantes (100,000 células por réplica) en cada tipo de célula (CD4, MO y NK). Se muestran cuatro loci ejemplares, dos con (locus IL2RA y PTPRC) y dos sin enriquecimiento para H3K27ac (CCR4 y MS4A1). (B) La correlación de Pearson entre los donantes y las gráficas de correlación correspondientes generadas utilizando una extensión de 300 pb y una ventana de 500 pb dentro del paquete MEDIPS para cada uno de los tipos de célula replica16. (C) Archivos de donantes fusionados para cada tipo de célula que muestran pistas H3K27ac (intensidad máxima del UCSC Genome Browser, función de suavizado de 4) en regiones específicas del tipo de célula (IL17R para CD4, CCR2 para MO y KLRC1 para NK) y el gen de mantenimiento doméstico B2M, presente en todos los tipos de células. Haga clic aquí para ver una versión más grande de esta figura.
Para el control de calidad visual, se prepararon pistas de enriquecimiento H3K27ac para su visualización en el navegador del genoma UCSC. Para cuatro loci genéticos, las pistas individuales para cada muestra mostraron una alta calidad de mapeo y una relación señal-ruido que refleja la alta consistencia y robustez de nuestro ensayo(Figura 3A). Los dos loci a la izquierda albergan genes bien expresados en estos tipos de células, mientras que los genes en los dos loci a la derecha no se expresan y sirven como controles de fondo13 (Figura 3A). Además, se utilizó el paquete de análisis MEDIPS como variable post-secuenciación para evaluar el índice de correlación entre réplicas técnicas(Figura 3B)5,16,17,estableciendo el grado de correlación para lecturas de nivel de enriquecimiento para contenedores de 500 pb16. Para la mayoría de las comparaciones por pares, los índices de correlaciones de Pearson mostraron una correlación de más del 90%, lo que sugiere un alto nivel de consistencia entre las réplicas biológicas(Figura 3B). Las réplicas con correlación aceptable se fusionaron para aumentar la relación señal-ruido. Mientras que los loci específicos del tipo celular mostraron un alto enriquecimiento en las células apropiadas, un gen de mantenimiento doméstico (B2M) mostró una modificación de histonas muy consistente(Figura 3C). Para el análisis, la fusión de pistas de réplicas aumentará el enriquecimiento, reforzará la señal específica, incluso para importantes potenciadores específicos del tipo de célula, y reducirá la variabilidad interindividual inherente a las muestras humanas5.
Aunque se utilizaron 100,000 células para este estudio, hubo una alta reproducibilidad para tan solo 10,000 células en una línea de células T cultivadas en humanos (HUT78). El análisis de correlación entre el conjunto de datos ChIP-Seq realizado a partir de muestras con menos de 100.000 células mostró una alta reproducibilidad y correlación hasta 10.000 células(Figura 4A).

Figura 4: Reproducibilidad de muestras de baja entrada. (A) Ejemplos de la consistencia de H3K27ac ChIP-Seq para células 100,000 hasta 10,000 en células HUT-78 (una línea celular de linfoma de células T). Las pistas (UCSC Genome Browser, intensidad máxima, función de suavizado de 4, todas con eje Y igualmente escalado) muestran el locus IL4. (B) Correlaciones de Pearson de las réplicas utilizando una extensión de 300 pb y una ventana de 500 pb dentro del paquete MEDIPS16. (C) Correlaciones de Pearson entre los diferentes grupos de números celulares (100.000, 50.000 y 10.000 células) utilizando los mismos parámetros MEDIPS que en (B)16. Haga clic aquí para ver una versión más grande de esta figura.
El análisis de correlación de Pearson mostró un alto índice de correlación (83% a 92%), lo que sugiere el mantenimiento de la señal en muestras de bajo número de células. Sin embargo, hubo un aumento de fondo a medida que se redujeron los números de células, así como una disminución de los coeficientes de correlación(Figura 4B). Para mantener señales de fondo bajas, se fusionaron duplicados técnicos y se probó la correlación entre grupos(Figura 4C).
| Búfer de fijación de celdas 10X | |
| Compuesto | Concentración final |
| Solución de formaldehído | 11% |
| NaCl | 100 metros |
| EDTA, pH 8.0 | 1 mM |
| EGTA, pH 8.0 | 0,5 mM |
| HEPES, pH 7.5 | 50 metros |
| Búfer de lisis completo | |
| Compuesto | Concentración final |
| Tris-HCI, pH 8.0 | 50 metros |
| EDTA, pH 8.0 | 10 mM |
| SDS | 0.25% |
| Butirato de sodio | 20 metros |
| Cóctel inhibidor de la proteasa | 1X |
| Tampón de lisis a corto plazo | |
| Compuesto | Concentración final |
| Tris-HCI, pH 8.0 | 50 metros |
| EDTA, pH 8.0 | 10 mM |
| SDS | 0.25% |
| Mezcla de tagmentación | |
| Compuesto | Concentración final |
| Tris-HCI, pH 8.0 | 10 mM |
| MgCl2 | 5 mM |
| N,N-dimetilformamida | 10% |
| Enzima de tagmentación Illumina | 1:24 vol:vol |
| Mezcla de CtD | |
| Compuesto | Por muestra (μL) |
| NextEra Índice Primer A (25 μM) | 0.275 |
| NextEra Índice Primer B (25 μM) | 0.275 |
| 2X KAPA HiFi HotStart Ready Mix | 2.75 |
| 1:1000 SYBR Tinte verde | 0.11 |
| Tinte pasivo ROX | 0.11 |
| Agua | Relleno a 4 μL |
| Mezcla amp | |
| Compuesto | Por muestra (μL) |
| 2X KAPA HiFi HotStart Ready Mix | 27.5 |
| Agua | Llenar a 31 μL |
Tabla suplementaria 1: Recetas tampón.
Tabla suplementaria 2: Correlaciones de muestras de Spearman y Pearson para los 6 donantes y cada tipo de célula. Haga clic aquí para descargar esta tabla.