Fine-mapping nicotine resistance loci in Drosophila using a multiparent advanced generation inter-cross population

Effects of cocaine, nicotine, and marijuana exposure in Drosophila Melanogaster development: A systematic review and meta-analysis

Abuse-related drug usage is a public health issue. Drosophila melanogaster has been used as an animal model to study the biological effects of these psychoactive substances in preclinical studies. Our objective in this review is to evaluate the adverse effects produced by cocaine, nicotine, and marijuana during the development of D. melanogaster. We searched experimental studies in which D. melanogaster was exposed to these three psychoactive drugs in seven online databases up to January 2023. Two reviewers independently extracted the data. Fifty-one studies met eligibility criteria and were included in the data extraction: nicotine (n = 26), cocaine (n = 20), and marijuana (n = 5). Fifteen studies were eligible for meta-analysis. Low doses (∼0.6 mM) of nicotine increased locomotor activity in fruit flies, while high doses (≥3 mM) led to a decrease. Similarly, exposure to cocaine increased locomotor activity, resulting in decreased climbing response in D. melanogaster. Studies with exposure to marijuana did not present a profile for our meta-analysis. However, this drug has been less associated with locomotor changes, but alterations in body weight and fat content and changes in cardiac function. Our analyses have shown that fruit flies exposed to drugs of abuse during different developmental stages, such as larvae and adults, exhibit molecular, morphological, behavioral, and survival changes that are dependent on the dosage. These phenotypes resemble the adverse effects of psychoactive substances in clinical medicine.

Differential gene expression and microRNA profile in corpora allata-corpora cardiaca of Aedes aegypti mosquitoes with weak juvenile hormone signalling

The corpora allata-corpora cardiaca (CA-CC) is an endocrine gland complex that regulates mosquito development and reproduction through the synthesis of juvenile hormone (JH). Epoxidase (Epox) is a key enzyme in the production of JH. We recently utilized CRISPR/Cas9 to establish an epoxidase-deficient (epox-/-) Aedes aegypti line. The CA from epox-/- mutants do not synthesize epoxidated JH III but methyl farneosate (MF), a weak agonist of the JH receptor, and therefore have reduced JH signalling. Illumina sequencing was used to examine the differences in gene expression between the CA-CC from wild type (WT) and epox-/- adult female mosquitoes. From 18,034 identified genes, 317 were significantly differentially expressed. These genes are involved in many biological processes, including the regulation of cell proliferation and apoptosis, energy metabolism, and nutritional uptake. In addition, the same CA-CC samples were also used to examine the microRNA (miRNA) profiles of epox-/- and WT mosquitoes. A total of 197 miRNAs were detected, 24 of which were differentially regulated in epox-/- mutants. miRNA binding sites for these particular miRNAs were identified using an in silico approach; they target a total of 101 differentially expressed genes. Our results suggest that a lack of epoxidase, besides affecting JH synthesis, results in the diminishing of JH signalling that have significant effects on Ae. aegypti CA-CC transcriptome profiles, as well as its miRNA repertoire.

Complexities of recapitulating polygenic effects in natural populations: replication of genetic effects on wing shape in artificially selected and wild-caught populations of Drosophila melanogaster

Identifying the genetic architecture of complex traits is important to many geneticists, including those interested in human disease, plant and animal breeding, and evolutionary genetics. Advances in sequencing technology and statistical methods for genome-wide association studies have allowed for the identification of more variants with smaller effect sizes, however, many of these identified polymorphisms fail to be replicated in subsequent studies. In addition to sampling variation, this failure to replicate reflects the complexities introduced by factors including environmental variation, genetic background, and differences in allele frequencies among populations. Using Drosophila melanogaster wing shape, we ask if we can replicate allelic effects of polymorphisms first identified in a genome-wide association studies in three genes: dachsous, extra-macrochaete, and neuralized, using artificial selection in the lab, and bulk segregant mapping in natural populations. We demonstrate that multivariate wing shape changes associated with these genes are aligned with major axes of phenotypic and genetic variation in natural populations. Following seven generations of artificial selection along the dachsous shape change vector, we observe genetic differentiation of variants in dachsous and genomic regions containing other genes in the hippo signaling pathway. This suggests a shared direction of effects within a developmental network. We also performed artificial selection with the extra-macrochaete shape change vector, which is not a part of the hippo signaling network, but showed a largely shared direction of effects. The response to selection along the emc vector was similar to that of dachsous, suggesting that the available genetic diversity of a population, summarized by the genetic (co)variance matrix (G), influenced alleles captured by selection. Despite the success with artificial selection, bulk segregant analysis using natural populations did not detect these same variants, likely due to the contribution of environmental variation and low minor allele frequencies, coupled with small effect sizes of the contributing variants.

Complex interactions between nicotine and resveratrol in the Drosophila melanogaster wing spot test

