Skip to main content

Advertisement

Optimized DNA extraction from neonatal dried blood spots: application in methylome profiling

Article metrics

Abstract

Background

Neonatal dried blood spots (DBS) represent an inexpensive method for long-term biobanking worldwide and are considered gold mines for research for several human diseases, including those of metabolic, infectious, genetic and epigenetic origin. However, the utility of DBS is restricted by the limited amount and quality of extractable biomolecules (including DNA), especially for genome wide profiling. Degradation of DNA in DBS often occurs during storage and extraction. Moreover, amplifying small quantities of DNA often leads to a bias in subsequent data, particularly in methylome profiles. Thus it is important to develop methodologies that maximize both the yield and quality of DNA from DBS for downstream analyses.

Results

Using combinations of in-house-derived and modified commercial extraction kits, we developed a robust and efficient protocol, compatible with methylome studies, many of which require stringent bisulfite conversion steps. Several parameters were tested in a step-wise manner, including blood extraction, cell lysis, protein digestion, and DNA precipitation, purification and elution. DNA quality was assessed based on spectrophotometric measurements, DNA detectability by PCR, and DNA integrity by gel electrophoresis and bioanalyzer analyses. Genome scale Infinium HumanMethylation450 and locus-specific pyrosequencing data generated using the refined DBS extraction protocol were of high quality, reproducible and consistent.

Conclusions

This study may prove useful to meet the increased demand for research on prenatal, particularly epigenetic, origins of human diseases and for newborn screening programs, all of which are often based on DNA extracted from DBS.

Background

Epigenetic mechanisms, such as DNA methylation, have been suggested as possible causal pathways linking environmental exposure to disease. Many of these studies depend on the epigenome-wide analysis of prospectively collected samples, in the context of large human cohorts. As epigenome-wide technologies are becoming available, the use of such cohort studies will provide large amounts of information in the coming years. Due to the general lack of biospecimen collection in observational human studies, many of these cohorts rely on the use of dried blood spots (DBS) obtained soon after birth as the main source of biological information [1].

The use of filter paper for blood collection and analysis was implemented as early as the 1960s by Guthrie et al. using dried-blood samples for newborn phenylketonuria detection [2]. “Guthrie cards” are widely used in many types of tests, including chemical, serological, and genetic applications [3]. More recently, Flinders Technology Associates chemically treated filter papers (FTA cards) were specifically developed for DNA/RNA analyses [4]. These chemically treated cards allow long-term storage of DNA at room temperature and are impregnated with denaturants that guard against oxidation, nuclease and ultraviolet damage, and both bacterial and fungal degradation.

Neonatal DBS are routinely collected in many countries and represent a cost-effective tool to store precious biological specimens for subsequent studies. However, reliable profiling the DNA methylome in DBS has proven to be technically challenging, particularly because such techniques require stringent bisulfite preprocessing that can degrade DNA [5]. Other limitations of their use include the variable degradation of DNA due to storage and extraction, the usually small amounts of DNA that can be obtained (typical blood spots are between 6 and 10 mm in diameter), and the identification of technical artifacts potentially associated with long term storage [6].

Recently, there has been increasing interest in the use of DBS in DNA methylome analyses, using Methylated DNA Immunoprecipitaion combined with sequencing (MeDIP-seq) [7], Methyl-CpG Binding Domain (MBD) protein-enrichment combined with sequencing (MBD-seq) [8], and Infinium (Illumina) bead arrays [7, 911]. The last version of Illumina’s bead array, Infinium HumanMethylation450 (HM450) Beadchip, is cost-effective, requires DNA amounts as low as 300 ng, enables the detection and quantitation of DNA methylation levels at 486,685 CpG sites across the genome and represents one of the most comprehensive microarray methods to date for investigating the methylome [12]. Three reports have addressed the utility of HM450 on DBS [7, 10, 11]. The first report validated the use and the high correlation between two different methylomic platforms on DBS DNA: HM450 and MeDIP-seq [7]. The second one generated good quality methylome-wide data from DBS, as compared to their matched frozen buffy coat [10]. The third one used DBS-based HM450 analyses to study the epigenetic effects of gestational age, as recently demonstrated by our group [11]. However, none of these studies analyzed methods for optimized DNA extraction and quality verification from DBS, both of which represent major upstream steps in the pipeline for DBS-based research, including epigenetics.

In this report, we tested and developed a range of DNA extraction methods from neonatal FTA cards, individually or in combinations. We incorporated or modified protocol steps that were crucial to increase the DNA yield and quality from DBS, and additionally tested their efficiency on Guthrie cards. Moreover, we suggest an optimal protocol for both, pyrosequencing- and HM450-based, methylation studies. This work could prove useful in meeting the increased demand for research on prenatal origins of human diseases and for newborn screening programs.

Results

Optimization of Phases I and II in the DNA extraction protocols

