The occurrence, distribution and biology of invasive fish species in fresh and brackish water bodies of NE Morocco

Abstract

The occurrence, distribution and biology of invasive fish species in fresh and brackish water bodies of NE Morocco

Monitoring the presence and expansion of alien species and upgrading their biological and ecological knowledge seems crucial to mitigate their possible impact on native communities. Within inland superficial waters, alien fish represent an important threat to the biodiversity and studies on their impact on native communities have increased around the world in the last years. However, little is known about their occurrence, biology and influences in North Africa in general, and more specifically in Morocco. In the present work we aimed to: 1) investigate the presence of any native Aphanius species, especially the Mediterranean killifish Aphanius fasciatus recorded from the lower basin of the Moulouya River (NE Morocco); 2) monitor the presence and expansion of two invasive species, the eastern mosquitofish Gambusia holbrooki and the mummichog Fundulus heteroclitus; and 3) contribute to the understanding of the ecological and abiotic affinities that govern the distribution of these alien fishes in North Africa. To achieve these goals, several field sampling campaigns were carried out between 2014 and 2018 across eastern Morocco, comprising the administrative Oriental Region and the Moulouya River Basin and covering an area of 119,268 km². No native Aphanius species were found. The eastern mosquitofish has invaded the freshwater hydrosystems of the northern part of Morocco, including the study area, while the mummichog is currently limited to the brackish and salty wetlands of Lower Moulouya. Our results show that the known ranges of these two alien species have expanded.

Data published through GBIF (Doi:10.15470/2qed9o)

Key words: Alien fishes, Mediterranean killifish, Eastern mosquitofish, Mummichog, Competition, Range expansion

Resumen

Incidencia, distribución y biología de las especies invasivas de peces en masas de agua dulce y salobre del NE de Marruecos

Monitorizar la presencia y expansión de especies exóticas y mejorar el conocimiento biológico y ecológico de las mismas parece clave para mitigar su posible impacto en las comunidades nativas. Las especies de peces exóticas constituyen una importante amenaza para la diversidad en las aguas superficiales continentales, por lo que los estudios sobre su impacto en las comunidades nativas se han incrementado en todo el mundo en los últimos años. No obstante, el conocimiento sobre la incidencia, biología e influencias de este fenómeno en el norte de África en general, y más específicamente en Marruecos, es escaso. En este trabajo nos propusimos: 1) investigar la presencia de alguna especie nativa de Aphanius, especialmente del fartet mediterráneo Aphanius fasciatus anteriormente registrado en la cuenca inferior del río Moulouya (NE de Marruecos); 2) monitorizar la presencia y expansión de dos especies invasoras, el pez mosquito Gambusia holbrooki y Fundulus heteroclitus, y 3) ayudar a comprender las afinidades ecológicas y abióticas que rigen la distribución de estos peces exóticos en el norte de África. Para alcanzar estos objetivos realizamos varias campañas de muestreo sobre el terreno en el este de Marruecos entre 2014 y 2018, incluidas la Región Oriental y la cuenca del río Moulouya, cubriendo una extensión de 119.268 km². No encontramos especies nativas de Aphanius. El pez mosquito oriental ha invadido los sistemas hídricos de agua dulce de la parte norte de Marruecos, incluida el área estudiada, mientras que Fundulus heteroclitus está limitado en la actualidad a las aguas salobres y saladas de los humedales del curso bajo del Moulouya. Como resultado de este estudio, las áreas de distribución conocidas de ambas especies exóticas se han ampliado.

Datos publicados en GBIF (Doi:10.15470/2qed9o)

Palabras clave: Peces exóticos, Fartet mediterráneo, Pez mosquito oriental, Fundulus heteroclitus, Competencia, Expansión potencial

Resum

Incidència, distribució i biologia de les espècies invasives de peixos en masses d’aigua dolça i salabrosa del NE del Marroc

Monitorar la presència i l’expansió d’espècies exòtiques i millorar-ne el coneixement biològic i ecològic sembla clau per mitigar l’impacte que poden exercir en les comunitats natives. Les espècies de peixos exòtiques constitueixen una amenaça important per a la diversitat a les aigües superficials continentals, per la qual cosa els estudis sobre el seu impacte en les comunitats natives s’han incrementat arreu del món els darrers anys. Tanmateix, el coneixement sobre la incidència, la biologia i les influències d’aquest fenomen al nord d’Àfrica en general, i més específicament al Marroc, és escàs. En aquest treball ens vam proposar: 1) investigar la presència d’alguna espècie nativa d’Aphanius, especialment del fartet mediterrani Aphanius fasciatus registrat anteriorment a la conca inferior del riu Moulouya (NE del Marroc); 2) monitorar la presència i l’expansió de dues espècies invasores, el peix mosquit Gambusia holbrooki i Fundulus heteroclitus, i 3) ajudar a entendre les afinitats ecològiques i abiòtiques que regeixen la distribució d’aquests peixos exòtics al nord d’Àfrica. Per assolir aquests objectius vam portar a terme diverses campanyes de mostreig sobre el terreny a l’est del Marroc entre 2014 i 2018, incloent-hi la Regió Oriental i la conca del riu Moulouya, cobrint una extensió de 119.268 km2. No hi vam trobar espècies natives d’Aphanius. El peix mosquit oriental ha envaït els sistemes hídrics d’aigua dolça de la part nord del Marroc, incloent-hi l’àrea estudiada, mentre que Fundulus heteroclitus està limitat actualment a les aigües salabroses i salades dels aiguamolls del curs baix del Moulouya. Com a resultat d’aquest estudi, les àrees de distribució conegudes de les dues espècies exòtiques s’han ampliat.

Dades publicades a GBIF (Doi:10.15470/2qed9o)

Paraules clau: Peixos exòtics, Fartet mediterrani, Peix mosquit oriental, Fundulus heteroclitus, Competència, Expansió potencial

Introduction

Biological invasions are one of the most important human impacts on a wide range of ecosystems and a major cause of global change (Ricciardi, 2006). The introduction and invasion of alien species are one of the main threats to biodiversity and ecosystem functioning worldwide, and represent the cause of huge ecological and economic costs around the world (Fausch and García-Berthou, 2013). Many invasive species have been implicated in species extinction, habitat degradation and ecosystem alteration (Cox, 1999; Pimentel et al., 2005; Coccia, 2015). Freshwater invasions have been less studied than terrestrial invasions, and the ecological and evolutionary consequences of most invasions remain unknown (Ricciardi and MacIsaac, 2011). Freshwaters are the most invaded and threatened ecosystem worldwide, with proportionally more invaders than terrestrial systems (Vitousek et al., 1997; Strayer, 2010).

Many members of the genus Aphanius are ‘victims’ of biological invasion, such as the Mediterranean killifish, A. fasciatus (Valenciennes, 1821), a cyprinodotid fish whose native range includes coastal brackish-waters mainly in the coastal zone of the central and eastern Mediterranean (Bianco, 1995). Nevertheless, it can also be found in various inland water bodies, such as lakes, inland streams and even in some North African oases (Hrbek and Meyer, 2003; Güçlü et al., 2013). Aphanius fasciatus is also present in various Mediterranean islands, such as Sardinia, Corsica and Cyprus (Bianco et al., 1996; Kottelat and Freyhof, 2007; Englezou et al., 2018). The Mediterranean killifish was reported in Morocco for the first time at the wetland of the Moulouya River mouth by Melhaoui (1994), but the species has not been reported since this record.

Aphanius fasciatus has disappeared from many sites where it once existed, and its populations have declined dramatically. In many cases it has even reached local extinction due to problems such as pollution of continental and coastal waters, brackish-water habitat degradation, destruction and reduction of salt-works and, most importantly, introduction of exotic fishes (Bianco, 1995). Competition with the alien species Gambusia holbrooki Girard, 1859 has strongly reduced the presence of A. fasciatus in many parts of its range (Kessabi et al., 2009; Valdesalici et al., 2015).

Commonly known as ‘Eastern mosquitofish’, Gambusia holbrooki is a freshwater poecilid native of America (Lloyd and Tomasov, 1985). It has been considered one of the world’ s one hundred worst invasive alien species by the GISP (Global Invasive Species Programme, http://www.issg.org/database/). It was introduced extensively for biological control against mosquito larvae and distributed in all continents except Antarctica, becoming a notorious pest worldwide (Cote et al., 2011; Pyke, 2005; Srean, 2015; Arnett, 2016) and causing serious problems for native fish such as A. fasciatus, Aphanius iberus (Valenciennes, 1846) or Valencia hispanica (Valenciennes, 1846) (Rincón et al., 2002; Caiola and de Sostoa, 2005; Valdesalici et al., 2015).

Another alien fish which has been reported to cause range reductions in native toothcarps is the mummichog Fundulus heteroclitus (Linnaeus, 1766), a small teleost fish, naturally occurring in saltmarshes of the Atlantic coast of North America from the Gulf of St Lawrence, southwards to north-eastern Florida (Shute, 1980). It was found in Iberia for the first time between 1970 and 1973 (Bernardi et al., 1995). Since then, F. heteroclitus has been found in many parts of Spain and Portugal (Morim et al., 2019). We found it for the first time in Morocco and Northern Africa at the Moulouya River mouth.

This study aimed to evaluate the presence of A. fasciatus or any native Aphanius species in NE Morocco, in parallel with the two invasive species G. holbrooki and F. heteroclitus. Furthermore, our goal was to study the interaction between the invasive fish, and their occurrence, distribution and biology in the fresh and brackish water bodies of the study area.

Material and methods

Study area

Morocco is currently divided into 12 administrative regions, including the Oriental Region (fig. 1), which occupies almost all the eastern side of the country and covers an area of 88,681 km² (see Mabrouki et al., 2018 for details). The Oriental Region includes the wilaya of Oujda (Oujda–Angad prefecture) and the provinces of Berkane, Driouch, Figuig, Guercif, Jerada, Nador and Taourirt. The watershed of the Moulouya (fig. 1), which includes nearly 43,412 km² of eastern Morocco, covers much of the Oriental Region. With a length of 600 km, the Moulouya is the longest Moroccan River flowing into the Mediterranean. Its main tributaries are the Oueds Ansegmir, Melloulou, Za and Msoun, all permanent. Other tributaries are presently intermittent (3–5 flash floods on average per year) (Bensaad et al., 2017; Mabrouki et al., 2017).

Fig. 1. Location of study area in Morocco and the prospected stations. Fig. 1. Localización del área de estudio en Marruecos y puntos de muestreo.

Surveys

Field investigations (often in the framework of various hydrobiological studies) have been carried out since 2014 at 45 stations throughout the Moulouya River basin, including its main affluents: Oued Anzegmir (side of the High Atlas), Oued Melloulou (Middle Atlas slope) and Oued Za (High Plateau), and at about 60 stations spread throughout eastern Morocco from the northern regions of Nador and Saïdia, to Figuig in the southeast and Talessint and Bouanane in the southwest (fig. 1). The stations that showed the occurrence of the mentioned fish species were re-examined between 2017 and 2018. Quantitative sampling of fish fauna was carried out using nets, searching in the most suitable places for the studied species. Sampling lasted an average of one hour over an area of 10 m² at each station, sufficient time to trap and catch virtually all the fish in each area. Invasive species caught were preserved in formalin solution, while native fish (especially fry) were returned to the water. Permission to perform these studied was granted by the authorities.

Abiotic factors

For this study, ten environmental parameters were selected: sulfate (SO₄^2–), biological oxygen demand after 5 days (BOD₅), phosphate (PO₄^3–) and nitrate (N-NO₃) were measured in the laboratory. Conductivity, salinity, pH, dissolved oxygen, temperature and the mean depth were measured in situ using a multiparametric measuring device (WTW, Multi-Line P4). Two replicas of water samples from each station were taken in 500 ml polyethylene bottles. The water samples were preserved with 2 ml of concentrated hydrochloric acid (pH = 2). According to standards ISO 5667-6 (1990), ISO 5667-2 (1991) and ISO 5667-3 (1994), water samples were transported in a cooler at a low temperature (± 4 ºC) to stop the metabolic activities of organisms in the water.

The current velocity, well known for its selective action on habitat and species distribution (Mabrouki et al., 2019a), was quantified by its mean value at three different locations of the same station. In the absence of a hydrometric reel to measure current velocity it was estimated using a stopwatch at various points of the watercourse by measuring the time it takes for a floating object to cross a given path. The average speed (converted to cm^–1) was semi-quantitatively estimated, followed by a transformation into three modalities: 1, very low to no current < 5 cm^–1; 2, low current 5 < 2 < 25 cm^–1; and 3, average current 25 < 3 < 50 cm^–1.

Statistical processing

Statistical analyses were carried out using software R package version 3.3.1. (R Core Team, 2019). In the modelling of counting processes, here the abundance of a species, two kinds of models are commonly implemented, namely the Poisson model and the negative binomial model (Hilbe, 2011). Frequently, counting data is characterized by overdispersion, whereby the sample variance is greater than its average. In the case where the anomalous dispersion is proved, the Poisson regression is no longer suitable for modelling this distribution, and the negative binomial regression model, allowing more flexibility in the dispersion, should be used. The relevance of Poisson regression or negative binomial models was evaluated by the Pearson Residue Test (Plackett, 1983).

Results

Occurrence, distribution and habitat types of the studied species

Despite our repeated sampling and the large number of prospected stations, we were unable to find A. fasciatus (or any Aphanius species). In contrast, its probable range was totally occupied by two exotic fish species F. heteroclitus and G. holbrooki. The species F. heteroclitus was found in six of 45 sampling sites while G. holbrooki appeared in 35.

In Morocco, the distribution of F. heteroclitus seems to be limited to the Oriental region so far, from the wetlands of Ain Chabbak to the mouth of the Moulouya River (fig. 2). Its habitat consists of coastal swamps and large sections of the Moulouya to its confluence with the Mediterranean, at an altitudinal ranging between 0 and 307 m a.s.l. In shallow habitats, where predators are less abundant, such as in the Ain Chebbak wetland, F. heteroclitus can be found swimming freely in pelagic and superficial parts along with G. holbrooki. On the other hand, in the main watercourse of the lower Moulouya, where the depth is greater and predatory fishes are more abundant, F. heteroclitus was only caught among the aquatic vegetation, always together with G. holbrooki, except near the river mouth in extreme lower Moulouya (stations S37, S78 and S39), where it was captured alone.

Fig. 2. Distribution of the studied species in Morocco (black dots represent new records, white circles represent bibliographic). Fig. 2. Distribución de las especies estudiadas en Marruecos (nuevos registros en negro, registros bibliográficos en blanco).

With regard to G. holbrooki, the species appears to be widely distributed through Morocco, including the Moulouya River Basin (fig. 2). It is highly confined to brackish waters near the coast (as in Nador lagoon and the Moulouya SIBE), and it is omnipresent along Za watercourses, Sebra and the complex Zeghzel-Cherraa rivers. It can go up to Melloulou and the Middle Moulouya during the summer period with the weakening current. The species also appears to be well established in other areas of the oriental region of Morocco, from Nador and Al Hoceima north and to Figuig south-east and the Talessint region south-west. The eastern mosquitofish seems able to cope remarkably well in a wide range of habitats, even those that are highly degraded as long as there is water. Its habitats range from large permanent streams of average and low altitudes, to dams, natural and managed springs, and small temporary tributaries (provided that they are connected to a permanent source or underground water). It can be also found in human constructions such as irrigation canals and cement tanks for water retention.

The physicochemical variables of water

Both G. holbrooki and F. heteroclitus species are well known for their tolerance to large variations in abiotic parameters (Kneib 1986, Lipcius and Subrahmanyam 1986; Pyke 2008). However, in the study area, G. holbrooki appears to excel F. heteroclitus in term of oscillation intervals of most of the studied parameters (see appendix). G. holbrooki showed great plasticity in terms of salinity; was recorded in low mineralized waters with a minimum of 0.587 g/l  (587.67 ppm), and also in the brackish waters of the coastal system with a maximum of 1,837g/l (1,839.09 ppm). The greater plasticity of the species was demonstrated mainly by the large intervals of the organic pollution indicators, showing that it can tolerate BOD5 concentrations between 1.49 and 20.3 mg/l, nitrate between 1.25 and 67 mg/l , sulfate between 33 and 403 mg/l phosphate between 0.016 and 2.82 mg/l and finally pH between 6.65 and 8.3. This same observation was noted for dissolved oxygen; the species seems to withstand great fluctuations of this vital element and it can be found not only in very well oxygenated waters 11.3 mg/l but also in weakly oxygenated waters with a minimum of 2.01 mg/l.

Appendix. Abundance and value of abiotic parameters measured at the sampling sites where G. holbrroki and F. heteroclitus were found: Gam, abundance of Gambusia; Fun, abundance of Fundulus; Depth, in cm; Temp, water temperature in ºC; Sali, salinity in g/l; Dis_O, dissolved oxygen in mg/l; NO3, nitrate in mg/l; SO, sulfate in mg/l; PO, phosphate in mg/l; BOD5 in mg/l; Vel, current velocity in cm/s. Apéndice. Abundancia y valor de los parámetros abióticos medidos en los puntos de muestreo donde se encontraron G. holbrroki y F. heteroclitus. (Para las abreviaturas, véase arriba).

Fundulus heteroclitus is instead more demanding in terms of the physicochemical quality of the water and showed weak ranges of variation for most of the studied parameters (see appendix). However, it can survive in a hypersaline waters compared to G. holbrooki, with a maximum of 2,757 g/l (2,760.15 ppm) recorded at the Moulouya mouth (S39).

Concerning temperature, it is difficult to accurately evaluate variation for these two species in the study area because we did not perform a monthly follow-up in all the stations because of the violent floods recorded during the cold period. However, water temperature is linked to local conditions (climate, duration of sunshine, flow and altitude) and to seasonality in the Moulouya River Basin (Mabrouki et al., 2016; Taybi, 2016). During the study period, the maximum temperature did not exceed 26 ºC for either species. Given the lack of sampling during winter, our study does not report the lowest temperature at which the fish remain alive or active.

Modelling results

Poisson regression and negative binomial regression showed the same explanatory variables of the abundance frequency, with similar effects. However, the Poisson regression model showed an abnormal over-dispersion, and its statistical relevance was rejected by the Pearson Residue Test. The negative binomial regression model was better adapted to these data, and its statistical relevance was verified.

Discussion

Growth in international trade and concurrent increases in transport capacity have accelerated the rate of introduction of alien species worldwide, with freshwater ecosystems, and their native fish communities, being particularly susceptible (Sala et al., 2000; Kolar and Lodge, 2002; Macdonald and Tonkin, 2008). The most important pathways of freshwater invasions are aquaculture, pet/aquarium trade and stocking activities, shipping, and inland canals, without forgetting that some species can be introduced intentionally (Coccia, 2015; Nunes et al., 2015). Invasive species are a major threat to global biodiversity, causing the decline and extinction of native freshwater fish species throughout the world (Strauss et al., 2006; García-Berthou et al., 2005; Srean, 2015). The piscivorous ‘Nile perch’ Lates niloticus (Linnaeus, 1758) is a classical example, as its introduction in the 1950s in Lake Victoria in East Africa led to the extinction of over 200 endemic fishes (Kitchell et al., 1997).

During our survey, we did not find A. fasciatus in Morocco in spite of our exhaustive sampling. This could be due to several reasons, such as pollution and destruction of the habitat on one hand, and the impact of alien fish on the other, as has occurred in many areas of its range (IUCN, 2006). We should take into account that the record of A. fasciatus from Morocco could be erroneous due to misidentification and the fact that scientific collections from Morocco contain no samples. Although healthy populations still exist in some parts of the Maghreb, occupying the same type of habitats in Algeria and Tunisia (Kraiem, 1983; Boumendjel et al., 2015), the palaeogeographic history of Morocco differs. Until 5.3 million years ago (Krijgsman et al., 1999a, 1999b), together with Spain, the north of Morocco formed the Betico-Rifian massif. As Spain and Morocco shared fauna not present in the other countries, in Morocco populations it would be more likely to find species related to Spanish aphaniids, such as A. baeticus Doadrio, Carmona and Fernández-Delgado, 2002 or A. iberus than A. fasciatus. On the other hand, the present distribution of A. iberus, A. baeticus, and even A. saourensis Blanco, Hrbek and Doadrio, 2006 occurs along the coastline that existed in the Miocene before the separation of the Betico-Rifian Massif at the Miocene-Pliocene boundary 5.3 million years ago (Doadrio, 1994; Krijgsman et al., 1999a, 1999b). It is therefore more likely to find Aphanius in inland Morocco than in coastal rivers.

In consensus with these observations, the previous record of A. fasciatus from Morocco would be the misidentification of a different species of the genus. Nevertheless, the idea that alien fishes have driven the Moroccan population of a native Aphanius species to extinction is also possible, especially in eastern Morocco, where G. holbrooki and F. heteroclitus have invaded all suitable hydrosystems, from the fresh continental waters of rivers and springs, to brackish and salty waters of coastal lagoons and wetlands.

The nature of G. holbrooki with its high reproductive potential, fast maturation rate, flexible behavior and broad environmental tolerances have contributed to its success as an invader, and the species is considered to pose a serious threat not only to native fishes worldwide but also to amphibians, macro-invertebrates, zooplankton and phytoplankton communities (Courtney and Meffe, 1989; Howe et al., 1997; Caiola and de Sostoa, 2005; Macdonald and Tonkin, 2008; Rowe et al., 2008; Mabrouki et al., 2019b). Usually, the species of the genus Aphanius can take refuge in estuarine habitats to escape competition with G. holbrookii due to its strong salt-tolerant character like all cyprinodontids, as has previously been reported (Nordlie and Mirandi, 1996; Doadrio, 2002; Alcaraz et al., 2008; Valdesalici et al., 2015). However, in our case these refuges are occupied by another alien fish, which is also a cyprinodontiform specialized in mesohalyne and hyperhalyne waters and which is possibly responsible for the population regression of A. baeticus in the Iberian Peninsula (Oliva-Paterna et al., 2006).

Regarding the ecological affinities of the two invasive species, the results of the negative binomial regression of mosquitofish abundance as a function of the environmental variables selected for this study (table 1) showed the existence of a highly significant static relationship between Gambusia abundance and salinity, and water velocity (p-value < 0.001), and a very significant relation with depth (p-value between [0.001, 0.01]). This correlation is negative for the three variables, i.e., if the salinity of the environment decreases by 1.4, the abundance of G. holbrooki will increase by one. G. holbrooki has extensive tolerance for several environmental factors as reported in other works (Pyke, 2008), but according to our results it seems to prefer shallow waters with average to low salinity and slow current velocity, in agreement with findings in previous studies (Stearns and Sage, 1980; Brown-Peterson and Peterson, 1990). Its peak of colonization will thus be carried out, most probably, in shallow lentic water bodies with low or moderate salinity. Brown-Peterson and Peterson (1990) reported that females from freshwater were in better somatic condition and matured later, since ovarian development was significantly more advanced in oligohalyne waters. The same model was used for F. heteroclitus; a positive correlation was detected between this species and salinity (table 2), and the relevance of this regression model was verified by the Pearson residue test; the model is well adapted to the data.

Table 1. Estimates of the parameters of the negative binomial regression model (for Gambusia holbrooki): SE, standard error; IC, interval of confidence. Tabla 1. Estimaciones de los parámetros del modelo de regresión binomial negativa (para Gambusia holbrooki): SE, error estándar; IC, intérvalo de confianza.

Table 2. Estimates of the parameters of the negative binomial regression model (for Fundulus heteroclitus): SE, standard error; IC, interval of confidence. Table 2. Estimaciones de los parámetros del modelo de regresión binomial negativa (para Fundulus heteroclitus): SE, error estándar; IC, intérvalo de confianza.

To conclude, both invasive species generally occur in water that is shallow, stagnant or slow moving, with dense aquatic and riparian vegetation. In water bodies that are moderately deep, they are generally found along the shallow edges, especially those that are well vegetated. However, F. heteroclitus seems to be specialized in brackish saline waters of coastal hydrosystems, while G. holbrooki seems specialized in freshwater and slightly brackish water, showing scarce selectivity in terms of habitat, as was observed in its native area (Moyle and Nichols, 1973). Indeed, this species can occur in environments that range from almost completely undisturbed swamps, lakes, and streams to highly disturbed water bodies including urban drains with various water quality problems such as pollution and dense growth of aquatic and surface weeds. G. holbrooki is often extremely abundant in disturbed habitats near urban areas (Arthington and Lloyd, 1989; Pyke, 2008). This could be explained by its high tolerance to fluctuations in abiotic factors (Odum and Caldwell, 1955; Keup and Bayliss, 1964; Cherry et al., 1976; Fastelli et al., 2012), this trait being one of the main causes behind its expansion success and outcompeting F. heteroclitus in the study area. Nevertheless, the apparently limited tolerance towards different environments observed in F. heteroclitus might also be due to its recent introduction in the study area, or to the presence of G. holbrooki, that could outcompete it in more internal water bodies.

Research on invasive species is an urgent research topic for the conservation of biodiversity (Dukes and Mooney, 2004; Rice and Sax, 2005; Taybi et al., 2020a). Studies on the impact of invasive fish on native fauna have flourished around the world in recent years. However, little is known about their occurrence, biology and influence in North Africa in general or in Morocco in particular. Since G. holbrooki was first introduced in Morocco in 1929 to fight against waterborne diseases, it has invaded all the northern Moroccan hydrosystems (Azeroual, 2003; Clavero et al., 2014). Our surveys in eastern Morocco have significantly increased its known range. Future assessments in other parts of the country could expand the species’ distribution area, especially in the Western Rif and the High and the Anti-Atlas valley, and even in the Sahara. On the other hand, F. heteroclitus seems to be limited to the wetland of Moulouya mouth for the time being, but thanks to its ‘salt-tolerant’ character, this newcomer could soon invade habitats of the same type. Indeed, considering the habitat destruction and deterioration of the quality of surface waters currently occurring in Morocco (Mabrouki et al., 2016; Taybi et al., 2016; Bensaad et al., 2017; Mabrouki et al., 2017; Taybi et al., 2020b), in addition to the climate change effects and intensification of saline stress (Moss et al., 2009), both alien fishes would be in a better position to invade North African ecosystems. It is also of note that freshwater ichtyofauna of Morocco is still little known and new species are being recorded (Ford et al., 2020). Several native fish species could co-occur with the two invaders in eastern Morocco and in the Moulouya River basin, such as the North African shad Alosa algeriensis Regan, 1916 and many endemic barbels (i.e. Luciobarbus guercifensis Doadrio, Casal-López and Yahyaoui, 2016; L. lanigarensis Brahimi, Libois, Henrard and Freyhof, 2018 and L. yahyaouii Doadrio, Casal-Lopez and Perea, 2016). Knowledge of the distribution of both native and non-native freshwater fishes in Morocco therefore merits urgent research.

References

Alcaraz, C, Bisazza, A., Garcıa-Berthou, E., 2008. Salinity mediates the competitive interactions between invasive mosquitofish and an endangered fish. Oecologia, 155: 205-213, Doi 10.1007/s00442-007-0899-4

Arnett, H., 2016. Sources of Ecologically Important Trait Variation in Mosquitofish (Gambusia affinis and Gambusia holbrooki). Electronic Theses and Dissertations, 2692, https://digitalcommons.library.umaine.edu/etd/2692

Arthington, A. H., Lloyd, L. N., 1989. Introduced Poeciliids in Australia and New Zealand. In: Ecology and evolution of livebearing fishes (Poeciliidae): 333-348 (G. K. Meffe, F. F. Snelson Jr., Eds.). Prentice Hall, Inc, New Jersey.

Azeroual, A., 2003. Monograph of the fishes of the continental waters of Morocco: systematics, distribution and ecology. PhD thesis, Mohammed V University, Rabat. [In French].

Bensaad, H., Mabrouki, Y., Taybi, A. F., Chafi, A., 2017. Assessment of wastewater discharges from Taourirt City on the water quality of the Oued Za (Eastern Morocco). Journal of Materials and Environmental Science, 8: 2365-2371.

Bernardi, G., Fernández‐Delgado, C., Gómez-Chiarri, M., Powers, D. A., 1995. Origin of a Spanish population of Fundulus heteroclitus inferred by cytochrome b sequence analysis. Journal of Fish Biology, 47: 737-740.

Bianco, P. G., 1995. Mediterranean endemic freshwater fishes of Italy. Biological Conservation, 72: 159-170.

Bianco, P. P., Ahnet, H., Economides, P. S., 1996. The freshwater fishes from eastern and large Mediterranean islands with comments on their safety status. Acta Universitatis Carolinae. Biologica, 40: 45-60.

Blanco, J. L., Hrbek, T., Doadrio, I., 2006. A new species of the genus Aphanius Nardo, 1832 (Actinopterygii, Cyprinodontidae from Algeria. Zootaxa, 1158(1): 39-53.

Boumendjel, M., Taibi, F., Hennouni, N., Faouzi Samar, M., Abdesselem, A., Slimane, B., Laabidi, Z., Laissani, A., Merdaci, D., Alleb, W., Djouadi, H., Menam, A., Allou, B., 2015. Ecologie et polymorphisme de l’Aphanius fasciatus Valenciennes 1821 Cyprinodontidae, dans le lac bleu, parc national d’El-Kala (nord-est Algérie). Sciences & Technologie, 41: 9-21.

Brown-Peterson, N., Peterson, M. S., 1990. Comparative life history of female mosquitofish, Gambusia affinis, in tidal freshwater and oligohaline habitats. Environmental Biology of Fishes, 27: 33-41.

Caiola, N., de Sostoa, A., 2005. Possible reasons for the decline of two native toothcarps in the Iberian Peninsula: Evidence of competition with the introduced Eastern mosquitofish. Journal of Applied Ichthyology, 21: 358-363, Doi:10.1111/j.1439-0426.2005.00684.x

Cherry, D. S., Guthrie, R. K., Rodgers, J. H., Cairns, J., Dickson, K. L., 1976. Responses of mosquitofish G. affinis to ash effluent and thermal stress. Transactions of the American Fisheries Society, 105: 686-94.

Clavero, M., Esquivias, J., Qninba, A., Riesco, M., Calzada, J., Ribeiro, F., Fernández, N., Delibes, M., 2014. Fish invading deserts: non-native species in arid Moroccan rivers. Aquatic Conservation: Marine and Freshwater Ecosystems, 25(1): 49-60, Doi:10.1002/aqc.2487

Coccia, C., 2015. Invasion biology of Trichocorixa verticalis in Doñana, SW Spain. PhD thesis, University of Sevilla.

Cote, J., Fogarty, S., Brodin, T., Weinersmith, K., Sih, A., 2010. Personality-dependent dispersal in the invasive mosquitofish: Group composition matters. Proceedings of the Royal Society B: Biological Sciences, 278: 1670-1678, Doi: 10.1098/rspb.2010.1892

Courtney, W. R. J., Meffe, G. K., 1989. Small fishes in strange places: a review of introduced poeciliids. In: Ecology and evolution of livebearing fishes (Poeciliidae): 319-332 (G. K. Meffe, F. F. Snelson Jr., Eds.). Prentice-Hall Inc., New Jersey.

Cox, G. W., 1999. Alien species in North America and Hawaii: impacts on natural ecosystems. Island Press, Washington DC.

Doadrio, I., 1994. Freshwater fish fauna of North Africa and its biogeography. Annales du Musée Royal de l’Afrique Centrale, 275: 21-34.

– 2002. Atlas y libro rojo de los peces continentales de España. Ministerio de Medio Ambiente, Madrid.

Dukes, J. S., Mooney, H. A., 2004. Disruption of ecosystem processes in western North America by invasive species. Revista Chilena de Historia Natural, 77: 411-437.

Englezou, C., Gücel, S., Zogaris, S., 2018. A new record of Aphanius fasciatus Valenciennes, 1821 on Cyprus; insights for conservation. Cahiers de Biologie Marine, 59: 379-384, Doi: 10.21411/CBM.A.171FB1DB

Fausch, K. D., García-Berthou, E., 2013. The problem of invasive species in river ecosystems. In: River Conservation: Challenges and Opportunities: 193-216 (S. Sabater, A. Elosegi, Eds.). Fundación BBVA, Bilbao.

Fastelli, P., Marcelli, M., Monti, F. 2012. Presenza e distribuzione di Aphanius fasciatus (Valenciennes, 1821) (Osteichthyes, Cyprinodontidae) e Gambusia holbrooki (Girard, 1859) (Osteichthyes, Poecilidae) nelle zone umide del parco regionale della Maremma (Toscana). Atti del Museo di Storia Naturale della Maremma, 23: 3-8.

Ford, M., Brahimi, A., Baikeche, L., Bergner, L., Clavero, M., Doadrio, I., Lopes-Lima, M., Perea, S., Yahyaoui, A., Freyhof, J., 2020. Freshwater fish distribution in the Maghreb: a call to contribute. OSF Preprints. Doi:10.31219/osf.io/kx4gc

García-Berthou, E., Alcaraz, C., Pou-Rovira, Q., Zamora, L., Coenders, G., Feo, C., 2005. Introduction pathways and establishment rates of invasive aquatic species in Europe. Canadian Journal of Fisheries and Aquatic Sciences, 62: 453-463.

Güçlü, S., Küçük, F., Osman, Ö., Güçlü, E., 2013. The fish fauna of the Büyük Menderes River (Turkey): taxonomic and zoogeographic Features. Turkish Journal of Fisheries and Aquatic Sciences, 13: 685-698.

Hilbe, J. M., 2011. Negative Binomial Regression. Cambridge University Press, Cambridge.

Howe, E., Howe, C. R. L., Burchett, M., 1997. Impact of the introduced poeciliid Gambusia holbrooki (Girard, 1859) on the growth and reproduction of Pseudomugil signifier (Kner, 1865) in Australia. Marine and Freshwater Research, 48: 425-434.

Hrbek, T., Meyer, A., 2003. Closing of the Tethys Sea and the phylogeny of Eurasian killifishes (Cyprinodontiformes: Cyprinodontidae). Journal of Evolutionary Biology, 16: 17-36.

ISO (International Organization for Standardization), 1990. Water Quality Sampling. Part 6: Guidelines for Sampling of Rivers and Streams. International Standard ISO 5667-6, Geneva, Switzerland.

– 1991. Water Quality Sampling. Part 2: Guidelines for Sampling of Rivers and Streams. International Standard ISO 5667-2, Geneva, Switzerland.

– 1994. Water Quality Sampling. Part 3: Guidance on the Preservation and Handling of Samples. International Standard ISO 5667-3, Geneva, Switzerland.

IUCN, 2006. Aphanius fasciatus (Mediterranean Killifish, South European Toothcarp). In: The IUCN Red List of Threatened Species 2006: e.T1847A8316811, Doi: 10.2305/IUCN.UK.2006.RLTS.T1847A8316811.en

Kessabi, K., Kerkeni, A., Saïd, K., Messaoudi, I., 2009. Involvement of Cd bioaccumulation in spinal deformities occurrence in natural populations of Mediterranean Killifish. Biological Trace Element Research, 128: 72-81,  Doi: 10.1007/s12011-008-8255-z

Keup, L., Bayliss, J., 1964. Fish distribution at varying salinities in Neuse River basin, North Carolina. Chesapeake Science, 5: 119-123.

Kitchell, J. F., Schindler, D. E., Oguto-Ohwayo, R., Reinthal, P., 1997. The Nile perch in lake Victoria: interactions betweeen predation and fisheries. Ecological Applications, 7: 653-664.

Kneib, R. T., 1986. The role of Fundulus heteroclitus in salt marsh trophic dynamics. American Zoologist, 26: 259–269.

Kolar, C. S., Lodge, D. M., 2002. Ecological predictions and risk assessment for alien fishes in North America. Science, 298: 1233-1236.

Kottelat, M., Freyhof, J., 2007. Handbook of European freshwater fishes. Publications Kottelat, Cornol, Switzerland.

Kraiem, M. M., 1983. Les poissons d’eau douce de Tunisie: inventaire comment et repartition géographique. Bulletin de l’institut National Scientifique et Technique d’océanographie et de Pêche de Salammbô, 10: 107-124.

Krijgsman, W., Hilgen, F. J., Marabini, S., Vai, G. B., 1999a. New paleomagnetic and cyclostratigraphic age constraints on the Messinian of the Northern Apennines (Vena del Gesso Basin, Italy). Memorie della Società Geologica Italiana, 54: 25-33.

Krijgsman, W., Hilgen, F. J., Raffi, I., Sierro, F. J., Wilson, D. S., 1999b. Chronology, causes and progression of the Messinian salinity crisis. Nature, 400: 652-655.

Lipcius, R. N., Subrahmanyam, C. B., 1986. Temporal factors influencing killifish abundance and recruitment in Gulf of Mexico salt marshes. Estuarine, Coastal, and Shelf Science, 22: 101-114.

Lloyd, L. N., Tomasov, J. F., 1985. Taxonomic status of mosquitofish, Gambusia affinis (Poeciliidae), in Australia. Australian Journal of Marine and Freshwater Research, 36: 447-451.

Mabrouki, Y., Taybi A. F., Bensaad, H., Berrahou, A., 2016. Variabilité spatio–temporelle de la qualité des eaux courantes de l’Oued Za (Maroc Oriental). Journal of Materials and Environmental Science, 7(1): 231-243.

Mabrouki, Y., Taybi, A. F., Berrahou, A., 2017. L’évolution spatio–temporelle de la qualité des eaux courantes de l’Oued Melloulou (Maroc). Revue des Sciences de l’Eau, 30(3): 213-225.

Mabrouki, Y., Taybi, A. F., Berrahou, A., Chaabane, K., Legssyer, B., 2018. Case study of the freshwater shrimp Atyaephyra desmarestii (Millet, 1831). Crustacea, Decapoda. in the watershed of Moulouya and Eastern Morocco. Life and Environment, 68(4): 175-183.

Mabrouki, Y., Taybi, A. F., El Alami, M., Berrahou, A., 2019a. Biotypology of stream macroinvertebrates from North African and semi arid catchment: Oued Za (Morocco). Knowledge and Management of Aquatic Ecosystems, 420: 17.

Mabrouki, Y., Taybi, A. F., Skalli, A., Sánchez-Vialas, A., 2019b. Amphibians of the Oriental Region and the Moulouya River Basin of Morocco: distribution and conservation notes. Basic and Applied Herpetology, 33:19-32.

Macdonald, J., Tonkin, Z., 2008. A review of the impact of eastern gambusia on native fishes of the Murray-Darling Basin. Arthur Rylah Institute for Environmental Research, Department of Sustainability and Environment, Heidelberg, Victoria.

Melhaoui, M., 1994. Ecologie des ressources halieutiques des eaux continentales à interêt économique: Anguillidae, Salomonidae. PhD thesis, University Mohamed First Oujda [in French].

Morim, T., Bigg, G. R., Madeira, P. M., Palma, J., Duvernell, D. D., Gisbert, E., Cunha, R. L., Castilho, R., 2019. Invasion genetics of the mummichog (Fundulus heteroclitus): recent anthropogenic introduction in Iberia. Peer J., 7: e6155, Doi: 10.7717/peerj.6155

Moss, B. D., Hering, A. J., Green, A., Aidoud, A., Becares, E., Beklioglu, M., Bennion, Boix, D., Brucet, S., Carvalho, L., Clement, B., Davidson, T., Declerck, S. H., Dodson, M., Van Donk, E., Dudley, B., Feuchtmayr, H., Friberg, N., Grenouillet, G., Hillebrand, H., Hobaek, A., Irvine, K., Jeppesen, E., Johnson, R., Jones, I., Kernan, M., Luridsen, T. L., Manca, M., Meerhoff, M., Olafsson, J., Ormerod, S., Papastergiadou, E., Ellis Penning, W., Ptacnik, R., Quintana, X., Sadin, L., Seferlis, M., Simpson, G., Triga, C., Verdoschot, P., Antonie, M., Weyhenmeyer, G. A., 2009. Climate change and the future of freshwater 340 biodiversity in Europe: A primer for policy-makers. Freshwater Reviews, 2: 103-130, Doi: 0.1608/FRJ-2.2.1

Moyle, P. B., Nichols, R. D., 1973. Ecology of some native and introduced fishes of the Sierra Nevada foothills in Central California. Copeia, 1973: 478-490.

Nordlie, F. G., Mirandi, A., 1996. Salinity relationships in a freshwater population of eastern mosquitofish. Journal of Fish Biology, 49(6): 1226-1232, Doi: 10.1111/j.1095-8649.1996.tb01791.x

Nunes, A. L., Tricarico, E., Panov, V. E., Cardoso, A. C., Katsanevakis, S., 2015. Pathways and gateways of freshwater invasions in Europe. Aquatic Invasions, 10(4): 359-370.

Odum, H. T., Caldwell, D. K., 1955. Fish respiration in the natural oxygen gradient of an anaerobic spring in Florida. Copeia, 104-106.

Oliva-Paterna, F. J., Doadrio, I., Fernández-Delgado, C., 2006. Threatened fishes of the world: Aphanius baeticus Doadrio, Carmona and Fernández Delgado, 2002 (Cyprinodontidae). Environmental Biology of Fishes, 75(4): 415-417.

Pimentel, D., Zuniga, R., Morrison, D., 2005. Update on the environmental and economic costs associated with alien-invasive species in the United States. Ecological Economics, 52: 273-288.

Plackett, R. L., 1983. Karl Pearson and the Chi-Squared Test. International Statistical Review, 51(1): 59-72.

Pyke, G. H., 2005. A review of the biology of Gambusia affinis and G. holbrooki. Reviews in Fish Biology and Fisheries, 15: 339-365, Doi: 10.1007/s11160-006-6394-x

– 2008. Plague Minnow or Mosquito Fish? A Review of the Biology and Impacts of Introduced Gambusia Species. Annual Review of Ecology, Evolution, and Systematics, 39: 171-191.

R Core Team, 2019. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. Url:https://www.R-project.org/

Ricciardi, A., 2006. Are modern biological invasions an unprecendent form of global change? Conservation Biology, 21: 329-336.

Ricciardi, A., MacIsaac, H. J., 2011. Impacts of biological invasions on freshwater ecosystems. In: Fifty Years of Invasion Ecology: The Legacy of Charles Elton: 211-224 (D. M. Richardson, Ed.). Wiley-Blackwell, Oxford, UK.

Rice, W. R., Sax, D., 2005. Testing fundamental evolutionary questions at large spatial and demographic scales: species invasions as an underappreciated tool. In: Species invasions: insights into ecology, evolution and biogeography: 291-308 (D. F. Sax, J. J. Stachowicz, S. D. Gaines, Eds.). Sinauer Associates, Sunderland, UK.

Rincón, P. A., Correas, A. M., Morcillo, F., Risueño, P., Lobón-Cerviá, J., 2002. Interaction between the introduced eastern mosquitofish and two autochthonous Spanish toothcarps. Journal of Fish Biology, 61: 1560-1585, Doi:10.1006/jfbi. 2002.2175

Rowe, D. K., Moore, A., Giorgetti, A., Maclean, C., Grace, P., Wadhwa, S., Cooke, J., 2008. Review of the impacts of gambusia, redfin perch, tench, roach, yellowfin goby and streaked goby in Australia. Draft report prepared for the Australian Government Department of the Environment, Water, Heritage and the Arts, Canberra.

Sala, O. E., Chapin, F. S., Armesto, J. J., Berlow, E., Bloomfield, J., Dirzo, R., Huber-Sanwald, E., Huenneke, L. F., Jackson, R. B., Kinzig, A., Leemans, R., Lodge, D. M., Mooney, H. A., Oesterheld, M., Poff, N. L., Sykes, M. T., Walker, B. H., Walker, M., Wall, D. H., 2000. Global biodiversity scenarios for the year 2100. Science, 287: 1770-1774.

Shute, P. W., 1980. Menidia extensa. In: An Atlas of North American Freshwater Fishes: 559 (D. S. Lee, C. R. Gilbert, C. H. Hocutt, R. E. Jenkins, D. E. McAllister, J. R. Stauffer, Jr. Eds.). North Carolina Museum of Natural History, Raleigh, NC.

Srean, P., 2015. Understanding the Ecological Success of Two Worldwide Fish Invaders (Gambusia holbrooki and Gambusia affinis). PhD thesis, University of Girona.

Stearns, S., Sage, R. D., 1980. Maladaptation in a marginal population of the mosquito fish. Evolution, 34: 65-75.

Strauss, S. Y., Lau, J. A., Carroll, S. P., 2006. Evolutionary responses of natives to introduced species: what do introductions tell us about natural communities? Ecology Letters, 9: 357-374.

Strayer, D. L., 2010. Alien species in fresh waters: ecological effects, interactions with other stressors, and prospects for the future. Freshwater Biology, 55: 152-174.

Taybi, A. F., 2016. Hydrobiological study of the Moulouya: Structure of the biodiversity and longitudinal zonation of benthic invertebrates. PhD thesis, University Mohamed First, Oujda, Morocco. [in French]

Taybi, A. F., Mabrouki, Y., Berrahou, A., Chaabane, K., 2016. Évolution spatiotemporelle des paramètres physico-chimiques de la Moulouya. Journal of Materials and Environmental Science, 7(1): 272-284.

Taybi, A. F, Mabrouki, Y., Chavanon, G., Millán, A., 2020a. The alien boatman Trichocorixa verticalis verticalis. Hemiptera: Corixidae. strongly increases its presence in Morocco. Limnetica, 39(1): 49-59, Doi: 10.23818/limn.39.04

Taybi, A. F., Mabrouki, Y., Legssyer, B., Berrahou, A., 2020b. Spatio-temporal typology of the physico-chemical parameters of a large North African River: the Moulouya and its main tributaries (Morocco). African Journal of Aquatic Sciences (in press).

Valdesalici, S., Langeneck, J., Barbieri, M., Castelli, A., Maltagliati, F., 2015. Distribution of natural populations of the killifish Aphanius fasciatus (Valenciennes, 1821). Teleostei: Cyprinodontidae. in Italy: past and current status, and future trends. Italian Journal of Zoology, 82: 212-223, Doi: 10.1080/ 11250003.2014.1003418

Vitousek, P. M., D’ Antonio, C. M., Loope, L., Rejmanek, M., Westbrooks, R. G., 1997. Introduced species: a significant component of human-caused global environmental change. New Zealand Journal of Ecology, 21:1-16.