Invasion success of a Lessepsian symbiont

Scientific Reports volume 13, Article number: 12578 (2023) Cite this article

1 Altmetric

Metrics details

Among the most successful Lessepsian invaders is the symbiont-bearing benthic foraminifera Amphistegina lobifera. In its newly conquered habitat, this prolific calcifier and ecosystem engineer is exposed to environmental conditions that exceed the range of its native habitat. To disentangle which processes facilitated the invasion success of A. lobifera into the Mediterranean Sea we analyzed a ~ 1400 bp sequence fragment covering the SSU and ITS gene markers to compare the populations from its native regions and along the invasion gradient. The genetic variability was studied at four levels: intra-genomic, population, regional and geographical. We observed that the invasion is not associated with genetic differentiation, but the invasive populations show a distinct suppression of intra-genomic variability among the multiple copies of the rRNA gene. A reduced genetic diversity compared to the Indopacific is observed already in the Red Sea populations and their high dispersal potential into the Mediterranean appears consistent with a bridgehead effect resulting from the postglacial expansion from the Indian Ocean into the Red Sea. We conclude that the genetic structure of the invasive populations reflects two processes: high dispersal ability of the Red Sea source population pre-adapted to Mediterranean conditions and a likely suppression of sexual reproduction in the invader. This discovery provides a new perspective on the cost of invasion in marine protists: The success of the invasive A. lobifera in the Mediterranean Sea comes at the cost of abandonment of sexual reproduction.

Biological invasions driven by climate change are currently profoundly modifying ecological landscapes1,2. Unlike normal range extensions, where species are largely tracking their climatic envelope, invasive species conquer entirely new spaces, with a higher probability of facing climatic (seasonality), physical (light), chemical (salinity) or biotic (microbiome and interactome) conditions that exceed the range they experienced within their native habitat. In this context, it is important to understand how a given species can become a successful invader. The challenge of exposure to foreign conditions in the newly conquered space could be counteracted by adaptations. In this scenario, successful invaders can be expected to display a high adaptive potential3,4. Alternatively, the native population could already possess the key adaptations, for example as a result of its evolutionary history5,6, past migration events7, or ecological filtering8.

A remarkable biological invasion phenomenon known as the Lessepsian invasion has taken place in the Mediterranean Sea since 1869. The opening of the Suez Canal in that year ignited a dramatic and largely uni-directional migration of Indo-Pacific marine species into the Mediterranean. So far, over 600 invasive marine species have been reported in the eastern Mediterranean9,10, with more new invaders appearing as the ongoing warming makes the Levantine Basin more tropics-like11. Among the particularly successful invaders are symbiont-bearing larger benthic foraminifera (LBF). LBFs inhabit shallow coastal waters, where they commonly live attached to algae or hard substrate12. Foraminifera have limited capacity for active movement during life, but adult specimens can be passively suspended and transported by currents13 and the passive mobility of the minute juveniles, propagules, or flagellated gametes produced by sexual reproduction is likely even larger14,15. In addition, foraminiferal dispersal mediated by travel in the digestive system of fish was documented in the Mediterranean14,14. This combination of multiple mechanisms of dispersal results in broad species ranges and lack of regional population differentiation17.

The most successful Lessepsian invader LBF is the diatom-bearing Amphistegina lobifera Larsen 1976, whose prolifically calcifying invasive populations modify the sediment substrate and displace native species18,19. After establishing itself in the Levantine Basin, its invasion towards the west accelerated in the last few decades18,20,21,22,23,24,25,26 and the species now expanded its range to Sicily27. The invasion appears to be sourced entirely from within the western Indian Ocean genotype Ia of the species17, but it remains unknown how the genetic diversity of the invading populations is structured along the invasion gradient.

Theoretically, an invading population could show a reduced genetic diversity compared to the source population because only a fraction of the source population participated in the invasion (founder effect28,29,30) or because the source population is highly structured and only one subpopulation possessed adaptations allowing it to invade31,32,33. The combination of a founder effect and exposure to the new environment could also lead to rapid emergence of new variation in the invasive population (e.g., as observed in a Lessepsian cornetfish34). In both cases, the reduced variability and increased divergence should show a gradient along the invasion progression35,36, with most severe effects visible at the invasion front. Alternatively, the source population could have already possessed the necessary adaptations to the new environment and in the presence of a large dispersal potential17, the opening of the Suez Canal could have merely removed a physical barrier after which the invader rapidly fills the free space without genetic differentiation.

Here we investigated the population structure of the invasive A. lobifera between the source population and populations representing different stages of the invasion. The analysis is based on a ~ 1400 bp long sequence fragment of the rRNA gene complex covering the end of the SSU coding region and the adjacent internal transcribed spacer (ITS) region. In foraminifera, the SSU rRNA gene contains fast-evolving variable regions, which provide a resolution within species37,38,39,40 and the ITS rRNA gene region provides an even higher resolution given its higher mutation rate41. Therefore, this marker should allow us to detect genetic differentiation (or a lack thereof) among the invasive populations of A. lobifera and thus constrain to what degree the invasion success of the species is due to novel adaptive change or preexisting adaptations.

To characterize the genetic variability of Amphistegina lobifera along the invasion gradient, we sampled populations of the species in Sicily, where the most recent invasion front has been identified by Guastella et al.27, and in Israel where populations representative of the source population (Red Sea) and pioneer invaders (Eastern Mediterranean Sea) could be collected (Fig. 1). During the sampling in Sicily in September 2019, we first carried out an exploratory survey to identify the position of the invasion front compared the observations of Guastella et al.27 that were done between 2015 and 2017. We focused our effort on the eastern coast of Sicily, where the invasion front was located by Guastella et al.27 between Capo Passero to the South and Brucoli to the North. We revisited these two locations and sampled two additional locations in between (Arenella and Plemmirio) and two additional locations North of Brucoli (Cannizzaro and Recanati), assuming that the invasion front may have moved further North since 2017. At each location, pebbles and macroalgae were collected from the depth of 0.5–5 m by snorkeling. The collected substrates were brushed, and the recovered material was sieved at 63–500 µm and transferred in 0.5 L jars that were filled with ambient seawater. We then examined the samples under a stereomicroscope and assessed qualitatively the presence of A. lobifera in the samples.

Occurrence of invasive Amphistegina in the Mediterranean Sea and sites sampled for genetic analysis (a, and e for zoom in the Sicilian sites) and progress of invasion in Sicily from 2015 until the present study assessment in 2019 (b)–(e). The map was created with the World Map Data from Natural Earth R package “rnaturalearth” version 0.1.0 available at https://CRAN.R-project.org/package=rnaturalearth.

