Genetic basis of natural variation in body pigmentation in Drosophila melanogaster

Body pigmentation in insects and other organisms is typically variable within and between species and is often associated with fitness. Regulatory variants with large effects at bab1, t and e affect variation in abdominal pigmentation in several populations of Drosophila melanogaster. Recently, we performed a genome wide association (GWA) analysis of variation in abdominal pigmentation using the inbred, sequenced lines of the Drosophila Genetic Reference Panel (DGRP). We confirmed the large effects of regulatory variants in bab1, t and e; identified 81 additional candidate genes; and validated 17 candidate genes (out of 28 tested) using RNAi knockdown of gene expression and mutant alleles. However, these analyses are imperfect proxies for the effects of segregating variants. Here, we describe the results of an extreme quantitative trait locus (xQTL) GWA analysis of female body pigmentation in an outbred population derived from light and dark DGRP lines. We replicated the effects on pigmentation of 28 genes implicated by the DGRP GWA study, including bab1, t and e and 7 genes previously validated by RNAi and/or mutant analyses. We also identified many additional loci. The genetic architecture of Drosophila pigmentation is complex, with a few major genes and many other loci with smaller effects.     

Identification of a genetic network for an ecologically relevant behavioural phenotype in Drosophila melanogaster

Pupation site choice of Drosophila third‐instar larvae is critical for the survival of individuals, as pupae are exposed to various biotic and abiotic dangers while immobilized during the 3–4 days of metamorphosis. This singular behavioural choice is sensitive to both environmental and genetic factors. Here, we developed a high‐throughput phenotyping approach to assay the variation in pupation height in Drosophila melanogaster, while controlling for possibly confounding factors. We find substantial variation of mean pupation height among sampled natural stocks and we show that the Drosophila Genetic Reference Panel (DGRP) reflects this variation. Using the DGRP stocks for genome‐wide association (GWA) mapping, 16 loci involved in determining pupation height could be resolved. The candidate genes in these loci are enriched for high expression in the larval central nervous system. A genetic network could be constructed from the candidate loci, which places scribble (scrib) at the centre, plus other genes known to be involved in nervous system development, such as Epidermal growth factor receptor (Egfr) and p53. Using gene disruption lines, we could functionally validate several of the initially identified loci, as well as additional loci predicted from network analysis. Our study shows that the combination of high‐throughput phenotyping with a genetic analysis of variation captured from the wild can be used to approach the genetic dissection of an environmentally relevant behavioural phenotype.     

Two genomic regions together cause dark abdominal pigmentation in Drosophila tenebrosa

Pigmentation is a rapidly evolving trait that is under both natural and sexual selection in many organisms. In the quinaria group of Drosophila, nearly all of the 30 species have an abdomen that is light in color with distinct markings; D. tenebrosa is the exception in that it has a completely melanic abdomen with no visible markings. In this study, we use a combination of quantitative genetic and candidate gene approaches to investigate the genetic basis of abdominal pigmentation in D. tenebrosa. We find that abdominal pigmentation is invariant across wild-caught lines of D. tenebrosa and is not sexually dimorphic. Quantitative genetic mapping utilizing crosses between D. tenebrosa and the light-colored D. suboccidentalis indicates that two genomic regions together underlie abdominal pigmentation, including the X-chromosome and an autosome (Muller Element C/E). Further support for their central importance in pigmentation is that experimental introgression of one phenotype into the other species, in either direction, results in introgression of these two genomic regions. Finally, the expression of the X-linked gene yellow in the pupae exactly foreshadows the adult melanization pattern in the abdomen of both species, suggesting that changes in the regulation of yellow are important for the phenotypic divergence of D. tenebrosa from the rest of the quinaria group. These results contribute to a body of work that demonstrates how changes in expression of highly conserved genes can cause substantial phenotypic differences even between closely related species.     

Epistatic partners of neurogenic genes modulate Drosophila olfactory behavior

