Genetic architecture of natural variation in visual senescence in Drosophila

The genetic basis of variation in Drosophila melanogaster mating behavior

Mating behavior is an essential fitness trait. We used the inbred, sequenced lines of the Drosophila Genetic Reference Panel (DGRP) to gain insights into the evolution of mating success and to evaluate the overlap in genetic architecture of mating behavior between the sexes. We found significant genetic variation for mating success when DGRP males and females from the same line were mated together, and when DGRP males and females were mated to an unrelated strain of the opposite sex. The mating success of DGRP males and females was not correlated when they were paired with the unrelated strain, suggesting independent genetic architecture of mating success in males and females that was confirmed by genome-wide association analyses. However, the mating success between pairs of the same or different DGRP lines was predicted accurately by the respective female and male mating success with the unrelated line.

Natural Variation in Age-Related Dopamine Neuron Degeneration is Glutathione-Dependent and Linked to Life Span

Aging is the biggest risk factor for Parkinson's disease (PD), suggesting that age-related changes in the brain promote dopamine neuron vulnerability. It is unclear, however, whether aging alone is sufficient to cause significant dopamine neuron loss and if so, how this intersects with PD-related neurodegeneration. Here, through examining a large collection of naturally varying Drosophila strains, we find a strong relationship between life span and age-related dopamine neuron loss. Naturally short-lived strains exhibit a loss of dopamine neurons but not generalized neurodegeneration, while long-lived strains retain dopamine neurons across age. Metabolomic profiling reveals lower glutathione levels in short-lived strains which is associated with elevated levels of reactive oxygen species (ROS), sensitivity to oxidative stress and vulnerability to silencing the familial PD gene parkin . Strikingly, boosting neuronal glutathione levels via glutamate-cysteine ligase (GCL) overexpression is sufficient to normalize ROS levels, extend life span and block dopamine neurons loss in short-lived backgrounds, demonstrating that glutathione deficiencies are central to neurodegenerative phenotypes associated with short longevity. These findings may be relevant to human PD pathogenesis, where glutathione depletion is frequently reported in idiopathic PD patient brain. Building on this evidence, we detect reduced levels of GCL catalytic and modulatory subunits in brain from PD patients harboring the LRRK2 G2019S mutation, implicating possible glutathione deficits in familial LRRK2-linked PD. Our study across Drosophila and human PD systems suggests that glutathione plays an important role in the influence of aging on PD neurodegeneration.

The genetic basis of variation in immune defense against Lysinibacillus fusiformis infection in Drosophila melanogaster

The genetic causes of phenotypic variation often differ depending on the population examined, particularly if the populations were founded by relatively small numbers of genotypes. Similarly, the genetic causes of phenotypic variation among similar traits (resistance to different xenobiotic compounds or pathogens) may also be completely different or only partially overlapping. Differences in genetic causes for variation in the same trait among populations suggests context dependence for how selection acts on those traits. Similarities in the genetic causes of variation for different traits, on the other hand, suggests pleiotropy which would also influence how natural selection shapes variation in a trait. We characterized immune defense against a natural Drosophila pathogen, the Gram-positive bacterium Lysinibacillus fusiformis, in three different populations and found almost no overlap in the genetic architecture of variation in survival post infection. However, when comparing our results to a similar experiment with the fungal pathogen, B. bassiana, we found a convincing shared QTL peak for both pathogens. This peak contains the Bomanin cluster of Drosophila immune effectors. Loss of function mutants and RNAi knockdown experiments confirms a role of some of these genes in immune defense against both pathogens. This suggests that natural selection may act on the entire cluster of Bomanin genes (and the linked region under the QTL) or specific peptides for specific pathogens.

Perspectives on the Drosophila melanogaster Model for Advances in Toxicological Science

The use of Drosophila melanogaster for studies of toxicology has grown considerably in the last decade. The Drosophila model has long been appreciated as a versatile and powerful model for developmental biology and genetics because of its ease of handling, short life cycle, low cost of maintenance, molecular genetic accessibility, and availability of a wide range of publicly available strains and data resources. These features, together with recent unique developments in genomics and metabolomics, make the fly model especially relevant and timely for the development of new approach methodologies and movements toward precision toxicology. Here, we offer a perspective on how flies can be leveraged to identify risk factors relevant to environmental exposures and human health. First, we review and discuss fundamental toxicologic principles for experimental design with Drosophila. Next, we describe quantitative and systems genetics approaches to resolve the genetic architecture and candidate pathways controlling susceptibility to toxicants. Finally, we summarize the current state and future promise of the emerging field of Drosophila metabolomics for elaborating toxic mechanisms. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC.

