Domoic Acid and Pseudo-nitzschia spp. Connected to Coastal Upwelling along Coastal Inhambane Province, Mozambique: A New Area of Concern

Abstract

Harmful algal blooms (HABs) are increasing globally in frequency, persistence, and geographic extent, posing a threat to ecosystem and human health. To date, no occurrences of marine phycotoxins have been recorded in Mozambique, which may be due to absence of a monitoring program and general awareness of potential threats. This study is the first documentation of neurotoxin, domoic acid (DA), produced by the diatom Pseudo-nitzschia along the east coast of Africa. Coastal Inhambane Province is a biodiversity hotspot where year-round Rhincodon typus (whale shark) sightings are among the highest globally and support an emerging ecotourism industry. Links between primary productivity and biodiversity in this area have not previously been considered or reported. During a pilot study, from January 2017 to April 2018, DA was identified year-round, peaking during Austral winter. During an intense study between May and August 2018, our research focused on identifying environmental factors influencing coastal productivity and DA concentration. Phytoplankton assemblage was diatom-dominated, with high abundances of Pseudo-nitzschia spp. Data suggest the system was influenced by nutrient pulses resulting from coastal upwelling. Continued and comprehensive monitoring along southern Mozambique would provide critical information to assess ecosystem and human health threats from marine toxins under challenges posed by global change.

1. Introduction

Wind-driven upwelling is a key driver of physical, biogeochemical, and ecological variability within coastal ocean environments, as nutrients delivered by coastal upwelling stimulate growth of phytoplankton, ultimately fueling diverse and productive marine ecosystems. The dynamics of wind-driven upwelling are fundamentally driven by large-scale atmospheric processes and have been extensively studied along the Eastern Boundary System (EBS) such as the Canary, California, Benguela, and Humbolt regions [1]. The Western Boundary System (WBS) such as the Gulf Stream, Kuroshio, and Agulhas regions are generally oligotrophic compared to EBS [2], and are known for mesoscale eddy activity and coastal upwelling [3,4,5]. Mesoscale eddies that upwell nutrients from deeper waters play an important role, increasing diversity within the food web due to the pulsed effect of introducing nutrients to the system [3,6,7].

Harmful algal blooms (HABs) are a natural phenomenon, which are increasing globally in frequency, persistence, and geographic extent [8,9,10,11,12]. These blooms can cause substantial harm to coastal societies through health impacts, increased costs associated with bloom management, and decreased revenue during fisheries and beach closures [13,14,15]. Domoic acid (DA), a toxin produced by some Pseudo-nitzschia spp., has been widely monitored since the first recorded outbreak resulted in the death of at least three people following the consumption of contaminated shellfish [16]. Domoic acid is a potent neurotoxin produced by several species of Pseudo-nitzschia, leading to both acute and chronic toxicity in mammals, fish, and birds [17,18,19,20,21]. Human consumption of highly contaminated shellfish can lead to Amnesic Shellfish Poisoning (ASP); symptoms include gastrointestinal distress, confusion, disorientation, seizures, short-term memory loss, and, in rare cases, death [16,22,23]. No human fatalities have been reported since the widespread implementation of monitoring programs [24,25,26].

Pseudo-nitzschia spp. is a cosmopolitan diatom, associated with upwelling regions, particularly the EBS, including the California, Humboldt, Canary, and Benguela currents [27,28], due to its ability to take advantage of disturbance or transitional ecological states. Blooms of Pseudo-nitzschia in EBS typically occur between late spring and summer when seasonal upwelling introduces key nutrients, particularly nitrogen [29,30,31]. Blooms of Pseudo-nitzschia along the WBS typically occur in shallow bays or estuaries where nutrients are terrestrially sourced, such as along the coast of China and eastern United States [32,33,34]. Incidences of DA have been less prevalent within the WBS than the EBS [24]. While the EBS has received attention for an increase in Pseudo-nitzschia abundance and DA events, coastal zones worldwide are being impacted by DA due to increased frequency of warm ocean anomalies [35], and the WBS could be a potential area of concern for future toxic blooms. In the southern hemisphere, DA emerged as a potential threat to fisheries within the WBS in New Zealand in the 1990s [36,37]. Along the east coast of South America, a bloom of Pseudo-nitzschia australis produced detectable levels of DA within the Argentine Sea in 2000 [38,39].

HABs within the Southwestern Indian Ocean (SWIO) are understudied due to a lack of monitoring infrastructure, despite the presence of potentially harmful species and documented marine toxins within the region [40]. Okadaic acid and ciguatoxin, produced by dinoflagellates, have been recorded along the east coast of Madagascar and nearby islands including Mauritius, the French Islands, Comoros, and Seychelles [40,41,42]. The consumption of a single shark contaminated with ciguatoxin along the coast of Madagascar led to the hospitalization of 188 people and the death of 50 people in 1993 [41]. South Africa has the most developed HAB monitoring program within the SWIO and has reported multiple derivatives of okadaic acid and yessotoxins within mussels along the western side of the country [43]. While Pseudo-nitzschia multiseries has been responsible for DA accumulation in mussels within Algoa Bay, South Africa [44], no marine toxins have caused shellfish closures along the eastern coast [45]. To date, there have been no reported HAB-related events in Mozambique [40], though Pseudo-nitzschia is among the coastal diatom assemblage [46].

The plankton communities of Delagoa Bight in Maputo Bay and Sofala Bank along Mozambique’s coastline [46,47] have high primary productivity that support higher trophic level diversity along southern Mozambique [48,49,50]. Among Mozambique’s coastal regions, Inhambane Province, located between Delagoa Bight and Sofala Bank, is known as a biodiversity hotspot [51,52], with whale sharks, manta rays, migrating humpback whales, bottlenose dolphins, humpback dolphins, dugongs, sea turtles, small-eye stingray, guitar sharks, corals, and many reef fishes utilizing the region [51,52,53,54,55,56]. This biodiversity supports a prominent tourism industry for SCUBA diving, snorkeling, surfing, and fishing [57]. In addition, roughly two-thirds of Mozambique’s population lives along the coast and relies on marine resources as a source of food and income [58,59].

Currently, there are limited studies contributing to the understanding of phytoplankton dynamics in coastal Mozambique [46,60,61,62,63]. Here we present the first known examination of the phytoplankton community composition within the coastal system of Inhambane Province and document the occurrence of DA from January 2017 to August 2018. We examine plankton abundance in relationship to inorganic dissolved nutrients (nitrogen (N) phosphorus (P) and silica (Si)) and other environmental parameters at four sites along the coast of Inhambane Province. We intend for this research to serve as a baseline to identify future changes to the coastal area, provide essential baseline data of the phytoplankton assemblage with the region, and initiate local monitoring of DA and other potential phycotoxins in Inhambane Province, Mozambique, to help ensure human and ecosystem health.

2. Results

2.1. Sites

During pilot DA monitoring from January 2017 to April 2018, all samples were collected at the soft substrate region. Unless otherwise noted, results refer to the intense study that occurred between May and August 2018. No significant difference in temperature (T; p = 0.81), dissolved oxygen (DO; p = 0.76), or salinity (p = 0.75) was observed between sites. Dissolved inorganic nutrient concentrations and particulate DA had no significant differences in concentrations among sites (N: p = 0.573, P: p = 0.300, Si: p = 0.479), which suggest that the coastal system from coastal Barra Beach to Praia de Jangamo (coastal Inhambane Province) is well-mixed and can be considered as a single region (Figure 1). Data from the four sites were combined into one region.

2.2. Upwelling Index

Localized upwelling index was used to estimate upwelling transport flux (m³ s⁻¹ km⁻¹) along the coast of Inhambane Province (Figure 2). Fluctuations lasting a few days to a few weeks occurred through May into June 2018. Prolonged increased upwelling flux occurred following the 14 July 2018 and peaked at the end of August (following the study period). This event (approximately the 14th of July through August) hereafter is referred to as a major upwelling event (MUE).

2.3. Environmental Conditions

Subsurface T (22.7–25.6 °C) and DO (6.36–7.04 mg·L⁻¹) decreased throughout the study period from May through August, while salinity remained relatively consistent (36.58 ± 0.4). A decrease in T was observed after the 14th of July 2018, which corresponded to MUE (Figure 2). Inorganic dissolved Si concentrations ranged from 0.99 mg·L⁻¹ to 2.27 mg·L⁻¹ with an average value of 1.41 ± 0.39 mg·L⁻¹. The MUE had no apparent effect on Si concentration (p = 0.82). Dissolved iron concentrations were below detection limits (0.007 mg·L⁻¹). Dissolved inorganic N within the system was dominated by NH₄⁺ with no NO₂ detected. Total dissolved N concentrations fluctuated within each region and ranged between 1.60 and 4.15 µM·L⁻¹ (Figure 3A). Though not statistically significant, the soft substrate site consistently had lower N concentrations (mean 2.15 ± 0.62 µM·L⁻¹) than the reef sites (mean 2.66 ± 0.78 µM·L⁻¹). Dissolved inorganic P concentration ranged from 0.161 µM·L⁻¹ to 0.36 µM·L⁻¹ throughout the study period with an average concentration of 0.21 ± 0.04 µM·L⁻¹ (Figure 3B). The soft substrate site tended to have a higher concentration of P with a mean of 0.25 ± 0.07 µM·L⁻¹, while reef sites had a collectively similar mean of 0.196 ± 0.02 µM·L⁻¹. Prior to the upwelling event, the mean N:P ratio was 13.4 ± 3.92, and following the MUE, this ratio decreased to 10.48 ± 2.44 with the exclusion of two samples from the central site, which were high in nitrogen (Figure 3C).

2.4. Particulate Organic Matter

The ratio of carbon (C) to N (C:N) within the particulate matter fluctuated throughout the study period but increased following the MUE (Figure 4). Prior to the MUE, the small fraction had lower average concentrations of both C (0.0087 ± 0.0047 mg·L⁻¹) and N (0.0024 ± 0.0011 mg·L⁻¹) than the large fraction (0.0138 ± 0.0050 mgC·L⁻¹, 0.0031 ± 0.0011 mgN·L⁻¹) (Figure 5A). Regressions between cellular C:N had coefficients of 0.144 (R² = 0.395) and 0.196 (R² = 0.836) for small and large fractions, respectively. Average C:N ratio was 4.07 ± 1.59 for the small fraction and 4.42 ± 0.61 for the large fraction.

