Highly Resolved Genomes of Two Closely Related Lineages of the Rodent Louse Polyplax serrata with Different Host Specificities.

Sucking lice of the parvorder Anoplura are permanent ectoparasites with specific lifestyle and highly derived features. Currently, genomic data are only available for a single species, the human louse Pediculus humanus. Here, we present genomes of two distinct lineages, with different host spectra, of a rodent louse Polyplax serrata. Genomes of these ecologically different lineages are closely similar in gene content and display a conserved order of genes, with the exception of a single translocation. Compared with P. humanus, the P. serrata genomes are noticeably larger (139 vs. 111 Mbp) and encode a higher number of genes. Similar to P. humanus, they are reduced in sensory-related categories such as vision and olfaction. Utilizing genome-wide data, we perform phylogenetic reconstruction and evolutionary dating of the P. serrata lineages. Obtained estimates reveal their relatively deep divergence (∼6.5 Mya), comparable with the split between the human and chimpanzee lice P. humanus and Pediculus schaeffi. This supports the view that the P. serrata lineages are likely to represent two cryptic species with different host spectra. Historical demographies show glaciation-related population size (Ne) reduction, but recent restoration of Ne was seen only in the less host-specific lineage. Together with the louse genomes, we analyze genomes of their bacterial symbiont Legionella polyplacis and evaluate their potential complementarity in synthesis of amino acids and B vitamins. We show that both systems, Polyplax/Legionella and Pediculus/Riesia, display almost identical patterns, with symbionts involved in synthesis of B vitamins but not amino acids.

Highly resolved genome assembly and comparative transcriptome profiling reveal genes related to developmental stages of tapeworm Ligula intestinalis.

Ligula intestinalis (Cestoda: Diphyllobothriidae) is an emerging model organism for studies on parasite population biology and host-parasite interactions. However, a well-resolved genome and catalogue of its gene content has not been previously developed. Here, we present the first genome assembly of L. intestinalis, based on Oxford Nanopore Technologies, Illumina and Omni-C sequencing methodologies. We use transcriptome profiling to compare plerocercoid larvae and adult worms and identify differentially expressed genes (DEGs) associated with these life stages. The genome assembly is 775.3 mega (M)bp in size, with scaffold N50 value of 118 Mbp and encodes 27 256 predicted protein-coding sequences. Over 60% of the genome consists of repetitive sequences. Synteny analyses showed that the 10 largest scaffolds representing 75% of the genome display high correspondence to full chromosomes of cyclophyllidean tapeworms. Mapping RNA-seq data to the new reference genome, we identified 3922 differentially expressed genes in adults compared with plerocercoids. Gene ontology analyses revealed over-represented genes involved in reproductive development of the adult stage (e.g. sperm production) and significantly enriched DEGs associated with immune evasion of plerocercoids in their fish host. This study provides the first insights into the molecular biology of L. intestinalis and provides the most highly contiguous assembly to date of a diphyllobothriid tapeworm useful for population and comparative genomic investigations of parasitic flatworms.

The history and organization of the Workshop on Population and Speciation Genomics.

With the advent of high-throughput genome sequencing, bioinformatics training has become essential for research in evolutionary biology and related fields. However, individual research groups are often not in the position to teach students about the most up-to-date methodology in the field. To fill this gap, extended bioinformatics courses have been developed by various institutions and provide intense training over the course of two or more weeks. Here, we describe our experience with the organization of a course in one of the longest-running extended bioinformatics series of workshops, the Evomics Workshop on Population and Speciation Genomics that takes place biennially in the UNESCO world heritage town of Ceský Krumlov, Czech Republic. We list the key ingredients that make this workshop successful in our view, explain the routine for workshop organization that we have optimized over the years, and describe the most important lessons that we have learned from it. We report the results of a survey conducted among past workshop participants that quantifies measures of effective teaching and provide examples of how the workshop setting has led to the cross-fertilisation of ideas and ultimately scientific progress. We expect that our account may be useful for other groups aiming to set up their own extended bioinformatics courses.

