SNP and mutation discovery using base-specific cleavage and MALDI-TOF mass spectrometry

MOTIVATION: Single Nucleotide Polymorphisms (SNPs) are believed to contribute strongly to the genetic variability in living beings, in particular their disease or drug side effect predispositions. Mutation-induced sequence variations are playing an important role in the development of cancer, among others. From this, it is clear that SNP and mutation discovery is of great interest in today's Life Sciences. Currently, such discovery is often performed utilizing electrophoresis-based Sanger Sequencing. Discovery of SNPs can also be performed by multiple sequence alignment of publicly available sequence data, but recent studies indicate that only a small percentage of SNPs can be discovered using this approach and, in particular, that SNPs with low frequency are often missed. Other SNP discovery methods only indicate the presence of a SNP in a sample region, but fail to resolve its characterization and localization. RESULTS: We present a method to discover mutations and SNPs using base-specific cleavage and mass spectrometry. An amplicon of known reference sequence with length usually between 100 and 1000 nt is amplified, transcribed, and cleaved using base-specific endonucleases such as RNAse A or T1. The resulting cleavage products (or fragments) are analyzed by MALDI-TOF mass spectrometry and, comparing the measured spectra with those predicted in-silico, the goal is to discover and pinpoint sequence variations of the sample sequence compared to the reference sequence. A time-efficient algorithm for discovering sequence variations is presented that enables fast analysis of such variations even if the sample sequence differs significantly from the reference sequence.     

Interactions of Escherichia coli RNA with bacteriophage MS2 coat protein: genomic SELEX

Genomic SELEX is a method for studying the network of nucleic acid–protein interactions within any organism. Here we report the discovery of several interesting and potentially biologically important interactions using genomic SELEX. We have found that bacteriophage MS2 coat protein binds several Escherichia coli mRNA fragments more tightly than it binds the natural, well-studied, phage mRNA site. MS2 coat protein binds mRNA fragments from rffG (involved in formation of lipopolysaccharide in the bacterial outer membrane), ebgR (lactose utilization repressor), as well as from several other genes. Genomic SELEX may yield experimentally induced artifacts, such as molecules in which the fixed sequences participate in binding. We describe several methods (annealing of oligonucleotides complementary to fixed sequences or switching fixed sequences) to eliminate some, or almost all, of these artifacts. Such methods may be useful tools for both randomized sequence SELEX and genomic SELEX.     

Sequencing from compomers: using mass spectrometry for DNA de novo sequencing of 200+ nt.

One of the main endeavors in today's life science remains the efficient sequencing of long DNA molecules. Today, most de novo sequencing of DNA is still performed using the electrophoresis-based Sanger concept of 1977, in spite of certain restrictions of this method. Methods using mass spectrometry to acquire the Sanger sequencing data are limited by short sequencing lengths of 15-25 nt. We propose a new method for DNA sequencing using base-specific cleavage and mass spectrometry that appears to be a promising alternative to classical DNA sequencing approaches. A single stranded DNA or RNA molecule is cleaved by a base-specific (bio-)chemical reaction using, for example, RNAses. The cleavage reaction is modified such that not all, but only a certain percentage of bases are cleaved. The resulting mixture of fragments is then analyzed using MALDI-TOF mass spectrometry, whereby we acquire the molecular masses of fragments. For every peak in the mass spectrum, we calculate those base compositions that will potentially create a peak of the observed mass and, repeating the cleavage reaction for all four bases, finally try to uniquely reconstruct the underlying sequence from these observed spectra. This leads us to the combinatorial problem of sequencing from compomers and, finally, to the graph-theoretical problem of finding a walk in a subgraph of the de Bruijn graph. Application of this method to simulated data indicates that it might be capable of sequencing DNA molecules with 200+ nt.     

Targeted Genome-Wide Enrichment of Functional Regions

