Speciation in parasites: a population genetics approach

Molecular phylogeny by Bayesian analysis and description of a new species of Henneguya, a parasitic cnidarian from Astyanax bimaculatus (Linnaeus, 1758) (Characiformes: Acestrorhamphidae) in Caatinga: A biome with neglected biodiversity data

Myxozoans are cnidarian parasites with a worldwide distribution and represent some of the most spectacular examples of parasite radiation. Based on an integrative taxonomic approach that includes spore morphology and molecular data, as well as other phenotypic characters such as host specificity, a new myxozoan, Henneguya caraubensis n. sp., is described from the gills of Astyanax bimaculatus (Linnaeus, 1758). Plasmodia were located within the gill lamellar epithelium, without connective tissue surrounding them. Mature myxospores were ellipsoidal in the frontal view, measuring 27.2 ± 2.6 μm in total length, 10.3 ± 0.7 μm in body length of spore, and 3.9 ± 0.2 μm in width. Two polar capsules elongated pyriform in shape, equal in size, measuring 3.4 ± 0.3 μm in length and 1.2 ± 0.1 μm in width. Polar tubules had 5 to 6 coils. Histopathological analysis showed large intralamellar polysporic plasmodia associated with displacement, deformation, and fusion of the nearest gill lamellae and with the adjacent lamellae being pushed laterally. Bayesian inference SSU rDNA-based phylogenetic analyses showed a strong tendency of Henneguya/Myxobolus species to form clusters mainly based on the order and/or family of the host, despite having different geographic origins. Phylogenetically, the novel species grouped within a clade containing Henneguya/Myxobolus spp. and appears in a well-supported subclade of Henneguya species parasites of Astyanax fishes. Our study reinforces highly endemic diversity still hidden that inhabits the Caatinga biome and points out the need to improve the strategies of protection and conservation of this threatened ecoregion with neglected biodiversity data.

Coral-associated Symbiodiniaceae exhibit host specificity but lack phylosymbiosis, with Cladocopium and Durusdinium showing different cophylogenetic patterns

Altering the composition of the Symbiodiniaceae community to adapt to anomalous sea water warming represents a potential survival mechanism for scleractinian corals. However, the processes of Symbiodiniaceae assembly and long-standing evolution of coral-Symbiodiniaceae interactions remain unclear. Here, we utilized ITS2 (internal transcribed spacer 2) amplicon sequencing and the SymPortal framework to investigate the diversity and specificity of Symbiodiniaceae across 39 scleractinian coral species. Furthermore, we tested phylosymbiosis and cophylogeny between coral hosts and their Symbiodiniaceae. In our study, environmental samples exhibited the highest Symbiodiniaceae diversity. Cladocopium and Durusdinium dominated the Symbiodiniaceae communities, with significant β-diversity differences among coral species. Additionally, host specificity was widespread in Symbiodiniaceae, especially Durusdinium spp., yet lacked a phylosymbiotic pattern. Moreover, Cladocopium spp. showed cophylogenetic congruence with their hosts, while there was no evidence for Durusdinium spp. Furthermore, host switching was the predominant evolutionary event, implying its contribution to Cladocopium diversification. These findings suggest that Symbiodiniaceae assembly does not recapitulate host phylogeny, and host specificity alone does not drive phylosymbiosis or cophylogeny. As environmental reservoirs, free-living Symbiodiniaceae may influence symbiotic communities. Additionally, Durusdinium-coral associations lack cophylogenetic signals, indicating more flexible partnerships than Cladocopium. Overall, our results enhance understanding of Symbiodiniaceae assembly and coral-Symbiodiniaceae evolutionary interactions.

Insights into the genetic diversity and species distribution of Oswaldocruzia nematodes (Trichostrongylida: Molineidae) in Europe: apparent absence of geographic and population structuring in amphibians

The genus Oswaldocruzia represents a taxonomically diverse group of nematodes with global distribution. Although Oswaldocruzia species are widespread and exhibit a remarkably wide host range in some species, their genetic diversity and biogeographic patterns remain poorly understood. This study investigated the genetic variability and distribution of Oswaldocruzia spp. in nine anuran species from the genera Bufo, Bufotes, Pelophylax, and Rana across Central Europe and the Balkans. Two species were identified: Oswaldocruzia filiformis and O. ukrainae, each exhibiting a different range of host associations. Phylogenetic analyses based on mitochondrial COI sequences revealed significant haplotype diversity in the generalist O. filiformis, with low geographic and host-associated genetic structuring. In contrast, O. ukrainae, which is closely associated with Bufotes viridis, exhibited only one genetic variant across all samples, highlighting its restricted genetic diversity. The findings emphasize contrasting genetic diversities among nematode parasites exhibiting different levels of host-specificity and expand the known distribution of O. filiformis into new regions of the Balkans. In addition, they highlight the need for additional studies on the ecological and evolutionary factors that influence the genetic diversity of parasites in amphibians.

Legacy parasite collections reveal species-specific population genetic patterns among three species of zoonotic schistosomes

Studies estimating genetic diversity and population structure in multi-host parasites are often constrained by temporally and spatially limited sampling. This study addresses these limitations by analyzing globally distributed samples of three congeneric avian schistosomes (Trematoda: Schistosomatidae: Trichobilharzia), including collections spanning 20 years archived at The Museum of Southwestern Biology, Parasites Division. The three species exhibited significant differences in population genetic parameters across one nuclear and two mitochondrial loci. Trichobilharzia querquedulae (TQ) maintained a well-connected, globally diverse metapopulation, with an effective population size approximately three times larger than that of the other two species, T. physellae (TP) and Trichobilharzia sp. A (TA). TP and TA had lower overall genetic diversity and greater population structure. These differences are likely shaped by the ecologies of the duck definitive hosts that disperse these parasites. This study highlights the value of natural history collections, particularly since Trichobilharzia is a key agent of zoonotic cercarial dermatitis, a disease whose etiology and epidemiology remain poorly understood. Within a comparative congeneric framework, population genetic data can provide insights into host-parasite natural history and its influence on microevolutionary patterns, including contributions to zoonotic disease.

Utility of 28S Ribosomal RNA Gene Domains for Molecular Classification and Phylogeny of Rhinonyssid Mites

The family Rhinonyssidae includes endoparasitic, blood-feeding mites that are parasitic on birds and that remain largely unstudied despite their potential role as vectors or reservoirs for various pathogens, like other Dermanyssoidea. Traditionally, the taxonomy of the group has been based on morphometric characteristics, which makes identification very difficult in many groups of closely related species. On the other hand, studies on the phylogenetic relationships within this group of mites have been neglected until the early years of the present century. In this study, twelve species belonging to five different species complexes were identified, and domains D1-D3 28S rRNA of each one were sequenced, for the first time, to investigate the sequence variation and its taxonomic implications for phylogenetic inference. Our data indicate that this molecular marker can effectively differentiate between species within the "motacillae", "sairae", "pari", and "hirsti" complexes of the genus Ptilonyssus and the "melloi" and "columbae" complexes of the genus Tinaminyssus. Furthermore, the phylogenetic tree that can be derived from the domain D1-D3 28S rRNA sequences presented in this study is congruent with the current taxonomy of the Rhinonyssidae. This research calls for a reassessment of the taxonomic status of some group of species.

Adaptive evolution of stress response genes in parasites aligns with host niche diversity

BACKGROUND

Stress responses are key the survival of parasites and, consequently, also the evolutionary success of these organisms. Despite this importance, our understanding of the evolution of molecular pathways dealing with environmental stressors in parasitic animals remains limited. Here, we tested the link between adaptive evolution of parasite stress response genes and their ecological diversity and species richness. We comparatively investigated antioxidant, heat shock, osmoregulatory, and behaviour-related genes (foraging) in two model parasitic flatworm lineages with contrasting ecological diversity, Cichlidogyrus and Kapentagyrus (Platyhelminthes: Monopisthocotyla), through whole-genome sequencing of 11 species followed by in silico exon bait capture as well as phylogenetic and codon analyses.

RESULTS

