A rapid phenotype change in the pathogen Perkinsus marinus was associated with a historically significant marine disease emergence in the eastern oyster

Co-infection is linked to infection prevalence and intensity in oysters amidst high environmental and spatial variation

Co-infecting parasites modify infection outcomes in the wild. However, it is unclear how multiple environmental factors influence co-infection. The Chesapeake Bay metapopulation of the eastern oyster, Crassostrea virginica, provides an opportunity to test the importance of co-infection across heterogeneous environments because multiple parasites infect oysters across a broad salinity gradient. This study leverages Maryland and Virginia oyster monitoring for a large-scale survey of four co-infecting organisms, including two tissue parasites and two shell bio-eroding parasites. We diagnosed infection in 440 oysters across 16 paired harvested and unharvested reefs and tested the importance of co-infecting organisms for each parasite relative to environmental conditions, host traits, and marine spatial management. Microscopic visual methods were used to diagnose prevalence and intensity of tissue infections with Perkinsus marinus (the causative agent of dermo disease) and Haplosporidium nelsoni (the causative agent of MSX disease). Macroscopic visual methods were used to diagnose prevalence and intensity of shell infections with Cliona boring sponges and blister-inducing Polydora worms. For the three oyster parasites that were detected [H. nelsoni infections were absent in all oysters], salinity was the overall strongest predictor, corresponding to bay-wide patterns of parasite prevalence and/or intensity. Despite high environmental and spatial variation, co-infections corresponded to altered prevalence and/or intensity for all three oyster parasites. The correlational patterns suggest that P. marinus acts as a lynchpin in co-infection, as its intensity increased with Cliona sponge prevalence and P. marinus co-infection predicted higher Polydora blister intensity. Oyster shell height, reef habitat, and harvest status also predicted parasite prevalence and intensity, further reflecting the multivariate drivers of infections in this system. Unharvested reefs had greater vertical habitat structure and higher intensities of Cliona sponge infections, but no differences in the prevalence of any of the three parasites. Spatial patterns unexpectedly show that reef-level predictors of parasite patterns were more important than differences between tributaries. This correlational survey provides novel insights through the statistical relationships between the three oyster parasites, environmental conditions, host traits, and human resource management. New and more detailed scenarios are needed to expand disease ecological theory to encompass co-infection in anthropogenically impacted wildlife populations.

High parasite prevalence in an ecosystem engineer correlated with both local- and landscape-level factors

Spatial variation in parasitic infection may have many physical and biological drivers. Uncovering these drivers may be especially important for parasites of ecosystem engineers because the engineers are foundational to their communities. Oysters are an important coastal ecosystem engineer that have declined drastically worldwide, in part due to enhanced cases of lethal oyster diseases, such as Dermo and MSX, caused by the protozoan parasites Perkinsus marinus and Haplosporidium nelsoni, respectively. Besides water quality and hydrodynamics, there is little information on how other variables influence the prevalence and intensity of these pathogens in oysters across a regional scale. To examine drivers of spatial variation in these oyster parasites-including host size, local reef properties, and landscape properties-we sampled 24 reefs systematically spread along the coast of Georgia, USA. Across sites, we found universally high prevalence of oysters with at least one of these parasites (91.02% ± 8.89, mean ± SD). Not only are high levels of parasite prevalence potentially problematic for a pivotal ecosystem engineer, but also low spatial variability may limit the explanatory power of variables across a regional scale. Our statistical models explained between 18 and 42% of the variation in spatial patterns of prevalence and intensity of these microparasites. Interestingly, landscape context was a positive predictor of P. marinus, but a negative predictor of H. nelsoni. Overall, our findings suggest that factors driving parasite prevalence and intensity operate across multiple spatial scales, and the same factor can both facilitate and hinder different parasites within the same host species.

Occurrence and seasonal variation of Perkinsus sp. Infection in wild mollusk populations from coastal waters of Qingdao, northern China

