Parametric alignment of Drosophila genomes

The classic algorithms of Needleman-Wunsch and Smith-Waterman find a maximum a posteriori probability alignment for a pair hidden Markov model (PHMM). To process large genomes that have undergone complex genome rearrangements, almost all existing whole genome alignment methods apply fast heuristics to divide genomes into small pieces that are suitable for Needleman-Wunsch alignment. In these alignment methods, it is standard practice to fix the parameters and to produce a single alignment for subsequent analysis by biologists. As the number of alignment programs applied on a whole genome scale continues to increase, so does the disagreement in their results. The alignments produced by different programs vary greatly, especially in non-coding regions of eukaryotic genomes where the biologically correct alignment is hard to find. Parametric alignment is one possible remedy. This methodology resolves the issue of robustness to changes in parameters by finding all optimal alignments for all possible parameters in a PHMM. Our main result is the construction of a whole genome parametric alignment of Drosophila melanogaster and Drosophila pseudoobscura. This alignment draws on existing heuristics for dividing whole genomes into small pieces for alignment, and it relies on advances we have made in computing convex polytopes that allow us to parametrically align non-coding regions using biologically realistic models. We demonstrate the utility of our parametric alignment for biological inference by showing that cis-regulatory elements are more conserved between Drosophila melanogaster and Drosophila pseudoobscura than previously thought. We also show how whole genome parametric alignment can be used to quantitatively assess the dependence of branch length estimates on alignment parameters.

A dicistronic gene pair within a cluster of "EF-hand" protein genes in the genomes of Drosophila species

Androcam is a Drosophila melanogaster calmodulin-related protein that functions specifically in the testis. We show that the Acam gene is part of a cluster of three intronless genes arranged in a head-to-tail manner. The additional genes also encode calmodulin-related proteins with testis-specific transcription. Acam and the 5'-most gene (gene1) generate monocistronic transcripts. Surprisingly, the central gene (gene2) is transcribed only as a dicistronic transcript with Acam. A similar cluster is found in D. yakuba. In D. pseudoobscura, the cluster contains four genes: two Acam-type genes downstream of a single gene related to both gene1 and gene2 and a fourth weakly related gene. Nevertheless, the D. pseudoobscura cluster also generates a dicistronic transcript from a gene pair analogous to the gene2-Acam pair. A cotranscribed gene1/2-Acam gene pair may be the founding feature of this locus. Although Acam protein is present in D. melanogaster and D. pseudoobscura testes, cognate proteins for the gene1/2-type ORFs are not detectable by immunoblotting and mass spectrometry techniques.

A Bayesian heterogeneous analysis of variance approach to inferring recent selective sweeps

The distribution of microsatellite allele sizes in populations aids in understanding the genetic diversity of species and the evolutionary history of recent selective sweeps. We propose a heterogeneous Bayesian analysis of variance model for inferring loci involved in recent selective sweeps by analyzing the distribution of allele sizes at multiple loci in multiple populations. Our model is shown to be consistent with a multilocus test statistic, ln RV, proposed for identifying microsatellite loci involved in recent selective sweeps. Our methodology differs in that it accepts original allele size data rather than summary statistics and allows the incorporation of prior knowledge about allele frequencies using a hierarchical prior distribution consisting of log normal and gamma probability distributions. Interesting features of the model are its ability to simultaneously analyze allele size data for any number of populations and to cope with the presence of any number of selected loci. The utility of the method is illustrated by application to two sets of microsatellite allele size data for a group of West African Anopheles gambiae populations. The results are consistent with the suppressed-recombination model of speciation, and additional candidate loci on chromosomes 2 (079 and 175) and 3 (088) are discovered that escaped former analysis.

DrosophilaPlexin B is a Sema-2a receptor required for axon guidance

