Degraded DNA Samples Research Articles

Development of a forensic DNA phenotyping panel using massive parallel sequencing

The HIrisPlex-S system, targeting a total of 41 SNPs, allows the simultaneous eye, hair and skin color prediction from DNA. In the present study, we developed a massive parallel sequencing (MPS) multiplex assay in order to genotype all the HIrisPlex-S markers in degraded casework samples. PCR amplicons sizes of target regions were kept below 180 bp, in order to allow analysis of degraded DNA samples. Individuals with known phenotype, artificially degraded DNA samples and a set of 2800M control DNA dilutions were sequenced on a Ion PGM System, in order to evaluate the concordance testing results and the forensic suitability of this 41-plex MPS assay. Full and reliable profiles could be obtained with 0.1 ng of input DNA. The increment of the number of PCR cycles results in improvement of sensitivity or in typing results but an increase of artifacts were also observed.

Targeted capture and sequencing of 1245 SNPs for forensic applications

Next generation sequencing (NGS) technologies have enabled the possibility of analyzing a large number of SNPs simultaneously from multiple samples in a single experiment, for complementing the shortcomings of STR based methods. To efficiently genotype the desired SNPs, it is critical to optimize the library construction procedures. In this study, we formulated a strategy combining the molecular inversion probe (MIP) based target region capture method and NGS for genotyping 1245 SNPs. All the SNPs we selected exhibited high heterozygosity (minor allele frequency (MAF) > 0.3) according to 1000 genomes data. We applied the method to genotype a population of 210 unrelated individuals from the Hubei province of China and assessed the allele frequencies, Hardy-Weinberg equilibrium and linkage disequilibrium. The MAFs of more than 95% of the SNPs were ≥0.2, and no significant deviation or strong linkage was observed for 98% of the SNPs. The data indicated that, even within a relatively confined region, our SNP panel is suitable for individual identifications. Furthermore, we performed paternity test for 7 trio families using low quality DNA samples. The conclusions are in total agreement with these of STR-based analyses, with higher confidence indexes. Finally, we evaluated the performance of the MIP-NGS method with mock degraded DNA samples. We were able to genotype most of the SNPs even when the genomic DNA was sonicated to ˜100 bp range. In summary, we established a highly accurate and cost-effective method of SNP genotyping, which is potentially capable of solving complex issues encountered in forensic practices.

Genotyping genetic markers from LCN and degraded DNA by HRM and their application in hair shaft

Degraded and low copy number (LCN) DNA samples are common challenging materials in forensic casework because they increase the difficulty of sample processing and reduce the possibility of obtaining genetic information from DNA. High-resolution melting (HRM) curve analysis is promising for genotyping genetic markers and has been applied to the detection of LCN and degraded DNA in the field of forensic science. However, the exact assessment based on HRM at multiple genetic markers from both degraded and LCN DNA materials has not been optimized. To explore the ability of HRM to genotype LCN and degraded DNA samples, we selected three genetic markers to genotype in experimental LCN and degraded DNA and practical hair shaft materials, which are often encountered as degraded and LCN DNA in forensic medicine. The results show that DNA samples of as low as 100pg and as short as 60bp were successfully genotyped by the HRM assay at all three genetic markers, whereas in hair shaft DNA, two loci were accurately genotyped. The HRM assay established in this study can be applied to LCN and degraded DNA analysis in forensic casework and can act as a reference point before genotyping short tandem repeat markers. Developing the HRM strategy for genotyping DNA genetic markers enriches detectable targets in hair shaft samples and provides valuable data for further exploration.

A single nucleotide polymorphism panel for individual identification and ancestry assignment in Caucasians and four East and Southeast Asian populations using a machine learning classifier.

