Search

With the increasing number of sequencing projects involving families, quality control tools optimized for family genome sequencing are needed. However, accurately quantifying contamination in a DNA mixture is particularly difficult when genetically related family members are the sources. We developed TrioMix, a maximum likelihood estimation (MLE) framework based on Mendel's law of inheritance, to quantify DNA mixture between family members in genome sequencing data of parent–offspring trios. TrioMix can accurately deconvolute any intrafamilial DNA contamination, including parent–offspring, sibling–sibling, parent–parent, and even multiple familial sources. In addition, TrioMix can be applied to detect genomic abnormalities that deviate from Mendelian inheritance patterns, such as uniparental disomy (UPD) and chimerism. A genome-wide depth and variant allele frequency plot generated by TrioMix facilitates tracing the origin of Mendelian inheritance deviations. We showed that TrioMix could accurately deconvolute genomes in both simulated and real data sets.

SummaryMyeloid malignancies, including acute myeloid leukemia (AML), arise from the proliferation and expansion of hematopoietic stem and progenitor cells which acquire somatic mutations. Bulk molecular profiling studies on patient samples have suggested that somatic mutations are obtained in a step-wise fashion, where mutant genes with high variant allele frequencies (VAFs) are proposed to occur early in disease development and mutations with lower VAFs are thought to be acquired later in disease progression1–3. Although bulk sequencing informs leukemia biology and prognostication, it cannot distinguish which mutations occur in the same clone(s), accurately measure clonal complexity and clone size, or offer definitive evidence of mutational order. To elucidate the clonal framework of myeloid malignancies, we performed single cell mutational profiling on 146 samples from 123 patients. We found AML is most commonly comprised of a small number of dominant clones, which in many cases harbor co-occurring mutations in epigenetic regulators. Conversely, mutations in signaling genes often occur more than once in distinct subclones consistent with increasing clonal diversity. We also used these data to map the clonal trajectory of each patient and found that specific mutation combinations (FLT3-ITD + NPM1c) synergize to promote clonal expansion and dominance. We combined cell surface protein expression with single cell mutational analysis to map somatic genotype and clonal architecture with immunophenotype. Our studies of clonal architecture at a single cell level provide novel insights into the pathogenesis of myeloid transformation and how clonal complexity contributes to disease progression.

Clonal hematopoiesis (CH) is frequent in cancer patients and associated with increased risk of therapy related myeloid neoplasms (tMN). To define the relationship between CH, oncologic therapy, and tMN progression, we studied 24,439 cancer patients. We show that previously treated patients have increased rates of CH, with enrichment of mutations in DNA Damage Response (DDR) genes (TP53, PPM1D, CHEK2). Exposure to radiation, platinum and topoisomerase II inhibitors have the strongest association with CH with evidence of dose-dependence and gene-treatment interactions. We validate these associations in serial sampling from 525 patients and show that exposure to cytotoxic and radiation therapy imparts a selective advantage specifically in hematopoietic cells with DDR mutations. In patients who progressed to tMN, the clone at CH demarcated the dominant clone at tMN diagnosis. CH mutational features predict risk of therapy-related myeloid neoplasm in solid tumor patients with clinical implications for early detection and treatment decisions.