Genome-wide mapping with biallelic markers in Arabidopsis thaliana

EST-derived single nucleotide polymorphism markers for assembling genetic and physical maps of the barley genome

In a panel of seven genotypes, 437 expressed sequence tag (EST)-derived DNA fragments were sequenced. Single nucleotide polymorphisms (SNPs) that were polymorphic between the parents of three mapping populations were mapped by heteroduplex analysis and a genome-wide consensus map comprising 216 EST-derived SNPs and 4 InDel (insertion/deletion) markers was constructed. The average frequency of SNPs amounted to 1/130 bp and 1/107.8 bp for a set of randomly selected and a set of mapped ESTs, respectively. The calculated nucleotide diversities (pi) ranged from 0 to 40.0 x 10(-3) (average 3.1 x 10(-3)) and 0.52 x 10(-3) to 39.51 x 10(-3) (average 4.37 x 10(-3)) for random and mapped ESTs, respectively. The polymorphism information content value for mapped SNPs ranged from 0.24 to 0.50 with an average of 0.34. As expected, combination of SNPs present in an amplicon (haplotype) exhibited a higher information content ranging from 0.24 to 0.85 with an average of 0.50. Cleaved amplified polymorphic sequence assays (including InDels) were designed for a total of 87 (39.5%) SNP markers. The high abundance of SNPs in the barley genome provides avenues for the systematic development of saturated genetic maps and their integration with physical maps.

Functional genomics in chickpea: an emerging frontier for molecular-assisted breeding

Chickpea is a valuable and important agricultural crop, but yield potential is limited by a series of biotic and abiotic stresses, including Ascochyta blight, Fusarium wilt, drought, cold and salinity. To accelerate molecular breeding efforts for the discovery and introgression of stress tolerance genes into cultivated chickpea, functional genomics approaches are rapidly growing. Recently a series of genetic tools for chickpea have become available that have allowed high-powered functional genomics studies to proceed, including a dense genetic map, large insert genome libraries, expressed sequence tag libraries, microarrays, serial analysis of gene expression, transgenics and reverse genetics. This review summarises the development of these genomic tools and the achievements made in initial and emerging functional genomics studies. Much of the initial research focused on Ascochyta blight resistance, and a resistance model has been synthesised based on the results of various studies. Use of the rich comparative genomics resources from the model legumes Medicago truncatula and Lotus japonicus is also discussed. Finally, perspectives on the future directions for chickpea functional genomics, with the goal of developing elite chickpea cultivars, are discussed.

Characterization of seed-specific benzoyloxyglucosinolate mutations in Arabidopsis thaliana

Glucosinolates are secondary metabolites involved in pathogen and insect defense of cruciferous plants. Although seeds and vegetative tissue often have very different glucosinolate profiles, few genetic factors that determine seed glucosinolate accumulation have been identified. An HPLC-based screen of 5500 mutagenized Arabidopsis thaliana lines produced 33 glucosinolate mutants, of which 21 have seed-specific changes. Five of these mutant lines, representing three genetic loci, are compromised in the biosynthesis of benzoyloxyglucosinolates, which are only found in seeds and young seedlings of A. thaliana. Genetic mapping and analysis of T-DNA insertions in candidate genes identified BZO1 (At1g65880), which encodes an enzyme with benzoyl-CoA ligase activity, as being required for the accumulation of benzoyloxyglucosinolates. Long-chain aliphatic glucosinolates are elevated in bzo1 mutants, suggesting substrate competition for the common short-chain aliphatic glucosinolate precursors. Whereas bzo1 mutations have seed-specific effects on benzoyloxyglucosinolate accumulation, the relative abundance of 3-benzoyloxypropyl- and 4-benzoyloxybutylglucosinolates depends on the maternal genotype.

Transcriptome profiling, sequence characterization, and SNP-based chromosomal assignment of the EXPANSIN genes in cotton

The knowledge of biological significance associated with DNA markers is very limited in cotton. SNPs are potential functional marker to tag genes of biological importance. Plant expansins are a group of extracellular proteins that directly modify the mechanical properties of cell walls, enable turgor-driven cell extension, and likely affect length and quality of cotton fibers. Here, we report the expression profiles of EXPANSIN transcripts during fiber elongation and the discovery of SNP markers, assess the SNP characteristics, and localize six EXPANSIN A genes to chromosomes. Transcriptome profiling of cotton fiber oligonucleotide microarrays revealed that seven EXPANSIN transcripts were differentially expressed when there was parallel polar elongation during morphogenesis at early stage of fiber development, suggesting that major and minor isoforms perform discrete functions during polar elongation and lateral expansion. Ancestral and homoeologous relationships of the six EXPANSIN A genes were revealed by phylogenetic grouping and comparison to extant A- and D-genome relatives of contemporary AD-genome cottons. The average rate of SNP per nucleotide was 2.35% (one SNP per 43 bp), with 1.74 and 3.99% occurring in coding and noncoding regions, respectively, in the selected genotypes. An unequal evolutionary rate of the EXPANSIN A genes at the subgenome level of tetraploid cotton was recorded. Chromosomal locations for each of the six EXPANSIN A genes were established by gene-specific SNP markers. Results revealed a strategy for discovering SNP markers in a polyploidy species like cotton. These markers could be useful to associate candidate genes with the complex fiber traits in MAS.

Single nucleotide polymorphisms in rye (Secale cereale L.): discovery, frequency, and applications for genome mapping and diversity studies

To elucidate the potential of single nucleotide polymorphism (SNP) markers in rye, a set of 48 barley EST (expressed sequence tag) primer pairs was employed to amplify from DNA prepared from five rye inbred lines. A total of 96 SNPs and 26 indels (insertion-deletions) were defined from the sequences of 14 of the resulting amplicons, giving an estimated frequency of 1 SNP per 58 bp and 1 indel per 214 bp in the rye transcriptome. A mean of 3.4 haplotypes per marker with a mean expected heterozygosity of 0.66 were observed. The nucleotide diversity index (pi) was estimated to be in the range 0.0059-0.0530. To improve assay cost-effectiveness, 12 of the 14 SNPs were converted to a cleaved amplified polymorphic sequence (CAPS) format. The resulting 12 SNP loci mapped to chromosomes 1R, 3R, 4R, 5R, 6R, and 7R, at locations consistent with their known map positions in barley. SNP genotypic data were compared with genomic simple sequence repeat (SSR) and EST-derived SSR genotypic data collected from the same templates. This showed a broad equivalence with respect to genetic diversity between these different data types.

