This Article OPEN ACCESS ARTICLE Full Text (PDF) Supplemental Research Data All Versions of this Article: gr.070169.107v1 gr.070169.107v2    most recent Alert me when this article is cited Alert me if a correction is posted Services Email this article to a friend Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Google Scholar Articles by Zeller, G. Articles by Rätsch, G. PubMed PubMed Citation Articles by Zeller, G. Articles by Rätsch, G. Social Bookmarking                   What's this?

Methods

Detecting polymorphic regions in Arabidopsis thaliana with resequencing microarrays

¹ Friedrich Miescher Laboratory of the Max Planck Society, Tübingen 72070, Germany; ² Max Planck Institute for Developmental Biology, Department of Molecular Biology, Tübingen 72070, Germany

Whole-genome oligonucleotide resequencing arrays have allowed the comprehensive discovery of single nucleotide polymorphisms (SNPs) in eukaryotic genomes of moderate to large size. With this technology, the detection rate for isolated SNPs is typically high. However, it is greatly reduced when other polymorphisms are located near a SNP as multiple mismatches inhibit hybridization to arrayed oligonucleotides. Contiguous tracts of suppressed hybridization therefore typify polymorphic regions (PRs) such as clusters of SNPs or deletions. We developed a machine learning method, designated margin-based prediction of polymorphic regions (mPPR), to predict PRs from resequencing array data. Conceptually similar to hidden Markov models, the method is trained with discriminative learning techniques related to support vector machines, and accurately identifies even very short polymorphic tracts (<10 bp). We applied this method to resequencing array data previously generated for the euchromatic genomes of 20 strains (accessions) of the best-characterized plant, Arabidopsis thaliana. Nonredundantly, 27% of the genome was included within the boundaries of PRs predicted at high specificity (97%). The resulting data set provides a fine-scale view of polymorphic sequences in A. thaliana; patterns of polymorphism not apparent in SNP data were readily detected, especially for noncoding regions. Our predictions provide a valuable resource for evolutionary genetic and functional studies in A. thaliana, and our method is applicable to similar data sets in other species. More broadly, our computational approach can be applied to other segmentation tasks related to the analysis of genomic variation.

⁴ Corresponding author.

E-mail Gunnar.Raetsch{at}tuebingen.mpg.de; fax 49-7071-601-801.

[Supplemental material is available online at www.genome.org.]

Article published online before print. Article and publication date are at http://www.genome.org/cgi/doi/10.1101/gr.070169.107

CiteULike    Connotea    Del.icio.us    Digg    Reddit    Technorati    What's this?