The extent to which epistasis affects the genetic architecture of complex traits is difficult to quantify, and identifying variants in natural populations with epistatic interactions is challenging. Previous studies in Drosophila implicated extensive epistasis between variants in genes that affect neural connectivity and contribute to natural variation in olfactory response to benzaldehyde. In this study, we implemented a powerful screen to quantify the extent of epistasis as well as identify candidate interacting variants using 203 inbred wild‐derived lines with sequenced genomes of the Drosophila melanogaster Genetic Reference Panel (DGRP). We crossed the DGRP lines to P[GT1]‐element insertion mutants in Sema‐5c and neuralized (neur), two neurodevelopmental loci which affect olfactory behavior, and to their coisogenic wild‐type control. We observed significant variation in olfactory responses to benzaldehyde among F₁ genotypes and for the DGRP line by mutant genotype interactions for both loci, showing extensive nonadditive genetic variation. We performed genome‐wide association analyses to identify the candidate modifier loci. None of these polymorphisms were in or near the focal genes; therefore, epistasis is the cause of the nonadditive genetic variance. Candidate genes could be placed in interaction networks. Several candidate modifiers are associated with neural development. Analyses of mutants of candidate epistatic partners with neur (merry‐go‐round (mgr), prospero (pros), CG10098, Alhambra (Alh) and CG12535) and Sema‐5c (CG42540 and bruchpilot (brp)) showed aberrant olfactory responses compared with coisogenic controls. Thus, integrating genome‐wide analyses of natural variants with mutations at defined genomic locations in a common coisogenic background can unmask specific epistatic modifiers of behavioral phenotypes.     

Genome-Wide Association Analysis of Tolerance to Methylmercury Toxicity in Drosophila Implicates Myogenic and Neuromuscular Developmental Pathways

Methylmercury (MeHg) is a persistent environmental toxin present in seafood that can compromise the developing nervous system in humans. The effects of MeHg toxicity varies among individuals, despite similar levels of exposure, indicating that genetic differences contribute to MeHg susceptibility. To examine how genetic variation impacts MeHg tolerance, we assessed developmental tolerance to MeHg using the sequenced, inbred lines of the Drosophila melanogaster Genetic Reference Panel (DGRP). We found significant genetic variation in the effects of MeHg on development, measured by eclosion rate, giving a broad sense heritability of 0.86. To investigate the influence of dietary factors, we measured MeHg toxicity with caffeine supplementation in the DGRP lines. We found that caffeine counteracts the deleterious effects of MeHg in the majority of lines, and there is significant genetic variance in the magnitude of this effect, with a broad sense heritability of 0.80. We performed genome-wide association (GWA) analysis for both traits, and identified candidate genes that fall into several gene ontology categories, with enrichment for genes involved in muscle and neuromuscular development. Overexpression of glutamate-cysteine ligase, a MeHg protective enzyme, in a muscle-specific manner leads to a robust rescue of eclosion of flies reared on MeHg food. Conversely, mutations in kirre, a pivotal myogenic gene identified in our GWA analyses, modulate tolerance to MeHg during development in accordance with kirre expression levels. Finally, we observe disruptions of indirect flight muscle morphogenesis in MeHg-exposed pupae. Since the pathways for muscle development are evolutionarily conserved, it is likely that the effects of MeHg observed in Drosophila can be generalized across phyla, implicating muscle as an additional hitherto unrecognized target for MeHg toxicity. Furthermore, our observations that caffeine can ameliorate the toxic effects of MeHg show that nutritional factors and dietary manipulations may offer protection against the deleterious effects of MeHg exposure.     

Dissecting the Biological Role of Mucin-type O-Glycosylation Using RNA Interference in Drosophila Cell Culture