The impact of geography and climate on the population structure and local adaptation in a wild bee

Deciphering processes that contribute to genetic differentiation and divergent selection of natural populations is useful for evaluating the adaptive potential and resilience of organisms faced with various anthropogenic stressors. Insect pollinator species, including wild bees, provide critical ecosystem services but are highly susceptible to biodiversity declines. Here, we use population genomics to infer the genetic structure and test for evidence of local adaptation in an economically important native pollinator, the small carpenter bee (Ceratina calcarata). Using genome-wide SNP data (n = 8302), collected from specimens across the species' entire distribution, we evaluated population differentiation and genetic diversity and identified putative signatures of selection in the context of geographic and environmental variation. Results of the analyses of principal component and Bayesian clustering were concordant with the presence of two to three genetic clusters, associated with landscape features and inferred phylogeography of the species. All populations examined in our study demonstrated a heterozygote deficit, along with significant levels of inbreeding. We identified 250 robust outlier SNPs, corresponding to 85 annotated genes with known functional relevance to thermoregulation, photoperiod, and responses to various abiotic and biotic stressors. Taken together, these data provide evidence for local adaptation in a wild bee and highlight genetic responses of native pollinators to landscape and climate features.

Genetic architecture of natural variation of cardiac performance from flies to humans

Deciphering the genetic architecture of human cardiac disorders is of fundamental importance but their underlying complexity is a major hurdle. We investigated the natural variation of cardiac performance in the sequenced inbred lines of the Drosophila Genetic Reference Panel (DGRP). Genome-wide associations studies (GWAS) identified genetic networks associated with natural variation of cardiac traits which were used to gain insights as to the molecular and cellular processes affected. Non-coding variants that we identified were used to map potential regulatory non-coding regions, which in turn were employed to predict transcription factors (TFs) binding sites. Cognate TFs, many of which themselves bear polymorphisms associated with variations of cardiac performance, were also validated by heart-specific knockdown. Additionally, we showed that the natural variations associated with variability in cardiac performance affect a set of genes overlapping those associated with average traits but through different variants in the same genes. Furthermore, we showed that phenotypic variability was also associated with natural variation of gene regulatory networks. More importantly, we documented correlations between genes associated with cardiac phenotypes in both flies and humans, which supports a conserved genetic architecture regulating adult cardiac function from arthropods to mammals. Specifically, roles for PAX9 and EGR2 in the regulation of the cardiac rhythm were established in both models, illustrating that the characteristics of natural variations in cardiac function identified in Drosophila can accelerate discovery in humans.

Proper control of R-loop homeostasis is required for maintenance of gene expression and neuronal function during aging

Age‐related loss of cellular function and increased cell death are characteristic hallmarks of aging. While defects in gene expression and RNA metabolism have been linked with age‐associated human neuropathies, it is not clear how the changes that occur in aging neurons contribute to loss of gene expression homeostasis. R‐loops are RNA–DNA hybrids that typically form co‐transcriptionally via annealing of the nascent RNA to the template DNA strand, displacing the non‐template DNA strand. Dysregulation of R‐loop homeostasis has been associated with both transcriptional impairment and genome instability. Importantly, a growing body of evidence links R‐loop accumulation with cellular dysfunction, increased cell death, and chronic disease onset. Here, we characterized the R‐loop landscape in aging Drosophila melanogaster photoreceptor neurons and showed that bulk R‐loop levels increased with age. Further, genome‐wide mapping of R‐loops revealed that transcribed genes accumulated R‐loops over gene bodies during aging, which correlated with decreased expression of long and highly expressed genes. Importantly, while photoreceptor‐specific down‐regulation of Top3β, a DNA/RNA topoisomerase associated with R‐loop resolution, lead to decreased visual function, over‐expression of Top3β or nuclear‐localized RNase H1, which resolves R‐loops, enhanced positive light response during aging. Together, our studies highlight the functional link between dysregulation of R‐loop homeostasis, gene expression, and visual function during aging.

Aging and Light Stress Result in Overlapping and Unique Gene Expression Changes in Photoreceptors

