Abstract
Zooplankton grazers, like copepods, can feed on toxic microalgae and live normally. We hypothesize that gut microbial communities (GMCs) may contribute to the detoxification of the host by changing their compositions and recruiting more beneficial bacteria. Here, we measured the physiological responses of two copepod species (Acartia sp. and Paracalanus sp.) fed with toxic (Alexandrium tamarense) and non-toxic (Alexandrium andersonii) dinoflagellates, respectively. Both copepods maintained consistently high survival rates but slightly reduced ingestion rates when feeding upon the toxic dinoflagellate (when compared to the non-toxic one), suggesting a compensatory mechanism. The compositional variation of copepod GMCs, at the amplicon sequence variant (ASV) level, was mostly significantly different among copepod host species (R = 0.83, by ANOSIM test), while diet type played minor but significant roles. Under the toxic diet, Acartia sp. enriched only five ASVs while Paracalanus sp. recruited a wide range of taxa (38 ASVs) mostly belonging to Alphaproteobacteria (e.g., Rhodobacteraceae) and Gammaproteobacteria (e.g., Alteromonadaceae). In contrast, when clustering GMCs by predicted functions, diet type was the key regulating factor, suggesting the functional convergence of copepod GMCs in response to algal toxins. This can be explained by the fact that most of the enriched bacteria under the toxic diet have similar functions on detoxification and maintaining the host homeostasis. This study deepens our understanding of the roles of GMC in the detoxification and adaptation mechanisms of copepods during harmful algal blooms.
1 Introduction
Harmful algal blooms (HABs) caused by dinoflagellates are a significant global issue that can have serious ecological and economic impacts. Ecologically, dinoflagellate HABs can lead to the depletion of oxygen levels in the water, the creation of dead zones, and the alteration of food webs (Turner, 2014). Economically, HABs can lead to the closure of fisheries, the loss of tourism revenue, and the costs associated with monitoring and mitigating HAB outbreaks (). For example, the high PSP toxins levels produced by toxic Alexandrium tamarense caused high mortality of caged salmon in Nova Scotia (). Additionally, the transfer of PSP toxins through the food chain has put the western North Atlantic right whale Eubalaena glacialis at risk of reduced reproduction rates (). Humans have also been impacted by this issue, with reports of poisonings from consuming crabs contaminated with PSP toxins (). Several marine dinoflagellate species such as Alexandrium spp., Pyrodinium spp., and Gymnodinium spp., can produce paralytic shellfish poisoning (PSP) toxins (; ). PSP is a foodborne illness that results from consuming seafood products contaminated with neurotoxins. These toxins include several derivatives, with gonyautoxins (GTX) 1–4, saxitoxin (STX), and neosaxitoxin (NEO) being the most potent (). These neurotoxins bind to voltage-gated sodium channels, preventing the passage of sodium ions and attenuating action potentials, progressively inhibit neurotransmission, relax smooth muscle, and lead to paralysis and respiratory failure in humans, as well as mussels, fishes, and seabirds (Tester et al., 2000; ; Turner, 2014), though several studies showed no evidence of these effects on some marine zooplankton (Van Dolah, 2000; ; ; ). Shellfish, zooplankton, and herbivorous fish consume these toxic microalgae and act as vectors to humans directly or transfer through the food web to higher trophic levels (Van Dolah, 2000).
Copepods are the dominant components of mesozooplankton in the ocean, playing a vital role in marine food webs by transferring energy from primary producers to higher trophic levels (; ; ). Interestingly, studies have shown that copepods can graze on toxic microalgae without experiencing any apparent negative effects (Turner and Anderson, 1983; Uye, 1986; Turner and Tester, 1989; Turner, 2014). For instance, studies have shown that Acartia bifilosa can consume toxic Dinophysis spp. which causes diarrhetic shellfish poisoning (DSP) without any food selectivity, regardless of the concentration of the toxic microalgae (; ). Both Acartia sp. and Temora sp. do not suffer any negative effects on survival, egg production, or hatching success when consuming toxic Pseudo-nitzschia spp., which is responsible for amnesic shellfish poisoning, and Acartia sp. even had a higher ingestion rate when consuming toxic P. multiseries (). Meanwhile, the okadaic acid, a diarrhetic toxin, excreted by copepods was lower than what copepods ingested, suggesting that copepods are efficient in detoxification and do not bring toxins to upper consumers (; ). However, mechanisms for the adaptation of marine copepods to toxic microalgae, especially toxic dinoflagellates, are still unclear.
Gut microbes play a crucial role in the detoxification of toxic diets for a wide range of animals including mammals and arthropods. For instance, thanks to their gut flora, the freshwater crustacean Daphnia magna can tolerate toxic cyanobacteria and the black soldier fly (Hermetia illucens) larvae can consume toxic seaweed Rugulopteryx okamurae without experiencing any adverse effects (; ). These benefits can be attributed to the ability of gut microbes to regulate the host’s expression of detoxifying enzymes, such as cytochrome P450 and glutathione S-transferase (GST), or directly digest toxic compounds before they diffuse from the gut into other tissues of the host (; ). However, little is known about the roles of gut microbes in helping copepods adapt to toxic algal diet and maintain homeostasis (i.e., stable physiology performance) during HAB events. One study reported that the abundance of the mlrA gene, a microcystin-degrading gene, was positively correlated with the density of toxic cyanobacteria, suggesting that gut bacteria carrying mlrA gene can contribute to the biodegradation of toxins ingested by the host and facilitate nutrition (). Still, the detailed changes of GMCs, both in terms of composition and function, were much less studied and require comprehensive analysis to understand the detoxification mechanisms.
The copepod gut creates a distinct acidic and suboxic-anoxic microhabitat that promotes the colonization of microbiota, supporting unique carbon and nutrient sources (Tang et al., 2011). Bacteria are introduced into the gut when copepods feed on detritus or plankton prey that have bacteria attached to them (; ). found that large amounts of active bacteria were isolated from fecal pellets and released into the surrounding seawater, potentially playing a significant role in marine biogeochemical cycles, sedimentation, and enhancing nutrient regeneration if these bacteria have specialized functions. Many kinds of bacteria, such as Bacteroidetes, Gammaproteobacteria, and Alphaproteobacteria, which are commonly associated with phytoplankton, can enter the copepod GMC through ingested food (Teeling et al., 2012). Due to the wide distribution of copepods, such as through vertical migration and passive dispersal by ocean currents, the microbial communities within their guts contribute greatly to the distribution of bacterioplankton in all oceans.
The impact of different copepod host species and ambient environmental conditions on copepod GMC are the main focus of recent studies. Researchers have shown that differences in copepod host species played the greatest role in the composition of copepod GMC (; ; Wäge et al., 2019), mainly regulated by the different physiology and life history of the hosts (Tang, 2005; ; ). However, other reports have suggested that environmental conditions, especially food supply, may be the dominant factor regulating zooplankton GMC (Tang et al., 2011; ; ). For instance, seasonal variations may bring different food sources, leading to dramatic alterations in copepod GMC (). The significant difference in the gut content of Calanus finmarchicus among the North Atlantic regions has been attributed to diverse environmental conditions (Yeh et al., 2020). Field diets can contain various toxins and these toxic diets can lead to significant changes of GMC as a detoxification strategy of hosts, which has been reported in terrestrial arthropods (). However, the changes of copepod GMC in response to toxic dinoflagellate diet, which is vital to the understanding of their detoxifying mechanisms, has received little attention and less is known about the regulating factors, such as species difference of copepods (as host) and dinoflagellates (as food).
In this study, we hypothesize that GMCs may help the host with the detoxification of toxic dinoflagellate as a diet by changing their compositions and recruiting more beneficial bacteria. To test this hypothesis, we ask the following questions: (1) Do copepod GMC change temporally and significantly in response to food supply? Further, can GMC be different or similar between toxic and non-toxic microalgae diets? (2) Do diet-induced GMC changes differ between copepod species? (3) How do GMC functionally help copepods detoxify? To answer these questions, we used two copepod species, Acartia sp. and Paracalanus sp., in our study and fed them two dinoflagellates species from the same genus, the toxic Alexandrium tamarense (as toxic microalgae) and the non-toxic Alexandrium andersonii (as control). The two copepods and dinoflagellates are typically abundant and frequently coexist in subtropical coastal waters (; Xiang et al., 2021). In particular, the two copepod species are from different taxonomic families and have different feeding modes (Acartia sp.: ambush and filter feeding, Paracalanus sp.: filter feeding) (Turner, 1984; ; ; ), providing insights into the host differences in feeding strategies may align with differences in host GMCs. The two dinoflagellate species are from the same genus, which minimizes the influence of other factors such as nutritional values. In addition to the measurement of the physiologic performance of copepods, we revealed the copepod GMCs using the high throughput sequencing of 16S rRNA gene at the individual level and performed comprehensive analyses on the compositional and functional changes of GMC in response to the toxic diet. We showed the significant changes in GMCs of copepods under different diets which confirmed our hypothesis and identified a number of enriched bacterial taxa which may contribute to the detoxification and homeostasis of the host.
2 Materials and methods
2.1 Copepod culturing
Copepods were collected from the Port Shelter (Hong Kong SAR, China) using a 200-µm mesh size zooplankton net, and immediately taken back to the laboratory. Two copepod species, Acartia sp. and Paracalanus sp. were identified under a dissection microscope (Olympus IX51, Japan), and adults (n = 500 for each species) with similar body lengths were picked and used for experiments immediately. Experimental copepods were divided into field, toxic, and non-toxic algal diet groups. The field groups were performed with gut dissection immediately after collection, while the toxic and non-toxic algal diet groups were fed with the corresponding dinoflagellates, and cultured in triplicated 5-L beakers filled with 0.2-μm filtered local coastal seawater with continuously pumped air to provide sufficient oxygen. Algal debris and copepod carcasses were removed by droppers to maintain good water quality in culturing flasks.
2.2 Dinoflagellate culture
Alexandrium tamarense (strain TIO855, isolated from Sishili Bay in Yantai, China) and Alexandrium andersonii (strain TIO940, isolated from the coastal water of Koh Mak, Thailand) were selected as the toxic and non-toxic algal diet for copepods, respectively. They were cultured in 800-mL culture flasks containing 0.2-μm filtered sterilized seawater supplemented with f/2 medium (). Cultures were maintained in a temperature-controlled walk-in chamber at 23.5°C with a 14h:10h light: dark cycle, 60 μmol m-2 s-1 light intensity and transferred to fresh medium every two weeks (Tillmann et al., 2009). To ensure consistent results and minimize the influence of external factors, such as toxin accumulation during algae cultivation (particularly in the post-exponential phase), we provided copepods with the algae consistently at the exponential growth stage (i.e., days 10 and 7 following transfer, for A. tamarense and A. andersonii, respectively) daily throughout the experiment period.
2.3 Grazing experiments
Grazing rate depends greatly on prey concentration and can be largely biased when the prey concentration is too high. Thus, preliminary experiments were conducted to determine the minimal dinoflagellate cell concentration leading to maximal copepod grazing rates to be used in the grazing experiment. Four algal concentrations (~ 500, 1000, 1500, 2000 cells mL-1) of A. tamarense and A. andersonii were added separately to experimental bottles (500 mL) with 10 copepods and control bottles (500 mL) without copepods, both in triplicates. After incubating at room temperature in the dark and on a rotating plankton wheel (0.5 rpm) for 8h, the number of surviving copepods was recorded and 500 μL algal samples from each bottle were collected and fixed with 2% (final concentration) acid Lugol’s solution. Cell counts were obtained in 24-well cell culture plates under an inverted microscope (Olympus IX51, Japan). Dinoflagellate cells were counted from top to bottom and from left to right, counting 5 fields of view for each direction, obtaining 10 fields of view. The ingestion rate of copepods under different algal concentrations was calculated following the formula ():
where I (cells cop-1 d-1) is the ingestion rate of the copepod, Cc1 (cell mL-1) and Ce1 (cell mL-1) are the algae concentrations of the control group and the experimental group at the end of the experiment separately, Ce0 (cell mL-1) is the algae concentration of the experimental group at the beginning of the experiment, V (mL) is the volume of the algae solution for foraging experiments, N is the number of surviving copepods, and t (Day) is the duration of the grazing experiment.
Grazing experiments were performed separately by feeding copepods toxic or non-toxic dinoflagellates on day 0, 1, 2, and 3 with the minimal cell concentrations achieving maximal ingestion rate (determined in the preliminary experiment) (Supplementary Table 1). The experimental settings such as the number of copepods, grouping, temperature, and time are the same as those in the preliminary experiment. The experiments were conducted in triplicates and grazing rates were estimated as explained above.
The copepod survival rate under different algal diets was measured daily by randomly picking out 10 individuals from the culture beaker of each treatment and culturing them in a 500 mL polycarbonate (PC) bottle with the same prey concentration for 24 h. The number of alive copepods was counted and then removed. Survival rate was calculated by:
where S (%) is the survival rate of copepods, and n is the number of surviving copepods. If there was no movement within 60 s after the mechanical stimulus, the copepod was defined as dead and excluded from the calculation of survival rate ().
2.4 Algal toxin measurement
Concurrently with the grazing experiments, 200 mL of exponentially grown algae were collected and filtered using a GF/F glass fiber membrane (Whatman). Filters were added with 0.5 mL of 0.1 mol L-1 acetic acid solution, extracted by ultrasonic breaking, and centrifuged at 12,000 rpm for 5 min. The supernatant was taken to another reagent bottle, by adding 0.5 mL of extractant to the residue, and repeated the above steps to combine the supernatant liquid. Then a C-18 solid-phase extraction cartridge was used for purification. The supernatant was eluted with 0.5 mL of 0.1 mol L-1 acetic acid solution after passing through the column. The liquid and eluent passing through the column were obtained and adjusted to 2 mL with 0.1 mol L-1 acetic acid solution and then passed through a 0.22-μm filter membrane to run on the LC-MS/MS machine (LCMS-8040, Japan) for toxin detection. Here, paralytic shellfish poisoning (PSP) toxins were measured, including gonyautoxins (GTX 1-4 and GTX 5-b: formerly B1), saxitoxin (STX), neosaxitoxin (NEO), N-sulfocarbamoyl gonyautoxins (C1 and C2).
2.5 Gut dissection of copepods
For each combination of copepod species and toxic/non-toxic dinoflagellates, 15 copepod individuals were collected from the culturing beakers on day 0 and day 2. Gut dissection was performed on each individual copepod, with sterilized (i.e., autoclaved and 70% ethanol steeped) dissection tweezers (Regine 5, Switzerland) in a sterile Petri dish under the stereomicroscope (Stemi 305, Germany) (Supplementary Figure 1) and stored in a 1.5 mL tube with small amount of 70% ethanol (~ 10 μL). All samples were kept at − 80°C for DNA extraction.
2.6 DNA extraction and PCR amplification
DNA was extracted from the tubes containing copepod gut using the DNeasy Blood & Tissue Kit (QIAGEN, Germany) following manufacturer protocol. Briefly, 180 μL proteinase K was added to each 1.5-mL microcentrifuge tube containing one copepod gut. Tubes were then vortexed and incubated at 56°C until complete lysis (around 1 h, vortexed every 20 mins), followed by the addition of 200 μL Buffer AL and vortexing. Copepod gut samples were incubated at 56°C for 10 min, then 200 μL absolute ethanol was added and vortexed. Afterward, the mixture was transferred to a DNeasy Mini spin column and placed into a 2 mL collection tube. These mixtures were centrifuged at 8,000 rpm for 2 min, then the flow-through and collection tube were discarded, the spin column was placed in a new 2 mL collection tube and added 500 μL Buffer AW1. The centrifugation procedure was repeated once but added 500 μL Buffer AW2. At last, the mixture was centrifuged at 14,000 rpm for 5 min and the spin column was transferred to a new 1.5 mL microcentrifuge tube. The extracted DNA was then eluted into 60 μL AE buffer for PCR amplification.
PCR was conducted to amplify the 16S rRNA gene of bacteria in the copepod gut with the forward primer 341F (5′-CCTACGGGRSGCAGCAG-3′) () and reverse primer 787R (5′-CTACNRGGGTATCTAA-3′) (). The PCR cycling program was performed as follows: pre-denaturing at 95 °C for 5 min; 30 cycles of denaturing at 95 °C for 45 s, annealing at 55 °C for 45 s, extension at 72 °C for 60 s; and a final extension at 72 °C for 10 min. The PCR reactions were conducted in triplicates to ensure the reliability of the results. The PCR products of 16S rDNA were pooled together and sequenced by a Hiseq 2500 System (Illumina) with 2 × 250 bp paired end read configurations.
2.7 Bioinformatic analysis
R raw reads were processed using QIIME2 (https://docs.qiime2.org/) with tutorial pipelines (). After removing the barcode and primer, the sequences (with paired ends for each sample) were imported into QIIME2 (v 2019.10), visually inspected for sequence quality after demultiplexing, and processed with DADA2 to remove contamination, trim reads, correct errors, merge read pairs and remove chimeras. Representative ASV (amplicon sequence variants) sequences and their abundances were extracted by feature table (file “table.qza”). Then, the representative ASV sequences were classified with detailed taxonomy information using the trained classifier. ASVs assigned to chloroplast and mitochondria (75 of 6773 ASVs) were removed. ASV abundance of each sample was obtained, annotated to different taxonomic levels, and rarefied to the same depth using the “rrarefy” function in the “vegan” package of R software () (https://www.r-project.org/). GMCs of copepods were constructed by the ASV abundance (number of sequences) shared by samples.
After rarefied to the same depth, the α-diversity (Richness and Shannon) of GMC was calculated using the “estimate” function in the “vegan” package of R. The values of α-diversity in different treatments were compared by the Wilcoxon test after the Shapiro–Wilk normality test. Community similarity was calculated by Bray-Curtis’s distance using the “vegdist” function. To show the effects of copepod species and diet type on the GMC, grouping the GMC based on community similarity was shown by Non-metric Multidimensional Scaling (nMDS) maps using the “monoMDS” function in the “vegan” package of R, with the difference and significance tested by analysis of similarity (ANOSIM) ().
A similarity percentage (SIMPER) analysis was conducted to identify the distinct ASVs between two experimental groups (i.e., toxic vs. non-toxic) of GMC with significance (p-value) by using the “vegdist” and “anosim” functions in the “vegan” package of R. The mean relative abundance of each ASV was calculated and compared between two groups, with largely and significantly changed ones (i.e., both |fold-change| > 2 and SIMPER p-value< 0.05) regarded as distinct ASVs (between two groups).
PICRUSt2 (phylogenetic investigation of communities by reconstruction of unobserved states) () was used for the functional prediction of GMCs. In brief, sequencing data of the 16S rRNA marker gene were annotated using the phylogenetic tree from the software and matched with the microbial metagenomic database to obtain the functional information. The abundance of predicted functions was normalized by the sequence abundance of each ASV. The functional annotation results were shown in both KEGG Orthology (KO) () and Enzyme Code (EC). ANOSIM tests were performed to analyze the effects of both copepod species and algal diet types on the functions of GMC.
3 Results
3.1 Dinoflagellate toxin production and copepod ingestion and survival rate
Saxitoxin, neosaxitoxin, and N-sulfocarbamoyl gonyautoxins were not detected in the two dinoflagellates, while gonyautoxins were detected only in A. tamarense including GTX-1, GTX-2, GTX-3, and GTX-4 with concentrations of 1.10 ± 0.54, 0.03 ± 0.01, 0.01 ± 0.01, and 0.34 ± 0.17 ng per cell, respectively.
Positive correlations between dinoflagellate cell concentration and copepod ingestion rate were observed in our preliminary experiment of grazing (Supplementary Figure 2). All the ingestion rates were saturated at the highest dinoflagellate cell concentration (2000 cells mL-1). For instance, when the two copepod species fed on the toxic A.tamarense, their feeding rates reached stable at the cell concentration of 1000 cells mL-1 (Supplementary Figures 2A, C) while the saturating point occurred at ~ 2000 cells mL-1 for the non-toxic A.andersonii (Supplementary Figures 2B, D). These cell concentrations were used for the later daily measurement of the grazing rate of copepods fed on different dinoflagellate diets (Supplementary Table 1).
In the daily-measured grazing experiment, both copepod species showed a stable ingestion rate on the non-toxic diet throughout the period (7.05 ± 1.06 × 104 cells copepod-1 day-1 for Acartia sp., 7.35 ± 0.64 ×104 cells copepod-1 day-1 for Paracalanus sp.) (Figures 1B, D). However, when fed the toxic diet, the ingestion rate of both copepod species showed a slight but significant decline after the second day (p< 0.05, by Kruskal-Wallis test) (Figures 1A, C).
Figure 1
Survival rate of the two copepod species was consistently high (> 80%) during the incubation period (Figure 2). In general, survival rates of Acartia sp. were relatively lower with larger variations (88% ± 10% fed on toxic diet, 86% ± 7% fed on non-toxic diet) when compared to Paracalanus sp. which showed consistently higher values (93% ± 4% fed on toxic A. tamarense, 98% ± 1% fed on non-toxic A. andersonii). Except for Acartia sp., which had a significant increase in survival rate when eating toxic dinoflagellate in day 2, the survival rates of other groups had no significant changes.
Figure 2
3.2 Copepod GMC composition and diversity
A total of 66 GMCs were constructed by 16S rDNA sequencing data (Figure 3). In general, at the class level, Alphaproteobacteria were the most dominant component across all communities (53.80%, on average). For Acartia sp. GMC, apart from Alphaproteobacteria, Actinobacteria (17.77%, on average) was the most abundant group in the field samples, while remarkable proportions of Campylobacteria (8.32%, on average) and Saccharimonadia (6.27%, on average) were found in the non-toxic and toxic diets, respectively. For Paracalanus sp. GMC, however, no apparently consistent change of bacterial class was observed between treatments, which was mainly due to the large heterogeneity of community composition among individuals, suggesting that the difference of GMC should be compared at finer taxonomic levels (i.e., lower than class level).
Figure 3
Alpha diversity of copepod GMC was compared between different groups using the indices of species richness (i.e., observed ASV number in a sample) and Shannon diversity (Supplementary Figure 3). In general, no significant difference was found between groups by two-sample Wilcoxon tests, except the significantly higher richness observed in the toxic diet treatment for Acartia sp. compared with the non-toxic diet (Supplementary Figure 3C).
Beta diversity of copepod GMC was calculated by the Bray-Curtis distance between any two samples. While there were no significant differences between any two groups by Wilcoxon test, the average value of community dissimilarity for each group is high (> 0.8), suggesting that the GMC varied greatly among individuals and should be revealed at the individual level for copepods (Supplementary Figure 4).
3.3 Impacts of copepod species and algal diet on the composition of GMC
Effects of copepod species and algal diet type on GMC were shown and compared on the nMDS maps which were constructed based on community similarity (Bray Curtis distance) (Figure 4). A clear separation of GMC between Acartia sp. and Paracalanus sp. was observed (ANOSIM-R = 0.83, p< 0.001) (Figure 4A), suggesting the strong effects of copepod species on the structure of GMC. Significant effects of three treatments on GMC were shown for each copepod species separately (Figures 4B, C).
Figure 4
3.4 Changes of copepod GMC at ASV level in response to toxic algal diet
Compared to non-toxic dinoflagellate diet, GMC of Acartia sp. fed on toxic dinoflagellate diet had a total of 5 ASVs significantly increased (Log2FoldChange > 1, p< 0.05) and 2 ASVs significantly decreased (Log2FoldChange< -1, p< 0.05) (Figure 5A). Notably, all the increased ASVs had a large fold-change value over 4, which were from different classes including Actinobacteria, Alphaproteobacteria, Campylobacteria and Parcubacteria (Supplementary file 2).
Figure 5
For Paracalanus sp., GMC under toxic diet, compared to non-toxic, had a large number of significantly increased ASVs (n = 38) which were mainly from Alphaproteobacteria (n = 23, mainly from the order Rhodobacterales), followed by Gammaproteobacteria (n = 8) (Figure 5B). Three ASVs belonged to Actinobacteria increased significantly with a high fold-change value over 6 (Supplementary file 2). Similar to Acartia sp., only 2 ASVs belonging to Alphaproteobacteria (order Rhodobacterales) were detected as significantly decreased in the toxic dinoflagellate diet treatment for the GMC of Paracalanus sp.
3.5 Effects of copepod species and algal diet on the functions of GMC
Functions of GMC were predicted in terms of KEGG Orthology (KO) and Enzyme Code (EC). Functional similarity of GMC was calculated by Bray-Curtis distance of the two indices between two communities and visualized by nMDS maps (Figure 6). Algal diet type had a more significant effect on copepod GMC (ANOSIM-R = 0.15 for both KO and EC, p< 0.001) than copepod host species (ANOSIM-R = 0.068 for KO, ANOSIM-R = 0.048 for EC, both p< 0.05) for both functional indices.
Figure 6
4 Discussion
4.1 High survival rate of copepods under the toxic dinoflagellate diet
We observed that the ingestion rate of both copepods significantly decreased after two days of feeding on a toxic dinoflagellate diet. This finding is in accord with the reports of Acartia hudsonica feeding upon toxic Alexandrium fundyense, Euterpina acutifrons and Acartia grani feeding upon toxic Gymnodinium catenatum (; ), but different from the study of Calanus finarchicus feeding on toxic Alexandrium fundyense which showed no significant change of ingestion rate (). Despite the reduced ingestion rate, both copepod species had consistently high survival rates (> 80% in all treatments), suggesting a compensatory mechanism in response to the algal toxins. It is supported by a previous study that showed copepod Acartia tonsa only consumed a few algae during an A. tamarense bloom, to maintain their basic metabolism and survive (Xu et al., 2018). The decrease in their grazing rate may be due to the reduced activity, as well as increased prey rejection and regurgitation of ingested toxic cells such as mouthpart twitching in Calanus pacificus, a passive filter-feeder (; Sykes and Huntley, 1987; ; Xu et al., 2018). This compensatory mechanism, by reducing the feeding rate of toxic algae to ensure survival, has also been reported in many other animals including arthropods (e.g., daphnids and insects) and mammals (e.g., rats) (; Wheeler and Isman, 2001; ).
In addition to the compensatory mechanism, detoxification by copepod GMCs may also play an important role in maintaining the high survival rate of copepods under the toxic diet. This is supported by the significant increase in the relative abundance of bacterial taxa at the ASV level in the toxic diet treatment when compared to the non-toxic group. In general, there was a relatively higher number of ASVs with an increased trend compared to those showing a decreased trend, suggesting the positive recruitment of beneficial bacterial taxa rather than the removal of harmful or functionally redundant ones. This finding is further supported by our results showing no significant change in bacterial richness among treatments, except for a significant increase of GMC richness in the toxic diet treatment of Acartia sp. (compared to the non-toxic diet, Supplementary Figure 3C).
4.2 Factors determining the GMC of copepods: host vs. diet
Our results showed that host difference had a much higher impact than diet on the composition (ASV relative abundance based) of copepod GMCs (Figure 4A). Previous studies reported that copepod GMCs may be affected by various host factors, such as differences in feeding behavior, feeding status, or sex (Wäge et al., 2019). Similarly, the two copepod species used in our experiment have different feeding modes, which could contribute to the variations in their GMCs. Acartia sp. acts as a feeding-current feeder or coarse-filter feeder when feeding on small phytoplankton cells and as an active ambush-feeder when feeding on motile prey, using its mechanosensory setae on the antennules to detect the prey (; ). In contrast, Paracalanus sp. is a fine-filter feeder that uses its antennules as open rakes to filter whole parcels of water containing food particles from the water, not just single food particles on the basis of size (Turner, 1984; ; ). Therefore, Paracalanus sp. may need to filter more ambient water than Acartia sp. to obtain adequate food, which may lead to variations in the GMC. Moreover, the two copepod species have different body sizes (Acartia sp.: ~ 1.2 mm, Paracalanus sp.: ~ 0.9 mm) and biomass, which may reflect their different metabolic rates, such as excretion rate of ammonium and phosphate, leading to variations in their gut GMC (). Additionally, studies have shown that many of the bacterial sequences were found only in one of the two copepod species, suggesting host-specific affiliation of bacteria ().
For each copepod species individually, we showed that diet type had significant effects on the composition of GMC, which is consistent with previous studies on Daphnia (; ; ). also demonstrated that food sources varied due to seasonal change which brought copepod GMC differences. In our experiment, to eliminate the potential effects from other factors influencing GMC, such as different nutritional values between algal taxa, we utilized two species of dinoflagellates from the same genus and at the same growth stage. We measured their toxin productions during the experiment and showed that A. tamarense had a remarkable production of gonyautoxins which was not detected in A. andersonii. Therefore, we argue that the significant change in copepod GMC under the toxic diet compared to the non-toxic diet or field samples could be attributed to the response to the algal toxin. Similarly, showed that when exposed to toxic cyanobacteria, the gut bacteria of both tolerant and susceptible genotypes of Daphnia tended to be more diverse, with some bacteria having the ability to degrade microcystin (; ). These findings suggest that changes in GMC in response to toxic diets may reflect copepods’ adaptive strategies to cope with the toxic effects of HAB species.
4.3 How can GMC help copepods deal with a toxic algal diet?
For Acartia sp., when fed with a toxic diet, remarkably lower numbers of ASVs (n = 7) were found to show significant changes in the GMC, compared to that of Paracalanus sp. (n = 40), suggesting the host differences in the recruitment of specific bacterial taxa for detoxification. These enriched bacterial taxa were previously reported to have special or strong abilities in detoxification, especially for algal toxins. For instance, Brevibacterium (Micrococcales, Actinobacteria) species have highly active and multiple proteolytic enzymes, with the ability to completely degrade neurotoxins (e.g., neurotoxic ochratoxin A, a mycotoxin produced by fungi species) (; ; ; ).
Notably, for Paracalanus sp. fed by a toxic diet, 20 ASVs belonging to the family Rhodobacteraceae (Rhodobacterales; Alphaproteobacteria) displayed significantly increased patterns, suggesting their crucial roles in copepods’ detoxification. While most of the Rhodobacteraceae ASVs detected in our study were previously reported to have detoxifying abilities on synthetic chemicals or heavy metals (i.e., methylmercury) (; ), their roles in detoxifying algal toxins remain poorly understood and require further investigation.
In addition to Rhodobacteraceae, other bacteria taxa within Actinobacteria, Alphaproteobacteria, and Gammaproteobacteria also displayed increased patterns in the GMC of Paracalanus sp. under a toxic dinoflagellate diet. Specially, three ASVs belonging to Actinobacteria were strongly enriched (with more than 6-fold change values), together with Brevibacterium sp. that is mentioned above for Acartia sp., further highlighting the key roles of Actinobacteria in copepod detoxification of HAB algae. Within Gammaproteobacteria, the genus Pseudomonas can produce an enzyme responsible for the transformation into less toxic products (). A similar function was also reported in the genus Stenotrophomonas (Gammaproteobacteria) which can strongly degrade Aflatoxin B1 and completely degrade two microcystin variants simultaneously (; Yang et al., 2014).
Our study also revealed that the unchanged survival rate of copepods under the toxic diet could be attributed to the increases of some bacterial groups, such as Rhodobacteraceae, Rhizobiaceae, and Parcubacteria, which have functions in maintaining host homeostasis through various mechanisms (; ; ). For instance, intestine bacteria from Rhodobacteraceae play an important role in the growth and health of arthropod hosts, such as shrimps (Wang et al., 2022). Species from Rhizobiaceae can supply auxiliary as well as essential metabolites to their host (). Parcubacteria species have also been shown to play an important role in celiac disease-related pathology (). Furthermore, researchers have demonstrated that the bacteria-free Daphnia were smaller, less fecund, and had higher mortality than those with microbiota (), and copepods with bacteria consumed more algae than those without bacteria (). Similarly, microbiota has been shown to increase the host resistance to toxic environments, as evidenced by the insecticide resistance in the bean bug Riptortus pedestris and facilitated intake of plant toxins in herbivores (; ). All these studies support our finding that gut bacteria play necessary roles in maintaining the homeostasis of copepods.
We acknowledge that functional prediction with 16S rDNA data using software, such as PICRUSt2 or Tax4Fun2, was widely used in the studies of GMCs (; ; ; ; Xu et al., 2023), however, the results may suffer from the incomplete database of bacterial genomes, especially for field samples with diverse bacteria taxa (). We argue that this bias would be minor for the comparison of effects between host and diet on copepod GMC functions in our results, as the same method and database were used. We showed that GMCs, constructed by predicted functional data (i.e., Enzyme codes and KEGG Orthologs), were highly regulated by diet type compared to host effects, which is opposite to the pattern observed in the communities with ASV relative abundance data. We regard this as a functional convergence of GMC when exposed to toxic algae. This could be explained by the fact that different bacteria taxa can perform the same or similar functions. This was supported by the above-mentioned species enriched in the toxic diet treatment, most of which have the similar function of detoxification or maintaining the homeostasis of hosts. Further study combined with metatranscriptomics data will be conducted to provide a more comprehensive understanding of the mechanisms, such as gene regulation, underlying the detoxification of algal toxins by copepods, and the roles of specific bacterial groups in this process.
4.4 Copepod GMC should be studied at the individual level
GMC plays an important role in the physiology performance (e.g., fitness, growth, survival, and reproduction) of zooplankton host and greatly influences the biogeochemical cycles in the ocean (; ; ). Therefore, disentangling the major factors regulating the spatiotemporal dynamics of GMC in zooplankton remains a key issue in recent studies.
Our results found high variations in GMC among copepod individuals in the same treatment, as evident from the differences in community composition (at the class level) and beta diversity (Bray-Curtis dissimilarity, at the ASV level). This finding is in line with a previous study that reported significant variation in the GMCs of single samples of identified bacterial taxa at the order level in Temora sp. and Acartia sp., which indicated that copepods should be considered as individuals rather than a group when studying their GMC (Wäge et al., 2019).
also reported that the most abundant operational taxonomic units (OTUs), such as Flavobacteriaceae, Oceanospirillaceae, and Pseudoalteromonadaceae, were not consistently present in the gut of Calanus finmarchicus. These OTUs showed patchy distributions with relative abundances varying greatly among individual copepods. Similar observations were made in studies of other arthropod species, including Daphnia and black tiger shrimp (Penaeus monodon) (; ; ; ). Therefore, we suggest that using individual gut to do PCR amplification and analyze GMC composition is more accurate than using the bulk community of copepods.
Despite their small body size and even smaller gut compared to other animals, we successfully dissected the entire gut of each copepod and performed PCR amplification. This demonstrates the feasibility of examining copepod GMCs at the individual level. Moreover, with a high number of replicates (maximum of 15) used in the statistical analyses (e.g., ANOSIM tests), we showed that variations in GMC were relatively higher between treatments than within a certain treatment, reflecting the significant effects from either host or diet.
Overall, we provided the first study of copepod GMC at the individual level to disentangle the roles and mechanisms of gut microbes on the detoxification of copepods in response to toxic dinoflagellates. We showed that copepod species difference had a much greater impact than algal diet type on shaping copepod GMC compositions. However, from the functional aspect, the pattern was the opposite (i.e., diet type predominated). Under a toxic diet, the number and taxonomy of enriched bacteria were largely different between two copepod species, and most of these taxa were reported to have the function of degrading algal toxins (e.g., species from Brevibacterium, Pseudomonas, Stenotrophomonas and Alteromonadaceae) and maintaining the homeostasis of hosts (e.g., species from Rhodobacteraceae, Rhizobiaceae, and Parcubacteria). This indicates the potentially functional convergence of GMC under algal toxins for different copepod species. Our study highlights the key role of gut microbes in the detoxification of copepods, offering new insights into the food webs, biogeochemical cycles, and management of HABs. Further studies, such as transcriptomic analysis, will be performed to investigate the detailed gene and pathway regulation of copepod GMCs.
Statements
Data availability statement
Raw sequencing reads for analyses in this study were deposited in the online open database of the National Center for Biotechnology Information Search database (NCBI) with the projection accession number PRJNA1032266. R scripts used for analysis can be found on the website: https://github.com/JingYang00/Copepod-gut-micorbes.git.
Ethics statement
The manuscript presents research on animals that do not require ethical approval for their study.
Author contributions
JY: Data curation, Formal analysis, Investigation, Methodology, Resources, Software, Validation, Visualization, Writing – original draft, Writing – review & editing. ZX: Conceptualization, Formal analysis, Methodology, Project administration, Visualization, Writing – review & editing. YC: Data curation, Investigation, Validation, Writing – review & editing. HX: Methodology, Resources, Validation, Writing – review & editing. ZG: Investigation, Methodology, Writing – review & editing. XZ: Methodology, Resources, Writing – review & editing. MP: Methodology, Validation, Writing – review & editing. SZ: Conceptualization, Methodology, Writing – review & editing. HL: Conceptualization, Funding acquisition, Methodology, Project administration, Supervision, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was supported by the National Key Research and Development Program of China (2022YFC3105301) and the Hong Kong Research Grants Council (16101622).
Acknowledgments
The authors sincerely thank Dr. Qian Wang and Dr. Meng Chen from Xiamen University for measuring algal toxins. We would like to acknowledge Prof. Haifeng Gu for providing the dinoflagellate cultures. We are grateful to Ms. Oi Ting KWOK, Ms. Yin Ki TAM, and Dr. Zhiyuan Shi from the Department of Ocean Science, HKUST, for sharing all kinds of experimental experience, equipment, and platforms.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2024.1368315/full#supplementary-material
References
1
AkbarS.GuL.SunY.ZhouQ.ZhangL.LyuK.et al. (2020). Changes in the life history traits of Daphnia magna are associated with the gut microbiota composition shaped by diet and antibiotics. Sci. Total Environ.705, 135827. doi: 10.1016/j.scitotenv.2019.135827
2
AllenY.CalowP.BairdD. J. (1995). A mechanistic model of contaminant-induced feeding inhibition in Daphnia magna. Environ. Toxicol. Chemistry: Int. J.14, 1625–1630. doi: 10.1002/etc.5620140923
3
AlmadaA. A.TarrantA. M. (2016). Vibrio elicits targeted transcriptional responses from copepod hosts. FEMS Microbiol. Ecol.92, fiw072. doi: 10.1093/femsec/fiw072
4
AndersonD. M.HoaglandP.KaoruY.WhiteA. W. (2000). Estimated annual economic impacts from harmful algal blooms (HABs) in the United States. Massachusetts, USA: Woods Hole Oceanographic Institution. doi: 10.1575/1912/96
5
AndersonD. M.KulisD. M.SullivanJ. J.HallS. (1990). Toxin composition variations in one isolate of the dinoflagellate Alexandrium fundyense. Toxicon28, 885–893. doi: 10.1016/0041-0101(90)90018-3
6
Band-SchmidtC. J.Duran-RiverollL. M.Bustillos-GuzmanJ. J.Leyva-ValenciaI.Lopez-CortesD. J.Nunez-VazquezE. J.et al. (2019). Paralytic toxin producing dinoflagellates in Latin America: Ecology and physiology. Front. Mar. Sci.6. doi: 10.3389/fmars.2019.00042
7
Barak-GavishN.FradaM. J.KuC.LeeP. A.DiTullioG. R.MalitskyS.et al. (2018). Bacterial virulence against an oceanic bloom-forming phytoplankter is mediated by algal DMSP. Sci. Adv.4, eaau5716. doi: 10.1126/sciadv.aau5716
8
BehrendsG.SchneiderG. (1995). Impact of Aurelia aurita medusae (Cnidaria, Scyphozoa) on the standing stock and community composition of mesozooplankton in the Kiel Bight (western Baltic Sea). Mar. Ecol. Prog. Ser.127, 39–45. doi: 10.3354/meps127039
9
BickelS. L.TangK. W.GrossartH.-P. (2014). Structure and function of zooplankton-associated bacterial communities in a temperate estuary change more with time than with zooplankton species. Aquat. Microbial Ecol.72, 1–15. doi: 10.3354/ame01676
10
BolyenE.RideoutJ. R.DillonM. R.BokulichN. A.AbnetC. C.Al-GhalithG. A.et al. (2019). Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol.37, 852–857. doi: 10.1038/s41587-019-0209-9
11
Buchholz-ClevenB. E. E.RattundeB.StraubK. L. (1997). Screening for genetic diversity of isolates of anaerobic Fe (II)-oxidizing bacteria using DGGE and whole-cell hybridization. Syst. Appl. Microbiol.20, 301–309. doi: 10.1016/S0723-2020(97)80077-X
12
CallensM.MackeE.MuylaertK.BossierP.LievensB.WaudM.et al. (2016). Food availability affects the strength of mutualistic host–microbiota interactions in Daphnia magna. ISME J.10, 911–920. doi: 10.1038/ismej.2015.166
13
CallensM.WatanabeH.KatoY.MiuraJ.DecaesteckerE. (2018). Microbiota inoculum composition affects holobiont assembly and host growth in Daphnia. Microbiome6, 1–12. doi: 10.1186/s40168-018-0444-1
14
CembellaA. D.QuilliamM. A.LewisN. I.BauderA. G.Dell’AversanoC.ThomasK.et al. (2002). The toxigenic marine dinoflagellate Alexandrium tamarense as the probable cause of mortality of caged salmon in Nova Scotia. Harmful Algae1, 313–325. doi: 10.1016/S1568-9883(02)00048-3
15
ChenM.YueY.BaoX.FengX.OuZ.QiuY.et al. (2022). Effects of polystyrene nanoplastics on oxidative stress, histopathology and intestinal microbiota in largemouth bass (Micropterus salmoides). Aquac Rep.27, 101423. doi: 10.1016/j.aqrep.2022.101423
16
ClarkeK. R. (1993). Non-parametric multivariate analyses of changes in community structure. Aust. J. Ecol.18, 117–143. doi: 10.1111/j.1442-9993.1993.tb00438.x
17
ColinS. P.DamH. G. (2003). Effects of the toxic dinoflagellate Alexandrium fundyense on the copepod Acartia hudsonica: a test of the mechanisms that reduce ingestion rates. Mar. Ecol. Prog. Ser.248, 55–65. doi: 10.3354/meps248055
18
da CostaR. M.PereiraL. C. C.FernándezF. (2008). Effects of toxic Alexandrium minutum strains on the feeding and survival rates of pelagic marine copepods Acartia grani and Euterpina acutifrons. Hydrobiologia614, 55–63. doi: 10.1007/s10750-008-9536-4
19
da CostaR. M.PereiraL. C. C.FerrnándezF. (2012). Deterrent effect of Gymnodinium catenatum Graham PSP-toxins on grazing performance of marine copepods. Harmful Algae17, 75–82. doi: 10.1016/j.hal.2012.03.002
20
DattaM. S.AlmadaA. A.BaumgartnerM. F.MincerT. J.TarrantA. M.PolzM. F. (2018). Inter-individual variability in copepod microbiomes reveals bacterial networks linked to host physiology. ISME J.12, 2103–2113. doi: 10.1038/s41396-018-0182-1
21
De CorteD.LekunberriI.SintesE.GarciaJ. A. L.GonzalesS.HerndlG. J. (2014). Linkage between copepods and bacteria in the North Atlantic Ocean. Aquat. Microbial Ecol.72, 215–225. doi: 10.3354/ame01696
22
De TrochM.BoeckxP.CnuddeC.Van GansbekeD.VanreuselA.VincxM.et al. (2012). Bioconversion of fatty acids at the basis of marine food webs: insights from a compound-specific stable isotope analysis. Mar. Ecol. Prog. Ser.465, 53–67. doi: 10.3354/meps09920
23
De WitP.YamadaK.PanovaM.AndréC.JohannessonK. (2018). Diet-dependent gene expression highlights the importance of Cytochrome P450 in detoxification of algal secondary metabolites in a marine isopod. Sci. Rep.8, 16824. doi: 10.1038/s41598-018-34937-z
24
DogsM.WemheuerB.WolterL.BergenN.DanielR.SimonM.et al. (2017). Rhodobacteraceae on the marine brown alga Fucus spiralis are abundant and show physiological adaptation to an epiphytic lifestyle. Syst. Appl. Microbiol.40, 370–382. doi: 10.1016/j.syapm.2017.05.006
25
DoucetteG. J.CembellaA. D.MartinJ. L.MichaudJ.ColeT. V. N.RollandR. M. (2006). Paralytic shellfish poisoning (PSP) toxins in North Atlantic right whales Eubalaena glacialis and their zooplankton prey in the Bay of Fundy, Canada. Mar. Ecol. Prog. Ser.306, 303–313. doi: 10.3354/meps306303
26
DouglasG. M.MaffeiV. J.ZaneveldJ. R.YurgelS. N.BrownJ. R.TaylorC. M.et al. (2020). PICRUSt2 for prediction of metagenome functions. Nat. Biotechnol.38, 685–688. doi: 10.1038/s41587-020-0548-6
27
DziallasC.GrossartH.-P.TangK. W.NielsenT. G. (2013). Distinct communities of free-living and copepod-associated microorganisms along a salinity gradient in Godthåbsfjord, West Greenland. Arct Antarct Alp Res.45, 471–480. doi: 10.1657/1938-4246.45.4.471
28
DzigaD.WasylewskiM.WladykaB.NybomS.MeriluotoJ. (2013). Microbial degradation of microcystins. Chem. Res. Toxicol.26, 841–852. doi: 10.1021/tx4000045
29
ErlacherA.CernavaT.CardinaleM.SohJ.SensenC. W.GrubeM.et al. (2015). Rhizobiales as functional and endosymbiontic members in the lichen symbiosis of Lobaria pulmonaria L. Front. Microbiol.6. doi: 10.3389/fmicb.2015.00053
30
EtheridgeS. M. (2010). Paralytic shellfish poisoning: seafood safety and human health perspectives. Toxicon56, 108–122. doi: 10.1016/j.toxicon.2009.12.013
31
FernándezJ.MohedanoA. F.GayaP.MedinaM.NuñezM. (2000). Purification and properties of two intracellular aminopeptidases produced by Brevibacterium linens SR3. Int. Dairy J.10, 241–248. doi: 10.1016/S0958-6946(00)00046-7
32
FiniguerraM.AveryD. E.DamH. G. (2014). No evidence for induction or selection of mutant sodium channel expression in the copepod Acartia husdsonica challenged with the toxic dinoflagellate Alexandrium fundyense. Ecol. Evol.4, 3470–3481. doi: 10.1002/ece3.1197
33
FreeseH. M.SchinkB. (2011). Composition and stability of the microbial community inside the digestive tract of the aquatic crustacean Daphnia magna. Microb. Ecol.62, 882–894. doi: 10.1007/s00248-011-9886-8
34
FrostB. W. (1972). Effects of size and concentration of food particles on the feeding behavior of the marine planktonic copepod Calanus pacificus 1. Limnol Oceanogr17, 805–815. doi: 10.4319/lo.1972.17.6.0805
35
GorokhovaE.El-ShehawyR.LehtiniemiM.GarbarasA. (2021). How copepods can eat toxins without getting sick: Gut bacteria help zooplankton to feed in cyanobacteria blooms. Front. Microbiol.11. doi: 10.3389/fmicb.2020.589816
36
GuanS.JiC.ZhouT.LiJ.MaQ.NiuT. (2008). Aflatoxin B1 degradation by Stenotrophomonas maltophilia and other microbes selected using coumarin medium. Int. J. Mol. Sci.9, 1489–1503. doi: 10.3390/ijms9081489
37
GuillardR. R. L. (1975). Culture of phytoplankton for feeding marine invertebrates. In: SmithWLChanleyMH, editors. Culture of marine invertebrate animals: proceedings—1st conference on culture of marine invertebrate animals greenport (Boston: Springer US), 29–60. doi: 10.1007/978-1-4615-8714-9_3
38
GuisandeC.FrangópulosM.CarotenutoY.ManeiroI.RiveiroI.VergaraA. R. (2002). Fate of paralytic shellfish poisoning toxins ingested by the copepod Acartia clausi. Mar. Ecol. Prog. Ser.240, 105–115. doi: 10.3354/meps240105
39
GuoX.LvM.LiJ.DingJ.WangY.FuL.et al. (2022). The distinct toxicity effects between commercial and realistic polystyrene microplastics on microbiome and histopathology of gut in zebrafish. J. Hazard Mater434, 128874. doi: 10.1016/j.jhazmat.2022.128874
40
GyllenbergG. (1980). Feeding of the copepod Eurytemora hirundoides (Crustacea) on different algal cultures. Annales Zoologici Fennici (JSTOR)17, 181–184. Available at: https://www.jstor.org/stable/23734049.
41
HallegraeffG. M. (2003). Harmful algal blooms: a global overview. Manual harmful Mar. microalgae33, 1–22.
42
HansenB.BjornsenP. K.HansenP. J. (1994). The size ratio between planktonic predators and their prey. Limnol Oceanogr39, 395–403. doi: 10.4319/lo.1994.39.2.0395
43
HathoutA. S.AlyS. E. (2014). Biological detoxification of mycotoxins: a review. Ann. Microbiol.64, 905–919. doi: 10.1007/s13213-014-0899-7
44
HeggA.RadersmaR.UllerT. (2021). A field experiment reveals seasonal variation in the Daphnia gut microbiome. Oikos130, 2191–2201. doi: 10.1111/oik.08530
45
IkedaT.KannoY.OzakiK.ShinadaA. (2001). Metabolic rates of epipelagic marine copepods as a function of body mass and temperature. Mar. Biol.139, 587–596. doi: 10.1007/s002270100608
46
KanehisaM.GotoS. (2000). KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res.28, 27–30. doi: 10.1093/nar/28.1.27
47
KikuchiY.HayatsuM.HosokawaT.NagayamaA.TagoK.FukatsuT. (2012). Symbiont-mediated insecticide resistance. Proc. Natl. Acad. Sci.109, 8618–8622. doi: 10.1073/pnas.1200231109
48
KohlK. D.WeissR. B.CoxJ.DaleC.Denise DearingM. (2014). Gut microbes of mammalian herbivores facilitate intake of plant toxins. Ecol. Lett.17, 1238–1246. doi: 10.1111/ele.12329
49
Kozlowsky-SuzukiB.CarlssonP.RühlA.GranéliE. (2006). Food selectivity and grazing impact on toxic Dinophysis spp. by copepods feeding on natural plankton assemblages. Harmful Algae5, 57–68. doi: 10.1016/j.hal.2005.05.002
50
LandryM. R. (1980). Detection of prey by Calanus pacificus: Implications of the first antennae 1. Limnol Oceanogr25, 545–549. doi: 10.4319/lo.1980.25.3.0545
51
LiY.XuZ.LiuH. (2021). Nutrient-imbalanced conditions shift the interplay between zooplankton and gut microbiota. BMC Genomics22, 1–18. doi: 10.1186/s12864-020-07333-z
52
LinW.HeY.LiR.MuC.WangC.ShiC.et al. (2023). Adaptive changes of swimming crab (Portunus trituberculatus) associated bacteria helping host against dibutyl phthalate toxification. Environ. pollut.324, 121328. doi: 10.1016/j.envpol.2023.121328
53
LincolnJ. A.TurnerJ. T.BatesS. S.LégerC.GauthierD. A. (2001). Feeding, egg production, and egg hatching success of the copepods Acartia tonsa and Temora longicornis on diets of the toxic diatom Pseudo-nitzschia multiseries and the non-toxic diatom Pseudo-nitzschia pungens. Hydrobiologia453, 107–120. doi: 10.1023/A:1013163816771
54
LiuY.ChenZ.GaoY.ZouJ.LuS.ZhangL. (2021). Identifying the source organisms producing paralytic shellfish toxins in a subtropical bay in the South China Sea. Environ. Sci. Technol.55, 3124–3135. doi: 10.1021/acs.est.0c06991
55
LuoM.AnR.FuJ.WanS.ZhuW.WangL.et al. (2022). Comparative analysis of the gut microbiota in bighead carp under different culture patterns. J. Appl. Microbiol.132, 1357–1369. doi: 10.1111/jam.15248
56
MackeE.CallensM.De MeesterL.DecaesteckerE. (2017). Host-genotype dependent gut microbiota drives zooplankton tolerance to toxic cyanobacteria. Nat. Commun.8, 1608. doi: 10.1038/s41467-017-01714-x
57
MaszczykP.KiersztynB.GozzoS.KowalczykG.Jimenez-LamanaJ.SzpunarJ.et al. (2022). Combined effects of polystyrene nanoplastics and enrofloxacin on the life histories and gut microbiota of Daphnia magna. Water (Basel)14, 3403. doi: 10.3390/w14213403
58
MoisanderP. H.SextonA. D.DaleyM. C. (2015). Stable associations masked by temporal variability in the marine copepod microbiome. PloS One10, e0138967. doi: 10.1371/journal.pone.0138967
59
MouX.LuX.JacobJ.SunS.HeathR. (2013). Metagenomic identification of bacterioplankton taxa and pathways involved in microcystin degradation in Lake Erie. PloS One8, e61890. doi: 10.1371/journal.pone.0061890
60
MuyzerG.De WaalE. C.UitterlindenA. (1993). Profiling of complex microbial populations by denaturing gradient gel electrophoresis analysis of polymerase chain reaction-amplified genes coding for 16S rRNA. Appl. Environ. Microbiol.59, 695–700. doi: 10.1128/aem.59.3.695-700.1993
61
NoguchiT.MiyazawaK.DaigoK.ArakawaO. (2011). Paralytic shellfish poisoning (PSP) toxin-and/or tetrodotoxin-contaminated crabs and food poisoning by them. Toxin Rev.30, 91–102. doi: 10.3109/15569543.2011.640743
62
OksanenJ.KindtR.LegendreP.O’HaraB.StevensM. H. H.OksanenM. J.et al. (2007). The vegan package. Community Ecol. Package10, 719.
63
PatónD.García-GómezJ. C.LoringJ.TorresA. (2023). Composting the invasive toxic seaweed Rugulopteryx okamurae using five invertebrate species, and a mini-review on composting macroalgae. Waste Biomass Valorization14, 167–184. doi: 10.1007/s12649-022-01849-z
64
PrimecM.KlemenakM.Di GioiaD.AloisioI.CionciN. B.QuagliarielloA.et al. (2019). Clinical intervention using Bifidobacterium strains in celiac disease children reveals novel microbial modulators of TNF-α and short-chain fatty acids. Clin. Nutr.38, 1373–1381. doi: 10.1016/j.clnu.2018.06.931
65
PujalteM. J.LucenaT.RuviraM. A.ArahalD. R.MaciánM. C. (2014). The family rhodobacteraceae. In Eds. RosenbergE.DeLongE. F.LoryS.StackebrandtE.ThompsonF.. The Prokaryotes: Alphaproteobacteria and Betaproteobacteria (Berlin, Heidelberg: Springer Berlin Heidelberg), 439–512. doi: 10.1007/978-3-642-30197-1_377
66
RajarajanA.WolinskaJ.WalserJ.-C.DennisS. R.SpaakP. (2023). Host-associated bacterial communities vary between Daphnia galeata genotypes but not by host genetic distance. Microb. Ecol.85, 1578–1589. doi: 10.1007/s00248-022-02011-x
67
RegueraB.RiobóP.RodríguezF.DíazP. A.PizarroG.PazB.et al. (2014). Dinophysis toxins: causative organisms, distribution and fate in shellfish. Mar. Drugs12, 394–461. doi: 10.3390/md12010394
68
RodriguezH.ReveronI.DoriaF.CostantiniA.De Las RivasB.MunozR.et al. (2011). Degradation of ochratoxin A by Brevibacterium species. J. Agric. Food Chem.59, 10755–10760. doi: 10.1021/jf203061p
69
RoncalliV.LenzP. H.CieslakM. C.HartlineD. K. (2017). Complementary mechanisms for neurotoxin resistance in a copepod. Sci. Rep.7, 14201. doi: 10.1038/s41598-017-14545-z
70
RoncalliV.TurnerJ. T.KulisD.AndersonD. M.LenzP. H. (2016). The effect of the toxic dinoflagellate Alexandrium fundyense on the fitness of the calanoid copepod Calanus finmarchicus. Harmful Algae51, 56–66. doi: 10.1016/j.hal.2015.11.003
71
RouraÁ.DoyleS. R.NandeM.StrugnellJ. M. (2017). You are what you eat: a genomic analysis of the gut microbiome of captive and wild Octopus vulgaris paralarvae and their zooplankton prey. Front. Physiol.8. doi: 10.3389/fphys.2017.00362
72
RungrassameeW.KlanchuiA.MaibunkaewS.ChaiyapecharaS.JiravanichpaisalP.KaroonuthaisiriN. (2014). Characterization of intestinal bacteria in wild and domesticated adult black tiger shrimp (Penaeus monodon). PloS One9, e91853. doi: 10.1371/journal.pone.0091853
73
SadaiappanB.PrasannaKumarC.NambiarV. U.SubramanianM.GaunsM. U. (2021). Meta-analysis cum machine learning approaches address the structure and biogeochemical potential of marine copepod associated bacteriobiomes. Sci. Rep.11, 3312. doi: 10.1038/s41598-021-82482-z
74
SaizE.KiørboeT. (1995). Predatory and suspension feeding of the copepod Acartia tonsa in turbulent environments. Mar. Ecol. Prog. Ser.122, 147–158. doi: 10.3354/meps122147
75
ScavottoR. E.DziallasC.Bentzon-TiliaM.RiemannL.MoisanderP. H. (2015). Nitrogen-fixing bacteria associated with copepods in coastal waters of the North Atlantic Ocean. Environ. Microbiol.17, 3754–3765. doi: 10.1111/1462-2920.12777
76
SchminkeH. K. (2007). Entomology for the copepodologist. J. Plankton Res.29, i149–i162. doi: 10.1093/plankt/fbl073
77
ShoemakerK. M.MoisanderP. H. (2017). Seasonal variation in the copepod gut microbiome in the subtropical North Atlantic Ocean. Environ. Microbiol.19, 3087–3097. doi: 10.1111/1462-2920.13780
78
ShuklaS. P.BeranF. (2020). Gut microbiota degrades toxic isothiocyanates in a flea beetle pest. Mol. Ecol.29, 4692–4705. doi: 10.1111/mec.15657
79
Sison-MangusM. P.MushegianA. A.EbertD. (2015). Water fleas require microbiota for survival, growth and reproduction. ISME J.9, 59–67. doi: 10.1038/ismej.2014.116
80
SkovgaardA.Castro-MejiaJ. L.HansenL. H.NielsenD. S. (2015). Host-specific and pH-dependent microbiomes of copepods in an extensive rearing system. PloS One10, e0132516. doi: 10.1371/journal.pone.0132516
81
ŠkrinjarM.RašićJ. L.StojičićV. (1996). Lowering of ochratoxin A level in milk by yoghurt bacteria and bifidobacteria. Folia Microbiol. (Praha)41, 26–28. doi: 10.1007/BF02816335
82
SommerU.StiborH. (2002). Copepoda–Cladocera–Tunicata: the role of three major mesozooplankton groups in pelagic food webs. Ecol. Res.17, 161–174. doi: 10.1046/j.1440-1703.2002.00476.x
83
SørhaugT. (1981). Comparison of peptide hydrolases from six strains of Brevibacterium linens. Milchwissenschaft36, 137–139. Available at: https://www.cabidigitallibrary.org/doi/full/10.5555/19810469449.
84
StromplC.HoldG. L.LunsdorfH.GrahamJ.GallacherS.AbrahamW.-R.et al. (2003). Oceanicaulis alexandrii gen. nov., sp. nov., a novel stalked bacterium isolated from a culture of the dinoflagellate Alexandrium tamarense (Lebour) Balech. Int. J. Syst. Evol. Microbiol.53, 1901–1906. doi: 10.1099/ijs.0.02635-0
85
SunZ.LiuY.XuH.YanC. (2022). Genome-wide identification of P450 genes in Chironomid Propsilocerus akamusi reveals candidate genes involved in gut microbiota-mediated detoxification of Chlorpyrifos. Insects13, 765. doi: 10.3390/insects13090765
86
SykesP. F.HuntleyM. E. (1987). Acute physiological reactions of Calanus pacificus to selected dinoflagellates: direct observations. Mar. Biol.94, 19–24. doi: 10.1007/BF00392895
87
TangK. W. (2005). Copepods as microbial hotspots in the ocean: effects of host feeding activities on attached bacteria. Aquat. Microbial Ecol.38, 31–40. doi: 10.3354/ame038031
88
TangK. W.GludR. N.GludA.RysgaardS.NielsenT. G. (2011). Copepod guts as biogeochemical hotspots in the sea: Evidence from microelectrode profiling of Calanus spp. Limnol Oceanogr56, 666–672. doi: 10.4319/lo.2011.56.2.0666
89
TeelingH.FuchsB. M.BecherD.KlockowC.GardebrechtA.BennkeC. M.et al. (2012). Substrate-controlled succession of marine bacterioplankton populations induced by a phytoplankton bloom. Sci. (1979)336, 608–611. doi: 10.1126/science.1218344
90
TesterP. A.TurnerJ. T.SheaD. (2000). Vectorial transport of toxins from the dinoflagellate Gymnodinium breve through copepods to fish. J. Plankton Res.22, 47–62. doi: 10.1093/plankt/22.1.47
91
TillmannU.AlpermannT. L.da PurificaçãoR. C.KrockB.CembellaA. (2009). Intra-population clonal variability in allelochemical potency of the toxigenic dinoflagellate Alexandrium tamarense. Harmful Algae8, 759–769. doi: 10.1016/j.hal.2009.03.005
92
TurnerJ. T. (1984). Zooplankton feeding ecology: Contents of fecal pellets of the copepods Eucalanus pileatus and Paracalanus quasimodo from continental shelf waters of the Gulf of Mexico. Marine ecology progress series. Oldendorf15, 27–46. doi: 10.3354/meps015027
93
TurnerJ. T. (2014). Planktonic marine copepods and harmful algae. Harmful Algae32, 81–93. doi: 10.1016/j.hal.2013.12.001
94
TurnerJ. T.AndersonD. M. (1983). Zooplankton grazing during dinoflagellate blooms in a cape cod embayment, with observations of predation upon tin-tinnids by copepods. Mar. Ecol.4, 359–374. doi: 10.1111/j.1439-0485.1983.tb00119.x
95
TurnerJ. T.TesterP. A. (1989). “Zooplankton feeding ecology: copepod grazing during an expatriate red tide,” in Phytoplankton Blooms, eds. CosperE. M.BriceljV. M.CarpenterE. J.. (Berlin: Springer-Verlag), 359–374. doi: 10.1007/978-3-642-75280-3
96
UyeS. (1986). Impact of copepod grazing on the red-tide flagellate Chattonella antiqua. Mar. Biol.92, 35–43. doi: 10.1007/BF00392743
97
Van DolahF. M. (2000). Marine algal toxins: origins, health effects, and their increased occurrence. Environ. Health Perspect.108, 133–141. doi: 10.1289/ehp.00108s1133
98
WägeJ.StrassertJ. F. H.LandsbergerA.Loick-WildeN.SchmaleO.StawiarskiB.et al. (2019). Microcapillary sampling of Baltic Sea copepod gut microbiomes indicates high variability among individuals and the potential for methane production. FEMS Microbiol. Ecol.95, fiz024. doi: 10.1093/femsec/fiz024
99
WangR.ChenH.ZhuY.Al-MasqariZ. A.YanM.WangG.et al. (2022). Survival status of Penaeus vannamei is associated with the homeostasis and assembly process of the intestinal bacterial community. Aquaculture558, 738398. doi: 10.1016/j.aquaculture.2022.738398
100
WheelerD. A.IsmanM. B. (2001). Antifeedant and toxic activity of Trichilia americana extract against the larvae of Spodoptera litura. Entomol Exp. Appl.98, 9–16. doi: 10.1046/j.1570-7458.2001.00751.x
101
XiangC.KeZ.LiK.LiuJ.ZhouL.LianX.et al. (2021). Effects of terrestrial inputs and seawater intrusion on zooplankton community structure in Daya Bay, South China Sea. Mar. pollut. Bull.167, 112331. doi: 10.1016/j.marpolbul.2021.112331
102
XuY.ChengH.MengJ.XuB.LiX.ShenP.et al. (2023). Study on the quality and symbiotic microbial composition of Artemia nauplii in three main producing areas. J. Freshw. Ecol.38, 1–12. doi: 10.1080/02705060.2022.2158140
103
XuJ.NielsenL. T.KiørboeT. (2018). Foraging response and acclimation of ambush feeding and feeding-current feeding copepods to toxic dinoflagellates. Limnol Oceanogr63, 1449–1461. doi: 10.1002/lno.10782
104
YangF.ZhouY.YinL.ZhuG.LiangG.PuY. (2014). Microcystin-degrading activity of an indigenous bacterial strain Stenotrophomonas acidaminiphila MC-LTH2 isolated from Lake Taihu. PloS One9, e86216. doi: 10.1371/journal.pone.0086216
105
YehH. D.QuestelJ. M.MaasK. R.BucklinA. (2020). Metabarcoding analysis of regional variation in gut contents of the copepod Calanus finmarchicus in the North Atlantic Ocean. Deep Sea Res. Part II: Topical Stud. Oceanography180, 104738. doi: 10.1016/j.dsr2.2020.104738
Summary
Keywords
harmful algal blooms, dinoflagellates, copepods, gut microbes, detoxification, toxic microalgae
Citation
Yang J, Xu Z, Chen Y, Xu H, Gao Z, Zhang X, Pang M, Zhang S and Liu H (2024) Community changes of gut microbes highlight their importance in the adaptation of copepods to toxic dinoflagellates. Front. Mar. Sci. 11:1368315. doi: 10.3389/fmars.2024.1368315
Received
10 January 2024
Accepted
22 March 2024
Published
05 April 2024
Volume
11 - 2024
Edited by
Paulo Sergio Salomon, Federal University of Rio de Janeiro, Brazil
Reviewed by
Catharina Alves-de-Souza, University of Concepcion, Chile
Ajit Kumar Mohanty, Indira Gandhi Centre for Atomic Research (IGCAR), India
Updates
Copyright
© 2024 Yang, Xu, Chen, Xu, Gao, Zhang, Pang, Zhang and Liu.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhimeng Xu, zxube@connect.ust.hk; Hongbin Liu, liuhb@ust.hk
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.