Limited quantity and quality are important drawbacks in the use of DNA obtained from DBS, particularly for epigenome-wide studies. Initially, we ruled out the possibility of using whole bisulfitome amplification (WGA) after confirming the introduction of biases, mostly in the middle range of DNA methylation levels (Additional file 1), further confirming the recently reported finding by Bundo et al. [13]. Then, to systematically optimize DNA extraction from DBS, we divided the different steps of this process into two phases (Figure 1). Critical steps in Phase I included blood extraction off the filter papers, cell lysis and protease digestion (Figure 1A, left panel). Phase II included DNA precipitation, purification and elution (Figure 1A, right panel).

Figure 1
figure1

Phases and classification of protocols used to extract DNA from DBS. Two sequential phases, each encompassing three steps, are outlined (A) and were optimized in the different protocols or method combinations (B) used to extract DNA from DBS. A spin basket is shown next to Phase I.A and consists of a tube with an embedded perforated basket used to separate blood solutions from the filter papers from which they were extracted. A silica-gel column with a funnel-shape design is shown next to Phase II.E and often used to elute small volumes (5–30 μl), as manufactured by Macherey-Nagel and supplied with the extra-small (XS) versions of NucleoSpin kits (B).

We have previously tested several genomic DNA extraction methods on DBS, including resin-based, lysis-based and magnetic bead-based [5]. Lysis- and bead-based methods were the best, but the latter is not suitable for beadchip methylation profiling, so it was not considered in this study [9]. Among lysis-based methods, several commercially available kits, including QIAamp DNA Micro Kit, GenSolve and NucleoSpin, have been shown to be efficient for DBS DNA extraction [10, 14, 15]. Therefore, we selected these three kits to optimize the two phases of DNA extraction. This optimization involved the combination of the different kits and modifications in several steps of the two protocol phases (Figure 1B) (described in detail in Additional file 2). A combination of GenSolve reagents in Phase I and QIAamp reagents in Phase II (referred to as GQ method) was set as a reference protocol to which other tested methods were compared (Tables 1, 2, 3, 4, 5 and 6). In all pairwise comparisons, two DBS punches from the same DBS were used, and assessment of quantity and quality were initially done using Nanodrop (with 260/280 and 260/230 spectrophotometric ratios as a measure of quality) (Tables 1, 2, 3, 4 and 5). These DBS were obtained from the National Children Study (NCS), USA, and were FTA-type, which preserves well DNA relative to other types of neonatal cards (Methods), hence, allowing comparisons across a wide range of DNA extraction protocols.

Table 1 Combinations of Gensolve and Qiagen protocols for DNA extraction from DBS: GQ, QQ and Qq methods
Table 2 Combinations of GenSolve, Qiagen and NucleoSpin protocols for DNA extraction from DBS: GQ versus GN and Gn methods
Table 3 Combinations of GenSolve, Qiagen and NucleoSpin protocols for DNA extraction from DBS: GQ versus GN-XS and Gn-XS methods
Table 4 Combinations of GenSolve, Qiagen and NucleoSpin protocols for DNA extraction from DBS: GQ versus NN and Nn methods
Table 5 Combinations of GenSolve, Qiagen and NucleoSpin protocols for DNA extraction from DBS: GQ versus NN-XS methods
Table 6 Cross-comparisons of DNA quantity and quality parameters among the different tested DNA extraction protocols

DNA yield and quality were consistently better for the reference GQ protocol when compared to Qiagen protocol (QQ) (Table 1, p < 0.05). Although DNA yield was drastically increased when ethanol was used in the Qiagen precipitation step (Qq, p < 0.001), DNA quality was still suboptimal compared to GQ, as assessed by nanodrop (Table 1). In contrast, the combination of Gensolve and NucleoSpin (GN and Gn) increased the DNA yield while preserving DNA quality, regardless of the use of ethanol in the precipitation step (Table 2, p < 0.05). A similar improvement was observed when using NucleoSpin kit in both phases of DNA extraction (NN and Nn) (Table 4, p < 0.001). The extra-small (XS) versions of NucleoSpin, with column designs specific for low elution volumes, did not consistently improve DNA quantity or quality, whether combined or not with other kits (GN-XS, Gn-XS and NN-XS; Tables 3, 5 and 6), the DNA precipitation buffer changed to ethanol, or the washing volume and frequency increased (Table 3, p > 0.05).

Cross-comparisons across the different tested DNA extraction protocols

DNA quality parameters assessed earlier were based on DNA 260/280 and 260/230 spectrophotometric ratios. Two other important quality parameters are DNA detectability by PCR and DNA integrity and size range, which can be assessed by gel electrophoresis and bioanalyzer analyses. DNA from all tested protocols exhibited detectable PCR bands of a housekeeping gene, GAPDH (Table 6 and data not shown), hence indicating that the DNA is amplifiable for specific short regions. DNA isolated by the GQ method exhibited a smear-like profile by gel electrophoresis, with peak intensity often greater than 1 kilo base pair (Kbp) (Figure 2 and Table 6). Bioanalyzer smear analyses confirmed the DNA average size peak to be greater than 1 Kbp, with an average size ranging across samples between 4.9-9.7 Kbp (Table 6 and data not shown). Compared to GQ, all tested protocols often showed more DNA degradation, except for protocol Qq which usually exhibited similar DNA smear profiles (Figure 2 and Table 6).

