Research Article |
Corresponding author: Aylin Ulman ( merseamed@gmail.com ) Academic editor: Cascade Sorte
© 2021 Aylin Ulman, Taner Yildiz, Nazli Demirel, Ozgur Canak, Emre Yemişken, Daniel Pauly.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Ulman A, Yildiz T, Demirel N, Canak O, Yemişken E, Pauly D (2021) The biology and ecology of the invasive silver-cheeked toadfish (Lagocephalus sceleratus), with emphasis on the Eastern Mediterranean. NeoBiota 68: 145-175. https://doi.org/10.3897/neobiota.68.71767
|
Invasive species pose threats to either human health or inflict ecological and/or economic damage. The silver-cheeked toadfish (Lagocephalus sceleratus), a Lessepsian species, is one of the most harmful species in the Mediterranean Sea, because of its potent neurotoxin, impacts on marine biodiversity, and the increased costs and labor they inflict on fishers. Since the catch and consumption of this pufferfish is prohibited by almost all countries bordering the Mediterranean, they have now expanded into the entire Mediterranean and Black Sea. We performed a comprehensive study of L. sceleratus covering ecological aspects, growth, reproduction, diet and trophic level based on samples from southwestern coasts of Turkey. The estimated growth parameters were L∞ = 88.7 cm, K = 0.27 year-1, C = 0.6 and WP = 0.1. Their sex-ratio was M/F = 1:0.69. Lagocephalus sceleratus appears to be a batch spawner with discontinuous oocyte recruitment and has different spawning seasons in the Eastern Mediterranean which seem to be based on temperature cues which get shorter in duration as one moves north from the Suez. We also report their first positive ecological trait, that they are controlling some other invasive species through their diets, such as lionfish, Red Sea goatfish, rabbitfish and longspine sea urchins, in addition to controlling themselves through cannibalism, which appears to be density-dependent. They are indeed a top predator in the region with a trophic level of 4.1. We suggest that targeted fishing using improved gear-types to reduce fishing gear damages are initiated, and that finding commercial markets for pufferfish could help to naturally fund ongoing control efforts.
Cannibalism, growth, Invasive Alien Species (IAS), pufferfish, reproduction, Tetraodontidae
Global biodiversity is currently being threatened by overfishing, pollution and invasive species (
Non-indigenous species (NIS) are called invasive when they cause either ecological, economical damage, or pose a threat to human health. Marine invasive species can pose major threats to biodiversity by altering community structure and function, and by modifying ecosystem processes, which can have long-lasting ecological and economic consequences (
There are 197 species of pufferfish globally, 112 of which live in marine environments, 48 in brackish environments and 37 in freshwater (
Pufferfish species found in the Mediterranean, their native region, first year of introduction, first locality introduced, established status in the Mediterranean, and reported in Turkey (Y= yes, N= no).
Common name | Scientific name | Native region | Year | Locality | Established | In Turkey |
---|---|---|---|---|---|---|
Prickly puffer | Ephippion guttiferum | E. Atlantic & W. Med. | NA | NA | Y | N |
Oceanic puffer | Lagocephalus lagocephalus | Subtropical | NA | NA | Y | Y |
Diamondback puffer | Lagocephalus guentheri | Indo-Pacific | 19501 | Egypt | Y | Y |
Suez puffer | Lagocephalus suezensis | W. Indian, Red Sea | 1977 | Lebanon | Y | Y |
Silver-cheeked toadfish | Lagocephalus sceleratus | Indo-Pacific | 2003 | Turkey | Y | Y |
Guinean puffer | Sphoeroides marmoratus | E. Atlantic | 1977 | Italy | Y | N |
Blunthead puffer | Sphoeroides pachygaster | Subtropical | 1979 | Spain | Y | Y |
Bandtail puffer | Sphoeroides spengleri | W. Atlantic | 2000 | Spain | Y | N |
Yellowspotted puffer | Torquigener flavimaculosus | W. Indian | 1987 | Israel | Y | Y |
Spiny blaasop | Tylerius spinosissimus | Indo-Pacific | 2004 | Greece | Y | Y |
Spotbase burrfish | Cyclichthys spilostylus | Indo-Pacific | 1993 | Israel | Y | Y |
Spotfin burrfish | Chilomycterus reticulatus | Subtropical | 2009 | Sardinia | N | N |
Spotfin porcupinefish | Diodon hysterix | Circumtropical | 1956, 2016 | Italy, Spain | N | N |
The silver-cheeked toadfish Lagocephalus sceleratus (Gmelin, 1789) A lateral view B ‘puffed’ lateral view (original drawings by Marc Dando).
Lagocephalus sceleratus (Fig.
An important part of the ecological sucess of L. sceleratus is due to their having one of the most advanced forms of teeth in the animal kingdom. The ‘first generation teeth’ are coated with recurring toothbands which are continously regenerated by stem cells (
The first record of Lagocephalus sceleratus in the Mediterranean was from Gökova Bay, southwestern Turkey in 2003 (
Distribution of Lagocephalus sceleratus in the Mediterranean, as documented by records in
In Lagocephalus sceleratus’s 18-year presence in the Mediterranean, only loggerhead turtles (Caretta caretta) have been documented preying on adult L. sceleratus, while garfish (Belone belone), common dolphinfish (Coryphaena hippurus), and cannibalism has been documented in juveniles (
Lagocephalus sceleratus combine two exceptional defense mechanisms which benefit them in evading predators, i.e., the ability to ‘puff’ themselves up, and their highly toxic tissues. The combination of these two factors contribute, in the Mediterranean, to a scarcity of predators.
The peculiar head, buccal cavity and pectoral girdle structures of these fishes facilitate their unique ability to ‘puff’ themselves up by rhythmic buccal pumping, swallowing and forcing water (or air if they are outside of water) into their stomach. While their ‘stomach’ can increase its size 50–100-fold depending on the species (
Lagocephalus sceleratus is the second most poisonous Mediterranean pufferfish species after Torquigener flavimaculosus Hardy & Randall, 1983 (Ayas 2017;
Tetrodotoxin (TTX) is an extremely potent neurotoxin found in L. sceleratus and other pufferfish (
Out of the 197 pufferfish species, only 55 (28%) are considered toxic (
Lagocephalus sceleratus has strong negative impacts on the livelihoods of small-scale fishers of the Eastern Mediterranean, most of whom are already marginalized due to declining catches and revenues (
Interestingly, around Turkey, this species is normally shy of humans and is not commonly encountered while snorkeling or scuba diving. In August 2019, a first human attack by L. sceleratus occurred in Kaledran, Turkey where L. sceleratus bit a child three times on the left hand, resulting in the amputation of her ring finger (
Lagocephalus sceleratus poisoning has caused dozens of human fatalities in the Mediterranean region, which is a severe underestimate given that most of these fatalities are not officially recorded (
Despite the multiple negative impacts of L. sceleratus, most Mediterranean studies have been dedicated to the high content of tetrodotoxin (TTX) in its tissues, with only a handful investigating its biology (Sabrah 2006;
This contribution is an attempt to correct this imbalance. Due to nearly a complete lack of control in the region, its negative impacts to marine biodiversity, human health and fishers’ livelihoods continue to worsen. This study presents new data on the species behaviour (eg., spawning, ecology and feeding) based on fishers’ knowledge, and from biological studies, presenting new data on their distribution, size, growth, spawning season and reproductive status, reproductive morphology and fecundity, and the taxonomic composition of their prey. This contribution aims to improve current knowledge about this invasive species, to help direct further research needs and management options.
Pufferfish samples were purchased from small-scale commercial fishers in southwestern Turkey primarily from Datça where they were targeted (36.726°N, 27.685°E) and about 15% of samples were caught as by-catch from Fethiye (36.659°N, 29.126°E), both Muğla Province, Turkey, from June 2019 to November 2020. This stock has not yet been studied and is understood to be a different stock from the neighbouring Antalya province, which has been somewhat studied. This area is very close to Gökova Bay, where the first Mediterranean L. sceleratus occurrence was reported (
A total of 1013 fish: 456 males, 270 females and 287 juveniles (where juveniles were generally < 25 cm and could not have their sex determined due a lack of gonadal development) were collected for this study from June 2019 to November 2020. Fishers were paid 10 Turkish Lira (≈ US$ 1.20; April 23/2021) per kg for L. sceleratus from June 2019 until mid-April 2020, and 20 Turkish Lira per kg from mid-April 2020 onwards. The fish were purchased from approximately 20 fishers from Fethiye and Datça, Muğla province, who all had special permissions to collect them for this study. Permission to collect pufferfish for the specified designated fishers for scientific research purposes was granted from the Turkish Ministry of Agriculture and Forestry and General Directorate of Water Products under Permission #67852565-140.03.03-E.1354602 & #6987137-663.08.
We formally surveyed 45 small-scale fishers face-to-face from the Muğla province (Fethiye to Bodrum) in April 2019 to help understand some of the behavior of this species and to inform them of this study. An initial structured survey consisting of 18 questions pertaining to their contact details, fisher characteristics, vessel and geartypes, average days fished, L. sceleratus catches, catch areas, caught depths, average sizes, maximum sizes, fishing gear losses in nets and longline hooks, and interest in catching pufferfish for this study was initially undertaken at the beginning of the study in April 2019. Twelve of those initially interviewed supplied fish afterwards for this study all using trammel nets, after permissions were granted for them to catch pufferfish. Any new information learnt as the study progressed was written down and transferred to the spreadsheets containing the other data. These data were then summarized for each topic. Their responses, aside from the new maximum depth record, should be viewed as anecdotal evidence.
Information on such basic biological parameters of species, such as growth, reproduction and fecundity are essential in understanding the basic life history traits of a species and are prerequisites needed to develop scientifically sound fisheries management policies. For all 1013 samples, the total length (L) and body weight (W) of fish were measured to the nearest 0.1 cm and the nearest 1 g, respectively, and gonads and livers were weighed to the nearest 0.01 g. The length-weight relationship yields authentic biological information about a species in a particular region and is of great importance in fishery assessments. The parameter of length-weight relationships (LWRs) of the form W = a·Lb were estimated through re-expression of the LWR equations in linearized form, i.e., log (W) = log (a) + b·log (L), where a is a scaling coefficient for the weight at length and b is a shape parameter; note that if b < 3, a fish become thinner as it grows, and plumper if b > 3.
The growth of water-breathing ectotherms such as fish can be conceived as the net result of two processes with opposing tendencies (
dW ∕ dt = HWd − kW (1)
where dW/dt is the growth rate, W is body weight (or mass), H and k are the coefficients of anabolism and catabolism, and d is the scaling exponents of anabolism, which depend on oxygen, and hence of the growth of gill surface area (Pauly 1984,
Lt = L∞ (1-e-K(t-t0)) (2)
where Lt is the length at age t, L∞ is the asymptotic length, i.e. the mean length the individuals of a given population would reach if they grew indefinitely, K is rate, or dimension time-1 (here: year-1) at which L∞ is approached, and t0 is the age at L = 0.
The mutual compatibility of the growth parameters L∞ and K can be evaluated by Ø’ = log(K)+2log(L∞) which should be roughly similar between populations of the same species (
Here, a seasonally oscillating variant of the von Bertalanffy growth function (VBGF) was used to estimate growth parameters from the length-frequency data available for L. sceleratus; this version of the VBGF has the form:
Lt = L∞{1-e-[K(t – t0)+S(t) – S(t0)]} (3)
where S (t) = (CK/2π)·sin (2π (t − ts)), S (t0) = (CK/2π)·sin (2π (t0 − ts), and L∞, K and t0 are defined as above; see
Equation (3) involves two parameters more than the standard VBGF: C and ts. Of these, the former is easiest to visualize, as it expresses the amplitude of the growth oscillations. When C = 0, Equation (3) reverts to Equation (2). When C = 0.5, the seasonal growth oscillations are such that growth rate increases by 50% at the peak of the ‘growth season’ (i.e., in ‘summer’), and, briefly, declines by 50% in ‘winter’. When C = 1, growth increases by 100%, doubling during ‘summer’, and becoming zero in the depth of ‘winter’. The other new parameter, ts expresses the time elapsed between t = 0 and the start of a sinusoid growth oscillation. However, visualization is facilitated if we define ts + 0.5 = WP (‘Winter Point’), which expresses, as a fraction of the year, the period when growth is slowest. WP is often close to 0.1 (i.e., early February) in the Northern Hemisphere and 0.6 (early August) in the Southern Hemisphere.
The parameters of Equation 3 were estimated through the ELEFAN method, which fits growth to the peaks of length-frequency (L/F) samples arranged in time (represented by black, positive histograms, and deemed to represent age classes) while avoiding the trough between peaks (represented by white, negative histograms). Peaks and troughs are identified by a simple high-pass filter, i.e., a running average which leads to definition of peaks as those parts of a length-frequency distribution that are above the corresponding running average and conversely for the troughs separating peaks. Then, hundreds of growth curves, each with a different set of growth parameters, are traced, and the growth curve (i.e., parameter set) is retained which has the highest score in linking the peaks of L/F distributions, whose ‘point’ values are positive, while avoiding troughs, whose point values are negative (
Variations in fish gonadal morphology explain important behavioral and ecological adaptations during reproduction. Particularly knowledge about the reproductive period is considered a major life-history trait and evaluating the changes in gonadal development, liver size and body weight can help to understand energy trade-offs in the development of reproductive strategies, notably in the inverse relationship between the gonadosomatic index (GSI) and the hepato-somatic index (HSI), while condition factor (CF) shows the relative health of the fish.
To estimate fecundity, the gonads were removed, weighed and preserved in formalin. To identify the reproductive season, temporal changes in the gonadosomatic index were assessed using the relation: GSI = 100·× [GW/(TW − GW)] where GW is the gonad weight and TW is the total weight. Also, the hepato-somatic index analyses was computed as an indicator of reserves in the liver, i.e., HSI = 100·× [HW/(TW − HW)] where HW and TW represent liver weight and total weight, respectively. Understanding changes in liver reserves, helps to better understand how energy is transferred from storage to reproduction. Finally, the overall plumpness of individuals was determined from their condition factor CF = 100·W/L3.
The size at first maturity (and spawning) was estimated by plotting the fraction of mature individual females and males against their lengths, and fitting a logistic curve. Mean length at first maturity (Lm) was the length at which, in a given population, 50% of individuals were mature. This was evaluated separately for fish sampled during the main spawning season (i.e., in June) and outside, to test if L. sceleratus reach maturity at smaller sizes within than outside the spawning season.
We also used the lengths of first maturity and maximum lengths in several population of L. sceleratus to indirectly estimate their ratios of metabolic rate at first maturity (Qm) to maintenance rates (Qmaint). These ratios were then used to test whether their mean value is compatible with earlier estimate ranging from 1.22 to 1.53 and suggesting that it is a declining relative oxygen supply which triggers maturation and spawning (
Knowledge on fecundity is used to calculate the reproductive potential of a stock and is another important factor for effective fish stock management. Ovary samples were collected in May and June 2020, to capture the peak GSI values. The oocyte size–frequency method (
Examination of oocyte development is evaluated to help identify reproductive strategies of species such as ovary organization, fecundity type and spawning patterns (
Knowledge on predator-prey interactions for species are essential to understanding their role in the ecosystem, impacts on biodiversity, and are essential in building accurate ecosystem models for a region. Two complementary studies were conducted on the diet of L. sceleratus. The prey/diet preferences were examined by a visual taxonomic examination of stomach content for 563 samples from Fethiye and Datça, Muğla province in Turkey. Food items were removed from the esophagus, stomach and intestine and identified to the lowest taxonomic level possible; fishing hooks and pieces of fishing net were also accounted for, as were sand and algae. The prey taxa were then grouped into three main categories: crustaceans, fish and cephalopods, and also identified as either indigenous or non-indigenous taxa where possible. A t-test was performed on the ratios of the three prey groups for juvenile (< 45 cm) and adult fish (> 45 cm) to determine if they target different taxonomic groups as they grow.
To better understand the role of L. sceleratus in the ecosystem, and to estimate their trophic level (TL), their mean fractional level of their prey for 34 stomachs, where the contribution of prey items in numbers (%N), weight (%W) and frequency of occurrence (%F) was recorded. These values were then used for calculating the Index of Relative Importance (IRI) of prey item (IRI = %F × (%N + %W)), which was then re-expressed using %IRI = (IRI/ SIRI) × 100 (Cortes 1997). SIRI is the percentage which a discrete prey taxon contributes to the sum of all IRI values in the prey spectrum. Based on the dietary composition (expressed as W%), the mean fractional trophic level (TP) of the L. sceleratus was estimated using the method of
The fishers who informed this study consisted of 12 using trammel nets, 12 using longlines, 21 using both trammel nets and longlines, and five occasionally using rods. The fishers who provided fish for this study used trammel nets, with three sometimes using fishing rods.
According to these fishers, when L. sceleratus first appeared along the southwestern Turkish coast, it was found mostly in rocky areas from depths of about 10 m, and never deeper than 100 m. However, over time L. sceleratus were increasingly found in deeper locations to a maximum of 220 m depth (recorded in April 2021 from Fethiye Bay). In June, i.e., during their spawning season, they aggregate in the shallows of bays, between 5–10 m depths; however, a few individuals have also been caught at the surface. Based on the accounts of 45 fishers in Muğla Province, Turkey, L. sceleratus regularly consumes bait from rods and longlines, severing many of the hooks and even steel lines in the process. Some fishers reported hook losses from 50–90% of their longlines in extreme cases, but the majority of long-line fishers claimed an average of about 10–20% of hooks lost. Hooks were found in 8% of L. sceleratus stomachs; nearly all samples were collected by net. Lagocephalus sceleratus uses its fused parrot-like teeth to bite holes in set trammel nets and consume the fish caught in the nets, as evidenced by nine pieces of fishing net between 3–20 cm in diameter in their stomachs, weighing up to 10 g. All fishers in the region are regularly affected by this and try to cast their nets away from L. sceleratus hotspots to minimize damages. One fisher from Fethiye (Mehmet Taniş, pers. comm.) explained that on several occasions, L. sceleratus bit through his trammel nets, and consumed the stingrays caught inside, leaving only the needle tail portion behind as evidence.
As one fisher, S. Taşkıran, was the main fisher in Datça that targeted L. sceleratus for this study, and thus has the most experience with this species, his observations are separately noted here. He estimated that in June 2020, there were approximately 10 tonnes of L. sceleratus spawning in InçiBurnu Bay near Datça. At this locality, during their spawning period, the fish were inactive at night, and actively fed at dusk and dawn. In July and August, they fed very little, but in September onwards for a few months, they again fed very aggressively.
In total, 1013 fish were examined, and of those, 456 were male, 270 were female and 287 were juveniles generally below 25 cm whose sex could not be determined. The overall sex-ratio was calculated as M:F = 1.0:0.69. Total length ranged from 13 to 77.2 cm. The mean lengths of females and males were not statistically different (p > 0.05, p = 0.71) but males were more abundant throughout the entire year. Suppl. material
Suppl. material
Illustrating the relationship between the multiplicative term (a) and the exponent (b) of length-weight relationship in Lagocephalus sceleratus. (Based on the results of 13 studies from data in Suppl. material
The close inverse relationship of log(a) vs b in Fig.
The best fit to the length-frequency data that we gathered (see Fig.
Seasonally oscillating growth curve fitted using ELEFAN to 14 length-frequency samples of Lagocephalus sceleratus (n = 1013) collected from June 2019 to November 2020; the estimated growth parameters were L∞ = 88.7 cm (TL), K = 0.27 year-1, C = 0.6 and WP = 0.1. The Roman numerals (& 0) refer to the 7 cohorts (= sequences of peak, i.e., black histograms) that were identified by ELEFAN (Based on data in Suppl. material
The ovarian organization of L. sceleratus appears to be based on synchronous development of groups of oocytes. Two concurrent populations of oocytes were found during spawning, i.e., larger oocytes and a more heterogeneous group of smaller oocytes (Fig.
Stages of oocyte development in Lagocephalus sceleratus. Whole oocytes on the slide A and histological sections B–E from nucleolus to vitellogenesis; and F hydration: Primary growth (pg), cortical alveoli (ca), nucleus (N), vitellogenic oocytes (Vit), oil droplets (od), atretic oocyte (A), and hydrated oocytes (H). Scale bars 1 mm (A); 200 μm (B–E); 400 μm (F).
Oocyte diameter during vitellogenesis were found to range between 0.42–0.58 mm, with an average oocyte size of 0.50 mm for the migratory nucleus stages. Oocyte counts were performed on 23 female ovaries from the peak reproductive period. Average fecundity was calculated as 134,000 oocytes for females of 55 cm and 2,000 g. The relationships between fecundity vs. length, and oocyte number vs. body weight are provided in Suppl. material
GSI starts to increase in April and May, peaks in June (9%), then declines sharply in July (see Fig.
Seasonal variation of the Gonadosomatic Index (GSI) of Lagocephalus sceleratus in the Mediterranean and the Suez region, based on data by Sabrah 2006, (1st location at the bottom of the figure), Syria-Leb/Khalaf 2014 (2nd location), Lebanon/Boustany 2015 (3rd location), S. Cyprus/Rousou, 2014 (4th location), N. Cyprus/Akbora 2020 (5th location), Antalya Bay/
Condition factors were similar between sexes, and its monthly variability (not shown) was not very pronounced; it exhibited a weak peak in June (during peak spawning season) and another in November. The baseline of the HSI index was around 3–4%. The HSI index started to increase in November to peak at 8% in April, thus suggesting that reserves were taken from the liver to be used for gonadal development.
As in other fish species, observed maturity stages in L. sceleratus were a function of size (Figure
Of the 563 fish that had their stomach contents examined, 48 (9%) of the stomachs were empty, 58 (10%) had food residues, 253 (45%) had stomachs less than half full, 170 (30%) were over half full, and 34 (6%) were full. A total of 34 specimens (Suppl. material
Stomach contents of Lagocephalus sceleratus presenting evidence of cannibalism (A) and predation on other invasive species, i.e., Torquigener flavimaculosus (B); and Parapeneus forsskali (C).
Crustaceans and fish made up the majority of diets being found in 26% and 24% of stomachs, respectively, with cephalopod remains in 11%. There was no statistical difference between the taxonomic prey composition between juvenile and adult L. sceleratus ratios of crustaceans, fish and cephalopods (p = 0.225). For crustaceans, small crabs were the major taxon, with only a few stomachs containing shrimp remains, as expected, since crab shells take longer to be digested and/or evacuated. Of the crabs, Carappa granulata was found in two stomachs, one Carcinus aestuarii, one Charybdis sp., and one Scyliarides latus. For fish, those that could be identified were three Pterois miles, Scorpaena spp., Epinephelus spp., Mugilidae spp., Atherina spp., Diplodus sp., Sparus aurata and Siganus spp. Of cephalopods, there were about a dozen cases each of common squid (Loligo vulgaris), common octopus (Octopus vulgaris), one violet blanket octopus (Tremoctopus violaceus) and unidentified cephalopod beaks and ink (Suppl. material
From the IRI examination, 34 additional L. sceleratus stomachs were analysed, and 91% (31) of those had food in their stomachs (coefficient of vacuity: 23.5). The IRI results are presented in Suppl. material
L. sceleratus offers a trifecta of highly negative impacts due to its high toxicity, economic losses to fishers, and negative effects on native marine biodiversity. Their unique ability to puff and high toxicity likely contribute to their invasive success in the Mediterranean. Due to a nearly complete lack of population control, this species has expanded to all corners of the basin, putting people, fishers, fisheries and the native ecology at risk. Its conquest of the Mediterranean is one of the most successful marine invasions in modern history, comparable with that of the invasive Western Atlantic lionfish Pterois volitans and Pterois miles, the latter having also established itself as a Lessepsian species in the Mediterranean in 2012 (
Their marked expansion benefits from both a lack of human control (as fishing and sales of pufferfish are prohibited in most countries, including all countries of the EU), and limited predatory control (due to their ability to puff, and high TTX content). Their success is also likely enhanced by the overfishing that characterizes the Mediterranean basin, which has lost its top predators (
Here we reveal the results of biological studies on their morphometrics and growth, reproduction, and diet before presenting some management advice and ideas for further directed research. The morphometric (LWRs) and growth studies conducted here produced results that were comparable to those of other authors. This also included the ratio LmaxD /LmD, which was statistically undistinguishable from estimates of this ratio in other teleosts (Pauly 1984;
Using traditional biological sampling combined with fishers’ knowledge improved the biological understanding of L. sceleratus, e.g, their spawning periodicity. We also found the HSI and GSI patterns to be asynchronous, which explains how its energy is stored and utilized (
From its ovarian organization, L. sceleratus was identified as a group synchronous batch spawner from the presence of both previtellogenic oocytes (in a range of sizes) and larger vitellogenic oocytes (of larger similar sizes) in the ovaries during the peak spawning period. The presence of these two clearly different size groups of oocytes is defined as group asynchronous ovarian organization, with a heterogeneous population of oocytes in their primary growth stage together with a synchronous population of larger oocytes in the yolked stage, indicating further recruitment into the oocyte stock at any time during the spawning season (
The most important finding of this study is that L. sceleratus appears to prey on a wide range of other invasive species, and its control of them is its first positive documented ecological trait. The spines of lionfish (P. miles) found inside three L. sceleratus suggest that pufferfish are preying on lionfish. L. sceleratus also provide some control on other invasive species such as Red Sea goatfish, rabbitfish, other pufferfish species, their own species, and even the longspine sea urchin (Diadema setosum). The finding that their target prey composition is nearly equally comprised of fish and crustaceans, and a lesser extent of cephalopods, did not differ between juveniles and adults which contrasts the earlier findings of
Lagocephalus sceleratus are now in direct competition with small-scale fishers in the Eastern Mediterranean, consuming their catches, revenue, time and thus much of their livelihoods. Their increasing damage to fishing gear also negatively impacts their livelihoods. In Cyprus, fishers often use newspaper articles mentioning pufferfish damage to lobby for financial support. In the Muğla province of Turkey, Ünal (2013) found over 90% of small-scale fishers were no longer generating a net income from their work. Small-scale fishers from southern Turkey, already highly marginalized, with many being forced out of the profession due to declining catches and incomes, have to completely replace their fishing nets every few months at an added cost of over $2000 US due to pufferfish damage, which previously lasted them several years. As both an incentive with the benefit of aiding fishers offset the increasing costs from pufferfish, a bounty program was recently initiated in Turkey.
The Turkish government recently completed a pilot bounty project collecting L. sceleratus tails from the Turkish Mediterranean coast in December 2020. A total of 46,000 tails were collected for a reward price of US $0.60 each. A second bounty program was established on June 27, 2021 for a duration of three years, during which L. sceleratus will be, this time, collected in its entirety (@ US $0.60 each) so that proper disposal can be ensured (Mahir Kanyılmaz, Fisheries Directorate, Ankara, Turkey, pers. comm.). Even if this initiative is not effective at reducing abundances, it will still add some positive economic benefit to some fishers. To predict how much the population of L. sceleratus should be reduced to negate its impacts to native biodiversity through predation, its biomass, the biomass or abundance of its prey and its feeding rates must first be known. However, only one stock assessment from a small area has been completed in Turkey (
After discussing the bounty program with twenty small-scale fishers from the Muğla province of Turkey, we strongly believe that this new bounty will not be effective at reducing their population enough to negate their effects. At present, small-scale fishers refuse to target this species due to the low reward and high costs of fishing gear damage; however, large-scale fishers may return specimens for reward if many are caught in a net at once. Alternatively, to control this species, we suggest a better solution would be to hire select commercial fishers, equip them with more resistant fishing nets, and have them specifically target L. sceleratus in their spawning season where they tend to aggregate. One Turkish fishing gear technology expert, Dr. Zafer Tosunoğlu, Ege University, who was contacted for advice on the most applicable net material to specifically target L. sceleratus suggested using Dyneema netting (used for catamaran trampoline netting), which is the strongest netting fibre currently available, 15 times stronger than steel, and should minimize fishing gear damage. Also, emerging ‘genetic biocontrol’ may be applied (
The development of commercial applications for invasive species such as L. sceleratus may financially support their ongoing removal (
This NIS top predator in the Mediterranean threatens local biodiversity, human health, fishing communities and potentially even tourism. Since it is currently lacking control on its population in most of the Mediterranean Sea, we suggest that removal through targeted fishing during its spawning period is the best control recommendation for decision-makers. Removals needs to be prioritized but can be expensive to fund, which is why commercialization of this species could help financially sustain their long-term control. The current commercial solutions, which would use a highly invasive species to benefit our teeth and skin, heal our wounds, alleviate our pain and protect our feet, if successful, could represent the largest turnaround in the history of marine invasions.
Prior to this study, most Mediterranean research on this species either studied its growth, or toxicity. This study added to this body of knowledge by determining the spawning strategy and reproductive ecology of L. sceleratus, factors relating to its growth, and its position within the trophic web, its density-dependent cannibalistic nature, and its potential ability to help control subsequent invasions. The diversity of their prey can be used as inputs for ecosystem modeling efforts, which, along with improved biomass estimates, can help to understand how much should be removed to help improve the state of native biodiversity. Further directed research needed to better understand and hence manage this invasion should involve mapping its various spawning habitats and seasons, its larval ecology and growth, its feeding rates, DNA stomach content analysis, a forum to update on their interactions with humans, and baseline stock assessments along the entire Levantine coast.
We thank the Turkish Ministry of Agriculture and Forestry for providing permissions to collect pufferfish from select fishers. We are grateful to Marc Dando for his incredible scientific illustrations, and we thank Dr. Alp Salman for aiding with the identification of the blanket octopus. We also thank friendly fishers for their assistance, especially Mr. Serhat Taşkiran for catching the majority of pufferfish we studied. Also, we thank Mr. Ioannis Giovos and Burak Çiçek for helping with the identification of stomach contents, and Ms. Elaine Chu for drafting Figures
Appendix 1, Figures S1–S4, Tables S1–S8
Data type: docx. file
Explanation note: Appendix 1. The Gill-Oxygen Limitation Theory. Figure S1. Plot of LmaxD vs. LmD (for D = 0.6) in different populations of L. sceleratus. Figure S2. Length-weight relationships for A) female and B) male Lagocephalus sceleratus. Figure S3. Relationships, in Lagocephalus sceleratus between fecundity and total length, total weight and gonad weight during peak spawning in June 2020. Figure S4. Stomach contents of Lagocephalus sceleratus showing freshly ingested cephalopods: Loligo vulgaris (A), Octopus vulgaris (B); and parts of Tremoctopus violaceus (C). Table S1. Sex ratios for L. sceleratus from Mediterranean and Suez Canal studies. Table S2. Length-weight relationships of Lagocephalus sceleratus from this study in southwestern Turkey. Table S3. Some length-weight parameters with sex, length range, length type and sample size (N) for L. sceleratus from Mediterranean and other studies. Table S4. Length-frequency data of Lagocephalus sceleratus collected in Southwestern Turkey from June 2019 to November 2020. Table S5. Growth parameters estimates for Lagocephalus sceleratus in the Mediterranean Sea (L in cm). Table S6. Mean length at first maturity (Lm) and maximum length (Lmax) for L. sceleratus at various locations of the Suez Canal and Mediterranean Sea. Table S7. Frequency of non-indigenous species (NIS) preyed upon by L. sceleratus in this study. Table S8. IRI Results of prey items for 34 sampled L. sceleratus.