Biofouling sponges as natural eDNA samplers for marine vertebrate biodiversity monitoring.

Environmental DNA (eDNA) analysis has now become a core approach in marine biodiversity research, which typically involves the collection of water or sediment samples. Yet, recently, filter-feeding organisms have received much attention for their potential role as natural eDNA samplers. While the indiscriminate use of living organisms as 'sampling tools' might in some cases raise conservation concerns, there are instances in which highly abundant sessile organisms may become a nuisance as biofouling on artificial marine structures. Here we demonstrate how a sea sponge species that colonizes the moorings of the world's largest curtain of hydroacoustic receivers can become a powerful natural collector of fish biodiversity information. By sequencing eDNA extracted from Vazella pourtalesii retrieved from moorings during routine biofouling maintenance, we detected 23 species of marine fish and mammals, compared to 19 and 15 species revealed by surface and bottom water eDNA respectively, and 28 species captured by groundfish survey in the surrounding area, which are more ecologically impactful and involve higher additional costs. Sponge-based species inventories proved at least as informative as those obtained by traditional survey methods, and are also able to detect seasonal differences in fish assemblages. We conclude that opportunistic sampling of marine sponge biofouling may become an efficient way to document and monitor biodiversity in our rapidly changing oceans.

Genome-wide DNA methylation predicts environmentally driven life history variation in a marine fish.

Epigenetic modifications are thought to be one of the molecular mechanisms involved in plastic adaptive responses to environmental variation. However, studies reporting associations between genome-wide epigenetic changes and habitat-specific variations in life history traits (e.g., lifespan, reproduction) are still scarce, likely due to the recent application of methylome resequencing methods to non-model species. In this study, we examined associations between whole genome DNA methylation and environmentally driven life history variation in 2 lineages of a marine fish, the capelin (Mallotus villosus), from North America and Europe. In both lineages, capelin harbor 2 contrasting life history tactics (demersal vs. beach-spawning). Performing whole genome and methylome sequencing, we showed that life history tactics are associated with epigenetic changes in both lineages, though the effect was stronger in European capelin. Genetic differentiation between the capelin harboring different life history tactics was negligible, but we found genome-wide methylation changes in both lineages. We identified 9,125 European and 199 North American differentially methylated regions (DMRs) due to life history. Gene ontology (GO) enrichment analysis for both lineages revealed an excess of terms related to neural function. Our results suggest that environmental variation causes important epigenetic changes that are associated with contrasting life history tactics in lineages with divergent genetic backgrounds, with variable importance of genetic variation in driving epigenetic variation. Our study emphasizes the potential role of genome-wide epigenetic variation in adaptation to environmental variation.

Biodiversity, environmental drivers, and sustainability of the global deep-sea sponge microbiome.

In the deep ocean symbioses between microbes and invertebrates are emerging as key drivers of ecosystem health and services. We present a large-scale analysis of microbial diversity in deep-sea sponges (Porifera) from scales of sponge individuals to ocean basins, covering 52 locations, 1077 host individuals translating into 169 sponge species (including understudied glass sponges), and 469 reference samples, collected anew during 21 ship-based expeditions. We demonstrate the impacts of the sponge microbial abundance status, geographic distance, sponge phylogeny, and the physical-biogeochemical environment as drivers of microbiome composition, in descending order of relevance. Our study further discloses that fundamental concepts of sponge microbiology apply robustly to sponges from the deep-sea across distances of >10,000 km. Deep-sea sponge microbiomes are less complex, yet more heterogeneous, than their shallow-water counterparts. Our analysis underscores the uniqueness of each deep-sea sponge ground based on which we provide critical knowledge for conservation of these vulnerable ecosystems.

Exploring ecosystem-based management in the North Atlantic.

