Name: Measuring Single-Cell Mitochondrial DNA Copy Number and Heteroplasmy Using Digital Droplet Polymerase Chain Reaction (Video) | JoVE Journal
Uploaded: 2022-07-12T13:30:00.000+00:00
Duration: 9 min 15 s
Description: Learn to measure mitochondrial DNA copy number and deletion heteroplasmy in single cells using digital droplet PCR. This protocol enables precise analysis for mitochondrial research.

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

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&#160;3C and Figure&#160;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&#160;a Wellcome Trust Principal Research Fellowship held by PFC (212219/Z/18/Z).

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 &amp; Δ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 &amp; 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.
<ol>
	<li>Design the primer and probe sequences according to the guidelines given by the system manufacturer. Validated primer and probe sequences suitable for measuring single-cell mtDNA copy number/deletion heteroplasmy in both human16,17 and mouse18 cells are provided in Table 1. When selecting alternative target sequences, consider the following points.
	<ol>
		<li>Multiplex two primer/probe sets in a single PCR assay but ensure that the two assays utilize one FAM-labeled probe and one HEX-labeled probe so that the target amplicons can be differentiated by the droplet reader. This protocol presents a duplex probe assay. With careful experimental design, multiplexing can be increased up to four targets, and alternative platforms can achieve higher-dimensional multiplexing.</li>
		<li>When multiplexing primer/probe sets targeting separate mtDNA sequences, ensure that the amplicons generated by the two primer sets do not overlap.</li>
		<li>When measuring heteroplasmy in samples carrying mtDNA deletions, ensure that one target amplicon falls within the expected deleted region and the other falls outside of the expected deleted region. 
		NOTE: Optimal PCR cycling conditions for alternative assays may differ from those quoted in Step 5.1; see the Discussion for further details on assay optimization.</li>
	</ol>
	</li>
</ol>
2. Isolation of DNA from single cells
NOTE: This method can also be used for small bulk samples of up to 100 cells.
<ol>
	<li>Collect single cells using an appropriate method, such as Fluorescence-Activated Cell Sorting (FACS)19 or laser capture microdissection20, in as small a volume as possible into a suitable receptacle (e.g., 96-well plate). Use of a cell viability dye prior to or during single-cell isolation will minimize the likelihood of isolating dead cells. Sort the single cells directly into lysis buffer if desired and store at -80 °C (up to 6 months) prior to analysis.</li>
	<li>Lyse the cells in a small volume (&lt;10 µL) of a suitable lysis buffer (the buffer used in this study is described in step 2.2.1.). 
	NOTE: All data obtained from cells in this study are from single-cell lysates or pools of 20 cells (sample sizes are stipulated in the figure legends). Efficient formation of an oil/droplet emulsion during the droplet generation step of the droplet generation PCR protocol can be significantly affected by the presence of detergents in the sample; therefore, it is recommended to validate the compatibility of all cell-lysis buffers with droplet generation at the intended input concentration before proceeding with experimental samples and using the smallest practical volume of lysis buffer. The following lysis protocol results in minimal impact on droplet generation efficiency.
	<ol>
		<li>Prepare the lysis buffer containing 50 mM Tris-HCl pH 8.3, 1% TWEEN-20, and 200 µg/mL Proteinase K.</li>
		<li>Add 2.5 µL of the lysis buffer to each sample, seal the plate with an adhesive plate seal, and then centrifuge the samples at 1,000 x g for 1 min at 4 °C.</li>
		<li>Incubate the samples at 37 °C for 30 min on the thermocycler with the heated lid set at 105 °C to prevent the condensation of liquid on the plate cover.</li>
		<li>Add 7.5 µL of nuclease free water to each sample to get a final sample volume of 10 µL, re-seal the plate, and then centrifuge the samples at 1,000 x g for 1 min at 4 °C.</li>
		<li>Incubate on the thermocycler at 80 °C for 15 min (heated lid set to 105 °C) to inactivate Proteinase K.</li>
		<li>Centrifuge the samples at 1,000 x g for 1 min at 4 °C and then keep on ice.</li>
	</ol>
	</li>
	<li>Proceed to step 3 below. 
	NOTE: Cell lysates can be stored at -20 °C prior to continuing with step 3; however, use freshly eluted DNA where possible to eliminate the risk of freeze-thaw cycles resulting in DNA degradation.</li>
