Intraspecific genetic lineages of a marine mussel show behavioural divergence and spatial segregation over a tropical/subtropical biogeographic transition

The effect of interspecific and intraspecific diversity on microplastic ingestion in two co-occurring mussel species in South Africa

Interspecific and intraspecific diversity are essential components of biodiversity with far-reaching implications for ecosystem function and service provision. Importantly, genotypic and phenotypic variation within a species can affect responses to anthropogenic pressures more than interspecific diversity. We investigated the effects of interspecific and intraspecific diversity on microplastic ingestion by two coexisting mussel species in South Africa, Mytilus galloprovincialis and Perna perna, the latter occurring as two genetic lineages. We found significantly higher microplastic abundance in M. galloprovincialis (0.54 ± 0.56 MP items g-1WW) than P. perna (0.16 ± 0.21 MP items g-1WW), but no difference between P. perna lineages. Microbeads were the predominant microplastic (76 % in P. perna, 99 % in M. galloprovincialis) and polyethylene the prevalent polymer. Interspecific differences in microplastic abundance varied across locations, suggesting diverse sources of contamination. We suggest that microplastic ingestion can be species-specific even in organisms that coexist and play similar functional roles within ecosystems.

Genomic phylogeography illuminates deep cyto-nuclear discordances in midwife toads (Alytes)

The advent of genomic methods allows us to revisit the evolutionary history of organismal groups for which robust phylogenies are still lacking, particularly in species complexes that frequently hybridize. In this study, we conduct RAD-sequencing (RAD-seq) analyses of midwife toads (genus Alytes), an iconic group of western Mediterranean amphibians famous for their parental care behavior, but equally infamous for the difficulties to reconstruct their evolutionary history. Through admixture and phylogenetic analyses of thousands of loci, we provide the most comprehensive phylogeographic framework for the A. obstetricans complex to date, as well as the first fully resolved phylogeny for the entire genus. As part of this effort, we carefully explore the influence of different sampling schemes and data filtering thresholds on tree reconstruction, showing that several, slightly different, yet robust topologies may be retrieved with small datasets obtained by stringent SNP calling parameters, especially when admixed individuals are included. In contrast, analyses of incomplete but larger datasets converged on the same phylogeny, irrespective of the reconstruction method used or the proportion of missing data. The Alytes tree features three Miocene-diverged clades corresponding to the proposed subgenera Ammoryctis (A. cisternasii), Baleaphryne (A. maurus, A. dickhilleni and A. muletensis), and Alytes (A. obstetricans complex). The latter consists of six evolutionary lineages, grouped into three clades of Pliocene origin, and currently delimited as two species: (1) A. almogavarii almogavarii and A. a. inigoi; (2) A. obstetricans obstetricans and A. o. pertinax; (3) A. o. boscai and an undescribed taxon (A. o. cf. boscai). These results contradict the mitochondrial tree, due to past mitochondrial captures in A. a. almogavarii (central Pyrenees) and A. o. boscai (central Iberia) by A. obstetricans ancestors during the Pleistocene. Patterns of admixture between subspecies appear far more extensive than previously assumed from microsatellites, causing nomenclatural uncertainties, and even underlying the reticulate evolution of one taxon (A. o. pertinax). All Ammoryctis and Baleaphryne species form shallow clades, so their taxonomy should remain stable. Amid the prevalence of cyto-nuclear discordance among terrestrial vertebrates and the usual lack of resolution of conventional nuclear markers, our study advocates for phylogeography based on next-generation sequencing, but also encourages properly exploring parameter space and sampling schemes when building and analyzing genomic datasets.

Transcriptome wide analyses reveal intraspecific diversity in thermal stress responses of a dominant habitat-forming species

