Genotype-environment interaction at quantitative trait loci affecting sensory bristle number in Drosophila melanogaster

Genomics of Developmental Plasticity in Animals

Developmental plasticity refers to the property by which the same genotype produces distinct phenotypes depending on the environmental conditions under which development takes place. By allowing organisms to produce phenotypes adjusted to the conditions that adults will experience, developmental plasticity can provide the means to cope with environmental heterogeneity. Developmental plasticity can be adaptive and its evolution can be shaped by natural selection. It has also been suggested that developmental plasticity can facilitate adaptation and promote diversification. Here, we summarize current knowledge on the evolution of plasticity and on the impact of plasticity on adaptive evolution, and we identify recent advances and important open questions about the genomics of developmental plasticity in animals. We give special attention to studies using transcriptomics to identify genes whose expression changes across developmental environments and studies using genetic mapping to identify loci that contribute to variation in plasticity and can fuel its evolution.

Dissection of Complex, Fitness-Related Traits in Multiple Drosophila Mapping Populations Offers Insight into the Genetic Control of Stress Resistance

We leverage two complementary Drosophila melanogaster mapping panels to genetically dissect starvation resistance-an important fitness trait. Using >1600 genotypes from the multiparental Drosophila Synthetic Population Resource (DSPR), we map numerous starvation stress QTL that collectively explain a substantial fraction of trait heritability. Mapped QTL effects allowed us to estimate DSPR founder phenotypes, predictions that were correlated with the actual phenotypes of these lines. We observe a modest phenotypic correlation between starvation resistance and triglyceride level, traits that have been linked in previous studies. However, overlap among QTL identified for each trait is low. Since we also show that DSPR strains with extreme starvation phenotypes differ in desiccation resistance and activity level, our data imply multiple physiological mechanisms contribute to starvation variability. We additionally exploited the Drosophila Genetic Reference Panel (DGRP) to identify sequence variants associated with starvation resistance. Consistent with prior work these sites rarely fall within QTL intervals mapped in the DSPR. We were offered a unique opportunity to directly compare association mapping results across laboratories since two other groups previously measured starvation resistance in the DGRP. We found strong phenotypic correlations among studies, but extremely low overlap in the sets of genomewide significant sites. Despite this, our analyses revealed that the most highly associated variants from each study typically showed the same additive effect sign in independent studies, in contrast to otherwise equivalent sets of random variants. This consistency provides evidence for reproducible trait-associated sites in a widely used mapping panel, and highlights the polygenic nature of starvation resistance.

Genetic basis of thermal plasticity variation in Drosophila melanogaster body size

Body size is a quantitative trait that is closely associated to fitness and under the control of both genetic and environmental factors. While developmental plasticity for this and other traits is heritable and under selection, little is known about the genetic basis for variation in plasticity that can provide the raw material for its evolution. We quantified genetic variation for body size plasticity in Drosophila melanogaster by measuring thorax and abdomen length of females reared at two temperatures from a panel representing naturally segregating alleles, the Drosophila Genetic Reference Panel (DGRP). We found variation between genotypes for the levels and direction of thermal plasticity in size of both body parts. We then used a Genome-Wide Association Study (GWAS) approach to unravel the genetic basis of inter-genotype variation in body size plasticity, and used different approaches to validate selected QTLs and to explore potential pleiotropic effects. We found mostly “private QTLs”, with little overlap between the candidate loci underlying variation in plasticity for thorax versus abdomen size, for different properties of the plastic response, and for size versus size plasticity. We also found that the putative functions of plasticity QTLs were diverse and that alleles for higher plasticity were found at lower frequencies in the target population. Importantly, a number of our plasticity QTLs have been targets of selection in other populations. Our data sheds light onto the genetic basis of inter-genotype variation in size plasticity that is necessary for its evolution.

Environmental conditions can influence development and lead to the production of phenotypes adjusted to the conditions adults will live in. This developmental plasticity, which can help organisms cope with environmental heterogeneity, is heritable and under selection. Its evolution will depend on available genetic variation. Using a panel of D. melanogaster flies representing naturally segregating alleles, we identified DNA sequence variants associated to variation in thermal plasticity for body size. We found that these variants correspond to a diverse set of functions and that their effects differ between body parts and properties of the thermal response. Our results shed new light onto the long discussed genes for plasticity.

Beyond mean allelic effects: A locus at the major color gene MC1R associates also with differing levels of phenotypic and genetic (co)variance for coloration in barn owls

The mean phenotypic effects of a discovered variant help to predict major aspects of the evolution and inheritance of a phenotype. However, differences in the phenotypic variance associated to distinct genotypes are often overlooked despite being suggestive of processes that largely influence phenotypic evolution, such as interactions between the genotypes with the environment or the genetic background. We present empirical evidence for a mutation at the melanocortin-1-receptor gene, a major vertebrate coloration gene, affecting phenotypic variance in the barn owl, Tyto alba. The white MC1R allele, which associates with whiter plumage coloration, also associates with a pronounced phenotypic and additive genetic variance for distinct color traits. Contrarily, the rufous allele, associated with a rufous coloration, relates to a lower phenotypic and additive genetic variance, suggesting that this allele may be epistatic over other color loci. Variance differences between genotypes entailed differences in the strength of phenotypic and genetic associations between color traits, suggesting that differences in variance also alter the level of integration between traits. This study highlights that addressing variance differences of genotypes in wild populations provides interesting new insights into the evolutionary mechanisms and the genetic architecture underlying the phenotype.

Accumulation, elimination, sequestration, and genetic variation of lead (Pb2+) loads within and between generations of Drosophila melanogaster

We examined accumulation, sequestration, elimination, and genetic variation for lead (Pb) loads within and between generations of Drosophila melanogaster. Flies were reared in control or leaded medium at various doses and tested for their Pb loads at different stages of development (larvae, eclosion, newly-eclosed adults, and mature adults). Pb loads were tested using Inductively Coupled Plasma Mass Spectrometry (ICP-MS). We found that D. melanogaster readily accumulated Pb throughout their lifespan and the levels of accumulation increased with Pb exposure in the medium. Wandering third-instar larvae accumulated more Pb than mature adults; this phenomenon may be due to elimination of Pb in the pupal cases during eclosion and/or depuration in adults post-eclosion. The accumulated Pb in mature adults was not transferred to F1 mature adult offspring. Using a set of recombinant inbred strains, we identified a quantitative trait locus for adult Pb loads and found that genetic variation accounted for 34% of the variance in Pb load. We concluded that D. melanogaster is a useful model organism for evaluating changes in Pb loads during development, as well as between generations. Furthermore, we found that genetic factors can influence Pb loads; this provides an essential foundation for evaluating phenotypic variation induced by the toxic effects of Pb.

