Quantifying the fitness effects of resistance alleles with and without anthelmintic selection pressure using Caenorhabditis elegans.

Albendazole (a benzimidazole) and ivermectin (a macrocyclic lactone) are the two most commonly co-administered anthelmintic drugs in mass-drug administration programs worldwide. Despite emerging resistance, we do not fully understand the mechanisms of resistance to these drugs nor the consequences of delivering them in combination. Albendazole resistance has primarily been attributed to variation in the drug target, a beta-tubulin gene. Ivermectin targets glutamate-gated chloride channels (GluCls), but it is unknown whether GluCl genes are involved in ivermectin resistance in nature. Using Caenorhabditis elegans, we defined the fitness costs associated with loss of the drug target genes singly or in combinations of the genes that encode GluCl subunits. We quantified the loss-of-function effects on three traits: (i) multi-generational competitive fitness, (ii) fecundity, and (iii) development. In competitive fitness and development assays, we found that a deletion of the beta-tubulin gene ben-1 conferred albendazole resistance, but ivermectin resistance required the loss of two GluCl genes (avr-14 and avr-15). The fecundity assays revealed that loss of ben-1 did not provide any fitness benefit in albendazole conditions and that no GluCl deletion mutants were resistant to ivermectin. Next, we searched for evidence of multi-drug resistance across the three traits. Loss of ben-1 did not confer resistance to ivermectin, nor did loss of any single GluCl subunit or combination confer resistance to albendazole. Finally, we assessed the development of 124 C. elegans wild strains across six benzimidazoles and seven macrocyclic lactones to identify evidence of multi-drug resistance between the two drug classes and found a strong phenotypic correlation within a drug class but not across drug classes. Because each gene affects various aspects of nematode physiology, these results suggest that it is necessary to assess multiple fitness traits to evaluate how each gene contributes to anthelmintic resistance.

Quantifying the fitness effects of resistance alleles with and without anthelmintic selection pressure using Caenorhabditis elegans.

Albendazole and ivermectin are the two most commonly co-administered anthelmintic drugs in mass-drug administration programs worldwide. Despite emerging resistance, we do not fully understand the mechanisms of resistance to these drugs nor the consequences of delivering them in combination. Albendazole resistance has primarily been attributed to variation in the drug target, a beta-tubulin gene. Ivermectin targets glutamate-gated chloride channel (GluCl) genes, but it is unknown whether these genes are involved in ivermectin resistance in nature. Using Caenorhabditis elegans, we defined the fitness costs associated with loss of the drug target genes singly or in combinations of the genes that encode GluCl subunits. We quantified the loss-of function effects on three traits: (i) multi-generational competitive fitness, (ii) fecundity, and (iii) development. In competitive fitness and development assays, we found that a deletion of the beta-tubulin gene ben-1 conferred albendazole resistance, but ivermectin resistance required loss of two GluCl genes (avr-14 and avr-15) or loss of three GluCl genes (avr-14, avr-15, and glc-1). The fecundity assays revealed that loss of ben-1 did not provide any fitness benefit in albendazole and that no GluCl deletion mutants were resistant to ivermectin. Next, we searched for evidence of multi-drug resistance across the three traits. Loss of ben-1 did not confer resistance to ivermectin, nor did loss of any single GluCl subunit or combination confer resistance to albendazole. Finally, we assessed the development of 124 C. elegans wild strains across six benzimidazoles and seven macrocyclic lactones to identify evidence of multi-drug resistance between the two drug classes and found a strong phenotypic correlation within a drug class but not across drug classes. Because each gene affects various aspects of nematode physiology, these results suggest that it is necessary to assess multiple fitness traits to evaluate how each gene contributes to anthelmintic resistance.

Praziquantel inhibits Caenorhabditis elegans development and species-wide differences might be cct-8-dependent.

