Host-specific Adaptation Research Articles

Genomic revolution in parasitology: Unlocking the mysteries of parasites

Whole-genome sequencing has revolutionized parasitology, offering insights into parasites' biology, evolution, and potential interventions. This discussion explores the profound impact of genomics on parasitology, from democratizing research to understanding parasites' intricacies and their relationships with hosts. High-throughput sequencing has made genomics accessible globally, allowing labs worldwide to study parasitic organisms. Whole-genome sequencing provides a comprehensive view of parasite genes and regulatory elements, facilitating in-depth biological investigations. It helps identify drug targets and reveals the genetic basis of drug resistance, which is vital for new antiparasitic drug development. Comparative genomics is a powerful tool for understanding the evolutionary relationships among parasites. Researchers create robust phylogenetic trees that clarify parasite evolution by comparing genomic data. Case studies of malaria parasites, trypanosomes, schistosomes, and Leishmania species illustrate the insights gained. Genomics also sheds light on parasite life cycles, uncovering gene expression, regulatory networks, key proteins, and molecular markers. This informs our understanding of host-specific adaptations. Practical applications encompass drug target discovery, vaccine development, diagnostics, intervention strategies, understanding drug resistance, and disease eradication efforts. In conclusion, genomics is transforming parasitology, especially whole-genome sequencing and comparative genomics. It offers a holistic view of parasites and their evolution, paving the way for innovative antiparasitic strategies and contributing to the control and potential eradication of parasitic diseases.

Comprehensive profile of the companion animal gut microbiome integrating reference-based and reference-free methods.

The gut microbiome of companion animals is relatively underexplored, despite its relevance to animal health, pet owner health, and basic microbial community biology. Here, we provide the most comprehensive analysis of the canine and feline gut microbiomes to date, incorporating 2639 stool shotgun metagenomes (2272 dog and 367 cat) spanning 14 publicly available datasets (n = 730) and 8 new study populations (n = 1909). These are compared with 238 and 112 baseline human gut metagenomes from the Human Microbiome Project 1-II and a traditionally living Malagasy cohort, respectively, processed in a manner identical to the animal metagenomes. All microbiomes were characterized using reference-based taxonomic and functional profiling, as well as de novo assembly yielding metagenomic assembled genomes clustered into species-level genome bins. Companion animals shared 184 species-level genome bins not found in humans, whereas 198 were found in all three hosts. We applied novel methodology to distinguish strains of these shared organisms either transferred or unique to host species, with phylogenetic patterns suggesting host-specific adaptation of microbial lineages. This corresponded with functional divergence of these lineages by host (e.g., differences in metabolic and antibiotic resistance genes) likely important to companion animal health. This study provides the largest resource to date of companion animal gut metagenomes and greatly contributes to our understanding of the "One Health" concept of a shared microbial environment among humans and companion animals, affecting infectious diseases, immune response, and specific genetic elements.

Pathogenic and Comparative Genomic Analysis of Ralstonia pseudosolanacearum Isolated from Casuarina.

Casuarina equisetifolia is crucial in protecting coastal regions of China against typhoon attacks but has faced a substantial challenge due to wilt disease caused by pathogens of the Ralstonia solanacearum species complex (RSSC). Although the initial outbreak of Casuarina wilt in the 1970s was effectively controlled by disease-resistant C. equisetifolia varieties, the disease has recently re-emerged in coastal regions of Guangdong. In this study, we report the isolation, characterization, and comparative genomic analysis of 11 RSSC strains from diseased C. equisetifolia at various locations along the coast of Guangdong. Phylogenomic analysis showed that the strains were closely related and clustered with phylotype I strains previously isolated from peanuts. Single-gene-based analysis further suggested these strains could be derived from strains present in Guangdong since the 1980s, indicating a historical context to their current pathogenicity. Casuarina-isolated strains exhibited notably higher virulence against C. equisetifolia and peanuts than the representative RSSC strains GMI1000 and EP1, suggesting host-specific adaptations that possibly contributed to the recent outbreak. Comparative genomic analysis among RSSC strains revealed a largely conserved genome structure and high levels of conservation in gene clusters encoding extracellular polysaccharide biosynthesis, secretion systems, and quorum sensing regulatory systems. However, we also found a number of unique genes in the Casuarina-isolated strains that were absent in GMI1000 and EP1, and vice versa, pointing to potential genetic factors underpinning their differential virulence. These unique genes offer promising targets for future functional studies. Overall, our findings provide crucial insights into the RSSC pathogens causing Casuarina wilt in Guangdong, guiding future efforts in disease control and prevention.

