Clonal Structure of Salmon ParasiteGyrodactylus salarison a Coevolutionary Gradient on Fennoscandian Salmon (Salmo salar)

Barriers to chimpanzee gene flow at the south-east edge of their distribution

Populations on the edge of a species' distribution may represent an important source of adaptive diversity, yet these populations tend to be more fragmented and are more likely to be geographically isolated. Lack of genetic exchanges between such populations, due to barriers to animal movement, can not only compromise adaptive potential but also lead to the fixation of deleterious alleles. The south-eastern edge of chimpanzee distribution is particularly fragmented, and conflicting hypotheses have been proposed about population connectivity and viability. To address this uncertainty, we generated both mitochondrial and MiSeq-based microsatellite genotypes for 290 individuals ranging across western Tanzania. While shared mitochondrial haplotypes confirmed historical gene flow, our microsatellite analyses revealed two distinct clusters, suggesting two populations currently isolated from one another. However, we found evidence of high levels of gene flow maintained within each of these clusters, one of which covers an 18,000 km2 ecosystem. Landscape genetic analyses confirmed the presence of barriers to gene flow with rivers and bare habitats highly restricting chimpanzee movement. Our study demonstrates how advances in sequencing technologies, combined with the development of landscape genetics approaches, can resolve ambiguities in the genetic history of critical populations and better inform conservation efforts of endangered species.

Geographical distribution of Gyrodactylus salaris Malmberg, 1957 (Monogenea, Gyrodactylidae)

BACKGROUND

Gyrodactylus salaris Malmberg, 1957 is an OIE (Office International des Epizooties)-listed parasitic pathogen and had until the current study been reported from 19 countries across Europe, although many of these records require confirmation. The last comprehensive evaluation regarding the distribution of G. salaris, however, was made in 2007, although some of the states identified as being G. salaris-positive were ascribed this status based on misidentifications, on partial data resulting from either morphological or molecular tests, or from records that have not been revisited since their early reporting. It is thus important to go through the reports on G. salaris to obtain a status for each country.

METHODS

To provide a revised update of the G. salaris distribution, a literature review was necessary. This literature, however, was not always readily accessible and, in certain cases, the article only made superficial reference to the parasite without providing details or data to support the identification. In most cases, the original specimens were not deposited in a national collection. Additional Gyrodactylus material for the current study was obtained from selected salmonid populations with the aim to contribute to current understanding regarding the distribution of G. salaris. Additional parasite material collected for this study was processed following standard procedures for species identification in Gyrodactylus [1].

RESULTS

From the work conducted in the current study, G. salaris is reported from a further three regions in Italy, alongside three other species, and appears to occur extensively throughout central Italy without causing significant mortalities to its rainbow trout, Oncorhynchus mykiss (Walbaum), host. The analysis of archive material from G. salaris-positive farms would suggest that G. salaris has been in this country since at least 2000. Material obtained from rainbow trout from Finland and Germany are confirmed as G. salaris, supporting existing data for these countries. No specimens of G. salaris, however, were found in the additional Gyrodactylus material obtained from rainbow trout reared in Portugal and Spain. A morphologically similar species, Gyrodactylus teuchis Lautraite, Blanc, Thiery, Daniel et Vigneulle, 1999, however, was found.

CONCLUSIONS

Following the present review, Gyrodactylus salaris is reported from 23 out of 50 recognised states throughout Europe; only records from 14 of these states have been confirmed by either morphology and/or by an appropriate molecular test and are considered valid, while only nine of these records have been confirmed by a combination of both methods.

A large wild salmon stock shows genetic and life history differentiation within, but not between, rivers

