Introduction
Animal hosts are exposed to complex pathogen assemblages that ultimately form a dynamic community within them (Johnson and Buller, Reference Johnson and Buller2011). However, research into host–pathogen interactions remains dominated by the study of one host–one pathogen systems despite the fact that hosts can simultaneously carry several infectious agents, with consequences for the dynamics of each agent and host health. Coinfections have been previously reported in shellfish hosts: including Vibrio tapetis, Perkinsus sp. and digenetic trematodes in cockles (Lassalle et al., Reference Lassalle, de Montaudouin, Soudant and Paillard2007); ostreid herpesvirus-1 (OsHV-1) and Vibrio spp. in oysters (Petton et al., Reference Petton, Bruto, James, Labreuche, Alunno-Bruscia and Le Roux2015; Alfaro et al., Reference Alfaro, Nguyen and Merien2019); and protozoans and bacteria in clams (Arzul et al., Reference Arzul, Chollet, Michel, Robert, Garcia, Joly, Francois and Miossec2012; Carella et al., Reference Carella, Elisabetta, Simone, Fulvio, Daniela, Prado, Rossella, Marino, Eleonora, Tobia and De Vico2020); and in some of the cases with a dramatic effect on disease susceptibility (Cox, Reference Cox2001; Lassalle et al., Reference Lassalle, de Montaudouin, Soudant and Paillard2007; Arzul et al., Reference Arzul, Chollet, Michel, Robert, Garcia, Joly, Francois and Miossec2012; Alfaro et al., Reference Alfaro, Nguyen and Merien2019). Processes at different scales of ecological organization from the within-host level (e.g. interactions of coinfecting parasites) to the ecosystem level (e.g. the influence of environmental variables) across space and over time have a key influence in complex natural systems (Hellard et al., Reference Hellard, Fouchet, Vavre and Pontier2015). Therefore, moving from one host–one pathogen systems towards an ecosystem view of host–pathogen interactions and their ecology is essential.
Detection of interspecific parasite/pathogen interactions in natural populations is not easy, due to complex networks of indirect effects making it difficult to infer underlying processes (Hellard et al., Reference Hellard, Pontier, Sauvage, Poulet and Fouchet2012). Interactions among coinfecting infectious agents can, in fact, alter host pathology, parasite transmission and virulence evolution, also influencing the spread of infections at a population level resulting in disease outbreaks, as described previously in bivalves (Lassalle et al., Reference Lassalle, de Montaudouin, Soudant and Paillard2007; Arzul et al., Reference Arzul, Chollet, Michel, Robert, Garcia, Joly, Francois and Miossec2012; Alfaro et al., Reference Alfaro, Nguyen and Merien2019). These effects can be the opposite, i.e. coinfections within the host may reduce infection success yet still enhance pathology or, in contrast, co-occurring parasites/pathogens may weaken host pathology increasing the infection success (Johnson and Hoverman, Reference Johnson and Hoverman2012). For example, mortality induced by one parasite can also eliminate the availability of hosts for other parasites (Jolles et al., Reference Jolles, Ezenwa, Etienne, Turner and Olff2008); likewise, morbidity induced by one parasite can increase exposure to a second, even if within-host interactions are antagonistic (Karvonen et al., Reference Karvonen, Seppala and Valtonen2009).
Simultaneous infections may therefore trigger a whole spectrum of outcomes within the hosts, both synergistic and antagonistic, i.e. the weight of one or both infectious agents may be suppressed (Fukami et al., Reference Fukami, Dickie, Wilkie, Paulus, Park, Roberts, Buchanan and Allen2010), one or both agents may be amplified resulting in coexistence (Thomas et al., Reference Thomas, Renaud and Poulin1998), or one may be amplified and the other suppressed (Johnson and Hoverman, Reference Johnson and Hoverman2012). Parasite/pathogen interactions can be not only direct, such as competition for attachment sites, competition for host resources and predation upon one another (intra-guild predation) (Hatcher et al., Reference Hatcher, Dick and Dunn2006, Reference Hatcher, Dick and Dunn2012; Johnson et al., Reference Johnson, Dobson, Lafferty, Marcogliese, Memmott, Orlofske, Poulin and Thieltges2010), but interactions can also be indirect or host-mediated often involving changes in immunity such as cross-immunity (apparent competition) and immune suppression (apparent facilitation) (Jolles et al., Reference Jolles, Ezenwa, Etienne, Turner and Olff2008). For instance, Johnson and Hoverman (Reference Johnson and Hoverman2012) suggested cross-reactive immunity as a cause of a decrease in parasite persistence in their experiments with a group of related larval trematodes in Pacific chorus frogs (Pseudacris regilla). Johansen and Sommer (Reference Johansen and Sommer2001) suggested that the immune suppression triggered by a primary infection could be the cause of multiple bacterial infections in Atlantic salmon (Salmo salar).
The co-occurrence of infectious agents may also result simply because the same risk factors promote their presence and not because the different pathogen groups are interacting synergistically (Vaumourin et al., Reference Vaumourin, Vourc'h, Gasqui and Vayssier-Taussat2015). Risk factors shared by two or more infectious agents can occur at the level of the host individual (e.g. host condition and age), population (e.g. host density) or landscape (e.g. climatic factors), increasing the probability of infection by both/all agents (Hellard et al., Reference Hellard, Fouchet, Vavre and Pontier2015). Such confounding factors may influence host exposure and/or host susceptibility promoting simultaneous infections (Vaumourin et al., Reference Vaumourin, Vourc'h, Gasqui and Vayssier-Taussat2015). Environmental parameters such as seawater temperature and salinity, host density or host physiological conditions and behaviour are common confounding factors that may influence the co-occurrence of different pathogen groups and the epidemiological and geographic patterns of infection and disease (Hellard et al., Reference Hellard, Pontier, Sauvage, Poulet and Fouchet2012; Vaumourin et al., Reference Vaumourin, Vourc'h, Gasqui and Vayssier-Taussat2015). Confounding factors must be considered in field studies since those factors can create statistical associations between infectious agents even if there is no true biological interaction between them (Kuris and Lafferty, Reference Kuris and Lafferty1994).
Generalist host species, such as the common cockle Cerastoderma edule, with wide spatial distributions, are exposed to a large range of environmental conditions, and a greater diversity of parasites and pathogens, and consequently are more likely to be coinfected (Vaumourin et al., Reference Vaumourin, Vourc'h, Gasqui and Vayssier-Taussat2015; Mahony et al., Reference Mahony, Lynch, Egerton, Cabral, de Montaudouin, Fitch, Magalhaes, Rocroy and Culloty2020; 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), directly affecting individual cockle health and cockle population dynamics (Longshaw and Malham, Reference Longshaw and Malham2013). Cockles may be host to multiple macro-parasite infections, particularly trematodes (de Montaudouin et al., Reference De Montaudouin, Paul-Pont, Lambert, Gonzalez, Raymond, Jude, Legeay, Baudrimont, Dang, Le Grand, Le Goic, Bourasseau and Paillard2010, Reference De Montaudouin, Arzul, Cao, Carballal, Bruno, Correia, Cuesta, Culloty, Daffe, Darriba, Díaz, Freitas, Garcia, Goedknegt, Grade, Groves, Iglesias, Jensen and Villalba2021), but may also be infected by a range of pathogens (Rowley et al., Reference Rowley, Cross, Culloty, Lynch, Mackenzie, Morgan, O'Riordan, Robins, Smith, Thrupp, Vogan, Wootton and Malham2014; 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), providing a microcosm for the parasite/pathogen community. Within this community, there are haplosporidian protists such as the genus Minchinia (Engelsma et al., Reference Engelsma, Roozenburg, Voorbergen-Laarman, Van den Brink, Troost, Ysebaert, Bateman and Longshaw2011; Longshaw and Malham, Reference Longshaw and Malham2013; Ramilo et al., Reference Ramilo, Abollo, Villalba and Carballal2018; Lynch et al., Reference Lynch, Lepee-Rivero, Kelly, Quinn, Coghlan, Bookelaar, Morgan, Finarelli, Carlsson and Culloty2020b; Albuixech-Martí et al., Reference Albuixech-Martí, Lynch and Culloty2020; 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), numerous members of the bacterial family Vibrio (Carballal et al., Reference Carballal, Iglesias, Santamarina, Ferro-Soto and Villalba2001; Lassalle et al., Reference Lassalle, de Montaudouin, Soudant and Paillard2007; 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), viruses such as ostreid herpesvirus-1 microVar (OsHV-1 μVar) (Carballal et al., Reference Carballal, Villalba, Iglesias and Hine2003; Evans et al., Reference Evans, Paul-Pont and Whittington2017; Bookelaar et al., Reference Bookelaar, Lynch and Culloty2020; 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) and microsporidian species reported in digeneans and paramyxids infecting cockle populations (Goater, Reference Goater1993; Fermer et al., Reference Fermer, Culloty, Kelly and O'Riordan2009, Reference Fermer, Culloty, Kelly and O'riordan2011; Villalba et al., Reference Villalba, Iglesias, Ramilo, Darriba, Parada, No, Abollo, Molares and Carballal2014; 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). Cockles are a keystone species responsible for different ecosystem services, i.e. carbon storage, energy cycling, food source for seabirds, etc., being also good sentinel and bioindicator species (Malham et al., Reference Malham, Hutchinson and Longshaw2012; Carss et al., Reference Carss, Brito, Chainho, Ciutat, de Montaudouin, Fernandez Otero, Incera Filgueira, Garbutt, Goedknegt, Lynch, Mahony, Maire, Malham, Orvain, Olivier and Jones2020). Moreover, the common cockle is one of the main non-cultured bivalve species harvested in western European waters, with a high economic value (Rowley et al., Reference Rowley, Cross, Culloty, Lynch, Mackenzie, Morgan, O'Riordan, Robins, Smith, Thrupp, Vogan, Wootton and Malham2014; Carss et al., Reference Carss, Brito, Chainho, Ciutat, de Montaudouin, Fernandez Otero, Incera Filgueira, Garbutt, Goedknegt, Lynch, Mahony, Maire, Malham, Orvain, Olivier and Jones2020). Consequently, cockles have been well-studied and are useful model species. Given the complex parasite/pathogen assemblages within the cockles and their intricate interactions and effects, a comprehensive and integrated approach that considers not only the synergistic impact of coinfections but also the role of ecological factors driving the infection across space and over time is needed.
Future climate scenarios predict a warming marine environment and increased precipitation events resulting in pulse events of reduced/fluctuating salinity in nearshore ecosystems (Allam and Raftos, Reference Allam and Raftos2015; Coen and Bishop, Reference Coen and Bishop2015). Consequently, climate change may be expected to have a significant effect on disease epidemiology (Coen and Bishop, Reference Coen and Bishop2015), and may even have significant consequences on zoonoses emergence (Hoarau et al., Reference Hoarau, Mavingui and Lebarbenchon2020). In this context, the need for the understanding of ecological processes involved in the transmission and dynamics of infectious agents in host reservoirs is critical (Hoarau et al., Reference Hoarau, Mavingui and Lebarbenchon2020). This study, therefore, aims to assess the interactions between infectious agents, their hosts and the environment.
For that purpose, this study examined the diversity and prevalence of significant pathogen groups, Haplosporidia, Vibrio spp., OsHV-1 μVar and Microsporidia, associated with bivalve mortalities globally as well as coinfections and associations of the different infectious agents at an individual and population level in wild and fished C. edule along the Irish and Celtic Seas. The role of abiotic (seawater temperature and salinity) and biotic (cockle size and age, and epiflora on their shells) factors in the risk of coinfection throughout the study was also assessed. Developing knowledge on pathogen diversity and co-occurrence in an economically and ecologically important species as is C. edule and integrating the role of the environment in the host–pathogen interplay is essential for a better understanding of the interaction networks that can affect disease outcomes and transmissibility in the current changing environment. Findings from this study highlight the complexity of pathogen interactions within a host and the effect of key biotic and environmental drivers shaping disease dynamics, which can also be applied to other coinfected animal host species.
Materials and methods
Sampling
Cockles were sampled from spring 2018 to spring 2019 at Cork Harbour (Ringaskiddy and Cuskinny, n = 257) on the south coast (Celtic Sea) of Ireland and from summer 2018 to spring 2019 at Youghal Bay (n = 69) and Dungarvan Harbour (n = 169), with an extensive farming of Pacific oysters (Crassostrea gigas), on the south coast (Celtic Sea) of Ireland, where there is no commercial fishing activity (Fig. 1). Additionally, 240 cockles were collected from summer 2018 to spring 2019 at the commercial fishery at Dundalk Bay (Annagassan and Cooley) on the northeast coast (Irish Sea) of Ireland (Fig. 1). Dundalk Bay is a classified Bivalve Mollusc Production Area and it has supported a commercial dredge cockle (C. edule) fishery since 2001 (Hervas et al., Reference Hervas, Tully, Hickey, O'Keefe and Kelly2008).
Map of Ireland highlighting the cockle sample sites with coordinates.
The four locations are bays with intertidal sand and mudflats, where a variety of bivalves inhabit. Cerastoderma edule densities in the selected sites vary from higher densities, 9.33 ± 3.5 ind m−2, in Dundalk Bay (Shellfish Stocks and Fisheries Review, Marine Institute and BIM, 2017) to lower densities in the no commercial fishing sites: 1.6 ± 2.1 ind m−2 in Youghal Bay; 0.8 ± 1.7 ind m−2 in Cork Harbour; and 0.4 ± 1.3 ind m−2 in Dungarvan Harbour (Fermer et al., Reference Fermer, Culloty, Kelly and O'riordan2011). Cockle mortalities associated with disseminated neoplasia have been previously reported in Cork Harbour (Morgan et al., Reference Morgan, O'Riordan, Kelly and Culloty2012). While a massive cockle mortality event associated with dinoflagellates occurred in Youghal Bay in the past (Ottway et al., Reference Ottway, Parker, McGrath and Crowley1979). Dungarvan Harbour has also a history of OsHV-1-associated mortality events (Lynch et al., Reference Lynch, Carlsson, Reilly, Cotter and Culloty2012; Prado-Alvarez et al., Reference Prado-Alvarez, Darmody, Hutton, O'Reilly, Lynch and Culloty2016). However, to the best of our knowledge, Dundalk Bay has not previously recorded mortality events.
Even though Cork Harbour is an important industrial area and is a busy ferry/shipping terminal, its trophic status (based on the nutrient levels such as phosphorus and nitrogen, growth of algae and dissolved oxygen concentration) was qualified as unpolluted by the EPA in 2017 along with Dungarvan Harbour, an extremely sheltered bay, while Youghal Bay, close to an urban hub, was classified as intermediate polluted area. In turn, Dundalk Bay showed a localized effect on the trophic status, with Annagassan classified as an unpolluted area and Cooley as an intermediate polluted area, probably due to the influence of water discharge from the local river.
Cockles were collected by hand raking from the intertidal area of each site at spring tide (0.5–1 m) and seasonal intervals. The sample size, outlined in Table 1, was dependant on the availability of cockles at each site and season. Cockles were kept at ambient temperature during the fieldwork and were processed either on the same day or were held out of water overnight at 4 °C and were processed the next day.
Sample sites, seasons and number of cockles collected
(–) No samples collected.
a In brackets the % representation of all cockles sampled.
Sample processing and diagnostic methods
Morphometrics
Measurements included (A) cockle body size [height (mm)] of each individual that was recorded using Vernier callipers. The (B) cockle's age was established by counting the annual winter growth rings. For verification and validation purposes, the number of growth rings was confirmed by a second researcher, as it is recognized that variations in growth during the same year lead to the possible formation of extra stria (de Montaudouin, Reference De Montaudouin1996). Ultimately, the presence or absence of epiflora on the shell of each individual was also recorded. Epiflora was mainly composed of epiphytic algae, although no specific species were identified.
Histology
A transverse section of each cockle was excised, including the digestive gland, gonads, gills and mantle tissues; as well as a section of the hinge ligament and a section of the tip of the foot were taken for histological examination. The three sections were placed in a labelled histocassette and were preserved in Davidson's fixative for 48 h (Shaw and Battle, Reference Shaw and Battle1957) prior to processing. Paraffin-embedded tissue blocks were sectioned at 5 μ m and stained using haematoxylin and eosin (Sigma Aldrich, USA) (Howard and Smith, Reference Howard and Smith1983). Tissue slides were screened at a magnification of 40x and under oil at 100x for visualization of the pathogens and tissue pathology.
Molecular techniques for cockle species identification and pathological screening
A small piece of gill tissue (2–5 mm2 of tissue) from each cockle was taken and stored at −20 °C for molecular assays, to determine the European cockle species [C. edule (Linnaeus, 1758) or Cerastoderma glaucum (Poiret, 1789)] and for subsequent pathological screening. Genomic DNA from the cockle gill tissue was extracted using the chelex-100 extraction method (Walsh et al., Reference Walsh, Metzger and Higuchi1991; Lynch et al., Reference Lynch, Mulcahy and Culloty2008).
Polymerase chain reaction (PCR) was carried out to amplify the nuclear DNA markers ITS-for/ITS Ce-R/ITS Cg-R to differentiate between the presence of C. edule, C. glaucum species or hybrids (Freire et al., Reference Freire, Arias, Mendez and Insua2011; Table 2). Amplification was conducted following the reaction mixture and thermocycling conditions, as well as the visualization of the product, described in Albuixech-Martí et al. (Reference Albuixech-Martí, Lynch and Culloty2020).
Description of PCR primer pairs showing sequences for each forward and reverse primer and expected product size
Standard PCR screening for haplosporidian detection was carried out in cockle gill tissue samples (n = 735) using generic haplosporidian HAP-F1 and HAP-R3 primers that amplify small regions of the SSU rDNA of most haplosporidians (Renault et al., Reference Renault, Stokes, Chollet, Cochennec, Berthe, Gerard and Burreson2000; Molloy et al., Reference Molloy, Giamberini, Stokes, Burreson and Ovcharenko2012; Table 2). Amplification was conducted following the reaction mixture and thermocycling conditions, as well as the visualization of the product, described in Albuixech-Martí et al. (Reference Albuixech-Martí, Lynch and Culloty2020).
Standard PCR screening for Minchinia spp. detection was carried out in the cockle samples that were positive for Haplosporidia by conventional PCR, using specific primers for Minchinia tapetis (TAP-For/Rev) and Minchinia mercenariae-like (MER-For/Rev) (Albuixech-Martí et al., Reference Albuixech-Martí, Lynch and Culloty2020; Table 2) in the appropriate pairings and separate PCR reactions. Amplification was conducted following the reaction mixture and thermocycling conditions, as well as the visualization of the product, described in Albuixech-Martí et al. (Reference Albuixech-Martí, Lynch and Culloty2020).
Standard PCR screening for Vibrio spp. in cockle gill tissues was conducted using generic primers Vib-F3/R3, developed in our lab for the detection of Vibrio genus [Kett et al., unpublished as per Notaro et al. (Reference Notaro, Culloty and Lynch2021), data confirmed to be specific for Vibrio spp. by Sanger sequencing; Table 2]. The screening was conducted in all the collected cockles (n = 735). A total of 2 μL of DNA per individual was used in a total volume of 27.5 μL of the reaction mixture containing: 5 μL of 5x green buffer, 5 μL of dNTPs (1.25 mm), 0.5 μL of MgCl2, 0.25 μL of forward and reverse primers, 0.1 μL of GoTaq polymerases, 1.5 μL of DMSO and 12.9 μL of ddH2O. Cycling conditions consisted of an initial denaturation of the sample at 95 °C for 1 min followed by 35 cycles of 94 °C for 20 s, 56 °C for 30 s, 72 °C for 30 s and a final elongation at 72 °C for 7 min. Electrophoresis of the amplification products was conducted in a 2% agarose gel.
Standard PCR screening for OsHV-1 and variants, including OsHV-1 μVar, was carried out in cockle gill tissue samples (n = 735) using specific primers OHVA/OHVB that amplify small products of the ORF4 gene in the OsHV-1 virus and its variants (Lynch et al., Reference Lynch, Dillane, Carlsson and Culloty2013; Table 2). Amplification was conducted following the reaction mixture and thermocycling conditions described in Lynch et al. (Reference Lynch, Dillane, Carlsson and Culloty2013). Electrophoresis of the amplification products was conducted in a 2% agarose gel.
Standard PCR using specific primers MicIF1/MicIR1 and MicIF2/MicIR2 (Lynch et al., unpublished data confirmed to be specific to Microsporidia spp. by Sanger sequencing; Table 2) to detect two unidentified microsporidian species was carried out in cockle gill tissue samples (n = 735). Both reaction mixture and thermocycling conditions were the same as used with the PCR for OsHV-1 and variants, as well as the visualization of the product.
Real-time quantitative PCR (qPCR) for detection and quantification of Vibrio aestuarianus in the cockle samples that were positive for Vibrio spp. was performed using specific primers (Table 3) for detection of the molecular chaperone dnaJ gene, following the protocol of McCleary and Henshilwood (Reference McCleary and Henshilwood2015). All samples were performed in triplicates using 5 μL of DNA. C t values were used to determine real-time PCR quantification and detection limits. A tested sample was considered positive if its mean C t value was below 37.
Description of qPCR primer pairs and probe
Viral load for OsHV-1 and variants was investigated by qPCR with a small number of individuals (n = 30, whose standard PCR results displayed some smear in the visualization of the product) to verify that there was no viral load by carrying out a more specific screening and increasing the amount of DNA processed. The protocol of Pepin et al. (Reference Pepin, Riou and Renault2008) was followed, using the specific primers HVDP-F/HVDP-R (Table 3). qPCR was conducted in duplicate using 5 μL DNA by a Thermo Hybaid PCR express thermal cycler.
Direct sequencing was carried out on representative PCR products (n = 10) amplified using generic Vibrio Vib-F3/R3 primers (Kett et al., unpublished as per Notaro et al., Reference Notaro, Culloty and Lynch2021) to identify the Vibrio species present in the samples. Likewise, direct sequencing was carried out on representative PCR products (n = 5) amplified using generic haplosporidian HAP-F1 and HAP-R3 primers (Renault et al., Reference Renault, Stokes, Chollet, Cochennec, Berthe, Gerard and Burreson2000; Molloy et al., Reference Molloy, Giamberini, Stokes, Burreson and Ovcharenko2012) to identify the species present in the samples beyond the Minchinia species screened. Genomic DNA from selected individuals was isolated and purified using the QIAquick Gel Extraction Kit (QIAgen) prior to direct sequencing. Both the forward and reverse strands of DNA samples were sequenced commercially (Source Bioscience). Each sequence was matched against the National Center for Biotechnology Information (NCBI) nucleotide database with BLASTn (Basic Local Alignment Search Tool), which finds regions of local similarity between sequences to identify and confirm the DNA being detected in the PCRs.
Environmental data
This study was conducted using E.U. Copernicus Marine Service Information (marine.copernicus.eu/services-portfolio/access-to-products). Two environmental parameters critical in shaping host–pathogen interactions in estuarine and marine environments (Coen and Bishop, Reference Coen and Bishop2015) were selected for this study as potential confounding factors: seawater temperature and salinity. Recorded data in seawater temperature (°C) and salinity (PSU), taken monthly at a depth of 0.5 m below sea level, were downloaded by ArcGIS Desktop 10.5.1, Redlands, CA, USA (Environmental Systems Research Institute, 2017) throughout the duration of the study (from April 2018 to May 2019) from the different sample sites (coordinate-based) along the Irish coast.
Statistical analysis
Statistical analysis was performed in R version 3.2.3. statistical software. The pathogen occurrence was modelled as pathogen presence/absence in each cockle for each pathogen group and their combinations, and it was used in the subsequent analysis to examine the associations between the different pathogen groups present through the samples.
Association screening approach (SCN), described in detail by Vaumourin et al. (Reference Vaumourin, Vourc'h, Telfer, Lambin, Salih, Seitzer, Morand, Charbonnel, Vayssier-Taussat and Gasqui2014), was conducted running 5000 simulations (NS) with the presence/absence of Haplosporidia and Vibrio spp. to determine a potential association between these two pathogen groups, under the null hypothesis of random pathogen associations [P(SCN) < 0.05]. Subsequently, a second SCN (NS = 5000) was performed with the presence/absence of the different pathogen species detected and quantified in the samples, i.e. M. tapetis, M. mercenariae-like and V. aestuarianus.
Corrected Pearson's χ 2 test, described in detail by Hellard et al. (Reference Hellard, Pontier, Sauvage, Poulet and Fouchet2012), was applied to the Haplosporidia and Vibrio spp. presence/absence data to determine if the coinfection was due to an underlying biological interaction or to confounding factors considered (seawater temperature, salinity, cockle size and age, and epiflora on their shells). Likewise, the test was applied to the M. tapetis and M. mercenariae-like presence/absence data. Biological interactions between two pathogens, therefore, will be suspected when the probability of coinfection is not random once confounding factors have been considered [P(Corr χ 2) < 0.05]. In the model, running 1000 bootstraps, the effects of the confounding factors are summarized as F1 + F2 + F3 + F4 + F5.
Moreover, Pearson's χ 2 tests were used to determine whether the prevalence observed of single infections and two-pathogen coinfections was significantly different [P(Chisq) < 0.05] among the sample sites and throughout the year. Fisher's exact test was conducted when the frequencies of the pathogen occurrence in a site or season were <4, to gain accuracy.
Results
In total, 735 cockles were screened to determine the prevalence of infection and coinfections of the target pathogens (Table 4). The cockle speciation PCR confirmed that all the individuals from all sample sites were C. edule.
Prevalence (%, positive individuals/screened individuals) of Haplosporidia spp., Minchinia tapetis, M. mercenariae-like, Vibrio spp. and Vibrio aestuarianus screening in Cerastoderma edule at each sample site by season
(–) No samples.
a The highest prevalence (%) for each pathogen group at each site/season is highlighted in bold.
b Potential Vibrio species were defined in Table 5.
Diversity and prevalence of single and multiple infections and the associations between them
No OsHV-1 μVar nor microsporidian species were detected in any of the cockles (n = 735) screened by PCR. The only pathogen groups detected during the study were Haplosporidia and Vibrio. Overall, Haplosporidia was the most common pathogen found through the samples, with a 37.7% (n = 277) occurrence (Fig. 2, detailed data in Table 4). Among these 277 infected cockles, two haplosporidian species were identified in the study: M. tapetis which was present in 29.6% (n = 82/277) of the tested individuals, and M. mercenariae-like which was present in 14.8% (n = 41/277) (Fig. 2, detailed data in Table 4). Nevertheless, 171 (61.7%) of the tested samples that were positive in the generic haplosporidian PCR did not amplify using the more specific PCR for Minchinia spp., which could indicate the presence of other haplosporidian species occurring in the samples. Thus, five of those samples amplified using generic HAP-F1/R3 primers were commercially sequenced; however, low-complexity sequences were obtained, and no significant similarity was found in GenBank. Therefore, the 277 positive PCR cases were defined as Haplosporidia spp. for statistical and analysis purposes.
Overall prevalence of infections and coinfections (%) with the different pathogen groups (Hap: Haplosporidia spp.; Vib: Vibrio spp.; M.tap: Minchinia tapetis; M.mer: Minchinia mercenariae-like; V.aes: Vibrio aestuarianus).
Vibrio was present in 25.3% (n = 186) of the total number of cockles screened (Fig. 2, detailed data in Table 4). Of the 186 samples that were positive for Vibrio, 19.9% (n = 37/186) were confirmed by qPCR to be V. aestuarianus (Fig. 2, detailed data in Table 4). Vibrio aestuarianus-positive controls (diluted purified DNA) have a C t = 30 (60 copies of V. aestuarianus DNA per μL). Based on that value, 13.5% (n = 5/37) of the individuals were classified as heavily infected individuals (C t ⩽ 30), while 86.5% (n = 32/37) were classified as low-intensity infections (37 > C t > 30).
Direct sequencing was carried out on representative PCR products (n = 10) amplified using generic Vibrio Vib-F3/R3 primers (Kett et al., unpublished as per Notaro et al., Reference Notaro, Culloty and Lynch2021) to identify other Vibrio species present in the samples apart from V. aestuarianus. Using BLASTn, three of the sequenced samples had a query length of 127–155 bp with an 83–91% query coverage and 92.6–97.9% identity to Vibrio splendidus deposited by Gao (2020) (MT445179.1; MT445177.1) and by Landreau et al. (2020) (MT345091.1) (Table 5). Three other PCR product samples that generated sequences had a query length of 82–179 bp, a query coverage of 78–83% and 92.6–98.6% identity to Vibrio kanaloae deposited by Zheng (2020) (MT505700.1) and by Barcia and Romalde (2020) (CP065151.1; CP065150.1) (Table 5). Another sample, with a query length of 158 bp, had 94% query coverage and 98% identity to V. aestuarianus deposited by Garcia (2018) (MK307696.1) and Vibrio mediterranei deposited by Gao et al. (2020) (MT269602.1) (Table 5). Two other samples were identified as Vibrio spp., as the blasted sequence was quite short (52–69 bp) to differentiate robustly between species (Table 5). The 186 positive PCR cases were defined as Vibrio spp. for statistical and analysis purposes.
Description of the BLASTn results obtained from the sequenced DNA of cockle samples using generic Vibrio primers Vib-F3/R3
Overall, 9.5% (n = 70/735) of the individuals displayed coinfection with the two pathogen groups present in the samples, Haplosporidia and Vibrio (Fig. 2). However, no significant association between them was displayed when testing by SCN (P value > 0.05), under the null hypothesis of random pathogen combinations (Table 6), i.e. the pathogen combination was considered not to occur differently than expected by chance.
Association screening analysis (SCN) results for all the pathogen combinations between Haplosporidia (HAP) and Vibrio (VIB) pathogen groups
The observed (Obs) frequency of each coinfection status is given along with the lower (LL) and upper (UL) limits of the 95% confidence envelope.
In order to identify the mechanisms driving the pathogen coinfection with Haplosporidia and Vibrio pathogen groups, corrected Pearson's χ2 test was also performed including seawater temperature, salinity, cockle height and age, as well as the presence/absence of epiflora on the shells (Tables 7 and 8), as confounding factors which may promote pathogen presence. Under the null hypothesis of independence between the two pathogens, the test showed that the observed frequencies for each pathogen combination were not significantly different (Corr χ 2 = 2.51; ĉ = 0.84; P value 1 = 0.08; P value 2 = 0.08) to the expected frequencies considering the confounding factors as drivers of the infection (Supplementary material). Therefore, the proportion of double infected individuals can be explained because the same factors promote the pathogen presence and not because the different pathogen groups are interacting synergistically. The differences in the considered environmental and biotic factors between sample sites and throughout the time, thus, might have been driving the risk of infection in C. edule populations.
Mean ± s.d. of abiotic (seawater temperature and salinity) and biotic factors (cockle height and age) and the overall percentage of the presence of epiflora on the cockle shells at each sample site
The largest value for each factor is highlighted in bold.
Mean ± s.d. of abiotic (seawater temperature and salinity) and biotic factors (cockle height and age) and the overall percentage of the presence of epiflora on the cockle shells by season
Regarding the specific species quantified in the samples, only 0.8% (n = 6/735) of the individuals were coinfected with V. aestuarianus and M. tapetis, while 0.4% (n = 3/735) were coinfected with V. aestuarianus and M. mercenariae-like (Fig. 2). Likewise, 13.8% of the Minchinia-infected cockles showed a coinfection with both Minchinia species (n = 17/123); however, the percentage fell to 2.3% when all the screened cockles were considered (n = 17/735) (Fig. 2). While coinfection with both Minchinia species and V. aestuarianus was only observed in one of the samples (0.1%, n = 1/735) (Fig. 2).
Associations between M. tapetis, M. mercenariae-like and V. aestuarianus and all the possible combinations were tested by SCN, revealing a significant association between the two species of Minchinia (Table 9). The observed number of coinfected individuals with both Minchinia species were significantly overrepresented compared to expected by random chance (P value < 0.001; Table 9), indicating a significant positive two-way association between these two species. While M. mercenariae-like was found to occur alone significantly less often than expected by chance (P value < 0.05; Table 9).
Association screening analysis (SCN) results for all the pathogen combinations between M. tapetis (MT), M. mercenariae-like (MM) and V. aestuarianus (VA)
The observed (Obs) frequency of each coinfection status is given along with the lower (LL) and upper (UL) limits of the 95% confidence envelope. Significant associations (P value < 0.05) are highlighted in bold.
The association between Minchinia species was further examined by corrected Pearson's χ 2 test in order to assess if the biotic and abiotic factors considered in the study (Tables 7 and 8) influenced this coinfection. The test showed that the observed frequencies for each pathogen combination were significantly different from the expected frequencies considering the confounding factors as drivers of the infection (Supplementary material). The significant outcome (Corr χ 2 = 3.04; ĉ = 0.59; P value 1 = 0.02; P value 2 = 0.02) showed that both pathogens are not independent, and the proportion of double infected individuals cannot be due only to shared confounding factors.
Site influence on single infections and coinfections
Spatial variability was assessed and significant differences [P(Chisq) < 0.05] were found between the sample sites when the prevalence of single infections and two-pathogen coinfections was considered (Fig. 3). Dundalk Bay was the site with the highest presence of Haploporidia spp. (52.5%, n = 126/240), showing a substantially higher prevalence of M. tapetis (23.3%, n = 56/240) than the other sites, while no presence of M. tapetis was observed in Youghal Bay. Dundalk Bay showed lower seawater salinity and a lower presence of epiflora on the cockle shells than the other sample sites (Table 7). Dungarvan Harbour, with the smallest and youngest cockles collected (Table 7), showed the highest presence of Vibrio spp. (44.4%, n = 75/169) and a peak of V. aestuarianus (13.6%, n = 23/169). Moreover, three out of the five individuals with the highest intensity of infection (C t ⩽ 30) with V. aestuarianus were observed at Dungarvan Harbour. Accordingly, the highest prevalence of coinfection with Haplosporidia and Vibrio spp. occurred in Dundalk Bay (14.6%, n = 35/240) and Dungarvan Harbour (11.8%, n = 20/169). However, coinfection with M. tapetis and M. mercenariae-like was absent in Dungarvan Harbour and Youghal Bay, with Cork Harbour being the site with the highest prevalence (4.3%, n = 11/257). The other two-pathogen coinfections with Minchinia species and V. aestuarianus showed no significant differences between the sites [P(Chisq) >0.05] since the total number of positive cases detected was too low [M. tapetis and V. aestuarianus (n = 6); M. mercenariae and V. aestuarianus (n = 3)]. Likewise, Cork Harbour was the only site that showed three-pathogen coinfection with both Minchinia species and V. aestuarianus (0.4%, n = 1/257). Cork Harbour, the site with the oldest cockles, had the highest presence of epiflora on the cockle shells compared to the other sample sites (Table 7).
Prevalence (%) of Haplosporidia spp. (Hap), Vibrio spp. (Vib), Minchinia tapetis (M.tap), M. mercenariae-like (M.mer) and V. aestuarianus (V.aes) single infections and coinfections (Hap_Vib; M.tap_M.mer; M.tap_M.mer_V.aes) at each sample site.
Seasonal impact on single infections and coinfections
Significant seasonal variability [P(Chisq) <0.001] was exhibited when the prevalence of single infections and two-pathogen coinfections were compared between seasons (Fig. 4). The absence of V. aestuarianus and, in general, a low prevalence of infections were observed in winter 2018/19 (Fig. 4), except for Haplosporidia spp. (38.6%, n = 61/158). Summer 2018 had the highest prevalence for all pathogens (Fig. 4), except for M. mercenariae-like that exhibited the highest prevalence in spring 2018 (23.3%, n = 14/60). Furthermore, the five individuals with the highest intensity of infection (C t ⩽ 30) with V. aestuarianus were observed during summer 2018. The highest prevalence of coinfection with Haplosporidia and Vibrio spp. occurred in summer 2018 (21.7%, n = 35/161), while the highest coinfection prevalence with M. tapetis and M. mercenariae-like occurred in spring 2018 (11.7%, n = 7/60), with no occurrence in winter 2018/19. The higher pathogen incidence during spring and summer corresponded with higher seawater temperature and a peak in the presence of epiflora on the cockle shells (Table 8). The other two-pathogen coinfections with Minchinia species and V. aestuarianus showed no significant differences throughout the year [P(Chisq) > 0.05] since the total number of positive cases detected was too low [M. tapetis and V. aestuarianus (n = 6); M. mercenariae and V. aestuarianus (n = 3)]. There was only one case of three-pathogen coinfection with M. tapetis, M. mercenariae-like and V. aestuarianus in spring 2019 (0.6%, n = 1/162).
(A) Prevalence (%) of Haplosporidia (Hap) and Vibrio (Vib) single infections and coinfection (Hap-Vib) by season. (B) Prevalence (%) of Minchinia tapetis (M.tap), M. mercenariae-like (M.mer) and V. aestuarianus (V.aes) single infections and coinfections (M.tap_M.mer and M.tap_M.mer_V.aes) by season.
Histopathology
The histological screening was used as a validation method of the PCR and qPCR results. Haplosporidian infection was observed in histological sections of 67 cockles out of a subsample of 70 samples (95.7% visual confirmation) that were positive for generic haplosporidian PCR. Mostly haplosporidians were observed as multinucleate plasmodia, multiple haplosporidia-like sporonts and developing spores similar to those described for Minchinia spp. (Ramilo et al., Reference Ramilo, Abollo, Villalba and Carballal2018; Lynch et al., Reference Lynch, Lepee-Rivero, Kelly, Quinn, Coghlan, Bookelaar, Morgan, Finarelli, Carlsson and Culloty2020b; Fig. 5), and as uninucleate cells stage in lower numbers. Among the sample sites, the haplosporidia-like sporonts were predominantly visualized in the infected individuals from Cork Harbour in spring 2018, autumn 2018 and spring 2019; while the developing spores were observed mainly in the infected individuals from Dundalk Bay in summer 2018 and winter 2018/19. A higher prevalence of haplosporidians was recorded by PCR at both sites. The multinucleate plasmodia were visualized in cockles throughout the year at the different sample sites. The multinucleate plasmodia and haplosporidia-like sporonts were mainly observed throughout the connective tissue of the gills and digestive gland, while the developing spores were principally found in the mantle of infected animals.
(A) Multiple haplosporidia-like sporonts (black arrow) along with haemocytes (red arrow) in the connective tissue of Cerastoderma edule; (B) a large number of haplosporidian spores (black arrow) in the mantle tissue of C. edule; and (C) red arrow shows a haemocyte accumulation in the connective tissue of C. edule.
Host pathology associated with Vibrio species included haemocyte accumulation in the connective tissues, as McCleary and Henshilwood (Reference McCleary and Henshilwood2015) described in C. gigas infected with V. aestuarianus. In this study, some of the sections examined (n = 28) showed haemocyte accumulation in the connective tissues of C. edule (Fig. 5). However, the presence of haemocyte infiltrations throughout the organs and tissues has also been previously associated with Minchinia infection (Longshaw and Malham, Reference Longshaw and Malham2013). In this study, only in a minor number of screened samples (n = 4), where Vibrio and Haplosporidia were detected by PCR, haemocyte infiltrations were observed, all of them at Cork Harbour.
Although OsHV-1 μVar is not visible using histology, cytopathic effects can be observed in infected hosts, such as hypertrophied cells with a low nucleus-cytoplasm ratio and nuclear abnormalities (Prado-Alvarez et al., Reference Prado-Alvarez, Darmody, Hutton, O'Reilly, Lynch and Culloty2016). Abnormal cell pathology associated with a viral infection was not observed in the histology (n = 70). Similarly, microsporidial spores may be visualized by microscopy in infected hosts (Garcia, Reference Garcia2002); however, it was not observed in the histological sections screened (n = 70).
Ten histological sections of samples that were negative by PCR were screened and no evident signal of infection was found in those, and tissues, organs and sinuses presented good integrity and a clear structure.
Discussion
The present study describes for the first time, to the best of our knowledge, simultaneous infections with haplosporidians and Vibrio spp. in wild and fished C. edule populations along the Irish coast. Coinfections with two and, especially, three pathogens, were rare in this study compared to single infected individuals. This finding may indicate that coinfected individuals succumb more easily to simultaneous infections and die. The lower number of coinfected individuals screened might indicate that the sampling was carried out on a much smaller cohort of survivors within the population. This survival bias is supported by previous studies on mortalities associated with coinfections of ostreid herpesvirus (OsHV-1 and variants) and Vibrio spp. in C. gigas in France (Petton et al., Reference Petton, Bruto, James, Labreuche, Alunno-Bruscia and Le Roux2015; Azéma et al., Reference Azéma, Travers, Benabdelmouna and Dégremont2016; Petton et al., Reference Petton, de Lorgeril, Mitta, Daigle, Pernet and Alunno-Bruscia2019), as well as with the co-occurrence of different protozoa and bacteria in the noble pen shell Pinna nobilis in the Mediterranean Sea (Carella et al., Reference Carella, Elisabetta, Simone, Fulvio, Daniela, Prado, Rossella, Marino, Eleonora, Tobia and De Vico2020). Only three individuals presented a three-pathogen coinfection with Minchinia spp. and Vibrio spp., all located at Cork Harbour. Cork Harbour was the site with the oldest cockles, and it is well known that older hosts have more time than younger hosts to accumulate disease-causing pathogens (Lafferty and Kuris, Reference Lafferty and Kuris2009; Breitburg et al., Reference Breitburg, Hondorp, Audemard, Carnegie, Burrell, Trice and Clark2015). Nevertheless, older bivalves have also more time available to develop immunotolerance and defence mechanisms compared to their younger cohorts (Renault and Novoa, Reference Renault and Novoa2004; Coen and Bishop, Reference Coen and Bishop2015). It may be the reason why cockles at Cork Harbour have not succumbed to the multiple infections.
Even though Haplosporidia and Vibrio were found simultaneously in double infected individuals, no significant biological interaction was statistically identified between those pathogen groups. Unlike other studies that have reported that multiple pathogens can interact synergistically, with the presence of one pathogen enhancing the abundance and/or virulence of the other (De Lorgeril et al., Reference De Lorgeril, Lucasson, Petton, Toulza, Montagnani, Clerissi, Vidal-Dupiol, Chaparro, Galinier, Escoubas, Haffner, Degremont, Charriere, Lafont, Delort, Vergnes, Chiarello, Faury, Rubio, Leroy, Perignon, Regler, Morga, Alunno-Bruscia, Boudry, Le Roux, Destoumieux-Garzon, Gueguen and Mitta2018; Carella et al., Reference Carella, Elisabetta, Simone, Fulvio, Daniela, Prado, Rossella, Marino, Eleonora, Tobia and De Vico2020). Concurrent infections by multiple pathogens species are more likely to occur in immune depressed individuals, which are more vulnerable to multiple opportunistic infections (De Lorgeril et al., Reference De Lorgeril, Lucasson, Petton, Toulza, Montagnani, Clerissi, Vidal-Dupiol, Chaparro, Galinier, Escoubas, Haffner, Degremont, Charriere, Lafont, Delort, Vergnes, Chiarello, Faury, Rubio, Leroy, Perignon, Regler, Morga, Alunno-Bruscia, Boudry, Le Roux, Destoumieux-Garzon, Gueguen and Mitta2018; Carella et al., Reference Carella, Elisabetta, Simone, Fulvio, Daniela, Prado, Rossella, Marino, Eleonora, Tobia and De Vico2020). These findings highlight that there can be variation amongst pathogen groups and the primary pathogen–opportunistic pathogen interplay.
It is well known that environmental and ecological factors potentially influence the mechanisms of disease transmission in marine systems, i.e. hydrodynamics, the biomass of infected animals, etc. (Petton et al., Reference Petton, Bruto, James, Labreuche, Alunno-Bruscia and Le Roux2015). Similarities in the host environment, behaviour or susceptibility can cause correlations in the risk of infection between two parasites (Vaumourin et al., Reference Vaumourin, Vourc'h, Telfer, Lambin, Salih, Seitzer, Morand, Charbonnel, Vayssier-Taussat and Gasqui2014). In this study, the proportion of double infected individuals with Haplosporidia and Vibrio was explained by risk factors (increased seawater temperature, reduced salinity and host condition) shared by the two pathogen groups, increasing the probability of infection by both infectious agents. Such confounding factors may influence host exposure and/or host susceptibility promoting simultaneous infections, even though involved pathogens do not interact biologically (Vaumourin et al., Reference Vaumourin, Vourc'h, Gasqui and Vayssier-Taussat2015). These factors may be naturally occurring or be outcomes of coastal development (e.g. pollution) or climate change (e.g. warming, increased precipitation events and subsequent freshwater loadings into bays and estuaries), and may impact the expression of diseases directly or indirectly (Coen and Bishop, Reference Coen and Bishop2015). The predicted warming, increased precipitation and reduced salinity of coastal waters in temperate regions will provide new areas for the natural occurrence of pathogenic strains and may also impact the resilience of bivalve hosts (Rowley et al., Reference Rowley, Cross, Culloty, Lynch, Mackenzie, Morgan, O'Riordan, Robins, Smith, Thrupp, Vogan, Wootton and Malham2014; Lynch et al., Reference Lynch, Coghlan, O'Leary, Morgan and Culloty2020a).
Consequently, spatial variability in the pathogen infection was detected driven by environmental and biological variables. The highest prevalence of coinfection with Haplosporidia and Vibrio spp. occurred in Dundalk Bay, which had a lower mean salinity throughout the year compared to the other sites. Gorrasi et al. (Reference Gorrasi, Pasqualetti, Franzetti, Pittino and Fenice2020) observed that the abundance of Vibrio genus in water decreased along the salinity gradient within a marine saltern hypersaline environment. In general, marine vibrios are considered halotolerant, with no extreme halophiles reported within this genus (Gorrasi et al., Reference Gorrasi, Pasqualetti, Franzetti, Pittino and Fenice2020). Previous studies have highlighted the intolerance of haplosporidians to low salinity in C. edule (Albuixech-Martí et al., Reference Albuixech-Martí, Lynch and Culloty2020), eastern oyster Crassostrea virginica (Ford and Haskin, Reference Ford and Haskin1982; Burreson and Ford, Reference Burreson and Ford2004) and European flat oyster Ostrea edulis (Flannery et al., Reference Flannery, Lynch, Longshaw, Stone, Martin, Ramilo, Villalba and Culloty2014). The high prevalence of haplosporidians at Dundalk Bay may be due to its commercial fishery activity. It is recognized that fishing and dredging may unbalance the interplay between infectious disease agents and their hosts in marine ecosystems (Rowley et al., Reference Rowley, Cross, Culloty, Lynch, Mackenzie, Morgan, O'Riordan, Robins, Smith, Thrupp, Vogan, Wootton and Malham2014; Coen and Bishop, Reference Coen and Bishop2015). Cranfield et al. (Reference Cranfield, Michael and Doonan1999, Reference Cranfield, Dunn, Doonan and Michael2005) hypothesized that extensive dredging in Foveaux Strait of New Zealand contributed to increased densities of Chilean oysters (Ostrea chilensis) and decreased densities of other filter-feeding organisms, that may have previously reduced the dispersal stages of the haplosporidian Bonamia exitiosa, enhancing disease transmission. The high cockle density found at Dundalk Bay may have facilitated the infection among the cockle population, due to the proximity of infected to non-infected cockles. The cultivated and harvested molluscs tend to have a much greater incidence of disease than wild molluscs, which may be because they are placed at high densities which enhance disease transmission and cause stress (competition for resources, space, etc.) among individuals (Ford, Reference Ford2002).
Dungarvan Harbour also displayed a high prevalence of coinfection with Haplosporidia and Vibrio spp. compared to the other sites. Cockles from Dungarvan Harbour were the smallest and the youngest ones. Although it is accepted that older hosts have more time to accumulate disease-causing pathogens (Lafferty and Kuris, Reference Lafferty and Kuris2009; Breitburg et al., Reference Breitburg, Hondorp, Audemard, Carnegie, Burrell, Trice and Clark2015), the susceptibility of small bivalves to viral and bacterial infections is generally greater than for larger adults (Renault and Novoa, Reference Renault and Novoa2004; Coen and Bishop, Reference Coen and Bishop2015). Moreover, the proximity of an extensive C. gigas culture site at Dungarvan Harbour, where Vibrio spp. are endemic (pers. comm. with Marine Institute Ireland), can have promoted the presence of Vibrio in the environment. Particularly, V. aestuarianus and V. splendidus, which have been identified in this study, have been described as increasingly harmful pathogens of C. gigas aquaculture (Le Roux et al., Reference Le Roux, Gay, Lambert, Waechter, Poubalanne, Chollet, Nicolas and Berthe2002; Garnier et al., Reference Garnier, Labreuche, Garcia, Robert and Nicolas2007; McCleary and Henshilwood, Reference McCleary and Henshilwood2015). Caballes et al. (Reference Caballes, Schupp, Pratchett and Rivera-Posada2012) demonstrated interspecific transmission of Vibrio rotiferianus between co-occurring echinoderm species.
Even though the natural host for OsHV-1 are C. gigas, this virus can behave opportunistically and infect other cohabiting bivalve species (Bookelaar et al., Reference Bookelaar, Lynch and Culloty2020). It has been previously reported in C. edule, in Ireland (Bookelaar et al., Reference Bookelaar, Lynch and Culloty2020) and in the Sydney cockle, Anadara trapezia, in Australia (Evans et al., Reference Evans, Paul-Pont and Whittington2017). Nevertheless, no presence of OsHV-1 μVar was detected in the study, despite the extensive farming of C. gigas at Dungarvan Harbour, which has a history of OsHV-1-associated mortality events (Lynch et al., Reference Lynch, Carlsson, Reilly, Cotter and Culloty2012; Prado-Alvarez et al., Reference Prado-Alvarez, Darmody, Hutton, O'Reilly, Lynch and Culloty2016). During summer 2018, the prevalence of OsHV-1 μVar in C. gigas at Dungarvan Harbour was confirmed by Kett et al. (unpublished data) to be 48.4% (n = 92/190), with low mortality rates (max. of 6% at the high shore in July). Moreover, Bookelaar et al. (Reference Bookelaar, Lynch and Culloty2020) observed OsHV-1μVar-infected C. edule (14.4%; n = 36/250) and infected C. gigas (range of 0–27% per month; n = 270) at Dungarvan Harbour from April 2015 to August 2015, with no recording of mortalities. Thus, the lack of OsHV-1 μVar detection in C. edule at Dungarvan Harbour in the current study may be because the virus remained within the oysters and was not shed from dying or dead oysters into the environment, as the high prevalence and low mortality recorded by Kett et al. (unpublished data) suggested. The smaller sample size (n = 81), compared to the above studies, taken during spring and summer 2018, when the OsHV-1 μVar can have a major emergence due to the higher temperatures (Renault et al., Reference Renault, Bouquet, Maurice, Lupo and Blachier2014), may have also influenced the results.
Seasonal variability was also observed in this study, with the highest prevalence of coinfection with Haplosporidia and Vibrio spp. occurring in summer 2018, which coincided with the highest intensity of infection with V. aestuarianus. Additionally, lower occurrences of the haplosporidian and Vibrio species quantified were detected in the cooler spring 2019 compared to spring 2018 in this study. Burreson and Ford (Reference Burreson and Ford2004) tested that low-temperature conditions caused a dramatic reduction in the prevalence of Haplosporidium nelsoni in C. virginica. Conversely, a higher prevalence of Bonamia ostreae in European flat oyster O. edulis was associated with cooler spring temperatures over a 30-year study period in Ireland (Engelsma et al., Reference Engelsma, Kerkhoff, Roozenburg, Haenen, van Gool, Sistermans, Wijnhoven and Hummel2010; Lynch et al., Reference Lynch, Morgan, Carlsson, Mackenzie, Wooton, Rowley, Malham and Culloty2014). These findings would suggest that a thermal tolerance range exists for haplosporidian species. It is also well-known that bacteria belonging to the genus Vibrio are strongly thermo-dependent and their occurrence is expected to rise with increasing temperature conditions (Garnier et al., Reference Garnier, Labreuche, Garcia, Robert and Nicolas2007; Cotter et al., Reference Cotter, Malham, O'Keeffe, Lynch, Latchford, King, Beaumont and Culloty2010; Vezzulli et al., Reference Vezzulli, Previati, Pruzzo, Marchese, Bourne, Cerrano and VibrioSea Consortium2010). Most of the reported Vibrio-associated mass mortalities in marine invertebrates have occurred during the summer when seawater temperatures were higher than 18°C (Le Roux et al., Reference Le Roux, Gay, Lambert, Waechter, Poubalanne, Chollet, Nicolas and Berthe2002; Paillard et al., Reference Paillard, Allam and Oubella2004; Garnier et al., Reference Garnier, Labreuche, Garcia, Robert and Nicolas2007; Vezzulli et al., Reference Vezzulli, Previati, Pruzzo, Marchese, Bourne, Cerrano and VibrioSea Consortium2010). Therefore, high temperatures in summer might have triggered a major pathogen propagation and, consequently, to a higher prevalence of two-pathogen coinfection. While lower temperatures in winter, with the lowest prevalence of two-pathogen coinfection, would have kept the pathogens at a sub-lethal level in the hosts reducing the transmission between cockles. It has been observed that Vibrio spp. may enter into a dormant state when the environmental conditions become unfavourable (i.e. suboptimal or reduced temperature, elevated salinity, low nutrient concentration and extreme pH) (Lin and Schwarz, Reference Lin and Schwarz2003; Nowakowska and Oliver, Reference Nowakowska and Oliver2013; Neogi et al., Reference Neogi, Lara, Alam, Harder, Yamasaki and Colwell2018). The low number of coinfected individuals in winter might also suggest that most C. edule had succumbed to infection in autumn after the peak in summer.
It was also observed that with the higher temperatures in summer and spring, the presence of epiflora on the cockle shells was more frequent, which may have consequences on cockle susceptibility to infectious agents. As seen in shellfish aquaculture, the presence of biofouling on the stock can have a direct impact, causing physical damage, biological competition, mechanical interference and environmental modification, while infrastructure is also impacted (Fitridge et al., Reference Fitridge, Dempster, Guenther and de Nys2012). Microalgae have been pointed as a repetitive and substantial source of bacterial inoculation into the bivalve larval culture systems (Salvesen et al., Reference Salvesen, Reitan, Skjermo and Oie2000). Microalgae may promote and/or inhibit bacterial growth by the production of organic exudates and toxic metabolites (Salvesen et al., Reference Salvesen, Reitan, Skjermo and Oie2000). For instance, Vibrio species such as V. splendidus or V. parahaemolyticus have been found to be associated with microalgae (Kumazawa et al., Reference Kumazawa, Nakagaki, Yonekawa, Ikura and Morimoto1991; Rehnstam-Holm et al., Reference Rehnstam-Holm, Godhe, Harnstrom, Raghunath, Saravanan, Collin, Karunasagar and Karunasagar2010; Notaro et al., Reference Notaro, Culloty and Lynch2021). Therefore, the presence of epiflora on the shells might have also exposed the cockles to pathogens, in particular Vibrio species.
Both species of Minchinia, M. tapetis and M. mercenariae-like, were found to occur together significantly more often than expected by random chance, indicating a positive two-way association between both species. However, previous studies examining pairwise coinfections indicate that antagonistic interactions between similar types of parasites may be common (Dobson, Reference Dobson1985; Johnson and Hoverman, Reference Johnson and Hoverman2012). A possible explanation of the established synergetic interaction is that, despite both Minchinia species being closely related, they exhibit differences in their pathogenic effect and use of resources within the host. Carballal et al. (Reference Carballal, Cao, Iglesias, Gonzalez and Villalba2020) observed that the distribution of infection with M. tapetis and M. mercenariae-like haplosporidian in cockles C. edule in Galicia was different throughout the cockle tissues. M. mercenariae-like haplosporidian was observed in the connective tissue of the digestive gland, gills and gonad, while M. tapetis was only observed in the digestive gland (Carballal et al., Reference Carballal, Cao, Iglesias, Gonzalez and Villalba2020). Further research is needed to reveal the underlying mechanism of this synergetic association in C. edule and the role that host exposure and host habitat play in it.
In brief, coinfected individuals may more easily succumb to pathogen load. The proportion of double infected individuals with haplosporidian and Vibrio was driven by environmental and biological factors, triggering a spatial and seasonal variability in the distribution of coinfection. The higher temperature in warmer months and the presence of epiflora, which can be associated with pathogens, on the cockle shells may have promoted the risk of coinfection, along with low salinity and the host condition in some of the sample sites. The close proximity to other infected host species may have also influenced the coinfection distribution in C. edule. Similarly, fished C. edule populations with high cockle densities might have been more exposed to the infection. Therefore, those factors may conform the ideal scenario for simultaneous infections, although the pathogen groups showed no biological interaction between them. Nevertheless, a positive two-way association between both Minchinia species was observed, which is suggested to be due to their different pathogenic effect and use of resources within the host. Our findings shed light on the complex interactions between multiple aetiological agents associated with host diseases in the frame of a natural system. This study highlights the consequences of the coinfections are not simply the sum of the effects caused by each infectious agent, but the result of a complex combination of direct and indirect pathogens effects, the host ecology and the environmental context.
These results are framed into the current scenario of climate change, which involves oceans warming on average by 0.6°C in the next 100 years (Wiltshire et al., Reference Wiltshire, Malzahn, Wirtz, Greve, Janisch, Mangelsdorf, Manly and Boersma2008), with European seas especially affected (Philippart et al., Reference Philippart, Anadon, Danovaro, Dippner, Drinkwater, Hawkins, Oguz, O'Sullivan and Reid2011). The sea surface temperature around the UK and Ireland has increased at six times the rate of the global average, particularly jeopardizing species that have small thermal tolerance windows (Frost et al., Reference Frost, Baxter, Buckley, Cox, Dye and Harvey2012). As a result, previously appropriate habitats for cockles may no longer be suitable due to changes in thermal regime, oxygen levels or altered rainfall patterns (Morgan et al., Reference Morgan, O'Riordan and Culloty2013; Coen and Bishop, Reference Coen and Bishop2015). Climate change may be also expected to have a large effect on disease occurrence, especially in host species showing strong relationships between temperature/salinity and diseases, as seen for the cockle populations studied, and for many other molluscan host–parasite systems (Allam and Raftos, Reference Allam and Raftos2015; Coen and Bishop, Reference Coen and Bishop2015; Lynch et al., Reference Lynch, Coghlan, O'Leary, Morgan and Culloty2020a). The physiological stress may also cause the host to become immune-compromised and unable to suppress an infection (Suttle, Reference Suttle2007). Therefore, climate change could be the driver of the spreading and range expansions of infectious agents, such as bacteria, haplosporidians and viruses, in the marine ecosystem and may lead to more frequent pathogen outbreaks (Burreson and Ford, Reference Burreson and Ford2004; Suttle, Reference Suttle2007; Rowley et al., Reference Rowley, Cross, Culloty, Lynch, Mackenzie, Morgan, O'Riordan, Robins, Smith, Thrupp, Vogan, Wootton and Malham2014; Coen and Bishop, Reference Coen and Bishop2015). As global climate change and anthropogenic impacts continue to modify the distribution and abundance of hosts and parasites and the ways in which they interact, integral studies like ours are important since concepts gained might be applied to other associations and to a changing environment.
A final important point to emphasize is that the relative potential of coinfections to cause the emergence of zoonoses in wild animals remains to be fully assessed (Hoarau et al., Reference Hoarau, Mavingui and Lebarbenchon2020). Given that two-thirds of emerging infectious diseases are zoonoses, with nearly 70% originating from wildlife (Hoarau et al., Reference Hoarau, Mavingui and Lebarbenchon2020), better knowledge of the interactions of infectious agents in wild reservoirs, such as cockles, will provide key insight for the understanding and management of spillover processes. Although challenging, findings from this study highlight the necessity to integrate data pertaining to host ecology, pathogen diversity, biogeography and seasonality of the disease and the environmental influence. This approach should be consistently applied to provide a more detailed picture of coinfection dynamics in individual hosts and communities over time.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/S0031182021001396.
Acknowledgements
We are grateful to Dr Eileen Dillane of the University College Cork (Ireland) (Research Support Officer and Genetic Laboratory Manager) for the provision of professional advice during our lab work. We also thank our colleagues Dr Kate Mahony, Dr Katie Costello and Gary Kett of University College Cork (Ireland), who assisted us with the transport, and fieldwork. For the permission to carry out the fieldwork in Dungarvan Harbour, we thank Dungarvan Shellfish Ltd (Ireland); we also thank Martin Hooey from Dundalk Fishery (Ireland) for showing us the best sampling areas in Dundalk Bay and to provide us samples.
Author contributions
Conceptualization and methodology: S.A.M., S.C.C., S.A.L.; formal analysis: S.A.M.; writing – original draft: S.A.M.; writing – review and editing: S.C.C., S.A.L.; supervision: S.C.C., S.A.L. All authors approved the final version of the manuscript before submission.
Financial support
This work was supported by the Bluefish (Grant Agreement No. 80991) Project, part-funded by the European Regional Development Fund (ERDF) through the Ireland Wales Co-operation Programme. A cross-border programme investing in the overall economic, environmental and social well-being of Ireland and Wales.
Conflict of interest
None.
Ethical standards
Not applicable.
Open access
