Research Article |
Corresponding author: Jesús Rosas-Guerrero ( jrosas@uma.es ) Academic editor: Michael McKinney
© 2025 Jesús Rosas-Guerrero, Raquel Carmona, Julio De la Rosa, Marianela Zanolla, María Altamirano.
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:
Rosas-Guerrero J, Carmona R, De la Rosa J, Zanolla M, Altamirano M (2025) A population dynamics approach to understand the invasiveness of the seaweed Rugulopteryx okamurae (Ochrophyta, Dictyotales). NeoBiota 97: 1-17. https://doi.org/10.3897/neobiota.97.137873
|
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.
Demography, density, generation, macroalga, Posidonia oceanica, recruitment
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 (
The study of seaweeds’ invasiveness has predominantly focused on a limited number of species (
Rugulopteryx okamurae (E.Y. Dawson) I.K. Hwang, W.J. Lee & H.S. Kim (Dictyotaceae, Phaeophyceae), a flattened brown seaweed with dichotomous branching (
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 (
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.
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.
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
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 (
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 (
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 (
Rugulopteryx okamurae was present throughout the whole year, displaying significant differences (DF: 6, F-value = 6.5, P = 0.002, Appendix
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
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).
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.
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.
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.
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.
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.
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 (
The lack of significant correlation between abiotic factors, density of individuals and biomass is useful for understanding the success and establishment of invasive species (
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 (
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 (
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.
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.
The authors have declared that no competing interests exist.
No ethical statement was reported.
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.
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.
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
All of the data that support the findings of this study are available in the main text.
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* |
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 |