We observed that the species was abundant at Capo Passero, Arenella and Plemmirio but absent or rare in the samples from Brucoli, Cannizzaro and Recanati (Fig. 1e). The collected samples for the exploratory phase were air-dried and transported to our laboratory in Bremen, Germany to quantitatively assess the progress of the invasion. The dried samples were weighed and split to obtain representative aliquots containing ~ 300 foraminifera per sample to determine the abundance of A. lobifera shells in the total assemblage (living + dead) and estimate its population density (individuals/g of sediment). The results are provided in Supplementary Table S1.

Following the exploratory survey, samples for genetic analyses were collected at Capo Passero and Plemmirio and at Avola located midway. The samples were collected as described above, but now living specimens of A. lobifera were isolated and individually transferred to micropaleontological slides, where they were air dried and stored at -20 °C (methods are detailed in Hallock et al.42; Schmidt et al.25; Stuhr et al.43). The isolated foraminifera were shipped on dry ice to our laboratory in Bremen, where they were stored at − 80 °C. The sampling in Israel was conducted in October 2019. We sampled at Shikmona, Haifa (Mediterranean Sea) and at Eilat, Gulf of Aqaba (Red Sea) where the presence of A. lobifera has already been documented25,44,45. The living specimens for genetic analyses were collected and isolated and transported as described above.

Between 10 and 24 specimens per site were selected for genetic analyses. We also analyzed nine specimens of A. lobifera collected at Okinawa, Japan (26.651819 N; 127.856243 E) in September 2015, that were present in our collection to serve as an outgroup. Each specimen was isolated in 50 µl of DOC buffer and the thick calcite shell was cracked with a sterilized crusher. Following the DOC protocol46, each specimen was incubated at 60 °C for 1 h followed by centrifugation at 10.000 rpm for 5 min and stored at 4 °C until further use.

Because we aimed at capturing a population dynamic process, we designed a protocol to access the most variable genomic region known for foraminifera, the Internal Transcribed Spacer (ITS)41 that is located between the ribosomal Small Sub-Unit (SSU) and the Large Sub-Unit (LSU). The genetic diversity of A. lobifera has been previously documented in the Indo-Pacific and the Mediterranean Sea based on a 600 bp fragment located at the end of the SSU (Prazeres et al. 2020). To make our results compatible with the existing data, we designed a protocol to amplify the same SSU fragment together with the entire ITS. We developed a semi-nested PCR protocol using the primers pairs S14F3 (5′-ACGCAMGTGTGAAACTTG-3′)—L5F (5′—TCGCCGTTACTAAGGRAATC—3′). and S14F1 (5′-AAGGGCACCACAAGAACGC-3′)—L5F46. The amplification was carried out using the GoTaq polymerase (Promega) with a PCR mix containing MilliQ water, 5× green buffer (final concentration: 1×), each primer (final concentration: 0.2 µmol/µl), MgCl2 (final concentration: 2.5 µmol/µl), dNTP mix (final concentration: 0.4 µmol/µl) and GoTaq polymerase (final concentration: 0.05 U/µl), and added DNA extract diluted 1:10 to reduce inhibition, within a final volume of 15 µl. The second PCR was carried out with the same mix but using 1 µl of the 1st PCR as the DNA template. The thermal cycling was as follows for both successive PCRs: initial denaturation at 95 °C for 2 min followed by 35 cycles of 30 s of denaturation at 95 °C, annealing for 30 s at 55 °C and extension at 72 °C for 45 s, followed by a final extension at 72 °C for 2 min.

The PCR product obtained was migrated on a 1.5% agarose gel and visually checked under UV light. The samples showing single bands were selected and purified using the QIAquick PCR purification kit (QIAGEN) following the manufacturer's instructions. Because of the presence of large intra-genomic variability among the multiple copies of the rRNA gene in foraminifera40, the purified PCR product was cloned using the TOPO® TA Cloning Kit (Invitrogen, USA). Amplicons were ligated to a pCR 2.1TOPO® vector, transformed into One Shot™ TOP10 chemically competent Escherichia coli cells, and grown overnight on LB-agar plates containing ampicillin (100 mg/ml). Eight to 16 clones per specimen were selected and placed in 1.5 ml tubes containing 30 µl of MilliQ water and a final PCR was performed and sent for Sanger sequencing with an external provider (LGC Genomics, Berlin). Due to the long fragment targeted (~ 1400 bp), each PCR product was sequenced from both ends using the primers S14F1 and L5F. The chromatograms were carefully checked and assembled, and the resulting sequences were deposited on NCBI under the accession numbers OP610171-OP610543. In addition to the new sequences, we also recovered all available A. lobifera sequences and associated metadata available on NCBI that were mostly produced by Prazeres et al.17 and Schmidt et al.25. We limited our query to sequences covering the entire ~ 600 bp fragment of the SSU, resulting in a total of 256 sequences. All newly generated sequences and publicly available sequences are provided with associated metadata in Supplementary Table S2.

To assess to what degree the invasive populations of A. lobifera differ from the source Indo-Pacific populations, we constructed phylogenetic trees and compared the distribution of patristic distances among sequences in the different genetic types of the species as identified by Prazeres et al.17. The patristic distances were analyzed for intra-genomic variability (genetic distances among sequences within the same specimen), for population-level variability (genetic distances among sequences from different specimens occurring at the same population), for regional variability (genetic distances among sequences from specimens occurring at different populations within the same oceanic basin) and for geographical variability (genetic distances among sequences from specimens occurring in different oceanic basins). To this end, we constructed two phylogenetic inferences, one including all the sequences of the dataset but covering only the SSU fragment, and one including all the sequences generated in this study and covering the SSU and ITS fragment. For each inference, the sequences were automatically aligned with MAFFT47 and a phylogenetic tree was inferred using RAxML-NG48 with 100 non-parametric bootstraps and using the substitution model TVM + I + G for the SSU alignment and the TIM2ef + I + G4 for the SSU + ITS alignment that were selected with Modeltest-NG49,50. Both trees are provided in Supplementary Fig. S1. After the inference, the patristic distances for the SSU tree and for the SSU + ITS tree were calculated and grouped according to the four categories of comparisons.

To identify the factors affecting intra-genomic distances within the invasive genotype, we used beta dispersion analysis based on Principle Coordinate analysis of square-root transformed patristic distances. We calculated the distance-to-centroid for each specimen for further statistical analyses (ANOVA and Wilcoxon test) to quantify the importance of the factors “oceanic basin” and “status of invasion” in structuring genetic diversity. The patristic distances and the statistical analyses were calculated using the packages ape51 and vegan52, respectively, from the software R 4.1.153.

