The Drosophila melanogaster Genetic Reference Panel.
Journal: 2012/March - Nature
ISSN: 1476-4687
Abstract:
A major challenge of biology is understanding the relationship between molecular genetic variation and variation in quantitative traits, including fitness. This relationship determines our ability to predict phenotypes from genotypes and to understand how evolutionary forces shape variation within and between species. Previous efforts to dissect the genotype-phenotype map were based on incomplete genotypic information. Here, we describe the Drosophila melanogaster Genetic Reference Panel (DGRP), a community resource for analysis of population genomics and quantitative traits. The DGRP consists of fully sequenced inbred lines derived from a natural population. Population genomic analyses reveal reduced polymorphism in centromeric autosomal regions and the X chromosome, evidence for positive and negative selection, and rapid evolution of the X chromosome. Many variants in novel genes, most at low frequency, are associated with quantitative traits and explain a large fraction of the phenotypic variance. The DGRP facilitates genotype-phenotype mapping using the power of Drosophila genetics.
Relations:
Content
Citations
(487)
References
(40)
Diseases
(1)
Genes
(54)
Organisms
(2)
Processes
(6)
Anatomy
(4)
Affiliates
(1)
Similar articles
Articles by the same authors
Discussion board
Nature. Feb/8/2012; 482(7384): 173-178
Published online Feb/7/2012

The Drosophila melanogaster Genetic Reference Panel

+43 authors

Abstract

A major challenge of biology is understanding the relationship between molecular genetic variation and variation in quantitative traits, including fitness. This relationship determines our ability to predict phenotypes from genotypes and to understand how evolutionary forces shape variation within and between species. Previous efforts to dissect the genotype-phenotype map were based on incomplete genotypic information. Here, we describe the Drosophila melanogaster Genetic Reference Panel (DGRP), a community resource for analysis of population genomics and quantitative traits. The DGRP consists of fully sequenced inbred lines derived from a natural population. Population genomic analyses reveal reduced polymorphism in centromeric autosomal regions and the X chromosome, evidence for positive and negative selection, and rapid evolution of the X chromosome. Many variants in novel genes, most at low frequency, are associated with quantitative traits and explain a large fraction of the phenotypic variance. The DGRP facilitates genotype-phenotype mapping using the power of Drosophila genetics.

Understanding how molecular variation maps to phenotypic variation for quantitative traits is central for understanding evolution, animal and plant breeding, and personalized medicine1,2. The principles of mapping quantitative trait loci (QTLs) by linkage to, or association with, marker loci are conceptually simple1,2. However, we have not yet achieved our goal of explaining genetic variation for quantitative traits in terms of the underlying genes; additive, epistatic, and pleiotropic effects as well as phenotypic plasticity of segregating alleles; and the molecular nature, population frequency and evolutionary dynamics of causal variants. Efforts to dissect the genotype-phenotype map in model organisms3,4 and humans 5-7 have revealed unexpected complexities, implicating many, novel loci; pervasive pleiotropy; and context-dependent effects.

Model organism reference populations of inbred strains that can be shared among laboratories studying diverse phenotypes, and for which environmental conditions can be controlled and manipulated, greatly facilitate efforts to dissect the genetic architecture of quantitative traits3,4. Measuring many individuals of the same homozygous genotype increases the accuracy of the estimates of genotypic value1 and the power to detect variants; and genotypes of molecular markers need only be obtained once. We constructed the Drosophila melanogaster Genetic Reference Panel (DGRP) as such a community resource. Unlike previous populations of recombinant inbred lines derived from limited samples of genetic variation, the DGRP consists of 192 inbred strains derived from a single outbred population. The DGRP contains a representative sample of naturally segregating genetic variation, has an ultra-fine grained recombination map suitable for precise localization of causal variants, and has almost complete euchromatic sequence information.

Here, we describe molecular and phenotypic variation in 168 re-sequenced lines comprising Freeze 1.0 of the DGRP, population genomic inferences of patterns of polymorphism and divergence and their correlation with genomic features, local recombination rate and selection acting on this population, genome wide association mapping analyses for three quantitative traits, and tools facilitating the use of this resource.

Molecular variation in the DGRP