Anthelmintic drugs are used to treat parasitic roundworm and flatworm infections in humans and other animals. Caenorhabditis elegans is an established model to investigate anthelmintics used to treat roundworms. In this study, we use C. elegans to examine the mode of action and the mechanisms of resistance against the flatworm anthelmintic drug praziquantel (PZQ), used to treat trematode and cestode infections. We found that PZQ inhibited development and that this developmental delay varies by genetic background. Interestingly, both enantiomers of PZQ are equally effective against C. elegans, but the right-handed PZQ (R-PZQ) is most effective against schistosome infections. We conducted a genome-wide association mapping with 74 wild C. elegans strains to identify a region on chromosome IV that is correlated with differential PZQ susceptibility. Five candidate genes in this region: cct-8, znf-782, Y104H12D.4, Y104H12D.2, and cox-18, might underlie this variation. The gene cct-8, a subunit of the protein folding complex TRiC, has variation that causes a putative protein coding change (G226V), which is correlated with reduced developmental delay. Gene expression analysis suggests that this variant correlates with slightly increased expression of both cct-8 and hsp-70. Acute exposure to PZQ caused increased expression of hsp-70, indicating that altered TRiC function might be involved in PZQ responses. To test if this variant affects development upon exposure to PZQ, we used CRISPR-Cas9 genome editing to introduce the V226 allele into the N2 genetic background (G226) and the G226 allele into the JU775 genetic background (V226). These experiments revealed that this variant was not sufficient to explain the effects of PZQ on development. Nevertheless, this study shows that C. elegans can be used to study PZQ mode of action and resistance mechanisms. Additionally, we show that the TRiC complex requires further evaluation for PZQ responses in C. elegans.

Variation in anthelmintic responses are driven by genetic differences among diverse C. elegans wild strains.

Treatment of parasitic nematode infections in humans and livestock relies on a limited arsenal of anthelmintic drugs that have historically reduced parasite burdens. However, anthelmintic resistance (AR) is increasing, and little is known about the molecular and genetic causes of resistance for most drugs. The free-living roundworm Caenorhabditis elegans has proven to be a tractable model to understand AR, where studies have led to the identification of molecular targets of all major anthelmintic drug classes. Here, we used genetically diverse C. elegans strains to perform dose-response analyses across 26 anthelmintic drugs that represent the three major anthelmintic drug classes (benzimidazoles, macrocyclic lactones, and nicotinic acetylcholine receptor agonists) in addition to seven other anthelmintic classes. First, we found that C. elegans strains displayed similar anthelmintic responses within drug classes and significant variation across drug classes. Next, we compared the effective concentration estimates to induce a 10% maximal response (EC10) and slope estimates of each dose-response curve of each strain to the laboratory reference strain, which enabled the identification of anthelmintics with population-wide differences to understand how genetics contribute to AR. Because genetically diverse strains displayed differential susceptibilities within and across anthelmintics, we show that C. elegans is a useful model for screening potential nematicides before applications to helminths. Third, we quantified the levels of anthelmintic response variation caused by genetic differences among individuals (heritability) to each drug and observed a significant correlation between exposure closest to the EC10 and the exposure that exhibited the most heritable responses. These results suggest drugs to prioritize in genome-wide association studies, which will enable the identification of AR genes.

Genomic signatures of selection associated with benzimidazole drug treatments in Haemonchus contortus field populations.

Genome-wide methods offer a powerful approach to detect signatures of drug selection. However, limited availability of suitable reference genomes and the difficulty of obtaining field populations with well-defined, distinct drug treatment histories mean there is little information on the signatures of selection in parasitic nematodes and on how best to detect them. This study addresses these knowledge gaps by using field populations of Haemonchus contortus with well-defined benzimidazole treatment histories, leveraging a recently completed chromosomal-scale reference genome assembly. We generated a panel of 49,393 genomic markers to genotype 20 individual adult worms from each of four H. contortus populations: two from closed sheep flocks with an approximate 20 year history of frequent benzimidazole treatment, and two populations with a history of little or no treatment. Sampling occurred in the same geographical region to limit genetic differentiation and maximise the detection sensitivity. A clear signature of selection was detected on chromosome I, centred on the isotype-1 ß-tubulin gene. Two additional, but weaker, signatures of selection were detected; one near the middle of chromosome I spanning 3.75 Mbp and 259 annotated genes, and one on chromosome II spanning a region of 3.3 Mbp and 206 annotated genes, including the isotype-2 ß-tubulin locus. We also assessed how sensitivity was impacted by sequencing depth, worm number, and pooled versus individual worm sequence data. This study provides the first known direct genome-wide evidence for any parasitic nematode, that the isotype-1 ß-tubulin gene is quantitatively the single most important benzimidazole resistance locus. It also identified two additional genomic regions that likely contain benzimidazole resistance loci of secondary importance. This study provides an experimental framework to maximise the power of genome-wide approaches to detect signatures of selection driven by anthelmintic drug treatments in field populations of parasitic nematodes.

