Drosophila Genetic Reference Panel Research Articles

Laboratory diet influences cold tolerance in a genotype-dependent manner in Drosophila melanogaster

Cold stress can reduce insect fitness and is an important determinant of species distributions and responses to climate change. Cold tolerance is influenced by genotype and environmental conditions, with factors such as day length and temperature having a particularly strong influence. Recent studies also indicate that diet impacts cold tolerance, but it is unclear whether diet-mediated shifts in cold tolerance are consistent across distinct genotypes. The goal of this study was to determine the extent to which commonly used artificial diets influence cold tolerance in Drosophila melanogaster, and whether these effects are consistent across genetically distinct lines. Specifically, we tested the impact of different fly diets on 1) ability to survive cold stress, 2) critical thermal minimum (CTmin), and 3) the ability to maintain reproduction after cold stress. Experiments were conducted across six isogenic lines from the Drosophila Genetic Reference Panel, and these lines were reared on different fly diets. Cold shock survival, CTmin, and reproductive output pre- and post-cold exposure varied considerably across diet and genotype combinations, suggesting strong genotype by environment interactions shape nutritionally mediated changes in cold tolerance. For example, in some lines cold shock survival remained consistently high or low across diets, while in others cold shock survival ranged from 5% to 75% depending on diet. Ultimately, these results add to a growing literature that cold tolerance is shaped by complex interactions between genotype and environment and inform practical considerations when selecting a laboratory diet for thermal tolerance experiments in Drosophila.

Genetic basis of susceptibility to low-dose paraquat and variation between the sexes in Drosophila melanogaster.

Toxicant resistance is a complex trait, affected both by genetics and the environment. Like most complex traits, it can exhibit sexual dimorphism, yet sex is often overlooked as a factor in studies of toxicant resistance. Paraquat, one such toxicant, is a commonly used herbicide and is known to produce mitochondrial oxidative stress, decrease dopaminergic neurons and dopamine (DA) levels, and decrease motor ability. While the main effects of paraquat are well-characterized, less is known about the naturally occurring variation in paraquat susceptibility. The purpose of this study was to map the genes contributing to low-dose paraquat susceptibility in Drosophila melanogaster, and to determine if susceptibility differs between the sexes. One hundred of the Drosophila Genetic Reference Panel (DGRP) lines were scored for susceptibility via climbing ability and used in a genome-wide association study (GWAS). Variation in seventeen genes in females and thirty-five genes in males associated with paraquat susceptibility. Only two candidate genes overlapped between the sexes despite a significant positive correlation between male and female susceptibilities. Many associated polymorphisms had significant interactions with sex, with most having conditionally neutral effects. Conditional neutrality between the sexes probably stems from sex-biased expression which may result from partial resolution of sexual conflict. Candidate genes were verified with RNAi knockdowns, gene expression analyses, and DA quantification. Several of these genes are novel associations with paraquat susceptibility. This research highlights the importance of assessing both sexes when studying toxicant susceptibility.

GWAS reveal a role for the central nervous system in regulating weight and weight change in response to exercise

Body size and weight show considerable variation both within and between species. This variation is controlled in part by genetics, but also strongly influenced by environmental factors including diet and the level of activity experienced by the individual. Due to the increasing obesity epidemic in much of the world, there is considerable interest in the genetic factors that control body weight and how weight changes in response to exercise treatments. Here, we address this question in the Drosophila model system, utilizing 38 strains of the Drosophila Genetics Reference Panel. We use GWAS to identify the molecular pathways that control weight and weight changes in response to exercise. We find that there is a complex set of molecular pathways controlling weight, with many genes linked to the central nervous system (CNS). The CNS also plays a role in the weight change with exercise, in particular, signaling from the CNS. Additional analyses revealed that weight in Drosophila is driven by two factors, animal size, and body composition, as the amount of fat mass versus lean mass impacts the density. Thus, while the CNS appears to be important for weight and exercise-induced weight change, signaling pathways are particularly important for determining how exercise impacts weight.

Accurate, ultra-low coverage genome reconstruction and association studies in Hybrid Swarm mapping populations.