</ol>
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).
<ol>
	<li>Defrost all reagents on ice, vortex, and spin down briefly before use.</li>
	<li>Prepare the master mix (minus sample DNA) on ice as described in (Table 2). 
	NOTE: If only one target amplicon is being measured, add an additional 2.55 µL of nuclease-free water per sample to achieve the final volume of 22 µL. Prepare sufficient master mix for the number of samples and NTCs to be analyzed, plus enough to account for any 'blank' wells required on the droplet generation PCR plate (see next step).</li>
	<li>Vortex the prepared master mix briefly to mix and aliquot the required volume (22 µL minus input DNA volume) into each well of a droplet generation PCR 96-well plate. Arrange the samples in full columns on the 96-well plate; fill any empty wells in incomplete columns with master mix and use them as NTCs.</li>
	<li>Add the input DNA to each sample well and nuclease-free water/elution buffer to NTC wells to bring the total volume of each up to 22 µL.</li>
	<li>Seal the plate with an adhesive plate seal and place on a shaker at 2,000 rpm for 1 min, and then centrifuge at 1,000 x g for one min at 4 °C.</li>
	<li>Place the plate on ice and proceed to Step 4. 
	NOTE: Prepared plates can be stored on ice and protected from light for short periods of time (e.g., 1-2 h) before proceeding with droplet generation.</li>
</ol>
4. Generation of droplets
NOTE: Refer to Figure 1A for the schematic of the droplet generator instrument deck referenced in this section.
<ol>
	<li>Power on the droplet generator.</li>
	<li>Check that a bottle of droplet generating oil is loaded to position E of the instrument deck. If prompted, follow the on-screen instructions to replace the bottle.</li>
	<li>Click Configure Sample Plate on the touch screen. Select all columns on the sample plate that contain samples/blanks. 
	NOTE: Entry of plate name and details of the experiment at this step is optional and not required to continue setting up the experiment.</li>
	<li>Click OK to confirm the plate configuration; orange lights will illuminate on the instrument deck to indicate positions where consumables need to be loaded.</li>
	<li>Load the droplet generator cartridges to position B on the instrument deck so that all lights turn green.</li>
	<li>Place an empty tip waste trough into position C on the instrument deck.</li>
	<li>Remove lids from filter tip boxes and load to position D on the instrument stage so that all lights turn green.</li>
	<li>Remove the adhesive plate seal from the sample plate and load to position F of the instrument deck so that the light turns green.</li>
	<li>Place an empty droplet collection 96-well plate into a cool block (pre-chilled at -20 °C) and load to position G of the instrument deck so that the light turns green.</li>
	<li>Click Start Droplet Generation on the instrument touch screen.</li>
	<li>Click Start Run on the instrument touch screen. The instrument lid will close automatically. After initializing, the time remaining until droplet generation is complete will display on the instrument touch screen.</li>
	<li>Once the droplet generation is complete, check the sample collection plate visually to confirm that a droplet emulsion layer is present above the oil phase in each well (Figure 1B). 
	NOTE: If a droplet emulsion layer is not visible, do not proceed further with the experiment and check for any errors that may have occurred during previous steps of the protocol.</li>
	<li>Clear the used consumables from the instrument deck and empty the tip-waste trough.</li>
	<li>Power on the plate sealer, set the temperature to 180 °C and seal time to 5 s, and allow the machine to reach the operating temperature.</li>
	<li>Put the sample collection plate into the plate sealer plate-holder and place a fresh foil seal on top of the plate with the red line facing upwards. Take care not to touch the underside of the foil seal. 
	NOTE: The plate holder should be stored at room temperature (RT) and only placed into the plate sealer drawer when a foil seal is being applied to prevent heat transfer to the sample collection plate.</li>
	<li>When the plate sealer has reached operating temperature, press Eject on the instrument touchscreen; the drawer will open automatically.</li>
	<li>Place the plate holder into the plate sealer drawer and press Seal. The drawer will close automatically and then open again once sealing is complete.</li>
	<li>Replace the sealed plate on the cold block, turn off the plate sealer, and proceed immediately to Step 5. 
	NOTE: At this stage, the droplets are unstable, so ensure that the PCR step is commenced within 1 h of droplet generation finishing.</li>
</ol>
5. PCR
<ol>
	<li>Place the sealed droplet collection plate onto a PCR cycler equipped with a 96-deep-well block and run the thermal cycling protocol described in Table 3. 
	NOTE: Set heated lid temperature to 105 °C and sample volume to 40 µL. Denaturation and annealing/extension steps must include a temperature ramp-rate of 2 °C/s to ensure that the oil has time to equilibrate to the correct temperature.</li>
	<li>Proceed to Step 6. 
	NOTE: Once the thermal cycling protocol is complete, the droplets are more stable, and the plate can be stored for up to 4 days at 4 °C before proceeding with the next step.</li>
