Research Article
Print
Research Article
A population dynamics approach to understand the invasiveness of the seaweed Rugulopteryx okamurae (Ochrophyta, Dictyotales)
expand article infoJesús Rosas-Guerrero, Raquel Carmona, Julio De la Rosa§, Marianela Zanolla, María Altamirano
‡ Universidad de Málaga, Málaga, Spain
§ Universidad de Granada, Granada, Spain
Open Access

Abstract

The success of invasive species can be measured by invasiveness, which depicts intrinsic characteristics that enable them to thrive in new environments. In invasive seaweeds, for example, the persistence of multiple overlapping cohorts throughout the year plays a key role in increasing plant cover and exerting unrelenting pressure on invaded areas. The marine brown macroalgae Rugulopteryx okamurae has recently established abundant populations in the Mediterranean Sea and Atlantic Ocean, negatively affecting both biodiversity and socioeconomic factors by unprecedently aggressive invasive behaviour. The objective of the study is to understand the invasiveness of R. okamurae through its popu lation dynamics. For this, a year-round study was conducted in a protected habitat of Posidonia oceanica in southern Spain, revealing that R. okamurae uses alternating mechanisms for population maintenance. It achieves high density of young individuals in late summer and autumn, peaking at 3285 individuals per square metre. In spring and early summer, the population shifts towards fewer – but larger – individuals, with densities dropping to 888 individuals per square metre and biomass reaching a peak of 170 g dry weight (DW) per square metre. Six overlapping cohorts were identified by Gaussian curves. They persisted throughout the year, but they were not related to environmental factors, which indicates adaptive physiological mechanisms that sustain dense monospecific populations. Additionally, the association between cohorts and different morphotypes suggests that R. okamurae phenotypic plasticity enables its persistence in introduced areas. These findings provide valuable insights into the biological traits underpinning its invasiveness in P. oceanica meadows, revealing temporal windows of invasiveness driven by different mechanisms. This knowledge is crucial for developing effective conservation and management strategies aimed at mitigating the impact of this invasive species.

Key words

Demography, density, generation, macroalga, Posidonia oceanica, recruitment

Introduction

Invasion success refers to an invasive species’ ability to progress through all stages of invasion, from introduction to dispersion, by overcoming ecological barriers and progressing through successive stages of establishment and expansion in a new environment (Gioria et al. 2023). Studies aiming to identify and understand the factors that contribute to the success of invasive species have led to the development of various invasive hypotheses (Lowry et al. 2013). Currently, there are more than thirty invasion hypotheses that try to explain the relationship between invasive species, invaded communities and their interactions (Enders et al. 2018). According to Gioria et al. (2023), invasion hypotheses can be grouped into three major categories, namely: (1) propagule pressure, which involves the introduction efforts by invasive species in a new area and is defined as the number and frequency of propagule release (Lockwood et al. 2005; Colautti et al. 2006), (2) invasibility, referring to the characteristics of the receiving ecosystem and its susceptibility to being invaded (Catford et al. 2012; Colautti et al. 2014) and (3) invasiveness, referring to intrinsic characteristics of the species, including functional traits that determine its invasive potential (Pyšek et al. 2014; Jehangir et al. 2024). In this context, identifying the traits that serve as good predictors of invasiveness involves several factors, including physical performance, growth rate, size, biomass allocation, physiology and phenotypic plasticity (Richardson and Pyšek 2006). The latter expands the range of responses to biological and abiotic factors encountered in introduced areas, allowing the species to adapt quickly to the new environment and facilitating their establishment and expansion (Richards et al. 2006). This trait has been used to evaluate invasiveness by comparing population abundances and survival rates with native populations (Ebeling et al. 2008; Colautti et al. 2014). Collecting censuses across all major life stages of invasive species and employing demographic models are essential tools for gaining deeper insights into the population dynamics of invasive species (Williams et al. 2010; Colautti et al. 2014).

The study of seaweeds’ invasiveness has predominantly focused on a limited number of species (Dalla Vecchia et al. 2020). The common traits that contribute to their invasiveness are morphology, physiologic plasticity, productivity and reproduction. These are the most studied characteristics (Mabey et al. 2023). Invasive seaweeds often show rapid growth rates and early maturation of reproductive structures, which allows them to quickly establish and spread in new environments (South et al. 2017; Zanolla et al. 2017). Other traits, including multiple overlapping cohorts maintained by both older individuals and new recruits, were described for Undaria pinnatifida (Harvey) Suringar (Schiel and Thompson 2012) and Asparagopsis taxiformis (Delile) Trevisan (Zanolla et al. 2019). This ability to persist year-round in introduced areas has been attributed to their morphological and physiological plasticity (Campbell 1999; Zanolla et al. 2015).