Evolution and host-specific adaptation of Pseudomonas aeruginosa.

The major human bacterial pathogen Pseudomonas aeruginosa causes multidrug-resistant infections in people with underlying immunodeficiencies or structural lung diseases such as cystic fibrosis (CF). We show that a few environmental isolates, driven by horizontal gene acquisition, have become dominant epidemic clones that have sequentially emerged and spread through global transmission networks over the past 200 years. These clones demonstrate varying intrinsic propensities for infecting CF or non-CF individuals (linked to specific transcriptional changes enabling survival within macrophages); have undergone multiple rounds of convergent, host-specific adaptation; and have eventually lost their ability to transmit between different patient groups. Our findings thus explain the pathogenic evolution of P. aeruginosa and highlight the importance of global surveillance and cross-infection prevention in averting the emergence of future epidemic clones.

Predicting hosts and cross-species transmission of Streptococcus agalactiae by interpretable machine learning

BackgroundStreptococcus agalactiae, commonly known as Group B Streptococcus (GBS), exhibits a broad host range, manifesting as both a beneficial commensal and an opportunistic pathogen across various species. In humans, it poses significant risks, causing neonatal sepsis and meningitis, along with severe infections in adults. Additionally, it impacts livestock by inducing mastitis in bovines and contributing to epidemic mortality in fish populations. Despite its wide host spectrum, the mechanisms enabling GBS to adapt to specific hosts remain inadequately elucidated. Therefore, the development of a rapid and accurate method differentiates GBS strains associated with particular animal hosts based on genome-wide information holds immense potential. Such a tool would not only bolster the identification and containment efforts during GBS outbreaks but also deepen our comprehension of the bacteria's host adaptations spanning humans, livestock, and other natural animal reservoirs. Methods and resultsHere, we developed three machine learning models—random forest (RF), logistic regression (LR), and support vector machine (SVM) based on genome-wide mutation data. These models enabled precise prediction of the host origin of GBS, accurately distinguishing between human, bovine, fish, and pig hosts. Moreover, we conducted an interpretable machine learning using SHapley Additive exPlanations (SHAP) and variant annotation to uncover the most influential genomic features and associated genes for each host. Additionally, by meticulously examining misclassified samples, we gained valuable insights into the dynamics of host transmission and the potential for zoonotic infections. ConclusionsOur study underscores the effectiveness of random forest (RF) and logistic regression (LR) models based on mutation data for accurately predicting GBS host origins. Additionally, we identify the key features associated with each GBS host, thereby enhancing our understanding of the bacteria's host-specific adaptations.

Pins Gene Table v2.0: An Online Genome Database of 37 Pythium insidiosum Strains for Gene Content Exploration and Phylogenomic Analysis.

Unlike most pathogenic oomycetes, Pythium insidiosum infects humans and animals instead of plants. P. insidiosum has three clinically relevant genotypes/clades that cause a severe disease called pythiosis. To develop strategies for infection control, it is necessary to understand the biology and pathogenesis of this pathogen. Investigating the evolutionary mechanisms behind the host-specific adaptation is vital, and comparative genomic analysis can help with this. To facilitate genomic analysis, an online bioinformatics tool called P. insidiosum (Pins) Gene Table v2.0 was developed. This tool includes genomic data from 37 genetically diverse P. insidiosum strains and four related species. The database contains 732,686 genes, grouped into 80,061 unique clusters and further divided into core and variable categories at genus, species, and genotype levels. A high-resolution phylogenomic relationship among P. insidiosum strains and other oomycetes was projected through hierarchical clustering and core gene analyses. 3156 P. insidiosum-specific genes were shared among all genotypes and may be responsible for causing disease in humans and animals. After comparing these species-specific genes to the MvirDB database, 112 had significant matches with 66 known virulence proteins, some of which might be involved in vascular occlusion, which is a pathological feature of pythiosis. The correlation of genotypes, geographic origins, and affected hosts of P. insidiosum suggests that clade-I strains are more specific to animals, while clade-II/III strains are more specific to humans. The clade-specific genes might link to host preference. In summary, Pins Gene Table v2.0 is a comprehensive genome database accessible to users with minimal bioinformatics experience for the analysis of P. insidiosum genomes.