Anadromous salmonid fishes frequently exhibit strong geographic population structuring. However, population genetic differentiation of Atlantic salmon (Salmo salar) at fine geographic scales differs across equivalent spatial extents in different regions. So far, fine-scale genetic differentiation has not been assessed in rivers of the Baltic Sea, a region that contains an evolutionarily distinct Atlantic salmon lineage. Thus, Baltic salmon are currently managed on the river level, without focus on potential genetic structure and diversity within rivers. Here, we used microsatellites to characterize the genetic structure of wild juvenile salmon sampled throughout the interconnected, northern Baltic Tornio and Kalix Rivers. We found genetic differentiation within the two rivers, but not between them: salmon in the upper reaches differed from individuals in the lower reaches, regardless of river system. Further, examining smolts migrating from the river to the sea and adults returning from the sea to spawn, we found an association between the genetic structure and seasonal migration timing. Out-migrating smolts genetically assigned to upper river reaches were older and tended to reach the sea later in the season than smolts from the lower reaches. In contrast, mature adults originating from the upper reaches returned to the river early in the season. Our observation of genetic population structuring between downstream and upstream reaches of the large Tornio and Kalix rivers, and its association with migration timing, implies that careful temporal management of the northern Baltic fisheries would help to preserve the diversity and sustainability of the wild salmon stocks of these rivers.

Revision of Gyrodactylus salaris phylogeny inspired by new evidence for Eemian crossing between lineages living on grayling in Baltic and White sea basins

In this research, grayling-specific Gyrodactylus salarisMalmberg, 1957 isolates from Baltic Sea basin were collected in Sweden for the first time. Samples were obtained in three drainage systems: Kalixälven (River Kaitum), Ljungan (River Sölvbacka strömmar), and Umeälven (River Juktån). Three molecular markers were analysed: nuclear ITS rDNA (Internal Transcribed Spacer) and ADNAM1 (Anonymous DNA Marker 1), and mitochondrial cox1 gene. As a result, four new mitochondrial haplotypes were identified (III-C1tt, III-C1tt_ht, IX-A1tt and X-A1tt). The ADNAM1 analyses resulted in revealing two new alleles (WS4 and BS9) and two new genotypes (T6 and T7). T7 seems to be an indicator of ancient crossing between Baltic and White Sea lineages of the parasite which happened during a first 3000-year period of Eemian interglacial about 130,000 years ago in the connection between Baltic and White Sea. Molecular clock estimates were adjusted, revealing the mean substitution rate and the divergence rate among branches of 3.6% (95% HPD: 2.2%–5.2%) and 7.2% per million years, respectively. As a result, cox1 phylogeny rooted with the introgressed haplotypes has been revised and altered in accordance to new data, revealing fourteen equidistant lineages five of which have been excluded from the study. Based on the new phylogenetic approach, including the molecular clock, this work suggests an overall revision of G. salaris phylogeny and attempts at precisely drawing the division of lineages within this polytypic species as well as proposes unification in nomenclature for its strains.

Genetic gradient of a host-parasite pair along a river persisted ten years against physical mobility: Baltic Salmo salar vs. Gyrodactylus salaris

The Atlantic salmon, Salmo salar L., in the Tornio River in the Northern Baltic Sea basin accommodates a monogenean ectoparasite, Gyrodactylus salaris. The aim of the study was to understand the population structure of apparently co-adapted host-parasite system: no parasite-associated mortality has been reported. The parasite burden among salmon juveniles (parr) was monitored along 460km of the river in 2000-2009. Among the parr, 33.0% were infected (n_fish=1913). The genetic structure of the parasite population was studied by sequencing an anonymous nuclear DNA marker (ADNAM1, three main genotypes) and mitochondrial CO1 (three clades, six haplotypes). During the ten years, the parasite population was strongly and stably genetically differentiated among up- and downstream nurseries (n_ADNAM1=411, F_ST=0.579; n_CO1=443, F_ST=0.534). Infection prevalence among the smolts migrating to sea was higher than in the sedentary parr populations (82.2%, n_fish=129). The spatial differentiation observed among the sedentary juveniles was reflected temporally in the smolt run: parasite genotypes dominating the upper part of the river arrived later than downstream dwellers (medians June 4 and June 2) to the trap 7km from the river mouth. The nuclear and mitochondrial markers were in stable disequilibrium which was not relaxed in the contact zone or among the smolts where the parasite clones often met on individual fish. Only five parasite specimens on smolts (n_worms=217) were putative recent sexual recombinants. The contribution of extant salmon hatcheries into the infection was negligible. The host salmon population in Tornio River is known to show significant spatial differentiation (F_ST=0.022). The stable spatial genetic structure of the parasite against the high physical mobility suggested a possibility of local co-adaptation of the host-parasite subpopulations.