We assembled the sequences of 48 stress-related genes and report the first foraging (For) gene orthologs in flatworms. We found duplications of heat shock (Hsp) and oxidative stress genes in Cichlidogyrus compared to Kapentagyrus. We also observed positive selection patterns in genes related to mitochondrial protein import (Hsp) and behaviour (For) in species of Cichlidogyrus infecting East African cichlids-a host lineage under adaptive radiation. These patterns are consistent with a potential adaptation linked to a co-radiation of these parasites and their hosts. Additionally, the absence of cytochrome P450 and kappa and sigma-class glutathione S-transferases in monogenean flatworms is reported, genes considered essential for metazoan life.

CONCLUSIONS

This study potentially identifies the first molecular function linked to a flatworm radiation. Furthermore, the observed gene duplications and positive selection indicate the potentially important role of stress responses for the ecological adaptation of parasite species.

Understanding the Influence of Host Radiation on Symbiont Speciation through Parasites of Species Flocks

(Adaptive) radiations have attracted evolutionary biologists for a long time as ideal model systems to study patterns and processes of often rapid speciation. However, whereas a wealth of (sometimes already genome-scale) data is available for host radiations, very few studies target the patterns of diversification in their symbionts, even though they would be excellent models to study symbiont speciation. Our review summarizes what little is known about general patterns of symbiont diversification in often iconic adaptive host radiations and to what extent these patterns are dependent on the evolutionary trajectories of their hosts. We identify research gaps that need to be addressed in the future and discuss the potential of approaches not yet typically used in these study systems, such as epidemiological disease modeling and new omics technologies, for significantly advancing our understanding of these complex eco-evolutionary relationships.

Evolutionary Analysis of Hypoderma Pantholopsum in Tibetan Antelopes on the Qinghai-Tibetan Plateau

PURPOSE

Hypoderma pantholopsum is a parasite that parasitizes Tibetan antelopes (Pantholops hodgsonii). This study aims was to reveal the genetic diversity within H. pantholopsum and contribute to the protection of Tibetan antelope.

METHODS

H. pantholopsum was subjected to mitochondrial genome sequencing, annotation, and phylogenetic analysis. PhyloSuite and BEAST were used to construct phylogenetic tree and divergence time for the parasite.

RESULTS

The complete H. pantholopsum genome was 16,265 bp in length. The complete mitochondrial genomes contained 37 typical genes, which included 13 protein-coding genes (PCGs), 22 tRNAs, and 2 rRNAs. Phylogenetic trees constructed based on the 18S rRNA, 28S rRNA, and mitochondrial genome sequences showed that H. pantholopsum clustered on the same branch as the Hypoderma species in the GenBank database. According to the divergence time for the COI gene, H. pantholopsum emerged and differentiated approximately 11.59 million years ago (Mya), which indicates that H. pantholopsum appeared much earlier than H. bovis and H. sinense in the genus Hypoderma.

CONCLUSION

The present study explored that the complete mitochondrial genome of H. pantholopsum, along with the phylogenetic evolution, and divergence time estimation, provide valuable data for future investigations into the phylogeny and differentiation of Hypoderma species on the Qinghai-Tibetan Plateau.

Gray foxes (Urocyon cinereoargenteus) as natural definitive hosts of Paragonimus mexicanus

Paragonimus mexicanus is a trematode that causes pulmonary and extrapulmonary infections in humans, characterized by chest pain, dyspnea, fever, and weight loss. The detection of Paragonimus spp. is primarily achieved through the microscopic observation of eggs in feces, sputum, and pleural fluid. Paragonimus mexicanus has been found in various wild animals, including dogs, cats, raccoons, and opossums. Although the reservoirs of P. mexicanus in Costa Rica are unknown, this study analyzed fecal samples from gray foxes (Urocyon cinereoargenteus) using microscopic and molecular methods. In the morphological analysis, characteristic eggs of the genus Paragonimus were identified. DNA was extracted from fecal samples, and a fragment of the ITS2 loci of trematodes was amplified, which showed a 100 % similarity with P. mexicanus metacercariae from crabs in Ecuador. Then, a Bayesian inference phylogenetic analysis was performed with the obtained data and pre-existing sequences of P. mexicanus found in America, showing that our sequence clustered firstly with others from Colima and Veracruz (Mexico), and Ecuador, while a second cluster contained sequences from Chiapas (Mexico), Ecuador, and Guatemala. These results provide evidence of the presence of P. mexicanus in the gray fox and suggest its role as a possible new wild reservoir, which could have zoonotic implications for the infection of other animal species and humans. Additionally, our phylogenetic analysis reveals low genetic differentiation among the compared P. mexicanus populations and the possibility of additional Paragonimus species currently classified as P. mexicanus. The finding of this parasite in our country, together with comparisons with previous studies, highlights the complex evolutionary history and population dynamics of P. mexicanus.

Testing the radiation cascade in postglacial radiations of whitefish and their parasites: founder events and host ecology drive parasite evolution

Reciprocal effects of adaptive radiations on the evolution of interspecific interactions, like parasitism, remain barely explored. We test whether the recent radiations of European whitefish (Coregonus spp.) across and within perialpine and subarctic lakes promote its parasite Proteocephalus fallax (Platyhelminthes: Cestoda) to undergo host repertoire expansion via opportunity and ecological fitting, or adaptive radiation by specialization. Using de novo genomic data, we examined P. fallax differentiation across lakes, within lakes across sympatric host species, and the contributions of host genetics versus host habitat use and trophic preferences. Whitefish intralake radiations prompted parasite host repertoire expansion in all lakes, whereas P. fallax differentiation remains incipient among sympatric fish hosts. Whitefish genetic differentiation per se did not explain the genetic differentiation among its parasite populations, ruling out codivergence with the host. Instead, incipient parasite differentiation was driven by whitefish phenotypic radiation in trophic preferences and habitat use in an arena of parasite opportunity and ecological fitting to utilize resources from emerging hosts. Whilst the whitefish radiation provides a substrate for the parasite to differentiate along the same water-depth ecological axis as Coregonus spp., the role of the intermediate hosts in parasite speciation may be overlooked. Parasite multiple-level ecological fitting to both fish and crustacean intermediate hosts resources may be responsible for parasite population substructure in Coregonus spp. We propose parasites' delayed arrival was key to the initial burst of postglacial intralake whitefish diversification, followed by opportunistic tapeworm host repertoire expansion and a delayed nonadaptive radiation cascade of incipient tapeworm differentiation. At the geographical scale, dispersal, founder events, and genetic drift following colonization of spatially heterogeneous landscapes drove strong parasite differentiation. We argue that these microevolutionary processes result in the mirroring of host-parasite phylogenies through phylogenetic tracking at macroevolutionary and geographical scales.

Species clarification of fairy inkcap ("Coprinellus disseminatus") in China

Coprinellus disseminatus and other morphologically similar species are widely dispersed worldwide and are commonly referred to as "fairy inkcap". Based on the molecular phylogenetic study and morphological observation, a thorough investigation was carried out utilising 74 collections of related species that were gathered from seventeen provinces and five Chinese fungaria between 1998 and 2023 and revealed 11 lineages of "fairy inkcap", nine of which were found in China, and which belonged to the two genera Coprinellus and Tulosesus. In sect. Disseminati, genetic diversities (π), and fixation index (Fst) amongst lineages were computed, and a haplotype-based network was established to ascertain the relationships amongst each clade. A new section of Coprinellus, sect. Aureodisseminati, were discovered. In addition, four new species (C. aureodisseminatus, C. austrodisseminatus, C. parcus, and C. velutipes), a new subspecies of C. disseminatus, a new combination (Tulosesus pseudodisseminatus), the first discovery of epigamous type of C. magnoliae and a new record to China (T. subdisseminatus) were also identified and thoroughly described with accompanying illustrations. Their differences in macro- and micro-features, as well as their character sequence, were discussed.

Genome-wide phylogenetic and genetic evolutionary analyses of mitochondria in Hypoderma bovis and H. sinense on the Qinghai-Tibetan Plateau