Effect of population size and selection on Toll-like receptor diversity in populations of Galápagos mockingbirds.

The interactions of evolutionary forces are difficult to analyse in free-living populations. However, when properly understood, they provide valuable insights into evolutionary biology and conservation genetics. This is particularly important for the interplay of genetic drift and natural selection in immune genes that confer resistance to disease. The Galápagos Islands are inhabited by four closely related species of mockingbirds (Mimus spp.). We used 12 different-sized populations of Galápagos mockingbirds and one population of their continental relative northern mockingbird (Mimus polyglottos) to study the effects of genetic drift on the molecular evolution of immune genes, the Toll-like receptors (TLRs: TLR1B, TLR4 and TLR15). We found that neutral genetic diversity was positively correlated with island size, indicating an important effect of genetic drift. However, for TLR1B and TLR4, there was little correlation between functional (e.g., protein) diversity and island size, and protein structural properties were largely conserved, indicating only a limited effect of genetic drift on molecular phenotype. By contrast, TLR15 was less conserved and even its putative functional polymorphism correlated with island size. The patterns observed for the three genes suggest that genetic drift does not necessarily dominate selection even in relatively small populations, but that the final outcome depends on the degree of selection constraint that is specific for each TLR locus.

Historical dispersal and host-switching formed the evolutionary history of a globally distributed multi-host parasite - The Ligula intestinalis species complex.

Studies on parasite biogeography and host spectrum provide insights into the processes driving parasite diversification. Global geographical distribution and a multi-host spectrum make the tapeworm Ligula intestinalis a promising model for studying both the vicariant and ecological modes of speciation in parasites. To understand the relative importance of host association and biogeography in the evolutionary history of this tapeworm, we analysed mtDNA and reduced-represented genomic SNP data for a total of 139 specimens collected from 18 fish-host genera across a distribution range representing 21 countries. Our results strongly supported the existence of at least 10 evolutionary lineages and estimated the deepest divergence at approximately 4.99-5.05 Mya, which is much younger than the diversification of the fish host genera and orders. Historical biogeography analyses revealed that the ancestor of the parasite diversified following multiple vicariance events and was widespread throughout the Palearctic, Afrotropical, and Nearctic between the late Miocene and early Pliocene. Cyprinoids were inferred as the ancestral hosts for the parasite. Later, from the late Pliocene to Pleistocene, new lineages emerged following a series of biogeographic dispersal and host-switching events. Although only a few of the current Ligula lineages show narrow host-specificity (to a single host genus), almost no host genera, even those that live in sympatry, overlapped between different Ligula lineages. Our analyses uncovered the impact of historical distribution shifts on host switching and the evolution of host specificity without parallel host-parasite co-speciation. Historical biogeography reconstructions also found that the parasite colonized several areas (Afrotropical and Australasian) much earlier than was suggested by only recent faunistic data.

Different hosts in different lakes: prevalence and population genetic structure of plerocercoids of Ligula intestinalis (Cestoda) in Czech water bodies.

Ligula intestinalis (Linnaeus, 1758) is a tapeworm parasite with a worldwide distribution that uses a wide variety of fish species as its second intermediate host. In the present study, we investigated the prevalence and population genetic structure of plerocercoids of L. intestinalis in five common cyprinoid species, roach Rutilus rutilus (Linnaeus), freshwater bream Abramis brama (Linnaeus), white bream Blicca bjoerkna (Linnaeus), bleak Alburnus alburnus (Linnaeus), and rudd Scardinius erythrophthalmus (Linnaeus), collected in six water bodies of the Czech Republic (Milada, Most, Medard, Jordán, Rímov and Lipno). Of the six study sites, the highest frequency of parasitism was recorded in Lake Medard (15%). The overall prevalence rate among the species was as follows: roach > rudd ≥ freshwater bream > bleak > white bream. Two mitochondrial genes (cytb and COI) were used to compare the population genetic structure of parasite populations using selected samples from the five fish species. The results of the phylogenetic analysis indicated that all populations of L. intestinalis were placed in Clade A, previously identified as the most common in Europe. At a finer scale, haplotype network and PCoA analyses indicated the possible emergence of host specificity of several mtDNA haplotypes to the freshwater bream. Moreover, pairwise Fixation indices (FST) revealed a significant genetic structure between the parasite population in freshwater bream and other host species. Parasite populations in roach not only showed the highest rate of prevalence but also depicted a maximum number of shared haplotypes with populations from bleak and rudd. Our results suggest that recent ecological differentiation might have influenced tapeworm populations at a fine evolutionary scale. Thus, the differences in prevalence between fish host species in different lakes might be influenced not only by the parasite's ecology, but also by its genetic diversity.