Nicotine (NIC) and resveratrol (RES) are chemicals in tobacco and wine, respectively, that are widely consumed concurrently worldwide. NIC is an alkaloid known to be toxic, addictive and to produce oxidative stress, while RES is thought of as an antioxidant with putative health benefits. Oxidative stress can induce genotoxic damage, yet few studies have examined whether NIC is genotoxic in vivo. In vitro studies have shown that RES can ameliorate deleterious effects of NIC. However, RES has been reported to have both antioxidant and pro-oxidant effects, and an in vivo study reported that 0.011 mM RES was genotoxic. We used the Drosophila melanogaster wing spot test to determine whether NIC and RES, first individually and then in combination, were genotoxic and/or altered the cell division. We hypothesized that RES would modulate NIC’s effects. NIC was genotoxic in the standard (ST) cross in a concentration-independent manner, but not genotoxic in the high bioactivation (HB) cross. RES was not genotoxic in either the ST or HB cross at the concentrations tested. We discovered a complex interaction between NIC and RES. Depending on concentration, RES was protective of NIC’s genotoxic damage, RES had no interaction with NIC, or RES had an additive or synergistic effect, increasing NIC’s genotoxic damage. Most NIC, RES, and NIC/RES combinations tested altered the cell division in the ST and HB crosses. Because we used the ST and HB crosses, we demonstrated that genotoxicity and cell division alterations were modulated by the xenobiotic metabolism. These results provide evidence of NIC’s genotoxicity in vivo at specific concentrations. Moreover, NIC’s genotoxicity can be modulated by its interaction with RES in a complex manner, in which their interaction can lead to either increasing NIC’s damage or protecting against it.

•

Nicotine was genotoxic at specific concentrations in the Drosophila wing spot test.

•

Resveratrol protected against nicotine’s genotoxic effects at some concentrations.

•

Resveratrol increased nicotine’s genotoxicity at specific concentrations.

•

Nicotine and resveratrol have a complex interaction in vivo.

•

Studying chemicals in combination in vivo may uncover unexpected interactions.

Trade-offs between cost of ingestion and rate of intake drive defensive toxin use

Animals that ingest toxins can become unpalatable and even toxic to predators and parasites through toxin sequestration. Because most animals rapidly eliminate toxins to survive their ingestion, it is unclear how populations transition from susceptibility and toxin elimination to tolerance and accumulation as chemical defence emerges. Studies of chemical defence have generally focused on species with active toxin sequestration and target-site insensitivity mutations or toxin-binding proteins that permit survival without necessitating toxin elimination. Here, we investigate whether animals that presumably rely on toxin elimination for survival can use ingested toxins for defence. We use the A4 and A3 Drosophila melanogaster fly strains from the Drosophila Synthetic Population Resource (DSPR), which respectively possess high and low metabolic nicotine resistance among DSPR fly lines. We find that ingesting nicotine increased A4 but not A3 fly survival against Leptopilina heterotoma wasp parasitism. Further, we find that despite possessing genetic variants that enhance toxin elimination, A4 flies accrued more nicotine than A3 individuals, likely by consuming more medium. Our results suggest that enhanced toxin metabolism can allow greater toxin intake by offsetting the cost of toxin ingestion. Passive toxin accumulation that accompanies increased toxin intake may underlie the early origins of chemical defence.

The UDP-Glycosyltransferase Family in Drosophila melanogaster: Nomenclature Update, Gene Expression and Phylogenetic Analysis

UDP-glycosyltransferases (UGTs) are important conjugation enzymes found in all kingdoms of life, catalyzing a sugar conjugation with small lipophilic compounds and playing a crucial role in detoxification and homeostasis. The UGT gene family is defined by a signature motif in the C-terminal domain where the uridine diphosphate (UDP)-sugar donor binds. UGTs have been identified in a number of insect genomes over the last decade and much progress has been achieved in characterizing their expression patterns and molecular functions. Here, we present an update of the complete repertoire of UGT genes in Drosophila melanogaster and provide a brief overview of the latest research in this model insect. A total of 35 UGT genes are found in the D. melanogaster genome, localized to chromosomes 2 and 3 with a high degree of gene duplications on the chromosome arm 3R. All D. melanogaster UGT genes have now been named in FlyBase according to the unified UGT nomenclature guidelines. A phylogenetic analysis of UGT genes shows lineage-specific gene duplications. Analysis of anatomical and induced gene expression patterns demonstrate that some UGT genes are differentially expressed in various tissues or after environmental treatments. Extended searches of UGT orthologs from 18 additional Drosophila species reveal a diversity of UGT gene numbers and composition. The roles of Drosophila UGTs identified to date are briefly reviewed, and include xenobiotic metabolism, nicotine resistance, olfaction, cold tolerance, sclerotization, pigmentation, and immunity. Together, the updated genomic information and research overview provided herein will aid further research in this developing field.

Characterizing the genetic basis of copper toxicity in Drosophila reveals a complex pattern of allelic, regulatory, and behavioral variation

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.