Hypoderma bovis (H. bovis) and Hypoderma sinense (H. sinense) are insects that cause hypodermosis in yaks and Bos taurus. Hypodermosis is a severe skin condition that not only impairs the development of local animal husbandry but also poses threats to human health as a zoonosis. The Qinghai-Tibetan Plateau (QTP) is known as the "Roof of the World." Its unique geographical environment and climate conditions have supported the growth of a wide range of mammals, providing favorable conditions for Hypoderma spp. to complete their life cycles. In this study, the whole mitochondrial genomes of H. bovis and H. sinense collected from the QTP were sequenced and phylogenetically analyzed. We found that the whole genomes of H. bovis and H. sinense are 16,283 bp and 16,300 bp in length, respectively. Both the H. bovis and H. sinense genomes have 37 mitochondrial genes, which include two rRNA genes (16S rRNA and 12S rRNA), 22 tRNA genes, the control region (D-loop region), the light chain replication initiation region, and 13 protein-coding genes (PCGs). The phylogenetic tree generated based on the 13 PCGs revealed close phylogenetic relationships between H. sinense, H. bovis, and Hypoderma lineatum. A similar result was also found in our phylogenetic analysis based on 18S rRNA and 28S rRNA. However, analysis of cytochrome oxidase subunit I (COI) showed cluster of H. bovis, H. sinense, and Cuterebra spp. on the same branch, all belonging to Oestridae. The differentiation time generated based on 13 PCGs indicates that H. bovis and H. sinense differentiated and formed ~4.69 million years ago (Mya) and ~4.06 Mya, respectively. This timing coincides with the differentiation and appearance of yak and Bos taurus in the Pliocene (~4.7 Mya), indicating that the parasites and mammals diverged in close temporal proximity. Of note, this period also witnessed a rapid uplift of the QTP, causing significant climate and environmental changes. Thus, we conjecture that the differentiation of Hypoderma spp. is potentially related to the differentiation of their host species, as well as climate changes caused by the uplift of the QTP. Overall, our study can provide valuable data to support further studies on the phylogeny and differentiation of Hypoderma spp. on the QTP.

POPULATION STRUCTURE AND BARTONELLA QUINTANA IN HEAD AND BODY LICE IN POKHARA, NEPAL (ANOPLURA: PEDICULIDAE)

This study examined the population structure of head and body lice infesting a random sample of people in Pokhara, Nepal during 2003, 2004, and 2005. A total of 106 participants (6 to 72 yr old, median = 12) volunteered to have lice collected from their heads and clothing. Most participants (70%) harbored only head lice, some (15%) had only body lice, and some (15%) had concurrent infestations of head and body lice (dual infestations). A total of 1,472 lice was collected. Significantly more nymphs were collected than adult lice. Louse populations were generally small (geometric mean = 8.8 lice per person) but a few participants harbored larger louse populations (maximum = 65 lice). People with dual infestations harbored significantly more lice than people with single infestations; however, there was no difference in the infestation intensities between people infested with head lice only vs. those infested with body lice only. Male participants harbored significantly more lice than did females. There were no significant differences in infestation intensity due to participant age or their socioeconomic level. The sex ratio of adult lice was increasingly female biased with increasing adult louse density. Infection of lice with Bartonella quintana was low (ca. 1.5%). Pediculosis is a common problem in urban Nepal.

Population Genomics of Pooled Samples: Unveiling Symbiont Infrapopulation Diversity and Host-Symbiont Coevolution

Microscopic symbionts represent crucial links in biological communities. However, they present technical challenges in high-throughput sequencing (HTS) studies due to their small size and minimal high-quality DNA yields, hindering our understanding of host-symbiont coevolution at microevolutionary and macroevolutionary scales. One approach to overcome those barriers is to pool multiple individuals from the same infrapopulation (i.e., individual host) and sequence them together (Pool-Seq), but individual-level information is then compromised. To simultaneously address both issues (i.e., minimal DNA yields and loss of individual-level information), we implemented a strategic Pool-Seq approach to assess variation in sequencing performance and categorize genetic diversity (single nucleotide polymorphisms (SNPs)) at both the individual-level and infrapopulation-level for microscopic feather mites. To do so, we collected feathers harboring mites (Proctophyllodidae: Amerodectes protonotaria) from four individual Prothonotary Warblers (Parulidae: Protonotaria citrea). From each of the four hosts (i.e., four mite infrapopulations), we conducted whole-genome sequencing on three extraction pools consisting of different numbers of mites (1 mite, 5 mites, and 20 mites). We found that samples containing pools of multiple mites had more sequencing reads map to the feather mite reference genome than did the samples containing only a single mite. Mite infrapopulations were primarily genetically structured by their associated individual hosts (not pool size) and the majority of SNPs were shared by all pools within an infrapopulation. Together, these results suggest that the patterns observed are driven by evolutionary processes occurring at the infrapopulation level and are not technical signals due to pool size. In total, despite the challenges presented by microscopic symbionts in HTS studies, this work highlights the value of both individual-level and infrapopulation-level sequencing toward our understanding of host-symbiont coevolution at multiple evolutionary scales.

Mitogenomic evolutionary rates in bilateria are influenced by parasitic lifestyle and locomotory capacity

The evidence that parasitic animals exhibit elevated mitogenomic evolutionary rates is inconsistent and limited to Arthropoda. Similarly, the evidence that mitogenomic evolution is faster in species with low locomotory capacity is limited to a handful of animal lineages. We hypothesised that these two variables are associated and that locomotory capacity is a major underlying factor driving the elevated rates in parasites. Here, we study the evolutionary rates of mitogenomes of 10,906 bilaterian species classified according to their locomotory capacity and parasitic/free-living life history. In Bilateria, evolutionary rates were by far the highest in endoparasites, much lower in ectoparasites with reduced locomotory capacity and free-living lineages with low locomotory capacity, followed by parasitoids, ectoparasites with high locomotory capacity, and finally micropredatory and free-living lineages. The life history categorisation (parasitism) explained ≈45%, locomotory capacity categorisation explained ≈39%, and together they explained ≈56% of the total variability in evolutionary rates of mitochondrial protein-coding genes in Bilateria. Our findings suggest that these two variables play major roles in calibrating the mitogenomic molecular clock in bilaterian animals.

First Insights into Population Structure and Genetic Diversity Versus Host Specificity in Trypanorhynch Tapeworms Using Multiplexed Shotgun Genotyping

Theory predicts relaxed host specificity and high host vagility should contribute to reduced genetic structure in parasites while strict host specificity and low host vagility should increase genetic structure. Though these predictions are intuitive, they have never been explicitly tested in a population genomic framework. Trypanorhynch tapeworms, which parasitize sharks and rays (elasmobranchs) as definitive hosts, are the only order of elasmobranch tapeworms that exhibit considerable variability in their definitive host specificity. This allows for unique combinations of host use and geographic range, making trypanorhynchs ideal candidates for studying how these traits influence population-level structure and genetic diversity. Multiplexed shotgun genotyping (MSG) data sets were generated to characterize component population structure and infrapopulation diversity for a representative of each trypanorhynch suborder: the ray-hosted Rhinoptericola megacantha (Trypanobatoida) and the shark-hosted Callitetrarhynchus gracilis (Trypanoselachoida). Adults of R. megacantha are more host-specific and less broadly distributed than adults of C. gracilis, allowing correlation between these factors and genetic structure. Replicate tapeworm specimens were sequenced from the same host individual, from multiple conspecific hosts within and across geographic regions, and from multiple definitive host species. For R. megacantha, population structure coincided with geography rather than host species. For C. gracilis, limited population structure was found, suggesting a potential link between degree of host specificity and structure. Conspecific trypanorhynchs from the same host individual were found to be as, or more, genetically divergent from one another as from conspecifics from different host individuals. For both species, high levels of homozygosity and positive FIS values were documented.

Loss of species and genetic diversity during colonization: Insights from acanthocephalan parasites in northern European seals

Studies on host-parasite systems that have experienced distributional shifts, range fragmentation, and population declines in the past can provide information regarding how parasite community richness and genetic diversity will change as a result of anthropogenic environmental changes in the future. Here, we studied how sequential postglacial colonization, shifts in habitat, and reduced host population sizes have influenced species richness and genetic diversity of Corynosoma (Acanthocephala: Polymorphidae) parasites in northern European marine, brackish, and freshwater seal populations. We collected Corynosoma population samples from Arctic, Baltic, Ladoga, and Saimaa ringed seal subspecies and Baltic gray seals, and then applied COI barcoding and triple-enzyme restriction-site associated DNA (3RAD) sequencing to delimit species, clarify their distributions and community structures, and elucidate patterns of intraspecific gene flow and genetic diversity. Our results showed that Corynosoma species diversity reflected host colonization histories and population sizes, with four species being present in the Arctic, three in the Baltic Sea, two in Lake Ladoga, and only one in Lake Saimaa. We found statistically significant population-genetic differentiation within all three Corynosoma species that occur in more than one seal (sub)species. Genetic diversity tended to be high in Corynosoma populations originating from Arctic ringed seals and low in the landlocked populations. Our results indicate that acanthocephalan communities in landlocked seal populations are impoverished with respect to both species and intraspecific genetic diversity. Interestingly, the loss of genetic diversity within Corynosoma species seems to have been less drastic than in their seal hosts, possibly due to their large local effective population sizes resulting from high infection intensities and effective intra-host population mixing. Our study highlights the utility of genomic methods in investigations of community composition and genetic diversity of understudied parasites.