Figure 2
figure2

DNA integrity and size range as assessed by agarose gel electrophoresis. (A) DNA size markers used to estimate size ranges are shown in addition to genomic DNA that was isolated from white blood cells (WBC) and used as a positive control. (B) Representative DBS from each of the tested protocols are shown, except for protocol QQ in which DNA amounts were insufficient to be analyzed by gel electrophoresis. Eight different gel sections are shown and are derived from either the same gel or different gels. In each section, two punches from the same NCS spot were run on the same gel, with the first punch, labeled ‘a’, representing protocol GQ and the second punch, labeled ‘b’, representing another unique protocol from the tested set. The two blue lines, representing the 100 and 1000 base pair (bp) size ranges, were set according to the molecular size marker used in each section. The 1000 base pair limit is a minimum size range with useful applications in many genetic and epigenetic studies, including Illumina’s HM450 Beadchip array. The results of other DBS analyzed by gel electrophoresis or bioanalyzer are summarized in Table 6.

Overall comparison of the tested DNA quantity and quality parameters across the protocols shows that the best two protocols that at least match GQ in most of the tested parameters are Qq and NN (or Nn) (Table 6). Relative to GQ, the only drawback of Qq is its 260/280 DNA ratios (Table 6). Qq 260/280 ratios were always out of range, indicating that the high Qq DNA quantity measurements recorded spectrophotometrically by nanodrop may not be accurate. However, quantification with Qubit, a fluorescent-based method, confirmed the yield in Qq to be 2.7× higher than in GQ (data not shown). The only drawback of NN (or Nn) relative to GQ was DNA integrity, with the DNA size ranges in NN or Nn being lower than in GQ (Table 6).

In conclusion, protocol GQ seems to be the most robust among the tested methods across all tested DNA quantity and quality parameters. Qq can be rather suitable for applications requiring larger DNA quantities from DBS, while maintaining large fragment sizes, but in which 260/280 ratios are not a necessity. On the other hand, protocols NN or Nn may be better suitable for applications requiring larger DNA quantities from DBS, relative to GQ, while maintaining optimal 260/280 and better 260/230 ratios, but in which large DNA fragment sizes are not a requirement.

Of note, when GQ was tested on DBS (Guthrie cards, Whatman 903) from the Tasmanian Infant Health Survey (TIHS), Australia, dating more than 20 years old, an average of 66 ± 15 ng (n = 3) of DNA could be extracted per two punches, each being 1 mm in diameter, with a mean 260/280 ratio of 1.66 ± 0.02; these DNA quantities are equivalent to 42.0 ng/mm2 for TIHS, compared to 12.7 ng/mm2 for NCS samples.

Performance of DNA extracted from DBS using methylome-wide analysis

Methylation probe call index

In order to assess the performance of DNA extracted from DBS in HM450 methylome-wide analyses, we used the DNA extracted by GQ, being the most robust protocol. DNA from two sample pairs, each representing two punches (serving as technical replicates) from the same DBS, were analyzed by HM450. In addition, DBS pairs were compared to reference DNA from neonatal blood or cell lines (Table 7). In all tested samples, whether originating from neonatal blood, DBS or cell lines, more than 99% of the 485,577 HM450 individual probes were detected, using the commonly accepted quality control detection p-value of 0.01, hence, indicating high quality data. The average beta-values were similar between the technical replicates NCS 37a and 37b (approximately 0.47 for either) and between the technical replicates NCS 38a and 38b (approximately 0.42 for either) (Table 7).

Table 7 Methylation quality control probe evaluation

Sample–dependent and –independent HM450 internal quality control probes

For sample and array quality, HM450 array includes 850 quality control (QC) probes. Fifteen QC probes are sample-independent and 835 QC probes are sample-dependent [12]. DNA from DBS, neonatal blood or cell lines passed the described HM450 QC (Additional file 3 and Figure 3). Background probes, wherever included, produced minimal signal (maximum limit is 1000 units, as recommended by Illumina Inc.), and the intended positive signals from the experimental QC probes were above background, for all of the three tested DNA sources (Additional file 3 and exemplified in Figure 3 using non-polymorphic probes, which are indicative of overall performance). In addition, performance of DBS samples was similar to that of subsets taken from reference neonatal blood and cell line samples (Additional file 4 and Figure 3). Bisulfite conversion efficiency for both, type I and II probes, was high for all tested samples (Additional file 4) and was confirmed by PCR using primers that are specific either to bisulfite converted or to non-modified GAPDH DNA regions (data not shown).

Figure 3
figure3

HM450 QC plot using Non-polymorphic probes which assess overall performance. In the green channel, background signals are shown in red and pink while positive signals in opaque and fluorescent green. In the red channel, background signals are shown in opaque and fluorescent green while positive signals in red and pink. One non-polymorphic control has been designed for each of the four nucleotides A, T, C, and G. Four DBS DNA samples are shown between four neonatal blood and four cell line DNA samples, in each of the two plots. The DBS samples represent two NCS spots, 37 and 38, each consisting of two tested punches labeled ‘a’ or ‘b’.