We constructed the DGRP by collecting mated females from the Raleigh, USA population, followed by 20 generations of full-sib inbreeding of their progeny. We sequenced 168 DGRP lines using a combination of Illumina and 454 sequencing technology: 29 of the lines were sequenced using both platforms, 129 lines have only Illumina sequence, and 10 lines have only 454 sequence. We mapped sequence reads to the D. melanogaster reference genome, re-calibrated base quality scores, and locally re-aligned Illumina reads. Mean sequence coverage was 21.4x per line for Illumina sequences and 12.1x per line for 454 sequences (Supplementary Table 1). On average, we assayed 113.5Mb (94.25%) of the euchromatic reference sequence with ~22,000 read mapping gaps per line (Supplementary Table 2). We called 4,672,297 single nucleotide polymorphisms (SNPs) using the Joint Genotyper for Inbred Lines (JGIL)8, which takes into account coverage and quality sequencing statistics, and expected allele frequencies after 20 generations of inbreeding from an outbred population initially in Hardy-Weinberg equilibrium. In cases where base calls were made by both technologies, concordance was 99.36% (Supplementary Table 3).

The SNP site frequency distribution (Fig. 1a) is characterized by a majority of low frequency variants. The numbers of SNPs vary by chromosome and site class (Fig. 1b). Linkage disequilibrium9 (LD) decays to r2 = 0.2 on average within 10 bp on autosomes and 30 bp on the X chromosome (Fig. 1c, Supplementary Fig. 1). This difference is expected because the population size of the X chromosome is ¾ that of autosomes, and the X chromosome can experience greater purifying selection due to exposure of deleterious recessive alleles in hemizygous males. There is little evidence of global population structure in the DGRP (Fig. 1d, Supplementary Fig. 2). The rapid decline in LD locally and lack of global population structure are favorable for genome wide association mapping.

Not all SNPs are fixed within individual DGRP lines (Supplementary Table 4). The expected inbreeding coefficient (F) after 20 generations of full-sib inbreeding1 is F = 0.986; therefore, we expect some SNPs to remain segregating by chance. Segregating SNPs can also arise from new mutations, or if natural selection opposes inbreeding, due to true overdominance for fitness at individual loci or associative overdominance due to complementary deleterious alleles that are closely linked or in segregating inversions.

We identified 390,873 microsatellite loci, 105,799 of which were polymorphic (Supplementary Table 5); 36,810 transposable element (TE) insertion sites and 197,402 total insertions (Supplementary Table 6). On average, each line contained 1,175 TE insertions (Supplementary Table 6), although most TE insertion sites (25,562) were present in only one line (Supplementary Table 7). We identified 149 TE families. The number of copies per family varied greatly from an average of 315.7 INE-1 elements per line to an average of 0.003 Gandalf-Dkoe-like elements per line (Supplementary Table 8).

Wolbachia pipientis is a maternally inherited bacterium found in insects, including Drosophila, and can affect reproduction10. We assessed Wolbachia infection status in the DGRP lines to account for it in analyses of genotype-phenotype associations, and found 51.2% of lines harboring sufficient Wolbachia DNA to imply infection (Supplementary Table 9).

Polymorphism and divergence

We used the DGRP Illumina sequence data and genome sequences from D. simulans and D. yakuba11 to perform genome wide analyses of polymorphism and divergence, assess the association of these parameters with genomic features and the recombination landscape, and infer the historical action of selection on a much larger scale than had been possible previously12-17. We computed polymorphism (π18 and θ19) and divergence (k20) for the whole genome, by chromosome arm (X, 2L, 2R, 3L, 3R), by chromosome region (three regions of equal size in Mb – telomeric, middle and centromeric), in 50 kbp non-overlapping windows, and by site class (synonymous and non-synonymous sites within coding sequences, and intronic, UTR and intergenic sites) (Supplementary Tables 10 and 11).

Averaged over the entire genome, π = 0.0056 and θ = 0.0067, similar to previous estimates from North American populations17,21. Average polymorphism on the X chromosome (πX = 0.0040) is reduced relative to the autosomes (πA = 0.0060) (X/A ratio = 0.67, Wilcoxon test p = 0), even after correcting for the X/A effective population size (X4/3 = 0.0054, Wilcoxon test p < 0.00002; Supplementary Table 10). Autosomal nucleotide diversity is reduced on average 2.4 fold in centromeric regions relative to non-centromeric regions, and at the telomeres (Fig. 2a, Supplementary Table 10), while diversity is relatively constant along the X chromosome. Thus, πX > πA in centromeric regions, but πA > πX in other chromosomal regions (Fig. 2a, Supplementary Table 10).