Novel insights into symbiont population structure: Globe-trotting avian feather mites contradict the specialist-generalist variation hypothesis

Researchers often examine symbiont host specificity as a species-level pattern, but it can also be key to understanding processes occurring at the population level, which are not as well understood. The specialist-generalist variation hypothesis (SGVH) attempts to explain how host specificity influences population-level processes, stating that single-host symbionts (specialists) exhibit stronger population genetic structure than multi-host symbionts (generalists) because of fewer opportunities for dispersal and more restricted gene flow between populations. However, this hypothesis has not been tested in systems with highly mobile hosts, in which population connectivity may vary temporally and spatially. To address this gap, we tested the SGVH on proctophyllodid feather mites found on migratory warblers (family Parulidae) with contrasting host specificities, Amerodectes protonotaria (a host specialist of Protonotaria citrea) and A. ischyros (a host generalist of 17 parulid species). We used a pooled-sequencing approach and a novel workflow to analyse genetic variants obtained from whole genome data. Both mite species exhibited fairly weak population structure overall, and contrary to predictions of the SGVH, the generalist was more strongly structured than the specialist. These results may suggest that specialists disperse more freely among conspecifics, whereas generalists sort according to geography. Furthermore, our results may reflect an unexpected period for mite transmission - during the nonbreeding season of migratory hosts - as mite population structure more closely reflects the distributions of hosts during the nonbreeding season. Our findings alter our current understanding of feather mite biology and highlight the potential for studies to explore factors driving symbiont diversification at multiple evolutionary scales.

Eimeria zuernii (Eimeriidae: Coccidia): mitochondrial genome and genetic diversity in the Chinese yak

BACKGROUND

Coccidiosis caused by Eimeria zuernii (Eimeriidae: Coccidia) represents a significant economic threat to the bovine industry. Understanding the evolutionary and genetic biology of E. zuernii can assist in new interaction developments for the prevention and control of this protozoosis.

METHODS

We defined the evolutionary and genetic characteristics of E. zuernii by sequencing the complete mitogenome and analyzing the genetic diversity and population structure of 51 isolates collected from eight yak breeding parks in China.

RESULTS

The 6176-bp mitogenome of E. zuernii was linear and encoded typical mitochondrial contents of apicomplexan parasites, including three protein-coding genes [PCGs; cytochrome c oxidase subunits I and III (cox1 and cox3), and cytochrome b (cytb)], seven fragmented small subunit (SSU) and 12 fragmented large subunit (LSU) rRNAs. Genome-wide comparative and evolutionary analyses showed cytb and cox3 to be the most and least conserved Eimeria PCGs, respectively, and placed E. zuernii more closely related to Eimeria mephitidis than other Eimeria species. Furthermore, cox1-based genetic structure defined 24 haplotypes of E. zuernii with high haplotype diversities and low nucleotide diversities across eight geographic populations, supporting a low genetic structure and rapid evolutionary rate as well as a previous expansion event among E. zuernii populations.

CONCLUSIONS

To our knowledge, this is the first study presenting the phylogeny, genetic diversity, and population structure of the yak E. zuernii, and such information, together with its mitogenomic data, should contribute to a better understanding of the genetic and evolutionary biological studies of apicomplexan parasites in bovines.

The maintenance of genetic diversity under host-parasite coevolution in finite, structured populations

As a corollary to the Red Queen hypothesis, host-parasite coevolution has been hypothesized to maintain genetic variation in both species. Recent theoretical work, however, suggests that reciprocal natural selection alone is insufficient to maintain variation at individual loci. As highlighted by our brief review of the theoretical literature, models of host-parasite coevolution often vary along multiple axes (e.g. inclusion of ecological feedbacks or abiotic selection mosaics), complicating a comprehensive understanding of the effects of interacting evolutionary processes on diversity. Here we develop a series of comparable models to explore the effect of interactions between spatial structures and antagonistic coevolution on genetic diversity. Using a matching alleles model in finite populations connected by migration, we find that, in contrast to panmictic populations, coevolution in a spatially structured environment can maintain genetic variation relative to neutral expectations with migration alone. These results demonstrate that geographic structure is essential for understanding the effect of coevolution on biological diversity.

Host body size, not host population size, predicts genome-wide effective population size of parasites

The effective population size (N_e) of an organism is expected to be generally proportional to the total number of individuals in a population. In parasites, we might expect the effective population size to be proportional to host population size and host body size, because both are expected to increase the number of parasite individuals. However, among other factors, parasite populations are sometimes so extremely subdivided that high levels of inbreeding may distort these predicted relationships. Here, we used whole-genome sequence data from dove parasites (71 feather louse species of the genus Columbicola) and phylogenetic comparative methods to study the relationship between parasite effective population size and host population size and body size. We found that parasite effective population size is largely explained by host body size but not host population size. These results suggest the potential local population size (infrapopulation or deme size) is more predictive of the long-term effective population size of parasites than is the total number of potential parasite infrapopulations (i.e., host individuals).

Morphometrical Identification and Phylogenetic Analysis of Rhinonyssidae (Acari: Mesostigmata) Parasitizing Avian Hosts: New Molecular Data

Members of the family Rhinonyssidae are tiny hematophagous endoparasitic mites that inhabit the nasal cavities of birds and can cause trauma to their hosts. Traditionally, identifying species in this group has relied on observing their morphometrical characteristics. Nevertheless, determining species within this particular group has become more challenging due to the rising number of newly discovered species. Moreover, the morphometrical traits vary depending on the specific genus or group of species being studied. In this study, the complete internal transcribed spacer ITS1, 5.8S rDNA, and ITS2 regions of the ribosomal DNA from eighteen species of rhinonyssid mites belonging to four genera were sequenced to assess the utility of this genomic region in resolving taxonomic questions in this group and to estimate the phylogenetic relationships among the species. Mites were collected by dissecting the nasal cavities of birds under a stereomicroscope. Specimens used for morphometrical analyses were cleared in 85% lactic acid for 1-48 h and mounted in Hoyer's medium. Other specimens were preserved at -20 °C for molecular studies. From the data obtained in this study, it can be concluded that a thorough review and an accurate morphometrical identification and determination of the discriminatory traits are needed in this group of mites. Moreover, although the ITS1-5.8S-ITS2 fragment solves different taxonomic and phylogenetic problems at the species level, it would be necessary to test new molecular markers, or even a combination of nuclear and mitochondrial markers or different domains of the nuclear 28S rDNA, to discover a reliable taxonomic situation for rhinonyssids.

All quiet on the western front? The evolutionary history of monogeneans (Dactylogyridae: Cichlidogyrus, Onchobdella) infecting a West and Central African tribe of cichlid fishes (Chromidotilapiini)

Owing to the largely unexplored diversity of metazoan parasites, their speciation mechanisms and the circumstances under which such speciation occurs - in allopatry or sympatry - remain vastly understudied. Cichlids and their monogenean flatworm parasites have previously served as a study system for macroevolutionary processes, e.g., for the role of East African host radiations on parasite communities. Here, we investigate the diversity and evolution of the poorly explored monogeneans infecting a West and Central African lineage of cichlid fishes: Chromidotilapiini, which is the most species-rich tribe of cichlids in this region. We screened gills of 149 host specimens (27 species) from natural history collections and measured systematically informative characters of the sclerotised attachment and reproductive organs of the parasites. Ten monogenean species (Dactylogyridae: Cichlidogyrus and Onchobdella) were found, eight of which are newly described and one redescribed herein. The phylogenetic positions of chromidotilapiines-infecting species of Cichlidogyrus were inferred through a parsimony analysis of the morphological characters. Furthermore, we employed machine learning algorithms to detect morphological features associated with the main lineages of Cichlidogyrus. Although the results of these experimental algorithms remain inconclusive, the parsimony analysis indicates that West and Central African lineages of Cichlidogyrus and Onchobdella are monophyletic, unlike the paraphyletic host lineages. Several instances of host sharing suggest occurrences of intra-host speciation (sympatry) and host switching (allopatry). Some morphological variation was recorded that may also indicate the presence of species complexes. We conclude that collection material can provide important insights on parasite evolution despite the lack of well-preserved DNA material.