Advanced age is one of the leading risk factors for vision loss and eye disease. Photoreceptors are the primary sensory neurons of the eye. The extended photoreceptor cell lifespan, in addition to its high metabolic needs due to phototransduction, makes it critical for these neurons to continually respond to the stresses associated with aging by mounting an appropriate gene expression response. Here, we sought to untangle the more general neuronal age-dependent transcriptional signature of photoreceptors with that induced by light stress. To do this, we aged flies or exposed them to various durations of blue light, followed by photoreceptor nuclei-specific transcriptome profiling. Using this approach, we identified genes that are both common and uniquely regulated by aging and light induced stress. Whereas both age and blue light induce expression of DNA repair genes and a neuronal-specific signature of death, both conditions result in downregulation of phototransduction. Interestingly, blue light uniquely induced genes that directly counteract the overactivation of the phototransduction signaling cascade. Lastly, unique gene expression changes in aging photoreceptors included the downregulation of genes involved in membrane potential homeostasis and mitochondrial function, as well as the upregulation of immune response genes. We propose that light stress contributes to the aging transcriptome of photoreceptors, but that there are also other environmental or intrinsic factors involved in age-associated photoreceptor gene expression signatures.

The Biology of Aging in Insects: From Drosophila to Other Insects and Back

An enormous amount of work has been done on aging in Drosophila melanogaster, a classical genetic and molecular model system, but also in numerous other insects. However, these two extensive bodies of work remain poorly integrated to date. Studies in Drosophila often explore genetic, developmental, physiological, and nutrition-related aspects of aging in the lab, while studies in other insects often explore ecological, social, and somatic aspects of aging in both lab and natural populations. Alongside exciting genomic and molecular research advances in aging in Drosophila, many new studies have also been published on aging in various other insects, including studies on aging in natural populations of diverse species. However, no broad synthesis of these largely separate bodies of work has been attempted. In this review, we endeavor to synthesize these two semi-independent literatures to facilitate collaboration and foster the exchange of ideas and research tools. While lab studies of Drosophila have illuminated many fundamental aspects of senescence, the stunning diversity of aging patterns among insects, especially in the context of their rich ecology, remains vastlyunderstudied. Coupled with field studies and novel, more easily applicable molecular methods, this represents a major opportunity for deepening our understanding of the biology of aging in insects and beyond.

Epistasis for head morphology in Drosophila melanogaster

Epistasis-gene-gene interaction-is common for mutations with large phenotypic effects in humans and model organisms. Epistasis impacts quantitative genetic models of speciation, response to natural and artificial selection, genetic mapping, and personalized medicine. However, the existence and magnitude of epistasis between alleles with small quantitative phenotypic effects are controversial and difficult to assess. Here, we use the Drosophila melanogaster Genetic Reference Panel of sequenced inbred lines to evaluate the magnitude of naturally occurring epistasis modifying the effects of mutations in jing and inv, two transcription factors that have subtle quantitative effects on head morphology as homozygotes. We find significant epistasis for both mutations and performed single marker genome-wide association analyses to map candidate modifier variants and loci affecting head morphology. A subset of these loci was significantly enriched for a known genetic interaction network, and mutations of the candidate epistatic modifier loci also affect head morphology.

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.

Natural variation in the regulation of neurodevelopmental genes modifies flight performance in Drosophila

The winged insects of the order Diptera are colloquially named for their most recognizable phenotype: flight. These insects rely on flight for a number of important life history traits, such as dispersal, foraging, and courtship. Despite the importance of flight, relatively little is known about the genetic architecture of flight performance. Accordingly, we sought to uncover the genetic modifiers of flight using a measure of flies’ reaction and response to an abrupt drop in a vertical flight column. We conducted a genome wide association study (GWAS) using 197 of the Drosophila Genetic Reference Panel (DGRP) lines, and identified a combination of additive and marginal variants, epistatic interactions, whole genes, and enrichment across interaction networks. Egfr, a highly pleiotropic developmental gene, was among the most significant additive variants identified. We functionally validated 13 of the additive candidate genes’ (Adgf-A/Adgf-A2/CG32181, bru1, CadN, flapper (CG11073), CG15236, flippy (CG9766), CREG, Dscam4, form3, fry, Lasp/CG9692, Pde6, Snoo), and introduce a novel approach to whole gene significance screens: PEGASUS_flies. Additionally, we identified ppk23, an Acid Sensing Ion Channel (ASIC) homolog, as an important hub for epistatic interactions. We propose a model that suggests genetic modifiers of wing and muscle morphology, nervous system development and function, BMP signaling, sexually dimorphic neural wiring, and gene regulation are all important for the observed differences flight performance in a natural population. Additionally, these results represent a snapshot of the genetic modifiers affecting drop-response flight performance in Drosophila, with implications for other insects.

