Genetic basis of variation in cocaine and methamphetamine consumption in outbred populations of Drosophila melanogaster.

We used Drosophila melanogaster to map the genetic basis of naturally occurring variation in voluntary consumption of cocaine and methamphetamine. We derived an outbred advanced intercross population (AIP) from 37 sequenced inbred wild-derived lines of the Drosophila melanogaster Genetic Reference Panel (DGRP), which are maximally genetically divergent, have minimal residual heterozygosity, are not segregating for common inversions, and are not infected with Wolbachia pipientis We assessed consumption of sucrose, methamphetamine-supplemented sucrose, and cocaine-supplemented sucrose and found considerable phenotypic variation for consumption of both drugs, in both sexes. We performed whole-genome sequencing and extreme quantitative trait locus (QTL) mapping on the top 10% of consumers for each replicate, sex, and condition and an equal number of randomly selected flies. We evaluated changes in allele frequencies among high consumers and control flies and identified 3,033 variants significantly (P < 1.9 × 10-8) associated with increased consumption, located in or near 1,962 genes. Many of these genes are associated with nervous system development and function, and 77 belong to a known gene-gene interaction subnetwork. We assessed the effects of RNA interference (RNAi) on drug consumption for 22 candidate genes; 17 had a significant effect in at least one sex. We constructed allele-specific AIPs that were homozygous for alternative candidate alleles for 10 single-nucleotide polymorphisms (SNPs) and measured average consumption for each population; 9 SNPs had significant effects in at least one sex. The genetic basis of voluntary drug consumption in Drosophila is polygenic and implicates genes with human orthologs and associated variants with sex- and drug-specific effects.

Gene expression networks in the Drosophila Genetic Reference Panel.

A major challenge in modern biology is to understand how naturally occurring variation in DNA sequences affects complex organismal traits through networks of intermediate molecular phenotypes. This question is best addressed in a genetic mapping population in which all molecular polymorphisms are known and for which molecular endophenotypes and complex traits are assessed on the same genotypes. Here, we performed deep RNA sequencing of 200 Drosophila Genetic Reference Panel inbred lines with complete genome sequences and for which phenotypes of many quantitative traits have been evaluated. We mapped expression quantitative trait loci for annotated genes, novel transcribed regions, transposable elements, and microbial species. We identified host variants that affect expression of transposable elements, independent of their copy number, as well as microbiome composition. We constructed sex-specific expression quantitative trait locus regulatory networks. These networks are enriched for novel transcribed regions and target genes in heterochromatin and euchromatic regions of reduced recombination, as well as genes regulating transposable element expression. This study provides new insights regarding the role of natural genetic variation in regulating gene expression and generates testable hypotheses for future functional analyses.

Context-dependent genetic architecture of Drosophila life span.

Understanding the genetic basis of variation in life span is a major challenge that is difficult to address in human populations. Evolutionary theory predicts that alleles affecting natural variation in life span will have properties that enable them to persist in populations at intermediate frequencies, such as late-life-specific deleterious effects, antagonistic pleiotropic effects on early and late-age fitness components, and/or sex- and environment-specific or antagonistic effects. Here, we quantified variation in life span in males and females reared in 3 thermal environments for the sequenced, inbred lines of the Drosophila melanogaster Genetic Reference Panel (DGRP) and an advanced intercross outbred population derived from a subset of DGRP lines. Quantitative genetic analyses of life span and the micro-environmental variance of life span in the DGRP revealed significant genetic variance for both traits within each sex and environment, as well as significant genotype-by-sex interaction (GSI) and genotype-by-environment interaction (GEI). Genome-wide association (GWA) mapping in both populations implicates over 2,000 candidate genes with sex- and environment-specific or antagonistic pleiotropic allelic effects. Over 1,000 of these genes are associated with variation in life span in other D. melanogaster populations. We functionally assessed the effects of 15 candidate genes using RNA interference (RNAi): all affected life span and/or micro-environmental variance of life span in at least one sex and environment and exhibited sex-and environment-specific effects. Our results implicate novel candidate genes affecting life span and suggest that variation for life span may be maintained by variable allelic effects in heterogeneous environments.

Genetic architecture of natural variation in visual senescence in Drosophila.