The United States, the EU and Canada established a trilateral working group on the ecosystem approach to ocean health and stressors under the Atlantic Ocean Research Alliance. Recognizing the Atlantic Ocean as a shared resource and responsibility, the working group sought to advance understanding of the Atlantic Ocean and its dynamic systems to improve ocean health, enhance ocean stewardship and promote the sustainable use and management of its resources. This included consideration of multiple ocean-use sectors such as fishing, shipping, tourism and offshore energy. The working group met for 4 years and worked through eight steps that covered the development of common language as a basis for collaboration, challenges of stakeholder engagement, review of the governance mandates, exploring the links between sectors and ecosystems effects, identifying gaps in knowledge and uptake of science, identification of tools for ecosystem-based management, customary best practice for tool development and communication of key research priorities. The key findings were that ecosystem-based management enables new benefits and opportunities, and that we need to make the business case. Further findings were that adequate mandates and effective tools exist for ecosystem-based management, and that ecosystem-based management urgently requires integration of human dimensions, so we must diversify the conversation. In addition, it was found that stakeholders do not see their stake in ecosystem-based management, so greater engagement with stakeholders and targeting of ocean literacy is required and a sustainable future requires a sustained investment in ecosystem-based management, so long-term commitment is key.

Oceanographic setting influences the prokaryotic community and metabolome in deep-sea sponges.

Marine sponges (phylum Porifera) are leading organisms for the discovery of bioactive compounds from nature. Their often rich and species-specific microbiota is hypothesised to be producing many of these compounds. Yet, environmental influences on the sponge-associated microbiota and bioactive compound production remain elusive. Here, we investigated the changes of microbiota and metabolomes in sponges along a depth range of 1232 m. Using 16S rRNA gene amplicon sequencing and untargeted metabolomics, we assessed prokaryotic and chemical diversities in three deep-sea sponge species: Geodia barretti, Stryphnus fortis, and Weberella bursa. Both prokaryotic communities and metabolome varied significantly with depth, which we hypothesized to be the effect of different water masses. Up to 35.5% of microbial ASVs (amplicon sequence variants) showed significant changes with depth while phylum-level composition of host microbiome remained unchanged. The metabolome varied with depth, with relative quantities of known bioactive compounds increasing or decreasing strongly. Other metabolites varying with depth were compatible solutes regulating osmolarity of the cells. Correlations between prokaryotic community and the bioactive compounds in G. barretti suggested members of Acidobacteria, Proteobacteria, Chloroflexi, or an unclassified prokaryote as potential producers.

Differential processing of dissolved and particulate organic matter by deep-sea sponges and their microbial symbionts.

Deep-sea sponges create hotspots of biodiversity and biological activity in the otherwise barren deep-sea. However, it remains elusive how sponge hosts and their microbial symbionts acquire and process food in these food-limited environments. Therefore, we traced the processing (i.e. assimilation and respiration) of 13C- and 15N-enriched dissolved organic matter (DOM) and bacteria by three dominant North Atlantic deep-sea sponges: the high microbial abundance (HMA) demosponge Geodia barretti, the low microbial abundance (LMA) demosponge Hymedesmia paupertas, and the LMA hexactinellid Vazella pourtalesii. We also assessed the assimilation of both food sources into sponge- and bacteria-specific phospholipid-derived fatty acid (PLFA) biomarkers. All sponges were capable of assimilating DOM as well as bacteria. However, processing of the two food sources differed considerably between the tested species: the DOM assimilation-to-respiration efficiency was highest for the HMA sponge, yet uptake rates were 4-5 times lower compared to LMA sponges. In contrast, bacteria were assimilated most efficiently and at the highest rate by the hexactinellid compared to the demosponges. Our results indicate that phylogeny and functional traits (e.g., abundance of microbial symbionts, morphology) influence food preferences and diet composition of sponges, which further helps to understand their role as key ecosystem engineers of deep-sea habitats.

Microbial Strategies for Survival in the Glass Sponge Vazella pourtalesii.