Different phylogenomic methods support monophyly of enigmatic 'Mesozoa' (Dicyemida + Orthonectida, Lophotrochozoa).

Dicyemids and orthonectids were traditionally classified in a group called Mesozoa, but their placement in a single clade has been contested and their position(s) within Metazoa is uncertain. Here, we assembled a comprehensive matrix of Lophotrochozoa (Metazoa) and investigated the position of Dicyemida (= Rhombozoa) and Orthonectida, employing multiple phylogenomic approaches. We sequenced seven new transcriptomes and one draft genome from dicyemids (Dicyema, Dicyemennea) and two transcriptomes from orthonectids (Rhopalura). Using these and published data, we assembled and analysed contamination-filtered datasets with up to 987 genes. Our results recover Mesozoa monophyletic and as a close relative of Platyhelminthes or Gnathifera. Because of the tendency of the long-branch mesozoans to group with other long-branch taxa in our analyses, we explored the impact of approaches purported to help alleviate long-branch attraction (e.g. taxon removal, coalescent inference, gene targeting). None of these were able to break the association of Orthonectida with Dicyemida in the maximum-likelihood trees. Contrastingly, the Bayesian analysis and site-specific frequency model in maximum-likelihood did not recover a monophyletic Mesozoa (but only when using a specific 50 gene matrix). The classic hypothesis on monophyletic Mesozoa is possibly reborn and should be further tested.

Unique genetic structure of the human tapeworm Dibothriocephalus latus from the Alpine lakes region - a successful adaptation?

Dibothriocephalus latus is the most frequent causative agent of fish-borne zoonosis (diphyllobothriosis) in Europe, where it is currently circulating mainly in the Alpine lakes region (ALR) and Russia. Three mitochondrial genes (cox1, cob and nad3) and 6 microsatellite loci were analysed to determine how is the recently detected triploidy/parthenogenesis in tapeworms from ALR displayed at the DNA level. A geographically distant population from the Krasnoyarsk Reservoir in Russia (RU-KR) was analysed as a comparative population. One or 2 alleles of each microsatellite locus was detected in plerocercoids from RU-KR, corresponding to the microsatellite pattern of a diploid organism. In contrast, 1­3 alleles were observed in tapeworms from ALR, in accordance with their triploidy. The high diversity of mitochondrial haplotypes in D. latus from RU-KR implied an original and relatively stable population, but the identical structure of mitochondrial genes of tapeworms from ALR was probably a consequence of a bottleneck typical of introduced populations. These results indicated that the diploid/sexually reproducing population from RU-KR was ancestral, located within the centre of the distribution of the species, and the triploid/parthenogenetically reproducing subalpine population was at the margin of the distribution. The current study revealed the allelic structure of the microsatellite loci in the triploid tapeworm for the first time.

Phylogeography of the parasitic mite Laelaps agilis in Western Palearctic shows lineages lacking host specificity but possessing different demographic histories.