Displaced phylogeographic signals from Gyrodactylus arcuatus, a parasite of the three-spined stickleback Gasterosteus aculeatus, suggest freshwater glacial refugia in Europe

We examined the global mitochondrial phylogeography of Gyrodactylus arcuatus, a flatworm ectoparasite of three-spined stickleback Gasterosteus aculeatus. In accordance with the suggested high divergence rate of 13%/million years, the genetic variation of the parasite was high: haplotype diversity h=0.985 and nucleotide diversity π=0.0161. The differentiation among the parasite populations was substantial (Φst=0.759), with two main allopatric clades (here termed Euro and North) accounting for 54% of the total genetic variation. The diversity center of the Euro clade was in the Baltic Sea, while the North clade was spread across the Barents and White Seas. A single haplotype within the North clade was found in the western and eastern Pacific Ocean. Divergence of main clades was estimated to be circa 200 thousand years ago. Each main clade was further divided into six distinct subclades, estimated to have diverged in isolation since 135 thousand years ago. This second division corresponds approximately to the Eemian interglacial predating the last glacial maximum. A demographic expansion of the subclades is associated with colonisation of northern Europe since the last glacial maximum, circa 15-40 thousand years ago. The parasite phylogeny is most likely explained by sequential isolated bottlenecks and expansions in numerous allopatric refugia. The postglacial intermingling and high variation in the marine parasite populations, separately in the Baltic and Barents Seas, suggest low competition of divergent parasite matrilines, coupled with a large population size and high rate of dispersal of hosts. The genetic contribution of the assumed refugial fish populations maintaining the parasite during the last glacial maximum was not detected among the marine sticklebacks, which perhaps were infected after range expansion.

Occurrence and morphogenetic characteristics of Gyrodactylus (Monogenea: Gyrodactylidae) from a rainbow trout farm (Lake Ladoga, Russia)

Gyrodactylus parasite infected juveniles on rainbow trout Oncorhynchus mykiss (Walbaum, 1792) from a fish farm in Lake Ladoga were investigated. The observed cases of infection in fish featured a high prevalence, when almost all of the fish were infected. However, if an outbreak of the monogenean infection is observed in spring, the intensity of the infection may be low, and when the infection occurs in the ice-covered period (late autumn - winter), the number of parasites on the fins of a single fish may exceed 3000 specimens. Molecular identification of the parasite demonstrated that the infecting clone was identical with rainbow trout specific strain of Gyrodactylus salaris RBT widely spread in Northern Europe, but a small proportion of the parasites were the hybrid clone Gyrodactylus pomeraniae x G. lavareti. Morphological variations of hooks and other opisthaptor parts in the monogenean Gyrodactylus depending on the intensity of infection in rainbow trout were demonstrated.

Genetic population structure of Gyrodactylus thymalli (Monogenea) in a large Norwegian river system

The extent of geographic genetic variation is the result of several processes such as mutation, gene flow, selection and drift. Processes that structure the populations of parasite species are often directly linked to the processes that influence the host. Here, we investigate the genetic population structure of the ectoparasite Gyrodactylus thymalli Žitňan, 1960 (Monogenea) collected from grayling (Thymallus thymallus L.) throughout the river Glomma, the largest watercourse in Norway. Parts of the mitochondrial dehydrogenase subunit 5 (NADH 5) and cytochrome oxidase I (COI) genes from 309 G. thymalli were analysed to study the genetic variation and investigated the geographical distribution of parasite haplotypes. Three main clusters of haplotypes dominated the three distinct geographic parts of the river system; one cluster dominated in the western main stem of the river, one in the eastern and one in the lower part. There was a positive correlation between pairwise genetic distance and hydrographic distance. The results indicate restricted gene flow between sub-populations of G. thymalli, most likely due to barriers that limit upstream migration of infected grayling. More than 80% of the populations had private haplotypes, also indicating long-time isolation of sub-populations. According to a molecular clock calibration, much of the haplotype diversity of G. thymalli in the river Glomma has developed after the last glaciation.