Only a small fraction of large genomes such as that of the human contains the functional regions such as the exons, promoters, and polyA sites. A platform technique for selective enrichment of functional genomic regions will enable several next-generation sequencing applications that include the discovery of causal mutations for disease and drug response. Here, we describe a powerful platform technique, termed “functional genomic fingerprinting” (FGF), for the multiplexed genomewide isolation and analysis of targeted regions such as the exome, promoterome, or exon splice enhancers. The technique employs a fixed part of a uniquely designed Fixed-Randomized primer, while the randomized part contains all the possible sequence permutations. The Fixed-Randomized primers bind with full sequence complementarity at multiple sites where the fixed sequence (such as the splice signals) occurs within the genome, and multiplex amplify many regions bounded by the fixed sequences (e.g., exons). Notably, validation of this technique using cardiac myosin binding protein-C (MYBPC3) gene as an example strongly supports the application and efficacy of this method. Further, assisted by genomewide computational analyses of such sequences, the FGF technique may provide a unique platform for high-throughput sample production and analysis of targeted genomic regions by the next-generation sequencing techniques, with powerful applications in discovering disease and drug response genes.     

High-resolution human papillomavirus genotyping by MALDI-TOF mass spectrometry

We describe a matrix-assisted laser desorption/ionization-time of flight (MALDI-TOF) mass spectrometry (MS)-based assay for human papillomavirus (HPV) genotyping--the restriction fragment mass polymorphism (RFMP) assay, which is based on mass measurement of genotype-specific oligonucleotide fragments generated by TypeIIS restriction endonuclease cleavage after recognition sites have been introduced by PCR amplification. The use of a TypeIIS restriction enzyme makes the RFMP assay independent of sequence and applicable to a wide variety of HPV genotypes, because these enzymes have cleavage sites at a fixed distance from their recognition sites. After PCR amplification, samples are subjected to restriction enzyme digestion with FokI and BtsCI and desalting using Oasis purification plates, followed by analysis by MALDI-TOF MS. Overall, the protocol is simple, takes approximately 4-4.5 h and can accurately detect and identify at least 74 different HPV genotypes.     

Target site selection for an RNA-cleaving catalytic DNA.

A small catalytic DNA, known as the 10-23 DNA enzyme or deoxyribozyme, has been shown to efficiently hydrolyze RNA at purine-pyrimidine (R-Y) junctions in vitro. Although these potentially cleavable junctions are ubiquitous, they are often protected from deoxyribozyme activity by RNA secondary structure. We have developed a multiplex cleavage assay for screening the entire length of a target RNA molecule for deoxyribozyme cleavage sites that are efficient, both in terms of kinetics and accessibility. This strategy allowed us to simultaneously compare the RNA cleaving activity of 80 deoxyribozymes for a model target gene (HPV16 E6), and an additional 60 deoxyribozymes against the rat c-myc target. The human papilloma virus (HPV) target was used primarily to characterize the multiplex system and determine its validity. The c-myc target, coupled with a smooth muscle cell proliferation assay, allowed us to assess the relationship between in vitro cleavage efficiency and c-myc gene suppression in cell culture. The multiplex reaction approach streamlines the process of revealing effective deoxyribozymes in a functional assay and provides accessibility data that may also be applicable to site selection for other hybridization-based agents.     

Structure of BamHI bound to nonspecific DNA: a model for DNA sliding

The central problem faced by DNA binding proteins is how to select the correct DNA sequence from the sea of nonspecific sequences in a cell. The problem is particularly acute for bacterial restriction enzymes because cleavage at an incorrect DNA site could be lethal. To understand the basis of this selectivity, we report here the crystal structure of endonuclease BamHI bound to noncognate DNA. We show that, despite only a single base pair change in the recognition sequence, the enzyme adopts an open configuration that is on the pathway between free and specifically bound forms of the enzyme. Surprisingly, the DNA drops out of the binding cleft with a total loss of base-specific and backbone contacts. Taken together, the structure provides a remarkable snapshot of an enzyme poised for linear diffusion (rather than cleavage) along the DNA.     

Association mapping with single-feature polymorphisms

We develop methods for exploiting "single-feature polymorphism" data, generated by hybridizing genomic DNA to oligonucleotide expression arrays. Our methods enable the use of such data, which can be regarded as very high density, but imperfect, polymorphism data, for genomewide association or linkage disequilibrium mapping. We use a simulation-based power study to conclude that our methods should have good power for organisms like Arabidopsis thaliana, in which linkage disequilibrium is extensive, the reason being that the noisiness of single-feature polymorphism data is more than compensated for by their great number. Finally, we show how power depends on the accuracy with which single-feature polymorphisms are called.     

Arrayed primer extension on in situ synthesized 5'-->3' oligonucleotides in microchannels