The impact of climate change on biodiversity has stimulated the need to understand environmental stress responses, particularly for ecosystem engineers whose responses to climate affect large numbers of associated organisms. Distinct species differ substantially in their resilience to thermal stress but there are also within-species variations in thermal tolerance for which the molecular mechanisms underpinning such variation remain largely unclear. Intertidal mussels are well-known for their role as ecosystem engineers. First, we exposed two genetic lineages of the intertidal mussel Perna perna to heat stress treatments in air and water. Next, we ran a high throughput RNA sequencing experiment to identify differences in gene expression between the thermally resilient eastern lineage and the thermally sensitive western lineage. We highlight different thermal tolerances that concord with their distributional ranges. Critically, we also identified lineage-specific patterns of gene expression under heat stress and revealed intraspecific differences in the underlying transcriptional pathways in response to warmer temperatures that are potentially linked to the within-species differences in thermal tolerance. Beyond the species, we show how unravelling within-species variability in mechanistic responses to heat stress promotes a better understanding of global evolutionary trajectories of the species as a whole in response to changing climate.

The relative effects of interspecific and intraspecific diversity on microplastic trapping in coastal biogenic habitats

Our understanding of how anthropogenic stressors such as climate change and plastic pollution interact with biodiversity is being widened to include diversity below the species level, i.e., intraspecific variation. The emerging appreciation of the key ecological importance of intraspecific diversity and its potential loss in the Anthropocene, further highlights the need to assess the relative importance of intraspecific versus interspecific diversity. One such issue is whether a species responds as a homogenous whole to plastic pollution. Using manipulative field transplant experiments and laboratory-controlled hydrodynamic simulations, we assessed the relative effects of intraspecific and interspecific diversity on microplastic trapping in coastal biogenic habitats dominated by two key bioengineers, the brown intertidal macroalgae Fucus vesiculosus and F. guiryi. At the individual level, northern morphotypes of F. guiryi trapped more microplastics than southern individuals, and F. vesiculosus trapped more microplastics than F. guiryi. Canopy density varied among species, however, leading to reversed patterns of microplastic accumulation, with F. guiryi canopies accumulating more microplastics than those of F. vesiculosus, while no differences were observed between the canopies of F. guiryi morphotypes. We emphasize the importance of assessing the effects of intraspecific variation which, along with other crucial factors such as canopy density, flow velocity and polymer composition, modulates the extent of microplastic accumulation in coastal biogenic habitats. Our findings indicate that a realistic estimation of plastic accumulation in biogenic habitats requires an understanding of within- and between-species traits at both the individual and population levels.

Effects of habitat quality on abundance, size and growth of mussel recruits

Recruitment of mussels is a complex process with the successful arrival of individuals hinging on the availability of suitable habitats. We examined the effects of adult mussels as settlement habitat and the degree to which the suitability of habitat they offer is species-specific by comparing the recruitment success of intertidal mussels. We hypothesised that mussel recruitment and early growth are dictated by the quality of habitat offered by conspecifics adults. We used a unique experimental arena on the south coast of South Africa, where Mytilus galloprovincialis and two lineages of Perna perna co-exist. Treatments were based on the translocation of individuals of M. galloprovincialis, western- and eastern lineage of P. perna to a single site, where artificial beds were created and sampled monthly over one year. Recruit's number, their sizes and growth were greater within beds of the western lineage of Perna than eastern lineage or Mytilus beds. The results clearly demonstrate that the quality of settlement habitat offered by adult beds differs among adult lineages/species and affects rates of settlement and the early growth of recruits. This effect extends to the intraspecific level; we found greater differences in density and growth of recruits between lineages of Perna than between either lineage and M. galloprovincialis.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s10750-022-04994-7.

Ectoparasites reduce scope for growth in a rocky-shore mussel (Perna perna) by raising maintenance costs