Genetic association studies seek to uncover the link between genotype and phenotype, and often utilize inbred reference panels as a replicable source of genetic variation. However, inbred reference panels can differ substantially from wild populations in their genotypic distribution, patterns of linkage-disequilibrium, and nucleotide diversity. As a result, associations discovered using inbred reference panels may not reflect the genetic basis of phenotypic variation in natural populations. To address this problem, we evaluated a mapping population design where dozens to hundreds of inbred lines are outbred for few generations, which we call the Hybrid Swarm. The Hybrid Swarm approach has likely remained underutilized relative to pre-sequenced inbred lines due to the costs of genome-wide genotyping. To reduce sequencing costs and make the Hybrid Swarm approach feasible, we developed a computational pipeline that reconstructs accurate whole genomes from ultra-low-coverage (0.05X) sequence data in Hybrid Swarm populations derived from ancestors with phased haplotypes. We evaluate reconstructions using genetic variation from the Drosophila Genetic Reference Panel as well as variation from neutral simulations. We compared the power and precision of Genome-Wide Association Studies using the Hybrid Swarm, inbred lines, recombinant inbred lines (RILs), and highly outbred populations across a range of allele frequencies, effect sizes, and genetic architectures. Our simulations show that these different mapping panels vary in their power and precision, largely depending on the architecture of the trait. The Hybrid Swam and RILs outperform inbred lines for quantitative traits, but not for monogenic ones. Taken together, our results demonstrate the feasibility of the Hybrid Swarm as a cost-effective method of fine-scale genetic mapping.

Prediction of complex phenotypes using the Drosophila melanogaster metabolome.

Understanding the genotype-phenotype map and how variation at different levels of biological organization is associated are central topics in modern biology. Fast developments in sequencing technologies and other molecular omic tools enable researchers to obtain detailed information on variation at DNA level and on intermediate endophenotypes, such as RNA, proteins and metabolites. This can facilitate our understanding of the link between genotypes and molecular and functional organismal phenotypes. Here, we use the Drosophila melanogaster Genetic Reference Panel and nuclear magnetic resonance (NMR) metabolomics to investigate the ability of the metabolome to predict organismal phenotypes. We performed NMR metabolomics on four replicate pools of male flies from each of 170 different isogenic lines. Our results show that metabolite profiles are variable among the investigated lines and that this variation is highly heritable. Second, we identify genes associated with metabolome variation. Third, using the metabolome gave better prediction accuracies than genomic information for four of five quantitative traits analyzed. Our comprehensive characterization of population-scale diversity of metabolomes and its genetic basis illustrates that metabolites have large potential as predictors of organismal phenotypes. This finding is of great importance, e.g., in human medicine, evolutionary biology and animal and plant breeding.

Natural genetic variation in Drosophila melanogaster reveals genes associated with Coxiella burnetii infection.

The gram-negative bacterium Coxiella burnetii is the causative agent of Query (Q) fever in humans and coxiellosis in livestock. Host genetics are associated with C. burnetii pathogenesis both in humans and animals; however, it remains unknown if specific genes are associated with severity of infection. We employed the Drosophila Genetics Reference Panel to perform a genome-wide association study to identify host genetic variants that affect host survival to C. burnetii infection. The genome-wide association study identified 64 unique variants (P < 10−5) associated with 25 candidate genes. We examined the role each candidate gene contributes to host survival during C. burnetii infection using flies carrying a null mutation or RNAi knockdown of each candidate. We validated 15 of the 25 candidate genes using at least one method. This is the first report establishing involvement of many of these genes or their homologs with C. burnetii susceptibility in any system. Among the validated genes, FER and tara play roles in the JAK/STAT, JNK, and decapentaplegic/TGF-β signaling pathways which are components of known innate immune responses to C. burnetii infection. CG42673 and DIP-ε play roles in bacterial infection and synaptic signaling but have no previous association with C. burnetii pathogenesis. Furthermore, since the mammalian ortholog of CG13404 (PLGRKT) is an important regulator of macrophage function, CG13404 could play a role in host susceptibility to C. burnetii through hemocyte regulation. These insights provide a foundation for further investigation regarding the genetics of C. burnetii susceptibility across a wide variety of hosts.