To investigate the phylogeographic relationships among the different populations, haplotype networks were constructed for both the SSU alignment and the newly assembled SSU + ITS alignment. We constructed Median-Joining Networks (MJN)54, following an algorithm analogous to that proposed by Excoffier and Smouse55 that first constructs Minimum Spanning Trees (MSTs) from a matrix of pairwise distances (absolute number of differences) among haplotypes and includes all possible MSTs using the parsimony criterion to infer and add missing node haplotypes to the MJN graph. We defined ε = 0 for a more stringent distance criterion to select the most parsimonious pathway56. To allow the comparison of the population structure within the different genotypes, an analysis of molecular variance (AMOVA) between and within oceanic basins and populations was conducted for each genotype.

In addition to the AMOVA, we also calculated the phi-statistics (Fst) that refers to relative contributions of between-population variations to the overall genetic variation in the whole dataset. The groups were tested for adherence to neutrality (random evolution) assumptions, with Tajima's D57. Negative values of Tajima's D indicate an excess of low-frequency polymorphisms, consistent with positive directional selection or recent population expansion, whereas positive values indicate an excess of intermediate frequency polymorphism potentially due to balancing selection or population contraction58. Nucleotide diversity, number of haplotypes, and number of segregating sites were also calculated to investigate the degree of polymorphism within the genotypes. The networks and AMOVA were performed in the software PopART 1.756, and the genetic diversity indexes and Tajima’s D were calculated in the packages haplotypes59, ape51, and pegas60, from the software R 4.1.153.

Our field sampling in 2019 revealed that the invasion front of A. lobifera along the eastern coast of Sicily reached at least Plemmirio, but not beyond Brucoli. At the same time, we observed that the abundance of Amphistegina in Sicily increased dramatically compared to the observations by Guastella et al.27. For instance, in Capo Passero the species represented 35% of the foraminifera in 2017 and in our sampling in 2019 the relative abundance has risen to 75%. In the other two sites where we found A. lobifera (Arenella and Plemmirio), the species represented 17 and 49% of the foraminifera (Fig. 1, Supplementary Table S1). Thus, the invasive population has now become a major component of the assemblages on the southern coast of Sicily, but the invasion seems to have halted at the 14 °C winter isotherm along the northern Sicilian coast, which is considered to represent the thermal limit for the species12. This means that the Sicily populations sampled for genetic analyses (Fig. 1) represent not only the invasion front but also an established invasive population in the South of the sampled coastal transect.

For the analysis of genetic variability of the invasive and source populations, we successfully amplified the SSU + ITS of four to 11 specimens per locality and sequenced between two to 20 clones per specimen, resulting in 373 SSU + ITS sequences from 37 specimens. For the SSU analysis, we combined the 259 sequences acquired from NCBI17 with the new data, resulting in a dataset with 632 SSU sequences from 88 specimens (Supplementary Table S2). The phylogenetic tree inference conducted on the SSU fragment showed that all the newly sampled specimens represent lineage Ia as defined by Prazeres et al.17, confirming that the invasive population is sourced exclusively from the Red Sea, where only that genotype occurs.

In the median-joining haplotype network, 247 haplotypes and 387 segregating (polymorphic) sites were observed across the 632 sequences in the SSU alignment (Fig. 2). The structure of the network was consistent with the phylogenetic tree with the four lineages being clearly separated. No further structure was observed within genotype Ia, where the invasive Mediterranean populations, the native population in the Red Sea and the Western-Indian Ocean populations share the same common haplotypes and there is no evidence for a systematic (geographical) divergence among them.

Median-joining network based on SSU rRNA gene sequences of Amphistegina lobifera populations, including the invasive genotype Ia and the invasion front in Sicily.

The AMOVA of the SSU alignment revealed that within all genotypes of A. lobifera, the largest part of the overall genetic variation was explained by variation within populations (specimens collected from the same locality) compared to variation among oceanic basins or among populations, including within the invasive genotype Ia. The genotype Ia showed the lowest nucleotide diversity, indicating a reduced degree of polymorphism and mutation (Table 1) and it was the only genotype that showed significant deviation from neutrality (i.e., evolving randomly) based on Tajima’s D (Tajima’s D = − 2.92, p = 0.003). The negative Tajima’s D value indicates fewer haplotypes than segregating sites, which is a sign of an excess of rare alleles and can be considered an indicator of population expansion after a recent bottleneck or recent selective sweep57,61,62.

To assess whether this pattern is due to genetic processes in the invasive populations of the genotype Ia, we carried out the analysis separately for sequences from the Mediterranean and from the Red Sea, and the Indian Ocean (Table 2). This revealed that the deviation from neutrality is not present in the entire genotype, but only in its Red Sea and Mediterranean populations. Finally, we carried out the analysis separately for sequences representing established invaders (Eastern Mediterranean) and the invasion front (Sicily). This revealed that both show deviation from neutrality, and it seems accentuated in the established invaders (Table 2).

As expected, given the known faster mutation rate in the ITS fragment41, the SSU + ITS haplotype network showed much higher polymorphism, with 130 haplotypes and 895 segregating sites across 373 sequences (Fig. 3). The network revealed two main haplogroups corresponding to the invasive genotypes Ia and the outgroup genotype IIc from Japan, with neither haplogroup possessing a shared central haplotype.

Median-joining network based on SSU + ITS regions of rRNA gene sequences of Amphistegina lobifera populations from invasive genotype Ia and out-group (genotype IIc).

Like for the SSU alignment, the AMOVA comparing the different oceanic basins within the invasive genotype Ia (Table 1) revealed that the genetic variation is higher within populations (86.5%) than among oceanic basins (− 12.5%) or among populations (26.0%). The negative value in the AMOVA is an artifact of the statistical approach and should be interpreted as 0%, which means that the genetic variation between the oceanic basins does not contribute to the overall total variation at all. This is also shown by the low values of Fst (0.135) and nucleotide diversity (0.019) within genotype Ia (Table 1). Like for the SSU, the Tajima’s D revealed a significant departure from neutrality (Table 1) and indicated population size expansion (e.g., after a bottleneck or a selective sweep) and this pattern is observed in both Red Sea and Mediterranean populations (Table 2). Within the Mediterranean population the deviation from neutrality is observed in both established and new invaders.

The lack of geographical structure in the haplotype networks is reflected by the distribution of patristic distances among sequences (Fig. 4). The patristic distances calculated to determine the amount of variability at the intra-genomic level (different clones from same specimen), at the population level (different specimens from same population), at the regional level (different populations from same oceanic basin) and at the geographical level (different oceanic basins) revealed little evidence for an increase in genetic divergence with geographical distance within the invasive genotype Ia. The pattern appeared both in the SSU and SSU + ITS (Fig. 4a) analyses, although in the latter the distances were lower.

