The genetic architecture of quantitative traits: lessons from Drosophila

Understanding the genetic architecture of quantitative traits begins with identifying the genes regulating these traits, mapping the subset of genetically varying quantitative trait loci (QTLs) in natural populations, and pinpointing the molecular polymorphisms defining QTL alleles. Studies in Drosophila have revealed large numbers of pleiotropic genes that interact epistatically to regulate quantitative traits, and large numbers of QTLs with sex-, environment- and genotype-specific effects. Multiple molecular polymorphisms in regulatory regions of candidate genes are often associated with variation for complex traits. These observations offer valuable lessons for understanding the genetic basis of variation for complex traits in other organisms, including humans.     

Complex genetic architecture of population differences in adult lifespan of a beetle: nonadditive inheritance, gender differences, body size and a large maternal effect

Evolutionary responses to selection can be complicated when there is substantial nonadditivity, which limits our ability to extrapolate from simple models of selection to population differentiation and speciation. Studies of Drosophila melanogaster indicate that lifespan and the rate of senescence are influenced by many genes that have environment- and sex-specific effects. These studies also demonstrate that interactions among alleles (dominance) and loci (epistasis) are common, with the degree of interaction differing between the sexes and among environments. However, little is known about the genetic architecture of lifespan or mortality rates for organisms other than D. melanogaster. We studied genetic architecture of differences in lifespan and shapes of mortality curves between two populations of the seed beetle, Callosobruchus maculatus (South India and Burkina Faso populations). These two populations differ in various traits (such as body size and adult lifespan) that have likely evolved via host-specific selection. We found that the genetic architecture of lifespan differences between populations differs substantially between males and females; there was a large maternal effect on male lifespan (but not on female lifespan), and substantial dominance of long-life alleles in females (but not males). The large maternal effect in males was genetically based (there was no significant cytoplasmic effect) likely due to population differences in maternal effects genes that influence lifespan of progeny. Rearing host did not affect the genetic architecture of lifespan, and there was no evidence that genes on the Y-chromosome influence the population differences in lifespan. Epistatic interactions among loci were detectable for the mortality rate of both males and females, but were detectable for lifespan only after controlling for body size variation among lines. The detection of epistasis, dominance, and sex-specific genetic effects on C. maculatus lifespan is consistent with results from line cross and quantitative trait locus studies of D. melanogaster.     

Genetic Load in Natural Populations: Is It Compatible with the Hypothesis That Many Polymorphisms Are Maintained by Natural Selection?

Recent studies of genetically controlled enzyme variation lead to an estimation that at least 30 to 60% of the structural genes are polymorphic in natural populations of many vertebrate and invertebrate species. Some authors have argued that a substantial proportion of these polymorphisms cannot be maintained by natural selection because this would result in an unbearable genetic load. If many polymorphisms are maintained by heterotic natural selection, individuals with much greater than average proportion of homozygous loci should have very low fitness. We have measured in Drosophila melanogaster the fitness of flies homozygous for a complete chromosome relative to normal wild flies. A total of 37 chromosomes from a natural population have been tested using 92 experimental populations. The mean fitness of homozygous flies is 0.12 for second chromosomes, and 0.13 for third chromosomes. These estimates are compatible with the hypothesis that many (more than one thousand) loci are maintained by heterotic selection in natural populations of D. melanogaster.     

Phenotypic variation and natural selection at catsup, a pleiotropic quantitative trait gene in Drosophila