Differential methylation and clustering analyses using HM450 data

Differential methylation of HM450 beta-values produced two major clusters (Figure 4). Cell line DNA samples formed one cluster while DBS and neonatal blood DNA formed another; this is expected because the cell lines used are of hepatic tissue origin while both, DBS and neonatal blood samples, are of blood tissue origin. Within the cluster of blood biospecimens, all four neonatal blood samples formed one sub-cluster, which was segregated away from the DBS sub-cluster. The two punches NCS 37a and 37b, representing the technical replicates from the same spot NCS 37, clustered together and away from the other two technical replicates, NCS 38a and 38b, which also clustered together (Figure 4). This further supports the observed higher correlations (p < 0.001, Steiger Z test) between technical duplicates punched from the same DBS (r2 range = 0.963-0.990) versus different DBS (r2 range = 0.949-0.951; Additional file 4A). When the analysis was limited to the top 1% of probes that showed the highest variance in M-values (transformed betas) across any of the four tested DBS, NCS_37a, NCS_37b, NCS_38a and NCS_38b, the correlation between replicates (r2 range = 0.850-0.972; p < 0.001, Pearson) become significantly higher and more predictive of replication (p < 0.001, Steiger Z test) than the correlation between non-replicates (r2 range = 0.365-0.371; p < 0.001, Pearson) (Additional file 4B). Similar observations were reported using Spearman correlations. Moreover, the frequency distributions of delta M-values (δM) between samples were centred at zero only between technical replicates (Additional file 4B). Hence, we can conclude that HM450 analyses using DBS DNA, extracted using the GQ protocol, is reproducible.

Figure 4
figure4

Differential methylation and unsupervised clustering analysis of HM450 data from neonatal blood, DBS and cell line DNA. Neonatal blood and cell line DNA samples are used as positive controls of good DNA quality for reference comparisons with DNA extracted from DBS. Neonatal blood and DBS are from different individuals. Four DBS DNA samples are shown between four different neonatal blood and four different cell line DNA samples. The DBS samples represent two NCS spots, 37 and 38, each consisting of two tested punches labeled ‘a’ or ‘b’. HM450 beta-values were clustered using Euclidean distance as the dissimilarity index. As shown in the color key, the red and blue signals represent relatively hypomethylated and hypermethylated regions, respectively.

Performance of DNA extracted from DBS using sequence-specific methylation analysis

The performance of DNA extracted from DBS by the GQ method was then tested using sequence-specific methylation analyses. For this purpose, we analyzed the methylation levels of several CpG sites in Line1 and AluYb8 sequences, both of which are proxy markers of global methylation, being transposable elements interspersed across the genome [16]. The technical replicates, NCS 37a and 37b, showed similar Line1 and AluYb8 methylation levels at each tested CpG site, and similar data were observed with the pair, NCS 38a and 38b (Figure 5A and B). These results show that inter-replicate variation in methylation levels is minimal using several CpG sites in two different loci, Line1 and AluYb8. Moreover, the observed difference in methylation levels between the pair NCS 37a and 37b versus NCS 38a and 38b was consistent at every single CpG tested and across both, Line1 and AluYb8 loci (p < 0.1 for CpG6 in Line1 and CpG3 in AluYb8 and p < 0.05 for all other CpGs; Mann–Whitney test) (Figure 5A and B). These results confirm that DNA extracted from DBS using GQ is suitable to detect small methylation differences in a consistent manner and with low inter-replicate variation.

Figure 5
figure5

Methylation analyses of Line1 and AluYb8 loci using bisulfite pyrosequencing. The methylation levels of six and four CpG sites were analyzed for Line1 (A) and AluYb8 (B), respectively, and are expressed as percent of the total number of CpGs analyzed for each individual CpG site. The DBS samples represent two NCS spots, 37 and 38, each consisting of two tested punches labeled ‘a’ or ‘b’.

Discussion

DBS have become an increasingly important tool for diagnostic purposes and for epigenetic, genetic and epidemiological research. We have previously tested a range of commercially available DNA extraction kits for purifying genomic DNA from fresh and dried blood for downstream PCR and DNA methylation applications [5]. We found that genomic DNA extraction, using the ChargeSwitch Forensic DNA Purification kit (Invitrogen), with subsequent bisulfite modification, using the MethylEasy kit (Human Genetic Signatures), was best in yielding bisulfite-converted DNA of sufficient quantity and quality for downstream candidate-gene DNA methylation analyses, such as SEQUENOM MassArray EpiTYPER analysis [5]. However, DNA extraction with ChargeSwitch was recently shown to be not suitable for beadchip methylation profiling, leading to up to 16% loss of detectable probes in Infinum HumanMethylation27 (Illumina Inc.) arrays analysis [9]. In relation to the limited amounts of DNA extracted from DBS, a recent report also pointed to the biases introduced by whole bisulfitome amplification and the need for careful data interpretation [13], as we have also observed in this study.

