Work

Revealing sequence variation patterns in rice with machine learning methods

Year: 2008

Type: article

Abstract: Based on experience gained in polymorphism studies for Arabidopsis thaliana [ 1 ] we developed a method for identifying single nucleotide polymorphisms (SNPs) from the array data using Support Vector Machines (SVMs). In a two-layered approach we trained SVMs to discriminate between SNP and non-SNP positions using information from each cultivar and, in a second step, across all cultivars. Wherever several SNPs or deletion/insertion polymorphisms occur in close vicinity, the hybridisation is suppressed and SNP calling in these regions becomes infeasible. We therefore adapted a machine learning method for sequence segmentation [ 2 , 3 ] to predict highly polymorphic regions in O. sativa (cf. Figure 1 ). These regions can then be analysed in more detail using alternative experimental techniques. Log 2 intensities for the maximally hybridising oligonucleotide at each tiled position are shown for the reference and a target cultivar together with the known and predicted polymorphisms for the target cultivar (based on data from the Arabidopsis project [ 1 , 2 ]). Most of the 21 known polymorphisms in the target cultivar could not be predicted with the SNP calling methods. For training and evaluation we compiled a set of reference polymorphisms obtained by dideoxy sequencing of more than 3,500 fragments from the 20 cultivars. Across all cultivars, we discovered 1,349,341 SNPs with the machine learning (ML) method at 316,373 non-redundant positions. In comparison to a model based (MB) SNP calling approach implemented by Perlegen Sciences, Inc. [ 4 ], the ML method was found to be much more sensitive by recovering 20.9% of all known SNPs at a precision of 91.7%, compared to 14.4% and 90.9%, respectively, for the MB approach (cf. Figure 2A ). The intersection of MB and ML predictions contained 761,606 SNPs predictions at 159,879 non-repetitive positions constituting a set of markedly higher quality with a precision of 97.1%. A . Comparison of the proposed SNP prediction (ML) method based on array intensities and additional information with a previously proposed one (MB). B . Accuracy of the polymorphic region predictions using the machine learning based segmentation algorithm. In addition to SNP predictions, our polymorphic region predictor discovered a substantial additional proportion of polymorphism regions, resulting in between ~65,000 and ~203,000 polymorphic regions per cultivar (cf. Figure 2B ). We identified hundreds of thousands polymorphisms on a genome-wide scale, providing the first whole genome set of polymorphisms for the world's most important crop plant. This polymorphism data represents a valuable resource for further functional studies and modern breeding of rice. Based on the SNP data, high-density genotyping arrays will be designed to investigate genomic variation in many more rice cultivars. The PR predictions will e.g. be helpful to constrain primer design to conserved regions and thus increase PCR success rates.

Source: BMC Bioinformatics

DOI: 10.1186/1471-2105-9-s10-o8