Mutations in an auxin receptor homolog AFB5 and in SGT1b confer resistance to synthetic picolinate auxins and not to 2,4-dichlorophenoxyacetic acid or indole-3-acetic acid in Arabidopsis

Although a wide range of structurally diverse small molecules can act as auxins, it is unclear whether all of these compounds act via the same mechanisms that have been characterized for 2,4-dichlorophenoxyacetic acid (2,4-D) and indole-3-acetic acid (IAA). To address this question, we used a novel member of the picolinate class of synthetic auxins that is structurally distinct from 2,4-D to screen for Arabidopsis (Arabidopsis thaliana) mutants that show chemically selective auxin resistance. We identified seven alleles at two distinct genetic loci that conferred significant resistance to picolinate auxins such as picloram, yet had minimal cross-resistance to 2,4-D or IAA. Double mutants had the same level and selectivity of resistance as single mutants. The sites of the mutations were identified by positional mapping as At4g11260 and At5g49980. At5g49980 is previously uncharacterized and encodes auxin signaling F-box protein 5, one of five homologs of TIR1 in the Arabidopsis genome. TIR1 is the recognition component of the Skp1-cullin-F-box complex associated with the ubiquitin-proteasome pathway involved in auxin signaling and has recently been shown to be a receptor for IAA and 2,4-D. At4g11260 encodes the tetratricopeptide protein SGT1b that has also been associated with Skp1-cullin-F-box-mediated ubiquitination in auxin signaling and other pathways. Complementation of mutant lines with their corresponding wild-type genes restored picolinate auxin sensitivity. These results show that chemical specificity in auxin signaling can be conferred by upstream components of the auxin response pathway. They also demonstrate the utility of genetic screens using structurally diverse chemistries to uncover novel pathway components.

A high-resolution map of Arabidopsis recombinant inbred lines by whole-genome exon array hybridization

Recombinant populations were the basis for Mendel's first genetic experiments and continue to be key to the study of genes, heredity, and genetic variation today. Genotyping several hundred thousand loci in a single assay by hybridizing genomic DNA to oligonucleotide arrays provides a powerful technique to improve precision linkage mapping. The genotypes of two accessions of Arabidopsis were compared by using a 400,000 feature exon-specific oligonucleotide array. Around 16,000 single feature polymorphisms (SFPs) were detected in ~8,000 of the ~26,000 genes represented on the array. Allelic variation at these loci was measured in a recombinant inbred line population, which defined the location of 815 recombination breakpoints. The genetic linkage map had a total length of 422.5 cM, with 676 informative SFP markers representing intervals of ~0.6 cM. One hundred fifteen single gene intervals were identified. Recombination rate, SFP distribution, and segregation in this population are not uniform. Many genomic regions show a clustering of recombination events including significant hot spots. The precise haplotype structure of the recombinant population was defined with unprecedented accuracy and resolution. The resulting linkage map allows further refinement of the hundreds of quantitative trait loci identified in this well-studied population. Highly variable recombination rates along each chromosome and extensive segregation distortion were observed in the population.

A goal of many genetic studies is to discover the underlying genetic condition (the genotype) of a specific physical manifestation in an organism (the phenotype), such as diabetes in humans or leaf rust in cultivated wheat. A limitation to making such discoveries is the ability to resolve genotype. Gene arrays carry representations of the genome, called features, at high-density on a surface the size of a thumbnail. In this study, microarrays designed to measure gene expression were used to detect DNA sequence polymorphisms. DNA from two different Arabidopsis strains was hybridized to arrays representing nearly the entire coding region of the genome. Differences in hybridization intensity indicated differences in DNA sequence. The sequence differences, termed single feature polymorphisms, were then assayed in a population of 100 plants derived through inbreeding the progeny from the two parental strains. The precise location of the genetic recombination breakpoints was defined for each line. As a result, Singer et al. were able to generate one of the first very high-resolution genotyping data sets in a multicellular organism that allowed the construction of a high-resolution genetic map of Arabidopsis. This map will greatly facilitate attempts to make definitive associations between genotypes and phenotypes.

Gene-associated single nucleotide polymorphism discovery in perennial ryegrass (Lolium perenne L.)

Molecular genetic marker development in perennial ryegrass has largely been dependent on anonymous sequence variation. The availability of a large-scale EST resource permits the development of functionally-associated genetic markers based on SNP variation in candidate genes. Genic SNP loci and associated haplotypes are suitable for implementation in molecular breeding of outbreeding forage species. Strategies for in vitro SNP discovery through amplicon cloning and sequencing have been designed and implemented. Putative SNPs were identified within and between the parents of the F(1)(NA(6) x AU(6)) genetic mapping family and were validated among progeny individuals. Proof-of-concept for the process was obtained using the drought tolerance-associated LpASRa2 gene. SNP haplotype structures were determined and correlated with predicted amino acid changes. Gene-length LD was evaluated across diverse germplasm collections. A survey of SNP variation across 100 candidate genes revealed a high frequency of SNP incidence (c. 1 per 54 bp), with similar proportions in exons and introns. A proportion (c. 50%) of the validated genic SNPs were assigned to the F(1)(NA(6) x AU(6)) genetic map, showing high levels of coincidence with previously mapped RFLP loci. The perennial ryegrass SNP resource will enable genetic map integration, detailed LD studies and selection of superior allele content during varietal development.

New Arabidopsis recombinant inbred line populations genotyped using SNPWave and their use for mapping flowering-time quantitative trait loci

The SNPWave marker system, based on SNPs between the reference accessions Colombia-0 and Landsberg erecta (Ler), was used to distinguish a set of 92 Arabidopsis accessions from various parts of the world. In addition, we used these markers to genotype three new recombinant inbred line populations for Arabidopsis, having Ler as a common parent that was crossed with the accessions Antwerp-1, Kashmir-2, and Kondara. The benefit of using multiple populations that contain many similar markers and the fact that all markers are linked to the physical map of Arabidopsis facilitates the quantitative comparison of maps. Flowering-time variation was analyzed in the three recombinant inbred line populations. Per population, four to eight quantitative trait loci (QTL) were detected. The comparison of the QTL positions related to the physical map allowed the estimate of 12 different QTL segregating for flowering time for which Ler has an allele different from one, two, or three of the other accessions.

Frequency of the canine leucocyte adhesion deficiency (CLAD) mutation among Irish red setters in Germany