Patristic distances calculated based on phylogenetic trees for the SSU and SSU + ITS rRNA gene sequences of Amphistegina lobifera including all genotypes (a) and beta dispersion analysis of intra-genomic distances within invasive genotype Ia: by oceanic basin (b, top panel) and by status of invasion in the Mediterranean Sea (b, bottom panel). Stars represent levels of significance in the Wilcoxon test shown in Table 3.

The most striking pattern revealed by the patristic distances for both the SSU and the SSU + ITS data is the large reduction of intra-genomic variability in the invasive genotype Ia (Fig. 4a). This observation is highlighted when the invasive genotype Ia is compared to the other genotypes, none of which shows such a large reduction in the intra-genomic variability, except for the genotype IIIa that has too few sequences. To characterize the nature of this apparent reduction in intra-genomic variability, we compared the intra-genomic distance to centroid in the beta dispersion analysis in different populations of the invasive genotype Ia (Fig. 4b and Table 3).

This analysis reveals that the large reduction in the intra-genomic distance is present in populations from the Red Sea and significantly accentuated in the Mediterranean, both for the SSU and SSU + ITS (Fig. 4b). The significant reduction in the intra-genomic distances in the Mediterranean is observed both among the established invaders and in populations at the invasion front in Sicily (Fig. 4b). The degree of reduction is similar for all invasive populations with the small differences indicated in the SSU dataset not confirmed by the more informative SSU + ITS analysis (Fig. 4b).

The observed lack of genetic difference between the invasive Mediterranean populations and the source population in the Red Sea (Figs. 2, 3) confirms the postulated large dispersal potential of A. lobifera17. It also indicates that the invasion must have involved many of the genotypes present in the Red Sea population, rather than a hypothetical pre-adapted subtype. This means that any adaptation facilitating the invasion success was already present in the source population and the opening of the Suez Canal represented an artificial removal of an obstacle for a population that would have otherwise been able to expand beyond the Red Sea. Simultaneously, there does not seem to be any evidence for genetic differentiation of the invasive populations. Theoretically, the time since the invasion may have been too short for unique mutations to accumulate, but this is unlikely considering the observed high within-population variability that demonstrates that the ITS is a suitable marker to capture such a process. Instead, it is more likely that genetic differentiation in the Mediterranean is counteracted by continuous re-seeding by populations from the Red Sea.

At the same time, we observe that the invasive and Red Sea populations show reduced genetic diversity compared to the Indian Ocean populations of the genotype Ia (Table 2). This pattern is consistent with a genetic bottleneck, which is an expected phenomenon associated with invasion (e.g.,63,64,65). Alternatively, the lower sequence diversity could also be a consequence of a selective sweep61,62. In this scenario, the reduced diversity in the observed marker could signal strong positive selection against another allele located in the proximity, indicating that the invasion could be associated with the presence of some particularly favorable traits. However, the pattern of low nucleotide diversity and high polymorphism is observed already in the source population from the Red Sea. This would imply that the observed bottleneck or selective sweep already affected the Red Sea population and was not associated with the Lessepsian invasion. Indeed, during each glacial sea level lowstand, the Red Sea becomes hypersaline and inhospitable to most marine organisms, and the basin is repopulated during each deglaciation from the Indian Ocean66, with the last such event dating to about 11,000 years ago (e.g.,67). The observed reduced genetic diversity could be the result of this last population expansion and seems consistent with a bridgehead effect68, where the source population for the invasion into the Mediterranean Sea was in fact a result of a previous expansion from the Indian Ocean into the Red Sea. This effect could have preconditioned the Red Sea population towards a greater invasiveness, thus facilitating the subsequent colonization of the Mediterranean, but further reducing the genetic variability of the invasive population.

This reduced genetic variability observed in the Mediterranean Sea is consistent with a recent invasion. Previous studies suggested that Amphistegina may have invaded the Mediterranean Sea on earlier occasions (e.g., pre-historic occurrence during the Pleistocene69,70), or even that A. lobifera is native to the Mediterranean Sea71. The genetic structure of all Mediterranean populations analyzed to date indicate that if Amphistegina colonized the Mediterranean Sea in a more distnat past, it must have failed to establish a permanent population, until the most recent invasion event. If the current A. lobifera populations occurring in the Mediterranean Sea would be descending from ancient native or earlier Pleistocene invasive populations, then we should have observed distinct genetic differences compared to the Red Sea population, which was not the case.

Surprisingly, rather than a signal of genetic differentiation within the invading populations, we observe a strong and significant reduction of gene-copy variability throughout the Mediterranean populations (Fig. 4; Table 3), which was clearly stronger than in the source population from the Red Sea. This suppression of variability between copy variants in the invasive populations is associated with the retention of high genetic variation among specimens of the same population, requiring an explanation that reduces variability within a genome but not among individuals. One might argue that this low genetic diversity between copy variants observed in Mediterranean populations could be the result of gene surfing, i.e., random loss of genetic diversity over space during range expansion72. However, the loss is likely not random, because it affected all the different types invading the Mediterranean. If it would occur in only one population, then we should see reduced intra-genomic variability and reduced variability among individuals as well. Therefore, the most plausible alternative that explains both, the reduced intra-genomic variability and a high genetic variation among specimens of the same population, would be a change in the reproductive strategy towards the suppression of meiotic recombination. This is because sexual reproduction with recombination would be expected to promote genetic variation both among and within individuals. Hence, the observed reduction in intra-genomic variability could be a consequence of a shift towards an asexual reproductive mode that favors gene conversion-like processes leading to a loss of heterozygosity (e.g., as observed in Daphnia by Tucker et al.73 and in Trypanosoma by Weir et al.74).

Such gene conversion processes during asexual reproduction could lead to homogenization within a genome. In the presence of continuous re-seeding of the invasive populations from the Red Sea, the large variability among individuals would be preserved. Therefore, we speculate that the genetic structure of the invasive population reflects two processes: high dispersal potential of a pre-adapted Red Sea source population and a suppression of sexual reproduction. It is important to highlight, however, that the observed homogenization could also be a result of concerted evolution, where multiple copies of a gene within a genome evolve in a coordinated manner75,76. However, such process should not be acting within the Mediterranean population only, and it is therefore less likely that this process would be responsible for the particular profile we observed.