Prevalence of Borrelia burgdorferi and diversity of its outer surface protein C (ospC) alleles in blacklegged ticks (Ixodes scapularis) in Delaware

Characterizing the diversity of genes associated with virulence and transmission of a pathogen across the pathogen's distribution can inform our understanding of host infection risk. Borrelia burgdorferi is a vector-borne bacterium that causes Lyme disease in humans and is common in the United States. The outer surface protein C (ospC) gene of B. burgdorferi exhibits substantial genetic variation across the pathogen's distribution and plays a critical role in virulence and transmission in vertebrate hosts. In fact, B. burgdorferi infections that disseminate across host tissues in humans are associated with only a subset of ospC alleles. Delaware has a high incidence of Lyme disease, but the diversity of ospC in B. burgdorferi in the state has not been evaluated. We used PCR to amplify ospC in B. burgdorferi-infected blacklegged ticks (Ixodes scapularis) in sites statewide and used short-read sequencing to identify ospC alleles. B. burgdorferi prevalence in blacklegged ticks varied across sites, but not significantly so. We identified 15 previously characterized ospC alleles accounting for nearly all of the expected diversity of alleles across the sites as estimated using the Chao1 index. Nearly 40% of sequenced infections (23/58) had more than one ospC allele present suggesting mixed strain infections and the relative frequencies of alleles in single infections were positively correlated with their relative frequencies in mixed infections. Turnover of ospC alleles was positively related to distance between sites with closer sites having more similar allele compositions than more distant sites. This suggests a degree of B. burgdorferi dispersal limitation or habitat specialization. OspC alleles known to cause disseminated infections in humans were found at the highest frequencies across sites, corresponding to Delaware's high incidence of Lyme disease.

A new species of Hedruris Nitzsch, 1821 (Nematoda: Hedruridae) parasitic in eel Synbranchus marmoratus Bloch (Synbranchiformes: Synbranchidae) from Argentina, Neotropical Region

During a parasitological survey, specimens belonging to the genus Hedruris (Nematoda: Hedruridae) were found in the stomach of Synbranchus marmoratus Bloch from the Paraná Guazú River, Paraná River basin (Argentina). Comparative morphological studies with its congeners showed that these nematodes belonged to a new species, herein described as Hedruris anguila sp. nov. The new species can be distinguished from remaining Hedruris species by the pseudolabia size, the morphology and size of spicules, the length of female caudal hook, the shape and size of eggs and the particular morphology of the male tail. Hedruris anguila sp. nov. closely resembles Hedruris bifida Rossin & Timi, 2016 and H. suttonae Brugni & Viozzi, 2010, both also described from fish hosts, in the arrangement of caudal papillae and morphology of eggs, but differ in the size of female hook and spicules, and the number of caudal papillae in males. The main morphological features, host range and geographical distribution of all species in the genus Hedruris species are analyzed and discussed. This is the first report of a representative of the family Hedruridae parasitizing a species of Synbranchidae and the eleventh report of the genus Hedruris for a fish host. The number of valid species of Hedruris described was increased to 26, including the new species; 11 for the Neotropical Region and represent the third species found in fish host from Argentina. Based on morphological analysis, the speciation process probably occurred first in Gondwana and the non-mammillated eggs could be an ancestral morphological character. The diet of Synbranchus marmoratus and values of prevalence and mean intensity could indicate that this host species is a true definitive host for this nematode.  

Population genomics of helminth parasites

Next generation sequencing technologies have facilitated a shift from a few targeted loci in population genetic studies to whole genome approaches. Here, we review the types of questions and inferences regarding the population biology and evolution of parasitic helminths being addressed within the field of population genomics. Topics include parabiome, hybridization, population structure, loci under selection and linkage mapping. We highlight various advances, and note the current trends in the field, particularly a focus on human-related parasites despite the inherent biodiversity of helminth species. We conclude by advocating for a broader application of population genomics to reflect the taxonomic and life history breadth displayed by helminth parasites. As such, our basic knowledge about helminth population biology and evolution would be enhanced while the diversity of helminths in itself would facilitate population genomic comparative studies to address broader ecological and evolutionary concepts.

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.

Low genetic heterogeneity of Leishmania major in different geographical regions of Iran

To examine the genetic diversity of Leishmania major, 100 Giemsa-stained positive slides were collected from endemic foci of Iran (Northeast, Central, and Southwest provinces) over two consecutive years during 2019-2021. The Leishmania ITS-rDNA gene was amplified and Leishmania sp. was recognized by PCR-RFLP and sequencing. In addition, 178 registered ITS-rDNA sequences from other geographical regions of Iran were retrieved from GenBank, including different host species (human, sandfly and rodent). A total of 40 new haplotypes were discovered using the ITS-rDNA sequence analysis. IR29 (20.6%) and IR34 (61%) were the two most common haplotypes, represented by a star-like feature in the overall population. Analysis of the molecular variance test revealed low genetic diversity of L. major in human cases (Haplotype diversity; 0.341), rodent (Hd; 0.387) and sandfly (Hd; 0.390) sequences. The lowest genetic diversity of L. major was observed in Southwest/Southeast Iran (Hd: 0.104-0.286). The statistically Fst value indicated that L. major is not genetically differentiated between geographic regions of Iran, except for the Northeast-Southwest (Fst: 0.29055) and Central-Southwest (Fst: 0.30294) population pairs. The current study as the first investigation discloses new perspectives for further evaluation in the identification local transmission paradigms and initiating effective prevention strategies.

Traveler Mites: Population Genetic Structure of the Wing Mites Periglischrus paracaligus (Acari: Mesostigmata: Spinturnicidae)

Wing mites of the genus Periglischrus are ectoparasites exclusively associated with phyllostomid bats. These mites show high host specificity and have been studied to understand the evolutionary history of their bat hosts mainly by using a morphological variation. Through a phylogeographic approach, we analyzed the genetic diversity and population genetic structure of the ectoparasite Periglischrus paracaligus Herrin and Tipton which parasitizes Leptonycteris yerbabuenae Martínez and Villa (lesser long-nosed bat) in Mexico. By the implementation of a multilocus approach, we found that P. paracaligus populations were diverse for haplotype diversity, and had values ranging from 0.5 to 1. No genetic structuring in the P. paracaligus parasites was observed along with the distribution of the host, L. yerbabuenae, in Mexico, nor when populations or regions were compared, but our results revealed a process of historical demographic expansion in all the analyzed markers. We discuss possible scenarios that could explain the lack of population structure in the light of the data analyzed for the parasites and the biology of L. yerbabuenae, such as the interplay between parasite and host traits being responsible for the genetic make-up of parasite populations. We also inferred its phylogenetic position among wing mites parasitizing the two other species of Leptonycteris bats. Long-nosed bats' monophyly helps to explain the observed presence of distinctive clades in the wing mite's phylogeny in specific association with each long-nosed bat host species.

Phylogenomics and Diversification of the Schistosomatidae Based on Targeted Sequence Capture of Ultra-Conserved Elements

Schistosomatidae Stiles and Hassall 1898 is a medically significant family of digenetic trematodes (Trematoda: Digenea), members of which infect mammals or birds as definitive hosts and aquatic or amphibious gastropods as intermediate hosts. Currently, there are 17 named genera, for many of which evolutionary interrelationships remain unresolved. The lack of a resolved phylogeny has encumbered our understanding of schistosomatid evolution, specifically patterns of host-use and the role of host-switching in diversification. Here, we used targeted sequence capture of ultra-conserved elements (UCEs) from representatives of 13 of the 17 named genera and 11 undescribed lineages that are presumed to represent either novel genera or species to generate a phylogenomic dataset for the estimation of schistosomatid interrelationships. This study represents the largest phylogenetic effort within the Schistosomatidae in both the number of loci and breadth of taxon sampling. We present a near-comprehensive family-level phylogeny providing resolution to several clades of long-standing uncertainty within Schistosomatidae, including resolution for the placement of the North American mammalian schistosomes, implying a second separate capture of mammalian hosts. Additionally, we present evidence for the placement of Macrobilharzia at the base of the Schistosoma + Bivitellobilharzia radiation. Patterns of definitive and intermediate host use and a strong role for intermediate host-switching are discussed relative to schistosomatid diversification.

High levels of inbreeding with spatial and host-associated structure in lice of an endangered freshwater seal

Host-specialist parasites of endangered large vertebrates are in many cases more endangered than their hosts. In particular, low host population densities and reduced among-host transmission rates are expected to lead to inbreeding within parasite infrapopulations living on single host individuals. Furthermore, spatial population structures of directly-transmitted parasites should be concordant with those of their hosts. Using population genomic approaches, we investigated inbreeding and population structure in a host-specialist seal louse (Echinophthirius horridus) infesting the Saimaa ringed seal (Phoca hispida saimensis), which is endemic to Lake Saimaa in Finland, and is one of the most endangered pinnipeds in the world. We conducted genome resequencing of pairs of lice collected from 18 individual Saimaa ringed seals throughout the Lake Saimaa complex. Our analyses showed high genetic similarity and inbreeding between lice inhabiting the same individual seal host, indicating low among-host transmission rates. Across the lake, genetic differentiation among individual lice was correlated with their geographic distance, and assignment analyses revealed a marked break in the genetic variation of the lice in the middle of the lake, indicating substantial population structure. These findings indicate that movements of Saimaa ringed seals across the main breeding areas of the fragmented Lake Saimaa complex may in fact be more restricted than suggested by previous population-genetic analyses of the seals themselves.

Population genetic structure and phenotypic diversity of Aspidodera raillieti (Nematoda: Heterakoidea), a parasite of Didelphini marsupials in Brazil's South and Southeast Atlantic Forest

Background

The population genetics of parasites may be influenced by host specificity, life cycle, host geographical range, evolutionary history, and host population structure. The nematode Aspidodera raillieti infects different marsupial and rodent hosts in the Nearctic and Neotropical regions, implying a gene flow among populations. However, niche diversification of the main hosts of A. raillieti in superimposed areas may provide conditions for population genetic structuring within this parasite species. We examined the genetic structuring of A. raillieti infecting three marsupial species co-occurring along the South and Southeast Brazilian Atlantic Forest, a hotspot of biodiversity.

Methods

We employed morphometric analyses and partial mitochondrial cytochrome c oxidase I gene sequences (MT-CO1) to characterize populations via phylogenetic and phylogeographic analyses.

Results

Among 175 A. raillieti specimens recovered from the marsupial hosts Didelphis aurita, D. albiventris, and Philander quica, we identified 99 MT-CO1 haplotypes forming four haplogroups and four clades in networks and phylogenetic trees, respectively. Clades I and II encompassed parasites of D. albiventris from the South region, clade III comprised parasites of D. aurita from the South and Southeast regions, and clade IV encompassed parasites of D. aurita and D. albiventris from the South and Southeast regions and parasites of P. quica from the South region. High genetic differentiation between clades, with a high fixation index and greater genetic variation in the analysis of molecular variance (AMOVA), indicated low gene flow between clades. Haplotypes shared among host species revealed a lack of host specificity. A significant correlation in the Mantel test suggested parasite isolation by distance, while there was no evidence of geographical structure between populations. Negative neutrality test values for clades III and IV suggested recent population expansion. Morphometric differentiation between A. raillieti specimens recovered from different host species, as well as from different localities, was more evident in males.

Conclusion

The genetic structure of A. raillieti populations in the South and Southeast Atlantic Forest resulted from historical events rather than from current geographical distribution or host specificity. We also demonstrate morphometric variation associated with host species and localities, suggesting phenotypic plasticity to host attributes and to spatial variables.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05288-6.

Why Plasmodium vivax and Plasmodium falciparum are so different? A tale of two clades and their species diversities

The global malaria burden sometimes obscures that the genus Plasmodium comprises diverse clades with lineages that independently gave origin to the extant human parasites. Indeed, the differences between the human malaria parasites were highlighted in the classical taxonomy by dividing them into two subgenera, the subgenus Plasmodium, which included all the human parasites but Plasmodium falciparum that was placed in its separate subgenus, Laverania. Here, the evolution of Plasmodium in primates will be discussed in terms of their species diversity and some of their distinct phenotypes, putative molecular adaptations, and host–parasite biocenosis. Thus, in addition to a current phylogeny using genome-level data, some specific molecular features will be discussed as examples of how these parasites have diverged. The two subgenera of malaria parasites found in primates, Plasmodium and Laverania, reflect extant monophyletic groups that originated in Africa. However, the subgenus Plasmodium involves species in Southeast Asia that were likely the result of adaptive radiation. Such events led to the Plasmodium vivax lineage. Although the Laverania species, including P. falciparum, has been considered to share “avian characteristics,” molecular traits that were likely in the common ancestor of primate and avian parasites are sometimes kept in the Plasmodium subgenus while being lost in Laverania. Assessing how molecular traits in the primate malaria clades originated is a fundamental science problem that will likely provide new targets for interventions. However, given that the genus Plasmodium is paraphyletic (some descendant groups are in other genera), understanding the evolution of malaria parasites will benefit from studying “non-Plasmodium” Haemosporida.

Somewhere I belong: phylogeny and morphological evolution in a species-rich lineage of ectoparasitic flatworms infecting cichlid fishes

A substantial portion of biodiversity has evolved through adaptive radiation. However, the effects of explosive speciation on species interactions remain poorly understood. Metazoan parasites infecting radiating host lineages could improve our knowledge because of their intimate host relationships. Yet limited molecular, phenotypic and ecological data discourage multivariate analyses of evolutionary patterns and encourage the use of discrete characters. Here, we assemble new molecular, morphological and host range data widely inferred from a species-rich lineage of parasites (Cichlidogyrus, Platyhelminthes: Monogenea) infecting cichlid fishes to address data scarcity. We infer a multimarker (28S/18S rDNA, ITS1, COI mtDNA) phylogeny of 58 of 137 species and characterize major lineages through synapomorphies inferred from mapping morphological characters. We predict the phylogenetic position of species without DNA data through shared character states, a morphological phylogenetic analysis, and a classification analysis with support vector machines. Based on these predictions and a cluster analysis, we assess the systematic informativeness of continuous characters, search for continuous equivalents for discrete characters, and suggest new characters for morphological traits not analysed to date. We also model the attachment/reproductive organ and host range evolution using the data for 136 of 137 described species and multivariate phylogenetic comparative methods (PCMs). We show that discrete characters not only can mask phylogenetic signals, but also are key for characterizing species groups. Regarding the attachment organ morphology, a divergent evolutionary regime for at least one lineage was detected and a limited morphological variation indicates host and environmental parameters affecting its evolution. However, moderate success in predicting phylogenetic positions, and a low systematic informativeness and high multicollinearity of morphological characters call for a revaluation of characters included in species characterizations.

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.

Expansion of Cyclophyllidea Biodiversity in Rodents of Qinghai-Tibet Plateau and the "Out of Qinghai-Tibet Plateau" Hypothesis of Cyclophyllideans

The Cyclophyllidea comprises the most species-rich order of tapeworms (Platyhelminthes, Cestoda) and includes species with some of the most severe health impact on wildlife, livestock, and humans. We collected seven Cyclophyllidea specimens from rodents in Qinghai-Tibet Plateau (QTP) and its surrounding mountain systems, of which four specimens in QTP were unsequenced, representing “putative new species.” Their complete mitochondrial (mt) genomes were sequenced and annotated. Phylogenetic reconstruction of partial 28S rDNA, cox1 and nad1 datasets provided high bootstrap frequency support for the categorization of three “putative new species,” assigning each, respectively, to the genera Mesocestoides, Paranoplocephala, and Mosgovoyia, and revealing that some species and families in these three datasets, which contain 291 species from nine families, may require taxonomic revision. The partial 18S rDNA phylogeny of 29 species from Taeniidae provided high bootstrap frequency support for the categorization of the “putative new species” in the genus Hydatigera. Combined with the current investigation, the other three known Taeniidae species found in this study were Taenia caixuepengi, T. crassiceps, and Versteria mustelae and may be widely distributed in western China. Estimates of divergence time based on cox1 + nad1 fragment and mt protein-coding genes (PCGs) showed that the differentiation rate of Cyclophyllidea species was strongly associated with the rate of change in the biogeographic scenarios, likely caused by the uplift of the QTP; i.e., species differentiation of Cyclophyllidea might be driven by host-parasite co-evolution caused by the uplift of QTP. We propose an “out of QTP” hypothesis for the radiation of these cyclophyllidean tapeworms.

