Next Article in Journal
The Baget Karstic System and the Interstitial Environment of Lachein, a Hotspot of Subterranean Biodiversity in the Pyrenees (France)
Next Article in Special Issue
Finding a Pied-à-Terre: Harbour Infrastructure Facilitates the Settlement of Non-Native Corals (Tubastraea spp.) in the Southern Caribbean
Previous Article in Journal
The First Attempt to Assess the Taxonomic Diversity of the Nepticulidae of Armenia Resulted in the Discovery of New Species and Cryptic Taxa in the Caucasus
Previous Article in Special Issue
Compilation, Revision, and Annotation of DNA Barcodes of Marine Invertebrate Non-Indigenous Species (NIS) Occurring in European Coastal Regions
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Diversity
Volume 16
Issue 1
10.3390/d16010061
diversity-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Coexistence of Two Copepods, Recorded for the First Time, in NW Iberian Shelf: The Case of Oithona atlantica and the Allochthonous O. davisae

by
Lara García-Alves
*,
Andrea Ramilo
,
Santiago Pascual
,
Ángel F. González
and
Elvira Abollo
*
ECOBIOMAR Research Group, Marine Research Institute IIM-CSIC, Spanish National Research Council, Eduardo Cabello 6, 36208 Vigo, Spain
*
Authors to whom correspondence should be addressed.
Diversity 2024, 16(1), 61; https://doi.org/10.3390/d16010061
Submission received: 13 December 2023 / Revised: 10 January 2024 / Accepted: 11 January 2024 / Published: 17 January 2024
(This article belongs to the Collection Marine Invasive Species)
Browse Figures
Versions Notes

Abstract

:
This paper compiles the data regarding the first occurrence of Oithona davisae and O. atlantica in NW Spain, which is supported by morphological and molecular analysis. Additionally, we investigated the seasonal dynamics of the invasive O. davisae, revealing that its abundance is conditioned by upwelling-downwelling patterns in the Rías Baixas of Galicia. Temperature was the most correlated factor, with higher abundances in upwelling relaxation-downwelling events. More studies in long-term zooplankton dynamics and molecular analysis are needed to determine if O. davisae is displacing other native species of the same genus, such as O. atlantica, in Galician waters.

1. Introduction

Copepods are crucial components in marine planktonic food webs, functioning as the main transfer node between primary producers and fish, and represent a crucial role in nutrient recycling and export [1,2]. Due to their short generation times, copepods can respond rapidly to environmental changes and have been identified as indicators of climatic variability [3,4,5].
The Galician coastal and shelf waters are highly productive ecosystems that support important commercial activities related to marine resources exploitation and aquaculture [6]. The particular hydrodynamic conditions in the estuary of the Ría de Arosa favor successful shellfish production, mainly represented by filtering mollusks that feed on the planktonic community [7,8]. Seasonal zooplankton dynamics in Galician waters have been described on a monthly basis by a zooplankton monitoring program (RADIALES), carried out by the Spanish Institute of Oceanography (IEO) in the nearshore waters of A Coruña and the proximity of the Ria of Vigo [9,10]. However, field works in the Ria de Arosa on the planktonic assemblages are scarce and mainly focused on the phytoplankton community [8,11,12,13]. In this area, the zooplankton taxonomic composition has been described only in 1978 [14], and more recent studies have been directed to Bivalvia veliger [15] and Cirripedia larvae [16,17,18].
Three species of the genus Oithona are considered native to the NW Iberian Peninsula: O. plumifera, O. nana, and O. similis [19]. All were recorded in the monitoring program RADIALES. Previously, Corral and Álvarez-Ossorio [14] described O. nana and O. helgolandica (unaccepted status, now accepted name O. similis) in the Ria of Arosa.
Oithona atlantica (Farran, 1908) (Copepoda, Cyclopoida) is considered an epi-mesopelagic and cosmopolitan species, present in the entire Atlantic Ocean from the Arctic to the Antarctic, including the Mediterranean Sea [20]. Based on its morphological characteristics, it is closely related to O. longispina and often confused with O. plumifera.
Oithona davisae (Ferrari F.D. and Orsi, 1984) is known to be indigenous to the coastal waters of Japan and the China Sea [21]. This species, sometimes initially misidentified as O. brevicornis, has colonized several American and European estuaries and bays in the last decades. O. davisae was originally described by Ferrari and Orsi [22] in Californian estuarine waters and also reported in the Bay of Varna [23], the Wadden Sea [24], Puget Sound [25], the Marmara Sea and the Golden Horn Estuary [26], and Venice Lagoon [27,28]. This species is known as an invasive species to the Black [29] and Mediterranean Seas [30].
In Iberian Peninsula waters, Oithona davisae has been recorded in Barcelona Port, in the north-western Mediterranean [31], and in the Estuary of Bilbao in the Cantabric Sea [32]. Currently, it is considered a non-indigenous species (NIS) and has been included in the European Alien Species Information Network (EASIN). Alien species can cause negative impacts on the ecosystems they colonize, displacing native species by increasing predation, competing with them for food and habitat resources [33] or introducing allochthonous pathogens in the colonized area. Increasing the knowledge of the zooplankton communities’ taxonomic composition and their seasonal variance could be useful to address the early detection of NIS and evaluate their influence on the ecosystems.
To resolve the status of this copepod in Galician waters, we aimed to identify the Oithona species collected in the Ria of Arosa, applying morphological and molecular techniques, and to describe the seasonal variability of the allochthonous species O. davisae.

2. Materials and Methods

2.1. Sampling

Oblique tows were performed to study the presence and variability of zooplankton species close to shellfish beds in three areas of the Ría de Arosa (NW Spain): Lombos do Ulla, Corón, and Sarrido (Figure 1), covering depths from 4 to 12 m. A total of 13 surveys were conducted from November 2015 to March 2017. Due to some logistic problems, there were no samplings in December 2015, January 2016, March 2016, or February 2017. Fifteen-minute oblique tows were carried out during the day using a 72 cm diameter double bongo plankton net of 200 µm mesh size and equipped with a flowmeter. The mesh size was chosen to capture mesozooplankton species that could potentially act as intermediate hosts for Bivalvia parasites. The mesh size used only collected adults of the genus Oithona.
Samples of each bongo were preserved differently: one was preserved in ethanol (96%) for molecular analyses, and the twin was preserved in 10% formaldehyde-seawater solution for 24 h and then preserved in ethanol (70%) for morphological studies.
Folsom splitter was used to obtain representative subsamples. Counts in the subsample were transformed into individuals per cubic meter for both total and relative abundance within all individual taxa. At least 500 individuals per sample were identified under a stereomicroscope. Identification was made to the lowest possible taxonomic level. Daily temperature, salinity, and fluorescence data were obtained from a network of oceanographic stations in the same study area and at a depth up to 5 m, provided by the Technological Institute for the Marine Environment Monitoring of Galicia. Upwelling index time series was provided by IEO (www.indicedeafloramiento.ieo.es) accessed on 18 May 2021. River discharge data were available at the River Gauging Network (www.meteogalicia.es) accessed on 19 May 2021. Monthly and daily average and Pearson correlations were calculated for each variable.

2.2. Morphological Analyses