Coevolution between Contracaecum (Nematoda, Anisakidae) and Austrolebias (Cyprinodontiformes, Rivulidae) host-parasite complex from SW Atlantic coastal basins

In recent years, molecular studies in host-parasite interactions in terms of coevolution have become important. Larvae (L3) of two species of Contracaecum were found parasitizing species of Rivulidae in the Atlantic coastal basins from Uruguay. The aim of this study is to determine the patterns of differentiation of this host-parasite complex in order to clarify possible coevolutionary events in such interaction throughout phylogeographic approach using both nuclear and mitochondrial molecular markers (internal transcribed spacers (ITS) and cytochrome oxidase subunit 1 (cox-1)). Based on both markers, intraspecific variation in Contracaecum species was lower than 2 %, while interspecific variation was greater than 10 %. Both species of Contracaecum constitute monophyletic groups. Contracaecum resulted in a paraphyletic genus when incorporating other Contracaecum species and closely related nematode sequences from GenBank. ITS regions showed that Contracaecum sp. 1 is more closely related to other species of the same genus than with their counterparts from Atlantic coastal basins in Uruguay. Haplotype network for both markers corroborate the existence of two distinct taxa. While ITS pairwise FST comparisons and the indirect estimate of gene flow confirm the existence of two distinct Contracaecum species, mitochondrial gene detected low levels of migrants between some of the populations from both species. Our results suggest that coevolution in this host-parasite complex species is plausible. Parasite cladogenetic events occur almost simultaneously with the separation of the hypothetical ancestors of each species complex of Austrolebias during Pliocene. Additionally, the two lineages of Contracaecum colonize differently the species within each of the Austrolebias complexes.

MicroRNA loci support conspecificity of Gyrodactylus salaris and Gyrodactylus thymalli (Platyhelminthes: Monogenea)

The monogenean flatworm Gyrodactylus salaris is a serious threat to wild and farmed Atlantic salmon stocks in Norway. Morphologically, the closely related but harmless Gyrodactylus thymalli on grayling can hardly be distinguished from G. salaris. Until now, molecular approaches could not resolve unambiguously whether G. salaris and G. thymalli represent just one polytypic species, two polytypic species or a complex of more than two species. In the first known genome-wide analysis utilizing 37 conserved microRNA loci, the genetic differentiation of seven populations of G. salaris and G. thymalli was assessed. The concatenated alignment spanned 21,742bp including 62 variable positions. A neighbor-joining cluster analysis did not support any host-based or mitochondrial haplotype-based grouping of strains. We conclude that a two species concept for G. salaris and G. thymalli does not reflect meaningful biological entities. Instead, G. salaris and G. thymalli are just one species comprising several pathogenic and non-pathogenic strains on various primary hosts. Following the International Code for Zoological Nomenclature, G. salaris Malmberg, 1957 is the valid species name with G. thymalli Žitňan, 1960 becoming the junior synonym. Accordingly, the range of G. salaris is significantly increased, given that formerly G. salaris-free countries such as e.g., Great Britain are now within the species' natural range. The synonymization of G. salaris and G. thymalli implies severe challenges to current disease management routines, which assume that G. salaris and G. thymalli are readily distinguishable. Protocols for reliable identification of pathogenic and non-pathogenic strains of G. salaris need to be developed.

Footprints of directional selection in wild Atlantic salmon populations: evidence for parasite-driven evolution?