Testing Implications of the Omnigenic Model for the Genetic Analysis of Loci Identified through Genome-wide Association

Organismal phenotypes usually have a quantitative distribution, and their genetic architecture can be studied by genome-wide association (GWA) mapping approaches. In most of such studies, it has become clear that many genes of moderate or small effects contribute to the phenotype.1-4 Hence, the attention has turned toward the loci falling below the GWA cut-off, which may contribute to the phenotype through modifier interactions with a set of core genes, as proposed in the omnigenic model.5 One can thus predict that both moderate effect GWA-derived candidate genes and randomly chosen genes should have a similar likelihood to affect a given phenotype when they are analyzed via gene disruption assays. We have tested this hypothesis by using an automated phenotyping system for Drosophila pupal phenotypes.6,7 We first identified candidate genes for pupal length in a GWA based on the Drosophila Genetic Reference Panel (DGRP)8,9 and showed that most of these candidate genes are indeed involved in the phenotype. We then randomly chose genes below a GWA significance threshold and found that three-quarters of them had also an effect on the trait with comparable effect sizes as the GWA candidate genes. We further tested the effects of these knockout lines on an independent behavioral pupal trait (pupation site choice) and found that a similar fraction had a significant effect as well. Our data thus confirm the implication that a large number of genes can influence independent quantitative traits.

Beta-blockers reduce intestinal permeability and early mortality following traumatic brain injury in Drosophila.

Traumatic brain injury (TBI) frequently leads to non-neurological consequences such as intestinal permeability. The beta-blocker drug labetalol, which inhibits binding of catecholamine neurotransmitters to adrenergic receptors, reduces intestinal permeability in a rat TBI model. Using a Drosophila melanogaster TBI model, we previously found a strong positive correlation between intestinal permeability and mortality within 24 hours of TBI in a standard laboratory line (w1118) and across genetically diverse inbred lines from the Drosophila Genetic Reference Panel (DGRP). Here, we report that feeding injured w1118flies the beta-blockers labetalol and metoprolol reduced intestinal permeability and mortality. Additionally, metoprolol reduced intestinal permeability when 18 DGRP fly lines were analyzed in aggregate, but neither beta-blocker affected mortality. These data indicate that the mechanism underlying disruption of the intestinal barrier by adrenergic signaling following TBI is conserved between humans and flies and that mortality following TBI in flies is not strictly dependent on disruption of the intestinal barrier. Thus, the fly TBI model is useful for shedding light on the mechanism and consequences of adrenergic signaling between the brain and intestine following TBI in humans.

Leveraging Multiple Layers of Data To Predict Drosophila Complex Traits.

The ability to accurately predict complex trait phenotypes from genetic and genomic data are critical for the implementation of personalized medicine and precision agriculture; however, prediction accuracy for most complex traits is currently low. Here, we used data on whole genome sequences, deep RNA sequencing, and high quality phenotypes for three quantitative traits in the ∼200 inbred lines of the Drosophila melanogaster Genetic Reference Panel (DGRP) to compare the prediction accuracies of gene expression and genotypes for three complex traits. We found that expression levels (r = 0.28 and 0.38, for females and males, respectively) provided higher prediction accuracy than genotypes (r = 0.07 and 0.15, for females and males, respectively) for starvation resistance, similar prediction accuracy for chill coma recovery (null for both models and sexes), and lower prediction accuracy for startle response (r = 0.15 and 0.14 for female and male genotypes, respectively; and r = 0.12 and 0.11, for females and male transcripts, respectively). Models including both genotype and expression levels did not outperform the best single component model. However, accuracy increased considerably for all the three traits when we included gene ontology (GO) category as an additional layer of information for both genomic variants and transcripts. We found strongly predictive GO terms for each of the three traits, some of which had a clear plausible biological interpretation. For example, for starvation resistance in females, GO:0033500 (r = 0.39 for transcripts) and GO:0032870 (r = 0.40 for transcripts), have been implicated in carbohydrate homeostasis and cellular response to hormone stimulus (including the insulin receptor signaling pathway), respectively. In summary, this study shows that integrating different sources of information improved prediction accuracy and helped elucidate the genetic architecture of three Drosophila complex phenotypes.

