Pleiotropy, epistasis and the genetic architecture of quantitative traits.

Pleiotropy (whereby one genetic polymorphism affects multiple traits) and epistasis (whereby non-linear interactions between genetic polymorphisms affect the same trait) are fundamental aspects of the genetic architecture of quantitative traits. Recent advances in the ability to characterize the effects of polymorphic variants on molecular and organismal phenotypes in human and model organism populations have revealed the prevalence of pleiotropy and unexpected shared molecular genetic bases among quantitative traits, including diseases. By contrast, epistasis is common between polymorphic loci associated with quantitative traits in model organisms, such that alleles at one locus have different effects in different genetic backgrounds, but is rarely observed for human quantitative traits and common diseases. Here, we review the concepts and recent inferences about pleiotropy and epistasis, and discuss factors that contribute to similarities and differences between the genetic architecture of quantitative traits in model organisms and humans.

The genetic architecture of behavioral canalization.

Behaviors are components of fitness and contribute to adaptive evolution. Behaviors represent the interactions of an organism with its environment, yet innate behaviors display robustness in the face of environmental change, which we refer to as 'behavioral canalization'. We hypothesize that positive selection of hub genes of genetic networks stabilizes the genetic architecture for innate behaviors by reducing variation in the expression of interconnected network genes. Robustness of these stabilized networks would be protected from deleterious mutations by purifying selection or suppressing epistasis. We propose that, together with newly emerging favorable mutations, epistatically suppressed mutations can generate a reservoir of cryptic genetic variation that could give rise to decanalization when genetic backgrounds or environmental conditions change to allow behavioral adaptation.

Gregor Mendel's legacy in quantitative genetics.

Gregor Mendel's discovery of the laws of segregation and independent assortment and his inference of the existence of non-mendelian interactions between loci remain at the heart of today's explorations of the genetic architecture of quantitative traits.

The Drosophila brain on cocaine at single-cell resolution.

Whereas the neurological effects of cocaine have been well documented, effects of acute cocaine consumption on genome-wide gene expression across the brain remain largely unexplored. This question cannot be readily addressed in humans but can be approached using the Drosophila melanogaster model, where gene expression in the entire brain can be surveyed at once. Flies exposed to cocaine show impaired locomotor activity, including climbing behavior and startle response (a measure of sensorimotor integration), and increased incidence of seizures and compulsive grooming. To identify specific cell populations that respond to acute cocaine exposure, we analyzed single-cell transcriptional responses in duplicate samples of flies that consumed fixed amounts of sucrose or sucrose supplemented with cocaine, in both sexes. Unsupervised clustering of the transcriptional profiles of a total of 86,224 cells yielded 36 distinct clusters. Annotation of clusters based on gene markers revealed that all major cell types (neuronal and glial) as well as neurotransmitter types from most brain regions were represented. The brain transcriptional responses to cocaine showed profound sexual dimorphism and were considerably more pronounced in males than females. Differential expression analysis within individual clusters indicated cluster-specific responses to cocaine. Clusters corresponding to Kenyon cells of the mushroom bodies and glia showed especially large transcriptional responses following cocaine exposure. Cluster specific coexpression networks and global interaction networks revealed a diverse array of cellular processes affected by acute cocaine exposure. These results provide an atlas of sexually dimorphic cocaine-modulated gene expression in a model brain.

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.

Evolution of Epistatic Networks and the Genetic Basis of Innate Behaviors.

Instinctive behaviors are genetically programmed behaviors that occur independent of experience. How genetic programs that give rise to the manifestation of such behaviors evolve remains an unresolved question. I propose that evolution of species-specific innate behaviors is accomplished through progressive modifications of pre-existing genetic networks composed of allelic variants. I hypothesize that changes in frequencies of one or more constituent allelic variants within the network leads to changes in gene network connectivity and the emergence of a reorganized network that can support the emergence of a novel behavioral phenotype and becomes stabilized when key allelic variants are driven to fixation.