Few studies have explored the microbiomes of glass sponges (Hexactinellida). The present study seeks to elucidate the composition of the microbiota associated with the glass sponge Vazella pourtalesii and the functional strategies of the main symbionts. We combined microscopic approaches with metagenome-guided microbial genome reconstruction and amplicon community profiling toward this goal. Microscopic imaging revealed that the host and microbial cells appeared within dense biomass patches that are presumably syncytial tissue aggregates. Based on abundances in amplicon libraries and metagenomic data, SAR324 bacteria, Crenarchaeota, Patescibacteria, and Nanoarchaeota were identified as abundant members of the V. pourtalesii microbiome; thus, their genomic potentials were analyzed in detail. A general pattern emerged in that the V. pourtalesii symbionts had very small genome sizes, in the range of 0.5 to 2.2 Mb, and low GC contents, even below those of seawater relatives. Based on functional analyses of metagenome-assembled genomes (MAGs), we propose two major microbial strategies: the "givers," namely, Crenarchaeota and SAR324, heterotrophs and facultative anaerobes, produce and partly secrete all required amino acids and vitamins. The "takers," Nanoarchaeota and Patescibacteria, are anaerobes with reduced genomes that tap into the microbial community for resources, e.g., lipids and DNA, likely using pilus-like structures. We posit that the existence of microbial cells in sponge syncytia together with the low-oxygen conditions in the seawater environment are factors that shape the unique compositional and functional properties of the microbial community associated with V. pourtalesii IMPORTANCE We investigated the microbial community of V. pourtalesii that forms globally unique, monospecific sponge grounds under low-oxygen conditions on the Scotian Shelf, where it plays a key role in its vulnerable ecosystem. The microbial community was found to be concentrated within biomass patches and is dominated by small cells (<1 µm). MAG analyses showed consistently small genome sizes and low GC contents, which is unusual compared to known sponge symbionts. These properties, as well as the (facultatively) anaerobic metabolism and a high degree of interdependence between the dominant symbionts regarding amino acid and vitamin synthesis, are likely adaptations to the unique conditions within the syncytial tissue of their hexactinellid host and the low-oxygen environment.

The BenBioDen database, a global database for meio-, macro- and megabenthic biomass and densities.

Benthic fauna refers to all fauna that live in or on the seafloor, which researchers typically divide into size classes meiobenthos (32/64 µm-0.5/1 mm), macrobenthos (250 µm-1 cm), and megabenthos (>1 cm). Benthic fauna play important roles in bioturbation activity, mineralization of organic matter, and in marine food webs. Evaluating their role in these ecosystem functions requires knowledge of their global distribution and biomass. We therefore established the BenBioDen database, the largest open-access database for marine benthic biomass and density data compiled so far. In total, it includes 11,792 georeferenced benthic biomass and 51,559 benthic density records from 384 and 600 studies, respectively. We selected all references following the procedure for systematic reviews and meta-analyses, and report biomass records as grams of wet mass, dry mass, or ash-free dry mass, or carbon per m2 and as abundance records as individuals per m2. This database provides a point of reference for future studies on the distribution and biomass of benthic fauna.

Climate-induced changes in the suitable habitat of cold-water corals and commercially important deep-sea fishes in the North Atlantic.

The deep sea plays a critical role in global climate regulation through uptake and storage of heat and carbon dioxide. However, this regulating service causes warming, acidification and deoxygenation of deep waters, leading to decreased food availability at the seafloor. These changes and their projections are likely to affect productivity, biodiversity and distributions of deep-sea fauna, thereby compromising key ecosystem services. Understanding how climate change can lead to shifts in deep-sea species distributions is critically important in developing management measures. We used environmental niche modelling along with the best available species occurrence data and environmental parameters to model habitat suitability for key cold-water coral and commercially important deep-sea fish species under present-day (1951-2000) environmental conditions and to project changes under severe, high emissions future (2081-2100) climate projections (RCP8.5 scenario) for the North Atlantic Ocean. Our models projected a decrease of 28%-100% in suitable habitat for cold-water corals and a shift in suitable habitat for deep-sea fishes of 2.0°-9.9° towards higher latitudes. The largest reductions in suitable habitat were projected for the scleractinian coral Lophelia pertusa and the octocoral Paragorgia arborea, with declines of at least 79% and 99% respectively. We projected the expansion of suitable habitat by 2100 only for the fishes Helicolenus dactylopterus and Sebastes mentella (20%-30%), mostly through northern latitudinal range expansion. Our results projected limited climate refugia locations in the North Atlantic by 2100 for scleractinian corals (30%-42% of present-day suitable habitat), even smaller refugia locations for the octocorals Acanella arbuscula and Acanthogorgia armata (6%-14%), and almost no refugia for P. arborea. Our results emphasize the need to understand how anticipated climate change will affect the distribution of deep-sea species including commercially important fishes and foundation species, and highlight the importance of identifying and preserving climate refugia for a range of area-based planning and management tools.