Arrays of oligonucleotides synthesized in the 5'-->3' direction have potential benefit in several areas of life sciences research because the free 3'-end can be modified by enzymatic reactions. A Geniom One instrument (febit biomed GmbH, Germany), with integrated chip fabrication, multiplex primer extension, fluorescence imaging, and data analysis, was evaluated for studies of genomic variations. Microchannels used for the array synthesis in Geniom One were not optimized before for the APEX method and, as preliminary experiments demonstrated in this study, the signals were strongly affected by the speed of the process inside reaction channels. Using the two-compartment model (TCM), target binding to feature were quantitatively analyzed, revealing profound mass-transport limitations in the observed kinetics and enabling us to draw a series of physicochemical conclusions of the optimal set-up for the APEX reaction. Some kinetically relevant parameters such as target concentration, reaction time, and temperature were comprehensively analyzed. Finally, we applied the arrays and methods in a proof-of-principle experiment where 36 individuals were typed with 900 oligonucleotide probes (sense and antisense primers for 450 markers), using the ABCR gene as a test system. A new DNA analysis method for studies of genomic variation was developed using this all-in-one platform.     

Mass spectrometry-based proteomics of human cannabinoid receptor 2: covalent cysteine 6.47(257)-ligand interaction affording megagonist receptor activation

The lack of experimental characterization of the structures and ligand-binding motifs of therapeutic G-protein coupled receptors (GPCRs) hampers rational drug discovery. The human cannabinoid receptor 2 (hCB2R) is a class-A GPCR and promising therapeutic target for small-molecule cannabinergic agonists as medicines. Prior mutational and modeling data constitute provisional evidence that AM-841, a high-affinity classical cannabinoid, interacts with cysteine C6.47(257) in hCB2R transmembrane helix 6 (TMH6) to afford improved hCB2R selectivity and unprecedented agonist potency. We now apply bottom-up mass spectrometry (MS)-based proteomics to define directly the hCB2R-AM-841 interaction at the amino-acid level. Recombinant hCB2R, overexpressed as an N-terminal FLAG-tagged/C-terminal 6His-tagged protein (FLAG-hCB2R-6His) with a baculovirus system, was solubilized and purified by immunochromatography as functional receptor. A multiplex multiple reaction monitoring (MRM)-MS method was developed that allowed us to observe unambiguously all seven discrete TMH peptides in the tryptic digest of purified FLAG-hCB2R-6His and demonstrate that AM-841 modifies hCB2R TMH6 exclusively. High-resolution mass spectra of the TMH6 tryptic peptide obtained by Q-TOF MS/MS analysis demonstrated that AM-841 covalently and selectively modifies hCB2R at TMH6 cysteine C6.47(257). These data demonstrate how integration of MS-based proteomics into a ligand-assisted protein structure (LAPS) experimental paradigm can offer guidance to structure-enabled GPCR agonist design.     

MARA: a novel approach for highly multiplexed locus-specific SNP genotyping using high-density DNA oligonucleotide arrays

We have developed a locus-specific DNA target preparation method for highly multiplexed single nucleotide polymorphism (SNP) genotyping called MARA (Multiplexed Anchored Runoff Amplification). The approach uses a single primer per SNP in conjunction with restriction enzyme digested, adapter-ligated human genomic DNA. Each primer is composed of common sequence at the 5′ end followed by locus-specific sequence at the 3′ end. Following a primary reaction in which locus-specific products are generated, a secondary universal amplification is carried out using a generic primer pair corresponding to the oligonucleotide and genomic DNA adapter sequences. Allele discrimination is achieved by hybridization to high-density DNA oligonucleotide arrays. Initial multiplex reactions containing either 250 primers or 750 primers across nine DNA samples demonstrated an average sample call rate of ~95% for 250- and 750-plex MARA. We have also evaluated >1000- and 4000-primer plex MARA to genotype SNPs from human chromosome 21. We have identified a subset of SNPs corresponding to a primer conversion rate of ~75%, which show an average call rate over 95% and concordance >99% across seven DNA samples. Thus, MARA may potentially improve the throughput of SNP genotyping when coupled with allele discrimination on high-density arrays by allowing levels of multiplexing during target generation that far exceed the capacity of traditional multiplex PCR.     

Quantitative multiplex real-time PCR for detection of PLP1 gene duplications in Pelizaeus-Merzbacher patients