Musashi expression in intestinal stem cells attenuates radiation-induced decline in intestinal permeability and survival in Drosophila

Exposure to genotoxic stress by environmental agents or treatments, such as radiation therapy, can diminish healthspan and accelerate aging. We have developed a Drosophila melanogaster model to study the molecular effects of radiation-induced damage and repair. Utilizing a quantitative intestinal permeability assay, we performed an unbiased GWAS screen (using 156 strains from the Drosophila Genetic Reference Panel) to search for natural genetic variants that regulate radiation-induced gut permeability in adult D. melanogaster. From this screen, we identified an RNA binding protein, Musashi (msi), as one of the possible genes associated with changes in intestinal permeability upon radiation. The overexpression of msi promoted intestinal stem cell proliferation, which increased survival after irradiation and rescued radiation-induced intestinal permeability. In summary, we have established D. melanogaster as an expedient model system to study the effects of radiation-induced damage to the intestine in adults and have identified msi as a potential therapeutic target.

Genetic architecture of natural variations of cardiac perfomances in flies

Drosophila, the simplest model system with a beating heart, has tremendous advantages for the identification and validation of genetic factors affecting heart development, function and aging. Our project aims to analyse the genetic architecture of the natural variation of heart function in Drosophila and to identify its main characteristics. In particular, we want to take advantage of this model system to disentangle the complex relationship between genetic and environmental factors influencing these complex traits. Cardiac performances were analysed in the Drosophila Genetic Reference Panel (DGRP), and Genome Wide Association Studies (GWAS) were conducted to identify variants associated to several cardiac traits (rate, rhythm, cardiac output, fractional shortening, etc). The compendium of variants and associated genes identified were studied using approaches based on the analysis of gene and protein interaction networks. Genetic and protein networks analyses pointed towards several signalling pathways involved, that were further validated in vivo by targeted gene knockdown. Phenotypic variance among genetically identical individuals was also analysed and associated variants were compared to those affecting variations of phenotypic mean, revealing that distinct though overlapping gene networks were affected. In addition, our results show that more than 90% of variants are located within non-coding regions, and we characterized and validated some aspects of the transcriptional regulatory network affected. We provide a comprehensive characterization of population scale diversity of cardiac traits and of their genetic basis. Our work may provide a framework to interpret similar data from highers eukaryotes, including humans.

Genetic Basis of Natural Variation in Spontaneous Grooming in Drosophila melanogaster.

Spontaneous grooming behavior is a component of insect fitness. We quantified spontaneous grooming behavior in 201 sequenced lines of the Drosophila melanogaster Genetic Reference Panel and observed significant genetic variation in spontaneous grooming, with broad-sense heritabilities of 0.25 and 0.24 in females and males, respectively. Although grooming behavior is highly correlated between males and females, we observed significant sex by genotype interactions, indicating that the genetic basis of spontaneous grooming is partially distinct in the two sexes. We performed genome-wide association analyses of grooming behavior, and mapped 107 molecular polymorphisms associated with spontaneous grooming behavior, of which 73 were in or near 70 genes and 34 were over 1 kilobase from the nearest gene. The candidate genes were associated with a wide variety of gene ontology terms, and several of the candidate genes were significantly enriched in a genetic interaction network. We performed functional assessments of 29 candidate genes using RNA interference, and found that 11 affected spontaneous grooming behavior. The genes associated with natural variation in Drosophila grooming are involved with glutamate metabolism (Gdh) and transport (Eaat); interact genetically with (CCKLR-17D1) or are in the same gene family as (PGRP-LA) genes previously implicated in grooming behavior; are involved in the development of the nervous system and other tissues; or regulate the Notch and Epidermal growth factor receptor signaling pathways. Several DGRP lines exhibited extreme grooming behavior. Excessive grooming behavior can serve as a model for repetitive behaviors diagnostic of several human neuropsychiatric diseases.

Genetic variation for resistance to the specific fly pathogen Entomophthora muscae