Quantitative traits are shaped by networks of pleiotropic genes . To understand the mechanisms that maintain genetic variation for quantitative traits in natural populations and to predict responses to artificial and natural selection, we must evaluate pleiotropic effects of underlying quantitative trait genes and define functional allelic variation at the level of quantitative trait nucleotides (QTNs). Catecholamines up (Catsup), which encodes a negative regulator of tyrosine hydroxylase , the rate-limiting step in the synthesis of the neurotransmitter dopamine, is a pleiotropic quantitative trait gene in Drosophila melanogaster. We used association mapping to determine whether the same or different QTNs at Catsup are associated with naturally occurring variation in multiple quantitative traits. We sequenced 169 Catsup alleles from a single population and detected 33 polymorphisms with little linkage disequilibrium (LD). Different molecular polymorphisms in Catsup are independently associated with variation in longevity, locomotor behavior, and sensory bristle number. Most of these polymorphisms are potentially functional variants in protein coding regions, have large effects, and are not common. Thus, Catsup is a pleiotropic quantitative trait gene, but individual QTNs do not have pleiotropic effects. Molecular population genetic analyses of Catsup sequences are consistent with balancing selection maintaining multiple functional polymorphisms.     

Genetic variation affecting host-parasite interactions: major-effect quantitative trait loci affect the transmission of sigma virus in Drosophila melanogaster

In natural populations, genetic variation affects resistance to disease. Whether that genetic variation comprises lots of small-effect polymorphisms or a small number of large-effect polymorphisms has implications for adaptation, selection and how genetic variation is maintained in populations. Furthermore, how much genetic variation there is, and the genes that underlie this variation, affects models of co-evolution between parasites and their hosts. We are studying the genetic variation that affects the resistance of Drosophila melanogaster to its natural pathogen — the vertically transmitted sigma virus. We have carried out three separate quantitative trait locus mapping analyses to map gene variants on the second chromosome that cause variation in the rate at which males transmit the infection to their offspring. All three crosses identified a locus in a similar chromosomal location that causes a large drop in the rate at which the virus is transmitted. We also found evidence for an additional smaller-effect quantitative trait locus elsewhere on the chromosome. Our data, together with previous experiments on the sigma virus and parasitoid wasps, indicate that the resistance of D. melanogaster to co-evolved pathogens is controlled by a limited number of major-effect polymorphisms.     

Genetic basis of adaptation in Arabidopsis thaliana: local adaptation at the seed dormancy QTL DOG1

Local adaptation provides an opportunity to study the genetic basis of adaptation and investigate the allelic architecture of adaptive genes. We study delay of germination 1 (DOG1), a gene controlling natural variation in seed dormancy in Arabidopsis thaliana and investigate evolution of dormancy in 41 populations distributed in four regions separated by natural barriers. Using F(ST) and Q(ST) comparisons, we compare variation at DOG1 with neutral markers and quantitative variation in seed dormancy. Patterns of genetic differentiation among populations suggest that the gene DOG1 contributes to local adaptation. Although Q(ST) for seed dormancy is not different from F(ST) for neutral markers, a correlation with variation in summer precipitation supports that seed dormancy is adaptive. We characterize dormancy variation in several F(2) -populations and show that a series of functionally distinct alleles segregate at the DOG1 locus. Theoretical models have shown that the number and effect of alleles segregatin at quantitative trait loci (QTL) have important consequences for adaptation. Our results provide support to models postulating a large number of alleles at quantitative trait loci involved in adaptation.     

Investigating natural variation in Drosophila courtship song by the evolve and resequence approach

A primary goal of population genetics is to determine the genetic basis of natural trait variation. We could significantly advance this goal by developing comprehensive genome-wide approaches to link genotype and phenotype in model organisms. Here we combine artificial selection with population-based resequencing to investigate the genetic basis of variation in the interpulse interval (IPI) of Drosophila melanogaster courtship song. We performed divergent selection on replicate populations for only 14 generations, but had considerable power to differentiate alleles that evolved due to selection from those that evolved stochastically. We identified a large number of variants that changed frequency in response to selection for this simple behavior, and they are highly underrepresented on the X chromosome. Though our power was adequate using this experimental technique, the ability to differentiate causal variants from those affected by linked selection requires further development.     

Four loci explain 83% of size variation in the horse