Particulate organic matter (POM) concentrations increased across all regions following the MUE (Figure 5B). The small fraction showed a decreased coefficient of 0.102 (R² = 0.881), whereas the large fraction maintained a similar relationship of 0.197 (R² = 0.90). Average cellular nutrient concentrations of small (0.0209 ± 0.0120 mgC·L⁻¹, 0.0038 ± 0.0013 mgN·L⁻¹) and large (0.0224 ± 0.0089 mgC·L⁻¹, 0.0046 ± 0.0019 mgN·L⁻¹) fractions increased, and average C:N ratios exhibited an increase of 0.24 for the small fraction and a 0.12-fold-change for the large fraction (Figure 5B). The high value measured within the small fraction corresponded to the high chlorophyll a seen on the 6th of August 2018 (Figure 6). One sample of the large fraction contained a crustacean larva (~10 mm), likely contributing to the high concentration of C and N.

2.5. Phytoplankton Composition

Chlorophyll a (chla) concentration was used as a proxy for productivity and was significantly affected by T (p < 0.05). Low chla concentrations were observed between the 24th of May through the 20th of July and increased following the MUE event (Figure 6). Concentrations ranged from 4 to 234 ng·L⁻¹ prior to the MUE and from 19 to 794 ng·L⁻¹ following the MUE.

The phytoplankton community was dominated by diatoms (Figure 7) based on photopigment signatures [64] and subsequent CHEMTAX analysis (Table S1) [65]. The diatom assemblage was composed of chain-forming genera including Pseudo-nitzschia spp., Bacteriastrum spp., Thalassiosira spp., Guinardia spp., and Chaetoceros spp., and the non-chain forming Rhizosolenia spp. Prevalent dinoflagellates included Protoperidinium spp., Ceratium spp., and Pyrophacus spp. Alexandrium spp., Dinophysis spp., Prorocentrum spp., and Akashiwo spp. were also observed in low abundance and have the ability to produce phycotoxins. The primary species identified within the “other” category were Phaeocystis spp., a prymnesiophyte, and Dictyocha spp., a unicellular flagellate.

Trichodesmium spp. tended to be more abundant prior to the MUE (prior to MUE = 4.62 ± 6%, post MUE = 0.66 ± 1.5%), whereas Pseudo-nitzschia spp. increased in community proportion across all regions following the MUE. Pseudo-nitzschia spp. made up 13.96 ± 15.6% of the community before the MUE, and 40.11 ± 17.5% after the MUE. The soft substrate site had the highest change in Pseudo-nitzschia spp. proportion, with a mean 25.2 ± 1.15% before the MUE, and became the dominant species with a mean of 54.5 ± 4.9% following the MUE.

2.6. Particulate DA

Surface samples collected as a pilot study between January 2017 and April 2018 combined with samples collected during the intensive sampling period of this study (May–August 2018) revealed the presence of pDA year-round (Figure 8). During pilot sampling, pDA was detected in 96% (n = 47) of samples with a range of 0.241 to 49.4 pg pDA·L⁻¹. The highest concentration (49.4 pg pDA·L⁻¹) was observed on the 11th of June 2017 at the soft substrate region. During the intensive sampling period, there was no statistical difference in pDA between the small fraction and large fraction samples (p = 0.62), so values reported are a combined average. During the intensive sampling period, pDA was detected in 80% of samples with a range of 0.295 to 11.71 pg pDA·L⁻¹. The highest concentration (11.7 pg pDA·L⁻¹) was observed at the soft substrate region on the 27th of July 2018. There was a greater than 2-fold increase of pDA following the MUE. Prior to the MUE, pDA concentrations averaged 0.937 ± 1.25 pg L⁻¹ (70% detected, n = 20). All samples following the MUE (n = 15) had detectable pDA with an average concentration of 4.81 ± 2.41 pg L⁻¹. Principle components analysis (PCA; Figure S1) suggested that pDA is negatively influenced by the N:P ratio; pDA is higher when the N:P ratio is low.

3. Discussion

Our results indicate that wind-driven upwelling likely reduced N limitation, which allowed for phytoplankton biomass to increase and supported a shift to a Pseudo-nitzschia spp.-dominated diatom assemblage. The relative increase in Pseudo-nitzschia spp. corresponded with an increase in pDA in the vicinity of Praia do Tofo, Mozambique. To our knowledge, this is the first identification of Pseudo-nitzschia spp. and DA in Inhambane Province, Mozambique. Preliminary data suggest a seasonal spike at the onset of Austral winter, with peaks in June of 2017 and July of 2018. Though pDA peaked after a coastal upwelling and decrease in T event, pDA was noted throughout the year. The Mozambique channel is known for having turbid waters due to the formation of eddies within the Mozambique current [66,67], primarily occurring during Austral winter. These mesoscale eddies play an important role in the interannual variability of phytoplankton blooms [63,66] as well as increased biodiversity of higher trophic levels [68]. During our study, satellite imagery identified the occurrence of a cyclonic eddy influencing the coastal waters of Inhambane Province, which corresponded to an increase in the calculated upwelling flux (Figure S2). It is possible that the increase in Pseudo-nitzschia spp. and pDA observed in our study were influenced by both wind-driven upwelling and the presence of the mesoscale eddy, and further research is necessary to determine if this trend occurs annually.

On a global scale, DA has been prevalent among EBS and in temperate climate regions where upwelling of nutrients promote toxicity [69,70]. Warming global climate conditions and increased nutrients from terrestrial sources create prime conditions for increases in community dominance by Pseudo-nitzschia spp. as well as increased toxicity [26,32,69,71,72,73]. There has been a recent emergence of DA contamination of shellfish fisheries within the WBS including Argentina [18,74], Uruguay [75], Brazil [76], New Zealand [36], and South Africa [44], supporting the increased range of this HAB species. In 2016, the Gulf of Maine recorded its first shellfish fishery closure due to DA concentrations exceeding regulatory limits [24]. Pseudo-nitzschia blooms in the Gulf of Maine occur in the fall and persist into the winter, similar to the Pseudo-nitzschia bloom and pDA identified in this study near Praia do Tofo, Mozambique. This study is the first documentation of DA along the eastern Africa within a subtropical WBS.

The upwelling event documented here correlated with changes in inorganic nutrient availability and changes to productivity as indicated by increases in chla and seston nutrient composition. Based on global averages of nutrient ratios within marine systems, the Redfield Ratio of 106_C:16_N:1_P:1_Si:0.01_Fe is used as a baseline for determining nutrient limitations [77]. Analysis of both the inorganic dissolved (N:P 13.4 ± 3.92) and seston (C:N 4.8 ± 1.8) nutrient composition suggested the system was slightly N limited prior to the MUE. Silica (Si) is another important nutrient for diatoms [78,79], with averaged ratios of 2:1 (N:Si) and 7:1 (Si:P). The low N:Si ratio further suggests that N is the primary limiting nutrient for diatom growth, and the high Si:P ratio suggests that Si is likely not a limiting nutrient. While low silica conditions have been suggested to provide Pseudo-nitzschia spp. an advantage over other diatom species [80], our data indicate high Si available during a transition to a Pseudo-nitzschia spp.-dominated phytoplankton assemblage, indicating a different driver was responsible for the shift. While the system experienced a decrease in dissolved nutrient ratio (N:P 10.4 ± 2.44, Figure 3C) following the MUE, we suspect that it was likely due to the rapid uptake of nutrients by phytoplankton as explained by the increase in particulate C and N (Figure 5B). The increase in C:N ratio after the MUE (Figure 4) resulted in a similar C:N ratio to the rest of the southern Indian Ocean (6.5–9.5) [81].

Energetically, high C:N ratios are more efficient at converting energy between trophic levels and therefore are considered more nutritious for primary consumers [82]. Increased seston C:N ratios following the MUE could be a crude indication of higher quality diet; however, Pseudo-nitzschia spp. became dominant in response to the MUE, which could have deleterious nutritional consequences on consumers due to decreased diet diversity and the potential for toxin consumption. Fatty acid analysis is suggested as a better indicator of nutritional quality for future research [83], as it measures the essential nutrients often not provided within a monospecific diet [82].

Pseudo-nitzschia spp. is well-adapted to shifts in environmental conditions, taking advantage of periotic disturbances such as upwelling events, eddies, or experimental manipulation [30,84]. Following the MUE, our data suggest Pseudo-nitzschia spp. was able to take advantage of the periotic disturbance and dominate the community assemblage within the system. Within our study, variability in environmental parameters such as T and nutrients was low between sampling sites; however, some trends require further investigation. Soft substrate sites tended to be low in dissolved N (2.51 ± 0.479 mg·L⁻¹) and high in P (0.255 ± 0.074 mg·L⁻¹), suggesting N limitation. These sites also consistently had the highest pDA (1.65–11.7 pg pDA L⁻¹). PCA analysis suggested a correlation between the low nutrient ratios and Pseudo-nitzschia spp. dominance (Figure S1). Another previously explored aspect to this relationship suggests a restriction of DA production under N limited conditions resulting from an insufficient pool of free N for synthesis of the nitrogen-rich toxin [85]. However, a study with multiple species of Pseudo-nitzschia revealed various growth responses to nutrient ratios [86], suggesting a more complex and potentially species-specific relationship with nutrients.

Overall, pDA concentrations within our study period were low (0.091–11.71 pg pDA L⁻¹) when detected (Figure 8). The solubility of the toxin in water prevents DA from accumulating within ecosystems without a source of production [87,88]; therefore, in areas where DA is known to occur with a seasonal bloom, monitoring is focused during blooms. However, along Inhambane Province, there could be concern for chronic exposure to the ecosystem and possibly human health [89] due to the year-round production of DA, which has been shown to cause declines in cognitive function to the central nervous system [90,91,92]. This study only focused on quantifying pDA within plankton samples, and further investigation into vector organisms is needed to assess the threat of chronic DA exposure.

