Quantification of human mitochondrial DNA using synthesized DNA standards
09765400 · 2017-09-19
Assignee
Inventors
Cpc classification
C12Q1/6883
CHEMISTRY; METALLURGY
International classification
Abstract
A real-time quantitative PCR assay that utilizes a duplex, synthetic DNA standard to ensure optimal quality assurance and quality control. One embodiment of the invention facilitates amplification of mtDNA by focusing on a 105-base pair target sequence that is minimally homologous to non-human mtDNA. The present invention can also be used to identify the presence of PCR inhibitors and thus indicate the need for sample repurification.
Claims
1. An oligonucleotide standard adapted for an assay for quantification of human mitochondrial DNA, said standard comprising: SEQ ID NO: 1 (forward) and SEQ ID NO: 2 (reverse).
Description
BRIEF DESCRIPTION OF THE DRAWINGS
(1) A more complete appreciation of the present invention, and many of the features and advantages of the present invention, will be readily apparent as the same becomes better understood by reference to the following detailed description when considered in conjunction with the accompanying drawings, wherein:
(2)
(3)
(4)
(5)
(6)
(7)
(8)
(9)
(10)
DETAILED DESCRIPTION OF THE INVENTION
(11) The invention is directed to real-time quantitative PCR (qPCR) assays to quantify human mtDNA and ensure sufficient quantity and quality exists for subsequent forensic analysis. However, the invention may also be used in the medical, anthropological, and other non-forensic fields. In one embodiment of the invention, the assay is used to quantify human mtDNA by detecting the presence of the NADH dehydrogenase subunit 5 gene in sample DNA preparations.
(12) In a preferred embodiment of the invention, the real-time qPCR assay uses a probe, primers, and a synthetic DNA standard. In a more preferred embodiment, the standard is a unique, linear, duplex, synthetic DNA standard. Benefits of using a synthetic standard include ease of production, high quality control, high purity, low cost (compared to the construction and production of traditional plasmid-based DNA standards), high yield, and rapid production. The inventive standard provides optimal quality assurance and quality control because of its purity and the absence of contaminating DNA, RNA, and protein, which could affect its absorbance of ultraviolet light (260 nm) and result in inaccurate quantification.
(13) In contrast to a linear synthetic standard, experiments involving a plasmid-based standard indicate that DNA supercoiling may result in inefficient amplification and significantly reduced Ct values for the plasmid standard, resulting in a ˜50-fold overestimation of copy numbers for samples, which are not supercoiled. This effect could be caused by ineffective annealing of primer and probe due to topological constraints. While this effect may be overcome by nicking or linearization of the plasmid standard followed by repurification, these steps are not required for linear synthetic standards.
(14) Target Sequence
(15) By nature, forensic DNA evidence may contain a mixture of human DNA and contaminating DNA from non-human sources. Furthermore, it can be difficult to determine whether small pieces of recovered skeletal remains are from a human or non-human source. qPCR primers and probes can potentially bind to and amplify non-human mtDNA, resulting in inaccurate yield determinations that may impact downstream analysis. Therefore, one criterion for selecting the qPCR target sequence was low homology with non-human mammalian DNA. Several potential qPCR target regions in the human mtDNA genome were compared to the mtDNA genome sequences of sixteen domesticated and wild mammalian species, which could pose as environmental contaminants, and were evaluated for sequence homology. A sequence within the NADH dehydrogenase subunit 5 gene exhibited relatively low homology with the non-human DNA sequences, and was selected as the qPCR target region for one embodiment of the present invention. Thus, in a preferred embodiment of the invention, the qPCR target region corresponds to positions 13,288-13,392 of the mtDNA revised Cambridge Reference Sequence (NCBI Accession No. NC_012920.1). The small amplicon size (105 bp) facilitates amplification even for degraded DNA.
(16) qPCR Synthetic Standard
(17) The target region contains known single nucleotide polymorphisms (SNP), any of which may be used in forensic analyses pursuant to the present invention. To address the possibility that the qPCR standard could be a potential contaminant in such cases, the qPCR standard, according to one embodiment of the present invention, contains a “signature” or marker sequence at positions 13,345-13,349 of the mtDNA revised Cambridge Reference Sequence (NCBI Accession No. NC_012920.1), representing a change from “GCCAT” to “TAACG” in the forward strand and a change from “atggc” to “cgtta” in the reverse strand. This signature sequence in the forward strand encodes a stop codon (Och, ochre) and a Met (methionine).fwdarw.Arg (arginine) substitution when translated (
(18) In one embodiment, the standard also contains five mtDNA base pairs beyond the primer binding sites, which eliminates qPCR quantification inaccuracies that may arise due to undesirable n−1, n−2, etc. oligo products that may be created during oligo synthesis and persist through oligo purification. The standard, containing the mtDNA region amplified (or target sequence) and five additional base pairs on each end, corresponds to positions 13,283-13,397 of the revised Cambridge Reference Sequence (NCBI Accession No. NC_012920.1), except that positions 13,345-13,349 are changed from “GCCAT” to “TAACG.”
(19) qPCR Assay
(20) Whereas the exemplary qPCR assay described herein employs a specific protocol and instrument, those of ordinary skill in the art will recognize that it may be easily adapted for any real-time qPCR protocol, instrument, and/or system that can detect 6FAM, VIC, TAMRA, ROX, and other compatible dyes. In addition, the exemplary assay described herein utilizes a “fast” qPCR amplification protocol, which results in a rapid run time of approximately 40 minutes; however, appropriate modifications known to those of ordinary skill in the art may be made to the protocol to carry out “standard” qPCR, e.g., modifications to, for example, the plate, sample block, master mix, and/or cycling parameters.
(21) The qPCR assay of the present invention is suitable for accurate and precise mtDNA quantification in DNA samples prior to forensic analysis and is reliable, robust, and highly sensitive. In addition, the invention may be utilized in biomedical and other applications that demand similar requirements.
EXAMPLES
(22) The examples below are intended to be purely exemplary and are not intended to limit the disclosure or any aspect of the invention.
(23) DNA Samples
(24) The primers and probe of the present invention were designed to specifically amplify human mtDNA. They hybridize to sequences that are minimally variable, maximizing amplification of a variety of polymorphic sequence types. To attain the highest quantification accuracy of human mtDNA in a sample, the non-specific amplification of contaminating non-human DNA should be minimal. However, the possibility does exist that the qPCR primers and probe may non-specifically bind to and amplify non-human DNA, resulting in non-specific amplification. Furthermore, in forensic DNA analysis, human DNA can be contaminated with DNA from domestic pets, farm animals, wildlife, bacteria, or fungi.
(25) To demonstrate specificity of the qPCR assay for human DNA, a standard quantity (100 ng) of DNA from each of twenty mammalian species and from ten other species, including non-mammalian vertebrates, bacteria, and fungi, was tested. The “apparent” copy numbers obtained via the mtDNA qPCR assay for each of the non-human species were compared to the copy number obtained for the same quantity (100 ng) of commercially prepared human DNA and were cited as a percentage of the human DNA copy number. For example, if the qPCR assay were to fully non-specifically amplify a given non-human sample, then the “apparent” copy number obtained for that sample would be equivalent to, or 100% of, the copy number obtained for the known human sample. This approach assumed that all non-human DNA was free of human DNA and vice versa and that the ratio of mtDNA to total DNA was equivalent for all species.
(26) The human (Zyagen GH-180) and other vertebrate DNA samples were purchased from Zyagen Laboratories (10225 Barnes Canyon Road, San Diego, Calif., 92121), and the bacterial and fungal DNA were purchased from the American Type Culture Collection (ATCC) (10801 University Boulevard, Manassas, Va., 20110). The non-human vertebrate DNA included DNA from bovine (Zyagen GB-110), cat (Zyagen GC-130), chicken (Zyagen GC-120), deer (Zyagen GD-140), dog (Zyagen GD-150), donkey (Zyagen GD-160), equine (Zyagen GE-170), ferret (Zyagen GF-180), fish (Zyagen GF-190), gerbil (Zyagen GG-120), goat (Zyagen GG-150), guinea pig (Zyagen GG-160), hamster (Zyagen GH-170), monkey (baboon) (Zyagen GM-140B), monkey (cynomolgus) (Zyagen GM 140C), monkey (rhesus) (Zyagen GM-140R), mouse (Zyagen GM-150), pig (Zyagen GP-160), rabbit (Zyagen GR-170), rat (Zyagen GR-180), sheep (Zyagen GS-190), and turkey (Zyagen GT-150). The bacterial DNA included DNA from Bacillus subtilis (ATCC 23857D-5), Clostridium perfringens (ATCC 13124D), Escherichia coli (ATCC 10798D), Pseudomonas aeruginosa (ATCC 17933D), Staphylococcus epidermidis (ATCC 12228D), and Streptococcus pneumoniae (ATCC BAA-334D). The fungal DNA included DNA from Aspergillus oryzae (ATCC 42149D-2) and Candida albicans (ATCC 14053D).
(27) The results showed that only slight non-specific amplification was apparent for 19 of the 30 non-human species examined with the greatest for cynomolgus monkey DNA, which was quantified at a level representing less than 1% of the human DNA copy number. At this level, an improbable amount of contaminating cynomolgus monkey DNA, 100-fold more, would be required to attain a signal equivalent to human DNA. Two additional primate species DNA, rhesus monkey and baboon, exhibited cross-reactivity approximately 30- and 90-fold lower than the cynomolgus monkey, respectively. Farm, wildlife, and pet species, often encountered with recovered human remains, were quantified at levels less than or equal to 0.4% that of the human DNA copy number, which demonstrated no appreciable non-specific amplification. Virtually no cross-reactivity was observed among the bacterial and fungal species examined. Given the high specificity of the qPCR assay according to the present invention, human mtDNA quantification should be virtually unaffected by the presence of mtDNA copies of the non-human DNA types listed here. See Table 4.
(28) TABLE-US-00001 TABLE 4 qPCR Cross-reactivity with Non-human DNA Apparent Copy Number of Human DNA Source mtDNA per 100 ng of DNA* % of Human Human 1,900,000 100% Cynomolgus Monkey 18,000 0.9% Cow 7,800 0.4% Guinea Pig 1,200 0.1% Rhesus Monkey 490 0.03% Baboon 140 0.01% *The following organisms were quantified at 35 copies or less per 100 ng DNA (<0.002% of human): chicken, deer, dog, horse, ferret, fish, gerbil, hamster, mouse, rabbit, rat, sheep, turkey, and C. albicans. mtDNA was not detected in DNA from cat, donkey, goat, pig, A. oryzae, B. subtilis, C. perfringens, E. coli, P. aeruginosa, S. epidermis, and S. pneumonia.
(29) To evaluate the present invention with a “forensic-type” specimen, buccal swab DNA was extracted using an organic extraction method. Briefly, a half buccal swab was incubated at 56° C. for 2 hours in a buffer containing 10 mM Tris, 100 mM NaCl, 39 mM DTT, 10 mM EDTA, 2% SDS, and 1.2 units of proteinase K. The swab was removed and the solution was extracted with an equal volume of phenol/chloroform/isoamyl alcohol (25:24:1, v/v). The aqueous layer was purified using a Microcon 100 centrifugal device (Millipore, Billerica, Mass.). DNA was recovered from the device with 60 μl of water. The samples were stored at −20° C. until experimentation and at 4° C. during the period of experimentation.
(30) qPCR Standard
(31) The dsT8sig qPCR standard consisted of two complementary, PAGE-purified, synthetic oligonucleotides (Ultramers™, Integrated DNA Technologies, Inc., Coralville, Iowa) (Table 1) that correspond to positions 13,283-13,397 of the mtDNA revised Cambridge Reference Sequence (NCBI Accession No. NC_012920.1), i.e., the target sequence (with the signature sequence change) plus five additional base pairs at both the 5′ and 3′ ends. For quality control, as noted, the forward strand of the standard contained a signature or marker sequence, TAACG, corresponding to positions 13,345-13,349 of the reference sequence and the reverse strand contained the reverse complement, CGTTA, at this position.
(32) TABLE-US-00002 TABLE 1 Sequences of DNA Standard, Primers, and Probe qPCR Standard (forward strand) Tfor8sig: 5′-CAA TCG GCA TCA ACC AAC CAC ACC TAG CAT TCC TGC ACA TCT GTA CCC ACG CCT TCT TCA AAT AAC GAC TAT TTA TGT GCT CCG GGT CCA TCA TCC ACA ACC TTA ACA ATG AAC A-3′ [SEQ ID NO: 1] qPCR Standard (reverse strand) Trev8sig: 5′-TGT TCA TTG TTA AGG TTG TGG ATG ATG GAC CCG GAG CAC ATA AAT AGT CGT TAT TTG AAG AAG GCG TGG GTA CAG ATG TGC AGG AAT GCT AGG TGT GGT TGG TTG ATG CCG ATT G-3′ [SEQ ID NO: 2] Forward Primer Qfor8: 5′-GGC ATC AAC CAA CCA CAC CTA-3′ [SEQ ID NO: 3] Reverse Primer Qrev8: 5′-ATT GTT AAG GTT GTG GAT GAT GGA-3′ [SEQ ID NO: 4] Probe (MGBNFQ = minor groove binder nonfluorescent quencher) QRL8: 5′-6FAM CAT TCC TGC ACA TCT G MGBNFQ-3′ [SEQ ID NO: 5]
(33) The forward and reverse qPCR standard oligonucleotides were separately reconstituted in Tris EDTA (TE) buffer (10 mM Tris-HCl, pH 8.0, 0.1 mM EDTA) and quantified by measuring their absorbance of ultraviolet light at a wavelength of 260 nm and applying the extinction coefficients 1,082,000 and 1,138,100 L/(mole.Math.cm), respectively. The concentration of each oligonucleotide was then adjusted to 2 μM and the adjusted oligonucleotides were mixed in equal proportions to generate a 1 μM double-stranded, primary qPCR standard stock. The primary stock was diluted with TE buffer to generate a secondary standard stock containing 10.sup.10 copies of standard per 2 μl. Primary and secondary standards were aliquoted (10 μl) and stored at −80° C. The assay dilution series was prepared fresh daily. The secondary stock was diluted serially with TE buffer to generate the following assay dilution series: 10.sup.8, 10.sup.7, 10.sup.6, 10.sup.5, 10.sup.4,10.sup.3, 10.sup.2, and 10.sup.1 copies per 2 μl (
(34) Real-time qPCR
(35) qPCR experiments were conducted with 2 μl of sample DNA, standard DNA, or control reagents (described below) per well of a 96-well plate in a 25 μl duplex reaction containing 12.5 μl of TaqMan 2× Fast Universal Master Mix (no UNG) (Applied Biosystems, Foster City, Calif.), 900 nM HPLC-purified forward and reverse primers (SEQ ID NOS: 3 and 4, respectively) (Integrated DNA Technologies, Inc.), and 250 nM HPLC-purified 5′ 6FAM-labeled probe containing a 3′ minor groove binder nonfluorescent quencher (SEQ ID NO: 5) (Applied Biosystems). The qPCR reaction also included the TaqMan Exogenous Internal Positive Control Reagents (Applied Biosystems) according to the manufacturer's recommendations, except for the blocking reagent, which, in the control wells, contained the Exogenous IPC Blocking Reagent at a final concentration of 0.8× instead of the recommended 1×. Additional control reactions consisted of water NTCs (no template controls) and TE NTCs. See
(36) Reactions were amplified in duplicate, unless otherwise noted, on a 7900HT FAST Sequence Detection System (Applied Biosystems) in “fast” mode: 20 seconds at 95° C. followed by 40 cycles of 1 second at 95° C. and 20 seconds at 60° C. Data were analyzed for both mtDNA- and IPC-specific probes with Sequence Detection Software (SDS) version 2.2.2 (Applied Biosystems). Cycle threshold (Ct) values were determined at 0.2 ΔRn using the automatic baseline algorithm. qPCR amplification efficiencies were calculated using the slope of the standard plot regression line: efficiency=[10.sup.(−1/slope)]−1.
(37) Exogenous IPC Analysis
(38) In another embodiment of the present invention, amplification of an IPC is used to determine whether PCR inhibitors are present in a DNA sample. The amplification involves inclusion of IPC-specific primers and probe in the qPCR reaction mixture as well as IPC DNA, which is expected to amplify in the absence of PCR inhibitors. The highest level of IPC amplification, and therefore the lowest IPC Ct values, is expected for the NTC controls that lack PCR inhibitors. In contrast, in the presence of inhibitors, IPC Ct values should increase relative to that for the NTCs or be “undetermined,” for low levels and high levels of inhibitor, respectively.
(39) Specifically, a controlled inhibition experiment was conducted using two-fold serial dilutions of EDTA varying in final concentration from 10 to 0.078 mM per well. The qualitative analysis with EDTA as inhibitor (
(40) IPC Ct values can also increase in the absence of inhibitors when high levels of sample mtDNA compete for polymerase and dNTPs during amplification, termed “inhibition by competition.” This phenomenon is readily observed with the inhibitor-free standards where the IPC Ct is slightly increased with 10.sup.4 mtDNA standard copies per well (
(41) mtDNA Hypervariable Region (HVR) PCR
(42) An important benefit of the assay of the present invention is the conservation of forensic DNA samples to allow for additional analysis. In this regard, the assay described herein can be used to determine the optimum number of mtDNA copies needed for HVR amplification.
(43) To determine the optimal quantity of human mtDNA needed for down-stream forensic analysis, ten-fold serial dilutions of human DNA (Zyagen Laboratories, GH-180) were subjected to HVR amplification in duplicate as follows. Hypervariable region 1 (HVR1) amplification: 10 μl of DNA in a reaction mixture containing 5 units AmpliTaq Gold DNA polymerase (Applied Biosystems), 1× GeneAmp PCR buffer (Applied Biosystems), 160 ng/μl BSA (#A3550, Sigma-Aldrich, St. Louis, Mo.), dNTP mixture (dATP, dCTP, dGTP, and dTTP, 250 μM each), 600 nM forward primer CAC CAT TAG CAC CCA AAG CT [SEQ ID NO: 6], and 600 nM reverse primer GAG GAT GGT GGT CAA GGG AC [SEQ ID NO: 7]. Hypervariable region 2 (HVR2) amplification reaction mixtures were similar, except for the forward primer, CTC ACG GGA GCT CTC CAT GC [SEQ ID NO: 8], and the reverse primer, CTG TTA AAA GTG CAT ACC GCC [SEQ ID NO: 9]. Reaction mixtures were amplified using a 95° C. hold for 9 minutes; 32 cycles consisting of 95° C. for 10 seconds, 60° C. for 30 seconds, and 72° C. for 30 seconds; and a 4° C. hold. Amplicons were quantified with the Agilent DNA 1000 kit (Agilent Technologies, Waldbronn, Germany).
(44) The results revealed that 10 pg of human genomic DNA (440 mtDNA copies) subjected to HVR1 and HVR2 amplification yielded 3 ng/μl and 2 ng/μl of amplicon, respectively, concentrations that are sufficient for sequence analysis (
(45) Assay Sensitivity and Reproducibility
(46) Data generated for the highly purified synthetic standard demonstrated consistently high assay sensitivity and reproducibility. Among fifteen individual assays, 100,000,000 (10.sup.8) copies of the highly purified standard per well exhibited an average Ct of 13.2 (range 12.5-13.7). The subsequent 10-fold dilution (10.sup.7/well) exhibited a Ct value of 16.6 (range 15.9-17.0). The difference in Ct values between each successive dilution of the standard averaged about 3.4 cycles. Ten copies of standard per well were always detected at an average Ct of 37.0 (range 35.9 - 38.2), well within the 40 cycles prescribed for the assay.
(47) High consistency of results and assay reproducibility were also demonstrated by compiling the standard curve results from the 15 separate assays (
(48) TABLE-US-00003 TABLE 2 Standard Ct Data for Fifteen Separate qPCR Assays Copies/well Average Minimum Maximum Range 10.sup.8 13.2 12.5 13.7 1.2 10.sup.7 16.6 15.9 17.0 1.1 10.sup.6 20.0 19.3 20.4 1.1 10.sup.5 23.5 22.9 23.9 0.9 10.sup.4 26.9 26.1 27.3 1.2 10.sup.3 30.4 29.5 30.9 1.5 10.sup.2 33.6 32.4 34.1 1.7 10.sup.1 37.0 35.9 38.2 2.3 10.sup.0* 40.5 39.5 41.0 1.5 *Represents Y-intercept data
(49) TABLE-US-00004 TABLE 3 Linear Regression Data for Fifteen Separate qPCR Assays Run Slope Y-Intercept R.sup.2 Efficiency Calc. Qty.* 1 −3.35 39.49 1.000 98.7% 6.49E+05 2 −3.41 39.98 1.000 96.4% 7.25E+05 3 −3.43 40.58 0.999 95.9% 1.02E+06 4 −3.40 40.33 0.999 96.9% 9.65E+05 5 −3.42 40.48 1.000 96.0% 9.70E+05 6 −3.32 39.77 0.999 100.2% 9.13E+05 7 −3.43 40.78 0.999 95.5% 1.13E+06 8 −3.39 40.54 0.997 97.2% 1.15E+06 9 −3.36 40.61 1.000 98.3% 1.34E+06 10 −3.44 40.79 0.999 95.1% 1.08E+06 11 −3.45 40.73 0.998 94.7% 9.99E+05 12 −3.41 40.68 1.000 96.3% 1.14E+06 13 −3.39 40.58 0.998 97.3% 1.19E+06 14 −3.46 40.92 1.000 94.7% 1.13E+06 15 −3.46 40.99 1.000 94.5% 1.15E+06 Average −3.41 40.48 0.999 96.5% 1.0E+06 Standard 0.04 0.42 0.001 1.6% 1.8E+05 Deviation Coefficient 0.012 0.010 0.001 0.017 0.171 of Variation Minimum −3.46 39.49 0.997 94.5% 6.49E+05 Maximum −3.32 40.99 1.000 100.2% 1.34E+06 Range 0.14 1.50 0.003 5.7% 6.87E+05 *mtDNA quantity calculated for a theoretical sample exhibiting a Ct value of 20.0, which corresponds to the average Ct for 10.sup.6 copies of standard.
(50) To examine whether similar assay results would be obtained using different lots of the qPCR standard, three additional sets of the synthetic oligos were obtained and new standard lots were prepared for each set. The four lots of qPCR standard, which included the original lot, were subjected to qPCR and a standard curve plot that consisted of all four lots combined was constructed (
(51) Finally, assay reproducibility was verified for an evidentiary type specimen. Various dilutions of a buccal swab DNA extract were assayed five times over a thirteen day period and the results demonstrated high consistency of data, both within assays and between assays, for most dilutions (