Senescence, i.e., functional decline with age, is a major determinant of health span in a rapidly aging population, but the genetic basis of interindividual variation in senescence remains largely unknown. Visual decline and age-related eye disorders are common manifestations of senescence, but disentangling age-dependent visual decline in human populations is challenging due to inability to control genetic background and variation in histories of environmental exposures. We assessed the genetic basis of natural variation in visual senescence by measuring age-dependent decline in phototaxis using Drosophila melanogaster as a genetic model system. We quantified phototaxis at 1, 2, and 4 wk of age in the sequenced, inbred lines of the Drosophila melanogaster Genetic Reference Panel (DGRP) and found an average decline in phototaxis with age. We observed significant genetic variation for phototaxis at each age and significant genetic variation in senescence of phototaxis that is only partly correlated with phototaxis. Genome-wide association analyses in the DGRP and a DGRP-derived outbred, advanced intercross population identified candidate genes and genetic networks associated with eye and nervous system development and function, including seven genes with human orthologs previously associated with eye diseases. Ninety percent of candidate genes were functionally validated with targeted RNAi-mediated suppression of gene expression. Absence of candidate genes previously implicated with longevity indicates physiological systems may undergo senescence independent of organismal life span. Furthermore, we show that genes that shape early developmental processes also contribute to senescence, demonstrating that senescence is part of a genetic continuum that acts throughout the life span.

The Drosophila melanogaster Genetic Reference Panel.

A major challenge of biology is understanding the relationship between molecular genetic variation and variation in quantitative traits, including fitness. This relationship determines our ability to predict phenotypes from genotypes and to understand how evolutionary forces shape variation within and between species. Previous efforts to dissect the genotype-phenotype map were based on incomplete genotypic information. Here, we describe the Drosophila melanogaster Genetic Reference Panel (DGRP), a community resource for analysis of population genomics and quantitative traits. The DGRP consists of fully sequenced inbred lines derived from a natural population. Population genomic analyses reveal reduced polymorphism in centromeric autosomal regions and the X chromosome, evidence for positive and negative selection, and rapid evolution of the X chromosome. Many variants in novel genes, most at low frequency, are associated with quantitative traits and explain a large fraction of the phenotypic variance. The DGRP facilitates genotype-phenotype mapping using the power of Drosophila genetics.

Genetic architecture of natural variation in Drosophila melanogaster aggressive behavior.

Aggression is an evolutionarily conserved complex behavior essential for survival and the organization of social hierarchies. With the exception of genetic variants associated with bioamine signaling, which have been implicated in aggression in many species, the genetic basis of natural variation in aggression is largely unknown. Drosophila melanogaster is a favorable model system for exploring the genetic basis of natural variation in aggression. Here, we performed genome-wide association analyses using the inbred, sequenced lines of the Drosophila melanogaster Genetic Reference Panel (DGRP) and replicate advanced intercross populations derived from the most and least aggressive DGRP lines. We identified genes that have been previously implicated in aggressive behavior as well as many novel loci, including gustatory receptor 63a (Gr63a), which encodes a subunit of the receptor for CO2, and genes associated with development and function of the nervous system. Although genes from the two association analyses were largely nonoverlapping, they mapped onto a genetic interaction network inferred from an analysis of pairwise epistasis in the DGRP. We used mutations and RNAi knock-down alleles to functionally validate 79% of the candidate genes and 75% of the candidate epistatic interactions tested. Epistasis for aggressive behavior causes cryptic genetic variation in the DGRP that is revealed by changing allele frequencies in the outbred populations derived from extreme DGRP lines. This phenomenon may pertain to other fitness traits and species, with implications for evolution, applied breeding, and human genetics.

Genetic basis of transcriptome diversity in Drosophila melanogaster.