BACKGROUND: Laelaps agilis C.L. Koch, 1836 is one the most abundant and widespread parasitic mite species in the Western Palearctic. It is a permanent ectoparasite associated with the Apodemus genus, which transmits Hepatozoon species via the host's blood. Phylogenetic relationships, genealogy and host specificity of the mite are uncertain in the Western Palearctic. Here, we investigated the population genetic structure of 132 individual mites across Europe from their Apodemus and Clethrionomys hosts. Phylogenetic relationships and genetic variation of the populations were analyzed using cytochrome c oxidase subunit I (COI) gene sequences. RESULTS: We recovered three main mtDNA lineages within L. agilis in the Western Palearctic, which differentiated between 1.02 and 1.79 million years ago during the Pleistocene period: (i) Lineage A, including structured populations from Western Europe and the Czech Republic, (ii) Lineage B, which included only a few individuals from Greece and the Czech Republic; and (iii) Lineage C, which comprised admixed populations from Western and Eastern Europe. Contrary to their population genetic differentiation, the lineages did not show signs of specificity to different hosts. Finally, we confirmed that the sympatric congener L. clethrionomydis is represented by a separated monophyletic lineage. CONCLUSION: Differences in the depth of population structure between L. agilis Lineages A and C, corroborated by the neutrality tests and demographic history analyses, suggested a stable population size in the structured Lineage A and a rapid range expansion for the geographically admixed Lineage C. We hypothesized that the two lineages were associated with hosts experiencing different glaciation histories. The lack of host specificity in L. agilis lineages was in contrast to the co-occurring highly host-specific lineages of Polyplax serrata lice, sharing Apodemus hosts. The incongruence was attributed to the differences in mobility between the parasites, allowing mites to switch hosts more often.

Comparative analysis of monozoic fish tapeworms Caryophyllaeus laticeps (Pallas, 1781) and recently described Caryophyllaeus chondrostomi Barcák, Oros, Hanzelová, Scholz, 2017, using microsatellite markers.

The monozoic tapeworm Caryophyllaeus laticeps has been characterized by five markedly different morphotypes largely corresponding to different fish hosts. Recently, the most distinct morphotype 4 from the common nase Chondrostoma nasus was studied in more details resulting in description of a new species Caryophyllaeus chondrostomi. The molecular study based on mitochondrial cox1 and ribosomal lsrDNA did not reveal any interspecific differences between C. laticeps and C. chondrostomi and did not provide any molecular support for recognition of these two species. In the current study, six polymorphic microsatellite markers were applied in order to detect molecular differences between the two species and to provide molecular evidence of validity of C. chondrostomi. While all six microsatellite loci were amplified in different geographic populations of C. laticeps, only two of them provided the amplification product in C. chondrostomi. Results on the Bayesian analysis assigned C. chondrostomi and all geographic populations of C. laticeps to distinct clusters. Neither any close relationships among C. laticeps populations nor specific position of C. chondrostomi were revealed. Contrary, the results of the principal coordinate analysis revealed striking genetic separation of C. chondrostomi with no overlaps with any of the C. laticeps population or morphotype. Caryophyllaeus chondrostomi very probably underwent morphological divergence as a result of ongoing speciation, but this process has not yet been accompanied by sufficient genetic divergence. In this particular case, microsatellites were proved to be better molecular discriminative markers than rDNA and mtDNA.

"Parasite turnover zone" at secondary contact: A new pattern in host-parasite population genetics.

We describe here a new pattern of population genetic structure in a host-parasite system that can arise after secondary contact of previously isolated populations. Due to different generation times, and therefore different tempos of molecular evolution, the host and parasite populations reach different degrees of genetic differentiation during their separation (e.g., in refugia). Consequently, upon secondary contact, the host populations are able to re-establish a single panmictic population across the area of contact, while the parasite populations stop their dispersal at the secondary contact zone and create a narrow hybrid zone. From the host's perspective, the parasite's hybrid zone functions on a microevolutionary scale as a "parasite turnover zone": while the hosts are passing from area A to area B, their parasites turn genetically from the area A genotypes to the area B genotypes. We demonstrate this novel pattern with a model composed of Apodemus mice and Polyplax lice by comparing maternally inherited markers (complete mitochondrial genomes, and complete genomes of the vertically transmitted symbiont Legionella polyplacis) with single nucleotide polymorphisms derived from louse genomic data. We discuss the circumstances that may lead to this pattern and possible reasons why it has been overlooked in studies of host-parasite population genetics.

Association between louse abundance and MHC II supertypes in Galápagos mockingbirds.

