Post-GWAS: where next? More samples, more SNPs or more biology?

The power of genome-wide association studies (GWAS) rests on several foundations: (i) there is a significant amount of additive genetic variation, (ii) individual causal polymorphisms often have sizable effects and (iii) they segregate at moderate-to-intermediate frequencies, or will be effectively 'tagged' by polymorphisms that do. Each of these assumptions has recently been questioned. (i) Why should genetic variation appear additive given that the underlying molecular networks are highly nonlinear? (ii) A new generation of relatedness-based analyses directs us back to the nearly infinitesimal model for effect sizes that quantitative genetics was long based upon. (iii) Larger effect causal polymorphisms are often low frequency, as selection might lead us to expect. Here, we review these issues and other findings that appear to question many of the foundations of the optimism GWAS prompted. We then present a roadmap emerging as one possible future for quantitative genetics. We argue that in future GWAS should move beyond purely statistical grounds. One promising approach is to build upon the combination of population genetic models and molecular biological knowledge. This combined treatment, however, requires fitting experimental data to models that are very complex, as well as accurate capturing of the uncertainty of resulting inference. This problem can be resolved through Bayesian analysis and tools such as approximate Bayesian computation-a method growing in popularity in population genetic analysis. We show a case example of anterior-posterior segmentation in Drosophila, and argue that similar approaches will be helpful as a GWAS augmentation, in human and agricultural research.

Whole-genome sequencing of two North American Drosophila melanogaster populations reveals genetic differentiation and positive selection.

The prevailing demographic model for Drosophila melanogaster suggests that the colonization of North America occurred very recently from a subset of European flies that rapidly expanded across the continent. This model implies a sudden population growth and range expansion consistent with very low or no population subdivision. As flies adapt to new environments, local adaptation events may be expected. To describe demographic and selective events during North American colonization, we have generated a data set of 35 individual whole-genome sequences from inbred lines of D. melanogaster from a west coast US population (Winters, California, USA) and compared them with a public genome data set from Raleigh (Raleigh, North Carolina, USA). We analysed nuclear and mitochondrial genomes and described levels of variation and divergence within and between these two North American D. melanogaster populations. Both populations exhibit negative values of Tajima's D across the genome, a common signature of demographic expansion. We also detected a low but significant level of genome-wide differentiation between the two populations, as well as multiple allele surfing events, which can be the result of gene drift in local subpopulations on the edge of an expansion wave. In contrast to this genome-wide pattern, we uncovered a 50-kilobase segment in chromosome arm 3L that showed all the hallmarks of a soft selective sweep in both populations. A comparison of allele frequencies within this divergent region among six populations from three continents allowed us to cluster these populations in two differentiated groups, providing evidence for the action of natural selection on a global scale.

Genetic variation in the Yolk protein expression network of Drosophila melanogaster: sex-biased negative correlations with longevity.

One of the persistent problems in biology is understanding how genetic variation contributes to phenotypic variation. Associations at many levels have been reported, and yet causal inference has remained elusive. We propose to rely on the knowledge of causal relationships established by molecular biology approaches. The existing molecular knowledge forms a firm backbone upon which hypotheses connecting genetic variation, transcriptional variation and phenotypic variation can be built. The sex determination pathway is a well-established molecular network, with the Yolk protein 1-3 (Yp) genes as the most downstream target. Our analyses reveal that genetic variation in expression for genes known to be upstream in the pathway explains variation in downstream targets. Relationships differ between the two sexes, and each Yp has a distinct transcriptional pattern. Yp expression is significantly negatively correlated with longevity, an important life history trait, for both males and females.

Allelic imbalance in Drosophila hybrid heads: exons, isoforms, and evolution.

Unraveling how regulatory divergence contributes to species differences and adaptation requires identifying functional variants from among millions of genetic differences. Analysis of allelic imbalance (AI) reveals functional genetic differences in cis regulation and has demonstrated differences in cis regulation within and between species. Regulatory mechanisms are often highly conserved, yet differences between species in gene expression are extensive. What evolutionary forces explain widespread divergence in cis regulation? AI was assessed in Drosophila melanogaster-Drosophila simulans hybrid female heads using RNA-seq technology. Mapping bias was virtually eliminated by using genotype-specific references. Allele representation in DNA sequencing was used as a prior in a novel Bayesian model for the estimation of AI in RNA. Cis regulatory divergence was common in the organs and tissues of the head with 41% of genes analyzed showing significant AI. Using existing population genomic data, the relationship between AI and patterns of sequence evolution was examined. Evidence of positive selection was found in 30% of cis regulatory divergent genes. Genes involved in defense, RNAi/RISC complex genes, and those that are sex regulated are enriched among adaptively evolving cis regulatory divergent genes. For genes in these groups, adaptive evolution may play a role in regulatory divergence between species. However, there is no evidence that adaptive evolution drives most of the cis regulatory divergence that is observed. The majority of genes showed patterns consistent with stabilizing selection and neutral evolutionary processes.

