Protocol for preparing Drosophila genomic DNA to create chromosome-level de novo genome assemblies.

De novo genome assemblies are common tools for examining novel biological phenomena in non-model organisms. Here, we present a protocol for preparing Drosophila genomic DNA to create chromosome-level de novo genome assemblies. We describe steps for high-molecular-weight DNA preparation with phenol or Genomic-tips, quality control, long-read nanopore sequencing, short-read DNA library preparation, and sequencing. We then detail procedures of genome assembly, annotation, and assessment that can be used for downstream comparison and functional analysis. For complete details on the use and execution of this protocol, please refer to Sperling et al.1.

A genetic basis for facultative parthenogenesis in Drosophila.

Facultative parthenogenesis enables sexually reproducing organisms to switch between sexual and asexual parthenogenetic reproduction. To gain insights into this phenomenon, we sequenced the genomes of sexually reproducing and parthenogenetic strains of Drosophila mercatorum and identified differences in the gene expression in their eggs. We then tested whether manipulating the expression of candidate gene homologs identified in Drosophila mercatorum could lead to facultative parthenogenesis in the non-parthenogenetic species Drosophila melanogaster. This identified a polygenic system whereby increased expression of the mitotic protein kinase polo and decreased expression of a desaturase, Desat2, caused facultative parthenogenesis in the non-parthenogenetic species that was enhanced by increased expression of Myc. The genetically induced parthenogenetic Drosophila melanogaster eggs exhibit de novo centrosome formation, fusion of the meiotic products, and the onset of development to generate predominantly triploid offspring. Thus, we demonstrate a genetic basis for sporadic facultative parthenogenesis in an animal.

A novel transposable element-mediated mechanism causes antiviral resistance in Drosophila through truncating the Veneno protein.

Hosts are continually selected to evolve new defenses against an ever-changing array of pathogens. To understand this process, we examined the genetic basis of resistance to the Drosophila A virus in Drosophila melanogaster. In a natural population, we identified a polymorphic transposable element (TE) insertion that was associated with an ∼19,000-fold reduction in viral titers, allowing flies to largely escape the harmful effects of infection by this virulent pathogen. The insertion occurs in the protein-coding sequence of the gene Veneno, which encodes a Tudor domain protein. By mutating Veneno with CRISPR-Cas9 in flies and expressing it in cultured cells, we show that the ancestral allele of the gene has no effect on viral replication. Instead, the TE insertion is a gain-of-function mutation that creates a gene encoding a novel resistance factor. Viral titers remained reduced when we deleted the TE sequence from the transcript, indicating that resistance results from the TE truncating the Veneno protein. This is a novel mechanism of virus resistance and a new way by which TEs can contribute to adaptation.

Drosophila Evolution over Space and Time (DEST): A New Population Genomics Resource.

Drosophila melanogaster is a leading model in population genetics and genomics, and a growing number of whole-genome data sets from natural populations of this species have been published over the last years. A major challenge is the integration of disparate data sets, often generated using different sequencing technologies and bioinformatic pipelines, which hampers our ability to address questions about the evolution of this species. Here we address these issues by developing a bioinformatics pipeline that maps pooled sequencing (Pool-Seq) reads from D. melanogaster to a hologenome consisting of fly and symbiont genomes and estimates allele frequencies using either a heuristic (PoolSNP) or a probabilistic variant caller (SNAPE-pooled). We use this pipeline to generate the largest data repository of genomic data available for D. melanogaster to date, encompassing 271 previously published and unpublished population samples from over 100 locations in >20 countries on four continents. Several of these locations have been sampled at different seasons across multiple years. This data set, which we call Drosophila Evolution over Space and Time (DEST), is coupled with sampling and environmental metadata. A web-based genome browser and web portal provide easy access to the SNP data set. We further provide guidelines on how to use Pool-Seq data for model-based demographic inference. Our aim is to provide this scalable platform as a community resource which can be easily extended via future efforts for an even more extensive cosmopolitan data set. Our resource will enable population geneticists to analyze spatiotemporal genetic patterns and evolutionary dynamics of D. melanogaster populations in unprecedented detail.

Broad geographic sampling reveals the shared basis and environmental correlates of seasonal adaptation in Drosophila.

To advance our understanding of adaptation to temporally varying selection pressures, we identified signatures of seasonal adaptation occurring in parallel among Drosophila melanogaster populations. Specifically, we estimated allele frequencies genome-wide from flies sampled early and late in the growing season from 20 widely dispersed populations. We identified parallel seasonal allele frequency shifts across North America and Europe, demonstrating that seasonal adaptation is a general phenomenon of temperate fly populations. Seasonally fluctuating polymorphisms are enriched in large chromosomal inversions, and we find a broad concordance between seasonal and spatial allele frequency change. The direction of allele frequency change at seasonally variable polymorphisms can be predicted by weather conditions in the weeks prior to sampling, linking the environment and the genomic response to selection. Our results suggest that fluctuating selection is an important evolutionary force affecting patterns of genetic variation in Drosophila.