Like many other foraminifera Amphistegina is known to have a trimorphic life cycle involving a sexual generation (agamont) and two asexual generations (gamont and schizont), with no necessity for strict alternation between asexual and sexual generations (e.g.,77,78,79,80,81,82,83). In symbiotic organisms, a change in reproduction strategy can also be linked to the process of obtaining and maintaining symbiosis. Like in many other symbiont-bearing organisms, in LBFs, asexual reproduction is associated with a vertical transfer of symbionts to the offspring. Offspring generated by multiple fission has a large size (~ 40–50 µm) and can receive the cytoplasm of the parent together with the symbionts78. On the other hand, the tiny gametes (2–3 µm) cannot carry the symbionts and therefore the zygotes must acquire them from the environment. As a result, in cyclic schizogony (i.e., no alternation with sexual reproduction), the symbiont culture is maintained without the need to receive new symbionts from the environment78, although such uptake of symbionts from the environment (horizontal transmission) remains a possibility. This could provide an explanation for the lack of sexual reproduction in the invasive A. lobifera. In this scenario, the ability to reproduce by cyclic schizogony would represent an advantage, or even a prerequisite, for populating regions with environmental conditions that do not allow the acquisition of new symbionts77. However, it is unclear whether enhanced cyclic schizogony in the Mediterranean individuals is an actual adaptation or a plastic response from A. lobifera facing extreme conditions. Therefore, it is equally possible that the invaders do not engage in sexual reproduction (gamete release), or they do, and their attempts are not successful.

Interestingly, the slightly but significantly reduced intra-genomic variability in the Red Sea population indicates that the same process, but to a lesser degree, may act already on the northern Red Sea populations of the species. Compared to the Indian Ocean range of the species, the northern Red Sea already represents thermal conditions close to the cold limit of the species range. However, the Red-Sea populations possess tolerance to high temperatures25. Like in northern Red Sea corals8, the retention of high tolerance in A. lobifera is likely the result of thermal filtering of Indian Ocean populations entering the Red Sea from the south. Therefore, it is possible that already the northern Red Sea populations of A. lobifera live under stressful environmental conditions that lead to partial suppression of sexual reproduction.

There are numerous observations of shifts in reproductive strategy in association with biological invasions and several well-known invasive species are asexual, with most examples known among plants (e.g.,84,85,86,87). For instance, the reed canary grass Phalaris arundinacea (Poaceae) rapid invasion in wetlands in North America has been attributed to clonal subsidy, morphological plasticity, and nutrient availability88. And the invasiveness of the aquatic plant Myriophyllum aquaticum (Haloragaceae) is regarded to be facilitated by clonal integration, which has been found to significantly improve the growth and photosynthetic performance of daughter ramets89.

That marine invasion may impede sexual reproduction among protists has, however, not been documented before. This is significant because sex has an obvious long-term advantage by creating new combinations of genes that allow adaptation to future changing conditions90. In contrast, asexual reproduction can be an effective strategy to rapidly increase population size during the colonization of new areas80 but has the disadvantage of decreasing the ability of the invasive population to react to future change by adaptation. Therefore, if the invasive Amphistegina is unable to reproduce sexually in the Mediterranean, its future proliferation may be affected by the loss of adaptive potential in the face of the projected continued environmental changes in their newly conquered space91. Finally, the hypothesis of the loss of sexual reproduction could be tested on other invasive foraminiferal species in the Mediterranean. The outcomes of this analysis may enhance our understanding of adaptive behavioral patterns within this group under diverse stressful conditions and contribute to predicting the species that are more likely to be successful invaders.

Our results revealed that the invasion of the symbiont-bearing foraminifera A. lobifera in the Mediterranean Sea is facilitated by the combination of preadaptation and a high dispersal ability, with sustained re-seeding of the Mediterranean from the Red Sea. The invasion involves many of the genotypes present in the Red Sea population, rather than a specific subtype, indicating that the preadaptation to invasion was widespread in the source population. At the same time, the invasive populations show reduced intragenomic variability associated with sustained high genetic variation among specimens, which can be explained by a lower average heterozygosity due to increased gene conversion during asexual reproduction. The invasion therefore appears to be associated with a sustained change in reproductive strategy towards the abandonment of sex. Either the sexual reproduction is not triggered or cannot be completed due to adverse environmental conditions in the new habitat or, alternatively, because the zygotes have difficulty in acquiring symbionts from the environment. Either way, this discovery provides a new perspective on the cost of invasion in marine protists. If the invasion is facilitated by or requires a shift towards cyclic schizogony, the short-term gain of invasion into new habitats may be offset by a long-term loss of adaptive potential.

The dataset with sequences and associated metadata generated during the current study are available in Supplementary Table S2 and in the NCBI repository, under the accession numbers OP610171-OP610543.

Simberloff, D. et al. Impacts of biological invasions: What’s what and the way forward. Trends Ecol. Evol. 28, 58–66 (2013).

Article PubMed Google Scholar

Bellard, C., Cassey, P. & Blackburn, T. M. Alien species as a driver of recent extinctions. Biol. Lett. https://doi.org/10.1098/rsbl.2015.0623 (2016).

Article PubMed PubMed Central Google Scholar

Parker, I. M., Rodriguez, J. & Loik, M. E. An evolutionary approach to understanding the biology of invasions: Local adaptation and general-purpose genotypes in the weed Verbascum thapsus. Conserv. Biol. 17, 59–72 (2003).

Article Google Scholar

Weinig, C. Plasticity versus canalization: Population differences in the timing of shade-avoidance responses. Author(s): Cynthia Weinig. Society for the Study of Evolution Stable. https://www.jstor.org/stable/2640851 Society for the Study of Ev. Evolution (N Y) 54, 441–451 (2000).

Davis, M. B. & Shaw, R. G. Range shifts and adaptive responses to quaternary climate change. Science 1979(292), 673–679 (2001).

Article ADS Google Scholar

de Lafontaine, G., Napier, J. D., Petit, R. J. & Hu, F. S. Invoking adaptation to decipher the genetic legacy of past climate change. Ecology 99, 1530–1546 (2018).

Article PubMed Google Scholar

Mimura, M. & Aitken, S. N. Adaptive gradients and isolation-by-distance with postglacial migration in Picea sitchensis. Heredity (Edinb) 99, 224–232 (2007).

Article CAS PubMed Google Scholar

Fine, M., Gildor, H. & Genin, A. A coral reef refuge in the Red Sea. Glob Chang Biol 19, 3640–3647 (2013).

Article ADS PubMed Google Scholar

Rilov, G. & Galil, B. Marine Bioinvasions in the Mediterranean Sea—History, Distribution and Ecology 549–575 (Springer, Berlin, 2009). https://doi.org/10.1007/978-3-540-79236-9_31.

Book Google Scholar

Zenetos, Α et al. Alien species in the Mediterranean Sea by 2012. A contribution to the application of European Union’s Marine Strategy Framework Directive (MSFD). Part 2. Introduction trends and pathways. Mediterr. Mar. Sci. 13, 328–352 (2012).

Article Google Scholar

Shaltout, M. & Omstedt, A. Recent sea surface temperature trends and future scenarios for the Mediterranean Sea. Oceanologia 56, 411–443 (2014).

Article Google Scholar

Langer, M. R. & Hottinger, L. Biogeography of selected ‘larger’ foraminifera. Micropaleontology 46, 105–126 (2000).