In this work, we have systematically compared different DNA extraction methods from DBS by dissecting different phases of extraction and optimizing several steps within each phase, using commercial and in-house extraction protocols. For these purposes, we used a homogenous set of DBS samples, spotted on the same day and stored in a similar manner, to provide a common platform for cross-protocol comparisons. Moreover, this study emphasizes DNA extraction protocols that have particular utility in a recent technology for studying methylome-wide methylation, Infinium HM450, and sequence-specific methylation, by pyrosequencing. The use of DBS for diagnostic and research purposes is not new, but there is a lack of quality standards for optimizing DNA extraction. This study suggests different DNA extraction protocols, each having specific advantages tailored for specific applications. Protocol GQ does not extract the highest DNA yield, but provides DNA in quantities and qualities sufficient for HM450 methylome-wide and sequence-specific methylation analyses. With GQ, the 260/280 ratios are consistently optimal and the extracted DNA is less fragmented relative to other protocols. Protocol Qq, on the other hand, produces twice as much DNA as GQ and with similar DNA integrity. However, 260/280 ratios in Qq are unreliable and cannot be used for sample selection, particularly for expensive downstream applications. As for protocols NN or Nn, they extract 1.7 folds more DNA than GQ and show optimal ranges for both, 260/280 and 260/230 ratios. However, these protocols lead to more DNA fragmentation relative to GQ. This may represent a limitation for bead array-based assays and other DNA methylation assays where DNA integrity is a requirement. It should be noted though that 260/280 and 260/230 ratios should be treated with caution; for example, different contaminants can compensate for each other’s’ deviations, resulting in misleading optimal 260/280 ratios.

Because the type of purification column was identical in GQ and Qq, but different from those used in all the other methods, column nature (Phase II) could be the reason for the better DNA integrity in the two protocols. This is further supported by our results showing that in GQ, changing Phase II (includes column type) while maintaining Phase I, as in GN/n or GN/n-XS, compromises DNA integrity (Table 6). Changing the DNA precipitation buffer to ethanol reduces the need in some protocols to vigorously vortex to dissolve resultant precipitates but does not seem to enhance DNA integrity (Figure 2; compare GN versus Gn, GN-XS versus Gn-XS, and NN versus Nn). However, ethanol, as a precipitating buffer, was essential in some protocols to increase the DNA yield, as in Qq versus QQ (Tables 1 and 6).

Other studies have compared DNA extraction protocols from DBS, but irrespective of epigenetic applications [1719]. It is difficult to compare DNA extraction protocols across different studies due to many reasons, such as differences in the structures of filter papers on which blood was soaked, the storage conditions and year-durations, and DNA quantifications methods used. However, these studies used QIAamp DNA Mini Kit (Qiagen) as a reference method, which is similar to the QQ protocol in this work, hence, allowing comparisons between our and their optimized methods. Sjoholm et al. reported that QIAamp DNA Mini Kit performed the best relative to EZNA (Omega Bio-Tek), Chelex 100 (AmershamBiosciences) and alkaline lysis (GenomiPhi DNA Amplification Kit, AmershamBiosciences) [19]. On the other hand, the in-house developed DNA extraction methods reported by Hue et al. and Hollegaard et al. yielded, by Nanodrop quantification, 3.3 and 2.5 fold more DNA, respectively than QIAamp DNA Mini Kit [18]. In comparison, the GQ protocol in our study yielded, also by Nanodrop quantification, on average 4.6 fold more DNA than matched QQ samples; in addition, protocols Qq, NN and Nn yielded at least 1.7 folds more DNA than matched GQ samples. Moreover, the in-house protocol by Hue et al. produced a low purity 260/280 average ratio (1.50) [18], while GQ, NN and Nn ratios were optimal in every tested sample. These findings support the good performance of our optimized methods relative to many other in-house and commercial DNA extraction protocols from DBS. Interestingly, one study reported a recent method suitable for performing scalable DNA extractions simultaneously from many DBS, but with less emphasis on DNA quality and yield comparisons across different methods [20]. The scale of our tested methods can be increased by implementing the QIAcube technology (Qiagen), and, with scalable designs, laser cutting of DBS punches would eliminate cross-contamination, as has been recently reported [21].

Conclusion

This study arises from an international effort across several cohorts and working groups aiming to fulfill the need to systematize quality standards in DNA extraction and to increase the DNA yield using DBS, particularly with the advent of high-throughput epigenomic technology that require high quality and quantity of DNA. Given the emerging appreciation of DBS collected at birth as a valuable resource for epigenetic analyses prior to phenotypic onset, our optimized methods for DNA extraction with application in methylation analyses have great potential for diagnostic and research purposes.

Methods

Sample overview