Insect flight is a widely recognizable phenotype of many winged insects, hence the name: flies. While fruit flies, or Drosophila melanogaster, are a genetically tractable model, flight performance is a highly integrative phenotype, and therefore challenging to identify comprehensively which genetic modifiers contribute to its genetic architecture. Accordingly, we screened 197 Drosophila Genetic Reference Panel lines for their ability to react and respond to an abrupt drop. Using several computational approaches, we identified additive, marginal, and epistatic variants, as well as whole genes and altered sub-networks of gene-gene and protein-protein interaction networks that contribute to variation in flight performance. More generally, we demonstrate the benefits of employing multiple methodologies to elucidate the genetic architecture of complex traits. Many variants and genes mapped to regions of the genome that affect neurodevelopment, wing and muscle development, and regulation of gene expression. We also introduce PEGASUS_flies, a Drosophila-adapted version of the PEGASUS platform first used in human studies, to infer gene-level significance of association based on the gene’s distribution of individual variant P-values. Our results contribute to the debate over the relative importance of individual, additive factors and epistatic, or higher order, interactions, in the mapping of genotype to phenotype.

Structural aspects of the aging invertebrate brain

Aging is characterized by a decline in neuronal function in all animal species investigated so far. Functional changes are accompanied by and may be in part caused by, structurally visible degenerative changes in neurons. In the mammalian brain, normal aging shows abnormalities in dendrites and axons, as well as ultrastructural changes in synapses, rather than global neuron loss. The analysis of the structural features of aging neurons, as well as their causal link to molecular mechanisms on the one hand, and the functional decline on the other hand is crucial in order to understand the aging process in the brain. Invertebrate model organisms like Drosophila and C. elegans offer the opportunity to apply a forward genetic approach to the analysis of aging. In the present review, we aim to summarize findings concerning abnormalities in morphology and ultrastructure in invertebrate brains during normal aging and compare them to what is known for the mammalian brain. It becomes clear that despite of their considerably shorter life span, invertebrates display several age-related changes very similar to the mammalian condition, including the retraction of dendritic and axonal branches at specific locations, changes in synaptic density and increased accumulation of presynaptic protein complexes. We anticipate that continued research efforts in invertebrate systems will significantly contribute to reveal (and possibly manipulate) the molecular/cellular pathways leading to neuronal aging in the mammalian brain.

Heat shock proteins and small nucleolar RNAs are dysregulated in a Drosophila model for feline hypertrophic cardiomyopathy

In cats, mutations in myosin binding protein C (encoded by the MYBPC3 gene) have been associated with hypertrophic cardiomyopathy (HCM). However, the molecular mechanisms linking these mutations to HCM remain unknown. Here, we establish Drosophila melanogaster as a model to understand this connection by generating flies harboring MYBPC3 missense mutations (A31P and R820W) associated with feline HCM. The A31P and R820W flies displayed cardiovascular defects in their heart rates and exercise endurance. We used RNA-seq to determine which processes are misregulated in the presence of mutant MYBPC3 alleles. Transcriptome analysis revealed significant downregulation of genes encoding small nucleolar RNA (snoRNAs) in exercised female flies harboring the mutant alleles compared to flies that harbor the wild-type allele. Other processes that were affected included the unfolded protein response and immune/defense responses. These data show that mutant MYBPC3 proteins have widespread effects on the transcriptome of co-regulated genes. Transcriptionally differentially expressed genes are also candidate genes for future evaluation as genetic modifiers of HCM as well as candidate genes for genotype by exercise environment interaction effects on the manifestation of HCM; in cats as well as humans.

Heat shock proteins and small nucleolar RNAs are dysregulated in a Drosophila model for feline hypertrophic cardiomyopathy

