Wider sampling reveals a non-sister relationship for geographically contiguous lineages of a marine mussel

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.

The Complete Mitochondrial Genome of Mytilisepta virgata (Mollusca: Bivalvia), Novel Gene Rearrangements, and the Phylogenetic Relationships of Mytilidae

The circular mitochondrial genome of Mytilisepta virgata spans 14,713 bp, which contains 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22 transfer RNA genes. Analysis of the 13 PCGs reveals that the mitochondrial gene arrangement of Mytilisepta is relatively conserved at the genus level. The location of the atp8 gene in Mytilisepta keenae differs from that of other species. However, compared with the putative molluscan ancestral gene order, M. virgata exhibits a high level of rearrangement. We constructed phylogenetic trees based on concatenated 12 PCGs from Mytilidae. As a result, we found that M. virgata is in the same clade as other Mytilisepta spp. The result of estimated divergence times revealed that M. virgata and M. keenae diverged around the early Paleogene period, although the oldest Mytilisepta fossil was from the late or upper Eocene period. Our results provide robust statistical evidence for a sister-group relationship within Mytilida. The findings not only confirm previous results, but also provide valuable insights into the evolutionary history of Mytilidae.

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 first use of LC-MS/MS proteomic approach in the brown mussel Perna perna after bacterial challenge: Searching for key proteins on immune response

The brown mussel Perna perna is a valuable fishing resource, primarily in tropical and subtropical coastal regions. Because of their filter-feeding habits, mussels are directly exposed to bacteria in the water column. Escherichia coli (EC) and Salmonella enterica (SE) inhabit human guts and reach the marine environment through anthropogenic sources, such as sewage. Vibrio parahaemolyticus (VP) is indigenous to coastal ecosystems but can be harmful to shellfish. In this study, we aimed to assess the protein profile of the hepatopancreas of P. perna mussel challenged by introduced - E. coli and S. enterica - and indigenous marine bacteria - V. parahaemolyticus. Bacterial-challenge groups were compared with non-injected (NC) and injected control (IC) - that consisted in mussels not challenged and mussels injected with sterile PBS-NaCl, respectively. Through LC-MS/MS proteomic analysis, 3805 proteins were found in the hepatopancreas of P. perna. From the total, 597 were significantly different among conditions. Mussels injected with VP presented 343 proteins downregulated compared with all the other conditions, suggesting that VP suppresses their immune response. Particularly, 31 altered proteins - upregulated or downregulated - for one or more challenge groups (EC, SE, and VP) compared with controls (NC and IC) are discussed in detail in the paper. For the three tested bacteria, significantly different proteins were found to perform critical roles in immune response at all levels, namely: recognition and signal transduction; transcription; RNA processing; translation and protein processing; secretion; and humoral effectors. This is the first shotgun proteomic study in P. perna mussel, therefore providing an overview of the protein profile of the mussel hepatopancreas, focused on the immune response against bacteria. Hence, it is possible to understand the immune-bacteria relationship at molecular levels better. This knowledge can support the development of strategies and tools to be applied to coastal marine resource management and contribute to the sustainability of coastal systems.

Comparative mitogenomic analyses and gene rearrangements reject the alleged polyphyly of a bivalve genus

Background

The order and orientation of genes encoded by animal mitogenomes are typically conserved, although there is increasing evidence of multiple rearrangements among mollusks. The mitogenome from a Brazilian brown mussel (hereafter named B1) classified as Perna perna Linnaeus, 1758 and assembled from Illumina short-length reads revealed an unusual gene order very different from other congeneric species. Previous mitogenomic analyses based on the Brazilian specimen and other Mytilidae suggested the polyphyly of the genus Perna.

Methods

To confirm the proposed gene rearrangements, we sequenced a second Brazilian P. perna specimen using the "primer-walking" method and performed the assembly using as reference Perna canaliculus. This time-consuming sequencing method is highly effective when assessing gene order because it relies on sequentially-determined, overlapping fragments. We also sequenced the mitogenomes of eastern and southwestern South African P. perna lineages to analyze the existence of putative intraspecific gene order changes as the two lineages show overlapping distributions but do not exhibit a sister relationship.

Results