For microscopic examination and body measurements, Oithona spp. individuals preserved in formaldehyde were selected and put in a glycerin drop in depression slides. Individuals of O. davisae were identified following the descriptions of Ferrari and Orsi [22] and Temnykh and Nishida [21]. O. atlantica (Farran, 1908) was identified following the description of Nishida [34]. Females were measured under a stereomicroscope (Nikon eclipse 80i, Tokyo, Japan) attached to NIS Elements D softwareVersion 4.13 Nikon, (Tokyo, Japan), using a magnification of 20×/0.4 objective lens for O. davisae and 10×/0.25 for O. atlantica. Total length of adult females was measured from the head to the posterior end to the furcal ramus, excluding setae. Confocal images were acquired using an upright Leica Microsystems (Wetzlar, Germany) SPE laser scanning microscope. Endogenous autofluorescence was detected by exciting the samples with a 488 nm laser, and the emission signal was between 500 and 633 nm. Plan-apochromatic 10× and 20× objectives were used. All images shown correspond to a maximal intensity projection of a confocal z-series. Images were processed with Adobe Photoshop CS (California, EEUU) for Figure 2. Adult females of both species were differentiated using various diagnostic characters: (1) dorsal head profile (the rostrum is dorsally not visible in O. davisae and pointed in O. atlantica) (Figure 2A); (2) lateral rostrum appearance (Figure 2B); (3) number of outer marginal spines on the exopod of swimming legs; and (4) the size of the adult (Figure 2C).

2.3. Molecular Analyses

Genomic DNA purification was performed on individual copepods, employing QIAamp DNA Micro Kit (QIAGEN, Hilden, Germany), according to the manufacturers’ protocol for Isolation of Genomic DNA from Tissues. DNA quality and quantity were checked in a spectrophotometer Nanodrop® ND-2000 (Thermo Fisher Scientific, Waltham, MA, USA). To identify Oithona species, a fragment of large subunit ribosomal RNA (28S rRNA) gene was amplified using the primers 28SF1: 5′-GCGGAGGAAAAGAAACTAAC-3′ and 28SR1: 5′-GCATAGTTTCACCATCTTTCGGG-3′, as described by Ortman [35]. PCRs were performed in a total volume of 25 mL containing 1 µL of genomic DNA, PCR buffer at 1x concentration, 0.3 µM primers, 0.2 mM nucleotides, and 0.025 U.µL−1 DreamTaq DNA Polymerase (Thermo Scientific, Waltham, MA, USA). A negative control (no DNA) was also used. The PCR assays were carried out in a Tgradient thermocycler (Biometra, Jena, Germany), under the following reaction conditions [36]: 4 min at 95 °C, 35 cycles of 40 s at 95 °C, 40 s at 50 °C, and 1.5 min at 72 °C, followed by 7 min at 72 °C. PCR products were separated on a 2% agarose gel in Tris acetate EDTA buffer, stained with Red Safe, and scanned in a GelDoc XR documentation system (Bio-Rad Laboratories, Hercules, CA, USA). Subsequently, they were cleaned for sequencing using ExoSAP-IT™ PCR Product Cleanup (Applied Biosystems, Waltham, MA, USA) for 15 min at 37 °C, followed by inactivation for 15 min at 80 °C. Sequencing was performed with a specialized service (StabVida, Caparica, Portugal), and the chromatograms were analyzed using ChromasPro 2.1.8. (Technelysium Pty Ltd., South Brisbane, Australia). All generated sequences were searched for identity using BLAST (Basic Local Alignment Search Tool) through web servers of the National Center for Biotechnology Information (NCBI, Bethesda, MD, USA). All the Oithona 28 rDNA sequences reported in this work have been deposited in the GenBank under the accession numbers: ON114103, ON114104, and ON114106–ON114118.
Generated sequences were aligned with 20 other Oithona sequences available in GenBank: O. atlantica (JF419541–JF419543, KP033173–KP033176), O. plumifera (KP033177, KP033178), O. frigida (KU982713), O. similis (KU982702–KU982703), O. dissimilis (AB604181, AB604184), O. brevicornis (AB604189, AB604190), O. davisae (KP033164), O. nana (KP033168, KP033172), and O. simplex (AF385458). A sequence of the Cyclopoida species Paracyclopina nana (KR048803) was also included as outgroup. Alignments were performed using Clustal W [37], included in MEGA 7 [38]. jModelTest Ver. 2.1.10 software under the Akaike information criterion (A.I.C.) allowed us to infer the best DNA substitution model for the maximum likelihood (ML) tree [39]. The tree was computed with 1000 bootstrap replications. Evolutionary divergence between the 15 28S rDNA sequences obtained from Oithona species from Galicia and the other 96 sequences of Oithona species deposited in GenBank (Table 1) was also calculated by pairwise distance estimation using the maximum composite likelihood method in MEGA 7 software.

3. Results

3.1. Morphological Analyses

Herein we report, for the first time, the coexistence of two species of Copepoda Cyclopoida of the family Oithonidae, O. davisae and O. atlantica, within the study area. The average total length was 0.57 ± 0.03 mm and 1.14 ± 0.09 mm for O. davisae females (N = 30) and O. atlantica females (N = 10), respectively. Male specimens could not be measured due to low abundance; they were almost absent in the study period.

3.2. Molecular Analyses

In a first attempt, we made pools of individuals, but after reading the electropherograms of O. atlantica, which revealed ambiguous sites’ C/T, we decided to conduct individual DNA extractions. Ten partial 28S rRNA gene sequences (727–835 bp), obtained from Oithona specimens collected from Lombos do Ulla (8) and Corón (2), were similar to O. davisae from the North Sea (KP033163–KP033167) and to O. davisae from South Korea (KR048808). Identity percentages ranged from 99.13% to 100% with query cover values between 60 and 100%. Five other partial sequences (811–823 bp), obtained from Oithona specimens collected from Sarrido (2), Corón (2), and Lombos do Ulla (1), shared 97-100% nucleotide identity with O. atlantica from the Atlantic Ocean (JF419541, JF419542, JF419543), 98.95–99.83% with O. atlantica sensu lato from the Norwegian Sea (KP033173), and 97.25–98.43% with O. atlantica sensu lato from the Azores (KP033174). The query cover value was 69% in all sequences.
A maximum likelihood (GTR + G model) phylogenetic tree was constructed with 412 sites and showed that the sequences of O. davisae from Ría of Arousa were included in the same clade with O. davisae sequences from the North Sea, with bootstrap values of 100% (Figure 3). However, the species O. atlantica constituted a polyphyletic group, including in the clade sequences of O. plumifera, O. frigida, and O. similis. Our sequences of O. atlantica were placed together, forming a subclade with other O. atlantica sequences from the Atlantic Ocean and Norwegian Sea with low bootstrap values (65%), but phylogenetically distant from two sequences from Azores (KP033174-5) and one sequence from the Atlantic Ocean (JF419543). Sequences of O. atlantica from Azores were grouped in a clade with a 93% bootstrap, being close to the clade of O. similis. The sequence JF419543 of O. atlantica did not cluster in the same clade as other sequences for this species and for the same geographic area, such as JF419541 and JF419542; rather, it was closer phylogenetically to the clade that included the species O. frigida, O similis, and O. atlantica from the Azores, but with a low bootstrap value. Of the eight species of Oithona included in the phylogenetic analysis, seven appeared to be monophyletic supported with high bootstrap values, and only O. atlantica appeared to be polyphyletic.
The pairwise distances calculated between Oithona species from Galicia and other Oithona species deposited in Genbank are shown in Table 2. The analysis involved 424 positions. Intraspecific distances between the five sequences obtained of O. atlantica from Galicia ranged from 0 to 0.012, which is in accordance with those for the O. atlantica sequences deposited (0–0.019). In the study area, the intraspecific distance of O. davisae were lower than those of O. atlantica (0–0.005), being in the same range as those deposited in GenBank (0–0.002).
Thus, the species with the highest intraspecific genetic distances were O. atlantica, O. dissimilis, and O. nana. Interspecific distance confirmed that O. atlantica from Galicia is more similar to O. atlantica, O. plumifera, O. frigida, and also to O. similis, with values lower than 0.043, while the genetic distances with other species were above 0.103. Likewise, O. davisae from Galicia showed interspecific distances as low as 0–0.005 with O. davisae deposited sequences, whereas the values obtained when other Oithona species are compared were higher than 0.117.