Single nucleotide polymorphism (SNP) profiling is an effective means of individual identification and ancestry inferences in forensic genetics. This study established a SNP panel for the simultaneous individual identification and ancestry assignment of Caucasian and four East and Southeast Asian populations. We analyzed 220 SNPs (125 autosomal, 17 X-chromosomal, 30 Y-chromosomal, and 48 mitochondrial SNPs) of the DNA samples from 563 unrelated individuals of five populations (89 Caucasian, 234 Taiwanese Han, 90 Filipino, 79 Indonesian and 71 Vietnamese) and 18 degraded DNA samples. Informativeness for assignment (In) was used to select ancestry informative SNPs (AISNPs). A machine learning classifier, support vector machine (SVM), was used for ancestry assignment. Of the 220 SNPs, 62 were individual identification SNPs (IISNPs) (51 autosomal and 11 X-chromosomal SNPs) and 191 were AISNPs (100 autosomal, 13 X-chromosomal, 30 Y-chromosomal, and 48 mitochondrial SNPs). The 51 autosomal IISNPs offered cumulative random match probabilities (cRMPs) ranging from 1.56 × 10-21 to 3.16 × 10-22 among these five populations. Using AISNPs with the SVM, the overall accuracy rate of ancestry inference achieved in the testing dataset between Caucasian, Taiwanese Han, and Filipino populations was 88.9%, whereas it was 70.0% between Caucasians and each of the four East and Southeast Asian populations. For the 18 degraded DNA samples with incomplete profiling, the accuracy rate of ancestry assignment was 94.4%. We have developed a 220-SNP panel for simultaneous individual identification and ethnic origin differentiation between Caucasian and the four East and Southeast Asian populations. This SNP panel may assist with DNA analysis of forensic casework.

Development of the MitoQ assay as a real-time quantification of mitochondrial DNA in degraded samples.

Mitochondrial DNA is a reliable genetic material for estimating maternally related haplogroups and ancestries. Exploring maternal DNA inheritance is particularly useful when nuclear DNA is degraded or limited, as the copy number of mitochondrial DNA is far greater than the copy number of nuclear DNA. Normal mitochondrial DNA copy number has been estimated to 100 copies per buccal epithelial cell, 4000 copies in skeletal cells and 7000 copies in myocardial cells. This estimation is usually performed via extrapolation from the nuclear DNA quantitation. It is essential to reduce this variability and accurately quantify the exact number of copies of mitochondrial DNA, especially in compromised samples of a forensic or ancient nature. While useful, the testing of mitochondrial DNA is often long and costly and comes with limited success. The accurate quantification of mitochondrial DNA using specific quantitative PCR assays can be used to make better decisions on the downstream testing and success of amplification. As a result, this study develops a real-time assay for the quantification of mitochondrial DNA copy number and assesses its performance on a set of degraded DNA samples. The developed MitoQ assay has been shown to be highly specific to the human mitochondrial genome with no amplification of nuclear pseudogenes being observed and outperformed a previously published concordant assay. Additionally, a high sensitivity was measured to 280 copies of mitochondrial DNA. Minimal variation was observed between each replication cycle, indicating the assay to be robust and repeatable. Overall, this study presents a real-time assay that is sensitive and robust to quantifying mitochondrial DNA copy number in degraded samples. Furthermore, there is potential to incorporate the assay as an additional target in current qPCR assays which use a six-dye chemistry and provide a complete overview of a sample's quality and quantity.

Molecular markers used in forensic genetics.

Forensic genetics is a field that has become subject to increasing interest in recent years. Both the technology and the markers used for forensic purposes have changed since the 1980s. The minisatellite sequences used in the famous Pitchfork case introduced genetics to the forensic sciences. Minisatellite sequences have now been replaced by more sensitive microsatellite markers, which have become the basis for the creation of genetic profile databases. Modern molecular methods also exploit single nucleotide polymorphisms, which are often the only way to identify degraded DNA samples. The same type of variation is taken into consideration in attempting to establish the ethnicity of a perpetrator and to determine phenotypic traits such as the eye or hair colour of the individual who is the source of the genetic material. This paper contains a review of the techniques and molecular markers used in human and animal forensic genetics, and also presents the potential trends in forensic genetics such as phenotyping.

Abstract 430: Detection of low-frequency variants from highly degraded DNA samples using a novel library preparation method