The three P. perna mitogenomes sequenced in this study exhibit the same gene order as the reference. CREx, a software that heuristically determines rearrangement scenarios, identified numerous gene order changes between B1 and our P. perna mitogenomes, rejecting the previously proposed gene order for the species. Our results validate the monophyly of the genus Perna and indicate a misidentification of B1.

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.

How does the brown mussel Perna perna respond to environmental pollution? A review on pollution biomarkers

The brown mussel Perna perna (Linnaeus, 1758) is a valuable resource for aquaculture in tropical and subtropical coastal regions. It presents desirable characteristics for biomonitoring, including being sessile, widely distributed and abundant, and is a filter-feeder able to accumulate several classes of pollutants (e.g., metals, hydrocarbons, among others). Mussels' biological responses to pollution exposure can be measured as biomarkers, which include alterations ranging from molecular to physiological levels, to estimate the degree of environmental contamination and its effects on biota. This full review compiles two decades (2000-2020) of literature concerning biological effects on P. perna mussel caused by environmental pollutants (i.e., metals, hydrocarbons, and emerging pollutants), considering environmental and farm-based biomonitoring. Biochemical markers related to mussels' oxidative status were efficient for the biomonitoring of metals (i.e., antioxidant enzymes associated with oxidative damage in biomolecules). Genotoxicity and cytotoxicity indicators (i.e., comet, micronucleus, and neutral red assays) provided a depiction of hydrocarbon contamination. The neutral red assay gave a time-concentration cytotoxic response to a wide range of pollutants, including emerging pollutants (e.g., pharmaceuticals and biocides) and hydrocarbons. Perna perna hemocyte parameters provided a useful approach for biocide biomonitoring. This paper summarizes useful biomarkers from molecular to physiological levels in this mussel species used to identify and quantify the degree of coastal pollution. An integrated biomarker analysis may provide a way to overcome possible biomarker variations and assess multi-polluted sites. Nevertheless, it is necessary to investigate biomarker variations according to natural factors (e.g., season and gonad maturation stage) to standardize them for trustworthy biomonitoring.

Coastal dunefields maintain pre-Holocene genetic structure in a rocky shore red alga

Most intertidal algae have limited dispersal potential, and areas that lack hard substratum suitable for attachment are thus expected to isolate regional populations from each other. Here, we used nuclear and mitochondrial genetic data to compare genetic structure in two co-distributed intertidal red algae with different dispersal potential along the South African coastline. Gelidium pristoides is divided into a south-eastern and a south-western evolutionary lineage separated by extensive, continuous sandy shoreline habitat adjacent to coastal dunefields. In contrast, Hypnea spicifera is genetically homogeneous throughout its range. In G. pristoides, the genetic breaks are associated with contemporary coastal dunefields. The age of the divergence event suggests that this may reflect the effect of older dispersal barriers, and that genetic structure was subsequently maintained by the formation of contemporary coastal dunefields.

Mussels get higher: A study on the occurrence of cocaine and benzoylecgonine in seawater, sediment and mussels from a subtropical ecosystem (Santos Bay, Brazil)

Data on the occurrence of cocaine (COC) and benzoylecgonine (BE) in marine environmental compartments are still limited, with few studies reporting superficial water contamination, mainly in tropical zones. In this sense, environmental data of these substances are essential to identify potential polluting sources, as well as their impact in costal ecosystems. The aim of this study was to evaluate the occurrence of COC and BE in seawater, sediment and mussels from a subtropical coastal zone (Santos Bay, São Paulo, Brazil), as well as to determine a field measured Bioaccumulation Factor (BAF). COC and BE were detected in all water samples in concentrations ranging from 1.91 ng·L^-1 to 12.52 ng·L^-1 and 9.88 ng·L^-1 to 28.53 ng·L^-1, respectively. In sediments, only COC was quantified in concentrations ranging from 0.94 ng·g^-1 to 46.85 ng·g^-1. Similarly, only COC was detected in tissues of mussels 0.914 μg·kg^-1 to 4.58 μg·kg^-1 (ww). The field-measured BAF ranged from 163 to 1454 (L·kg^-1). Our results pointed out a widespread contamination by cocaine and its main human metabolite benzoylecgonine in Santos Bay. Mussels were able to accumulate COC in areas used by residents and tourists for bathing, fishing, and harvest, denoting concern to human health. Therefore, our data can be considered a preliminary assessment, which indicates the need to evaluate drugs (including illicit as COC) in environmental and seafood monitoring programs, in order to understand their risks on the ecosystem and human health.

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.