Metabarcoding in two isolated populations of wild roe deer (Capreolus capreolus) reveals variation in gastrointestinal nematode community composition between regions and among age classes.

BACKGROUND: Gastrointestinal nematodes are ubiquitous for both domestic and wild ungulates and have varying consequences for health and fitness. They exist as complex communities of multiple co-infecting species, and we have a limited understanding of how these communities vary in different hosts, regions and circumstances or of how this affects their impacts. METHODS: We have undertaken ITS2 rDNA nemabiome metabarcoding with next-generation sequencing on populations of nematode larvae isolated from 149 fecal samples of roe deer of different sex and age classes in the two isolated populations of Chizé and Trois Fontaines in France not co-grazing with any domestic ungulate species. RESULTS: We identified 100 amplified sequence variants (ASVs) that were assigned to 14 gastrointestinal nematode taxa overall at either genus (29%) or species (71%) level. These taxa were dominated by parasites classically found in cervids-e.g. Ostertagia leptospicularis, Spiculopteragia spp. Higher parasite species diversity was present in the Trois Fontaines population than in the Chizé population including the presence of species more typically seen in domestic livestock (Haemonchus contortus, Bunostomum sp., Cooperia punctata, Teladorsagia circumcincta). No differences in parasite species diversity or community composition were seen in the samples collected from three zones of differing habitat quality within the Chizé study area. Young roe deer hosted the highest diversity of gastrointestinal nematodes, with more pronounced effects of age apparent in Trois Fontaines. The effect of host age differed between gastrointestinal nematode species, e.g. there was little effect on O. leptospicularis but a large effect on Trichostrongylus spp. No effect of host sex was detected in either site. CONCLUSIONS: The presence of some livestock parasite species in the Trois Fontaines roe deer population was unexpected given the isolation of this population away from grazing domestic livestock since decades. Overall, our results illustrate the influence of host traits and the local environment on roe deer nemabiome and demonstrate the power of the nemabiome metabarcoding approach to elucidate the composition of gastrointestinal nematode communities in wildlife.

Natural variation in Caenorhabditis elegans responses to the anthelmintic emodepside.

Treatment of parasitic nematode infections depends primarily on the use of anthelmintics. However, this drug arsenal is limited, and resistance against most anthelmintics is widespread. Emodepside is a new anthelmintic drug effective against gastrointestinal and filarial nematodes. Nematodes that are resistant to other anthelmintic drug classes are susceptible to emodepside, indicating that the emodepside mode of action is distinct from previous anthelmintics. The laboratory-adapted Caenorhabditis elegans strain N2 is sensitive to emodepside, and genetic selection and in vitro experiments implicated slo-1, a large K+ conductance (BK) channel gene, in emodepside mode of action. In an effort to understand how natural populations will respond to emodepside, we measured brood sizes and developmental rates of wild C. elegans strains after exposure to the drug and found natural variation across the species. Some of the observed variation in C. elegans emodepside responses correlates with amino acid substitutions in slo-1, but genetic mechanisms other than slo-1 coding variants likely underlie emodepside resistance in wild C. elegans strains. Additionally, the assayed strains have higher offspring production in low concentrations of emodepside (a hormetic effect). We find that natural variation affects emodepside sensitivity, supporting the suitability of C. elegans as a model system to study emodepside responses across natural nematode populations.

Two novel loci underlie natural differences in Caenorhabditis elegans abamectin responses.