DBS that were used to perform the comparisons across the DNA extraction protocols were obtained from NCS, USA, and have been spotted on Flinders Technology Associates (FTA) mini cards on the same day and dried in air-sealed containers for approximately two years at room temperature. Guthrie cards (Whatman 903) from TIHS, Australia, dating more than 20 years old, were also used to test the efficiency of the robust protocol GQ. Both NCS and TIHS samples were heel-prick without anticoagulants added. Permissions from the ethical committees of the International Agency for Research on Cancer (IARC), as well as both, NCS and TIHS, were obtained. TIHS is one of the founder cohorts of the International Childhood Cancer Cohort Consortium (I4C) [1].

DNA extraction protocols

Combinations of different commercially available DNA extraction kits were used, including QIAamp DNA Micro Kit (Qiagen 56304), GenSolve (Gen Vault, GVR110), NucleoSpin (gDNA clean-up, Macherey-Nagel 740230), and the extra-small (XS) version of NucleoSpin (gDNA clean-up XS, Macherey-Nagel 740904). Reported quantifications were done using Nanodrop, unless indicated otherwise using Qubit™ dsDNA High Sensitivity Assay, Invitrogen Q32851. Detailed protocols are included in Additional Methods (Additional File 2).

Gel electrophoresis and bioanalyzer analysis

Samples were run on a 0.8% agarose gel (Eurobio GEPAGA07-65) in 1 × Tris Acetate-EDTA buffer and stained with GelRed. 300 ng of DBS DNA were utilized per sample in electrophoresis analyses of DNA integrity and size range. The following DNA size markers were used: 80–10,000 bp ladder (Thermo SM0403), 500–5,000 bp ladder (Takara 3411A), and 100–1,000 bp ladder (Thermo SM0243). As for bioanalyzer analysis, 500 pg of DNA per sample was loaded on the chip and analyzed on Agilent 2100 Bioanalyzer, as per manufacturer’s instructions (Agilent Technology, High Sensitivity DNA Kit, 5067–4626).

Bisulfite conversion and PCR

DBS DNA (300 ng) samples were bisulfite converted using EZ DNA Methylation Kit (Zymo Research D5001) according to manufacturer’s instructions. To assess the efficiency of bisulfite conversion, DNA was amplified using PCR primers that were specific either to bisulfite-converted or to non-modified DNA, and spanning the region of the housekeeping gene, GAPDH. The primer pairs spanning bisulfite-converted GAPDH regions are termed GAPDH-bc and consisted of the following forward and reverse primers, respectively: 5’-GTATTTGTTGATGGGTTAAGG-3’ and 5’-ATAAAAACAAATCCCCTACCC-3’. The primer pairs spanning non-modified GAPDH regions are termed GAPDH-nm and consisted of the following forward and reverse primers, respectively: 5’-CTCTTGCTACTCTGCTCTGG-3’ and 5’-GCTAAGTTTAGCCTGCCTGG-3’. Efficient conversion is observed when PCR bands are detected with GAPDH-bc but not with GAPDH-nm for a given sample. The PCR conditions used were: 95°C 15 min, [95°C 30 s, 57°C 30 s, 72°C 30 s] × 50 cycles, 72°C 10 min, and pause 4°C.

Pyrosequencing

DNA methylation analysis by pyrosequencing of bisulfie-converted DNA was performed as described [22]. Briefly, the region of interest was amplified with forward and reverse primers, one of which is biotinylated (btn), and then the methylation levels of the amplified region were analyzed using a sequencing primer. The forward, reverse and sequencing primers used for Line1 were the following, respectively, and adopted from Daskalos et al. [16]: 5’-btn-TAGGGAGTGTTAGATAGTGG-3’, 5’-AACTCCCTAACCCCTTAC-3’ and 5’-CAAATAAAACAATACCTC-3’. The forward, reverse and sequencing primers used for AluYb8 were the following, respectively: 5’-AGATTATTTTGGTTAATAAG-3’, 5’-btn-AACTACRAACTACAATAAC-3’ and 5’-GTTTGTAGTTTTAGTTATT- 3’, as previously described [23].

Illumina Infinium HM450 array and data processing