Genes on the X chromosome evolve faster (kX = 0.140) than autosomal genes (kA = 0.126) (X/A ratio = 1.131, Wilcoxon test p =0) (Figure 2b, Supplementary Table 10). Divergence is more uniform (CVk = 0.2841) across chromosome arms than is polymorphism (CVπ = 0.4265). The peaks of divergence near the centromeres could be attributable to the reduced quality of alignments in these regions. Patterns of divergence are similar regardless of the outgroup species used (Fig. 2b, Supplementary Table 11).

The pattern of polymorphism and divergence by site class is consistent within and among chromosomes (π(k)Synonymous > π(k)Intron > π(k)Intergenic > π(k)UTR > π(k)Non-synonymous), in agreement with previous studies on smaller data sets13,16 (Supplementary Figs. 3 and 4, Supplementary Table 11). Polymorphism levels between synonymous and non-synonymous sites differ by an order of magnitude. Variation and divergence patterns within the site classes generally follow the same patterns observed overall, with reduced polymorphism for all site classes on the X chromosome relative to autosomes, increased X chromosome divergence relative to autosomes for all but synonymous sites, decreased polymorphism in centromeric regions, and greater variation among regions and arms for polymorphism than for divergence. Other diversity measures and more detailed patterns at different window-sizes for each chromosome arm can be accessed from the Population Drosophila Browser (PopDrowser) (Table 1, Supplementary Methods).

Recombination landscape

Evolutionary models of hitchhiking and background selection22,23 predict a positive correlation between polymorphism and recombination rate. This expectation is realized in regions where recombination is < 2 cM/Mb (Spearman’s ρ = 0.471, p = 0), but recombination and polymorphism are independent in regions where recombination exceeds 2 cM/Mb (Spearman’s ρ = −0.0044, p = 0.987) (Fig. 2a, Supplementary Table 12). The average rate of recombination of the X chromosome (2.9 cM/Mb) is greater than that of autosomes (2.1 cM/Mb), which may account for the low overall X-linked correlation between recombination rate and π. The lack of correlation between recombination and divergence (Supplementary Table 12) excludes mutation associated with recombination as the cause of the correlation. We assessed the independent effects of recombination rate, divergence, chromosome region and gene density on nucleotide variation of autosomes and the X chromosome (Supplementary Table 13). Recombination is the major predictor of polymorphism on the X chromosome and autosomes; however, the significant effect of autosomal chromosome region remains after accounting for variation in recombination rates between centromeric and non-centromeric regions.

Selection regimes

We used the standard24 and generalized13,25,26 McDonald Kreitman tests (MKT) to scan the genome for evidence of selection. These tests compare the ratio of polymorphism at a selected site with that of a neutral site to the ratio of divergence at a selected site to divergence at a neutral site. The standard MKT is applied to coding sequences, and synonymous and non-synonymous sites are used as putative neutral and selected sites, respectively. The generalized MKT is applied to non-coding sequences and uses four-fold degenerate sites as neutral sites. Using polymorphism and divergence data avoids confounding inference of selection with mutation rate differences, and restricting the tests to closely linked sites controls for shared evolutionary history27-29. We infer adaptive divergence when there is an excess of divergence relative to polymorphism, and segregation of slightly deleterious mutations when there is an excess of polymorphism over divergence. Estimates of α, the proportion of adaptive divergence, are biased downwards by low frequency, slightly deleterious mutations30,31. Rather than eliminate low frequency variants32, we incorporated information on the site frequency distribution to the MKT test framework to obtain estimates of the proportion of sites that are strongly deleterious (d), weakly deleterious (b), neutral (f), and recently neutral ( ) at segregating sites, as well as unbiased estimates of α (Supplementary Methods).

Deleterious and neutral sites

Averaged over the entire genome, we infer that 58.5% of the segregating sites are neutral or nearly neutral, 1.9% are weakly deleterious, and 39.6% are strongly deleterious. However, these proportions vary between the X chromosome and autosomes, site classes, and chromosome regions (Supplementary Tables 14-16, Fig. 3). Non-synonymous sites are the most constrained (d = 77.6%), while in non-coding sites d ranges from 29.1% in 5′UTRs to 41.3% in 3′intergenic regions. The inferred pattern of selection differs between autosomal centromeric and non-centromeric regions: d is reduced and f is increased in centromeric regions for all site categories (Fig. 3). We observe an excess of polymorphism relative to divergence in autosomal centromeric regions, even after correcting for weakly deleterious mutations, implying a relaxation of selection from the time of separation of D. melanogaster and D. yakuba. Since selection coefficients depend on the effective population size33 (Ne), this could occur if the recombination rate has specifically diminished in centromeric regions during the divergence between D. melanogaster and D. yakuba; or with an overall reduction of Ne associated with the colonization of North American habitats34,35. In the latter case, we expect a genome-wide signature of an excess of low frequency polymorphisms and of polymorphism relative to divergence, exacerbated in regions of low recombination. We indeed find an excess of low frequency polymorphism relative to neutral expectation as indicated by the negative estimates of Tajima’s D statistic36 ( D = −0.686 averaged over the whole genome and D = −0.997 in autosomal centromeric regions). In contrast, the X chromosome does not show a differential pattern of selection in the centromeric region, has a lower fraction of relaxation of selection, fewer neutral alleles, and a higher percentage of strongly deleterious alleles for all site classes and regions (Fig. 3, Supplementary Tables 14-16).