Common genetic associations between age-related diseases.

Age is a common risk factor in many diseases, but the molecular basis for this relationship is elusive. In this study we identified 4 disease clusters from 116 diseases in the UK Biobank data, defined by their age-of-onset profiles, and found that diseases with the same onset profile are genetically more similar, suggesting a common etiology. This similarity was not explained by disease categories, co-occurrences or disease cause-effect relationships. Two of the four disease clusters had an increased risk of occurrence from age 20 and 40 years respectively. They both showed an association with known aging-related genes, yet differed in functional enrichment and evolutionary profiles. Moreover, they both had age-related expression and methylation changes. We also tested mutation accumulation and antagonistic pleiotropy theories of aging and found support for both.

Functional conservation in genes and pathways linking ageing and immunity.

At first glance, longevity and immunity appear to be different traits that have not much in common except the fact that the immune system promotes survival upon pathogenic infection. Substantial evidence however points to a molecularly intertwined relationship between the immune system and ageing. Although this link is well-known throughout the animal kingdom, its genetic basis is complex and still poorly understood. To address this question, we here provide a compilation of all genes concomitantly known to be involved in immunity and ageing in humans and three well-studied model organisms, the nematode worm Caenorhabditis elegans, the fruit fly Drosophila melanogaster, and the house mouse Mus musculus. By analysing human orthologs among these species, we identified 7 evolutionarily conserved signalling cascades, the insulin/TOR network, three MAPK (ERK, p38, JNK), JAK/STAT, TGF-ß, and Nf-κB pathways that act pleiotropically on ageing and immunity. We review current evidence for these pathways linking immunity and lifespan, and their role in the detrimental dysregulation of the immune system with age, known as immunosenescence. We argue that the phenotypic effects of these pathways are often context-dependent and vary, for example, between tissues, sexes, and types of pathogenic infection. Future research therefore needs to explore a higher temporal, spatial and environmental resolution to fully comprehend the connection between ageing and immunity.

Allelic polymorphism at foxo contributes to local adaptation in Drosophila melanogaster.

The insulin/insulin-like growth factor signalling pathway has been hypothesized as a major determinant of life-history profiles that vary adaptively in natural populations. In Drosophila melanogaster, multiple components of this pathway vary predictably with latitude; this includes foxo, a conserved gene that regulates insulin signalling and has pleiotropic effects on a variety of fitness-associated traits. We hypothesized that allelic variation at foxo contributes to genetic variance for size-related traits that vary adaptively with latitude. We first examined patterns of variation among natural populations along a latitudinal transect in the eastern United States and show that thorax length, wing area, wing loading, and starvation tolerance exhibit significant latitudinal clines for both males and females but that development time does not vary predictably with latitude. We then generated recombinant outbred populations and show that naturally occurring allelic variation at foxo, which exhibits stronger clinality than expected, is associated with the same traits that vary with latitude in the natural populations. Our results suggest that allelic variation at foxo contributes to adaptive patterns of life-history variation in natural populations of this genetic model.

Transposable Element Landscape in Drosophila Populations Selected for Longevity.

Transposable elements (TEs) inflict numerous negative effects on health and fitness as they replicate by integrating into new regions of the host genome. Even though organisms employ powerful mechanisms to demobilize TEs, transposons gradually lose repression during aging. The rising TE activity causes genomic instability and was implicated in age-dependent neurodegenerative diseases, inflammation, and the determination of lifespan. It is therefore conceivable that long-lived individuals have improved TE silencing mechanisms resulting in reduced TE expression relative to their shorter-lived counterparts and fewer genomic insertions. Here, we test this hypothesis by performing the first genome-wide analysis of TE insertions and expression in populations of Drosophila melanogaster selected for longevity through late-life reproduction for 50-170 generations from four independent studies. Contrary to our expectation, TE families were generally more abundant in long-lived populations compared with nonselected controls. Although simulations showed that this was not expected under neutrality, we found little evidence for selection driving TE abundance differences. Additional RNA-seq analysis revealed a tendency for reducing TE expression in selected populations, which might be more important for lifespan than regulating genomic insertions. We further find limited evidence of parallel selection on genes related to TE regulation and transposition. However, telomeric TEs were genomically and transcriptionally more abundant in long-lived flies, suggesting improved telomere maintenance as a promising TE-mediated mechanism for prolonging lifespan. Our results provide a novel viewpoint indicating that reproduction at old age increases the opportunity of TEs to be passed on to the next generation with little impact on longevity.

Transcriptomic profiling of long- and short-lived mutant mice implicates mitochondrial metabolism in ageing and shows signatures of normal ageing in progeroid mice.