Parasitic nematodes cause a massive worldwide burden on human health along with a loss of livestock and agriculture productivity. Anthelmintics have been widely successful in treating parasitic nematodes. However, resistance is increasing, and little is known about the molecular and genetic causes of resistance for most of these drugs. The free-living roundworm Caenorhabditis elegans provides a tractable model to identify genes that underlie resistance. Unlike parasitic nematodes, C. elegans is easy to maintain in the laboratory, has a complete and well annotated genome, and has many genetic tools. Using a combination of wild isolates and a panel of recombinant inbred lines constructed from crosses of two genetically and phenotypically divergent strains, we identified three genomic regions on chromosome V that underlie natural differences in response to the macrocyclic lactone (ML) abamectin. One locus was identified previously and encodes an alpha subunit of a glutamate-gated chloride channel (glc-1). Here, we validate and narrow two novel loci using near-isogenic lines. Additionally, we generate a list of prioritized candidate genes identified in C. elegans and in the parasite Haemonchus contortus by comparison of ML resistance loci. These genes could represent previously unidentified resistance genes shared across nematode species and should be evaluated in the future. Our work highlights the advantages of using C. elegans as a model to better understand ML resistance in parasitic nematodes.

Complementary Approaches with Free-living and Parasitic Nematodes to Understanding Anthelmintic Resistance.

Anthelmintic drugs are the major line of defense against parasitic nematode infections, but the arsenal is limited and resistance threatens sustained efficacy of the available drugs. Discoveries of the modes of action of these drugs and mechanisms of resistance have predominantly come from studies of a related nonparasitic nematode species, Caenorhabditis elegans, and the parasitic nematode Haemonchus contortus. Here, we discuss how our understanding of anthelmintic resistance and modes of action came from the interplay of results from each of these species. We argue that this 'cycle of discovery', where results from one species inform the design of experiments in the other, can use the complementary strengths of both to understand anthelmintic modes of action and mechanisms of resistance.

Genomic and transcriptomic variation defines the chromosome-scale assembly of Haemonchus contortus, a model gastrointestinal worm.

Haemonchus contortus is a globally distributed and economically important gastrointestinal pathogen of small ruminants and has become a key nematode model for studying anthelmintic resistance and other parasite-specific traits among a wider group of parasites including major human pathogens. Here, we report using PacBio long-read and OpGen and 10X Genomics long-molecule methods to generate a highly contiguous 283.4 Mbp chromosome-scale genome assembly including a resolved sex chromosome for the MHco3(ISE).N1 isolate. We show a remarkable pattern of conservation of chromosome content with Caenorhabditis elegans, but almost no conservation of gene order. Short and long-read transcriptome sequencing allowed us to define coordinated transcriptional regulation throughout the parasite's life cycle and refine our understanding of cis- and trans-splicing. Finally, we provide a comprehensive picture of chromosome-wide genetic diversity both within a single isolate and globally. These data provide a high-quality comparison for understanding the evolution and genomics of Caenorhabditis and other nematodes and extend the experimental tractability of this model parasitic nematode in understanding helminth biology, drug discovery and vaccine development, as well as important adaptive traits such as drug resistance.

High levels of third-stage larvae (L3) overwinter survival for multiple cattle gastrointestinal nematode species on western Canadian pastures as revealed by ITS2 rDNA metabarcoding.

BACKGROUND: The ability of infective larvae of cattle gastrointestinal nematode (GIN) species to overwinter on pastures in northerly climatic zones with very cold dry winters is poorly understood. This is an important knowledge gap with critical implications for parasite risk assessment and control. METHODS: Infective third-stage larvae (L3) were quantified in samples of fecal pats, together with adjacent grass and soil, before and after winter on three farms in southern, central and northern Alberta. Nemabiome ITS2 metabarcoding was then performed on the harvested L3 populations to determine the species composition. Finally, parasite-free tracer calves were used to investigate if the L3 surviving the winter could infect calves and develop to adult worms in spring. RESULTS: Farm level monitoring, using solar powered weather stations, revealed that ground temperatures were consistently higher, and less variable, than the air temperatures; minimum winter air and ground temperatures were - 32.5 °C and - 24.7 °C respectively. In spite of the extremely low minimum temperatures reached, L3 were recovered from fecal pats and grass before and after winter with only a 38% and 61% overall reduction over the winter, respectively. Nemabiome ITS2 metabarcoding assay revealed that the proportion of L3 surviving the winter was high for both Cooperia oncophora and Ostertagia ostertagi although survival of the former species was statistically significantly higher than the latter. Nematodirus helvetinaus and Trichostrongylus axei could be detected after winter whereas Haemonchus placei L3 could not overwinter at all. Adult C. oncophora, O. ostertagi and N. helvetianus could be recovered from tracer calves grazing after the winter. CONCLUSIONS: The largest proportion of L3 were recovered from fecal pats suggesting this is important refuge for L3 survival. Results also show that L3 of several GIN parasite species can survive relatively efficiently on pastures even in the extreme winter conditions in western Canada. Tracer calf experiments confirmed that overwintered L3 of both C. oncophora and O. ostertagi were capable of establishing a patent infection in the following spring. These results have important implications for the epidemiology, risk of production impact and the design of effective control strategies. The work also illustrates the value of applying ITS2 nemabiome metabarcoding to environmental samples.