Google Scholar

Alve, E. Colonization of new habitats by benthic foraminifera: A review. Earth-Sci. Rev. 46, 167–185 (1999).

Article ADS Google Scholar

Alve, E. & Goldstein, S. T. Propagule transport as a key method of dispersal in benthic foraminifera (Protista). Limnol. Oceanogr. 48, 2163–2170 (2003).

Article ADS Google Scholar

Hallock, P. & Reymond, C. E. Contributions of trimorphic life cycles to dispersal and evolutionary trends in large benthic foraminifers. J. Earth Sci. 33, 1425–1433 (2022).

Article CAS Google Scholar

Guy-Haim, T., Hyams-Kaphzan, O., Yeruham, E., Almogi-Labin, A. & Carlton, J. T. A novel marine bioinvasion vector: Ichthyochory, live passage through fish. Limnol. Oceanogr. Lett. 2, 81–90 (2017).

Article Google Scholar

Prazeres, M. et al. High dispersal capacity and biogeographic breaks shape the genetic diversity of a globally distributed reef-dwelling calcifier. Ecol. Evol. 10, 5976–5989 (2020).

Article PubMed PubMed Central Google Scholar

Langer, M. R., Weinman, A. E., Lotters, S. & Rodder, D. ‘“Strangers”’ in paradise: Modeling the biogeographic range. J. Foraminifer. Res. 42, 234–244 (2012).

Article Google Scholar

Mouanga, G. H. & Langer, M. R. At the front of expanding ranges: Shifting community structures at amphisteginid species range margins in the mediterranean sea. Neues Jahrb. Geol. Palaontol. Abh. 271, 141–150 (2014).

Article Google Scholar

Caruso, A. & Cosentino, C. The first colonization of the genus Amphistegina and other exotic benthic foraminifera of the Pelagian Islands and south-eastern Sicily (central Mediterranean Sea). Mar. Micropaleontol. 111, 38–52 (2014).

Article ADS Google Scholar

El Kateb, A., Stalder, C., Stainbank, S., Fentimen, R. & Spezzaferri, S. The genus Amphistegina (benthic foraminifera): Distribution along the southern Tunisian coast. Bioinvasions. Rec. 7, 391–398 (2018).

Article Google Scholar

Guastella, R. et al. Reconstructing bioinvasion dynamics through micropaleontologic analysis highlights the role of temperature change as a driver of alien foraminifera invasion. Front. Mar. Sci. 8, 1–15 (2021).

Article ADS Google Scholar

Langer, M. R. & Mouanga, G. H. Invasion of amphisteginid foraminifera in the Adriatic Sea. Biol. Invasions 18, 1335–1349 (2016).

Article Google Scholar

Meriç, E., Avşar, N. & Yokeş, M. B. Some alien foraminifers along the Aegean and southwestern coasts of Turkey. Micropaleontology 54, 307–349 (2008).

Article Google Scholar

Schmidt, C., Morard, R., Prazeres, M., Barak, H. & Kucera, M. Retention of high thermal tolerance in the invasive foraminifera Amphistegina lobifera from the Eastern Mediterranean and the Gulf of Aqaba. Mar. Biol. https://doi.org/10.1007/s00227-016-2998-4 (2016).

Article Google Scholar

Yokes, M. B., Meric, E. & Avsar, N. On the presence of alien foraminifera Amphistegina lobifera Larsen on the coasts of the Maltese Islands. Aquat. Invasions 2, 439–441 (2007).

Article Google Scholar

Guastella, R. et al. “Hidden invaders” conquer the Sicily Channel and knock on the door of the Western Mediterranean sea. Estuar. Coast. Shelf Sci. 225, 106234 (2019).

Article Google Scholar

Dlugosch, K. M. & Parker, I. M. Founding events in species invasions: Genetic variation, adaptive evolution, and the role of multiple introductions. Mol. Ecol. 17, 431–449 (2008).

Article CAS PubMed Google Scholar

Lee, C. E. Evolutionary genetics of invasive species. Trends Ecol. Evol. 17, 386–391 (2002).

Article Google Scholar

Sakai, A. K. et al. The population biology of invasive species. Annu. Rev. Ecol. Syst. 32, 305–332 (2001).

Article Google Scholar

Leger, E. A. & Rice, K. J. Assessing the speed and predictability of local adaptation in invasive California poppies (Eschscholzia californica). J. Evol. Biol. 20, 1090–1103 (2007).

Article CAS PubMed Google Scholar

Maron, J. L., Vilà, M., Bommarco, R., Elmendorf, S. & Beardsley, P. Rapid evolution of an invasive plant. Ecol. Monogr. 74, 261–280 (2004).

Article Google Scholar

Ward, S. M., Gaskin, J. F. & Wilson, L. M. Ecological genetics of plant invasion: What do we know?. Invasive Plant Sci. Manag. 1, 98–109 (2008).

Article Google Scholar

Bernardi, G., Azzurro, E., Golani, D. & Miller, M. R. Genomic signatures of rapid adaptive evolution in the bluespotted cornetfish, a Mediterranean Lessepsian invader. Mol. Ecol. 25, 3384–3396 (2016).

Article CAS PubMed Google Scholar

Quinn, T. P., Unwin, M. J. & Kinnison, M. T. Evolution of temporal isolation in the wild: Genetic divergence in timing of migration and breeding by introduced chinook salmon populations. Evolution (N Y) 54, 1372–1385 (2000).

CAS Google Scholar

Yue, G. H., Li, J., Bai, Z., Wang, C. M. & Feng, F. Genetic diversity and population structure of the invasive alien red swamp crayfish. Biol. Invasions 12, 2697–2706 (2010).

Article Google Scholar

André, A. et al. SSU rDNA divergence in planktonic foraminifera: Molecular taxonomy and biogeographic implications. PLoS One 9, e104641 (2014).

Article ADS PubMed PubMed Central Google Scholar

Pawlowski, J. & Lecroq, B. Short rDNA barcodes for species identification in foraminifera. J. Eukaryot. Microbiol. 57, 197–205 (2010).

Article CAS PubMed Google Scholar

Pillet, L., Fontaine, D. & Pawlowski, J. Intra-genomic ribosomal RNA polymorphism and morphological variation in Elphidium macellum suggests inter-specific hybridization in foraminifera. PLoS One 7, e32373 (2012).

Article ADS CAS PubMed PubMed Central Google Scholar

Weber, A. A. T. & Pawlowski, J. Wide occurrence of SSU rDNA intragenomic polymorphism in foraminifera and its implications for molecular species identification. Protist 165, 645–661 (2014).

Article CAS PubMed Google Scholar