Pelizaeus-Merzbacher disease (PMD) is an X-linked recessive disorder of central nervous system (CNS) myelination typically affecting males. A genomic duplication of variable size at Xq22.2, containing the entire proteolipid protein 1 gene (PLP1), is responsible for approximately 60-70% of PMD cases. The aim of this study was to develop a rapid and robust method for determination of PLP1 gene dosage. We optimized two multiplex real-time quantitative PCR (Q-PCR) assays targeting exons 3 and 6 of the PLP1 gene, and then validated these assays by retrospective analysis of a set of genomic DNAs from 67 previously tested patients and 43 normal controls. Samples were analyzed in multiplex PCR reactions using TaqMan chemistry and the ABI Prism 7000 Sequence Detection System. PLP1 dosage was determined by the relative quantitative comparative threshold cycle method (DeltaDeltaCt) using the human serum albumin gene as the endogenous reference gene. Three clearly non-overlapping ranges of results, corresponding to the presence of one, two, or three PLP1 copies, were detected in both assays. The results were completely concordant with gender and previous PLP1 gene dosage testing based on quantitative fluorescent multiplex PCR and analysis of a dinucleotide polymorphism in the first intron of the PLP1 gene. We conclude that multiplex real-time Q-PCR represents a fast and reliable tool for PLP1 duplication testing in PMD families.     

Signatures of Sex-Antagonistic Selection on Recombining Sex Chromosomes

Sex-antagonistic (SA) selection has major evolutionary consequences: it can drive genomic change, constrain adaptation, and maintain genetic variation for fitness. The recombining (or pseudoautosomal) regions of sex chromosomes are a promising setting in which to study SA selection because they tend to accumulate SA polymorphisms and because recombination allows us to deploy the tools of molecular evolution to locate targets of SA selection and quantify evolutionary forces. Here we use coalescent models to characterize the patterns of polymorphism expected within and divergence between recombining X and Y (or Z and W) sex chromosomes. SA selection generates peaks of divergence between X and Y that can extend substantial distances away from the targets of selection. Linkage disequilibrium between neutral sites is also inflated. We show how the pattern of divergence is altered when the SA polymorphism or the sex-determining region was recently established. We use data from the flowering plant Silene latifolia to illustrate how the strength of SA selection might be quantified using molecular data from recombining sex chromosomes.     

Cot-based sampling of genomes for polymorphic low-copy DNA

DNA polymorphisms are powerful tools for many evolutionary and genomic studies in plants including molecular breeding. Single nucleotide polymorphisms (SNPs) are the most elemental DNA marker for genomic studies, but even with advances in DNA sequencing technology, SNP discovery remains costly and computationally demanding, especially in large genomes that are rich in repetitive DNA such as those of many plants. Here we report a method using DNA renaturation kinetics (Cot techniques), sequencing, and BLAST-based screening to identify low-copy, non-coding DNA sequences that were subsequently found to be relatively rich in polymorphisms. A total of of 63 such fragments isolated from a diploid D genome cotton species (Gossypium raimondii) revealed a higher frequency of polymorphisms than that observed for cotton expressed sequence tags or hypomethylated (PstI-susceptible) genomic DNA. While microsatellite-derived loci show still higher polymorphism rates, they often fall in repetitive elements and their sequence analysis is often complicated by alignment difficulties. The potential applications of Cot-filtered noncoding (CFNC) DNA in development of DNA markers are discussed.     

BuildSummary: using a group-based approach to improve the sensitivity of peptide/protein identification in shotgun proteomics

The target-decoy database search strategy is widely accepted as a standard method for estimating the false discovery rate (FDR) of peptide identification, based on which peptide-spectrum matches (PSMs) from the target database are filtered. To improve the sensitivity of protein identification given a fixed accuracy (frequently defined by a protein FDR threshold), a postprocessing procedure is often used that integrates results from different peptide search engines that had assayed the same data set. In this work, we show that PSMs that are grouped by the precursor charge, the number of missed internal cleavage sites, the modification state, and the numbers of protease termini and that the proteins grouped by their unique peptide count should be filtered separately according to the given FDR. We also develop an iterative procedure to filter the PSMs and proteins simultaneously, according to the given FDR. Finally, we present a general framework to integrate the results from different peptide search engines using the same FDR threshold. Our method was tested with several shotgun proteomics data sets that were acquired by multiple LC/MS instruments from two different biological samples. The results showed a satisfactory performance. We implemented the method in a user-friendly software package called BuildSummary, which can be downloaded for free from http://www.proteomics.ac.cn/software/proteomicstools/index.htm as part of the software suite ProteomicsTools.     