Horse body size varies greatly due to intense selection within each breed. American Miniatures are less than one meter tall at the withers while Shires and Percherons can exceed two meters. The genetic basis for this variation is not known. We hypothesize that the breed population structure of the horse should simplify efforts to identify genes controlling size. In support of this, here we show with genome-wide association scans (GWAS) that genetic variation at just four loci can explain the great majority of horse size variation. Unlike humans, which are naturally reproducing and possess many genetic variants with weak effects on size, we show that horses, like other domestic mammals, carry just a small number of size loci with alleles of large effect. Furthermore, three of our horse size loci contain the LCORL, HMGA2 and ZFAT genes that have previously been found to control human height. The LCORL/NCAPG locus is also implicated in cattle growth and HMGA2 is associated with dog size. Extreme size diversification is a hallmark of domestication. Our results in the horse, complemented by the prior work in cattle and dog, serve to pinpoint those very few genes that have played major roles in the rapid evolution of size during domestication.     

Climatic selection on genes and traits after a 100 year-old invasion: a critical look at the temperate-tropical clines in Drosophila melanogaster from eastern Australia

Drosophila melanogaster invaded Australia around 100 years ago, most likely through a northern invasion. The wide range of climatic conditions in eastern Australia across which D. melanogaster is now found provides an opportunity for researchers to identify traits and genes that are associated with climatic adaptation. Allozyme studies indicate clinal patterns for at least four loci including a strong linear cline in Adh and a non-linear cline in alpha-Gpdh. Inversion clines were initially established from cytological studies but have now been validated with larger sample sizes using molecular markers for breakpoints. Recent collections indicate that some genetic markers (Adh and In(3R)Payne) have changed over the last 20 years reflecting continuing evolution. Heritable clines have been established for quantitative traits including wing length/area, thorax length and cold and heat resistance. A cline in egg size independent of body size and a weak cline in competitive ability have also been established. Postulated clinal patterns for resistance to desiccation and starvation have not been supported by extensive sampling. Experiments under laboratory and semi-natural conditions have suggested selective factors generating clinal patterns, particularly for reproductive patterns over winter. Attempts are being made to link clinal variation in traits to specific genes using QTL analysis and the candidate locus approach, and to identify the genetic architecture of trait variation along the cline. This is proving difficult because of inversion polymorphisms that generate disequilibrium among genes. Substantial gaps still remain in linking clines to field selection and understanding the genetic and physiological basis of the adaptive shifts. However D. melanogaster populations in eastern Australia remain an excellent resource for understanding past and future evolutionary responses to climate change.     

Beyond pathways: genetic dissection of tocopherol content in maize kernels by combining linkage and association analyses

Although tocopherols play an important role in plants and animals, the genetic architecture of tocopherol content in maize kernels has remained largely unknown. In this study, linkage and association analyses were conducted to examine the genetic architecture of tocopherol content in maize kernels. Forty‐one unique quantitative trait loci (QTLs) were identified by linkage mapping in six populations of recombinant inbred lines (RILs). In addition, 32 significant loci were detected via genome‐wide association study (GWAS), 18 of which colocalized with the QTLs identified by linkage mapping. Fine mapping of a major QTL validated the accuracy of GWAS and QTL mapping results and suggested a role for nontocopherol pathway genes in the modulation of natural tocopherol variation. We provided genome‐wide evidence that genes involved in fatty acid metabolism, chlorophyll metabolism and chloroplast function may affect natural variation in tocopherols. These findings were confirmed through mutant analysis of a particular gene from the fatty acid pathway. In addition, the favourable alleles for many of the significant SNPs/QTLs represented rare alleles in natural populations. Together, our results revealed many novel genes that are potentially involved in the variation of tocopherol content in maize kernels. Pyramiding of the favourable alleles of the newly elucidated genes and the well‐known tocopherol pathway genes would greatly improve tocopherol content in maize.     