Tsuchiya, M., Tazume, M. & Kitazato, H. Molecular characterization of the non-costate morphotypes of buliminid foraminifers based on internal transcribed region of ribosomal DNA (ITS rDNA) sequence data. Mar. Micropaleontol. 69, 212–224 (2008).

Article ADS Google Scholar

Hallock, P., Williams, S. K., Toler, E. M., Fisher & Talge, H. K. Bleaching in reef-dwelling foraminifers: Implications for reef decline. In 10th International Coral Reef Symposium, Vol. 737, 729–737 (2006).

Stuhr, M. et al. Variable thermal stress tolerance of the reef-associated symbiont-bearing foraminifera Amphistegina linked to differences in symbiont type. Coral Reefs 37, 811–824 (2018).

Article ADS Google Scholar

Hottinger, L., Halicz, E. & Reiss, Z. Recent Foraminiferida from the Gulf of Aqaba, Red Sea. (Dela SAZU, 1993).

Titelboim, D. et al. Thermal tolerance and range expansion of invasive foraminifera under climate changes. Sci. Rep. 9, 1–5 (2019).

Article CAS Google Scholar

Weiner, A. K. M. et al. Methodology for single-cell genetic analysis of planktonic foraminifera for studies of protist diversity and evolution. Front. Mar. Sci. 3, 1–15 (2016).

Article MathSciNet Google Scholar

Katoh, K. & Standley, D. M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 30, 772–780 (2013).

Article CAS PubMed PubMed Central Google Scholar

Kozlov, A. M., Darriba, D., Flouri, T., Morel, B. & Stamatakis, A. RAxML-NG: A fast, scalable and user-friendly tool for maximum likelihood phylogenetic inference. Bioinformatics 35, 4453–4455 (2019).

Article CAS PubMed PubMed Central Google Scholar

Darriba, Di. et al. ModelTest-NG: A new and scalable tool for the selection of DNA and protein evolutionary models. Mol. Biol. Evol. 37, 291–294 (2020).

Article MathSciNet CAS PubMed Google Scholar

Flouri, T. et al. The phylogenetic likelihood library. Syst. Biol. 64, 356–362 (2015).

Article CAS PubMed Google Scholar

Paradis, E. & Schliep, K. Ape 5.0: An environment for modern phylogenetics and evolutionary analyses in R.. Bioinformatics 35, 526–528 (2019).

Article CAS PubMed Google Scholar

Oksanen, J. et al. vegan: Community Ecology Package_. R package version 2.6-2, https://CRAN.R-project.org/package=vegan. Preprint at https://cran.r-project.org/package=vegan (2022).

R Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org/. Preprint at https://www.r-project.org/ (2022).

Bandelt, H. J., Forster, P. & Röhl, A. Median-joining networks for inferring intraspecific phylogenies. Mol. Biol. Evol. 16, 37–48 (1999).

Article CAS PubMed Google Scholar

Excoffier, L. & Smouse, P. E. Using allele frequencies and geographic subdivision to reconstruct gene trees within a species: Molecular variance parsimony. Genetics 136, 343–359 (1994).

Article CAS PubMed PubMed Central Google Scholar

Leigh, J. W. & Bryant, D. POPART: Full-feature software for haplotype network construction. Methods Ecol. Evol. 6, 1110–1116 (2015).

Article Google Scholar

Tajima, F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 123, 585–595 (1989).

Article CAS PubMed PubMed Central Google Scholar

Wachowiak, W., Salmela, M. J., Ennos, R. A., Iason, G. & Cavers, S. High genetic diversity at the extreme range edge: Nucleotide variation at nuclear loci in Scots pine (Pinus sylvestris L.) in Scotland. Heredity 106, 775–787. https://doi.org/10.1038/hdy.2010.118 (2011).

Article CAS PubMed Google Scholar

Aktas, C. Haplotypes: Manipulating DNA Sequences and Estimating Unambiguous Haplotype Network with Statistical Parsimony. R package version 1.1.2. URL https://cran.r-project.org/package=haplotypes. Preprint at https://cran.r-project.org/package=haplotypes (2020).

Paradis, E. Pegas: An R package for population genetics with an integrated-modular approach. Bioinformatics 26, 419–420 (2010).

Article CAS PubMed Google Scholar

Nielsen, R. et al. Genomic scans for selective sweeps using SNP data. Genome Res 15, 1566–1575 (2005).

Article CAS PubMed PubMed Central Google Scholar

Smith, J. M. & Haigh, J. The hitch-hiking effect of a favourable gene. Genet. Res. (Camb.) 23, 23–35 (1974).

Article CAS Google Scholar

Meimberg, H. et al. Molecular evidence for an extreme genetic bottleneck during introduction of an invading grass to California. Biol. Invasions 8, 1355–1366 (2006).

Article Google Scholar

Xue, D. X. et al. Successful worldwide invasion of the veined rapa whelk, Rapana venosa, despite a dramatic genetic bottleneck. Biol. Invasions 20, 3297–3314 (2018).

Article Google Scholar

Zepeda-Paulo, F. et al. Signatures of genetic bottleneck and differentiation after the introduction of an exotic parasitoid for classical biological control. Biol. Invasions 18, 565–581 (2016).

Article Google Scholar

Fenton, M., Geiselhart, S., Rohling, E. J. & Hemleben, C. Aplanktonic zones in the Red Sea. Mar. Micropaleontol. 40, 277–294 (2000).

Article ADS Google Scholar

Hemleben, C. et al. Three hundred eighty thousand year long stable isotope and faunal records from the Red Sea: Influence of global sea level change on hydrography. Paleoceanography 11, 147–156 (1996).

Article ADS Google Scholar

Lombaert, E. et al. Bridgehead effect in the worldwide invasion of the biocontrol Harlequin ladybird. PLoS One 5, e9743 (2010).

Article ADS PubMed PubMed Central Google Scholar

Meriç, E. et al. Did Amphistegina lobifera Larsen reach the Mediterranean via the Suez Canal?. Quat. Int. 401, 91–98 (2016).

Article ADS Google Scholar

Stulpinaite, R., Hyams-Kaphzan, O. & Langer, M. R. Alien and cryptogenic Foraminifera in the Mediterranean Sea: A revision of taxa as part of the EU 2020 marine strategy framework directive. Mediterr. Mar. Sci. 21, 719–758 (2020).

Google Scholar

Albano, P. G. et al. Alleged Lessepsian foraminifera prove native and suggest Pleistocene range expansions into the Mediterranean Sea. Mar. Ecol. Prog. Ser. 700, 65–78 (2022).

Article ADS CAS Google Scholar

Hallatschek, O. & Nelson, D. R. Gene surfing in expanding populations. Theor. Popul. Biol. 73, 158–170 (2008).

Article PubMed MATH Google Scholar

