Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2014 Aug 19;9(8):e105585.
doi: 10.1371/journal.pone.0105585. eCollection 2014.

A quantitative comparison of single-cell whole genome amplification methods

Affiliations

A quantitative comparison of single-cell whole genome amplification methods

Charles F A de Bourcy et al. PLoS One. .

Abstract

Single-cell sequencing is emerging as an important tool for studies of genomic heterogeneity. Whole genome amplification (WGA) is a key step in single-cell sequencing workflows and a multitude of methods have been introduced. Here, we compare three state-of-the-art methods on both bulk and single-cell samples of E. coli DNA: Multiple Displacement Amplification (MDA), Multiple Annealing and Looping Based Amplification Cycles (MALBAC), and the PicoPLEX single-cell WGA kit (NEB-WGA). We considered the effects of reaction gain on coverage uniformity, error rates and the level of background contamination. We compared the suitability of the different WGA methods for the detection of copy-number variations, for the detection of single-nucleotide polymorphisms and for de-novo genome assembly. No single method performed best across all criteria and significant differences in characteristics were observed; the choice of which amplifier to use will depend strongly on the details of the type of question being asked in any given experiment.

PubMed Disclaimer

Conflict of interest statement

Competing Interests: The authors have declared that no competing interests exist.

Figures

Figure 1
Figure 1. Design of experiments.
Overview of experiments, where n denotes the number of experiments of a given type. In the “E. coli single cells” column, the box that straddles both the “Microfluidic” and the “Tube” fields corresponds to the method of carrying out a first round of amplification in a microfluidic chamber and then a second round of amplification in a test tube. This method will be denoted by “microfluidic+tube” in subsequent figure captions.
Figure 2
Figure 2. Sequence read classification.
(A) Breakdown of read pairs in each experiment according to type of mapping achieved. (B) Breakdown of unmapped reads by organism of origin, expressed as a fraction of the total number of reads.
Figure 3
Figure 3. Amplification bias and uniformity.
(A) Local mapping density from properly mapped reads (at fixed 5x sampling depth), normalized to average 1, as a function of position along the reference sequence, for single-cell MDA (microfluidic in red, microfluidic+tube in orange). (B) Same as panel A, but for single-cell MALBAC and bulk NEB-WGA. (C) Fractional genome coverage from properly mapped read pairs, plotted as a function of gain. Here, each set of properly mapped read pairs was randomly down-sampled to 5x depth. Experiments that did not generate this many properly mapped reads were not included in the figure. (D) Fraction of the genome covered by mapped read pairs when the set of raw read pairs was down-sampled to a fixed depth of 20x, plotted as a function of gain. Filled symbols signify bulk experiments, open symbols single-cell experiments.
Figure 4
Figure 4. CNV resolution.
Size of resolvable duplications (W, minimum width of a sliding window average filter that gives rise to a relative genome mapping density smaller than 2 across all positions in the genome) versus gain. Filled symbols signify bulk experiments, open symbols single-cell experiments.
Figure 5
Figure 5. Combined single-nucleotide error rates.
Main panel: Experimental error rates D versus gain G, for low gains. Here D is the fraction of bases differing from the reference in the mapped reads. Linear fits for D as a function of formula image are also shown: their slope approximately indicates the per-base per-cycle replication error rate. Inset: D versus G over the entire gain range. Filled symbols signify bulk experiments, open symbols single-cell experiments.
Figure 6
Figure 6. De-novo assemblies.
LG50, the minimal number of assembled contigs (≥500 bp) needed to cover 50% of the E. coli reference genome, versus reaction gain (at fixed raw sequencing depth 30x). Assemblies that failed to cover 50% of the reference sequence were symbolically assigned the maximum value that LG50 can take in this scenario ((50%•genome length/500) = 4686). Filled symbols signify bulk experiments, open symbols single-cell experiments.

References

    1. Zong C, Lu S, Chapman AR, Xie XS (2012) Genome-wide detection of single-nucleotide and copy-number variations of a single human cell. Science 338: 1622–1626. - PMC - PubMed
    1. Navin N, Kendall J, Troge J, Andrews P, Rodgers L, et al. (2011) Tumour evolution inferred by single-cell sequencing. Nature 472: 90–94. - PMC - PubMed
    1. Jan M, Snyder TM, Corces-Zimmerman MR, Vyas P, Weissman IL, et al. (2012) Clonal Evolution of Preleukemic Hematopoietic Stem Cells Precedes Human Acute Myeloid Leukemia. Science Translational Medicine 4: 149ra118. - PMC - PubMed
    1. Wells D, Delhanty JD (2000) Comprehensive chromosomal analysis of human preimplantation embryos using whole genome amplification and single cell comparative genomic hybridization. Mol Hum Reprod 6: 1055–1062. - PubMed
    1. McConnell MJ, Lindberg MR, Brennand KJ, Piper JC, Voet T, et al. (2013) Mosaic Copy Number Variation in Human Neurons. Science 342: 632–637. - PMC - PubMed

Publication types

MeSH terms