Understanding how DNA sequence variation is translated into variation for complex phenotypes has remained elusive but is essential for predicting adaptive evolution, for selecting agriculturally important animals and crops, and for personalized medicine. Gene expression may provide a link between variation in DNA sequence and organismal phenotypes, and its abundance can be measured efficiently and accurately. Here we quantified genome-wide variation in gene expression in the sequenced inbred lines of the Drosophila melanogaster Genetic Reference Panel (DGRP), increasing the annotated Drosophila transcriptome by 11%, including thousands of novel transcribed regions (NTRs). We found that 42% of the Drosophila transcriptome is genetically variable in males and females, including the NTRs, and is organized into modules of genetically correlated transcripts. We found that NTRs often were negatively correlated with the expression of protein-coding genes, which we exploited to annotate NTRs functionally. We identified regulatory variants for the mean and variance of gene expression, which have largely independent genetic control. Expression quantitative trait loci (eQTLs) for the mean, but not for the variance, of gene expression were concentrated near genes. Notably, the variance eQTLs often interacted epistatically with local variants in these genes to regulate gene expression. This comprehensive characterization of population-scale diversity of transcriptomes and its genetic basis in the DGRP is critically important for a systems understanding of quantitative trait variation.

Natural variation in genome architecture among 205 Drosophila melanogaster Genetic Reference Panel lines.

The Drosophila melanogaster Genetic Reference Panel (DGRP) is a community resource of 205 sequenced inbred lines, derived to improve our understanding of the effects of naturally occurring genetic variation on molecular and organismal phenotypes. We used an integrated genotyping strategy to identify 4,853,802 single nucleotide polymorphisms (SNPs) and 1,296,080 non-SNP variants. Our molecular population genomic analyses show higher deletion than insertion mutation rates and stronger purifying selection on deletions. Weaker selection on insertions than deletions is consistent with our observed distribution of genome size determined by flow cytometry, which is skewed toward larger genomes. Insertion/deletion and single nucleotide polymorphisms are positively correlated with each other and with local recombination, suggesting that their nonrandom distributions are due to hitchhiking and background selection. Our cytogenetic analysis identified 16 polymorphic inversions in the DGRP. Common inverted and standard karyotypes are genetically divergent and account for most of the variation in relatedness among the DGRP lines. Intriguingly, variation in genome size and many quantitative traits are significantly associated with inversions. Approximately 50% of the DGRP lines are infected with Wolbachia, and four lines have germline insertions of Wolbachia sequences, but effects of Wolbachia infection on quantitative traits are rarely significant. The DGRP complements ongoing efforts to functionally annotate the Drosophila genome. Indeed, 15% of all D. melanogaster genes segregate for potentially damaged proteins in the DGRP, and genome-wide analyses of quantitative traits identify novel candidate genes. The DGRP lines, sequence data, genotypes, quality scores, phenotypes, and analysis and visualization tools are publicly available.

Epistasis dominates the genetic architecture of Drosophila quantitative traits.

Epistasis-nonlinear genetic interactions between polymorphic loci-is the genetic basis of canalization and speciation, and epistatic interactions can be used to infer genetic networks affecting quantitative traits. However, the role that epistasis plays in the genetic architecture of quantitative traits is controversial. Here, we compared the genetic architecture of three Drosophila life history traits in the sequenced inbred lines of the Drosophila melanogaster Genetic Reference Panel (DGRP) and a large outbred, advanced intercross population derived from 40 DGRP lines (Flyland). We assessed allele frequency changes between pools of individuals at the extremes of the distribution for each trait in the Flyland population by deep DNA sequencing. The genetic architecture of all traits was highly polygenic in both analyses. Surprisingly, none of the SNPs associated with the traits in Flyland replicated in the DGRP and vice versa. However, the majority of these SNPs participated in at least one epistatic interaction in the DGRP. Despite apparent additive effects at largely distinct loci in the two populations, the epistatic interactions perturbed common, biologically plausible, and highly connected genetic networks. Our analysis underscores the importance of epistasis as a principal factor that determines variation for quantitative traits and provides a means to uncover genetic networks affecting these traits. Knowledge of epistatic networks will contribute to our understanding of the genetic basis of evolutionarily and clinically important traits and enhance predictive ability at an individualized level in medicine and agriculture.

Complex genetic architecture of Drosophila aggressive behavior.