Functional host-specific adaptation of the intestinal microbiome in hominids

Fine-scale knowledge of the changes in composition and function of the human gut microbiome compared that of our closest relatives is critical for understanding the evolutionary processes underlying its developmental trajectory. To infer taxonomic and functional changes in the gut microbiome across hominids at different timescales, we perform high-resolution metagenomic-based analyzes of the fecal microbiome from over two hundred samples including diverse human populations, as well as wild-living chimpanzees, bonobos, and gorillas. We find human-associated taxa depleted within non-human apes and patterns of host-specific gut microbiota, suggesting the widespread acquisition of novel microbial clades along the evolutionary divergence of hosts. In contrast, we reveal multiple lines of evidence for a pervasive loss of diversity in human populations in correlation with a high Human Development Index, including evolutionarily conserved clades. Similarly, patterns of co-phylogeny between microbes and hosts are found to be disrupted in humans. Together with identifying individual microbial taxa and functional adaptations that correlate to host phylogeny, these findings offer insights into specific candidates playing a role in the diverging trajectories of the gut microbiome of hominids. We find that repeated horizontal gene transfer and gene loss, as well as the adaptation to transient microaerobic conditions appear to have played a role in the evolution of the human gut microbiome.

Transcriptional landscape of the pMP7017 megaplasmid and its impact on the Bifidobacterium breve UCC2003 transcriptome.

The 190 kb megaplasmid pMP7017 of Bifidobacterium breve JCM7017 represents the first conjugative and largest plasmid characterised within this genus to date. In the current study, we adopted an integrated approach combining transcriptomics, whole genome comparative analysis and metagenomic data mining to understand the biology of pMP7017 and related megaplasmids, and to assess the impact of plasmid-carriage on the host strain. The data generated revealed variations within basic features of promoter elements which correlate with a high level of transcription on the plasmid and highlight the transcriptional activity of genes encoding both offensive and defensive adaptations, including a Type IIL restriction-modification system, an anti-restriction system and four Type II toxin-antitoxin systems. Furthermore, a highly transcribed tmRNA, which likely provides translational support to the host strain, was identified, making pMP7017 the first plasmid of the Bifidobacterium genus and the smallest plasmid known to express a tmRNA. Analyses of synteny and variability among pMP7017 and related plasmids indicate substantial diversity in gene organisation and accessory gene cargo highlighting diverse (co-)evolution and potential host-specific rearrangements and adaptations. Systematic analysis of the codon usage profile of transcriptionally active pMP7017-encoded genes suggests that pMP7017 originated from (sub)species of Bifidobacterium longum. Furthermore, mining of metagenomic data suggests the presence of pMP7017-homologues in ~10% of microbiome samples, mostly infants and/or mothers from various geographical locations. Comparative transcriptome analysis of the B. breve UCC2003 chromosome in the presence or absence of pMP7017 revealed differential expression of genes representing 8% of the total gene pool. Genes involved in genetic information processing were exclusively upregulated, while altered expression of genes involved in biofilm production and polysaccharide biosynthesis was also observed.

Ecomorphological convergence in the walking leg dactyli of two clades of ascidian- and mollusc-associated shrimps (Decapoda: Caridea: Palaemonidae).

Symbiotic species, living within or on the surface of host organisms, may evolve a wide range of adaptations as a result of various selection pressures, host specificity of the symbiont and the nature of the symbiosis. In tropical marine coral reef ecosystems, palaemonid shrimps (Crustacea: Decapoda: Caridea) live in association with at least five different invertebrate phyla. Host switches between (distantly) related host groups, and the thereby associated selection pressures were found to play a major role in the diversification of these shrimp lineages, giving rise to various host-specific adaptations. Two lineages of palaemonid shrimp, which have switched from an ectosymbiotic association towards endosymbiosis, are studied for their morphological diversification and possible convergence. Special attention is given to the between-phyla host switches involving ascidian and bivalve hosts, which are characteristic for these lineages. Using landmark-based (phylo)morphospace analyses and Scanning Electron Microscopy, the walking leg dactylus shape and the microstructures on these dactyli are studied. No specific bivalve- or ascidian-associated morphotypes were found, but morphological convergence in dactylus morphology was found in various species within the two studied clades with similar host groups. In addition, multiple lineages of bivalve-associated species appear to be morphologically diverging more than their ascidian-associated relatives, with 'intermediate' morphotypes found near host-switching events.