Perkinsosis has been recognized as one of the major threats to natural and farmed bivalve populations, many of which are of commercial as well as environmental significance. Three Perkinsus species have been identified in China, and the Manila clam (Ruditapes philippinarum) was the most frequently infected species in northern China. Although the occurrence and seasonal variation of Perkinsus spp. have previously been examined, the pathological characteristics of these infections in wild Manila clams and sympatric species in China have seldom been reported. In the present study, the prevalence and intensity of Perkinsus infection in wild populations of Manila clams and 10 sympatric species from three sites were investigated by Ray's fluid thioglycolate medium (RFTM) assay seasonally across a single year. Perkinsus infection was only identified in Manila clams, with a high prevalence (274/284 = 96.48 %) and low intensity (89.8 % with a Mackin value ≤ 2, suggesting generally low-intensity infections) throughout the year. Heavily infected clams were mainly identified in Tianheng in January, which displayed no macroscopic signs of disease. An overview of the whole visceral mass section showed that the trophozoites mostly aggregated in gills and connective tissue of the digestive tract, to a lesser extent in the mantle and foot, and even less frequently in adductor muscle and connective tissues of the gonad. PCR and ITS-5.8S rRNA sequencing of 93 representative RFTM-positive samples revealed a 99.69 to 100 % DNA sequence identity to Perkinsus olseni. Unexpectedly, significantly higher infection intensities were usually identified in January and April when the Condition Index (CI) was relatively high. We propose that factors associated with the anthropogenic harvesting pressure and irregular disturbances should be responsible for the uncommon seasonal infection dynamics of perkinsosis observed in the present study.

Into the danger zone: How the within-host distribution of parasites controls virulence

Despite the importance of virulence in epidemiological theory, the relative contributions of host and parasite to virulence outcomes remain poorly understood. Here, we use reciprocal cross experiments to disentangle the influence of host and parasite on core virulence components-infection and pathology-and understand dramatic differences in parasite-induced malformations in California amphibians. Surveys across 319 populations revealed that amphibians' malformation risk was 2.7× greater in low-elevation ponds, even while controlling for trematode infection load. Factorial experiments revealed that parasites from low-elevation sites induced higher per-parasite pathology (reduced host survival and growth), whereas there were no effects of host source on resistance or tolerance. Parasite populations also exhibited marked differences in within-host distribution: ~90% of low-elevation cysts aggregated around the hind limbs, relative to <60% from high-elevation. This offers a novel, mechanistic basis for regional variation in parasite-induced malformations while promoting a framework for partitioning host and parasite contributions to virulence.

Morphological and functional characterization of the oyster Crassostrea gasar circulating hemocytes: Cell types and phagocytosis activity

Hemocytes are the circulating cells of the hemolymph of oysters and are responsible for numerous physiological functions, including immune defense. The oyster Crassostrea gasar is a native species inhabiting mangrove habitat and is of great commercial interest, cultured throughout the Brazilian coast, mainly in the north and northeast. Despite its commercial importance, little is known about its immunological aspects and defense cells, the hemocytes. This work aimed to morphologically characterize hemocytes of the oyster C. gasar and to study one of the main cellular defense response, phagocytosis, using light microscopy and flow cytometry. The results showed the presence of six hemocyte populations in C. gasar hemolymph. These comprise of large and small granulocytes, large and small hyalinocytes, blast-like cells and a rare type classified as vesicular or serous hemocytes. Hyalinocytes were highly abundant and the most heterogeneous cell population, while small granulocytes, along with vesicular hemocytes were the less abundant population. Hemocytes of C. gasar oysters demonstrated capabilities to phagocytose three different types of particles tested: zymosan A, latex particles and Escherichia coli, indicating a broad defense capacity. The zymosan A were the most engulfed particles, followed by beads, mainly phagocytized by granulocytes, the most phagocytic cells, and finally E. coli, which were the least phagocytized. This study is the first characterization of C. gasar oyster hemocytes and will support future studies that aim to understand the participation of different hemocyte types in defense responses against pathogens and/or environmental changes.

Decline of Marteilia cochillia in Ría de Arousa may be due to increased resistance in host Cerastoderma edule

A huge, unprecedented mortality of cockle Cerastoderma edule caused by the protist Marteilia cochillia, which had never before been detected in Galicia (NW Spain), brought on a cockle fishery collapse in the Ría de Arousa (Galicia) in 2012. Since then, the disease dynamic pattern in the shellfish bed of Lombos do Ulla (at the inner area of that ria) involved an overwhelming annual wave of infections and subsequent cockle mass mortality that caused the near extinction of every cohort recruited to that bed. However, a pattern shift was detected among wild cohorts recruiting since 2016, with progressive declines of marteiliosis prevalence and increments in cockle survival. This suggested 2 non-exclusive hypotheses: increasing marteiliosis resistance through natural selection, and reduced abundance and/or virulence of the parasite. A field experiment was performed to assess these hypotheses by comparing marteiliosis prevalence and severity, as well as mortality, in cockles that naturally recruited to this bed in 2017 and 2018 with those of naïve cockles collected from a marteiliosis-free area and transplanted into Lombos do Ulla in 2017 and 2018. Marteiliosis prevalence and cumulative cockle mortality quickly reached very high values among the transplanted cockles, demonstrating that the parasite remained present and virulent in the area. Conversely, marteiliosis prevalence and cockle mortality were much lower in the cockles that recruited to Lombos do Ulla, suggesting increased resistance that may have been driven by natural selection. The young age at which cockles start reproduction and the very high mortality caused by marteiliosis may have enhanced natural selection.