Genomic and transcriptomic analyses in Drosophila suggest that the ecdysteroid kinase-like (EcKL) gene family encodes the 'detoxification-by-phosphorylation' enzymes of insects

Phosphorylation is a phase II detoxification reaction that, among animals, occurs near exclusively in insects, but the enzymes responsible have never been cloned or otherwise identified. We propose the hypothesis that members of the arthropod-specific ecdysteroid kinase-like (EcKL) gene family encode detoxicative kinases. To test this hypothesis, we annotated the EcKL gene family in 12 species of Drosophila and explored their evolution within the genus. Many ancestral EcKL clades are evolutionarily unstable and have experienced repeated gene gain and loss events, while others are conserved as single-copy orthologs. Leveraging multiple published gene expression datasets from D. melanogaster, and using the cytochrome P450s-a classical detoxification family-as a test case, we demonstrate relationships between xenobiotic induction, detoxification tissue-enriched expression and evolutionary instability in the EcKLs and the P450s. We devised a systematic method for identifying candidate detoxification genes in large gene families that is concordant with experimentally determined functions of P450 genes in D. melanogaster. Applying this method to the EcKLs suggested a significant proportion of these genes play roles in detoxification, and that the EcKLs may constitute a detoxification gene family in insects. Additionally, we estimate that between 11 and 16 uncharacterised D. melanogaster P450s are strong detoxification candidates. Lastly, we also found previously unreported genomic and transcriptomic variation in a number of EcKLs and P450s associated with toxic stress phenotypes using a targeted phenome-wide association study (PheWAS) approach in D. melanogaster, presenting multiple future avenues of research for detoxification genetics in this species.

A naturally-occurring 22-bp coding deletion in Ugt86Dd reduces nicotine resistance in Drosophila melanogaster

Objective

Segregating genetic variants contribute to the response to toxic, xenobiotic compounds, and identifying these causative sites can help describe the mechanisms underlying metabolism of toxic compounds. In previous work we implicated the detoxification gene Ugt86Dd in the genetic control of larval nicotine resistance in Drosophila melanogaster. Furthermore, we suggested that a naturally-occurring 22-bp deletion that leads to a stop codon in exon 2 of the gene markedly reduces resistance. Here we use homology directed CRISPR/Cas9 gene editing to specifically test this hypothesis.

Results

We edited chromosome three from an inbred strain named A4 which carries the insertion allele at Ugt86Dd, successfully generated four alleles carrying the 22-bp Ugt86Dd deletion, and substituted edited chromosomes back into the A4 background. The original A4 strain, and an un-edited control strain in the same A4 background, show no significant difference in egg-to-adult or larva-to-adult viability on either control media or nicotine-supplemented media, and only slightly delayed development in nicotine media. However, strains carrying the 22-bp deletion showed reduced viability in nicotine conditions, and significantly longer development. Our data strongly suggest that the naturally-occurring 22-bp insertion/deletion event in Ugt86Dd directly impacts variation in nicotine resistance in D. melanogaster.

Genome-Wide Association Study in Two Cohorts from a Multi-generational Mouse Advanced Intercross Line Highlights the Difficulty of Replication Due to Study-Specific Heterogeneity

There has been extensive discussion of the "Replication Crisis" in many fields, including genome-wide association studies (GWAS). We explored replication in a mouse model using an advanced intercross line (AIL), which is a multigenerational intercross between two inbred strains. We re-genotyped a previously published cohort of LG/J x SM/J AIL mice (F34; n = 428) using a denser marker set and genotyped a new cohort of AIL mice (F39-43; n = 600) for the first time. We identified 36 novel genome-wide significant loci in the F34 and 25 novel loci in the F39-43 cohort. The subset of traits that were measured in both cohorts (locomotor activity, body weight, and coat color) showed high genetic correlations, although the SNP heritabilities were slightly lower in the F39-43 cohort. For this subset of traits, we attempted to replicate loci identified in either F34 or F39-43 in the other cohort. Coat color was robustly replicated; locomotor activity and body weight were only partially replicated, which was inconsistent with our power simulations. We used a random effects model to show that the partial replications could not be explained by Winner's Curse but could be explained by study-specific heterogeneity. Despite this heterogeneity, we performed a mega-analysis by combining F34 and F39-43 cohorts (n = 1,028), which identified four novel loci associated with locomotor activity and body weight. These results illustrate that even with the high degree of genetic and environmental control possible in our experimental system, replication was hindered by study-specific heterogeneity, which has broad implications for ongoing concerns about reproducibility.

Structural variants exhibit widespread allelic heterogeneity and shape variation in complex traits