A database for ITS2 sequences from nematodes.

BACKGROUND: Marker gene surveys have a wide variety of applications in species identification, population genetics, and molecular epidemiology. As these methods expand to new types of organisms and additional markers beyond 16S and 18S rRNA genes, comprehensive databases are a critical requirement for proper analysis of these data. RESULTS: Here we present an ITS2 rDNA database for marker gene surveys of both free-living and parasitic nematode populations and the software used to build the database. This is currently the most complete and up-to-date ITS2 database for nematodes and is able to reproduce previous analysis that used a smaller database. CONCLUSIONS: This database is an important resource for researchers working on nematodes and also provides a tool to create ITS2 databases for any given taxonomy.

Wild ruminants as reservoirs of domestic livestock gastrointestinal nematodes.

Gastrointestinal nematode (GIN) infections in cattle cause appetite suppression which leads to poor feed conversion, reduced weight gain and reduced milk production. Overuse and exclusive reliance on anthelmintic drugs has resulted in widespread resistance in many parasitic nematode species infecting livestock making control increasingly difficult. Wild ruminants are competent hosts of a number of nematode species that typically infect and are best adapted for cattle, sheep, and goats. Thus, the potential exists for wild ruminants to act as reservoirs in the translocation of domestic GIN, including those carrying anthelmintic resistance mutations as well as susceptible genotypes. The potential for parasite exchange is heightened by interfaces or ecotones between managed and wild rangelands, and by perturbations linked to climate warming that can increasingly alter the distributions of wild ungulates and their interactions with domestic and free-ranging ruminants. To investigate the extent to which wild ruminants harbour parasites capable of infecting domestic ruminants we first performed an epidemiological study of feces from wildlife hosts that spanned 16 states and included white-tailed deer (85 % of the samples), pronghorn, elk, mule deer, bighorn sheep, moose, cattle, and caribou across the United States. All samples were cultured to third stage larvae and nematode DNA was isolated and PCR amplified. Among the 548 wild ruminant samples received, 33 % (181 samples) were positive for nematode DNA, among which half (84 samples) contained DNA from GIN species commonly found in cattle. DNA from cattle GIN species was detected in 46 % of samples from the Northeast, 42 % from the Southeast, 10 % from the Midwest, 0 % from the Southwest and 11 % from the West. Deep amplicon sequencing of the ITS-2 rDNA indicated that Ostertagia and Trichostrongylus were present in 90 % and 69 % of the nematode DNA positive samples, respectively, whereas Haemonchus, Cooperia and Oesophagostomum were present in 26 %, 2 % and 10 % of the samples, respectively. These data clearly show that wild ruminants commonly harbour multiple parasite species whose primary hosts are domestic cattle, and suggest that further work is warranted to investigate their specific roles in the management of anthelmintic resistance.

Benchmarking hybrid assemblies of Giardia and prediction of widespread intra-isolate structural variation.