Tucker, A. E., Ackerman, M. S., Eads, B. D., Xu, S. & Lynch, M. Population-genomic insights into the evolutionary origin and fate of obligately asexual Daphnia pulex. Proc. Natl. Acad. Sci. U. S. A. 110, 15740–15745 (2013).

Article ADS CAS PubMed PubMed Central Google Scholar

Weir, W. et al. Population genomics reveals the origin and asexual evolution of human infective trypanosomes. Elife 5, 1–14 (2016).

Article Google Scholar

Liao, D. Concerted evolution: Molecular mechanism and biological implications. Am. J. Hum. Genet. 64, 24–30 (1999).

Article CAS PubMed PubMed Central Google Scholar

Dover, G. Molecular drive: A cohesive mode of species evolution. Nature 299, 111–117 (1982).

Article ADS CAS PubMed Google Scholar

Beavington-Penney, S. J. & Racey, A. Ecology of extant nummulitids and other larger benthic foraminifera: Applications in palaeoenvironmental analysis. Earth Sci. Rev. 67, 219–265 (2004).

Article ADS Google Scholar

Dettmering, C., Rottger, R., Hohenegger, J. & Schmaljohann, R. The trimorphic life cycle in foraminifera: Observations from cultures allow new evaluation. Eur. J. Protistol. 34, 363–368 (1998).

Article Google Scholar

Hallock, P. & Seddighi, M. Why did some larger benthic foraminifera become so large and flat?. Sedimentology https://doi.org/10.1111/sed.12837 (2020).

Article Google Scholar

Harney, J. N., Hallock, P. & Talge, H. K. Observations on a trimorphic life cycle in Amphistegina gibbosa populations from the Florida keys. J Foraminifer. Res. 28, 141–147 (1998).

Article Google Scholar

Hollaus, S. S. & Hottinger, L. Temperature dependance of endosymbiontic relationships? Evidence from the depth range of mediterranean Amphistegina lessonii (Foraminiferida) truncated by the thermocline. Eclogae Geol. Helv. 90, 591–597 (1997).

Google Scholar

Hottinger, L. Larger foraminifera, giant cells with a historical background. Naturwissenschaften 69, 361–371 (1982).

Article ADS Google Scholar

Reiss, Z. & Hottinger, L. The Gulf of Aqaba (Springer, Berlin, 1984). https://doi.org/10.1007/978-3-642-69787-6.

Book Google Scholar

Liu, J. et al. Invasive alien plants in China: Role of clonality and geographical origin. Biol. Invasions 8, 1461–1470 (2006).

Article Google Scholar

Poulin, J., Weller, S. G. & Sakai, A. K. Genetic diversity does not affect the invasiveness of fountain grass (Pennisetum setaceum) in Arizona, California and Hawaii. Divers. Distrib. 11, 241–247 (2005).

Article Google Scholar

Raybould, A. F., Gray, A. J., Lawrence, M. J. & Marshall, D. F. The evolution of Spartina anglica C.E. Hubbard (Gramineae): Origin and genetic variability. Biol. J. Linn. Soc. 43, 111–126 (1991).

Article Google Scholar

Lyman, J. C. & Ellstrand, N. C. Clonal diversity in taraxacum officinale (compositae), an apomict. Heredity (Edinb.) 53, 1–10 (1984).

Article Google Scholar

Maurer, D. A. & Zedler, J. B. Differential invasion of a wetland grass explained by tests of nutrients and light availability on establishment and clonal growth. Oecologia 131, 279–288 (2002).

Article ADS PubMed Google Scholar

You, W., Yu, D., Liu, C., Xie, D. & Xiong, W. Clonal integration facilitates invasiveness of the alien aquatic plant Myriophyllum aquaticum L. under heterogeneous water availability. Hydrobiologia 718, 27–39 (2013).

Article Google Scholar

De Meeûs, T., Prugnolle, F. & Agnew, P. Asexual reproduction: Genetics and evolutionary aspects. Cell. Mol. Life Sci. 64, 1355–1372 (2007).

Article PubMed Google Scholar

Lejeusne, C., Chevaldonné, P., Pergent-Martini, C., Boudouresque, C. F. & Pérez, T. Climate change effects on a miniature ocean: The highly diverse, highly impacted Mediterranean Sea. Trends Ecol. Evol. 25, 250–260 (2010).

Article PubMed Google Scholar

Download references

We express our deepest gratitude to the four anonymous reviewers for their constructive feedback, which significantly improved the quality of this manuscript. We also express our gratitude to Dr. Giulia Visconti, director of the Plemmirio marine reserve. Sampling was possible thanks to the agreement with the Dipartimento di Scienze della Terra e del Mare (DiSTeM) of Palermo University. We thank Israel Nature and Parks Authority to have supported both sampling campaigns that were conducted in Israel under the permit 42055/2018, and Dr. Miguel Frada who kindly provided lab access at the IUI biological station. This study was funded by the BMBF-MOST cooperation in Marine Sciences Grant No. 03F0820A “ForaInva”, and by the Cluster of Excellence “The Ocean Floor—Earth’s Uncharted Interface” funded by the German Research Foundation (DFG).

Open Access funding enabled and organized by Projekt DEAL.

Center for Marine Environmental Sciences, MARUM, Universität Bremen, Bremen, Germany

Débora S. Raposo, Christiane Hassenrück, Michal Kucera & Raphaël Morard

Department of Biology and Biochemistry, University of Houston, Houston, USA

Rebecca A. Zufall

Dipartimento di Scienze e Tecnologie Biologiche Chimiche e Farmaceutiche, Università degli Studi di Palermo, Palermo, Italy

Antonio Caruso

Department of Earth Sciences, University of Oxford, Oxford, UK

Danna Titelboim

Department of Earth and Environmental Sciences, Ben Gurion University of the Negev, Beer Sheva, Israel

Sigal Abramovich

Department of Biological Oceanography, Leibniz Institute for Baltic Sea Research Warnemünde (IOW), Rostock, Warnemünde, Germany

Christiane Hassenrück

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

All authors contributed to the writing and discussion of the study. Material preparation and data collection were performed by D.S.R., D.T., S.A., and R.M. Data analysis was performed by D.S.R., C.H., M.K., and R.M. All authors read and approved the final manuscript.

Correspondence to Débora S. Raposo.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

S. Raposo, D., A. Zufall, R., Caruso, A. et al. Invasion success of a Lessepsian symbiont-bearing foraminifera linked to high dispersal ability, preadaptation and suppression of sexual reproduction. Sci Rep 13, 12578 (2023). https://doi.org/10.1038/s41598-023-39652-y

Download citation

Received: 04 October 2022

Accepted: 28 July 2023

Published: 03 August 2023

DOI: https://doi.org/10.1038/s41598-023-39652-y

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.