It has been hypothesized that individually-rare hidden structural variants (SVs) could account for a significant fraction of variation in complex traits. Here we identified more than 20,000 euchromatic SVs from 14 Drosophila melanogaster genome assemblies, of which ~40% are invisible to high specificity short-read genotyping approaches. SVs are common, with 31.5% of diploid individuals harboring a SV in genes larger than 5kb, and 24% harboring multiple SVs in genes larger than 10kb. SV minor allele frequencies are rarer than amino acid polymorphisms, suggesting that SVs are more deleterious. We show that a number of functionally important genes harbor previously hidden structural variants likely to affect complex phenotypes. Furthermore, SVs are overrepresented in candidate genes associated with quantitative trait loci mapped using the Drosophila Synthetic Population Resource. We conclude that SVs are ubiquitous, frequently constitute a heterogeneous allelic series, and can act as rare alleles of large effect.

UDP-glucosyltransferases potentially contribute to imidacloprid resistance in Aphis gossypii glover based on transcriptomic and proteomic analyses

The cotton aphid, Aphis gossypii Glover, is a destructive global crop pest. Control of A. gossypii has relied heavily on the application of chemical insecticides. The cotton aphid has developed resistance to numerous insecticides, including imidacloprid, which has been widely used to control cotton pests in China since the 1990s. Our objective was to investigate the potential role of UDP-glycosyltransferases (UGTs) in imidacloprid resistance based on transcriptomic and proteomic analyses of field-originated imidacloprid-resistant (IMI_R) and -susceptible (IMI_S) A. gossypii clones. The transcriptomic and proteomic analyses revealed that 12 out of 512 differentially expressed genes and three out of 510 differentially expressed proteins were predicted as UDP-glycosyltransferase (UGT). Based on quantitative real-time PCR analysis, nine UGT genes, UGT343A4, UGT344A15, UGT344A16, UGT344B4, UGT344C7, UGT344C9, UGT344N4, UGT 24541, and UGT7630, were up-regulated in the IMI_R clone compared to the IMI_S clone. Meanwhile, UGT344A16, UGT344B4, UGT344C7, and UGT344N4 were overexpressed at the protein level based on western blot analysis. Furthermore, knockdown of UGT344B4 or UGT344C7 using RNA interference (RNAi) significantly increased sensitivity to imidacloprid in the IMI_R clone. In conclusion, UGTs potentially contributed to imidacloprid resistance in A. gossypii originating from cotton-growing regions of China. These results provide insights into the way we study insecticide resistance in cotton aphids.

Multiple genetic loci affect place learning and memory performance in Drosophila melanogaster

Learning and memory are critical functions for all animals, giving individuals the ability to respond to changes in their environment. Within populations, individuals vary, however the mechanisms underlying this variation in performance are largely unknown. Thus, it remains to be determined what genetic factors cause an individual to have high learning ability and what factors determine how well an individual will remember what they have learned. To genetically dissect learning and memory performance, we used the Drosophila synthetic population resource (DSPR), a multiparent mapping resource in the model system Drosophila melanogaster, consisting of a large set of recombinant inbred lines (RILs) that naturally vary in these and other traits. Fruit flies can be trained in a "heat box" to learn to remain on one side of a chamber (place learning) and can remember this (place memory) over short timescales. Using this paradigm, we measured place learning and memory for ~49 000 individual flies from over 700 DSPR RILs. We identified 16 different loci across the genome that significantly affect place learning and/or memory performance, with 5 of these loci affecting both traits. To identify transcriptomic differences associated with performance, we performed RNA-Seq on pooled samples of seven high performing and seven low performing RILs for both learning and memory and identified hundreds of genes with differences in expression in the two sets. Integrating our transcriptomic results with the mapping results allowed us to identify nine promising candidate genes, advancing our understanding of the genetic basis underlying natural variation in learning and memory performance.

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

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

What can genetic association panels tell us about evolutionary processes in insects?

If we are to fully comprehend the evolution of insect diversity at a genomic level we need to understand how natural selection can alter genetically encoded characters within populations. Genetic association panels have the potential to be standard bearers in this endeavour. They enable the mapping of phenotypes to genotypes at unprecedented resolution while simultaneously providing population genomic samples that can be interrogated for the tell-tale signs of selection. Analyses of these panels promise to elucidate the entanglement of gene ontologies, pathways, developmental processes and evolutionary constraints, and inform how these are shaped by adaptation.

Shared Genomic Regions Underlie Natural Variation in Diverse Toxin Responses

Phenotypic complexity is caused by the contributions of environmental factors and multiple genetic loci, interacting or acting independently. Studies of yeast and Arabidopsis often find that the majority of natural variation across phenotypes is attributable to independent additive quantitative trait loci (QTL). Detected loci in these organisms explain most of the estimated heritable variation. By contrast, many heritable components underlying phenotypic variation in metazoan models remain undetected. Before the relative impacts of additive and interactive variance components on metazoan phenotypic variation can be dissected, high replication and precise phenotypic measurements are required to obtain sufficient statistical power to detect loci contributing to this missing heritability. Here, we used a panel of 296 recombinant inbred advanced intercross lines of Caenorhabditis elegans and a high-throughput fitness assay to detect loci underlying responses to 16 different toxins, including heavy metals, chemotherapeutic drugs, pesticides, and neuropharmaceuticals. Using linkage mapping, we identified 82 QTL that underlie variation in responses to these toxins, and predicted the relative contributions of additive loci and genetic interactions across various growth parameters. Additionally, we identified three genomic regions that impact responses to multiple classes of toxins. These QTL hotspots could represent common factors impacting toxin responses. We went further to generate near-isogenic lines and chromosome substitution strains, and then experimentally validated these QTL hotspots, implicating additive and interactive loci that underlie toxin-response variation.