TE insertions are thought to be largely deleterious. There are more singleton insertions in regions of high recombination (≥ 2cM/Mb) and more insertions shared in multiple lines in regions of low recombination (< 2cM/Mb) (Fisher’s exact test p = 0), and comparison of observed and expected site occupancy spectra reveals an excess of singleton insertions (p = 0, Supplementary Fig. 5).

Adaptive fixation

We find substantial evidence for positive selection in autosomal non-centromeric regions and the X chromosome (Fig. 2c, Supplementary Tables 15, 17). We estimated α by aggregating all sites in each region analyzed to avoid underestimation by averaging across genes37 in comparisons of chromosomes, regions and site classes. We also computed the Direction of Selection, DoS38, which is positive with adaptive selection, zero under neutrality, and negative when weakly deleterious or new nearly neutral mutations are segregating. Estimates of α from the standard and generalized MKT indicate that on average 25.2% of the fixed sites between D. melanogaster and D. yakuba are adaptive, ranging from 30% in introns to 7% in UTR sites (Supplementary Fig. 6). Estimates of DoS and α are negative for non-synonymous and UTR sites in the autosomal centromeres, consistent with underestimating the fraction of adaptive substitutions in regions of low recombination because weakly deleterious or nearly neutral mutations are more common than adaptive fixations. The majority of adaptive fixation on autosomes occurs in non-centromeric regions (Fig. 2c). We find over four times as many adaptive fixations on the X chromosome relative to autosomes. The pattern holds for all site classes, in particular non-synonymous sites and UTRs, as well as individual genes; and is not solely due to the autosomal centromeric effect (Supplementary Table 15; Supplementary Figs. 6 and 7). Finally, when we consider DoS in recombination environments above and below 2 cM/Mb, we find greater adaptive propensity in genes whose recombination context is ≥ 2 cM/Mb (Wilcoxon test, p = 0) (Supplementary Fig. 8).

To understand the global patterns of divergence and constraint across functional classes of genes, we examined the distributions of ω (dN/dS, the ratio of non-synonymous to synonymous divergence) and DoS across Gene Ontology (GO) categories. The 4.9% GO categories with significantly elevated DoS include the biological process categories of behavior, developmental process involved in reproduction, reproduction and ion transport (Supplementary Table 18). Recombination context is the major determinant of variation in DoS (Supplementary Table 19) while GO category is as important as recombinational context for predicting variation in ω (Supplementary Table 19).

GO categories enriched for positive DoS values differ from those associated with high values of ω (Supplementary Table 18), indicating that positive selection does not occur necessarily on genes with high ω values. If adaptive substitutions are common, high values of ω reflect the joint contributions of neutral and adaptive substitutions. Further, equating high constraint (low ω) with functional importance overlooks the functional role of adaptive changes16. Unlike ω, DoS takes into account the constraints inferred from the current polymorphism, distinguishing negative, neutral and adaptive selection.

Genome wide association of phenotypes with molecular variation

We measured resistance to starvation stress, chill coma recovery time, and startle response39 in the DGRP. We found considerable genetic variation for all traits, with high broad sense heritabilities. We also found variation in sex dimorphism for starvation resistance and chill coma recovery with cross-sex genetic correlations significantly different from unity (Supplementary Tables 20-22).

We performed genome wide association (GWA) analyses for these traits, using the 2,490,165 SNPs and 77,756 microsatellites for which the minor allele was represented in four or more lines, using single locus analyses pooled across sexes and separately for males and females. At p < 10−5 (p < 10−6), we find 203 (32) SNPs and 2 (0) microsatellites associated with starvation resistance; 90 (7) SNPs and 4 (2) microsatellites associated with startle response; and 235 (45) SNPs and 5 (3) microsatellites associated with chill coma recovery time (Fig. 4a, Supplementary Fig. 9, Supplementary Tables 23 and 24). The minor allele frequencies for most of the associated SNPs are low, and there is an inverse relationship between effect sizes and minor allele frequency (Supplementary Fig. 10).