BACKGROUND: Currently available short read genome assemblies of the tetraploid protozoan parasite Giardia intestinalis are highly fragmented, highlighting the need for improved genome assemblies at a reasonable cost. Long nanopore reads are well suited to resolve repetitive genomic regions resulting in better quality assemblies of eukaryotic genomes. Subsequent addition of highly accurate short reads to long-read assemblies further improves assembly quality. Using this hybrid approach, we assembled genomes for three Giardia isolates, two with published assemblies and one novel, to evaluate the improvement in genome quality gained from long reads. We then used the long reads to predict structural variants to examine this previously unexplored source of genetic variation in Giardia. METHODS: With MinION reads for each isolate, we assembled genomes using several assemblers specializing in long reads. Assembly metrics, gene finding, and whole genome alignments to the reference genomes enabled direct comparison to evaluate the performance of the nanopore reads. Further improvements from adding Illumina reads to the long-read assemblies were evaluated using gene finding. Structural variants were predicted from alignments of the long reads to the best hybrid genome for each isolate and enrichment of key genes was analyzed using random genome sampling and calculation of percentiles to find thresholds of significance. RESULTS: Our hybrid assembly method generated reference quality genomes for each isolate. Consistent with previous findings based on SNPs, examination of heterozygosity using the structural variants found that Giardia BGS was considerably more heterozygous than the other isolates that are from Assemblage A. Further, each isolate was shown to contain structural variant regions enriched for variant-specific surface proteins, a key class of virulence factor in Giardia. CONCLUSIONS: The ability to generate reference quality genomes from a single MinION run and a multiplexed MiSeq run enables future large-scale comparative genomic studies within the genus Giardia. Further, prediction of structural variants from long reads allows for more in-depth analyses of major sources of genetic variation within and between Giardia isolates that could have effects on both pathogenicity and host range.

Deep amplicon sequencing as a powerful new tool to screen for sequence polymorphisms associated with anthelmintic resistance in parasitic nematode populations.

Parasitic gastrointestinal nematodes contribute to significant human morbidity and cause billions of dollars per year in lost agricultural production. Control is dependent on the use of anthelmintic drugs which, in the case of livestock parasites, is severely compromised by the widespread development of drug resistance. There are now concerns regarding the emergence of anthelmintic resistance in parasitic nematodes of humans in response to the selection pressure resulting from mass drug administration programs. Consequently, there is an urgent need for sensitive, scalable and accurate diagnostic tools to detect the emergence of anthelmintic resistance. Detecting and measuring the frequency of resistance-associated mutations in parasite populations has the potential to provide sensitive and quantitative assessment of resistance emergence from an early stage. Here, we describe the development and validation of deep amplicon sequencing as a powerful new approach to detect and quantify the frequency of single nucleotide polymorphisms associated with benzimidazole resistance. We have used parasite communities in sheep to undertake a proof-of-concept study of this approach. Sheep provide an excellent host system, as there are multiple co-infecting trichostrongylid nematode species, each likely with a varying prevalence of benzimidazole resistance. We demonstrate that the approach provides an accurate measure of resistance allele frequencies, and can reliably detect resistance alleles down to a frequency of 0.1%, making it particularly valuable for screening mutations in the early stages of resistance. We illustrate the power of the technique by screening UK sheep flocks for benzimidazole resistance-associated single nucleotide polymorphisms at three different codons of the ß-tubulin gene in seven different parasite species from 164 populations (95 from ewes and 69 from lambs) in a single MiSeq sequencing run. This approach provides a powerful new tool to screen for the emergence of anthelmintic resistance mutations in parasitic nematode populations of both animals and humans.

Resequencing Helminth Genomes for Population and Genetic Studies.

Next-generation sequencing has become increasingly accessible and economical, making genome-wide studies routine for many species, including humans, model organisms, and domestic livestock. However, in the case of helminth parasites, there are still major practical challenges to the application of these approaches for genetic and population studies. Dozens to hundreds of individual parasites from multiple populations may need to be re-sequenced which, together with the relatively large size of helminth genomes, can still make whole-genome resequencing of individual parasites unfeasible for many studies. Fortunately, there are alternative approaches to the complete sequencing of genomes when conducting genome-wide studies. Here we review some strategies, including genome subsampling and pooling, that enable genome-wide analysis of large numbers of parasites in populations.

Parallel Gene Expression Differences between Low and High Latitude Populations of Drosophila melanogaster and D. simulans.