A Fly's Eye View of Natural and Drug Reward

Animals encounter multiple stimuli each day. Some of these stimuli are innately appetitive or aversive, while others are assigned valence based on experience. Drugs like ethanol can elicit aversion in the short term and attraction in the long term. The reward system encodes the predictive value for different stimuli, mediating anticipation for attractive or punishing stimuli and driving animal behavior to approach or avoid conditioned stimuli. The neurochemistry and neurocircuitry of the reward system is partly evolutionarily conserved. In both vertebrates and invertebrates, including Drosophila melanogaster, dopamine is at the center of a network of neurotransmitters and neuromodulators acting in concert to encode rewards. Behavioral assays in D. melanogaster have become increasingly sophisticated, allowing more direct comparison with mammalian research. Moreover, recent evidence has established the functional modularity of the reward neural circuits in Drosophila. This functional modularity resembles the organization of reward circuits in mammals. The powerful genetic and molecular tools for D. melanogaster allow characterization and manipulation at the single-cell level. These tools are being used to construct a detailed map of the neural circuits mediating specific rewarding stimuli and have allowed for the identification of multiple genes and molecular pathways that mediate the effects of reinforcing stimuli, including their rewarding effects. This report provides an overview of the research on natural and drug reward in D. melanogaster, including natural rewards such as sugar and other food nutrients, and drug rewards including ethanol, cocaine, amphetamine, methamphetamine, and nicotine. We focused mainly on the known genetic and neural mechanisms underlying appetitive reward for sugar and reward for ethanol. We also include genes, molecular pathways, and neural circuits that have been identified using assays that test the palatability of the rewarding stimulus, the preference for the rewarding stimulus, or other effects of the stimulus that indicate how it can modify behavior. Commonalities between mechanisms of natural and drug reward are highlighted and future directions are presented, putting forward questions best suited for research using D. melanogaster as a model organism.

A Genomic Reference Panel for Drosophila serrata

Here we describe a collection of re-sequenced inbred lines of Drosophila serrata, sampled from a natural population situated deep within the species endemic distribution in Brisbane, Australia. D. serrata is a member of the speciose montium group whose members inhabit much of south east Asia and has been well studied for aspects of climatic adaptation, sexual selection, sexual dimorphism, and mate recognition. We sequenced 110 lines that were inbred via 17-20 generations of full-sib mating at an average coverage of 23.5x with paired-end Illumina reads. 15,228,692 biallelic SNPs passed quality control after being called using the Joint Genotyper for Inbred Lines (JGIL). Inbreeding was highly effective and the average levels of residual heterozygosity (0.86%) were well below theoretical expectations. As expected, linkage disequilibrium decayed rapidly, with r² dropping below 0.1 within 100 base pairs. With the exception of four closely related pairs of lines which may have been due to technical errors, there was no statistical support for population substructure. Consistent with other endemic populations of other Drosophila species, preliminary population genetic analyses revealed high nucleotide diversity and, on average, negative Tajima’s D values. A preliminary GWAS was performed on a cuticular hydrocarbon trait, 2-Me-C₂₈ revealing 4 SNPs passing Bonferroni significance residing in or near genes. One gene Cht9 may be involved in the transport of CHCs from the site of production (oenocytes) to the cuticle. Our panel will facilitate broader population genomic and quantitative genetic studies of this species and serve as an important complement to existing D. melanogaster panels that can be used to test for the conservation of genetic architectures across the Drosophila genus.

Hidden genetic variation shapes the structure of functional elements in Drosophila

Mutations that add, subtract, rearrange, or otherwise refashion genome structure often affect phenotypes, although the fragmented nature of most contemporary assemblies obscures them. To discover such mutations, we assembled the first new reference-quality genome of Drosophila melanogaster since its initial sequencing. By comparing this new genome to the existing D. melanogaster assembly, we created a structural variant map of unprecedented resolution and identified extensive genetic variation that has remained hidden until now. Many of these variants constitute candidates underlying phenotypic variation, including tandem duplications and a transposable element insertion that amplifies the expression of detoxification-related genes associated with nicotine resistance. The abundance of important genetic variation that still evades discovery highlights how crucial high-quality reference genomes are to deciphering phenotypes.

Polygenicity and Epistasis Underlie Fitness-Proximal Traits in the Caenorhabditis elegans Multiparental Experimental Evolution (CeMEE) Panel

