Introduction
In coastal ecosystems, trematodes are the most abundant and common metazoan parasites (Lauckner, Reference Lauckner and Kinne1983; Sousa, Reference Sousa1991; Mouritsen and Poulin, Reference Mouritsen and Poulin2002). These macroparasites are exclusively endoparasites and have a complex and heteroxenous life cycle, generally involving three hosts and exhibiting alternation between asexual multiplication and sexual reproduction phases (Esch, Reference Esch2002; Bartoli and Gibson, Reference Bartoli and Gibson2007). The adult stage of these parasites reproduces sexually in the final host, which is a vertebrate (generally a fish or a shorebird). Eggs are released into the environment (through final host feces) and either evolve into miracidium, a free-living stage, to infect the first intermediate host (usually a mollusc), or hatches in miracidium after they have been ingested by the first intermediate host. Each larva develops into a sac-like sporocyst or a redia, depending on the trematode species, which will asexually produce cercariae, a second free-living stage. These cercariae emerge from the first intermediate host and swim actively to penetrate the second intermediate host (a vertebrate or an invertebrate) and settle as a metacercaria, a latent stage. When the second host is predated by the final host, metacercaria transforms into the adult stage, achieving the life cycle.
Molluscs are the common first and second intermediate hosts of trematode parasites and almost all known bivalves are parasitized, with predominant infection by metacercariae compared to infection by sporocysts (Lauckner, Reference Lauckner and Kinne1983; Sousa, Reference Sousa1991; Galaktionov and Dobrovolskij, Reference Galaktionov and Dobrovolskij2003). More particularly, the edible cockle Cerastoderma edule (Linnaeus, 1758) is one of the most widespread and abundant bivalves in soft bottom shallow coastal ecosystems along the northeast Atlantic (Malham et al., Reference Malham, Hutchinson and Longshaw2012). They are suitable hosts for harbouring one of the highest diversities of trematode species (Krakau et al., Reference Krakau, Thieltges and Reise2006; Thieltges et al., Reference Thieltges, Krakau, Andresen, Fottner and Reise2006), with 16 known species (de Montaudouin et al., Reference de Montaudouin, Arzul, Cao, Carballal, Bruno, Correia, Cuesta, Culloty, Daffe, Darriba, Díaz, Freitas, Garcia, Goedknegt, Grade, Groves, Iglesias, Jensen and Villalba2021). Living buried a few centimetres into the sediment, this bivalve is a key species in coastal ecosystem functioning (Carss et al., Reference Carss, Brito, Chainho, Ciutat, de Montaudouin, Fernández Otero, Incera Filgueira, Grabutt, Goedknegt, Lynch, Mahony, Maire, Malham, Orvain, van der Schatte Olivier and Jones2020). In particular, due to their bioturbation (i.e. biomixing of the sediment) and biodeposition activities, they modulate the physical properties and biogeochemical dynamics of the sediment (Ciutat et al., Reference Ciutat, Widdows and Pope2007; Eriksson et al., Reference Eriksson, Westra, van Gerwen, Weerman, van der Zee, van der Heide, van de Koppel, Olff, Piersma and Donadi2017; Dairain et al., Reference Dairain, Maire, Meynard, Richard, Rodolfo-Damiano and Orvain2020b), and they have an important role in connection between trophic levels (Rakotomalala et al., Reference Rakotomalala, Grangeré, Ubertini, Forêt and Orvain2015). Finally, through their filtration activity, cockles can regulate phytoplankton biomass and turbidity (Cloern, Reference Cloern1982; Newell, Reference Newell2004). The effect of trematodes is not restricted to cockle activity, but may significantly alter their fitness, with a subsequent impact at the population scale. Indeed, parasites can significantly contribute to C. edule mortality and population decline (Burdon et al., Reference Burdon, Callaway, Elliott, Smith and Wither2014). The influence of trematode parasites on cockle survival is highly dependent on species, infection intensity (i.e. number of parasites per infected cockle) or abundance (i.e. number of parasites per cockle, infected or not), prevalence (i.e. percentage of infected cockles) and parasitic stage (Lauckner, Reference Lauckner and Kinne1983). When cockles are the first intermediate hosts, the effects are particularly deleterious. For example, Bucephalus minimus sporocysts invade most of the tissues, including the gonads, digestive gland and gills (Dubois et al., Reference Dubois, Savoye, Sauriau, Billy, Martinez and de Montaudouin2009). This invasion leads to castration, starvation, reduction of the cockle growth rate and condition index, as well as modulation of their impact on sediment erodibility (de Montaudouin et al., Reference de Montaudouin, Thieltges, Gam, Krakau, Pina, Bazairi, Dabouineau, Russell-Pinto and Jensen2009, Reference de Montaudouin, Arzul, Cao, Carballal, Bruno, Correia, Cuesta, Culloty, Daffe, Darriba, Díaz, Freitas, Garcia, Goedknegt, Grade, Groves, Iglesias, Jensen and Villalba2021; Magalhães et al., Reference Magalhães, Freitas and de Montaudouin2015; Dairain et al., Reference Dairain, Maire, Meynard and Orvain2020a). In addition, Gymnophallus choledochus, also using C. edule as first (and second) intermediate host, occupies the entire mantle cavity, causing gonad structure loss and mass mortality (Thieltges, Reference Thieltges2006b; Magalhães et al., Reference Magalhães, Daffe, Freitas and de Montaudouin2020a). Nevertheless, prevalence is usually low and the effect on the cockle population scale is considered moderate (Thieltges et al., Reference Thieltges, de Montaudouin, Fredensborg, Jensen, Koprivnikar and Poulin2008; de Montaudouin et al., Reference de Montaudouin, Thieltges, Gam, Krakau, Pina, Bazairi, Dabouineau, Russell-Pinto and Jensen2009, Reference de Montaudouin, Arzul, Cao, Carballal, Bruno, Correia, Cuesta, Culloty, Daffe, Darriba, Díaz, Freitas, Garcia, Goedknegt, Grade, Groves, Iglesias, Jensen and Villalba2021). Metacercariae have a more limited impact since they do not multiply in their second intermediate host tissues. For instance, Renicola roscovitus is one of the dominant metacercariae encysting in cockle palps (Krakau et al., Reference Krakau, Thieltges and Reise2006; Lassalle et al., Reference Lassalle, de Montaudouin, Soudant and Paillard2007). Its impact on cockles has been reported as moderate, with a slight reduction in oxygen consumption, no impact on the condition index, low antioxidant defence activation and an intermediate level of cellular damage (Magalhães et al., Reference Magalhães, Freitas and de Montaudouin2020b). However, other species of trematodes infecting cockles as metacercariae can have more deleterious impacts on their host population, especially when the prevalence and intensity become high. For example, Gymnophallus minutus (de Montaudouin et al., Reference de Montaudouin, Kisielewski, Bachelet and Desclaux2000; Thieltges and Reise, Reference Thieltges and Reise2006a; Gam et al., Reference Gam, Bazaïri, Jensen and de Montaudouin2008; Fermer et al., Reference Fermer, Culloty, Kelly and O'Riordan2010) causes pathology in cockles, modifies their behaviour (emerging at the sediment–water interface) and provokes significant mortality (Bowers et al., Reference Bowers, Bartoli, Russell-Pinto and James1996; Gam et al., Reference Gam, de Montaudouin and Bazairi2009b; Fermer et al., Reference Fermer, Culloty, Kelly and O'Riordan2011a).
The Himasthla genus occurs along the northeastern Atlantic coasts (James, Reference James1968; Blakeslee and Byers, Reference Blakeslee and Byers2008, Galaktionov et al., Reference Galaktionov, Solovyeva and Miroliubov2021), and four species constitute the subject of this review (i.e. Himasthla continua, Himasthla elongata, Himasthla interrupta and Himasthla quissetensis). This is one of the most prevalent, abundant and widespread trematode genera infecting cockles as the second intermediate host (Thieltges and Reise, Reference Thieltges and Reise2006a; Gam et al., Reference Gam, Bazaïri, Jensen and de Montaudouin2008; de Montaudouin et al., Reference de Montaudouin, Arzul, Cao, Carballal, Bruno, Correia, Cuesta, Culloty, Daffe, Darriba, Díaz, Freitas, Garcia, Goedknegt, Grade, Groves, Iglesias, Jensen and Villalba2021). The aims of this review were: (1) to compile literature concerning these parasites in cockles and to summarize the main findings; (2) to provide a molecular signature with the potential to accompany stereomicroscope morphological identification and (3) to analyse a 20-year long-term database concerning a cockle population and its associated trematode species in Banc d'Arguin, France, in order to describe the infection patterns of cockles by Himasthla spp. In the latter case, the tested hypotheses were: (1) infection increases with age and with seasonal modulation; (2) infection success may be limited by cockle abundance (dilution effects) and (3) Himasthla species occupy different ecological niches (i.e. different organs in the cockle) and do not compete inside their individual host.
Materials and methods
Literature review
The references gathered in this review were found in Scopus using relevant terms such as ‘Cerastoderma (or Cardium) edule’ and ‘Himasthla’, published before March 2021. The list of articles was restricted to those studies that clearly identified the occurrence of H. continua, H. elongata, H. interrupta and H. quissetensis in C. edule. A reference list of relevant papers was provided and their main findings were summarized. Thus, a total of 46 publications was examined.
Long-term monitoring
Sampling and trematode identification
From November 1997 to October 2018, cockle monitoring was performed in Banc d'Arguin (44°40′N; 1°10′W), a National Nature Reserve in France. The sampled station is an intertidal semi-sheltered sandflat. The sediment is composed of medium sands (grain-size median = 330 μm) (de Montaudouin and Lanceleur, Reference de Montaudouin and Lanceleur2011), the temperature of water fluctuates seasonally between 9.5 and 21.5°C and the salinity is constant (34–35). The tide is semidiurnal (Gassiat, Reference Gassiat1989). Cockles were collected monthly by sampling six 0.25 m2 quadrates sieved with a 1 mm mesh. Cockle shell length was measured to the nearest millimetre with a digital calliper. Cohorts were identified by the analysis of length frequency histograms (Bhattacharya, Reference Bhattacharya1967). Ten cockles per cohort were dissected and squeezed between two glass slides for trematode observation under a stereomicroscope. All trematodes were identified to the species level using morphological criteria (de Montaudouin et al., Reference de Montaudouin, Thieltges, Gam, Krakau, Pina, Bazairi, Dabouineau, Russell-Pinto and Jensen2009, Reference de Montaudouin, Arzul, Cao, Carballal, Bruno, Correia, Cuesta, Culloty, Daffe, Darriba, Díaz, Freitas, Garcia, Goedknegt, Grade, Groves, Iglesias, Jensen and Villalba2021). However, the different species of the Himasthla genus remain difficult to distinguish using morphological analysis and light microscopy. Therefore, several metacercariae were punctually dissected in different cockle tissues, identified morphologically under the microscope (size, number of spines) and then molecularly characterized (see the molecular biology section). Four species were identified in Banc d'Arguin: H. continua, H. elongata, H. interrupta and H. quissetensis. This long-term monthly survey was not based on the dissection of all metacercariae, so that our present analysis was restricted to H. interrupta and H. quissetensis. Indeed, we noticed frequent mistakes concerning stereomicroscope identification between the two other species: H. continua and H. elongata.
Data analysis
During the 20 years, 5820 cockles were analysed (Fig. 1), with shell lengths ranging between 2 and 38 mm. A Spearman test was conducted to investigate the relationship between the prevalence (i.e. the percentage of infected cockles) and the cockle shell length for both Himasthla species. Then, the seasonality of the intensity of infection (i.e. number of metacercariae per infected cockle) was studied. Firstly, cockles’ shell lengths were transformed into relative age using a local Von Bertalanffy growth function (Gam et al., Reference Gam, de Montaudouin and Bazairi2009b):
where t is the relative cockle age (years), Lt is the cockle length at age t (mm), k = 1.5 year−1 and L ∞ = 38.3 mm. Absolute age was deduced from a probable recruitment date in May (de Montaudouin et al., Reference de Montaudouin, Arzul, Cao, Carballal, Bruno, Correia, Cuesta, Culloty, Daffe, Darriba, Díaz, Freitas, Garcia, Goedknegt, Grade, Groves, Iglesias, Jensen and Villalba2021). Then, for each Himasthla species, parasite intensity was compared between months (i.e. between cockle age), using a Wilcoxon non-parametric test, in order to detect significant infection (i.e. increase of the parasite intensity) or parasite-dependent mortality (i.e. decrease of the parasite intensity) processes.
For the following tests, all 1-year-old cockles (corresponding to 15–25 mm) were pooled. This cockle range was selected in order to exclude younger cockles, which are always poorly infected (whatever the environmental conditions are) and older cockles, which do not occur every year. For each Himasthla species, the associated trematode community was compared at the cockle specimen scale with and without this parasite species, through relative abundance (χ 2 test), species richness and metacercariae abundance (Wilcoxon test). The effect of cockle density on Himasthla infection was tested with a Spearman correlation test. Trematodes using cockles as first intermediate hosts were not considered (Bucephalus minutus, G. choledochus, Monorchis parvus), as they have recently been studied in detail (Magalhães et al., Reference Magalhães, Freitas and de Montaudouin2015, Reference Magalhães, Daffe, Freitas and de Montaudouin2020a). All statistical analyses were performed using the open-source program R (v3.6.1) in R studio (v1.3.1056) (www.R-project.org, accessed on 1 August 2020).
Molecular identification: DNA isolation, amplification and sequencing
The cockles were dissected to extract metacercariae from the four Himasthla species that can occur in cockles. Prior to molecular biology analysis, identification was performed based on morphology (metacercariae diameter and number of oral spines) and tissue location (mantle, foot, digestive gland). Metacercariae were assigned to H. quissetensis when 31 oral spines were present (Stunkard, Reference Stunkard1938). Metacercariae were confidently assigned to H. interrupta when they presented 29 oral spines with a diameter <140 μm, and occurred in the mantle margin (de Montaudouin et al., Reference de Montaudouin, Thieltges, Gam, Krakau, Pina, Bazairi, Dabouineau, Russell-Pinto and Jensen2009). A mismatch between H. continua and H. elongata was possible when metacercariae had 29 spines with a diameter >150 μm and occurred in the foot. Periwinkles (Littorina littorea) were also collected as the first intermediate host of H. elongata. These periwinkles were disposed in small dishes at ambient temperature in order to stimulate cercariae emission (Wegeberg et al., Reference Wegeberg, de Montaudouin and Jensen1999). Then, cercariae were sampled with a micropipette for molecular analysis and subsequent comparison with large 29-spine metacercariae found in the cockle foot.
Metacercariae and cercariae were sampled under a stereomicroscope for DNA analysis. For all species, three replicates (i.e. metacercariae) were collected. They were placed in microtubes and immediately frozen at −20°C. DNA extraction was performed using the QIAamp DNA micro kit (QIAGEN, Hilden, Germany), following the protocol supplied by the manufacturer. Using primers Bb18S and Bb18AS for small subunit ribosomal RNA gene (18S) (de Montaudouin et al., Reference de Montaudouin, Bazairi, Mlik and Gonzalez2014), BbITS and BbITAS for internal transcribed spacer 1 (ITS1) (de Montaudouin et al., Reference de Montaudouin, Bazairi, Mlik and Gonzalez2014) and TremCOIS2 and TremCOIAS2 for cytochrome c oxidase subunit I (COI) (Magalhães et al., Reference Magalhães, Daffe, Freitas and de Montaudouin2020a), with sequences given in Table 1, about 530 bp of 18S, 600 bp for ITS1 and 300 bp of COI genes were amplified. The polymerase chain reaction (PCR) was performed with Gotaq G2 Flexi DNA polymerase (PROMEGA, Madison, Wisconsin, USA), with 50 μL mixtures containing: 10 μL of 5× Colorless GoTaq® reaction buffer (final concentration of 1×), 1.5 μL of MgCl2 solution (final concentration of 1.5 mm), 1 μL of PCR nucleotide mix (final concentration of 0.2 mm each dNTP), 0.5 μL of each primer (final concentration of 1 μ m), 0.2 μL of GoTaq® G2 Flexi DNA polymerase (5 U μL−1), 1 μL template DNA and 33.8 μL of nuclease-free water. The temperature profile was 94°C/10 min–(94°C/60 s–59°C/30 s–72°C/90 s) × 40 cycles–72°C/10 min–4°C for 18S and ITS1, and 95°C/10 min–(95°C/60 s–43°C/30 s–72°C/60 s) × 40 cycles–72°C/10 min–4°C for COI. The amplified PCR products were analysed by electrophoresis in a 1% p/v agarose gel stained with ethidium bromide. They were then sent to Eurofins Company for complete double-strain sequencing, using the same set of primers as used for the PCR. Overlapping sequence (forward and reverse) fragments were merged into consensus sequences and aligned using Clustal Omega. For COI, the sequences were translated into amino acid alignments, and checked for stop codons to avoid pseudogenes. All sequences obtained in this study were deposited in GenBank (Table 2).
Nucleotide sequences of specific primer pairs
Accession numbers when DNA sequences were deposited in GenBank, for each gene (18S, ITS1 and COI) and the four Himasthla species
Results
Literature review
Description and life cycle
Himasthla continua (Loos-Frank, 1967), H. elongata (Mehlis, 1831) Dietz, 1909, H. interrupta (Loos-Frank, 1967) and H. quissetensis (Miller & Northup, 1926) Stunkard, Reference Stunkard1938, belong to the Platyhelminthes phylum, Trematoda class, Digenea subclass and Himasthlidae family (Table 3). Himasthla continua, H. elongata and H. interrupta are considered to be native parasites of C. edule, whereas H. quissetensis could have been introduced to Europe from North America (de Montaudouin et al., Reference de Montaudouin, Jensen, Desclaux, Wegeberg and Sajus2005; Longshaw and Malham, Reference Longshaw and Malham2013), and was not reported before 1990 on the eastern Atlantic coastline (Russell-Pinto, Reference Russell-Pinto1993). Himasthla species can be differentiated by morphometric measures, number of spines and to a lesser extent by their location in cockle organs (Russell-Pinto et al., Reference Russell-Pinto, Gonçalves and Bowers2006; de Montaudouin et al., Reference de Montaudouin, Thieltges, Gam, Krakau, Pina, Bazairi, Dabouineau, Russell-Pinto and Jensen2009). Himasthla quissetensis is the only one of these four species with 31 oral spines. Himasthla interrupta displays the smallest cyst diameter (80–140 μm) and is mainly located in the cockle mantle margin (in the anterior edge, at the opposite side of siphons). In contrast, H. elongata displays the largest cysts (210–270 μm). The size of H. continua metacercariae ranges between 150 and 210 μm.
Characteristics of the four studied Himasthla species in terms of target organs, number of oral spines, metacercariae mean diameter and different host species within their life cycle
Geographic distribution, abundance and prevalence
Himasthla spp. metacercariae infect C. edule from Norway to North Africa, with different successes according to the species (Table 4). Himasthla continua and H. interrupta are ubiquitous species, present from Denmark to Morocco. They use the same first intermediate host, Peringia ulvae (Bordalo et al., Reference Bordalo, Ferreira, Jensen and Pardal2011), which is also widely distributed along the east Atlantic shoreline and could explain that they follow similar distribution patterns. However, the intensity of infection and the prevalence are generally higher for H. interrupta than those for H. continua (de Montaudouin et al., Reference de Montaudouin, Thieltges, Gam, Krakau, Pina, Bazairi, Dabouineau, Russell-Pinto and Jensen2009). This can be attributed to the fact that H. continua cercariae have more difficulties penetrating through the cockle inhalant siphon due to their larger dimension (Wegeberg et al., Reference Wegeberg, de Montaudouin and Jensen1999). Himasthla elongata is absent South to Portugal, and the highest abundance occurs in the northern European countries (e.g. Norway) (de Montaudouin et al., Reference de Montaudouin, Thieltges, Gam, Krakau, Pina, Bazairi, Dabouineau, Russell-Pinto and Jensen2009). This range could be related to the distribution of the first intermediate host, the periwinkle L. littorea which is present from Portugal to Norway and Russia (Johannesson, Reference Johannesson1988). In contrast, H. quissetensis is mainly reported in the southern part of the C. edule geographical distribution, with the highest rate of infection in France, Portugal and Morocco. In this case, the first intermediate host is Tritia reticulata, which is widespread from the north to south of Europe (Russell-Pinto et al., Reference Russell-Pinto, Gonçalves and Bowers2006). This lack of direct relationship between host and parasite distribution shows that the abundance of metacercariae also depends on other factors, such as cockle density, size, age and fitness, as well as the ambient benthic community (Gam et al., Reference Gam, de Montaudouin and Bazairi2009b; de Montaudouin and Lanceleur, Reference de Montaudouin and Lanceleur2011; Magalhães et al., Reference Magalhães, Freitas, Dairain and De Montaudouin2017; Welsh et al., Reference Welsh, Hempel, Markovic, van der Meer and Thieltges2019; Correia et al., Reference Correia, Magalhães, Freitas, Bazairi, Gam and de Montaudouin2020a). Moreover, H. quissetensis has also been recorded in different Mediterranean lagoons, infecting a close-related cockle, Cerastoderma glaucum (Prévot, Reference Prévot1974; Bartoli and Gibson, Reference Bartoli and Gibson2007).
Review of the literature regarding H. continua, H. elongata, H. interrupta and H. quissetensis metacercariae infection in C. edule
Abundance was the number of metacercariae per infected or uninfected cockle, intensity was the number of metacercariae per infected cockle and prevalence was the percentage of infected cockles.
Effects of second intermediate hosts
In most studies, the pathogenicity of Himasthla metacercariae is reported as low in C. edule, as this stage is considered energetically inert (Lauckner, Reference Lauckner and Kinne1983). Indeed, laboratory and field experiments have demonstrated that, under moderate infection and normal environmental conditions, H. continua and H. interrupta do not increase cockle mortality (Jensen et al., Reference Jensen, Castro and Bachelet1999; Wegeberg and Jensen, Reference Wegeberg and Jensen2003). Similarly, H. interrupta and H. quissetensis do not impair C. edule shell growth and production (Wegeberg and Jensen, Reference Wegeberg and Jensen2003; Gam et al., Reference Gam, de Montaudouin and Bazairi2009b), and H. elongata has no significant effect on cockle bioturbation activity (sediment reworking and bioirrigation rates) (Richard et al., Reference Richard, de Montaudouin, Rubiello and Maire2021). Nevertheless, when cercariae encyst in the cockle foot, they can induce damages such as muscle fibre destruction (Jensen et al., Reference Jensen, Castro and Bachelet1999) through mechanical pressure and tissue lysis related to the secretion of enzymes by the cercariae (Lauckner, Reference Lauckner and Kinne1983). In addition, when the abundance of Himasthla spp. metacercariae exceeds a certain threshold (the value could depend on environmental conditions), cockle survival is reduced, as exemplified for the 4 species: (1) H. elongata induces mechanical obstruction in the cockle foot, increasing their burrowing time and making them more vulnerable to predators (Lauckner, Reference Lauckner and Kinne1983). Infection also induces a strong cockle immune response (Paul-Pont et al., Reference Paul-Pont, Gonzalez, Baudrimont, Jude, Raymond, Bourrasseau, Le Goïc, Haynes, Legeay, Paillard and de Montaudouin2010). Moreover, it modulates cockle biochemical performance and physiology by reducing their oxygen consumption, increasing antioxidant enzyme activity and modifying their energy allocation (Magalhães et al., Reference Magalhães, de Montaudouin, Figueira and Freitas2018b, Reference Magalhães, de Montaudouin, Figueira and Freitas2018c, Reference Magalhães, Freitas and de Montaudouin2020b). Finally, infection can significantly reduce (around 40%) cockle survival compared to non-infected cockles after 30 h under hypoxic conditions (Wegeberg and Jensen, Reference Wegeberg and Jensen1999). (2) Himasthla quissetensis promotes cockle emergence at the sediment surface, exposing them to other threats, like predation (Desclaux et al., Reference Desclaux, de Montaudouin and Bachelet2002) and can contribute to up to 46% of cockle population mortality (Desclaux et al., Reference Desclaux, de Montaudouin and Bachelet2004). (3) Himasthla interrupta moderately significantly reduces the cockle growth rate (de Montaudouin et al., Reference de Montaudouin, Binias and Lassalle2012b), and a marginal but significant loss of infected cockle flesh weight and body condition was observed by Wegeberg and Jensen (Reference Wegeberg and Jensen2003). (4) In contrast, no effect on cockles was reported concerning H. continua, with the exception of cockle burrowing time increasing at the sediment surface (Jensen et al., Reference Jensen, Castro and Bachelet1999).
Long-term monitoring
The dataset included cockles from 2 to 38 mm, corresponding to 0+ to 3+ year old cockles. Globally, the parasite community was dominated by G. minutus (mean of 62.8% of the total number of metacercariae per cockle), H. interrupta and H. quissetensis (16.2 and 5.5%, respectively). The other species were Curtuteria arguinae, Diphterostomum brusinae, H. continua, H. elongata, Psilostomum brevicolle and R. roscovitus. The following results aimed to obtain a mean Himasthla-host phenology calculated from our 20-year monthly monitoring.
Himasthla interrupta
Infection by H. interrupta started with 2 mm cockles, and prevalence regularly increased with shell length (ρ = 0.88, P < 0.001), to attain a median asymptotic prevalence of 80% (Fig. 1).
Prevalence of Himasthla interrupta (black line) and Himasthla quissetensis (grey line) by shell length class and number of dissected cockles (bars).
Mean intensity of infection increased significantly between recruitment in May [4 metacercariae per cockle, standard deviation (s.d.) = 7)] and December (50 metacercariae per cockle, s.d. = 70) (Wilcoxon test, P < 0.001) (Fig. 2A). Then, a significant decrease was observed until February (28 metacercariae per cockle, s.d. = 44) (P < 0.001), before a sharp increase during the second summer, reaching 59 metacercariae per cockle in August (s.d. = 62) (P < 0.001). Another decrease was observed in September (22 metacercariae per cockle, s.d. = 24) (P < 0.001), with a stabilization of the parasite intensity until December, which was the limit for assigning a shell length to an age, and thus identifying a cohort.
Boxplot of H. interrupta (A) and H. quissetensis (B) intensity per cockle shell length and corresponding age and seasons. Absolute age was deduced from a recruitment date in May. The box (25–75% of the data) contains a black line (median) and a red line (mean). Whiskers represent the lower and upper values in the range of ±1.5 interquartile range, with outliers as black circles. Grey arrows indicate significant variation between successive months (Wilcoxon test, P < 0.01). For example, in the case of H. interrupta, the first value that is significantly different from May 0+ intensity is in December 0+.
When excluding the smallest cockles (below 15 mm shell length, corresponding to 3-month old) which are rarely infected regardless of environmental conditions, and excluding the largest cockles (over 25 mm), for which age is uncertain due to slow growth, there was a weak correlation between cockle metacercariae infection and cockle density (Spearman test, ρ = −0.20, P = 0.03, not shown).
The structure of the trematode community was similar between cockles with and without H. interrupta (χ 2 test, P = 0.647). Himasthla quissetensis, G. minutus and C. arguinae co-dominated, representing 81% of the total abundance (Fig. 3A and B). However, species richness in cockles without H. interrupta (2.0 species, s.d. = 1.2) was lower than in those with H. interrupta (3.4 species, s.d. = 1.1) (Wilcoxon test, P < 0.001). Similarly, the mean number of metacercariae in cockles without H. interrupta (28 metacercariae per cockle, s.d. = 35) was lower than those with H. interrupta (66 metacercariae per cockle, s.d. = 107) (Wilcoxon test, P < 0.001).
Percentage of metacercariae per species (Curtuteria arguinae, Gymnophallus minutus, Psilostomum brevicolle, Renicola roscovitus, Diphterostomum brusinae and H. quissetensis) in Cerastoderma edule without (A) and with (B) H. interrupta.
Himasthla quissetensis
Infection by H. quissetensis started for 2 mm cockles, and prevalence regularly increased with shell length (ρ = 0.78, P < 0.001), to reach a median asymptotic prevalence of 80% (Fig. 1). Mean intensity of infection increased significantly between recruitment in May (1 metacercariae per cockle, s.d. = 1) and February (11 metacercariae per cockle, s.d. = 14) (Wilcoxon test, P < 0.001) (Fig. 2B). Then, a significant decrease was observed until July (7 metacercariae per cockle, s.d. = 9) (P < 0.001), followed by a stagnation in summer before another decrease in October of the second year (P < 0.001). A late infection was observed between October and November (6 metacercariae per cockle, s.d. = 12) (P = 0.005). Beyond December of the second year, the correspondence between cockle shell length and cockle age determination was no longer reliable due to slower growth.
For cockles with a shell length ranging between 15 and 25 mm (see section on H. interrupta) there was a weak correlation between cockle metacercariae infection and cockle density (Spearman test, ρ = −0.29, P < 0.001). The structure of the trematode community was similar between cockles with and cockles without H. quissetensis (χ 2 test, P = 0.183). Himasthla interrupta, G. minutus and C. arguinae co-dominated (Fig. 4A and B). However, species richness in cockles without H. quissetensis (1.9 species, s.d. = 1.4) was lower than those with H. quissetensis (2.9 species, s.d. = 1.3) (Wilcoxon test, P < 0.001). Similarly, the mean number of metacercariae in cockles without H. quissetensis (37 metacercariae per cockle, s.d. = 78) was lower than those with H. quissetensis (62 metacercariae per cockle, s.d. = 101) (Wilcoxon test, P < 0.001).
Percentage of metacercariae per species (C. arguinae, G. minutus, P. brevicolle, R. roscovitus, D. brusinae and H. interrupta) in C. edule without (A) and with (B) H. quissetensis.
Molecular identification
All metacercariae from cockles were first identified under a stereomicroscope based on morphological characteristics. Then, these metacercariae were compared with molecular tools. Himasthla continua metacercariae are morphologically very similar to H. elongata (the latter are slightly larger, and they both have 29 oral spines) and both occupy the same niche (i.e. the cockle foot). This difficulty of identification explains why they were not considered in the monitoring. Additionally, no specific type sequence was determined for H. continua due to a large variability of gene sequences within the specimens considered as H. continua. Conversely, sequences of 18S, ITS and COI were obtained for H. elongata, and sequences of 18S and COI were confirmed from cercariae emitted by L. littorea. Himasthla interrupta metacercariae have very small and light metacercariae settled in the mantle margin, at the opposite side to the siphon. For this species, only sequences of 18S and COI were obtained, with 100% similarity between samples and 100% agreement with morphological identification. Himasthla quissetensis is the only species with 31 oral spines. Sequences of 18S, ITS and COI were obtained with 100% similarity between samples and 100% agreement with morphological identification. The amplified products of 18S, ITS1 and COI presented 549, 798 and 273 bp, respectively, for H. elongata and 516, 529 and 284 bp for H. quissetensis. The amplified products of 18S and COI for H. interrupta presented, respectively, 541 and 259 bp, and 535 and 281 bp for cercariae collected from L. littorea.
Discussion
Size- and density-dependent infection and seasonality
For both Himasthla species, the long-term data analysis showed that infestation started rapidly after recruitment for cockles with a 2 mm shell length. This early infection in cockles was experimentally observed for all Himasthla species (Jensen et al., Reference Jensen, Castro and Bachelet1999; Wegeberg et al., Reference Wegeberg, de Montaudouin and Jensen1999; de Montaudouin et al., Reference de Montaudouin, Jensen, Desclaux, Wegeberg and Sajus2005) and is consistent with previous field studies (de Montaudouin et al., Reference de Montaudouin, Kisielewski, Bachelet and Desclaux2000; Desclaux et al., Reference Desclaux, de Montaudouin and Bachelet2004). The positive relationship between parasite prevalence and cockle shell length was also documented and ascribed to the higher filtration rate and longer exposure time of older/larger individuals (André et al., Reference André, Jonsson and Lindegarth1993; Riisgård, Reference Riisgård2001), resulting in a higher exposure to infective stages and thus greater parasite accumulation (de Montaudouin et al., Reference de Montaudouin, Wegeberg, Jensen and Sauriau1998; Mouritsen et al., Reference Mouritsen, McKechnie, Meenken, Toynbee and Poulin2003; Thieltges and Reise, Reference Thieltges and Reise2006a).
A moderate negative correlation between cockle density and intensity of H. interrupta and H. quissetensis was highlighted in this study, suggesting a dilution effect. Indeed, dense cockle populations can filter a high volume of water and thus eliminate parasitic cercariae, with subsequent lower metacercariae infection in cockles (Mouritsen et al., Reference Mouritsen, McKechnie, Meenken, Toynbee and Poulin2003; Thieltges and Reise, Reference Thieltges and Reise2006b; Buck and Lutterschmidt, Reference Buck and Lutterschmidt2017; Magalhães et al., Reference Magalhães, Freitas, Dairain and De Montaudouin2017; Correia et al., Reference Correia, Magalhães, Freitas, Bazairi, Gam and de Montaudouin2020a). However, the density of cockles only explained 4–8% of metacercariae intensity, implying that other factors, such as host condition, first intermediate host density and environmental parameters modulated the infection of the cockles (Wilson et al., Reference Wilson, Bjørnstad, Dobson, Merler, Poglayen, Read, Skorping, Hudson, Rizzoli, Grenfell, Heesterbeek and Dobson2001; Mouritsen et al., Reference Mouritsen, McKechnie, Meenken, Toynbee and Poulin2003; Thieltges and Reise, Reference Thieltges and Reise2006b; Welsh et al., Reference Welsh, Hempel, Markovic, van der Meer and Thieltges2019).
The phenology's pattern of infection was similar during the first year in both H. interrupta and H. quissetensis. The infection occurred during the warmer season, as already reported for H. quissetensis (Prévot, Reference Prévot1974; Desclaux et al., Reference Desclaux, de Montaudouin and Bachelet2004), but also for closely related species such as Himasthla littorinae (Nikolaev et al., Reference Nikolaev, Levakin and Galaktionov2020), H. elongata (Nikolaev et al., Reference Nikolaev, Levakin and Galaktionov2021), C. arguinae (Desclaux et al., Reference Desclaux, Russell-Pinto, de Montaudouin and Bachelet2006) and other trematode families such as gymnophallids (Gam et al., Reference Gam, de Montaudouin and Bazairi2009b), renicolids (Thieltges and Rick, Reference Thieltges and Rick2006) or microphallids (Meißner, Reference Meißner2001). Temperature of water is an important trigger to stimulate infection by cercariae (Lo and Lee, Reference Lo and Lee1996; Mouritsen and Jensen, Reference Mouritsen and Jensen1997; Mouritsen, Reference Mouritsen2002; Koprivnikar et al., Reference Koprivnikar, Ellis, Shim and Forbes2014). In particular, in situ experiments showed that the thermal window for cockle infection by H. quissetensis was between 15 and 23°C, with maximum infection being at 19–20°C (de Montaudouin et al., Reference de Montaudouin, Blanchet, Desclaux-Marchand, Lavesque and Bachelet2016). During the first winter, the infection intensity per cockle decreased significantly, by 44 and 38% for H. interrupta and H. quissetensis, respectively. The decrease in the mean parasite intensity could be due to immigration of the less heavily infected cockles or emigration of highly parasitized cockles. This hypothesis seems irrelevant regarding the low locomotive capacity of adult cockles (Richardson et al., Reference Richardson, Ibarrola and Ingham1993). It could also be explained by the death of metacercariae in cockles. There are very few studies exploring the dynamics of parasite infrapopulations (i.e. populations at the scale of a host individual). Mortality of parasites was observed in cockles for the non-encysted metacercariae of G. minutus (de Montaudouin et al., Reference de Montaudouin, Binias and Lassalle2012b). In this case, the authors had transplanted cockles, and the new site could have been deleterious to parasites, but in other cases G. minutus can suffer from hyperparasitism (Fermer et al., Reference Fermer, Culloty, Kelly and O'Riordan2010). However, in the case of Himasthla spp. and their encysted metacercariae, empty cysts that suggest parasite death have been observed and registered at a very low intensity (Desclaux et al., Reference Desclaux, de Montaudouin and Bachelet2004), leading to the exclusion of this conjecture as well. Finally, the death of the most heavily infected cockles could explain the reduction of Himasthla metacercariae intensity in winter. This third hypothesis is the most likely, and has been mentioned in several studies concerning trematodes in their second intermediate hosts (Kennedy, Reference Kennedy1984; Desclaux et al., Reference Desclaux, de Montaudouin and Bachelet2004, Reference Desclaux, Russell-Pinto, de Montaudouin and Bachelet2006; Gam et al., Reference Gam, de Montaudouin and Bazairi2009b) or first intermediate hosts (Bowers, Reference Bowers1969; Schmidt and Fried, Reference Schmidt and Fried1997; Rantanen et al., Reference Rantanen, Valtonen and Holopainen1998; Watters, Reference Watters1998). During the second summer, the infection pattern was less obvious. It is noteworthy that during summer infections, a stable parasite intensity in cockles can result from a balance between parasite infection and parasite-dependent mortality processes.
Parasite co-occurrence
The trematode species richness presented in this study was similar to what has been reported in similar ecosystems along the northeast Atlantic coast (Krakau et al., Reference Krakau, Thieltges and Reise2006; Thieltges et al., Reference Thieltges, Krakau, Andresen, Fottner and Reise2006; Gam et al., Reference Gam, de Montaudouin and Bazairi2009b; Magalhães et al., Reference Magalhães, Correia, de Montaudouin and Freitas2018a; Correia et al., Reference Correia, Picado, de Montaudouin, Freitas, Rocha, Dias and Magalhães2020b). Negative interactions among parasites within their host have been poorly documented, and in particular few studies have investigated the effect of invasive sporocyst stages on the global diversity of trematodes. Neither M. parvus nor G. choledochus sporocysts influence the prevalence or abundance of other trematode species (Magalhães et al., Reference Magalhães, Daffe, Freitas and de Montaudouin2020a), contrary to those observed concerning B. minimus whose presence is linked to a higher abundance of other trematode species (Magalhães et al., Reference Magalhães, Freitas and de Montaudouin2015). Magalhães et al. (Reference Magalhães, Freitas and de Montaudouin2015) suggested that B. minimus infection could impair cockle resistance to metacercariae infection, or that high metacercariae infection could facilitate B. minimus infestation. However, a second hypothesis is that all parasites co-occur independently of one another, and infect cockles because all conditions are favourable to infection by all parasite species (environmental factors, cockle fitness, other host presence). In the present study, the fact that H. interrupta (or H. quissetensis) occurrence is associated with higher trematode species richness and abundance, with similar community structure, favours the second hypothesis. Indeed, the relatively low metacercariae intensity values observed, combined with the occupation of specific organs by most trematode species (de Montaudouin et al., Reference de Montaudouin, Thieltges, Gam, Krakau, Pina, Bazairi, Dabouineau, Russell-Pinto and Jensen2009), are weak arguments supporting an interspecific metacercarial competition, as observed by Thieltges and Reise (Reference Thieltges and Reise2006b) and Lassalle et al. (Reference Lassalle, de Montaudouin, Soudant and Paillard2007). In addition, an interspecific competition between Himasthla species was not expected, as their metacercariae do not grow inside their second intermediate host (de Montaudouin et al., Reference de Montaudouin, Jensen, Desclaux, Wegeberg and Sajus2005).
Molecular identity
For H. elongata, H. quissetensis and H. interrupta, the metacercariae molecular identification was performed using 18S and COI sequences. Concerning H. elongata, all analysed sequences matched each other, and also sequences that were isolated from L. littorea cercariae. These results validate the molecular identification of H. elongata since L. littorea is the first intermediate host of only this Himasthla species. For H. quissetensis and H. interrupta, the sequences matched each other, and thus provided a good molecular identification. All sequenced H. quissetensis came from samples extracted from the cockles' foot, while samples corresponding to H. interrupta were extracted from the mantle. However, a mismatch occurred for H. continua, with no match among the analysed sequences (high variability). We cannot rule out that this high variability of the obtained sequences was associated with the presence of some larvae of another Himasthla species, e.g. Himasthla leptosoma. The larvae of these two species are hardly distinguishable by microscopic methods, being very similar in the size of their cysts as well as in the number of spines on the collar (Galaktionov et al., Reference Galaktionov, Solovyeva and Miroliubov2021). Finally, our results confirm the identity of three species of Himasthla metacercariae, which can be difficult to distinguish under a stereomicroscope based on morphological identification.
Conclusion
Trematodes of the Himasthla genus are very common parasites of cockles. Their effect on the cockle individuals or populations is usually reported as low. From an evolution point of view, the metacercarial stage is an opportunity to accumulate diverse parasite genotypes in order to contribute to the genetic diversity of trematode populations (Leung et al., Reference Leung, Poulin and Keeney2009). Thus, the main objective of the parasite in this parasitic stage would not be to consume the host energy, as occurs in the first and final hosts. However, a literature review and analysis of a 20-year database revealed that some Himasthla-dependent negative effects occur when the metacercariae infection reaches high levels. Considering that this infection is often related to temperature, this parasite dynamics should be monitored according to different climate change scenarios. While morphological identification is particularly difficult concerning Himasthla genus, new molecular sequences provided in this study may be helpful for an accurate identification of some species, although uncertainties still remain concerning H. continua.
Data
The data presented in this study are available on request from the corresponding author.
Acknowledgements
The authors acknowledge the SEPANSO (Société d'Etude pour la Protection et l'Aménagement de la Nature dans le Sud-Ouest) which manages the National Reserve of Arguin. The authors also acknowledge MDPI for correcting the English of the manuscript. The authors are grateful to the referee for his help in improving the paper.
Author contributions
X. M. conducted data gathering. A. R. produced the literature review. G. D. realized molecular analyses. A. R., X. M. and L. M. performed statistical analyses. A. R. and X. M. wrote the original draft paper. A. R., G. D., L. M., O. M. and X. M. reviewed and corrected the article.
Financial support
This research was partly funded by the INTERREG-ATLANTIC programme through the research project COCKLES (EAPA_458/2016 COCKLES: Cooperation for restoring cockle shellfisheries and their ecosystem services in the Atlantic Area). This work is part of Anaïs Richard's doctoral thesis (University of Bordeaux – 2018-SG-D-13) financed by a doctoral grant of the French ‘Ministère de l'Enseignement Supérieur et de la Recherche’. Sampling was performed thanks to Planula 4 vessels (CNRS-INSU, Flotte Océanographique Française). Luísa Magalhães acknowledges the financial support of CESAM by FCT/MCTES (UIDP/50017/2020 + UIDB/50017/2020 + LA/P/0094/2020), through national funds.
Conflict of interest
The authors declare there are no conflicts of interest.
Ethical Standards
Not applicable.
Open access