Active buoyancy adjustment increases dispersal potential in benthic marine animals.

While the study of dispersal and connectivity in the ocean typically centres on pelagic species and planktonic larval stages of benthic species, the present work explores an overlooked locomotor means in post-settlement benthic stages that redefines their dispersal potential. Members of the echinoderm class Holothuroidea colonize a diversity of marine environments world-wide, where they play key ecological and economical roles, making their conservation a priority. Holothuroids are commonly called sea cucumbers or sea slugs to reflect their slow movements and are assumed to disperse chiefly through pelagic larvae. The present study documents and explores their unexpected ability to actively modify their buoyancy, leading them to tumble or float at speeds orders of magnitudes faster than through benthic crawling. Two focal species representing different taxonomic orders, geographic distributions and reproductive strategies were studied over several years. Active buoyancy adjustment (ABA) was achieved through a rapid increase in water-to-flesh ratio by up to 740%, leading to bloating, and simultaneously detachment from the substrate. It occurred as early as 6 months post settlement in juveniles and was recorded in wild adult populations. In experimental trials, ABA was triggered by high conspecific density, decreasing salinity and increasing water turbidity. Based on field video footage, ABA-assisted movements generated speeds of up to 90 km/day. These findings imply that displacement during planktonic larval stages may not supersede the locomotor capacity of benthic stages, challenging the notion of sedentarity. Combining the present results and anecdotal reports, ABA emerges as a generalized means of dispersal among benthic animals, with critical implications for world-wide management and conservation of commercially and ecologically significant species.

Predicted distribution of the glass sponge Vazella pourtalesi on the Scotian Shelf and its persistence in the face of climatic variability.

Emerald Basin on the Scotian Shelf off Nova Scotia, Canada, is home to a globally unique aggregation of the glass sponge Vazella pourtalesi, first documented in the region in 1889. In 2009, Fisheries and Oceans Canada (DFO) implemented two Sponge Conservation Areas to protect these sponge grounds from bottom fishing activities. Together, the two conservation areas encompass 259 km2. In order to ascertain the degree to which the sponge grounds remain unprotected, we modelled the presence probability and predicted range distribution of V. pourtalesi on the Scotian Shelf using random forest modelling on presence-absence records. With a high degree of accuracy the random forest model predicted the highest probability of occurrence of V. pourtalesi in the inner basins on the central Scotian Shelf, with lower probabilities at the shelf break and in the Fundian and Northeast Channels. Bottom temperature was the most important determinant of its distribution in the model. Although the two DFO Sponge Conservation Areas protect some of the more significant concentrations of V. pourtalesi, much of its predicted distribution remains unprotected (over 99%). Examination of the hydrographic conditions in Emerald Basin revealed that the V. pourtalesi sponge grounds are associated with a warmer and more saline water mass compared to the surrounding shelf. Reconstruction of historical bottom temperature and salinity in Emerald Basin revealed strong multi-decadal variability, with average bottom temperatures varying by 8°C. We show that this species has persisted in the face of this climatic variability, possibly indicating how it will respond to future climate change.

Barcoding Atlantic Canada's mesopelagic and upper bathypelagic marine fishes.