Plexin receptors play a crucial role in the transduction of axonal guidance events elicited by semaphorin proteins. In Drosophila, Plexin A(PlexA) is a receptor for the transmembrane semaphorin semaphorin-1a (Sema-1a)and is required for motor and central nervous system (CNS) axon guidance in the developing embryonic nervous system. However, it remains unknown how PlexB functions during neural development and which ligands serve to activate this receptor. Here, we show that plexB, like plexA, is robustly expressed in the developing CNS and is required for motor and CNS axon pathfinding. PlexB and PlexA serve both distinct and shared neuronal guidance functions. We observe a physical association between these two plexin receptors in vivo and find that they can utilize common downstream signaling mechanisms. PlexB does not directly bind to the cytosolic semaphorin signaling component MICAL (molecule that interacts with CasL), but requires MICAL for certain axonal guidance functions. Ligand binding and genetic analyses demonstrate that PlexB is a receptor for the secreted semaphorin Sema-2a,suggesting that secreted and transmembrane semaphorins in Drosophilause PlexB and PlexA, respectively, for axon pathfinding during neural development. These results establish roles for PlexB in central and peripheral axon pathfinding, define a functional ligand for PlexB, and implicate common signaling events in plexin-mediated axonal guidance.

A genome-wide Drosophila RNAi screen identifies DYRK-family kinases as regulators of NFAT

Precise regulation of the NFAT (nuclear factor of activated T cells) family of transcription factors (NFAT1-4) is essential for vertebrate development and function. In resting cells, NFAT proteins are heavily phosphorylated and reside in the cytoplasm; in cells exposed to stimuli that raise intracellular free Ca2+ levels, they are dephosphorylated by the calmodulin-dependent phosphatase calcineurin and translocate to the nucleus. NFAT dephosphorylation by calcineurin is countered by distinct NFAT kinases, among them casein kinase 1 (CK1) and glycogen synthase kinase 3 (GSK3). Here we have used a genome-wide RNA interference (RNAi) screen in Drosophila to identify additional regulators of the signalling pathway leading from Ca2+-calcineurin to NFAT. This screen was successful because the pathways regulating NFAT subcellular localization (Ca2+ influx, Ca2+-calmodulin-calcineurin signalling and NFAT kinases) are conserved across species, even though Ca2+-regulated NFAT proteins are not themselves represented in invertebrates. Using the screen, we have identified DYRKs (dual-specificity tyrosine-phosphorylation regulated kinases) as novel regulators of NFAT. DYRK1A and DYRK2 counter calcineurin-mediated dephosphorylation of NFAT1 by directly phosphorylating the conserved serine-proline repeat 3 (SP-3) motif of the NFAT regulatory domain, thus priming further phosphorylation of the SP-2 and serine-rich region 1 (SRR-1) motifs by GSK3 and CK1, respectively. Thus, genetic screening in Drosophila can be successfully applied to cross evolutionary boundaries and identify new regulators of a transcription factor that is expressed only in vertebrates.

Genomewide comparative analysis of the highly abundant transposable element DINE-1 suggests a recent transpositional burst in Drosophila yakuba

DINE-1 (Drosophila interspersed element) is the most abundant repetitive sequence in the Drosophila genome derived from transposable elements. It comprises >1% of the Drosophila melanogaster genome (DMG) and is believed to be a relic from an ancient transpositional burst that occurred approximately 5-10 MYA. We performed a genomewide comparison of the abundance, sequence variation, and chromosomal distribution of DINE-1 in D. melanogaster and D. yakuba. Unlike the highly diverged copies in the DMG (pairwise distance approximately 15%), DINE-1's in the Drosophila yakuba genome (DYG) have diverged by only 3.4%. Moreover, the chromosomal distribution of DINE-1 in the two species is very different, with a significant number of euchromatic insertions found only in D. yakuba. We propose that these different patterns are caused by a second transpositional burst of DINE-1's in the D. yakuba genome approximately 1.5 MYA. On the basis of the sequence of these recently transposed copies, we conclude that DINE-1 is likely to be a family of nonautomomous DNA transposons. Analysis of the chromosomal distribution of two age groups of DINE-1's in D. yakuba indicates that (1) there is a negative correlation between recombination rates and the density of DINE-1's and (2) younger copies are more evenly distributed in the chromosome arms, while older copies are mostly located near the centromere regions. Our results fit the predictions of a selection-transposition balance model. Our data on whole-genome comparison of a highly abundant TE among Drosophila sibling species demonstrate the unexpectedly dynamic nature of TE activity in different host genomes.

Polycomb silencing mechanisms and genomic programming

Polycomb complexes, best known for their role in the epigenetic silencing of homeotic genes, are now known to regulate a large number of functions in organisms from flies to man. They control transcription activators, pattern-forming genes, maintenance of stem cells and are implicated in cell proliferation and oncogenesis. Our understanding of Polycomb mechanisms derives principally from the study of homeotic genes in Drosophila, where they act in an all-or-none fashion to silence expression in inappropriate parts of the organism. This review summarizes what has been learned from homeotic genes and examines the possible extensions of Polycomb mechanisms to allow for dynamic regulatory behavior and the reprogramming of silenced chromatin states.