Social effects for locomotion vary between environments in Drosophila melanogaster females

Despite strong purifying or directional selection, variation is ubiquitous in populations. One mechanism for the maintenance of variation is indirect genetic effects (IGEs), as the fitness of a given genotype will depend somewhat on the genes of its social partners. IGEs describe the effect of genes in social partners on the expression of the phenotype of a focal individual. Here, we ask what effect IGEs, and variation in IGEs between abiotic environments, has on locomotion in Drosophila. This trait is known to be subject to intralocus sexually antagonistic selection. We estimate the coefficient of interaction, Ψ, using six inbred lines of Drosophila. We found that Ψ varied between abiotic environments, and that it may vary across among male genotypes in an abiotic environment specific manner. We also found evidence that social effects of males alter the value of a sexually dimorphic trait in females, highlighting an interesting avenue for future research into sexual antagonism. We conclude that IGEs are an important component of social and sexual interactions and that they vary between individuals and abiotic environments in complex ways, with the potential to promote the maintenance of phenotypic variation.

Genotypic variation in biomass allocation in response to field drought has a greater affect on yield than gas exchange or phenology

BACKGROUND

Plant performance in agricultural and natural settings varies with moisture availability, and understanding the range of potential drought responses and the underlying genetic architecture is important for understanding how plants will respond to both natural and artificial selection in various water regimes. Here, we raised genotypes of Brassica rapa under well-watered and drought treatments in the field. Our primary goal was to understand the genetic architecture and yield effects of different drought-escape and dehydration-avoidance strategies.

RESULTS

Drought treatments reduced soil moisture by 62 % of field capacity. Drought decreased biomass accumulation and fruit production by as much as 48 %, whereas instantaneous water-use efficiency and root:shoot ratio increased. Genotypes differed in the mean value of all traits and in the sensitivity of biomass accumulation, root:shoot ratio, and fruit production to drought. Bivariate correlations involving gas-exchange and phenology were largely constant across environments, whereas those involving root:shoot varied across treatments. Although root:shoot was typically unrelated to gas-exchange or yield under well-watered conditions, genotypes with low to moderate increases in root:shoot allocation in response to drought survived the growing season, maintained maximum photosynthesis levels, and produced more fruit than genotypes with the greatest root allocation under drought. QTL for gas-exchange and yield components (total biomass or fruit production) had common effects across environments while those for root:shoot were often environment-specific.

CONCLUSIONS

Increases in root allocation beyond those needed to survive and maintain favorable water relations came at the cost of fruit production. The environment-specific effects of root:shoot ratio on yield and the differential expression of QTL for this trait across water regimes have important implications for efforts to improve crops for drought resistance.

Natural variability in Drosophila larval and pupal NaCl tolerance

The regulation of NaCl is essential for the maintenance of cellular tonicity and functionality, and excessive salt exposure has many adverse effects. The fruit fly, Drosophila melanogaster, is a good osmoregulator and some strains can survive on media with very low or high NaCl content. Previous analyses of mutant alleles have implicated various stress signaling cascades in NaCl sensitivity or tolerance; however, the genes influencing natural variability of NaCl tolerance remain for the most part unknown. Here, we use two approaches to investigate natural variation in D. melanogaster NaCl tolerance. We describe four D. melanogaster lines that were selected for different degrees of NaCl tolerance, and present data on their survival, development, and pupation position when raised on varying NaCl concentrations. After finding evidence for natural variation in salt tolerance, we present the results of Quantitative Trait Loci (QTL) mapping of natural variation in larval and pupal NaCl tolerance, and identify different genomic regions associated with NaCl tolerance during larval and pupal development.

Novel genetic capacitors and potentiators for the natural genetic variation of sensory bristles and their trait specificity in Drosophila melanogaster

Cryptic genetic variation (CGV) is defined as the genetic variation that has little effect on phenotypic variation under a normal condition, but contributes to heritable variation under environmental or genetic perturbations. Genetic buffering systems that suppress the expression of CGV and store it in a population are called genetic capacitors, and the opposite systems are called genetic potentiators. One of the best-known candidates for a genetic capacitor and potentiator is the molecular chaperone protein, HSP90, and one of its characteristics is that it affects the genetic variation in various morphological traits. However, it remains unclear whether the wide-ranging effects of HSP90 on a broad range of traits are a general feature of genetic capacitors and potentiators. In the current study, I searched for novel genetic capacitors and potentiators for quantitative bristle traits of Drosophila melanogaster and then investigated the trait specificity of their genetic buffering effect. Three bristle traits of D. melanogaster were used as the target traits, and the genomic regions with genetic buffering effects were screened using the 61 genomic deficiencies examined previously for genetic buffering effects in wing shape. As a result, four and six deficiencies with significant effects on increasing and decreasing the broad-sense heritability of the bristle traits were identified, respectively. Of the 18 deficiencies with significant effects detected in the current study and/or by the previous study, 14 showed trait-specific effects, and four affected the genetic buffering of both bristle traits and wing shape. This suggests that most genetic capacitors and potentiators exert trait-specific effects, but that general capacitors and potentiators with effects on multiple traits also exist.

QTL architecture of reproductive fitness characters in Brassica rapa

BACKGROUND

Reproductive output is critical to both agronomists seeking to increase seed yield and to evolutionary biologists interested in understanding natural selection. We examine the genetic architecture of diverse reproductive fitness traits in recombinant inbred lines (RILs) developed from a crop (seed oil) × wild-like (rapid cycling) genotype of Brassica rapa in field and greenhouse environments.

RESULTS

Several fitness traits showed strong correlations and QTL-colocalization across environments (days to bolting, fruit length and seed color). Total fruit number was uncorrelated across environments and most QTL affecting this trait were correspondingly environment-specific. Most fitness components were positively correlated, consistent with life-history theory that genotypic variation in resource acquisition masks tradeoffs. Finally, we detected evidence of transgenerational pleiotropy, that is, maternal days to bolting was negatively correlated with days to offspring germination. A QTL for this transgenerational correlation was mapped to a genomic region harboring one copy of FLOWERING LOCUS C, a genetic locus known to affect both days to flowering as well as germination phenotypes.