Major histocompatibility complex class II (MHC II) is an essential molecule triggering the adaptive immune response by the presentation of pathogens to helper T cells. The association between individual MHC II variants and various parasites has become a frequent finding in studies of vertebrate populations. However, although bird ectoparasites have a significant effect on their host's fitness, and the host's immune system can regulate ectoparasitic infections, no study has yet investigated the association between MHC II polymorphism and ectoparasite infection in the populations of free-living birds. Here, we test whether an association exists between the abundance of a chewing louse (Myrsidea nesomimi) and MHC II polymorphism of its hosts, the Galápagos mockingbirds (Mimus). We have found that the presence of two MHC II supertypes (functionally differentiated clusters) was significantly associated with louse abundance. This pattern supports the theory that a co-evolutionary interaction stands behind the maintenance of MHC polymorphism. Moreover, we have found a positive correlation between louse abundance and heterophil/lymphocyte ratio (an indicator of immunological stress) that serves as an additional piece of evidence that ectoparasite burden is affected by immunological state of Galápagos mockingbirds.

Author Correction: From taxonomic deflation to newly detected cryptic species: Hidden diversity in a widespread African squeaker catfish.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

From taxonomic deflation to newly detected cryptic species: Hidden diversity in a widespread African squeaker catfish.

Cryptic genetic diversity and erroneous morphological species determination represent frequent problems in biodiversity research. Here, examination of 138 specimens of Synodontis (Mochokidae, Siluriformes) from the Nile River and Lake Turkana revealed the presence of both S. schall-like and S. frontosus-like morphotypes, with a phenotypic gradient between them. We concluded phylogenetic and population genetic analyses based on two mitochondrial and one nuclear marker including 131 coxI (565 bp), 96 cytb (973 bp) and 19 RAG2 (896 bp) sequences from the Nile-Turkana population, plus additional GenBank data of Synodontis spp. Whilst nuclear data were inconclusive, mitochondrial sequences suggested that both morphotypes and intermediate forms are conspecific. The results imply probable synonymy of S. frontosus with S. schall. Conversely, a strong biogeographical signal was revealed among widely distributed and supposedly conspecific S. schall-like catfish of the Nilo-Sudanian ichthyological province. Synodontis schall sensu stricto (=Eastern clade), as defined by type locality in the Nile, is apparently restricted to the eastern part of the Nilo-Sudanian ichthyological province (e.g. Nile, Turkana, Chad). Synodontis schall Western clade (Senegambia, Niger, Chad) most probably represents a cryptic taxon, unrecognized thus far due to the absence of distinctive morphological differences.

Population co-divergence in common cuttlefish (Sepia officinalis) and its dicyemid parasite in the Mediterranean Sea.

Population structure and biogeography of marine organisms are formed by different drivers than in terrestrial organisms. Yet, very little information is available even for common marine organisms and even less for their associated parasites. Here we report the first analysis of population structure of both a cephalopod host (Sepia officinalis) and its dicyemid parasite, based on a homologous molecular marker (cytochrome oxidase I). We show that the population of common cuttlefish in the Mediterranean area is fragmented into subpopulations, with some areas featuring restricted level of gene flow. Amongst the studied areas, Sardinia was genetically the most diverse and Cyprus the most isolated. At a larger scale, across the Mediterranean, the population structure of the parasite shows co-diversification pattern with its host, but a slower rate of diversification. Differences between the two counterparts are more obvious at a finer scale, where parasite populations show increased level of fragmentation and lower local diversities. This discrepancy can be caused by local extinctions and replacements taking place more frequently in the dicyemid populations, due to their parasitic lifestyle.

Host specificity driving genetic structure and diversity in ectoparasite populations: Coevolutionary patterns in Apodemus mice and their lice.