Canine leucocyte adhesion deficiency (CLAD) is an autosomal recessive hereditary disease occurring among Irish red setters. The genetic defect causative for this disorder was recently identified as a missense mutation in the ITGB2 gene. Irish red setters with one copy of the mutant gene appear normal, while dogs with two copies of the mutant gene manifest the disease. The present report describes the analysis of the single nucleotide polymorphism in 289 Irish red setters by DNA sequencing. The frequency of CLAD carriers in Germany is 11%.

Genome-wide single-nucleotide polymorphism map for Candida albicans

Single-nucleotide polymorphisms (SNPs) are essential tools for studying a variety of organismal properties and processes, such as recombination, chromosomal dynamics, and genome rearrangement. This paper describes the development of a genome-wide SNP map for Candida albicans to study mitotic recombination and chromosome loss. C. albicans is a diploid yeast which propagates primarily by clonal mitotic division. It is the leading fungal pathogen that causes infections in humans, ranging from mild superficial lesions in healthy individuals to severe, life-threatening diseases in patients with suppressed immune systems. The SNP map contains 150 marker sequences comprising 561 SNPs and 9 insertions-deletions. Of the 561 SNPs, 437 were transition events while 126 were transversion events, yielding a transition-to-transversion ratio of 3:1, as expected for a neutral accumulation of mutations. The average SNP frequency for our data set was 1 SNP per 83 bp. The map has one marker placed every 111 kb, on average, across the 16-Mb genome. For marker sequences located partially or completely within coding regions, most contained one or more nonsynonymous substitutions. Using the SNP markers, we identified a loss of heterozygosity over large chromosomal fragments in strains of C. albicans that are frequently used for gene manipulation experiments. The SNP map will be useful for understanding the role of heterozygosity and genome rearrangement in the response of C. albicans to host environments.

Identification and mapping of SNPs from ESTs in sunflower

More than 67,000 expressed sequence tags (ESTs) have recently been generated for sunflower (Helianthus), including 44,000 from cultivated confectionery (RHA280) and oilseed (RHA801) lines of Helianthus annuus and 23,000 from drought- and salt-tolerant wild sunflowers, H. argophyllus and H. paradoxus, respectively. To create a transcript map for sunflower, we identified 605 ESTs that displayed small insertion-deletion polymorphism (SNP) variation in silico, had apparent tissue-specific expression patterns, and/or were ESTs with candidate functions in traits such as development, cell transport, metabolism, plant defense, and tolerance to abiotic stress. Primer pairs for 535 of the loci were designed from the ESTs and screened for polymorphism in recombinant inbred lines derived from a cross between the same cultivars (RHA280 x RHA801) employed for sequencing. In total, 273 of the loci amplified polymorphic products, of which 243 mapped to the 17 linkage groups previously identified for sunflower. Comparisons with previously mapped QTL revealed some cases where ESTs with putatively related functions mapped near QTLs identified in other crosses for salt tolerance and for domestication traits such as stem diameter, shattering, flowering time, and achene size.

Microarray methods in Drosophila neurobiology

In the relatively short period since their development, DNA microarrays have been used increasingly in the study of genetic and cellular processes, thereby offering a genome-wide approach to gene expression studies. With the advent of genome sequencing programs for organisms from yeast to man, the number of organisms which now have ready-made commercial arrays continues to increase. Here, the principle of DNA microarrays is introduced, with particular attention being given to the role of this technology in studies of the nervous system of the fruitfly Drosophila melanogaster. The importance of experimental design and sample preparation, in line with minimum information about microarray experiment (MIAME) compliance, is emphasised. The technical platforms available to the Drosophila neurobiologist have been illustrated and a brief number of data analysis tools that are readily available reviewed.

Discovery of SNPs in soybean genotypes frequently used as the parents of mapping populations in the United States and Korea

Single nucleotide polymorphisms (SNPs) including insertion/deletions (indels) serve as useful and informative genetic markers. The availability of high-throughput and inexpensive SNP typing systems has increased interest in the development of SNP markers. After fragments of genes were amplified with primers derived from 110 soybean GenBank ESTs, sequencing data of PCR products from 15 soybean genotypes from Korea and the United States were analyzed by SeqScape software to find SNPs. Among 35 gene fragments with at least one SNP among the 15 genotypes, SNPs occurred at a frequency of 1 per 2,038 bp in 16,302 bp of coding sequence and 1 per 191 bp in 16,960 bp of noncoding regions. This corresponds to a nucleotide diversity (theta) of 0.00017 and 0.00186, respectively. Of the 97 SNPs discovered, 78 or 80.4% were present in the six North American soybean mapping parents. The addition of "Hwaeomputkong," which originated from Japan, increased the number to 92, or 94.8% of the total number of SNPs present among the 15 genotypes. Thus, Hwaeomputkong and the six North American mapping parents provide a diverse set of soybean genotypes that can be successfully used for SNP discovery in coding DNA and closely associated introns and untranslated regions.

Demonstration of loss of heterozygosity by single-nucleotide polymorphism microarray analysis and alterations in strain morphology in Candida albicans strains during infection

Candida albicans is a diploid yeast with a predominantly clonal mode of reproduction, and no complete sexual cycle is known. As a commensal organism, it inhabits a variety of niches in humans. It becomes an opportunistic pathogen in immunocompromised patients and can cause both superficial and disseminated infections. It has been demonstrated that genome rearrangement and genetic variation in isolates of C. albicans are quite common. One possible mechanism for generating genome-level variation among individuals of this primarily clonal fungus is mutation and mitotic recombination leading to loss of heterozygosity (LOH). Taking advantage of a recently published genome-wide single-nucleotide polymorphism (SNP) map (A. Forche, P. T. Magee, B. B. Magee, and G. May, Eukaryot. Cell 3:705-714, 2004), an SNP microarray was developed for 23 SNP loci residing on chromosomes 5, 6, and 7. It was used to examine 21 strains previously shown to have undergone mitotic recombination at the GAL1 locus on chromosome 1 during infection in mice. In addition, karyotypes and morphological properties of these strains were evaluated. Our results show that during in vivo passaging, LOH events occur at observable frequencies, that such mitotic recombination events occur independently in different loci across the genome, and that changes in karyotypes and alterations of phenotypic characteristics can be observed alone, in combination, or together with LOH.

Single-nucleotide polymorphisms detected in expressed sequence tags of melon (Cucumis melo L.)