A de novo transcriptome assembly of Lucilia sericata (Diptera: Calliphoridae) with predicted alternative splices, single nucleotide polymorphisms and transcript expression estimates.

The blow fly Lucilia sericata (Diptera: Calliphoridae) (Meigen) is a nonmodel organism with no reference genome that is associated with numerous areas of research spanning the ecological, evolutionary, medical, veterinary and forensic sciences. To facilitate scientific discovery in this species, the transcriptome was assembled from more than six billion bases of Illumina and twenty-one million bases of 454 sequence derived from embryonic, larval, pupal, adult and larval salivary gland libraries. The assembly was carried out in a manner that enabled identification of putative single nucleotide polymorphisms (SNPs) and alternative splices, and that provided expression estimates for various life history stages and for salivary tissue. The assembled transcriptome was also used to identify transcribed transposable elements in L. sericata. The results of this study will enable blow fly biologists, dipterists and comparative genomicists to more rapidly develop and test molecular and genetic hypotheses, especially those regarding blow fly development and salivary gland biology.

Genome-enabled hitchhiking mapping identifies QTLs for stress resistance in natural Drosophila.

Identification of genes underlying complex traits is an important problem. Quantitative trait loci (QTL) are mapped using marker-trait co-segregation in large panels of recombinant genotypes. Most frequently, recombinant inbred lines derived from two isogenic parents are used. Segregation patterns are also studied in pedigrees from multiple families. Great advances have been made through creative use of these techniques, but narrow sampling and inadequate power represent strong limitations. Here, we propose an approach combining the strengths of both techniques. We established a mapping population from a sample of natural genotypes, and applied artificial selection for a complex character. Selection changed the frequencies of alleles in QTLs contributing to the selection response. We infer QTLs with dense genotyping microarrays by identifying blocks of linked markers undergoing selective changes in allele frequency. We demonstrated this approach with an experimental population composed from 20 isogenic strains. Selection for starvation survival was executed in three replicated populations with three control non-selected populations. Three individuals per population were genotyped using Affymetrix GeneChips. Two regions of the genome, one each on the left arms of the second and third chromosomes, showed significant divergence between control and selected populations. For the former region, we inferred allele frequencies in selected and control populations by pyrosequencing. We conclude that the allele frequency difference, averaging approximately 40% between selected and control lines, contributed to selection response. Our approach can contribute to the fine scale decomposition of the genetics of direct and indirect selection responses, and genotype by environment interactions.

Additivity and trans-acting effects on gene expression in male Drosophila simulans.

Understanding how genetic variation is maintained begins with a comprehensive description of what types of genetic variation exist, the extent and magnitude of the variation, and patterns discernable in that variation. However, such studies have focused primarily on DNA sequence data and have ignored genetic variation at other hierarchical levels of genetic information. Microarray technology permits an examination of genetic variation at the level of mRNA abundance. Utilizing a round-robin design, we present a quantitative description of variation in mRNA abundance in terms of GCA (general combining ability or additive variance). We test whether genes significant for GCA are randomly distributed across chromosomes and use a nonparametric approach to demonstrate that the magnitude of the variation is not random for GCA. We find that there is a paucity of genes significant for GCA on the X relative to the autosomes. The overall magnitude of the effects for GCA on the X tends to be lower than that on the autosomes and is contributed by rare alleles of larger effect. Due to male hemizygosity, GCA for X-linked phenotypes must be due to trans-acting factors, while GCA for autosomal phenotypes may be due to cis- or trans-acting factors. The contrast in the amount of variation between the X and the autosomes suggests that both cis and trans factors contribute to variation for expression in D. simulans with the preponderance of effects being trans. This nonrandom patterning of genetic variation in gene expression data with respect to chromosomal context may be due to hemizygosity in the male.