The gut microbiome mediates adaptation to scarce food in Coleoptera.

Beetles are ubiquitous cave invertebrates worldwide that adapted to scarce subterranean resources when they colonized caves. Here, we investigated the potential role of gut microbiota in the adaptation of beetles to caves from different climatic regions of the Carpathians. The beetles' microbiota was host-specific, reflecting phylogenetic and nutritional adaptation. The microbial community structure further resolved conspecific beetles by caves suggesting microbiota-host coevolution and influences by local environmental factors. The detritivore species hosted a variety of bacteria known to decompose and ferment organic matter, suggesting turnover and host cooperative digestion of the sedimentary microbiota and allochthonous-derived nutrients. The cave Carabidae, with strong mandibula, adapted to predation and scavenging of animal and plant remains, had distinct microbiota dominated by symbiotic lineages Spiroplasma or Wolbachia. All beetles had relatively high levels of fermentative Carnobacterium and Vagococcus involved in lipid accumulation and a reduction of metabolic activity, and both features characterize adaptation to caves.

Plant Responses of Maize to Two formae speciales of Sporisorium reilianum Support Recent Fungal Host Jump.

Host jumps are a major factor for the emergence of new fungal pathogens. In the evolution of smut fungi, a putative host jump occurred in Sporisorium reilianum that today exists in two host-adapted formae speciales, the sorghum-pathogenic S. reilianum f. sp. reilianum and maize-pathogenic S. reilianum f. sp. zeae. To understand the molecular host-specific adaptation to maize, we compared the transcriptomes of maize leaves colonized by both formae speciales. We found that both varieties induce many common defense response-associated genes, indicating that both are recognized by the plant as pathogens. S. reilianum f. sp. reilianum additionally induced genes involved in systemic acquired resistance. In contrast, only S. reilianum f. sp. zeae induced expression of chorismate mutases that function in reducing the level of precursors for generation of the defense compound salicylic acid (SA), as well as oxylipin biosynthesis enzymes necessary for generation of the SA antagonist jasmonic acid (JA). In accordance, we found reduced SA levels as well as elevated JA and JA-Ile levels in maize leaves inoculated with the maize-adapted variety. These findings support a model of the emergence of the maize-pathogenic variety from a sorghum-specific ancestor following a recent host jump.

Genomic Heterogeneity of Cryptosporidium parvum Isolates From Children in Bangladesh: Implications for Parasite Biology and Human Infection.

Cryptosporidium species are a major cause of diarrhea and associated with growth failure. There is currently only limited knowledge of the parasite's genomic variability. We report a genomic analysis of Cryptosporidium parvum isolated from Bangladeshi infants and reanalysis of sequences from the United Kingdom. Human isolates from both locations shared 154 variants not present in the cattle-derived reference genome, suggesting host-specific adaptation of the parasite. Remarkably 34.6% of single-nucleotide polymorphisms unique to human isolates were nonsynonymous and 8.2% of these were in secreted proteins. Linkage disequilibrium decay indicated frequent recombination. The genetic diversity of C. parvum has potential implications for vaccine and therapeutic design. Clinical Trials Registration. NCT02764918.

Host-specific adaptation drove the coevolution of leek yellow stripe virus and Allium plants.