Mucin type O-glycosylation is a highly conserved form of post-translational modification initiated by the family of enzymes known as the polypeptide α-N-acetylgalactosaminyltransferases (ppGalNAcTs in mammals and PGANTs in Drosophila). To address the cellular functions of the many PGANT family members, RNA interference (RNAi) to each pgant gene was performed in two independent Drosophila cell culture lines. We demonstrate that RNAi to individual pgant genes results in specific reduction in gene expression without affecting the expression of other family members. Cells with reduced expression of individual pgant genes were then examined for changes in viability, morphology, adhesion, and secretion to assess the contribution of each family member to these cellular functions. Here we find that RNAi to pgant3, pgant6, or pgant7 resulted in reduced secretion, further supporting a role for O-glycosylation in proper secretion. Additionally, RNAi to pgant3 or pgant6 resulted in altered Golgi organization, suggesting a role for each in establishing or maintaining proper secretory apparatus structure. Other subcellular effects observed included multinucleated cells seen after RNAi to either pgant2 or pgant35A, suggesting a role for these genes in the completion of cytokinesis. These studies demonstrate the efficient and specific knockdown of pgant gene expression in two Drosophila cell culture systems, resulting in specific morphological and functional effects. Our work provides new information regarding the biological roles of O-glycosylation and illustrates a new platform for interrogating the cellular and subcellular effects of this form of post-translational modification.     

Quantitative Genetics of Food Intake in Drosophila melanogaster

Food intake is an essential animal activity, regulated by neural circuits that motivate food localization, evaluate nutritional content and acceptance or rejection responses through the gustatory system, and regulate neuroendocrine feedback loops that maintain energy homeostasis. Excess food consumption in people is associated with obesity and metabolic and cardiovascular disorders. However, little is known about the genetic basis of natural variation in food consumption. To gain insights in evolutionarily conserved genetic principles that regulate food intake, we took advantage of a model system, Drosophila melanogaster, in which food intake, environmental conditions and genetic background can be controlled precisely. We quantified variation in food intake among 182 inbred, sequenced lines of the Drosophila melanogaster Genetic Reference Panel (DGRP). We found significant genetic variation in the mean and within-line environmental variance of food consumption and observed sexual dimorphism and genetic variation in sexual dimorphism for both food intake traits (mean and variance). We performed genome wide association (GWA) analyses for mean food intake and environmental variance of food intake (using the coefficient of environmental variation, CV _E, as the metric for environmental variance) and identified molecular polymorphisms associated with both traits. Validation experiments using RNAi-knockdown confirmed 24 of 31 (77%) candidate genes affecting food intake and/or variance of food intake, and a test cross between selected DGRP lines confirmed a SNP affecting mean food intake identified in the GWA analysis. The majority of the validated candidate genes were novel with respect to feeding behavior, and many had mammalian orthologs implicated in metabolic diseases.     

Identifying Loci Contributing to Natural Variation in Xenobiotic Resistance in Drosophila

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

Natural genetic variation in Drosophila melanogaster reveals genes associated with Coxiella burnetii infection

The gram-negative bacterium Coxiella burnetii is the causative agent of Query (Q) fever in humans and coxiellosis in livestock. Host genetics are associated with C. burnetii pathogenesis both in humans and animals; however, it remains unknown if specific genes are associated with severity of infection. We employed the Drosophila Genetics Reference Panel to perform a genome-wide association study to identify host genetic variants that affect host survival to C. burnetii infection. The genome-wide association study identified 64 unique variants (P < 10⁻⁵) associated with 25 candidate genes. We examined the role each candidate gene contributes to host survival during C. burnetii infection using flies carrying a null mutation or RNAi knockdown of each candidate. We validated 15 of the 25 candidate genes using at least one method. This is the first report establishing involvement of many of these genes or their homologs with C. burnetii susceptibility in any system. Among the validated genes, FER and tara play roles in the JAK/STAT, JNK, and decapentaplegic/TGF-β signaling pathways which are components of known innate immune responses to C. burnetii infection. CG42673 and DIP-ε play roles in bacterial infection and synaptic signaling but have no previous association with C. burnetii pathogenesis. Furthermore, since the mammalian ortholog of CG13404 (PLGRKT) is an important regulator of macrophage function, CG13404 could play a role in host susceptibility to C. burnetii through hemocyte regulation. These insights provide a foundation for further investigation regarding the genetics of C. burnetii susceptibility across a wide variety of hosts.     