Accumulation of transposable elements in the genome of Drosophila melanogaster is associated with a decrease in fitness.

Replicates of the two isogenic laboratory strains of Drosophila melanogaster, 2b and Harwich, contain different average transposable element (TE) copy numbers in the same genetic background. These lines were used to analyze the correlation between TE copy number and fitness. Assuming a weak deleterious effect of each TE insertion, a decrease in fitness is expected with an increase in genomic TE copy number. Higher rates of ectopic exchanges and, consequently, chromosomal rearrangements resulting in early embryonic death are also predicted from an increase in TE copy number. Therefore egg hatchability is expected to decrease as the genomic TE copy number increases. In 2b, where replicate lines have diverged up by 90 TE copies per haploid genome, a negative correlation between the number of TE insertions and both fitness and egg hatchability were found. Neither correlation was significant for the Harwich replicates, which have only diverged by 30 TE copies. The average deleterious effect of a TE insertion on fitness and its components was estimated as 0.004. Both homozygous and heterozygous TE insertions were shown to have deleterious effects on fitness and its components.

[Heterosis of quantitative trait loci modifying lifespan in Drosophila melanogaster]. / Geterozis poligennykh lokusov, modifitsiruiushchikh prodolzhitel'nost' zhizni u Drosophila melanogaster.

Knowledge of genes responsible for aging and death is a prerequisite for determining the relative contributions of the different evolutionary factors responsible for the limited duration of life. Polymorphism of these genes probably accounts for the variation in lifespan. Previously, quantitative trait loci (QTLs) controlling this variation were mapped with the use of 98 recombinant inbred (RI) lines originating from two parental isogenic Drosophila melanogaster stocks. In each RI line, lifespan was measured for 25 males and 25 females, and alleles were established for 93 marker genes segregating between the parental lines. Significant correlation between marker segregation and lifespan was revealed for several chromosome regions. The lifespan genes had sex-specific effects and late age onset. In the present work, the effects of the QTLs were compared for homozygous and heterozygous flies. In Six out of the eight detected QTLs alleles that decreased lifespan were recessive. Heterosis was observed for a of QTL at 33E-38A. Thus, heterosis might contribute to maintaining variation in lifespan in natural populations.

Quantitative trait loci affecting a courtship signal in Drosophila melanogaster.

Courtship plays a major role in the sexual isolation of species, yet the genetics underlying courtship behaviour are poorly understood. Here we analyse quantitative trait loci (QTL) for a major component of courtship song in recombinant inbred lines derived from two laboratory strains of Drosophila melanogaster. The total variance among lines exceeds that between parental strains, and is broadly similar to that seen among geographic strains of the Cosmopolitan form of this species. Previous studies of the quantitative genetics of fly song have implied a polygenic additive inheritance with numerous genes spread throughout the genome. We find evidence for only three significant QTLs explaining 54% of the genetic variance in total. Thus there is evidence for a few large effect genes contributing to the genetic variance among lines. Interestingly, almost all of the candidate song genes previously described for D. melanogaster do not coincide with our QTLs.

Quantitative trait locus mapping of fitness-related traits in Drosophila melanogaster.

We examined the genetic architecture of four fitness-related traits (reproductive success, ovariole number, body size and early fecundity) in a panel of 98 Oregon-R x 2b3 recombinant inbred lines (RILs). Highly significant genetic variation was observed in this population for female, but not male, reproductive success. The cross-sex genetic correlation for reproductive success was 0.20, which is not significantly different from zero. There was significant genetic variation segregating in this cross for ovariole number, but not for body size or early fecundity. The RILs were genotyped for cytological insertion sites of roo transposable elements, yielding 76 informative markers with an average spacing of 3.2 cM. Quantitative trait loci (QTL) affecting female reproductive success and ovariole number were mapped using a composite interval mapping procedure. QTL for female reproductive success were located at the tip of the X chromosome between markers at cytological locations 1B and 3E; and on the left arm of chromosome 2 in the 30D-38A cytological region. Ovariole number QTL mapped to cytological intervals 62D-69D and 98A-98E, both on the third chromosome. The regions harbouring QTL for female reproductive success and ovariole number were also identified as QTL for longevity in previous studies with these lines.

Are the same genes responsible for intra- and interspecific variability for sex comb tooth number in Drosophila?