A Reanalysis of Protein Polymorphism in Drosophila Melanogaster, D. Simulans, D. Sechellia and D. Mauritiana: Effects of Population Size and Selection

Comparison of synonymous and nonsynonymous variation/substitution within and between species at individual genes has become a widely used general approach to detect the effect of selection versus drift. The sibling species group comprised of two cosmopolitan (Drosophila melanogaster and Drosophila simulans) and two island (Drosophila mauritiana and Drosophila sechellia) species has become a model system for such studies. In the present study we reanalyzed the pattern of protein variation in these species, and the results were compared against the patterns of nucleotide variation obtained from the literature, mostly available for melanogaster and simulans. We have mainly focused on the contrasting patterns of variation between the cosmopolitan pair. The results can be summarized as follows: (1) As expected the island species D. mauritiana and D. sechellia showed much less variation than the cosmopolitan species D. melanogaster and D. simulans. (2) The chromosome 2 showed significantly less variation than chromosome 3 and X in all four species which may indicate effects of past selective sweeps. (3) In contrast to its overall low variation, D. mauritiana showed highest variation for X-linked loci which may indicate introgression from its sibling, D. simulans. (4) An average population of D. simulans was as heterozygous as that of D. melanogaster (14.4% v.s. 13.9%) but the difference was large and significant when considering only polymorphic loci (37.2% v.s. 26.1%). (5) The species-wise pooled populations of these two species showed similar results (all loci = 18.3% v.s. 20.0%, polymorphic loci = 47.2% v.s. 37.6%). (6) An average population of D. simulans had more low-frequency alleles than D. melanogaster, and the D. simulans alleles were found widely distributed in all populations whereas the D. melanogaster alleles were limited to local populations. As a results of this, pooled populations of D. melanogaster showed more polymorphic loci than those of D. simulans (48.0% v.s. 32.0%) but the difference was reduced when the comparison was made on the basis of an average population (29.1% v.s. 21.4%). (7) While the allele frequency distributions within populations were nonsignificant in both D. melanogaster and D. simulans, melanogaster had fewer than simulans, but more than expected from the neutral theory, low frequency alleles. (8) Diallelic loci with the second allele with a frequency less than 20% had similar frequencies in all four species but those with the second allele with a frequency higher than 20% were limited to only melanogaster the latter group of loci have clinal (latitudinal) patterns of variation indicative of balancing selection. (9) The comparison of D. simulans/D. melanogaster protein variation gave a ratio of 1.04 for all loci and 1.42 for polymorphic loci, against a ratio of approximately 2-fold difference for silent nucleotide sites. This suggests that the species ratios of protein and silent nucleotide polymorphism are too close to call for selective difference between silent and allozyme variation in D. simulans. In conclusion, the contrasting levels of allozyme polymorphism, distribution of rare alleles, number of diallelic loci and the patterns of geographic differentiation between the two species suggest the role of natural selection in D. melanogaster, and of possibly ancient population structure and recent worldwide migration in D. simulans. Population size differences alone are insufficient as an explanation for the patterns of variation between these two species.     

The genetic architecture of resistance to virus infection in Drosophila

Variation in susceptibility to infection has a substantial genetic component in natural populations, and it has been argued that selection by pathogens may result in it having a simpler genetic architecture than many other quantitative traits. This is important as models of host–pathogen co‐evolution typically assume resistance is controlled by a small number of genes. Using the Drosophila melanogaster multiparent advanced intercross, we investigated the genetic architecture of resistance to two naturally occurring viruses, the sigma virus and DCV (Drosophila C virus). We found extensive genetic variation in resistance to both viruses. For DCV resistance, this variation is largely caused by two major‐effect loci. Sigma virus resistance involves more genes – we mapped five loci, and together these explained less than half the genetic variance. Nonetheless, several of these had a large effect on resistance. Models of co‐evolution typically assume strong epistatic interactions between polymorphisms controlling resistance, but we were only able to detect one locus that altered the effect of the main effect loci we had mapped. Most of the loci we mapped were probably at an intermediate frequency in natural populations. Overall, our results are consistent with major‐effect genes commonly affecting susceptibility to infectious diseases, with DCV resistance being a near‐Mendelian trait.     