DNA barcode sequences were developed from 557 mesopelagic and upper bathypelagic teleost specimens collected in waters off Atlantic Canada. Confident morphological identifications were available for 366 specimens, of 118 species and 93 genera, which yielded 328 haplotypes. Five of the species were novel to the Barcode of Life Database (BOLD). Most of the 118 species conformed to expectations of monophyly and the presence of a "barcode gap", though some known weaknesses in existing taxonomy were confirmed and a deficiency in published keys was revealed. Of the specimens for which no firm morphological identification was available, 156 were successfully identified to species, and a further 11 to genus, using their barcode sequences and a combination of distance- and character-based methods. The remaining 24 specimens were from species for which no reference barcode is yet available or else ones confused by apparent misidentification of publicly available sequences in BOLD. Addition of the new sequences to those previously in BOLD contributed support to recent taxonomic revisions of Chiasmodon and Poromitra, while it also revealed 18 cases of potential cryptic speciation. Most of the latter appear to result from genetic divergence among populations in different ocean basins, while the general lack of strong horizontal environmental gradients within the deep sea has allowed morphology to be conserved. Other examples of divergence appear to distinguish individuals living under the sub-tropical gyre of the North Atlantic from those under that ocean's sub-polar gyre. In contrast, the available sequences for two myctophid species, Benthosema glaciale and Notoscopelus elongatus, showed genetic structuring on finer geographic scales. The observed structure was not consistent with recent suggestions that "resident" populations of myctophids can maintain allopatry despite the mixing of ocean waters. Rather, it indicates that the very rapid speciation characteristic of the Myctophidae is both on-going and detectable using barcodes.

Influence of Benthic Macrofauna as a Spatial Structuring Agent for Juvenile Haddock (Melanogrammus aeglefinus) on the Eastern Scotian Shelf, Atlantic Canada.

We examined the habitat of juvenile haddock on the eastern Scotian Shelf (off Nova Scotia, Canada) in relation to grab-sampled benthic macrofaunal invertebrate species assemblages in order to determine whether there were significant differences in benthic macrofauna between areas of historically persistent high and low juvenile haddock abundance. Our analyses were conducted over two spatial scales in each of two years: among banks (Emerald, Western and Sable Island), approximately 60 km distant from each other, and between areas of high and low juvenile haddock abundance at distances of 10 to 30 km-all in an area that had not experienced groundfishing in the decade prior to sampling. We also examined fine-scale (10s of metres) within-site variability in the macrofauna and used surficial sediment characteristics, along with hydrographic variables, to identify environmental correlates. PERMANOVA identified statistically significant differences in biomass, density and composition of the benthos associated with juvenile haddock abundance; however it was difficult to determine whether the results had biological relevance. Post hoc tests showed that these differences occurred only on Sable Island Bank where both fish and benthos may have been independently responding to sediment type which was most different there (100% sand in the area of low haddock abundance vs. 22% gravel in the area of high haddock abundance). In total, 383 benthic taxa representing 13 phyla were identified. Annelida was the most specious phylum (36.29% of taxa, representing 33 families), followed by Arthropoda (with Crustaceans, mostly Amphipoda, accounting for 25.07% of the total number of taxa). The strongest pattern in the macrofauna was expressed at the largest scale, between banks, accounting for approximately 25% of the variation in the data. Emerald Bank, deeper, warmer and saltier than the Western and Sable Island Banks, had a distinctive fauna.

Morphological and genetic variation in North Atlantic giant file clams, Acesta spp. (Bivalvia: Limidae), with description of a new cryptic species in the northwest Atlantic.

We analyze the morphological and genetic variability within and between seven species of Acesta and specimens recently collected in the northwest Atlantic using traditional morphological measurements, landmark-based geometric morphometrics, and the cytochrome oxidase subunit I (COI) gene sequences, with particular emphasis on North Atlantic species. Shell morphology and external shell appearance do not allow reliable distinction between the widely recognized northeastern Atlantic A. excavata and other northwest Atlantic species or populations of Acesta, with the exception of A. oophaga. Similarly, shape analysis reveals a wide variability within northeastern Atlantic A. excavata, and significant morphological overlap with A. bullisi from the Gulf of Mexico and A. rathbuni from the southwestern Pacific and South China Sea. Specimens from the northwestern and Mid-Atlantic display shell shapes marginally similar to that of A. excavata. These differences are at least partly related to anterior or posterior shifting of the shell body and to the opposite shifting of the hinge line/dorsal region and upper lunule. These morphological variations, along with the midline-width-ratio, explain much of the variability extracted by principal component analysis. Results from a mitochondrial DNA barcode approach (COI), however, suggest that the northwest Atlantic specimens belong to a new species for which we propose the name Acesta cryptadelphe sp. nov. Differences in larval shell sizes between northeastern and northwestern Atlantic specimens are consistent with this result.

