Mutations and quantitative genetic variation: lessons from Drosophila

Within-population plastic responses to combined thermal-nutritional stress differ from those in response to single stressors, and are genetically independent across traits in both males and females

Phenotypic plasticity helps animals to buffer the effects of increasing thermal and nutritional stress created by climate change. Plastic responses to single and combined stressors can vary among genetically diverged populations. However, less is known about how plasticity in response to combined stress varies among individuals within a population or whether such variation changes across life-history traits. This is important because individual variation within populations shapes population-level responses to environmental change. Here, we used isogenic lines of Drosophila melanogaster to assess the plasticity of egg-to-adult viability and sex-specific body size for combinations of 2 temperatures (25 °C or 28 °C) and 3 diets (standard diet, low caloric diet, or low protein:carbohydrate ratio diet). Our results reveal substantial within-population genetic variation in plasticity for egg-to-adult viability and wing size in response to combined thermal-nutritional stress. This genetic variation in plasticity was a result of cross-environment genetic correlations that were often < 1 for both traits, as well as changes in the expression of genetic variation across environments for egg-to-adult viability. Cross-sex genetic correlations for body size were weaker when the sexes were reared in different conditions, suggesting that the genetic basis of traits may change with the environment. Furthermore, our results suggest that plasticity in egg-to-adult viability is genetically independent from plasticity in body size. Importantly, plasticity in response to diet and temperature individually differed from plastic shifts in response to diet and temperature in combination. By quantifying plasticity and the expression of genetic variance in response to combined stress across traits, our study reveals the complexity of animal responses to environmental change, and the need for a more nuanced understanding of the potential for populations to adapt to ongoing climate change.

The condition-dependence of male genital size and shape

The male genitals of internal fertilisers evolve rapidly and divergently, and sexual selection is generally responsible for this. Many sexually selected traits are condition-dependent-with their expression dependent upon the resources available to be allocated to them-as revealed by genetic or environmental manipulations of condition. However, it is not clear whether male genitals are also condition-dependent. Here we manipulate condition in two ways (via inbreeding and diet) to test the condition-dependence of the genital arch of Drosophila simulans. We found that genital size but not genital shape suffered from inbreeding depression, whereas genital size and shape were affected by dietary manipulation of condition. The differential effects of these treatments likely reflect underlying genetic architecture that has been shaped by past selection: inbreeding depression is only expected when traits have a history of directional selection, while diet impacts traits regardless of historical selection. Nonetheless, our results suggest genitals can be condition-dependent like other sexually selected traits.

Pervasive G × E interactions shape adaptive trajectories and the exploration of the phenotypic space in artificial selection experiments

Quantitative genetics models have shown that long-term selection responses depend on initial variance and mutational influx. Understanding limits of selection requires quantifying the role of mutational variance. However, correlative responses to selection on nonfocal traits can perturb the selection response on the focal trait; and generations are often confounded with selection environments so that genotype by environment (G×E) interactions are ignored. The Saclay divergent selection experiments (DSEs) on maize flowering time were used to track the fate of individual mutations combining genotyping data and phenotyping data from yearly measurements (DSEYM) and common garden experiments (DSECG) with four objectives: (1) to quantify the relative contribution of standing and mutational variance to the selection response, (2) to estimate genotypic mutation effects, (3) to study the impact of G×E interactions in the selection response, and (4) to analyze how trait correlations modulate the exploration of the phenotypic space. We validated experimentally the expected enrichment of fixed beneficial mutations with an average effect of +0.278 and +0.299 days to flowering, depending on the genetic background. Fixation of unfavorable mutations reached up to 25% of incoming mutations, a genetic load possibly due to antagonistic pleiotropy, whereby mutations fixed in the selection environment (DSEYM) turned to be unfavorable in the evaluation environment (DSECG). Global patterns of trait correlations were conserved across genetic backgrounds but exhibited temporal patterns. Traits weakly or uncorrelated with flowering time triggered stochastic exploration of the phenotypic space, owing to microenvironment-specific fixation of standing variants and pleiotropic mutational input.

A possible genomic footprint of polygenic adaptation on population divergence in seed beetles?

Efforts to unravel the genomic basis of incipient speciation are hampered by a mismatch between our toolkit and our understanding of the ecology and genetics of adaptation. While the former is focused on detecting selective sweeps involving few independently acting or linked speciation genes, the latter states that divergence typically occurs in polygenic traits under stabilizing selection. Here, we ask whether a role of stabilizing selection on polygenic traits in population divergence may be unveiled by using a phenotypically informed integrative approach, based on genome‐wide variation segregating in divergent populations. We compare three divergent populations of seed beetles (Callosobruchus maculatus) where previous work has demonstrated a prominent role for stabilizing selection on, and population divergence in, key life history traits that reflect rate‐dependent metabolic processes. We derive and assess predictions regarding the expected pattern of covariation between genetic variation segregating within populations and genetic differentiation between populations. Population differentiation was considerable (mean F_ST = 0.23–0.26) and was primarily built by genes showing high selective constraints and an imbalance in inferred selection in different populations (positive Tajima's D_NS in one and negative in one), and this set of genes was enriched with genes with a metabolic function. Repeatability of relative population differentiation was low at the level of individual genes but higher at the level of broad functional classes, again spotlighting metabolic genes. Absolute differentiation (d_XY) showed a very different general pattern at this scale of divergence, more consistent with an important role for genetic drift. Although our exploration is consistent with stabilizing selection on polygenic metabolic phenotypes as an important engine of genome‐wide relative population divergence and incipient speciation in our study system, we note that it is exceedingly difficult to firmly exclude other scenarios.

Inter- and intraspecific variation in mycotoxin tolerance: A study of four Drosophila species

Many mycophagous Drosophila species have adapted to tolerate high concentrations of mycotoxins, an ability not reported in any other eukaryotes. Although an association between mycophagy and mycotoxin tolerance has been established in many Drosophila species, the genetic mechanisms of the tolerance are unknown. This study presents the inter‐ and intraspecific variation in the mycotoxin tolerance trait. We studied the mycotoxin tolerance in four Drosophila species from four separate clades within the immigrans‐tripunctata radiation from two distinct locations. The effect of mycotoxin treatment on 20 isofemale lines per species was studied using seven gross phenotypes: survival to pupation, survival to eclosion, development time to pupation and eclosion, thorax length, fecundity, and longevity. We observed interspecific variation among four species, with D. falleni being the most tolerant, followed by D. recens, D. neotestacea, and D. tripunctata, in that order. The results also revealed geographical variation and intraspecific genetic variation in mycotoxin tolerance. This report provides the foundation for further delineating the genetic mechanisms of the mycotoxin tolerance trait.

Natural variation in Drosophila shows weak pleiotropic effects

Background

Pleiotropy describes the phenomenon in which a gene affects multiple phenotypes. The extent of pleiotropy is still disputed, mainly because of issues of inadequate power of analyses. A further challenge is that empirical tests of pleiotropy are restricted to a small subset of all possible phenotypes. To overcome these limitations, we propose a new measurement of pleiotropy that integrates across many phenotypes and multiple generations to improve power.

Results

We infer pleiotropy from the fitness cost imposed by frequency changes of pleiotropic loci. Mixing Drosophila simulans populations, which adapted independently to the same new environment using different sets of genes, we show that the adaptive frequency changes have been accompanied by measurable fitness costs.

Conclusions

Unlike previous studies characterizing the molecular basis of pleiotropy, we show that many loci, each of weak effect, contribute to genome-wide pleiotropy. We propose that the costs of pleiotropy are reduced by the modular architecture of gene expression, which facilitates adaptive gene expression changes with low impact on other functions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02680-4.

Radiobiological features in offspring of natural populations of Drosophila melanogaster after Chernobyl accident

In their natural habitats, populations of organisms are faced with different levels of chronic low-intensity radiation, causing a wide range of radiobiological effects (from radiosensitivity to radioadaptive response and hormesis). In this study, specimens of Drosophila melanogaster were selected from territories of the Chernobyl nuclear power plant with different levels of radioactive contamination. The isogenic stocks derived from these specimens represent the genetic systems of current populations and make it possible to study radioresistance and its mechanisms in future generations under controlled laboratory conditions. Previous studies have shown that transgenerational radiation effects at the level of lethal mutations and survival rate are unstable and depend not only on the level of chronic low-intensity irradiation, but also on other factors. A single acute irradiation exposure of offspring whose parents inhabited a site with a higher level of chronic irradiation made it possible to reveal pronounced radioresistant features in the offspring. And the offspring whose parents were exposed to radiation levels close to the natural radiation background, on the contrary, acquired radiosensitive features. Their response to acute exposure includes a high-frequency of lethal mutations and a short lifespan. The differential response to different levels of chronic parental exposure is caused by differences in the activities of certain transposons that destabilize the genome. Our data contribute to the understanding of genetic and epigenetic mechanisms (via transposon activity) of the effect of parental radiation exposure on the health and adaptive potential of populations affected by the technogenically increased radiation background.

Interplay between extreme drift and selection intensities favors the fixation of beneficial mutations in selfing maize populations

Population and quantitative genetic models provide useful approximations to predict long-term selection responses sustaining phenotypic shifts, and underlying multilocus adaptive dynamics. Valid across a broad range of parameters, their use for understanding the adaptive dynamics of small selfing populations undergoing strong selection intensity (thereafter High Drift-High selection regime, HDHS) remains to be explored. Saclay Divergent Selection Experiments (DSEs) on maize flowering time provide an interesting example of populations evolving under HDHS, with significant selection responses over 20 generations in two directions. We combined experimental data from Saclay DSEs, forward individual-based simulations, and theoretical predictions to dissect the evolutionary mechanisms at play in the observed selection responses. We asked two main questions: How do mutations arise, spread, and reach fixation in populations evolving under HDHS? How does the interplay between drift and selection influence observed phenotypic shifts? We showed that the long-lasting response to selection in small populations is due to the rapid fixation of mutations occurring during the generations of selection. Among fixed mutations, we also found a clear signal of enrichment for beneficial mutations revealing a limited cost of selection. Both environmental stochasticity and variation in selection coefficients likely contributed to exacerbate mutational effects, thereby facilitating selection grasp and fixation of small-effect mutations. Together our results highlight that despite a small number of polymorphic loci expected under HDHS, adaptive variation is continuously fueled by a vast mutational target. We discuss our results in the context of breeding and long-term survival of small selfing populations.

The effects of the sex chromosomes on the inheritance of species-specific traits of the copulatory organ shape in Drosophila virilis and Drosophila lummei

The shape of the male genitalia in many taxa is the most rapidly evolving morphological structure, often driving reproductive isolation, and is therefore widely used in systematics as a key character to distinguish between sibling species. However, only a few studies have used the genital arch of the male copulatory organ as a model to study the genetic basis of species-specific differences in the Drosophila copulatory system. Moreover, almost nothing is known about the effects of the sex chromosomes on the shape of the male mating organ. In our study, we used a set of crosses between D. virilis and D. lummei and applied the methods of quantitative genetics to assess the variability of the shape of the male copulatory organ and the effects of the sex chromosomes and autosomes on its variance. Our results showed that the male genital shape depends on the species composition of the sex chromosomes and autosomes. Epistatic interactions of the sex chromosomes with autosomes and the species origin of the Y-chromosome in a male in interspecific crosses also influenced the expression of species-specific traits in the shape of the male copulatory system. Overall, the effects of sex chromosomes were comparable to the effects of autosomes despite the great differences in gene numbers between them. It may be reasonably considered that sexual selection for specific genes associated with the shape of the male mating organ prevents the demasculinization of the X chromosome.

Microevolutionary Dynamics of Chicken Genomes under Divergent Selection for Adiposity

Decades of artificial selection have significantly improved performance and efficiency of animal production systems. However, little is known about the microevolution of genomes due to intensive breeding. Using whole-genome sequencing, we document dynamic changes of chicken genomes under divergent selection on adiposity over 19 generations. Directional selection reduced within-line but increased between-line genomic differences. We observed that artificial selection tended to result in recruitment of preexisting variations of genes related to adipose tissue growth. In addition, novel mutations contributed to divergence of phenotypes under selection but contributed significantly less than preexisting genomic variants. Integration of 15 generations genome sequencing, genome-wide association study, and multi-omics data further identified that genes involved in signaling pathways important to adipogenesis, such as autophagy and lysosome (URI1, MBL2), neural system (CHAT), and endocrine (PCSK1) pathways, were under strong selection. Our study provides insights into the microevolutionary dynamics of domestic animal genomes under artificial selection.

Life-History Evolution and the Genetics of Fitness Components in Drosophila melanogaster

Life-history traits or "fitness components"-such as age and size at maturity, fecundity and fertility, age-specific rates of survival, and life span-are the major phenotypic determinants of Darwinian fitness. Analyzing the evolution and genetics of these phenotypic targets of selection is central to our understanding of adaptation. Due to its simple and rapid life cycle, cosmopolitan distribution, ease of maintenance in the laboratory, well-understood evolutionary genetics, and its versatile genetic toolbox, the "vinegar fly" Drosophila melanogaster is one of the most powerful, experimentally tractable model systems for studying "life-history evolution." Here, I review what has been learned about the evolution and genetics of life-history variation in D. melanogaster by drawing on numerous sources spanning population and quantitative genetics, genomics, experimental evolution, evolutionary ecology, and physiology. This body of work has contributed greatly to our knowledge of several fundamental problems in evolutionary biology, including the amount and maintenance of genetic variation, the evolution of body size, clines and climate adaptation, the evolution of senescence, phenotypic plasticity, the nature of life-history trade-offs, and so forth. While major progress has been made, important facets of these and other questions remain open, and the D. melanogaster system will undoubtedly continue to deliver key insights into central issues of life-history evolution and the genetics of adaptation.

Unpredictable Effects of the Genetic Background of Transgenic Lines in Physiological Quantitative Traits

Physiology, fitness and disease phenotypes are complex traits exhibiting continuous variation in natural populations. To understand complex trait gene functions transgenic lines of undefined genetic background are often combined to assess quantitative phenotypes ignoring the impact of genetic polymorphisms. Here, we used inbred wild-type strains of the Drosophila Genetics Reference Panel to assess the phenotypic variation of six physiological and fitness traits, namely, female fecundity, survival and intestinal mitosis upon oral infection, defecation rate and fecal pH upon oral infection, and terminal tracheal cell branching in hypoxia. We found continuous variation in the approximately 150 strains tested for each trait, with extreme values differing by more than four standard deviations for all traits. In addition, we assessed the effects of commonly used Drosophila UAS-RNAi transgenic strains and their backcrossed isogenized counterparts, in the same traits plus baseline intestinal mitosis and tracheal branching in normoxia, in heterozygous conditions, when only half of the genetic background was different among strains. We tested 20 non-isogenic strains (10 KK and 10 GD) from the Vienna Drosophila Resource Center and their isogenized counterparts without Gal4 induction. Survival upon infection and female fecundity exhibited differences in 50% and 40% of the tested isogenic vs. non-isogenic pairs, respectively, whereas all other traits were affected in only 10–25% of the cases. When 11 isogenic and their corresponding non-isogenic UAS-RNAi lines were expressed ubiquitously with Gal4, 4 isogenic vs. non-isogenic pairs exhibited differences in survival to infection. Furthermore, when a single UAS-RNAi line was crossed with the same Gal4 transgene inserted in different genetic backgrounds, the quantitative variations observed were unpredictable on the basis of pure line performance. Thus, irrespective of the trait of interest, the genetic background of commonly used transgenic strains needs to be considered carefully during experimentation.

Multi-trait Improvement by Predicting Genetic Correlations in Breeding Crosses

The many quantitative traits of interest to plant breeders are often genetically correlated, which can complicate progress from selection. Improving multiple traits may be enhanced by identifying parent combinations - an important breeding step - that will deliver more favorable genetic correlations (rG ). Modeling the segregation of genomewide markers with estimated effects may be one method of predicting rG in a cross, but this approach remains untested. Our objectives were to: (i) use simulations to assess the accuracy of genomewide predictions of rG and the long-term response to selection when selecting crosses on the basis of such predictions; and (ii) empirically measure the ability to predict genetic correlations using data from a barley (Hordeum vulgare L.) breeding program. Using simulations, we found that the accuracy to predict rG was generally moderate and influenced by trait heritability, population size, and genetic correlation architecture (i.e., pleiotropy or linkage disequilibrium). Among 26 barley breeding populations, the empirical prediction accuracy of rG was low (-0.012) to moderate (0.42), depending on trait complexity. Within a simulated plant breeding program employing indirect selection, choosing crosses based on predicted rG increased multi-trait genetic gain by 11-27% compared to selection on the predicted cross mean. Importantly, when the starting genetic correlation was negative, such cross selection mitigated or prevented an unfavorable response in the trait under indirect selection. Prioritizing crosses based on predicted genetic correlation can be a feasible and effective method of improving unfavorably correlated traits in breeding programs.

Genome-Wide Association Analysis of Anoxia Tolerance in Drosophila melanogaster

As the genetic bases to variation in anoxia tolerance are poorly understood, we used the Drosophila Genetics Reference Panel (DGRP) to conduct a genome-wide association study (GWAS) of anoxia tolerance in adult and larval Drosophila melanogaster Survival ranged from 0-100% in adults exposed to 6 h of anoxia and from 20-98% for larvae exposed to 1 h of anoxia. Anoxia tolerance had a broad-sense heritability of 0.552 in adults and 0.433 in larvae. Larval and adult phenotypes were weakly correlated but the anoxia tolerance of adult males and females were strongly correlated. The GWA identified 180 SNPs in adults and 32 SNPs in larvae associated with anoxia tolerance. Gene ontology enrichment analysis indicated that many of the 119 polymorphic genes associated with adult anoxia-tolerance were associated with ionic transport or immune function. In contrast, the 22 polymorphic genes associated with larval anoxia-tolerance were mostly associated with regulation of transcription and DNA replication. RNAi of mapped genes generally supported the hypothesis that disruption of these genes reduces anoxia tolerance. For two ion transport genes, we tested predicted directional and sex-specific effects of SNP alleles on adult anoxia tolerance and found strong support in one case but not the other. Correlating our phenotype to prior DGRP studies suggests that genes affecting anoxia tolerance also influence stress-resistance, immune function and ionic balance. Overall, our results provide evidence for multiple new potential genetic influences on anoxia tolerance and provide additional support for important roles of ion balance and immune processes in determining variation in anoxia tolerance.

Stress response, behavior, and development are shaped by transposable element-induced mutations in Drosophila

Most of the current knowledge on the genetic basis of adaptive evolution is based on the analysis of single nucleotide polymorphisms (SNPs). Despite increasing evidence for their causal role, the contribution of structural variants to adaptive evolution remains largely unexplored. In this work, we analyzed the population frequencies of 1,615 Transposable Element (TE) insertions annotated in the reference genome of Drosophila melanogaster, in 91 samples from 60 worldwide natural populations. We identified a set of 300 polymorphic TEs that are present at high population frequencies, and located in genomic regions with high recombination rate, where the efficiency of natural selection is high. The age and the length of these 300 TEs are consistent with relatively young and long insertions reaching high frequencies due to the action of positive selection. Besides, we identified a set of 21 fixed TEs also likely to be adaptive. Indeed, we, and others, found evidence of selection for 84 of these reference TE insertions. The analysis of the genes located nearby these 84 candidate adaptive insertions suggested that the functional response to selection is related with the GO categories of response to stimulus, behavior, and development. We further showed that a subset of the candidate adaptive TEs affects expression of nearby genes, and five of them have already been linked to an ecologically relevant phenotypic effect. Our results provide a more complete understanding of the genetic variation and the fitness-related traits relevant for adaptive evolution. Similar studies should help uncover the importance of TE-induced adaptive mutations in other species as well.

Influence of a growth hormone transgene on the genetic architecture of growth-related traits: A comparative analysis between transgenic and wild-type coho salmon

Genetic engineering has been increasingly applied to many commercially important plant and animal species, generating phenotypic changes that are not observed in natural populations and creating genetic interactions that have not experienced natural selection. The degree to and way in which such human-induced genetic variation interacts with the rest of the genome is currently largely unknown. Integrating such information into ecological and risk assessment frameworks is crucial to understand the potential effects of genetically modified organisms in natural environments. Here, we performed QTL mapping to investigate the genetic architecture of growth-related traits in nontransgenic (NT) and growth hormone transgenic (T) coho salmon with large changes in growth and related physiology, with the aim of identifying how an inserted transgene might influence the opportunity for selection. These fish shared the same parental genetic background, thus allowing us to determine whether the same or different loci influence these traits within the two groups. The use of over 1,700 loci, derived from restriction site-associated DNA sequencing, revealed that different genomic regions were linked with growth over time between the two groups. Additionally, the effect sizes of detected QTL appear to have been influenced by the transgene. Direct comparison of QTL between the T and NT fish during two size-matched periods identified little overlap in their location. Taken together, the results showed that the transgene altered the genetic basis of growth-related traits in this species. The study has important implications for effective conservation and management of wild populations experiencing introduction of transgenes. Evolutionary changes and their ecological consequences may occur at different rates and in different directions in NT versus T individuals in response to selection. Thus, assessments of phenotypic change, and hence ecological risk, should be determined periodically to evaluate whether initial estimates made with founder strains remain valid.

Dissection of the molecular bases of genotype x environment interactions: a study of phenotypic plasticity of Saccharomyces cerevisiae in grape juices

Background

The ability of a genotype to produce different phenotypes according to its surrounding environment is known as phenotypic plasticity. Within different individuals of the same species, phenotypic plasticity can vary greatly. This contrasting response is caused by gene-by-environment interactions (GxE). Understanding GxE interactions is particularly important in agronomy, since selected breeds and varieties may have divergent phenotypes according to their growing environment. Industrial microbes such as Saccharomyces cerevisiae are also faced with a large range of fermentation conditions that affect their technological properties. Finding the molecular determinism of such variations is a critical task for better understanding the genetic bases of phenotypic plasticity and can also be helpful in order to improve breeding methods.

Results

In this study we implemented a QTL mapping program using two independent cross (~ 100 progeny) in order to investigate the molecular basis of yeast phenotypic response in a wine fermentation context. Thanks to whole genome sequencing approaches, both crosses were genotyped, providing saturated genetic maps of thousands of markers. Linkage analyses allowed the detection of 78 QTLs including 21 with significant interaction with the environmental conditions. Molecular dissection of a major QTL demonstrated that the sulfite pump Ssu1p has a pleiotropic effect and impacts the phenotypic plasticity of several traits.

Conclusions

The detection of QTLs and their interactions with environment emphasizes the complexity of yeast industrial traits. The validation of the interaction of SSU1 allelic variants with the nature of the fermented juice increases knowledge about the impact of the sulfite pump during fermentation. All together these results pave the way for exploiting and deciphering the genetic determinism of phenotypic plasticity.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5145-4) contains supplementary material, which is available to authorized users.

Biocontrol characteristics of the fruit fly pupal parasitoid Trichopria drosophilae (Hymenoptera: Diapriidae) emerging from different hosts

Trichopria drosophilae (Hymenoptera: Diapriidae) is an important pupal endoparasitoid of Drosophila melanogaster Meigen (Diptera: Drosophilidae) and some other fruit fly species, such as D. suzukii, a very important invasive and economic pest. Studies of T. drosophilae suggest that this could be a good biological control agent for fruit fly pests. In this research, we compared the parasitic characteristics of T. drosophilae reared in D. melanogaster (TD_m) with those reared in D. hydei (TD_h). TD_h had a larger size than TD_m. The number of maximum mature eggs of a female TD_h was 133.6 ± 6.9, compared with the significantly lower value of 104.8 ± 11.4 for TD_m. Mated TD_h female wasp continuously produced female offspring up to 6 days after mating, compared with only 3 days for TD_m. In addition, the offspring female ratio of TD_h, i.e., 82.32%, was significantly higher than that of TD_m, i.e., 61.37%. Under starvation treatment, TD_h survived longer than TD_m. TD_h also survived longer than TD^m at high temperatures, such as 37 °C, although they both survived well at low temperatures, such as 18 °C and 4 °C. Old-age TD_h females maintained a high parasitism rate and offspring female ratio, while they were declined in old-age TD_m. Overall, TD_h had an advantage in terms of body size, fecundity, stress resistance ability and the parasitism rate compared with TD_m. Therefore, T. drosophilae from D. hydei could improve biocontrol efficacy with enormous economic benefits in the field, especially in the control of many frugivorous Drosophilidae species worldwide.

Quantitative trait locus mapping and analysis of heritable variation in affiliative social behavior and co-occurring traits

Humans exhibit broad heterogeneity in affiliative social behavior. Twin and family studies show that individual differences in core dimensions of social behavior are heritable, yet there are knowledge gaps in understanding the underlying genetic and neurobiological mechanisms. Animal genetic reference panels (GRPs) provide a tractable strategy for examining the behavioral and genetic architecture of complex traits. Here, using males from 50 mouse strains from the BXD GRP, 4 domains of affiliative social behavior-social approach, social recognition, direct social interaction (DSI) (partner sniffing) and vocal communication-were examined in 2 widely used behavioral tasks-the 3-chamber and DSI tasks. There was continuous and broad variation in social and nonsocial traits, with moderate to high heritability of social approach sniff preference (0.31), ultrasonic vocalization (USV) count (0.39), partner sniffing (0.51), locomotor activity (0.54-0.66) and anxiety-like behavior (0.36). Principal component analysis shows that variation in social and nonsocial traits are attributable to 5 independent factors. Genome-wide mapping identified significant quantitative trait loci for USV count on chromosome (Chr) 18 and locomotor activity on Chr X, with suggestive loci and candidate quantitative trait genes identified for all traits with one notable exception-partner sniffing in the DSI task. The results show heritable variation in sociability, which is independent of variation in activity and anxiety-like traits. In addition, a highly heritable and ethological domain of affiliative sociability-partner sniffing-appears highly polygenic. These findings establish a basis for identifying functional natural variants, leading to a new understanding typical and atypical sociability.

Antagonistic pleiotropy in species with separate sexes, and the maintenance of genetic variation in life-history traits and fitness

Antagonistic pleiotropy (AP)-where alleles of a gene increase some components of fitness at a cost to others-can generate balancing selection, and contribute to the maintenance of genetic variation in fitness traits, such as survival, fecundity, fertility, and mate competition. Previous theory suggests that AP is unlikely to maintain variation unless antagonistic selection is strong, or AP alleles exhibit pronounced differences in genetic dominance between the affected traits. We show that conditions for balancing selection under AP expand under the likely scenario that the strength of selection on each fitness component differs between the sexes. Our model also predicts that the vast majority of balanced polymorphisms have sexually antagonistic effects on total fitness, despite the absence of sexual antagonism for individual fitness components. We conclude that AP polymorphisms are less difficult to maintain than predicted by prior theory, even under our conservative assumption that selection on components of fitness is universally sexually concordant. We discuss implications for the maintenance of genetic variation, and for inferences of sexual antagonism that are based on sex-specific phenotypic selection estimates-many of which are based on single fitness components.

Dissecting genetic architecture of startle response in Drosophila melanogaster using multi-omics information

Startle behavior is important for survival, and abnormal startle responses are related to several neurological diseases. Drosophila melanogaster provides a powerful system to investigate the genetic underpinnings of variation in startle behavior. Since mechanically induced, startle responses and environmental conditions can be readily quantified and precisely controlled. The 156 wild-derived fully sequenced lines of the Drosophila Genetic Reference Panel (DGRP) were used to identify SNPs and transcripts associated with variation in startle behavior. The results validated highly significant effects of 33 quantitative trait SNPs (QTSs) and 81 quantitative trait transcripts (QTTs) directly associated with phenotypic variation of startle response. We also detected QTT variation controlled by 20 QTSs (tQTSs) and 73 transcripts (tQTTs). Association mapping based on genomic and transcriptomic data enabled us to construct a complex genetic network that underlies variation in startle behavior. Based on principles of evolutionary conservation, human orthologous genes could be superimposed on this network. This study provided both genetic and biological insights into the variation of startle response behavior of Drosophila melanogaster, and highlighted the importance of genetic network to understand the genetic architecture of complex traits.

Natural and laboratory mutations in kuzbanian are associated with zinc stress phenotypes in Drosophila melanogaster

Organisms must cope with altered environmental conditions such as high concentrations of heavy metals. Stress response to heavy metals is mediated by the metal-responsive transcription factor 1 (MTF-1), which is conserved from Drosophila to humans. MTF-1 binds to metal response elements (MREs) and changes the expression of target genes. kuzbanian (kuz), a metalloendopeptidase that activates the evolutionary conserved Notch signaling pathway, has been identified as an MTF-1 target gene. We have previously identified a putatively adaptive transposable element in the Drosophila melanogaster genome, named FBti0019170, inserted in a kuz intron. In this work, we investigated whether a laboratory mutant stock overexpressing kuz is associated with zinc stress phenotypes. We found that both embryos and adult flies overexpressing kuz are more tolerant to zinc compared with wild-type flies. On the other hand, we found that the effect of FBti0019170 on zinc stress tolerance depends on developmental stage and genetic background. Moreover, in the majority of the genetic backgrounds analyzed, FBti0019170 has a deleterious effect in unpolluted environments in pre-adult stages. These results highlight the complexity of natural mutations and suggest that besides laboratory mutations, natural mutations should be studied in order to accurately characterize gene function and evolution.

Natural quantitative genetic variance in plant growth differs in response to ecologically relevant temperature heterogeneity

Adaptation to large‐scale spatial heterogeneity in the environment accounts for a major proportion of genetic diversity within species. Theory predicts the erosion of adaptive genetic variation on a within‐population level, but considerable genetic diversity is often found locally. Genetic diversity could be expected to be maintained within populations in temporally or spatially variable conditions if genotypic rank orders vary across contrasting microenvironmental settings. Taking advantage of fine‐resolution environmental data, we tested the hypothesis that temperature heterogeneity among years could be one factor maintaining quantitative genetic diversity within a natural and genetically diverse plant population. We sampled maternal families of Boechera stricta, an Arabidopsis thaliana relative, at one location in the central Rocky Mountains and grew them in three treatments that, based on records from an adjacent weather station, simulated hourly temperature changes at the native site during three summers with differing mean temperatures. Treatment had a significant effect on all traits, with 2–3‐fold increase in above‐ and belowground biomass and the highest allocation to roots observed in the treatment simulating the warmest summer on record at the site. Treatment affected bivariate associations between traits, with the weakest correlation between above‐ and belowground biomass in the warmest treatment. The magnitude of quantitative genetic variation for all traits differed across treatments: Genetic variance of biomass was 0 in the warmest treatment, while highly significant diversity was found in average conditions, resulting in broad‐sense heritability of 0.31. Significant genotype × environment interactions across all treatments were found only in root‐to‐shoot ratio. Therefore, temperature variation among summers appears unlikely to account for the observed levels of local genetic variation in size in this perennial species, but may influence family rank order in growth allocation. Our results indicate that natural environmental fluctuations can have a large impact on the magnitude of within‐population quantitative genetic variance.

Drosophila wing modularity revisited through a quantitative genetic approach

To predict the response of complex morphological structures to selection it is necessary to know how the covariation among its different parts is organized. Two key features of covariation are modularity and integration. The Drosophila wing is currently considered a fully integrated structure. Here, we study the patterns of integration of the Drosophila wing and test the hypothesis of the wing being divided into two modules along the proximo-distal axis, as suggested by developmental, biomechanical, and evolutionary evidence. To achieve these goals we perform a multilevel analysis of covariation combining the techniques of geometric morphometrics and quantitative genetics. Our results indicate that the Drosophila wing is indeed organized into two main modules, the wing base and the wing blade. The patterns of integration and modularity were highly concordant at the phenotypic, genetic, environmental, and developmental levels. Besides, we found that modularity at the developmental level was considerably higher than modularity at other levels, suggesting that in the Drosophila wing direct developmental interactions are major contributors to total phenotypic shape variation. We propose that the precise time at which covariance-generating developmental processes occur and/or the magnitude of variation that they produce favor proximo-distal, rather than anterior-posterior, modularity in the Drosophila wing.

Quantitative perturbation-based analysis of gene expression predicts enhancer activity in early Drosophila embryo

Enhancers constitute one of the major components of regulatory machinery of metazoans. Although several genome-wide studies have focused on finding and locating enhancers in the genomes, the fundamental principles governing their internal architecture and cis-regulatory grammar remain elusive. Here, we describe an extensive, quantitative perturbation analysis targeting the dorsal-ventral patterning gene regulatory network (GRN) controlled by Drosophila NF-κB homolog Dorsal. To understand transcription factor interactions on enhancers, we employed an ensemble of mathematical models, testing effects of cooperativity, repression, and factor potency. Models trained on the dataset correctly predict activity of evolutionarily divergent regulatory regions, providing insights into spatial relationships between repressor and activator binding sites. Importantly, the collective predictions of sets of models were effective at novel enhancer identification and characterization. Our study demonstrates how experimental dataset and modeling can be effectively combined to provide quantitative insights into cis-regulatory information on a genome-wide scale.

DOI:http://dx.doi.org/10.7554/eLife.08445.001

DNA contains regions known as genes, which may be “transcribed” to produce the RNA molecules that act as templates for building proteins and regulate cell activity. Proteins called transcription factors can bind to specific sequences of DNA to influence whether nearby genes are transcribed. For example, so-called enhancer regions of DNA contain several binding sites for transcription factors, and this binding activates gene transcription. Little is known about how the transcription factor binding sites are organized in enhancer regions, which makes it difficult to use DNA sequence information alone to predict the regulation of genes.

A transcription factor called Dorsal controls the activity of a network of genes that plays a crucial role in the development of fruit fly embryos. Dorsal binds to the enhancer region of a gene called rhomboid, which has been well studied and is known to be a fairly typical example of an enhancer region.

To understand the regulatory information encoded in the DNA sequences of enhancers, Sayal, Dresch et al. have now used a technique called perturbation analysis to investigate the interactions that are likely to occur between Dorsal and other transcription factors as they bind to the rhomboid enhancer. This technique involves systematically mutating the enhancer to remove different combinations of transcription factor binding sites and quantitatively investigating the effect this has on gene activity. A large set of mathematical models were then trained using this data and shown to correctly predict the activity of a range of other gene regulatory regions. The collective predictions of the models identified new enhancer regions and revealed details about how different types of transcription factor binding sites are arranged within enhancers.

As we enter an era where the DNA sequences of entire human populations are increasingly accessible, we would like to know the functional significance of changes in gene regulatory regions. Sayal, Dresch et al. show that the regulatory properties of specific control proteins are accessible by employing quantitative experiments and mathematical models. Similar studies will be required to learn how mutations found across the genome may alter gene expression, leading to better diagnosis and treatment of disease.

DOI:http://dx.doi.org/10.7554/eLife.08445.002

Armour GE, Mackay TFC, Anholt RRH. Epistatic partners of neurogenic genes modulate Drosophila olfactory behavior

The extent to which epistasis affects the genetic architecture of complex traits is difficult to quantify, and identifying variants in natural populations with epistatic interactions is challenging. Previous studies in Drosophila implicated extensive epistasis between variants in genes that affect neural connectivity and contribute to natural variation in olfactory response to benzaldehyde. In this study, we implemented a powerful screen to quantify the extent of epistasis as well as identify candidate interacting variants using 203 inbred wild-derived lines with sequenced genomes of the Drosophila melanogaster Genetic Reference Panel (DGRP). We crossed the DGRP lines to P[GT1]-element insertion mutants in Sema-5c and neuralized (neur), two neurodevelopmental loci which affect olfactory behavior, and to their coisogenic wild-type control. We observed significant variation in olfactory responses to benzaldehyde among F1 genotypes and for the DGRP line by mutant genotype interactions for both loci, showing extensive nonadditive genetic variation. We performed genome-wide association analyses to identify the candidate modifier loci. None of these polymorphisms were in or near the focal genes; therefore, epistasis is the cause of the nonadditive genetic variance. Candidate genes could be placed in interaction networks. Several candidate modifiers are associated with neural development. Analyses of mutants of candidate epistatic partners with neur (merry-go-round (mgr), prospero (pros), CG10098, Alhambra (Alh) and CG12535) and Sema-5c (CG42540 and bruchpilot (brp)) showed aberrant olfactory responses compared with coisogenic controls. Thus, integrating genome-wide analyses of natural variants with mutations at defined genomic locations in a common coisogenic background can unmask specific epistatic modifiers of behavioral phenotypes.

Quantitative Genetics of Food Intake in Drosophila melanogaster

Food intake is an essential animal activity, regulated by neural circuits that motivate food localization, evaluate nutritional content and acceptance or rejection responses through the gustatory system, and regulate neuroendocrine feedback loops that maintain energy homeostasis. Excess food consumption in people is associated with obesity and metabolic and cardiovascular disorders. However, little is known about the genetic basis of natural variation in food consumption. To gain insights in evolutionarily conserved genetic principles that regulate food intake, we took advantage of a model system, Drosophila melanogaster, in which food intake, environmental conditions and genetic background can be controlled precisely. We quantified variation in food intake among 182 inbred, sequenced lines of the Drosophila melanogaster Genetic Reference Panel (DGRP). We found significant genetic variation in the mean and within-line environmental variance of food consumption and observed sexual dimorphism and genetic variation in sexual dimorphism for both food intake traits (mean and variance). We performed genome wide association (GWA) analyses for mean food intake and environmental variance of food intake (using the coefficient of environmental variation, CV_E, as the metric for environmental variance) and identified molecular polymorphisms associated with both traits. Validation experiments using RNAi-knockdown confirmed 24 of 31 (77%) candidate genes affecting food intake and/or variance of food intake, and a test cross between selected DGRP lines confirmed a SNP affecting mean food intake identified in the GWA analysis. The majority of the validated candidate genes were novel with respect to feeding behavior, and many had mammalian orthologs implicated in metabolic diseases.

Genetic approaches in comparative and evolutionary physiology

Whole animal physiological performance is highly polygenic and highly plastic, and the same is generally true for the many subordinate traits that underlie performance capacities. Quantitative genetics, therefore, provides an appropriate framework for the analysis of physiological phenotypes and can be used to infer the microevolutionary processes that have shaped patterns of trait variation within and among species. In cases where specific genes are known to contribute to variation in physiological traits, analyses of intraspecific polymorphism and interspecific divergence can reveal molecular mechanisms of functional evolution and can provide insights into the possible adaptive significance of observed sequence changes. In this review, we explain how the tools and theory of quantitative genetics, population genetics, and molecular evolution can inform our understanding of mechanism and process in physiological evolution. For example, lab-based studies of polygenic inheritance can be integrated with field-based studies of trait variation and survivorship to measure selection in the wild, thereby providing direct insights into the adaptive significance of physiological variation. Analyses of quantitative genetic variation in selection experiments can be used to probe interrelationships among traits and the genetic basis of physiological trade-offs and constraints. We review approaches for characterizing the genetic architecture of physiological traits, including linkage mapping and association mapping, and systems approaches for dissecting intermediary steps in the chain of causation between genotype and phenotype. We also discuss the promise and limitations of population genomic approaches for inferring adaptation at specific loci. We end by highlighting the role of organismal physiology in the functional synthesis of evolutionary biology.

Exploring the Phenotypic Space and the Evolutionary History of a Natural Mutation in Drosophila melanogaster

A major challenge of modern Biology is elucidating the functional consequences of natural mutations. Although we have a good understanding of the effects of laboratory-induced mutations on the molecular- and organismal-level phenotypes, the study of natural mutations has lagged behind. In this work, we explore the phenotypic space and the evolutionary history of a previously identified adaptive transposable element insertion. We first combined several tests that capture different signatures of selection to show that there is evidence of positive selection in the regions flanking FBti0019386 insertion. We then explored several phenotypes related to known phenotypic effects of nearby genes, and having plausible connections to fitness variation in nature. We found that flies with FBti0019386 insertion had a shorter developmental time and were more sensitive to stress, which are likely to be the adaptive effect and the cost of selection of this mutation, respectively. Interestingly, these phenotypic effects are not consistent with a role of FBti0019386 in temperate adaptation as has been previously suggested. Indeed, a global analysis of the population frequency of FBti0019386 showed that climatic variables explain well the FBti0019386 frequency patterns only in Australia. Finally, although FBti0019386 insertion could be inducing the formation of heterochromatin by recruiting HP1a (Heterochromatin Protein 1a) protein, the insertion is associated with upregulation of sra in adult females. Overall, our integrative approach allowed us to shed light on the evolutionary history, the relevant fitness effects, and the likely molecular mechanisms of an adaptive mutation and highlights the complexity of natural genetic variants.

Power analysis of artificial selection experiments using efficient whole genome simulation of quantitative traits

Evolve and resequence studies combine artificial selection experiments with massively parallel sequencing technology to study the genetic basis for complex traits. In these experiments, individuals are selected for extreme values of a trait, causing alleles at quantitative trait loci (QTL) to increase or decrease in frequency in the experimental population. We present a new analysis of the power of artificial selection experiments to detect and localize quantitative trait loci. This analysis uses a simulation framework that explicitly models whole genomes of individuals, quantitative traits, and selection based on individual trait values. We find that explicitly modeling QTL provides qualitatively different insights than considering independent loci with constant selection coefficients. Specifically, we observe how interference between QTL under selection affects the trajectories and lengthens the fixation times of selected alleles. We also show that a substantial portion of the genetic variance of the trait (50–100%) can be explained by detected QTL in as little as 20 generations of selection, depending on the trait architecture and experimental design. Furthermore, we show that power depends crucially on the opportunity for recombination during the experiment. Finally, we show that an increase in power is obtained by leveraging founder haplotype information to obtain allele frequency estimates.

Gene expression variation in Drosophila melanogaster due to rare transposable element insertion alleles of large effect

Transposable elements are a common source of genetic variation that may play a substantial role in contributing to gene expression variation. However, the contribution of transposable elements to expression variation thus far consists of a handful of examples. We used previously published gene expression data from 37 inbred Drosophila melanogaster lines from the Drosophila Genetic Reference Panel to perform a genome-wide assessment of the effects of transposable elements on gene expression. We found thousands of transcripts with transposable element insertions in or near the transcript and that the presence of a transposable element in or near a transcript is significantly associated with reductions in expression. We estimate that within this example population, ∼2.2% of transcripts have a transposable element insertion, which significantly reduces expression in the line containing the transposable element. We also find that transcripts with insertions within 500 bp of the transcript show on average a 0.67 standard deviation decrease in expression level. These large decreases in expression level are most pronounced for transposable element insertions close to transcripts and the effect diminishes for more distant insertions. This work represents the first genome-wide analysis of gene expression variation due to transposable elements and suggests that transposable elements are an important class of mutation underlying expression variation in Drosophila and likely in other systems, given the ubiquity of these mobile elements in eukaryotic genomes.

The genotype-phenotype map of yeast complex traits: basic parameters and the role of natural selection

Most phenotypic traits are controlled by many genes, but a global picture of the genotype-phenotype map (GPM) is lacking. For example, in no species do we know generally how many genes affect a trait and how large these effects are. It is also unclear to what extent GPMs are shaped by natural selection. Here we address these fundamental questions using the reverse genetic data of 220 morphological traits in 4,718 budding yeast strains, each of which lacks a nonessential gene. We show that 1) the proportion of genes affecting a trait varies from <1% to >30%, averaging 6%, 2) most traits are impacted by many more small-effect genes than large-effect genes, and 3) the mean effect of all nonessential genes on a trait decreases precipitously as the estimated importance of the trait to fitness increases. An analysis of 3,116 yeast gene expression traits in 754 gene-deletion strains reveals a similar pattern. These findings illustrate the power of genome-wide reverse genetics in genotype-phenotype mapping, uncover an enormous range of genetic complexity of phenotypic traits, and suggest that the GPM of cellular organisms has been shaped by natural selection for mutational robustness.

The action of stabilizing selection, mutation, and drift on epistatic quantitative traits

For a quantitative trait under stabilizing selection, the effect of epistasis on its genetic architecture and on the changes of genetic variance caused by bottlenecking were investigated using theory and simulation. Assuming empirical estimates of the rate and effects of mutations and the intensity of selection, we assessed the impact of two-locus epistasis (synergistic/antagonistic) among linked or unlinked loci on the distribution of effects and frequencies of segregating loci in populations at the mutation-selection-drift balance. Strong pervasive epistasis did not modify substantially the genetic properties of the trait and, therefore, the most likely explanation for the low amount of variation usually accounted by the loci detected in genome-wide association analyses is that many causal loci will pass undetected. We investigated the impact of epistasis on the changes in genetic variance components when large populations were subjected to successive bottlenecks of different sizes, considering the action of genetic drift, operating singly (D), or jointly with mutation (MD) and selection (MSD). An initial increase of the different components of the genetic variance, as well as a dramatic acceleration of the between-line divergence, were always associated with synergistic epistasis but were strongly constrained by selection.

Stability and response of polygenic traits to stabilizing selection and mutation

When polygenic traits are under stabilizing selection, many different combinations of alleles allow close adaptation to the optimum. If alleles have equal effects, all combinations that result in the same deviation from the optimum are equivalent. Furthermore, the genetic variance that is maintained by mutation-selection balance is [Formula: see text] per locus, where μ is the mutation rate and S the strength of stabilizing selection. In reality, alleles vary in their effects, making the fitness landscape asymmetric and complicating analysis of the equilibria. We show that that the resulting genetic variance depends on the fraction of alleles near fixation, which contribute by [Formula: see text], and on the total mutational effects of alleles that are at intermediate frequency. The interplay between stabilizing selection and mutation leads to a sharp transition: alleles with effects smaller than a threshold value of [Formula: see text] remain polymorphic, whereas those with larger effects are fixed. The genetic load in equilibrium is less than for traits of equal effects, and the fitness equilibria are more similar. We find that if the optimum is displaced, alleles with effects close to the threshold value sweep first, and their rate of increase is bounded by [Formula: see text] Long-term response leads in general to well-adapted traits, unlike the case of equal effects that often end up at a suboptimal fitness peak. However, the particular peaks to which the populations converge are extremely sensitive to the initial states and to the speed of the shift of the optimum trait value.

Allelic diversity and its implications for the rate of adaptation

Genetic variation is usually estimated empirically from statistics based on population gene frequencies, but alternative statistics based on allelic diversity (number of allelic types) can provide complementary information. There is a lack of knowledge, however, on the evolutionary implications attached to allelic-diversity measures, particularly in structured populations. In this article we simulated multiple scenarios of single and structured populations in which a quantitative trait subject to stabilizing selection is adapted to different fitness optima. By forcing a global change in the optima we evaluated which diversity variables are more strongly correlated with both short- and long-term adaptation to the new optima. We found that quantitative genetic variance components for the trait and gene-frequency-diversity measures are generally more strongly correlated with short-term response to selection, whereas allelic-diversity measures are more correlated with long-term and total response to selection. Thus, allelic-diversity variables are better predictors of long-term adaptation than gene-frequency variables. This observation is also extended to unlinked neutral markers as a result of the information they convey on the demographic population history. Diffusion approximations for the allelic-diversity measures in a finite island model under the infinite-allele neutral mutation model are also provided.

Allelic diversity and its implications for the rate of adaptation

Genetic variation is usually estimated empirically from statistics based on population gene frequencies, but alternative statistics based on allelic diversity (number of allelic types) can provide complementary information. There is a lack of knowledge, however, on the evolutionary implications attached to allelic-diversity measures, particularly in structured populations. In this article we simulated multiple scenarios of single and structured populations in which a quantitative trait subject to stabilizing selection is adapted to different fitness optima. By forcing a global change in the optima we evaluated which diversity variables are more strongly correlated with both short- and long-term adaptation to the new optima. We found that quantitative genetic variance components for the trait and gene-frequency-diversity measures are generally more strongly correlated with short-term response to selection, whereas allelic-diversity measures are more correlated with long-term and total response to selection. Thus, allelic-diversity variables are better predictors of long-term adaptation than gene-frequency variables. This observation is also extended to unlinked neutral markers as a result of the information they convey on the demographic population history. Diffusion approximations for the allelic-diversity measures in a finite island model under the infinite-allele neutral mutation model are also provided.

Promises and limitations of hitchhiking mapping

Building the connection between genetic and phenotypic variation is an important 'work in progress', and one that will enable proactive diagnosis and treatment in medicine, promote development of environment-targeted varieties in agriculture, and clarify the limits of species adaptation to changing environments in conservation. Quantitative trait loci (QTL) mapping and genome wide association (GWA) studies have recently been allied to an additional focus on 'hitchhiking' (HH) mapping--using changes in allele frequency due to artificial or natural selection. This older technique has been popularized by the falling costs of high throughput sequencing. Initial HH-resequensing experiments seem to have found many thousands of polymorphisms responding to selection. We argue that this interpretation appears too optimistic, and that the data might in fact be more consistent with dozens, rather than thousands, of loci under selection. We propose several developments required for sensible data analyses that will fully realize the great power of the HH technique, and outline ways of moving forward.

Genetic complexity in a Drosophila model of diabetes-associated misfolded human proinsulin

Drosophila melanogaster has been widely used as a model of human Mendelian disease, but its value in modeling complex disease has received little attention. Fly models of complex disease would enable high-resolution mapping of disease-modifying loci and the identification of novel targets for therapeutic intervention. Here, we describe a fly model of permanent neonatal diabetes mellitus and explore the complexity of this model. The approach involves the transgenic expression of a misfolded mutant of human preproinsulin, hINS^C96Y, which is a cause of permanent neonatal diabetes. When expressed in fly imaginal discs, hINS^C96Y causes a reduction of adult structures, including the eye, wing, and notum. Eye imaginal discs exhibit defects in both the structure and the arrangement of ommatidia. In the wing, expression of hINS^C96Y leads to ectopic expression of veins and mechano-sensory organs, indicating disruption of wild-type signaling processes regulating cell fates. These readily measurable “disease” phenotypes are sensitive to temperature, gene dose, and sex. Mutant (but not wild-type) proinsulin expression in the eye imaginal disc induces IRE1-mediated XBP1 alternative splicing, a signal for endoplasmic reticulum stress response activation, and produces global change in gene expression. Mutant hINS transgene tester strains, when crossed to stocks from the Drosophila Genetic Reference Panel, produce F₁ adults with a continuous range of disease phenotypes and large broad-sense heritability. Surprisingly, the severity of mutant hINS-induced disease in the eye is not correlated with that in the notum in these crosses, nor with eye reduction phenotypes caused by the expression of two dominant eye mutants acting in two different eye development pathways, Drop (Dr) or Lobe (L), when crossed into the same genetic backgrounds. The tissue specificity of genetic variability for mutant hINS-induced disease has, therefore, its own distinct signature. The genetic dominance of disease-specific phenotypic variability in our model of misfolded human proinsulin makes this approach amenable to genome-wide association study in a simple F₁ screen of natural variation.