Epistasis and pleiotropy feature prominently in the genetic architecture of quantitative traits but are difficult to assess in outbred populations. We performed a diallel cross among coisogenic Drosophila P-element mutations associated with hyperaggressive behavior and showed extensive epistatic and pleiotropic effects on aggression, brain morphology, and genome-wide transcript abundance in head tissues. Epistatic interactions were often of greater magnitude than homozygous effects, and the topology of epistatic networks varied among these phenotypes. The transcriptional signatures of homozygous and double heterozygous genotypes derived from the six mutations imply a large mutational target for aggressive behavior and point to evolutionarily conserved genetic mechanisms and neural signaling pathways affecting this universal fitness trait.

Quantitative and molecular genetic analyses of mutations increasing Drosophila life span.

Understanding the genetic and environmental factors that affect variation in life span and senescence is of major interest for human health and evolutionary biology. Multiple mechanisms affect longevity, many of which are conserved across species, but the genetic networks underlying each mechanism and cross-talk between networks are unknown. We report the results of a screen for mutations affecting Drosophila life span. One third of the 1,332 homozygous P-element insertion lines assessed had quantitative effects on life span; mutations reducing life span were twice as common as mutations increasing life span. We confirmed 58 mutations with increased longevity, only one of which is in a gene previously associated with life span. The effects of the mutations increasing life span were highly sex-specific, with a trend towards opposite effects in males and females. Mutations in the same gene were associated with both increased and decreased life span, depending on the location and orientation of the P-element insertion, and genetic background. We observed substantial--and sex-specific--epistasis among a sample of ten mutations with increased life span. All mutations increasing life span had at least one deleterious pleiotropic effect on stress resistance or general health, with different patterns of pleiotropy for males and females. Whole-genome transcript profiles of seven of the mutant lines and the wild type revealed 4,488 differentially expressed transcripts, 553 of which were common to four or more of the mutant lines, which include genes previously associated with life span and novel genes implicated by this study. Therefore longevity has a large mutational target size; genes affecting life span have variable allelic effects; alleles affecting life span exhibit antagonistic pleiotropy and form epistatic networks; and sex-specific mutational effects are ubiquitous. Comparison of transcript profiles of long-lived mutations and the control line reveals a transcriptional signature of increased life span.

Asymmetrical lineage introgression and recombination in populations of Aspergillus flavus: Implications for biological control.

Aspergillus flavus is an agriculturally important fungus that causes ear rot of maize and produces aflatoxins, of which B1 is the most carcinogenic naturally-produced compound. In the US, the management of aflatoxins includes the deployment of biological control agents that comprise two nonaflatoxigenic A. flavus strains, either Afla-Guard (member of lineage IB) or AF36 (lineage IC). We used genotyping-by-sequencing to examine the influence of both biocontrol agents on native populations of A. flavus in cornfields in Texas, North Carolina, Arkansas, and Indiana. This study examined up to 27,529 single-nucleotide polymorphisms (SNPs) in a total of 815 A. flavus isolates, and 353 genome-wide haplotypes sampled before biocontrol application, three months after biocontrol application, and up to three years after initial application. Here, we report that the two distinct A. flavus evolutionary lineages IB and IC differ significantly in their frequency distributions across states. We provide evidence of increased unidirectional gene flow from lineage IB into IC, inferred to be due to the applied Afla-Guard biocontrol strain. Genetic exchange and recombination of biocontrol strains with native strains was detected in as little as three months after biocontrol application and up to one and three years later. There was limited inter-lineage migration in the untreated fields. These findings suggest that biocontrol products that include strains from lineage IB offer the greatest potential for sustained reductions in aflatoxin levels over several years. This knowledge has important implications for developing new biocontrol strategies.

Genotype by environment interaction for gene expression in Drosophila melanogaster.

The genetics of phenotypic responses to changing environments remains elusive. Using whole-genome quantitative gene expression as a model, here we study how the genetic architecture of regulatory variation in gene expression changed in a population of fully sequenced inbred Drosophila melanogaster strains when flies developed in different environments (25 °C and 18 °C). We find a substantial fraction of the transcriptome exhibited genotype by environment interaction, implicating environmentally plastic genetic architecture of gene expression. Genetic variance in expression increases at 18 °C relative to 25 °C for most genes that have a change in genetic variance. Although the majority of expression quantitative trait loci (eQTLs) for the gene expression traits in the two environments are shared and have similar effects, analysis of the environment-specific eQTLs reveals enrichment of binding sites for two transcription factors. Finally, although genotype by environment interaction in gene expression could potentially disrupt genetic networks, the co-expression networks are highly conserved across environments. Genes with higher network connectivity are under stronger stabilizing selection, suggesting that stabilizing selection on expression plays an important role in promoting network robustness.

