Origin matters: Using a local reference genome improves measures in population genomics.

Genome sequencing enables answering fundamental questions about the genetic basis of adaptation, population structure and epigenetic mechanisms. Yet, we usually need a suitable reference genome for mapping population-level resequencing data. In some model systems, multiple reference genomes are available, giving the challenging task of determining which reference genome best suits the data. Here, we compared the use of two different reference genomes for the three-spined stickleback (Gasterosteus aculeatus), one novel genome derived from a European gynogenetic individual and the published reference genome of a North American individual. Specifically, we investigated the impact of using a local reference versus one generated from a distinct lineage on several common population genomics analyses. Through mapping genome resequencing data of 60 sticklebacks from across Europe and North America, we demonstrate that genetic distance among samples and the reference genomes impacts downstream analyses. Using a local reference genome increased mapping efficiency and genotyping accuracy, effectively retaining more and better data. Despite comparable distributions of the metrics generated across the genome using SNP data (i.e. π, Tajima's D and FST ), window-based statistics using different references resulted in different outlier genes and enriched gene functions. A marker-based analysis of DNA methylation distributions had a comparably high overlap in outlier genes and functions, yet with distinct differences depending on the reference genome. Overall, our results highlight how using a local reference genome decreases reference bias to increase confidence in downstream analyses of the data. Such results have significant implications in all reference-genome-based population genomic analyses.

In vitro leukocyte response of three-spined sticklebacks (Gasterosteus aculeatus) to helminth parasite antigens.

Helminth parasites of teleost fish have evolved strategies to evade and manipulate the immune responses of their hosts. Responsiveness of fish host immunity to helminth antigens may therefore vary depending on the degree of host-parasite counter-adaptation. Generalist parasites, infective for a number of host species, might be unable to adapt optimally to the immune system of a certain host species, while specialist parasites might display high levels of adaptation to a particular host species. The degree of adaptations may further differ between sympatric and allopatric host-parasite combinations. Here, we test these hypotheses by in vitro exposure of head kidney leukocytes from three-spined sticklebacks (Gasterosteus aculeatus) to antigens from parasites with a broad fish host range (Diplostomum pseudospathaceum, Triaenophorus nodulosus), a specific fish parasite of cyprinids (Ligula intestinalis) and parasites highly specific only to a single fish species as second intermediate host (Schistocephalus pungitii, which does not infect G. aculeatus, and Schistocephalus solidus, infecting G. aculeatus). In vitro responses of stickleback leukocytes to S. solidus antigens from six European populations, with S. solidus prevalence from <1% to 66% were tested in a fully crossed experimental design. Leukocyte cultures were analysed by means of flow cytometry and a chemiluminescence assay to quantify respiratory burst activity. We detected decreasing magnitudes of in vitro responses to antigens from generalist to specialist parasites and among specialists, from parasites that do not infect G. aculeatus to a G. aculeatus-infecting species. Generalist parasites seem to maintain their ability to infect different host species at the costs of relatively higher immunogenicity compared to specialist parasites. In a comparison of sympatric and allopatric combinations of stickleback leukocytes and antigens from S. solidus, magnitudes of in vitro responses were dependent on the prevalence of the parasite in the population of origin, rather than on sympatry. Antigens from Norwegian (prevalence 30-50%) and Spanish (40-66%) S. solidus induced generally higher in vitro responses compared to S. solidus from two German (<1%) populations. Likewise, leukocytes from stickleback populations with a high S. solidus prevalence showed higher in vitro responses to S. solidus antigens compared to populations with low S. solidus prevalence. This suggests a rather low degree of local adaptation in S. solidus populations, which might be due to high gene flow among populations because of their extremely mobile final hosts, fish-eating birds.

Induction of diploid gynogenesis in an evolutionary model organism, the three-spined stickleback (Gasterosteus aculeatus).

BACKGROUND: Rapid advances in genomics have provided nearly complete genome sequences for many different species. However, no matter how the sequencing technology has improved, natural genetic polymorphism complicates the production of high quality reference genomes. To address this problem, researchers have tried using artificial modes of genome manipulation such as gynogenesis for fast production of inbred lines. RESULTS: Here, we present the first successful induction of diploid gynogenesis in an evolutionary model system, the three-spined sticklebacks (Gasterosteus aculeatus), using a combination of UV-irradiation of the sperm and heat shock (HS) of the resulting embryo to inhibit the second meiotic division. Optimal UV irradiation of the sperm was established by exposing stickleback sperm to a UV- light source at various times. Heat shock parameters like temperature, duration, and time of initiation were tested by subjecting eggs fertilized with UV inactivated sperm 5, 10, 15, 20, 25, or 30 minutes post fertilization (mpf) to 30°C, 34°C, or 38°C for 2, 4, 6 or 8 minutes. Gynogen yield was highest when stickleback eggs were activated with 2 minutes UV-irradiated sperm and received HS 5 mpf at 34°C for 4 minutes. CONCLUSIONS: Diploid gynogenesis has been successfully performed in three-spined stickleback. This has been confirmed by microsatellite DNA analysis which revealed exclusively maternal inheritance in all gynogenetic fry tested. Ploidy verification by flow cytometry showed that gynogenetic embryos/larvae exhibiting abnormalities were haploids and those that developed normally were diploids, i.e., double haploids that can be raised until adult size.