The DGRP is a powerful tool for rapidly reducing the search space for molecular variants affecting quantitative traits from the entire genome to candidate polymorphisms and genes. While we cannot infer which of these polymorphisms are causal due to LD between SNPs in close physical proximity as well as occasional spurious long range LD (Fig. 4a, Supplementary Fig. 9), the candidate gene lists are likely to be enriched for causal variants. The majority of associations are in computationally predicted genes or genes with annotated functions not obviously associated with the three traits. However, genes previously associated with startle response40 (Sema-1a and Eip-75B) and starvation resistance41 (pnt) were identified in this study; and a SNP in CG3213, previously identified in a Drosophila obesity screen42, is associated with variation in starvation resistance. Several genes associated with quantitative traits are rapidly evolving (psq, Egfr; Supplementary Tables 17, 23) or are plausible candidates based on SNP or Gene Ontology annotations (Supplementary Table 23).

Predicting phenotypes from genotypes

We used regression models to predict trait phenotypes from SNP genotypes and estimate the total variance explained by SNPs. The latter cannot be done by summing the individual contributions of the single marker effects because markers are not completely independent, and estimates of effects of single markers are biased when more than one locus affecting the trait segregates in the population. We derived gene-centered multiple regression models to estimate the effects of multiple SNPs simultaneously. In all cases 6-10 SNPs explain from 51-72% of the phenotypic variance and 65-90% of the genetic variance (Supplementary Tables 25 and 26; Supplementary Figs. 11-13). We also derived partial least square regression models using all SNPs for which the single marker effect was significant at p < 10−5. These models explain 72-85% of the phenotypic variance (Fig. 4b-c, Supplementary Fig. 14).

Discussion

The DGRP lines, sequences, variant calls, phenotypes, and web tools for molecular population genomics and GWA analysis are publicly available (Table 1). The DGRP lines contain at least 4,672,297 SNPs, 105,799 polymorphic microsatellites and 36,810 TEs, as well as insertion/deletion events and copy number variants and are a valuable resource for understanding the genetic architecture of quantitative traits of ecological and evolutionary relevance as well as Drosophila models of human quantitative traits. These novel mutations have survived the sieve of natural selection and will enhance the functional annotation of the Drosophila genome, complementing the Drosophila Gene Disruption Project43 and the Drosophila modENCODE project44. .

Genome-wide molecular population genetic analyses show that patterns of polymorphism, but not divergence, differ by autosomal chromosome region, and between the X chromosome and autosomes. Polymorphism is lower in autosomal centromeric than non-centromeric regions, but not for the X chromosome. We hypothesize that the correlation of polymorphism with recombination in regions where recombination is < 2 cM/Mb is due to the reduced effective population size in regions of low recombination9. Selection is less efficient in regions of low recombination33, consistent with our observation that the fraction of strongly deleterious mutations and positively selected sites are reduced in these regions.

All molecular population genomic analyses support the ‘faster X’ hypothesis45. Relative to the autosomes, the X chromosome exhibits lower polymorphism, faster rates of molecular evolution, a higher percentage of gene regions undergoing adaptive evolution, a higher fraction of strongly deleterious sites, and a lower level of weak negative selection and relaxation of selection. New X-linked mutations are directly exposed to selection each generation in hemizygous males, and the X chromosome has greater recombination than autosomes45; both of these factors could contribute to this observation.

GWA analyses of three fitness-related quantitative traits reveal hundreds of novel candidate genes, highlighting our ignorance of the genetic basis of complex traits. Most variants associated with the traits are at low frequency, and there is an inverse relationship between frequency and effect. Given that low frequency alleles are likely to be deleterious for traits under directional or stabilizing selection, these results are consistent with the mutation-selection balance hypothesis1 for the maintenance of quantitative genetic variation. Regression models incorporating significant SNPs explain most of the phenotypic variance of the traits, in contrast with human association studies, where significant SNPs have tiny effects and together explain a small fraction of the total phenotypic variance7. If the genetic architecture of human complex traits is also dominated by low frequency causal alleles, we expect estimates of effect size based on LD with common variants to be strongly biased downwards.

In the future, the full power of Drosophila genetics can be applied to validating marker-trait associations: mutations, RNAi constructs and QTL mapping populations. The DGRP is an ideal resource for systems genetics analyses of the relationship between molecular variation, causal molecular networks, and genetic variation for complex traits4,39,46, and will anchor evolutionary studies in comparison with sequenced Drosophila species to assess to what extent variation within a species corresponds to variation among species.