In cats, mutations in myosin binding protein C (encoded by the MYBPC3 gene) have been associated with hypertrophic cardiomyopathy (HCM). However, the molecular mechanisms linking these mutations to HCM remain unknown. Here, we establish Drosophila melanogaster as a model to understand this connection by generating flies harboring MYBPC3 missense mutations (A31P and R820W) associated with feline HCM. The A31P and R820W flies displayed cardiovascular defects in their heart rates and exercise endurance. We used RNA-seq to determine which processes are misregulated in the presence of mutant MYBPC3 alleles. Transcriptome analysis revealed significant downregulation of genes encoding small nucleolar RNA (snoRNAs) in exercised female flies harboring the mutant alleles compared to flies that harbor the wild-type allele. Other processes that were affected included the unfolded protein response and immune/defense responses. These data show that mutant MYBPC3 proteins have widespread effects on the transcriptome of co-regulated genes. Transcriptionally differentially expressed genes are also candidate genes for future evaluation as genetic modifiers of HCM as well as candidate genes for genotype by exercise environment interaction effects on the manifestation of HCM; in cats as well as humans.

Leveraging Multiple Layers of Data To Predict Drosophila Complex Traits

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

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

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

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

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

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.

Evolution of Reproductive Behavior

Behaviors associated with reproduction are major contributors to the evolutionary success of organisms and are subject to many evolutionary forces, including natural and sexual selection, and sexual conflict. Successful reproduction involves a range of behaviors, from finding an appropriate mate, courting, and copulation, to the successful production and (in oviparous animals) deposition of eggs following mating. As a consequence, behaviors and genes associated with reproduction are often under strong selection and evolve rapidly. Courtship rituals in flies follow a multimodal pattern, mediated through visual, chemical, tactile, and auditory signals. Premating behaviors allow males and females to assess the species identity, reproductive state, and condition of their partners. Conflicts between the "interests" of individual males, and/or between the reproductive strategies of males and females, often drive the evolution of reproductive behaviors. For example, seminal proteins transmitted by males often show evidence of rapid evolution, mediated by positive selection. Postmating behaviors, including the selection of oviposition sites, are highly variable and Drosophila species span the spectrum from generalists to obligate specialists. Chemical recognition features prominently in adaptation to host plants for feeding and oviposition. Selection acting on variation in pre-, peri-, and postmating behaviors can lead to reproductive isolation and incipient speciation. Response to selection at the genetic level can include the expansion of gene families, such as those for detecting pheromonal cues for mating, or changes in the expression of genes leading to visual cues such as wing spots that are assessed during mating. Here, we consider the evolution of reproductive behavior in Drosophila at two distinct, yet complementary, scales. Some studies take a microevolutionary approach, identifying genes and networks involved in reproduction, and then dissecting the genetics underlying complex behaviors in D. melanogaster Other studies take a macroevolutionary approach, comparing reproductive behaviors across the genus Drosophila and how these might correlate with environmental cues. A full synthesis of this field will require unification across these levels.

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.

Circadian gene variants and the skeletal muscle circadian clock contribute to the evolutionary divergence in longevity across Drosophila populations

Organisms use endogenous clocks to adapt to the rhythmicity of the environment and to synchronize social activities. Although the circadian cycle is implicated in aging, it is unknown whether natural variation in its function contributes to differences in lifespan between populations and whether the circadian clock of specific tissues is key for longevity. We have sequenced the genomes of Drosophila melanogaster strains with exceptional longevity that were obtained via multiple rounds of selection from a parental strain. Comparison of genomic, transcriptomic, and proteomic data revealed that changes in gene expression due to intergenic polymorphisms are associated with longevity and preservation of skeletal muscle function with aging in these strains. Analysis of transcription factors differentially modulated in long-lived versus parental strains indicates a possible role of circadian clock core components. Specifically, there is higher period and timeless and lower cycle expression in the muscle of strains with delayed aging compared to the parental strain. These changes in the levels of circadian clock transcription factors lead to changes in the muscle circadian transcriptome, which includes genes involved in metabolism, proteolysis, and xenobiotic detoxification. Moreover, a skeletal muscle-specific increase in timeless expression extends lifespan and recapitulates some of the transcriptional and circadian changes that differentiate the long-lived from the parental strains. Altogether, these findings indicate that the muscle circadian clock is important for longevity and that circadian gene variants contribute to the evolutionary divergence in longevity across populations.

Evolutionary insights from large scale resequencing datasets in Drosophila melanogaster

Drosophila melanogaster has long been used as an evolutionary model system. Its small genome size, well-annotated genome, and ease of sampling, also makes it a choice species for genome resequencing studies. Hundreds of genomic samples from populations worldwide are available and are currently being used to tackle a wide range of evolutionary questions. In this review, we focused on three insights that have increased our understanding of the evolutionary history of this species, and that have implications for the study of evolutionary processes in other species as well. Because of technical limitations, most of the studies so far have focused on SNP variants. However, long-read sequencing techniques should allow us in the near future to include other type of genomic variants that also influence genome evolution.