A search was performed for single-nucleotide polymorphisms (SNP) and short insertions-deletions (indels) in 34 melon (Cucumis melo L.) expressed sequence tag (EST) fragments between two distantly related melon genotypes, a group Inodorus 'Piel de sapo' market class breeding line T111 and the Korean accession PI 161375. In total, we studied 15 kb of melon sequence. The average frequency of SNPs between the two genotypes was one every 441 bp. One indel was also found every 1666 bp. Seventy-five percent of the polymorphisms were located in introns and the 3'untranslated regions. On average, there were 1.26 SNPs plus indels per amplicon. We explored three different SNP detection systems to position five of the SNPs in a melon genetic map. Three of the SNPs were mapped using cleaved amplified polymorphic sequence (CAPS) markers, one SNP was mapped using the single primer extension reaction with fluorescent-labelled dideoxynucleotides, and one indel was mapped using polyacrilamide gel electrophoresis separation. The discovery of SNPs based on ESTs and a suitable system for SNP detection has broad potential utility in melon genome mapping.

Using cDNA and genomic sequences as tools to develop SNP strategies in cassava (Manihot esculenta Crantz)

Single nucleotide polymorphisms (SNP) are the most abundant type of DNA polymorphism found in animal and plant genomes. They provide an important new source of molecular markers that are useful in genetic mapping, map-based positional cloning, quantitative trait locus mapping and the assessment of genetic distances between individuals. Very little is known on the frequency of SNPs in cassava. We have exploited the recently-developed collection of cassava expressed sequence tags (ESTs) to detect SNPs in the five cultivars of cassava used to generate the sequences. The frequency of intra-cultivar and inter-cultivar SNPs after analysis of 111 contigs was one polymorphism per 905 and one per 1,032 bp, respectively; totaling 1 each 509 bp. We have obtained further information on the frequency of SNPs in six cassava cultivars by analysis of 33 amplicons obtained from 3' EST and BAC end sequences. Overall, about 11 kb of DNA sequence was obtained for each cultivar. A total of 186 SNPs (136 and 50 from ESTs and BAC ends, respectively) were identified. Among these, 146 were intra-cultivar polymorphisms, while 80 were inter-cultivar polymorphisms. Thus the total frequency of SNPs was one per 62 bp. This information will help to develop new strategies for EST mapping as well as their association with phenotypic characteristics.

Reduced variation on the chicken Z chromosome

Understanding the population genetic factors that shape genome variability is pivotal to the design and interpretation of studies using large-scale polymorphism data. We analyzed patterns of polymorphism and divergence at Z-linked and autosomal loci in the domestic chicken (Gallus gallus) to study the influence of mutation, effective population size, selection, and demography on levels of genetic diversity. A total of 14 autosomal introns (8316 bp) and 13 Z-linked introns (6856 bp) were sequenced in 50 chicken chromosomes from 10 highly divergent breeds. Genetic variation was significantly lower at Z-linked than at autosomal loci, with one segregating site every 39 bp at autosomal loci (theta(W) = 5.8 +/- 0.8 x 10(-3)) and one every 156 bp on the Z chromosome (theta(W) = 1.4 +/- 0.4 x 10(-3)). This difference may in part be due to a low male effective population size arising from skewed reproductive success among males, evident both in the wild ancestor-the red jungle fowl-and in poultry breeding. However, this effect cannot entirely explain the observed three- to fourfold reduction in Z chromosome diversity. Selection, in particular selective sweeps, may therefore have had an impact on reducing variation on the Z chromosome, a hypothesis supported by the observation of heterogeneity in diversity levels among loci on the Z chromosome and the lower recombination rate on Z than on autosomes. Selection on sex-linked genes may be particularly important in organisms with female heterogamety since the heritability of sex-linked sexually antagonistic alleles advantageous to males is improved when fathers pass a Z chromosome to their sons.

Ligation detection reaction-TaqMan procedure for single nucleotide polymorphism detection on genomic DNA

In this article, we describe a genotyping approach applicable to both individual and multiplexed single nucleotide polymorphism (SNP) analysis, based on a ligation detection reaction (LDR) performed directly on genomic DNA. During the ligation, the biallelic state of the SNP locus is converted into a bimarker state of ligated detector oligonucleotides. The state of the markers is then determined by a 5'-nuclease assay (TaqMan) with universal fluorescent probes. The LDR-TaqMan method was successfully applied for the genotyping of 30 SNP loci of Arabidopsis thaliana. The technology is cost-effective, needs no locus-specific optimization, requires minimal manipulations, and has very good potential for automation.

Influence of microarrays experiments missing values on the stability of gene groups by hierarchical clustering

Background

Microarray technologies produced large amount of data. The hierarchical clustering is commonly used to identify clusters of co-expressed genes. However, microarray datasets often contain missing values (MVs) representing a major drawback for the use of the clustering methods. Usually the MVs are not treated, or replaced by zero or estimated by the k-Nearest Neighbor (kNN) approach. The topic of the paper is to study the stability of gene clusters, defined by various hierarchical clustering algorithms, of microarrays experiments including or not MVs.

Results

In this study, we show that the MVs have important effects on the stability of the gene clusters. Moreover, the magnitude of the gene misallocations is depending on the aggregation algorithm. The most appropriate aggregation methods (e.g. complete-linkage and Ward) are highly sensitive to MVs, and surprisingly, for a very tiny proportion of MVs (e.g. 1%). In most of the case, the MVs must be replaced by expected values. The MVs replacement by the kNN approach clearly improves the identification of co-expressed gene clusters. Nevertheless, we observe that kNN approach is less suitable for the extreme values of gene expression.

Conclusion

The presence of MVs (even at a low rate) is a major factor of gene cluster instability. In addition, the impact depends on the hierarchical clustering algorithm used. Some methods should be used carefully. Nevertheless, the kNN approach constitutes one efficient method for restoring the missing expression gene values, with a low error level. Our study highlights the need of statistical treatments in microarray data to avoid misinterpretation.

Mining single nucleotide polymorphisms from EST data of silkworm, Bombyx mori, inbred strain Dazao

We made use of 81,635 expressed sequence tags (ESTs) derived from 12 different cDNA libraries of the silkworm, Bombyx mori, inbred strain Dazao (P50), to identify high-quality candidate single nucleotide polymorphisms (SNPs). By PHRAP assembling, 12,980 contigs containing 11,537 contigs assembled by more than one read were obtained, and 101 candidate SNPs and 27 single base insertions/deletions were identified from 117 contigs assembled from 1576 high-quality reads base-called with PHRED and screened on the basis of the neighborhood quality standard (NQS). Simultaneously, we also predicted 40 SNPs in coding regions (cSNPs), of which 26 were predicted to lead to amino acid non-synonymous variations and 14 synonymous substitutions. Also, the 1.66:1 ratio of transition/transversion is different from that of other insects. As the first SNP analysis of a Lepidoptera, B. mori, the single nucleotide polymorphic density is estimated to be 1.3 x 10(-3) by sequence diversity. This analysis shows that expressed sequences from multiple libraries may provide an abundant source of comparative reads to mine for cSNPs from the silkworm genome.