Systems genetics of the Drosophila metabolome.

How effects of DNA sequence variants are transmitted through intermediate endophenotypes to modulate organismal traits remains a central question in quantitative genetics. This problem can be addressed through a systems approach in a population in which genetic polymorphisms, gene expression traits, metabolites, and complex phenotypes can be evaluated on the same genotypes. Here, we focused on the metabolome, which represents the most proximal link between genetic variation and organismal phenotype, and quantified metabolite levels in 40 lines of the Drosophila melanogaster Genetic Reference Panel. We identified sex-specific modules of genetically correlated metabolites and constructed networks that integrate DNA sequence variation and variation in gene expression with variation in metabolites and organismal traits, including starvation stress resistance and male aggression. Finally, we asked to what extent SNPs and metabolites can predict trait phenotypes and generated trait- and sex-specific prediction models that provide novel insights about the metabolomic underpinnings of complex phenotypes.

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.

Natural genetic variation in a dopamine receptor is associated with variation in female fertility in Drosophila melanogaster.

Fertility is a major component of fitness but its genetic architecture remains poorly understood. Using a full diallel cross of 50 Drosophila Genetic Reference Panel inbred lines with whole genome sequences, we found substantial genetic variation in fertility largely attributable to females. We mapped genes associated with variation in female fertility by genome-wide association analysis of common variants in the fly genome. Validation of candidate genes by RNAi knockdown confirmed the role of the dopamine 2-like receptor (Dop2R) in promoting egg laying. We replicated the Dop2R effect in an independently collected productivity dataset and showed that the effect of the Dop2R variant was mediated in part by regulatory gene expression variation. This study demonstrates the strong potential of genome-wide association analysis in this diverse panel of inbred strains and subsequent functional analyses for understanding the genetic architecture of fitness traits.

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.

Modulation of the Drosophila transcriptome by developmental exposure to alcohol.

BACKGROUND: Prenatal exposure to ethanol can cause fetal alcohol spectrum disorder (FASD), a prevalent, preventable pediatric disorder. Identifying genetic risk alleles for FASD is challenging since time, dose, and frequency of exposure are often unknown, and manifestations of FASD are diverse and evident long after exposure. Drosophila melanogaster is an excellent model to study the genetic basis of the effects of developmental alcohol exposure since many individuals of the same genotype can be reared under controlled environmental conditions. RESULTS: We used 96 sequenced, wild-derived inbred lines from the Drosophila melanogaster Genetic Reference Panel (DGRP) to profile genome-wide transcript abundances in young adult flies that developed on ethanol-supplemented medium or standard culture medium. We found substantial genetic variation in gene expression in response to ethanol with extensive sexual dimorphism. We constructed sex-specific genetic networks associated with alcohol-dependent modulation of gene expression that include protein-coding genes, Novel Transcribed Regions (NTRs, postulated to encode long non-coding RNAs) and female-specific coordinated regulation of snoRNAs that regulate pseudouridylation of ribosomal RNA. We reared DGRP lines which showed extreme upregulation or downregulation of snoRNA expression during developmental alcohol exposure on standard or ethanol supplemented medium and demonstrated that developmental exposure to ethanol has genotype-specific effects on adult locomotor activity and sleep. CONCLUSIONS: There is significant and sex-specific natural genetic variation in the transcriptional response to developmental exposure to ethanol in Drosophila that comprises networks of genes affecting nervous system development and ethanol metabolism as well as networks of regulatory non-coding RNAs.

Pleiotropic fitness effects of the lncRNA Uhg4 in Drosophila melanogaster.

