# Trematode genetic patterns at host individual and population scales provide insights about infection mechanisms

Simão Correia, Sergio Fernández-Boo, Luísa Magalhães, Xavier de Montaudouin, Guillemine Daffe, Robert Poulin, Manuel Vera

► To cite this version:

Simão Correia, Sergio Fernández-Boo, Luísa Magalhães, Xavier de Montaudouin, Guillemine Daffe, et al. Trematode genetic patterns at host individual and population scales provide insights about infection mechanisms. Parasitology, 2023, 150 (13), pp.1207-1220. 10.1017/s0031182023000987. hal-04513408

HAL Id: hal-04513408

https://hal.science/hal-04513408v1

Submitted on 20 Mar 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Trematode genetic patterns at host individual and population scales provide insights about infection mechanisms

## Research Article

**Cite this article:** Correia S, Fernández-Boo S, Magalhães L, de Montaudouin X, Daffe G, Poulin R, Vera M (2023). Trematode genetic patterns at host individual and population scales provide insights about infection mechanisms. *Parasitology* **150**, 1207–1220.
https://doi.org/10.1017/S0031182023000987

Received: 22 May 2023
Revised: 2 October 2023
Accepted: 7 October 2023
First published online: 20 October 2023

**Keywords:**
*Bucephalus minimus*; *Cerastoderma edule*; clonal diversity; COI; host–parasite interactions; parasite; population genetics

**Corresponding authors:**
Simão Correia;
Email: simaocorreia@ua.pt;
simaopedro.correia@usc.es;
Manuel Vera;
Email: manuel.vera@usc.es

Simão Correia<fnref n="1" />,²,³,⁴ ID, Sergio Fernández-Boo² ID, Luísa Magalhães<fnref n="1" /> ID, Xavier de Montaudouin⁵ ID, Guillemine Daffe⁶, Robert Poulin⁴ ID and Manuel Vera³ ID

## Abstract

Multiple parasites can infect a single host, creating a dynamic environment where each parasite must compete over host resources. Such interactions can cause greater harm to the host than single infections and can also have negative consequences for the parasites themselves. In their first intermediate hosts, trematodes multiply asexually and can eventually reach up to 20% of the host's biomass. In most species, it is unclear whether this biomass results from a single infection or co-infection by 2 or more infective stages (miracidia), the latter being more likely *a priori* in areas where prevalence of infection is high. Using as model system the trematode *Bucephalus minimus* and its first intermediate host cockles, we examined the genetic diversity of the cytochrome c oxidase subunit I region in *B. minimus* from 3 distinct geographical areas and performed a phylogeographic study of *B. minimus* populations along the Northeast Atlantic coast. Within localities, the high genetic variability found across trematodes infecting different individual cockles, compared to the absence of variability within the same host, suggests that infections could be generally originating from a single miracidium. On a large spatial scale, we uncovered significant population structure of *B. minimus*, specifically between the north and south of Bay of Biscay. Although other explanations are possible, we suggest this pattern may be driven by the population structure of the final host.

## Introduction

With about 45 000 species described in a wide range of ecosystems, trematodes are one of the most common and widespread group of parasites (Carlson *et al.*, 2020). They can be found at almost across all trophic levels of dynamic food chains (Bartoli and Gibson, 2007). Trematodes are an important component of ecosystem biodiversity with significant impacts at the host individuals (Schulte-Oehlmann *et al.*, 1997; Curtis *et al.*, 2000; Thieltges, 2004), host populations (Fredensborg *et al.*, 2005) and ecosystem communities (Poulin, 1999; Mouritsen and Poulin, 2002; Goedknegt *et al.*, 2016). Besides, by contributing to the nutrient cycle, acting as indicators of environmental changes or as proxy of environmental diversity (due to multi-host life cycles), trematode presence may indicate a healthy and resilient ecosystem (Johnson *et al.*, 2010; Hatcher and Dunn, 2011).

Trematodes have a complex life cycle that alternates between free-living and parasitic stages. The miracidium, the trematode larva hatched from an egg, infects the first intermediate host (usually a mollusc) and transforms into sporocysts or rediae (parasitic stage). At this stage, sporocysts or rediae, through asexual multiplication, produce cercariae (free-living stage) that emerge from the first host to infect the second intermediate host (a vertebrate or invertebrate) where they settle as metacercariae. In the trematode's final host (a vertebrate), after ingestion of the second host, metacercariae develop into adult flukes, reproduce sexually, and complete the life cycle (Cribb *et al.*, 2003; Bartoli and Gibson, 2007). In the first intermediate host, sporocyst stages are overtly destructive, replacing host tissue and reaching up to 20% of the host's biomass (Dubois *et al.*, 2009; Preston *et al.*, 2013), with direct consequences for host reproduction (Carballal *et al.*, 2001), growth (Bowers, 1969) and energy demand (Jokela *et al.*, 1993), leading to eventual host death (Thieltges, 2006). In most species, it is currently unknown whether this sporocyst biomass results from a single miracidium, that excludes other miracidia by predation or intraspecific competition, or from co-infection. If co-infection is the rule, trematode invasion may result in a burden that the host might not be able to bear (Fredensborg and Poulin, 2005; Mideo, 2009). On the other hand, co-infection can, occasionally, benefit the host by lessening the overall burden of infection, by reducing

© The Author(s), 2023. Published by Cambridge University Press. This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.

<figure><img src="image_2.png" /></figure>

CAMBRIDGE
UNIVERSITY PRESS

<footnote n="1">¹Department of Biology, CESAM, University of Aveiro, 3810-193 Aveiro, Portugal; ²Aquatic and Animal Health Group, CIIMAR, University of Porto, 4450-208 Matosinhos, Portugal; ³Department of Zoology, Genetics and Physical Anthropology, Campus Terra, University of Santiago de Compostela, 27002 Lugo, Spain; ⁴Department of Zoology, University of Otago, 9054 Dunedin, New Zealand; ⁵CNRS, UMR EPOC, Station Marine, University of Bordeaux, F-33120 Arcachon, France and ⁶Université de Bordeaux, CNRS, Observatoire Aquitain des Sciences de l'Univers, F-33615 Pessac, France</footnote>
parasite infection success, or by strengthening the host's immune response, leading to higher resistance to infection (Dumont *et al.*, 2007; Balmer *et al.*, 2009).

Determining the genetic diversity of trematode sporocysts within and among individual first intermediate hosts is therefore important to understand the biology, behaviour and evolutionary patterns of these parasites. Nonetheless, current knowledge regarding biology of these parasites, and particularly about the sporocyst life stage, is still very scarce, despite trematodes' wide distribution and importance to the ecosystem. Few studies have focused on the conspecific diversity of trematode sporocysts within the first intermediate host (Rauch *et al.*, 2005; Keeney *et al.*, 2007; Lagrue *et al.*, 2007), showing that the likelihood of infection by conspecifics increased with the prevalence of the parasite in the community (Keeney *et al.*, 2008; Louhi *et al.*, 2013). Moreover, at the population level, the host has a significant impact on genetic diversity and population structure of parasites, with substantial gene flow occurring in parasite species with efficient dispersion mechanisms (Agola *et al.*, 2009; Feis *et al.*, 2015). However, most studies of trematode genetic diversity have focused on intermediate and final hosts, with particular emphasis on trematodes with harmful impacts on human (Theron *et al.*, 2004; Bell *et al.*, 2006; Balmer *et al.*, 2009) or socio-economically important species, namely fish (Vilas *et al.*, 2003; Criscione and Blouin, 2006).

Bucephalus minimus is a marine trematode parasite that occurs in several aquatic systems along the Northeast Atlantic coast and Mediterranean Sea (Magalhães *et al.*, 2015). In the Atlantic area, this parasite infects the European edible cockle, *Cerastoderma edule*, which serves as the first intermediate host when a miracidium penetrates its tissue. In cockles, the prevalence of this parasite varies greatly among coastal systems and season; depending on the time since infection, the parasite's dry mass in infected cockles can range from 1 to 20% of the total living tissue within the cockle shell (Magalhães *et al.*, 2015; de Montaudouin *et al.*, 2021). *Bucephalus minimus* initially infects the cockle's gonad and digestive gland but promptly spreads to other parts of the host, eventually invading the entire body (Desclaux *et al.*, 2002; de Montaudouin *et al.*, 2009). Infection by *B. minimus* results in castration (Carballal *et al.*, 2001) and energy consumption (Dubois *et al.*, 2009), leading to starvation and autolysis of the host's digestive tract. This parasite is considered as one of the most harmful trematode parasites infecting *C. edule* (Magalhães *et al.*, 2015; de Montaudouin *et al.*, 2021). Inside cockles, *B. minimus* produces sporocysts and cercariae through asexual multiplication, which emerge and infect the goby *Pomatoschistus* spp. (second intermediate host), where they encyst and develop into metacercariae. The final host, *Dicentrarchus labrax*, the European seabass, is infected after consumption of parasitized gobies. Metacercariae develop into adult flukes and produce eggs, through sexual reproduction, to complete the cycle (Pina *et al.*, 2009; Magalhães *et al.*, 2015).

Due to the limited information regarding this trematode's sporocyst stages, the primary goal of the present study was to assess the genetic variability of the cytochrome c oxidase subunit I (COI) region of *B. minimus* sporocyst DNA within and among the first intermediate host, *C. edule*. We tested the hypothesis that higher genetic variability at the host individual scale (resulting in 2 or more haplotypes among sporocysts in a single cockle) are more common in localities with high prevalence of infection, where joint infections should be more frequent by chance alone. As our second goal, a phylogeographic study of *B. minimus* was also carried out combining information available in the literature and from samples taken on cockle beds that were examined in this study for the first time (i.e. Iberian Peninsula and Great Britain).

## Materials and methods

### Bucephalus minimus samples

The genetic variability of *B. minimus* sporocyst haplotypes within the same host was studied by collecting specimens present in the first intermediate host, the edible cockle, at 3 different beds with different prevalence along the European Atlantic coast: Ria de Aveiro, Aveiro, Portugal (lowest prevalence [Magalhães *et al.*, 2018]); de la Ramallosa Lagoon, Baiona, Spain (moderate prevalence [Intecmar, 2021]); and Île aux Oiseaux, Arcachon, France (with high levels of prevalence [Magalhães *et al.*, 2015]) (Fig. 1). Adult cockles (between 20 and 30 mm shell length) were haphazardly collected at low tide and dissected in the laboratory to morphologically identify *B. minimus* infection. The flesh was then transferred and observed under a stereomicroscope by carefully compressing it between 2 sterilized glass slides. Four sporocyst replicates were extracted per cockle (in a total of 5 cockles per sampling site) using forceps, and preserved separately in 100% ethanol at -20°C. All material was sterilized between samples. To enhance the possibility of different clones, sporocysts were taken from different infected tissues of the same cockle (i.e. 1 from the foot, 1 from the gills and 2 from the digestive gland).

### DNA isolation, amplification and sequencing

Genomic DNA extraction from *B. minimus* specimens was performed using E.Z.N.A Mollusc DNA kit (Omega Bio-Tek, Norcross, GA, USA) in accordance with the manufacturer's instructions. Nanodrop was used to assess DNA concentration, and, if needed, aliquots were created to dilute DNA to approximately 30 ng µL⁻¹.

The mitochondrial COI fragment was amplified using the MplatCOX1-dF (5'-TTW CIT TRG ATC ATA AG-3') and MplatCOX1-dR (5'-TGA AAY AAY AII GGA TCI CCA CC-3') primers (Moszczynska *et al.*, 2009), resulting in sequences of 587 bp. The polymerase chain reaction (PCR) was carried out in a final volume of 20 µL composed of 1X reaction buffer, 2.5 mM MgCl₂, 100 µM deoxynucleotide triphosphates (dNTPs), 0.5 µM of forward and reverse primers, 0.65 units of ThermoFisher AmpliTaq Gold DNA polymerase (ThermoFisher, Waltham, MA, USA) and 60 ng of DNA. The PCR programme employed had an initial denaturation step for 10 min at 95°C, followed by 35 cycles of 30 s at 94°C, 30 s at 50°C and 1 min at 72°C and a final extension of 10 min at 72°C. Following this initial PCR, a second PCR was carried out using the same conditions previously described but using 2 µL of PCR product instead of DNA. The amplified PCR products were analysed by electrophoresis through a 1% agarose gel dyed with SYBR safe DNA Gel Stain (Invitrogen, Carlsbad, CA, USA) and visualized under UV light.

PCR products were enzymatically purified with ExoSAP mix (10 µL of PCR product, 0.6 units of EXO I [DNA nuclease] and 0.3 units of shrimp alkaline phosphatase [SAP] for a final volume of 12 µL) under the following conditions: 60 min at 37°C and 15 min at 85°C. Purified PCR products were sequenced using the ABI Prism BigDye™ Terminator v3.1 Cycle Sequencing Kit protocol on an ABI Prism 3730 xl automatic sequencer (Applied Biosystems, Foster City, CA, USA). All sequences obtained in this study were deposited in GenBank (accession numbers: OQ625925-OQ625936; Table 1). Variable sites were manually checked using the SEQSCAPE 2.5 program (Applied Biosystems) and aligned using ClustalW algorithm implemented on BioEdit v.7.2.5 (Hall, 1999) with LaB haplotype (GenBank accession number: KF880429.1) as reference. Identification of the different haplotypes in the *B. minimus* specimens analysed was carried out using the software DNAsp v5.10 (Librado and Rozas, 2009). Finally, to check for the presence of premature
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Geographical location of the *Cerastoderma edule* cockle beds sampled for the study of *Bucephalus minimus* genetic variability at host level (Aveiro, Baiona and Arcachon [in italics]) and for phylogeographic analysis. \*Cockle beds sampled for the first time in this study.</figcaption></figure>

stop codons, COI haplotypes identified were translated to amino acid sequences using the flatworm mtDNA code in the online software EMBOSS Transeq (Rice *et al.*, 2000; Goujon *et al.*, 2010).

## Data analysis

### Bucephalus minimus genetic variability at host level

To determine the genetic variability of *B. minimus* specimens within the same cockle and among different cockles, haplotypes found in the different tissues of each analysed cockle were compared using BioEdit. A haplotype network was built using the obtained data, identifying haplotypes for each cockle with a different colour. The haplotype network was constructed by calculating the distance (based on number of base pair differences) between DNA sequences and determining the number of mutations between haplotypes using ‘pegas’ and ‘ape’ packages of R Statistical Software v.4.2.2 (Paradis, 2010; Paradis and Schliep, 2019).

### Phylogeographic analysis