Approaches for testing hypotheses for the hypometric scaling of aerobic metabolic rate in animals

Hypometric scaling of aerobic metabolism [larger organisms have lower mass-specific metabolic rates (MR/g)] is nearly universal for interspecific comparisons among animals, yet we lack an agreed upon explanation for this pattern. If physiological constraints on the function of larger animals occur and limit MR/g, these should be observable as direct constraints on animals of extant species and/or as evolved responses to compensate for the proposed constraint. There is evidence for direct constraints and compensatory responses to O2 supply constraint in skin-breathing animals, but not in vertebrates with gas-exchange organs. The duration of food retention in the gut is longer for larger birds and mammals, consistent with a direct constraint on nutrient uptake across the gut wall, but there is little evidence for evolving compensatory responses to gut transport constraints in larger animals. Larger placental mammals (but not marsupials or birds) show evidence of greater challenges with heat dissipation, but there is little evidence for compensatory adaptations to enhance heat loss in larger endotherms, suggesting that metabolic rate (MR) more generally balances heat loss for thermoregulation in endotherms. Size-dependent patterns in many molecular, physiological, and morphological properties are consistent with size-dependent natural selection, such as stronger selection for neurolocomotor performance and growth rate in smaller animals and stronger selection for safety and longevity in larger animals. Hypometric scaling of MR very likely arises from different mechanisms in different taxa and conditions, consistent with the diversity of scaling slopes for MR.

Proper splicing contributes to visual function in the aging Drosophila eye

Changes in splicing patterns are a characteristic of the aging transcriptome; however, it is unclear whether these age-related changes in splicing facilitate the progressive functional decline that defines aging. In Drosophila, visual behavior declines with age and correlates with altered gene expression in photoreceptors, including downregulation of genes encoding splicing factors. Here, we characterized the significance of these age-regulated splicing-associated genes in both splicing and visual function. To do this, we identified differential splicing events in either the entire eye or photoreceptors of young and old flies. Intriguingly, aging photoreceptors show differential splicing of a large number of visual function genes. In addition, as shown previously for aging photoreceptors, aging eyes showed increased accumulation of circular RNAs, which result from noncanonical splicing events. To test whether proper splicing was necessary for visual behavior, we knocked down age-regulated splicing factors in photoreceptors in young flies and examined phototaxis. Notably, many of the age-regulated splicing factors tested were necessary for proper visual behavior. In addition, knockdown of individual splicing factors resulted in changes in both alternative splicing at age-spliced genes and increased accumulation of circular RNAs. Together, these data suggest that cumulative decreases in splicing factor expression could contribute to the differential splicing, circular RNA accumulation, and defective visual behavior observed in aging photoreceptors.

A Cyclin E Centered Genetic Network Contributes to Alcohol-Induced Variation in Drosophila Development

Prenatal exposure to ethanol causes a wide range of adverse physiological, behavioral and cognitive consequences. However, identifying allelic variants and genetic networks associated with variation in susceptibility to prenatal alcohol exposure is challenging in human populations, since time and frequency of exposure and effective dose cannot be determined quantitatively and phenotypic manifestations are diverse. Here, we harnessed the power of natural variation in the Drosophila melanogaster Genetic Reference Panel (DGRP) to identify genes and genetic networks associated with variation in sensitivity to developmental alcohol exposure. We measured development time from egg to adult and viability of 201 DGRP lines reared on regular or ethanol- supplemented medium and identified polymorphisms associated with variation in susceptibility to developmental ethanol exposure. We also documented genotype-dependent variation in sensorimotor behavior after developmental exposure to ethanol using the startle response assay in a subset of 39 DGRP lines. Genes associated with development, including development of the nervous system, featured prominently among genes that harbored variants associated with differential sensitivity to developmental ethanol exposure. Many of them have human orthologs and mutational analyses and RNAi targeting functionally validated a high percentage of candidate genes. Analysis of genetic interaction networks identified Cyclin E (CycE) as a central, highly interconnected hub gene. Cyclin E encodes a protein kinase associated with cell cycle regulation and is prominently expressed in ovaries. Thus, exposure to ethanol during development of Drosophila melanogaster might serve as a genetic model for translational studies on fetal alcohol spectrum disorder.

Charting the genotype-phenotype map: lessons from the Drosophila melanogaster Genetic Reference Panel