Abstract Diagnostic tools based on next generation sequencing are fundamentally transforming clinical oncology. However, there is a lack of adequate library preparation strategies for highly degraded, clinically relevant samples, such as cell-free DNA (cfDNA) and FFPE DNA. Furthermore, clinical samples are often available in limited quantity, making the detection of pathologic variants particularly challenging. Due to the extreme heterogeneity of these sample types, targeted sequencing is often used to achieve deep coverage of genomic loci and enable detection of low-frequency variants. Commercially available protocols for library preparation require stringent size-selection to remove adapter-dimers, which reduces library complexity. Achieving high specificity can be challenging because low-frequency artifacts arise from a variety of sources, including DNA extraction, library construction, PCR, hybrid selection, and sequencing. These artifacts can be identified by “duplex sequencing”, where strand-specific unique molecular identifiers (UMIs) are used to confirm the presence of an alteration on both strands of an input molecule. However, duplex sequencing typically delivers low conversion rates with degraded samples due to poor ligation efficiency and template loss during size-selection. We present a novel library construction chemistry that eliminates adapter dimers and utilizes duplexed UMIs to increase sensitivity. The method employs a two-step ligation procedure that does not require size-selection or adjustment of AMPure bead ratios during cleanup. We achieve maximal library conversion using a unique, mutant DNA ligase and proprietary sequencing adapters that increase ligation efficiency and suppress chimera formation. We demonstrate performance using three sample types: sheared genomic DNA, cfDNA, and FFPE DNA, across a wide input range (1-1000 ng). To mimic performance with degraded DNA, we created libraries from mixtures of NA12878 and NA24385 cell line DNA sheared to 150 bp with inputs of 1-25 ng and mutant allele fractions (MAFs) down to 0.25%. Libraries were enriched using a 75 kb panel of xGen Lockdown Probes, followed by ultra-deep sequencing and variant calling. When compared to commercially available methods, our approach yielded a 50-100% increase in library complexity with significantly improved sensitivity to <1% variants. We also obtained 100% specificity using duplexed UMI correction, which removed all false-positive calls. To highlight clinical utility, we extended our study to cfDNA samples with inputs of 5-25 ng and MAFs down to 1%, and FFPE DNA samples with inputs of 25-100 ng and MAFs down to 1%. Citation Format: Ushati Das Chakravarty, Zac Zwirko, Yu Zheng, Madelyn Light, Kevin Lai, Keith Bryan, Scott Rose, Yun Bao, Mirna Jarosz, Caifu Chen. Detection of low-frequency variants from highly degraded DNA samples using a novel library preparation method [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 430.

Systematic evaluation of the early access applied biosystems precision ID Globalfiler mixture ID and Globalfiler NGS STR panels for the ion S5 system

Systematic evaluation of the early access applied biosystems precision ID Globalfiler mixture ID and Globalfiler NGS STR panels for the ion S5 system

The HIrisPlex-S system for eye, hair and skin colour prediction from DNA: Introduction and forensic developmental validation

Forensic DNA Phenotyping (FDP), i.e. the prediction of human externally visible traits from DNA, has become a fast growing subfield within forensic genetics due to the intelligence information it can provide from DNA traces. FDP outcomes can help focus police investigations in search of unknown perpetrators, who are generally unidentifiable with standard DNA profiling. Therefore, we previously developed and forensically validated the IrisPlex DNA test system for eye colour prediction and the HIrisPlex system for combined eye and hair colour prediction from DNA traces. Here we introduce and forensically validate the HIrisPlex-S DNA test system (S for skin) for the simultaneous prediction of eye, hair, and skin colour from trace DNA. This FDP system consists of two SNaPshot-based multiplex assays targeting a total of 41 SNPs via a novel multiplex assay for 17 skin colour predictive SNPs and the previous HIrisPlex assay for 24 eye and hair colour predictive SNPs, 19 of which also contribute to skin colour prediction. The HIrisPlex-S system further comprises three statistical prediction models, the previously developed IrisPlex model for eye colour prediction based on 6 SNPs, the previous HIrisPlex model for hair colour prediction based on 22 SNPs, and the recently introduced HIrisPlex-S model for skin colour prediction based on 36 SNPs. In the forensic developmental validation testing, the novel 17-plex assay performed in full agreement with the Scientific Working Group on DNA Analysis Methods (SWGDAM) guidelines, as previously shown for the 24-plex assay. Sensitivity testing of the 17-plex assay revealed complete SNP profiles from as little as 63 pg of input DNA, equalling the previously demonstrated sensitivity threshold of the 24-plex HIrisPlex assay. Testing of simulated forensic casework samples such as blood, semen, saliva stains, of inhibited DNA samples, of low quantity touch (trace) DNA samples, and of artificially degraded DNA samples as well as concordance testing, demonstrated the robustness, efficiency, and forensic suitability of the new 17-plex assay, as previously shown for the 24-plex assay. Finally, we provide an update to the publically available HIrisPlex website https://hirisplex.erasmusmc.nl/, now allowing the estimation of individual probabilities for 3 eye, 4 hair, and 5 skin colour categories from HIrisPlex-S input genotypes. The HIrisPlex-S DNA test represents the first forensically validated tool for skin colour prediction, and reflects the first forensically validated tool for simultaneous eye, hair and skin colour prediction from DNA.

A simple strategy for improving bird sexing from highly degraded DNA samples

A simple methodology was developed to select new sex-specific primers for bird sexing from degraded and low-quantity DNA sources. The strategy was validated using highly degraded DNA extracted from Giemsa-stained blood smears of common cranes (Grus grus). The new primers allowed the accurate molecular sexing using (i) a classic approach of PCR followed by agarose gel electrophoresis and (ii) an advanced real-time PCR method. The simplicity, speed and low cost make this methodology a versatile molecular tool for selection of novel markers/primers for bird sex differentiation from complex DNA sources, which can be used as basis or complement in several fields of ornithological research.

Performance of the Early Access AmpliSeq™ Mitochondrial Panel with degraded DNA samples using the Ion Torrent™ platform.

The Early Access AmpliSeq™ Mitochondrial Panel amplifies whole mitochondrial genomes for phylogenetic and kinship identifications, using Ion Torrent™ technology. There is currently limited information on its performance with degraded DNA, a common occurrence in forensic samples. This study evaluated the performance of the Panel with DNA samples degraded in vitro, to mimic conditions commonly found in forensic investigations. Purified DNA from five individuals was heat-treated at five time points each (125°C for 0, 30, 60, 120, and 240 min; total n = 25). The quality of DNA was assessed via a real-time DNA assay of genomic DNA and prepared for massively parallel sequencing on the Ion Torrent™ platform. Mitochondrial sequences were obtained for all samples and had an amplicon coverage averaging between 66X to 2803X. Most amplicons (157/162) displayed high coverages (452 ± 333X), while reads with less than 100X coverage were recorded in five amplicons only (90 ± 5X). Amplicon coverage was decreased with prolonged heating. At 72% strand balance, reads were well balanced between forward and reverse strands. Using a coverage threshold of ten reads per SNP, complete sequences were recovered in all samples and resolved kinship and, haplogroup relations. Additionally, the HV1 and HV2 regions of the reference and 240-min heat-treated samples (n = 10) were Sanger-sequenced for concordance. Overall, this study demonstrates the efficacy of a novel forensic Panel that recovers high quality mitochondrial sequences from degraded DNA samples.

Analysis of 30 INDEL Polymorphic Markers in the Panamanian Population: Gene Admixture Estimates, Population Structure and Forensic Parameters

Short tandem repeat (STR) loci have been successfully employed for forensic genetic, paternity and anthropological analyses in the admixed and Amerindian populations of Panamanian for several years. Nevertheless, reports indicate that the use of STRs might be limited in cases involving degraded DNA samples due to their PCR amplicon size and in paternity because of their relatively high mutation rates. Therefore, as a complement to STRs, markers with higher PCR efficiency, such as insertion/deletion (INDEL), have been developed. However, the genetic polymorphisms and distribution of INDELs in the Panamanian population was unknown. Using the Investigator® DIPplex kit (Qiagen), we report here for the first time the genetic profile of 30 INDEL markers in the Panamanian population. Gene admixture estimates and population structure indicated that the Panamanian population is differentially admixed highly polymorphic. Interestingly, admixture estimates where highly similar to our previous report using STR indicating Amerindian, African and European gene contributions. INDELs showed three gene clusters in these proportions: 0.46 (cluster 1), 0.24 (cluster 2) and 0.30 (cluster 3) versus 0.51 (Amerindian), 0.24 (African) and 0.25 (European) for STRs. We also found that both, INDELs and STRs indicated that ancestral gene distribution is heterogeneous among the provinces, across the country. Furthermore, allele frequency variation showed that all loci accomplished Hardy-Weinberg expectations and display high diversity as indicated by the average of observed heterozygosity (0.4477), expected heterozygosity (0.4553) and high similarity with reference U.S. and world populations. Additionally, forensic statistic parameters showed a strong significance in the combined power of discrimination (0.9999999999987) and the combined matching probability (1.27 × 10-12) but relatively lower significance in the combined power of exclusion (0.992655332). The forensic parameters calculated indicate that INDEL markers can be affectively applicable to the Panamanian population for forensic uses but should be complemented with additional markers, such as STRs for paternity analyses.

Development and validation of open-source software for DNA mixture interpretation based on a quantitative continuous model.

In criminal investigations, forensic scientists need to evaluate DNA mixtures. The estimation of the number of contributors and evaluation of the contribution of a person of interest (POI) from these samples are challenging. In this study, we developed a new open-source software “Kongoh” for interpreting DNA mixture based on a quantitative continuous model. The model uses quantitative information of peak heights in the DNA profile and considers the effect of artifacts and allelic drop-out. By using this software, the likelihoods of 1–4 persons’ contributions are calculated, and the most optimal number of contributors is automatically determined; this differs from other open-source software. Therefore, we can eliminate the need to manually determine the number of contributors before the analysis. Kongoh also considers allele- or locus-specific effects of biological parameters based on the experimental data. We then validated Kongoh by calculating the likelihood ratio (LR) of a POI’s contribution in true contributors and non-contributors by using 2–4 person mixtures analyzed through a 15 short tandem repeat typing system. Most LR values obtained from Kongoh during true-contributor testing strongly supported the POI’s contribution even for small amounts or degraded DNA samples. Kongoh correctly rejected a false hypothesis in the non-contributor testing, generated reproducible LR values, and demonstrated higher accuracy of the estimated number of contributors than another software based on the quantitative continuous model. Therefore, Kongoh is useful in accurately interpreting DNA evidence like mixtures and small amounts or degraded DNA samples.

Using long ssDNA polynucleotides to amplify STRs loci in degraded DNA samples.

Obtaining informative short tandem repeat (STR) profiles from degraded DNA samples is a challenging task usually undermined by locus or allele dropouts and peak-high imbalances observed in capillary electrophoresis (CE) electropherograms, especially for those markers with large amplicon sizes. We hereby show that the current STR assays may be greatly improved for the detection of genetic markers in degraded DNA samples by using long single stranded DNA polynucleotides (ssDNA polynucleotides) as surrogates for PCR primers. These long primers allow a closer annealing to the repeat sequences, thereby reducing the length of the template required for the amplification in fragmented DNA samples, while at the same time rendering amplicons of larger sizes suitable for multiplex assays. We also demonstrate that the annealing of long ssDNA polynucleotides does not need to be fully complementary in the 5’ region of the primers, thus allowing for the design of practically any long primer sequence for developing new multiplex assays. Furthermore, genotyping of intact DNA samples could also benefit from utilizing long primers since their close annealing to the target STR sequences may overcome wrong profiling generated by insertions/deletions present between the STR region and the annealing site of the primers. Additionally, long ssDNA polynucleotides might be utilized in multiplex PCR assays for other types of degraded or fragmented DNA, e.g. circulating, cell-free DNA (ccfDNA).

Measuring human DNA degradation and gender detection in forensic DNA samples by q-PCR/HRM analysis

Human DNA quantification, DNA degradation assessment and gender determination are key aspects in most field of human DNA analysis. The assay reported here is a tri-plex Real Time quantitative PCR reaction followed by high resolution melting (HRM) using Syto9 as intercalating dye. The system produces three amplicons: 1- transducin (beta)-like 1, Y-linked −TBL1Y (84bp), 2- DeGraded small target DNA–DGst- (152bp) and 3-DeGraded large target DNA-DGlt- (244bp). DNA quantitation is based on total fluorescence; TBL1Y amplicon allows detecting male DNA and the ratio DGst/DGlt to assess DNA degradation level. q-PCR quantitation proved good linearity in triplicates among 3.2pg/ul–50ng/ul DNA concentration range. Amplification efficiency (E) and reaction slope (m) mean values were 1.04 and 3.23 respectively. Upon HRM analysis, three melting peaks are detected in a male DNA sample and two if only female DNA is present. We define the parameter D as the ratio DGst/DGlt that reflects the extent of DNA degradation in a given sample. A direct correlation has been demonstrated between DNA damage and increased value of parameter D. This q-PCR approach is rapid, sensitive, and a cost-effective method suitable for detecting degraded DNA samples and applicable to any field where human DNA quantitation-qualification is required.

Degraded DNA samples made informative by using superprimers

The use of long ssDNA polynucleotides in PCR of degraded DNA samples can discriminate Short Tandem Repeat (STR) genotypes otherwise missed by current commercial kits.

Analysis of degraded casework DNA by redesigning a mini Y-STR multiplex

As partial genetic profile from degraded DNA samples of a case limiting the application of Y-STRs in forensic sciences, a new multiplex mini Y-STR typing was developed to obtain more useful information and complement current commercially available multiplex kits. The multiplex mini Y-STR amplification including three Y-STR loci (DYS385a/b, DYS389 II, and DYS392) was performed to generate amplicons smaller than those of PowerPlex® Y system over 140bp by redesigning primers. Allelic ladders for three mini Y-STR loci were created using a dilution of allelic ladders from the PowerPlex® Y system. A sensitivity test using serially diluted standard 9948 male DNA showed that three mini Y-STR loci were reliable at template concentrations as low as 0.25ng. By comparing of the multiplex mini Y-STR amplification and PowerPlex® Y system with referenced DNA samples, the results indicated stable concordance in genotyping. And then the successful analysis for the degraded DNA samples from the forensic case produced complete genetic profiles with these three loci. So the multiplex mini Y-STR typing combining the aforementioned commercially available kit increased chance of gaining of 16 Y-STR loci profiles within some core loci in Y Chromosome Haplotype Reference Database from forensic casework sample. It is an alternative strategy for recovering locus-specific information from degraded to redesign the suitable mini Y-STR complement besides new commercial Y-STR kits.

Sequence analysis of ABO and its homologues is valid for species identification.

ABO and its paralogues, such as A3GALT2 and GGTA1, encoding α1,3-Gal(NAc) transferases, belong to the glycosyltransferase 6 (GT6) gene family. We have developed an alternative method for the identification of species based on sequence variations within the GT6 gene family, which is applicable to degraded DNA. DNA samples prepared from control mammalian species, together with an unknown sample, were polymerase chain reaction (PCR)-amplified using one universal primer pair targeting the sequences in the last coding exons of the GT6 gene family, yielding 141-bp products derived from those multiple loci. After cloning, sequence determination and Basic Local Alignment Search Tool analysis, phylogenetic trees were constructed. Comparison of the sequences obtained with those references showed good concordance with each of the starting species of mammals. This system was able to identify 'mouse' or 'rodent' as the origin of the unknown sample. For the identification of species, genotyping of ABO and its homologues would be applicable for the analysis of degraded DNA samples. Although the method employed in this study is likely valid for mammals, it would not be suitable for birds, fish and reptiles. It may be possible to improve the present method for use with other species by employing an alternative universal primer set.

A preliminary evaluation study of new generation multiplex STR kits comprising of the CODIS core loci and the European Standard Set loci

The GlobalFiler™ (Life Technologies), Investigator® 24plex QS (Qiagen), and PowerPlex® Fusion 6C (Promega) kits are the latest generation 6-dye fluorescent chemistry STR-PCR amplification kits. These kits allow for the simultaneous amplification of the CODIS core loci and the European Standard Set loci, as well as a few Y-STR loci in addition to the standard sex-determining marker Amelogenin. The present study was designed to be a preliminary evaluation of the three STR-PCR kits in terms of sensitivity, profile recovery from degraded DNA samples, tolerance to PCR inhibitors, and detection of minor components in DNA mixtures. The results showed that the three STR-PCR kits had relatively similar performance with each kit faring better for the different aspects studied. The PowerPlex® Fusion 6C and the Investigator® 24plex QS kits were shown to tolerate inhibitors better, while the GlobalFiler™ kit appeared to have a higher mean percentage recovery of alleles from low template DNA samples and for minor components in DNA mixtures.

Rapid and affordable genome-wide bisulfite DNA sequencing by XmaI-reduced representation bisulfite sequencing.

To develop a reduced representation bisulfite sequencing (RRBS) approach for rapid and affordable genome-wide DNA methylation analysis. We have selected restriction endonuclease XmaI to produce RRBS library fragments. After digestion and partial fill-in DNA fragments were ligated to barcoded adapters, bisulfite converted, size-selected, and sequenced on the Ion Torrent Personal Genome Machine. XmaI-RRBS results were compared with the previously published RRBS data. We have developed an XmaI-RRBS method for rapid and affordable genome-wide DNA methylation analysis, with library preparation taking only 4 days and sequencing possible within 4 h. We have also addressed several challenges in order to further improve the RRBS technology. XmaI-RRBS may be performed on degraded DNA samples and is compatible with the bench-top next-generation sequencing machines.