Host adaptation plays a crucial role in virus evolution and is a consequence of long-term interactions between virus and host in a complex arms race between host RNA silencing and viral RNA silencing suppressor (RSS) as counterdefense. Leek yellow stripe virus (LYSV), a potyvirus causing yield loss of garlic, infects several species of Allium plants. The unexpected discovery of an interspecific hybrid of garlic, leek, and great-headed (GH) garlic motivated us to explore the host-adaptive evolution of LYSV. Here, using Bayesian phylogenetic comparative methods and a functional assay of viral RSS activity, we show that the evolutionary context of LYSV has been shaped by the host adaptation of the virus during its coevolution with Allium plants. Our phylogenetic analysis revealed that LYSV isolates from leek and their taxonomic relatives (Allium ampeloprasum complex; AAC) formed a distinct monophyletic clade separate from garlic isolates and are likely to be uniquely adapted to AAC. Our comparative studies on viral accumulation indicated that LYSV accumulated at a low level in leek, whereas LYSVs were abundant in other Allium species such as garlic and its relatives. When RSS activity of the viral P1 and HC-Pro of leek LYSV isolate was analyzed, significant synergism in RSS activity between the two proteins was observed in leek but not in other species, suggesting that viral RSS activity may be important for the viral host-specific adaptation. We thus consider that LYSV may have undergone host-specific evolution at least in leek, which must be driven by speciation of its Allium hosts. IMPORTANCE Potyviruses are the most abundant plant RNA viruses and are extremely diversified in terms of their wide host range. Due to frequent host switching during their evolution, host-specific adaptation of potyviruses may have been shaped by numerous host factors. However, any critical determinants for viral host range remain largely unknown, possibly because of the repeated gain and loss of virus infectivity of plants. Leek yellow stripe virus (LYSV) is a species of the genus Potyvirus, which has a relatively narrow host range, generally limited to hosts in the genus Allium. Our investigations on leek and leek relatives (Allium ampeloprasum complex), which must have been generated through interspecies hybridization, revealed that LYSV accumulation remained low in leek as a result of viral host adaptation in competition with host resistance such as RNA silencing. This study presents LYSV as an ideal model to study the process of host-adaptive evolution and virus-host coevolution.

The evolving role of methanogenic archaea in mammalian microbiomes.

Methanogenic archaea (methanogens) represent a diverse group of microorganisms that inhabit various environmental and host-associated microbiomes. These organisms play an essential role in global carbon cycling given their ability to produce methane, a potent greenhouse gas, as a by-product of their energy production. Recent advances in culture-independent and -dependent studies have highlighted an increased prevalence of methanogens in the host-associated microbiome of diverse animal species. Moreover, there is increasing evidence that methanogens, and/or the methane they produce, may play a substantial role in human health and disease. This review addresses the expanding host-range and the emerging view of host-specific adaptations in methanogen biology and ecology, and the implications for host health and disease.

Feeding experience on hairy vetch induces a loss of host-specific adaptation to maize and changes in the gut microbiota communities of the fall armyworm

Feeding experience on hairy vetch induces a loss of host-specific adaptation to maize and changes in the gut microbiota communities of the fall armyworm

Host Derivation of Sindbis Virus Influences Mammalian Type I Interferon Response to Infection.

Arboviruses are defined by their ability to replicate in both mosquito vectors and mammalian hosts. There is good evidence that arboviruses "prime" their progeny for infection of the next host, such as via differential glycosylation of their outer glycoproteins or packaging of host ribosomal subunits. We and others have previously shown that mosquito-derived viruses more efficiently infect mammalian cells than mammalian-derived viruses. These observations are consistent with arboviruses acquiring host-specific adaptations, and we hypothesized that a virus derived from either the mammalian host or mosquito vector elicits different responses when infecting the mammalian host. Here, we perform an RNA-sequencing analysis of the transcriptional response of Human Embryonic Kidney 293 (HEK-293) cells to infection with either mosquito (Aedes albopictus, C7/10)- or mammalian (Baby Hamster Kidney, BHK-21)-derived Sindbis virus (SINV). We show that the C7/10-derived virus infection leads to a more robust transcriptional response in HEK-293s compared to infection with the BHK-derived virus. Surprisingly, despite more efficient infection, we found an increase in interferon-β (IFN-β) and interferon-stimulated gene (ISG) transcripts in response to the C7/10-derived virus infection versus the BHK-derived virus infection. However, translation of interferon-stimulated genes was lower in HEK-293s infected with the C7/10-derived virus, starkly contrasting with the transcriptional response. This inhibition of ISG translation is reflective of a more rapid overall shut-off of host cell translation following infection with the C7/10-derived virus. Finally, we show that the C7/10-derived virus infection of HEK-293 cells leads to elevated levels of phosphorylated eukaryotic translation elongation factor-2 (eEF2), identifying a potential mechanism leading to the more rapid shut-off of host translation. We postulate that the rapid shut-off of host translation in mammalian cells infected with the mosquito-derived virus acts to counter the IFN-β-stimulated transcriptional response.

Maintenance of a host-specific minority mutation in the West Nile virus NS3