Genetic structure of capelin (Mallotus villosus) in the northwest Atlantic Ocean.

Capelin (Mallotus villosus) is a commercially exploited, key forage-fish species found in the boreal waters of the North Pacific and North Atlantic Oceans. We examined the population structure of capelin throughout their range in the Canadian northwest Atlantic Ocean using genetic-based methods. Capelin collected at ten beach and five demersal spawning locations over the period 2002 through 2008 (N = 3,433 fish) were genotyped using six polymorphic microsatellite loci. Temporally distinct samples were identified at three beach spawning locations: Chance Cove, Little Lawn and Straitsview, Newfoundland. Four capelin stocks are assumed for fisheries management in the northwest Atlantic Ocean based on meristics, morphometrics, tag returns, and seasonal distribution patterns. Our results suggested groupings that were somewhat different than the assumed structure, and indicate at least seven genetically defined populations arising from two ancestral populations. The spatial mosaic of capelin from each of the two basal cluster groups explains much of the observed geographic variability amongst neighbouring samples. The genetic-defined populations were resolved at Jost's Dest ≥ 0.01 and were composed of fish collected 1) in the Gulf of St. Lawrence, 2) along the south and east coasts of Newfoundland, 3) along coastal northern Newfoundland and southern Labrador, 4) along coastal northern Labrador, 5) near the Saguenay River, and at two nearshore demersal spawning sites, 6) one at Grebes Nest off Bellevue Beach on the east coast of Newfoundland, and 7) one off the coast of Labrador at Domino Run. Moreover, the offshore demersal spawners on the Scotian Shelf and Southeast Shoal appeared to be related to the inshore demersal spawners at Grebes Nest and in Domino Run and to beach spawners from the Gulf of St. Lawrence.

Kernel density surface modelling as a means to identify significant concentrations of vulnerable marine ecosystem indicators.

The United Nations General Assembly Resolution 61/105, concerning sustainable fisheries in the marine ecosystem, calls for the protection of vulnerable marine ecosystems (VME) from destructive fishing practices. Subsequently, the Food and Agriculture Organization (FAO) produced guidelines for identification of VME indicator species/taxa to assist in the implementation of the resolution, but recommended the development of case-specific operational definitions for their application. We applied kernel density estimation (KDE) to research vessel trawl survey data from inside the fishing footprint of the Northwest Atlantic Fisheries Organization (NAFO) Regulatory Area in the high seas of the northwest Atlantic to create biomass density surfaces for four VME indicator taxa: large-sized sponges, sea pens, small and large gorgonian corals. These VME indicator taxa were identified previously by NAFO using the fragility, life history characteristics and structural complexity criteria presented by FAO, along with an evaluation of their recovery trajectories. KDE, a non-parametric neighbour-based smoothing function, has been used previously in ecology to identify hotspots, that is, areas of relatively high biomass/abundance. We present a novel approach of examining relative changes in area under polygons created from encircling successive biomass categories on the KDE surface to identify "significant concentrations" of biomass, which we equate to VMEs. This allows identification of the VMEs from the broader distribution of the species in the study area. We provide independent assessments of the VMEs so identified using underwater images, benthic sampling with other gear types (dredges, cores), and/or published species distribution models of probability of occurrence, as available. For each VME indicator taxon we provide a brief review of their ecological function which will be important in future assessments of significant adverse impact on these habitats here and elsewhere.

Modeling the distribution of Geodia sponges and sponge grounds in the Northwest Atlantic.