Genetically modified mouse models of ageing are the living proof that lifespan and healthspan can be lengthened or shortened, and provide a powerful context in which to unravel the molecular mechanisms at work. In this study, we analysed and compared gene expression data from 10 long-lived and 8 short-lived mouse models of ageing. Transcriptome-wide correlation analysis revealed that mutations with equivalent effects on lifespan induce more similar transcriptomic changes, especially if they target the same pathway. Using functional enrichment analysis, we identified 58 gene sets with consistent changes in long- and short-lived mice, 55 of which were up-regulated in long-lived mice and down-regulated in short-lived mice. Half of these sets represented genes involved in energy and lipid metabolism, among which Ppargc1a, Mif, Aldh5a1 and Idh1 were frequently observed. Based on the gene sets with consistent changes, and also the whole transcriptome, the gene expression changes during normal ageing resembled the transcriptome of short-lived models, suggesting that accelerated ageing models reproduce partially the molecular changes of ageing. Finally, we identified new genetic interventions that may ameliorate ageing, by comparing the transcriptomes of 51 mouse mutants not previously associated with ageing to expression signatures of long- and short-lived mice and ageing-related changes.

A clinal polymorphism in the insulin signaling transcription factor foxo contributes to life-history adaptation in Drosophila.

A fundamental aim of adaptation genomics is to identify polymorphisms that underpin variation in fitness traits. In Drosophila melanogaster, latitudinal life-history clines exist on multiple continents and make an excellent system for dissecting the genetics of adaptation. We have previously identified numerous clinal single-nucleotide polymorphism in insulin/insulin-like growth factor signaling (IIS), a pathway known from mutant studies to affect life history. However, the effects of natural variants in this pathway remain poorly understood. Here we investigate how two clinal alternative alleles at foxo, a transcriptional effector of IIS, affect fitness components (viability, size, starvation resistance, fat content). We assessed this polymorphism from the North American cline by reconstituting outbred populations, fixed for either the low- or high-latitude allele, from inbred DGRP lines. Because diet and temperature modulate IIS, we phenotyped alleles across two temperatures (18°C, 25°C) and two diets differing in sugar source and content. Consistent with clinal expectations, the high-latitude allele conferred larger body size and reduced wing loading. Alleles also differed in starvation resistance and expression of insulin-like receptor, a transcriptional target of FOXO. Allelic reaction norms were mostly parallel, with few GxE interactions. Together, our results suggest that variation in IIS makes a major contribution to clinal life-history adaptation.

Evolution of longevity improves immunity in Drosophila.

Much has been learned about the genetics of aging from studies in model organisms, but still little is known about naturally occurring alleles that contribute to variation in longevity. For example, analysis of mutants and transgenes has identified insulin signaling as a major regulator of longevity, yet whether standing variation in this pathway underlies microevolutionary changes in lifespan and correlated fitness traits remains largely unclear. Here, we have analyzed the genomes of a set of Drosophila melanogaster lines that have been maintained under direct selection for postponed reproduction and indirect selection for longevity, relative to unselected control lines, for over 35 years. We identified many candidate loci shaped by selection for longevity and late-life fertility, but - contrary to expectation - we did not find overrepresentation of canonical longevity genes. Instead, we found an enrichment of immunity genes, particularly in the Toll pathway, suggesting that evolutionary changes in immune function might underpin - in part - the evolution of late-life fertility and longevity. To test whether this genomic signature is causative, we performed functional experiments. In contrast to control flies, long-lived flies tended to downregulate the expression of antimicrobial peptides upon infection with age yet survived fungal, bacterial, and viral infections significantly better, consistent with alleviated immunosenescence. To examine whether genes of the Toll pathway directly affect longevity, we employed conditional knockdown using in vivo RNAi. In adults, RNAi against the Toll receptor extended lifespan, whereas silencing the pathway antagonist cactus--causing immune hyperactivation - dramatically shortened lifespan. Together, our results suggest that genetic changes in the age-dependent regulation of immune homeostasis might contribute to the evolution of longer life.

Genomic Evidence for Adaptive Inversion Clines in Drosophila melanogaster.