For phylogeographic analyses, together with the specimens collected in the present study, DNA extractions of 13 specimens from 7 cockle beds sampled as part of the COCKLES Interreg project (http://cockles-project.eu/) were sequenced as previously described. Samples were available for Aveiro, Portugal (1 sample), Noia, Spain (1 sample), Arcachon, France (1 sample), Bay of Somme, France (4 samples), The Dee, Wales (2 samples), Burry Inlet, Wales (2 samples) and Wadden Sea, the Netherlands (2 samples). Additionally, the analysis included another 54 COI

sequences available in the GenBank database retrieved in January 2023 (see Feis *et al.*, 2015 and Table 1). In total, sequences represented specimens from 11 different cockle beds located in 7 countries, covering a large part of the natural distributional range of cockles in the Atlantic area (Fig. 1).

Phylogenetic relationships were studied using different and complementary approaches. First, a haplotype network was computed for the full dataset as previously described, which included information regarding the cockle bed in which the haplotypes were found and the frequency of occurrence. Identification of unique haplotypes present in the dataset was carried out with DnaSP v.5.10 (Librado and Rozas, 2009). Using MEGA X software (Kumar *et al.*, 2018), the Hasegawa-Kishino-Yano nucleotide substitution rate (HKY) with a $\gamma$ value of 0.655 and invariable sites of 0.728 was identified as the most probable nucleotide substitution model for our data. Phylogenetic trees were constructed using exclusively the different haplotypes identified and the nucleotide substitution model described above. Maximum likelihood (ML) and neighbour-joining rooted and unrooted trees were constructed using the R Statistical Software v.4.2.2. The rooted tree was created using the COI sequences of *Rhipidocotyle* sp. (a trematode from the same family as *B. minimus*, Bucephalidae, GenBank accession number: KM538111.1) and *Himasthla quissetensis* (a trematode that infects cockles as second intermediate host but from a different family, Himasthliidae; GenBank accession number: MN272732.1) as outgroups. The sequences were trimmed to 540 bp to remove missing data. The unrooted trees were constructed using the full 587 bp sequences of *B. minimus*. The robustness of the branches for
<table><caption>**Table 1.** Accession number for *B. minimus* COI gene DNA sequences downloaded (with reference) and deposited (in bold) in GenBank</caption><thead><tr><th>Haplotype name</th><th>GenBank accession number</th><th>Reference</th></tr></thead><tbody><tr><td>LaA</td><td>KF880428.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaB</td><td>KF880429.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaC</td><td>KF880430.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaD</td><td>KF880431.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaE</td><td>KF880432.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaF</td><td>KF880433.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaG</td><td>KF880434.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaH</td><td>KF880435.1</td><td>Feis et al. (2015)</td></tr><tr><td>Lal</td><td>KF880436.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaJ</td><td>KF880437.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaK</td><td>KF880438.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaL</td><td>KF880439.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaM</td><td>KF880440.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaN</td><td>KF880441.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaO</td><td>KF880442.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaP</td><td>KF880443.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaQ</td><td>KF880444.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaR</td><td>KF880445.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaS</td><td>KF880446.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaT</td><td>KF880447.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaU</td><td>KF880448.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaV</td><td>KF880449.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaW</td><td>KF880450.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaX</td><td>KF880451.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaY</td><td>KF880452.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaZ</td><td>KF880453.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaAA</td><td>KF880454.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaAB</td><td>KF880455.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaAC</td><td>KF880456.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaAD</td><td>KF880457.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaAE</td><td>KF880458.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaAF</td><td>KF880459.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaAG</td><td>KF880460.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaAH</td><td>KF880461.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaAI</td><td>KF880462.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaAJ</td><td>KF880463.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaAK</td><td>KF880464.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaAL</td><td>KF880465.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaAM</td><td>KF880466.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaAN</td><td>KF880467.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaAO</td><td>KF880468.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaAP</td><td>KF880469.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaAQ</td><td>KF880470.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaAR</td><td>KF880471.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaAS</td><td>KF880472.1</td><td>Feis et al. (2015)</td></tr></tbody></table>

(Continued)

<table><caption>**Table 1.** (Continued.)</caption><thead><tr><th>Haplotype name</th><th>GenBank accession number</th><th>Reference</th></tr></thead><tbody><tr><td>LaAT</td><td>KF880473.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaAV</td><td>KF880474.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaAW</td><td>KF880475.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaAX</td><td>KF880476.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaAY</td><td>KF880477.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaAZ</td><td>KF880478.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaBA</td><td>KF880479.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaBB</td><td>KF880480.1</td><td>Feis et al. (2015)</td></tr><tr><td>LaBC</td><td>KF880481.1</td><td>Feis et al. (2015)</td></tr><tr><td>BmA</td><td>OQ625925</td><td>This study</td></tr><tr><td>BmB</td><td>OQ625926</td><td>This study</td></tr><tr><td>BmC</td><td>OQ625927</td><td>This study</td></tr><tr><td>BmD</td><td>OQ625928</td><td>This study</td></tr><tr><td>BmE</td><td>OQ625929</td><td>This study</td></tr><tr><td>BmF</td><td>OQ625930</td><td>This study</td></tr><tr><td>BmG</td><td>OQ625931</td><td>This study</td></tr><tr><td>BmH</td><td>OQ625932</td><td>This study</td></tr><tr><td>Bml</td><td>OQ625933</td><td>This study</td></tr><tr><td>BmJ</td><td>OQ625934</td><td>This study</td></tr><tr><td>BmK</td><td>OQ625935</td><td>This study</td></tr><tr><td>BmL</td><td>OQ625936</td><td>This study</td></tr></tbody></table>

the phylogenetic trees was estimated with 1000 bootstrap replicates and a likelihood ratio test was performed based on the minimum Akaike information criterion values for ML. All phylogeographic analyses were performed with the ‘ape’, ‘pegas’, ‘ggtree’ and ‘phangorn’ packages of R Statistical Software v.4.2.2 (Paradis, 2010; Schliep, 2011; Yu et al., 2017; Paradis and Schliep, 2019).

Genetic diversity parameters, calculated as haplotype diversity ($h$) and nucleotide diversity ($\pi$), were estimated within each cockle bed studied using Arlequin v.3.5.1.3 (Excoffier et al., 2005). The HKY model is not available in this software. For this reason, the Tamura–Nei (TN) model with a $\gamma$ value of 0.639 (similar to the HKY and identified as the third best option for MEGA X) was used for $\pi$ estimations. Genetic structure and population differentiation were assessed with global and pairwise coefficients of population differentiation applying the TN with $\gamma$ value of 0.639 substitution rate ($\phi_{ST}$ values) with Arlequin. Analysis of molecular variance (AMOVA) applying different models of a priori clustering (based on host population genetics [Souche et al., 2015; Vera et al., 2022] and observed data – see results) was carried out to study the distribution of genetic variation within ($\phi_{SC}$) and among ($\phi_{CT}$) bed groups using Arlequin. The significance for all the $\phi$ statistics was evaluated with 10 000 permutations.

## Results

### Bucephalus minimus genetic variability at host level

During this study, a total of 210 cockles were analysed, of which 17 were found to be infected with *B. minimus* (5 each in Aveiro and Baiona and 7 in Arcachon). The prevalence of *B. minimus* in the cockle beds sampled varied from 3.3% in Aveiro
(Portugal) to 23.3% in Arcachon (France). In Baiona (Spain), *B. minimus* was present in 16.7% of the sampled cockles. Five infected cockles per bed were used to extract 4 sporocysts per cockle, yielding a total of 60 samples. Fifty-six samples were successfully sequenced, while 4 samples, from a single cockle from Arcachon (France), were not successfully sequenced due to DNA extraction problems. From the 56 sequenced sporocysts, belonging to 14 infected cockles, 12 different haplotypes were identified, with 5, 3 and 4 haplotypes found in Aveiro, Baiona and Arcachon, respectively (Table 2). Six of the identified haplotypes were characterized for the first time (named as BmA-BmF; see Table 1).

All sporocysts of *B. minimus* from the same cockle had identical haplotype, however *B. minimus* haplotypes identified in different cockles from the same bed were different, except in Baiona where the haplotype found in 3 different cockles was identical (haplotype LaE; see Table 2 and Fig. 2). Moreover, haplotypes were not shared among the 3 beds (Fig. 2).

## Phylogeographic analysis

From the 69 DNA sequenced samples (56 for the *B. minimus* haplotype genetic variability at host level study, and 13 from different European sites selected from the COCKLES project), 12 resulted in novel haplotypes (Table 1). No premature STOP codons were identified in these sequences (data not shown). Thus, when GenBank resources were included, a total of 162 COI gene sequences of *B. minimus* specimens from 11 cockle beds were analysed. From these available sequences, 66 represented unique haplotype sequences. Shared haplotypes (i.e. those found in more than 1 cockle bed) accounted for 17% of the total. The LaAQ haplotype was the most prevalent and abundant haplotype, occurring 42 times across 6 different beds. On the other hand, 83% of the haplotypes were exclusively found in a single bed, with several reported only once (i.e. singletons). Arcachon presented the highest number of different *B. minimus* haplotypes detected (19, see Table 3).

Due to the high variability found in the COI region, the analysis of *B. minimus* haplotypes across the various cockle beds produced a complex network made up of several closely connected haplotypes and associated mutational steps, with no more than 5 mutations separating any 2 successive haplotypes identified. A common haplotype (LaAQ), observed in several beds situated north of Arcachon (44°N), was located in the centre of the network, from which numerous other haplotypes diverged in a star-like pattern. These haplotypes were exclusively found in a single bed or shared between relatively close beds (Fig. 3). Nonetheless, haplotype clusters (i.e. haplogroups) were identified in specific geographic areas, and beds from the South (Merja Zerga, Aveiro, Baiona and Noia) and North (Bay of Somme, English Channel, Celtic Sea, Burry Inlet, The Dee and Wadden Sea) did not share any haplotype, with the

exception of the LaE haplotype identified in Baiona, Arcachon and the English Channel. Phylogenetic relationships observed in the network were also confirmed with the phylogenetic trees (Fig. 4).

The Arcachon bed, located in the centre of cockle's distributional range between the northern and southern geographic areas, exhibited the highest number of detected haplotypes (19), sharing haplotypes with locations from both regions (Table 3). Excluding Noia, where only 1 individual was analysed, haplotype diversity ranged from 0.3846 in Merja Zerga to 1.0000 in beds where all individuals analysed had a distinct haplotype (Aveiro, Burry Inlet and The Dee). Nucleotide diversity ranged from 0.0007 in Merja Zerga to 0.0076 in Baiona (Table 3). The high diversity and the haplotype distribution among locations were also reflected in the $\phi_{ST}$ values. Global $\phi_{ST}$ for the whole region was 0.2922 ($P$ value < 0.001). Many pairwise $\phi_{ST}$ values resulted in significant differences, although many comparisons between close locations were non-significant, mainly among those involving northern beds (Supplementary Table S1). Moreover, all pairwise $\phi_{ST}$ values involving Merja Zerga were high and highly significant ($P$ value < 0.001), suggesting the singularity of this bed (global $\phi_{ST}$ in the whole region excluding Merja Zerga = 0.1706, $P$ value < 0.001). These results suggest the presence of one northern group (composed by Bay of Somme, English Channel, Celtic Sea, Burry Inlet, The Dee and Wadden Sea) more homogeneous genetically ($\phi_{ST}$ = 0.0295, $P$ value = 0.073) than the southern one (composed by Merja Zerga, Aveiro, Baiona and Noia; $\phi_{ST}$ = 0.6694, $P$ value < 0.001), with Arcachon representing a potential contact region between both geographic areas. The $\phi_{ST}$ value in the northern group increased up to 0.0336 ($P$ value = 0.018) when Arcachon was included, while this value decreased in the southern group when this location was included although it remained quite high ($\phi_{ST}$ = 0.4184, $P$ value < 0.001). Hence, these results suggest a closer relationship of Arcachon with the northern group. AMOVA analysis assigned 34.58% of the genetic differentiation to differences between northern and southern groups ($\phi_{CT}$ = 0.3458, $P$ value = 0.008), this percentage being three times higher than those assigned to differences among beds within groups ($\phi_{SC}$ = 0.2109, $P$ value < 0.001, percentage of genetic differentiation = 13.80%). The AMOVA model including Arcachon in the northern group yielded similar values ($\phi_{CT}$ = 0.3100, $P$ value = 0.006, percentage of genetic differentiation = 31.00%; $\phi_{SC}$ = 0.1567, $P$ value < 0.001, percentage of genetic differentiation = 10.81%). This model assigned a higher percentage of genetic differentiation among groups and a lower percentage to differences among beds within groups than the model including Arcachon in the southern group ($\phi_{CT}$ = 0.1327, $P$ value = 0.048, percentage of genetic differentiation = 13.27%; $\phi_{SC}$ = 0.2299, $P$ value < 0.001, percentage of genetic differentiation = 19.13%), suggesting a more coherent grouping of the beds in the former model.

<table><caption>**Table 2.** Coordinates of each cockle bed, the number of analysed cockles ($N_{cockles}$), number of parasites sequenced ($N_{B. minimus}$), prevalence of *B. minimus*, number of haplotypes ($k$), number of polymorphic sites (PS) and haplotype composition (between parentheses the number of individuals bearing the same haplotype when different from one)</caption><thead><tr><th rowspan="2">Cockle bed</th><th rowspan="2">Coordinates</th><th colspan="2">N</th><th rowspan="2">Prevalence</th><th rowspan="2">k</th><th rowspan="2">PS</th><th rowspan="2">Haplotype composition</th></tr><tr><th>Cockles</th><th><i>B. minimus</i></th></tr></thead><tbody><tr><td>Aveiro</td><td>40.710123, -8.704596</td><td>5</td><td>20</td><td>3.3%</td><td>5</td><td>8</td><td>BmA, BmB, BmC, BmD, BmE</td></tr><tr><td>Baiona</td><td>42.117020, -8.820283</td><td>5</td><td>20</td><td>16.7%</td><td>3</td><td>3</td><td>BmF, LaD, LaE (3)</td></tr><tr><td>Arcachon</td><td>44.690111, -1.182944</td><td>4</td><td>16</td><td>23.3%</td><td>4</td><td>5</td><td>BmG, LaAE, LaAQ, LaAW</td></tr><tr><td><strong>General</strong></td><td></td><td><strong>14</strong></td><td><strong>56</strong></td><td><strong>8.1%</strong></td><td><strong>12</strong></td><td><strong>14</strong></td><td></td></tr></tbody></table>
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Haplotype network of *Bucephalus minimus* samples from Aveiro (Portugal), Baiona (Spain) and Arcachon (France) based on genetic distance (number of base pair differences) of cytochrome c oxidase subunit 1 (COI) gene sequences. Different haplotypes with respective names are represented by circles, with circle size proportional to observed frequency. Inferred mutation steps are shown by black dots. Colours depict samples taken from the same cockle.</figcaption></figure>

Discussion

*Bucephalus minimus* genetic variability at host level

Co-infection by multiple parasites, from the same or different species, within the same host is a well-recognized phenomenon in the parasitological literature (Poulin, 2001; Read and Taylor, 2001). This pattern has been extensively studied for several parasite species, namely with an impact on human health (Theron et al., 2004; Bell et al., 2006). For example, in the case of malaria, more than 5 strains have been found to be infecting the same host (Bell et al., 2006). The same trend was observed for the trematode parasite *Schistosoma mansoni* within their second and final host (Theron et al., 2004). Similar to what is observed for metacercarial or adult stages of trematode parasites, it would

be anticipated that different clones would infect the same first intermediate host when thousands of eggs per infected definitive host are shed into the water column, i.e. thousands of miracidia hatching within metres of each other. This was observed for some trematode species (Rauch *et al.*, 2005; Keeney *et al.*, 2007; Lagrue *et al.*, 2007). In the present study, only 1 COI haplotype was found inside each infected cockle (regardless of the samples' origin), in contrast with what has been previously recorded. Nevertheless, it should be noted that in the present study only the COI region (maternally inherited) was sequenced, while for previous studies, microsatellite markers were used to identify individual variability. In fact, microsatellite markers are more accurate for population structure analysis and individual identification since they are highly variable polymorphic regions
<table><caption>**Table 3.** Haplotype frequency for *B. minimus* COI gene per cockle bed including haplotype diversity ($h \pm$ standard deviation) and nucleotide diversity ($\pi \pm$ standard deviation)</caption><thead><tr><th>Haplotype</th><th>Merja Zerga</th><th>Aveiro</th><th>Baiona</th><th>Noia</th><th>Arcachon</th><th>Bay of Somme</th><th>English Channel</th><th>Celtic Sea</th><th>Burry Inlet</th><th>The Dee</th><th>Wadden Sea</th><th>Sum</th></tr></thead><tbody><tr><td>LaA</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>1</td></tr><tr><td>LaB</td><td>20</td><td>0</td><td>0</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>22</td></tr><tr><td>LaC</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>LaD</td><td>0</td><td>0</td><td>1</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>LaE</td><td>0</td><td>0</td><td>3</td><td>0</td><td>2</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>6</td></tr><tr><td>LaF</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>LaG</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>LaH</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>1</td></tr><tr><td>Lal</td><td>5</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>5</td></tr><tr><td>LaJ</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>8</td><td>1</td><td>0</td><td>0</td><td>9</td></tr><tr><td>LaK</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>LaL</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>LaM</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>LaN</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>LaO</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>LaP</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>1</td></tr><tr><td>LaQ</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>1</td></tr><tr><td>LaR</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>1</td><td>0</td><td>1</td><td>1</td><td>4</td></tr><tr><td>LaS</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>LaT</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>LaU</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>1</td></tr><tr><td>LaV</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>LaW</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>LaX</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>LaY</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>LaZ</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>LaAA</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>LaAB</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>LaAC</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>LaAD</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>1</td></tr><tr><td>LaAE</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td></tr></tbody></table>

(Continued)
<table><caption>**Table 3. (Continued.)**</caption><thead><tr><th>Haplotype</th><th>Merja Zerga</th><th>Aveiro</th><th>Baiona</th><th>Noia</th><th>Arcachon</th><th>Bay of Somme</th><th>English Channel</th><th>Celtic Sea</th><th>Burry Inlet</th><th>The Dee</th><th>Wadden Sea</th><th>Sum</th></tr></thead><tbody><tr><td>LaAF</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>4</td><td>1</td><td>0</td><td>0</td><td>0</td><td>5</td></tr><tr><td>LaAG</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>LaAH</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>LaAI</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>LaAJ</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>1</td></tr><tr><td>LaAK</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>1</td></tr><tr><td>LaAL</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>1</td><td>2</td></tr><tr><td>LaAM</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>LaAN</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>0</td><td>3</td></tr><tr><td>LaAO</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>1</td></tr><tr><td>LaAP</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>1</td></tr><tr><td>LaAQ</td><td>0</td><td>0</td><td>0</td><td>0</td><td>7</td><td>2</td><td>12</td><td>4</td><td>1</td><td>0</td><td>16</td><td>42</td></tr><tr><td>LaAR</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>LaAS</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>LaAT</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>1</td><td>2</td></tr><tr><td>LaAV</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>1</td></tr><tr><td>LaAW</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>LaAX</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>LaAY</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>LaAZ</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>LaBA</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>LaBB</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>LaBC</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>BmA</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>BmB</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>BmC</td><td>0</td><td>1</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>BmD</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>BmE</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>BmF</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>BmG</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>BmH</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>Bml</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td></tr></tbody></table>
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Haplotype network of *Bucephalus minimus* based on genetic distance (number of base pair differences) of cytochrome c oxidase subunit 1 (COI) gene sequences across 11 beds of *Cerastoderma edule* where *Bucephalus minimus* specimens were collected. Different haplotypes with respective names are represented by circles, with circle size proportional to observed frequency. Inferred mutation steps are shown by black dots. Colours depict sample location.</figcaption></figure>

may also result in a rise in co-infections, which may heighten the host's susceptibility and mortality. Theoretical and empirical studies provide conflicting findings, and it is therefore unclear whether susceptibility and subsequent host death are directly connected to co-infection by conspecific parasites (Read and Taylor, 2001; Davies et al., 2002; Alizon and van Baalen, 2008). Regardless, a study conducted using different strains of the trematode *S. mansoni* in their host showed an increase in overall pathogenicity (Davies et al., 2002).

As an alternative and excluding the above scenarios, it is possible that after parasite settlement, co-infection becomes unlikely due to alterations in host characteristics. For instance, the nematode *Acanthocheilonema viteae* induces an immune response by the host to further infection by free-living larvae, with no effect on the nematode's persistence within the host (Rajakumar et al., 2006). The extensive use of host tissues, particularly the gonad and digestive gland (Dubois et al., 2009), could be another

reason for the absence of co-infection. This exhaustive use of resources might prevent the settlement of other parasite individuals. Interestingly, co-infection with another trematode species at the same sporocyst stage is exceedingly rare and less frequent than expected (in terms of probability), without predominance by either species (Magalhães et al., 2020). The exhaustive use of tissue preventing the establishment of further miracidia seems to be the most likely explanation for our findings. Thus, the present results cannot fully demonstrate whether co-infection by conspecifics is present or not due to the limited sample size and the genetic marker used. They do, however, suggest that co-infections might be rare and not the rule, and that, for a co-infection to occur, different miracidia must infect simultaneously or within a short period of time. Further studies using codominant nuclear DNA molecular markers, such as microsatellites or SNPs (unavailable to date for *B. minimus*), or attempts at experimental infection of previously infected cockles, would help
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Maximum likelihood phylogenetic tree based on cytochrome c oxidase subunit I (COI) gene sequences of *Bucephalus minimus* haplotypes from this study and retrieved from GenBank database. *Rhipodocotyle* sp. and *Himasthla quissetensis* haplotypes were used as outgroups. Numbers at nodes represent the percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 replicates) using neighbour-joining and maximum likelihood methods (NJ/ML). The scale bar indicates the distances in substitutions per nucleotide.</figcaption></figure>

to determine whether co-infections are possible and confirm our
results.

Bucephalus minimus phylogeography

The COI sequences amplified in this study, as well as those avail-
able from the literature (Feis *et al.*, 2015), were all grouped
together, with short branches separating each sequence, indicating
that the sequences all belong to the same species throughout the
entire cockle distributional range, with no evidence of any
subspecies.

The haplotype network constructed resulted in a complex pat-
tern with several different haplotypes, generally with few muta-
tional steps between them (at most 5 steps between 2 adjacent
haplotypes) in accordance with the previous study by Feis *et al.*
(2015). These results are typical of a population with relatively
stable age structure and a fixed ratio of natural growth (i.e. stable
demography) (Loewe and Hill, 2010). The absence of shared hap-
lotypes among the various beds (most haplotypes were present in
a single bed) supports this conclusion. These results were
expected and consistent with previous reports on this species
(Feis *et al.*, 2015). Trematodes have a complex life cycle that
involves a vertebrate, usually a fish or a bird, as definitive host
(Cribb *et al.*, 2003). In the case of *B. minimus*, its final host is
the European seabass, which has limited migratory capacity com-
pared to birds (de Pontual *et al.*, 2023). This limited mobility
results in a higher degree of isolation between populations of
*B. minimus*.

In our study, 2 haplogroups were identified, one located south of the Bay of Biscay (including Merja Zerga, Aveiro, Baiona and Noia) and the other to the north (including Bay of Somme, English Channel, Celtic Sea, Burry Inlet, The Dee and Wadden Sea). These groups did not share haplotypes except one (LaE). However, this division also coincides with the cockle beds where a lower number of individuals were assessed (see Table 3). Besides, it is noteworthy that there is a core haplotype present in 6 different beds from which most of the other haplotypes diverge. This haplotype has only been identified in the northernmost beds. However, it may represent an ancestral haplotype that has spread across the different regions, and then undergone evolutionary mutations that are specific to each bed in which it occurs. For example, the haplotype network created for the *B. minimus* genetic variability at host level study (Fig. 2) revealed higher similarity of Aveiro (south group) and Arcachon (central/north group) haplotypes than Aveiro and Baiona (south group) haplotypes due to a parsimonious position (position 50 of our alignment) in the COI sequences where a cytosine replaces a thymine in all Baiona haplotypes. Consequently, the presence or absence of haplotypes in a particular bed should be interpreted cautiously, especially for southern beds.

Despite the limitations previously described, the population
genetics of *B. minimus* exhibited similar geographic clusters as
those of its first intermediate host, *C. edule* [structured as north-
wards and southwards of French Brittany (Vera *et al.*, 2022)].
Cockle dispersal mainly occurs during their larval stages
(Martel and Chia, 1991), when *B. minimus* cannot infect them
(Magalhães et al., 2015). It may also occur through human-mediated movements, although no information is available regarding this possibility. Therefore, it is more logical to expect a correlation with the population genetic structure of seabass, the most mobile host in the trematode life cycle (Zemmer et al., 2020). Two distinct genetic population units, one in the Atlantic area and another in the Mediterranean, have been previously identified for seabass (Souche et al., 2015), with a slight genetic differentiation in the Atlantic area observed south of the Strait of Gibraltar (Morocco) attributed to a hybrid zone between 2 evolutionary lineages (Lemaire et al., 2005; Vandeputte et al., 2019). A genomic survey conducted with more than 2700 molecular markers (i.e. SNPs) throughout the Atlantic area revealed the presence of 3 different groups weakly differentiated and geographically distributed. Thus, all Atlantic wild fish belonged to a single group, except specimens from the northern North Sea (i.e. Norway) and the Strait of Gibraltar (AQUATRACE, 2017), matching the results of Souche et al. (2015). This structure (see Fig. S1) could explain the high B. minimus population differentiation found when Merja Zerga (located on the Moroccan coast) was compared with the remaining beds, although more beds south of the Strait of Gibraltar need to be studied, along with greater sample sizes from southern beds, especially in the Iberian Peninsula, to corroborate this explanation. Following seabass genetic structure, the other cockle beds analysed may be included within the same genetic group. Despite this homogeneity, Quéré et al. (2010) identified genetic differentiation between seabass from the Bay of Biscay and North Sea using a molecular marker (microsatellite) under selection associated to the somatolactin gene, although this pattern was not confirmed either by other candidate genes or adaptive variation screenings (Souche et al., 2015; AQUATRACE, 2017). Moreover, an electronic tagging study in seabass has identified the Bay of Biscay as a potential hybridization zone for several subpopulations of seabass from various Atlantic coast locations (de Pontual et al., 2023), which could account for the region's high haplotype diversity and differentiation in the population genetics of B. minimus. In any event, the effect of selective processes on seabass population structure cannot be ruled out and could explain the current genetic structure in B. minimus, with the Bay of Biscay (here represented by Arcachon bed) being a possible hybrid region between northern and southern groups. Thus, the genetic structure among B. minimus population uncovered in our study may reflect local adaptation of the parasite to the most common host genotypes occurring in each of the regions sampled (e.g. Sasal et al., 2000). One possible approach to further investigate this possibility would be to conduct co-phylogeographic analyses between the genetic structure of the parasite and that of their cockle and seabass hosts (e.g. Nieberding et al., 2004; Nieberding and Olivieri, 2007).

In conclusion, with a single COI haplotype observed per host, B. minimus genetic variability at the host level was very limited. This suggests that a single miracidium may be infecting C. edule, possibly indicating the existence of strong mechanisms operating in the background to reduce the likelihood of multiple infections. However, the processes involved are currently poorly understood, with a long way to go to fully comprehend trematode host-parasite interactions. Particularly, when discussing the early stage of the parasite's life cycle (sporocysts), laboratory studies are required to confirm the likelihood of co-infections and the defence mechanisms involved. Similar to the population structure of its first host, C. edule, 2 B. minimus groups (north and south of Bay of Biscay) were found. It is probable that this pattern derives from D. labrax, this parasite's final host, when specific genomic areas under selection are investigated. However, this structure may also be the result of sampling limitations mostly in the

southern beds. The full parasite population connectivity will only be revealed by additional research, specifically by increasing the number of samples of understudied beds, such as in the Mediterranean.

**Supplementary material.** The supplementary material for this article can be found at https://doi.org/10.1017/S0031182023000987

**Data availability statement.** The data that support the findings of this manuscript are provided in the text and are also available from the corresponding authors (S.C. and M.V.) upon reasonable request.

**Acknowledgements.** The authors are grateful to Txetxu Santiago and Leslie Stout for their help in cockle sampling and Cristina Alonso for her technical support. We are also grateful to Professor J. Russell Stothard and the anonymous reviewers for constructive comments on an earlier version. Simão Correia benefited from PhD grant (2020.04688.BD) funded by National Funds through the Portuguese Science Foundation (Fundação para a Ciência e a Tecnologia, FCT). Luísa Magalhães is funded by national funds (OE), through FCT (grant 2021.01858.CEECIND). Thanks are also due to CESAM (UIDP/50017/2020 + UIDB/50017/2020 + LA/P/0094/2020). Sampling in Arcachon (France) was performed thanks to Planula 4 vessel (CNRS-INSU, Flotte Océanographique Française) and its staff.

**Author contributions.** S.C., S.F.-B., L.M. and M.V. conceived and designed the study. S.C. conducted data gathering. X.d.M. and G.D. contributed with data. S.C. and M.V. performed data analyses. S.C. wrote the original draft. S.F.-B., L.M., X.d.M., R.P. and M.V. reviewed and edited the draft. S.F.-B., L.M., X.d.M., G.D. and M.V. provided resources. S.F.-B., L.M., R.P. and M.V. supervised the work.

**Financial support.** The COCKLES project (EAPA_458/2016 COCKLES Co-Operation for Restoring Cockle Shellfisheries and its Ecosystem Services in the Atlantic Area) contributed with samples for the phylogeographic study. The project was funded by the Interreg Atlantic Area Programme through the European Regional Development Fund.

**Competing interests.** None.

**Ethical standards.** Not applicable.

## References

* **Abdul-Muneer PM** (2014) Application of microsatellite markers in conservation genetics and fisheries management: recent advances in population structure analysis and conservation strategies. *Genetics Research International* **2014**, 691759. doi: 10.1155/2014/691759

* **Agola LE, Steinauer ML, Mburu DN, Mungai BN, Mwangi IN, Magoma GN, Loker ES and Mijoji GM** (2009) Genetic diversity and population structure of *Schistosoma mansoni* within human infrapopulations in Mwea, central Kenya assessed by microsatellite markers. *Acta Tropica* **111**, 219–225.

* **Alizon S and van Baalen M** (2008) Multiple infections, immune dynamics, and the evolution of virulence. *American Naturalist* **172**, E150–E168.

* **AQUATRACE** (2017) The development of tools for tracing and evaluating the genetic impact of fish from aquaculture: ‘AquaTrace’. CORDIS EU research accessed 20 March 2023. Available at https://cordis.europa.eu/project/id/311920/reporting

* **Baiocchi T, Lee G, Choe DH and Dillman AR** (2017) Host seeking parasitic nematodes use specific odors to assess host resources. *Scientific Reports* **7**, 6270.

* **Balmer O, Stearns SC, Schotzau A and Brun R** (2009) Intraspecific competition between co-infecting parasite strains enhances host survival in African trypanosomes. *Ecology* **90**, 3367–3378.

* **Bartoli P and Gibson DI** (2007) Synopsis of the life cycles of Digenea (Platyhelminthes) from lagoons of the northern coast of the western Mediterranean. *Journal of Natural History* **41**, 1553–1570.

* **Bell AS, De Roode JC, Sim D and Read AF** (2006) Within-host competition in genetically diverse malaria infections: parasite virulence and competitive success. *Evolution* **60**, 1358–1371.

* **Bowers EA** (1969) *Cercariae Bucephalosis haimeana* (Lacaze-Duthiers, 1854) (Digenea – Bucephalidae) in cockle, *Cardium edule* L. in south Wales. *Journal of Natural History* **3**, 409–422.
Burman M (1982) *Neoaplectana carpocapsae*: toxin production by axenic insect parasitic nematodes. *Nematologica* **28**, 62–70.

Carballal MJ, Iglesias D, Santamarina J, Ferro-Soto B and Villalba A (2001) Parasites and pathologic conditions of the cockle *Cerastoderma edule* populations of the coast of Galicia (NW Spain). *Journal of Invertebrate Pathology* **78**, 87–97.

Carlson CJ, Dallas TA, Alexander LW, Phelan AL and Phillips AJ (2020) What would it take to describe the global diversity of parasites? *Proceedings of the Royal Society B-Biological Sciences* **287**, 20201841. doi: 10.1098/rspb.2020.1841

Correia S, Magalhães L, Freitas R, Bazairi H, Gam M and de Montaudouin X (2020) Large scale patterns of trematode parasite communities infecting *Cerastoderma edule* along the Atlantic coast from Portugal to Morocco. *Estuarine Coastal and Shelf Science* **233**, 106546. doi: 10.1016/j.ecss.2019.106546

Cribb TH, Bray RA, Olson PD and Littlewood DTJ (2003) Life cycle evolution in the Digenea: a new perspective from phylogeny. *Advances in Parasitology* **54**, 197–254.

Criscione CD and Blouin MS (2006) Minimal selfing, few clones, and no among-host genetic structure in a hermaphroditic parasite with asexual larval propagation. *Evolution* **60**, 553–562.

Curtis LA, Kinley JL and Tanner NL (2000) Longevity of oversized individuals: growth, parasitism, and history in an estuarine snail population. *Journal of the Marine Biological Association of the United Kingdom* **80**, 811–820.

Davies CM, Fairbrother E and Webster JP (2002) Mixed strain schistosome infections of snails and the evolution of parasite virulence. *Parasitology* **124**, 31–38.

de Montaudouin X, Thieltgès DW, Gam M, Krakau M, Pina S, Bazairi H, Dabouineau L, Russell-Pinto F and Jensen KT (2009) Digenean trematode species in the cockle *Cerastoderma edule*: identification key and distribution along the north-eastern Atlantic shoreline. *Journal of the Marine Biological Association of the United Kingdom* **89**, 543–556.

de Montaudouin X, Blanchet H, Desclaux-Marchand C, Lavesque N and Bachelet G (2016) Cockle infection by *Himasthla quissetensis* - I. From certa-ria emergence to metacercariae infection. *Journal of Sea Research* **113**, 99–107.

de Montaudouin X, Arzul I, Cao A, Carballal MJ, Chollet B, Correia S, Cuesta J, Culloty S, Daffe G, Darriba S, Díaz S, Engelsma M, Freitas R, Garcia C, Goedknegt A, Gonzalez P, Grade A, Groves E, Iglesias D, Jensen KT, Joaquim S, Lynch S, Magalhães L, Mahony K, Maia F, Malham S, Matias D, Nowaczyk A, Ruano F, Thieltgès D and Villalba A (2021) *Parasites and Diseases of the Common Cockle *Cerastoderma edule**, 1st Edn. Aveiro: UA Editora-Universidade de Aveiro.

de Pontual H, Heerah K, Goossens J, Garren F, Martin S, Le Ru L, Le Roy D and Woillez M (2023) Seasonal migration, site fidelity, and population structure of European seabass (*Dicentrarchus labrax*). *ICES Journal of Marine Science* **80**, 1606–1618.

Desclaux C, de Montaudouin X and Bachelet G (2002) Cockle emergence at the sediment surface: ‘favourization’ mechanism by digenean parasites? *Diseases of Aquatic Organisms* **52**, 137–149.

Dubois SY, Savoye N, Sauriau PG, Billy I, Martinez P and de Montaudouin X (2009) Digenean trematodes-marine mollusc relationships: a stable isotope study. *Diseases of Aquatic Organisms* **84**, 65–77.

Dumont M, Mone H, Mouahid G, Idris MA, Shaban M and Boissier J (2007) Influence of pattern of exposure, parasite genetic diversity and sex on the degree of protection against reinfection with *Schistosoma mansoni*. *Parasitology Research* **101**, 247–252.

Excoffier L, Laval G and Schneider S (2005) Arlequin (version 3.0): an integrated software package for population genetics data analysis. *Evolutionary Bioinformatics Online* **1**, 47–50.

Feis ME, Thieltgès DW, Olsen JL, de Montaudouin X, Jensen KT, Bazairi H, Culloty SC and Luttikhuizen PC (2015) The most vague host as the main determinant of population connectivity in marine macroparasites. *Marine Ecology Progress Series* **520**, 85–99.

Fredensborg BL and Poulin R (2005) Larval helminths in intermediate hosts: does competition early in life determine the fitness of adult parasites? *International Journal for Parasitology* **35**, 1061–1070.

Fredensborg BL, Mouritsen KN and Poulin R (2005) Impact of trematodes on host survival and population density in the intertidal gastropod *Zeacumantus subcarinatus*. *Marine Ecology Progress Series* **290**, 109–117.

Goedknegt MA, Feis ME, Wegner KM, Luttikhuizen PC, Buschbaum C, Camphuysen K, van der Meer J and Thieltgès DW (2016) Parasites and marine invasions: ecological and evolutionary perspectives. *Journal of Sea Research* **113**, 11–27.

Goujon M, McWilliam H, Li W, Valentin F, Squizzato S, Paern J and Lopez R (2010) A new bioinformatics analysis tools framework at EMBL-EBI. *Nucleic Acids Research* **38**, W695–W699.

Hall TA (1999) *BioEdit*: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. *Nucleic Acids Symposium Series* **41**, 95–98.

Hatcher MJ and Dunn AM (2011) *Parasites in Ecological Communities: From Interactions to Ecosystems*. Cambridge, UK: Cambridge University Press.

Intecmar (2021) *Informe Epidemiológico de moluscos bivalvos de Galicia, ano 2021*. Intecmar - Instituto Tecnolóxico para o controlo do medio marino de Galicia, Unidade de Patoloxía, Galicia, Spain. Available at http://www.intecmar.gal/Informacion/Patoloxia/Default.aspx

Johnson PTJ, Dobson A, Lafferty KD, Marcogliese DJ, Memmott J, Orlofske SA, Poulin R and Thieltgès DW (2010) When parasites become prey: ecological and epidemiological significance of eating parasites. *Trends in Ecology & Evolution* **25**, 362–371.

Jokela J, Uotila L and Taskinen J (1993) Effect of the castrating trematode parasite *Rhipidocotyle fennica* on energy allocation of fresh-water clam *Anodonta piscinalis*. *Functional Ecology* **7**, 332–338.

Keeney DB, Waters JM and Poulin R (2007) Clonal diversity of the marine trematode *Maritrema novaezealandensis* within intermediate hosts: the molecular ecology of parasite life cycles. *Molecular Ecology* **16**, 431–439.

Keeney DB, Boessenkool S, King TM, Leung TLF and Poulin R (2008) Effects of interspecific competition on asexual proliferation and genetic diversity in larval trematode infections of snails. *Parasitology* **135**, 741–747.

Kumar S, Stecher G, Li M, Knyaz C and Tamura K (2018) MEGA X: molecular evolutionary genetics analysis across computing platforms. *Molecular Biology and Evolution* **35**, 1547–1549.

Lagrue C, McEwan J, Poulin R and Keeney DB (2007) Co-occurrences of parasite clones and altered host phenotype in a snail–trematode system. *International Journal for Parasitology* **37**, 1459–1467.

Lemaire C, Versini JJ and Bonhomme F (2005) Maintenance of genetic differentiation across a transition zone in the sea: discordance between nuclear and cytoplasmic markers. *Journal of Evolutionary Biology* **18**, 70–80.

Librado P and Rozas J (2009) DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. *Bioinformatics* **25**, 1451–1452.

Loewe L and Hill WG (2010) The population genetics of mutations: good, bad and indifferent. *Philosophical Transactions of the Royal Society B-Biological Sciences* **365**, 1153–1167.

Louhi KR, Karvonen A, Rellstab C, Louhi R and Jokela J (2013) Prevalence of infection as a predictor of multiple genotype infection frequency in parasites with multiple-host life cycle. *Journal of Animal Ecology* **82**, 191–200.

Magalhães L, Freitas R and de Montaudouin X (2015) Review: *Bucephalus minimus*, a deleterious trematode parasite of cockles *Cerastoderma* spp. *Parasitology Research* **114**, 1263–1278.

Magalhães L, Freitas R, Dairain A and de Montaudouin X (2017) Can host density attenuate parasitism? *Journal of the Marine Biological Association of the United Kingdom* **97**, 497–505.

Magalhães L, Correia S, de Montaudouin X and Freitas R (2018) Spatio-temporal variation of trematode parasites community in *Cerastoderma edule* cockles from Ria de Aveiro (Portugal). *Environmental Research* **164**, 114–123.

Magalhães L, Daffe G, Freitas R and de Montaudouin X (2020) *Monorchis parvus* and *Gymnophallus choledochus*: two trematode species infecting cockles as first and second intermediate host. *Parasitology* **147**, 643–658.

Martel A and Chia FS (1991) Drifting and dispersal of small bivalves and gastropods with direct development. *Journal of Experimental Marine Biology and Ecology* **150**, 131–147.

Mideo N (2009) Parasite adaptations to within-host competition. *Trends in Parasitology* **25**, 261–268.

Morley NJ (2012) Thermodynamics of miracidial survival and metabolism. *Parasitology* **139**, 1640–1651.

Morley NJ and Lewis JW (2017) Thermodynamics of egg production, development and hatching in trematodes. *Journal of Helminthology* **91**, 284–294.

Moszczynska A, Locke SA, McLaughlin JD, Marcogliese DJ and Crease TJ (2009) Development of primers for the mitochondrial cytochrome c oxidase I gene in digenetic trematodes (Platyhelminthes) illustrates the challenge of barcoding parasitic helminths. *Molecular Ecology Resources* **9**, 75–82.
Mouritsen KN and Poulin R (2002) Parasitism, community structure and biodiversity in intertidal ecosystems. *Parasitology* **124**, S101–S117.

Nieberding C and Olivieri I (2007) Parasites: proxies for host genealogy and ecology? *Trends in Ecology & Evolution* **22**, 156–165.

Nieberding C, Morand S, Libois R and Michaux JR (2004) A parasite reveals cryptic phylogeographic history of its host. *Proceeding of the Royal Society B* **271**, 2559–2568.

Paradis E (2010) pegas: an R package for population genetics with an integrated-modular approach. *Bioinformatics* **26**, 419–420.

Paradis E and Schliep K (2019) ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R. *Bioinformatics* **35**, 526–528.

Pina S, Barandela T, Santos MJ, Russell-Pinto F and Rodrigues P (2009) Identification and description of *Bucephalus minimus* (Digenea: Bucephalidae) life cycle in Portugal: morphological, histopathological and molecular data. *Journal of Parasitology* **95**, 353–359.

Poulin R (1999) The functional importance of parasites in animal communities: many roles at many levels? *International Journal for Parasitology* **29**, 903–914.

Poulin R (2001) Interactions between species and the structure of helminth communities. *Parasitology* **122**, S3–S11.

Poulin R (2006) Global warming and temperature-mediated increases in cer- carial emergence in trematode parasites. *Parasitology* **132**, 143–151.

Preston DL, Orlofske SA, Lambden JP and Johnson PTJ (2013) Biomass and productivity of trematode parasites in pond ecosystems. *Journal of Animal Ecology* **82**, 509–517.

Queiroz N, Humphries NE, Noble LR, Santos AM and Sims DW (2012) Spatial dynamics and expanded vertical niche of blue sharks in oceano- graphic fronts reveal habitat targets for conservation. *PLoS ONE* **7**, e32374. doi: 10.1371/journal.pone.0032374

Quéré N, Guinand B, Kuhl H, Reinhardt R, Bonhomme F and Desmarais E (2010) Genomic sequences and genetic differentiation at associated tandem repeat markers in growth hormone, somatolactin and insulin-like growth factor-1 genes of the sea bass, *Dicentrarchus labrax*. *Aquatic Living Resources* **23**, 285–296.

Rajakumar S, Bleiss W, Hartmann S, Schierack P, Marko A and Lucius R (2006) Concomitant immunity in a rodent model of filariasis: the infection of *Meriones unguiculatus* with *Acanthocheilonema viteae*. *Journal of Parasitology* **92**, 41–45.

Rauch G, Kalbe M and Reusch TBH (2005) How a complex life cycle can improve a parasite’s sex life. *Journal of Evolutionary Biology* **18**, 1069–1075.

Read AF and Taylor LH (2001) The ecology of genetically diverse infections. *Science* **292**, 1099–1102.

Rice P, Longden I and Bleasby A (2000) EMBOSS: the European molecular biology open software suite. *Trends in Genetics* **16**, 276–277.

Sasal S, Durand P, Faliex E and Morand S (2000) Experimental approach to the importance of parasitism in biological conservation. *Marine Ecology Progress Series* **198**, 293–302.

Schliep K (2011) phangorn: phylogenetic analysis in R. *Bioinformatics* **27**, 592–593.

Schulte-Oehlmann U, Oehlmann J, Fioroni P and Bauer B (1997) Imposedx and reproductive failure in *Hydrobia ulvae* (Gastropoda: Prosobranchia). *Marine Biology* **128**, 257–266.

Selbach C and Poulin R (2020) Some like it hotter: trematode transmission under changing temperature conditions. *Oecologia* **194**, 745–755.

Selva L, Viana D, Regev-Yochay G, Trzcinski K, Corpa JM, Lasa I, Novick RP and Penades JR (2009) Killing niche competitors by remote-control bacteriophage induction. *Proceedings of the National Academy of Sciences of the USA* **106**, 1234–1238.

Souche EL, Hellemans B, Babbucci M, MacAoidh E, Guinard B, Bargelloni L, Chistiakov DA, Patarnello T, Bonhomme F, Martinsohn JT and Volckaert FAM (2015) Range-wide population structure of European sea bass *Dicentrarchus labrax*. *Biological Journal of the Linnean Society* **116**, 86–105.

Theron A, Sire C, Rognon A, Prugnolle F and Durand P (2004) Molecular ecology of *Schistosoma mansoni* transmission inferred from the genetic composition of larval and adult infrapopulations within intermediate and definitive hosts. *Parasitology* **129**, 57–585.

Thieltges DW (2004) Effects of trematode parasites on mortality in the common cockle *Cerastoderma edule*. *International Journal of Medical Microbiology* **293**, 57–57.

Thieltges DW (2006) Parasite induced summer mortality in the cockle *Cerastoderma edule* by the trematode *Gymnophallus choledochus*. *Hydrobiologia* **559**, 455–461.

Thieltges DW and Rick J (2006) Effect of temperature on emergence, survival and infectivity of cercariae of the marine trematode *Renicola roscovita* (Digenea: Renicolidae). *Diseases of Aquatic Organisms* **73**, 63–68.

Vandeputte M, Gagnaire PA and Allal F (2019) The European sea bass: a key marine fish model in the wild and in aquaculture. *Animal Genetics* **50**, 195–206.

Vera M, Maroso F, Wilmes SB, Hermida M, Blanco A, Fernández C, Groves E, Malham SK, Bouza C, Cockle’s Consortium, Robins PE and Martínez P (2022) Genomic survey of edible cockle (*Cerastoderma edule*) in the Northeast Atlantic: a baseline for sustainable management of its wild resources. *Evolutionary Applications* **15**, 262–285.

Vilas R, Paniagua E and Sanmartin ML (2003) Genetic variation within and among infrapopulations of the marine digenetic trematode *Lecithochirium fusiforme*. *Parasitology* **126**, 465–472.

Yu G, Smith D, Zhu H, Guan Y and Lam TT (2017) ggtree: and R package for visualization and annotation of phylogenetic trees with their covariates and other associated data. *Methods in Ecology and Evolution* **8**, 28–36.

Zemmer SA, Detwiler JT, Sokol ER, Da Silva Neto JG, Wyderko J, Potts K, Gajewski ZJ, Sarment LV, Benfield EF and Belden LK (2020) Spatial scale and structure of complex life cycle trematode parasite communities in streams. *PLOS ONE* **15**, e0241973. doi: 10.1371/journal.pone.0241973