Deep-sea sponge grounds provide structurally complex habitat for fish and invertebrates and enhance local biodiversity. They are also vulnerable to bottom-contact fisheries and prime candidates for Vulnerable Marine Ecosystem designation and related conservation action. This study uses species distribution modeling, based on presence and absence observations of Geodia spp. and sponge grounds derived from research trawl catches, as well as spatially continuous data on the physical and biological ocean environment derived from satellite data and oceanographic models, to model the distribution of Geodia sponges and sponge grounds in the Northwest Atlantic. Most models produce excellent fits with validation data although fits are reduced when models are extrapolated to new areas, especially when oceanographic regimes differ between areas. Depth and minimum bottom salinity were important predictors in most models, and a Geodia spp. minimum bottom salinity tolerance threshold in the 34.3-34.8 psu range was hypothesized on the basis of model structure. The models indicated two currently unsampled regions within the study area, the deeper parts of Baffin Bay and the Newfoundland and Labrador slopes, where future sponge grounds are most likely to be found.

Proteomic analysis of eggs from Mytilus edulis females differing in mitochondrial DNA transmission mode.

Many bivalves have an unusual mechanism of mitochondrial DNA (mtDNA) inheritance called doubly uniparental inheritance (DUI) in which distinctly different genomes are inherited through the female (F genome) and male (M genome) lineages. In fertilized eggs that will develop into male embryos, the sperm mitochondria remain in an aggregation, which is believed to be delivered to the primordial germ cells and passed to the next generation through the sperm. In fertilized eggs that will develop into female embryos, the sperm mitochondria are dispersed throughout the developing embryo and make little if any contribution to the next generation. The frequency of embryos with the aggregated or dispersed mitochondrial type varies among females. Previous models of DUI have predicted that maternal nuclear factors cause molecular differences among unfertilized eggs from females producing embryos with predominantly dispersed or aggregated mitochondria. We test this hypothesis using females of each of the two types from a natural population. We have found small, yet detectable, differences of the predicted type at the proteome level. We also provide evidence that eggs of females giving the dispersed pattern have consistently lower expression for different proteasome subunits than eggs of females giving the aggregated pattern. These results, combined with those of an earlier study in which we used hatchery lines of Mytilus, and with a transcriptomic study in a clam that has the DUI system of mtDNA transmission, reinforce the hypothesis that the ubiquitin-proteasome system plays a key role in the mechanism of DUI and sex determination in bivalves. We also report that eggs of females giving the dispersed pattern have higher expression for arginine kinase and enolase, enzymes involved in energy production, whereas ferritin, which is involved in iron homeostasis, has lower expression. We discuss these results in the context of genetic models for DUI and suggest experimental methods for further understanding the role of these proteins in DUI.

Expression of M6 and M7 lysin in Mytilus edulis is not restricted to sperm, but occurs also in oocytes and somatic tissue of males and females.

Sperm proteins of marine sessile invertebrates have been extensively studied to understand the molecular basis of reproductive isolation. Apart from molecules such as bindin of sea urchins or lysin of abalone species, the acrosomal protein M7 lysin of Mytilus edulis has been analyzed. M7 lysin was found to be under positive selection, but mechanisms driving the evolution of this protein are not fully understood. To explore functional aspects, this study investigated the protein expression pattern of M7 and M6 lysin in gametes and somatic tissue of male and female M. edulis. The study employs a previously published monoclonal antibody (G26-AG8) to investigate M6 and M7 lysin protein expression, and explores expression of both genes. It is shown that these proteins and their encoding genes are expressed in gametes and somatic tissue of both sexes. This is in contrast to sea urchin bindin and abalone lysin, in which gene expression is strictly limited to males. Although future studies need to clarify the functional importance of both acrosomal proteins in male and female somatic tissue, new insights into the evolution of sperm proteins in marine sessile invertebrates are possible. This is because proteins with male-specific expression (bindin, lysin) might evolve differently than proteins with expression in both sexes (M6/M7 lysin), and the putative function of both proteins in females opens the possibility that the evolution of M6/M7 lysin is under sexual antagonistic selection, for example, mutations beneficial to the acrosomal function that are less beneficial the function in somatic tissue of females.