3.3. Seasonal Abundance Variation

The analysis of seasonal variation was focused on the inner part of the Ría de Arosa (Lombos do Ulla shellfish bed), due to the presence of oceanographic stations in the same area. We examined the relationship between abiotic factors and the abundance of copepods, trying to untangle the limiting factors.
O. atlantica was detected only in October in this area (Figure 4) at very low abundances (10.39 ind.m−3). In this area, this species did not present a clear seasonal variation.
Abundance of O. davisae was mainly concentrated from May to October 2016, coinciding with higher temperatures in the study area. Its abundance peaked in June (120 ind.m−3) and September (108 ind.m−3), with an average temperature over 17 °C. On the other hand, density of O. davisae was halved from July, and its individuals disappeared in August. The abundance reached another maximum in September, and then decreased slightly in October. From that month on, the population decreased sharply, reaching minimum abundances, coinciding with low temperatures in the area (below 14 °C). In august, when O. davisae disappeared from the area, the fluorescence was minimal (0.15 mg.m−3) (Figure 5b). Salinity fluctuated greatly due to rainfall and river discharge in this inner zone of the estuary (Figure 5c–e).
The greatest presence of O. davisae occurred during the upwelling favorable season in this region (March–April to September–October) (Figure 6a). The upwelling index showed that the periods of maximum abundance of O. davisae were characterized by periods of upwelling relaxation (June, July and October) or downwelling (May and September) in days prior to sampling (Figure 6b). O. davisae did not appear in August, when the upwelling index peaked (5.30 × 103 m3 s−1 km−1). Regarding rainfall regime (Figure 6c), when the presence of this species began to be relevant in May, there were several events of intense rains, reaching maximums of 25.2 L.m−2 a few days before sampling. In September, with previous downwelling and high abundance of O. davisae, the highest rainfall reached 27.2 L.m−2.
The Pearson correlation of abundance to each parameter showed that the presence of O. davisae is significantly positively correlated with temperature (r = 0.749). Our data revealed positive correlations with the oceanographic parameters analyzed, but there was an almost null correlation with salinity (r = 0.076). In relation to fluorescence, the correlation was at the limits of significance (r = 0.442).

3.4. Expansion of the Two Copepod Species in the Ría de Arosa

Surveys in the Ría de Arosa showed the presence of O. davisae in shellfish beds of Corón (depth 13.8 ± 0.43 m), an area located in the center of the estuary, and in Sarrido (depth 2.97 ± 0.40 m), an outer area of that Ría. Abundance in these shellfish beds was lower than in the inner zone of the estuary (Lombos do Ulla) at depths of 7.15 ± 1.03 m.
Both species (O. davisae and O. atlantica) presented different spatial distributions. While O. davisae was the dominant species of Oithonidae in the innermost zone of the estuary, O. atlantica was more prevalent in the outermost zones. O. atlantica was detected in the inner part of estuary (Lombos do Ulla) only in October (10.39 ind.m−3) and in other shellfish beds beginning in May. Its greatest abundance was recorded in Corón in May (41.04 ind.m−3) with another peak of abundance in October (29.30 ind.m−3). In the outermost area (Sarrido), both species presented low abundances (<6 ind.m−3).

4. Discussion