SNPWave: a flexible multiplexed SNP genotyping technology

Scalable multiplexed amplification technologies are needed for cost-effective large-scale genotyping of genetic markers such as single nucleotide polymorphisms (SNPs). We present SNPWave, a novel SNP genotyping technology to detect various subsets of sequences in a flexible fashion in a fixed detection format. SNPWave is based on highly multiplexed ligation, followed by amplification of up to 20 ligated probes in a single PCR. Depending on the multiplexing level of the ligation reaction, the latter employs selective amplification using the amplified fragment length polymorphism (AFLP) technology. Detection of SNPWave reaction products is based on size separation on a sequencing instrument with multiple fluorescence labels and short run times. The SNPWave technique is illustrated by a 100-plex genotyping assay for Arabidopsis, a 40-plex assay for tomato and a 10-plex assay for Caenorhabditis elegans, detected on the MegaBACE 1000 capillary sequencer.

An AFLP-based genome-wide mapping strategy

To efficiently determine the chromosomal location of phenotypic mutants, we designed a genome-wide mapping strategy that can be used in any crop for which a dense AFLP (Amplified Fragment Length Polymorphism) map is available or can be made. The AFLP technique is particularly suitable to initiate map-based cloning projects because it detects many markers per reaction. First a standard set of AFLP primer combinations that results in a framework of AFLP markers well dispersed over the genome is selected. These primer combinations are applied to a limited number of mutant individuals from a segregating population to register linkage and non-linkage of the AFLP markers to the gene-of-interest. Further delineation of the area of interest is accomplished by analyzing the remaining recombinants and additional mutant individuals with AFLP markers that lie within the identified region. We illustrate the usefulness of the method by mapping three rotunda ( ron) leaf-form mutant loci of Arabidopsis thaliana and show that in the initial phase of map-based cloning projects a 400-600 kb interval can be identified for the average mutant locus within a few weeks. Once such an area is identified and before initiating the more time-consuming fine-mapping procedure, it is essential to examine publicly available databases for candidate genes and known mutants in the identified region. The 390-kb interval on chromosome 4 that harbors the ron2 mutation, also carries a known flower mutant, leunig ( lug); upon crossing, the two mutants appeared to be allelic. When no such candidates are found, the mapping procedure should be continued. We present a strategy to efficiently select recombinants that can be used for fine mapping.

Discrimination of homoeologous gene expression in hexaploid wheat by SNP analysis of contigs grouped from a large number of expressed sequence tags

Single-nucleotide polymorphisms (SNPs) are useful markers for gene diagnosis and mapping of genes on chromosomes. However, polyploidy, which is characteristic of the evolution of higher plants, complicates the analysis of SNPs in the duplicated genes. We have developed a new method for SNP analysis in hexaploid wheat. First, we classified a large number of expressed sequence tags (ESTs) from wheat in silico. Those grouped into contigs were anticipated to correspond to transcripts from homoeologous loci. We then selected relatively abundant ESTs, and assigned these contigs to each of the homoeologous chromosomes using a nullisomic/tetrasomic series of Chinese Spring wheat strains in combination with pyrosequencing. The ninety genes assigned were almost evenly distributed into seven homologous chromosomes. We then created a virtual display of the relative expression of these genes. Expression patterns of genes from the three genomes in hexaploid wheat were classified into two major groups: (1) genes almost equally expressed from all three genomes; and (2) genes expressed with a significant preference, which changed from tissue to tissue, from certain genomes. In 11 cases, one of the three genes in the allopolyploid was found to be silenced. No preference for gene-silencing in particular genomes or chromosomes was observed, suggesting that gene-silencing occurred after polyploidization, and at the gene level, not at the chromosome or genome level. Thus, the use of this SNP method to distinguish the expression profiles of three homoeologous genes may help to elucidate the molecular basis of heterosis in polyploid plants.

First-generation SNP/InDel markers tagging loci for pathogen resistance in the potato genome

A panel of 17 tetraploid and 11 diploid potato genotypes was screened by comparative sequence analysis of polymerase chain reaction (PCR) products for single nucleotide polymorphisms (SNPs) and insertion-deletion polymorphisms (InDels), in regions of the potato genome where genes for qualitative and/or quantitative resistance to different pathogens have been localized. Most SNP and InDel markers were derived from bacterial artificial chromosome (BAC) insertions that contain sequences similar to the family of plant genes for pathogen resistance having nucleotide-binding-site and leucine-rich-repeat domains (NBS-LRR-type genes). Forty-four such NBS-LRR-type genes containing BAC-insertions were mapped to 14 loci, which tag most known resistance quantitative trait loci (QTL) in potato. Resistance QTL not linked to known resistance-gene-like (RGL) sequences were tagged with other markers. In total, 78 genomic DNA fragments with an overall length of 31 kb were comparatively sequenced in the panel of 28 genotypes. 1498 SNPs and 127 InDels were identified, which corresponded, on average, to one SNP every 21 base pairs and one InDel every 243 base pairs. The nucleotide diversity of the tetraploid genotypes (pi = 0.72 x 10(-3)) was lower when compared with diploid genotypes (pi = 2.31 x 10(-3)). RGL sequences showed higher nucleotide diversity when compared with other sequences, suggesting evolution by divergent selection. Information on sequences, sequence similarities, SNPs and InDels is provided in a database that can be queried via the Internet.

Snipping polymorphisms from large EST collections in barley (Hordeum vulgare L.)

