Whole-genome Sequencing of Reed-sternberg Cells Identifies Outcome Determinants in Hodgkin Lymphoma

Resumen Técnico: La secuenciación del genoma completo de las células de Reed-Sternberg identifica determinantes de los resultados en el linfoma de Hodgkin

Planteamiento del problema
El linfoma clásico de Hodgkin (LCH) es una neoplasia derivada de células B caracterizada por la presencia de raras células malignas de Hodgkin y Reed-Sternberg (HRS) dentro de un denso microambiente inflamatorio. Aunque la mayoría de los pacientes logran una remisión duradera con la terapia de primera línea, un subgrupo significativo experimenta recaída o progresión refractaria. Las herramientas pronósticas actuales, como la imagen PET intervencionista, no logran capturar la complejidad molecular de la enfermedad y carecen de sensibilidad para la detección de la enfermedad residual mínima. El valor pronóstico de las características genómicas globales —específicamente las variaciones en el número de copias (CNV), las variantes estructurales (SV) y las firmas mutacionales— en las células HRS purificadas aún no está suficientemente definido. Este estudio aborda la necesidad de mejorar las estrategias para identificar a pacientes con mayor riesgo de fracaso terapéutico mediante la caracterización del panorama genómico completo de las células HRS.

Metodología
El estudio empleó un diseño prospectivo multicéntrico que involucró a 76 pacientes con LCH confirmado.

• Procesamiento de muestras: Biopsias diagnósticas en bloques de parafina embebidos en formalina (FFPE) se disociaron en suspensiones de células individuales. Las células HRS se purificaron mediante la tecnología DEPArray™ basada en la expresión de CD30 y la morfología de célula grande, obteniendo grupos masivos (bulk) de células HRS (media de 24.5 células por muestra).
• Secuenciación: Las células purificadas se sometieron a amplificación del genoma completo (PicoPLEX® Gold) y fueron secuenciadas en una plataforma Illumina NovaSeq™ X (2 × 150 bp), alcanzando una cobertura media de ~15× en las muestras tumorales.
• Bioinformática:
  • CNVs y SVs: Las CNVs somáticas se infirieron utilizando el flujo de trabajo GATK CNV y GISTIC2.0 para picos recurrentes. Las variantes estructurales se identificaron mediante svABA y GaTSV.
  • Firmas mutacionales: Las variantes de nucleótido único (SNV) se identificaron mediante Mutect2. Las firmas mutacionales (SBS y SV) se analizaron utilizando la suite SigProfiler contra los catálogos de referencia COSMIC v3.3.
  • Modelado de riesgo: Se derivó una Puntuación de Inestabilidad Genómica (GIS) a partir de los perfiles de CNV (combinando el porcentaje de gen alterado y la ploidía). Se desarrolló un modelo de riesgo compuesto, TRACER-HL, mediante regresión logística multivariable, integrando la GIS y una puntuación de riesgo de CNV binaria basada en tres loci específicos asociados con la recaída.
• Análisis estadístico: Las asociaciones entre las características genómicas y los resultados clínicos (supervivencia libre de progresión, respuesta al tratamiento) se evaluaron mediante curvas de Kaplan–Meier, pruebas de log-rank y regresión logística. El rendimiento del modelo se validó mediante remuestreo de bootstrap.

Contribuciones clave

1. Perfil genómico de alta pureza: El estudio generó con éxito datos de secuenciación del genoma completo a partir de grupos de células HRS purificadas aisladas de biopsias FFPE, superando el desafío de la baja celularidad tumoral en el LCH.
2. Panorama genómico integral: Los autores mapearon el panorama de CNV, SV y firmas mutacionales somáticas de las células HRS, identificando alteraciones recurrentes consistentes con la biología conocida del Loh (NF-κB, JAK/STAT, evasión inmunitaria) al tiempo que revelaron nuevos marcadores pronósticos.
3. Identificación de determinantes de resultados: El estudio identificó características genómicas específicas vinculadas a una menor supervivencia libre de progresión (PFS) y respuesta al tratamiento, yendo más allá de la genómica descriptiva hacia la utilidad pronóstica.
4. Desarrollo de TRACER-HL: Se creó un nuevo modelo de riesgo compuesto que combina la inestabilidad cromosómica global con picos de CNV focales para estratificar el riesgo de recaída al inicio del tratamiento.