Patterns of deer ked (Diptera: Hippoboscidae) and tick (Ixodida: Ixodidae) infestation on white-tailed deer (Odocoileus virginianus) in the eastern United States

Background

White-tailed deer (Odocoileus virginianus) host numerous ectoparasitic species in the eastern USA, most notably various species of ticks and two species of deer keds. Several pathogens transmitted by ticks to humans and other animal hosts have also been found in deer keds. Little is known about the acquisition and potential for transmission of these pathogens by deer keds; however, tick-deer ked co-feeding transmission is one possible scenario. On-host localization of ticks and deer keds on white-tailed deer was evaluated across several geographical regions of the eastern US to define tick-deer ked spatial relationships on host deer, which may impact the vector-borne disease ecology of these ectoparasites.

Methods

Ticks and deer keds were collected from hunter-harvested white-tailed deer from six states in the eastern US. Each deer was divided into three body sections, and each section was checked for 4 person-minutes. Differences in ectoparasite counts across body sections and/or states were evaluated using a Bayesian generalized mixed model.

Results

A total of 168 white-tailed deer were inspected for ticks and deer keds across the study sites. Ticks (n = 1636) were collected from all surveyed states, with Ixodes scapularis (n = 1427) being the predominant species. Counts of I. scapularis from the head and front sections were greater than from the rear section. Neotropical deer keds (Lipoptena mazamae) from Alabama and Tennessee (n = 247) were more often found on the rear body section. European deer keds from Pennsylvania (all Lipoptena cervi, n = 314) were found on all body sections of deer.

Conclusions

The distributions of ticks and deer keds on white-tailed deer were significantly different from each other, providing the first evidence of possible on-host niche partitioning of ticks and two geographically distinct deer ked species (L. cervi in the northeast and L. mazamae in the southeast). These differences in spatial distributions may have implications for acquisition and/or transmission of vector-borne pathogens and therefore warrant further study over a wider geographic range and longer time frame.

Graphical Abstract

Gene migration of giardiasis in Iran; a microevolutionary scale for reflecting transmission patterns of Giardia lamblia assemblages in symptomatic patients

In the microevolutionary scale of Giardia lamblia, the gene migration indicates how G. lamblia assemblages have transmitted between adjacent counties. 33 positive fecal samples were taken from patients suffering gastrointestinal disorders (nausea, bloating, burping constipation and fatty diarrhea) at Tabriz and Ardabil cities, where located in the cold regions of northwest Iran. Following parasitological examinations, DNA samples were extracted, amplified and digested by single-step PCR-RFLP assay, targeting the glutamate dehydrogenase (gdh) locus to distinguish within and between assemblages A and B. PCR products were directly sequenced to reconfirm their heterogeneity traits and phylogenetic analysis. Of the 33 isolates, 81.9% (n: 27), 9% (n: 3) and 9% (n: 3) were successfully identified as assemblages A (genotype AII), B (genotype BIII) and the mixed of genotypes AII and B, respectively. Despite the presence of heterogeneous clinical backgrounds, a low genetic diversity of sub-assemblage AII was identified among symptomatic cases. A low value of pairwise fixation index showed that G. lamblia sub-assemblage AII is not genetically differentiated among northwest regions of Iran. The occurrence of haplotypes TAB-1/ARD-1 between two regional populations indicates that there is a dawn of G. lamblia gene flow due to transfer of alleles through host mobility and/or ecological alterations. To assess the hypothetical evolutionary scenario, further studies are essential for multilocus genotyping of G. lamblia in tropical regions of Iran and neighboring countries.

Comparative mitogenomics of Spirocerca lupi from South Africa and China: Variation and possible heteroplasmy

The complete mitochondrial genome of Spirocerca lupi isolated from a dog in South Africa was sequenced using next generation sequencing (NGS) technology and the 12 protein coding genes along with the two rRNA genes were compared to 18 other nematode species as well as S. lupi from China. The mitochondrial genome of S. lupi South Africa had a mean genetic diversity of 6.1 % compared to S. lupi China with some variation in nucleotide composition, gene positioning and size. Pairwise distance results indicated slightly higher variation when compared to the pairwise distances of other closely related species, however, this variation was not high enough for it to be considered a cryptic species. Phylogenetic analysis indicated that S. lupi from the two continents are very similar. In addition, single nucleotide polymorphisms were detected in the nad2 gene with ten sequence variants identified from 10 clones from a single nematode, suggesting possible heteroplasmy. The origin of the heteroplasmy is currently unknown but it is speculated to have arisen from accumulated mutations in the mitochondria during somatic replication.

Slow crabs - fast genomes: Locomotory capacity predicts skew magnitude in crustacean mitogenomes

Base composition skews (G-C/G+C) of mitochondrial genomes are believed to be primarily driven by mutational pressure, which is positively correlated with metabolic rate. In marine animals, metabolic rate is also positively correlated with locomotory capacity. Given the central role of mitochondria in energy metabolism, we hypothesised that selection for locomotory capacity should be positively correlated with the strength of purifying selection (dN/dS), and thus be negatively correlated with the skew magnitude. Therefore, these two models assume diametrically opposite associations between the metabolic rate and skew magnitude: positive correlation in the prevailing paradigm, and negative in our working hypothesis. We examined correlations between the skew magnitude, metabolic rate, locomotory capacity, and several other variables previously associated with mitochondrial evolution on 287 crustacean mitogenomes. Weakly locomotory taxa had higher skew magnitude and ω (dN/dS) values, but not the gene order rearrangement rate. Skew and ω magnitudes were correlated. Multilevel regression analyses indicated that three competing variables, body size, gene order rearrangement rate, and effective population size, had negligible impacts on the skew magnitude. In most crustacean lineages selection for locomotory capacity appears to be the primary factor determining the skew magnitude. Contrary to the prevailing paradigm, this implies that adaptive selection outweighs nonadaptive selection (mutation pressure) in crustaceans. However, we found indications that effective population size (nonadaptive factor) may outweigh the impact of locomotory capacity in sessile crustaceans (Thecostraca). In conclusion, skew magnitude is a product of the interplay between adaptive and nonadaptive factors, the balance of which varies among lineages.

Comparative population genetics of swimming crab host (Portunus pelagicus) and common symbiotic barnacle (Octolasmis angulata) in Vietnam

Background

By comparing spatial geographical structures of host populations with that of their symbionts light can be shed on their biological interactions, and the degree of congruence between host and symbiont phylogeographies should reflect their life histories and especially dispersal mechanisms.

Methods

Here, we analyzed the genetic diversity and structure of a host, the blue swimming crab, Portunus pelagicus, and its symbiotic pedunculate barnacle Octolasmis angulata from six location sites representing three geographic regions (north, central and south) along the Vietnam coastline. High levels of congruence in their phylogeographic patterns were expected as they both undergo planktonic larval stages.

Results

Based on the COI mtDNA markers, O. angulata populations showed higher genetic diversity in comparison with their host P. pelagicus (number of haplotype/individuals, haplotype and nucleotide diversity are 119/192, 0.991 ± 0.002 and 0.02; and 89/160, 0.913 ± 0.02 and 0.015, respectively). Pairwise Fst and AMOVA analyses showed a more pronounced population structure in the symbiotic barnacle than in its crab host. The DAPC analyses identified three genetic clusters. However, both haplotype networks and scatter plots supported connectivity of the host and the symbiotic barnacle throughout their distribution range, except for low subdivision of southern population. Isolation by distance were detected only for the symbiont O. angulata (R² = 0.332, P = 0.05), while dbMEM supported spatial structure of both partners, but only at MEM-1 (Obs. 0.2686, P < 0.01 and Obs. 0.2096, P < 0.01, respectively).

Association of putatively adaptive genetic variation with climatic variables differs between a parasite and its host