A degree of host specificity, manifested by the processes of host-parasite cospeciations and host switches, is assumed to be a major determinant of parasites' evolution. To understand these patterns and formulate appropriate ecological hypotheses, we need better insight into the coevolutionary processes at the intraspecific level, including the maintenance of genetic diversity and population structure of parasites and their hosts. Here, we address these questions by analyzing large-scale molecular data on the louse Polyplax serrata and its hosts, mice of the genus Apodemus, across a broad range of European localities. Using mitochondrial DNA sequences and microsatellite data, we demonstrate the general genetic correspondence of the Apodemus/Polyplax system to the scenario of the postglacial recolonization of Europe, but we also show several striking discrepancies. Among the most interesting are the evolution of different degrees of host specificity in closely related louse lineages in sympatry, or decoupled population structures of the host and parasites in central Europe. We also find strong support for the prediction that parasites with narrower host specificity possess a lower level of genetic diversity and a deeper pattern of interpopulation structure as a result of limited dispersal and smaller effective population size.

Development of microsatellite loci in zoonotic tapeworm Dibothriocephalus latus (Linnaeus, 1758), Lühe, 1899 (syn. Diphyllobothrium latum) using microsatellite library screening.

The broad fish tapeworm Dibothriocephalus latus is a causative agent of human food-borne disease called diphyllobothriosis. Medical importance, scattered geographical distribution and unknown origin of D. latus in Europe and North America make this species to be an interesting model for population genetics. Microsatellite markers were originally designed by library screening using NGS approach and validated as tools for future studies on population genetics of D. latus. Out of 122 candidates selected after NGS analysis, 110 yielded PCR products of the expected size, and in 78 of them, a declared repetitive motif was confirmed by Sanger sequencing. After the fragment analysis, six loci were proved to be polymorphic and tested for observed (Ho) and expected (He) heterozygosity, and deviations from Hardy-Weinberg equilibrium (HWE). They promise future application in studies on genetic interrelationships, origin and migratory routes of this medically important emerging tapeworm.

Tour around the globe: The case of invasive tapeworm Atractolytocestus huronensis (Cestoda: Caryophyllidea), a parasite of common carp.

The monozoic tapeworm Atractolytocestus huronensis Anthony, 1958 (Cestoda: Caryophyllidea), an intestinal parasite of the common carp, is characterized by its invasive character and potential to colonize new territories. It was initially described from North America and has also been found in several European countries. The most recent findings of A. huronensis originated from China and South Africa; however, no data on genetic relationships of these populations were available. The current study provides the first molecular characterisation of A. huronensis from South Africa and China using a partial sequence of mitochondrial cytochrome c oxidase subunit 1 (cox1) and a complete ribosomal ITS2 spacer. Ribosomal and mitochondrial data were applied for phylogenetic analyses in order to assess the genetic interrelationships among global A. huronensis populations. Divergent intragenomic copies of ribosomal ITS2 were detected in all analysed specimens; the structure and frequency of the ITS2 variants of tapeworms from China and South Africa corresponded with the data on ITS2 paralogues observed previously in A. huronensis from Slovakia, the United States and the United Kingdom. The phylogenetic analysis of cox1 indicated that A. huronensis exist in two slightly differentiated clusters; one cluster was supported by all phylogenetic approaches (NJ, ML, BI) and was represented by samples from China, the USA and the UK. A second cluster was represented by tapeworms from continental Europe (Slovakia, Hungary, Romania, Croatia) and South Africa. Haplotype network analysis revealed that the highest population diversity occurs in China. The results provide useful pilot information about the interrelationships of A. huronensis on four continents and indicate that China, or the eastern Palaearctic, served as the original source population for the global expansion of this invasive tapeworm. Data on the origin and distribution of the common carp, the only specific host of A. huronensis, are also discussed.

Development of polymorphic microsatellites for the invasive Asian fish tapeworm Schyzocotyle acheilognathi.

We describe the development of ready-to-use set of fifteen polymorphic microsatellite markers to benefit future population biology and phylogeographic studies on the invasive Asian fish tapeworm Schyzocotyle acheilognathi. The microsatellite loci were selected from partial Illumina shotgun genome sequences of three parasite specimens and their universality tested on a set of 12 geographically distant populations of the parasite. Particularly low levels of heterozygosity have been detected in the Chinese population pointing towards possible hidden population structure that deserves further attention in future population genetic studies.