Historical and contemporary range expansion of an invasive mussel, Semimytlius algosus, in Angola and Namibia despite data scarcity in an infrequently surveyed region

Understanding the spread of invasive species in many regions is difficult because surveys are rare. Here, historical records of the invasive marine mussel, Semimytilus algosus, on the shores of Angola and Namibia are synthesised to re-construct its invasive history. Since this mussel was first discovered in Namibia about 90 years ago, it has spread throughout the western coast of southern Africa. By the late 1960s, the species was well established across a range of 1005 km of coastline in southern Angola and northern Namibia. Although only coarse spatial resolution data are available since the 1990s, the distribution of S. algosus clearly increased substantially over the subsequent decades. Today, the species is distributed over 2785 km of coastline, appearing in southern Namibia in 2014, whence it spread across the border to northern South Africa in 2017, and in northern Angola in 2015. Conspicuously, its current range appears to be relatively contiguous across at least 810 km of shore in southern Angola and throughout Namibia, with isolated, spatially disjunct occurrences towards the southern and northern limits of its distribution. Despite there being few occurrence records that are unevenly distributed spatially and temporally, data for the distributional patterns of S. algosus in Angola and Namibia provide invaluable insights into how marine invasive species spread in developing regions that are infrequently monitored.

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.

Comparative analysis of the ribosomal DNA repeat unit (rDNA) of Perna viridis (Linnaeus, 1758) and Perna canaliculus (Gmelin, 1791)

Perna viridis and P. canaliculus are economically and ecologically important species of shellfish. In this study, the complete ribosomal DNA (rDNA) unit sequences of these species were determined for the first time. The gene order, 18S rRNA-internal transcribed spacer (ITS) 1-5.8S rRNA-ITS2-28S rRNA-intergenic spacer (IGS), was similar to that observed in other eukaryotes. The lengths of the P. viridis and P. canaliculus rDNA sequences ranged from 8,432 to 8,616 bp and from 7,597 to 7,610 bp, respectively, this variability was mainly attributable to the IGS region. The putative transcription termination site and initiation site were confirmed. Perna viridis and P. canaliculus rDNA contained two (length: 93 and 40 bp) and one (length: 131 bp) repeat motifs, respectively. Individual intra-species differences mainly involved the copy number of repeat units. In P. viridis, three cytosine-guanine (CpG) sites with sizes of 440, 1,075 and 537 bp were found to cover nearly the entire IGS sequence, whereas in P. canaliculus, two CpG islands with sizes of 361 and 484 bp were identified. The phylogenetic trees constructed with maximum likelihood and neighbour-joining methods and based on ITS sequences were identical and included three major clusters. Species of the same genus were easily clustered together.

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.

Genetic relationships among freshwater mussel species from fifteen Amazonian rivers and inferences on the evolution of the Hyriidae (Mollusca: Bivalvia: Unionida)

The current phylogenetic framework for the South American Hyriidae is solely based on morphological data. However, freshwater bivalve morphology is highly variable due to both genetic and environmental factors. The present study used both mitochondrial (COI and 16S) and nuclear (18S-ITS1) sequences in molecular phylogenetic analyses of nine Neotropical species of Hyriidae, collected from 15 South American rivers, and sequences of hyriids from Australia and New Zealand obtained from GenBank. The present molecular findings support traditional taxonomic proposals, based on morphology, for the South American subfamily Hyriinae, currently divided in three tribes: Hyriini, Castaliini and Rhipidodontini. Phylogenetic trees based on COI nucleotide sequences revealed at least four geographical groups of Castalia ambigua: northeast Amazon (Piriá, Tocantins and Caeté rivers), central Amazon, including C. quadrata (Amazon and Aripuanã rivers), north (Trombetas river), and C. ambigua from Peru. Genetic distances suggest that some specimens may be cryptic species. Among the Hyriini, a total evidence data set generated phylogenetic trees indicating that Paxyodon syrmatophorus and Prisodon obliquus are more closely related, followed by Triplodon corrugatus. The molecular clock, based on COI, agreed with the fossil record of Neotropical hyriids. The ancestor of both Australasian and Neotropical Hyriidae is estimated to have lived around 225million years ago.