Clines in chromosomal inversion polymorphisms-presumably driven by climatic gradients-are common but there is surprisingly little evidence for selection acting on them. Here we address this long-standing issue in Drosophila melanogaster by using diagnostic single nucleotide polymorphism (SNP) markers to estimate inversion frequencies from 28 whole-genome Pool-seq samples collected from 10 populations along the North American east coast. Inversions In(3L)P, In(3R)Mo, and In(3R)Payne showed clear latitudinal clines, and for In(2L)t, In(2R)NS, and In(3R)Payne the steepness of the clinal slopes changed between summer and fall. Consistent with an effect of seasonality on inversion frequencies, we detected small but stable seasonal fluctuations of In(2R)NS and In(3R)Payne in a temperate Pennsylvanian population over 4 years. In support of spatially varying selection, we observed that the cline in In(3R)Payne has remained stable for >40 years and that the frequencies of In(2L)t and In(3R)Payne are strongly correlated with climatic factors that vary latitudinally, independent of population structure. To test whether these patterns are adaptive, we compared the amount of genetic differentiation of inversions versus neutral SNPs and found that the clines in In(2L)t and In(3R)Payne are maintained nonneutrally and independent of admixture. We also identified numerous clinal inversion-associated SNPs, many of which exhibit parallel differentiation along the Australian cline and reside in genes known to affect fitness-related traits. Together, our results provide strong evidence that inversion clines are maintained by spatially-and perhaps also temporally-varying selection. We interpret our data in light of current hypotheses about how inversions are established and maintained.

Genome-wide association studies reveal a simple genetic basis of resistance to naturally coevolving viruses in Drosophila melanogaster.

Variation in susceptibility to infectious disease often has a substantial genetic component in animal and plant populations. We have used genome-wide association studies (GWAS) in Drosophila melanogaster to identify the genetic basis of variation in susceptibility to viral infection. We found that there is substantially more genetic variation in susceptibility to two viruses that naturally infect D. melanogaster (DCV and DMelSV) than to two viruses isolated from other insects (FHV and DAffSV). Furthermore, this increased variation is caused by a small number of common polymorphisms that have a major effect on resistance and can individually explain up to 47% of the heritability in disease susceptibility. For two of these polymorphisms, it has previously been shown that they have been driven to a high frequency by natural selection. An advantage of GWAS in Drosophila is that the results can be confirmed experimentally. We verified that a gene called pastrel--which was previously not known to have an antiviral function--is associated with DCV-resistance by knocking down its expression by RNAi. Our data suggest that selection for resistance to infectious disease can increase genetic variation by increasing the frequency of major-effect alleles, and this has resulted in a simple genetic basis to variation in virus resistance.

Genome-wide patterns of latitudinal differentiation among populations of Drosophila melanogaster from North America.

Understanding the genetic underpinnings of adaptive change is a fundamental but largely unresolved problem in evolutionary biology. Drosophila melanogaster, an ancestrally tropical insect that has spread to temperate regions and become cosmopolitan, offers a powerful opportunity for identifying the molecular polymorphisms underlying clinal adaptation. Here, we use genome-wide next-generation sequencing of DNA pools ('pool-seq') from three populations collected along the North American east coast to examine patterns of latitudinal differentiation. Comparing the genomes of these populations is particularly interesting since they exhibit clinal variation in a number of important life history traits. We find extensive latitudinal differentiation, with many of the most strongly differentiated genes involved in major functional pathways such as the insulin/TOR, ecdysone, torso, EGFR, TGFß/BMP, JAK/STAT, immunity and circadian rhythm pathways. We observe particularly strong differentiation on chromosome 3R, especially within the cosmopolitan inversion In(3R)Payne, which contains a large number of clinally varying genes. While much of the differentiation might be driven by clinal differences in the frequency of In(3R)P, we also identify genes that are likely independent of this inversion. Our results provide genome-wide evidence consistent with pervasive spatially variable selection acting on numerous loci and pathways along the well-known North American cline, with many candidates implicated in life history regulation and exhibiting parallel differentiation along the previously investigated Australian cline.

Male-killing Wolbachia do not protect Drosophila bifasciata against viral infection.

BACKGROUND: Insect symbionts employ multiple strategies to enhance their spread through populations, and some play a dual role as both a mutualist and a reproductive manipulator. It has recently been found that this is the case for some strains of Wolbachia, which both cause cytoplasmic incompatibility and protect their hosts against viruses. Here, we carry out the first test as to whether a male-killing strain of Wolbachia also provides a direct benefit to its host by providing antiviral protection to its host Drosophila bifasciata. We infected flies with two positive sense RNA viruses known to replicate in a range of Drosophila species (Drosophila C virus and Flock House virus) and measure the rate of death in Wolbachia positive and negative host lines with the same genetic background. RESULTS: Both viruses caused considerable mortality to D. bifasciata flies, with Drosophila C virus killing 43% more flies than the uninfected controls and Flock House virus killing 78% more flies than the uninfected controls. However, viral induced mortality was unaffected by the presence of Wolbachia. CONCLUSION: In the first male-killing Wolbachia strain tested for antiviral effects, we found no evidence that it conferred protection against two RNA viruses. We show that although antiviral resistance is widespread across the Wolbachia phylogeny, the trait seems to have been lost or gained along some lineages. We discuss the potential mechanisms of this, and can seemingly discount protection against these viruses as a reason why this symbiont has spread through Drosophila populations.