A Comprehensive Study of Genic Variation in Natural Populations of Drosophila melanogaster. I. Estimates of Gene Flow from Rare Alleles

In order to assess the evolutionary significance of molecular variation in natural populations of Drosophila melanogaster, we have started a comprehensive genetic variation study program employing a relatively large number of gene-protein loci and an array of populations obtained from various geographic locations throughout the world. In this first report we provide estimates of gene flow based on the spatial distributions of rare alleles at 117 gene loci in 15 worldwide populations of D. melanogaster . Estimates of Nm (number of migrants exchanged per generation among populations) range from 1.09 in East-Asian populations (Taiwan, Vietnam and Australia) to 2.66 in West-Coast populations of North America. These estimates, among geographic populations separated by hundreds or even thousands of miles, suggest that gene flow among neighboring populations of D. melanogaster is quite extensive. This means that, for selectively neutral genes, we should expect little differentiation among neighboring populations. A survey of eight West-Coast populations of D. melanogaster (geographically comparable to Drosophila pseudoobscura) showed that in spite of extensive gene flow, populations of D. melanogaster show much more geographic differentiation than comparable populations of D. pseudoobscura. From this we conclude that migration in combination with natural selection rather than migration alone is responsible for the geographic uniformity of molecular polymorphisms in D. pseudoobscura.     

QTL mapping reveals a striking coincidence in the positions of genomic regions associated with adaptive variation in body size in parallel clines of Drosophila melanogaster on different continents

Latitudinal genetic clines in body size are common in many ectotherm species and are attributed to climatic adaptation. Here, we use Quantitative Trait Loci (QTL) mapping to identify genomic regions associated with adaptive variation in body size in natural populations of Drosophila melanogaster from extreme ends of a cline in South America. Our results show that there is a significant association between the positions of QTL with strong effects on wing area in South America and those previously reported in a QTL mapping study of Australian cline end populations (P < 0.05). In both continents, the right arm of the third chromosome is associated with QTL with the strongest effect on wing area. We also show that QTL peaks for wing area and thorax length are associated with the same genomic regions, indicating that the clinal variation in the body size traits may have a similar genetic basis. The consistency of the results found for the South American and Australian cline end populations indicate that the genetic basis of the two clines may be similar and future efforts to identify the genes producing the response to selection should be focused on the genomic regions highlighted by the present work.     

Whole‐genome resequencing of extreme phenotypes in collared flycatchers highlights the difficulty of detecting quantitative trait loci in natural populations

Dissecting the genetic basis of phenotypic variation in natural populations is a long‐standing goal in evolutionary biology. One open question is whether quantitative traits are determined only by large numbers of genes with small effects, or whether variation also exists in large‐effect loci. We conducted genomewide association analyses of forehead patch size (a sexually selected trait) on 81 whole‐genome‐resequenced male collared flycatchers with extreme phenotypes, and on 415 males sampled independent of patch size and genotyped with a 50K SNP chip. No SNPs were genomewide statistically significantly associated with patch size. Simulation‐based power analyses suggest that the power to detect large‐effect loci responsible for 10% of phenotypic variance was <0.5 in the genome resequencing analysis, and <0.1 in the SNP chip analysis. Reducing the recombination by two‐thirds relative to collared flycatchers modestly increased power. Tripling sample size increased power to >0.8 for resequencing of extreme phenotypes (N = 243), but power remained <0.2 for the 50K SNP chip analysis (N = 1245). At least 1 million SNPs were necessary to achieve power >0.8 when analysing 415 randomly sampled phenotypes. However, power of the 50K SNP chip to detect large‐effect loci was nearly 0.8 in simulations with a small effective population size of 1500. These results suggest that reliably detecting large‐effect trait loci in large natural populations will often require thousands of individuals and near complete sampling of the genome. Encouragingly, far fewer individuals and loci will often be sufficient to reliably detect large‐effect loci in small populations with widespread strong linkage disequilibrium.     