Understanding the genetic basis of complex traits remains a major challenge in biology. Polygenicity, phenotypic plasticity, and epistasis contribute to phenotypic variance in ways that are rarely clear. This uncertainty can be problematic for estimating heritability, for predicting individual phenotypes from genomic data, and for parameterizing models of phenotypic evolution. Here, we report an advanced recombinant inbred line (RIL) quantitative trait locus mapping panel for the hermaphroditic nematode Caenorhabditis elegans, the C. elegans multiparental experimental evolution (CeMEE) panel. The CeMEE panel, comprising 507 RILs at present, was created by hybridization of 16 wild isolates, experimental evolution for 140-190 generations, and inbreeding by selfing for 13-16 generations. The panel contains 22% of single-nucleotide polymorphisms known to segregate in natural populations, and complements existing C. elegans mapping resources by providing fine resolution and high nucleotide diversity across > 95% of the genome. We apply it to study the genetic basis of two fitness components, fertility and hermaphrodite body size at time of reproduction, with high broad-sense heritability in the CeMEE. While simulations show that we should detect common alleles with additive effects as small as 5%, at gene-level resolution, the genetic architectures of these traits do not feature such alleles. We instead find that a significant fraction of trait variance, approaching 40% for fertility, can be explained by sign epistasis with main effects below the detection limit. In congruence, phenotype prediction from genomic similarity, while generally poor ([Formula: see text]), requires modeling epistasis for optimal accuracy, with most variance attributed to the rapidly evolving chromosome arms.

Naturally Segregating Variation at Ugt86Dd Contributes to Nicotine Resistance in Drosophila melanogaster

Identifying the sequence polymorphisms underlying complex trait variation is a key goal of genetics research, since knowing the precise causative molecular events allows insight into the pathways governing trait variation. Genetic analysis of complex traits in model systems regularly starts by constructing QTL maps, but generally fails to identify causative sequence polymorphisms. Previously we mapped a series of QTL contributing to resistance to nicotine in a Drosophila melanogaster multiparental mapping resource and here use a battery of functional tests to resolve QTL to the molecular level. One large-effect QTL resided over a cluster of UDP-glucuronosyltransferases, and quantitative complementation tests using deficiencies eliminating subsets of these detoxification genes revealed allelic variation impacting resistance. RNAseq showed that Ugt86Dd had significantly higher expression in genotypes that are more resistant to nicotine, and anterior midgut-specific RNA interference (RNAi) of this gene reduced resistance. We discovered a segregating 22-bp frameshift deletion in Ugt86Dd, and accounting for the InDel during mapping largely eliminates the QTL, implying the event explains the bulk of the effect of the mapped locus. CRISPR/Cas9 editing of a relatively resistant genotype to generate lesions in Ugt86Dd that recapitulate the naturally occurring putative loss-of-function allele, leads to a large reduction in resistance. Despite this major effect of the deletion, the allele appears to be very rare in wild-caught populations and likely explains only a small fraction of the natural variation for the trait. Nonetheless, this putatively causative coding InDel can be a launchpad for future mechanistic exploration of xenobiotic detoxification.

Loci Contributing to Boric Acid Toxicity in Two Reference Populations of Drosophila melanogaster

Populations maintain considerable segregating variation in the response to toxic, xenobiotic compounds. To identify variants associated with resistance to boric acid, a commonly-used household insecticide with a poorly understood mechanism of action, we assayed thousands of individuals from hundreds of strains. Using the Drosophila Synthetic Population Resource (DSPR), a multi-parental population (MPP) of inbred genotypes, we mapped six QTL to short genomic regions containing few protein-coding genes (3–188), allowing us to identify plausible candidate genes underlying resistance to boric acid toxicity. One interval contains multiple genes from the cytochrome P450 family, and we show that ubiquitous RNAi of one of these genes, Cyp9b2, markedly reduces resistance to the toxin. Resistance to boric acid is positively correlated with caffeine resistance. The two phenotypes additionally share a pair of QTL, potentially suggesting a degree of pleiotropy in the genetic control of resistance to these two distinct xenobiotics. Finally, we screened the Drosophila Genetic Reference Panel (DGRP) in an attempt to identify sequence variants within mapped QTL that are associated with boric acid resistance. The approach was largely unsuccessful, with only one QTL showing any associations at QTL-specific 20% False Discovery Rate (FDR) thresholds. Nonetheless, these associations point to a potential candidate gene that can be targeted in future validation efforts. Although the mapping data resulting from the two reference populations do not clearly overlap, our work provides a starting point for further genetic dissection of the processes underlying boric acid toxicity in insects.

α-amanitin resistance in Drosophila melanogaster: A genome-wide association approach