Synergistic Interactions between Drosophila Orthologues of Genes Spanned by De Novo Human CNVs Support Multiple-Hit Models of Autism

Autism spectrum disorders (ASDs) are highly heritable and characterised by deficits in social interaction and communication, as well as restricted and repetitive behaviours. Although a number of highly penetrant ASD gene variants have been identified, there is growing evidence to support a causal role for combinatorial effects arising from the contributions of multiple loci. By examining synaptic and circadian neurological phenotypes resulting from the dosage variants of unique human:fly orthologues in Drosophila, we observe numerous synergistic interactions between pairs of informatically-identified candidate genes whose orthologues are jointly affected by large de novo copy number variants (CNVs). These CNVs were found in the genomes of individuals with autism, including a patient carrying a 22q11.2 deletion. We first demonstrate that dosage alterations of the unique Drosophila orthologues of candidate genes from de novo CNVs that harbour only a single candidate gene display neurological defects similar to those previously reported in Drosophila models of ASD-associated variants. We then considered pairwise dosage changes within the set of orthologues of candidate genes that were affected by the same single human de novo CNV. For three of four CNVs with complete orthologous relationships, we observed significant synergistic effects following the simultaneous dosage change of gene pairs drawn from a single CNV. The phenotypic variation observed at the Drosophila synapse that results from these interacting genetic variants supports a concordant phenotypic outcome across all interacting gene pairs following the direction of human gene copy number change. We observe both specificity and transitivity between interactors, both within and between CNV candidate gene sets, supporting shared and distinct genetic aetiologies. We then show that different interactions affect divergent synaptic processes, demonstrating distinct molecular aetiologies. Our study illustrates mechanisms through which synergistic effects resulting from large structural variation can contribute to human disease.     

Large Neurological Component to Genetic Differences Underlying Biased Sperm Use in Drosophila

Sperm competition arises as a result of complex interactions among male and female factors. While the roles of some male factors are known, little is known of the molecules or mechanisms that underlie the female contribution to sperm competition. The genetic tools available for Drosophila allow us to identify, in an unbiased manner, candidate female genes that are critical for mediating sperm competition outcomes. We first screened for differences in female sperm storage and use patterns by characterizing the natural variation in sperm competition in a set of 39 lines from the sequenced Drosophila Genetic Reference Panel (DGRP) of wild-derived inbred lines. We found extensive female variation in sperm competition outcomes. To generate a list of candidate female genes for functional studies, we performed a genome-wide association mapping, utilizing the common single-nucleotide polymorphisms (SNPs) segregating in the DGRP lines. Surprisingly, SNPs within ion channel genes and other genes with roles in the nervous system were among the top associated SNPs. Knockdown studies of three candidate genes (para, Rab2, and Rim) in sensory neurons innervating the female reproductive tract indicate that some of these candidate female genes may affect sperm competition by modulating the neural input of these sensory neurons to the female reproductive tract. More extensive functional studies are needed to elucidate the exact role of all these candidate female genes in sperm competition. Nevertheless, the female nervous system appears to have a previously unappreciated role in sperm competition. Our results indicate that the study of female control of sperm competition should not be limited to female reproductive tract-specific genes, but should focus also on diverse biological pathways.     

Genome-Wide Analysis Reveals Novel Regulators of Growth in Drosophila melanogaster

Organismal size depends on the interplay between genetic and environmental factors. Genome-wide association (GWA) analyses in humans have implied many genes in the control of height but suffer from the inability to control the environment. Genetic analyses in Drosophila have identified conserved signaling pathways controlling size; however, how these pathways control phenotypic diversity is unclear. We performed GWA of size traits using the Drosophila Genetic Reference Panel of inbred, sequenced lines. We find that the top associated variants differ between traits and sexes; do not map to canonical growth pathway genes, but can be linked to these by epistasis analysis; and are enriched for genes and putative enhancers. Performing GWA on well-studied developmental traits under controlled conditions expands our understanding of developmental processes underlying phenotypic diversity.     