Differential expression ofHR38in the mushroom bodies of the honeybee brain depends on the caste and division of labor

We used a cDNA microarray to identify genes expressed in a caste (worker)- and division of labor (nurse bees or foragers)-dependent manner in the honeybee brain. Among the identified genes, one encoded a putative orphan receptor (HR38) homologue that mediates ecdysteroid-signaling. Real-time reverse transcription-polymerase chain reaction indicated that expression of this gene is higher in forager brains, as compared to nurse bees and queens. In the forager brain, expression was concentrated in a subset of the mushroom body neurons, suggesting that ecdysteroid-signaling in the mushroom bodies might be involved in the division of labor of the workers.

Whitefly (Bemisia tabaci) genome project: analysis of sequenced clones from egg, instar, and adult (viruliferous and non-viruliferous) cDNA libraries

Background

The past three decades have witnessed a dramatic increase in interest in the whitefly Bemisia tabaci, owing to its nature as a taxonomically cryptic species, the damage it causes to a large number of herbaceous plants because of its specialized feeding in the phloem, and to its ability to serve as a vector of plant viruses. Among the most important plant viruses to be transmitted by B. tabaci are those in the genus Begomovirus (family, Geminiviridae). Surprisingly, little is known about the genome of this whitefly. The haploid genome size for male B. tabaci has been estimated to be approximately one billion bp by flow cytometry analysis, about five times the size of the fruitfly Drosophila melanogaster. The genes involved in whitefly development, in host range plasticity, and in begomovirus vector specificity and competency, are unknown.

Results

To address this general shortage of genomic sequence information, we have constructed three cDNA libraries from non-viruliferous whiteflies (eggs, immature instars, and adults) and two from adult insects that fed on tomato plants infected by two geminiviruses: Tomato yellow leaf curl virus (TYLCV) and Tomato mottle virus (ToMoV). In total, the sequence of 18,976 clones was determined. After quality control, and removal of 5,542 clones of mitochondrial origin 9,110 sequences remained which included 3,843 singletons and 1,017 contigs. Comparisons with public databases indicated that the libraries contained genes involved in cellular and developmental processes. In addition, approximately 1,000 bases aligned with the genome of the B. tabaci endosymbiotic bacterium Candidatus Portiera aleyrodidarum, originating primarily from the egg and instar libraries. Apart from the mitochondrial sequences, the longest and most abundant sequence encodes vitellogenin, which originated from whitefly adult libraries, indicating that much of the gene expression in this insect is directed toward the production of eggs.

Conclusion

This is the first functional genomics project involving a hemipteran (Homopteran) insect from the subtropics/tropics. The B. tabaci sequence database now provides an important tool to initiate identification of whitefly genes involved in development, behaviour, and B. tabaci-mediated begomovirus transmission.

Genome-wide identification of mRNAs associated with the translational regulator PUMILIO in Drosophila melanogaster

Genome-wide identification of RNAs associated with RNA-binding proteins is crucial for deciphering posttranscriptional regulatory systems. PUMILIO is a member of the evolutionary conserved Puf-family of RNA-binding proteins that repress gene expression posttranscriptionally. We generated transgenic flies expressing affinity-tagged PUMILIO under the control of an ovary-specific promoter, and we purified PUMILIO from whole adult flies and embryos and analyzed associated mRNAs by using DNA microarrays. Distinct sets comprising hundreds of mRNAs were associated with PUMILIO at the two developmental stages. Many of these mRNAs encode functionally related proteins, supporting a model for coordinated regulation of posttranscriptional modules by specific RNA-binding proteins. We identified a characteristic sequence motif in the 3'-untranslated regions of mRNAs associated with PUMILIO, and the sufficiency of this motif for interaction with PUMILIO was confirmed by RNA pull-down experiments with biotinylated synthetic RNAs. The RNA motif strikingly resembles the one previously identified for Puf3p, one of five Saccharomyces cerevisiae Puf proteins; however, proteins encoded by the associated mRNAs in yeast and Drosophila do not appear to be related. The results suggest extensive posttranscriptional regulation by PUMILIO and uncover evolutionary features of this conserved family of RNA-binding proteins.