Understanding the genetic architecture (causal molecular variants, their effects, and frequencies) of quantitative traits is important for precision agriculture and medicine and predicting adaptive evolution, but is challenging in most species. The Drosophila melanogaster Genetic Reference Panel (DGRP) is a collection of 205 inbred strains with whole genome sequences derived from a single wild population in Raleigh, NC, USA. The large amount of quantitative genetic variation, lack of population structure, and rapid local decay of linkage disequilibrium in the DGRP and outbred populations derived from DGRP lines present a favorable scenario for performing genome-wide association (GWA) mapping analyses to identify candidate causal genes, polymorphisms, and pathways affecting quantitative traits. The many GWA studies utilizing the DGRP have revealed substantial natural genetic variation for all reported traits, little evidence for variants with large effects but enrichment for variants with low P-values, and a tendency for lower frequency variants to have larger effects than more common variants. The variants detected in the GWA analyses rarely overlap those discovered using mutagenesis, and often are the first functional annotations of computationally predicted genes. Variants implicated in GWA analyses typically have sex-specific and genetic background-specific (epistatic) effects, as well as pleiotropic effects on other quantitative traits. Studies in the DGRP reveal substantial genetic control of environmental variation. Taking account of genetic architecture can greatly improve genomic prediction in the DGRP. These features of the genetic architecture of quantitative traits are likely to apply to other species, including humans. WIREs Dev Biol 2018, 7:e289. doi: 10.1002/wdev.289 This article is categorized under: Invertebrate Organogenesis > Flies.

Charting the genotype-phenotype map: lessons from the Drosophila melanogaster Genetic Reference Panel

Understanding the genetic architecture (causal molecular variants, their effects, and frequencies) of quantitative traits is important for precision agriculture and medicine and predicting adaptive evolution, but is challenging in most species. The Drosophila melanogaster Genetic Reference Panel (DGRP) is a collection of 205 inbred strains with whole genome sequences derived from a single wild population in Raleigh, NC, USA. The large amount of quantitative genetic variation, lack of population structure, and rapid local decay of linkage disequilibrium in the DGRP and outbred populations derived from DGRP lines present a favorable scenario for performing genome-wide association (GWA) mapping analyses to identify candidate causal genes, polymorphisms, and pathways affecting quantitative traits. The many GWA studies utilizing the DGRP have revealed substantial natural genetic variation for all reported traits, little evidence for variants with large effects but enrichment for variants with low P-values, and a tendency for lower frequency variants to have larger effects than more common variants. The variants detected in the GWA analyses rarely overlap those discovered using mutagenesis, and often are the first functional annotations of computationally predicted genes. Variants implicated in GWA analyses typically have sex-specific and genetic background-specific (epistatic) effects, as well as pleiotropic effects on other quantitative traits. Studies in the DGRP reveal substantial genetic control of environmental variation. Taking account of genetic architecture can greatly improve genomic prediction in the DGRP. These features of the genetic architecture of quantitative traits are likely to apply to other species, including humans. WIREs Dev Biol 2018, 7:e289. doi: 10.1002/wdev.289 This article is categorized under: Invertebrate Organogenesis > Flies.

Transcriptome profiling of aging Drosophila photoreceptors reveals gene expression trends that correlate with visual senescence

BACKGROUND

Aging is associated with functional decline of neurons and increased incidence of both neurodegenerative and ocular disease. Photoreceptor neurons in Drosophila melanogaster provide a powerful model for studying the molecular changes involved in functional senescence of neurons since decreased visual behavior precedes retinal degeneration. Here, we sought to identify gene expression changes and the genomic features of differentially regulated genes in photoreceptors that contribute to visual senescence.

RESULTS

To identify gene expression changes that could lead to visual senescence, we characterized the aging transcriptome of Drosophila sensory neurons highly enriched for photoreceptors. We profiled the nuclear transcriptome of genetically-labeled photoreceptors over a 40 day time course and identified increased expression of genes involved in stress and DNA damage response, and decreased expression of genes required for neuronal function. We further show that combinations of promoter motifs robustly identify age-regulated genes, suggesting that transcription factors are important in driving expression changes in aging photoreceptors. However, long, highly expressed and heavily spliced genes are also more likely to be downregulated with age, indicating that other mechanisms could contribute to expression changes at these genes. Lastly, we identify that circular RNAs (circRNAs) strongly increase during aging in photoreceptors.

