Single-cell mutation profiles.

<p>Variant profiles across targeted somatic mutations in single-cell samples (sAML bone marrow) in (<b>A</b>) UPN461282, (<b>B</b>) UPN182896, and (<b>C</b>) UPN288033. Rows display positive and negative variant calls color-coded by mutation cluster for each single-cell sample, and columns indicate specific SNVs somatic at sAML diagnosis. Variants are grouped and color-coded by cluster as predicted from sequencing unfractionated material (uppermost track in each panel). Each cell is grouped by the clone it is inferred to represent. Outlier SNVs (purple) were those which could not be confidently clustered based on bulk sequencing. Here, many of these are merged into predicted clusters based upon their presence/absence in single-cell libraries (i.e., harboring the same pattern as well-defined clones). Positions where reference calls were made are colored grey; positions where no call was made (<25× coverage) are colored white. Pairs of variants that always travel in the same state (reference or variant) likely arose in the same clonal expansion. Pairs of variants that are called together in some cells but not others are likely related by linear evolution. Pairs of variants that are mutually exclusive suggest evolutionary branch points, and were rare. This suggested that variants in subclone 5 in UPN461282 (<b>A</b>), and subclone 1 in UPN288033 (<b>C</b>) were divided among additional subclones (now 5A/5B, 2A/2B). See <b><a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1004462#pgen.1004462.s007" target="_blank">Figure S7</a></b> for data presentation with unmodified clone and cluster definitions (derived from bulk sequencing).</p

Depth and distribution of coverage for each sequencing library (n = 56).

<p>(<b>A</b>) Cumulative coverage represented as the proportion of the capture target (y-axis) with read depth greater than or equal to specific coverage thresholds (x-axis). Coverage values are derived from quality-filtered data (de-duplicated, phred-scaled alignment quality ≥10, phred- scaled base quality ≥13). The intersection of each curve with y = 0.5 identifies the median coverage. Higher coverage was obtained for the unsorted samples (median 228×), compared to the single- or two-cell samples (median 28×). (<b>B</b>) Lorenz curve detailing uniformity of coverage as proportion of targeted bases versus proportion of sequenced bases. Dashed line (y = x) represents a perfectly uniform distribution of read depth across the capture target. Libraries prepared from WGA samples (single- and two-cells) exhibit significantly less uniform representation, compared to libraries derived from unfractionated material. See <a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1004462#pgen-1004462-t001" target="_blank"><b>Table 1</b></a> and <b><a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1004462#pgen.1004462.s011" target="_blank">Table S3</a></b> for additional details.</p

Performance of variant calling at germline SNPs.

<p>True Positive (TP): ≥1 non-reference allele called by Affymetrix array, ≥1 non-reference allele called by sequencing.</p><p>True Negative (TN): 0 non-reference alleles called by Affymetrix array, 0 non-reference alleles called by sequencing.</p><p>False Positive (FP): 0 non-reference alleles called by Affymetrix array, ≥1 non-reference allele called by sequencing.</p><p>False Negative (FN): ≥1 non-reference allele called by Affymetrix array, ≥0 non-reference alleles called by sequencing.</p><p>True Positive Rate (TPR): TP/(TP+FN) = sensitivity = power.</p><p>False Positive Rate (FPR): FP/(FP+TN) = 1-specificity.</p><p>False Negative Rate (FNR): FN/(TP+FN) = 1 - sensitivity = type II error.</p

Performance of variant calling.

<p>The specificity (<b>A</b>) and sensitivity (<b>B, C</b>) of three separate variant callers—SAMtools, VarScan2, and GATK—were evaluated by analyzing single-cell variant calls at germline SNPs previously ascertained by Affymetrix SNP arrays <a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1004462#pgen.1004462-Walter2" target="_blank">[14]</a>. As we have defined true positive and true negative, sensitivity is undefined at homozygous reference positions (there are no true positives) and specificity is undefined at heterozygous and homozygous variant positions (there are no true negatives). Sensitivity and specificity were similar among all three callers at homozygous positions, but GATK demonstrated greater sensitivity at heterozygous sites. Variants were called jointly across all single-cell libraries with the GATK Unified Genotyper utility, whereas variants were called independently for each sample using SAMtools and VarScan2. See <a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1004462#pgen-1004462-t002" target="_blank"><b>Table 2</b></a> and <b><a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1004462#pgen.1004462.s014" target="_blank">Table S6</a></b> for additional details. TPR: true positive rate. FPR: false positive rate.</p

Model of tumor clonality predicted from unfractionated samples.

<p>SciClone analysis and cluster assignment of previously validated somatic mutations present in unsorted MDS and sAML bone marrow cells <a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1004462#pgen.1004462-Walter2" target="_blank">[14]</a>, <a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1004462#pgen.1004462-Miller1" target="_blank">[21]</a>. Variant allele fractions of variants at MDS (x-axis) and sAML (y-axis) predict specific tumor substructure and evolution over time. Color-coded mutation clusters are completely non-overlapping between subjects. (<b>A</b>) At sAML diagnosis, UPN461282 was predicted to harbor 10% non-tumor cells in addition to 5 subclones: 1) 6% of cells harboring cluster 1 variants (clone 1), 2) 4% of cells harboring cluster 1 and cluster 2 variants (clone 2), 3) 33% of cells harboring cluster 1, cluster 2, and cluster 3 variants (clone 3), 4) 33% of cells harboring cluster 1, cluster 2, cluster 3, and cluster 4 variants (clone 4), and 5) 14% of cells harboring either cluster 1, cluster 2, cluster 3, and cluster 5 variants or cluster 1, cluster 2, cluster 3, cluster 4, and cluster 5 variants (clone 5). (<b>B</b>) At diagnosis with sAML, UPN182896 was predicted to harbor 35% non-tumor cells in addition to 2 subclones: 1) 13% of cells harboring cluster 1 variants (clone 1), and 2) 52% of cells harboring cluster 1 and cluster 2 variants (clone 2). (<b>C</b>) At diagnosis with sAML, UPN288033 was predicted to harbor 7% non-tumor cells in addition to 2 subclones: 1) 62% of cells harboring cluster 1 variants (clone 1), and 2) 31% of cells harboring cluster 1 and cluster 2 variants (clone 2). (<b>D</b>) Schematic summarizing our initial models of clonal evolution inferred from SciClone analysis—the question mark denotes ambiguity of clone 5 origin in UPN461282.</p