West Nile virus (WNV), the most prevalent arthropod-borne virus (arbovirus) in the United States, is maintained in a cycle between Culex spp. mosquitoes and birds. Arboviruses exist within hosts and vectors as a diverse set of closely related genotypes. In theory, this genetic diversity can facilitate adaptation to distinct environments during host cycling, yet host-specific fitness of minority genotypes has not been assessed. Utilizing WNV deep-sequencing data, we previously identified a naturally occurring, mosquito-biased substitution, NS3 P319L. Using both cell culture and experimental infection in natural hosts, we demonstrated that this substitution confers attenuation in vertebrate hosts and increased transmissibility by mosquitoes. Biochemical assays demonstrated temperature-sensitive ATPase activity consistent with host-specific phenotypes. Together these data confirm the maintenance of host-specific minority variants in arbovirus mutant swarms, suggest a unique role for NS3 in viral fitness, and demonstrate that intrahost sequence data can inform mechanisms of host-specific adaptation.

Multiple genotypes infection and molecular characterization of Torque teno neovison virus: A novel Anelloviridae of mink in China

A novel Torque teno neovison virus (TTVs) was identified in specimens collected from dead mink during an outbreak of the Aleutian mink disease virus. Eighteen complete genomic sequences were obtained, ranging from 2109 to 2158 nucleotides in length and consisting of an untranslated region and three open reading frames. The genomic organization of mink TTVs is similar to previously reported anelloviruses. However, the deduced amino acid sequence of its ORF1 protein shows genetic diversity compared to related anelloviruses, suggesting that it represents a putative new species within the Anelloviridae family. This study provides a detailed molecular characterization of the novel mink anelloviruses, including its codon usage pattern, origin, and evolution. Analysis of the viral genomic sequences reveals the existence of multiple genotypes of co-infection. Principal component analysis and phylogenetic trees confirm the coexistence of multiple genotypes. Furthermore, the codon usage analyses indicate that mink TTVs have a genotype-specific codon usage pattern and show a low codon usage bias. Host-specific adaptation analysis suggests that TTVs are less adapted to mink. The possible origin and evolutionary history of mink TTVs were elucidated. Mink TTVs was genetically closely related to giant panda anellovirus, representing a new species. The observed incongruence between the phylogenetic history of TTVs and that of their hosts suggests that the evolution of anellovirus is largely determined by cross-species transmission. The study provides insights into the co-infection and genetic evolution of anellovirus in China.

Human-to-Animal Transmission of SARS-CoV-2, South Korea, 2021

To investigate SARS-CoV-2 transmission from humans to animals in Seoul, South Korea, we submitted samples from companion animals owned by persons with confirmed COVID-19. Real-time PCR indicated higher SARS-CoV-2 viral infection rates for dogs and cats than previously reported from the United States and Europe. Host-specific adaptations could introduce mutant SARS-CoV-2 to humans.

Phylogeny of Shrew- and Mole-Borne Hantaviruses in Poland and Ukraine.

Earlier, we demonstrated the co-circulation of genetically distinct non-rodent-borne hantaviruses, including Boginia virus (BOGV) in the Eurasian water shrew (Neomys fodiens), Seewis virus (SWSV) in the Eurasian common shrew (Sorex araneus) and Nova virus (NVAV) in the European mole (Talpa europaea), in central Poland. To further investigate the phylogeny of hantaviruses harbored by soricid and talpid reservoir hosts, we analyzed RNAlater®-preserved lung tissues from 320 shrews and 26 moles, both captured during 1990-2017 across Poland, and 10 European moles from Ukraine for hantavirus RNA through RT-PCR and DNA sequencing. SWSV and Altai virus (ALTV) were detected in Sorex araneus and Sorex minutus in Boginia and the Białowieża Forest, respectively, and NVAV was detected in Talpa europaea in Huta Dłutowska, Poland, and in Lviv, Ukraine. Phylogenetic analyses using maximum-likelihood and Bayesian methods showed geography-specific lineages of SWSV in Poland and elsewhere in Eurasia and of NVAV in Poland and Ukraine. The ATLV strain in Sorex minutus from the Białowieża Forest on the Polish-Belarusian border was distantly related to the ATLV strain previously reported in Sorex minutus from Chmiel in southeastern Poland. Overall, the gene phylogenies found support long-standing host-specific adaptation.