Quantitative trait loci for locomotor behavior in Drosophila melanogaster

Locomotion is an integral component of most animal behaviors and many human diseases and disorders are associated with locomotor deficits, but little is known about the genetic basis of natural variation in locomotor behavior. Locomotion is a complex trait, with variation attributable to the joint segregation of multiple interacting quantitative trait loci (QTL), with effects that are sensitive to the environment. We assessed variation in a component of locomotor behavior (locomotor reactivity) in a population of 98 recombinant inbred lines of Drosophila melanogaster and mapped four QTL affecting locomotor reactivity by linkage to polymorphic roo transposable element insertion sites. We used complementation tests of deficiencies to fine map these QTL to 12 chromosomal regions and complementation tests of mutations to identify 13 positional candidate genes affecting locomotor reactivity, including Dopa decarboxylase (Ddc), which catalyzes the final step in the synthesis of serotonin and dopamine. Linkage disequilibrium mapping in a population of 164 second chromosome substitution lines derived from a single natural population showed that polymorphisms at Ddc were associated with naturally occurring genetic variation in locomotor behavior. These data implicate variation in the synthesis of bioamines as a factor contributing to natural variation in locomotor reactivity.     

Enzyme and Chromosome Polymorphisms in Japanese Natural Populations of DROSOPHILA MELANOGASTER

Collections of D. melanogaster from Japanese populations were analyzed for enzyme and chromosomal polymorphisms. Allelic frequencies at the Adh and alphaGpd loci were compared with polymorphic inversion (In(2L)B, In(2R)C) frequencies in the second chromosome. There was a significant positive correlation between the frequencies of AdhS and In(2L)B, caused by linkage. On the other hand, inversion-free cage populations maintained in the laboratory for a long time showed considerably larger variation in the frequencies of these enzyme alleles, which seem very likely to be a consequence of random drift. Two fitness components of these enzyme and chromosomal variants were measured in two different environmental conditions; neither of the two loci showed heterozygote superiority in viability or productivity, while the inversion heterozygotes showed a superior productivity compared to the corresponding homozygotes in the fluctuating environment. These findings are compatible with the hypothesis that polymorphic isozyme genes are maintained by random drift of neutral genes in natural populations, and that association with linked inversions is a historical accident.     

The adaptive genomic landscape of beak morphology in Darwin's finches

Beak shape in Darwin's ground finches (Geospiza) is emblematic of natural selection and adaptive radiation, yet our understanding of the genetic basis of beak shape variation, and thus the genetic target of natural selection, is still evolving. Here we reveal the genomic architecture of beak shape variation using genomewide comparisons of four closely related and hybridizing species across 13 islands subject to parallel natural selection. Pairwise contrasts among species were used to identify a large number of genomic loci that are consistently related to species differences across a complex landscape. These loci are associated with hundreds of genes that have enriched GO categories significantly associated with development. One genomic region of particular interest is a section of Chromosome 1A with many candidate genes and increased linkage. The distinct, pointed beak shape of the cactus finch is linked to an excess of intermediate frequency alleles and increased heterozygosity in significant SNPs, but not across the rest of the genome. Alleles associated with pointier beaks among species were associated with pointier‐beaked populations within each species, thus establishing a common basis for natural selection, species divergence and adaptive radiation. The adaptive genomic landscape for Darwin's finches mirrors theoretical expectations based on morphological variation. The implication that a large number of genes are actively maintained to facilitate beak variation across parallel populations with documented interspecies admixture challenges our understanding of evolutionary processes in the wild.