Abstract:Background: Sequencing studies have advanced our understanding of the genomic underpinnings of Multiple Myeloma (MM). While whole genome sequencing (WGS) has emerged as a comprehensive approach for profiling the genomic landscape of MM, WGS has limited sensitivity in defining the intrapatient genomic heterogeneity. Single cell RNA approaches can capture phenotypic heterogeneity but fail to recover most of the key genomic drivers. To overcome these historical limitations, we utilized a newer single cell whole genome direct library preparation (scWGS DLP+; Lacks et al. Cell 2019) sequencing technology able to detect early changes and study clonal population dynamics, providing a more precise understanding of the emergence and clinical impact of genomic events during MM evolution. Methods: scWGS DLP+ was applied to CD138 + Plasma Cells (PC) isolated from bone marrow samples of MM patients. Using newly developed computational approaches, we characterized the Single Nucleotide Variants (SNVs), copy number changes and complex Structural Variants (SVs) at a single cell level and use this information to reconstruct the evolutionary histories of each patient's myeloma. Results: We applied DLP to 7 MM patients obtaining 800-1600 single cell whole genomes per patient. All canonical IgH translocations identified by FISH were recovered in the scWGS data. Single cell copy number heterogeneity was observed in all but 1 patient, including heterogeneity affecting known drivers. The number of subclones ranged from 1 to 15, with an average of 6 subclones detected per patient (≥10 cells, ≥5MB distinct segmental difference). This is in striking contrast with bulk WGS where we are usually able to define no more than one or two subclones. Gain of 1q (3 copies) was observed in all patients but was either rare (1-3%, 3 patients) or clonal (>99%, 4 patients), with two clonal 1q patients exhibiting convergent subclonal 1q amplification (>3 copies). Loss of previously gained chromosomes was observed including loss of ancestrally gained chromosome 15 in hypodiploid Patient 2, and sporadic loss of 1q in patients with ubiquitous 1q gain. For 3 of 7 patients, hypodiploidy (HRD) had arisen before the most recent clonal expansion. HRD clones continued to diversify, accruing multiple additional events on HRD chromosomes including large or focal gains, or deletions subsequent to a gain. Interestingly, one case had 1 small hypodiploid clone co-existing with a completely different major clone with 1q and t(14;20) IgH translocation, suggesting the existence of two distinct clonal plasma cell populations (one malignant and one likely bening). Two of the 7 patients harbored both small (<20 cells) whole genome doubled (WGD) subclones, and small numbers of highly aneuploid likely produced by multi-polar cell division post WGD. We explored two patient's evolutionary histories in detail, highlighting the unique resolution of scWGS DLP+ to detect CNV but also SNV, mutational signatures and SVs. Patient 3 exhibited 4 distinct subclones with distant evolutionary origin as evidenced by the relatively small number of shared SNV. APOBEC signature mutations were detected in all 4 subclones suggesting APOBEC mutagenesis was still ongoing across the cell populations in this patient. In Patient 1, a complex SV on chromosome 17 resulted in loss of most of 17q and an amplification of 17q segments. The heterogeneous structure of chromosome 17 across clones suggests a late onset of this complex event and ongoing instability. Chromosome 16 exhibited a pattern of cell specific arm loss likely resulting from telomere loss. Interestingly, a minor clone showed complex amplification of 16p affecting TNFRSF17, suggesting that chromosome 16 is at an earlier stage of a process of instability that starts with cell specific arm loss, and ends with clonal expansion of a complex amplification SV targeting a gene with fitness benefit. Discussion: Applied to Multiple Myeloma, scWGS shows higher clonal heterogeneity than previously described with multiple clones exhibiting diverse genetic alterations. The findings underscore the genomic complexity and dynamic nature of MM, highlighting the need for sensitive techniques like scWGS to capture the full spectrum and chronological order of genomic events.

Genome‐wide analysis of single nucleotide variants allows for robust and accurate assessment of clonal derivation in cell lines used to produce biologics

Next-generation sequencing of immunoglobulin gene rearrangements for clonality assessment: a technical feasibility study by EuroClonality-NGS

Methods for Estimating the Probability of Clonality in Cell Line Development

Assurance of monoclonality in one round of cloning through cell sorting for single cell deposition coupled with high resolution cell imaging

Assuring Clonality on the Beacon Digital Cell Line Development Platform

Genomic barcoding for clonal diversity monitoring and control in cell-based complex antibody production

Measuring the aggregation of CHO cells prior to single cell cloning allows a more accurate determination of the probability of clonality

Validation of Next Generation Sequencing-Based Clonality Assays for Diagnosis and Minimal Residual Disease Tracking in Lymphoid Malignancies

Single Cell Multi-Omic Correlation of Single Nucleotide Variants, Copy Number Variation and Surface Epitopes for Clonal Profiling of Myeloma

Clonal expansion dictates the efficacy of mitochondrial lineage tracing in single cells

Abstract 3479: Evaluation of subclonal deconvolution pipelines using reference cell-lines and patient plasma samples

A single-step FACS sorting strategy in conjunction with fluorescent vital dye imaging efficiently assures clonality of biopharmaceutical production cell lines.

A unique molecular identifier-based and clonal hematopoiesis-aware approach for accurate mutation calling in cell-free DNA assays.

Exploring the impact of clonal definition on B-cell diversity: implications for the analysis of immune repertoires

Single Cell Whole Genome Measurement of Genomic Heterogeneity and Clonal Evolution in Multiple Myeloma

Identifying cancer cells from calling single-nucleotide variants in scRNA-seq data

Coupling picodroplet microfluidics with plate imaging for the rapid creation of biomanufacturing suitable cell lines with high probability and improved multi-step assurance of monoclonality

Lymphoma Cell Sorting to Enrich Low Degree Infiltration Clonal B-Cells for Diagnostic Delineation of Molecular Genetic Background By NGS

Abstract 2326: Integrating single nucleotide variants (SNVs), copy number alterations (CNAs), and structural variants (SVs) into single-cell clonal lineage inference

Single-cell sequencing of genomic DNA resolves sub-clonal heterogeneity in a melanoma cell line

High‐throughput and automation advances for accelerating single‐cell cloning, monoclonality and early phase clone screening steps in mammalian cell line development for biologics production