</ol>
6. Droplet reading
<ol>
	<li>Power on the droplet reader and connected computer.</li>
	<li>Check the fluid levels in the droplet reader oil and droplet reader waste bottles and fill/empty these, respectively, if required.</li>
	<li>Load the plate containing the post-PCR samples into the droplet reader plate holder, place the cover on top of the holder and secure it in place with the black locking clips. Do not remove the foil lid from the sample plate.</li>
	<li>Press the Open/Close button on the lid of the droplet reader to open it.</li>
	<li>Load the droplet reader plate holder containing the sample plate into the reader chamber and locate it securely onto the magnetic base.</li>
	<li>Press the Open/Close button on the lid of the droplet reader to close it.</li>
	<li>Open the analysis software interface, select the Add Plate tab and prepare a new sample plate as follows:
	<ol>
		<li>Click Add Plate and then Configure Plate. </li>
		<li>In the Plate Information tab enter a Plate Name, select Supermix for Probes (no dUTP) from the Supermix drop-down menu and enter a file name under Save Data File As.</li>
		<li>In the Well Selection tab highlight the wells to be analyzed and click Include Selected Wells. 
		NOTE: Steps 6.7.4 – 6.7.7 are optional.</li>
		<li>In the Well Information tab, select the wells to be annotated. 
		NOTE: Refer to Figure 1D for an example of the Well Information tab interface.</li>
		<li>Select Direct Quantification (DQ) from the Experiment Type drop-down menu, populate the Sample Description, Sample Type &amp; Target Name boxes, and add Well Notes and/or Plate Notes as required.</li>
		<li>Click Apply. The information in the fields will be applied to the selected wells.</li>
		<li>Repeat steps 6.7.4–6.7.6 until all sample wells have been annotated. </li>
	</ol>
	</li>
	<li>Once the plate configuration is complete, click Start Run.</li>
	<li>When the run is complete, discard the empty sample plate and empty the droplet reader waste bottle if necessary. </li>
	<li>Proceed to Step 7.</li>
</ol>
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 &gt;10,000 droplets and has two clearly separate populations of low fluorescence (negative) and high fluorescence (positive) droplets in each channel (Figure 2A).
<ol>
	<li>Select the Data Analysis tab and open the file to be analyzed from the My Datafiles menu.</li>
	<li>Click the 1D Amplitude tab for single-color experiments or the 2D Amplitude tab for two-color experiments.</li>
	<li>Apply a fluorescence threshold to each channel to differentiate the positive and negative droplets. The analysis software will attempt to apply an automatic threshold on a sample-by-sample and channel-by-channel basis, but if this appears to have failed, or looks inaccurate, then manually apply the threshold as follows:
	<ol>
		<li>Select the wells that require thresholding.</li>
		<li>Manually apply the threshold in the clear space between the positive and negative populations on the amplitude plot using the Threshold Line Mode tool. When setting manual thresholds in 2D view, ensure that the crosshair is placed so that the four droplet populations (double negative, single positive FAM, single positive HEX, and double positive) are clearly separated. 
		NOTE: Thresholding can either be done well-by-well or applied to multiple wells at once.</li>
	</ol>
	</li>
	<li>Once all samples have had a threshold applied, export the results as a .csv file for further analysis by clicking on the Data Table tab, clicking Import/Export and selecting Export Visible Data to CSV. Enter a suitable file name and click Save.</li>
	<li>Calculate the mtDNA copy number in the initial sample as follows: 
	Copy Number = mtDNA target concentration × 22 × 1/fraction of total lysate input 
	NOTE: If two mitochondrial probes are used, then the concentration of the two targets is averaged before making the calculation. For samples containing more than one cell, the copy number per cell is calculated by dividing by the number of cells in the sample. It is important to use the 'Concentration' value (i.e., number of copies per microliter) multiplied by the initial sample volume of 22 µL, rather than the 'Copies per 20 µL' value calculated in the .csv file, since the mtDNA copies left in the 2 µL overage required by the droplet generator must be taken into account when calculating the absolute copy number of the initial sample.</li>
	<li>For cells carrying heteroplasmic mtDNA deletions, calculate the copy number using the results from the mtDNA target outside the deleted region (i.e., the target present in both deleted and WT mtDNA molecules). The deletion heteroplasmy is calculated using the concentrations of the target inside the deleted region and the target outside the deleted region as follows: 
	<img alt="Percentage deleted mtDNA molecules formula; concentration ratio equation in genetic analysis." src="/files/ftp_upload/63870/63870eq01.jpg" style="margin: auto;" data-altupdatedat="1788282455000" data-alt="Equation 1"></li>
</ol>