A sensitive mass spectrometry method for simultaneous quantification of DNA methylation and hydroxymethylation levels in biological samples

The recent discovery of 5-hydroxymethyl-cytosine (5hmC) in embryonic stem cells and postmitotic neurons has triggered the need for quantitative measurements of both 5-methyl-cytosine (5mC) and 5hmC in the same sample. We have developed a method using liquid chromatography electrospray ionization tandem mass spectrometry with multiple reaction monitoring (LC–ESI–MS/MS–MRM) to simultaneously measure levels of 5mC and 5hmC in digested genomic DNA. This method is fast, robust, and accurate, and it is more sensitive than the current 5hmC quantitation methods such as end labeling with thin layer chromatography and radiolabeling by glycosylation. Only 50 ng of digested genomic DNA is required to measure the presence of 0.1% 5hmC in DNA from mouse embryonic stem cells. Using this procedure, we show that human induced pluripotent stem cells exhibit a dramatic increase in 5mC and 5hmC levels compared with parental fibroblast cells, suggesting a dynamic regulation of DNA methylation and hydroxymethylation during cellular reprogramming.     

Sequence specific detection of DNA using nicking endonuclease signal amplification (NESA)

We have developed a new method for identifying specific single- or double-stranded DNA sequences called nicking endonuclease signal amplification (NESA). A probe and target DNA anneal to create a restriction site that is recognized by a strand-specific endonuclease that cleaves the probe into two pieces leaving the target DNA intact. The target DNA can then act as a template for fresh probe and the process of hybridization, cleavage and dissociation repeats. Laser-induced fluorescence coupled with capillary electrophoresis was used to measure the probe cleavage products. The reaction is rapid; full cleavage of probe occurs within one minute under ideal conditions. The reaction is specific since it requires complete complementarity between the oligonucleotide and the template at the restriction site and sufficient complementarity overall to allow hybridization. We show that both Bacillus subtilis and B. anthracis genomic DNA can be detected and specifically differentiated from DNA of other Bacillus species. When combined with multiple displacement amplification, detection of a single copy target from less than 30 cfu is possible. This method should be applicable whenever there is a requirement to detect a specific DNA sequence. Other applications include SNP analysis and genotyping. The reaction is inherently simple to multiplex and is amenable to automation.     

Global substrate specificity profiling of post‐translational modifying enzymes

Enzymes that modify the proteome, referred to as post‐translational modifying (PTM) enzymes, are central regulators of cellular signaling. Determining the substrate specificity of PTM enzymes is a critical step in unraveling their biological functions both in normal physiological processes and in disease states. Advances in peptide chemistry over the last century have enabled the rapid generation of peptide libraries for querying substrate recognition by PTM enzymes. In this article, we highlight various peptide‐based approaches for analysis of PTM enzyme substrate specificity. We focus on the application of these technologies to proteases and also discuss specific examples in which they have been used to uncover the substrate specificity of other types of PTM enzymes, such as kinases. In particular, we highlight our multiplex substrate profiling by mass spectrometry (MSP‐MS) assay, which uses a rationally designed, physicochemically diverse library of tetradecapeptides. We show how this method has been applied to PTM enzymes to uncover biological function, and guide substrate and inhibitor design. We also briefly discuss how this technique can be combined with other methods to gain a systems‐level understanding of PTM enzyme regulation and function.     

Automated band mapping in electrophoretic gel images using background information

Some popular methods for polymorphism and mutation discovery involve ascertainment of novel bands by the examination of electrophoretic gel images. Although existing strategies for mapping bands work well for specific applications, such as DNA sequencing, these strategies are not well suited for novel band detection. Here, we describe a general strategy for band mapping that uses background banding patterns to facilitate lane calling and size calibration. We have implemented this strategy in GelBuddy, a user-friendly Java-based program for PC and Macintosh computers, which includes several utilities to assist discovery of mutations and polymorphisms. We demonstrate the use of GelBuddy in applications based on single-base mismatch cleavage of heteroduplexed PCR products. Use of software designed to facilitate novel band detection can significantly shorten the time needed for image analysis and data entry in a high-throughput setting. Furthermore, the interactive strategy implemented in GelBuddy has been successfully applied to DNA fingerprinting applications, such as AFLP. GelBuddy promises to make electrophoretic gel analysis a viable alternative to DNA resequencing for discovery of mutations and polymorphisms.