High concentrations of DA did not correspond to high densities of Pseudo-nitzschia spp., corroborating previous studies [24,93,94]. High variance in cell abundance suggests high day-to-day growth fluctuations in the region. Light microscopy identification suggested at least three Pseudo-nitzschia spp. present simultaneously along the coast of Inhambane Province, namely P. multiseries (Figure S3A,B), P. australis (Figure S3C–E), and P. turgidula (Figure S3F,G) [38,44,70], which are all known to produce DA [24]. Prior reports have also identified P. multiseries, P. pungens, P. cuspidate, and P. seriata along the eastern coast of Africa [44,95]. Each individual species of Pseudo-nitzschia exhibits specific growth and toxin production dynamics in response to temperature, salinity, irradiance, photoperiod, and variation in nutrients [24], which suggests both cell abundance and toxin concentration must be included in monitoring efforts.

Along the coast of Inhambane Province, the highest toxin concentration was measured in the central, soft substrate region (Figure 1), a common feeding area for aggregations of Rhincodon typus, Manta alfredi, and M. birostris [55]. Though not significantly different (p = > 0.05), pDA concentrations within the large plankton fraction at the soft substrate region were repeatedly higher than those in the small plankton fraction. Not enough data were collected to support the threat of bioaccumulation; however, there is a pressing need to better understand the risks of toxin bioaccumulation in order to advise local residents on this potential public health issue.

As climate changes, the Praia do Tofo region would benefit from a better understanding of the effects of seasonal influences on plankton assemblage and community dynamics to predict how the system might be affected in the upcoming decades. The current rate of temperature increase in the Indian Ocean will likely have profound effects on surrounding coastal communities [96]. The fluctuation of Trichodesmium spp. abundance within the community composition was not significant within the study period but is noted to increase seasonally during the monsoon season in nearby Tanzania [97], promoting the recycling of limited N through nitrogen fixation properties. Trichodesmium spp. is also recognized as a harmful species in this area [40] and is predicted to increase in prevalence under elevated surface temperature conditions. Changes to mesoscale processes affecting environmental conditions and nutrient transport are predicted to occur within the Southwestern Indian Ocean [6,98], possibly impacting the factors influencing Pseudo-nitzschia spp. assemblage dominance and pDA concentration. Monsoon events, which can affect Mozambique annually from January through March, are predicted to intensify and become more frequent. These events have been shown to increase productivity in the Northern Indian Ocean [99] and south of Madagascar [100]. In addition, increased nitrate during the monsoon season is documented to promote Pseudo-nitzschia blooms in Tanzania [101]. Increased rain intensity has also promoted harmful cyanobacteria along coastlines in Tanzania [97]. However, cyclones are a relatively recent threat to southern coastal systems and have not been studied along the coast of Mozambique. Warming is predicted to have delayed effects on the Indian Ocean Dipole, which accounts for interannual variability in rainfall and is likely to decrease total rainfall [102]. Overall, rainfall in southern Africa is expected to become less frequent, but to increase in intensity [103]. These processes are likely to affect coastal systems with decreased terrestrial runoff, such as Inhambane Province, Mozambique.

Anthropogenic nutrient loading can increase growth of harmful algae, including Pseudo-nitzschia spp. [28]. Agriculture [104] and aquaculture [59,105] industries are forecasted to continue to grow in Mozambique. As coastal communities develop and increase agriculture and aquaculture production, nutrient levels of freshwater runoff will likely change the affective nutrient input to the coastal phytoplankton communities. Data from this research suggest that the coastal system along Inhambane Province was not highly impacted by freshwater input during the dry season. However, additional research during the wet season (Austral summer) would be necessary to assess if changes in land use will impact this region. Coastal Inhambane Province is a biological hot spot reliant on local fisheries as a food source and economic development through tourism due to the year-round presence of actively feeding planktivores [52,57], and changes to the coastal system can have compounding affects to the rural community. Further monitoring and education of HABs is needed in the region to ensure ecosystem, human, and wildlife health.

4. Methods

4.1. Study Site

The coastal waters of Southern Mozambique feed the Agulhas Current, one of the strongest WBSs in the world [106,107,108]. This current system is influenced by the East Madagascar Current and eddies within the Mozambique Channel [108]. Water flowing south within the Mozambique Channel along the steep slope and narrow shelf [109] forms both cyclonic and anticyclonic eddies [108]. Westward cyclonic (cold core) eddies are common along the south of Madagascar during Austral winter with a lifetime of a weeks to months [60,66]. These waters are primarily oligotrophic; however, mesoscale features such as coastal upwelling caused by wind stress and cyclonic, cold-core eddies promote increased primary production along the southeastern coast of Africa [61,107].

The study site (Figure 1) was located along the coast of Inhambane Province, Mozambique, stretching approximately 40 km from Barra Beach (23°47′35.00 S, 35°31′05.00 E) to Praia de Jangamo (24°03′52 S, 35°29′33 E). The region is defined by a narrow continental shelf with a steep slope [109], warm waters (21–29 °C), and wind-driven coastal upwelling [107], providing for a diverse coral reef complex. This semi-tropical system is located slightly south of the Tropic of Capricorn and exhibits many characteristics of tropical systems, including abundant light availability and the occurrence of tropical species such as nitrogen-fixing Trichodesmium [110] and hard corals [111]; however, the system experiences annual fluctuations in sea surface temperature similar to temperate areas. Based on the geophysical and benthic characteristics, four regions were chosen for detailed study: northern reefs, central reefs, soft substrate, and southern reefs. The northern reefs are dominated by hard plated corals such as Acropora spp. (plate coral), Porites spp. (porous coral), and Favia spp. (false honeycomb coral). The site is marked by strong currents typically flowing from the northwest to the southeast and is also the closest location to Baia de Inhambane and the estuary system. Reefs within the north site range between depths of 23 and 31 m. While the area provides habitat for large aggregations of schooling fish, due to the variable currents, this region is less frequented by recreational divers.

The central reef sites are sheltered within Praia do Tofo coastal bay. Moderate currents flow from north to south along the coast, allowing for the development of branching hard corals such as Tubastrea micranthus (green tree coral), as well as soft corals such as Dendronephthya spp. (branching soft coral). Depth at reefs ranged between 18 and 30 m. Based on its close proximity to Praia do Tofo and depth, this location is heavily trafficked by recreational divers and fishermen.

The central soft substrate region ranging from 14 m to 18 m depth is located on the leeward side of Tofino Point, approximately 5 km south of Praia do Tofo. The site is an established year-round feeding location of R. typus, M. birostris, and M. alfredi, and it is frequently visited by recreational snorkelers.

The southern reefs are dominated by rocky substrate and soft corals, such as Lobophytum spp. (leather coral), Sarcophyton spp. (fleshy soft coral), and Dendronephthya spp. (branching soft coral). The reefs are located approximately 20 km south of Praia do Tofo and at a depth of 18–27 m.

4.2. Satellite Imagery and Analysis

MODIS-Aqua satellite Level 2 images were obtained from the NASA Goddard Space Flight Center (GSFC) Ocean Color data (https://oceancolor.gsfc.nasa.gov, accessed on the 26 September 2018). All images within the sampling period were visually assessed for cloud contamination. Data were processed using SeaDAS v7.4 and re-projected with a 1 km spatial resolution.

4.3. Upwelling Index Calculations

Upwelling index (UI), characterized by the volume of water transported through an along-shore transect by a time and distance unit, was calculated using the methodology proposed by Bakun [112] as indicated in technical reports published by the Instituto Español de Oceanografía [113,114]. Sea level pressure field over the ocean is required to calculate UI; however, this study used meteorological data from four coastal World Meteorological Organization (WMO) stations (Vilanculos, #67315; Inhambane, #67323; Panda, #67327; and Xai-Xai, #67335) to calculate wind force components, which were subsequently used to approximate UI at −23.5 °S, 35.3 °E. Pressure data used for the UI calculation were daily average values (hourly). The meridional (v) and zonal (u) wind components were calculated using the following equations:

(1)$v=\frac{R·T}{p_a}·\frac{\Delta p_x}{L·f}$

(2)$u=\frac{R·T}{p_a}·\frac{\Delta p_y}{L·f}$

where R is the ideal gas constant (R = 287 j Kg⁻¹ K⁻¹), T is the standard atmosphere temperature (T = 288 K), p_a is standard atmospheric pressure (p_a = 1013.5 hPa), Dp_x and Dp_y are, respectively, the zonal and meridional pressure difference (hPa), L (km) is the distance between the locations used to calculate pressure gradient, and f is the Coriolis parameter defined as twice the vertical component of the Earth’s angular velocity (Ω = 7.3$\times$10⁻⁵ rad s⁻¹) about the local vertical or f = 2Ωsin (θ) at latitude θ.

Ekman transport (Q, m³ s⁻¹ km⁻¹) was calculated using the components of wind speed (u and v), seawater density (ρ_w = 1025 kg m⁻³), a dimensionless empirical drag coefficient (C_d = 1.4 × 10⁻³), and air density (ρ_a = 1.2 kg m⁻³, normal conditions) by means of

(3)$Q_x=\frac{{\mathsf{\tau }}_y}{f·{\mathsf{\rho }}_w}$

(4)${\mathsf{\tau }}_y={\mathsf{\rho }}_a·C_d·\sqrt{u^2+v^2}·v$

where the x subscript corresponds to the zonal component and the y subscript to the meridional component. UI is defined as –Q_x (m³ s⁻¹ km⁻¹), which is the volume transported per distance unit of an alongshore section. The sign of Ekman transport is changed to define positive (negative) values of UI as response of upwelling (downwelling) favorable winds.

4.4. Field Site Monitoring and Sample Collection

Plankton surveys along the coast of Praia do Tofo were designed to investigate the temporal and spatial patterns of DA-producing Pseudo-nitzschia abundance. Ships of opportunity, through Peri-Peri Divers, were used to conduct environmental and plankton sampling between January 2017 and August 2018. A pilot monitoring period consisted of collection of whole water samples at least monthly in the soft substrate region pDA analyses from January 2017 to April 2018. Sea surface temperatures were also measured at this time. Measurements and sample collections were a part of a citizen scientist initiative, in conjunction with All Out Africa and Peri-Peri Divers, to explore the plankton community in a known R. typus feeding location. An intense study took place from the 20 May 2018 to 10 August 2018. During this time, sampling was planned in each region approximately once per week within an 11 week period; however, weather and access limited sampling events. There was a total of nine sites within the four regions. Samples were collected as the boat drifted above the reef site. Physicochemical parameters (T, conductivity, and DO) were recorded at GPS waypoint markings at the beginning of each sample collection using a YSI ProDSS multiparameter water quality meter (model #626870-1, #627150-4). Calibrations for conductivity (using certified calibration solution at 50,000 µS/cm) and DO were completed prior to each use. Weather activity, wind direction, percent cloud cover, and surface activity (Beaufort scale) were also recorded. Wind data were gathered from World Meteorological Organization station #67323 located at the Inhambane Airport (FQIN).

Approximately 30 L of water was collected at a depth of 5 m via Niskin (General Oceanics Model 1010) during each sampling event. Roughly 1 L of whole water was transferred to an acid-washed container for nutrient analysis. Samples for dissolved inorganic N (NO₂, NO₃, and NH₄), P (PO₄), Si, and iron (Fe) analysis were filtered through glass fiber filters using gentle vacuum and collected in water quality sterile sample bags (VWR # 89085-546) and frozen. The remaining water was gently poured through consecutive 80 µm mesh (i.e., large fraction) and 36 µm mesh (i.e., small fraction). The larger 80 µm size mesh was used to separate zooplankton from Pseudo-nitzschia spp. cells based on the assumption that Pseudo-nitzschia cells would pass through, which was confirmed by microscopic analysis. The smaller 36 µm mesh was used to target Pseudo-nitzschia cells for pDA analysis. The particulate matter collected on the 80 µm and 36 µm mesh (hereafter referred to as “large” and “small” fraction, respectively) was transported on ice and in the dark to All Out Africa Marine Research Station for processing and preservation. Subsamples of small (duplicate) and large (single) fractions for pDA analysis were filtered onto 25 mm glass fiber filters and frozen. A subsample of the small fraction was preserved with acidified Lugols and used to determine Pseudo-nitzschia spp. abundance using the Utermöhl method [115]. Duplicate subsamples from both small and large fractions were used for organic C and N analysis; samples were filtered through combusted glass fiber filters under gentle vacuum and frozen for later analysis. Analyses were conducted following standard procedures by Bianchi and Bauer [116]. All samples were frozen or preserved and transported to Louisiana State University or University of Minnesota Duluth Large Lakes Observatory for analysis.

4.5. Dissolved Nutrient Analysis

Ammonium, nitrate + nitrite, and orthophosphate concentrations were determined using standard AutoAnalyzer techniques (EPA 353.4, EPA 350.1, and EPA 365.5). Dissolved inorganic N was calculated as the sum of nitrate, nitrite, and ammonium concentrations, and dissolved inorganic P was considered equal to the orthophosphate concentration. Fe and Si were analyzed by ICP-OES using standard methods (EPA #200.7). Quantification limits were 0.1 µM for nitrate, 0.2 µM for nitrite, 0.1 µM for ammonia, 0.04 µM for P, 0.012 mg L⁻¹ for Si, and 0.007 mg L⁻¹ for Fe.

4.6. Pigment Analysis

Duplicate samples for chla were filtered through glass fiber filters under gentle vacuum for the small fraction. Filters were frozen and stored for photopigment analysis by high performance liquid chromatography (HPLC) following the method from Pinckney et al. [64]. Filters were sonicated in 3 mL of 100% acetone for 30 s and extracted in the dark for 20 to 24 h at −20 °C. Extracts were filtered through 0.2 µm filter, and 300 µL was injected into an HPLC system equipped with reverse-phase C18 columns in series (Rainin Microsorb-MV, 0.46 × 10 cm, 3 mm, Vydac 201TP, 0.46 × 25 cm, 5 mm). A nonlinear binary gradient adapted from Van Heukelem et al. [117] was used for pigment separations. Solvent A consisted of 80% methanol and 20% ammonium acetate (0.5 m adjusted to pH 7.2), and solvent B consisted of 80% methanol and 20% acetone. Absorption spectra and chromatograms were acquired using a Shimadzu SPD-M10av photodiode array detector, where pigment peaks were quantified at 440 nm.

Photopigment concentrations were examined with CHEMTAX (v.1.95) [65] to determine the absolute abundance of major phytoplankton groups (µg chl a⁻¹) using the convergence procedure outlined by Latasa [118]. The initial pigment matrix was derived from Barlow et al. [60] and modified using Higgins et al. [119] (Table S1). Measured pigments include chlorophyll a (chla), chlorophyll b (Chl b), chlorophyll c₁ (Chlc1), chlorophyll c₂ (Chlc2), chlorophyll c₃ (Chlc3), peridinin (Per), 19′-butanoyloxyfucoxanthin (But), fucoxanthin (Fuc), neoxanthin (Neo), violaxanthin (Viol), prasinoxanthin (Pras), 19′-hexanoylfucoxanthin (Hex), alloxanthin (Allo), zeaxanthin (Zea), antheraxanthin (Anth), and lutein (Lut). Previous research suggests that the majority of the phytoplankton community is represented within the 20–200 µm size range [120]; phytoplankton groups whose size were <36 µm were excluded from the matrix. Pigment ratio values for cyanobacteria were derived from Higgins et al. [119] and reported as Trichodesmium. Pseudo-nitzschia spp. was distinguished from other diatoms through the presence of Chlc1 [119].

A portion of the small fraction samples (13%) were analyzed via microscopy to verify CHEMTAX results. A settling procedure [121] was used to concentrate the samples to 5 mL, and cell counts were performed using a Sedgwick–Rafter counting cell on an inverted microscope at 10× objective (ZIESS Axis Observer-A1 axiovert 135) and identified to the genus level. A subsample of the large fraction was collected and preserved with 5% acidified Lugols, and organisms were identified to the genus level using a dissecting scope (ZIESS SteREO Discovery.V8). Chaetocerous spp. were noted within several samples among the large fraction; however, due to preservation effects, biomass could not be accurately quantified, so presence/absence was noted.

4.7. Particulate Domoic Acid Analysis

Concentration of pDA was assessed within three months of collection by Biosense ASP ELISA (Amnesic Shellfish Poison, enzyme-linked immunosorbent assay) (Biosense #A31300401) using the protocol described in Baustian et al. [122] expressed as pg pDA L⁻¹. Samples were extracted using 20% methanol, sonicated, then filtered to remove cell debris. Plates were read using a Thermo Varioskan Flash spectral scanning multimode reader and SkanIt Software version 2.4.1. Toxin concentrations (pDA) were calculated using a non-linear 4-parameter logistic curve fit provided by Abraxis (Abraxis, Inc, Warminster, PA, USA). Detection limit was 9.33 ± 1.1 pg DA mL⁻¹ across the analyzed plates; concentrations below detection limit were reported as below detection limit.

4.8. Statistical Analysis

Multivariate analyses were performed using SigmaPlot v14 (Inpixon, Palo Alto, CA, USA). Distinct one-way analysis of variance (ANOVA) for temperature, dissolved oxygen, and salinity parameters were used to compare regions. Pearson correlations (Table S2) for physicochemical and biological variables in surface water throughout the study were tested. Significance was defined as a p value < 0.05. PCA analyses were performed to characterize the relationships between physicochemical parameters (T, salinity, and nutrients) and biological parameters (chla, Pseudo-nitzschia cell counts, and pDA).

Author Contributions

Conceptualization, H.K. and R.M.E.; methodology, H.K. and R.M.E.; validation, R.M.E., K.M.S. and S.B.; formal analysis, H.K. and R.M.E.; investigation, H.K. and R.M.E.; resources, K.G.R.; data curation, H.K. and K.E.R.-A.; writing-original draft preparation, H.K.; writing-review and editing, R.M.E., K.M.S., S.B. and K.E.R.-A.; supervision, R.M.E. and K.E.R.-A.; project administration, R.M.E.; funding acquisition, R.M.E. and K.G.R. All authors have read and agreed to the published version of the manuscript.

Funding

Funding for this project was provided by the Louisiana State University AgCenter. A sincere thanks to Peri Peri Divers and staff for their collaboration and help collecting data for this project. This project would not have been possible without the support of All Out Africa Marine Research staff and volunteers.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Role of Coastal Environmental Conditions During Austral Winter on Plankton Community Dynamics and the Occurrence of Pseudo-nitzschia spp. and Domoic Acid in Inhambane Province, Mozambique. Available online: https://digitalcommons.lsu.edu/gradschool_theses/5056 (accessed on 15 November 2021).

Conflicts of Interest

The authors declare no conflict of interest.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figures

Figure 1. Study site: Inhambane Province, Mozambique. Four areas within the region were identified between Barra and Praia de Jangamo (40 km). During the pilot monitoring period between January 2017 and April 2018, only the soft substrate region was sampled. During the intensive study period, all four regions were sampled, but no statistical significance was noted.

Figure 2. Calculated upwelling index. Atmospheric pressure data were used to estimate upwelling transport flux (m3 s−1 km−1) along the coast of Inhambane Province. The dotted line shows the weekly average of transport flux. The box represents the timeframe of the intensive study period. A prolonged increase in upwelling flux occurred in late July through August, corresponding to multiple biological responses, referred to as a major upwelling event (MUE).

Figure 3. Inorganic dissolved nutrients. (A) Total dissolved nitrogen (N), (B) total dissolved phosphorus (P), and (C) N:P ratio over the intense study period (May–August 2018) indicated the system is primarily controlled by nitrogen limitation. The dotted line indicates the MUE that occurred following 19 July 2018.

Figure 4. Particulate carbon (C) to N ratio (C:N) from May to August 2018. Total (combined fractions) ratios fluctuated throughout the study period, but average ratio increased following the MUE (dotted line).

Figure 5. Particulate organic matter (POM). (A) Prior to the MUE, regressions between cellular C:N had coefficients of 0.144 (R2 = 0.395) and 0.196 (R2 = 0.836) for small and large fractions, respectively. (B) Following the MUE, the small fraction showed a decreased coefficient of 0.102 (R2 = 0.881), whereas the large fraction maintained a similar relationship of 0.197 (R2 = 0.90).

Figure 6. Chlorophyll a (chla) concentrations and upwelling index for May to August 2018. Chla, used as a proxy for primary productivity, increased in late July into August, corresponding to MUE.

Figure 7. Phytoplankton assemblage composition. CHEMTAX analysis [65] revealed the assemblage was diatom-dominated (shades of green) throughout the study period. Following the MUE on 19 July 2018 (indicated on timeline), there was an increase in assemblage dominance by Pseudo-nitzschia spp. (light green).

Figure 8. Particulate domoic acid (pDA) concentration in coastal Inhambane Province, Mozambique. Samples collected from January 2017 through August 2018 detected toxin year-round. The shaded box highlights samples collected during the intensive study during May–August 2018. Data suggest a possible spike in pDA production during the onset of Austral winter.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/toxins13120903/s1, Table S1. Pigment: chlorophyll a starting and output ratios in the CHEMTAX analysis of HPLC pigments; Table S2. Pearson correlation for physiochemical and biological variables over study period (May–August 2018); Figure S1 Principle Components Analysis. Variable factors from May to August 2018 defined by the two first axes showing the link between physiochemical parameters, chla, pDA and Pseudo-nitzschia spp. abundance; Figure S2 Satellite Imagery. Sea surface temperature (SST) and chlorophyll a (chla) images from 28 July 2018 at 13:15 LST reveal the occurrence of a cyclonic eddy along the coast of Inhambane Province; Figure S3 Pseudo-nitzschia spp. light microscopy analysis revealed at least three species of Pseudo-nitzschia spp. present in samples from coastal Inhambane Province based on width and shape.

References

1. Bakun, A.; Nelson, C.S. The seasonal cycle of wind-stress curl in subtropical Eastern Boundary Current regions. J. Phys. Oceanogr.; 1991; 21, pp. 1815-1834. [DOI: https://dx.doi.org/10.1175/1520-0485(1991)021<1815:TSCOWS>2.0.CO;2]

2. Carr, M.-E.; Kearns, E.J. Production regimes in four Eastern Boundary Current systems. Deep Sea Res. Pt. II Top. Stud. Oceanogr.; 2003; 50, pp. 3199-3221. [DOI: https://dx.doi.org/10.1016/j.dsr2.2003.07.015]

3. Doblin, M.A.; Petrou, K.; Sinutok, S.; Seymour, J.R.; Messer, L.F.; Brown, M.V.; Norman, L.; Everett, J.D.; McInnes, A.S.; Ralph, P.J. et al. Nutrient uplift in a cyclonic eddy increases diversity, primary productivity and iron demand of microbial communities relative to a western boundary current. Peer J.; 2016; 4, e1973. [DOI: https://dx.doi.org/10.7717/peerj.1973] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27168982]

4. Everett, J.D.; Doblin, M.A. Characterising primary productivity measurements across a dynamic western boundary current region. Deep Sea Res. Pt. I Oceanogr. Res. Pap.; 2015; 100, pp. 105-116. [DOI: https://dx.doi.org/10.1016/j.dsr.2015.02.010]

5. Kimura, S.; Kasai, A.; Nakata, H.; Sugimoto, T.; Simpson, J.H.; Cheok, J.V.S. Biological productivity of meso-scale eddies caused by frontal disturbances in the Kuroshio. ICES J. Mar. Sci.; 1997; 54, pp. 179-192. [DOI: https://dx.doi.org/10.1006/jmsc.1996.0209]

6. Behagle, N.; du Buisson, L.; Josse, E.; Lebourges-Dhaussy, A.; Roudaut, G.; Menard, F. Mesoscale features and micronekton in the Mozambique Channel: An acoustic approach. Deep Sea Res. Pt. Ii; 2014; 100, pp. 164-173. [DOI: https://dx.doi.org/10.1016/j.dsr2.2013.10.024]

7. Wells, R.J.D.; Rooker, J.R.; Quigg, A.; Wissel, B. Influence of mesoscale oceanographic features on pelagic food webs in the Gulf of Mexico. Mar. Biol.; 2017; 164, 92. [DOI: https://dx.doi.org/10.1007/s00227-017-3122-0]

8. Davidson, K.; Gowen, R.J.; Harrison, P.J.; Fleming, L.E.; Hoagland, P.; Moschonas, G. Anthropogenic nutrients and harmful algae in coastal waters. J. Environ. Manag.; 2014; 146, pp. 206-216. [DOI: https://dx.doi.org/10.1016/j.jenvman.2014.07.002] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25173729]

9. Graneli, E.; Turner, J.T. Ecology of Harmful Algae; Springer: Berlin/Heidelberg, Germany, 2006.

10. Paerl, H.W.; Paul, V.J. Climate change: Links to global expansion of harmful cyanobacteria. Water Res.; 2012; 46, pp. 1349-1363. [DOI: https://dx.doi.org/10.1016/j.watres.2011.08.002]

11. Smayda, T.J. Novel and nuisance phytoplankton blooms in the sea: Evidence for a global epidemic. Toxic Marine Phytoplankton; Graneli, E.; Sundstorm, B.; Edler, I.; Anderson, D.M. Elsevier: New York, NY, USA, 1990; pp. 29-40.

12. Van Dolah, F.M. Marine algal toxins: Origins, health effects, and their increased occurrence. Environ. Health Persp.; 2000; 108, pp. 133-141. [DOI: https://dx.doi.org/10.1289/ehp.00108s1133] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/10698729]

13. Anderson, D.M.; Hoagland, P.; Kaoru, Y.; White, A.W. Estimated Annual Economic Impacts from Harmful Algal Blooms (HABs) in the United States; Woods Hole Oceanographic Institution: Falmouth, MA, USA, 2000.

14. Grattan, L.M.; Holobaugh, S.; Morris, J.G. Harmful algal blooms and public health. Harmful Algae; 2016; 57, pp. 2-8. [DOI: https://dx.doi.org/10.1016/j.hal.2016.05.003] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27616971]

15. Willis, C.; Papathanasopoulou, E.; Russel, D.; Artioli, Y. Harmful algal blooms: The impacts on cultural ecosystem services and human well-being in a case study setting, Cornwall, UK. Mar. Policy; 2018; 97, pp. 232-238. [DOI: https://dx.doi.org/10.1016/j.marpol.2018.06.002]

16. Bates, S.S.; Bird, C.J.; Defreitas, A.S.W.; Foxall, R.; Gilgan, M.; Hanic, L.A.; Johnson, G.R.; McCulloch, A.W.; Odense, P.; Pocklington, R. et al. Pennate diatom Nitzschia-pungens as the primary source of domoic acid, a toxin in shellfish from eastern Prince Edward Island, Canada. Can. J. Fish. Aquat. Sci.; 1989; 46, pp. 1203-1215. [DOI: https://dx.doi.org/10.1139/f89-156]

17. Bargu, S.; Smith, E.; Ozhan, K. Toxic diatom Pseudo-nitzschia and its primary consumers (vectors). The Diatom World; Cellular Origin, Life in Extreme Habitats and Astrobiology; Springer: Berlin/Heidelberg, Germany, 2011; pp. 491-512.

18. D’Agostino, V.C.; Degrati, M.; Sastre, V.; Santinelli, N.; Krock, B.; Krohn, T.; Dans, S.L.; Hoffmeyer, M.S. Domoic acid in a marine pelagic food web: Exposure of southern right whales Eubalaena australis to domoic acid on the Peninsula Valdes calving ground, Argentina. Harmful Algae; 2017; 68, pp. 248-257. [DOI: https://dx.doi.org/10.1016/j.hal.2017.09.001]

19. Lefebvre, K.A.; Bargu, S.; Kieckhefer, T.; Silver, M.W. From sanddabs to blue whales: The pervasiveness of domoic acid. Toxicon; 2002; 40, pp. 971-977. [DOI: https://dx.doi.org/10.1016/S0041-0101(02)00093-4]

20. Work, T.M.; Barr, B.; Beale, A.M.; Fritz, L.; Quilliam, M.A.; Wright, J.L.C. Epidemiology of domoic acid poisoning in brown pelicans (Pelecanus occidentalis) and Brandt’s cormorants (Phalacrocorax penicillatus) in California. J. Zoo. Wildl. Med.; 1993; 24, pp. 54-62.

21. Zabaglo, K.; Chrapusta, E.; Bober, B.; Kaminski, A.; Adamski, M.; Bialczyk, J. Environmental roles and biological activity of domoic acid: A review. Algal Res.; 2016; 13, pp. 94-101. [DOI: https://dx.doi.org/10.1016/j.algal.2015.11.020]

22. Jeffery, B.; Barlow, T.; Moizer, K.; Paul, S.; Boyle, C. Amnesic shellfish poison. Food Chem. Toxicol.; 2004; 42, pp. 545-557. [DOI: https://dx.doi.org/10.1016/j.fct.2003.11.010]

23. Lefebvre, K.A.; Robertson, A. Domoic acid and human exposure risks: A review. Toxicon; 2010; 56, pp. 218-230. [DOI: https://dx.doi.org/10.1016/j.toxicon.2009.05.034]

24. Bates, S.S.; Hubbard, K.A.; Lundholm, N.; Montresor, M.; Leaw, C.P. Pseudo-nitzschia, Nitzschia, and domoic acid: New research since 2011. Harmful Algae; 2018; 79, pp. 3-43. [DOI: https://dx.doi.org/10.1016/j.hal.2018.06.001] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30420013]

25. Saeed, A.F.; Awan, S.A.; Ling, S.; Wang, R.; Wang, S. Domoic acid: Attributes, exposure risks, innovative detection techniques and therapeutics. Algal Res.; 2017; 24, pp. 97-110. [DOI: https://dx.doi.org/10.1016/j.algal.2017.02.007]

26. Trainer, V.L.; Bates, S.S.; Lundholm, N.; Thessen, A.E.; Cochlan, W.P.; Adams, N.G.; Trick, C.G. Pseudo-nitzschia physiological ecology, phylogeny, toxicity, monitoring and impacts on ecosystem health. Harmful Algae; 2012; 14, pp. 271-300. [DOI: https://dx.doi.org/10.1016/j.hal.2011.10.025]

27. Pitcher, G.; Pillar, S. Harmful algal blooms in eastern boundary upwelling systems. Prog. Oceanogr.; 2010; 85, pp. 1-4. [DOI: https://dx.doi.org/10.1016/j.pocean.2010.02.001]

28. Trainer, V.L.; Pitcher, G.C.; Reguera, B.; Smayda, T.J. The distribution and impacts of harmful algal bloom species in eastern boundary upwelling systems. Prog. Oceanogr.; 2010; 85, pp. 33-52. [DOI: https://dx.doi.org/10.1016/j.pocean.2010.02.003]

29. Fawcett, A.; Pitcher, G.C.; Bernard, S.; Cembella, A.D.; Kudela, R.M. Contrasting wind patterns and toxigenic phytoplankton in the southern Benguela upwelling system. Mar. Ecol. Prog. Ser.; 2007; 348, pp. 19-31. [DOI: https://dx.doi.org/10.3354/meps07027]

30. Trainer, V.L.; Adams, N.G.; Bill, B.D.; Stehr, C.M.; Wekell, J.C.; Moeller, P.; Busman, M.; Woodruff, D. Domoic acid production near California coastal upwelling zones, June 1998. Limnol. Oceanogr.; 2000; 45, pp. 1818-1833. [DOI: https://dx.doi.org/10.4319/lo.2000.45.8.1818]

31. Trainer, V.L.; Hickey, B.M.; Homer, R.A. Biological and physical dynamics of domoic acid production off the Washington coast. Limnol. Oceanogr.; 2002; 47, pp. 1438-1446. [DOI: https://dx.doi.org/10.4319/lo.2002.47.5.1438]

32. Anderson, C.R.; Sapiano, M.R.P.; Prasad, M.B.K.; Long, W.; Tango, P.J.; Brown, C.W.; Murtugudde, R. Predicting potentially toxigenic Pseudo-nitzschia blooms in the Chesapeake Bay. J. Mar. Syst.; 2010; 83, pp. 127-140. [DOI: https://dx.doi.org/10.1016/j.jmarsys.2010.04.003]

33. Jiang, R.; Wang, Y.-S. Modeling the ecosystem response to summer coastal upwelling in the northern South China Sea. Oceanologia; 2018; 60, pp. 32-51. [DOI: https://dx.doi.org/10.1016/j.oceano.2017.05.004]

34. Van Meerssche, E.; Pinckney, J.L. The influence of salinity in the domoic acid effect on estuarine phytoplankton communities. Harmful Algae; 2017; 69, pp. 65-74. [DOI: https://dx.doi.org/10.1016/j.hal.2017.10.003]

35. McKibben, S.M.; Peterson, W.; Wood, A.M.; Trainer, V.L.; Hunter, M.; White, A.E. Climatic regulation of the neurotoxin domoic acid. Proc. Natl. Acad. Sci. USA; 2017; 114, pp. 239-244. [DOI: https://dx.doi.org/10.1073/pnas.1606798114] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28069959]

36. Rhodes, L.; Jiang, W.; Knight, B.; Adamson, J.; Smith, K.; Langi, V.; Edgar, M. The genus Pseudo-nitzschia (Bacillariophyceae) in New Zealand: Analysis of the last decade’s monitoring data. N. Z. J. Mar. Fresh.; 2013; 47, pp. 490-503. [DOI: https://dx.doi.org/10.1080/00288330.2013.803489]

37. Rhodes, L.; Scholin, C.; Garthwaite, I. Pseudo-nitzschia in New Zealand and the role of DNA probes and immunoassays in refining marine biotoxin monitoring programmes. Nat. Toxins; 1998; 6, pp. 105-111. [DOI: https://dx.doi.org/10.1002/(SICI)1522-7189(199805/08)6:3/4<105::AID-NT13>3.0.CO;2-9]

38. Almandoz, G.O.; Fabro, E.; Ferrario, M.; Tillmann, U.; Cembella, A.; Krock, B. Species occurrence of the potentially toxigenic diatom genus Pseudo-nitzschia and the associated neurotoxin domoic acid in the Argentine Sea. Harmful Algae; 2017; 63, pp. 45-55. [DOI: https://dx.doi.org/10.1016/j.hal.2017.01.007] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28366399]

39. Negri, R.M.; Montoya, N.G.; Carreto, J.I.; Akselman, R.; Inza, D. Pseudo-nitzschia australis, Mytilus edulis, Engraulis anchoita, and domoic acid in the Argentine Sea. Harmful Algae; Reguera, B.; Blanco, J.; Fernandes, M.L.; Wyatt, T. Intergovernmental Oceanographic Commission of UNESCO: Paris, France, 2004.

40. Tamele, I.J.; Silva, M.; Vasconcelos, V. The incidence of marine toxins and the associated seafood poisoning episodes in the African countries of the Indian Ocean and the Red Sea. Toxins; 2019; 11, 58. [DOI: https://dx.doi.org/10.3390/toxins11010058]

41. Boisier, P.; Ranaivoson, G.; Rasolofonirina, N.; Andriamahefazafy, B.; Roux, J.; Chanteau, S.; Satake, M.; Yasumoto, T. Fatal mass poisoning in Madagascar following ingestion of a shark (Carcharhinus leucas): Clinical and epidemiological aspects and isolation of toxins. Toxicon; 1995; 33, pp. 1359-1364. [DOI: https://dx.doi.org/10.1016/0041-0101(95)00051-M]

42. Bouaïcha, N.; Chézeau, A.; Turquet, J.; Quod, J.P.; Puiseux-Dao, S. Morphological and toxicological variability of Prorocentrum lima clones isolated from four locations in the south-west Indian Ocean. Toxicon; 2001; 39, pp. 1195-1202. [DOI: https://dx.doi.org/10.1016/S0041-0101(00)00258-0]

43. Pitcher, G.C.; Calder, D. Harmful algal blooms of the southern Benguela Current: A review and appraisal of monitoring from 1989 to 1997. S. Afr. J. Mar. Sci.; 2000; 22, pp. 255-271. [DOI: https://dx.doi.org/10.2989/025776100784125681]

44. Pitcher, G.C.; Cembella, A.D.; Krock, B.; Macey, B.M.; Mansfield, L.; Probyn, T.A. Identification of the marine diatom Pseudo-nitzschia multiseries (Bacillariophyceae) as a source of the toxin domoic acid in Algoa Bay, South Africa. Afr. J. Mar. Sci.; 2014; 36, pp. 523-528. [DOI: https://dx.doi.org/10.2989/1814232X.2014.983971]

45. DAFF. South African Live Molluscan Shellfish Monitoring and Control Programme; Republic of South Africa Department of Agriculture, Forestry and Fisheries: Cape Town, South Africa, 2016.

46. Sá, C.; Leal, M.C.; Silva, A.; Nordez, S.; André, E.; Paula, J.; Brotas, V. Variation of phytoplankton assemblages along the Mozambique coast as revealed by HPLC and microscopy. J. Sea Res.; 2013; 79, pp. 1-11. [DOI: https://dx.doi.org/10.1016/j.seares.2013.01.001]

47. Barlow, R.; Stuart, V.; Lutz, V.; Sessions, H.; Sathyendranath, S.; Platt, T.; Kyewalyanga, M.; Clementson, L.; Fukasawa, M.; Watanabe, S. et al. Seasonal pigment patterns of surface phytoplankton in the subtropical southern hemisphere. Deep Sea Res. Pt. I; 2007; 54, pp. 1687-1703. [DOI: https://dx.doi.org/10.1016/j.dsr.2007.06.010]

48. Ternon, J.F.; Bach, P.; Barlow, R.; Huggett, J.; Jaquemet, S.; Marsac, F.; Ménard, F.; Penven, P.; Potier, M.; Roberts, M.J. The Mozambique Channel: From physics to upper trophic levels. Deep Sea Res. Pt. II Top. Stud. Oceanogr.; 2014; 100, pp. 1-9. [DOI: https://dx.doi.org/10.1016/j.dsr2.2013.10.012]

49. Tew Kai, E.; Marsac, F. Influence of mesoscale eddies on spatial structuring of top predators’ communities in the Mozambique Channel. Prog. Oceanogr.; 2010; 86, pp. 214-223. [DOI: https://dx.doi.org/10.1016/j.pocean.2010.04.010]

50. Tew-Kai, E.; Marsac, F. Patterns of variability of sea surface chlorophyll in the Mozambique Channel: A quantitative approach. J. Mar. Syst.; 2009; 77, pp. 77-88. [DOI: https://dx.doi.org/10.1016/j.jmarsys.2008.11.007]

51. Fordyce, A. Reef fishes of Praia do Tofo and Praia da Barra, Inhambane, Mozambique. West. Indian O. J. Mar. Sci.; 2018; 17, pp. 71-91.

52. Rohner, C.A.; Richardson, A.J.; Jaine, F.R.A.; Bennett, M.B.; Weeks, S.J.; Cliff, G.; Robinson, D.P.; Reeve-Arnold, K.E.; Pierce, S.J. Satellite tagging highlights the importance of productive Mozambican coastal waters to the ecology and conservation of whale sharks. Peer J.; 2018; 6, e4161. [DOI: https://dx.doi.org/10.7717/peerj.4161]

53. Pereira, M.A.M.; Litulo, C.; Santos, R.; Leal, M.C.; Fernandes, R.S.; Tibiriçá, Y.; Williams, J.; Atanassov, B.; Carreira, F.; Massingue, A. et al. Mozambique Marine Ecosystems Review; Biodinâmica/CTV: Maputo, Mozambique, 2014; 139.

54. Pierce, S.J.; White, W.T.; Marshall, A.D. New record of the smalleye stingray, Dasyatis microps (Myliobatiformes: Dasyatidae), from the western Indian Ocean. Zootaxa; 2008; 1734, pp. 65-68. [DOI: https://dx.doi.org/10.11646/zootaxa.1734.1.5]

55. Rohner, C.A.; Pierce, S.J.; Marshall, A.D.; Weeks, S.J.; Bennett, M.B.; Richardson, A.J. Trends in sightings and environmental influences on a coastal aggregation of manta rays and whale sharks. Mar. Ecol. Prog. Ser.; 2013; 482, pp. 153-168. [DOI: https://dx.doi.org/10.3354/meps10290]

56. Williams, J.L.; Pierce, S.J.; Rohner, C.A.; Fuentes, M.M.P.B.; Hamann, M. Spatial distribution and residency of green and loggerhead sea turtles using coastal reef habitats in Southern Mozambique. Front. Mar. Sci.; 2017; 3, 228. [DOI: https://dx.doi.org/10.3389/fmars.2016.00288]

57. Tibiriçá, Y.; Birtles, A.; Valentine, P.; Miller, D.K. Diving tourism in Mozambique: An opportunity at risk?. Tour. Mar. Environ.; 2011; 7, pp. 141-151. [DOI: https://dx.doi.org/10.3727/154427311X13195453162732]

58. FAO. The Future of Food and Agriculture—Trends and Challenges; FAO: Rome, Italy, 2017.

59. World Bank. Communities Livelihoods Fisheries: Fisheries Governance and Shared Growth in Mozambique; World Bank: Washington, DC, USA, 2018.

60. Barlow, R.; Lamont, T.; Gibberd, M.J.; Airs, R.; Jacobs, L.; Britz, K. Phytoplankton communities and acclimation in a cyclonic eddy in the southwest Indian Ocean. Deep Sea Res. Pt. I; 2017; 124, pp. 18-30. [DOI: https://dx.doi.org/10.1016/j.dsr.2017.03.013]

61. Barlow, R.; Lamont, T.; Morris, T.; Sessions, H.; van den Berg, M. Adaptation of phytoplankton communities to mesoscale eddies in the Mozambique Channel. Deep Sea Res. Pt. I; 2014; 100, pp. 106-118. [DOI: https://dx.doi.org/10.1016/j.dsr2.2013.10.020]

62. Olofsson, M.; Karlberg, M.; Lage, S.; Ploug, H. Phytoplankton community composition and primary production in the tropical tidal ecosystem, Maputo Bay (the Indian Ocean). J. Sea Res.; 2017; 125, pp. 18-25. [DOI: https://dx.doi.org/10.1016/j.seares.2017.05.007]

63. Raj, R.P.; Peter, B.N.; Pushpadas, D. Oceanic and atmospheric influences on the variability of phytoplankton bloom in the Southwestern Indian Ocean. J. Mar. Syst.; 2010; 82, pp. 217-229. [DOI: https://dx.doi.org/10.1016/j.jmarsys.2010.05.009]

64. Pinckney, J.L.; Millie, D.F.; Howe, K.E.; Paerl, H.W.; Hurley, J.P. Flow scintillation counting of 14C-labeled microalgal photosynthetic pigments. J. Plankton Res.; 1996; 18, pp. 1867-1880. [DOI: https://dx.doi.org/10.1093/plankt/18.10.1867]

65. Mackey, M.D.; Mackey, D.J.; Higgens, H.W.; Wright, S.W. CHEMTAX—A program for estimating class abundances from chemical markers: Application to HPLC measurements of phytoplankton. Mar. Ecol. Prog. Ser.; 1996; 144, pp. 265-283. [DOI: https://dx.doi.org/10.3354/meps144265]

66. Lamont, T.; Barlow, R.G.; Morris, T.; van den Berg, M.A. Characterisation of mesoscale features and phytoplankton variability in the Mozambique Channel. Deep Sea Res. Pt. II Top. Stud. Oceanogr.; 2014; 100, pp. 94-105. [DOI: https://dx.doi.org/10.1016/j.dsr2.2013.10.019]

67. Quartly, G.D.; Srokosz, M.A. Eddies in the southern Mozambique Channel. Deep-Sea Res. Pt. Ii; 2004; 51, pp. 69-83. [DOI: https://dx.doi.org/10.1016/j.dsr2.2003.03.001]

68. Weimerskirch, H.; Le Corre, M.; Jaquemet, S.; Potier, M.; Marsac, F. Foraging strategy of a top predator in tropical waters: Great frigatebirds in the Mozambique Channel. Mar. Ecol. Prog. Ser.; 2004; 275, pp. 297-308. [DOI: https://dx.doi.org/10.3354/meps275297]

69. McCabe, R.M.; Hickey, B.M.; Kudela, R.; Lefebvre, K.; Adams, N.G.; Bill, B.D.; Fulland, F.M.D.; Thomson, R.E.; Cochlan, W.; Trainer, V. An unprecedented coastwide toxic algal bloom linked to anomalous ocean conditions. Geophys. Res. Lett.; 2016; 43, pp. 366-376. [DOI: https://dx.doi.org/10.1002/2016GL070023] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27917011]

70. Trainer, V.; Hickey, B.M.; Bates, S.S. Toxic diatoms. Oceans and Human Health; Walsh, P.; Smith, S.; Fleming, L.E.; Solo-Gabriele, H.; Gerwick, W. Academic Press: Cambridge, MA, USA, 2008.

71. Lundholm, N.; Clarke, A.; Ellegaard, M. A 100-year record of changing Pseudo-nitzschia species in a sill-fjord in Denmark related to nitrogen loading and temperature. Harmful Algae; 2010; 9, pp. 449-457. [DOI: https://dx.doi.org/10.1016/j.hal.2010.03.001]

72. Rhodes, L.; Smith, K.; Moisan, C. Shifts and stasis in marine HAB monitoring in New Zealand. Environ. Sci. Pollut. Res. Int.; 2013; 20, pp. 6872-6877. [DOI: https://dx.doi.org/10.1007/s11356-012-0898-9]

73. Smith, J.; Connell, P.; Evans, R.H.; Gellene, A.G.; Howard, M.D.A.; Jones, B.H.; Kaveggia, S.; Palmer, L.; Schnetzer, A.; Seegers, B.N. et al. A decade and a half of Pseudo-nitzschia spp. and domoic acid along the coast of southern California. Harmful Algae; 2018; 79, pp. 87-104. [DOI: https://dx.doi.org/10.1016/j.hal.2018.07.007]

74. Almandoz, G.O.; Ferrario, M.E.; Ferreyra, G.A.; Schloss, I.R.; Esteves, J.L.; Paparazzo, F.E. The genus Pseudo-nitzschia (Bacillariophyceae) in continental shelf waters of Argentina (Southwestern Atlantic Ocean, 38–55° S). Harmful Algae; 2007; 6, pp. 93-103. [DOI: https://dx.doi.org/10.1016/j.hal.2006.07.003]

75. Méndez, S.M.; Ferrario, M.; Cefarelli, A.O. Description of toxigenic species of the genus Pseudo-nitzschia in coastal waters of Uruguay: Morphology and distribution. Harmful Algae; 2012; 19, pp. 53-60. [DOI: https://dx.doi.org/10.1016/j.hal.2012.05.007]

76. Fernandes, L.F.; Brandini, F. The potentially toxic diatom Pseudo-nitzschia H. Peragallo in the Paraná and Santa Catarina States, Southern Brazil. Ineringia Ser. Bot.; 2010; 65, pp. 47-62.

77. Redfield, A.C. Biological control of chemical factors in the environment. Am. Sci.; 1958; 46, pp. 204-221.

78. Bristow, L.A.; Mohr, W.; Ahmerkamp, S.; Kuypers, M.M.M. Nutrients that limit growth in the ocean. Curr. Biol.; 2017; 27, pp. R474-R478. [DOI: https://dx.doi.org/10.1016/j.cub.2017.03.030]

79. Reynolds, C.S. Ecology of Phytoplankton; Cambridge University Press: New York, NY, USA, 2006.

80. Parsons, M.L.; Dortch, Q.; Turner, R.E. Sedimentological evidence of an increase in Pseudo-nitzschia (Bacillariophyceae) abundance in response to coastal eutrophication. Limnol. Oceanogr.; 2002; 47, pp. 551-558. [DOI: https://dx.doi.org/10.4319/lo.2002.47.2.0551]

81. Martiny, A.C.; Pham, C.T.A.; Primeau, F.W.; Vrugt, J.A.; Moore, J.K.; Levin, S.A.; Lomas, M.W. Strong latitudinal patterns in the elemental ratios of marine plankton and organic matter. Nat. Geosci.; 2013; 6, pp. 279-283. [DOI: https://dx.doi.org/10.1038/ngeo1757]

82. Jones, R.H.; Flynn, K.J. Nutritional status and diet composition affect the value of diatoms as copepod prey. Science; 2005; 307, pp. 1457-1459. [DOI: https://dx.doi.org/10.1126/science.1107767] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15746424]

83. Sterner, R.W.; Schultz, K.L. Zooplankton nutrition: Recent progress and a reality check. Aquat. Ecol.; 1998; 32, pp. 261-279. [DOI: https://dx.doi.org/10.1023/A:1009949400573]

84. Cohen, N.R.; Gong, W.; Moran, D.M.; McIlvin, M.R.; Saito, M.A.; Marchetti, A. Transcriptomic and proteomic responses of the oceanic diatom Pseudo-nitzschia granii to iron limitation. Environ. Microbiol.; 2018; 20, pp. 3109-3126. [DOI: https://dx.doi.org/10.1111/1462-2920.14386]

85. Pan, Y.; Bates, S.S.; Cembella, A.D. Environmental stress and domoic acid production by Pseudo-nitzschia: A physiological perspective. Nat. Toxins; 1998; 6, pp. 127-135. [DOI: https://dx.doi.org/10.1002/(SICI)1522-7189(199805/08)6:3/4<127::AID-NT9>3.0.CO;2-2]

86. Hu, Z.-X.; Xu, N.; Aifen, L.; Duan, S. Effects of different N:P ratios on the growth of Pseudo-nitzschia pungens, Prorocentrum donghaiense and Phaeocystis globosa. Acta Hydrobiol. Sin.; 2008; 32, pp. 482-487. [DOI: https://dx.doi.org/10.3724/SP.J.1035.2008.00482]

87. Quilliam, M.A.; Sim, P.G.; McCulloch, A.W.; McInnes, A.G. High-performance liquid chromatography of domoic acid, a marine neurotoxin, with application to shellfish and plankton. Int. J. Environ. Anal. Chem.; 2006; 36, pp. 139-154. [DOI: https://dx.doi.org/10.1080/03067318908026867]

88. Wohlgeschaffen, G.D.; Mann, K.H.; Subba Rao, D.V.; Pocklington, R. Dynamics of the phycotoxin domoic acid: Accumulation and excretion in two commercially important bivalves. J. Appl. Phycol.; 1992; 4, pp. 297-310. [DOI: https://dx.doi.org/10.1007/BF02185786]

89. Petroff, R.; Hendrix, A.; Shum, S.; Grant, K.S.; Lefebvre, K.A.; Burbacher, T.M. Public health risks associated with chronic, low-level domoic acid exposure: A review of the evidence. Pharmacol. Ther.; 2021; 227, 107865. [DOI: https://dx.doi.org/10.1016/j.pharmthera.2021.107865]

90. Ferriss, B.E.; Marcinek, D.J.; Ayres, D.; Borchert, J.; Lefebvre, K.A. Acute and chronic dietary exposure to domoic acid in recreational harvesters: A survey of shellfish consumption behavior. Environ. Int.; 2017; 101, pp. 70-79. [DOI: https://dx.doi.org/10.1016/j.envint.2017.01.006]

91. Hiolski, E.M.; Kendrick, P.S.; Frame, E.R.; Myers, M.S.; Bammler, T.K.; Beyer, R.P.; Farin, F.M.; Wilkerson, H.W.; Smith, D.R.; Marcinek, D.J. et al. Chronic low-level domoic acid exposure alters gene transcription and impairs mitochondrial function in the CNS. Aquat. Toxicol.; 2014; 155, pp. 151-159. [DOI: https://dx.doi.org/10.1016/j.aquatox.2014.06.006]

92. Lefebvre, K.A.; Kendrick, P.S.; Ladiges, W.; Hiolski, E.M.; Ferriss, B.E.; Smith, D.R.; Marcinek, D.J. Chronic low-level exposure to the common seafood toxin domoic acid causes cognitive deficits in mice. Harmful Algae; 2017; 64, pp. 20-29. [DOI: https://dx.doi.org/10.1016/j.hal.2017.03.003] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28427569]

93. Husson, B.; Hernandez-Farinas, T.; Le Gendre, R.; Schapira, M.; Chapelle, A. Two decades of Pseudo-nitzschia spp. blooms and king scallop (Pecten maximus) contamination by domoic acid along the French Atlantic and English Channel coasts: Seasonal dynamics, spatial heterogeneity and interannual variability. Harmful Algae; 2016; 51, pp. 26-39. [DOI: https://dx.doi.org/10.1016/j.hal.2015.10.017] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28003060]

94. Sekula-Wood, E.; Benitez-Nelson, C.; Morton, S.; Anderson, C.; Burrell, C.; Thunell, R. Pseudo-nitzschia and domoic acid fluxes in Santa Barbara Basin (CA) from 1993 to 2008. Harmful Algae; 2011; 10, pp. 567-575. [DOI: https://dx.doi.org/10.1016/j.hal.2011.04.009]

95. Lugomela, C. Population dynamics of Pseudo-nitzschia species (Bacillariophyceae) in the near shore waters of Dar Es Salaam, Tanzania. Tanz. J. Sci.; 2013; 39, pp. 38-48.

96. Hermes, J.C.; Masumoto, Y.; Beal, L.M.; Roxy, M.K.; Vialard, J.; Andres, M.; Annamalai, H.; Behera, S.; D’Adamo, N.; Doi, T. et al. A sustained ocean observing system in the Indian Ocean for climate related scientific knowledge and societal needs. Front. Mar. Sci.; 2019; 6, 355. [DOI: https://dx.doi.org/10.3389/fmars.2019.00355]

97. Lugomela, C.; Lyimo, T.J.; Bryceson, I.; Semesi, A.K.; Bergman, B. Trichodesmium in coastal waters of Tanzania: Diversity, seasonality, nitrogen and carbon fixation. Hydrobiologia; 2002; 477, pp. 1-13. [DOI: https://dx.doi.org/10.1023/A:1021017125376]

98. Backeberg, B.C.; Reason, C.J.C. A connection between the South Equatorial Current north of Madagascar and Mozambique Channel Eddies. Geophys. Res. Lett.; 2010; 37, L04604. [DOI: https://dx.doi.org/10.1029/2009GL041950]

99. Madhu, N.V.; Maheswaran, P.A.; Jyothibabu, R.; Sunil, V.; Revichandran, C.; Balasubramanian, T.; Gopalakrishnan, T.C.; Nair, K.K.C. Enhanced biological production off Chennai triggered by October 1999 super cyclone (Orissa). Curr. Sci.; 2002; 82, pp. 1472-1479.

100. Uz, B.M. What causes the sporadic phytoplankton bloom southeast of Madagascar?. J. Geophys. Res.; 2007; 112, 1129. [DOI: https://dx.doi.org/10.1029/2006JC003685]

101. Hamisi, M.I.; Mamboya, F.A. Nutrient and phytoplankton dynamics along the Ocean Road Sewage Discharge Channel, Dar es Salaam, Tanzania. J. Ecosys.; 2014; 2014, 271456. [DOI: https://dx.doi.org/10.1155/2014/271456]

102. Saji, N.H.; Goswami, B.N.; Vinayachandran, P.N.; Yamagata, T. A dipole mode in the tropical Indian Ocean. Nature; 1999; 401, pp. 360-363. [DOI: https://dx.doi.org/10.1038/43854]

103. Usman, M.T.; Reason, C.J.C. Dry spell frequencies and their variability over southern Africa. Climate Res.; 2004; 26, pp. 199-211. [DOI: https://dx.doi.org/10.3354/cr026199]

104. Porter, A.; Bohl, D.; Kwasi, S.; Donnenfeld, Z.; Cilliers, J. Prospects and Challenges: Mozambique’s Growth and Human Development Outlook to 2040; Institute for Security Studies,: Pretoria, South Africa, 2017.

105. FAO. National Aquaculture Sector Overview: Mozambique; FAO: Rome, Italy, 2005.

106. Bryden, H.L.; Beal, L.M.; Duncan, L.M. Structure and transport of the Agulhas Current and its temporal variability. J. Oceanogr.; 2005; 61, pp. 479-492. [DOI: https://dx.doi.org/10.1007/s10872-005-0057-8]

107. Lamont, T.; Roberts, M.J.; Barlow, R.G.; Morris, T.; van den Berg, M.A. Circulation patterns in the Delagoa Bight, Mozambique, and the influence of deep ocean eddies. Afr. J. Mar. Sci.; 2010; 32, pp. 553-562. [DOI: https://dx.doi.org/10.2989/1814232X.2010.538147]

108. Lutjeharms, J.R.E. The ocean environment off southeastern Africa: A review. S. Afr. J. Sci.; 2006; 102, pp. 419-426.

109. Swart, N.C.; Lutjeharms, J.R.E.; Ridderinkhof, H.; de Ruijter, W.P.M. Observed characteristics of Mozambique Channel eddies. J. Geophys. Res.; 2010; 115, C09006. [DOI: https://dx.doi.org/10.1029/2009JC005875]

110. Westberry, T.K.; Siegel, D.A. Spatial and temporal distribution of Trichodesmium blooms in the world’s oceans. Glob. Biogeochem. Cy.; 2006; 20, GB4016. [DOI: https://dx.doi.org/10.1029/2005GB002673]

111. Goreau, T.F.; Goreau, N.I.; Goreau, T.J. Corals and coral reefs. Sci. Am.; 1979; 241, pp. 124-137. [DOI: https://dx.doi.org/10.1038/scientificamerican0879-124]

112. Bakun, A. Coastal upwelling indices, west coast of North America. Spec. Sci. Rep. Fish. Ser.; 1973; 231, pp. 1946-1971.

113. Gonzalez-Nuevo, G.; Gago, J.; Cabanas, J.M. Upwelling index: A powerful tool for marine research in the NW Iberian upwelling system. J. Oper. Oceano.; 2014; 7, pp. 47-57. [DOI: https://dx.doi.org/10.1080/1755876X.2014.11020152]

114. Lavín, A.; Díaz del Río, G.; Cabanas, J.M.; Casas, G. Afloramiento en el Noroeste de la Peninsula Iberica Indices de Afloramiento Para El Punto 43.0N–11.0W Periodo 1966-1989; Instituto Español de Oceanografia: Madrid, Spain, 1991.

115. Utermöhl, H. Toward the improvement of the quantitative phytoplankton method. Mitt. Int. Verieningung Fur Limnol.; 1958; 9, pp. 1-39.

116. Bianchi, T.S.; Bauer, J.E. Particulate organic carbon cycling and transformation. Treatise on Estuarine and Coastal Science; Academic Press: Cambridge, MA, USA, 2011; pp. 69-117.

117. Van Heukelem, L.; Lewitus, A.J.; Kana, T.M.; Craft, N.E. Improved separations of phytoplankton pigments using temperature-controlled high performance liquid chromatography. Mar. Ecol. Prog. Ser.; 1994; 114, pp. 303-313. [DOI: https://dx.doi.org/10.3354/meps114303]

118. Latasa, M. Improving estimations of phytoplankton class abundances using CHEMTAX. Mar. Ecol. Prog. Ser.; 2007; 329, pp. 13-21. [DOI: https://dx.doi.org/10.3354/meps329013]

119. Higgins, H.W.; Wright, S.W.; Schluter, L. Quantitative interpretation of chemotaxonomic pigment data. Phytoplankton Pigments: Characterization, Chemotaxonomy, and Applications in Oceanography; Roy, S.; Llewellyn, C.A. Cambridge University Press: New York, NY, USA, 2011; pp. 257-313.

120. Van der Molen, J.S.; Scharler, U.M.; Muir, D. Species composition, abundance and biomass of microphytoplankton in the KwaZulu-Natal Bight on the east coast of South Africa. Afr. J. Mar. Sci.; 2016; 38, pp. S139-S153. [DOI: https://dx.doi.org/10.2989/1814232X.2016.1145597]

121. Harris, R.P.; Wiebe, P.H.; Lenz, J.; Skjoldal, H.R.; Huntly, M. Zooplankton Methodology Manual; Academic Press: Cambridge, MA, USA, 2000.

122. Baustian, M.M.; Bargu, S.; Morrison, W.; Sexton, C.; Rabalais, N.N. The polychaete, Paraprionospio pinnata, is a likely vector of domoic acid to the benthic food web in the northern Gulf of Mexico. Harmful Algae; 2018; 79, pp. 44-49. [DOI: https://dx.doi.org/10.1016/j.hal.2018.06.002] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30420014]