Infinium HM450 arrays were processed according to manufacturer’s instructions. GenomeStudio was used to analyze quality controls. The raw colour channels were corrected using the internal control probes and converted, without background subtraction and normalization, into absolute methylation levels (beta-values). Data were then imported into R (3.0.0), using the minfi package version 1.2.0 (http://www.bioconductor.org). Subset-quantile Within Array Normalisation (SWAN) normalisation was performed to correct for technical discrepancies between Type I and Type II [24]. Probes with detection p-values above 0.01 were considered as background noise and omitted from further analysis. Sex chromosome-specific probes were eliminated to minimize gender-specific variation of the X versus Y chromosomes [25]. Logarithmic transformation of the beta-values into M values was done, as previously described [26]. Statistical tests were performed using M-values. Clustering plots were generated using the lumi R package and based on coefficient of variation, which is calculated by standard deviation divided by mean across samples [27]. Differentially methylated CpGs were identified using an F-test in minfi.

Abbreviations

bp:

base pair

DBS:

dried blood spots

gDNA:

genomic DNA

GQ:

GenSolve-QIAamp

GN:

GenSolve-NucleoSpin

Gn:

GN with NucleoSpin’s precipitation buffer replaced by ethanol

GN-XS:

GN extra small

Gn-XS:

GN-XS with NucleoSpin’s precipitation buffer replaced by ethanol

HM450:

Infinium HumanMethylation450 Beadchip

Kbp:

kilo base pair

NCS:

National Children Study - USA

NN:

NucleoSpin-NucleoSpin

Nn:

NN with NucleoSpin’s precipitation buffer replaced by ethanol

NN-XS:

NN extra small

QC:

quality control

QQ:

QIAamp-QIAamp

Qq:

QQ with QIAamp’s precipitation buffer replaced by ethanol

SWAN:

Subset-quantile Within Array Normalisation

TIHS:

Tasmanian Infant Health Survey - Australia.

References

  1. 1.

    Brown RC, Dwyer T, Kasten C, Krotoski D, Li Z, Linet MS, Olsen J, Scheidt P, Winn DM: Cohort profile: the International Childhood Cancer Cohort Consortium (I4C). Int J Epidemiol. 2007, 36: 724-730.

  2. 2.

    Guthrie R, Susi A: A simple phenylalanine method for detecting phenylketonuria in large populations of newborn infants. Pediatrics. 1963, 32: 338-343.

  3. 3.

    Mei JV, Alexander JR, Adam BW, Hannon WH: Use of filter paper for the collection and analysis of human whole blood specimens. J Nutr. 2001, 131: 1631S-1636S.

  4. 4.

    He H, Argiro L, Dessein H, Chevillard C: Improved technique that allows the performance of large-scale SNP genotyping on DNA immobilized by FTA technology. Infect Genet Evol. 2007, 7: 128-132.

  5. 5.

    Wong N, Morley R, Saffery R, Craig J: Archived Guthrie blood spots as a novel source for quantitative DNA methylation analysis. Biotechniques. 2008, 45: 423-424. 426, 428 passim

  6. 6.

    Cruickshank MN, Pitt J, Craig JM: Going back to the future with Guthrie-powered epigenome-wide association studies. Genome Med. 2012, 4: 83-

  7. 7.

    Beyan H, Down TA, Ramagopalan SV, Uvebrant K, Nilsson A, Holland ML, Gemma C, Giovannoni G, Boehm BO, Ebers GC, Lernmark Å, Cilio CM, Leslie RD, Rakyan VK: Guthrie card methylomics identifies temporally stable epialleles that are present at birth in humans. Genome Res. 2012, 22: 2138-2145.

  8. 8.

    Aberg KA, Xie LY, Nerella S, Copeland WE, Costello EJ, van den Oord EJCG: High quality methylome-wide investigations through next-generation sequencing of DNA from a single archived dry blood spot. Epigenetics. 2013, 8: 542-547.

  9. 9.

    Hollegaard MV, Grauholm J, Nørgaard-Pedersen B, Hougaard DM: DNA methylome profiling using neonatal dried blood spot samples: a proof-of-principle study. Mol Genet Metab. 2013, 108: 225-231.

  10. 10.

    Joo JE, Wong EM, Baglietto L, Jung C-H, Tsimiklis H, Park DJ, Wong NC, English DR, Hopper JL, Severi G, Giles GG, Southey MC: The use of DNA from archival dried blood spots with the Infinium HumanMethylation450 array. BMC Biotechnol. 2013, 13: 23-

  11. 11.

    Cruickshank MN, Oshlack A, Theda C, Davis PG, Martino D, Sheehan P, Dai Y, Saffery R, Doyle LW, Craig JM: Analysis of epigenetic changes in survivors of preterm birth reveals the effect of gestational age and evidence for a long term legacy. Genome Med. 2013, 5: 96-

  12. 12.

    Bibikova M, Barnes B, Tsan C, Ho V, Klotzle B, Le JM, Delano D, Zhang L, Schroth GP, Gunderson KL, Fan J-B, Shen R: High density DNA methylation array with single CpG site resolution. Genomics. 2011, 98: 288-295.

  13. 13.

    Bundo M, Sunaga F, Ueda J, Kasai K, Kato T, Iwamoto K: A systematic evaluation of whole genome amplification of bisulfite-modified DNA. Clin Epigenetics. 2012, 4: 22-

  14. 14.

    Hardin J, Finnell RH, Wong D, Hogan ME, Horovitz J, Shu J, Shaw GM: Whole genome microarray analysis, from neonatal blood cards. BMC Genet. 2009, 10: 38-

  15. 15.

    Bank S, Nexø BA, Andersen V, Vogel U, Andersen PS: High-quality and -quantity DNA extraction from frozen archival blood clots for genotyping of single-nucleotide polymorphisms. Genet Test Mol Biomarkers. 2013, 17: 501-503.

  16. 16.

    Daskalos A, Nikolaidis G, Xinarianos G, Savvari P, Cassidy A, Zakopoulou R, Kotsinas A, Gorgoulis V, Field JK, Liloglou T: Hypomethylation of retrotransposable elements correlates with genomic instability in non-small cell lung cancer. Int J cancer. 2009, 124: 81-87.

  17. 17.

    Hollegaard MV, Thorsen P, Norgaard-Pedersen B, Hougaard DM: Genotyping whole-genome-amplified DNA from 3- to 25-year-old neonatal dried blood spot samples with reference to fresh genomic DNA. Electrophoresis. 2009, 30: 2532-2535.

  18. 18.

    Hue NT, Chan NDH, Phong PT, Linh NTT, Giang ND: Extraction of human genomic DNA from dried blood spots and hair roots. Int J Biosci Biochem Bioinforma. 2012, 2: 21-26.

  19. 19.

    Sjöholm MIL, Dillner J, Carlson J: Assessing quality and functionality of DNA from fresh and archival dried blood spots and recommendations for quality control guidelines. Clin Chem. 2007, 53: 1401-1407.

  20. 20.

    Saavedra-Matiz CA, Isabelle JT, Biski CK, Duva SJ, Sweeney ML, Parker AL, Young AJ, Diantonio LL, Krein LM, Nichols MJ, Caggana M: Cost-effective and scalable DNA extraction method from dried blood spots. Clin Chem. 2013, 59: 1045-1051.

  21. 21.

    Murphy SC, Daza G, Chang M, Coombs R: Laser cutting eliminates nucleic acid cross-contamination in dried-blood-spot processing. J Clin Microbiol. 2012, 50: 4128-4130.

  22. 22.

    Paliwal A, Vaissière T, Krais A, Cuenin C, Cros M-P, Zaridze D, Moukeria A, Boffetta P, Hainaut P, Brennan P, Herceg Z: Aberrant DNA methylation links cancer susceptibility locus 15q25.1 to apoptotic regulation and lung cancer. Cancer Res. 2010, 70: 2779-2788.

  23. 23.

    Choi SH, Worswick S, Byun H-M, Shear T, Soussa JC, Wolff EM, Douer D, Garcia-Manero G, Liang G, Yang AS: Changes in DNA methylation of tandem DNA repeats are different from interspersed repeats in cancer. Int J cancer J. 2009, 125: 723-729.

  24. 24.

    Maksimovic J, Gordon L, Oshlack A: SWAN: Subset-quantile within array normalization for illumina infinium HumanMethylation450 BeadChips. Genome Biol. 2012, 13: R44-

  25. 25.

    Morris T, Lowe R: Report on the Infinium 450 k methylation array analysis workshop: April 20, 2012 UCL, London, UK. Epigenetics. 2012, 7: 961-962.

  26. 26.

    Du P, Zhang X, Huang C-C, Jafari N, Kibbe WA, Hou L, Lin SM: Comparison of Beta-value and M-value methods for quantifying methylation levels by microarray analysis. BMC Bioinformatics. 2010, 11: 587-

  27. 27.

    Du P, Kibbe WA, Lin SM: lumi: a pipeline for processing Illumina microarray. Bioinformatics. 2008, 24: 1547-1548.

Download references

Acknowledgements

We thank Dr. Florence Le Calvez-Kelm and Mr. Geoffroy Durand for helping in the HM450 experimental process and bioanalyzer measurements. This work was partially supported by the Association pour le Recherche Contre le Cancer (l’ARC). AG is supported by the IARC Postdoctoral Fellowship and Marie Curie Actions-People-COFUND. We thank the International Childhood Cancer Cohort Consortium for contributing to the coordination between participating cohorts.

The findings and conclusions in this report are those of the authors and do not necessarily represent the views of the National Children’s Study, National Institutes of Health, or U.S. Department of Health and Human Services. Supported in part by NICHD Contracts HHSN275200503414C, HHSN275200503411C, HHSN275200603416C, HHSN275200503415C, HHSN275200503413C, HHSN275200503410C, HHSN275200503396C, HHSN275201000121U, and HHSN275200900010C.

This manuscript is dedicated to the participating children and their families.

Author information

Correspondence to Hector Hernandez-Vargas.

Additional information

Competing interests

No conflicts of interest are disclosed.

Authors’ contributions

AG carried out all experiments and drafted the manuscript. MPC participated in the optimization of the protocols. ALP participated in the collection of TIHS samples. SH and KC participated in the collection of NCS samples. ZH, TD, RS and HH conceived the study and participated in its design and coordination and helped to draft the manuscript. AG and HH performed the statistical and bioinformatics analyses. All authors read and approved the final manuscript.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Ghantous, A., Saffery, R., Cros, M. et al. Optimized DNA extraction from neonatal dried blood spots: application in methylome profiling. BMC Biotechnol 14, 60 (2014) doi:10.1186/1472-6750-14-60

Download citation

Keywords

  • Blood spot
  • DNA extraction
  • Epigenetics
  • Methylome
  • HM450
  • Pyrosequencing
  • Whole bisulfitome amplification
  • QIAamp
  • GenSolve
  • NucleoSpin