The public EST (expressed sequence tag) databases represent an enormous but heterogeneous repository of sequences, including many from a broad selection of plant species and a wide range of distinct varieties. The significant redundancy within large EST collections makes them an attractive resource for rapid pre-selection of candidate sequence polymorphisms. Here we present a strategy that allows rapid identification of candidate SNPs in barley (Hordeum vulgare L.) using publicly available EST databases. Analysis of 271,630 EST sequences from different cDNA libraries, representing 23 different barley varieties, resulted in the generation of 56,302 tentative consensus sequences. In all, 8171 of these unigene sequences are members of clusters with six or more ESTs. By applying a novel SNP detection algorithm (SNiPpER) to these sequences, we identified 3069 candidate inter-varietal SNPs. In order to verify these candidate SNPs, we selected a small subset of 63 present in 36 ESTs. Of the 63 SNPs selected, we were able to validate 54 (86%) using a direct sequencing approach. For further verification, 28 ESTs were mapped to distinct loci within the barley genome. The polymorphism information content (PIC) and nucleotide diversity (pi) values of the SNPs identified by the SNiPpER algorithm are significantly higher than those that were obtained by random sequencing. This demonstrates the efficiency of our strategy for SNP identification and the cost-efficient development of EST-based SNP-markers.

Forward genetics and map-based cloning approaches

Whereas reverse genetics strategies seek to identify and select mutations in a known sequence, forward genetics requires the cloning of sequences underlying a particular mutant phenotype. Map-based cloning is tedious, hampering the quick identification of candidate genes. With the unprecedented progress in the sequencing of whole genomes, and perhaps even more with the development of saturating marker technologies, map-based cloning can now be performed so efficiently that, at least for some plant model systems, it has become feasible to identify some candidate genes within a few months. This, in turn, will boost the use of forward genetics approaches, as applied (for example) to isolating genes involved in natural variation and genes causing phenotypic mutations as derived from (second-site) mutagenesis screens.

Establishment of a high-efficiency SNP-based framework marker set for Arabidopsis

The major goal of this project was the establishment of a tool for rapid mapping of new mutations and genotyping in Arabidopsis consisting of at least 100 evenly spaced framework markers. We assembled a single nucleotide polymorphism (SNP)-based marker set consisting of 112 polymorphic sites with average spacing of 1.15 Mbp derived from an SNP database that we recently developed. This information was used to set up efficient SNP detection reactions based on multiplexed primer extension assays. The 112 Columbia (Col-0)/C24 framework markers were used to assemble 18 multiplexed SNaPshot assays with which up to eight separate loci can be genotyped in a single-tube/single-capillary format. In addition, for 110 framework markers matrix-assisted laser desorption/ionization time of flight (MALDI-ToF) assays have been established for high throughput analyses. We demonstrated the usefulness and the robustness of both procedures of this tool by genotyping 48 BC3F1 individuals created between the accessions Col-0 and C24. Subsets of 10-62 of the established markers discriminate between various combinations of the accessions Col-0, C24, Landsberg erecta (Ler), Cape Verdi Islands (Cvi) and Niederzenz (Nd). Using a subset of 17 evenly distributed and established SNP markers that are also polymorphic between Ler and Col-0, we were able to rapidly map a mutant gene (tbr1) to an interval of 2.3 Mbp in an Ler (tbr1) x Col-0 cross.

Combined flash- and glow-type chemiluminescent reactions for high-throughput genotyping of biallelic polymorphisms

The difference in light-emission kinetics between the Ca(2+)-triggered bioluminescent reaction of the photoprotein aequorin (AEQ) and the alkaline phosphatase (ALP)-catalyzed chemiluminescent hydrolysis of dioxetane aryl phosphate substrates was exploited for the analysis of both alleles of biallelic polymorphisms in a single microtiter well. The genotyping of the IVS-1-110 locus of the human beta-globin gene was chosen as a model. Genomic DNA, isolated from whole blood, was first subjected to polymerase chain reaction using primers flanking the polymorphic site. A single oligonucleotide-ligation reaction employing two allele-specific probes, labeled with biotin and digoxigenin, and a common probe carrying a characteristic tail was then performed. The ligation products were captured in a microtiter well through hybridization of the tail with an immobilized complementary oligonucleotide. The products were detected by adding a mixture of streptavidin-aequorin complex and antidigoxigenin-alkaline phosphatase conjugate. AEQ was measured first by adding Ca(2+) and integrating the signal for 3s followed by the addition of the substrate for ALP. The ratio of the luminescence signals obtained from ALP and AEQ gives the genotype of each sample. The coefficient of variation of the dual assay ranged from 7 to 11% for each allele. The reproducibility of the ALP/AEQ signal ratio was about 14%. The proposed assay allows for many samples to be screened in parallel in a single microtiter plate, for single-nucleotide polymorphisms.

Heterosis for biomass yield and related traits in five hybrids of Arabidopsis thaliana L. Heynh

Heterosis is of utmost economic importance in plant breeding. However, its underlying molecular causes are still unknown. Given the numerous advantages of Arabidopsis thaliana as a model species in plant genetics and genomics, we assessed the extent of heterosis in this species using five hybrids derived from five ecotypes. Parents, F(1) and F(2), generations in both reciprocal forms were grown in a greenhouse experiment with four replications. Mid-parent heterosis (MPH) and best-parent heterosis (BPH) averaged across hybrids were surprisingly high for biomass yield (MPH: 60.3%; BPH: 32.9%) and rosette diameter (MPH: 49.4%; BPH: 34.8%), but smaller for flowering date (MPH: 27.5%; BPH: 18.5%), seed yield (MPH: 18.9%; BPH: 1.7%), and yield components. Individual hybrids varied considerably in their MPH and BPH values for all traits, one cross displaying 140.1% MPH for biomass yield. MPH was not associated with parental genetic distance determined from molecular markers. Reciprocal effects were significant only in a few cases. With a proper choice of hybrids, our results encourage the use of Arabidopsis as a model species for investigating the molecular causes of heterosis.

Large-scale identification and analysis of genome-wide single-nucleotide polymorphisms for mapping in Arabidopsis thaliana

Genetic markers such as single nucleotide polymorphisms (SNPs) are essential tools for positional cloning, association, or quantitative trait locus mapping and the determination of genetic relationships between individuals. We identified and characterized a genome-wide set of SNP markers by generating 10,706 expressed sequence tags (ESTs) from cDNA libraries derived from 6 different accessions, and by analysis of 606 sequence tagged sites (STS) from up to 12 accessions of the model flowering plant Arabidopsis thaliana. The cDNA libraries for EST sequencing were made from individuals that were stressed by various means to enrich for transcripts from genes expressed under such conditions. SNPs discovered in these sequences may be useful markers for mapping genes involved in interactions with the biotic and abiotic environment. The STS loci are distributed randomly over the genome. By comparison with the Col-0 genome sequence, we identified a total of 8051 SNPs and 637 insertion/deletion polymorphisms (InDel). Analysis of STS-derived SNPs shows that most SNPs are rare, but that it is possible to identify intermediate frequency framework markers that can be used for genetic mapping in many different combinations of accessions. A substantial proportion of SNPs located in ORFs caused a change of the encoded amino acid. A comparison of the density of our SNP markers among accessions in both the EST and STS datasets, revealed that Cvi-0 is the most divergent accession from Col-0 among the 12 accessions studied. All of these markers are freely available via the internet.