<ol><li>Taanman, J. W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+mitochondrial+genome%3A+structure%2C+transcription%2C+translation+and+replication%5BTitle%5D)+AND+%22Biochimica+Biophysica+Acta%22%5BJournal%5D)">The mitochondrial genome: structure, transcription, translation and replication.</a> Biochimica Biophysica Acta. 1410 (2), 103-123 (1999).</li>
<li>Wai, T., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+role+of+mitochondrial+DNA+copy+number+in+mammalian+fertility%5BTitle%5D)+AND+%22Biology+of+Reproduction%22%5BJournal%5D)">The role of mitochondrial DNA copy number in mammalian fertility.</a> Biology of Reproduction. 83 (1), 52-62 (2010).</li>
<li>D'Erchia, A. M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Tissue-specific+mtDNA+abundance+from+exome+data+and+its+correlation+with+mitochondrial+transcription%2C+mass+and+respiratory+activity%5BTitle%5D)+AND+%22Mitochondrion%22%5BJournal%5D)">Tissue-specific mtDNA abundance from exome data and its correlation with mitochondrial transcription, mass and respiratory activity.</a> Mitochondrion. 20, 13-21 (2015).</li>
<li>Stewart, J. B., Chinnery, P. F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+dynamics+of+mitochondrial+DNA+heteroplasmy%3A+implications+for+human+health+and+disease%5BTitle%5D)+AND+%22Nature+Reviews%3A+Genetics%22%5BJournal%5D)">The dynamics of mitochondrial DNA heteroplasmy: implications for human health and disease.</a> Nature Reviews: Genetics. 16 (9), 530-542 (2015).</li>
<li>Durham, S. E., Samuels, D. C., Cree, L. M., Chinnery, P. F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Normal+levels+of+wild-type+mitochondrial+DNA+maintain+cytochrome+c+oxidase+activity+for+two+pathogenic+mitochondrial+DNA+mutations+but+not+for+m.3243A--%3EG%5BTitle%5D)+AND+%22American+Journal+of+Human+Genetics%22%5BJournal%5D)">Normal levels of wild-type mitochondrial DNA maintain cytochrome c oxidase activity for two pathogenic mitochondrial DNA mutations but not for m.3243A-->G.</a> American Journal of Human Genetics. 81 (1), 189-195 (2007).</li>
<li>Liu, H., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Wild-type+mitochondrial+DNA+copy+number+in+urinary+cells+as+a+useful+marker+for+diagnosing+severity+of+the+mitochondrial+diseases%5BTitle%5D)+AND+%22PloS+One%22%5BJournal%5D)">Wild-type mitochondrial DNA copy number in urinary cells as a useful marker for diagnosing severity of the mitochondrial diseases.</a> PloS One. 8 (6), 67146(2013).</li>
<li>Filograna, R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Modulation+of+mtDNA+copy+number+ameliorates+the+pathological+consequences+of+a+heteroplasmic+mtDNA+mutation+in+the+mouse%5BTitle%5D)+AND+%22Science+Advances%22%5BJournal%5D)">Modulation of mtDNA copy number ameliorates the pathological consequences of a heteroplasmic mtDNA mutation in the mouse.</a> Science Advances. 5 (4), (2019).</li>
<li>Wang, Y., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+increase+of+mitochondrial+DNA+content+in+endometrial+adenocarcinoma+cells%3A+a+quantitative+study+using+laser-captured+microdissected+tissues%5BTitle%5D)+AND+%22Gynecologic+Oncology%22%5BJournal%5D)">The increase of mitochondrial DNA content in endometrial adenocarcinoma cells: a quantitative study using laser-captured microdissected tissues.</a> Gynecologic Oncology. 98 (1), 104-110 (2005).</li>
<li>Boulet, L., Karpati, G., Shoubridge, E. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Distribution+and+threshold+expression+of+the+tRNA%28Lys%29+mutation+in+skeletal+muscle+of+patients+with+myoclonic+epilepsy+and+ragged-red+fibers+%28MERRF%29%5BTitle%5D)+AND+%22American+Journal+of+Human+Genetics%22%5BJournal%5D)">Distribution and threshold expression of the tRNA(Lys) mutation in skeletal muscle of patients with myoclonic epilepsy and ragged-red fibers (MERRF).</a> American Journal of Human Genetics. 51 (6), 1187-1200 (1992).</li>
<li>Lee, J., Hyeon, D. Y., Hwang, D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Single-cell+multiomics%3A+technologies+and+data+analysis+methods%5BTitle%5D)+AND+%22Experimental+and+Molecular+Medicine%22%5BJournal%5D)">Single-cell multiomics: technologies and data analysis methods.</a> Experimental and Molecular Medicine. 52 (9), 1428-1442 (2020).</li>
<li>Taylor, S. C., Laperriere, G., Germain, H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Droplet+Digital+PCR+versus+qPCR+for+gene+expression+analysis+with+low+abundant+targets%3A+from+variable+nonsense+to+publication+quality+data%5BTitle%5D)+AND+%22Scientific+Reports%22%5BJournal%5D)">Droplet Digital PCR versus qPCR for gene expression analysis with low abundant targets: from variable nonsense to publication quality data.</a> Scientific Reports. 7 (1), 2409(2017).</li>
<li>Hindson, C. M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Absolute+quantification+by+droplet+digital+PCR+versus+analog+real-time+PCR%5BTitle%5D)+AND+%22Nature+Methods%22%5BJournal%5D)">Absolute quantification by droplet digital PCR versus analog real-time PCR.</a> Nature Methods. 10 (10), 1003-1005 (2013).</li>
<li>Herbst, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Digital+PCR+quantitation+of+muscle+mitochondrial+DNA%3A+age%2C+fiber+type%2C+and+mutation-induced+changes%5BTitle%5D)+AND+%22Journals+of+Gerontology.+Series+A%3A+Biological+Sciences+and+Medical+Sciences%22%5BJournal%5D)">Digital PCR quantitation of muscle mitochondrial DNA: age, fiber type, and mutation-induced changes.</a> Journals of Gerontology. Series A: Biological Sciences and Medical Sciences. 72 (10), 1327-1333 (2017).</li>
<li>O'Hara, R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Quantitative+mitochondrial+DNA+copy+number+determination+using+droplet+digital+PCR+with+single-cell+resolution%5BTitle%5D)+AND+%22Genome+Research%22%5BJournal%5D)">Quantitative mitochondrial DNA copy number determination using droplet digital PCR with single-cell resolution.</a> Genome Research. 29 (11), 1878-1888 (2019).</li>
<li>Diaz, F., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Human+mitochondrial+DNA+with+large+deletions+repopulates+organelles+faster+than+full-length+genomes+under+relaxed+copy+number+control%5BTitle%5D)+AND+%22Nucleic+Acids+Research%22%5BJournal%5D)">Human mitochondrial DNA with large deletions repopulates organelles faster than full-length genomes under relaxed copy number control.</a> Nucleic Acids Research. 30 (21), 4626-4633 (2002).</li>
<li>Krishnan, K. J., Bender, A., Taylor, R. W., Turnbull, D. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+multiplex+real-time+PCR+method+to+detect+and+quantify+mitochondrial+DNA+deletions+in+individual+cells%5BTitle%5D)+AND+%22Analytical+Biochemistry%22%5BJournal%5D)">A multiplex real-time PCR method to detect and quantify mitochondrial DNA deletions in individual cells.</a> Analytical Biochemistry. 370 (1), 127-129 (2007).</li>
<li>Lowes, H., Pyle, A., Duddy, M., Hudson, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Cell-free+mitochondrial+DNA+in+progressive+multiple+sclerosis%5BTitle%5D)+AND+%22Mitochondrion%22%5BJournal%5D)">Cell-free mitochondrial DNA in progressive multiple sclerosis.</a> Mitochondrion. 46, 307-312 (2019).</li>
<li>Perier, C., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Accumulation+of+mitochondrial+DNA+deletions+within+dopaminergic+neurons+triggers+neuroprotective+mechanisms%5BTitle%5D)+AND+%22Brain%22%5BJournal%5D)">Accumulation of mitochondrial DNA deletions within dopaminergic neurons triggers neuroprotective mechanisms.</a> Brain. 136, Pt 8 2369-2378 (2013).</li>
<li>Cossarizza, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Guidelines+for+the+use+of+flow+cytometry+and+cell+sorting+in+immunological+studies+%28second+edition%29%5BTitle%5D)+AND+%22European+Journal+of+Immunology%22%5BJournal%5D)">Guidelines for the use of flow cytometry and cell sorting in immunological studies (second edition).</a> European Journal of Immunology. 49 (10), 1457(2019).</li>
<li>Espina, V., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Laser-capture+microdissection%5BTitle%5D)+AND+%22Nature+Protocols%22%5BJournal%5D)">Laser-capture microdissection.</a> Nature Protocols. 1 (2), 586-603 (2006).</li>
<li>Cree, L. M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+reduction+of+mitochondrial+DNA+molecules+during+embryogenesis+explains+the+rapid+segregation+of+genotypes%5BTitle%5D)+AND+%22Nature+Genetics%22%5BJournal%5D)">A reduction of mitochondrial DNA molecules during embryogenesis explains the rapid segregation of genotypes.</a> Nature Genetics. 40 (2), 249-254 (2008).</li>
<li>Belmonte, F. R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Digital+PCR+methods+improve+detection+sensitivity+and+measurement+precision+of+low+abundance+mtDNA+deletions%5BTitle%5D)+AND+%22Scientific+Reports%22%5BJournal%5D)">Digital PCR methods improve detection sensitivity and measurement precision of low abundance mtDNA deletions.</a> Scientific Reports. 6, 25186(2016).</li>
<li>Samuels, D. C., Schon, E. A., Chinnery, P. F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Two+direct+repeats+cause+most+human+mtDNA+deletions%5BTitle%5D)+AND+%22Trends+in+Genetics%22%5BJournal%5D)">Two direct repeats cause most human mtDNA deletions.</a> Trends in Genetics. 20 (9), 393-398 (2004).</li>
<li>Nissanka, N., Minczuk, M., Moraes, C. T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Mechanisms+of+mitochondrial+DNA+deletion+formation%5BTitle%5D)+AND+%22Trends+in+Genetics%22%5BJournal%5D)">Mechanisms of mitochondrial DNA deletion formation.</a> Trends in Genetics. 35 (3), 235-244 (2019).</li>
<li>Macaulay, I. C., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Separation+and+parallel+sequencing+of+the+genomes+and+transcriptomes+of+single+cells+using+G%26T-seq%5BTitle%5D)+AND+%22Nature+Protocols%22%5BJournal%5D)">Separation and parallel sequencing of the genomes and transcriptomes of single cells using G&T-seq.</a> Nature Protocols. 11 (11), 2081-2103 (2016).</li>
<li>Ludwig, L. S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Lineage+tracing+in+humans+enabled+by+mitochondrial+mutations+and+single-cell+genomics%5BTitle%5D)+AND+%22Cell%22%5BJournal%5D)">Lineage tracing in humans enabled by mitochondrial mutations and single-cell genomics.</a> Cell. 176 (6), 1325-1339 (2019).</li>
<li>Rooney, J. P., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((PCR+based+determination+of+mitochondrial+DNA+copy+number+in+multiple+species%5BTitle%5D)+AND+%22Methods+in+Molecular+Biology%22%5BJournal%5D)">PCR based determination of mitochondrial DNA copy number in multiple species.</a> Methods in Molecular Biology. 1241, 23-38 (2015).</li>
<li>Kamitaki, N., Usher, C. L., McCarroll, S. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Using+droplet+digital+PCR+to+analyze+allele-specific+RNA+expression%5BTitle%5D)+AND+%22Methods+in+Molecular+Biology%22%5BJournal%5D)">Using droplet digital PCR to analyze allele-specific RNA expression.</a> Methods in Molecular Biology. 1768, 401-422 (2018).</li>
<li>Maeda, R., Kami, D., Maeda, H., Shikuma, A., Gojo, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((High+throughput+single+cell+analysis+of+mitochondrial+heteroplasmy+in+mitochondrial+diseases%5BTitle%5D)+AND+%22Scientific+Reports%22%5BJournal%5D)">High throughput single cell analysis of mitochondrial heteroplasmy in mitochondrial diseases.</a> Scientific Reports. 10 (1), 10821(2020).</li>
<li>Quan, P. L., Sauzade, M., Brouzes, E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((dPCR%3A+A+Technology+Review%5BTitle%5D)+AND+%22Sensors+%28Basel%29%22%5BJournal%5D)">dPCR: A Technology Review.</a> Sensors (Basel). 18 (4), (2018).</li>
<li>Lin, X., Huang, X., Urmann, K., Xie, X., Hoffmann, M. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Digital+loop-mediated+isothermal+amplification+on+a+commercial+membrane%5BTitle%5D)+AND+%22ACS+Sensors%22%5BJournal%5D)">Digital loop-mediated isothermal amplification on a commercial membrane.</a> ACS Sensors. 4 (1), 242-249 (2019).</li>
<li>Li, Z., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Fully+integrated+microfluidic+devices+for+qualitative%2C+quantitative+and+digital+nucleic+acids+testing+at+point+of+care%5BTitle%5D)+AND+%22Biosensors+and+Bioelectronics%22%5BJournal%5D)">Fully integrated microfluidic devices for qualitative, quantitative and digital nucleic acids testing at point of care.</a> Biosensors and Bioelectronics. 177, 112952(2021).</li>
</ol>