Mechanisms of host-parasite co-adaptation have long been of interest in evolutionary biology; however, determining the genetic basis of parasite resistance has been challenging. Current advances in genome technologies provide new opportunities for obtaining a genome-scale view of the action of parasite-driven natural selection in wild populations and thus facilitate the search for specific genomic regions underlying inter-population differences in pathogen response. European populations of Atlantic salmon (Salmo salar L.) exhibit natural variance in susceptibility levels to the ectoparasite Gyrodactylus salaris Malmberg 1957, ranging from resistance to extreme susceptibility, and are therefore a good model for studying the evolution of virulence and resistance. However, distinguishing the molecular signatures of genetic drift and environment-associated selection in small populations such as land-locked Atlantic salmon populations presents a challenge, specifically in the search for pathogen-driven selection. We used a novel genome-scan analysis approach that enabled us to i) identify signals of selection in salmon populations affected by varying levels of genetic drift and ii) separate potentially selected loci into the categories of pathogen (G. salaris)-driven selection and selection acting upon other environmental characteristics. A total of 4631 single nucleotide polymorphisms (SNPs) were screened in Atlantic salmon from 12 different northern European populations. We identified three genomic regions potentially affected by parasite-driven selection, as well as three regions presumably affected by salinity-driven directional selection. Functional annotation of candidate SNPs is consistent with the role of the detected genomic regions in immune defence and, implicitly, in osmoregulation. These results provide new insights into the genetic basis of pathogen susceptibility in Atlantic salmon and will enable future searches for the specific genes involved.

Beyond MHC: signals of elevated selection pressure on Atlantic salmon (Salmo salar) immune-relevant loci

Using Atlantic salmon (Salmo salar) as a model system, we investigated whether 18 microsatellites tightly linked to immune-relevant genes have experienced different selection pressures than 76 loci with no obvious association with immune function. Immune-relevant loci were identified as outliers by two outlier tests significantly more often than nonimmune linked loci (22% vs. 1.6%). In addition, the allele frequencies of immune relevant markers were more often correlated with latitude and temperature. Combined, these results support the hypothesis that immune-relevant loci more frequently exhibit footprints of selection than other loci. They also indicate that the correlation between immune-relevant loci and latitude may be due to temperature-induced differences in pathogen-driven selection or some other environmental factor correlated with latitude.

Introgression of distant mitochondria into the genome of Gyrodactylus salaris: Nuclear and mitochondrial markers are necessary to identify parasite strains

Novel combinations of mitochondrial DNA (CO1) and internal transcribed spacers of nuclear ribosomal DNA (ITS) were detected among Gyrodactylus parasites on brown trout (Salmo trutta L.), rainbow trout (Oncorhynchus mykiss (Walbaum)), and Ohrid trout (Salmo letnica (Karaman)) from salmonid farms in Poland and Macedonia. Some clones differed from standard ITS only by ≤ 4 nucleotides, but they belonged to a mtDNA clade that differed from the Northern European lineages of G. salaris by d MCL = 0.266 ± 0.108 (maximum composite likelihood distance). The divergence of d MCL = 0.013 ± 0.005 within the alien mtDNA clade suggested introgression from an unknown maternal ancestor into the G. salaris Malmberg genome 137 to 57 kyr ago (or, less probably, repeated introgression). A comparable modern hybrid was detected based on permanently heterozygous ITS (28 bp/1264 = 2.2%) in a clone that is widespread throughout Finnish rainbow trout farms. This was a F1 hybrid of maternal G. pomeraniae Kuusela, Ziętara et Lumme (on roach, Rutilus rutilus (L.)) and G. lavareti Malmberg (on whitefish, Coregonus lavaretus (L.)). The mtDNA of the parental species differed by d MCL = 0.290 ± 0.130. The observations emphasize that both nuclear and maternally-inherited DNA markers are necessary to characterize sexually-produced lineages or clones of Gyrodactylus. The hybridization events demonstrated predict incongruence of mitochondrial vs. nuclear gene trees, i.e., reticulate evolution in species trees.