Effector-Triggered Immune Response in Arabidopsis thaliana Is a Quantitative Trait

We identified loci responsible for natural variation in Arabidopsis thaliana (Arabidopsis) responses to a bacterial pathogen virulence factor, HopAM1. HopAM1 is a type III effector protein secreted by the virulent Pseudomonas syringae strain Pto DC3000. Delivery of HopAM1 from disarmed Pseudomonas strains leads to local cell death, meristem chlorosis, or both, with varying intensities in different Arabidopsis accessions. These phenotypes are not associated with differences in bacterial growth restriction. We treated the two phenotypes as quantitative traits to identify host loci controlling responses to HopAM1. Genome-wide association (GWA) of 64 Arabidopsis accessions identified independent variants highly correlated with response to each phenotype. Quantitative trait locus (QTL) mapping in a recombinant inbred population between Bur-0 and Col-0 accessions revealed genetic linkage to regions distinct from the top GWA hits. Two major QTL associated with HopAM1-induced cell death were also associated with HopAM1-induced chlorosis. HopAM1-induced changes in Arabidopsis gene expression showed that rapid HopAM1-dependent cell death in Bur-0 is correlated with effector-triggered immune responses. Studies of the effect of mutations in known plant immune system genes showed, surprisingly, that both cell death and chlorosis phenotypes are enhanced by loss of EDS1, a regulatory hub in the plant immune-signaling network. Our results reveal complex genetic architecture for response to this particular type III virulence effector, in contrast to the typical monogenic control of cell death and disease resistance triggered by most type III effectors.     

UP-TORR: Online Tool for Accurate and Up-to-Date Annotation of RNAi Reagents

RNA interference (RNAi) is a widely adopted tool for loss-of-function studies but RNAi results only have biological relevance if the reagents are appropriately mapped to genes. Several groups have designed and generated RNAi reagent libraries for studies in cells or in vivo for Drosophila and other species. At first glance, matching RNAi reagents to genes appears to be a simple problem, as each reagent is typically designed to target a single gene. In practice, however, the reagent–gene relationship is complex. Although the sequences of oligonucleotides used to generate most types of RNAi reagents are static, the reference genome and gene annotations are regularly updated. Thus, at the time a researcher chooses an RNAi reagent or analyzes RNAi data, the most current interpretation of the RNAi reagent–gene relationship, as well as related information regarding specificity (e.g., predicted off-target effects), can be different from the original interpretation. Here, we describe a set of strategies and an accompanying online tool, UP-TORR (for Updated Targets of RNAi Reagents; www.flyrnai.org/up-torr), useful for accurate and up-to-date annotation of cell-based and in vivo RNAi reagents. Importantly, UP-TORR automatically synchronizes with gene annotations daily, retrieving the most current information available, and for Drosophila, also synchronizes with the major reagent collections. Thus, UP-TORR allows users to choose the most appropriate RNAi reagents at the onset of a study, as well as to perform the most appropriate analyses of results of RNAi-based studies.     

The Tribolium ortholog of knirps and knirps-related is crucial for head segmentation but plays a minor role during abdominal patterning

Segment formation in the long germ insect Drosophila is dominated by overlapping gap gene domains in the syncytial blastoderm. In the short germ beetle Tribolium castaneum abdominal segments arise from a cellular growth zone, implying different patterning mechanisms. We describe here the single Tribolium ortholog of the Drosophila genes knirps and knirps-related (called Tc-knirps). Tc-knirps expression is conserved during head patterning and at later stages. However, posterior Tc-knirps expression in the ectoderm is limited to a stripe in A1, instead of a broad abdominal domain covering segment primordia A2-A5 as in Drosophila. Tc-knirps RNAi yields only mild defects in the abdomen, at a position posterior to the abdominal Tc-knirps domain. In addition, Tc-knirps RNAi larvae lack the antennal and mandibular segments. These defects are much more severe than the head defects caused by combined inactivation of Dm-knirps and Dm-knirps-related. Our findings support the notion that the role of gap gene homologs in abdominal segmentation differs fundamentally in long and short germ insects. Moreover, the pivotal role of Tc-knirps in the head suggests an ancestral role for knirps as head patterning gene. Based on this RNAi analysis, Tc-knirps functions neither in the head nor the abdomen as a canonical gap gene.     