The identification of genes contributing to speciation has the potential to provide important insights into the mechanisms of evolution. One of the most interesting unresolved puzzles is the relationship between intraspecific variability in morphological traits and their interspecific divergence. Intraspecific polymorphisms without major detrimental side-effects may serve as a substrate for selection response during speciation. The same quantitative trait loci (QTLs) may, then, account for the trait variability both within and between species. In contrast, the vast majority of intraspecific variants could be deleterious mutations that have not yet been selected out. In this case intraspecific variation would not ultimately generate interspecific trait differences. In previous work, QTLs responsible for morphological differences between animal species, including those for the secondary sexual trait sex comb tooth number, have been mapped with the resolution of chromosome segments. Here, we mapped QTLs for which alleles segregated between two laboratory lines of Drosophila melanogaster. The two QTLs identified mapped to the X chromosome and accounted for only 8% of the between-line differences, implying that a large number of small-effect genes modify sex combs. One intraspecific QTL mapped to the same interval as the QTL for interspecific differences between D. simulans and D. mauritiana. Whether or not these effects result from the same genes requires further examination.

The genetic architecture of selection response. Inferences from fine-scale mapping of bristle number quantitative trait loci in Drosophila melanogaster.

Quantitative trait loci (QTL) affecting responses and correlated responses to selection for abdominal and sternopleural bristle number have been mapped with high resolution to the X and third chromosomes. Advanced intercross recombinant isogenic chromosomes were constructed from high and low selection lines in an unselected inbred background, and QTL were detected using composite interval mapping and high density transposable element marker maps. We mapped a total of 26 bristle number QTL with large effects, which were in or immediately adjacent to intervals previously inferred to contain bristle number QTL on these chromosomes. The QTL contributing to response to selection for high bristle number were not the same as those contributing to response to selection for low bristle number, suggesting that distributions of allelic effects per locus may be asymmetrical. Correlated responses were more often attributable to loose linkage than pleiotropy or close linkage. Bristle number QTL mapping to the same locations have been inferred in studies with different parental strains. Of the 26 QTL, 20 mapped to locations consistent with candidate genes affecting peripheral nervous system development and/or bristle number. This facilitates determining the molecular basis of quantitative variation and allele frequencies by associating molecular variation at the candidate genes with phenotypic variation in bristle number in samples of alleles from nature.

High-resolution mapping of quantitative trait loci for sternopleural bristle number in Drosophila melanogaster.

We have mapped quantitative trait loci (QTL) harboring naturally occurring allelic variation for Drosophila bristle number. Lines with high (H) and low (L) sternopleural bristle number were derived by artificial selection from a large base population. Isogenic H and L sublines were extracted from the selection lines, and populations of X and third chromosome H/L recombinant isogenic lines were constructed in the homozygous low line background. The polymorphic cytological locations of roo transposable elements provided a dense molecular marker map with an average intermarker distance of 4.5 cM. Two X chromosome and six chromosome 3 QTL affecting response to selection for sternopleural bristle number and three X chromosome and three chromosome 3 QTL affecting correlated response in abdominal bristle number were detected using a composite interval mapping method. The average effects of bristle number QTL were moderately large, and some had sex-specific effects. Epistasis between QTL affecting sternopleural bristle number was common, and interaction effects were large. Many of the intervals containing bristle number QTL coincided with those mapped in previous studies. However, resolution of bristle number QTL to the level of genetic loci is not trivial, because the genomic regions containing bristle number QTL often did not contain obvious candidate loci, and results of quantitative complementation tests to mutations at candidate loci affecting adult bristle number were ambiguous.

New estimates of the rates and effects of mildly deleterious mutation in Drosophila melanogaster.

The genomic rate and distribution of effects of deleterious mutations are important parameters in evolutionary theory. The most detailed information comes from the work of Mukai and Ohnishi, who allowed mutations to accumulate on Drosophila melanogaster second chromosomes, shielded from selection and recombination by being maintained heterozygous in males. Averaged over studies, the estimated rate of nonlethal viability mutations per second chromosome per generation under an equal-effects model, UBM, was 0. 12, suggesting a high genomic mutation rate. We have performed a mutation-accumulation experiment similar to those of Mukai and Ohnishi, except that three large homozygous control populations were maintained. Egg-to-adult viability of 72 nonlethal mutation-accumulation (MA) lines and the controls was assayed after 27-33 generations of mutation accumulation. The rate of decline in mean viability was significantly lower than observed by Mukai, and the rate of increase in among-line variance was significantly higher. Our UBM estimate of 0.02 is much lower than the previous estimates. Our results suggest that the rate of mutations that detectably reduce viability may not be much greater than the lethal mutation rate (0.01 in these lines), but the results also are consistent with models that include many mutations with very small effects.