Parasitism is a pervasive phenomenon in nature with the relationship between species driving evolution in both parasite and host. Due to their host-dependent lifestyle, parasites may adapt to the abiotic environment in ways that differ from their hosts or from free-living relatives; yet rarely has this been assessed. Here, we test two competing hypotheses related to whether putatively adaptive genetic variation in a specialist mistletoe associates with the same, or different, climatic variables as its host species. We sampled 11 populations of the specialist mistletoe Amyema gibberula var. tatei (n = 154) and 10 populations of its associated host Hakea recurva subsp. recurva (n = 160). Reduced-representation sequencing was used to obtain genome-wide markers and putatively adaptive variation detected using genome scan methods. Climate associations were identified using generalized dissimilarity modelling, and these were mapped geographically to visualize the spatial patterns of genetic composition. Our results supported the hypothesis of parasites and host species responding differently to climatic variables. Temperature was relatively more important in predicting allelic turnover in the specialist mistletoe while precipitation was more important for the host. This suggests that parasitic plants and host species may respond differently to selective pressures, potentially as a result of differing nutrient acquisition strategies. Specifically, mistletoes acquire water from hosts (rather than the abiotic environment), which may provide a buffer to precipitation as a selective pressure. This work deepens and complements the physiological and other ecological studies of adaptation and provides a window into the evolutionary processes that underlie previously observed phenomena. Applying these methods to a comparative study in a host-parasite system has also highlighted factors that affect the study of selection pressure on nonmodel organisms, such as differing adaptation rates and lack of reference genomes.

Host tropism determination by convergent evolution of immunological evasion in the Lyme disease system

Pathogens possess the ability to adapt and survive in some host species but not in others-an ecological trait known as host tropism. Transmitted through ticks and carried mainly by mammals and birds, the Lyme disease (LD) bacterium is a well-suited model to study such tropism. Three main causative agents of LD, Borrelia burgdorferi, B. afzelii, and B. garinii, vary in host ranges through mechanisms eluding characterization. By feeding ticks infected with different Borrelia species, utilizing feeding chambers and live mice and quail, we found species-level differences in bacterial transmission. These differences localize on the tick blood meal, and specifically complement, a defense in vertebrate blood, and a polymorphic bacterial protein, CspA, which inactivates complement by binding to a host complement inhibitor, Factor H (FH). CspA selectively confers bacterial transmission to vertebrates that produce FH capable of allele-specific recognition. CspA is the only member of the Pfam54 gene family to exhibit host-specific FH-binding. Phylogenetic analyses revealed convergent evolution as the driver of such uniqueness, and that FH-binding likely emerged during the last glacial maximum. Our results identify a determinant of host tropism in Lyme disease infection, thus defining an evolutionary mechanism that shapes host-pathogen associations.

Genome reduction and relaxed selection is associated with the transition to symbiosis in the basidiomycete genus Podaxis

Insights into the genomic consequences of symbiosis for basidiomycete fungi associated with social insects remain sparse. Capitalizing on viability of spores from centuries-old herbarium specimens of free-living, facultative, and specialist termite-associated Podaxis fungi, we obtained genomes of 10 specimens, including two type species described by Linnaeus >240 years ago. We document that the transition to termite association was accompanied by significant reductions in genome size and gene content, accelerated evolution in protein-coding genes, and reduced functional capacities for oxidative stress responses and lignin degradation. Functional testing confirmed that termite specialists perform worse under oxidative stress, while all lineages retained some capacity to cleave lignin. Mitochondrial genomes of termite associates were significantly larger; possibly driven by smaller population sizes or reduced competition, supported by apparent loss of certain biosynthetic gene clusters. Our findings point to relaxed selection that mirrors genome traits observed among obligate endosymbiotic bacteria of many insects.

Molecular approaches reveal weak sibship aggregation and a high dispersal propensity in a non-native fish parasite

Inferring parameters related to the aggregation pattern of parasites and to their dispersal propensity are important for predicting their ecological consequences and evolutionary potential. Nonetheless, it is notoriously difficult to infer these parameters from wildlife parasites given the difficulty in tracking these organisms. Molecular-based inferences constitute a promising approach that has yet rarely been applied in the wild. Here, we combined several population genetic analyses including sibship reconstruction to document the genetic structure, patterns of sibship aggregation, and the dispersal dynamics of a non-native parasite of fish, the freshwater copepod ectoparasite Tracheliastes polycolpus. We collected parasites according to a hierarchical sampling design, with the sampling of all parasites from all host individuals captured in eight sites spread along an upstream-downstream river gradient. Individual multilocus genotypes were obtained from 14 microsatellite markers, and used to assign parasites to full-sib families and to investigate the genetic structure of T. polycolpus among both hosts and sampling sites. The distribution of full-sibs obtained among the sampling sites was used to estimate individual dispersal distances within families. Our results showed that T. polycolpus sibs tend to be aggregated within sites but not within host individuals. We detected important upstream-to-downstream dispersal events of T. polycolpus between sites (modal distance: 25.4 km; 95% CI [22.9, 27.7]), becoming scarcer as the geographic distance from their family core location increases. Such a dispersal pattern likely contributes to the strong isolation-by-distance observed at the river scale. We also detected some downstream-to-upstream dispersal events (modal distance: 2.6 km; 95% CI [2.2-23.3]) that likely result from movements of infected hosts. Within each site, the dispersal of free-living infective larvae among hosts likely contributes to increasing genetic diversity on hosts, possibly fostering the evolutionary potential of T. polycolpus.

Genetic Evolution and Implications of the Mitochondrial Genomes of Two Newly Identified Taenia spp. in Rodents From Qinghai-Tibet Plateau

The larva of Taeniidae species can infect a wide range of mammals, causing major public health and food safety hazards worldwide. The Qinghai-Tibet Plateau (QTP), a biodiversity hotspot, is home to many species of rodents, which act as the critical intermediate hosts of many Taeniidae species. In this study, we identified two new larvae of Taenia spp., named T. caixuepengi and T. tianguangfui, collected from the plateau pika (Ochotona curzoniae) and the Qinghai vole (Neodon fuscus), respectively, in QTP, and their mitochondrial genomes were sequenced and annotated. Phylogenetic trees based on the mitochondrial genome showed that T. caixuepengi has the closest genetic relationship with T. pisiformis, while T. tianguangfui was contained in a monophyletic group with T. crassiceps, T. twitchelli, and T. martis. Biogeographic scenarios analysis based on split time speculated that the speciation of T. caixuepengi (∼5.49 Mya) is due to host switching caused by the evolution of its intermediate host. Although the reason for T. tianguangfui (∼13.11 Mya) speciation is not clear, the analysis suggests that it should be infective to a variety of other rodents following the evolutionary divergence time of its intermediate host and the range of intermediate hosts of its genetically close species. This study confirms the species diversity of Taeniidae in the QTP, and speculates that the uplift of the QTP has not only a profound impact on the biodiversity of plants and animals, but also that of parasites.

Elucidation of Himasthla leptosoma (Creplin, 1829) Dietz, 1909 (Digenea, Himasthlidae) life cycle with insights into species composition of the north Atlantic Himasthla associated with periwinkles Littorina spp

Trematodes of the genus Himasthla are usual parasites of coastal birds in nearshore ecosystems of northern European seas and the Atlantic coast of North America. Their first intermediate hosts are marine and brackish-water gastropods, while second intermediate hosts are various invertebrates. We analysed sequences of partial 28S rRNA and nad1 genes and the morphology of intramolluscan stages, particularly cercariae of Himasthla spp. parasitizing intertidal molluscs Littorina spp. in the White Sea, the Barents Sea and coasts of North Norway and Iceland. We showed that only three Himasthla spp. are associated with periwinkles in these regions. Intramolluscan stages of H. elongata were found in Littorina littorea, of H. littorinae, in both L. saxatilis and L. obtusata, and of Cercaria littorinae obtusatae, predominantly, in L. obtusata. Other Himasthla spp. previously reported from Littorina spp. in North Atlantic are either synonymous with one of these species or described erroneously. Based on a comparison of newly generated 28S rDNA sequences with GenBank data, rediae and cercariae of C. littorinae obtusatae were identified as belonging to H. leptosoma. Some previously unknown morphological features of young and mature rediae and cercariae of the three Himasthla spp. are described. We provide a key to the rediae and highlight characters important for identification of cercariae. Genetic diversity within the studied species was only partially determined by their specificity to the molluscan host. The nad1 network constructed for H. leptosoma lacked geographical structure, which is explained by a high gene flow owing to highly vagile definitive hosts, shorebirds.