Abd-B suppresses lepidopteran proleg development in posterior abdomen

Pterygotes lack abdominal appendages except for pleuropods and prolegs. The larvae of some holometabolous insects develop prolegs, which are used for locomotion. We analyzed the role of the homeotic genes abd-A and Abd-B in lepidopteran proleg development using mutant analysis and embryonic RNAi in the silkworm Bombyx mori. The E^Mu mutant developed extra prolegs in its posterior abdomen and showed the misexpression of both genes, suggesting their involvement in proleg formation. The depletion of Abd-B by embryonic RNAi caused the development of extra prolegs on all segments posterior to A6, indicating the suppressive function of Abd-B. The abd-A RNAi animals failed to develop prolegs. These results indicate that abd-A and Abd-B are involved in proleg development in B. mori.     

Genetic Architecture and Functional Characterization of Genes Underlying the Rapid Diversification of Male External Genitalia Between Drosophila simulans and Drosophila mauritiana

Male sexual characters are often among the first traits to diverge between closely related species and identifying the genetic basis of such changes can contribute to our understanding of their evolutionary history. However, little is known about the genetic architecture or the specific genes underlying the evolution of male genitalia. The morphology of the claspers, posterior lobes, and anal plates exhibit striking differences between Drosophila mauritiana and D. simulans. Using QTL and introgression-based high-resolution mapping, we identified several small regions on chromosome arms 3L and 3R that contribute to differences in these traits. However, we found that the loci underlying the evolution of clasper differences between these two species are independent from those that contribute to posterior lobe and anal plate divergence. Furthermore, while most of the loci affect each trait in the same direction and act additively, we also found evidence for epistasis between loci for clasper bristle number. In addition, we conducted an RNAi screen in D. melanogaster to investigate if positional and expression candidate genes located on chromosome 3L, are also involved in genital development. We found that six of these genes, including components of Wnt signaling and male-specific lethal 3 (msl3), regulate the development of genital traits consistent with the effects of the introgressed regions where they are located and that thus represent promising candidate genes for the evolution these traits.     

Divergent and conserved roles of extradenticle in body segmentation and appendage formation, respectively, in the cricket Gryllus bimaculatus

The cricket Gryllus bimaculatus is a typical hemimetabolous intermediate germ insect, in which the processes of segmentation and appendage formation differ from those in Drosophila, a holometabolous long germ insect. In order to compare their developmental mechanisms, we have focused on Gryllus orthologs of the Drosophila developmental regulatory genes and studied their functions. Here, we report a functional analysis of the Gryllus ortholog of extradenticle (Gb′exd) using embryonic and parental RNA interference (RNAi) techniques. We found the following: (1) RNAi suppression of Gb′exd results in the deletion or fusion of body segments. Especially the head was often very severely affected. This gap-like phenotype may be related to reduced expression of the gap genes hunchback and Krüppel in early RNAi germbands. (2) In the appendages, several segments (podomeres) were fused. (3) Head appendages including the antenna were transformed to a leg-like structure consisting of at least one proximal podomere as well as several tarsomeres. The defects in appendages are reminiscent of the phenotype caused by large exd clones in Drosophila antennal discs. These findings led us to the conclusion that (1) Gb′exd is required for segment patterning in the gnathal to abdominal region, acting in a gap gene-like manner in the anterior region. (2) Gb′exd plays important roles in formation of the appendages and the determination of their identities, acting as a regulatory switch that chooses between the fates of head appendages versus the appendage ground state. Although functions of Gb′exd in appendage patterning appear fundamentally conserved between Gryllus and Drosophila, its role in body segmentation may differ from that of Drosophila exd.