CONCLUSIONS

Overall, we identified changes in gene expression in aging Drosophila photoreceptors that could account for visual senescence. Further, we show that genomic features predict these age-related changes, suggesting potential mechanisms that could be targeted to slow the rate of age-associated visual decline.

Charting the genotype-phenotype map: lessons from the Drosophila melanogaster Genetic Reference Panel

Understanding the genetic architecture (causal molecular variants, their effects, and frequencies) of quantitative traits is important for precision agriculture and medicine and predicting adaptive evolution, but is challenging in most species. The Drosophila melanogaster Genetic Reference Panel (DGRP) is a collection of 205 inbred strains with whole genome sequences derived from a single wild population in Raleigh, NC, USA. The large amount of quantitative genetic variation, lack of population structure, and rapid local decay of linkage disequilibrium in the DGRP and outbred populations derived from DGRP lines present a favorable scenario for performing genome-wide association (GWA) mapping analyses to identify candidate causal genes, polymorphisms, and pathways affecting quantitative traits. The many GWA studies utilizing the DGRP have revealed substantial natural genetic variation for all reported traits, little evidence for variants with large effects but enrichment for variants with low P-values, and a tendency for lower frequency variants to have larger effects than more common variants. The variants detected in the GWA analyses rarely overlap those discovered using mutagenesis, and often are the first functional annotations of computationally predicted genes. Variants implicated in GWA analyses typically have sex-specific and genetic background-specific (epistatic) effects, as well as pleiotropic effects on other quantitative traits. Studies in the DGRP reveal substantial genetic control of environmental variation. Taking account of genetic architecture can greatly improve genomic prediction in the DGRP. These features of the genetic architecture of quantitative traits are likely to apply to other species, including humans. WIREs Dev Biol 2018, 7:e289. doi: 10.1002/wdev.289 This article is categorized under: Invertebrate Organogenesis > Flies.

Genetics of alcohol consumption in Drosophila melanogaster

Individual variation in alcohol consumption in human populations is determined by genetic, environmental, social and cultural factors. In contrast to humans, genetic contributions to complex behavioral phenotypes can be readily dissected in Drosophila, where both the genetic background and environment can be controlled and behaviors quantified through simple high-throughput assays. Here, we measured voluntary consumption of ethanol in ∼3000 individuals of each sex from an advanced intercross population derived from 37 lines of the Drosophila melanogaster Genetic Reference Panel. Extreme quantitative trait loci mapping identified 385 differentially segregating allelic variants located in or near 291 genes at P < 10^-8 . The effects of single nucleotide polymorphisms associated with voluntary ethanol consumption are sex-specific, as found for other alcohol-related phenotypes. To assess causality, we used RNA interference knockdown or P{MiET1} mutants and their corresponding controls and functionally validated 86% of candidate genes in at least one sex. We constructed a genetic network comprised of 23 genes along with a separate trio and a pair of connected genes. Gene ontology analyses showed enrichment of developmental genes, including development of the nervous system. Furthermore, a network of human orthologs showed enrichment for signal transduction processes, protein metabolism and developmental processes, including nervous system development. Our results show that the genetic architecture that underlies variation in voluntary ethanol consumption is sexually dimorphic and partially overlaps with genetic factors that control variation in feeding behavior and alcohol sensitivity. This integrative genetic architecture is rooted in evolutionarily conserved features that can be extrapolated to human genetic interaction networks.

Armour GE, Thakkar E, Mackay TFC, Anholt RRH. 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.

Although physiological effects of environmental toxins are well documented, we know little about the genetic factors that determine individual variation in susceptibility to toxins. Such information is difficult to obtain in human populations due to heterogeneity in genetic background and environmental exposure, and the diversity of symptoms and time lag with which they appear after toxic exposure. Here, we show that the fruit fly, Drosophila, can serve as a powerful genetic model system to elucidate the genetic underpinnings that contribute to individual variation in resistance to toxicity, using lead and cadmium exposure as an experimental paradigm. We identified genes that harbor genetic variants that contribute to individual variation in resistance to heavy metal exposure. Furthermore, we constructed genetic networks on which we could superimpose human counterparts of Drosophila genes. We were able to place human genes previously implicated in heavy metal toxicity in biological context and identify novel targets for heavy metal toxicity. Thus, we demonstrate that based on evolutionary conservation of fundamental biological processes, we can use Drosophila as a powerful translational model for toxicogenomics studies.

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.

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.