# Matrix-Assisted Laser Desorption/Ionisation, Time-of-Flight Mass Spectrometry in Genomics Research

Jiannis Ragoussis, Gareth P. Elvidge, Kulvinder Kaur, Stefano Colella

► To cite this version:

Jiannis Ragoussis, Gareth P. Elvidge, Kulvinder Kaur, Stefano Colella. Matrix-Assisted Laser Desorption/Ionisation, Time-of-Flight Mass Spectrometry in Genomics Research. PLoS Genetics, 2006, 2 (7), 10 p. 10.1371/journal.pgen.0020100. [https://doi.org/10.1371/journal.pgen.0020100](https://doi.org/10.1371/journal.pgen.0020100)

HAL Id: hal-02664834

https://hal.inrae.fr/hal-02664834

Submitted on 31 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Review

# Matrix-Assisted Laser Desorption/Ionisation, Time-of-Flight Mass Spectrometry in Genomics Research

Jiannis Ragoussis*, Gareth P. Elvidge, Kulvinder Kaur, Stefano Colella

## ABSTRACT

The beginning of this millennium has seen dramatic advances in genomic research. Milestones such as the complete sequencing of the human genome and of many other species were achieved and complemented by the systematic discovery of variation at the single nucleotide (SNP) and whole segment (copy number polymorphism) level. Currently most genomics research efforts are concentrated on the production of whole genome functional annotations, as well as on mapping the epigenome by identifying the methylation status of CpGs, mainly in CpG islands, in different tissues. These recent advances have a major impact on the way genetic research is conducted and have accelerated the discovery of genetic factors contributing to disease. Technology was the critical driving force behind genomics projects: both the combination of Sanger sequencing with high-throughput capillary electrophoresis and the rapid advances in microarray technologies were keys to success. MALDI-TOF MS-based genome analysis represents a relative newcomer in this field. Can it establish itself as a long-term contributor to genetics research, or is it only suitable for niche areas and for laboratories with a passion for mass spectrometry? In this review, we will highlight the potential of MALDI-TOF MS-based tools for resequencing and for epigenetics research applications, as well as for classical complex genetic studies, allele quantification, and quantitative gene expression analysis. We will also identify the current limitations of this approach and attempt to place it in the context of other genome analysis technologies.

## Introduction

It is now commonly accepted that genetics plays a key role in the occurrence of common diseases, from diabetes and cancer to psychiatric disorders. At the same time, genetic approaches have been mostly successful in identifying the genetic component of Mendelian diseases. The construction of physical and genetic maps for the human genome [1-3] was followed by the complete human genome sequence and analysis [4]. This has allowed the experimental [5] and computational identification of transcribed sequences, and driven the execution of cost-effective positional cloning approaches, leading to the identification of about 1,000 genes causing Mendelian diseases. The completion of the human sequence was rapidly followed by the identification of the most common variations in the human genome, namely single nucleotide polymorphisms (SNPs), estimated at 10 million [6].

These major advances are being complemented by efforts aimed at the detection of the methylation status of CpG islands (Epigenome Project) [7] and by large-scale gene expression studies aiming to link gene expression patterns to genome variation [8].

The underlying force for these major advances is a collection of enabling technologies, spearheaded by Sanger sequencing combined with capillary electrophoresis and microarray-based technologies. Current efforts are concentrated towards dramatically improving the cost efficiency of whole genome sequencing approaches to achieve the goal of sequencing one genome a day at $100,000 each, and potentially for $10,000 or less [9].

One of the technologies that has been put forward as an alternative to the Sanger sequencing/capillary electrophoresis combination is matrix-assisted laser desorption/ionisation, time-of-flight mass spectrometry (MALDI-TOF MS). MALDI-TOF MS was introduced, and has established itself, as the tool of choice in proteomics applications, while the full potential for DNA analysis was demonstrated in 1995 [10] and for RNA in 1998 [11].

In brief, for MALDI-TOF MS analysis single-stranded nucleic acid molecules of 3-29 bp in length (1,000-8,600 Dalton range) need to be generated and deposited on a matrix (e.g., 3-hydroxy picolinic acid). The analyte/matrix molecules are then irradiated by a laser inducing their desorption and ionisation, upon which the molecules pass through a flight tube connected to a detector on the other

**Editor:** Elizabeth M. C. Fisher, University College London, United Kingdom

**Citation:** Ragoussis J, Elvidge GP, Kaur K, Colella S (2006) Matrix-assisted laser desorption/ionisation, time-of-flight mass spectrometry in genomics research. PLoS Genet (2)7: e100. DOI: 10.1371/journal.pgen.0020100

**DOI:** 10.1371/journal.pgen.0020100

**Copyright:** © 2006 Ragoussis et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Abbreviations:** hMC, homogeneous MassCLEAVE; hME, homogeneous MassEXTEND reaction (or protocol); indels, insertions and deletions; MALDI-TOF MS, matrix-assisted laser desorption/ionisation time-of-flight mass spectrometry; MSI, microsatellite instability; RC-PCR, real competitive PCR; SNP, single nucleotide polymorphism

Jiannis Ragoussis, Gareth P. Elvidge, Kulvinder Kaur, and Stefano Colella are at the Genomics Laboratory, Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford, United Kingdom.

<figure><img src="image_2.png" /></figure>

<footnote>* To whom correspondence should be addressed. E-mail: Ioannis.Ragoussis@well.ox.ac.uk</footnote>
<figure><img src="image_3.png" /></figure>

DOI: 10.1371/journal.pgen.0020100.g001

<caption>**Figure 1.** MALDI-TOF Mass Spectrometry</caption>

The analyte is deposited on a matrix crystal containing spot. A laser is fired at the matrix leading to the desorption of the analytes and their ionisation. For the majority of DNA analyses, the matrix is designed to produce positive ions. The molecules accelerate into the flight tube where they fly towards the detector on the right. Low mass molecules arrive in a shorter time than heavier ones and molecules of different mass are thus separated (1). After data processing, a spectrum is produced with relative intensity on the y-axis and mass in Daltons on the x-axis (2).

end (see Figure 1). Separation occurs by the time of flight, which is proportional to the mass of the individual molecules. The main advantage of the method is that it directly measures an intrinsic physical property of the molecules, namely their mass, and at a very high speed (about 100 µs). Limitations lie in the size of the DNA molecules that can be detected intact to less than 100 bp (due to size-dependent fragmentation during the MALDI process); and that the analytes must be free from ion adducts which lead to mass distortion. Efforts in the last ten years concentrated on developing simple, robust, homogeneous, and automatable assays suitable for a wide spectrum of genomic applications.

## SNP Genotyping

One of the most effective uses of MALDI-TOF MS in genome analysis is for SNP genotyping. Several methods have been developed for this application based on either a locus-specific PCR step followed by a primer extension, invader reaction, or hybridization steps with or without exonuclease digestion. The products are then detected by MALDI-TOF MS (reviewed in [12]). From all these possibilities, the most widely used are the GOOD assay [13], the PinPoint [14], the GenoSNIP (Bruker Daltonics, Billerica, Massachusetts, United States) and the MassEXTEND (Sequenom, San Diego, California, United States) [15]. The homogeneous MassEXTEND reaction (hME) represented an important step forward since it was the only truly homogeneous assay that is performed in one tube (well) using unmodified primers, and therefore best suited to high throughput analyses. The assay is based on generating PCR products, which are subsequently treated with shrimp alkaline phosphatase, followed by primer extension reactions using Thermosequenase (GE Healthcare, Bucks, United Kingdom) and stop mixes containing cocktails

of ddNTPs and dNTPs. Adduct-forming ions are removed by adding an anion exchange resin. MALDI-TOF analysis is performed by piezo electrically dispensing 15 nl of the reaction products into silica-based chip arrays containing matrix spots (3-hydroxy picolinic acid) in a 384 format. A batch of ten of these chips can be loaded into a typical MALDI instrument (e.g., Bruker Biflex II/III or Autoflex I/II), allowing the analysis of 3,840 reactions in a few hours.

To be cost effective, the SNP assays must be multiplexed at a high level. Although 12-plex reactions have been demonstrated [16], designing multiplex PCR assays involves dealing with the competing parameters of specificity and uniform cycling conditions, and minimising interactions between primer sets.

A step towards optimising multiplexed PCR design was achieved recently through the generation of a set of assay design algorithms that screen out primer interactions and perform simultaneous BLAT searches on the human genome for specificity optimisations, thus allowing automated design of primers for 20–30 plex reactions [17]. Sequenom optimised the primer extension reactions for MALDI-TOF MS by using acyclic mass modified base terminators that achieve at least 16 Dalton gaps between the four possible bases incorporated into a single base extension. Finally, noncomplementary bases are added to the 5' of the extension primers to uniformly fill the 4,500–9,000 Dalton range of the spectrum. In parallel, the PCR and primer extension reaction conditions were optimised. As a result, this type of assay, called iPLEX, features a high design efficiency, low assay dropout rates, and high pass rates (about 95% in our hands) at typically the 26–29 plex level (Figure 2). At the same time, iPLEX is cost efficient, able to offer a very competitive cost in the “few cent” per genotype category, while the ability to process high sample numbers is retained. Due to the high efficiency PCR conditions, extra care is needed to avoid contaminations.

hME is suitable for large-scale collaborative projects requiring a mid-range number of SNPs (10–400) and high numbers of samples [18] and has also contributed to the first phase of the HapMap project. It is also suitable for projects requiring validation steps of ultra-high-throughput SNP-based disease association studies, as well as for genetic studies involving model organisms. A recent example is the development of SNP assays within the framework of the domestic dog genome sequencing project [19]. The ability to design and obtain primers within a very short time enabled us to perform rapid candidate gene analysis by hME using as few as two SNPs and 800 samples up to ten SNPs and 2,215 samples for association studies [20,21]. In addition, 200 SNPs have been analysed in 360 individuals to validate LD data [22]. The hME assay is the cost-effective choice where 2–9 plexes are sufficient on 384 samples or more, while iPLEX is very cost-effective when the assays involve multiples of 25–29 SNPs on 384 samples or more. In both cases, 7,680 samples can be processed within a few days.

Comparable alternative SNP genotyping technologies in terms of throughput include SNPStream from Beckman Coulter (Fullerton, California, United States) [23] and SNPlex from ABI (Foster City, California, United States) [24], both of which have the capability to multiplex 12–48 SNP assays in a 384-well format. SNPStream is based on primer extension and product capture on an array, while SNPlex is based on oligonucleotide ligation and PCR, followed by hybridization

<figure><img src="image_4.png" /></figure>
<figure><img src="image_5.png" /></figure>

DOI: 10.1371/journal.pgen.0020100.g002

<caption>**Figure 2.** Example of a 29plex iPLEX Reaction</caption>

For the iPLEX reactions, the primers are split into three groups according to the position of their respective mass peaks to use the whole spectrum. The concentrations of the three groups are adjusted so that peaks detected at the low mass end of the spectrum have a similar peak height to peaks at the high end of the spectrum. Three examples are highlighted in different colours in different areas of the spectrum. The unextended primer peak is marked by an asterisk and dotted arrow, while the other solid arrows indicate the two different alleles. In red, a homozygous genotype example is shown. In green and blue, two heterozygous alleles are shown. Note the similar peak heights at which extension products corresponding to the two alleles of an SNP are detected and their relative proximity. This proximity enables a higher accuracy to be achieved in the assay by applying the same level of local background subtraction and by eliminating misinterpretations resulting from small peaks, detected just above background, at the expected mass.

to fluorescent probes (ZipChute, Applied Biosystems) and separation by capillary electrophoresis. Equipment costs depend on the level of automation required, but are similar between Sequenom and ABI, while the Beckman Coulter equipment is less expensive. In terms of consumable cost/ genotype, the iPLEX reaction is the most cost-efficient due to the combination of low SNP assay specific and general reagents costs. In terms of pass rate and accuracy, both MassEXTEND and SNPstream were used within the framework of the GenomEUtwin project [18], and comparable results were obtained. The MALDI-TOF-based methods have the advantage over other assays in that the user is in total control of the assay design process and is able to modify procedures and to develop his/her own protocols. The hME/iPLEX type assays have the additional advantage of requiring only unmodified primers, which allows for a minimal gap between SNP selection, design, and project start.

Determining SNP haplotypes is important for many research and clinical applications ranging from association studies to pharmacogenetics. Haplotypes are constructed either using pedigrees and statistical methods, or allele-specific PCR, or somatic cell hybrids [25]. Ivo Gut and colleagues combined allele-specific PCR with MALDI-TOF MS-based genotyping to determine haplotypes within a 2-4 kb range at high throughput [26]. Ding and Cantor took advantage of the sensitivity offered by the MALDI-TOF MS analysis and combined it with high efficiency PCR on highly diluted (3 pg) DNA samples to achieve molecular haplotyping [27]. The method, called M1-PCR, allows the direct determination of haplotypes within up to 24 kb of DNA. Although this latter method is able to generate haplotypes of longer distances than allele-specific PCR-based techniques, it is very sensitive to contamination.

Furthermore, MALDI-TOF MS-based genotyping has also found its way into many clinical genetic applications for SNP

or known mutation detection such as in cystic fibrosis [28-31], and in other clinical applications involving the typing of bacteria and viruses [32,33].

## Quantitative DNA Analysis Applications

In MALDI-TOF MS analyses, the amount of analyte that is co-crystallized with the matrix and then ionized by laser excitation is proportional to the areas under the detected peaks. This allows the relative quantification of the molecules present in the analyte by calculating the peak area ratios. In this setting, variation in peak ratios in a large number of samples can be easily determined, and this quantitative data is one of the strengths of this method. The relative quantitative power of the approach has been used successfully to determine SNP allele frequencies in DNA pools [34-36] in allelotyping experiments. The area under the peaks corresponding to the two alleles is measured and used to determine a ratio, which is normalised to correct the skew due to different extension efficiencies using the ratio obtained from nonpooled genomic DNA (where we expect a 1:1 ratio between the two alleles). For accuracy, each measurement involves at least three technical replicates. MALDI-TOF MS emerged as one of the most accurate allelotyping techniques when compared with Pyrosequencing (Biotage, Uppsala, Sweden), TaqMan, and SNaPshot (Applied Biosystems) [37]. This method has been applied to disease association studies, particularly at a time when genotyping costs were relatively high [38,39]. Even if DNA pooling achieves considerable cost savings and has now been applied on high density Affymetrix GeneChip mapping arrays [40], generating good quality DNA pools requires careful adjustment of each individual DNA, particularly for large numbers of samples (thousands), and therefore becomes labour-intensive. This limitation may be overcome by introducing a whole genome amplification step, whereby the

<figure><img src="image_6.png" /></figure>
end DNA concentration of all samples after amplification is quite homogeneous [41]. The MALDI-TOF MS-based methods utilising DNA pools would be useful for allelotyping studies involving either candidate functional SNPs as shown in Butcher et al. [42] or to validate results from high-density SNP arrays [43].

## Detection of Indels Associated with MSI and Short Tandem Repeats

In addition to SNPs, another important polymorphism is represented by insertions and deletions (indels), some of which can be disease-causing mutations or associated with microsatellite repeat instability (MSI) in the form of length variation of short repetitive sequences. Repeat instability as a form of mutation is linked to a number of disease conditions [44] and cancer [45]. To date, indels and MSIs are detected using electrophoresis and DNA sequencing, which apart from being time-consuming and expensive, are not particularly sensitive in detecting quantitative differences. To solve this problem, MALDI-TOF MS has been applied in combination with assays utilising primer extension [46] or the generation of RNA fragments [47]. Bonk et al. [46] investigated MSIs linked to neoplastic lesions (colorectal tumours, cell lines, and normal control samples) and used primers designed to anneal directly at the 5' of short mononucleotide repeats. Analysis of the primer extension products resulted in specific peak patterns, allowing the identification of unstable repeats in the tumour samples. Another approach was used by Sasayama et al. [47]. Briefly, in vitro transcription of PCR products was followed by cleavage at 5' and 3' of the indel site using molecular scissors and analysis of the resulting RNA fragment by MALDI-TOF, as demonstrated in the analysis of the *ApoE3* gene.

A related approach involving ribozymes was used by Krebs et al. [48] to type short tandem repeat sequences. The method was subsequently refined to utilize base specific cleavage of RNA molecules using RNase T1 [49] (see sequencing below).

## Gene Expression Analysis

Now that the HapMap project has produced working assays for about one million SNPs [50], future efforts will concentrate more on the identification and characterization of functional variants linked to phenotypes [51]. There is great interest in the identification of genetic elements influencing gene expression levels [52]: measuring relative expression levels of specific alleles carrying coding polymorphisms (with potential functional effects) and identifying the effect of DNA polymorphisms on the way hnRNA is spliced to produce mRNA. Once functional variants have been identified, it will be possible to link them to particular phenotypes, covering a broad spectrum from common disease to infection resistance and drug responses. The allelotyping approach described in the quantitative analysis applications can be adapted for allele-specific expression analysis involving coding SNPs or SNPs in UTRs. Briefly, the RNA is converted to cDNA by reverse transcription, followed by PCR of the segment containing the SNP and the interrogation of the SNP using a primer extension-based approach. In heterozygous individuals, we would expect two peaks in MALDI-TOF MS analysis. The area under the peaks corresponding to the two alleles is measured

and used to determine a ratio similar to allelotyping [53]. The method has been used successfully to determine allele-specific expression at the *IL-8* locus [54]. In the absence of a transcribed SNP, it is possible to use the HaploChip technique [55] to assay for SNP alleles present near the transcription start site (within 1 kb) using Polymerase II loading as a marker for transcriptional activity. Antibodies against phosphorylated Pol II are used to immunoprecipitate cross-linked Pol II-DNA complexes. Upon reversal of the cross-link, the segment containing the transcription start site is amplified by PCR and the SNPs within the PCR product are interrogated by primer extension and MALDI-TOF MS. A distortion in the relative ratio of the two alleles indicates an imbalance in their respective levels of transcription.

Total levels of expression can be determined by the real competitive PCR technique (RC-PCR) [56]. This also involves converting RNA to cDNA and amplification of a segment by PCR (usually short amplicons 70-90 bp long). In addition, a competitor oligo of known concentration with the same sequence as the PCR region, apart from one single base close to the middle of the amplicon, is amplified in the same reaction. This type of sequence eliminates differences in the amplification rates between target and competitor and can be used to titrate the amount of target present in the total cDNA (see Figure 3). Typically 7-15 different concentrations of competitor are used and the results plotted in a titration curve. We have applied linear regression analysis using the software package TITAN, and this can be used to analyze data from any type of competitor PCR-based approach [57]. We have compared RC-PCR to real time PCR and found that the main advantage of RC-PCR was the relatively little optimisation required for each assay as compared with SYBR green-based real time PCR. Since the assay uses an internal standard, both competitor and cDNA are amplified under the same conditions and therefore the results are relatively insensitive to nonspecific PCR products or primer-dimers. In addition, the two rounds of amplification (PCR followed by primer extension) allow for two levels of specificity, ensuring that only specific products are detected by MALDI-TOF MS analysis. Disadvantages include the need to prepare dilution series; the addition of reaction components to each well can be complicated and laborious, ideally requiring the use of liquid handling systems and further post-PCR processing. The RC-PCR method is suited to the simultaneous analysis of 70-700 genes rather than one (or small number of) gene(s) in one (or few) sample(s) [57], particularly when combined with iPLEX. The method is also suitable for detecting copy number changes in the genome and is particularly suitable for the validation of expression-profiling results, or multiple signature genes [57]. Finally, primer extension-based techniques can be used effectively in multiplexed assays to identify RNA splice variants by interrogating diagnostic bases at the splice junctions [58].

## DNA and RNA Sequencing

The advantages of MALDI-TOF MS over gel-based separation methods can be summarised as high resolution, high speed, and an absence of compression zones. Considerable amount of work has gone into establishing this technology as a viable alternative to capillary electrophoresis of Sanger sequencing or enzymatic digestion reactions.

<figure><img src="image_7.png" /></figure>
A

<figure><img src="image_8.png" /></figure>

DOI: 10.1371/journal.pgen.0020100.g003

<caption>**Figure 3.** RC-PCR Experiment</caption>

(A) A triplex assay for the genes *BNIP3*, *CA9*, and *NDRG1* is shown. Using a high competitor concentration (10<sup>-12</sup>M), competitor is predominantly PCR-amplified and only the primer extension products resulting from the competitor oligos (blue arrows, top graph) are detected. Using low competitor concentrations (10<sup>-16</sup>M), cDNA is predominantly PCR-amplified and only the primer extension products resulting from the cDNA are detected (blue arrows, bottom graph). The equivalence point, whereby an equal area of competitor and cDNA peaks is derived, would lie at a competitor concentration of between 10<sup>-12</sup> M and 10<sup>-16</sup> M.

(B) Output of the TITAN software (http://www.well.ox.ac.uk/~tprice/titan) showing linear regression analysis of two samples (red and black) for BNIP3 expression. x-axis: competitor concentrations; y-axis (left) is the log 10 ratio of test (cDNA) to competitor signal. The equivalence point is given by the concentration at which y = 0, and is calculated by interpolation (black/red arrows). An approximated difference in expression is detected between the two samples.

B

<figure><img src="image_9.png" /></figure>

Promising developments involving the use of solid phase
capturable dideoxynucleotides, and the generation of RNA
ladders from a DNA template have been reviewed extensively
elsewhere [59,60]. Despite these developments, the read
lengths obtained at a maximum of approximately 100 bp are
not competitive when compared with Sanger sequencing/
capillary electrophoresis [60,61]. As a result, alternative
strategies are being pursued to increase throughput and cut
costs. However, MALDI-TOF MS is well-suited to analysing
RNA [62,63], which is more stable than DNA under MALDI
conditions. In RNA analysis applications, this can be applied
to identify post-transcriptional modifications in ribosomal
RNA. This has been investigated in the context of resistance
towards antibiotics mediated through post-transcriptional
modification of rRNA using RNase digestion-based protocols
[64]. Similar RNase digestion-based protocols have also been
used to investigate RNA secondary structure to optimise the
sequence of antisense oligonucleotides [65]. The assay
involves using oligodeoxynucleotides that are used to form
RNA-oligodeoxynucleotides duplexes which are subsequently
recognised and digested by RNaseH. MALDI-TOF MS is the
method of choice to analyse the resulting fragments.

For DNA analysis, methods have been developed that
involve in vitro transcription of DNA into RNA. RNA
fragments are generated by base specific cleavage and
analysed by MALDI-TOF MS, with the aim of identifying
basepair changes in the sequence as compared with a
computationally predicted outcome [66-69] (Figure 4A).
These methods (for example MassCLEAVE (hMC) by
Sequenom) are suited to resequencing projects for SNP and
frameshift mutations detection (see Figure 4B) as well as SNP
discovery using PCR products with a length of 200-1500 bp.
More complex mutations involving longer indels or

combinations thereof are difficult to determine correctly
without further improvements in the algorithms (K. Kaur and
J. Ragoussis, unpublished data). Although further
bioinformatic optimisations have been incorporated which
allow reactions to be multiplexed, the wide use of capillary
electrophoresis to analyse Sanger-sequencing reaction
products combined with the wealth of algorithms developed
to interpret electropherograms and to detect mutations
mean that it is likely that MALDI-TOF MS-based sequencing
will remain a niche method. Considerable investment will also
be required to compete against microarray-based
resequencing, proven to be effective in detecting point
mutations [70]. However, there is one exciting application
that is a result of the efforts to develop MALDI-TOF MS-
based sequencing approaches and that is methylation
detection.

DNA Methylation Analysis

DNA methylation at CpG sites in the genome contributes to several inherited human genetic diseases [71], and it has been postulated to contribute in a significant way to complex disease susceptibility [72]. The involvement of DNA methylation in gene expression regulation, genome stability, and imprinting has prompted the development and implementation of the Human Epigenome Project. This project will generate a genome-wide DNA methylation profile of the human genome with the analysis of different individuals and multiple tissues [7]. The project is under way and the profile of the human major histocompatibility complex, the target of the Human Epigenome Project pilot phase, has been recently completed [73]. The identification of the epigenetic factors contributing to disease etiology needs

<figure><img src="image_10.png" /></figure>
<figure><img src="image_11.png" /></figure>

DOI: 10.1371/journal.pgen.0020100.g004

<caption>**Figure 4.** Base Specific Cleavage for Mutation Detection</caption>

(A) Illustration of the MassCLEAVE procedure.

(A1) PCR using one primer with a T7 RNA polymerase tag and one with a universal tag.

(A2) In vitro transcription using T7 RNA polymerase and a mixture of dNTPs and NTPs.

(A3) Cleavage by RNAse A. RNAse A cleaves at pyrimidines. By adding one of the pyrimidines as dCTP or dTTP only, the positions containing CTP or TTP are cleaved by RNAse A. Four cleavage reactions are set up in total.

Box to the right of (A): Bioanalyser trace showing a full-length 660 bp RNA produced by IVT (IVT) compared with negative control (CTRL). The final steps involve sample conditioning and analysis by MALDI-TOF MS. Other possibilities include RNAseT1, which cleaves at a G. For methylation analysis, only the two reverse reactions are required.

(B) Spectra and interpretation: A 396-bp fragment containing exon 6 of the *EXT1* gene was analysed by MassCLEAVE to identify a heterozygous deletion of a T at 160 bp. The informative T-forward and C-reverse reaction spectra are shown, whereby the spectrum produced by the sample under investigation is at the top and the control underneath. The numbers at the top indicate expected peak positions along the spectrum and are indicators of reaction quality. In the T-forward reaction, a new peak corresponding to the A1C6G3T1-composed primer appears at 3,425 Dalton (indicated by a red arrow). In the C-reverse reaction, a new peak at 2,699 Daltons is observed, corresponding to the T1G6C1 fragment, while the height of the 3,028 Dalton peak corresponding to the normal A1T1G6C1 fragment is reduced. The relevant normal and mutation-containing fragments are indicated and the deleted T (or A in reverse) is underlined in red. In total, 11/15 mutations in *EXT1* exon 6 and *EXT2* exon 5 (associated with hereditary multiple exostosis) were detectable with this method.

high-throughput technologies to detect DNA methylation. Many technologies have been developed [74], but only some of them lend themselves to high-throughput application. Among these, several are based on the detection of C/T SNPs that are chemically induced at CpG sites in the genome sequence by bisulfite treatment (where a C will be detected if methylation is present) [75]. Direct sequencing has been used, and several technologies developed for SNP detection have been adapted to bisulfite treatment-based methylation detection: including single nucleotide primer extension (SNuPE) combined with dot-blot [76] or with ion pair reverse

phase high performance liquid chromatography (IP RP HPLC) [77], capillary electrophoresis [78,79], pyrosequencing technology [78,80,81], probe-based [82], and melting curve-based real-time PCR methods [83].

Whole genome approaches have been used to globally profile methylation levels including restriction landmark genomic scanning [84] and microarray-based approaches [85]. Methylation has been successfully detected as a chemically induced SNP after bisulfite treatment with specific oligo arrays [86] or using a methylated DNA positive (¹⁰C antibody) or negative (digestion with methylation-sensitive enzymes)

<figure><img src="image_12.png" /></figure>
<figure><img src="image_13.png" /></figure>

DOI: 10.1371/journal.pgen.0020100.g005

<caption>**Figure 5.** Methylation Detection Using MALDI-TOF MS</caption>

MALDI-TOF MS spectra for the detection of different DNA methylation levels (examples shown: 100% methylated DNA; 50:50 mix of methylated and unmethylated DNA; 10% methylated DNA samples mixed with 90% unmethylated DNA) after bisulfite treatment in specific CpG positions in the promoter of *INK4A* (p16) gene (specific amplicon mapping to Chromosome 9: 21964960–21965128 on the Human Genome May 2004 assembly). (A) Primer extension approach using iPLEX chemistry to detect DNA methylation at CpG 2 (Chromosome 9: 21965105–21965106), C call corresponding to the methylated and T call to the nonmethylated alleles.

(B) MassCLEAVE assay using T-reverse specific cleavage, detection of the methylation levels at CpG 10–11 (Chromosome 9: 21964989–21964992 on the Human Genome May 2004 assembly). Fragments detected on the reverse strand are as follows: 4484.85 Dalton 5OH-AAAAAAAAACCACAAT-3p for the unmethylated (Unmeth) and 4516.84 Dalton 5OH-AAAAAAAAACC GCGAT-3p with both CpG sites methylated (Meth).

selection approach on BAC arrays and/or tiling arrays for fine methylation mapping [87,88]. More recently, the Golden Gate genotyping assay on the BeadArray platform (Illumina, San Diego, California, United States) has been successfully adapted to detect methylation [89]. All microarray-based approaches allow the global mapping of overall methylation variation, but they do not offer high-resolution detection at the level of individual CpG sites. Furthermore, these approaches work on the assumption of extended methylation in the region interrogated, which could limit the detection sensitivity when complex unknown patterns of DNA methylation are present.

MALDI-TOF MS approaches, originally developed for SNP genotyping and SNP discovery, have been adapted to the detection of the chemically induced SNP present in the DNA after bisulfite treatment. A primer extension-based approach in the form of the GOOD assay has been successfully developed and used for the quantification of methylation levels [90]. Based on the same principle, the protocol for allelotyping analysis as implemented in the MassARRAY system (Sequenom) (see "Quantitative DNA analysis applications") can be applied to primer extension reactions. We have used both the MassEXTEND protocol (hME) and the iPLEX protocol, to successfully detect and quantify methylation at specific CpG positions (Figure 5A; S. Colella and J. Ragoussis, unpublished data). The major disadvantage

of primer extension-based approaches, irrespective of the detection platform, is the assay design limitation within CpG rich regions, as it is better to avoid having a CpG site (and therefore a potential SNP) in the sequence complementary to the extension primer. This problem can be circumvented now by using peptide nucleic acid probe hybridization, instead of primer extension, to interrogate methylated CpGs. MALDI-TOF analysis of successfully hybridized peptide nucleic acids was shown to produce quantitative data at high sensitivity [91].

An alternative approach is based on sequencing by base specific cleavage as described in the previous section. This method generates different cleavage patterns for the methylated and unmethylated CpG positions present in the fragment analysed. This approach for methylation detection was developed by Schatz et al. using a guanosine-specific cleavage by RNase T1 [92]. More recently, Ehrich et al. have developed an assay using RNase A [69], based on modification of the hMC SNP discovery assay. In this method only the reverse strand is used as a template for the in vitro transcription, followed by the two base specific fragmentations (T-reaction examples in Figure 5B). The spectra obtained from the two reactions are combined and allow the quantification of methylation at CpG sites present in the amplicon analyzed. The final CpG coverage achieved depends on the amplicon size (a maximum of ~800 bp) and the specific sequence content. The average number of CpG

<figure><img src="image_14.png" /></figure>
sites detected in a single amplicon is comparable to sequencing-based approaches.

The RNA cleavage-based approach, e.g., hMC, is quite unique and competitive in rapidly generating a comprehensive CpG island quantitative methylation map. The method is comparable to direct sequencing in terms of throughput, but superior in terms of quantification. The quantification in turn is comparable to a combined cloning and sequencing approach, with the advantage that RNA cleavage is much easier and faster to perform. The RNA cleavage offers the possibility of detecting more CpG sites in a single amplicon, when compared with pyrosequencing, and, furthermore, it is independent from the extension primer design issues. One possible limitation is that some CpG sites cannot be analysed if they are in fragments that have a molecular weight outside the mass detection range. Array-based methods offer the highest throughput in detecting multiple promoter regions in a single experiment, but at present they do not compare to the detailed profile that can be generated with the RNA cleavage-based MALDI-TOF MS approaches. In conclusion, the MALDI-TOF MS resequencing-based method for methylation detection offers a level of resolution that enables high-throughput quantitative analysis of full CpG islands. Data generated with this approach can lead to a better definition of complex methylation patterns of specific regions and provide a better understanding of the effect of specific methylation events in the human genome.

## Discussion

It is clear that MALDI-TOF MS is a versatile tool for genomic research, allowing the analysis of DNA as well as RNA fragments. Due to the high sensitivity and ability to quantify, it is also suitable for clinical diagnostic use. In terms of number of applications, it is comparable to a DNA sequencer. The main advantages are that the technology, by measuring the mass of DNA molecules, allows a high degree of freedom for the individual investigator to conceive assays where the products can be analysed by MALDI-TOF MS. Even as a complete commercially available genomic applications system, as offered by Sequenom or Bruker Daltonics, it is relatively straightforward to apply custom assay design. The other advantage is the high sample number that can be analysed in each run, up to 3,840 samples at a time. This makes applications such as allele-specific gene expression and high-resolution methylation analysis very attractive for the study of functional polymorphisms and epigenetic modifications. Since all applications so far are based on a PCR step, developing new tools to assess PCR primers (and where applicable extension primers) interactions could lead to further improvements in efficient assay design and increase the level of multiplexing, particularly in genotyping assays.

For genotyping applications, the level of multiplexing that can be achieved with current methodology has limits. For example, the currently used part of the spectrum is in the 4,500-9,000 Dalton range. If we leave the PCR multiplexing issues aside, it is theoretically possible to accommodate 93 peak trios of an iPLEX reaction (extension primer and two alleles) at 16 Dalton intervals, thus increasing the number of reactions that can be multiplexed at current resolution levels. This level of multiplexing would result in significant genotyping cost reductions.

An additional problem being addressed is that the

instruments (mass spectrometers such as Bruker Biflex II/III or Autoflex and spotting devices) themselves are expensive and maintenance-intensive, while high throughput applications require the use of liquid handling systems. Both Bruker and Sequenom offer compact, benchtop instruments in the form of the Microflex (Bruker Daltonics) and the Compact (Sequenom) to address the cost (and maintenance) issue.

One important advantage of MALDI-TOF MS that has not been fully utilized is the sensitivity of the technology. For example, no more than 15 nl of analyte is used on a MassARRAY chip. This low level of reaction volume would achieve high cost efficiency and could be realized with miniaturized nanoliter scale devices (for an example see [93], for a review see [94]). For example, the iPLEX assay takes place in 5 µl or 10 µl volume/sample, where the costs, in broad terms, are equally divided among PCR, extension reaction, and chip. It is easy to calculate that such an approach will result in significant savings by cutting two-thirds of the cost by a factor of ten or more. All methods described so far, such as resequencing and expression profiling, are amenable to the same cost reduction principle, and this would allow compact MALDI-TOF MS instruments to be integrated with miniaturized liquid handling systems.

So which applications are going to be most widely used? We think that the genotyping assays (such as iPLEX) are going to be the main applications. This is because current whole genome association studies using hundreds of thousands of SNPs will produce genomic segments or genes with disease associations requiring further rounds of replication in additional sample collections. Furthermore, additional SNPs must be investigated in order to identify the ones directly associated with disease. These downstream studies will require genotyping technologies that offer high-cost efficiency at medium levels of multiplexing and high sample throughput, currently a role fulfilled by the MALDI-TOF MS-based techniques. If further instrument development produces compact, simple-to-use, and lower-maintenance equipment (for example the benchtop MALDI microMX from Waters, [Milford, Massachusetts, United States]), then the added benefits of being able to perform a wide spectrum of additional assays will be a clear advantage. The ability to use the same platform to perform genomic studies beyond SNP genotyping, such as methylation detection, quantitative gene expression, and allele-specific expression, will be a bonus when equipment decisions have to be made for genomics and functional genomics laboratories.

Furthermore, development of compact systems such as the Microflex or MALDI microMX that can be adapted for use in genomic or proteomic applications, will allow MALDI-TOF MS to expand its position in the future as a versatile platform for basic research and clinical diagnostic applications alike. ■

## Supporting Information

### Accession Numbers

The accession numbers for the following genes from the EntrezGene databank (http://srs.pku.edu.cn/srs/srs) are: *ApoE3* (348), *BNIP3* (664), *CA9* (768), *EXT1* (2131), *IL-8* (3576), *INK4A* (1029), and *NDRG1* (10397).

### Acknowledgments

We would like to thank Helen Butler for useful discussions and Aaron Abbott, Lucy Marcelline, and Lorne Lonie for technical assistance.

<figure><img src="image_15.png" /></figure>
**Author contributions.** JR, GPE, and SC conceived and designed the experiments. GPE, KK, and SC performed the experiments. GPE, KK, and SC analyzed the data. JR and SC wrote the paper.

**Funding.** This work was supported by the Wellcome Trust.

**Competing interests.** The authors have declared that no competing interests exist.

## References

1. Deloukas P, Schuler GD, Gyapay G, Beasley EM, Soderlund C, et al. (1998) A physical map of 30,000 human genes. Science 282: 744–746.

2. McPherson JD, Marra M, Hillier L, Waterston RH, Chinwalla A, et al. (2001) A physical map of the human genome. The International Human Genome Mapping Consortium. Nature 409: 934–941.

3. Bentley DR, Deloukas P, Dunham A, French L, Gregory SG, et al. (2001) The physical maps for sequencing human chromosomes 1, 6, 9, 10, 13, 20, and X. Nature 409: 942–943.

4. Lander ES, Linton LM, Birren B, Nusbaum C, Zody MC, et al. (2001) Initial sequencing and analysis of the human genome. International Human Genome Sequencing Consortium. Nature 409: 860–921.

5. Shoemaker DD, Schadt EE, Armour CD, He YD, Garrett-Engle P, et al. (2001) Experimental annotation of the human genome using microarray technology. Nature 409: 922–927.

6. Sachidanandam R, Weissman D, Schmidt SC, Kakol JM, Stein LD, et al. (2001) A map of human genome sequence variation containing 1.42 million single nucleotide polymorphisms. The International SNP Map Working Group. Nature 409: 928–933.

7. Novik KL, Nimmrich I, Genc B, Maier S, Piepenbrock C, et al. (2002) Epigenomics: Genome-wide study of methylation phenomena. Curr Issues Mol Biol 4: 111–128.

8. Cheung VG, Spielman RS, Ewens KG, Weber TM, Morley M, et al. (2005) Mapping determinants of human gene expression by regional and genome-wide association. Nature 437: 1365–1369.

9. Chan EY (2005) Advances in sequencing technology. Mutat Res 573: 13–40.

10. Tang K, Fu D, Kotter S, Cotter RJ, Cantor CR, et al. (1995) Matrix-assisted laser desorption/ionization mass spectrometry of immobilized duplex DNA probes. Nucleic Acids Res 23: 3126–3131.

11. Tolson DA, Nicholson NH (1998) Sequencing RNA by a combination of exonuclease digestion and urine specific chemical cleavage using MALDI-TOF. Nucleic Acids Res 26: 446–451.

12. Tost J, Gut IG (2002) Genotyping single nucleotide polymorphisms by mass spectrometry. Mass Spectrom Rev 21: 388–418.

13. Sauer S, Lechner D, Berlin K, Plancon C, Heuermann A, et al. (2000) Full flexibility genotyping of single nucleotide polymorphisms by the GOOD assay. Nucleic Acids Res 28: E100.

14. Paracchini S, Arredi B, Chalk R, Tyler-Smith C (2002) Hierarchical high-throughput SNP genotyping of the human Y chromosome using MALDI-TOF mass spectrometry. Nucleic Acids Res 30: e27.

15. Buettow KH, Edmonson M, MacDonald R, Clifford R, Yip P, et al. (2001) High-throughput development and characterization of a genomewide collection of gene-based single nucleotide polymorphism markers by chip-based matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. Proc Natl Acad Sci U S A 98: 581–584.

16. van den Boom D, Beaulieu M, Oeth P, Roth R, Honisch C, et al. (2004) MALDI-TOF MS: A platform technology for genetic discovery. Int J Mass Spectrom 238: 173–188.

17. Rachlin J, Ding C, Cantor C, Kasif S (2005) MuPlex: Multi-objective multiplex PCR assay design. Nucleic Acids Res 33: W544–W547.

18. Silander K, Axelsson T, Widen E, Dahlgren A, Palotie A, et al. (2003) Analysis of genetic variation in the GenomEUtwin project. Twin Res 6: 391–398.

19. Lindblad-Toh K, Wade CM, Mikkelsen TS, Karlsson EK, Jaffe DB, et al. (2005) Genome sequence, comparative analysis, and haplotype structure of the domestic dog. Nature 438: 803–819.

20. Bonora E, Beyer KS, Lamb JA, Parr JR, Klauck SM, et al. (2003) Analysis of reelin as a candidate gene for autism. Mol Psychiatry 8: 885–892.

21. Yamazaki K, McGovern D, Ragoussis J, Paolucci M, Butler H, et al. (2005) Single nucleotide polymorphisms in TNFSF15 confer susceptibility to Crohn's disease. Hum Mol Genet 14: 3499–3506.

22. Lawrence R, Evans DM, Morris AP, Ke X, Hunt S, et al. (2005) Genetically indistinguishable SNPs and their influence on inferring the location of disease-associated variants. Genome Res 15: 1503–1510.

23. Bell PA, Chaturvedi S, Gelfand CA, Huang CY, Kochersperger M, et al. (2002) SNPstream UHT: Ultra-high throughput SNP genotyping for pharmacogenomics and drug discovery. Biotechniques 33 (Supplement): 70–72, 74, 76–77.

24. Tobler AR, Short S, Andersen MR, Paner TM, Briggs JC, et al. (2005) The SNplex genotyping system: A flexible and scalable platform for SNP genotyping. J Biomol Tech 16: 398–406.

25. Crawford DC, Nickerson DA (2005) Definition and clinical importance of haplotypes. Annu Rev Med 56: 303–320.

26. Tost J, Brandt O, Boussicault F, Derbala D, Caloustian C, et al. (2002) Molecular haplotyping at high throughput. Nucleic Acids Res 30: e96.

27. Ding C, Cantor CR (2003) Direct molecular haplotyping of long-range genomic DNA with M1-PCR. Proc Natl Acad Sci U S A 100: 7449–7453.

28. Buyse IM, McCarthy SE, Lurix P, Pace RP, Vo D, et al. (2004) Use of MALDI-TOF mass spectrometry in a 51-mutation test for cystic fibrosis: Evidence that 3199del6 is a disease-causing mutation. Genet Med 6: 426–430.

29. Jaremko M, Justenhoven C, Abraham BK, Schroth W, Fritz P, et al. (2005) MALDI-TOF MS and TaqMan assisted SNP genotyping of DNA isolated from formalin-fixed and paraffin-embedded tissues (FFPET). Hum Mutat 25: 232–238.

30. Liao HK, Su YN, Kao HY, Hung CC, Wang HT, et al. (2005) Parallel minisequencing followed by multiplex matrix-assisted laser desorption/ionization mass spectrometry assay for beta-thalassemia mutations. J Hum Genet 50: 139–150.

31. Tost J, Gut IG (2005) Genotyping single nucleotide polymorphisms by MALDI mass spectrometry in clinical applications. Clin Biochem 38: 335–350.

32. Govorun VM, Lokhov PG, Moshkovskii SA, Momynaliev KT, Selesnyova OV, et al. (2004) Comparative analysis of different typing methods for *Helicobacter pylori* clinical isolates. Biochemistry (Mosc) 69: 536–541.

33. Kim KH, Yoon MS, Na YJ, Park CS, Oh MR, et al. (2005) Development and evaluation of a highly sensitive human papillomavirus genotyping DNA chip. Gynecol Oncol 100: 38–43.

34. Mohlke KL, Erdos MR, Scott LJ, Fingerlin TE, Jackson AU, et al. (2002) High-throughput screening for evidence of association by using mass spectrometry genotyping on DNA pools. Proc Natl Acad Sci U S A 99: 16928–16933.

35. Werner M, Sych M, Herbon N, Illig T, König IR, et al. (2002) Large-scale determination of SNP allele frequencies in DNA pools using MALDI-TOF mass spectrometry. Hum Mutat 20: 57–64.

36. Downes K, Barratt BJ, Akan P, Bumpstead SJ, Taylor SD, et al. (2004) SNP allele frequency estimation in DNA pools and variance components analysis. Biotechniques 36: 840–845.

37. Shifman S, Pisante-Shalom A, Yakir B, Darvasi A (2002) Quantitative technologies for allele frequency estimation of SNPs in DNA pools. Mol Cell Probes 16: 429–434.

38. Bansal A, van den Boom D, Kammerer S, Honisch C, Adam G, et al. (2002) Association testing by DNA pooling: An effective initial screen. Proc Natl Acad Sci U S A 99: 16871–16874.

39. Shifman S, Bronstein M, Sternfeld M, Pisante-Shalom A, Lev-Lehman E, et al. (2002) A highly significant association between a COMT haplotype and schizophrenia. Am J Hum Genet 71: 1296–1302.

40. Craig DW, Huentelman MJ, Hu-Lince D, Zismann VL, Kruer MC, et al. (2005) Identification of disease causing loci using an array-based genotyping approach on pooled DNA. BMC Genomics 6: 138.

41. Zhao ZZ, Nyholt DR, James MR, Mayne R, Treloar SA, et al. (2005) A comparison of DNA pools constructed following whole genome amplification for two-stage SNP genotyping designs. Twin Res Hum Genet 8: 353–361.

42. Butcher LM, Meaburn E, Dale PS, Sham P, Schalkwyk LC, et al. (2005) Association analysis of mild mental impairment using DNA pooling to screen 432 brain-expressed single-nucleotide polymorphisms. Mol Psychiatry 10: 384–392.

43. Butcher LM, Meaburn E, Liu L, Fernandes C, Hill L, et al. (2004) Genotyping pooled DNA on microarrays: A systematic genome screen of thousands of SNPs in large samples to detect QTLs for complex traits. Behav Genet 34: 549–555.

44. Pearson CE, Edamura KN, Cleary JD (2005) Repeat instability: Mechanisms of dynamic mutations. Nat Rev Genet 6: 729–742.

45. Ionov Y, Peinado MA, Malkhosyan S, Shibata D, Perucho M (1993) Ubiquitous somatic mutations in simple repeated sequences reveal a new mechanism for colonic carcinogenesis. Nature 363: 558–561.

46. Bonk T, Humeny A, Gebert J, Sutter C, von Knebel Doeberitz M, et al. (2003) Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry-based detection of microsatellite instabilities in coding DNA sequences: A novel approach to identify DNA-mismatch repair-deficient cancer cells. Clin Chem 49: 552–561.

47. Sasayama T, Kato M, Aburatani H, Kuzuya A, Komiyama M (2006) Simultaneous genotyping of indels and SNPs by mass spectroscopy. J Am Soc Mass Spectrom 17: 3–8.

48. Krebs S, Seichter D, Forster M (2001) Genotyping of dinucleotide tandem repeats by MALDI mass spectrometry of ribozyme-cleaved RNA transcripts. Nat Biotechnol 19: 877–880.

49. Seichter D, Krebs S, Forster M (2004) Rapid and accurate characterisation of short tandem repeats by MALDI-TOF analysis of endonuclease cleaved RNA transcripts. Nucleic Acids Res 32: e16.

50. Altshuler D, Brooks LD, Chakravarti A, Collins FS, Daly MJ, et al. (2005) A haplotype map of the human genome. Nature 437: 1299–1320.

51. Gibbs R (2005) Deeper into the genome. Nature 437: 1233–1234.

52. Stranger BE, Forrest MS, Clark AG, Minichiello MJ, Deutsch S, et al. (2005) Genome-wide associations of gene expression variation in humans. PLoS Genet 1: e78.

53. Tang K, Oeth P, Kammerer S, Denissenko MF, Ekblom J, et al. (2004) Mining disease susceptibility genes through SNP analyses and expression profiling using MALDI-TOF mass spectrometry. J Proteome Res 3: 218–227.

54. Hacking D, Knight JC, Rockett K, Brown H, Frampton J, et al. (2004)

<figure><img src="image_16.png" /></figure>
Increased in vivo transcription of an IL-8 haplotype associated with respiratory syncytial virus-disease susceptibility. *Genes Immun* 5: 274-282.

55. Knight JC, Keating BJ, Rockett KA, Kwiatkowski DP (2003) In vivo characterization of regulatory polymorphisms by allele-specific quantification of RNA polymerase loading. *Nat Genet* 33: 469-475.

56. Ding C, Maier E, Roscher AA, Braun A, Cantor CR (2004) Simultaneous quantitative and allele-specific expression analysis with real competitive PCR. *BMC Genet* 5: 8.

57. Elvidge GP, Price TS, Glenny L, Ragoussis J (2005) Development and evaluation of real competitive PCR for high-throughput quantitative applications. *Anal Biochem* 339: 231-241.

58. McCullough RM, Cantor CR, Ding C (2005) High-throughput alternative splicing quantification by primer extension and matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. *Nucleic Acids Res* 33: e99.

59. Gut IG (2004) DNA analysis by MALDI-TOF mass spectrometry. *Hum Mutat* 23: 437-441.

60. Edwards JR, Ruparel H, Ju J (2005) Mass-spectrometry DNA sequencing. *Mutat Res* 573: 3-12.

61. Fu DJ, Broude NE, Koster H, Smith CL, Cantor CR (1996) Efficient preparation of short DNA sequence ladders potentially suitable for MALDI-TOF DNA sequencing. *Genet Anal* 12: 137-142.

62. Kirpekar F, Nordhoff E, Kristiansen K, Roepstorff P, Lezius A, et al. (1994) Matrix assisted laser desorption/ionization mass spectrometry of enzymatically synthesized RNA up to 150 kDa. *Nucleic Acids Res* 22: 3866-3870.

63. Amexis G, Oeth P, Abel K, Ivshina A, Pelloquin F, et al. (2001) Quantitative mutant analysis of viral quasispecies by chip-based matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. *Proc Natl Acad Sci U S A* 98: 12097-12102.

64. Kirpekar F, Douthwaite S, Roepstorff P (2000) Mapping posttranscriptional modifications in 5S ribosomal RNA by MALDI mass spectrometry. *RNA* 6: 296-306.

65. Gabler A, Krebs S, Seichter D, Forster M (2003) Fast and accurate determination of sites along the FUT2 in vitro transcript that are accessible to antisense oligonucleotides by application of secondary structure predictions and RNase H in combination with MALDI-TOF mass spectrometry. *Nucleic Acids Res* 31: e79.

66. Stanssens P, Zabeau M, Meersseman G, Remes G, Gansemans Y, et al. (2004) High-throughput MALDI-TOF discovery of genomic sequence polymorphisms. *Genome Res* 14: 126-133.

67. von Wintzingerode F, Bocker S, Schlotelburg C, Chiu NH, Storm N, et al. (2002) Base-specific fragmentation of amplified 16S rRNA genes analyzed by mass spectrometry: A tool for rapid bacterial identification. *Proc Natl Acad Sci U S A* 99: 7039-7044.

68. Bocker S (2004) Sequencing from compomers: Using mass spectrometry for DNA de novo sequencing of 200+ nt. *J Comput Biol* 11: 1110-1134.

69. Ehrich M, Nelson MR, Stanssens P, Zabeau M, Liloglou T, et al. (2005) Quantitative high-throughput analysis of DNA methylation patterns by base-specific cleavage and mass spectrometry. *Proc Natl Acad Sci U S A* 102: 15785-15790.

70. Maitra A, Cohen Y, Gillespie SE, Mambo E, Fukushima N, et al. (2004) The human MitoChip: A high-throughput sequencing microarray for mitochondrial mutation detection. *Genome Res* 14: 812-819.

71. Robertson KD (2005) DNA methylation and human disease. *Nat Rev Genet* 6: 597-610.

72. Petronis A (2001) Human morbid genetics revisited: Relevance of epigenetics. *Trends Genet* 17: 142-146.

73. Rakyan VK, Hildmann T, Novik KL, Lewin J, Tost J, et al. (2004) DNA methylation profiling of the human major histocompatibility complex: A

pilot study for the human epigenome project. *PLoS Biol* 2: DOI: 10.1371/journal.pbio.0020405

74. Fraga MF, Esteller M (2002). Fraga MF, Esteller M (2002) DNA methylation: A profile of methods and applications. *Biotechniques* 33: 632, 634, 636-649.

75. Clark SJ, Harrison J, Paul CL, Frommer M (1994) High sensitivity mapping of methylated cytosines. *Nucleic Acids Res* 22: 2990-2997.

76. Gonzalgo ML, Jones PA (2002) Quantitative methylation analysis using methylation-sensitive single-nucleotide primer extension (Ms-SNuPE). *Methods* 27: 128-133.

77. El-Maarri O, Herbiniaux U, Walter J, Oldenburg J (2002) A rapid, quantitative, non-radioactive bisulfite-SNuPE-IP RP HPLC assay for methylation analysis at specific CpG sites. *Nucleic Acids Res* 30: e25.

78. Uhlmann K, Brinckmann A, Toliat MR, Ritter H, Nurnberg P (2002) Evaluation of a potential epigenetic biomarker by quantitative methyl- single nucleotide polymorphism analysis. *Electrophoresis* 23: 4072-4079.

79. Kaminsky ZA, Assadzadeh A, Flanagan J, Petronis A (2005) Single nucleotide extension technology for quantitative site-specific evaluation of metC/C in GC-rich regions. *Nucleic Acids Res* 33: e95.

80. Colella S, Shen L, Baggerly KA, Issa JP, Krahe R (2003) Sensitive and quantitative universal Pyrosequencing methylation analysis of CpG sites. *Biotechniques* 35: 146-150.

81. Tost J, Dunker J, Gut IG (2003) Analysis and quantification of multiple methylation variable positions in CpG islands by Pyrosequencing. *Biotechniques* 35: 152-156.

82. Eads CA, Danenberg KD, Kawakami K, Saltz LB, Blake C, et al. (2000) *MethyLight*: A high-throughput assay to measure DNA methylation. *Nucleic Acids Res* 28: E32.

83. Akey DT, Akey JM, Zhang K, Jin L (2002) Assaying DNA methylation based on high-throughput melting curve approaches. *Genomics* 80: 376-384.

84. Hatada I, Hayashizaki Y, Hirotsune S, Komatsubara H, Mukai T (1991) A genomic scanning method for higher organisms using restriction sites as landmarks. *Proc Natl Acad Sci U S A* 88: 9523-9527.

85. van Steensel B, Henikoff S (2003) Epigenomic profiling using microarrays. *Biotechniques* 35: 346-350, 352-354, 356-347.

86. Shi H, Maier S, Nimmrich I, Yan PS, Caldwell CW, et al. (2003) Oligonucleotide-based microarray for DNA methylation analysis: Principles and applications. *J Cell Biochem* 88: 138-143.

87. Mockler TC, Chan S, Sundaresan A, Chen H, Jacobsen SE, et al. (2005) Applications of DNA tiling arrays for whole-genome analysis. *Genomics* 85: 1-15.

88. Weber M, Davies JJ, Wittig D, Oakeley EJ, Haase M, et al. (2005) Chromosome-wide and promoter-specific analyses identify sites of differential DNA methylation in normal and transformed human cells. *Nat Genet* 37: 853-862.

89. Bibikova M, Lin Z, Zhou L, Chudin E, Garcia EW, et al. (2006) High-throughput DNA methylation profiling using universal bead arrays. *Genome Res* 16: 383-393.

90. Tost J, Schatz P, Schuster M, Berlin K, Gut IG (2003) Analysis and accurate quantification of CpG methylation by MALDI mass spectrometry. *Nucleic Acids Res* 31: e50.

91. Schatz P, Distler J, Berlin K, Schuster M (2006) Novel method for high throughput DNA methylation marker evaluation using PNA-probe library hybridization and MALDI-TOF detection. *Nucleic Acids Res* 34: e59.

92. Schatz P, Dietrich D, Schuster M (2004) Rapid analysis of CpG methylation patterns using RNase T1 cleavage and MALDI-TOF. *Nucleic Acids Res* 32: e167.

93. Liu J, Enzelberger M, Quake S (2002) A nanoliter rotary device for polymerase chain reaction. *Electrophoresis* 23: 1531-1536.

94. Zhang C, Xu J, Ma W, Zheng W (2005) PCR microfluidic devices for DNA amplification. *Biotechnol Adv* 24: 243-284.

<figure><img src="image_17.png" /></figure>