Sure facts, speculations, and open questions about the evolution of transposable element copy number.

Transposable elements (TEs) are sequences capable of multiplying in their host's genome. They survive by increasing copy numbers due to transpositions, and natural selection washes them out because hosts with heavier loads of TEs have lower fitness. The available phylogenetic evidence supports the view that TEs have existed in living organisms for hundreds of millions of years. A fundamental question facing the field is how can an equilibrium be attained between transposition and selection which allows these parasitic genetic elements to persist for such a long time period? To answer this question, it is necessary to understand how the rate of TE transposition is controlled and to describe the mechanisms with which natural selection opposes TE accumulation. Perhaps the best models for such a study are copia and gypsy retrotransposons in Drosophila. Their average rate of transposition in nature is between 10(-5) - 10(-4) transpositions per copy per generation. Unlike nature, transposition rates vary widely, from zero to 10(-2), between laboratory lines. This variability in transposition rate is controlled by host genes. It is probable that in nature TE site heterogeneity is caused by frequent transpositions in rare flies with permissive alleles, and no transpositions happen in the rest of flies. The average rate of TE transposition in nature may thus depend on the frequency of permissive alleles, which is a function of the rate of mutation from restrictive to permissive alleles, the mechanism and the strength of selection opposing TE multiplication, and population size. Thus, evolution of the frequency of permissive alleles of genes controlling transposition must be accounted for to understand evolution of TE copy numbers.

The relationship between the rate of transposition and transposable element copy number for copia and Doc retrotransposons of Drosophila melanogaster.

We present data on the relationship between the rate of transposition and copy number in the genome for the copia and Doc retrotransposons of Drosophila melanogaster. copia and Doc transposition rates were directly measured in sublines of the isogenic 2b line using individual males or females, respectively, with a range of copia copy numbers from 49 to 103 and Doc copy numbers from 112 to 235 per genome. Transposition rates varied from 3 x 10(-4) to 2 x 10(-2) for copia and from 2 x 10(-4) to 2 x 10(-3) for Doc. A positive relationship between transposition rate and copy number was found both for copia and for Doc when the data were analysed across all the 2b individuals; no significant correlation was found when the data were analysed across the subline means for both retrotransposons tested. Overall, correlation between copia and Doc transposition rate and their copy number in the genome, if any, was not negative, which would be expected if transposable elements (TEs) self-regulate their copy number. Thus, for copia and Doc no evidence for self-regulation was provided, and at least for these two TEs this hypothesis is not favoured for explaining the maintenance of the stable copy number that is characteristic for natural populations. The transposition rate of copia was measured twice, and a strong positive correlation between copy number and transposition rate both across individuals and subline means was found in 1994, while in 1995 no correlation was found. This fact is in agreement with the hypothesis that a positive correlation between the rate of transposition and TE copy number may be a default starting point for future host-TE coevolution.

Quantitative genetic analysis of copia retrotransposon activity in inbred Drosophila melanogaster lines.

The rates of transcription and transposition of retrotransposons vary between lines of Drosophila melanogaster. We have studied the genetics of differences in copia retrotransposon activity by quantitative trait loci (QTL) mapping. Ninety-eight recombinant inbred lines were constructed from two parental lines exhibiting a 10-fold difference in copia transcript level and a 100-fold difference in transposition rate. The lines were scored for 126 molecular markers, copia transcript level, and rate of copia transposition. Transcript level correlated with copia copy number, and the difference in copia copy number between parental lines accounted for 45.1% of copia transcript-level difference. Most of the remaining difference was accounted for by two transcript-level QTL mapping to cytological positions 27B-30D and 50F-57C on the second chromosome, which accounted for 11.5 and 30.4%, respectively. copia transposition rate was controlled by interacting QTL mapping to the region 27B-48D on the second and 61A-65A and 97D-100A on the third chromosome. The genes controlling copia transcript level are thus not necessarily those involved in controlling copia transposition rate. Segregation of modifying genes, rather than mutations, might explain the variability in copia retrotransposon activity between lines.