Phenotypic variation and natural selection at catsup, a pleiotropic quantitative trait gene in Drosophila.

Quantitative traits are shaped by networks of pleiotropic genes . To understand the mechanisms that maintain genetic variation for quantitative traits in natural populations and to predict responses to artificial and natural selection, we must evaluate pleiotropic effects of underlying quantitative trait genes and define functional allelic variation at the level of quantitative trait nucleotides (QTNs). Catecholamines up (Catsup), which encodes a negative regulator of tyrosine hydroxylase , the rate-limiting step in the synthesis of the neurotransmitter dopamine, is a pleiotropic quantitative trait gene in Drosophila melanogaster. We used association mapping to determine whether the same or different QTNs at Catsup are associated with naturally occurring variation in multiple quantitative traits. We sequenced 169 Catsup alleles from a single population and detected 33 polymorphisms with little linkage disequilibrium (LD). Different molecular polymorphisms in Catsup are independently associated with variation in longevity, locomotor behavior, and sensory bristle number. Most of these polymorphisms are potentially functional variants in protein coding regions, have large effects, and are not common. Thus, Catsup is a pleiotropic quantitative trait gene, but individual QTNs do not have pleiotropic effects. Molecular population genetic analyses of Catsup sequences are consistent with balancing selection maintaining multiple functional polymorphisms.

Phenotypic plasticity and genotype by environment interaction for olfactory behavior in Drosophila melanogaster.

Genotype by environment interactions (GEI) play a major part in shaping the genetic architecture of quantitative traits and are confounding factors in genetic studies, for example, in attempts to associate genetic variation with disease susceptibility. It is generally not known what proportion of phenotypic variation is due to GEI and how many and which genes contribute to GEI. Behaviors are complex traits that mediate interactions with the environment and, thus, are ideally suited for studies of GEI. Olfactory behavior in Drosophila melanogaster presents an opportunity to systematically dissect GEI, since large numbers of genetically identical individuals can be reared under defined environmental conditions and the olfactory system of Drosophila and its behavioral response to odorants have been well characterized. We assessed variation in olfactory behavior in a population of 41 wild-derived inbred lines and asked to what extent different larval-rearing environments would influence adult olfactory behavior and whether GEI is a minor or major contributing source of phenotypic variation. We found that approximately 50% of phenotypic variation in adult olfactory behavior is attributable to GEI. In contrast, transcriptional analysis revealed that only 20 genes show GEI at the level of gene expression [false discovery rate (FDR) < 0.05], some of which are associated with physiological responses to environmental chemicals. Quantitative complementation tests with piggyBac-tagged mutants for 2 of these genes (CG9664 and Transferrin 1) demonstrate that genes that show transcriptional GEI are candidate genes for olfactory behavior and that GEI at the level of gene expression is correlated with GEI at the level of phenotype.

Lisinopril Preserves Physical Resilience and Extends Life Span in a Genotype-Specific Manner in Drosophila melanogaster.

Physical resiliency declines with age and comorbid conditions. In humans, angiotensin-converting enzyme (ACE) has been associated with attenuation of the decline in physical performance with age. ACE-inhibitor compounds, commonly prescribed for hypertension, often have beneficial effects on physical performance however the generality of these effects are unclear. Here, we tested the effects of the ACE-inhibitor Lisinopril on life span, and age-specific speed, endurance, and strength using three genotypes of the Drosophila melanogaster Genetic Reference Panel. We show that age-related decline in physical performance and survivorship varies with genetic background. Lisinopril treatment increased mean life span in all Drosophila Genetic Reference Panel lines, but its effects on life span, speed, endurance, and strength depended on genotype. We show that genotypes with increased physical performance on Lisinopril treatment experienced reduced age-related protein aggregation in muscle. Knockdown of skeletal muscle-specific Ance, the Drosophila ortholog of ACE, abolished the effects of Lisinopril on life span, implying a role for skeletal muscle Ance in survivorship. Using transcriptome profiling, we identified genes involved in stress response that showed expression changes associated with genotype and age-dependent responsiveness to Lisinopril. Our results demonstrate that Ance is involved in physical decline and demonstrate genetic variation in phenotypic responses to an ACE inhibitor.