No conflicts of interest to disclose.

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 ensu...

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.

<table><tbody><tr><td>50% Tween-20 solution</td><td>Novex</td><td>3005</td><td></td></tr><tr><td>Automated droplet-generating oil&amp;nbsp;</td><td>Bio Rad</td><td>1864110</td><td>Commercial oil formulation used to generate the oil/droplet emulsion (used in Protocol Step 4.1)</td></tr><tr><td>C1000 PCR machine with deep-well block</td><td>Bio Rad</td><td>1851197</td><td>PCR thermocycler equipped with a deep-well heating block, used for cell lysis (Protocol Step 2.1.2.) and PCR cycling (Protocol Step 5)</td></tr><tr><td>Collection plate cooling block</td><td>Bio Rad</td><td>12002819</td><td>Cooling block that keeps samples chilled during droplet generation (used in Protocol step 4.3)</td></tr><tr><td>ddPCR 96-well plates&amp;nbsp;</td><td>Bio Rad</td><td>12001925</td><td>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)</td></tr><tr><td>ddPCR droplet reader oil&amp;nbsp;</td><td>Bio Rad</td><td>1863004</td><td>Commercial oil formulation used by the droplet reader (used in Protocol step 6.1)</td></tr><tr><td>ddPCR Supermix for Probes (no dUTP)</td><td>Bio Rad</td><td>1863023</td><td>Commercial supermix for use in ddPCR experiments utilising probes (used in Protocol Step 3.3)</td></tr><tr><td>DG32 automated droplet generator cartridges&amp;nbsp;</td><td>Bio Rad</td><td>1864108</td><td>Microfluidic cartridges used in the QX200 AutoDG droplet generator to generate the oil/droplet emulsion (used in Protocol Step 4.3)</td></tr><tr><td>Fetal bovine serum</td><td>Gibco</td><td>10270-106</td><td>Qualified fetal bovine serum</td></tr><tr><td>Foil plate covers&amp;nbsp;</td><td>Bio Rad</td><td>1814040</td><td>Foil plate covers used to seal droplet collection plates after droplet generation (used in Protocol step 4.6)</td></tr><tr><td>HEK 293T cells</td><td>Takara</td><td>632180</td><td>Commercial subclone of the transformed human embryonic kidney cell line, HEK 293, expressing the SV40 Large-T antigen</td></tr><tr><td>HeLa cells</td><td>ECACC</td><td>93021013</td><td>Human cervix epitheloid carcinoma cells</td></tr><tr><td>High glucose DMEM</td><td>Gibco</td><td>13345364</td><td>4.5g/L D-Glucose, with L-glutamine and sodium pyruvate</td></tr><tr><td>Human cybrids</td><td>University of Miami</td><td></td><td></td></tr><tr><td>Mouse embryonic fibroblasts&amp;nbsp;</td><td>Newcastle University</td><td>Immortalized from C57Bl/6 mice</td><td></td></tr><tr><td>Nuclease-free water</td><td>Ambion</td><td>AM9937</td><td></td></tr><tr><td>PCR plate seals</td><td>Pierce</td><td>SP-0027</td><td>Clear adhesive plate seals, only used pre-droplet generation (foil seal must be used in step 4.6)</td></tr><tr><td>Pipet Tip Waste Bins</td><td>Bio Rad</td><td>1864125</td><td>Disposable collection bin used to collect discarded tips in the QX200 AutoDG droplet generator (used in Protocol step 4.3)</td></tr><tr><td>Pipet tips for AutoDG system&amp;nbsp;</td><td>Bio Rad</td><td>1864120</td><td>Filtered pipet tips designed for use in the QX200 AutoDG droplet generator (used in Protocol step 4.3)</td></tr><tr><td>Primary human dermal fibroblast cells</td><td>Newcastle Biobank</td><td></td><td></td></tr><tr><td>Primers/Probes</td><td>IDT</td><td>N/A</td><td>Exact primer/probe sequences will be assay dependent. Primers and probes used in this study are given in Table 1</td></tr><tr><td>Proteinase K 20 mg/mL solution</td><td>Ambion</td><td>AM2546</td><td></td></tr><tr><td>PX1 PCR plate sealer</td><td>Bio Rad</td><td>1814000</td><td>Applies foil seals to ddPCR sample plates after droplet generation (used in Protocol Step 4.6)</td></tr><tr><td>QX Manager software</td><td>Bio Rad</td><td>12012172</td><td>Droplet reader set up &amp;amp; analysis software (used in Protocol Steps 6 &amp;amp; 7)</td></tr><tr><td>QX200 AutoDG droplet generator</td><td>Bio Rad</td><td>1864101</td><td>Automated microfluidic droplet generator (used in Protocol Step 4)</td></tr><tr><td>QX200 droplet reader</td><td>Bio Rad</td><td>1864003</td><td>Droplet reader (used in Protocol Step 6)</td></tr><tr><td>Trizma pre-set crystals pH 8.3</td><td>Sigma</td><td>T8943-100G</td><td></td></tr></tbody></table>