We found substantial variation in resistance to the fly-specific pathogen Entomophthora muscae 'Berkeley' (Entomophthoromycota), in 20 lines from the Drosophila melanogaster Genetic Reference Panel (DGRP). Resistance to E. muscae is positively (r = 0.55) correlated with resistance to the broad host range ascomycete entomopathogen Metarhizium anisopliae (Ma549), indicative of generalist (non-specific) defenses. Most of the lines showing above average resistance to Ma549 showed cross-resistance to E. muscae. However, lines that succumbed quickly to Ma549 exhibited the full range of resistance to E. muscae. This suggests fly populations differ in E. muscae-specific resistance mechanisms as well as generic defences effective against both Ma549 and E. muscae. We looked for trade-offs that could account for inter-line variation, but increases (decreases) in disease resistance to E. muscae are not consistently associated with increases (decreases) of resistance to oxidative stress, starvation stress and sleep indices. That these pathogens are dynamic agents of selection on hosts is reflected in this genetic variation for resistance in lines derived from wild populations.

Integrating GWAS and Transcriptomics to Identify the Molecular Underpinnings of Thermal Stress Responses in Drosophila melanogaster.

Thermal tolerance of an organism depends on both the ability to dynamically adjust to a thermal stress and preparatory developmental processes that enhance thermal resistance. However, the extent to which standing genetic variation in thermal tolerance alleles influence dynamic stress responses vs. preparatory processes is unknown. Here, using the model species Drosophila melanogaster, we used a combination of Genome Wide Association mapping (GWAS) and transcriptomic profiling to characterize whether genes associated with thermal tolerance are primarily involved in dynamic stress responses or preparatory processes that influence physiological condition at the time of thermal stress. To test our hypotheses, we measured the critical thermal minimum (CTmin) and critical thermal maximum (CTmax) of 100 lines of the Drosophila Genetic Reference Panel (DGRP) and used GWAS to identify loci that explain variation in thermal limits. We observed greater variation in lower thermal limits, with CTmin ranging from 1.81 to 8.60°C, while CTmax ranged from 38.74 to 40.64°C. We identified 151 and 99 distinct genes associated with CTmin and CTmax, respectively, and there was strong support that these genes are involved in both dynamic responses to thermal stress and preparatory processes that increase thermal resistance. Many of the genes identified by GWAS were involved in the direct transcriptional response to thermal stress (72/151 for cold; 59/99 for heat), and overall GWAS candidates were more likely to be differentially expressed than other genes. Further, several GWAS candidates were regulatory genes that may participate in the regulation of stress responses, and gene ontologies related to development and morphogenesis were enriched, suggesting many of these genes influence thermal tolerance through effects on development and physiological status. Overall, our results suggest that thermal tolerance alleles can influence both dynamic plastic responses to thermal stress and preparatory processes that improve thermal resistance. These results also have utility for directly comparing GWAS and transcriptomic approaches for identifying candidate genes associated with thermal tolerance.

GWAS for Lifespan and Decline in Climbing Ability in Flies upon Dietary Restriction Reveal decima as a Mediator of Insulin-like Peptide Production

Dietary restriction (DR) is the most robust means to extend lifespan and delay age-related diseases across species. An underlying assumption in the aging field is that DR enhances both lifespan and physical activity through similar mechanisms, but this has not been rigorously tested in different genetic backgrounds. Furthermore, nutrient response genes responsible for lifespan extension or age-related decline in functionality remain underexplored in natural populations. To address this, we measured nutrient-dependent changes in lifespan and age-related decline in climbing ability in the Drosophila Genetic Reference Panel fly strains. On average, DR extended lifespan and delayed decline in climbing ability, but there was a lack of correlation between these traits across individual strains, suggesting that distinct genetic factors modulate these traits independently and that genotype determines response to diet. Only 50% of strains showed positive response to DR for both lifespan and climbing ability, 14% showed a negative response for one trait but not both, and 35% showed no change in one or both traits. Through GWAS, we uncovered a number of genes previously not known to be diet responsive nor to influence lifespan or climbing ability. We validated decima as a gene that alters lifespan and daedalus as one that influences age-related decline in climbing ability. We found that decima influences insulin-like peptide transcription in the GABA receptor neurons downstream of short neuropeptide F precursor (sNPF) signaling. Modulating these genes produced independent effects on lifespan and physical activity decline, which suggests that these age-related traits can be regulated through distinct mechanisms.

Phenotypic Responses to and Genetic Architecture of Sterility Following Exposure to Sub-Lethal Temperature During Development.

Thermal tolerance range, based on temperatures that result in incapacitating effects, influences species’ distributions and has been used to predict species’ response to increasing temperature. Reproductive performance may also be negatively affected at less extreme temperatures, but such sublethal heat-induced sterility has been relatively ignored in studies addressing the potential effects of, and ability of species’ to respond to, predicted climate warming. The few studies examining the link between increased temperature and reproductive performance typically focus on adults, although effects can vary between life history stages. Here we assessed how sublethal heat stress during development impacted subsequent adult fertility and its plasticity, both of which can provide the raw material for evolutionary responses to increased temperature. We quantified phenotypic and genetic variation in fertility of Drosophila melanogaster reared at standardized densities in three temperatures (25, 27, and 29°C) from a set of lines of the Drosophila Genetic Reference Panel (DGRP). We found little phenotypic variation at the two lower temperatures with more variation at the highest temperature and for plasticity. Males were more affected than females. Despite reasonably large broad-sense heritabilities, a genome-wide association study found little evidence for additive genetic variance and no genetic variants were robustly linked with reproductive performance at specific temperatures or for phenotypic plasticity. We compared results on heat-induced male sterility with other DGRP results on relevant fitness traits measured after abiotic stress and found an association between male susceptibility to sterility and male lifespan reduction following oxidative stress. Our results suggest that sublethal stress during development has profound negative consequences on male adult reproduction, but despite phenotypic variation in a population for this response, there is limited evolutionary potential, either through adaptation to a specific developmental temperature or plasticity in response to developmental heat-induced sterility.

Ageing and genetic background influence anaesthetic effects in a D. melanogaster model of blunt trauma with brain injury†

BackgroundGeneral anaesthetics interact with the pathophysiological mechanisms of traumatic brain injury (TBI). We used a Drosophila melanogaster (fruit fly) model to test the hypothesis that ageing and genetic background modulate the effect of anaesthetics and hyperoxia on TBI-induced mortality in the context of blunt trauma. MethodsWe exposed flies to isoflurane or sevoflurane under normoxic or hyperoxic conditions and TBI, and subsequently quantified the effect on mortality 24 h after injury. To determine the effect of age on anaesthetic-induced mortality, we analysed flies at 1–8 and 43–50 days old. To determine the effect of genetic background, we performed a genome-wide association study (GWAS) analysis on a collection of young inbred, fully sequenced lines. ResultsExposure to anaesthetics and hyperoxia differentially affected mortality in young and old flies. Pre-exposure of young but not old flies to anaesthetics reduced mortality. Post-exposure selectively increased mortality. For old but not young flies, hyperoxia enhanced the effect on mortality of post-exposure to isoflurane but not to sevoflurane. Post-exposure to isoflurane in hyperoxia increased the mortality of young fly lines in the Drosophila Genetic Reference Panel collection to different extents. GWAS analysis of these data identified single nucleotide polymorphisms in genes involved in cell water regulation and oxygen sensing as being associated with the post-exposure effect on mortality. ConclusionsAgeing and genetic background influence the effects of volatile general anaesthetics and hyperoxia on mortality in the context of traumatic brain injury. Polymorphisms in specific genes are identified as potential causes of ageing and genetic effects.

The metabolome as a link in the genotype-phenotype map for peroxide resistance in the fruit fly, Drosophila melanogaster

BackgroundGenetic association studies that seek to explain the inheritance of complex traits typically fail to explain a majority of the heritability of the trait under study. Thus, we are left with a gap in the map from genotype to phenotype. Several approaches have been used to fill this gap, including those that attempt to map endophenotype such as the transcriptome, proteome or metabolome, that underlie complex traits. Here we used metabolomics to explore the nature of genetic variation for hydrogen peroxide (H2O2) resistance in the sequenced inbred Drosophila Genetic Reference Panel (DGRP).ResultsWe first studied genetic variation for H2O2 resistance in 179 DGRP lines and along with identifying the insulin signaling modulator u-shaped and several regulators of feeding behavior, we estimate that a substantial amount of phenotypic variation can be explained by a polygenic model of genetic variation. We then profiled a portion of the aqueous metabolome in subsets of eight ‘high resistance’ lines and eight ‘low resistance’ lines. We used these lines to represent collections of genotypes that were either resistant or sensitive to the stressor, effectively modeling a discrete trait. Across the range of genotypes in both populations, flies exhibited surprising consistency in their metabolomic signature of resistance. Importantly, the resistance phenotype of these flies was more easily distinguished by their metabolome profiles than by their genotypes. Furthermore, we found a metabolic response to H2O2 in sensitive, but not in resistant genotypes. Metabolomic data further implicated at least two pathways, glycogen and folate metabolism, as determinants of sensitivity to H2O2. We also discovered a confounding effect of feeding behavior on assays involving supplemented food.ConclusionsThis work suggests that the metabolome can be a point of convergence for genetic variation influencing complex traits, and can efficiently elucidate mechanisms underlying trait variation.

Rapid evolution of piRNA-mediated silencing of an invading transposable element was driven by abundant de novo mutations.

The regulation of transposable element (TE) activity by small RNAs is a ubiquitous feature of germlines. However, despite the obvious benefits to the host in terms of ensuring the production of viable gametes and maintaining the integrity of the genomes they carry, it remains controversial whether TE regulation evolves adaptively. We examined the emergence and evolutionary dynamics of repressor alleles after P-elements invaded the Drosophila melanogaster genome in the mid-twentieth century. In many animals including Drosophila, repressor alleles are produced by transpositional insertions into piRNA clusters, genomic regions encoding the Piwi-interacting RNAs (piRNAs) that regulate TEs. We discovered that ∼94% of recently collected isofemale lines in the Drosophila melanogaster Genetic Reference Panel (DGRP) contain at least one P-element insertion in a piRNA cluster, indicating that repressor alleles are produced by de novo insertion at an exceptional rate. Furthermore, in our sample of approximately 200 genomes, we uncovered no fewer than 80 unique P-element insertion alleles in at least 15 different piRNA clusters. Finally, we observe no footprint of positive selection on P-element insertions in piRNA clusters, suggesting that the rapid evolution of piRNA-mediated repression in D. melanogaster was driven primarily by mutation. Our results reveal for the first time how the unique genetic architecture of piRNA production, in which numerous piRNA clusters can encode regulatory small RNAs upon transpositional insertion, facilitates the nonadaptive rapid evolution of repression.

Genetic Variation and Potential for Resistance Development to the tTA Overexpression Lethal System in Insects.

Release of insect pests carrying the dominant lethal tetracycline transactivator (tTA) overexpression system has been proposed as a means for population suppression. High levels of the tTA transcription factor are thought to be toxic due to either transcriptional squelching or interference with protein ubiquitination. Here we utilized the Drosophila melanogaster Genetic Reference Panel (DGRP) to examine the influence of genetic variation on the efficacy of a female-specific tTA overexpression system. The level of female lethality between DGRP lines varied from 11 to 97% with a broad sense heritability of 0.89. A genome-wide association analysis identified 192 allelic variants associated with high or low lethality (P < 10−5), although none were significant when corrected for multiple testing. 151 of the variants fell within 108 genes that were associated with several biological processes including transcription and protein ubiquitination. In four lines with high female lethality, tTA RNA levels were similar or higher than in the parental tTA overexpression strain. In two lines with low lethality, tTA levels were about two fold lower than in the parental strain. However, in two other lines with low lethality, tTA levels were similar or approximately 30% lower. RNAseq analysis identified genes that were up or downregulated in the four low female lethal lines compared to the four high lethal lines. For example, genes associated with RNA processing and rRNA maturation were significantly upregulated in low lethal lines. Our data suggest that standing genetic variation in an insect population could provide multiple mechanisms for resistance to the tTA overexpression system.