BACKGROUND: Long noncoding RNAs (lncRNAs) are a diverse class of RNAs that are critical for gene regulation, DNA repair, and splicing, and have been implicated in development, stress response, and cancer. However, the functions of many lncRNAs remain unknown. In Drosophila melanogaster, U snoRNA host gene 4 (Uhg4) encodes an antisense long noncoding RNA that is host to seven small nucleolar RNAs (snoRNAs). Uhg4 is expressed ubiquitously during development and in all adult tissues, with maximal expression in ovaries; however, it has no annotated function(s). RESULTS: We used CRISPR-Cas9 germline gene editing to generate multiple deletions spanning the promoter region and first exon of Uhg4. Females showed arrested egg development and both males and females were sterile. In addition, Uhg4 deletion mutants showed delayed development and decreased viability, and changes in sleep and responses to stress. Whole-genome RNA sequencing of Uhg4 deletion flies and their controls identified co-regulated genes and genetic interaction networks associated with Uhg4. Gene ontology analyses highlighted a broad spectrum of biological processes, including regulation of transcription and translation, morphogenesis, and stress response. CONCLUSION: Uhg4 is a lncRNA essential for reproduction with pleiotropic effects on multiple fitness traits.

Functional Diversification, Redundancy, and Epistasis among Paralogs of the Drosophila melanogaster Obp50a-d Gene Cluster.

Large multigene families, such as the insect odorant-binding proteins (OBPs), are thought to arise through functional diversification after repeated gene duplications. Whereas many OBPs function in chemoreception, members of this family are also expressed in tissues outside chemosensory organs. Paralogs of the Obp50 gene cluster are expressed in metabolic and male reproductive tissues, but their functions and interrelationships remain unknown. Here, we report the genetic dissection of four members of the Obp50 cluster, which are in close physical proximity without intervening genes. We used CRISPR technology to excise the entire cluster while introducing a PhiC31 reintegration site to reinsert constructs in which different combinations of the constituent Obp genes were either intact or rendered inactive. We performed whole transcriptome sequencing and assessed sexually dimorphic changes in transcript abundances (transcriptional niches) associated with each gene-edited genotype. Using this approach, we were able to estimate redundancy, additivity, diversification, and epistasis among Obp50 paralogs. We analyzed the effects of gene editing of this cluster on organismal phenotypes and found a significant skewing of sex ratios attributable to Obp50a, and sex-specific effects on starvation stress resistance attributable to Obp50d. Thus, there is functional diversification within the Obp50 cluster with Obp50a contributing to development and Obp50d to stress resistance. The deletion-reinsertion approach we applied to the Obp50 cluster provides a general paradigm for the genetic dissection of paralogs of multigene families.

Genetics of cocaine and methamphetamine consumption and preference in Drosophila melanogaster.

Illicit use of psychostimulants, such as cocaine and methamphetamine, constitutes a significant public health problem. Whereas neural mechanisms that mediate the effects of these drugs are well-characterized, genetic factors that account for individual variation in susceptibility to substance abuse and addiction remain largely unknown. Drosophila melanogaster can serve as a translational model for studies on substance abuse, since flies have a dopamine transporter that can bind cocaine and methamphetamine, and exposure to these compounds elicits effects similar to those observed in people, suggesting conserved evolutionary mechanisms underlying drug responses. Here, we used the D. melanogaster Genetic Reference Panel to investigate the genetic basis for variation in psychostimulant drug consumption, to determine whether similar or distinct genetic networks underlie variation in consumption of cocaine and methamphetamine, and to assess the extent of sexual dimorphism and effect of genetic context on variation in voluntary drug consumption. Quantification of natural genetic variation in voluntary consumption, preference, and change in consumption and preference over time for cocaine and methamphetamine uncovered significant genetic variation for all traits, including sex-, exposure- and drug-specific genetic variation. Genome wide association analyses identified both shared and drug-specific candidate genes, which could be integrated in genetic interaction networks. We assessed the effects of ubiquitous RNA interference (RNAi) on consumption behaviors for 34 candidate genes: all affected at least one behavior. Finally, we utilized RNAi knockdown in the nervous system to implicate dopaminergic neurons and the mushroom bodies as part of the neural circuitry underlying experience-dependent development of drug preference.

Genetics of aggression.