Genome-wide identification of direct targets of the Drosophila retinal determination protein Eyeless

The discovery of direct downstream targets of transcription factors (TFs) is necessary for understanding the genetic mechanisms underlying complex, highly regulated processes such as development. In this report, we have used a combinatorial strategy to conduct a genome-wide search for novel direct targets of Eyeless (Ey), a key transcription factor controlling early eye development in Drosophila. To overcome the lack of high-quality consensus binding site sequences, phylogenetic shadowing of known Ey binding sites in sine oculis (so) was used to construct a position weight matrix (PWM) of the Ey protein. This PWM was then used for in silico prediction of potential binding sites in the Drosophila melanogaster genome. To reduce the false positive rate, conservation of these potential binding sites was assessed by comparing the genomic sequences from seven Drosophila species. In parallel, microarray analysis of wild-type versus ectopic ey-expressing tissue, followed by microarray-based epistasis experiments in an atonal (ato) mutant background, identified 188 genes induced by ey. Intersection of in silico predicted conserved Ey binding sites with the candidate gene list produced through expression profiling yields a list of 20 putative ey-induced, eye-enriched, ato-independent, direct targets of Ey. The accuracy of this list of genes was confirmed using both in vitro and in vivo methods. Initial analysis reveals three genes, eyes absent, shifted, and Optix, as novel direct targets of Ey. These results suggest that the integrated strategy of computational biology, genomics, and genetics is a powerful approach to identify direct downstream targets for any transcription factor genome-wide.

Repetitive DNA in tropical tasar silkworm Antheraea mylitta

Antheraea mylitta is an endemic insect species producing the world famous tasar silk. Its populations occupying different ecological and geographical regions show certain degree of phenotypic variability for which they are known as 'eco-races'. In order to understand the genetic variability and phylogenetic relationship among the different eco-races we characterized a repetitive TaqI genomic DNA fragment as a genetic marker. The sequence analysis and Southern hybridization show the repetitive nature of TaqI DNA fragment, designated as A. mylitta TaqI family repeat, AmTFR. The PCR amplification of AmTFR reveals its presence in all the tested eco-races of tasar silkworm and some other silk producing insects. The AmTFR is evenly distributed in all 31 meiotic metaphase I chromosomes as observed by fluorescent in situ hybridization. The AmTFR based phylogenetic analysis of the eco-races is not congruent with the morphological variations and their geographical distribution.

The nicotinic acetylcholine receptor Dalpha7 is required for an escape behavior in Drosophila

Acetylcholine is the major excitatory neurotransmitter in the central nervous system of insects. Mutant analysis of the Dα7 nicotinic acetylcholine receptor (nAChR) ofDrosophila shows that it is required for the giant fiber-mediated escape behavior. The Dα7 protein is enriched in the dendrites of the giant fiber, and electrophysiological analysis of the giant fiber circuit showed that sensory input to the giant fiber is disrupted, as is transmission at an identified cholinergic synapse between the peripherally synapsing interneuron and the dorsal lateral muscle motor neuron. Moreover, we found thatgfA¹, a mutation identified in a screen for giant fiber defects more than twenty years ago, is an allele ofDα7. Therefore, a combination of behavioral, electrophysiological, anatomical, and genetic data indicate an essential role for the Dα7 nAChR in giant fiber-mediated escape inDrosophila.

The authors integrate behavioral, electrophysiological, anatomical, and genetic data to establish an essential role for the Dα7 nAChR in giant fiber-mediated escape inDrosophila.

Identifying neuropeptide and protein hormone receptors in Drosophila melanogaster by exploiting genomic data

Most neuropeptide and protein hormone receptors belong to the large superfamily of G-protein-coupled receptors (GPCRs). These cell membrane proteins steer many important processes such as development, reproduction, homeostasis and behaviour when activated by their corresponding ligands. The first insect genome, that of the fruitfly Drosophila melanogaster, was sequenced in 2000, and about 200 GPCRs have been annnotated in this model insect. About 50 of these receptors were predicted to have neuropeptides or protein hormones as their ligands. Since 2000, the cDNAs of most of these candidate receptors have been cloned and for many receptors the endogenous ligand has been identified. In this review, we will give an update about the current knowledge of all Drosophila neuropeptide and protein hormone receptors, and discuss their phylogenetic relationships.