Large-scale identification of single-feature polymorphisms in complex genomes

We have developed a high-throughput genotyping platform by hybridizing genomic DNA from Arabidopsis thaliana accessions to an RNA expression GeneChip (AtGenome1). Using newly developed analytical tools, a large number of single-feature polymorphisms (SFPs) were identified. A comparison of two accessions, the reference strain Columbia (Col) and the strain Landsberg erecta (Ler), identified nearly 4000 SFPs, which could be reliably scored at a 5% error rate. Ler sequence was used to confirm 117 of 121 SFPs and to determine the sensitivity of array hybridization. Features containing sequence repeats, as well as those from high copy genes, showed greater polymorphism rates. A linear clustering algorithm was developed to identify clusters of SFPs representing potential deletions in 111 genes at a 5% false discovery rate (FDR). Among the potential deletions were transposons, disease resistance genes, and genes involved in secondary metabolism. The applicability of this technique was demonstrated by genotyping a recombinant inbred line. Recombination break points could be clearly defined, and in one case delimited to an interval of 29 kb. We further demonstrate that array hybridization can be combined with bulk segregant analysis to quickly map mutations. The extension of these tools to organisms with complex genomes, such as Arabidopsis, will greatly increase our ability to map and clone quantitative trait loci (QTL).

Arabidopsis gene expression changes during cyst nematode parasitism revealed by statistical analyses of microarray expression profiles

With the availability of microarray technology, the expression profiles of thousands of genes can be monitored simultaneously to help determine the mechanisms of these biological processes. We conducted Affymetrix GeneChip microarray analyses of the Arabidopsis-cyst nematode interaction and employed a statistical procedure to analyze the resultant data, which allowed us to identify significant gene expression changes. Quantitative real-time RT-PCR assays were used to confirm the microarray analyses. The results of the expression profiling revealed 128 genes with altered steady-state mRNA levels following infection by the sugar beet cyst nematode (Heterodera schachtii; BCN), in contrast to only 12 genes that had altered expression following infection by the soybean cyst nematode (H. glycines; SCN). The expression of these 12 genes also changed following infection by BCN, i.e. we did not identify any genes regulated exclusively by SCN. The identification of 116 genes whose expression changes during successful cyst nematode parasitism by BCN suggests a potential involvement of these genes in the infection events starting with successful syncytium induction. Further characterization of these genes will permit the formulation of testable hypotheses to explain successful cyst nematode parasitism.

A decade of high-resolution liquid chromatography of nucleic acids on styrene-divinylbenzene copolymers

The introduction of alkylated, nonporous poly-(styrene-divinylbenzene) microparticles in 1992 enabled the subsequent development of denaturing HPLC that has emerged as the most sensitive screening method for mutations to date. Denaturing HPLC has provided unprecedented insight into human origins and prehistoric migrations, accelerated the cloning of genes involved in mono- and polygenic traits, and facilitated the mutational analysis of more than a hundred candidate genes of human disease. A significant step toward increased sample-throughput and information content was accomplished by the recent introduction of monolithic poly(styrene-divinylbenzene) capillary columns. They have enabled the construction of capillary arrays amenable to multiplex analysis of fluorescent dye-labeled nucleic acids by laser-induced fluorescence detection. Hyphenation of denaturing HPLC with electrospray ionization mass spectrometry, on the other hand, has allowed the direct elucidation of the chemical nature of DNA variation and determination of phase of multiple alleles on a chromosome.

Sequence variation and the biological function of genes: methodological and biological considerations

Single nucleotide polymorphisms (SNPs) are expected to facilitate the chromosomal mapping and eventual cloning of genetic determinants of complex quantitative phenotypes. To date, more than 2.5 million non-redundant human SNPs have been reported in the public domain, of which approximately 100000 have been validated by either independent investigators or by independent methods. Equally impressive is the myriad of methods developed for allelic discrimination. Nevertheless, reports of successful applications of SNPs to genome-wide linkage analysis of both mono- and polygenic traits are rare and limited to a few model organisms, that provide affordable platforms to test both novel methodological and biological concepts at a whole-genome scale under conditions that can be reasonably controlled. Progress in the analysis of SNPs needs to be complemented by methods that allow the systematic elucidation of both primary and secondary phenotypes of genes. Importantly, observations made in one species may very well be of immediate applicability to other species including human. This is particularly true for conserved biological processes such as mitochondrial respiration and DNA repair.

Rapid mapping of zebrafish mutations with SNPs and oligonucleotide microarrays

Large-scale genetic screens in zebrafish have identified thousands of mutations in hundreds of essential genes. The genetic mapping of these mutations is necessary to link DNA sequences to the gene functions defined by mutant phenotypes. Here, we report two advances that will accelerate the mapping of zebrafish mutations: (1) The construction of a first generation single nucleotide polymorphism (SNP) map of the zebrafish genome comprising 2035 SNPs and 178 small insertions/deletions, and (2) the development of a method for mapping mutations in which hundreds of SNPs can be scored in parallel with an oligonucleotide microarray. We have demonstrated the utility of the microarray technique in crosses with haploid and diploid embryos by mapping two known mutations to their previously identified locations. We have also used this approach to localize four previously unmapped mutations. We expect that mapping with SNPs and oligonucleotide microarrays will accelerate the molecular analysis of zebrafish mutations.

SNP typing by apyrase-mediated allele-specific primer extension on DNA microarrays

This study reports the development of a microarray-based allele-specific extension method for typing of single nucleotide polymorphisms (SNPs). The use of allele-specific primers has been employed previously to identify single base variations but it is acknowledged that certain mismatches are not refractory to extension. Here we have overcome this limitation by introducing apyrase, a nucleotide-degrading enzyme, to the extension reaction. We have shown previously that DNA polymerases exhibit slower reaction kinetics when extending a mismatched primer compared with a matched primer. This kinetic difference is exploited in the apyrase-mediated allele-specific extension (AMASE) assay, allowing incorporation of nucleotides when the reaction kinetics are fast but degrading the nucleotides before extension when the reaction kinetics are slow. Here we show that five homozygous variants (14% of the total number of variants) that were incorrectly scored in the absence of apyrase were correctly typed when apyrase was included in the extension reaction. AMASE was performed in situ on the oligonucleotide microarrays using fluorescent nucleotides to type 10 SNPs and two indels in 17 individuals generating approximately 200 genotypes. Cluster analysis of these data shows three distinct clusters with clear-cut boundaries. We conclude that SNP typing on oligonucleotide microarrays by AMASE is an efficient, rapid and accurate technique for large-scale genotyping.