Rugulopteryx okamurae (E.Y. Dawson) I.K. Hwang, W.J. Lee & H.S. Kim (Dictyotaceae, Phaeophyceae), a flattened brown seaweed with dichotomous branching (Hwang et al. 2009), is native to the western Pacific Ocean (Verlaque et al. 2009). Starting as an unprecedented cryptic invasion in the western Mediterranean in 2015, it now forms abundant populations in Spain, including the Canary Islands and the Chafarinas Archipelago (Altamirano et al. 2016; REDEXOS 2022); France (Ruitton et al. 2021); Italy (Bellisimo et al. 2024); Morocco (El Aamri et al. 2018); and Portugal, including the Azores and Madeira (Faria et al. 2021; Bernal-Ibáñez et al. 2022; Liulea et al. 2023). Currently, R. okamurae continues its expansion both in the Mediterranean Sea and the Atlantic Ocean, being most recently recorded in northern Spain, which indicates a worrying upward trend in its spread (Díaz-Tapia et al. 2024). Predictive models suggest that the European and north African coasts are highly favourable habitats for this species, so it is likely that its expansion will continue in the coming years (Muñoz et al. 2019). Like other Dictyotaceae species, R. okamurae exhibits a digenetic isomorphic life cycle, alternating haploid gametophytes with diploid sporophytes. At a morphological level, it presents up to three different morphotypes throughout the year in both native and introduced areas (Sun et al. 2006; Salido and Altamirano 2020). The species can reproduce by propagules, asexually by mitotic monospores and sexually by gametes and meiotic tetraspores, but gametangia have only recently been observed in northern Spain (Díaz-Tapia et al. 2024).

Rugulopteryx okamurae can settle on both horizontal and vertical rocky surfaces within a bathymetric range from eulittoral zones to depths of more than 50 m (García-Gómez et al. 2018; Altamirano et al. 2019). In southern Spain, R. okamurae has seriously impacted native flora, including the seagrass Posidonia oceanica (Linnaeus) Delille (García-Gómez et al. 2018; Junta de Andalucía 2019). Posidonia oceanica is a native seagrass listed in the Spanish List of Wild Species under the Special Protection Regime (Listado Español de Especies Silvestres en Régimen de Protección Especial y Catálogo Español de Especies Amenazadas) (Real Decreto 139/2011). It is also considered to form a priority habitat by the Habitats Directive (Council Directive 92/43/EEC). Similar impacts on native flora have been observed in the Azores Islands (Faria et al. 2022), where R. okamurae is replacing dominant species and altering the structure of shallow-water benthic communities. The fast and abundant proliferation of R. okamurae biomass not only produces high environmental impacts, but also socioeconomic ones (Altamirano et al. 2016; MITECO 2022). In particular, the fishing industry has been suffering economic losses estimated in millions of euros because of a decrease in fish stocks (MITECO 2022). Moreover, the extensive efforts to remove castaway biomass accumulated from beaches in tourist areas is also affecting local administrations (MITECO 2022). All these impacts have led to its inclusion in the List of Invasive Alien Species of Union Concern, being the first seaweed to be added (Commission Implementing Regulation (EU) 2022/1203 of 12 July 2022 amending Implementing Regulation (EU) 2016/1141). However, despite numerous studies about the species’ potential biotechnological uses (e.g. Romero-Vargas et al. (2024)), only few studies focused on its basic biology, which is key to our understanding of the risks associated with its invasion and to lay the groundwork for implementing effective conservation and management strategies for affected marine ecosystems.

In this context, the objective of the present study is to understand the invasiveness of R. okamurae by analysing its population dynamics on an invaded P. oceanica meadow in southern Spain and exploring the relationship of several demographic descriptors with abiotic environmental factors throughout the year.

Methods

Study area and sampling procedure

The study was conducted on a population of R. okamurae established on a P. oceanica meadow located in front of Cambriles Cliff, Granada, Spain (36°44.0033'N, 3°20.6767'W), at a depth of 10 metres (Fig. 1). This P. oceanica meadow is 48 hectares in extent, ranging from 7 metres to 13 metres in depth (Portal Ambiental de Andalucía 2024).

Figure 1.

Sampling site of R. okamurae in Cambriles cliff, Granada (Spain).

Sampling for the study was carried out by scuba every two months from July 2021 to July 2022. Due to the protected status of the P. oceanica habitat (Real Decreto 139/2011) and the status of the invasive species R. okamurae, appropriate permits were obtained from the relevant authorities. Prior to the study, the minimal sampling area of R. okamurae was estimated following Cain and Castro (1959). All algal material in three one-metre by one-metre quadrats, each subdivided into 20 cm × 20 cm squares and spaced five metres apart, was carefully scraped and bagged separately in plastic bags. Samples were transported to the laboratory for weighing and for estimating the minimal sampling area through iterative measurements of contiguous squares. The process involved analysing the abundance data to estimate the smallest square size that showed no statistical differences in abundance between contiguous squares and through the estimation of the aggregation index, based on the ratio of variance to mean biomass, as described by Blackmann (1942). The minimum sampling area was 30 cm × 30 cm. Thus, within three areas with a homogeneous biomass density of R. okamurae, four replicates of the minimal area were taken and collected at the same depth, separated by 5 m between sampling areas. All algal material was thoroughly removed avoiding damaging P. oceanica plants and placed in plastic bags, preventing accidental dispersal of R. okamurae. Subsequently, the samples were transported to the laboratory under cold and dark conditions for further analysis.

Population dynamics