Resultados

• Panorama de CNV: Los genomas de las células HRS exhibieron una marcada inestabilidad cromosómica. Se observaron amplificaciones recurrentes en 19q13.32 (BCL3/SPIB), 2p23.3 (DNMT3A), 2p16.1 (REL/XPO1) y 9p24.1 (JAK2/PD-L1/PD-L2). Las deleciones recurrentes incluyeron 1p36.32 (TNFRSF14).
• CNVs pronósticas: La amplificación en 2p11.2 y las deleciones en 11q25 y 2p13.2 se asociaron significativamente con una menor PFS. Un clasificador compuesto que requiere ≥2 de estos eventos identificó un subgrupo de alto riesgo con una PFS notablemente más corta (log-rank p = 0.00033).
• Firmas mutacionales: El espectro mutacional estuvo dominado por cinco firmas COSMIC: SBS5, SBS8, SBS84, SBS9 y SBS1.
  • Se observó una relación inversa entre la actividad de SBS5 (tipo reloj) y SBS9 (tipo AID).
  • Un perfil "tipo AID" (SBS9-alto/SBS5-bajo) se enriqueció significativamente en el subtipo histológico de esclerosis nodular (p = 0.0111).
  • Las firmas de SV fueron mayoritariamente homogéneas, dominadas por SV2. Una mayor exposición a SV2 se asoció con toxicidad relacionada con el tratamiento de grado 3–4.
• Rendimiento del modelo TRACER-HL: El modelo integrado alcanzó un Área Bajo la Curva (AUC) de 0.81. Con un umbral optimizado, el modelo demostró un Valor Predictivo Negativo (VPN) del 95%, identificando correctamente a 10 de los 12 pacientes con recaída (sensibilidad 83.3%) y clasificando a 38 de los 52 pacientes sin recaída como de bajo riesgo (especificidad 73.1%). La validación por bootstrap confirmó la robustez del modelo (media de AUC 0.80).

Significancia y afirmaciones
El artículo afirma que la secuenciación del genoma completo de las células HRS purificadas revela un panorama genómico complejo dominado por CNV recurrentes que afectan las vías de regulación inmunitaria y de señalización. El estudio sostiene que ciertas CNV focales (amplificación en 2p11.2, deleciones en 11q25/2p13.2) y patrones mutacionales (eje SBS5–SBS9) están asociados con los resultados clínicos y la respuesta al tratamiento.

Crucialmente, los autores afirman que el modelo TRACER-HL proporciona un marco genómico para la estratificación del riesgo basal. Al integrar la inestabilidad cromosómica global con las CNV asociadas con la recaída, el modelo respalda la identificación de un subgrupo de pacientes con una probabilidad muy baja de recaída (VPN alto). Los autores concluyen que estos hallazgos apoyan la evaluación de estrategias de manejo adaptadas al riesgo e informadas genómicamente en el LCH, posicionando la caracterización genómica como un enfoque complementario a las evaluaciones clínicas e de imagen actuales.

El estudio reconoce limitaciones, incluyendo la dependencia de una cohorte prospectiva única con un número limitado de eventos de recaída, la potencial incapacidad de los grupos masivos de células HRS para capturar la heterogeneidad subclonal rara y las restricciones técnicas de la secuenciación del genoma completo de célula única basada en FFPE. En consecuencia, los autores enfatizan que se requiere la validación externa en series independientes más grandes antes de su implementación clínica.