Genomics and the Evolution of Phenotypic Traits

Evolutionary genetics has entered an unprecedented era of discovery, catalyzed in large part by the development of technologies that provide information about genome sequence and function. An important benefit is the ability to move beyond a handful of model organisms in lab settings to identify the genetic basis for evolutionarily interesting traits in many organisms in natural settings. Other benefits are the abilities to identify causal mutations and validate their phenotypic consequences more readily and in many more species. Genomic technologies have reinvigorated interest in some of the most fundamental and persistent questions in evolutionary genetics, revealed previously unsuspected evolutionary phenomena, and opened the door to a wide range of new questions.

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.

Assumptions and properties of limiting pathway models for analysis of epistasis in complex traits

For most complex traits, results from genome-wide association studies show that the proportion of the phenotypic variance attributable to the additive effects of individual SNPs, that is, the heritability explained by the SNPs, is substantially less than the estimate of heritability obtained by standard methods using correlations between relatives. This difference has been called the "missing heritability". One explanation is that heritability estimates from family (including twin) studies are biased upwards. Zuk et al. revisited overestimation of narrow sense heritability from twin studies as a result of confounding with non-additive genetic variance. They propose a limiting pathway (LP) model that generates significant epistatic variation and its simple parametrization provides a convenient way to explore implications of epistasis. They conclude that over-estimation of narrow sense heritability from family data ('phantom heritability') may explain an important proportion of missing heritability. We show that for highly heritable quantitative traits large phantom heritability estimates from twin studies are possible only if a large contribution of common environment is assumed. The LP model is underpinned by strong assumptions that are unlikely to hold, including that all contributing pathways have the same mean and variance and are uncorrelated. Here, we relax the assumptions that underlie the LP model to be more biologically plausible. Together with theoretical, empirical, and pragmatic arguments we conclude that in outbred populations the contribution of additive genetic variance is likely to be much more important than the contribution of non-additive variance.

Abundance and distribution of transposable elements in two Drosophila QTL mapping resources

Here we present computational machinery to efficiently and accurately identify transposable element (TE) insertions in 146 next-generation sequenced inbred strains of Drosophila melanogaster. The panel of lines we use in our study is composed of strains from a pair of genetic mapping resources: the Drosophila Genetic Reference Panel (DGRP) and the Drosophila Synthetic Population Resource (DSPR). We identified 23,087 TE insertions in these lines, of which 83.3% are found in only one line. There are marked differences in the distribution of elements over the genome, with TEs found at higher densities on the X chromosome, and in regions of low recombination. We also identified many more TEs per base pair of intronic sequence and fewer TEs per base pair of exonic sequence than expected if TEs are located at random locations in the euchromatic genome. There was substantial variation in TE load across genes. For example, the paralogs derailed and derailed-2 show a significant difference in the number of TE insertions, potentially reflecting differences in the selection acting on these loci. When considering TE families, we find a very weak effect of gene family size on TE insertions per gene, indicating that as gene family size increases the number of TE insertions in a given gene within that family also increases. TEs are known to be associated with certain phenotypes, and our data will allow investigators using the DGRP and DSPR to assess the functional role of TE insertions in complex trait variation more generally. Notably, because most TEs are very rare and often private to a single line, causative TEs resulting in phenotypic differences among individuals may typically fail to replicate across mapping panels since individual elements are unlikely to segregate in both panels. Our data suggest that “burden tests” that test for the effect of TEs as a class may be more fruitful.

Exploiting gene expression variation to capture gene-environment interactions for disease

Gene-environment interactions have long been recognized as a fundamental concept in evolutionary, quantitative, and medical genetics. In the genomics era, study of how environment and genome interact to shape gene expression variation is relevant to understanding the genetic architecture of complex phenotypes. While genetic analysis of gene expression variation focused on main effects, little is known about the extent of interaction effects implicating regulatory variants and their consequences on transcriptional variation. Here we survey the current state of the concept of transcriptional gene-environment interactions and discuss its utility for mapping disease phenotypes in light of the insights gained from genome-wide association studies of gene expression.

Proteomic characterization of inbreeding-related cold sensitivity in Drosophila melanogaster

Inbreeding depression is a widespread phenomenon of central importance to agriculture, medicine, conservation biology and evolutionary biology. Although the population genetic principles of inbreeding depression are well understood, we know little about its functional genomic causes. To provide insight into the molecular interplay between intrinsic stress responses, inbreeding depression and temperature tolerance, we performed a proteomic characterization of a well-defined conditional inbreeding effect in a single line of Drosophila melanogaster, which suffers from extreme cold sensitivity and lethality. We identified 48 differentially expressed proteins in a conditional lethal line as compared to two control lines. These proteins were enriched for proteins involved in hexose metabolism, in particular pyruvate metabolism, and many were found to be associated with lipid particles. These processes can be linked to known cold tolerance mechanisms, such as the production of cryoprotectants, membrane remodeling and the build-up of energy reserves. We checked mRNA-expression of seven genes with large differential protein expression. Although protein expression poorly correlated with gene expression, we found a single gene (CG18067) that, after cold shock, was upregulated in the conditional lethal line both at the mRNA and protein level. Expression of CG18067 also increased in control flies after cold shock, and has previously been linked to cold exposure and chill coma recovery time. Many differentially expressed proteins in our study appear to be involved in cold tolerance in non-inbred individuals. This suggest the conditional inbreeding effect to be caused by misregulation of physiological cold tolerance mechanisms.

Mapping Wolbachia distributions in the adult Drosophila brain

The maternally inherited bacterium Wolbachia infects the germline of most arthropod species. Using Drosophila simulans and D. melanogaster, we demonstrate that localization of Wolbachia to the fat bodies and adult brain is likely also a conserved feature of Wolbachia infection. Examination of three Wolbachia strains (WMel , WRiv , WPop ) revealed that the bacteria preferentially concentrate in the central brain with low titres in the optic lobes. Distribution within regions of the central brain is largely determined by the Wolbachia strain, while the titre is influenced by both, the host species and the bacteria strain. In neurons of the central brain and ventral nerve cord, Wolbachia preferentially localizes to the neuronal cell bodies but not to axons. All examined Wolbachia strains are present intracellularly or in extracellular clusters, with the pathogenic WPop strain exhibiting the largest and most abundant clusters. We also discovered that 16 of 40 lines from the Drosophila Genetic Reference Panel are Wolbachia infected. Direct comparison of Wolbachia infected and cured lines from this panel reveals that differences in physiological traits (chill coma recovery, starvation, longevity) are partially due to host line influences. In addition, a tetracycline-induced increase in Drosophila longevity was detected many generations after treatment.