Gene expression variation within species is relatively common, however, the role of natural selection in the maintenance of this variation is poorly understood. Here we investigate low and high latitude populations of Drosophila melanogaster and its sister species, D. simulans, to determine whether the two species show similar patterns of population differentiation, consistent with a role for spatially varying selection in maintaining gene expression variation. We compared at two temperatures the whole male transcriptome of D. melanogaster and D. simulans sampled from Panama City (Panama) and Maine (USA). We observed a significant excess of genes exhibiting differential expression in both species, consistent with parallel adaptation to heterogeneous environments. Moreover, the majority of genes showing parallel expression differentiation showed the same direction of differential expression in the two species and the magnitudes of expression differences between high and low latitude populations were correlated across species, further bolstering the conclusion that parallelism for expression phenotypes results from spatially varying selection. However, the species also exhibited important differences in expression phenotypes. For example, the genomic extent of genotype × environment interaction was much more common in D. melanogaster. Highly differentiated SNPs between low and high latitudes were enriched in the 3' UTRs and CDS of the geographically differently expressed genes in both species, consistent with an important role for cis-acting variants in driving local adaptation for expression-related phenotypes.

Laboratory selection for increased longevity in Drosophila melanogaster reduces field performance.

Drosophila melanogaster is frequently used in ageing studies to elucidate which mechanisms determine the onset and progress of senescence. Lines selected for increased longevity have often been shown to perform as well as or superior to control lines in life history, stress resistance and behavioural traits when tested in the laboratory. Functional senescence in longevity selected lines has also been shown to occur at a slower rate. However, it is known that performance in a controlled laboratory setting is not necessarily representative of performance in nature. In this study the effect of ageing, environmental temperature and longevity selection on performance in the field was tested. Flies from longevity selected and control lines of different ages (2, 5, 10 and 15 days) were released in an environment free of natural food sources. Control flies were tested at low, intermediate and high temperatures, while longevity selected flies were tested at the intermediate temperature only. The ability of flies to locate and reach a food source was tested. Flies of intermediate age were generally better at locating resources than both younger and older flies, where hot and cold environments accelerate the senescent decline in performance. Control lines were better able to locate a resource compared to longevity selected lines of the same age, suggesting that longevity comes at a cost in early life field fitness, supporting the antagonistic pleiotropy theory of ageing.

Longevity for free? Increased reproduction with limited trade-offs in Drosophila melanogaster selected for increased life span.

Selection for increased life span in Drosophila melanogaster has been shown to correlate with decreased early fecundity and increased fecundity later in life. This phenomenon has been ascribed to the existence of trade-offs in which limited resources can be invested in either somatic maintenance or reproduction. In our longevity selection lines, we did not find such a trade-off. Rather, we find that females have similar or higher fecundity throughout life compared to non-selected controls. To determine whether increased longevity affects responses in other traits, we looked at several stress resistance traits (chill coma recovery, heat knockdown, desiccation and starvation), geotactic behaviour, egg-to-adult viability, body size, developmental time as well as metabolic rate. Longevity selected flies were more starvation resistant. However, in females longevity and fecundity were not negatively correlated with the other traits assayed. Males from longevity selected lines were slower at recovering from a chill induced coma and resting metabolic rate increased with age, but did not correlate with life span.

Confirming candidate genes for longevity by RT-qPCR using two different genetic backgrounds and selection methods.

Elucidating genes that affect life span or that can be used as biomarkers for ageing has received attention in diverse studies in recent years. Using model organisms and various approaches several genes have been linked to the longevity phenotype. For Drosophila melanogaster those studies have usually focussed on one sex and on flies originating from one genetic background, and results from different studies often do not overlap. Using D. melanogaster selected for increased longevity we aimed to find robust longevity related genes by examining gene expression in both sexes of flies originating from different genetic backgrounds. Further, we compared expression changes across three ages, when flies were young, middle aged or old, to examine how candidate gene expression changes with the onset of ageing. We selected 10 genes based on their expression differences in prior microarray studies. For about 50% of these we confirmed their potential as a candidate longevity gene. We found one robust candidate gene for longevity, CG32638. Three other genes, CG8934, mRpS10 and Spn43Ad, showed a tendency to be involved in life span determination in both backgrounds tested.