We investigated the mechanisms of mushroom toxin resistance in the Drosophila Genetic Reference Panel (DGRP) fly lines, using genome-wide association studies (GWAS). While Drosophila melanogaster avoids mushrooms in nature, some lines are surprisingly resistant to α-amanitin—a toxin found solely in mushrooms. This resistance may represent a pre-adaptation, which might enable this species to invade the mushroom niche in the future. Although our previous microarray study had strongly suggested that pesticide-metabolizing detoxification genes confer α-amanitin resistance in a Taiwanese D. melanogaster line Ama-KTT, none of the traditional detoxification genes were among the top candidate genes resulting from the GWAS in the current study. Instead, we identified Megalin, Tequila, and widerborst as candidate genes underlying the α-amanitin resistance phenotype in the North American DGRP lines, all three of which are connected to the Target of Rapamycin (TOR) pathway. Both widerborst and Tequila are upstream regulators of TOR, and TOR is a key regulator of autophagy and Megalin-mediated endocytosis. We suggest that endocytosis and autophagy of α-amanitin, followed by lysosomal degradation of the toxin, is one of the mechanisms that confer α-amanitin resistance in the DGRP lines.

Genetic dissection of intraspecific variation in a male-specific sexual trait in Drosophila melanogaster

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.

Characterization of Three Rice Multiparent Advanced Generation Intercross (MAGIC) Populations for Quantitative Trait Loci Identification

Three new rice ( L.) multiparent advanced generation intercross (MAGIC) populations were developed using eight elite rice varieties from different breeding programs. These three populations were two recombinant inbred line (RIL) populations derived from two 4-way crosses, DC1 and DC2, and one RIL population derived from an 8-way cross. These populations were genotyped using an Illumina Infinium rice 6K SNP chip. The potential of the three MAGIC populations in identifying marker-trait associations was demonstrated using the plant height (PH) and heading date (HD) measured in 2014. A population of 248 IRRI breeding lines and a population of 323 Chinese breeding lines were also included to compare genetic diversity and linkage disequilibrium (LD) pattern. Our study discovered that (i) the 8-way population had a higher gene diversity than the DC1, DC2, and IRRI populations; (ii) all three MAGIC populations showed no clear population structure; (iii) LD decayed to < 0.2 at about 2.5, 2.5, 1.25, 1.75, and 4.0 Mb for the DC1, DC2, 8-way, IRRI, and Chinese populations, respectively; and (iv) the 8-way population was more powerful than the DC1, DC2, and IRRI populations on QTL identification. The association analysis identified two and three QTL for PH and HD, respectively. Four of the five QTL had peak markers close to known genes. A novel QTL for PH was identified on chromosome 12 using the 8-way population. Therefore, our study suggests that the three new MAGIC populations are valuable resources for QTL identification.

Using Drosophila to study the evolution of herbivory and diet specialization

Herbivory evolved many times independently across the insect phylogeny, and its evolution is linked with increased rates of diversification. Plants present many barriers to potential herbivores, among them are the so-called secondary chemicals and other molecular defenses such as protease inhibitors that deter herbivores. To understand the mechanisms behind the emergence of herbivory and subsequent species radiations of insects driven largely by diet specialization, it is important to identify the molecular basis associated with these evolutionary transitions. However, most herbivore species lack the genomic information and genetic tools required to identify functionally important genes. The notable exception is the genus Drosophila in which herbivory evolved at least three times independently, and for which abundant genomic data are available. Furthermore, contained within the family Drosophilidae is Drosophila melanogaster, the first genetic model animal. Here, we provide a synthesis of the salient tools that the D. melanogaster system provides to identify functionally important genes required for herbivory and subsequent diet specialization across insects.

The Wiggle Index: An Open Source Bioassay to Assess Sub-Lethal Insecticide Response in Drosophila melanogaster

Toxicological assays measuring mortality are routinely used to describe insecticide response, but sub-lethal exposures to insecticides can select for resistance and yield additional biological information describing the ways in which an insecticide impacts the insect. Here we present the Wiggle Index (WI), a high-throughput method to quantify insecticide response by measuring the reduction in motility during sub-lethal exposures in larvae of the vinegar fly Drosophila melanogaster. A susceptible wild type strain was exposed to the insecticides chlorantraniliprole, imidacloprid, spinosad, and ivermectin. Each insecticide reduced larval motility, but response times and profiles differed among insecticides. Two sets of target site mutants previously identified in mortality studies on the basis of imidacloprid or spinosad resistance phenotypes were tested. In each case the resistant mutant responded significantly less than the control. The WI was also able to detect a spinosad response in the absence of the primary spinosad target site. This response was not detected in mortality assays suggesting that spinosad, like many other insecticides, may have secondary targets affecting behaviour. The ability of the WI to detect changes in insecticide metabolism was confirmed by overexpressing the imidacloprid metabolizing Cyp6g1 gene in digestive tissues or the central nervous system. The data presented here validate the WI as an inexpensive, generic, sub-lethal assay that can complement information gained from mortality assays, extending our understanding of the genetic basis of insecticide response in D. melanogaster.

Identifying Loci Contributing to Natural Variation in Xenobiotic Resistance in Drosophila