Aggression mediates competition for food, mating partners, and habitats and, among social animals, establishes stable dominance hierarchies. In humans, abnormal aggression is a hallmark of neuropsychiatric disorders and can be elicited by environmental factors acting on an underlying genetic susceptibility. Identifying the genetic architecture that predisposes to aggressive behavior in people is challenging because of difficulties in quantifying the phenotype, genetic heterogeneity, and uncontrolled environmental conditions. Studies on mice have identified single-gene mutations that result in hyperaggression, contingent on genetic background. These studies can be complemented by systems genetics approaches in Drosophila melanogaster, in which mutational analyses together with genome-wide transcript analyses, artificial selection studies, and genome-wide analysis of epistasis have revealed that a large segment of the genome contributes to the manifestation of aggressive behavior with widespread epistatic interactions. Comparative genomic analyses based on the principle of evolutionary conservation are needed to enable a complete dissection of the neurogenetic underpinnings of this universal fitness trait.

A Drosophila model for toxicogenomics: Genetic variation in susceptibility to heavy metal exposure.

The genetic factors that give rise to variation in susceptibility to environmental toxins remain largely unexplored. Studies on genetic variation in susceptibility to environmental toxins are challenging in human populations, due to the variety of clinical symptoms and difficulty in determining which symptoms causally result from toxic exposure; uncontrolled environments, often with exposure to multiple toxicants; and difficulty in relating phenotypic effect size to toxic dose, especially when symptoms become manifest with a substantial time lag. Drosophila melanogaster is a powerful model that enables genome-wide studies for the identification of allelic variants that contribute to variation in susceptibility to environmental toxins, since the genetic background, environmental rearing conditions and toxic exposure can be precisely controlled. Here, we used extreme QTL mapping in an outbred population derived from the D. melanogaster Genetic Reference Panel to identify alleles associated with resistance to lead and/or cadmium, two ubiquitous environmental toxins that present serious health risks. We identified single nucleotide polymorphisms (SNPs) associated with variation in resistance to both heavy metals as well as SNPs associated with resistance specific to each of them. The effects of these SNPs were largely sex-specific. We applied mutational and RNAi analyses to 33 candidate genes and functionally validated 28 of them. We constructed networks of candidate genes as blueprints for orthologous networks of human genes. The latter not only provided functional contexts for known human targets of heavy metal toxicity, but also implicated novel candidate susceptibility genes. These studies validate Drosophila as a translational toxicogenomics gene discovery system.

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 road less traveled: from genotype to phenotype in flies and humans.

Understanding how genomic variation gives rise to phenotypic variation is essential for elucidating mechanisms of adaptive evolution, plant and animal breeding, and precision medicine. However, identifying causal links between DNA sequence variants and variation in phenotypes is challenging in human populations, due to large blocks of linkage disequilibrium in the genome and heterogeneous developmental histories, lifestyles, and social and physical environments. Drosophila melanogaster presents a powerful genetic model, since linkage disequilibrium decays rapidly, facilitating assignment of causality to polymorphisms associated with phenotypic variation, and large numbers of individuals can be reared under defined environmental conditions, economically, and without regulatory restrictions. The D. melanogaster Genetic Reference Panel (DGRP), a population of 205 sequenced, inbred wild-derived flies, has enabled genome-wide association studies of morphological, physiological, behavioral, and life history traits, and demonstrated that genetic architectures of complex traits are highly polygenic, sexually dimorphic, and context dependent with extensive sex-, environment-, and genetic background (epistatic) effects. These features together with a modular organization of the transcriptome illustrate a dynamic integrative genetic architecture for complex traits. The complexity of the genetic architectures for complex traits in Drosophila provides important caveats for the interpretation of genetic studies in human populations. Aspects of the genetic underpinnings of complex traits can be represented as simplified gene networks on which human orthologues can be superimposed to provide blueprints for subsequent studies on analogous traits in human populations. Fundamental principles of the genetic architectures of Drosophila complex traits are likely applicable across phyla, from the DGRP to human populations.

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.