Host genetic identity determines parasite community structure across time and space in oyster restoration

Intraspecific variation in host susceptibility to individual parasite species is common, yet how these effects scale to mediate the structure of diverse parasite communities in nature is less well understood. To address this knowledge gap, we tested how host genetic identity affects parasite communities on restored reefs seeded with juvenile oysters from different sources-a regional commercial hatchery or one of two wild progenitor lines. We assessed prevalence and intensity of three micro- and two macroparasite species for 4 years following restoration. Despite the spatial proximity of restored reefs, oyster source identity strongly predicted parasite community prevalence across all years, with sources varying in their relative susceptibility to different parasites. Oyster seed source also predicted reef-level parasite intensities across space and through time. Our results highlight that host intraspecific variation can shape parasite community structure in natural systems, and reinforce the importance of considering source identity and diversity in restoration design.

A reconstruction of parasite burden reveals one century of climate-associated parasite decline

Long-term data allow ecologists to assess trajectories of population abundance. Without this context, it is impossible to know whether a taxon is thriving or declining to extinction. For parasites of wildlife, there are few long-term data-a gap that creates an impediment to managing parasite biodiversity and infectious threats in a changing world. We produced a century-scale time series of metazoan parasite abundance and used it to test whether parasitism is changing in Puget Sound, United States, and, if so, why. We performed parasitological dissection of fluid-preserved specimens held in natural history collections for eight fish species collected between 1880 and 2019. We found that parasite taxa using three or more obligately required host species-a group that comprised 52% of the parasite taxa we detected-declined in abundance at a rate of 10.9% per decade, whereas no change in abundance was detected for parasites using one or two obligately required host species. We tested several potential mechanisms for the decline in 3+-host parasites and found that parasite abundance was negatively correlated with sea surface temperature, diminishing at a rate of 38% for every 1 °C increase. Although the temperature effect was strong, it did not explain all variability in parasite burden, suggesting that other factors may also have contributed to the long-term declines we observed. These data document one century of climate-associated parasite decline in Puget Sound-a massive loss of biodiversity, undetected until now.

Development of a novel electroporation method for the oyster parasite Perkinsus marinus

Gene manipulation techniques are fundamental to molecular biology and are continuously being improved. However, gene transfection methods are not established for many unicellular eukaryotes (protists), thereby hindering molecular biological investigations. The oyster parasite Perkinisus marinus is one of the few protists with established gene transfection and drug selection. Nevertheless, the present protocols are tedious, requiring a specific electroporator and pulse conditions which limits the accessibility of this technique across different research groups. Here, we present alternative buffer and electroporation conditions that make the protocol less restrictive. We revealed the pulse condition that enables the introduction of plasmids into P. marinus cell using Ingenio electroporation buffer and NEPA21 electroporator. We found that number of cells and plasmid concentration were critical parameters for the electroporation system. We also constructed a simpler expression plasmid that is removed needless regions for gene expression in the parasite. Our findings resolved the equipment restriction in electroporation of P. marinus and would be a good reference for electroporation in other protists, in particular other Perkinsozoa parasites and core dinoflagellates.

Development of a multiplex qPCR for the quantification of three protozoan parasites of the eastern oyster Crassostrea virginica

A multiplex quantitative PCR (qPCR) assay for the simultaneous detection of 3 eastern oyster Crassostrea virginica parasites, Perkinsus marinus, Haplosporidium nelsoni, and H. costale, was developed using 3 different fluorescently labeled hydrolysis probes. The primers and probe from a previously validated singleplex qPCR for P. marinus detection were combined with newly designed primers and probes specific for H. nelsoni and H. costale. The functionality of the multiplex assay was demonstrated on 2 different platforms by the linear relationship of the standard curves and similar cycle threshold (CT) values between parasites. Efficiency of the multiplex qPCR assay on the Roche and BioRad platforms ranged between 93 and 101%. The sensitivity of detection ranged between 10 and 100 copies of plasmid DNA for P. marinus and Haplosporidium spp., respectively. The concordance between the Roche and BioRad platforms in the identification of the parasites P. marinus, H. nelsoni, and H. costale was 91, 97, and 97%, respectively, with a 10-fold increase in the sensitivity of detection of Haplosporidium spp. on the BioRad thermocycler. The concordance between multiplex qPCR and histology for P. marinus, H. nelsoni, and H. costale was 54, 57, and 87%, respectively. Discordances between detection methods were largely related to localized or low levels of infections in oyster tissues, and qPCR was the more sensitive diagnostic. The multiplex qPCR developed here is a sensitive diagnostic tool for the quantification and surveillance of single and mixed infections in the eastern oyster.