All individuals within each sample were measured from the base of the thallus to the most distal dichotomy and categorised into eight different size classes in two-centimetre increments, except the first class (four centimetres), which corresponds to young individuals or recruits. The abundance of R. okamurae was estimated as the number of individuals and their dry weight (DW) biomass per square metre for each sampling event. Dry weight biomass was quantified after drying at 60 °C for 48 hours in an oven. Distribution frequency of size classes plots was constructed using the percentage of individuals in each size class against the total number of individuals in each replicate and the mean values calculated. Cohorts, defined as individuals sharing a particular event during their lifespan (Crisp 1971), were identified by tracking the displacement of Gaussian curves in the frequency versus time histograms, based on distributions of size classes along the year (Aranda et al. 1984; Zanolla et al. 2019). If more than one cohort overlapped within a given sampling time, a separation was achieved through a numerical fit using Gaussian curves (Petersen 1912; Zanolla et al. 2019). Each cohort was represented using all replicates by a Gaussian curve characterised by its mean and standard deviation (SD) through the Distr. Norm function. The adjustment of these parameters was performed using the SOLVER application (Excel, Windows Office 365). Finally, a statistical comparison of each mean and SD from the curves was carried out using a t-test statistical analysis. If no statistical differences were found, we assumed the existence of a single cohort. The model’s significance was validated through R2.

Environmental factors

Daily data on average temperature, maximum and minimum temperature and salinity from a buoy situated near the study area were obtained from Puertos del Estado website (Ministerio de Transportes y Agenda Urbana), (SIMAR 2044080; 36.67°N, 3.5°E). The photoperiod was provided by the Observatorio Astronómico Nacional website, located in Granada (Instituto Geográfico Nacional, Ministerio de Fomento). In addition, three seawater samples were collected on each sampling date at the same depth as the R. okamurae population and analysed using an automated nutrient analyser QuAAtro AQ2 AACE (Seal Analytical Ltd. Fareham, UK) for ammonium (Slawyk and MacIsaac 1972), nitrate and nitrite (Shinn 1941; Wood et al. 1967) and phosphate (Murphy and Riley 1962) following standard procedures.

Statistical analysis

Each biological variable was analysed using a one-way model ANOVA (P < 0.05), with time as a fixed factor. Homoscedasticity and normality were tested prior to the ANOVA by Levene’s and Saphiro-Wilks test, respectively. When significant differences were found for a given biological variable, the Student-Newman-Keuls test (SNK) was applied for post hoc comparisons. Statistical analyses of ANOVA were carried out by SigmaPlot 11.0 software (Systat Software Inc., Chicago, IL, USA). To assess the relationship between environmental factors and biological variables, a principal component analysis (PCA) was conducted. An environmental space was constructed with axes derived from PCA using environmental factors and plotted in InfoStat version 2008 (Di Rienzo et al. 2008). Subsequently, Pearson correlations were used to investigate the relationship between biotic variables (density of individuals and biomass of R. okamurae) and the two principal axes derived from PCA (PC 1 and PC 2). These environmental axes produced in PCA were computed with the software PAST (Hammer et al. 2001).

Results

Population dynamics

Rugulopteryx okamurae was present throughout the whole year, displaying significant differences (DF: 6, F-value = 6.5, P = 0.002, Appendix 1) in density of individuals for different months (Fig. 2A). The highest values were recorded in September and November, with an average value of 3285 individuals per square metre, whereas the lowest values were observed in January, March and May, with a decrease of more than one-third (888 individuals per square metre) (Fig. 2A). There were also significant monthly differences in biomass (Fig. 2B, DF: 6, F-value = 60.6, P < 0.001, Appendix 1). These values increased from November to July 2022, when they were 14 times higher.

Figure 2.

Density of R. okamurae throughout the study period (July 2021–July 2022) referred as A individuals (103 individuals per square metre) and B biomass (g DW per square metre). Data are expressed as mean ± SD (n = 3). Different letters denote significant differences among months following ANOVA results (P < 0.05).

The distribution of size classes varied throughout the year, reaching 18 cm in May and July 2022, representing less than 3% of the total individuals in those months. In contrast, in November, individuals did not exceed six centimetres (Figs 3, 4). Over 70% of the total population belonged to the size class of 0–4 cm in September and November (Figs 3, 4). However, there was an 85% reduction in the density of the smallest size class from November to January and it continued to decrease until May, when the lowest numbers were recorded (Figs 3, 4).

Figure 3.

Density (individuals per square metre) of each size class throughout the study period (July 2021 to July 2022). Data are expressed as mean ± SD (n = 3).

Figure 4.

Frequency of individuals of each size class (%) related to the total of individuals of R. okamurae in each sampling event. Discontinuous lines indicate the six cohorts (labelled I to VI) found throughout the study (July 2021–July 2022). Values are expressed as means ± SD (n = 3).

Cohorts

The frequency distribution of thallus size classes followed a normal distribution, revealing the presence of six cohorts (named using Roman numerals I to VI) with different longevity and temporal distribution (Fig. 4). The longest cohort started in November and persisted until May (cohort IV), while cohorts I, III, V and VI had a lifespan of one month each (Fig. 4). Cohort II was present in July and September 2021. The composition of these cohorts also showed variations (Fig. 4). Cohort III was composed solely of young individuals, whereas cohorts IV and VI included individuals from different size classes, including individuals up to 18 cm (Fig. 4). However, it took six months for plants in cohort IV to achieve the largest size class, but only two months for those in cohort VI (Fig. 4). Additionally, overlapping cohorts were observed in July 2021 (cohorts I–II), September (II–III) and May (V–VI) (Fig. 4).