Herein, we confirm the coexistence of invasive copepod Oithona davisae and O. atlantica in the Ría of Arosa (NW Spain), based on morphological and molecular features. Previous studies of the zooplankton community in Galician waters did not record these species [19]. The absence of studies of the mesozooplankton community in the Ria of Arosa does not allow us to define when they first appeared and whether O. atlantica is a native species. O. plumifera, morphologically a very similar species to O. atlantica, is included in the Project Checklist of marine species from Galician waters [19]. Thus, Oithona plumifera atlantica (Farran, 1913) was even considered a synonymized name [34] of O. atlantica. Additionally, considering that the Galician waters are located within the distribution range of O. atlantica, we could infer a potential confusion between O. plumifera and O. atlantica, since most of the studies in these waters were based on morphological characters, with an absence of molecular studies.
The species O. nana, O. similis, and O. plumifera, described so far in the Ría of Vigo and the waters of A Coruña (by the monitoring program RADIALES), were not detected in our area and time of study. Due to the sampling design in shallow shellfish beds, the mesh size used heavily underestimates small Oithona, and the short study period prevents us from concluding that they have been displaced by the presence of the invasive species O. davisae.
DNA-based analyses have become consolidated as part of biodiversity studies and have proven to be powerful tools for assessing the diversity of marine zooplankton species. For molecular analysis, we decide to use the large-subunit (28S) rRNA gene, because it has proven to be a useful marker to identify and discriminate among Oithona species [36,40,41] and to identify invasive species [24]. DNA extraction from Oithona species was difficult due to the small size of the individuals. Maximum likelihood analysis of the 28S rRNA gene points out that O. davisae from Galicia formed a well-supported monophyletic clade together with sequences of O. davisae from the North Sea, as described in previous studies [24,40]. Likewise, the phylogenetic tree showed that the sequences of O. atlantica obtained in this study were grouped with those of O. atlantica from the Norwegian Sea and Atlantic Ocean in a polyphyletic group, which also included O. plumifera, O. similis, and O. frigida species. Our tree revealed several well-supported clades for Oithona species but clearly did not solve the phylogeny of O. atlantica. Cornils et al. [40] described a similar phylogenetic topology, and it was supported by the interspecific genetic distances, which were significantly lower between the species of the O. atlantica group than between O. atlantica and the other Oithona species. Intraspecific genetic distances obtained for both O. davisae and O. atlantica sequences also agree with those calculated in other studies [36,40].
The existence of cosmopolitan species, living In at least two oceans, is widely accepted in marine plankton ecosystems due to the lack of apparent geographic barriers and potentially broad oceanic dispersal. In recent years, it has been shown that some of the Oithona species are circumglobally distributed and constitute a complex of cryptic species [40], demonstrating that O. similis is a complex of at least seven distinct mitochondrial lineages. Ueda et al. [41] also suggested that O. dissimilis is a species complex containing more than two cryptic species. For O. atlantica, the absence of a monophyletic clade and its mixture with other species seem to suggest that a species complex could be housed under its denomination, so it would be interesting to carry out population genetic studies for this species worldwide.
The main objective of the project that involved this zooplankton sampling was to detect which zooplankton organisms acted as intermediate hosts for the Bivalvia parasite. For this reason, the sampling design involved the use of 200 µm mesh size, necessary for the collection of the mesozooplankton targeted (0.2 to 20 mm in size). The optimal mesh size for the estimation of O. davisae would be 100 µm [42], and thus the mesh size used in this work was not completely appropriate for a correct and complete study of the life cycle of O. davisae. Nauplii, juveniles, and adult males are not retained in a 200 µm mesh size, which would explain why these smaller individuals were absent in the samples. However, we consider that its presence throughout the study period and its temporal variability are of interest since it is an invasive species that has not been described in this study area.
Regarding the temporal distribution, O. davisae was present throughout the study period (November 2015–March 2017). The highest abundances were observed in the warm season (May–October). The maximum abundance in June coincides with previous studies by O. davisae in other geographic areas such Venice Lagoon [28]. Other studies, in Sevastopol Bay [29] and Estuary of Bilbao [32], showed maximum values in autumn, and in the Aegean Sea [43], higher abundances were observed during summer months (July, September) and in winter (February). A second peak of abundance was recorded in late summer (September, 108.9 ind.m−3). In our study, the coldest season showed a quasi-absence of this species. We have to consider that we cannot compare the abundances obtained in our study with other works ([44] in southern Black Sea with maximum of 49,761 ind.m−3 at the coastal station; Pansera et al. [28] in Venice Lagoon with records of 6490 ind.m−3 at the inner station), because the mesh size used is different and our abundances are underestimated. Seawater temperature seems to be the main oceanographic parameter that influences the abundance of O. davisae (r = 0.749), with peaks of abundance in periods of higher temperatures and minimums when the temperature was less than 14 degrees. This is consistent with the thermophilic character of O. davisae. An increase in seawater temperature favors the appearance and settlement of thermophilic invasive species.
Fluorescence also appears to be related to increasing population density. During the present study, fluorescence reached two maximums, in June (9.47 mg.m−3) and in October (6.36 mg.m−3), coinciding with its peaks of abundance. These results suggest that there is a correlation between fluorescence and maximum abundance of O. davisae (r = 0.442). This coincidence seems to be related to the omnivorous feeding of this species. Previous studies [45,46,47,48] showed that Oithona species feed mainly on ciliates and dinoflagellates and that O. davisae avoids diatoms [49].
Monitored oceanographic conditions allow us to shed light on the seasonal dynamics of O. davisae. Although the greatest presence of this invasive species is concentrated in the warm season, it is necessary to reduce the time scale (monthly to daily) to untangle its dynamics in the study area, especially in species with such short life cycles and rapid growth. It was observed that higher densities during the warm period correspond with upwelling relaxation-downwelling in the Ría of Arosa coupling temperature and upwelling index data in the days prior to sampling. A large number of studies have analyzed the hydrodynamic characterization of the Galician Rías and this particular coastal upwelling system [50,51]. Additionally, previous studies have investigated the influence of upwelling-downwelling events in seasonal plankton dynamics [12,52,53,54,55]. The presence of O. davisae presented herein seems to be directly related to upwelling relaxation or downwelling events, and consequently with seawater temperature. Considering various factors (upwelling index, temperature, river discharge, and rainfall), the data suggest that retention and accumulation of seawater in the innermost area is favored during rainfall events (May, June, and September; Figure 6c). This event slows down the seawater renewal, and therefore favors warming up the surface water. These conditions of high seawater temperatures are optimal for the growth and reproduction of O. davisae [56]. Furthermore, upwelling relaxation and downwelling events favor the dominance of dinoflagellates in phytoplankton communities in Rías of Galicia [57,58]. Conversely, upwelling events favor the prevalence of diatoms. As we have already mentioned, dinoflagellates constitute the main food of O. davisae. Other species of the same genus, such as O. plumifera, are considered an upwelling-indicative species [52]. Our results suggest that O. davisae, despite being a species adapted to a wide range of temperatures [29,59,60], prefers upwelling relaxation and downwelling events. Furthermore, in recent years the average intensity of upwelling and the length of the upwelling season have been reduced [61], which may favor the settlement of these invasive species from warmer areas.
Our results showed that the presence of O. davisae in terms of abundance is limited by physical (upwelling-downwelling, seawater temperature) and biological variables such as food availability (dinoflagellates) and predation by gelatinous plankton [56] and by fish larvae [62]. According to the available literature, O. davisae is more abundant in inner areas [28,42]. In our series, we observed the expansion of this invasive species from the inner zone of the estuary (Lombos do Ulla) to other more external areas of the Ría of Arosa (Corón and Sarrido), but to a lesser extent. Consequently, the detection of O. davisae in these outermost areas seems to indicate that this species is adapting to different hydrographic conditions, according to the high tolerance of this species to variability of environmental conditions [29]. This high tolerance can become a danger to the ecosystem because this species seems to be able to establish itself successfully in a new environment [24]. Besides this, the short sampling period and the mesh size used do not make it possible to determine whether this invasive species has displaced native species of the same genus Oithona.
In this study, O. davisae and O. atlantica presented an opposite pattern. O. atlantica was abundant in more external areas with higher densities in the middle area (41.04 ind. m−3 in May and 29.30 ind.m−3 in October). The reason why the abundances are higher in the middle zone (Corón) than in the outermost zone (Sarrido) of Ría de Arosa, may be that the outermost zone where the samplings were carried out was shallow (2.97 ± 0.40 m). If we observe its geographical distribution, O. atlantica could be considered a typical species of these waters. Furthermore, other studies [63] have shown that the copepod size composition can decrease with eutrophication and in confined coastal waters [64]. This fact could also explain the greater abundance of O. davisae in the inner part of the Ría de Arosa.
In the sampling area, ballast water and bivalve transplants are probably the main routes of entry of invasive species. Organisms included in this ballast water and originating from a distant port of origin are represented by bacteria, protozoa, and invertebrates, mainly as resting stages able to survive travel conditions [65] and establish a new reproductive population in the new ecosystem. Additionally, bivalve transplants activities are also frequent, since the three sampling areas are important cockle and clam shellfish beds in Galicia. In fact, bivalve seed is imported from countries such as France or Italy, Mediterranean countries where O. davisae is settled. Moreover, the introduction of a new species entails additional implications if we assimilate the holobiont concept, which includes the host and all its symbionts. It is well known that copepods act as intermediate hosts in the biological cycle of numerous parasitic species [66,67,68,69,70,71,72,73]. The co-introduction of parasitic pathogens in wild shellfish production systems can cause unexplained mortality phenomena in populations of bivalve mollusks, with its consequent impact on native species, ecosystems, and ecosystem services [74,75,76]. Undertaking interdisciplinary studies that include the assessment of the impacts of the invasion by any particular given holobiont is a very challenging topic in socio-ecological system management.

Author Contributions

All authors contributed to the study conceptualization and design. Data collection, Á.F.G. and L.G.-A.; methodology and formal analysis, L.G.-A. and A.R.; writing—original draft, L.G.-A. and A.R.; supervision, E.A.; writing—review and editing, E.A., Á.F.G. and S.P.; funding acquisition and project administration, Á.F.G. All authors have read and agreed to the published version of the manuscript.

Funding

We thank the ECOSUMA project (PID2019-110088RB-I00) and the Marine Science programme (ThinkInAzul), supported by the Ministerio de Ciencia e Innovación and Xunta de Galicia, with funding from the European Union NextGenerationEU (PRTR-C17.I1) and European Maritime and Fisheries Fund.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Acknowledgments