Supplementary Material

Fig. 1

SNP variation in the DGRP lines. (a) Site frequency spectrum. (b) Numbers of SNPs per site class. (c) Decay of LD (r2) with physical distance for the five major chromosome arms. (d) Lack of population structure. The red curve depicts the ranked eigenvalues of the genetic covariance matrix in decreasing order with respect to the marginal variance explained; the blue curve shows their cumulative sum as a fraction of the total with respect to cumulative variance explained. The partitioning of total genetic variance is balanced among the eigenvectors. The principal eigenvector explains <1.1% of the total genetic variance.

Fig. 2

Pattern of polymorphism, divergence, α and recombination rate along chromosome arms in non-overlapping 50 kbp windows. (a) Nucleotide polymorphism (π). The solid curves give the recombination rate (cM/Mb). (b) Divergence (k) for D. simulans (light green) and D. yakuba (dark green). (c) Polymorphism/Divergence, estimated as 1- [(π0-fold/π4-fold)/(k0-fold/k4-fold)]. An excess of 0-fold divergence relative to polymorphism (k0-fold/k4-fold) > (π0-fold/π4-fold) is interpreted as adaptive fixation while an excess of 0-fold polymorphism relative to divergence (π0-fold/π4-fold) > (k0-fold/k4-fold) indicates weakly deleterious or nearly neutral mutations are segregating in the population. Tel: telomere; Cent: centromere.

Fig. 3

The fraction of alleles segregating under different selection regimes by site class and chromosome region, for the autosomes (A) and the X chromosome (X). The selection regimes are strongly deleterious (d, dark blue), weakly deleterious (b, blue), recently neutral (γ, white) and old neutral (f-γ, light blue) Each chromosome arm has been divided in three regions of equal size (in Mb): centromere, middle and telomere.

Fig. 4

Genotype-phenotype associations for starvation resistance. (a) GWA results for significant SNPs. The lower triangle depicts LD (r2) among SNPs, with the five major chromosome arms demarcated by black lines. The upper panels give the significance threshold (-log10p, uncorrected for multiple tests), the effect in phenotypic standard deviation units, and the minor allele frequency (MAF). (b-c) Partial Least Squares regressions of phenotypes predicted using SNP data on observed phenotypes. The blue dots represent the predicted and observed phenotypes of lines that were not included in the initial study. (b) Females (r2 = 0.81). (c) Males (r2 = 0.85).

Table 1
Community resources
ResourceLocation
DGRP LinesBloomington Drosophila Stock Center
http://flystocks.bio.indiana.edu/Browse/RAL.php
SequencesBaylor College of Medicine Human Genome Sequencing Center
http://www.hgsc.bcm.tmc.edu/project-species-i-DGRP_lines.hgsc
National Center for Biotechnology Information Short Read Archive
http://www.ncbi.nlm.nih.gov/sra?term=DGRP
Mackay Laboratory
http://dgrp.gnets.ncsu.edu/
Read AlignmentsBaylor College of Medicine Human Genome Sequencing Center
http://www.hgsc.bcm.tmc.edu/projects/dgrp/
SNPsBaylor College of Medicine Human Genome Sequencing Center
http://www.hgsc.bcm.tmc.edu/projects/dgrp/freeze1_July_2010/snp_calls/
National Center for Biotechnology Information dbSNP
http://www.ncbi.nlm.nih.gov/SNP/snp_viewBatch.cgi?sbid=1052186
Mackay Laboratory
http://dgrp.gnets.ncsu.edu/
MicrosatellitesBaylor College of Medicine Human Genome Sequencing Center
http://www.hgsc.bcm.tmc.edu/projects/dgrp/freeze1_July_2010/microsat_calls/
Mittelman Laboratory
http://genome.vbi.vt.edu/public/DGRP/
TEsMackay Laboratory
http://dgrp.gnets.ncsu.edu/
Molecular
Population
Genomics
PopDrowser
http://PopDrowser.uab.cat
PhenotypesMackay Laboratory
http://dgrp.gnets.ncsu.edu/
GWA AnalysisMackay Laboratory
http://dgrp.gnets.ncsu.edu/

Footnotes

Author Contributions Conceived the project: T. F. C. M., S. R., R. A. G.

Wrote the main manuscript: T. F. C. M., S. R., A. B., E. A. S.