Environmental factors

Average temperature ranged from 15 °C to 23 °C, with the maximum temperature registered in August 2021 (25.1 °C) and the minimum (14.8 °C) in January (Fig. 5A). Low variability (> 1 °C) was registered in the winter months but, in summer 2021, it reached a difference of 8 °C (Fig. 5A). The photoperiod varied within a range of 10 to 14 hours (Fig. 5B). Salinity remained relatively stable during the study period, ranging from 36.8 to 37.3 (Fig. 5B).

Figure 5.

Environmental factors during the study period (July 2021–July 2022) in the study area (buoy SIMAR 2044080) A average data of mean, maximum and minimum temperature B monthly average photoperiod (dotted line) and salinity levels (solid line). Data expressed as mean ± SD of daily data.

Nitrate concentration during the study period ranged from < 0.2 µM (the detection limit of the analytical method) to 1.4 µM (Fig. 6). Ammonium concentration did not display significant seasonal variations, oscillating between 0.5 and 1.6 µM (Fig. 6). Concentrations of nitrite and phosphate were too low to be detected (< 0.2 µM).

Figure 6.

Nitrate (solid line) and ammonium (dotted line) concentration during the study period (July 2021–July 2022) in the study area. Data expressed as mean ± SD (n = 3).

Principal component analysis (PCA) of environmental factors indicated that specific abiotic factors contributed differently along the two main axes (PC 1 and PC 2), explaining 79.8% of the total variance (Fig. 7, Appendix 2). Temperature, photoperiod and salinity contributed mainly for the PC 1, while nitrate and ammonium concentrations were the primary factors contributing for the PC 2 (Fig. 7). Pearson analysis between biological variables (density of individuals and biomass) and the two principal components derived from the PCA showed no significant correlations (P values > 0.05, Appendix 3).

Figure 7.

Environmental space derived from PCA constructed by environmental factors measured in the study area. Cohorts were labelled as I to VI. T, average temperature; Max T, maximum average temperature; Min T, minimum average temperature; S, salinity; LH, Photoperiod; NO3-, concentration of nitrate; NH4+, concentration of ammonium; PC, principal component.

Discussion

Population dynamics of R. okamurae on a P. oceanica meadow during the year explains its invasive behaviour, characterised by the constant presence of a high density of individuals. This is achieved through its ability to continuously produce recruits and the presence of year-round short-lived overlapping cohorts that are not affected by the environmental parameters of the area.

The seaweeds’ continuous presence is a characteristic of invasiveness shared amongst invasive seaweed species in invaded areas. Examples include the red seaweeds Womersleyella setacea (Hollenberg) R.E. Norris (Cebrian and Rodríguez-Prieto 2012), Asparagopsis armata Harvey (Aranda et al. 1984) and A. taxiformis (Zanolla et al. 2019), the brown seaweeds Sargassum muticum (Yendo) Fensholt (Thomsen et al. 2006) and U. pinnatifida (Schiel and Thompson 2012) and the green seaweed genus Caulerpa, specifically C. taxifolia (M. Vahl) C. Agardh (de Villèle and Verlaque 1995) and C. cylindracea Sonder (Ruitton et al. 2005). All these species showed the ability to maintain stable populations throughout the year, suggesting the existence of adaptive physiological mechanisms to the newly-invaded area and posing a continuous impact over native communities. Likewise, R. okamurae demonstrated a remarkable ability to establish dense monospecific populations, maintaining densities exceeding 3000 individuals per square metre and reaching biomass up to 170 g DW per square metre. These values greatly surpass other documented invasive seaweeds. For instance, A. taxiformis in the southern Iberian Peninsula, reached the highest biomass in the summer (Zanolla et al. 2017), but these values were 40% lower than those recorded for R. okamurae in this study. Another example was U. pinnatifida with a maximum density of about 100 individuals per square metre in affected areas in New Zealand (Schiel and Thompson 2012) or S. muticum with a maximum of 1000 individuals per square metre in Spain (Arenas and Fernández 2000). These are 97% and 67% lower than for R. okamurae, respectively, but this can be explained by species size differences. Density peaks of individuals were recorded in late summer and autumn, while biomass peaks were observed in late spring and summer. This pattern suggests a seasonal variation in the species’ invasiveness in P. oceanica meadows, with space occupation influenced either by the high influx of new recruits or the substantial biomass of larger adult individuals. A similar pattern in temporal fluctuation of invasive behaviour was observed in different invasive seaweeds, such as S. muticum, showing a population dynamic characterised by alternating periods dominated by either few large individuals or numerous small individuals (Arenas and Fernández 2000).

The lack of significant correlation between abiotic factors, density of individuals and biomass is useful for understanding the success and establishment of invasive species (Jehangir et al. 2024). This result suggests that R. okamurae can thrive under diverse physical-chemical conditions, tolerating a wide range of abiotic environmental factors, which enhances invasiveness (Alpert et al. 2000). A similar pattern was documented in A. taxiformis, which maintains a high productive succession of cohorts throughout the year, regardless of changing environmental conditions (Zanolla et al. 2019). During the study period, R. okamurae exhibited continuous recruitment, with peaks of young individuals observed in summer and autumn. Although the specific mechanisms involved were not detailed, the presence of propagules – diminute proliferous branches developed from cortical cells of the parent blade (Kajimura 1992) – was noticeable throughout the year, except in January (pers. obs.). This ability, coupled with its capacity for clonal multiplication through fragmentation (Rosas-Guerrero et al. 2020), facilitates ongoing recruitment regardless of abiotic factors. This mechanism not only adds to its continuous presence in the P. oceanica meadow, but also contributes to the release of new individuals into the water column. These unattached thalli can be detected at a depth of 1141 metres, where they remain photosynthetically active after light exposure (Mateo-Ramírez et al. 2023). This agrees with the propagule pressure invasion hypothesis, which emphasises the continuous input of new individuals into an invaded area (Simberloff 2009).

The identification of six successive cohorts throughout the year revealed the specific ecological strategies of R. okamurae for maintaining its population in the introduced area. The presence of different cohorts under varying environmental conditions suggests distinct environmental requirements for each cohort, a trait shared with other invasive species, such as A. taxiformis (Zanolla et al. 2019) and U. pinnatifida (Schiel and Thompson 2012). The population’s prevalence is driven by constant recruitment and continuous generational succession. Most cohorts had a brief duration of about two months, except cohort IV, which lasted for six months. This finding suggests that the species is unable to complete the typical haplodiplontic digenetic life cycle exhibited in its native area, which usually spans more than two years and involves an annual alternation between gametophytic and sporophytic phases (Agatsuma et al. 2005). Although reproductive characterisation is still being processed, the continuous presence of propagules and monospores (pers. obs.) may indicate that the population is sustained primarily through asexual reproduction and vegetative multiplication, with sporophytes playing a key role in population maintenance. Furthermore, overlapping cohorts observed in some months allows for a continuous succession of individuals and sustained occupation of space resulting in high coverage throughout the year. This pattern, like that observed for U. pinnatifida in New Zealand, enhances the invasion success of the species and increases propagule pressure in invaded areas (Schiel and Thompson 2012; South et al. 2017). This continuous occupation can have more impact on the structure and function of the ecosystem than intermittent periods of high and low coverage (Schiel and Thompson 2012) and is likely to disrupt native species and alter ecosystem dynamics, contributing to its invasive potential.

The identified cohorts might be linked to the occurrence of the species’ morphotypes described in both native and introduced areas, characterised by variations in thallus thickness and width and the number of dichotomies (Sun et al. 2006; Salido and Altamirano 2020). These morphotypes appear to follow a seasonal pattern. For example, in introduced areas, the thick morphotype arises in winter and the thinner one in summer, alternating with the intermediate morphotype (Salido and Altamirano 2020). Therefore, cohort IV could align with the thicker morphotype, cohorts II and VI with the thinner morphotype and cohorts I, III and V with the intermediate morphotype. This potential association between morphotypes and cohorts provides valuable insights into the adaptability of the different morphotypes of R. okamurae in response to changing environmental conditions, which contributes to understanding the species’ dynamics in marine ecosystems. Although this study did not explore the relationship between morphotypes and abiotic factors, identifying the environmental drivers behind morphotype development would be essential for understanding the species’ invasion success and its impact on ecosystems. Theoretical frameworks exploring the connection between phenotypic plasticity and invasion success propose that invasive species could maintain high fitness across diverse environmental conditions owing to their physiological plasticity and can even thrive under unfavourable conditions (Richards et al. 2006; Gioria et al. 2023), which may explain the observed patterns in R. okamurae. Understanding this behaviour could deepen knowledge about the species’ invasiveness and the factors influencing morphotype succession throughout the year.

Conclusion

Overall, R. okamurae exhibited a notable ability to persist in the P. oceanica meadow despite fluctuating environmental conditions. Its capacity to endure under a wide range of abiotic factors highlights its invasiveness, which is facilitated through population dynamics. This invasiveness observed in this specific habitat is further enhanced by continuous recruitment and a succession of distinct cohorts. This study also reveals temporal windows of invasiveness for R. okamurae, driven by an intense density of new small-sized individuals during summer and early autumn, while late spring and summer are characterised by high biomass accumulation of larger-sized individuals. These findings contribute to a deeper understanding of the invasiveness of R. okamurae in the P. oceanica meadows through population dynamics, underscoring its ability to dominate space and persist in a variety of environmental conditions. Identifying these patterns offers crucial insights into the success of R. okamurae and can guide effective management strategies.

Acknowledgements

The authors thank to Gianluca Nania and the other volunteers for their support in the measurements of thalli. We further thank comments provided by reviewers.

Additional information

Conflict of interest

The authors have declared that no competing interests exist.

Ethical statement

No ethical statement was reported.

Funding

Proyecto RUGULOPTERYX (Fundación Biodiversidad, 2020), FEDERJA-006 (FEDER-Junta de Andalucía); TED2021-130080B-I00 (Ministerio de Ciencia e Innovación, NextGenerationEU); RNM 262 Biogeography, Diversity and Conservation Research Team, University of Málaga.

Author contributions

Conceptualization: JRG, JDR, RC, MA. Data curation: JRG. Formal analysis: RC, MZ, JRG. Funding acquisition: JDR, RC, MA. Investigation: MA, RC, JRG, JDR, MZ. Methodology: JRG, MZ, JDR, MA, RC. Project administration: MA. Resources: MA. Supervision: MA. Validation: MA. Visualization: JRG. Writing – original draft: JRG. Writing – review and editing: JDR, RC, MZ, JRG, MA.

Author ORCIDs

Jesús Rosas-Guerrero https://orcid.org/0000-0001-6042-8031

Raquel Carmona https://orcid.org/0000-0002-9656-3195

Julio De la Rosa https://orcid.org/0000-0002-4402-4405

Marianela Zanolla https://orcid.org/0000-0001-9585-5906

María Altamirano https://orcid.org/0000-0003-0912-3704

Data availability

All of the data that support the findings of this study are available in the main text.

References

  • Alpert P, Bone E, Holzapfel C (2000) Invasiveness, invasibility and the role of environmental stress in the spread of non-native plants. Perspectives in Plant Ecology, Evolution and Systematics 3(1): 52–66. https://doi.org/10.1078/1433-8319-00004
  • Altamirano M, De La Rosa J, Martínez FJ (2016) Arribazones de la especie exótica Rugulopteryx okamurae (E.Y. Dawson) I.K. Wang, W. J. Lee and H.S. Kim (Dictyotales, Ochrophyta) en el Estrecho de Gibraltar: Primera cita para el Atlántico y España. Algas 52: 20.
  • Altamirano M, De La Rosa J, Carmona R, Zanolla M, Muñoz AR (2019) Macroalgas invasoras en las costas andaluzas. Algas 55e: 10–13.
  • Aranda J, Niell FX, Fernández JA (1984) Production of Asparagopsis armata (Harvey) in a thermally-stressed intertidal system of low tidal amplitude. Journal of Experimental Marine Biology and Ecology 84(3): 285–295. https://doi.org/10.1016/0022-0981(84)90186-2
  • Bellisimo G, Altamirano M, Munoz AR, De la Rosa J, Hung TH et al. (2024) The invasive brown seaweed Rugulopteryx okamurae (Dictyotales, Ochrophyta) continues to expand: First record in Italy. BioInvasions Records 13(2): 385–401. https://doi.org/10.3391/bir.2024.13.2.08
  • Bernal-Ibáñez A, Chebaane S, Sempere-Valverde J, Faria J, Ramalhosa P, Kaufmann M, Florido M, Albert-Fonseca A, Canning-Clode J, Gestoso I, Cacabelos E (2022) A worrying arrival: The first record of brown macroalga Rugulopteryx okamurae in Madeira Island and its invasive risk. BioInvasions Records 11(4): 912–924. https://doi.org/10.3391/bir.2022.11.4.10
  • Blackmann GE (1942) Statistical and ecological studies in the distribution of species in plant communities. I Dispersion as a factor in the study of changes in plant populations. Annals of Botany 6(2): 351–370. https://doi.org/10.1093/oxfordjournals.aob.a088411
  • Campbell SJ (1999) Uptake of ammonium by four species of macroalgae in Port Phillip Bay, Victoria, Australia. Marine and Freshwater Research 50(6): 515–522. https://doi.org/10.1071/MF98058
  • Catford JA, Vesk PA, Richardson DM, Pyšek P (2012) Quantifying levels of biological invasion: Towards the objective classification of invaded and invasible ecosystems. Global Change Biology 18(1): 44–62. https://doi.org/10.1111/j.1365-2486.2011.02549.x
  • Crisp DJ (1971) Energy flow measurements. In: Holme A, McIntyre AD (Eds) Methods for the Study of Marine Benthos. Blackwell, Oxford, 197–279.
  • de Villèle X, Verlaque M (1995) Changes and degradation in a Posidonia oceanica bed invaded by the introduced tropical alga Caulerpa taxifolia in the northwestern Mediterranean. Botanica Marina 38(1–6): 79–87. https://doi.org/10.1515/botm.1995.38.1-6.79
  • Di Rienzo JA, Casanoves F, Balzarini MG, Gonzalez L, Tablada MRCW, Robledo CW (2008) InfoStat, version 2008. Group Infostat, fca, Universidad Nacional de Córdoba, Argentina, 115.
  • Díaz-Tapia P, Alvite N, Bañón R, Barreiro R, Barrientos S et al. (2024) Multiple Introduction Events Expand the Range of the Invasive Brown Alga Rugulopteryx okamurae to Northern Spain. https://doi.org/10.2139/ssrn.4884141
  • El Aamri F, Idhalla M, Tamsouri MN (2018) Ocurrence of the invasive brown seaweed Rugulopteryx okamurae (EY Dawson) IK Wang, WJ Lee and HS Kim (Dictyotales, Phaeophyta) in Morocco (Mediterranean Sea). Mediterranean Fisheries and Aquaculture Research 1(2): 92–96.
  • Enders M, Hütt MT, Jeschke JM (2018) Drawing a map of invasion biology based on a network of hypotheses. Ecosphere 9(3): e02146. https://doi.org/10.1002/ecs2.2146
  • Faria J, Prestes AC, Moreu I, Martins GM, Neto AI, Cacabelos E (2021) Arrival and proliferation of the invasive seaweed Rugulopteryx okamurae in NE Atlantic islands. Botanica Marina 65(1): 45–50. https://doi.org/10.1515/bot-2021-0060
  • Faria J, Prestes AC, Moreu I, Cacabelos E, Martins GM (2022) Dramatic changes in the structure of shallow-water marine benthic communities following the invasion by Rugulopteryx okamurae (Dictyotales, Ochrophyta) in Azores (NE Atlantic). Marine Pollution Bulletin 175: 113358. https://doi.org/10.1016/j.marpolbul.2022.113358
  • García-Gómez JC, Sempere-Valverde J, Ostalé-Valriberas E, Martínez M, Olaya L et al. (2018) Rugulopteryx okamurae (EY Dawson) IK Hwang, WJ Lee & HS Kim (Dictyotales, Ochrophyta), alga exótica "explosiva" en el estrecho de Gibraltar. Observaciones preliminares de su distribución e impacto. Almoraima : Revista de Estudios Campogibraltareños 49: 97–113.
  • Hammer Ø, Harper DAT, Ryan P (2001) PAST: Paleontological statistics software package for education and data analysis. Palaeontologia Electronica 4: 9.
  • Hwang IK, Lee WJ, Kim HS, De Clerck O (2009) Taxonomic Reappraisal of Dilophus Okamurae (Dictyotales, Phaeophyta) from The Western Pacific Ocean. Phycologia 48(1): 1–12. https://doi.org/10.2216/07-68.1
  • Jehangir S, Khan SM, Ahmad Z, Ejaz U, Ain QU et al. (2024) Distribution of the Cannabis sativa L. in the Western Himalayas: A tale of the ecological factors behind its continuous invasiveness. Global Ecology and Conservation 49: e02779. https://doi.org/10.1016/j.gecco.2023.e02779
  • Kajimura M (1992) Vegetative propagation in Dilophus okamurae and Zonaria flabellata (Dictyotaceae, Phaeophyta). Shimane University 26: 95–106.
  • Liulea S, Serrão EA, Santos R (2023) Spread and impact of the invasive brown algae Rugulopteryx Okamurae on the Algarve Coast, Southern Portugal (Ne Atlantic). https://doi.org/10.2139/ssrn.4446622
  • Lowry E, Rollinson EJ, Laybourn AJ, Scott TE, Aiello‐Lammens ME, Gray SM, Mickley J, Gurevitch J (2013) Biological invasions: A field synopsis, systematic review, and database of the literature. Ecology and Evolution 3(1): 182–196. https://doi.org/10.1002/ece3.431
  • Mateo-Ramírez Á, Iñiguez C, Fernández-Salas LM, Sánchez-Leal RF, Farias C, Bellanco MJ, Gil J, Rueda JL (2023) Healthy thalli of the invasive seaweed Rugulopteryx okamurae (Phaeophyceae) being massively dragged into deep-sea bottoms by the Mediterranean Outflow Water. Phycologia 62(2): 99–108. https://doi.org/10.1080/00318884.2023.2177057
  • Muñoz AR, Martín-Taboada A, De La Rosa J, Carmona R, Zanolla M, Altamirano M (2019) La modelación de la distribución de especies como herramienta en la gestión de invasiones biológicas en el medio marino: El caso de Rugulopteryx okamurae (Dictyotaceae, Ochrophyta) en el mediterráneo. Algas 55e: 37–41.
  • Petersen CJ (1912) Über Menge und Jahresproduktion der Benthospflanzen an den westeuropäischen Küsten. Internationale Revue der Gesamten Hydrobiologie 5(1): 47–52. https://doi.org/10.1002/iroh.19120050104
  • Pyšek P, Jarošík V, Pergl J, Moravcová L, Chytrý M, Kühn I (2014) Temperate trees and shrubs as global invaders: The relationship between invasiveness and native distribution depends on biological traits. Biological Invasions 16(3): 577–589. https://doi.org/10.1007/s10530-013-0600-2
  • Richards CL, Bossdorf O, Muth NZ, Gurevitch J, Pigliucci M (2006) Jack of all trades, master of some? On the role of phenotypic plasticity in plant invasions. Ecology Letters 9(8): 981–993. https://doi.org/10.1111/j.1461-0248.2006.00950.x
  • Richardson DM, Pyšek P (2006) Plant invasions: Merging the concepts of species invasiveness and community invasibility. Progress in Physical Geography 30(3): 409–431. https://doi.org/10.1191/0309133306pp490pr
  • Romero-Vargas A, Cala K, Blandino A, Díaz AB (2024) Bioconversion of the invasive seaweed Rugulopteryx okamurae into enzymes and polyhydroxyalkanoates. Algal Research 81: 103587. https://doi.org/10.1016/j.algal.2024.103587
  • Rosas-Guerrero J, Carmona R, Altamirano M (2020) Efecto de la temperatura y la irradiancia sobre el crecimiento, la propagación vegetativa y la actividad fotosintética del alga invasora Rugulopteryx okamurae (Dictyotales, Ochrophyta). Algas 56: 106.
  • Ruitton S, Verlaque M, Boudouresque CF (2005) Seasonal changes of the introduced Caulerpa racemosa var. cylindracea (Caulerpales, Chlorophyta) at the northwest limit of its Mediterranean range. Aquatic Botany 82(1): 55–70. https://doi.org/10.1016/j.aquabot.2005.02.008
  • Ruitton S, Blanfune A, Boudouresque CF, Guillemain D, Michotey V, Roblet S, Thibault D, Thibaut T, Verlaque M (2021) Rapid spread of the invasive brown alga Rugulopteryx okamurae in a National Park in Provence (France, Mediterranean Sea). Water 13(16): 2306. https://doi.org/10.3390/w13162306
  • Salido M, Altamirano M (2020) Variabilidad temporal de la morfología e invasividad de Rugulopteryx okamurae (Dictyotales, Ochrophyta) en el Estrecho de Gibraltar. Algas 56: 101.
  • Schiel DR, Thompson GA (2012) Demography and population biology of the invasive kelp Undaria pinnatifida on shallow reefs in southern New Zealand. Journal of Experimental Marine Biology and Ecology 434: 25–33. https://doi.org/10.1016/j.jembe.2012.07.023
  • South PM, Floerl O, Forrest BM, Thomsen MS (2017) A review of three decades of research on the invasive kelp Undaria pinnatifida in Australasia: An assessment of its success, impacts and status as one of the world’s worst invaders. Marine Environmental Research 131: 243–257. https://doi.org/10.1016/j.marenvres.2017.09.015
  • Sun Z, Hasegawa K, Tanaka J (2006) A Morphological Study of Dilophus okamurae (Dictyotales, Phaeophyceae) in Japan. Shokubutsu Kenkyu Zasshi 81: 235–244.
  • Thomsen MS, Wernberg T, Stæhr PA, Pedersen MF (2006) Spatio-temporal distribution patterns of the invasive macroalga Sargassum muticum within a Danish Sargassum-bed. Helgoland Marine Research 60(1): 50–58. https://doi.org/10.1007/s10152-005-0016-1
  • Verlaque M, Steen F, De Clerck O (2009) Rugulopteryx (Dictyotales, Phaeophyceae), a genus recently introduced to the Mediterranean. Phycologia 48(6): 536–542. https://doi.org/10.2216/08-103.1
  • Williams JL, Auge H, Maron JL (2010) Testing hypotheses for exotic plant success: Parallel experiments in the native and introduced ranges. Ecology 91(5): 1355–1366. https://doi.org/10.1890/08-2142.1
  • Wood ED, Armstrong FAJ, Richards FA (1967) Determination of nitrate in sea water by Cadmium-copper reduction to nitrtite. Journal of the Marine Biological Association of the United Kingdom 47(1): 23–31. https://doi.org/10.1017/S002531540003352X
  • Zanolla M, Altamirano M, Carmona R, De La Rosa J, Sherwood A, Andreakis N (2015) Photosynthetic plasticity of the genus Asparagopsis (Bonnemaisoniales, Rhodophyta) in response to temperature: Implications for invasiveness. Biological Invasions 17(5): 1341–1353. https://doi.org/10.1007/s10530-014-0797-8
  • Zanolla M, Carmona R, Altamirano M (2017) Reproductive ecology of an invasive lineage 2 population of Asparagopsis taxiformis (Bonnemaisoniales, Rhodophyta) in the Alboran Sea (western Mediterranean Sea). Botanica Marina 60(6): 627–638. https://doi.org/10.1515/bot-2017-0056
  • Zanolla M, Altamirano M, Niell FX, Carmona R (2019) There is more than meets the eye: Primary production of the invasive seaweed Asparagopsis taxiformis (Bonnemaisoniaceae, Rhodophyta) is provided by six cohorts with distinctive characteristics. Aquatic Botany 153: 24–28. https://doi.org/10.1016/j.aquabot.2018.11.007

Appendix 1

Table A1.

One-way (month) ANOVA for density (individuals x 103 per square metre) of R. okamurae and biomass (g DW per square metre). Asterisks represent significant differences (P < 0.05).

variable df MS F-value P value
One-way ANOVA Density 6 3761825 6.49 0.002*
Biomass 6 9291 60.61 < 0.001*

Appendix 2

Table A2.

Loading factors and percentage of the variance of the two first PCA.

PC 1 PC 2
Percentage of variance 53.2% 26.6%
Cumulative percentage 53.2% 79.8%
Variable PC 1 PC 2
Loading factors
T 0.51 -0.05
Max T 0.49 0.10
Min T 0.49 -0.14
S -0.30 0.01
LH 0.41 0.05
NO3 -0.01 0.70
NH4 0.04 0.69

Appendix 3

Table A3.

Correlation analysis between two principal components and biological variables measured (density and biomass) in the population of R. okamurae (n = 7).

PC 1 PC 2
Variable r P value r P value
Density 0.57 0.18 0.24 0.61
Biomass 0.39 0.36 0.09 0.85
login to comment