Genome-wide integration of chromatin accessibility and endogenous DNA break signals across cancer cell lines

Resumen Técnico: Integración en todo el genoma de la accesibilidad de la cromatina y las señales de roturas de ADN endógenas

Planteamiento del problema
Las roturas de la cadena de ADN endógenas son inherentes a procesos celulares normales como la replicación, la transcripción y el metabolismo. Cuando no se reparan, estas lesiones impulsan la inestabilidad genómica, un sello distintivo del cáncer. Si bien se sabe que los paisajes de mutaciones somáticas en genomas de cáncer se correlacionan con los estados de la cromatina (por ejemplo, la cromatina abierta suele mostrar tasas de mutación más bajas), los perfiles de mutación representan el resultado acumulativo de la formación, reparación y selección clonal del daño. Por consiguiente, no reflejan directamente la distribución inicial de las roturas de ADN endógenas. Además, la relación específica entre la accesibilidad de la cromatina y el paisaje de las roturas de ADN endógeno en diversas líneas celulares de cáncer no tratadas sigue estando insuficientemente caracterizada, particularmente en lo que respecta a si la cromatina abierta es intrínsecamente más susceptible a la formación de roturas o si simplemente es más accesible para la maquinaria de reparación.

Metodología
El estudio integró dos ensayos de secuenciación de alto rendimiento en 20 muestras biológicas derivadas de 10 líneas celulares de cáncer humano (que representan orígenes hematopoyéticos, de pulmón, próstata y riñón).

• mDEtail-seq (Multiplex DNA End tailing and sequencing): Un método adaptado de DEtail-seq para perfilar las roturas de ADN endógenas. A diferencia del DEtail-seq original, el mDEtail-seq utiliza ADN genómico purificado como entrada, lo que permite el procesamiento multiplexado de muestras. El protocolo implica el tratamiento con RNasa A, desnaturalización, adición de colas con códigos de barras (barcoded tailing) y ligación de splint, sonicación y secuenciación de alto rendimiento para capturar las señales de rotura de ADN en estado estacionario.
• ATAC-seq: Realizado en alícuotas de las mismas preparaciones celulares para mapear la accesibilidad de la cromatina.
• Integración de datos: Los autores correlacionaron las señales de mDEtail-seq con los picos de ATAC-seq, los sitios de inicio de la transcripción (TSS) y los cuerpos de los genes. También integraron conjuntos de datos públicos de RNA-seq para cinco líneas celulares para evaluar la relación entre la abundancia de ARNm y las roturas de ADN.
• Análisis: Los análisis estadísticos incluyeron correlaciones de Pearson y Spearman, Análisis de Componentes Principales (PCA) para el agrupamiento de muestras y cálculos de enriquecimiento comparando las fracciones de señal con las fracciones de cobertura genómica. Se realizó un enriquecimiento de la Ontología Génica (GO) en los genes con las señales de rotura de ADN más altas.

Resultados clave

1. Reproducibilidad y especificidad de la línea celular: Tanto mDEtail-seq como ATAC-seq demostraron una alta reproducibilidad entre réplicas biológicas (Pearson r = 0.93–0.98 para mDEtail-seq; 0.96–0.98 para ATAC-seq). El análisis de PCA confirmó que los patrones de distribución de las roturas de ADN son distintos y reproducibles entre diferentes líneas celulares.
2. Enriquecimiento en la cromatina accesible: Las señales de rotura de ADN se enriquecieron consistentemente dentro de los picos de ATAC-seq en las 10 líneas celulares. La relación de enriquecimiento (fracción de señal dividida por la fracción de cobertura genómica) osciló entre 1.37 y 2.66, lo que indica que las regiones de cromatina accesible contienen una fracción desproporcionadamente alta de roturas de ADN en comparación con su tamaño genómico.
3. Correlación con la intensidad de la accesibilidad: Se observó una asociación positiva entre la intensidad de la accesibilidad de la cromatina y los niveles de rotura de ADN. En las células RPMI-8226, el 1% superior de los picos de ATAC-seq más accesibles exhibió las señales de rotura de ADN más altas, con niveles de señal que disminuyen progresivamente a medida que la accesibilidad disminuye. Esta tendencia se mantuvo en las 20 muestras.
4. Enriquecimiento proximal al TSS: Las señales de rotura de ADN fueron significativamente más altas en regiones accesibles cerca de los Sitios de Inicio de la Transcripción (TSS) en comparación con las regiones inaccesibles dentro de las ventanas de ±200 a ±1,000 pb. El enriquecimiento más fuerte ocurrió inmediatamente en el TSS.
5. Distribución de características genómicas: Los promotores (definidos como TSS ± 2 kb) mostraron el mayor enriquecimiento de señales de rotura de ADN (valores de enriquecimiento 1.37–1.89), seguidos por los cuerpos de los genes (enriquecimiento modesto), mientras que otras regiones genómicas estuvieron relativamente depletadas.
6. Independencia de los niveles de expresión: A pesar del enriquecimiento de las roturas en regiones asociadas a genes, no se detectó una correlación positiva entre la abundancia de ARNm (RNA-seq) y las señales de rotura en el cuerpo de los genes en las cinco líneas celulares analizadas. Los genes con altos niveles de ARN no necesariamente exhibieron mayores señales de rotura.
7. Asimetría de la cadena: Se observó un sesgo débil pero consistente hacia la cadena molde en los cuerpos de los genes, con un enriquecimiento en la cadena molde significativamente mayor que en la no molde en las diez líneas celulares (P = 0.0020).
8. Enriquecimiento funcional: Los genes clasificados en el 1%, 5% y 10% superior por señales de rotura de ADN estuvieron significativamente enriquecidos para términos de la Ontología Génica relacionados con la estructura de la cromatina, el ensamblaje de nucleosomas y la organización de complejos proteína-ADN.

Significado y afirmaciones
Los autores afirman que este estudio proporciona una visión de todo el genoma de la distribución de las roturas de ADN endógenas en células cancerosas, demostrando que estas roturas se distribuyen de forma no aleatoria y se asocian preferentemente con la cromatina accesible.

• Desacoplamiento de las roturas y las mutaciones: Los hallazgos sugieren que los paisajes de mutaciones somáticas no pueden servir como indicadores directos de la distribución inicial del daño al ADN. Las regiones con cargas de mutación bajas (a menudo cromatina abierta) pueden, de hecho, experimentar niveles altos de roturas de ADN endógenas, probablemente debido a una reparación eficiente o a un procesamiento de daño activo en estas regiones.
• Contexto de la cromatina: El estudio establece que la cromatina abierta representa un entorno genómico asociado con mayores roturas de ADN endógenas. Esto respalda la hipótesis de que el ADN con depleción de nucleosomas es más vulnerable a la formación de roturas, mientras que simultáneamente es más accesible para los factores de reparación.
• Información mecanística: Los resultados indican que la accesibilidad local de la cromatina, más que la producción transcripcional por sí sola, es el factor principal asociado con los niveles elevados de rotura de ADN cerca de los promotores.
• Limitaciones: Los autores declaran explícitamente que su análisis es correlativo. No pretenden haber establecido un vínculo causal directo donde la accesibilidad de la cromatina cause el aumento de la formación de roturas. Señalan que mDEtail-seq captura roturas en estado estacionario, lo que significa que las señales altas podrían resultar de un aumento en la formación, un retraso en la reparación, o ambos. Se requieren experimentos de perturbación futuros para distinguir estos mecanismos.

En conclusión, el artículo argumenta que el perfilado directo de las roturas de ADN ofrece un marco complementario al análisis de mutaciones para comprender la inestabilidad genómica, revelando que la cromatina accesible es un contexto local clave para el enriquecimiento de las roturas de ADN en las células cancerosas.