Molecular methods for assessing insect parasitism

Determining insect parasitism rates is problematic due to the small size and lack of useful distinguishing morphological characters of many parasitoid taxa. To solve this problem, entomologists have employed one of four general methods to detect parasitoid protein or nucleic acid markers: serological assay; random amplified polymorphic DNA-polymerase chain reaction (RAPD-PCR); allozyme electrophoresis; or specific PCR. Serological methods, especially with monoclonal antibodies, are unrivalled for specificity, enabling discrimination at the stage as well as species level. However, they have not found favour with many workers, possibly due to complexity and expense. RAPD-PCR has been widely used, but can only be recommended for restricted applications because of its poor reproducibility. Allozyme electrophoresis provides reproducible detection and discrimination of closely related species. Specific-PCR is highly specific and reproducible, and also has the shortest latency for detection, usually 24 h or less after parasitization. The substantial existing literature on allozyme electrophoresis and specific PCR is used to support recommendations on what are apt to be fruitful enzyme systems or genomic regions for detecting and discriminating parasitoids in untried parasitoid-host assemblages.

[Transposition of the hobo element in Drosophila melanogaster somatic cells]

Somatic mutation and recombination test on wing cells of Drosophila melanogaster showed that the recombination frequency in the somatic tissues of strains studied correlated with the presence of a full-length copy of the hobo transposable element in the genome. Transposition of hobo in somatic tissue cells at a frequency 3.5 x 10-2 per site per X chromosome was shown by fluorescence in situ hybridization with salivary gland polytene chromosomes of larvae of one of the D. melanogaster strains having a full-length hobo copy.

Cloning and structural characterization of the 6-phosphogluconate dehydrogenase locus of the medfly Ceratitis capitata and the olive fruit fly Bactrocera oleae

The pentose phosphate cycle is considered as a major source of NADPH and pentose needed for nucleic acid biosynthesis. 6-Phosphogluconate dehydrogenase (6PGD), an enzyme participating in this cycle, catalyzes the oxidative decarboxylation of 6PGD to ribulose 5-phosphate with the subsequent release of CO(2) and the reduction of NADP. We have determined the genomic sequences of 6PGD of two species of Tephritidae, the medfly Ceratitis capitata and olive fruit fly Bactrocera oleae, and constructed a three-dimensional model of 6PGD of C. capitata based on the homologous known sheep structure. In a comparative study of 6PGD sequences from seven species, all the conserved and variable sites of the enzyme were analyzed and the regions of functional importance were localized, an attempt promoted also by the direct involvement of the enzyme in various human diseases. The enzymes between the two species of Tephritidae have a very high homology and further examination of the variable positions with respect to the highly conserved binding site residues enabled their grouping in three distinct categories, with possible association to dimer formation, functional specificity, and antigenicity. Moreover, placement of sequence differences on the 3-D model suggests probable sites accommodating variations appearing at the allozymic variants of both species.

[18] Mapping Drosophila Genomic Aberration Breakpoints with Comparative Genome Hybridization on Microarrays

Chromosomal aberrations are genetic "reagents" that are commonly used in Drosophila research. Stocks containing chromosomes carrying large deletions of DNA (deficiency stocks, designated Df) as well as stocks carrying an extra copy of a chromosomal region (duplication stocks, designated Dp) are essential for a variety of genetic analyses. The extent of what is deleted or duplicated has typically been determined cytologically by salivary gland polytene chromosome squashes, which identify the edges of the aberration (so-called breakpoints) of each Df or Dp at low resolution. The margin of error for this technique can be quite high, however, because it is dependent on the quality of the squash and the experience of the scientist interpreting the data. Comparative genome hybridization on microarrays provides a precise molecular method to identify which regions of the genome are deleted or duplicated in these stocks by examining a change in chromosomal ploidy across the whole genome. Furthermore, this technique allows genetic data obtained with these strains to be placed in a molecular genomic context.

Genomic structure of the luciferase gene from the bioluminescent beetle, Nyctophila cf. caucasica

