Here we present a protocol for measuring absolute mitochondrial (mt)DNA copy number and mtDNA deletion heteroplasmy levels in single cells.
The mammalian mitochondrial (mt)DNA is a small, circular, double-stranded, intra-mitochondrial DNA molecule, encoding 13 subunits of the electron transport chain. Unlike the diploid nuclear genome, most cells contain many more copies of mtDNA, ranging from less than 100 to over 200,000 copies depending on cell type. MtDNA copy number is increasingly used as a biomarker for a number of age-related degenerative conditions and diseases, and thus, accurate measurement of the mtDNA copy number is becoming a key tool in both research and diagnostic settings. Mutations in the mtDNA, often occurring as single nucleotide polymorphisms (SNPs) or deletions, can either exist in all copies of the mtDNA within the cell (termed homoplasmy) or as a mixture of mutated and WT mtDNA copies (termed heteroplasmy). Heteroplasmic mtDNA mutations are a major cause of clinical mitochondrial pathology, either in rare diseases or in a growing number of common late-onset diseases such as Parkinson’s disease. Determining the level of heteroplasmy present in cells is a critical step in the diagnosis of rare mitochondrial diseases and in research aimed at understanding common late-onset disorders where mitochondria may play a role. MtDNA copy number and heteroplasmy have traditionally been measured by quantitative (q)PCR-based assays or deep sequencing. However, the recent introduction of ddPCR technology has provided an alternative method for measuring both parameters. It offers several advantages over existing methods, including the ability to measure absolute mtDNA copy number and sufficient sensitivity to make accurate measurements from single cells even at low copy numbers. Presented here is a detailed protocol describing the measurement of mtDNA copy number in single cells using ddPCR, referred to as droplet generation PCR henceforth, with the option for simultaneous measurement of heteroplasmy in cells with mtDNA deletions. The possibility of expanding this method to measure heteroplasmy in cells with mtDNA SNPs is also discussed.
Mammalian mitochondrial (mt)DNA is a small (approx. 16.5 Kb), circular DNA genome residing in the mitochondrial matrix that encodes 37 genes, comprising two rRNAs, 22 tRNAs, and 13 protein-coding genes1. Unlike the nuclear genome, which contains one (haploid) or two (diploid) copies of each gene per cell, mtDNA is present in multiple copies in the mitochondria of each cell, ranging from tens of copies (e.g., mature spermatocytes) to hundreds of thousands of copies (e.g., oocytes)2,3. A consequence of this multi-copy nature is that mutations in the mtDNA genome, which may exist as single-nucleotide polymorphisms (SNPs), deletions, or duplications, can be present at varying levels in any given cell, making up anywhere from 0% to 100% of the cell’s total mtDNA population. The existence of wild-type and mutant mtDNA genomes in the same cell is termed heteroplasmy, and pathogenic heteroplasmic mtDNA mutations are a major cause of mitochondrial disease, with several common neurological syndromes linked to underlying heteroplasmic mtDNA mutations4.
Two key parameters that contribute to the likelihood of a heteroplasmic mtDNA mutation causing clinical disease are the heteroplasmy level and the mtDNA copy number. Many heteroplasmic mutations display a threshold effect, with biochemical and clinical phenotypes only becoming apparent above a certain heteroplasmy level, typically around 80%5, and subsequently worsening as heteroplasmy increases further4. However, it is also important to consider the number of copies of mtDNA present in the cell, as this will influence the number of wild-type (i.e., ‘healthy’) mtDNA genomes that are present at a given heteroplasmy level. Studies in mitochondrial disease patients have highlighted the importance of this interplay between heteroplasmy and copy number5,6, and Filograna et al. recently reported an increase in mtDNA copy number alleviating symptoms despite unchanged heteroplasmy in a mouse model of mitochondrial disease7.
Whilst much has been done in recent years to improve the understanding of the pathogenesis and transmission of diseases caused by mtDNA heteroplasmy, most of this work has been conducted at the level of tissues rather than cells, comparing mean tissue heteroplasmy levels and copy number measurements acquired from bulk tissue biopsies and blood samples. Some well-established techniques, such as laser capture microdissection, allow for such measurements at the cellular level8,9; however, the recent explosion of high-throughput single-cell analysis methods, or so-called “Single-cell Omics”10, has created a requirement for methods that can accurately measure these crucial mtDNA parameters at the single-cell level.
The droplet generation PCR method takes advantage of recent advances in microfluidic technology to improve upon existing methods of quantitating unknown DNA concentrations via PCR amplification of specific target amplicons in the sample DNA11. Unlike qPCR, where the sample DNA is amplified in a single reaction, with the relative rate of accumulation of PCR product acting as the readout, this method splits the initial sample into thousands of individual droplets, thus compartmentalizing the sample DNA molecules into spatially separate reactions12. The subsequent PCR reaction proceeds in each individual droplet, with the PCR product only accumulating in droplets that contain the target DNA. The result of this reaction is a digital output in the form of a pool of droplets that either contain a copy of the target DNA and amplified PCR product, or do not contain target DNA. Using fluorescent DNA probes or a double-stranded (ds)DNA binding dye, the droplets containing amplified product can be counted, and the ratio of ‘positive’ to ‘negative’ droplets can be used to calculate the absolute number of DNA copies that were present in the initial sample. This absolute measurement, as opposed to the relative measurement acquired from a qPCR reaction, is the key factor that allows this methodology to be used for accurate measurement of mtDNA copy number and heteroplasmy in single cells11, and several recent studies have already utilized droplet generation PCR technology for this purpose13,14. This article presents a method for measuring mtDNA copy number and deletion heteroplasmy in single cells from both human and mouse tissue.
All experiments followed the ARRIVE guidelines and were approved by the University of Cambridge Animal Welfare Ethical Review Body (AWERB).
NOTE: All sample preparation steps prior to droplet generation must be performed in a clean pre-PCR work area, ideally in a UV-sterilized cabinet where possible. The protocol described here uses specific droplet generation PCR equipment (see Table of Materials), and while the general method should be applicable to other systems, it is recommended to consult the manufacturer's guidelines regarding primer/probe concentrations, PCR cycling conditions, etc., as they may differ from those described here. Cell lines/primary cells used in this study were as follows: Human HeLa cells (commercially obtained), human HEK 293T cells (commercially obtained), primary human dermal fibroblast cells (obtained from the Newcastle Biobank), human cybrids (WT & ΔH2.1 deletion15, obtained from C. Moraes, University of Miami), mouse embryonic fibroblasts (immortalized, from C57Bl/6 mice, obtained from J. Stewart, Newcastle University), mouse primordial germ cells (obtained from C57Bl/6 mouse embryos), and mouse MII oocytes (obtained from adult female C57Bl/6 mice). All cultured cells were maintained in High Glucose (4.5g/L) DMEM supplemented with 10% fetal bovine serum at 37 °C with 5% CO2. Primary mouse cells used in this study were isolated from animals kept according to the Animal (Scientific Procedures) Act 1986 under Home Office Project License P6C97520A.
1. Design and synthesize primer & probe sets targeting DNA sequences of interest
NOTE: Droplet generation PCR can also be performed using a dsDNA binding dye in place of the amplicon-specific probes.
2. Isolation of DNA from single cells
NOTE: This method can also be used for small bulk samples of up to 100 cells.
3. Preparation of samples
NOTE: Ensure that technical replicates of samples and non-template controls (NTCs) are included on each assay plate to ensure the accuracy of results. The dynamic range of the droplet generation PCR assay is up to 120,000 copies of the target amplicon per reaction. For single-cell or small bulk-cell sample lysates, dilution is unlikely to be necessary, and the lysate mixture can be input directly to the reaction for the measurement of mtDNA copy number (e.g., a lysate sample of 10 µL can be split and run in triplicate with 3 µL input directly into each assay). However, using undiluted lysate from cells with very high copy mtDNA numbers (e.g., oocytes) or containing larger numbers of cells (e.g., 50-100) may exceed this dynamic range. In such cases, an initial serial dilution is required to identify an appropriate dilution factor that avoids saturation of the assay for that specific cell type/cell number (see Representative Results for further details).
4. Generation of droplets
NOTE: Refer to Figure 1A for the schematic of the droplet generator instrument deck referenced in this section.
5. PCR
6. Droplet reading
7. Analysis of results
NOTE: The droplet reader measures the fluorescence intensity in the FAM and HEX channel for each droplet in a sample. In a successful assay, droplets fall into one of two categories for each probe: Negative (meaning the target was not present in the droplet) or Positive (meaning the target was present in the droplet). Prior to analyzing samples, ensure that each well contains >10,000 droplets and has two clearly separate populations of low fluorescence (negative) and high fluorescence (positive) droplets in each channel (Figure 2A).
Following droplet generation, a clear layer of opaque droplets is visible floating on top of the oil phase in each well (Figure 1B). Droplet formation can be adversely affected by the presence of detergents in the input lysate when performing experiments on single cells. Using the lysis protocol described in 2.1.2., droplet yields above the recommended level of 10,000 are routinely achieved, despite the presence of a small residual amount of TWEEN-20 in the final sample (Figure 1C). However, with other lysis buffers, such as buffer RLT Plus, a dramatic reduction is seen in the number of droplets formed (Figure 1C), which is likely to adversely affect the accuracy of the copy number measurement.
Following PCR and droplet reading, clearly separated populations of positive and negative droplets will be seen in each of the channels used for each sample well (Figure 2A), with NTC wells containing only negative droplets. Potential issues that may be encountered when analyzing results include:
The lack of clear separation between positive and negative populations, and the presence of droplets in the gap between the two populations, termed 'rain' (Figure 2B). This can occur if a non-optimal annealing/extension temperature is used, especially in assays where the difference in fluorescence between positive and negative droplets is relatively small. Running a gradient of annealing/extension temperatures may help identify an optimal temperature to use. Alternatively, primers/probes may need to be redesigned to improve assay performance.
A 'double band' of droplets in the negative population when using the '1D Amplitude' view. In some cases, a double band of negative droplets may be seen in the HEX channel (Figure 2Bii); this is the result of differing levels of background fluorescence between the double-negative and single-negative (i.e., FAM-positive, HEX-negative) droplet populations, as seen in (Figure 2Biii), which shows FAM single-positive population (top left quadrant) with lower mean fluorescence in the HEX channel than the double-negative population (bottom left quadrant). This fluorescence artifact does not affect the results of the assay, and use of the '2D Amplitude' allows the four droplet populations to be clearly identified.
Absence of positive droplets (Figure 2C). For validated assays, this is most likely due to an error during sample preparation (e.g., incomplete cell lysis, incorrect primers/probe added to the master mix, failure of single-cell sorting, etc.).
Absence of negative droplets (Figure 2D). This occurs when the assay is saturated with a high target DNA input, resulting in every droplet containing one or more target sequences. This is most likely to be seen when measuring high copy number samples (e.g., oocytes) and can be resolved by diluting the lysate prior to step 3.5. A serial dilution of the input DNA (Figure 2D) can be used to identify a suitable dilution factor when this problem is encountered.
When using two independent mtDNA probes, the integrity of the sample DNA can be assessed by selecting the '2D Amplitude' view: Intact mtDNA molecules carry both target sequences; therefore, droplets containing complete mtDNA genomes will be positive in both channels (Figure 3A). When analyzing samples carrying mtDNA deletions, a mixture of double-positive droplets and droplets positive only for the non-deleted target sequence will be present, even in non-degraded samples (Figure 3B). When mtDNA is degraded, some target sequences are likely to exist on separate DNA fragments, and therefore droplets are more likely to be positive for just one probe, and fewer droplets will be positive for both (Figure 3C). These differences are accounted for during data analysis, as concentrations are always estimated independently for each target using the total population of all positive droplets in that fluorescence channel, and then either averaged to estimate copy number or used to calculate heteroplasmy (see steps 7.3.-7.4.).
This method can be used to measure the copy number in a wide variety of cell types from both mouse and human tissues, including laboratory cell lines such as HeLa, HEK 293T, and immortalized mouse embryonic fibroblasts (MEFs) and primary cells such as human dermal fibroblasts and mouse primordial germ cells (PGCs) (Figure 4A). In the somatic cell types tested, a range of copy numbers is seen from a few hundred (e.g., dermal fibroblasts) to several thousand (e.g., HeLa cells). To quantify assay precision, mature mouse oocytes, which have a very high copy number (>150,000 copies per cell)2, were lysed, and the resulting lysate was split equally across 22 separate reactions per cell. After running the droplet generation PCR protocol, very little difference was seen in the number of copies detected in each reaction, with two representative cells (Figure 4B) yielding a total of 229,647 copies (mean 10,439 per well, standard deviation 380) and 330,121 copies (mean 15,006 per well, standard deviation 340) per oocyte. The coefficients of variation (CV) for these two cells were 3.64% and 2.26%, respectively. To test reproducibility across biological specimens, PGCs were harvested from the gonads of four embryonic mice at day 13.5 post-conception, and the resulting single-cell copy number of these cell populations was compared across the different embryos (Figure 4C). Whilst some biological variation is to be expected, the mean copy number and standard deviation remained comparable across the four embryos tested, and these results were in agreement with copy number data previously reported in this cell type at a similar stage of embryonic development21.
In addition to measuring copy number, this method can also be utilized to measure heteroplasmy level in cybrid cells carrying a heteroplasmic mtDNA deletion. This requires the use of two mtDNA probes, one falling inside the deleted region of the mtDNA and one falling outside (Figure 5A). After analysis, the level of heteroplasmy can be found by comparing the concentration of each probe to calculate the proportion of copies that contain both targets (i.e., WT mtDNA genomes) and the proportion that only contain the probe which falls outside the deleted region (i.e., deleted mtDNA genomes) as described in step 7.4. Since the target located outside the deleted region is present in all mtDNA copies, this can also be used to calculate the overall copy number of the cell. In cells with high deletion heteroplasmy, relatively few copies of the probe inside the deleted region are detected (Figure 5Bi), while in cells with low deletion heteroplasmy, the concentrations of the two probes are much more similar (Figure 5Bii). Use of 'in-out' probe assays to measure deletion heteroplasmy by droplet generation PCR in bulk-extracted DNA samples has shown improved precision over qPCR-based methods, but lacks sensitivity when dealing with very low heteroplasmy levels (<5%)22, and therefore is not recommended in cells with very small deletion burdens. The primer/probe sets for both mouse and human mtDNA described in Table 1 comprise one target within the mtDNA major arc and one target outside, and therefore are suitable for most common mtDNA deletions23,24.
Species | Target | Primer/Probe Name | Sequence | 5’ Probe | 3’ Quencher | ||
Human | ND4 | huND4_F Primer | AGTGCGATGAGTAGGGGAAGG | ||||
huND4_R Primer | ACCTTGGCTATCATCACCCGAT | ||||||
huND4_FAM Probe | CAACCAGCCAGAACGCCTGAACGCA | FAM | BHQ 2 | ||||
ND1 | huND1_F Primer | GGGTTCATAGTAGAAGAGCGATGG | |||||
huND1_R Primer | ACGCCATAAAACTCTTCACCAAAG | ||||||
huND1_HEX Probe | ACCCGCCACATCTACCATCACCCTC | HEX | BHQ 1 | ||||
Mouse | ND1 | musND1_F Primer | GAGCCTCAAACTCCAAATACTCACT | ||||
musND1_R Primer | GAACTGATAAAAGGATAATA GCTATGGTTACTTCA |
||||||
musND1_FAM Probe | CCGTAGCCCAAACAAT | FAM | BHQ 1 | ||||
COX3 | musCOX3_F Primer | CCTCGTACCAACACATGATCTAGG | |||||
musCOX3_R Primer | AGTGGGACTTCTAGAGGGTTAAGTG | ||||||
musCOX3_HEX Probe | ACCTCCAACAGGAATTTCA | HEX | BHQ 1 |
Table 1: Validated droplet generation PCR primer and probe sets for use in human and mouse cells. Primer and probe sets for previously validated assays amplifying two independent mtDNA targets for human16,17 and mouse cells18.
Component | Volume per sample | Final Concentration |
2X ddPCR Supermix for Probes (no dUTP) | 11 µl | 1X |
Target 1 Primer F (20 µM) | 1 µl | 900 nM |
Target 1 Primer R (20 µM) | 1 µl | 900 nM |
Target 2 Primer F (20 µM) | 1 µl | 900 nM |
Target 2 Primer R (20 µM) | 1 µl | 900 nM |
Target 1 Probe (FAM, 10 µM) | 0.55 µl | 250 nM |
Target 2 Probe (HEX, 10 µM) | 0.55 µl | 250 nM |
Sample DNA (added in step 3.4) | Variable | Variable |
Nuclease-free water | Variable | |
Total Volume | 22 µl |
Table 2: Droplet generation PCR reaction mix recipe. Reaction mix required for a single-well droplet generation PCR assay. Recommended primer/probe sequences can be found in Table 1; the recipe can be scaled up when running multiple samples.
Step | Temp (°C) | Time | No. Cycles |
Enzyme Activation | 95 | 10 min | 1 |
Denaturation | 94 | 30 sec | 40 |
Annealing/Extension | 58 | 1 min | |
Enzyme Deactivation | 98 | 10 min | 1 |
Hold | 4 | Infinite | 1 |
Table 3: Droplet generation PCR thermal cycling protocol. Cycling protocol for amplification of target DNA following droplet generation.
Figure 1: Droplet generation. (A) Schematic layout of the droplet generation PCR droplet generator instrument deck, positions as follows: A = automated liquid handling module, B = droplet generator cartridge (x3), C = tip waste bucket, D = filter tip box (x2), E = droplet generator oil holder, F = sample plate, G = collection plate (housed in cool block). (B) Representative image of droplet emulsion floating on top of the oil phase immediately after droplet generation. (C) Impact of different lysis buffers on droplet number; with 1% TWEEN-20 there is a small decrease in droplet number, whilst RLT Plus buffer results in an almost threefold reduction in droplet numbers, well below the recommended number of 10,000 (results shown as mean of three technical repeats per condition with standard deviations). (D) Screenshot of the software Well Information tab showing completed sample information applied to well A1 of a 96-well template. Please click here to view a larger version of this figure.
Figure 2: Representative droplet generation PCR assay analyses showing expected results and common issues. (A) Dual mtDNA probe assay results from three HeLa single-cell lysates showing clearly separated positive and negative droplet populations in the (i) FAM channel, (ii) HEX channel, and (iii) combined 2D amplitude view. (B) Assay results from four bulk HeLa 20-cell samples showing 'rain' between the positive and negative populations in the (i) FAM channel, (ii) poor separation between positive and negative droplets in the HEX channel, and (iii) a double droplet band in the negative droplets of the HEX channel (ii) due to different background fluorescence levels in the HEX-negative populations, as seen in 2D amplitude view. (C) An example of a WT cybrid single-cell lysate with no positive droplets (i) in either fluorescence channel, with a successful WT cybrid assay from the same experiment (ii) for comparison. (D) Representative results from a 10-fold serial dilution of PCR-amplified human ND1 target amplicon with known concentration, showing saturation of the droplet generation PCR above the assay upper limit of 120,000 copies. RFU = Relative Fluorescence Units. Please click here to view a larger version of this figure.
Figure 3: Assessing mtDNA integrity using the '2D Amplitude' view. (A) Representative droplet generation PCR results from a WT cybrid single-cell lysate containing mostly intact mtDNA genomes, indicated by the presence of mainly FAM and HEX double-positive droplets. (B) Droplet generation PCR results from a heteroplasmic deletion cybrid single-cell lysate, showing a mixture of HEX single-positive droplets (containing deleted mtDNA molecules carrying only the HEX target) and FAM and HEX double-positive droplets (containing complete mtDNA molecules carrying both targets), with very few FAM single-positive droplets, indicating no major mtDNA degradation. (C) Droplet generation PCR results from a mouse oocyte single-cell lysate containing a mixture of FAM single-positive, HEX single-positive, and FAM and HEX double-positive droplets, indicating that some mtDNA degradation is likely to have occurred in this sample. RFU = Relative Fluorescence Units. Please click here to view a larger version of this figure.
Figure 4: Single-cell copy number measurements reveal variation between cell types but show robust intra- and inter-experimental repeatability. (A) Representative single-cell copy number measurements from a range of cell lines from both human (squares) and mouse (triangles) tissues (n = 12 cells per tissue). (B) Serial measurement of copy number in lysate from two individual mouse oocytes; total lysate volume was split across 22 separate droplet generation PCR reactions; each point on the box and whisker plots represents the copy number measured in a single well. Mean, standard deviation (SD), and coefficient of variation (CV) are reported for each oocyte. (C) Single-cell copy number of primary mouse primordial germ cells collected from four separate embryos at 13.5 days post-fertilization (n = 94 single cells per embryo). Please click here to view a larger version of this figure.
Figure 5: Measuring single-cell heteroplasmy in a deletion cybrid line. (A) Schematic showing the region of mtDNA deleted in the cybrid line used for this study and the location of the ND1 and ND4 mtDNA probes. (B) Representative results from (i) a high heteroplasmy single-cell lysate and (ii) a low heteroplasmy single-cell lysate showing the differing ratios of ND4:ND1 probe and the resulting copy number and heteroplasmy measurements. RFU = Relative Fluorescence Units. Please click here to view a larger version of this figure.
The protocol described here is applicable across a wide range of cell types and species in addition to those discussed above, although careful optimization of new assay designs will be key to ensure that accuracy and repeatability of the method are maintained when moving away from previously validated primer/probe combinations. When working with single cells it is vital to ensure that sample collection is performed as accurately as possible (e.g., using stringent single-cell parameters when sorting cells by FACS) to ensure that the results obtained truly reflect those of individual cells, and also that all preparatory work is carried out in a clean pre-PCR area to reduce the risk of contamination occurring. Efficient cell lysis is also critical to ensure that mtDNA copies can be successfully compartmentalized into separate droplets; however, this needs to be balanced against the risk of disrupting droplet formation when using higher concentrations of detergent.
Due to the assay design and lysis steps being the most critical in the protocol, these are also the steps most likely to require troubleshooting. Primer/probe sets that do not initially produce well-separated populations of positive and negative droplets may show improved performance if the annealing/extension temperature of the PCR protocol is optimized. The suggested temperature of 58 °C given in step 5.1. has been validated for the primer/probe sets in (Table 1); however, the optimal temperature can vary from 55-65 °C, and a gradient PCR across this range can be run to find the most suitable temperature for new primer/probe sets. If this is unsuccessful, then it may be necessary to redesign primers and/or probes to improve performance, referring to manufacturer guidelines on assay design (see Table of Materials).
If different lysis conditions to those described in step 2.1.2. are required, then the impact of the alternative lysis buffer on droplet generation must be assessed. It may be possible to improve droplet generation by reducing the concentration of detergent in the lysis buffer or diluting the lysate in a larger volume of nuclease-free water, although we have found that some buffers (e.g., RLT Plus buffer) still have a significant effect on droplet formation even at very low concentrations. One possible solution to this problem is to perform a post-lysis cleanup of the DNA, for example using SPRI beads to bind the DNA and allow the detergents to be washed away. Whilst SPRI beads are used to purify both genomic DNA25 and mtDNA26 from single cell lysates, any loss of mtDNA genomes during this process will affect the subsequent copy number measurement, and there is currently no data that quantifies this potential source of error when using commercially available SPRI bead products on single-cell lysates.
Traditionally, mtDNA copy number and deletion heteroplasmy has been measured using qPCR assays27. Whilst the underlying concept of the qPCR approach is very similar to the droplet generation PCR technique described here, one major drawback of using qPCR is its inability to directly measure absolute copy number11. For each sample, qPCR measures a cycle threshold (ct)value – the number of cycles of PCR required to reach a certain amount of product from the initial input DNA. Since this threshold is arbitrary, it is necessary to generate a standard curve from inputs containing a range of known target concentrations before absolute copy numbers can be assigned to samples. At the bulk DNA level, this process is very accurate and comparable to droplet generation PCR in its ability to measure mtDNA copy number. However, at very low copy numbers, such as those typically found in single cells, it becomes increasingly difficult to accurately generate standard curves for qPCR, and therefore, harder to obtain accurate measures of absolute copy number11. It should be noted that droplet generation PCR is currently more expensive to perform than qPCR, with consumables and reagents costing roughly 2.5x as much on a per-sample basis, in addition to higher initial equipment costs. This means that the cost of droplet generation PCR may be prohibitive for certain high-throughput applications, in which case qPCR still represents a viable alternative.
Since droplet generation PCR directly measures the number of copies of a given target present in a sample, there is no requirement for a standard curve when measuring absolute copy number, and the droplet generation PCR method is therefore extremely sensitive even with very small amounts of input DNA11,12. Given the emerging importance of single-cell technologies in the field of mitochondrial biology, this droplet generation PCR-based approach represents a key advance over previous methodologies and allows for continued measurement of key mitochondrial parameters in the single-cell setting. It is also possible to use droplet generation PCR to measure mtDNA copy number in DNA extracted from bulk tissue samples; this requires measurement of nuclear DNA copy number in addition to mtDNA copy number to allow for normalization to cell number. Therefore, the inclusion of a primer/probe set targeting a nuclear gene is required alongside a primer/probe set targeting mtDNA. As indicated above, at the bulk tissue level, droplet generation PCR and standard qPCR methods are of comparable accuracy, and therefore, droplet generation PCR does not offer significant advantages over existing technologies in this setting.
One major limitation of the method presented here is the inability to measure heteroplasmy levels in cells carrying heteroplasmic SNPs rather than mtDNA deletions, thus requiring an additional assay (e.g., pyrosequencing or qPCR) to measure both copy number and heteroplasmy from a single cell. Successful detection of SNPs demonstrating Allele-specific expression (ASE) using this platform has been reported28, and a recent report described the use of this method to measure heteroplasmy of the common m.3243A>G MELAS mutation29. Therefore, whilst it seems that this limitation can be overcome using carefully designed probe pairs, it requires a separate primer/probe set to be designed/validated for each individual SNP, and the design parameters are much more restricted by the need for probes to target a specific base in the mtDNA sequence. As such, incorporating SNP detection into the droplet generation PCR protocol described here has not yet been attempted.
It is important to note that droplet-based digital PCR platforms, such as the type used in this study, are not the only technologies capable of absolute quantification of DNA targets, and commercial products utilizing alternative methods such as well plate-based partitioning are also available30. Furthermore, technologies using alternatives to PCR, such as loop-mediated isothermal amplification (LAMP)31 and recombinase polymerase amplification (RPA)32, may prove to be more suitable than droplet generation PCR for use in clinical settings.
Overall, this protocol presents a highly sensitive method for measuring mtDNA copy number and deletion heteroplasmy in single cells, which can be applied to a wide range of different cell types and species.
The authors have nothing to disclose.
Thank you to Dr. L Bozhilova for advice on the statistical analysis of droplet generation PCR data. Thank you to Dr. H Zhang for providing the oocytes used to generate data in Figure 3C and Figure 4B. This work was carried out by SPB at the Medical Research Council Mitochondrial Biology Unit (MC_UU_00015/9), University of Cambridge, and funded by a Wellcome Trust Principal Research Fellowship held by PFC (212219/Z/18/Z).
50% Tween-20 solution | Novex | 3005 | |
Automated droplet-generating oil | Bio Rad | 1864110 | Commercial oil formulation used to generate the oil/droplet emulsion (used in Protocol Step 4.1) |
C1000 PCR machine with deep-well block | Bio Rad | 1851197 | PCR thermocycler equipped with a deep-well heating block, used for cell lysis (Protocol Step 2.1.2.) and PCR cycling (Protocol Step 5) |
Collection plate cooling block | Bio Rad | 12002819 | Cooling block that keeps samples chilled during droplet generation (used in Protocol step 4.3) |
ddPCR 96-well plates | Bio Rad | 12001925 | 96-well plates pipet tips designed for use in the QX200 AutoDG droplet generator, used for sample preparation (Protocol step 3.4) and droplet collection (Protocol step 4.3) |
ddPCR droplet reader oil | Bio Rad | 1863004 | Commercial oil formulation used by the droplet reader (used in Protocol step 6.1) |
ddPCR Supermix for Probes (no dUTP) | Bio Rad | 1863023 | Commercial supermix for use in ddPCR experiments utilising probes (used in Protocol Step 3.3) |
DG32 automated droplet generator cartridges | Bio Rad | 1864108 | Microfluidic cartridges used in the QX200 AutoDG droplet generator to generate the oil/droplet emulsion (used in Protocol Step 4.3) |
Fetal bovine serum | Gibco | 10270-106 | Qualified fetal bovine serum |
Foil plate covers | Bio Rad | 1814040 | Foil plate covers used to seal droplet collection plates after droplet generation (used in Protocol step 4.6) |
HEK 293T cells | Takara | 632180 | Commercial subclone of the transformed human embryonic kidney cell line, HEK 293, expressing the SV40 Large-T antigen |
HeLa cells | ECACC | 93021013 | Human cervix epitheloid carcinoma cells |
High glucose DMEM | Gibco | 13345364 | 4.5g/L D-Glucose, with L-glutamine and sodium pyruvate |
Human cybrids | University of Miami | ||
Mouse embryonic fibroblasts | Newcastle University | Immortalized from C57Bl/6 mice | |
Nuclease-free water | Ambion | AM9937 | |
PCR plate seals | Pierce | SP-0027 | Clear adhesive plate seals, only used pre-droplet generation (foil seal must be used in step 4.6) |
Pipet Tip Waste Bins | Bio Rad | 1864125 | Disposable collection bin used to collect discarded tips in the QX200 AutoDG droplet generator (used in Protocol step 4.3) |
Pipet tips for AutoDG system | Bio Rad | 1864120 | Filtered pipet tips designed for use in the QX200 AutoDG droplet generator (used in Protocol step 4.3) |
Primary human dermal fibroblast cells | Newcastle Biobank | ||
Primers/Probes | IDT | N/A | Exact primer/probe sequences will be assay dependent. Primers and probes used in this study are given in Table 1 |
Proteinase K 20 mg/mL solution | Ambion | AM2546 | |
PX1 PCR plate sealer | Bio Rad | 1814000 | Applies foil seals to ddPCR sample plates after droplet generation (used in Protocol Step 4.6) |
QX Manager software | Bio Rad | 12012172 | Droplet reader set up & analysis software (used in Protocol Steps 6 & 7) |
QX200 AutoDG droplet generator | Bio Rad | 1864101 | Automated microfluidic droplet generator (used in Protocol Step 4) |
QX200 droplet reader | Bio Rad | 1864003 | Droplet reader (used in Protocol Step 6) |
Trizma pre-set crystals pH 8.3 | Sigma | T8943-100G |