Gene-environment contributions to the development of infant vagal reactivity: the interaction of dopamine and maternal sensitivity.

This study investigated dopamine receptor genes (DRD2 and DRD4) and maternal sensitivity as predictors of infant respiratory sinus arrhythmia (RSA) and RSA reactivity, purported indices of vagal tone and vagal regulation, in a challenge task at 3, 6, and 12 months in 173 infant-mother dyads. Hierarchical linear modeling (HLM) revealed that at 3 and 6 months, RSA withdrawal in response to maternal separation was greater (suggesting expected physiological regulation) in infants without the DRD2 risk allele than those with the risk allele. At 12 months, infants with the risk allele who were also exposed to maternal sensitivity showed levels of RSA withdrawal comparable to infants who were not at genetic risk. Findings demonstrate the importance of developmental analysis of gene-environment interaction.

Regulation of Drosophila Lifespan by bellwether Promoter Alleles.

Longevity varies among individuals, but how natural genetic variation contributes to variation in lifespan is poorly understood. Drosophila melanogaster presents an advantageous model system to explore the genetic underpinnings of longevity, since its generation time is brief and both the genetic background and rearing environment can be precisely controlled. The bellwether (blw) gene encodes the α subunit of mitochondrial ATP synthase. Since metabolic rate may influence lifespan, we investigated whether alternative haplotypes in the blw promoter affect lifespan when expressed in a co-isogenic background. We amplified 521 bp upstream promoter sequences containing alternative haplotypes and assessed promoter activity both in vitro and in vivo using a luciferase reporter system. The AG haplotype showed significantly greater expression of luciferase than the GT haplotype. We then overexpressed a blw cDNA construct driven by either the AG or GT haplotype promoter in transgenic flies and showed that the AG haplotype also results in greater blw cDNA expression and a significant decrease in lifespan relative to the GT promoter haplotype, in male flies only. Thus, our results show that naturally occurring regulatory variants of blw affect lifespan in a sex-specific manner.

Quantitative trait loci affecting the difference in pigmentation between Drosophila yakuba and D. santomea.

Using quantitative trait locus (QTL) mapping, we studied the genetic basis of the difference in pigmentation between two sister species of Drosophila: Drosophila yakuba, which, like other members of the D. melanogaster subgroup, shows heavy black pigmentation on the abdomen of males and females, and D. santomea, an endemic to the African island of São Tomé, which has virtually no pigmentation. Here we mapped four QTL with large effects on this interspecific difference in pigmentation: two on the X chromosome and one each on the second and third chromosomes. The same four QTL were detected in male hybrids in the backcrosses to both D. santomea and D. yakuba and in the female D. yakuba backcross hybrids. All four QTL exhibited strong epistatic interactions in male backcross hybrids, but only one pair of QTL interacted in females from the backcross to D. yabuka. All QTL from each species affected pigmentation in the same direction, consistent with adaptive evolution driven by directional natural selection. The regions delimited by the QTL included many positional candidate loci in the pigmentation pathway, including genes affecting catecholamine biosynthesis, melanization of the cuticle, and many additional pleiotropic effects.

Spontaneous mutations and the origin and maintenance of quantitative genetic variation.

Mutation and natural selection shape the genetic variation in natural populations. Here, we directly estimated the spontaneous mutation rate by sequencing new Drosophila mutation accumulation lines maintained with minimal natural selection. We inferred strong stabilizing natural selection on quantitative traits because genetic variation among wild-derived inbred lines was much lower than predicted from a neutral model and the mutational effects were much larger than allelic effects of standing polymorphisms. Stabilizing selection could act directly on the traits, or indirectly from pleiotropic effects on fitness. However, our data are not consistent with simple models of mutation-stabilizing selection balance; therefore, further empirical work is needed to assess the balance of evolutionary forces responsible for quantitative genetic variation.