Endolithic cyanobacteria are ubiquitous colonisers of organic and inorganic carbonate substrata that frequently attack the shells of mussels, eroding the shell to extract carbon, often with population infestation rates of >80%. This reduces host physiological condition and ultimately leads to shell collapse and mortality, compromising the services provided by these important ecosystem engineers. While the ecological implications of this and similar interactions have been examined, our understanding of the underlying mechanisms driving the physiological responses of infested hosts remains limited. Using field and laboratory experiments, we assessed the energetic costs of cyanobacterial infestation to the intertidal brown mussel (Perna perna). In the field we found that growth (measured as both increase in shell length and rate of biomineralization) and reproductive potential of clean mussels are greater than those of infested individuals. To explore the mechanisms behind these effects, we compared the energy allocation of parasite-free and infested mussels using the scope for growth (SFG) framework. This revealed a lower SFG in parasitized mussels attributed to an energetic imbalance caused by increased standard metabolic rates, without compensation through increased feeding or reduced excretion of ammonia. Separate laboratory assays showed no differences in calcium uptake rates, indicating that infested mussels do not compensate for shell erosion through increased mineralization. This suggests that the increased maintenance costs detected reflect repair of the organic component of the inner nacreous layer of the shell, an energetically more demanding process than mineralization. Thus, parasite-inflicted damage reduces SFG directly through the need for increased basal metabolic rate to drive shell repair without compensatory increases in energy intake. This study provides a first perspective of the physiological mechanisms underlying this parasite-host interaction, a critical step towards a comprehensive understanding of the ecological processes driving dynamics of this intertidal ecosystem engineer.

Intraspecific diversity in an ecological engineer functionally trumps interspecific diversity in shaping community structure

Can intraspecific diversity functionally supersede interspecific diversity? Recent studies have established the ecological effects of intraspecific variation on a number of ecosystem dynamics including resilience and productivity and we hypothesised that they may functionally exceed those of species diversity. We focused on a coastal ecosystem dominated by two coexisting bioengineering mussel species, one of which, Perna perna, displays two distinct phylogeographic lineages. A manipulative field experiment revealed greater habitat structural complexity and a more benign microscale environment within beds of the eastern lineage than those of the western lineage or the second species (Mytilus galloprovincialis); the latter two did not differ. Similarly, while infaunal species abundance and biomass differed significantly between the two lineages of Perna, there was no such difference between Mytilus and the western Perna lineage. The evenness and diversity of associated infaunal assemblages responded differently. Diversity differed relatively weakly between species, while evenness showed a very strong difference between conspecific lineages. Our results show that variation within a species can functionally supersede diversity between species. As the two P. perna lineages have different physiological tolerances, we expect them to react differently to environmental change. Our findings indicate that predicting the ecosystem-level consequences of climate change requires an understanding of the relative strengths of within- and between-species differences in functionality.

Rejection of the genetic implications of the "Abundant Centre Hypothesis" in marine mussels

The ‘Abundant-Centre Hypothesis’ is a well-established but controversial hypothesis stating that the abundance of a species is highest at the centre of its range and decreases towards the edges, where conditions are unfavourable. As genetic diversity depends on population size, edge populations are expected to show lower intra-population genetic diversity than core populations, while showing high inter-population genetic divergence. Here, the genetic implications of the Abundant-Centre Hypothesis were tested on two coastal mussels from South Africa that disperse by means of planktonic larvae, the native Perna perna and the invasive Mytilus galloprovincialis. Genetic structure was found within P. perna, which, together with evidence from Lagrangian particle simulations, points to significant reductions in gene flow between sites. Despite this, the expected diversity pattern between centre and edge populations was not found for either species. We conclude that the genetic predictions of the Abundant-Centre Hypothesis are unlikely to be met by high-dispersal species with large population sizes, and may only become evident in species with much lower levels of connectivity.

Biogeographical Patterns of Endolithic Infestation in an Invasive and an Indigenous Intertidal Marine Ecosystem Engineer

By altering the phenotypic properties of their hosts, endolithic parasites can modulate the engineering processes of marine ecosystem engineers. Here, we assessed the biogeographical patterns of species assemblages, prevalence and impact of endolithic parasitism in two mussel species that act as important ecosystem engineers in the southern African intertidal habitat, Perna perna and Mytilus galloprovincialis. We conducted large-scale surveys across three biogeographic regions along the South African coast: the subtropical east coast, dominated by the indigenous mussel, P. perna, the warm temperate south coast, where this species coexists with the invasive Mediterranean mussel, M. galloprovincialis, and the cool temperate west coast dominated by M. galloprovincialis. Infestation increased with mussel size, and in the case of M. galloprovincialis we found a significantly higher infestation in the cool temperate bioregion than the warm temperate region. For P. perna, the prevalence of infestation was higher on the warm temperate than the subtropical region, though the difference was marginally non-significant. On the south coast, there was no significant difference in infestation prevalence between species. Endolith-induced mortality rates through shell collapse mirrored the patterns for prevalence. For P. perna, endolith species assemblages revealed clear grouping by bioregions. Our findings indicate that biogeography affects cyanobacteria species composition, but differences between biogeographic regions in their effects are driven by environmental conditions.

