Eunhye Lee,
Hyun Jeong Lim* and
Ahjeong Son*
Department of Environmental Science and Engineering, Ewha Womans University, 52 Ewhayeodae-gil, Seodaemun-gu, Seoul 03760, Republic of Korea. E-mail: hj.lim77@ewha.ac.kr; ason@ewha.ac.kr; ahjeong.son@gmail.com; Tel: +82-2-3277-3773 Tel: +82-2-3277-3339
First published on 11th June 2025
The complexity and heterogeneity of soil samples necessitate the inclusion of extensive purification steps prior to genomic assays, such as quantitative PCR (qPCR). Although conventional DNA extraction kits have notably enhanced the convenience of the process, those designed for soil vary considerably in terms of reagents, time, and equipment. Therefore, purified gDNA quality varies depending on the DNA extraction kits used, which leads to discrepancies in gene quantification using qPCR. This issue can be amplified considerably when more complicated (or contaminated) soils are analyzed, even if extensive DNA extraction is employed. Here, we evaluated the influences of the DNA extraction method to the gene quantification using qPCR across soil types. Further Mg2+ ion spiking experiments were performed to observe multiple inhibitory effects on qPCR analysis performance. The results suggest that discrepancies in gene quantification are evident in the presence of qPCR inhibitors in soil samples. Furthermore, discrepancies in quantification results are exacerbated by gDNA template quality, which is attributed to DNA extraction. The observed multiple inhibitory effects underscore the importance of careful consideration of both DNA template quality and soil type to ensure more accurate gene quantification in soils.
The adaptation of qPCR to soil microbiology was initially challenging owing to soil sample complexity and heterogeneity, and physical and chemical property variations. Soils often contain substances that inhibit qPCR, including humic acids, polysaccharides, urea, phenolic compounds, cations, and heavy metals. These substances can remain in genomic DNA via co-precipitation during DNA extraction and inhibit qPCR enzymatic reactions.2–5 For instance, qPCR is inhibited by low levels of humic acid through interactions with template DNA and Taq polymerase6,7 as well as in the presence of Mg2+ ions, even after extensive purification.8 The Mg2+ ion is a part of the qPCR reagent and functions as a cofactor of Taq polymerase during PCR. However, excess Mg2+ can interact with DNA and inhibit coagulation-based inhibition.9,10 Therefore, DNA purification is necessary to remove these substances from soil samples and circumvent this potential inhibition.11–15 Current DNA extraction kits for soil comprise multiple process steps as well as the use of multiple reagents, and are suitable even for the most delicate samples.
Although conventional DNA extraction has notably enhanced convenience for decades, methods designed for soil vary considerably in terms of reagents, time, and equipment. Therefore, the quality of purified gDNA varies depending on the DNA extraction kit used, leading to discrepancies in gene quantification using qPCR. This issue can be amplified when complicated or contaminated soils are analyzed, even if extensive DNA extraction is employed. Such discrepancies compromise the performance of qPCR analysis and are not easily circumvented, as each case varies. If the PCR assay discrepancy appears to be a false negative due to inhibition, the issue can become more serious (e.g., detection of pathogens). Therefore, we revisited and evaluated this issue to identify any potential assay discrepancies and developed remedial actions for qPCR assays in soils.
This study aimed to evaluate the impact of soil type and DNA extraction kit on the bacterial gene quantification using qPCR. The soils were sampled to assess various physicochemical properties, including soil texture, organic content, and Mg2+ ion content. We then investigated the quantification of the target gene from DNA templates extracted from Pseudomonas putida–seeded soils—a model bacterium commonly found in contaminated environments—using different DNA extraction kits. Further Mg2+ ion-spiking experiments were performed to observe multiple inhibitory effects on the performance of the qPCR assay.
Soil properties—including moisture content, total organic content, pH, texture, and cation-exchange capacity—were determined. The soil moisture content (5 g) was measured using a protocol suggested by the soil pollution process test standards.16 The total organic matter content of the soil was measured using the loss of ignition method in a furnace (SH Scientific, Sejong, Korea) for combustion at 550 °C for 4 h. The total organic content (%) was calculated by dividing the weight loss upon ignition by initial dry weight. All sieved soils (5 g) were subjected to pH measurements after adding 25 mL of deionized water and stirring at ambient temperature for 1 h. The pH was measured using a PB-10 Basic Benchtop pH meter (Sartorius, Göttingen, Germany). The Korean Society of Forest Environment Research (Namyangju, Gyeonggi-do, Korea) measured cation exchange capacity and soil texture.
The Mg2+ ions were extracted from the soil samples, and their concentration was measured using inductively coupled plasma mass spectrometry (ICP-MS, 7850, Agilent Technologies, Santa Clara, CA, USA). The soil sample (5 g) and 10 mL of 0.2 M hydrochloric acid (ACS reagent, 37%, Sigma-Aldrich, St. Louis, MO, USA) were placed in 40 mL amber glass bottles with a Teflon™ lining. The samples were then incubated overnight at ambient temperature (25 ± 1 °C) in an orbital shaker (SHO-2D; Daihan Scientific, Wonju, Kangwon-do, Korea).17 Impurities were removed from the extract using a 0.2 μm syringe filter (PTFE, ADVENTEC Co., Tokyo, Japan). The Mg2+ ion standard (1000 mg L−1, Merck, Darmstadt, Germany) solution was prepared at 0, 0.5, 1, 1.5, 2, 2.5, 5, and 10 mg L−1 by serial dilution with 0.2 M hydrochloric acid. The detailed parameters for the ICP-MS analysis are listed in Table S2.†
DNA extraction from cell-seeded soil samples was conducted using three commercial kits—Kit A: Power soil pro kit (QIAGEN, Venlo, Netherlands), Kit B: FastDNA™ SPIN Kit (MP Biomedicals, Santa Ana, CA, USA), and Kit C: Nusoil kit (MACHEREY-NAGEL, Duren, Germany). Subsequent DNA extraction was performed according to the protocols provided by the three DNA extraction kits which included cell lysis, DNA binding, washing, and elution steps. Kit A is a DNA extraction kit that follows the most basic protocol, which includes two washing steps using one spin filter. Kit B allows for relatively fast DNA extraction compared to that of the other kits and requires only one washing step using one spin filter. In addition, during the elution step, an incubation step is performed at a temperature higher than ambient temperature (25 ± 1 °C). Conversely, Kit C includes a filter to remove additional inhibitors following cell and before DNA binding. A singular spin filter involves four washing steps. The detailed features of each kit are listed in Table 1.
DNA extraction process | Kit A | Kit B | Kit C |
---|---|---|---|
Cell lysis step | Chemical and physical lysis using beads | Chemical and physical lysis using beads | Chemical and physical lysis using beads |
Inhibitor removal step | Remove inhibition factors using chemical precipitation | — | Remove inhibition factors using an inhibitor column filter |
DNA binding step | Chemical DNA binding | Chemical DNA binding | Chemical DNA binding |
Washing step | Two-step washing with ethanol-based washing reagent | One-step washing with ethanol-based washing reagent | Four-step washing with ethanol-based washing reagent |
DNA elution step | Elution at room temperature (18–25 °C) | Elution at high temperature (55 °C) | Elution at room temperature (18–25 °C) |
![]() | (1) |
The qPCR was performed using the Applied Biosystems QuantStudio 1 Real-Time PCR System (Thermo Fisher Scientific). Primers targeting a 306-bp region of the PAH-RHDα gene were synthesized by Bioneer (Daejeon, Korea; Table S3†). The qPCR was performed via a denaturation process at 95 °C for 15 min, followed by annealing at 94 °C for 15 s, 57 °C for 30 s, and 72 °C for 30 s for 40 cycles. The qPCR reaction (total 25 μL) composition was as follows: SYBR Green Master mix (2X; Thermo Fisher Scientific) 12.5 μL, forward primer (10 μM) 1 μL, reverse primer (10 μM) 1 μL, UltraPure™ DNase/RNase-free distilled water (Invitrogen) 5.5 μL, and template gDNA 5 μL. Afterwards, the elongation process was carried out at 80 °C for 10 s and 72 °C for 7 min.
To determine the inhibitory effect of Mg2+ ion content on the qPCR calibration curve, Mg2+ ions were added to the qPCR standard solution. Five points five microliter of diluted Mg2+ ion solutions (22 and 87 mg L−1 Mg2+ ion content) were added to the qPCR reaction, instead of deionized water, to achieve 5 and 20 mg L−1 Mg2+ ion content, respectively, resulting in a total volume of 25 μL reaction (equivalent to 0.1 and 0.5 μg Mg2+ in a reaction, respectively). The same volume (5.5 μL) of deionized water was subjected to the qPCR reaction as a negative control (0 mg L−1 Mg2+ ion content). The qPCR temperature program was as described above. Amplification and calibration curves were obtained and replotted using SigmaPlot software (12.5 version, Systat Software, Inc., San Jose, CA, US).
To examine the effect of a high Mg2+ ion concentration solely on the qPCR analysis, 2.4 mL of Mg2+ ion stock solution (120 mg L−1) was added to the pure cell culture (9.6 mL, OD ∼0.1) to achieve a final concentration of 24 mg L−1 (equivalent to 0.288 mg Mg2+ in a sample). Furthermore, the effect of a high Mg2+ concentration in the soil matrix was assessed by choosing sandy soil as the sample matrix among the four soil samples. A high Mg2+ ion concentration (325 μL of 120 mg L−1 stock solution, 78 mg kg−1 as the final concentration, equivalent to 0.039 mg Mg2+ in a sample) was spiked to 500 mg of the sandy soil to increase the Mg2+ ion content, while ensuring that the final concentration remained below 500 mg kg−1, consistent with relevant environmental conditions.18–20 The calculated final Mg2+ ion concentrations in the sandy soil are plotted in Fig. S6.† Notably, the same number of bacterial cells (9.6 mL, OD ∼0.1) were seeded in the sandy soil. For the negative control, deionized water (2.4 mL for the pure cell sample and 325 μL for the soil sample) was added instead of Mg2+ ions.
The mixture was subjected to three different DNA extraction kits (Section 2.3). Subsequently, a qPCR analysis was performed as described in Section 2.4. The qPCR results of the Mg2+ ion-spiked pure cell cultures and Mg2+ ion-spiked sandy soils were compared.
The total organic content of the soil samples was measured to determine the amount of organic matter, which is a representative inhibitor of qPCR analyses (Fig. 1b).4,6,7,21 The results showed that Mountain contained the highest total organic content (9.46 ± 0.16%) followed by 4.57 ± 0.15% in paddy, 1.18 ± 0.05% in river, and 2.01 ± 0.81% in sandy soil.
Mg2+ ion concentrations in the soil samples were measured using ICP-MS and quantified using a calibration curve (Fig. S5†). The results are presented in Fig. 1c and Table S1.† River contained the highest Mg2+ ion concentration of 287.52 ± 2.98 mg kg−1, followed by Paddy (272.47 ± 2.62 mg kg−1), with mountain and sandy soils showing a relatively lower Mg2+ ion content at 138.12 ± 1.79 and 95.75 ± 0.77 mg kg−1, respectively. Additionally, we assessed other physicochemical parameters to understand soil properties (Fig. 1d). All four soil samples had slightly acidic pHs ranging from 4.50 to 5.59, with moisture contents ranging from 2.24 to 16.30%. They had varying soil cation exchange capacity values of ∼<13 cmolc kg−1 (12.54 cmolc kg−1 in mountain; 3.85 cmolc kg−1 in paddy; 2.86 cmolc kg−1 in river; and 0.77 cmolc kg−1 in sandy soil).
The three-dimensional plot of Fig. 2 presents the qPCR result of the target gene (PAH-RHDα) of P. putida-seeded (all with the same amount of bacterial cells: 9.6 mL of OD 0.1) soil samples (i.e., Mountain, Paddy, and River), which were extracted using three different DNA extraction kits. We observed a substantial discrepancy in the qPCR results depending on soil type with different organic matter content. Moreover, the discrepancy in gene quantification was pronounced depending on the type of template DNA extracted. The target gene in the DNA extracted by Kit A from Paddy and River was successfully quantified (2.5 × 107 and 7.9 × 107 gene copies, respectively) via qPCR, whereas that in Mountain was not detected (depicted as N.D. in Fig. 2). The gene in the extracted DNA from Mountain was also not determined (depicted as N.D. in Fig. 2) when it was treated by Kit B. Moreover, the target gene in the extracted DNA using Kit B from Paddy was close to N.D. (in the order of 10−1 gene copy number). Those from River were underestimated by two orders of magnitude compared with those of Kits A and C. In contrast, the target gene in the extracted DNA from the three soils treated with Kit C could be quantified by qPCR, even though Mountain presented a slight underestimation (three orders of magnitude). This observation suggests that substantial discrepancies were found in qPCR results depending on soil type and DNA extraction kit used. Further investigation is required to understand the influence of template DNA quality on qPCR results.
To elucidate the relationship between the discrepancy in qPCR results and type of DNA extraction (Fig. 2), it is important to understand the differences in DNA template qualities resulting from the selected DNA extraction kits. Therefore, we first examined DNA recovery, purity, and Mg2+ ion removal in bacterial cell-seeded soils using these kits (Fig. 3). The overall DNA concentrations were observed within the range of 7–24, 7–125, and 10–24 ng μL−1 for Kits A, B, and C, respectively (Fig. 3a). Kits A and C exhibited similar DNA yields in the Mountain (t-test, P-value > 0.770) and River (t-test, P-value > 0.734) samples, whereas Kit B presented an erroneously large variation in DNA extraction for all soil types (t-test, P-values ≪ 0.001 when comparing Kit B with Kits A and C, indicating statistically significant differences). This result was consistent with previous studies that reported significant variations in DNA yield depending on the extraction kit used.22–25
To further evaluate the variation in the three kits, the purity of the extracted DNA was examined based on 260/280 and 260/230 ratios, which generally indicate the level of contamination of extracted DNA by protein and organic matter, respectively. According to the NanoDrop™ 2000 instructions and,26 it is recommended to obtain the 260/280 ratio in a range similar to 1.8 and 260/230 ratio in a range of 2.0–2.2 in well-purified DNA samples. As shown in Fig. 3b, DNA templates were successfully purified from Mountain, Paddy, and River by Kits A and C while avoiding protein contamination (overall ratio of 1.9–2.0), whereas the extracted DNA by Kit B presented a slightly lower 260/280 ratio of 1.4–1.7 in the three soils. In comparison, the range of 260/230 ratios differed remarkably for each kit (Fig. 3c, ANOVA test, P-value ≪ 0.001 for all soil type), with the lowest values for Kit B (<0.3). Although the DNA template extracted using Kit C exhibited higher 260/230 ratios (>0.9), all template DNA samples seemed to contain residual organic matter, resulting in low 260/230 ratios. These results indicate that all three kits struggled to remove organic matter in the soil samples to a level sufficient for qPCR, despite the common belief that they should have shown some level of bias among the kits.
The Mg2+ ion removal efficiency also varied according to the kit type (Fig. 3d), but we did not observe any notable differences between soil samples using the Kit A and C (ANOVA test, P-value = 0.258 and 0.991, respectively). This indicates that Mg2+ ion contents in the soil samples were within the removable ranges determined through the purification steps suggested by each kit. Kit B exhibited relatively lower and varied Mg2+ ion removal efficiency (ANOVA test, P-value ≪ 0.001) with residual Mg2+ ion concentrations of ∼5.6–6.1 mg L−1 in the extracted DNA from the three soils, whereas Kit C presented marginal residual Mg2+ ion concentrations of <1.1 mg L−1 in all samples.
This can be explained by comparing the kit protocols (Table 1). The variation in the purity results of the 260/230 ratio and Mg2+ ion removal may have occurred because of the purification differences. The absence of the inhibitor removal and simplified washing steps (i.e., one step) in Kit B resulted in a large deviation in the DNA recovery and purification of bacterial debris, organic matter, and Mg2+ ions in the soil matrix. In contrast, inhibitor column filtration and intensive purification steps (i.e., four steps) in Kit C resulted in a more thoroughly purified DNA template from the soil samples. Referring to the qPCR results (Fig. 2), differences in the quality of the DNA template influenced by its extraction caused a considerable discrepancy in quantifying the target genes.
The Mg2+ ion causing qPCR inhibition was first determined by observing the amplification plots and Ct values constructed by using the PAH-RHDα gene standards. Fig. 4a–c presents the amplification plots of the PAH-RHDα gene with varying amounts of Mg2+ (i.e., 0, 0.1, and 0.5 μg Mg2+ per reaction, equivalent to 0, 5, and 10 mg L−1 Mg2+, respectively) in the standard plasmid for qPCR. The signal-interfered amplification plot was observed in the presence of 0.5 μg of Mg2+ ions (Fig. 4c) in the plasmid template solution, exhibiting non-replicability and low ΔRn. In comparison, the PAH-RHDα gene was successfully amplified in the absence of Mg2+ ions (Fig. 4a) and presence of 0.1 μg of Mg2+ ions (Fig. 4b), and was well-replicated and had high ΔRn values. Since the poor quality of amplification plots with low ΔRn values refers to the occurrence of inhibition during the PCR reaction, this result infers that the gene quantification can be hindered by an excess amount of Mg2+ ions (e.g., 0.5 μg Mg2+ ions per qPCR reaction solution).
We compared the Ct values from each amplification plot with the gene copy number (Fig. 4d). In the presence of 0.1 μg Mg2+ ions, the qPCR assay showed similar Ct values to that without Mg2+ ions across the entire gene copy number range assessed (from 100 to 108). This result also suggests the negligible inhibitory effects on the gene quantification caused by 0.1 μg of Mg2+ ions. However, in the presence of 0.5 μg of Mg2+ ions, the overall Ct values were higher than those with or without 0.1 μg of Mg2+ ions. Moreover, a higher Ct value was observed when the gene copy number was <104. Thus, the lower limit of detection became two orders of magnitude higher than that without Mg2+ ion inhibition when the qPCR reaction solution contained 0.5 μg of Mg2+ ions. Consequently, we confirmed that the inhibition was caused by an excessive amount of Mg2+ ions (0.5 μg of Mg2+ ions per reaction solution, equivalent to 20 mg L−1 Mg2+ ions in the solution) disrupting the calibration curve of the qPCR assay under the reaction conditions in this study, which was far less than the optimum Mg2+ ion concentration for PCR (i.e., 2–5 mM in a reaction, equivalent to 49–122 mg L−1).29
![]() | ||
Fig. 5 Comparison of gene copy number of extracted DNA templates from Mg2+ ion-spiked (a) pure cells and (b) sandy soil sample using different DNA extraction kits. |
The target gene in sandy soil with a high Mg2+ ion content was hardly quantified using the DNA extraction kits (Fig. 5b), even though the spiked amount of Mg2+ in sandy soil was 10 times lower than that in the pure cell sample. The DNA template prepared using Kit A exhibited a slightly lower gene copy number than the mean gene copy number obtained from the Mg2+ ion-spiked pure cell culture when extracted from the soil matrix. Moreover, the target gene in the DNA extracted from sandy soil using Kit B could not be detected by qPCR. However, the DNA extracted using Kit C only successfully demonstrated target gene quantification using the soil matrix, even with a high Mg2+ ion content, presenting a gene copy number similar to that of the Mg2+ ion-spiked pure cell culture. Although sandy soil has a simple composition without critical inhibitors, the matrix effect of the soil may notably affect the quality of the extracted DNA templates. In sand containing less than 100 mg kg−1 of Mg2+ ion (equivalent to 0.1%) as a single inhibitor, a soil-specific DNA extraction kit was sufficient to obtain DNA acceptable quality for qPCR.8 However, in actual sandy soil with a more complex composition, discrepancies in qPCR results appeared to be magnified depending on the type of DNA extraction kit used.
Overall, an apparent discrepancy in the qPCR results was observed when the sample matrix contained a high inhibitor content in soils compared to that of the pure cell culture. In addition, the different DNA template extraction procedures further exacerbated the inconsistency in the qPCR results of the soil samples, consistent with the results shown in Fig. 2 and 3.
An exacerbated discrepancy in target gene quantification was observed in environmental soil samples (i.e., Mountain, Paddy, and River) with high Mg2+ ion-spiking. As shown in Fig. 6, the target gene in the DNA templates extracted from Paddy and River were successfully quantified despite the high Mg2+ ion concentration, which reached final levels of 350.47 and 365.52 mg kg−1, respectively, due to Mg2+ addition (Fig. S6†). However, in Mountain with Mg2+ spiking, the target gene was quantified at a 2.3 × 102 gene copy number, which is reduced by two orders of magnitude compared to that in Mountain without Mg2+ ion-spiking (Fig. 2). This reduction may have resulted from an antagonistic interaction between the high organic matter content of Mountain and elevated Mg2+ ion concentration, which reached a final level of 216.12 mg kg−1 due to spiking (Fig. S6†).
![]() | ||
Fig. 6 Comparison of the gene copy number of the extracted DNA template by Kit C from soil samples with and without Mg2+ ion-spiking. Note that the qPCR results from the soils without Mg2+ ion-spiking are replicated in Fig. 2. |
These results suggest that discrepancies in gene quantification are evident in the presence of qPCR inhibitors. Furthermore, discrepancies in the quantification results are exacerbated by the DNA template quality, which is attributed to DNA extraction. The possibility of multiple inhibitory effects caused by soil properties and DNA extraction kits was also observed.
This study highlights the importance of selecting appropriate DNA extraction kits based on soil characteristics. For fine-textured soils with high organic content, kits equipped with effective metal ion inhibitor removal steps such as chelation are recommended. In cases where soils have relatively low or unpredictable levels of inhibitors, kits that offer a balanced inhibitor removal and DNA yield, along with versatile extraction protocols, are preferable to ensure reliable qPCR results. If significant inhibition persists despite using a high-performance kit, further dilution or post-purification of the DNA template should be considered to improve qPCR performance.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d5ra02689j |
This journal is © The Royal Society of Chemistry 2025 |