Natural populations exhibit a great deal of interindividual genetic variation in the response to toxins, exemplified by the variable clinical efficacy of pharmaceutical drugs in humans, and the evolution of pesticide resistant insects. Such variation can result from several phenomena, including variable metabolic detoxification of the xenobiotic, and differential sensitivity of the molecular target of the toxin. Our goal is to genetically dissect variation in the response to xenobiotics, and characterize naturally-segregating polymorphisms that modulate toxicity. Here, we use the Drosophila Synthetic Population Resource (DSPR), a multiparent advanced intercross panel of recombinant inbred lines, to identify QTL (Quantitative Trait Loci) underlying xenobiotic resistance, and employ caffeine as a model toxic compound. Phenotyping over 1,700 genotypes led to the identification of ten QTL, each explaining 4.5-14.4% of the broad-sense heritability for caffeine resistance. Four QTL harbor members of the cytochrome P450 family of detoxification enzymes, which represent strong a priori candidate genes. The case is especially strong for Cyp12d1, with multiple lines of evidence indicating the gene causally impacts caffeine resistance. Cyp12d1 is implicated by QTL mapped in both panels of DSPR RILs, is significantly upregulated in the presence of caffeine, and RNAi knockdown robustly decreases caffeine tolerance. Furthermore, copy number variation at Cyp12d1 is strongly associated with phenotype in the DSPR, with a trend in the same direction observed in the DGRP (Drosophila Genetic Reference Panel). No additional plausible causative polymorphisms were observed in a full genomewide association study in the DGRP, or in analyses restricted to QTL regions mapped in the DSPR. Just as in human populations, replicating modest-effect, naturally-segregating causative variants in an association study framework in flies will likely require very large sample sizes.

The plastic response of Manduca sexta to host and non-host plants

Specialist insect herbivores have evolved efficient ways to adapt to the major defenses of their host plants. Although Manduca sexta, specialized on Solanaceous plants, has become a model organism for insect molecular biology, little is known about its adaptive responses to the chemical defenses of its hosts. To study larval performance and transcriptomic responses to host and non-host plants, we conducted developmental assays and replicated RNAseq experiments with Manduca larvae fed on different Solanaceous plants as well as on a Brassicaceous non-host plant, Brassica napus. Manduca larvae developed fastest on Nicotiana attenuata, but no significant differences in performance were found on larvae fed on other Solanaceae or the non-host B. napus. The RNAseq experiments revealed that Manduca larvae display plastic responses at the gene expression level, and transcriptional signatures specific to the challenges of each host- and non-host plant. Our observations are not consistent with expectations that specialist herbivores would perform poorly on non-host plants. Instead, our findings demonstrate the ability of this specialized insect herbivore to efficiently use a larger repertoire of host plants than it utilizes in the field.

Bayesian modeling of haplotype effects in multiparent populations

A general Bayesian model, Diploffect, is described for estimating the effects of founder haplotypes at quantitative trait loci (QTL) detected in multiparental genetic populations; such populations include the Collaborative Cross (CC), Heterogeneous Socks (HS), and many others for which local genetic variation is well described by an underlying, usually probabilistically inferred, haplotype mosaic. Our aim is to provide a framework for coherent estimation of haplotype and diplotype (haplotype pair) effects that takes into account the following: uncertainty in haplotype composition for each individual; uncertainty arising from small sample sizes and infrequently observed haplotype combinations; possible effects of dominance (for noninbred subjects); genetic background; and that provides a means to incorporate data that may be incomplete or has a hierarchical structure. Using the results of a probabilistic haplotype reconstruction as prior information, we obtain posterior distributions at the QTL for both haplotype effects and haplotype composition. Two alternative computational approaches are supplied: a Markov chain Monte Carlo sampler and a procedure based on importance sampling of integrated nested Laplace approximations. Using simulations of QTL in the incipient CC (pre-CC) and Northport HS populations, we compare the accuracy of Diploffect, approximations to it, and more commonly used approaches based on Haley–Knott regression, describing trade-offs between these methods. We also estimate effects for three QTL previously identified in those populations, obtaining posterior intervals that describe how the phenotype might be affected by diplotype substitutions at the modeled locus.

Dissecting complex traits using the Drosophila Synthetic Population Resource

For most complex traits we have a poor understanding of the positions, phenotypic effects, and population frequencies of the underlying genetic variants contributing to their variation. Recently, several groups have developed multi-parent advanced intercross mapping panels in different model organisms in an attempt to improve our ability to characterize causative genetic variants. These panels are powerful and are particularly well suited to the dissection of phenotypic variation generated by rare alleles and loci segregating multiple functional alleles. We describe studies using one such panel, the Drosophila Synthetic Population Resource (DSPR), and the implications for our understanding of the genetic basis of complex traits. In particular, we note that many loci of large effect appear to be multiallelic. If multiallelism is a general rule, analytical approaches designed to identify multiallelic variants should be a priority for both genome-wide association studies (GWASs) and multi-parental panels.