Evidence for rangewide panmixia despite multiple barriers to dispersal in a marine mussel

Oceanographic features shape the distributional and genetic patterns of marine species by interrupting or promoting connections among populations. Although general patterns commonly arise, distributional ranges and genetic structure are species-specific and do not always comply with the expected trends. By applying a multimarker genetic approach combined with Lagrangian particle simulations (LPS) we tested the hypothesis that oceanographic features along northeastern Atlantic and Mediterranean shores influence dispersal potential and genetic structure of the intertidal mussel Perna perna. Additionally, by performing environmental niche modelling we assessed the potential and realized niche of P. perna along its entire native distributional range and the environmental factors that best explain its realized distribution. Perna perna showed evidence of panmixia across >4,000 km despite several oceanographic breaking points detected by LPS. This is probably the result of a combination of life history traits, continuous habitat availability and stepping-stone dynamics. Moreover, the niche modelling framework depicted minimum sea surface temperatures (SST) as the major factor shaping P. perna distributional range limits along its native areas. Forthcoming warming SST is expected to further change these limits and allow the species to expand its range polewards though this may be accompanied by retreat from warmer areas.

Hypoxia effects on oxidative stress and immunocompetence biomarkers in the mussel Perna perna (Mytilidae, Bivalvia)

This study investigated the effects of hypoxia on oxidative stress response and immune function in mussels Perna perna exposed to air for 6, 12, 24 and 48 h. In air-exposed mussels, the antioxidant enzymes superoxide dismutase (SOD), catalase, and glutathione reductase (GR) activities were lower in gill tissues (24-48 h) and digestive gland (12 h), while the glutathione peroxidase and GR activities were increased in the digestive gland (48 h). In both tissues, aerial exposure promoted a rapid (6 h) and persistent (up to 48 h) increase of glutathione levels. Decreased hemocyte count and viability, as well as increased phagocytic activity and cellular adhesion capacity were detected after prolonged aerial exposure (>12 h). In summary, induction of thiol pools, altered antioxidant enzyme activities, and activation of immune responses were detected in hypoxia exposed brown mussels, indicating hypoxia induced tissue-specific responses in both antioxidant and immune systems.

Oceanographic Conditions Limit the Spread of a Marine Invader along Southern African Shores

Invasive species can affect the function and structure of natural ecological communities, hence understanding and predicting their potential for spreading is a major ecological challenge. Once established in a new region, the spread of invasive species is largely controlled by their dispersal capacity, local environmental conditions and species interactions. The mussel Mytilus galloprovincialis is native to the Mediterranean and is the most successful marine invader in southern Africa. Its distribution there has expanded rapidly and extensively since the 1970s, however, over the last decade its spread has ceased. In this study, we coupled broad scale field surveys, Ecological Niche Modelling (ENM) and Lagrangian Particle Simulations (LPS) to assess the current invaded distribution of M. galloprovincialis in southern Africa and to evaluate what prevents further spread of this species. Results showed that all environmentally suitable habitats in southern Africa have been occupied by the species. This includes rocky shores between Rocky Point in Namibia and East London in South Africa (approx. 2800 km) and these limits coincide with the steep transitions between cool-temperate and subtropical-warmer climates, on both west and southeast African coasts. On the west coast, simulations of drifting larvae almost entirely followed the northward and offshore direction of the Benguela current, creating a clear dispersal barrier by advecting larvae away from the coast. On the southeast coast, nearshore currents give larvae the potential to move eastwards, against the prevalent Agulhas current and beyond the present distributional limit, however environmental conditions prevent the establishment of the species. The transition between the cooler and warmer water regimes is therefore the main factor limiting the northern spread on the southeast coast; however, biotic interactions with native fauna may also play an important role.