Wrote the supplementary methods: T. F. C. M., S. R., A. B., E. A. S., J. F. A., K. R. T., J. M. C., C. M. B., D. M.

Performed experiments: M. M. M., C. B., K. P. B., M. A. C., L. C., L. D., Y. H., M. J., J. C. J., S. N. J., K. W. J., F. L., F. L., S. L. L., R. F. L., M. M., D. M. M., L. N., I. M., L. P., L. L. P., C. Q., J. G. R., S. M. R., L. T., K. C. W., Y.-Q. W., A. Y., Y. Z.

Bioinformatics and data analysis: T. F. C. M., A. B., J. F. A., D. Z., S. C., M. M. M., J. M. C., M. F. R., M. B., D. C., R. S. L., A. M., C. M. B., K. R. T., D. M., E. A. S.

Methods and web site development: J. F. A., S. C., M. M. M., Z. H., P. L., M. R., J. R., E. A. S., Contributed resources: R. R. H. A.

Sequences have been deposited at the National Center for Biotechnology Information Short Read Archives: http://www.ncbi.nlm.nih.gov/sra?term=DGRP

Reprints and permissions information is available at www.nature.com/reprints

The authors declare they have no completing financial interests

Acknowledgements

This work was supported by National Institutes of Health grant GM 45146 to T. F. C. M., E. A. S and R. R. H. A; R01 GM 059469 to R. R. H. A., MCI BFU 2009-09504 to A.B., R01 GM 085183 to K. R. T., NHGRI U54 HG003273 to R.A.G.; and an award through the NVIDIA Foundation’s “Compute the Cure” program to D. M.