We would like to thank Javier Tamame González for assisting with laboratory work and also Noelia Rivero Chaver for her assistance in collecting zooplankton samples and for helping to sort the individuals in microscopic studies. We are grateful to Lucía Sánchez Ruiloba, responsible for the microscopy service IIM, for obtaining and processing the confocal images.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Turner, J.F. The importance of small planktonic copepods and their roles in pelagic marine food webs. Zool. Stud. 2004, 43, 255–266. [Google Scholar]
  2. Valdés, V.; Carlotti, F.; Escribano, R.; Donoso, K.; Pagano, M.; Molina, V.; Fernandez, C. Nitrogen and phosphorus recycling mediated by copepods in Western Tropical South Pacific. Biogeosciences 2018, 15, 6019–6032. [Google Scholar] [CrossRef]
  3. Taylor, A.H.; Allen, J.I.; Clark, P.A. Extraction of a weak climatic signal by an ecosystem. Nature 2002, 416, 629–632. [Google Scholar] [CrossRef] [PubMed]
  4. Perry, R.I.; Batchelder, H.P.; Mackas, D.L.; Chiba, S.; Durbin, E.; Greve, W.; Verheye, H.M. Identifying global synchronies in marine zooplankton populations: Issues and opportunities. ICES J. Mar. Sci. 2004, 61, 445–456. [Google Scholar] [CrossRef]
  5. Hays, G.C.; Richardson, A.J.; Robinson, C. Climate change and marine plankton. Trends Ecol. Evol. 2005, 20, 337–344. [Google Scholar] [CrossRef]
  6. Alvarez-Salgado, X.A.; Rosón, G.; Pérez, F.F.; Figueiras, F.G.; Pazos, Y. Hydrographic variability off the Rías Baixas (NW Spain) during the upwelling season. J. Geohys. Res. 1993, 98, 14447–14455. [Google Scholar] [CrossRef]
  7. Blanton, J.O.; Tenore, K.R.; Castillejo, F.F.; Atkinson, L.P.; Schwing, F.B.; Lavin, A. The relationship of upwelling to mussel production in the rias on the western coast of Spain. J. Mar. Res. 1987, 45, 497–511. [Google Scholar] [CrossRef]
  8. Figueiras, F.G.; Pazos, Y. Microplankton assemblages in three Rías Baixas (Vigo, Arosa and Muros, Spain) with a subsurface chlorophyll maximum: Their relationships to hydrography. Mar. Ecol. Prog. Ser. 1991, 76, 219–233. Available online: http://www.jstor.org/stable/24825583 (accessed on 6 July 2021). [CrossRef]
  9. Bode, A.; Álvarez-Ossorio, M.T.; González, N.; Stephens, J.; Jordan, M.; Taylor, A.; Proctor, R. Estimations of mesozooplankton biomass in a coastal upwelling area off NW Spain. J. Plankton Res. 1998, 20, 1005–1014. [Google Scholar] [CrossRef]
  10. Bode, A.; Álvarez-Ossorio, M.T.; González, N.; Lorenzo, J.; Rodríguez, C.; Varela, M.; Varela, M.M. Seasonal variability of plankton blooms in the Ria de Ferrol (NW Spain): II. Plankton abundance, composition and biomass. Estuar. Coast. Shelf Sci. 2005, 63, 285–300. [Google Scholar] [CrossRef]
  11. Varela, M.; Penas, E. Primary production of benthic microalgae in an intertidal sand flat of the Ria de Arosa, NW Spain. Mar. Ecol. Prog. Ser. 1985, 25, 111–119. [Google Scholar] [CrossRef]
  12. Hanson, R.B.; Alvarez-Ossorio, M.T.; Cal, R.; Campos, M.J.; Roman, M.; Santiago, G.; Varela, M.; Yoder, J.A. Plankton response following a spring upwelling event in the Ria de Arosa, Spain. Mar. Ecol. Prog. Ser. 1986, 32, 101–113. [Google Scholar] [CrossRef]
  13. Barciela, R.M.; García, E.; Fernández, E. Modelling primary production in a coastal embayment affected by upwelling using dynamic ecosystem models and artificial neural networks. Ecol. Model. 1999, 120, 199–211. [Google Scholar] [CrossRef]
  14. Corral, J.; Alvarez-Ossorio, M. El zooplancton de la Ria de Arosa (NW de España): Composicion y distribucion de las comunidades en un ciclo anual. Bol. Inst. Esp. Oceanogr. 1978, 265, 133–163. [Google Scholar]
  15. Molares, J.; Fuentes, J. Recruitment of the mussel Mytilus galloprovincialis on collectors situated on the intertidal zone in the Riá de Arousa (NW Spain). Aquaculture 1995, 138, 131–137. [Google Scholar] [CrossRef]
  16. Macho, G.; Molares, J.; Vázquez, E. Timing of larval release by three barnacles from the NW Iberian Peninsula. Mar. Ecol. Prog. Ser. 2005, 298, 251–260. [Google Scholar] [CrossRef]
  17. Macho, G. Ecología Reproductiva y Larvaria del Percebe y otros Cirrípedos en Galicia. Ph.D. Thesis, Universidade de Vigo, Vigo, Spain, 2006. [Google Scholar]
  18. Macho, G.; Vázquez, E.; Giráldez, R.; Molares, J. Spatial and temporal distribution of barnacle larvae in the partially mixed estuary of the Ría de Arousa (Spain). J. Exp. Mar. Biol. Ecol. 2010, 392, 129–139. [Google Scholar] [CrossRef]
  19. Miranda, A.; Alvarez-Ossorio, M.; Valdés, L.; Casas, G. Copepoda. In Inventario de la Biodiversidad Marina de Galicia: Proyecto LEMGAL; Bañón, R., Ed.; Consellería do Mar, Xunta de Galicia: Santiago de Compostela, Spain, 2017; pp. 392–401. [Google Scholar]
  20. Mazzocchi, M.G. Oithona Baird, 1843; ICES Identification Leaflets for Plankton No. 188; ICES: Copenhagen, Denmark, 2019. [Google Scholar] [CrossRef]
  21. Temnykh, A.; Nishida, S. New record of the planktonic copepod Oithona davisae Ferrari and Orsi in the Black Sea with notes on the identity of “Oithona brevicornis”. Aquat. Invasions 2012, 7, 425–431. [Google Scholar] [CrossRef]
  22. Ferrari, F.D.; Orsi, J. Oithona Davisae, New Species, and Limnoithona sinensis (Burckhardt, 1912) (Copepoda: Oithonidae) from the Sacramento-San Joaquin Estuary, California. J. Crustac. Biol. 1984, 4, 106–126. [Google Scholar] [CrossRef]
  23. Mihneva, V.; Stefanova, K. The non-native copepod Oithona davisae (Ferrari F.D. and Orsi, 1984) in the Western Black Sea: Seasonal and annual abundance variability. BioInvasions Rec. 2013, 2, 119–124. [Google Scholar] [CrossRef]
  24. Cornils, A.; Wend-Heckmann, B. First report of the planktonic copepod Oithona davisae in the northern Wadden Sea (North Sea): Evidence for recent invasion? Helgol. Mar. Res. 2015, 69, 243–248. [Google Scholar] [CrossRef]
  25. Cordell, J.R.; Lawrence, D.J.; Ferm, N.C.; Tear, L.M.; Smith, S.S.; Herwig, R.P. Factors influencing densities of non-indigenous species in the ballast water of ships arriving at ports in Puget Sound, Washington, United States. Aquat. Conserv. Mar. Freshw. Ecosyst. 2008, 19, 322–343. [Google Scholar] [CrossRef]
  26. Isinibilir, M.; Svetlichny, L.S.; Hubareva, E.S. Small oithonid copepods in the northeastern Marmara Sea. In Proceedings of the 41st CIESM Congress, Kiel, Germany, 12–16 September 2016; Volume 41, p. 483. [Google Scholar]
  27. Vidjak, O.; Bojanić, N.; de Olazabal, A.; Benzi, M.; Brautović, I.; Camatti, E.; Hure, M.; Lipej, L.; Lučić, D.; Pansera, M.; et al. Zooplankton in Adriatic port environments: Indigenous communities and non-indigenous species. Mar. Pollut. Bull. 2018, 147, 133–149. [Google Scholar] [CrossRef] [PubMed]
  28. Pansera, M.; Camatti, E.; Schroeder, A.; Zagami, G.; Bergamasco, A. The non-indigenous Oithona davisae in a Mediterranean transitional environment: Coexistence patterns with competing species. Sci. Rep. 2021, 11, 8341. [Google Scholar] [CrossRef]
  29. Altukhov, D.A.; Gubanova, A.D.; Mukhanov, V.S. New invasive copepod Oithona davisae Ferrari and Orsi, 1984: Seasonal dynamics in Sevastopol Bay and expansion along the Black Sea coasts. Mar. Ecol. 2014, 35, 28–34. [Google Scholar] [CrossRef]
  30. Zagami, G.; Brugnano, C.; Granata, A.; Guglielmo, L.; Minutoli, R.; Aloise, A. Biogeographical distribution and ecology of the planktonic copepod Oithona davisae: Rapid invasion in Lakes Faro and Ganzirri (Central Mediterranean Sea). In Trends in Copepod Studies—Distribution, Biology and Ecology; Uteri, M., Ed.; Nova Science Publishers: New York, NY, USA, 2019; pp. 59–82. [Google Scholar]
  31. Saiz, E.; Calbet, A.; Broglio, E. Effects of small-scale turbulence on copepods: The case of Oithona davisae. Limnol. Oceanogr. 2003, 48, 1304–1311. [Google Scholar] [CrossRef]
  32. Uriarte, I.; Villate, F.; Iriarte, A. Zooplankton recolonization of the inner estuary of Bilbao: Influence of pollution abatement, climate and non-indigenous species. J. Plankton Res. 2015, 38, 718–731. [Google Scholar] [CrossRef]
  33. Cinar, M.E.; Bilecenoglu, M.; Ozturk, B.; Katagan, T.; Aysel, V. Alien species on the coasts of Turkey. Mediterr. Mar. Sci. 2005, 6, 119. [Google Scholar] [CrossRef]
  34. Nishida, S. Taxonomy and distribution of the family Oithonidae (Copepoda, Cyclopoida) in the Pacific and Indian Oceans. Bull. Ocean. Res. Inst. Univ. Tokyo 1985, 20, 167. [Google Scholar]
  35. Ortman, B.D. DNA Barcoding the Medusozoa and Ctenophora. Ph.D. Thesis, University of Connecticut, Storrs, CT, USA, 2008. [Google Scholar]
  36. Cepeda, G.D.; Blanco-Bercial, L.; Bucklin, A.; Berón, C.M.; Viñas, M.D. Molecular Systematic of Three Species of Oithona (Copepoda, Cyclopoida) from the Atlantic Ocean: Comparative Analysis Using 28S rDNA. PLoS ONE 2012, 7, e35861. [Google Scholar] [CrossRef]
  37. Thompson, J.D.; Higgins, D.G.; Gibson, T.J. CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22, 4673–4680. [Google Scholar] [CrossRef] [PubMed]
  38. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [PubMed]
  39. Darriba, D.; Taboada, G.L.; Doallo, R.; Posada, D. jModelTest 2: More models, new heuristics and parallel computing. Nat. Methods 2012, 9, 772. [Google Scholar] [CrossRef]
  40. Cornils, A.; Wend-Heckmann, B.; Held, C. Global phylogeography of Oithona similis s.l. (Crustacea, Copepoda, Oithonidae)—A cosmopolitan plankton species or a complex of cryptic lineages? Mol. Phylogenetics Evol. 2017, 107, 473–485. [Google Scholar] [CrossRef] [PubMed]
  41. Ueda, H.; Yamaguchi, A.; Saitoh, S.-I.; Sakaguchi, S.O.; Tachihara, K. Speciation of two salinity-associated size forms of Oithona dissimilis (Copepoda: Cyclopoida) in estuaries. J. Nat. Hist. 2011, 45, 2069–2079. [Google Scholar] [CrossRef]
  42. Pansera, M.; Granata, A.; Guglielmo, L.; Minutoli, R.; Zagami, G.; Brugnano, C. How does mesh-size selection reshape the description of zooplankton community structure in coastal lakes? Estuar. Coast. Shelf Sci. 2014, 151, 221–235. [Google Scholar] [CrossRef]
  43. Kurt, T.T.; Beşiktepe, S. First distribution record of the invasive copepod Oithona davisae Ferrari and Orsi, 1984, in the coastal waters of the Aegean Sea. Mar. Ecol. 2019, 40, e12548. [Google Scholar] [CrossRef]
  44. Yıldız, I.; Feyzioglu, A.M.; Besiktepe, S. First observation and seasonal dynamics of the new invasive planktonic copepod Oithona davisae Ferrari and Orsi, 1984 along the southern Black Sea (Anatolian Coast). J. Nat. Hist. 2017, 51, 127–139. [Google Scholar] [CrossRef]
  45. Atkinson, A. Omnivory and feeding selectivity in five copepod species during spring in the Bellingshausen Sea, Antarctica. ICES J. Mar. Sci. 1995, 52, 385–396. [Google Scholar] [CrossRef]
  46. Nakamura, Y.; Turner, J.T. Predation and respiration by the small cyclopoid copepod Oithona similisr: How important is feeding on ciliates and heterotrophic flagellates? J. Plankton Res. 1997, 19, 1275–1288. [Google Scholar] [CrossRef]
  47. Broglio, E.; Saiz, E.; Calbet, A.; Trepat, I.; Alcaraz, M. Trophic impact and prey selection by crustacean zooplankton on the microbial communities of an oligotrophic coastal area (NW Mediterranean Sea). Aquat. Microb. Ecol. 2004, 35, 65–78. [Google Scholar] [CrossRef]
  48. Saiz, E.; Griffell, K.; Calbet, A.; Isari, S. Feeding rates and prey: Predator size ratios of the nauplii and adult females of the marine cyclopoid copepod Oithona davisae. Limnol. Oceanogr. 2014, 59, 2077–2088. [Google Scholar] [CrossRef]
  49. Uchima, M.; Hirano, R. Food of Oithona davisae (Copepoda: Cyclopoida) and the effect of food concentration at first feeding on the larval growth. Bull. Plankton Soc. Jpn. 1986, 33, 21–28. [Google Scholar]
  50. Alvarez-Salgado, X.A.; Rosón, G.; Pérez, F.F.; Figueiras, F.G.; Pazos, Y. Nitrogen cycling in an estuarine upwelling system, the Ría de Arousa (NW Spain). I. Short-time-scale patterns of hydrodynamic and biogeochemical circulation. Mar. Ecol. Prog. Ser. 1996, 135, 259–273. [Google Scholar] [CrossRef]
  51. Villacieros-Robineau, N.; Herrera, J.L.; Castro, C.G.; Piedracoba, S.; Roson, G. Hydrodynamic characterization of the bottom boundary layer in a coastal upwelling system (Ría de Vigo, NW Spain). Cont. Shelf Res. 2013, 68, 67–79. [Google Scholar] [CrossRef]
  52. Blanco-Bercial, L.; Álvarez-Marqués, F.; Cabal, J.A. Changes in the mesozooplankton community associated with the hidrographic off the northwestern Iberian Peninsula. ICES J. Mar. Sci. 2006, 63, 799–810. [Google Scholar] [CrossRef]
  53. Bode, A.; Alvarez-Ossorio, M.T.; Cabanas, J.M.; Miranda, A.; Varela, M. Recent trends in plankton and upwelling intensity off Galicia (NW Spain). Prog. Oceanogr. 2009, 83, 342–350. [Google Scholar] [CrossRef]
  54. Buttay, L.; Miranda, A.; Casas, G.; González-Quirós, R.; Nogueira, E. Long-term and seasonal zooplankton dynamics in the northwest Iberian shelf and its relationship with meteo-climatic and hydrographic variability. J. Plankton Res. 2015, 38, 106–121. [Google Scholar] [CrossRef]
  55. Olmos, L. Ecology of Cephalopod Paralarvae in a Seasonal Upwelling System. Ph.D. Thesis, Universidade de Vigo, Vigo, Spain, 2018. [Google Scholar]
  56. Uye, S.; Sano, K. Seasonal reproductive biology of the small cyclopoid copepod Oithona davisae in a temperate eutrophic inlet. Mar. Ecol. Prog. Ser. 1995, 118, 121–128. [Google Scholar] [CrossRef]
  57. Fraga, S.; Anderson, D.M.; Bravo, I.; Reguera, B.; Steidinger, K.A.; Yentsch, C.M. Influence of upwelling relaxation on dinoflagellates and shellfish toxicity in Ria de Vigo, Spain. Estuar. Coast. Shelf Sci. 1988, 27, 349–361. [Google Scholar] [CrossRef]
  58. Crespo, B.G.; Figueiras, F.G.; Porras, P.; Teixeira, I. Downwelling and dominance of autochthonous dinoflagellates in the NW Iberian margin: The example of the Ría de Vigo. Harmful Algae 2006, 5, 770–781. [Google Scholar] [CrossRef]
  59. Uye, S.; Sano, K. Seasonal variations in biomass, growth rate and production rate of the small cyclopoid copepod Oithona davisae in a temperate eutrophic inlet. Estuar. Coast. Shelf Sci. 1998, 163, 37–44. [Google Scholar]
  60. Martens, P.; Beusekom, J.E.E. Zooplankton response to a warmer northern Wadden Sea. Helgol. Mar. Res. 2008, 62, 67–75. [Google Scholar] [CrossRef]
  61. Álvarez-Salgado, X.A.; Labarta, U.; Fernández-Reiriz, M.J.; Figueiras, F.; Rosón, G.; Piedracoba, S.; Filgueira, R.; Cabanas, J.M. Renewal time and the impact of harmful algal blooms on the extensive mussel raft culture of the Iberian coastal upwelling system (SW Europe). Harmful Algae 2008, 7, 849–855. [Google Scholar] [CrossRef]
  62. Sanchez-Velasco, L. Diet composition and feeding habits of fish larvae of two co-occurring species (Pisces: Callionymidae and Bothidae) in the North-western Mediterranean. ICES J. Mar. Sci. 1998, 55, 299–308. [Google Scholar] [CrossRef]
  63. Uye, S.-I. Replacement of large copepods by small ones with eutrophication of embayments: Cause and consequence. Hydrobiologia 1994, 292, 513–519. [Google Scholar] [CrossRef]
  64. Belmonte, G.; Cavallo, A. Body size and its variability in the copepod Acartia margalefi (Calanoida) from Lake Acquatina (SE Italy). Ital. J. Zoól. 1997, 64, 377–382. [Google Scholar] [CrossRef]
  65. Belmonte, G.; Rubino, F. Cysts and Resting Eggs from Marine Zooplankton: Dimension of the Phenomenon, Physiology of Rest, and Ecological and Biogeographic Implications. In Dormancy in Aquatic Organisms. Theory, Human Use and Modeling; Editors Taylor & Francis: Abingdon, UK, 2019; pp. 71–94. [Google Scholar] [CrossRef]
  66. Arzul, I.; Chollet, B.; Boyer, S.; Bonnet, D.; Gaillard, J.; Baldi, Y.; Robert, M.; Joly, J.P.; Garcia, C.; Bouchoucha, M. Contribution to the understanding of the cycle of the protozoan parasite Marteilia refringens. Parasitology 2014, 141, 227–240. [Google Scholar] [CrossRef]
  67. Carrasco, N.; López-Flores, I.; Alcaraz, M.; Furones, M.D.; Berthe, F.C.J.; Arzul, I. Dynamics of the parasite Marteilia refringens (Paramyxea) in Mytilus galloprovincialis and zooplankton populations in Alfacs Bay (Catalonia, Spain). Parasitology 2007, 134, 1541–1550. [Google Scholar] [CrossRef]
  68. Foster, R.; Peeler, E.; Bojko, J.; Clark, P.F.; Morritt, D.; Roy, H.E.; Stebbing, P.; Tidbury, H.J.; Wood, L.E.; Bass, D. Pathogens co-transported with invasive non-native aquatic species: Implications for risk analysis and legislation. NeoBiota 2021, 69, 79–102. [Google Scholar] [CrossRef]
  69. Gregori, M.; Aznar, F.J.; Abollo, E.; Roura, Á.; González, Á.F.; Pascual, S. Nyctiphanes couchii as intermediate host for the acanthocephalan Bolbosoma balaenae in temperate waters of the NE Atlantic. Dis. Aquat. Org. 2012, 99, 37–47. [Google Scholar] [CrossRef] [PubMed]
  70. Gregori, M.; Aznar, F.J.; Abollo, E.; Roura, A.; González, A.F.; Pascual, S. Nyctiphanes couchii as intermediate host for Rhadinorhynchus sp. (Acanthocephala, Echinorhynchidae) from NW Iberian Peninsula waters. Dis. Aquat. Org. 2013, 105, 9–20. [Google Scholar] [CrossRef] [PubMed]
  71. Gregori, M.; Roura, Á.; Abollo, E.; González, Á.F.; Pascual, S. Anisakis simplex complex (Nematoda: Anisakidae) in zooplankton communities from temperate NE Atlantic waters. J. Nat. Hist. 2015, 49, 755–773. [Google Scholar] [CrossRef]
  72. Gregori, M.; Fernández-Leborans, G.; Roura, A.; Abollo, E.; González, Á.F.; Pascual, S. Description of a new epibiotic relationship (Suctorian-Copepoda) in NE Atlantic waters: From morphological to phylogenetic analyses. Acta Zoól. 2016, 97, 165–176. [Google Scholar] [CrossRef]
  73. Peeler, E.J.; Oidtmann, B.C.; Midtlyng, P.J.; Miossec, L.; Gozlan, R.E. Non-native aquatic animals introductions have driven disease emergence in Europe. Biol. Invasions 2011, 13, 1291–1303. [Google Scholar] [CrossRef]
  74. Carnegie, R.B.; Arzul, I.; Bushek, D. Managing marine mollusc diseases in the context of regional and international commerce: Policy issues and emerging concerns. Philos. Trans. R. Soc. B Biol. Sci. 2016, 371, 20150215. [Google Scholar] [CrossRef]
  75. Roy, H.E.; Hesketh, H.; Purse, B.V.; Eilenberg, J.; Santini, A.; Scalera, R.; Stentiford, G.D.; Adriaens, T.; Bacela-Spychalska, K.; Bass, D.; et al. Alien Pathogens on the Horizon: Opportunities for Predicting their Threat to Wildlife. Conserv. Lett. 2016, 10, 477–484. [Google Scholar] [CrossRef]
  76. Villalba, A.; Iglesias, D.; Ramilo, A.; Darriba, S.; Parada, J.B.; No, E.; Abollo, E.; Molares, J.; Carballal, M.J. Cockle Cerastoderma edule fishery collapse in the Ría de Arousa (Galicia, NW Spain) associated with the protistan parasite Marteilia cochillia. Dis. Aquat. Org. 2014, 109, 55–80. [Google Scholar] [CrossRef]
