Identifying IGH disease clones for MRD monitoring in childhood B-cell acute lymphoblastic leukemia using RNA-Seq

Identifying patient-specific clonal IGH/TCR junctional sequences is critical for minimal residual disease (MRD) monitoring. Conventionally these junctional sequences are identified using laborious Sanger sequencing of excised heteroduplex bands. We found that the IGH is highly expressed in our diagnostic B-cell acute lymphoblastic leukemia (B-ALL) samples using RNA-Seq. Therefore, we used RNA-Seq to identify IGH disease clone sequences in 258 childhood B-ALL samples for MRD monitoring. The amount of background IGH rearrangements uncovered by RNA-Seq followed the Zipf's law with IGH disease clones easily identified as outliers. Four hundred and ninety-seven IGH disease clones (median 2, range 0-7 clones/patient) are identified in 90.3% of patients. High hyperdiploid patients have the most IGH disease clones (median 3) while DUX4 subtype has the least (median 1) due to the rearrangements involving the IGH locus. In all, 90.8% of IGH disease clones found by Sanger sequencing are also identified by RNA-Seq. In addition, RNA-Seq identified 43% more IGH disease clones. In 69 patients lacking sensitive IGH targets, targeted NGS IGH MRD showed high correlation (R = 0.93; P = 1.3 × 10^-14), better relapse prediction than conventional RQ-PCR MRD using non-IGH targets. In conclusion, RNA-Seq can identify patient-specific clonal IGH junctional sequences for MRD monitoring, adding to its usefulness for molecular diagnosis in childhood B-ALL.