References

  • 1. FalconerDSMackayTFCIntroduction to Quantitative Genetics1996Addison Wesley LongmanHarlow
  • 2. LynchMWalshBGenetics and Analysis of Quantitative Traits1998Sinauer AssociatesSunderland, Massachusetts
  • 3. FlintJMackayTFCGenetic architecture of quantitative traits in flies, mice and humansGenome Res200919723733[PubMed][Google Scholar]
  • 4. MackayTFCStoneEAAyrolesJFThe genetics of quantitative traits: Challenges and prospectsNat. Rev. Genet200910565577[PubMed][Google Scholar]
  • 5. AltshulerDDalyMJLanderESGenetic mapping in human diseaseScience2008322881888[PubMed][Google Scholar]
  • 6. DonnellyPProgress and challenges in genome-wide association studies in humansNature2008456728731[PubMed][Google Scholar]
  • 7. ManolioTAFinding the missing heritability of complex diseasesNature2009461747753[PubMed][Google Scholar]
  • 8. StoneEAJoint genotyping on the fly: Identifying variation among a sequenced panel of inbred linesGenome Res2011submitted[Google Scholar]
  • 9. HillWGRobertsonAThe effect of linkage on limits to artificial selectionGenet. Res19668269294[PubMed][Google Scholar]
  • 10. WerrenJHBiology of WolbachiaAnn. Rev. Entomol199742587609[PubMed][Google Scholar]
  • 11. ClarkAGEvolution of genes and genomes on the Drosophila phylogenyNature2007450203218[PubMed][Google Scholar]
  • 12. SmithNGEyre-WalkerAAdaptive protein evolution in DrosophilaNature200241510221024[PubMed][Google Scholar]
  • 13. AndolfattoPAdaptive evolution of non-coding DNA in DrosophilaNature200543711491152[PubMed][Google Scholar]
  • 14. PresgravesDCRecombination enhances protein adaptation in Drosophila melanogasterCurr. Biol20051516511656[PubMed][Google Scholar]
  • 15. CasillasSBarbadillaABergmanCPurifying selection maintains highly conserved noncoding sequences in DrosophilaMol. Evol. Biol20072422222234[Google Scholar]
  • 16. SellaGPervasive natural selection in the Drosophila genome?PLoS Genet20095e1000495[PubMed][Google Scholar]
  • 17. SacktonTBPopulation genomic inferences from sparse high-throughput sequencing of two populations of Drosophila melanogasterGenome Biol. Evol20091449465[PubMed][Google Scholar]
  • 18. NeiMMolecular Evolutionary Genetics1987Columbia University Press
  • 19. WattersonGAOn the number of segregating sites in genetical models without recombinationTheor. Popul. Biol19757256276[PubMed][Google Scholar]
  • 20. JukesTHCantorCRMunroHNEvolution of protein moleculesMammalian Protein Metabolism III196921132Academic Press[Google Scholar]
  • 21. AndolfattoPPrzeworskiMRegions of lower crossing over harbor more rare variants in African Drosophila melanogasterGenetics2001158657665[PubMed][Google Scholar]
  • 22. BegunDJAquadroCFLevels of naturally occurring DNA polymorphism correlate with recombination rates in D. melanogasterNature1992356519520[PubMed][Google Scholar]
  • 23. CharlesworthBMorganMTCharlesworthDThe effect of deleterious mutations on neutral molecular variationGenetics199313412891303[PubMed][Google Scholar]
  • 24. McDonaldJHKreitmanMAdaptive protein evolution at the Adh locus in DrosophilaNature1991351652654[PubMed][Google Scholar]
  • 25. JenkinsDLOrtoriCABrookfieldJFA test for adaptive change in DNA sequences controlling transcriptionProc. Biol. Sci1995261203207[PubMed][Google Scholar]
  • 26. EgeaRCasillasSBarbadillaAStandard and generalized McDonald-Kreitman test: a website to detect selection by comparing different classes of DNA sitesNuc. Acids Res200836W157162[Google Scholar]
  • 27. SawyerSAHartlDLPopulation genetics of polymorphism and divergenceGenetics199213211611176[PubMed][Google Scholar]
  • 28. NielsenRStatistical tests of selective neutrality in the age of genomicsHeredity200186641647[PubMed][Google Scholar]
  • 29. Eyre-WalkerAChanging effective population size and the McDonald-Kreitman testGenetics200216220172024[PubMed][Google Scholar]
  • 30. CharlesworthJEyre-WalkerAThe McDonald-Kreitman test and slightly deleterious mutationsMol. Biol. Evol20082510071015[PubMed][Google Scholar]
  • 31. Eyre-WalkerAKeightleyPDEstimating the rate of adaptive molecular evolution in the presence of slightly deleterious mutations and population size changeMol. Evol. Biol20092620972108[Google Scholar]
  • 32. FayJCWyckoffGJWuCITesting the neutral theory of molecular evolution with genomic data from DrosophilaNature200241510241026[PubMed][Google Scholar]
  • 33. OhtaTSlightly deleterious mutant substitutions in evolutionNature19732469698[PubMed][Google Scholar]
  • 34. DavidJHCapyPGenetic variation of Drosophila melanogaster natural populationsTrends Genet19884106111[PubMed][Google Scholar]
  • 35. BegunDJAquadroCFAfrican and North American populations of Drosophila melanogaster are very different at the DNA levelNature1993365548550[PubMed][Google Scholar]
  • 36. TajimaFStatistical methods to test for nucleotide mutation hypothesis by DNA polymorphismGenetics1989123585595[PubMed][Google Scholar]
  • 37. SmithNGEyre-Walker, AAdaptive protein evolution in Drosophila. Nature200241510221024[Google Scholar]
  • 38. StoletzkiNEyre-WalkerAEstimation of the Neutrality IndexMol. Biol. Evol2011286370[PubMed][Google Scholar]
  • 39. AyrolesJFSystems genetics of complex traits in Drosophila melanogasterNat. Genet200941299307[PubMed][Google Scholar]
  • 40. YamamotoAZwartsLCallaertsPNorgaKMackayTFCAnholtRRHNeurogenetic networks for startle-induced locomotion in Drosophila melanogasterProc. Natl. Acad. Sci. USA20081051239312398[PubMed][Google Scholar]
  • 41. HarbisonSTYamamotoAHFanaraJJNorgaKKMackayTFCQuantitative trait loci affecting starvation resistance in Drosophila melanogasterGenetics200416618071823[PubMed][Google Scholar]
  • 42. PospisilikJADrosophila genome-wide obesity screen reveals hedgehog as a determinant of brown versus white adipose cell fateCell2010140148160[PubMed][Google Scholar]
  • 43. BellenHJThe BDGP gene disruption project: single transposon insertions associated with 40% of Drosophila genesGenetics2004167761781[PubMed][Google Scholar]
  • 44. ModENCODE ConsortiumIdentification of functional elements and regulatory circuits by Drosophila modENCODEScience201033017871797[PubMed][Google Scholar]
  • 45. CharlesworthBCoyneJABartonNHThe relative rates of evolution of sex chromosomes and autosomesAm. Nat1987130113146[Google Scholar]
  • 46. SiebertsSKSchadtEEMoving toward a system genetics view of diseaseMamm. Genome200718389401[PubMed][Google Scholar]
Collaboration tool especially designed for Life Science professionals.Drag-and-drop any entity to your messages.