Arabidopsis map-based cloning in the post-genome era

Map-based cloning is an iterative approach that identifies the underlying genetic cause of a mutant phenotype. The major strength of this approach is the ability to tap into a nearly unlimited resource of natural and induced genetic variation without prior assumptions or knowledge of specific genes. One begins with an interesting mutant and allows plant biology to reveal what gene or genes are involved. Three major advances in the past 2 years have made map-based cloning in Arabidopsis fairly routine: sequencing of the Arabidopsis genome, the availability of more than 50,000 markers in the Cereon Arabidopsis Polymorphism Collection, and improvements in the methods used for detecting DNA polymorphisms. Here, we describe the Cereon Collection and show how it can be used in a generic approach to mutation mapping in Arabidopsis. We present the map-based cloning of the VTC2 gene as a specific example of this approach.

A draft sequence of the rice genome (Oryza sativa L. ssp. indica)

We have produced a draft sequence of the rice genome for the most widely cultivated subspecies in China, Oryza sativa L. ssp. indica, by whole-genome shotgun sequencing. The genome was 466 megabases in size, with an estimated 46,022 to 55,615 genes. Functional coverage in the assembled sequences was 92.0%. About 42.2% of the genome was in exact 20-nucleotide oligomer repeats, and most of the transposons were in the intergenic regions between genes. Although 80.6% of predicted Arabidopsis thaliana genes had a homolog in rice, only 49.4% of predicted rice genes had a homolog in A. thaliana. The large proportion of rice genes with no recognizable homologs is due to a gradient in the GC content of rice coding sequences.

Annotation and BAC/PAC localization of nonredundant ESTs from drought-stressed seedlings of an indica rice

To decipher the genes associated with drought stress response and to identify novel genes in rice, we utilized 1540 high-quality expressed sequence tags (ESTs) for functional annotation and mapping to rice genomic sequences. These ESTs were generated earlier by 3'-end single-pass sequencing of 2000 cDNA clones from normalized cDNA libraries constructed form drought-stressed seedlings of an indica rice. A rice UniGene set of 1025 transcripts was constructed from this collection through the BLASTN algorithm. Putative functions of 559 nonredundant ESTs were identified by BLAST similarity search against public databases. Putative functions were assigned at a stringency E value of 10(-6) in BLASTN and BLASTX algorithms. To understand the gene structure and function further, we have utilized the publicly available finished and unfinished rice BAC/PAC (BAC, bacterial artificial chromosome; PAC, P1 artificial chromosome) sequences for similarity search using the BLASTN algorithm. Further, 603 nonredundant ESTs have been mapped to BAC/PAC clones. BAC clones were assigned by a homology of above 95% identity along 90% of EST sequence length in the aligned region. In all, 700 ESTs showed rice EST hits in GenBank. Of the 325 novel ESTs, 128 were localized to BAC clones. In addition, 127 ESTs with identified putative functions but with no homology in IRGSP (International Rice Genome Sequencing Program) BAC/PAC sequences were mapped to the Chinese WGS (whole genome shotgun contigs) draft sequence of the rice genome. Functional annotation uncovered about a hundred candidate ESTs associated with abiotic stress in rice and Arabidopsis that were previously reported based on microarray analysis and other studies. This study is a major effort in identifying genes associated with drought stress response and will serve as a resource to rice geneticists and molecular biologists.

Sequence variation in genes and genomic DNA: methods for large-scale analysis

The large-scale typing of sequence variation in genes and genomic DNA presents new challenges for which it is not clear that current technologies are sufficiently sensitive, robust, or scalable. This review surveys the current platform technologies: separation-based approaches, which include mass spectrometry; homogeneous assays; and solid-phase/array-based assays. We assess techniques for discovering and typing variation on a large scale, especially that of single-nucleotide polymorphisms. The in-depth focus is the DNA chip/array platform, and some of the published large-scale studies are closely examined. The problem of large-scale amplification is addressed, and emerging technologies for present and future needs are indicated.

The extent of linkage disequilibrium in Arabidopsis thaliana

Linkage disequilibrium (LD), the nonrandom occurrence of alleles in haplotypes, has long been of interest to population geneticists. Recently, the rapidly increasing availability of genomic polymorphism data has fueled interest in LD as a tool for fine-scale mapping, in particular for human disease loci. The chromosomal extent of LD is crucial in this context, because it determines how dense a map must be for associations to be detected and, conversely, limits how finely loci may be mapped. Arabidopsis thaliana is expected to harbor unusually extensive LD because of its high degree of selfing. Several polymorphism studies have found very strong LD within individual loci, but also evidence of some recombination. Here we investigate the pattern of LD on a genomic scale and show that in global samples, LD decays within approximately 1 cM, or 250 kb. We also show that LD in local populations may be much stronger than that of global populations, presumably as a result of founder events. The combination of a relatively high level of polymorphism and extensive haplotype structure bodes well for developing a genome-wide LD map in A. thaliana.

Genetic mapping with SNP markers in Drosophila

Map-based positional cloning of Drosophila melanogaster genes is hampered by both the time-consuming, error-prone nature of traditional methods for genetic mapping and the difficulties in aligning the genetic and cytological maps with the genome sequence. The identification of sequence polymorphisms in the Drosophila genome will make it possible to map mutations directly to the genome sequence with high accuracy and resolution. Here we report the identification of 7,223 single-nucleotide polymorphisms (SNPs) and 1,392 insertions/deletions (InDels) in common laboratory strains of Drosophila. These sequence polymorphisms define a map of 787 autosomal marker loci with a resolution of 114 kb. We have established PCR product-length polymorphism (PLP) or restriction fragment-length polymorphism (RFLP) assays for 215 of these markers. We demonstrate the use of this map by delimiting two mutations to intervals of 169 kb and 307 kb, respectively. Using a local high-density SNP map, we also mapped a third mutation to a resolution of approximately 2 kb, sufficient to localize the mutation within a single gene. These methods should accelerate the rate of positional cloning in Drosophila.