Field Measurements Indicate Unexpected, Serious Underestimation of Mussel Heart Rates and Thermal Tolerance by Laboratory Studies

Attempts to predict the response of species to long-term environmental change are generally based on extrapolations from laboratory experiments that inevitably simplify the complex interacting effects that occur in the field. We recorded heart rates of two genetic lineages of the brown mussel Perna perna over a full tidal cycle in-situ at two different sites in order to evaluate the cardiac responses of the two genetic lineages present on the South African coast to temperature and the immersion/emersion cycle. "Robomussel" temperature loggers were used to monitor thermal conditions at the two sites over one year. Comparison with live animals showed that robomussels provided a good estimate of mussel body temperatures. A significant difference in estimated body temperatures was observed between the sites and the results showed that, under natural conditions, temperatures regularly approach or exceed the thermal limits of P. perna identified in the laboratory. The two P. perna lineages showed similar tidal and diel patterns of heart rate, with higher cardiac activity during daytime immersion and minimal values during daytime emersion. Comparison of the heart rates measured in the field with data previously measured in the laboratory indicates that laboratory results seriously underestimate heart rate activity, by as much as 75%, especially during immersion. Unexpectedly, field estimates of body temperatures indicated an ability to tolerate temperatures considered lethal on the basis of laboratory measurements. This suggests that the interaction of abiotic conditions in the field does not necessarily raise vulnerability to high temperatures.

The complete mitochondrial genome of the golden mussel Limnoperna fortunei and comparative mitogenomics of Mytilidae

Here we describe the mitochondrial genome of the golden mussel Limnoperna fortunei, an Asian bivalve which has become one of the most aggressive invasive species in Japan and South America. The mitochondrial genome of L. fortunei does not present conserved gene arrangement when compared to the other Mytilidae species suggesting a high degree of gene recombination in the mitochondria of this clade. In addition, the golden mussel mitogenome encodes two copies of tRNA(Lys) and presents a putative pseudogene for the atp8 gene sequence that encodes a 27 amino acid peptide containing an in-frame stop codon. The presence of this pseudogene raises the question as to whether atp8 is encoded in some bivalve mitochondrial genomes or not. The phylogenetic analysis of all complete mitochondrial genomes available from Mytilidae mussels confirmed the close evolutionary relationships among bivalves from the genus Mytilys and placed L. fortunei coming from a more ancestral branch on the family. The supermatrix phylogeny described used the concatenation of all 12 genes from the mitochondria and disputed the monophyly of the genus Perna, as Perna perna was shown to be more closely related to Brachidontes exustus than to Perna viridis. The comparative analysis of mitogenome synteny also confirmed the polyphyly of the genus Perna. The complete and annotated mitogenome has been published in GenBank under the accession number KP756905.

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

Background

Intraspecific variability is seen as a central component of biodiversity. We investigated genetic differentiation, contemporary patterns of demographic connectivity and intraspecific variation of adaptive behavioural traits in two lineages of an intertidal mussel (Perna perna) across a tropical/subtropical biogeographic transition.

Results

Microsatellite analyses revealed clear genetic differentiation between western (temperate) and eastern (subtropical/tropical) populations, confirming divergence previously detected with mitochondrial (COI) and nuclear (ITS) markers.

Gene flow between regions was predominantly east-to-west and was only moderate, with higher heterozygote deficiency where the two lineages co-occur. This can be explained by differential selection and/or oceanographic dynamics acting as a barrier to larval dispersal.

Common garden experiments showed that gaping (periodic closure and opening of the shell) and attachment to the substratum differed significantly between the two lineages. Western individuals gaped more and attached less strongly to the substratum than eastern ones.

Conclusions

These behavioural differences are consistent with the geographic and intertidal distributions of each lineage along sharp environmental clines, indicating their strong adaptive significance. We highlight the functional role of diversity below the species level in evolutionary trends and the need to understand this when predicting biodiversity responses to environmental change.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0366-5) contains supplementary material, which is available to authorized users.