|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Molecular Diagnostics and Genetics |
1 BioAnalytical Innovation Team, LGC Ltd., Teddington, United Kingdom.
2 Experimental Fetal Medicine Group, Institute of Reproductive & Developmental Biology, Imperial College London, and Centre for Fetal Care, Queen Charlottes & Chelsea Hospital, Hammersmith Campus, London, United Kingdom.
aAddress correspondence to this author at: BioAnalytical Innovation Team, LGC Ltd., Queens Road, Teddington TW11 0LY, United Kingdom. Fax 44-20-8943-2767; e-mail jacquie.keer{at}lgc.co.uk.
| Abstract |
|---|
|
|
|---|
Methods: Blood samples were collected from 201 women between 5 and 41 weeks of pregnancy. Quantitative PCR was used to assess total and fetal DNA concentrations, and allelic discrimination analysis was investigated as a route to detecting specifically fetal DNA.
Results: Male fetuses were detectable from 5 weeks amenorrhea with increasing fetal DNA concentrations across gestation. The sensitivity of fetal male gender determination in pregnancies with live birth confirmation was 99%, with 100% specificity. Total DNA concentrations did not correlate with gestational age, but appeared slightly higher in the first and third trimesters than in mid-pregnancy. Analysis of short tandem repeats demonstrated that significant improvements in the detection limit are required for specific detection of fetal DNA.
Conclusions: The high sensitivity of PCR-based detection, together with quantification provided by real-time DNA analysis, has clear potential for clinical application in noninvasive prenatal diagnosis. However, accurate quantification using best-fit data analysis, standardization of methods, and performance control indicators are necessary for robust routine noninvasive diagnostics.
| Introduction |
|---|
|
|
|---|
Lo et al. first quantified fetal DNA in women tested during both the first and third trimesters of pregnancy, the earliest detection of fetal DNA being at 11 weeks of gestation. Most subsequent investigations have focused on a particular trimester (15)(16), and more recently, fetal DNA has been detected as early as 5 weeks of gestation (17)(18). However, few groups have monitored the concentrations of both fetal and total DNA in maternal plasma across all three trimesters or included the steps taken to ensure accurate quantification.
The way in which detected DNA is quantified has generally not been detailed in these reports, although there are brief references to methodology (15)(19), typically describing the use of serially diluted human genomic DNA as a calibrant for unknown samples. Some studies have provided Pearson correlation coefficients (R2) for calibration curves (15)(19), but few specified the methods used to prepare calibration curves, the number of replicate points in each assay, or how consistency was monitored between assays. Because quantitative analysis by real-time quantitative PCR (Q-PCR)1 is only as accurate as the calibration, the method used for DNA dilution series preparation is central to assay reliability. This is especially important at low concentrations, such as first trimester analyses, where sampling effects relating to the stochastic distribution of the molecules may noticeably affect the accuracy of quantification (20)(21). Here we address the accuracy of quantification and the variation in DNA concentrations throughout gestation.
| Materials and Methods |
|---|
|
|
|---|
processing of blood samples
EDTA blood was processed at the Institute for Reproductive and Developmental Biology immediately after collection (with a few exceptions, which were stored overnight at 4 °C). In most cases, the exact volume of blood collected was recorded. Plasma was prepared by centrifugation at 1600g for 10 min (Sorvall RC 5C Pluss). The plasma fraction was carefully transferred to a clean tube (volume again recorded) and frozen at 80 °C; it was then sent frozen by courier to LGC, where it was stored at 80 °C until further processing. On the day of analysis, the plasma was thawed and subjected to further centrifugation at 16 000g for 10 min (Jouan MR22i) to pellet any remaining cellular debris. The plasma was then carefully transferred to a clean polypropylene tube. All blood processing was performed in a class II safety cabinet (MDH).
dna extraction from maternal plasma
DNA was prepared from 1.5 mL of plasma by use of a QIAamp Blood Mini Kit (Qiagen) according to the "blood and body fluid protocol". The extraction method was performed as recommended by the manufacturer, except that all reagent volumes were scaled up to accommodate the increased volume of plasma used in each extraction. The sample was loaded on the spin-column by repeated centrifugation, and a vacuum line was used to aspirate the filtrate, minimizing contamination. A final elution volume of 50 µL was used. The extraction efficiency was determined by adding known quantities of male human genomic DNA to nonpregnant female plasma samples and calculating the recovery by Q-PCR. The efficiency ranged from 45% to 75%, with some losses possibly caused by reduced elution volume, leading to diminished yields as stated in the manufacturers guidelines.
generation of calibration curves
A human male control DNA (Promega) was used to prepare calibration curves for the chosen amplification targets. A PicoGreen® (Molecular Probes) quantification analysis was performed with a purified DNA calibrator (Cambio), prepared gravimetrically, to determine the concentration of the stock DNA. This value was converted to genome-equivalents/mL by use of the conversion factor of one diploid genome being equivalent to 6.6 pg of DNA. High-concentration (5000 genome-equivalents/µL) DNA stocks were prepared after the initial PicoGreen analysis, aliquoted into small volumes, and stored at 4 °C. For each run, one aliquot of stock DNA was serially diluted with molecular biology grade water (Eppendorf) to construct a nine-point calibration curve in the range 0.52500 genome-equivalents/µL.
real-time q-pcr
Real-time Q-PCR was performed with the 5'-nuclease assay. Primers and fluorogenic 5' exonuclease (TaqMan) probes were designed according to the guidelines recommended by Applied Biosystems. A coding region of the male-specific SRY gene was chosen to monitor the presence of fetal DNA, and a coding region of the housekeeping gene GAPDH was used to quantify total DNA. Both targets were detected independently in singleplex reactions to prevent outcompetition of the SRY amplification (data not shown). Reaction conditions for both targets were as follows: 450 nM forward primer, 450 nM reverse primer, 225 nM TaqMan probe, 1x PCR mixture (BioGene), and 5 µL of genomic DNA extract in a final reaction volume of 25 µL. PCR was performed for 55 cycles (uracil N-glycosylase activation at 50 °C for 2 min and denaturation at 95 °C for 10 min, followed by 55 cycles of 95 °C for 5 s and 60 °C for 1 min). The gene-specific primer and probe sequences were as follows:
SRY: Forward primer, 5'-CGATCAGAGGCGCAAGATG-3'; reverse primer, 5'-TGGTATCCCAGCTGCTTGCT-3'; probe, 5'-VIC-TCTAGAGAATCCCAGAATGCGAAACTCAGAGA-TAMRA-3' (where TAMRA is 6-carboxytetramethylrhodamine)
GAPDH: Forward primer, 5'-AGGTTTACATGTTCCAATATGATTCCA-3'; reverse primer, 5'-ATGGGATTTCCATTGATGACAAG-3'; probe, 5'-FAM-CCGTTCTCAGCCTTGACGGTGC-TAMRA-3' (where FAM is 6-carboxyfluorescein)
All probes and primers were obtained from Sigma-Genosys with the exception of the SRY VIC-labeled probe, which was supplied by Applied Biosystems.
Amplification data were collected and analyzed with an ABI Prism 7700 Sequence Detector (Applied Biosystems. The cycle threshold value (Ct) was the measurand in all cases. Each sample was analyzed in triplicate, and multiple negative reaction blanks were included in every analysis for both sample extraction and amplification stages. A calibration curve (triplicate samples) was analyzed on the same reaction plate for each run. Fetal gender of each sample was confirmed phenotypically at birth or by fluorescence in situ hybridization analysis of fetal material. The study was performed blind at LGC; none of the confirmed fetal genders were disclosed by the Institute for Reproductive and Developmental Biology until after the analysis was completed and the results were recorded.
quality control
All reagents for DNA extraction and Q-PCR were prealiquoted into one-shot volumes. Filter tips were used throughout, and in addition 1-mL Reach filter tips (Molecular BioProducts) were used when pipetting from 15-mL tubes. All equipment and work areas were swabbed with fresh 10 mL/L hypochlorite, and DNA extraction was performed in a separated area in a class II hood. A vacuum line was used during the extraction process to minimize generation of aerosols. Ensuring that all lids were tightly closed during centrifugation significantly reduced contamination. During Q-PCR, carryover contamination was prevented by the use of uracil N-glycosylase. Reaction mixtures were prepared in a separated area in a hood (CleneCab; Herolab), and DNA was added to the reaction in a dedicated template addition area.
Both extraction and amplification negative controls were checked for positive signals, which if present would indicate contamination occurring during the process. In the event of contamination, either the extraction and/or the amplification was repeated for the potentially affected batch of samples. The reliability of quantification was assessed by monitoring of the calibration curve for each individual amplification run.
analysis of short tandem repeats
We analyzed five pairs of age-matched samples (gestational ages, 29 weeks/6 days to 41 weeks) from confirmed male and female gender pregnancies as well as artificial mixtures of male and female DNA in ratios of 1:2 to 1:660, respectively, using an ampFLSTR® SGM PlusTM Kit (Applied Biosystems). With the ampFLSTR SGM Plus, we could amplify 11 different loci, but only results for the Amelogenin amplicons are presented here. Amplification was performed according to the manufacturers protocol with the following modifications: 5 µL of DNA extract plus 15 µL of 1x Tris-EDTA buffer was used in each reaction to make the volume up to 20 µL. Positive and negative controls were included with every run. Results were analyzed with GeneScan® 3.1.2 and Genotyper® 2.5 software (Applied Biosystems).
statistical analysis
We used Microsoft Excel to perform a factorial ANOVA to test for significant interrun differences in the calibration curve data from all experiments. We used an empirical method to analyze the data, using quasi-Newton exponential weighted regression because this was determined to be the most appropriate model describing the relationship between the target concentration and Ct value for both GAPDH and SRY calibration curves (English et al., Improving the accuracy and reliability of low level DNA analysis using real-time PCR, submitted for publication). We used this model to quantify the amount of each target in the unknown samples, performing quasi-Newton weighted analysis with the software package STATISTICA 6 (StatSoft Ltd). Because the volume of blood taken and plasma prepared was recorded for almost every sample, we were able accurately to calculate the amount of each target gene present in the samples in terms of genome-equivalents/mL of blood, using the following equation:
![]() |
We analyzed the relationships of SRY and GAPDH concentrations with gestational age and GAPDH concentration with maternal age by linear regression using the least-squares method; we also used a Pearson test to establish any correlation (Microsoft Excel). We used STATISTICA 6 to calculate the 95% prediction intervals for SRY concentrations against gestational age from the linear regression equation. Reference intervals were derived for both circulating fetal and total DNA concentrations in each trimester based on 95% prediction intervals obtained by use of the residual standard deviation from the regression analysis (22)(23).
| Results |
|---|
|
|
|---|
|
correlation of fetal dna (sry) copy numbers in maternal plasma with gestational age
The analysis gave a range of quantifications for SRY from 0 to 507.5 genome-equivalents/mL of blood, with a clear difference between positive signals for SRY, which indicated the presence of a male fetus, and the absence of amplifiable SRY, which signified a female fetus. No contamination events were detected, and the cutoff between the two populations was 1 genome-equivalent/mL. Either live birth examination or fluorescence in situ hybridization analysis on invasively sampled fetal material was used to confirm fetal gender. The sensitivity of the SRY assay was 99% in the 183 samples for which confirmation of fetal sex could be obtained. A single SRY false negative, at 7 weeks of gestation, was the male of a dizygotic twin pregnancy. The range of SRY concentrations in each trimester is represented as a box and whisker plot in Fig. 2A
, and the log mean value of each set of triplicate results per sample is plotted against gestational age in Fig. 3A
. In the male samples, there was a strong positive correlation (R2 = 0.52; y = 0.0051x + 0.6009) between log SRY concentration and gestational age, indicating that >50% of the variation was attributable to this factor.
|
|
We calculated the 95% prediction intervals for SRY (Fig. 3A
); these intervals predict with 95% probability the upper and lower values for an individual sample at a given time-point. We also calculated indicative time-specific reference intervals for both SRY and GAPDH for each trimester, using the 95% prediction intervals derived from the residual standard deviation in the regression analysis; these reference intervals are shown in Table 1
together with the mean and median values for each trimester.
|
effect of gestational age and maternal age on total dna in maternal plasma
A composite calibration curve, including the calibrators from each GAPDH assay, was used to calculate the amount of DNA present in each sample. During this study, one batch of nine samples was reextracted because the extraction negative controls showed low positive signals for GAPDH. In addition, one PCR run of 18 samples was repeated because positive GAPDH signals in the PCR negative controls indicated that contamination had occurred. The repeat results were comparable to the initial data in both cases. Box and whisker plots of the data, shown in Fig. 2B
, illustrated a noticeably increased proportion of samples with increased GAPDH concentrations in the first and third trimesters compared with mid-pregnancy.
A plot of the log GAPDH concentrations against gestational age (Fig. 3B
) showed that there was no significant relationship (R2 = 0.06). Similarly, there was no significant correlation between maternal age and log GAPDH concentrations (R2 = 0.02). The mean and median trimester values are displayed in Table 1
together with the indicative reference intervals for total DNA concentrations, which were calculated using 95% prediction intervals based on the residual standard deviation from the regression analysis.
positive control for fetal dna
We investigated the potential of using short tandem repeat (STR) analysis to provide positive confirmation of sufficient extracted fetal nucleic acid to be detectable by fluorogenic PCR analysis. STR profiling of plasma samples from known male and female pregnancies, matched by gestational age, and of artificial male/female DNA mixtures was performed as described. Electropherograms of the products of amelogenin-specific PCR of the STR multiplex revealed peaks for the 106-bp amplicon from the X chromosome and a 112-bp amplicon from the Y chromosome, if present (Fig. 4
). The heights of the fluorescent peaks indicated the amount of PCR product produced and thus of the amount of target material present in the sample. The in-house mixtures showed a detectable male signal only when the male DNA was present in a 1:2 ratio with female DNA. In mixtures containing lower proportions of male target, the 112-bp male amplicon was not detectable. In the profiles from very late pregnancy samples (
38 weeks), a small minority peak detectable in the male samples at 112 bp was produced from the paternally inherited Y-chromosome amelogenin locus (indicated by arrows in Fig. 4
). In the later pregnancy samples, both the total amount of SRY present and the ratio of fetal to total DNA in the sample, as assessed by Q-PCR, were increased considerably.
|
| Discussion |
|---|
|
|
|---|
Performing a log transformation of the data (Fig. 3A
) allows the application of normal statistics to heteroscedastic data sets. The increase in SRY concentration with gestational age up to 41 weeks (R2 = 0.52 by linear regression analysis) is clearly demonstrated through the complete coverage across the entirety of gestation presented here, in contrast to more limited studies (15)(16). The broad linear relationship between SRY concentrations and gestational age is not as strong in the very early and late stages of pregnancy. Biological variation in individuals relating to placental size, amounts of DNA shed, fetal cell trafficking, and apoptosis may all contribute to the variation observed (24)(25). The lowest concentration detected was 1.2 genome-equivalents/mL of blood at 37 days and increased to 507.5 genome-equivalents/mL of blood at 291 days, which are within the calculated reference intervals shown in Table 1
. Toward the very late stages of pregnancy, sharp increases in circulating fetal DNA concentrations can occur as a result of placental breakdown (26), which may explain the very high SRY concentration observed in the single sample taken at 41 weeks of gestation. The earliest detection was at 37 days, which demonstrates a sensitivity equivalent to that described previously (17)(18). Overall the fetal DNA concentrations reported here are in broad agreement with previous studies for which the trimester data can be directly compared (2)(17)(27). Maternal age did not correlate significantly with SRY concentrations, concurring with previous findings (28). We also did not observe a significant relationship between total and fetal DNA concentrations (R2 = 0.16), in contrast to an earlier report (16).
Recent concerns that persistence of fetal cells from previous pregnancies may lead to false positives (29) are not supported by our results. Twenty-seven of the samples in this study group (13.4%) were female fetuses from mothers who had previously carried at least one male fetus. No SRY signal was detected in any of these samples, indicating that no detectable male DNA persisted from the previous pregnancies. Subjecting samples to centrifugation at 16 000g before analysis, as performed in the present study, has been reported to remove any potential cellular material that could contribute to the free-floating DNA signal after extraction (30)(31), explaining the lack of persistence observed.
We used weighted quasi-Newton exponential fit to the pooled calibration curve data to quantify GAPDH in the unknowns. We observed no strong correlation between variations in the concentration of total DNA and gestational or maternal age, or fetal DNA concentrations. However, the highest concentrations appear to occur at the beginning and toward the later stages of pregnancy (Fig. 2B
and Table 1
), although both the highest and lowest GAPDH values were observed in samples taken in the third trimester. Mean and median values for the second trimester showed a decrease in total concentrations and a tighter distribution of results compared with both the first and third trimesters. Differences in total DNA concentrations are primarily determined by the concentration of maternal DNA and may be attributable to inherent variation between individuals, reflecting lower cell shedding or apoptosis in mid-pregnancy (32). This effect has not been reported previously, but earlier studies either have focused on a single trimester or have not been sufficiently large to observe such trends. Technical causes for this observation are unlikely because the GAPDH quantification performed consistently and samples were analyzed in batches of mixed gestational ages, ensuring that the performance of any individual run would not bias the results of one age group or trimester specifically. Other recent studies are in agreement with the present work, including that of Chan et al. (16), who concluded that ß-globin concentrations (measuring total DNA) and gestational age were not significantly correlated.
A potential route to improving the reliability of noninvasive testing is the use of robust and traceable quantitative methods. Ultimately, standardization of techniques will improve comparability among laboratories and increase confidence in analytical results, and other investigators have addressed the importance of standardization in noninvasive diagnosis. A recent interlaboratory comparison to assess noninvasive male fetal detection (33) highlighted the need for optimized protocols, particularly in the extraction process. The importance of sample preparation, including methods for blood processing, has been the subject of several studies (30)(31). Recently, formaldehyde treatment of blood samples at collection has been advocated to reduce the concentrations of circulating maternally derived nucleic acids by reducing cell lysis (34). Extending clear methodologic guidelines to encompass the entire analytical procedure will improve the consistency and reliability of routine diagnostic noninvasive methods.
Development of a universally applicable positive control for fetally derived DNA would allow better use to be made of negative noninvasive diagnostic results, but in the absence of a positive control, such data must be interpreted cautiously (35). Approaches to providing a positive control include use of distinguishable tracer DNA (18) or detection of fetally expressed genes to confirm the presence of fetally derived nucleic acid in plasma extracts (36). Placenta-specific mRNA expression profiles have recently been reported (37)(38); alternatively, a combination of methylation patterns from imprinting and known inherited polymorphisms can be exploited to identify fetal material specifically (39).
To investigate another potential route, we carried out STR analysis on extracted plasma to detect fetal DNA in the maternal samples by identification of paternally inherited alleles. For most loci, either one or two majority peaks should be seen, corresponding to homozygous or heterozygous maternal targets at each locus. When fetal DNA is present in detectable quantities, an additional minority peak corresponding to the paternally inherited fetal allele may also be seen, if distinguishable by size from the two maternally derived targets. This approach could also be used in detection of expanded repeat mutations. In practice, however, we found that the minority fetal DNA was detectable only when present at high concentrations, as determined by use of artificial mixtures of known amounts of female and male genomic DNA (Fig. 4
). The 112-bp amelogenin male-specific amplicon was detectable only in the 2:1 mixture of female to male DNA, whereas the female 106-bp amplicon was clearly seen in all samples. Informed by the results from the artificial mixtures, we tested third-trimester plasma samples because earlier samples did not have sufficient fetal DNA to allow detection by STR analysis. The samples analyzed ranged from 29 to 41 weeks of gestation, but only a single weak 112-bp amplicon was detected, in the male 41-week sample shown (Fig. 4
).
To exploit the STR approach for noninvasive prenatal analysis, an increase in the sensitivity of the method is required. Potentially this could be achieved either by preferential amplification of fetal DNA or through enrichment strategies during the DNA extraction process (34)(40).
Accurate measurement is underpinned by effective methods for target quantification. There are numerous methods for quantifying DNA, many of which are imprecise. The performance of quantification approaches is ill-defined because of the lack of certified DNA reference materials.
In an attempt to overcome these problems, we used a calibration process traceable to a commercially available quantified human DNA calibration solution (Cambio). This approach has also been reported recently in the quantification of cDNA targets by real-time PCR (41). To minimize interrun variation, we prepared one stock of the human DNA and divided it into aliquots of a high-concentration DNA solution (5000 genome-equivalents/µL), which were stored at 4 °C. All subsequent calibration curves came from a single quantification, and the same analyst was used throughout to reduce variation. The volume used in preparing the dilution series affected the results of the analysis; smaller volumes (2050 µL) led to greater variability in the Ct values than when larger volumes were used (300 µL to 1 mL), with the greatest effect observed at lower target concentrations. In addition, we observed during initial method development that low-concentration dilutions gave a gradually decreasing signal with length of storage (days) at 4 °C, in agreement with other recent findings (42). To minimize variation through this effect, we made fresh dilution series on the day of each run. The calibration curve for each run was used as an effective performance control indicator, providing an effective means for assessing interrun variation and identifying problems with the assay.
A previous study (43) demonstrated that increasing the number of sample replicates from three to six in a 50-cycle real-time PCR improved the rate of detection of male DNA in maternal plasma from 92.6% to 100%. In the present study, the PCR was run for 55 cycles to increase sensitivity sufficiently to detect fetal material in early first-trimester samples. Because of limitations in the volume of DNA extracted, samples were only run in triplicate for both the GAPDH and SRY assays, but this approach was sufficient to achieve 99% sensitivity of fetal detection. Increasing the cycle number from the 3840 cycles recommended by the instrument manufacturers improved analytical sensitivity both by detecting positive signals that were measurable only after >40 cycles of amplification and by increasing our confidence in the distinction between low-concentration positive and true negative results.
In conclusion, these data represent the most complete analysis to date of the concentrations of fetal and total DNA across pregnancy, with coverage from 5 to 41 weeks of gestation. The large sample size permitted clear demonstration of the positive correlation between fetal, but not total, DNA concentrations and gestational age. Although STR analysis has potential as a positive control for fetal DNA, as shown in later pregnancy, the current approach is not sensitive enough for such analyses. An empirical approach using weighted quasi-Newton exponential analysis to increase the accuracy of low-level quantification combined with modification of the Q-PCR method to maximize analytical sensitivity increases confidence in the reliability of quantitative trace molecular measurements. Incorporating modifications such as increasing the number of replicate analyses and amplification cycles as well as anchoring quantitative analysis by use of traceable calibrators, where available, has the potential to dramatically increase confidence in measurements of low-concentration nucleic acids and thus remove barriers to routine clinical application.
| Acknowledgments |
|---|
| Footnotes |
|---|
| References |
|---|
|
|
|---|
The following articles in journals at HighWire Press have cited this article:
![]() |
C. F. Wright and H. Burton The use of cell-free fetal nucleic acids in maternal blood for non-invasive prenatal diagnosis Hum. Reprod. Update, January 1, 2009; 15(1): 139 - 151. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Yin, E.H.Y. Ng, X. Zhang, Y. He, J. Wu, and K.Y. Leung Correlation of maternal plasma total cell-free DNA and fetal DNA levels with short term outcome of first-trimester vaginal bleeding Hum. Reprod., June 1, 2007; 22(6): 1736 - 1743. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. Stanghellini, R. Bertorelli, L. Capone, V. Mazza, C. Neri, A. Percesepe, and A. Forabosco Quantitation of fetal DNA in maternal serum during the first trimester of pregnancy by the use of a DAZ repetitive probe Mol. Hum. Reprod., September 1, 2006; 12(9): 587 - 591. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |