Polygenic mutation in Drosophila melanogaster: Mapping spontaneous mutations affecting sensory bristle number

Empirical measures of mutational effects define neutral models of regulatory evolution in Saccharomyces cerevisiae

Understanding how phenotypes evolve requires disentangling the effects of mutation generating new variation from the effects of selection filtering it. Tests for selection frequently assume that mutation introduces phenotypic variation symmetrically around the population mean, yet few studies have tested this assumption by deeply sampling the distributions of mutational effects for particular traits. Here, we examine distributions of mutational effects for gene expression in the budding yeast Saccharomyces cerevisiae by measuring the effects of thousands of point mutations introduced randomly throughout the genome. We find that the distributions of mutational effects differ for the 10 genes surveyed and are inconsistent with normality. For example, all 10 distributions of mutational effects included more mutations with large effects than expected for normally distributed phenotypes. In addition, some genes also showed asymmetries in their distribution of mutational effects, with new mutations more likely to increase than decrease the gene's expression or vice versa. Neutral models of regulatory evolution that take these empirically determined distributions into account suggest that neutral processes may explain more expression variation within natural populations than currently appreciated.

Causes and consequences of genetic background effects illuminated by integrative genomic analysis

The phenotypic consequences of individual mutations are modulated by the wild-type genetic background in which they occur. Although such background dependence is widely observed, we do not know whether general patterns across species and traits exist or about the mechanisms underlying it. We also lack knowledge on how mutations interact with genetic background to influence gene expression and how this in turn mediates mutant phenotypes. Furthermore, how genetic background influences patterns of epistasis remains unclear. To investigate the genetic basis and genomic consequences of genetic background dependence of the scalloped(E3) allele on the Drosophila melanogaster wing, we generated multiple novel genome-level datasets from a mapping-by-introgression experiment and a tagged RNA gene expression dataset. In addition we used whole genome resequencing of the parental lines-two commonly used laboratory strains-to predict polymorphic transcription factor binding sites for SD. We integrated these data with previously published genomic datasets from expression microarrays and a modifier mutation screen. By searching for genes showing a congruent signal across multiple datasets, we were able to identify a robust set of candidate loci contributing to the background-dependent effects of mutations in sd. We also show that the majority of background-dependent modifiers previously reported are caused by higher-order epistasis, not quantitative noncomplementation. These findings provide a useful foundation for more detailed investigations of genetic background dependence in this system, and this approach is likely to prove useful in exploring the genetic basis of other traits as well.

Lines of least resistance and genetic architecture of house mouse (Mus musculus) mandible shape

Evolution along "lines of least resistance (LLR)" is a well-established phenomenon, which is ultimately based on variational properties of the developmental system. However, it has remained unclear which aspects of development are responsible for observed distributions of variation. This article suggests that house mouse mandible shape may present an opportunity for future research in this field. The genetic architecture of mouse mandible shape has been investigated in the laboratory, and its evolution has been surveyed in studies on natural variation. Here, I ask whether evolutionary diversification of mandible shape follows "LLR" by comparing principal directions of genetic and evolutionary variation, and I assess the potential contribution of specific genomic regions to evolutionary divergence along these directions. The role of two aspects of development, QTL number (presumably reflecting mutational target size) and canalization, is comparatively assessed. Overall, our results demonstrate a strong tendency of genetic and evolutionary systems to vary along similar directions, suggesting that mandible shape evolves along LLR at the level of populations, subspecies and species. At the level of genetic factors, effects bearing similarity to directions of evolution are significantly overrepresented, providing support for the mutational target size hypothesis. A role for canalization is not supported by a clear correlation between directions of evolution and size of genetic effects; however, the evidence for canalization remains ambiguous. These results provide some insights into how the developmental system may shape the variational properties of genetic systems and thus influence the direction of evolution.

Flowering time in maize: linkage and epistasis at a major effect locus

In a previous study, we identified a candidate fragment length polymorphism associated with flowering time variation after seven generations of selection for flowering time, starting from the maize inbred line F252. Here, we characterized the candidate region and identified underlying polymorphisms. Then, we combined QTL mapping, association mapping, and developmental characterization to dissect the genetic mechanisms responsible for the phenotypic variation. The candidate region contained the Eukaryotic Initiation Factor (eIF-4A) and revealed a high level of sequence and structural variation beyond the 3'-UTR of eIF-4A, including several insertions of truncated transposable elements. Using a biallelic single-nucleotide polymorphism (SNP) (C/T) in the candidate region, we confirmed its association with flowering time variation in a panel of 317 maize inbred lines. However, while the T allele was correlated with late flowering time within the F252 genetic background, it was correlated with early flowering time in the association panel with pervasive interactions between allelic variation and the genetic background, pointing to underlying epistasis. We also detected pleiotropic effects of the candidate polymorphism on various traits including flowering time, plant height, and leaf number. Finally, we were able to break down the correlation between flowering time and leaf number in the progeny of a heterozygote (C/T) within the F252 background consistent with causal loci in linkage disequilibrium. We therefore propose that both a cluster of tightly linked genes and epistasis contribute to the phenotypic variation for flowering time.

Dissecting genetic networks underlying complex phenotypes: the theoretical framework

Great progress has been made in genetic dissection of quantitative trait variation during the past two decades, but many studies still reveal only a small fraction of quantitative trait loci (QTLs), and epistasis remains elusive. We integrate contemporary knowledge of signal transduction pathways with principles of quantitative and population genetics to characterize genetic networks underlying complex traits, using a model founded upon one-way functional dependency of downstream genes on upstream regulators (the principle of hierarchy) and mutual functional dependency among related genes (functional genetic units, FGU). Both simulated and real data suggest that complementary epistasis contributes greatly to quantitative trait variation, and obscures the phenotypic effects of many 'downstream' loci in pathways. The mathematical relationships between the main effects and epistatic effects of genes acting at different levels of signaling pathways were established using the quantitative and population genetic parameters. Both loss of function and "co-adapted" gene complexes formed by multiple alleles with differentiated functions (effects) are predicted to be frequent types of allelic diversity at loci that contribute to the genetic variation of complex traits in populations. Downstream FGUs appear to be more vulnerable to loss of function than their upstream regulators, but this vulnerability is apparently compensated by different FGUs of similar functions. Other predictions from the model may account for puzzling results regarding responses to selection, genotype by environment interaction, and the genetic basis of heterosis.

The population genetics of mutations: good, bad and indifferent

Population genetics is fundamental to our understanding of evolution, and mutations are essential raw materials for evolution. In this introduction to more detailed papers that follow, we aim to provide an oversight of the field. We review current knowledge on mutation rates and their harmful and beneficial effects on fitness and then consider theories that predict the fate of individual mutations or the consequences of mutation accumulation for quantitative traits. Many advances in the past built on models that treat the evolution of mutations at each DNA site independently, neglecting linkage of sites on chromosomes and interactions of effects between sites (epistasis). We review work that addresses these limitations, to predict how mutations interfere with each other. An understanding of the population genetics of mutations of individual loci and of traits affected by many loci helps in addressing many fundamental and applied questions: for example, how do organisms adapt to changing environments, how did sex evolve, which DNA sequences are medically important, why do we age, which genetic processes can generate new species or drive endangered species to extinction, and how should policy on levels of potentially harmful mutagens introduced into the environment by humans be determined?

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.

Standing variation and new mutations both contribute to a fast response to selection for flowering time in maize inbreds

BACKGROUND

In order to investigate the rate and limits of the response to selection from highly inbred genetic material and evaluate the respective contribution of standing variation and new mutations, we conducted a divergent selection experiment from maize inbred lines in open-field conditions during 7 years. Two maize commercial seed lots considered as inbred lines, F252 and MBS847, constituted two biological replicates of the experiment. In each replicate, we derived an Early and a Late population by selecting and selfing the earliest and the latest individuals, respectively, to produce the next generation.

RESULTS

All populations, except the Early MBS847, responded to selection despite a short number of generations and a small effective population size. Part of the response can be attributed to standing genetic variation in the initial seed lot. Indeed, we identified one polymorphism initially segregating in the F252 seed lot at a candidate locus for flowering time, which explained 35% of the trait variation within the Late F252 population. However, the model that best explained our data takes into account both residual polymorphism in the initial seed lots and a constant input of heritable genetic variation by new (epi)mutations. Under this model, values of mutational heritability range from 0.013 to 0.025, and stand as an upper bound compare to what is reported in other species.

CONCLUSIONS

Our study reports a long-term divergent selection experiment for a complex trait, flowering time, conducted on maize in open-field conditions. Starting from a highly inbred material, we created within a few generations populations that strikingly differ from the initial seed lot for flowering time while preserving most of the phenotypic characteristics of the initial inbred. Such material is unique for studying the dynamics of the response to selection and its determinants. In addition to the fixation of a standing beneficial mutation associated with a large phenotypic effect, a constant input of genetic variance by new mutations has likely contributed to the response. We discuss our results in the context of the evolution and mutational dynamics of populations characterized by a small effective population size.

Epistatic interactions attenuate mutations affecting startle behaviour in Drosophila melanogaster

Epistasis is an important feature of the genetic architecture of quantitative traits. Previously, we showed that startle-induced locomotor behaviour of Drosophila melanogaster, a critical survival trait, is highly polygenic and exhibits epistasis. Here, we characterize epistatic interactions among homozygous P-element mutations affecting startle-induced locomotion in the Canton-S isogenic background and in 21 wild-derived inbred genetic backgrounds. We find pervasive epistasis for pairwise combinations of homozygous P-element insertional mutations as well as for mutations in wild-derived backgrounds. In all cases, the direction of the epistatic effects is to suppress the mutant phenotypes. The magnitude of the epistatic interactions in wild-derived backgrounds is highly correlated with the magnitude of the main effects of mutations, leading to phenotypic robustness of the startle response in the face of deleterious mutations. There is variation in the magnitude of epistasis among the wild-derived genetic backgrounds, indicating evolutionary potential for enhancing or suppressing effects of single mutations. These results provide a partial glimpse of the complex genetic network underlying the genetic architecture of startle behaviour and provide empirical support for the hypothesis that suppressing epistasis is the mechanism underlying genetic canalization of traits under strong stabilizing selection. Widespread suppressing epistasis will lead to underestimates of the main effects of quantitative trait loci (QTLs) in mapping experiments when not explicitly accounted for. In addition, suppressing epistasis could lead to underestimates of mutational variation for quantitative traits and overestimates of the strength of stabilizing selection, which has implications for maintenance of variation of complex traits by mutation-selection balance.

Mesosternal bristle number in a cosmopolitan drosophilid: an X-linked variable trait independent of sternopleural bristles

Mesosternal (MS) bristles in Drosophila are a pair of machrochaetae found at the sternal end of the sternopleural (STP) microchaetae, and are thought to be invariable. In a closely related drosophilid genus, Zaprionus, their number is four and, in contrast to Drosophila, they show interspecific and intraspecific variability. The genetic basis of MS bristle number variability was studied in Z. indianus, the only cosmopolitan species of the genus. The trait responded rapidly to selection and two lines were obtained, one lacking any bristles (0-0) and the other bearing the normal phenotype (2-2). Other symmetrical phenotypes, (1-1) and (3-3), could also be selected for, but with lesser success. By contrast, STP bristle number did not vary significantly between the two lines (0-0) and (2-2), revealing its genetic independence from MS bristle number. Reciprocal crosses between these two lines showed that MS bristle number is mainly influenced by a major gene on the X chromosome (i.e. F(1) males always resembled their mothers) with codominant expression (i.e. heterozygous F(1) females harboured an average phenotype of 2 bristles). However, trait penetrance was incomplete and backcrosses revealed that this variability was partly due to genetic modifiers, most likely autosomal. The canalization of MS bristle number was investigated under different temperatures, and the increased appearance of abnormal phenotypes mainly occurred at extreme temperatures. There was a bias, however, towards bristle loss, as shown by a liability (developmental map) analysis. Finally, when ancestral and introduced populations were compared, the latter were far less stable, suggesting that genetic bottlenecks may perturb the MS bristle number canalization system. MS bristle number, thus, appears to be an excellent model for investigating developmental canalization at both the quantitative and the molecular level.

Associations between sperm competition and natural variation in male reproductive genes on the third chromosome of Drosophila melanogaster

We applied association analysis to elucidate the genetic basis for variation in phenotypes affecting postcopulatory sexual selection in a natural population of Drosophila melanogaster. We scored 96 third chromosome substitution lines for nine phenotypes affecting sperm competitive ability and genotyped them at 72 polymorphisms in 13 male reproductive genes. Significant heterogeneity among lines (P < 0.01) was detected for all phenotypes except male-induced refractoriness (P = 0.053). We identified 24 associations (8 single-marker associations, 12 three-marker haplotype associations, and 4 cases of epistasis revealed by single-marker interactions). Fewer than 9 of these associations are likely to be false positives. Several associations were consistent with previous findings [Acp70A with the male's influence on the female's refractoriness to remating (refractory), Esterase-6 with a male's remating probability (remating) and a measure of female offspring production (fecundity)], but many are novel associations with uncharacterized seminal fluid proteins. Four genes showed evidence for pleiotropic effects [CG6168 with a measure of sperm competition (P2') and refractory, CG14560 with a defensive measure of sperm competition (P1') and a measure of female fecundity, Acp62F with P2' and a measure of female fecundity, and Esterase-6 with remating and a measure of female fecundity]. Our findings provide evidence that pleiotropy and epistasis are important factors in the genetic architecture of male reproductive success and show that haplotype analyses can identify associations missed in the single-marker approach.

High-resolution mapping of quantitative trait loci affecting increased life span in Drosophila melanogaster

Limited life span and senescence are near-universal characteristics of eukaryotic organisms, controlled by many interacting quantitative trait loci (QTL) with individually small effects, whose expression is sensitive to the environment. Analyses of mutations in model organisms have shown that genes affecting stress resistance and metabolism affect life span across diverse taxa. However, there is considerable segregating variation for life span in nature, and relatively little is known about the genetic basis of this variation. Replicated lines of Drosophila that have evolved increased longevity as a correlated response to selection for postponed senescence are valuable resources for identifying QTL affecting naturally occurring variation in life span. Here, we used deficiency complementation mapping to identify at least 11 QTL on chromosome 3 that affect variation in life span between five old (O) lines selected for postponed senescence and their five base (B) population control lines. Most QTL were sex specific, and all but one affected multiple O lines. The latter observation is consistent with alleles at intermediate frequency in the base population contributing to the response to selection for postponed senescence. The QTL were mapped with high resolution and contained from 12 to 170 positional candidate genes.