measuring single-cell mitochondrial dna copy number and heteroplasmy using digital droplet polymerase chain reaction

Measuring Single-Cell Mitochondrial DNA Copy Number and Heteroplasmy Using Digital Droplet Polymerase Chain Reaction

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.

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 (<strong class="xfig"...

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 ...

measuring-single-cell-mitochondrial-dna-copy-number-heteroplasmy

Research

JoVE Journal

Medicine

In this video, we demonstrate the duplex digital PCR (ddPCR) technique — a modification of the traditional PCR technique that is useful in detecting two different genetic markers simultaneously. A single PCR reaction is partitioned into nanoliter-sized emulsified droplets that are independently amplified, and the detection of differently colored fluorescence amplification signals from the fraction of droplets is used to compute the initial concentration of the target sequences.

Duplex Digital PCR for Simultaneous Quantification of Dual Genetic Markers

This video demonstrates the multiplex droplet polymerase chain reaction technique for the simultaneous quantification of methylated DNA regions in different white blood cell populations as an early biomarker for clinical diagnosis and prognosis.

Multiplex Droplet Polymerase Chain Reaction to Detect DNA Methylation in Leukocytes

A sensitive and accurate method for cell-free microRNAs quantification using a dye-based chemistry and droplet digital PCR technology is described.

Circulating MicroRNA Quantification Using DNA-binding Dye Chemistry and Droplet Digital PCR

Here we describe a method amenable to simultaneously quantitate and genome-wide map ribonucleotides in highly intact DNA at single-nucleotide resolution, combining enzymatic cleavage of genomic DNA with its alkaline hydrolysis and subsequent 5´-end sequencing.

Simultaneous Mapping and Quantitation of Ribonucleotides in Human Mitochondrial DNA

Here, we present a protocol to allow accurate quantification of mitochondrial DNA (mtDNA) methylation. In this protocol, we describe an enzymatic digestion of DNA with BamHI coupled with a bioinformatic analysis pipeline which can be used to avoid overestimation of mtDNA methylation levels caused by the secondary structure of mtDNA.

Methodology for Accurate Detection of Mitochondrial DNA Methylation

Here, we present a protocol to accurately quantify multiple genetic alterations of a target region in a single reaction using drop-off ddPCR and a unique pair of hydrolysis probes.

Single Droplet Digital Polymerase Chain Reaction for Comprehensive and Simultaneous Detection of Mutations in Hotspot Regions

Epigenetic markers are used for white blood cell (WBC) subtyping through the quantification of DNA methylation patterns. This protocol presents a multiplex droplet polymerase chain reaction (mdPCR) method using a thermoplastic elastomer (TPE)-based microfluidic device for droplet generation allowing for precise and multiplex methylation-specific target quantification of WBC differential counts.

Methylation Specific Multiplex Droplet PCR using Polymer Droplet Generator Device for Hematological Diagnostics

Pyrosequencing assays enable the robust and rapid genotyping of mitochondrial DNA single nucleotide polymorphisms in heteroplasmic cells or tissues.

Genotyping Single Nucleotide Polymorphisms in the Mitochondrial Genome by Pyrosequencing

This protocol describes the detailed method of digital droplet PCR (dd-PCR) for precise quantification of circular RNA (circRNA) levels in cells using divergent primers.

Quantification of Circular RNAs Using Digital Droplet PCR

A procedure for studying the dynamics of mitochondrial DNA (mtDNA) metabolism in cells using a multi-well plate format and automated immunofluorescence imaging to detect and quantify mtDNA synthesis and distribution is described. This can be further used to investigate the effects of various inhibitors, cellular stresses, and gene silencing on mtDNA metabolism.

High-Throughput Image-Based Quantification of Mitochondrial DNA Synthesis and Distribution

This study presents a rapid ddPCR method for quantifying soil microbial biomass and fungal-to-bacterial ratios. Although ddPCR-derived microbial biomass carbon shows only weak correlation with chloroform fumigation estimates, the method provides a sensitive, scalable molecular approach for assessing microbial abundance without additional soil sampling or the use of harsh chemicals.

Determination of Microbial Biomass in Soil using Digital Droplet PCR (ddPCR)

Learn to measure mitochondrial DNA copy number and deletion heteroplasmy in single cells using digital droplet PCR. This protocol enables precise analysis for mitochondrial research.

Measuring Single-Cell Mitochondrial DNA Copy Number and Heteroplasmy Using Digital Droplet Polymerase Chain Reaction

Explore More Videos