CONCLUSIONS

This study characterizes the genetic structure of important fitness/yield traits within and between generations in B. rapa. Several identified QTL are suitable candidates for fine-mapping for the improvement of yield in crop Brassicas. Specifically, brFLC1, warrants further investigation as a potential regulator of phenology between generations.

Quantitative trait loci × environment interactions for plant morphology vary over ontogeny in Brassica rapa

Growth in plants occurs via the addition of repeating modules, suggesting that the genetic architecture of similar subunits may vary between earlier- and later-developing modules. These complex environment × ontogeny interactions are not well elucidated, as studies examining quantitative trait loci (QTLs) expression over ontogeny have not included multiple environments. Here, we characterized the genetic architecture of vegetative traits and onset of reproduction over ontogeny in recombinant inbred lines of Brassica rapa in the field and glasshouse. The magnitude of genetic variation in plasticity of seedling internodes was greater than in those produced later in ontogeny. We correspondingly detected that QTLs for seedling internode length were environment-specific, whereas later in ontogeny the majority of QTLs affected internode lengths in all treatments. The relationship between internode traits and onset of reproduction varied with environment and ontogenetic stage. This relationship was observed only in the glasshouse environment and was largely attributable to one environment-specific QTL. Our results provide the first evidence of a QTL × environment × ontogeny interaction, and provide QTL resolution for differences between early- and later-stage plasticity for stem elongation. These results also suggest potential constraints on morphological evolution in early vs later modules as a result of associations with reproductive timing.

No boundaries: genomes, organisms, and ecological interactions responsible for divergence and reproductive isolation

Revealing the genetic basis of traits that cause reproductive isolation, particularly premating or sexual isolation, usually involves the same challenges as most attempts at genotype-phenotype mapping and so requires knowledge of how these traits are expressed in different individuals, populations, and environments, particularly under natural conditions. Genetic dissection of speciation phenotypes thus requires understanding of the internal and external contexts in which underlying genetic elements are expressed. Gene expression is a product of complex interacting factors internal and external to the organism including developmental programs, the genetic background including nuclear-cytotype interactions, epistatic relationships, interactions among individuals or social effects, stochasticity, and prevailing variation in ecological conditions. Understanding of genomic divergence associated with reproductive isolation will be facilitated by functional expression analysis of annotated genomes in organisms with well-studied evolutionary histories, phylogenetic affinities, and known patterns of ecological variation throughout their life cycles. I review progress and prospects for understanding the pervasive role of host plant use on genetic and phenotypic expression of reproductive isolating mechanisms in cactophilic Drosophila mojavensis and suggest how this system can be used as a model for revealing the genetic basis for species formation in organisms where speciation phenotypes are under the joint influences of genetic and environmental factors.

Deficiency screening for genomic regions with effects on environmental sensitivity of the sensory bristles of Drosophila melanogaster

Environmental canalization is defined as a reduction in the effect of external environmental perturbations on a phenotype, while phenotypic plasticity is defined as the production of different phenotypes in alternative environments. These terms describe different aspects of the same phenomenon, that is, the sensitivity of the phenotype to the environment. Genetic regulation of the environmental sensitivity has been a central topic in the field of evolutionary biology. In this study, we performed deficiency screening to detect genomic regions with effects on the environmental sensitivity of Drosophila melanogaster sensory bristles. We used a collection of isogenic deficiency strains established by the DrosDel Project for screening. We screened 423 genomic deficiencies that encompassed approximately 63.6% of the entire D. melanogaster genome. We identified 29 genomic deficiencies showing significant effects on environmental sensitivity, suggesting that multiple genomic regions may influence phenotypic variation. We also found significant correlations among the effects of deficiencies on environmental sensitivity for different bristle traits, suggesting that the same genetic mechanism can regulate environmental sensitivity of multiple traits. Current high-resolution mapping will facilitate the examination of individual candidate genes using mutations or RNAi approaches in future studies.

Quantitative variation in water-use efficiency across water regimes and its relationship with circadian, vegetative, reproductive, and leaf gas-exchange traits

Drought limits light harvesting, resulting in lower plant growth and reproduction. One trait important for plant drought response is water-use efficiency (WUE). We investigated (1) how the joint genetic architecture of WUE, reproductive characters, and vegetative traits changed across drought and well-watered conditions, (2) whether traits with distinct developmental bases (e.g. leaf gas exchange versus reproduction) differed in the environmental sensitivity of their genetic architecture, and (3) whether quantitative variation in circadian period was related to drought response in Brassica rapa. Overall, WUE increased in drought, primarily because stomatal conductance, and thus water loss, declined more than carbon fixation. Genotypes with the highest WUE in drought expressed the lowest WUE in well-watered conditions, and had the largest vegetative and floral organs in both treatments. Thus, large changes in WUE enabled some genotypes to approach vegetative and reproductive trait optima across environments. The genetic architecture differed for gas-exchange and vegetative traits across drought and well-watered conditions, but not for floral traits. Correlations between circadian and leaf gas-exchange traits were significant but did not vary across treatments, indicating that circadian period affects physiological function regardless of water availability. These results suggest that WUE is important for drought tolerance in Brassica rapa and that artificial selection for increased WUE in drought will not result in maladaptive expression of other traits that are correlated with WUE.

The genetics of phenotypic plasticity. X. Variation versus uncertainty

Despite the apparent advantages of adaptive plasticity, it is not common. We examined the effects of variation and uncertainty on selection for plasticity using an individual-based computer simulation model. In the model, the environment consisted of a linear gradient of 50 demes with dispersal occurring either before or after selection. Individuals consisted of multiple loci whose phenotypic expression either are affected (plastic) or are not affected (nonplastic) by the environment. Typically, evolution occurred first as genetic differentiation, which was then replaced by the evolution of adaptive plasticity, opposite to the evolutionary trend that is often assumed. Increasing dispersal rates selected for plasticity, if selection occurred before dispersal. If selection occurred after dispersal, the highest plasticity was at intermediate dispersal rates. Temporal variation in the environment occurring after development, but before selection, favored the evolution of plasticity. With dispersal before selection, such temporal variation resulted in hyperplasticity, with a reaction norm much steeper than the optimum. This effect was enhanced with negative temporal autocorrelation and can be interpreted as representing a form of bet hedging. As the number of nonplastic loci increased, plasticity was disfavored due to an increase in the uncertainty of the genomic environment. This effect was reversed with temporal variation. Thus, variation and uncertainty affect whether or not plasticity is favored with different sources of variation-arising from the amount and timing of dispersal, from temporal variation, and even from the genetic architecture underlying the phenotype-having contrasting, interacting, and at times unexpected effects.

Genetic architecture of life history traits and environment-specific trade-offs

Life history theory predicts the evolution of trait combinations that enhance fitness, and the occurrence of trade-offs depends in part on the magnitude of variation in growth rate or acquisition. Using recombinant inbred lines, we examined the genetic architecture of age and size at reproduction across abiotic conditions encountered by cultivars and naturalized populations of Brassica rapa. We found that genotypes are plastic to seasonal setting, such that reproduction was accelerated under conditions encountered by summer annual populations and genetic variances for age at reproduction varied across simulated seasonal settings. Using an acquisition-allocation model, we predicted the likelihood of trade-offs. Consistent with predicted relationships, we observed a trade-off where early maturity is associated with small size at maturity under simulated summer and fall annual conditions but not under winter annual conditions. The trade-off in the summer annual setting was observed despite significant genotypic variation in growth rate, which is often expected to decouple age and size at reproduction because rapidly growing genotypes could mature early and attain a larger size relative to slowly growing genotypes that mature later. The absence of a trade-off in the winter setting is presumably attributable to the absence of genotypic differences in age at reproduction. We observed QTL for age at reproduction that jointly regulated size at reproduction in both the summer and fall annual settings, but these QTL were environment-specific (i.e. different QTL contributed to the trade-off in the fall vs. summer annual settings). Thus, at least some of the genetic mechanisms underlying observed trade-offs differed across environments.

The genetic architecture of ecophysiological and circadian traits in Brassica rapa

Developmental mechanisms that enable perception of and response to the environment may enhance fitness. Ecophysiological traits typically vary depending on local conditions and contribute to resource acquisition and allocation, yet correlations may limit adaptive trait expression. Notably, photosynthesis and stomatal conductance vary diurnally, and the circadian clock, which is an internal estimate of time that anticipates diurnal light/dark cycles, may synchronize physiological behaviors with environmental conditions. Using recombinant inbred lines of Brassica rapa, we examined the quantitative-genetic architecture of ecophysiological and phenological traits and tested their association with the circadian clock. We also investigated how trait expression differed across treatments that simulated seasonal settings encountered by crops and naturalized populations. Many ecophysiological traits were correlated, and some correlations were consistent with expected biophysical constraints; for example, stomata jointly regulate photosynthesis and transpiration by affecting carbon dioxide and water vapor diffusion across leaf surfaces, and these traits were correlated. Interestingly, some genotypes had unusual combinations of ecophysiological traits, such as high photosynthesis in combination with low stomatal conductance or leaf nitrogen, and selection on these genotypes could provide a mechanism for crop improvement. At the genotypic and QTL level, circadian period was correlated with leaf nitrogen, instantaneous measures of photosynthesis, and stomatal conductance as well as with a long-term proxy (carbon isotope discrimination) for gas exchange, suggesting that gas exchange is partly regulated by the clock and thus synchronized with daily light cycles. The association between circadian rhythms and ecophysiological traits is relevant to crop improvement and adaptive evolution.

Quantitative trait loci mapping of phenotypic plasticity and genotype-environment interactions in plant and insect performance

Community genetic studies generally ignore the plasticity of the functional traits through which the effect is passed from individuals to the associated community. However, the ability of organisms to be phenotypically plastic allows them to rapidly adapt to changing environments and plasticity is commonly observed across all taxa. Owing to the fitness benefits of phenotypic plasticity, evolutionary biologists are interested in its genetic basis, which could explain how phenotypic plasticity is involved in the evolution of species interactions. Two current ideas exist: (i) phenotypic plasticity is caused by environmentally sensitive loci associated with a phenotype; (ii) phenotypic plasticity is caused by regulatory genes that simply influence the plasticity of a phenotype. Here, we designed a quantitative trait loci (QTL) mapping experiment to locate QTL on the barley genome associated with barley performance when the environment varies in the presence of aphids, and the composition of the rhizosphere. We simultaneously mapped aphid performance across variable rhizosphere environments. We mapped main effects, QTL × environment interaction (QTL×E), and phenotypic plasticity (measured as the difference in mean trait values) for barley and aphid performance onto the barley genome using an interval mapping procedure. We found that QTL associated with phenotypic plasticity were co-located with main effect QTL and QTL×E. We also located phenotypic plasticity QTL that were located separately from main effect QTL. These results support both of the current ideas of how phenotypic plasticity is genetically based and provide an initial insight into the functional genetic basis of how phenotypically plastic traits may still be important sources of community genetic effects.

Gene-environment interactions at nucleotide resolution

Interactions among genes and the environment are a common source of phenotypic variation. To characterize the interplay between genetics and the environment at single nucleotide resolution, we quantified the genetic and environmental interactions of four quantitative trait nucleotides (QTN) that govern yeast sporulation efficiency. We first constructed a panel of strains that together carry all 32 possible combinations of the 4 QTN genotypes in 2 distinct genetic backgrounds. We then measured the sporulation efficiencies of these 32 strains across 8 controlled environments. This dataset shows that variation in sporulation efficiency is shaped largely by genetic and environmental interactions. We find clear examples of QTN:environment, QTN: background, and environment:background interactions. However, we find no QTN:QTN interactions that occur consistently across the entire dataset. Instead, interactions between QTN only occur under specific combinations of environment and genetic background. Thus, what might appear to be a QTN:QTN interaction in one background and environment becomes a more complex QTN:QTN:environment:background interaction when we consider the entire dataset as a whole. As a result, the phenotypic impact of a set of QTN alleles cannot be predicted from genotype alone. Our results instead demonstrate that the effects of QTN and their interactions are inextricably linked both to genetic background and to environmental variation.

Phenotypic variation among individuals is caused by naturally occurring genetic differences, or alleles. The relationship between an allele and the phenotype is extremely complex; for example, the effect of an allele often depends upon both the environment and the individual's genetic background. To better understand these complex relationships, we examined the effects of four quantitative trait nucleotides (QTN) in three genes that cause variation in sporulation efficiency between vineyard and oak tree strains of yeast. We measured the effects of the QTN while varying both the genetic makeup of the strains and their growth environments. We found that the effects of each of the four QTN alleles depended upon the genotypes at the other QTN, the growth environment, and whether the strain carried the oak or vineyard parent genome. There were no simple rules that describe the effects of the alleles across all environments; instead, detailed models were needed to account for environmental and genetic variation in order to predict the effects of alleles in specific individuals.

The quantitative-genetic and QTL architecture of trait integration and modularity in Brassica rapa across simulated seasonal settings

Within organisms, groups of traits with different functions are frequently modular, such that variation among modules is independent and variation within modules is tightly integrated, or correlated. Here, we investigated patterns of trait integration and modularity in Brassica rapa in response to three simulated seasonal temperature/photoperiod conditions. The goals of this research were to use trait correlations to understand patterns of trait integration and modularity within and among floral, vegetative and phenological traits of B. rapa in each of three treatments, to examine the QTL architecture underlying patterns of trait integration and modularity, and to quantify how variation in temperature and photoperiod affects the correlation structure and QTL architecture of traits. All floral organs of B. rapa were strongly correlated, and contrary to expectations, floral and vegetative traits were also correlated. Extensive QTL co-localization suggests that covariation of these traits is likely due to pleiotropy, although physically linked loci that independently affect individual traits cannot be ruled out. Across treatments, the structure of genotypic and QTL correlations was generally conserved. Any observed variation in genetic architecture arose from genotype × environment interactions (GEIs) and attendant QTL × E in response to temperature but not photoperiod.

Mutations and quantitative genetic variation: lessons from Drosophila

A central issue in evolutionary quantitative genetics is to understand how genetic variation for quantitative traits is maintained in natural populations. Estimates of genetic variation and of genetic correlations and pleiotropy among multiple traits, inbreeding depression, mutation rates for fitness and quantitative traits and of the strength and nature of selection are all required to evaluate theoretical models of the maintenance of genetic variation. Studies in Drosophila melanogaster have shown that a substantial fraction of segregating variation for fitness-related traits in Drosophila is due to rare deleterious alleles maintained by mutation-selection balance, with a smaller but significant fraction attributable to intermediate frequency alleles maintained by alleles with antagonistic pleiotropic effects, and late-age-specific effects. However, the nature of segregating variation for traits under stabilizing selection is less clear and requires more detailed knowledge of the loci, mutation rates, allelic effects and frequencies of molecular polymorphisms affecting variation in suites of pleiotropically connected traits. Recent studies in D. melanogaster have revealed unexpectedly complex genetic architectures of many quantitative traits, with large numbers of pleiotropic genes and alleles with sex-, environment- and genetic background-specific effects. Future genome wide association analyses of many quantitative traits on a common panel of fully sequenced Drosophila strains will provide much needed empirical data on the molecular genetic basis of quantitative traits.

Cryptic intraspecific variation in sex determination in Caenorhabditis elegans revealed by mutations

Sex determination mechanisms (SDMs) show striking diversity and appear to evolve rapidly. Although interspecific comparisons and studies of ongoing major transitions in sex determination (such as the establishment of new sex chromosomes) have shed light on how SDMs evolve, comparatively little attention has been paid to intraspecific variation with less drastic effects. In this study, I used mutant strains carrying a temperature-sensitive sex determination mutation, along with a second null mutation, in different wild genetic backgrounds to uncover hidden variation in the SDM of the model nematode Caenorhabditis elegans. I then used quantitative trait locus (QTL) mapping to begin to investigate its genetic basis. I identified several QTLs, and although this variation apparently involved genotype-by-temperature interactions, QTL effects were generally consistent across temperatures. These QTLs collectively and individually explained a relatively large fraction of the variance in tail morphology (a sexually dimorphic trait), and two QTLs contained no genes known to be involved in somatic sex determination. These results show the existence of within-species variation in sex determination in this species, and underscore the potential for microevolutionary change in this important developmental pathway.

Genetical toxicogenomics in Drosophila identifies master-modulatory loci that are regulated by developmental exposure to lead

The genetics of gene expression in recombinant inbred lines (RILs) can be mapped as expression quantitative trait loci (eQTLs). So-called "genetical genomics" studies have identified locally acting eQTLs (cis-eQTLs) for genes that show differences in steady-state RNA levels. These studies have also identified distantly acting master-modulatory trans-eQTLs that regulate tens or hundreds of transcripts (hotspots or transbands). We expand on these studies by performing genetical genomics experiments in two environments in order to identify trans-eQTL that might be regulated by developmental exposure to the neurotoxin lead. Flies from each of 75 RIL were raised from eggs to adults on either control food (made with 250 microM sodium acetate), or lead-treated food (made with 250 microM lead acetate, PbAc). RNA expression analyses of whole adult male flies (5-10 days old) were performed with Affymetrix DrosII whole genome arrays (18,952 probesets). Among the 1389 genes with cis-eQTL, there were 405 genes unique to control flies and 544 genes unique to lead-treated ones (440 genes had the same cis-eQTLs in both samples). There are 2396 genes with trans-eQTL which mapped to 12 major transbands with greater than 95 genes. Permutation analyses of the strain labels but not the expression data suggests that the total number of eQTL and the number of transbands are more important criteria for validation than the size of the transband. Two transbands, one located on the 2nd chromosome and one on the 3rd chromosome, co-regulate 33 lead-induced genes, many of which are involved in neurodevelopmental processes. For these 33 genes, rather than allelic variation at one locus exerting differential effects in two environments, we found that variation at two different loci are required for optimal effects on lead-induced expression.

Pleiotropic effects of environment-specific adaptation in Arabidopsis thaliana

Local adaptation may be important for the preservation of genetic diversity and the promotion of speciation. However, local adaptation may also constrain establishment in different environments. The consequences of local adaptation depend strongly on the pleiotropic effects of the genes involved in adaptation. Here, we investigated the pleiotropic effects of the genetic response to selection in outbred lines of Arabidopsis artificially selected to flower earlier under both winter- and spring-annual simulated conditions. The consequences of adaptation were evaluated by reciprocally transplanting selected and control lines between the two conditions. Selected lines always flower earlier than their controls, independent of growing conditions. However, selected lines, growing in the same condition in which they were selected, flower earlier than plants selected in the alternative environment. Plants selected to flower earlier in spring produce more fruits than controls when growing in the spring, and less fruits when growing in the winter; indicating that local adaptation has negative pleiotropic effects in another environment. Our results indicate that local adaptation can arise even when selection targets the same trait in the same direction. Furthermore, it suggests that adaptation under the two different environments can generate fitness trade-offs that can maintain genetic variation for flowering time.

Quantitative trait loci for aggressive behavior in Drosophila melanogaster

Aggressive behavior is observed across animal taxa and is likely to be evolutionarily conserved. Although potentially advantageous, aggression can have social and health consequences in humans, and is a component of a number of psychiatric disorders. As a complex genetic trait, it is modulated by numerous quantitative trait loci (QTL) with allelic effects that can vary in direction and magnitude and that are sensitive to environmental perturbations. Assays to quantify aggressive behavior in Drosophila melanogaster have been developed, making this an ideal model system in which to dissect the genomic architecture underlying manifestation of and variation in aggressive behavior. Here, we map QTL affecting variation in aggression between two wild-type Drosophila strains. We identified a minimum of five QTL in a genomewide scan: two on chromosome 2 and three on chromosome 3. At least three and possibly all five of these QTL interact epistatically. We used deficiency complementation mapping to subdivide two linked, epistatically interacting QTL of large effect on chromosome 3 into at least six QTL, and complementation tests to mutations identified four candidate quantitative trait genes. Extensive epistasis poses a serious challenge for understanding the genetic basis of complex traits.

Genetic architecture of quantitative traits in mice, flies, and humans

We compare and contrast the genetic architecture of quantitative phenotypes in two genetically well-characterized model organisms, the laboratory mouse, Mus musculus, and the fruit fly, Drosophila melanogaster, with that found in our own species from recent successes in genome-wide association studies. We show that the current model of large numbers of loci, each of small effect, is true for all species examined, and that discrepancies can be largely explained by differences in the experimental designs used. We argue that the distribution of effect size of common variants is the same for all phenotypes regardless of species, and we discuss the importance of epistasis, pleiotropy, and gene by environment interactions. Despite substantial advances in mapping quantitative trait loci, the identification of the quantitative trait genes and ultimately the sequence variants has proved more difficult, so that our information on the molecular basis of quantitative variation remains limited. Nevertheless, available data indicate that many variants lie outside genes, presumably in regulatory regions of the genome, where they act by altering gene expression. As yet there are very few instances where homologous quantitative trait loci, or quantitative trait genes, have been identified in multiple species, but the availability of high-resolution mapping data will soon make it possible to test the degree of overlap between species.

Variations at a quantitative trait locus (QTL) affect development of behavior in lead-exposed Drosophila melanogaster

We developed Drosophila melanogaster as a model to study correlated behavioral, neuronal and genetic effects of the neurotoxin lead, known to affect cognitive and behavioral development in children. We showed that, as in vertebrates, lead affects both synaptic development and complex behaviors (courtship, fecundity, locomotor activity) in Drosophila. By assessing differential behavioral responses to developmental lead exposure among recombinant inbred Drosophila lines (RI), derived from parental lines Oregon R and Russian 2b, we have now identified a genotype by environment interaction (GEI) for a behavioral trait affected by lead. Drosophila Activity Monitors (TriKinetics, Waltham, MA), which measure activity by counting the number of times a single fly in a small glass tube walks through an infrared beam aimed at the middle of the tube, were used to measure activity of flies, reared from eggs to 4 days of adult age on either control or lead-contaminated medium, from each of 75 RI lines. We observed a significant statistical association between the effect of lead on Average Daytime Activity (ADA) across lines and one marker locus, 30AB, on chromosome 2; we define this as a Quantitative Trait Locus (QTL) associated with behavioral effects of developmental lead exposure. When 30AB was from Russian 2b, lead significantly increased locomotor activity, whereas, when 30AB was from Oregon R, lead decreased it. 30AB contains about 125 genes among which are likely "candidate genes" for the observed lead-dependent behavioral changes. Drosophila are thus a useful, underutilized model for studying behavioral, synaptic and genetic changes following chronic exposure to lead or other neurotoxins during development.

Quantitative trait loci affecting phenotypic plasticity and the allometric relationship of ovariole number and thorax length in Drosophila melanogaster

Environmental factors during juvenile growth such as temperature and nutrition have major effects on adult morphology and life-history traits. In Drosophila melanogaster, ovary size, measured as ovariole number, and body size, measured as thorax length, are developmentally plastic traits with respect to larval nutrition. Herein we investigated the genetic basis for plasticity of ovariole number and body size, as well the genetic basis for their allometric relationship using recombinant inbred lines (RILs) derived from a natural population in Winters, California. We reared 196 RILs in four yeast concentrations and measured ovariole number and body size. The genetic correlation between ovariole number and thorax length was positive, but the strength of this correlation decreased with increasing yeast concentration. Genetic variation and genotype-by-environment (G x E) interactions were observed for both traits. We identified quantitative trait loci (QTL), epistatic, QTL-by-environment, and epistatic-by-environment interactions for both traits and their scaling relationships. The results are discussed in the context of multivariate trait evolution.

Phenotypic plasticity and genotype by environment interaction for olfactory behavior in Drosophila melanogaster

Genotype by environment interactions (GEI) play a major part in shaping the genetic architecture of quantitative traits and are confounding factors in genetic studies, for example, in attempts to associate genetic variation with disease susceptibility. It is generally not known what proportion of phenotypic variation is due to GEI and how many and which genes contribute to GEI. Behaviors are complex traits that mediate interactions with the environment and, thus, are ideally suited for studies of GEI. Olfactory behavior in Drosophila melanogaster presents an opportunity to systematically dissect GEI, since large numbers of genetically identical individuals can be reared under defined environmental conditions and the olfactory system of Drosophila and its behavioral response to odorants have been well characterized. We assessed variation in olfactory behavior in a population of 41 wild-derived inbred lines and asked to what extent different larval-rearing environments would influence adult olfactory behavior and whether GEI is a minor or major contributing source of phenotypic variation. We found that approximately 50% of phenotypic variation in adult olfactory behavior is attributable to GEI. In contrast, transcriptional analysis revealed that only 20 genes show GEI at the level of gene expression [false discovery rate (FDR) < 0.05], some of which are associated with physiological responses to environmental chemicals. Quantitative complementation tests with piggyBac-tagged mutants for 2 of these genes (CG9664 and Transferrin 1) demonstrate that genes that show transcriptional GEI are candidate genes for olfactory behavior and that GEI at the level of gene expression is correlated with GEI at the level of phenotype.