The gene coding for beetle luciferase, the enzyme responsible for bioluminescence in over two thousand coleopteran species has, to date, only been characterized from one Palearctic species of Lampyridae. Here we report the characterization of the luciferase gene from a female beetle of an Iranian lampyrid species, Nyctophila cf. caucasica (Coleoptera:Lampyridae). The luciferase gene was composed of seven exons, coding for 547 amino acids, separated by six introns spanning 1976 bp of genomic DNA. The deduced amino acid sequences of the luciferase gene of N. caucasica showed 98.9% homology to that of the Palearctic species Lampyris noctiluca. Analysis of the 810 bp upstream region of the luciferase gene revealed three TATA boxes and several other consensus transcriptional factor recognition sequences presenting evidence for a putative core promoter region conserved in Lampyrinae from -190 through to -155 upstream of the luciferase start codon. Along with the core promoter region the luciferase gene was compared with orthologous sequences from other lampyrid species and found to have greatest identity to Lampyris turkistanicus and Lampyris noctiluca. The significant sequence identity to the former is discussed in relation to taxonomic issues of Iranian lampyrids.

Sequencing a new target genome: the Boophilus microplus (Acari: Ixodidae) genome project

The southern cattle tick, Boophilus microplus (Canestrini), causes annual economic losses in the hundreds of millions of dollars to cattle producers throughout the world, and ranks as the most economically important tick from a global perspective. Control failures attributable to the development of pesticide resistance have become commonplace, and novel control technologies are needed. The availability of the genome sequence will facilitate the development of these new technologies, and we are proposing sequencing to a 4-6X draft coverage. Many existing biological resources are available to facilitate a genome sequencing project, including several inbred laboratory tick strains, a database of approximately 45,000 expressed sequence tags compiled into a B. microplus Gene Index, a bacterial artificial chromosome (BAC) library, an established B. microplus cell line, and genomic DNA suitable for library synthesis. Collaborative projects are underway to map BACs and cDNAs to specific chromosomes and to sequence selected BAC clones. When completed, the genome sequences from the cow, B. microplus, and the B. microplus-borne pathogens Babesia bovis and Anaplasma marginale will enhance studies of host-vector-pathogen systems. Genes involved in the regeneration of amputated tick limbs and transitions through developmental stages are largely unknown. Studies of these and other interesting biological questions will be advanced by tick genome sequence data. Comparative genomics offers the prospect of new insight into many, perhaps all, aspects of the biology of ticks and the pathogens they transmit to farm animals and people. The B. microplus genome sequence will fill a major gap in comparative genomics: a sequence from the Metastriata lineage of ticks. The purpose of the article is to synergize interest in and provide rationales for sequencing the genome of B. microplus and for publicizing currently available genomic resources for this tick.

Learning to fly—getting the best out of microarrays

Microarrays are becoming pervasive in biology and can offer unparalleled insights into biological processes. This article suggests that a more considered approach to the design of array probe sequences can substantially improve the technology.

Identification of transposable elements using multiple alignments of related genomes

Accurate genome-wide cataloging of transposable elements (TEs) will facilitate our understanding of mobile DNA evolution, expose the genomic effects of TEs on the host genome, and improve the quality of assembled genomes. Using the availability of several nearly complete Drosophila genomes and developments in whole genome alignment methods, we introduce a large-scale comparative method for identifying repetitive mobile DNA regions. These regions are highly enriched for transposable elements. Our method has two main features distinguishing it from other repeat-finding methods. First, rather than relying on sequence similarity to determine the location of repeats, the genomic artifacts of the transposition mechanism itself are systematically tracked in the context of multiple alignments. Second, we can derive bounds on the age of each repeat instance based on the phylogenetic species tree. We report results obtained using both complete and draft sequences of four closely related Drosophila genomes and validate our results with manually curated TE annotations in the Drosophila melanogaster euchromatin. We show the utility of our findings in exploring both transposable elements and their host genomes: In the study of TEs, we offer predictions for novel families, annotate new insertions of known families, and show data that support the hypothesis that all known TE families in D. melanogaster were recently active; in the study of the host, we show how our findings can be used to determine shifts in the eu-heterochromatin junction in the pericentric chromosome regions.

A low-cost open-source SNP genotyping platform for association mapping applications

Association mapping aimed at identifying DNA polymorphisms that contribute to variation in complex traits entails genotyping a large number of single-nucleotide polymorphisms (SNPs) in a very large panel of individuals. Few technologies, however, provide inexpensive high-throughput genotyping. Here, we present an efficient approach developed specifically for genotyping large fixed panels of diploid individuals. The cost-effective, open-source nature of our methodology may make it particularly attractive to those working in nonmodel systems.