Diversity 16 00061 g001
Figure 1. Study area. Shellfish beds: 1: Sarrido; 2: Lombos do Ulla; 3: Corón.
Figure 1. Study area. Shellfish beds: 1: Sarrido; 2: Lombos do Ulla; 3: Corón.
Diversity 16 00061 g001
Diversity 16 00061 g002
Figure 2. Confocal images: (A) dorsal view, (B) rostrum, (C) lateral view. Scale bar = 100 µm.
Figure 2. Confocal images: (A) dorsal view, (B) rostrum, (C) lateral view. Scale bar = 100 µm.
Diversity 16 00061 g002
Diversity 16 00061 g003
Figure 3. Maximum likelihood (GTR + G model) analysis of 28S rRNA gene showing the taxonomic position of Oithona species from Ría de Arousa. Numbers at branch nodes indicate bootstrap confidence values in percent. Sequences obtained in this study: ON114103–ON114104, ON114106–ON114118.
Figure 3. Maximum likelihood (GTR + G model) analysis of 28S rRNA gene showing the taxonomic position of Oithona species from Ría de Arousa. Numbers at branch nodes indicate bootstrap confidence values in percent. Sequences obtained in this study: ON114103–ON114104, ON114106–ON114118.
Diversity 16 00061 g003
Diversity 16 00061 g004
Figure 4. Seasonal variations of O. davisae abundance and water temperature.
Figure 4. Seasonal variations of O. davisae abundance and water temperature.
Diversity 16 00061 g004
Diversity 16 00061 g005
Figure 5. Oceanographic conditions in the study area: (a) temperature, (b) fluorescence, (c) salinity, (d) rainfall, (e) river discharge.
Figure 5. Oceanographic conditions in the study area: (a) temperature, (b) fluorescence, (c) salinity, (d) rainfall, (e) river discharge.
Diversity 16 00061 g005
Diversity 16 00061 g006
Figure 6. Upwelling index: (a) interannual variability, (b) UI monthly variation (upwelling index) during the period of maximum abundance and average daily temperature (temperature records), (c) rainfall.
Figure 6. Upwelling index: (a) interannual variability, (b) UI monthly variation (upwelling index) during the period of maximum abundance and average daily temperature (temperature records), (c) rainfall.
Diversity 16 00061 g006
Table 1. 28S rDNA sequences of Oithona species from Galicia and others deposited in GenBank used to estimate the genetic divergence. N = number of sequences.
Table 1. 28S rDNA sequences of Oithona species from Galicia and others deposited in GenBank used to estimate the genetic divergence. N = number of sequences.
SpeciesNAccession Number
O. atlantica12KP033173–KP033176; JF419541–JF419543; ON114103; ON114104; ON114106–ON114108
O. plumifera2KP033177–KP0331778
O. frigida1KU982713
O. similis55JF419529–JF419540; FM991723; EF460779; KF153700; KX495626; KP033160–KP033162; KU982646–KU982648; KU982650–KU982651; KU982655; KU982657; KU982659-KU982660; KU982663; KU982665; KU982671–KU982674; KU982676; KU982678–KU982686; KU982702–KU982703; KU982694; KU982699; KU982705–KU982711
O. dissimilis12AB604177–AB604188
O. brevicornis2AB604189–AB604190
O. davisae15KP033163–KP033167; ON114109–ON114118
O. nana11FM991725–FM991727; JF419544; JF419546–JF419547; KP033168–KP033172
O. simplex1AF385458
Table 2. Genetic divergence between 28S rDNA sequences of Oithona atlantica/O. davisae from Galicia (Gal) and sequences of other Oithona species deposited in GenBank. Intraspecific distances are also shown.
Table 2. Genetic divergence between 28S rDNA sequences of Oithona atlantica/O. davisae from Galicia (Gal) and sequences of other Oithona species deposited in GenBank. Intraspecific distances are also shown.
Intraspecific DistancesInterspecific Distances
O. atlantica GalO. davisae Gal
O. davisae Gal0–0.005--
O. davisae0–0.0020.120–0.1300–0.005
O. atlantica Gal0–0.012-0.120–0.133
O. atlantica0–0.0190–0.0220.117–0.126
O. plumifera0.0000.014–0.0270.126–0.129
O. frigida-0.012–0.0240.120–0.146
O. similis0–0.0100.025–0.0430.138–0.146
O. dissimilis0–0.0200.103–0.1250.134–0.145
O. brevicornis0.0020.130–0.1470.156–0.168
O. nana0–0.0270.180–0.2200.195–0.230
O. simplex-0.130–0.1470.156–0.163
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