Of flies and man: Drosophila as a model for human complex traits

Understanding the genetic and environmental factors affecting human complex genetic traits and diseases is a major challenge because of many interacting genes with individually small effects, whose expression is sensitive to the environment. Dissection of complex traits using the powerful genetic approaches available with Drosophila melanogaster has provided important lessons that should be considered when studying human complex traits. In Drosophila, large numbers of pleiotropic genes affect complex traits; quantitative trait locus alleles often have sex-, environment-, and genetic background-specific effects, and variants associated with different phenotypic are in noncoding as well as coding regions of candidate genes. Such insights, in conjunction with the strong evolutionary conservation of key genes and pathways between flies and humans, make Drosophila an excellent model system for elucidating the genetic mechanisms that affect clinically relevant human complex traits, such as alcohol dependence, sleep, and neurodegenerative diseases.

Single-locus polymorphism in a heterogeneous two-deme model

Environmental heterogeneity has long been considered a likely explanation for the high levels of genetic variation found in most natural populations: selection in a spatially heterogeneous environment can maintain more variation. While this theoretical result has been extensively studied in models with limited parameters (e.g., two alleles, fixed gene flow, and particular selection schemes), the effect of spatial heterogeneity is poorly understood for models with a wider range of parameters (e.g., multiple alleles, different levels of gene flow, and more general selection schemes). We have compared the volume of fitness space that maintains variation in a single-deme model to the volume in a two-deme model for multiple alleles, random selection schemes, and various levels of migration. Furthermore, equilibrium allele-frequency vectors were examined to see if particular patterns of variation are more prevalent than first expected. The two-deme model maintains variation for substantially larger volumes of fitness space with lower heterozygote fitness than the single-deme model. This result implies that selection schemes in the two-deme model can have a wider range of fitness patterns while still maintaining variation. The equilibrium allele-frequency patterns emerging from the two-deme model are more variable and strongly influenced by gene flow.

Joint estimates of quantitative trait locus effect and frequency using synthetic recombinant populations of Drosophila melanogaster

We develop and implement a strategy to map QTL in two synthetic populations of Drosophila melanogaster each initiated with eight inbred founder strains. These recombinant populations allow simultaneous estimates of QTL location, effect, and frequency. Five X-linked QTL influencing bristle number were resolved to intervals of approximately 1.3 cM. We confirm previous observations of bristle number QTL distal to 4A at the tip of the chromosome and identify two novel QTL in 7F-8C, an interval that does not include any classic bristle number candidate genes. If QTL at the tip of the X are biallelic they appear to be intermediate in frequency, although there is evidence that these QTL may reside in multiallelic haplotypes. Conversely, the two QTL mapping to the middle of the X chromosome are likely rare: in each case the minor allele is observed in only 1 of the 16 founders. Assuming additivity and biallelism we estimate that identified QTL contribute 1.0 and 8.7%, respectively, to total phenotypic variation in male abdominal and sternopleural bristle number in nature. Models that seek to explain the maintenance of genetic variation make different predictions about the population frequency of QTL alleles. Thus, mapping QTL in eight-way recombinant populations can distinguish between these models.

Ain't misbehavin' ? Genotype-environment interactions and the genetics of behavior

Although the multiple interacting genes affecting complex traits can readily be dissected, how much genotype-environment interactions contribute to variation in complex traits remains elusive. A recent study that quantified several behavioral phenotypes on the same mouse strains in different laboratories decades apart, shows that some behavioral differences between laboratories remain greatly replicable, whereas others are less robust over time. This report, together with studies from Drosophila, stresses the importance of understanding genotype-environment interactions.

How repeatable are associations between polymorphisms in achaete-scute and bristle number variation in Drosophila?

Currently, the relevance of common genetic variants--particularly those significantly associated with phenotypic variation in laboratory studies--to standing phenotypic variation in the wild is poorly understood. To address this, we quantified the relationship between achaete-scute complex (ASC) polymorphisms and Drosophila bristle number phenotypes in several new population samples. MC22 is a biallelic, nonrepetitive-length polymorphism 97 bp downstream of the scute transcript. It has been previously shown to be associated with sternopleural bristle number variation in both sexes in a set of isogenic lines. We replicated this association in a large cohort of wild-caught Drosophila melanogaster. We also detected a significant association at MC22 in an outbred population maintained under laboratory conditions for approximately 25 years, but the phenotypic effects in this sample were opposite from the direction estimated in the initial study. Finally, no significant associations were detected in a second large wild-caught cohort or in a set of 134 nearly isogenic lines. Our ability to repeat the initial association in wild samples suggests that it was not spurious. Nevertheless, inconsistent results from the other three panels suggest that the relationship between polymorphic genetic markers and loci contributing to continuous variation is not a simple one.

Mapping phenotypic plasticity and genotype-environment interactions affecting life-history traits in Caenorhabditis elegans

Phenotypic plasticity and genotype-environment interactions (GEI) play an important role in the evolution of life histories. Knowledge of the molecular genetic basis of plasticity and GEI provides insight into the underlying mechanisms of life-history changes in different environments. We used a genomewide single-nucleotide polymorphism map in a recombinant N2 x CB4856 inbred panel of the nematode Caenorhabditis elegans to study the genetic control of phenotypic plasticity to temperature in four fitness-related traits, that is, age at maturity, fertility, egg size and growth rate. We mapped quantitative trait loci (QTL) for the respective traits at 12 and 24 degrees C, as well as their plasticities. We found genetic variation and GEI for age at maturity, fertility, egg size and growth rate. GEI in fertility and egg size was attributed to changes in rank order of reaction norms. In case of age at maturity and growth rate, GEI was caused mainly by differences in the among-line variance. In total, 11 QTLs were detected, five QTL at 12 degrees C and six QTL at 24 degrees C, which were associated with life-history traits. Five QTL associated with age at maturity, fertility and growth rate showed QTL x environment interaction. These colocalized with plasticity QTL for the respective traits suggesting allelic sensitivity to temperature. Further fine mapping, complementation analyses and gene silencing are planned to identify candidate genes underlying phenotypic plasticity for age at maturity, fertility and growth.

Genotype by environment interactions in viability and developmental time in populations of cactophilic Drosophila

The genetic and ecological basis of viability and developmental time differences between Drosophila buzzatii and D. koepferae were analysed using the isofemale line technique. Several isofemale lines were sampled from pairs of allopatric/sympatric populations of each species. Flies were reared in media prepared with decaying tissues of two of the main natural cactus hosts of each species. This experimental design enabled us to evaluate the relative contribution of phenotypic plasticity, genetic variation and genotype by environment interaction (G x E) to total phenotypic variation for two fitness traits, viability and developmental time. Our results revealed significant G x E in both traits, suggesting that the maintenance of genetic variation can be explained, at least in part, by diversifying selection in different patches of a heterogeneous environment in both species. However, the relative importance of the factors involved in the G x E varied between traits and populations within species. For viability, the G x E can be mainly attributed to changes in the rank order of lines across cacti. However, the pattern was different for developmental time. In D. buzzatii the G x E can be mainly accounted for by changes in among line variance across cacti, whereas changes in the rank order of lines across cacti was the main component in D. koepferae. These dissimilar patterns of variation between traits and species suggest that the evolutionary forces shaping genetic variation for developmental time and viability vary between populations within species and between species.

Epistasis for fitness-related quantitative traits in Arabidopsis thaliana grown in the field and in the greenhouse

The extent to which epistasis contributes to adaptation, population differentiation, and speciation is a long-standing and important problem in evolutionary genetics. Using recombinant inbred (RI) lines of Arabidopsis thaliana grown under natural field conditions, we have examined the genetic architecture of fitness-correlated traits with respect to epistasis; we identified both single-locus additive and two-locus epistatic QTL for natural variation in fruit number, germination, and seed length and width. For fruit number, we found seven significant epistatic interactions, but only two additive QTL. For seed germination, length, and width, there were from two to four additive QTL and from five to eight epistatic interactions. The epistatic interactions were both positive and negative. In each case, the magnitude of the epistatic effects was roughly double that of the effects of the additive QTL, varying from -41% to +29% for fruit number and from -5% to +4% for seed germination, length, and width. A number of the QTL that we describe participate in more than one epistatic interaction, and some loci identified as additive also may participate in an epistatic interaction; the genetic architecture for fitness traits may be a network of additive and epistatic effects. We compared the map positions of the additive and epistatic QTL for germination, seed width, and seed length from plants grown in both the field and the greenhouse. While the total number of significant additive and epistatic QTL was similar under the two growth conditions, the map locations were largely different. We found a small number of significant epistatic QTL x environment effects when we tested directly for them. Our results support the idea that epistatic interactions are an important part of natural genetic variation and reinforce the need for caution in comparing results from greenhouse-grown and field-grown plants.

Drosophila bristles and the nature of quantitative genetic variation

Numbers of Drosophila sensory bristles present an ideal model system to elucidate the genetic basis of variation for quantitative traits. Here, we review recent evidence that the genetic architecture of variation for bristle numbers is surprisingly complex. A substantial fraction of the Drosophila genome affects bristle number, indicating pervasive pleiotropy of genes that affect quantitative traits. Further, a large number of loci, often with sex- and environment-specific effects that are also conditional on background genotype, affect natural variation in bristle number. Despite this complexity, an understanding of the molecular basis of natural variation in bristle number is emerging from linkage disequilibrium mapping studies of individual candidate genes that affect the development of sensory bristles. We show that there is naturally segregating genetic variance for environmental plasticity of abdominal and sternopleural bristle number. For abdominal bristle number this variance can be attributed in part to an abnormal abdomen-like phenotype that resembles the phenotype of mutants defective in catecholamine biosynthesis. Dopa decarboxylase (Ddc) encodes the enzyme that catalyses the final step in the synthesis of dopamine, a major Drosophila catecholamine and neurotransmitter. We found that molecular polymorphisms at Ddc are indeed associated with variation in environmental plasticity of abdominal bristle number.

A major bristle QTL from a selected population of Drosophila uncovers the zinc-finger transcription factor poils-au-dos, a repressor of achaete-scute

Traditional screens aiming at identifying genes regulating development have relied on mutagenesis. Here, we describe a new gene involved in bristle development, identified through the use of natural variation and selection. Drosophila melanogaster bears a pattern of 11 macrochaetes per heminotum. From a population initially sampled in Marrakech, a strain was selected for an increased number of thoracic macrochaetes. Using recombination and single nucleotide polymorphisms, the factor responsible was mapped to a single locus on the third chromosome, poils au dos, that encodes a zinc-finger-ZAD protein. The original, as well as new, presumed null, alleles of poils au dos, is associated with ectopic achaete-scute expression that results in the additional bristles. This suggests a possible role for Poils au dos as a repressor of achaete and scute. Ectopic expression appears to be independent of the activity of known cis-regulatory enhancer sequences at the achaete-scute complex that mediate activation at specific sites on the notum. The target sequences for Poils au dos activity were mapped to a 14 kb region around scute. In addition, we show that pad interacts synergistically with the repressor hairy and with Dpp signaling in posterior and anterior regions of the notum, respectively.

Quantitative trait loci and interaction effects responsible for variation in female postmating mortality in Drosophila simulans and D. sechellia introgression lines

Mating appears to inflict a cost to Drosophila females, resulting in a reduction of their lifespan shortly after mating. Males from different chromosome extracted lines differ significantly in their detrimental effects on postmating female survival, and seminal fluid proteins produced in the male accessory glands are at least partially responsible for the effect. This suggests that there is a genetic basis underlying the male inflicted effect on female's postmating mortality. However, the genes responsible for this effect remain elusive. Using males from introgression lines between D. simulans and D. sechellia genomes and a quantitative trait locus (QTL) mapping approach, we identified chromosomal regions that affect postmating mortality of females. We found a second chromosome QTL with an effect on average female lifespan after mating and a third chromosome QTL with an effect on postmating female mortality rate. Under the general observation of a faster divergence of sex-related genes among closely related species, it is predicted that genes for reproductive traits other than hybrid sterility will show evidence of epistatic effects when brought into a heterospecific background. We detected a significant epistatic genetic effect on postmating female mortality rate that supports this prediction.