Emerging Parasitic Protists: The Case of Perkinsea

The last century has witnessed an increasing rate of new disease emergence across the world leading to permanent loss of biodiversity. Perkinsea is a microeukaryotic parasitic phylum composed of four main lineages of parasitic protists with broad host ranges. Some of them represent major ecological and economical threats because of their geographically invasive ability and pathogenicity (leading to mortality events). In marine environments, three lineages are currently described, the Parviluciferaceae, the Perkinsidae, and the Xcellidae, infecting, respectively, dinoflagellates, mollusks, and fish. In contrast, only one lineage is officially described in freshwater environments: the severe Perkinsea infectious agent infecting frog tadpoles. The advent of high-throughput sequencing methods, mainly based on 18S rRNA assays, showed that Perkinsea is far more diverse than the previously four described lineages especially in freshwater environments. Indeed, some lineages could be parasites of green microalgae, but a formal nature of the interaction needs to be explored. Hence, to date, most of the newly described aquatic clusters are only defined by their environmental sequences and are still not (yet) associated with any host. The unveiling of this microbial black box presents a multitude of research challenges to understand their ecological roles and ultimately to prevent their most negative impacts. This review summarizes the biological and ecological traits of Perkinsea-their diversity, life cycle, host preferences, pathogenicity, and highlights their diversity and ubiquity in association with a wide range of hosts.

Transcriptomic Response to Perkinsus marinus in Two Crassostrea Oysters Reveals Evolutionary Dynamics of Host-Parasite Interactions

Infectious disease outbreaks are causing widespread declines of marine invertebrates including corals, sea stars, shrimps, and molluscs. Dermo is a lethal infectious disease of the eastern oyster Crassostrea virginica caused by the protist Perkinsus marinus. The Pacific oyster Crassostrea gigas is resistant to Dermo due to differences in the host-parasite interaction that is not well understood. We compared transcriptomic responses to P. marinus challenge in the two oysters at early and late infection stages. Dynamic and orchestrated regulation of large sets of innate immune response genes were observed in both species with remarkably similar patterns for most orthologs, although responses in C. virginica were stronger, suggesting strong or over-reacting immune response could be a cause of host mortality. Between the two species, several key immune response gene families differed in their expansion, sequence variation and/or transcriptional response to P. marinus, reflecting evolutionary divergence in host-parasite interaction. Of note, significant upregulation of inhibitors of apoptosis (IAPs) was observed in resistant C. gigas but not in susceptible C. virginica, suggesting upregulation of IAPs is an active defense mechanism, not a passive response orchestrated by P. marinus. Compared with C. gigas, C. virginica exhibited greater expansion of toll-like receptors (TLRs) and positive selection in P. marinus responsive TLRs. The C1q domain containing proteins (C1qDCs) with the galactose-binding lectin domain that is involved in P. marinus recognition, were only present and significantly upregulated in C. virginica. These results point to previously undescribed differences in host defense genes between the two oyster species that may account for the difference in susceptibility, providing an expanded portrait of the evolutionary dynamics of host-parasite interaction in lophotrochozoans that lack adaptive immunity. Our findings suggest that C. virginica and P. marinus have a history of coevolution and the recent outbreaks may be due to increased virulence of the parasite.

Disease-driven mass mortality event leads to widespread extirpation and variable recovery potential of a marine predator across the eastern Pacific

The prevalence of disease-driven mass mortality events is increasing, but our understanding of spatial variation in their magnitude, timing and triggers are often poorly resolved. Here, we use a novel range-wide dataset comprised 48 810 surveys to quantify how sea star wasting disease affected Pycnopodia helianthoides, the sunflower sea star, across its range from Baja California, Mexico to the Aleutian Islands, USA. We found that the outbreak occurred more rapidly, killed a greater percentage of the population and left fewer survivors in the southern half of the species's range. Pycnopodia now appears to be functionally extinct (greater than 99.2% declines) from Baja California, Mexico to Cape Flattery, Washington, USA and exhibited severe declines (greater than 87.8%) from the Salish Sea to the Gulf of Alaska. The importance of temperature in predicting Pycnopodia distribution rose more than fourfold after the outbreak, suggesting latitudinal variation in outbreak severity may stem from an interaction between disease severity and warmer waters. We found no evidence of population recovery in the years since the outbreak. Natural recovery in the southern half of the range is unlikely over the short term. Thus, assisted recovery will probably be required to restore the functional role of this predator on ecologically relevant time scales.