García-Alves, L.; Ramilo, A.; Pascual, S.; González, Á.F.; Abollo, E. Coexistence of Two Copepods, Recorded for the First Time, in NW Iberian Shelf: The Case of Oithona atlantica and the Allochthonous O. davisae. Diversity 2024, 16, 61. https://doi.org/10.3390/d16010061

AMA Style

García-Alves L, Ramilo A, Pascual S, González ÁF, Abollo E. Coexistence of Two Copepods, Recorded for the First Time, in NW Iberian Shelf: The Case of Oithona atlantica and the Allochthonous O. davisae. Diversity. 2024; 16(1):61. https://doi.org/10.3390/d16010061

Chicago/Turabian Style

García-Alves, Lara, Andrea Ramilo, Santiago Pascual, Ángel F. González, and Elvira Abollo. 2024. "Coexistence of Two Copepods, Recorded for the First Time, in NW Iberian Shelf: The Case of Oithona atlantica and the Allochthonous O. davisae" Diversity 16, no. 1: 61. https://doi.org/10.3390/d16010061

APA Style

García-Alves, L., Ramilo, A., Pascual, S., González, Á. F., & Abollo, E. (2024). Coexistence of Two Copepods, Recorded for the First Time, in NW Iberian Shelf: The Case of Oithona atlantica and the Allochthonous O. davisae. Diversity, 16(1), 61. https://doi.org/10.3390/d16010061

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop