A range of heavy metals are required for normal cell function and homeostasis. However, the anthropogenic release of metal compounds into soil and water sources presents a pervasive health threat. Copper is one of many heavy metals that negatively impacts diverse organisms at a global scale. Using a combination of quantitative trait locus (QTL) mapping and RNA sequencing in the Drosophila Synthetic Population Resource, we demonstrate that resistance to the toxic effects of ingested copper in D. melanogaster is genetically complex and influenced by allelic and expression variation at multiple loci. QTL mapping identified several QTL that account for a substantial fraction of heritability. Additionally, we find that copper resistance is impacted by variation in behavioral avoidance of copper and may be subject to life-stage specific regulation. Gene expression analysis further demonstrated that resistant and sensitive strains are characterized by unique expression patterns. Several of the candidate genes identified via QTL mapping and RNAseq have known copper-specific functions (e.g., Ccs, Sod3, CG11825), and others are involved in the regulation of other heavy metals (e.g., Catsup, whd). We validated several of these candidate genes with RNAi suggesting they contribute to variation in adult copper resistance. Our study illuminates the interconnected roles that allelic and expression variation, organism life stage, and behavior play in copper resistance, allowing a deeper understanding of the diverse mechanisms through which metal pollution can negatively impact organisms.

Despite the value of Recombinant Inbred Lines (RILs) for the dissection of complex traits, large panels can be difficult to maintain, distribute, and phenotype. An attractive alternative to RILs for many traits leverages selecting phenotypically extreme individuals from a segregating population, and subjecting pools of selected and control individuals to sequencing. Under a bulked or extreme segregant analysis paradigm, genomic regions contributing to trait variation are revealed as frequency differences between pools. Here we describe such an extreme quantitative trait locus, or X-QTL, mapping strategy that builds on an existing multiparental population, the DSPR (Drosophila Synthetic Population Resource), and involves phenotyping and genotyping a population derived by mixing hundreds of DSPR RILs. Simulations demonstrate that challenging, yet experimentally tractable X-QTL designs ( > =4 replicates, > =5000 individuals/replicate, and selecting the 5-10% most extreme animals) yield at least the same power as traditional RIL-based QTL mapping and can localize variants with sub-centimorgan resolution. We empirically demonstrate the effectiveness of the approach using a 4-fold replicated X-QTL experiment that identifies 7 QTL for caffeine resistance. Two mapped X-QTL factors replicate loci previously identified in RILs, 6/7 are associated with excellent candidate genes, and RNAi knock-downs support the involvement of 4 genes in the genetic control of trait variation. For many traits of interest to drosophilists, a bulked phenotyping/genotyping X-QTL design has considerable advantages.

1A range of heavy metals are required for normal cell function and homeostasis. Equally, 2 the anthropogenic release of heavy metals into soil and water sources presents a pervasive 3 health threat. Copper is one such metal; it functions as a critical enzymatic cofactor, but at high 4 concentrations is toxic, and can lead to the production of reactive oxygen species. Using a 5 combination of quantitative trait locus (QTL) mapping and RNA sequencing in the Drosophila 6Synthetic Population Resource (DSPR), we demonstrate that resistance to the toxic effects of 7 ingested copper in D. melanogaster is genetically complex, and influenced by allelic and 8 expression variation at multiple loci. Additionally, we find that copper resistance is impacted by 9 variation in behavioral avoidance of copper and may be subject to life-stage specific regulation. 10Multiple genes with known copper-specific functions, as well as genes that are involved in the 11 regulation of other heavy metals were identified as potential candidates to contribute to 12 variation in adult copper resistance. We demonstrate that nine of 16 candidates tested by RNAi 13 knockdown influence adult copper resistance, a number of which may have pleiotropic effects 14 since they have previously been shown to impact the response to other metals. Our work 15 provides new understanding of the genetic complexity of copper resistance, highlighting the 16 diverse mechanisms through which copper pollution can negatively impact organisms. 17Additionally, we further support the similarities between copper metabolism and that of other 18 essential and nonessential heavy metals. 19 20 Rearing and assay conditions 126Strains from the DSPR were maintained, reared, and tested in the same incubator under 127 a 12:12hr light:dark photoperiod at 25°C and 50% humidity. To obtain female flies for the adult 128 copper response assay, RNA sequencing, and RNAi validation, adults were transferred to 129 cornmeal-molasses-yeast food, allowed to oviposit for two days, then discarded. Experimental 130 female, presumably mated, flies from the following generation were sorted over CO2 and were 131 placed into vials with new cornmeal-molasses-yeast media for 24 hours prior to the start of each 132 assay before they were transferred to copper-supplemented food. All adult assays were 133 performed on 3-5 day old individuals. 134

An open question in evolutionary biology is the relationship between standing variation for a trait and the variation that leads to interspecific divergence. By identifying loci underlying phenotypic variation in intra- and interspecific crosses we can determine the extent to which polymorphism and divergence are controlled by the same genomic regions. Sexual traits provide abundant examples of morphological and behavioral diversity within and among species, and here we leverage variation in the Drosophila sex comb to address this question. The sex comb is an array of modified bristles or ‘teeth' present on the male forelegs of several Drosophilid species. Males use the comb to grasp females during copulation, and ablation experiments have shown that males lacking comb teeth typically fail to mate. We measured tooth number in >700 genotypes derived from a multiparental advanced-intercross population, mapping three moderate-effect loci contributing to trait heritability. Two quantitative trait loci (QTLs) coincide with previously identified intra- and interspecific sex comb QTL, but such overlap can be explained by chance alone, in part because of the broad swathes of the genome implicated by earlier, low-resolution QTL scans. Our mapped QTL regions encompass 70–124 genes, but do not include those genes known to be involved in developmental specification of the comb. Nonetheless, we identified plausible candidates within all QTL intervals, and used RNA interference to validate effects at four loci. Notably, TweedleS expression knockdown substantially reduces tooth number. The genes we highlight are strong candidates to harbor segregating, functional variants contributing to sex comb tooth number.

Despite the value of Recombinant Inbred Lines (RILs) for the dissection of complex traits, large panels can be difficult to maintain, distribute, and phenotype. An attractive alternative to RILs for many traits leverages selecting phenotypically-extreme individuals from a segregating population, and subjecting pools of selected and control individuals to sequencing. Under a bulked or extreme segregant analysis paradigm, genomic regions contributing to trait variation are revealed as frequency differences between pools. Here we describe such an extreme quantitative trait locus, or X-QTL mapping strategy that builds on an existing multiparental population, the DSPR (Drosophila Synthetic Population Resource), and involves phenotyping and genotyping a population derived by mixing hundreds of DSPR RILs. Simulations demonstrate that challenging, yet experimentally tractable X-QTL designs (>=4 replicates, >=5000 individuals/replicate, and a selection intensity of 5-10%) yield at least the same power as traditional RIL-based QTL mapping, and can localize variants with sub-centimorgan resolution. We empirically demonstrate the effectiveness of the approach using a 4-fold replicated X-QTL experiment that identifies 7 QTL for caffeine resistance. Two mapped X-QTL factors replicate loci previously identified in RILs, 6/7 are associated with excellent candidate genes, and RNAi knock-downs support the involvement of 4 genes in the genetic control of trait variation. For many traits of interest to drosophilists a bulked phenotyping/genotyping X-QTL design has considerable advantages.