Quantification of marine benthic communities with metabarcoding.

DNA metabarcoding methods have been implemented in studies aimed at detecting and quantifying marine benthic biodiversity. In such surveys, universal barcodes are amplified and sequenced from environmental DNA. To quantify biodiversity with DNA metabarcoding, a relation between the number of DNA sequences of a species and its biomass and/or the abundance is required. However, this relationship is complicated by many factors, and it is often unknown. In this study, we validate estimates of biomass and abundance from molecular approaches with those from the traditional morphological approach. Abundance and biomass were quantified from 126 samples of benthic intertidal mudflat using traditional morphological approaches and compared with frequency of occurrence and relative read abundance estimates from a molecular approach. A relationship between biomass and relative read abundance was found for two widely dispersed annelid taxa (Pygospio and Scoloplos). None of the other taxons, however, showed such a relationship. We discuss how quantification of abundance and biomass using molecular approaches are hampered by the ecology of DNA i.e. all the processes that determine the amount of DNA in the environment, including the ecology of the benthic species as well as the compositional nature of sequencing data.

Marine stepping-stones: Connectivity of Mytilus edulis populations between offshore energy installations.

Recent papers have suggested that epifaunal organisms use artificial structures as stepping-stones to spread to areas that are too distant to reach in a single generation. With thousands of artificial structures present in the North Sea, we test the hypothesis that these structures are connected by water currents and act as an interconnected reef. Population genetic structure of the blue mussel, Mytilus edulis, was expected to follow a pattern predicted by a particle tracking model (PTM). Correlation between population genetic differentiation, based on microsatellite markers, and particle exchange was tested. Specimens of M. edulis were found at each location, although the PTM indicated that locations >85 km offshore were isolated from coastal subpopulations. The fixation coefficient FST correlated with the number of arrivals in the PTM. However, the number of effective migrants per generation as inferred from coalescent simulations did not show a strong correlation with the arriving particles. Isolation by distance analysis showed no increase in isolation with increasing distance and we did not find clear structure among the populations. The marine stepping-stone effect is obviously important for the distribution of M. edulis in the North Sea and it may influence ecologically comparable species in a similar way. In the absence of artificial shallow hard substrates, M. edulis would be unlikely to survive in offshore North Sea waters.

Global invasion genetics of two parasitic copepods infecting marine bivalves.

Invasive species, and especially invasive parasites, represent excellent models to study ecological and evolutionary mechanisms in the wild. To understand these processes, it is crucial to obtain more knowledge on the native range, invasion routes and invasion history of invasive parasites. We investigated the consecutive invasions of two parasitic copepods (Mytilicola intestinalis and Mytilicola orientalis) by combining an extensive literature survey covering the reported putative native regions and the present-day invaded regions with a global phylogeography of both species. The population genetic analyses based on partial COI sequences revealed significant population differentiation for M. orientalis within the native region in Japan, while introduced populations in North America and Europe could not be distinguished from the native ones. Thus, M. orientalis' invasion history resembles the genetic structure and recent spread of its principal host, the Pacific oyster, Crassostrea gigas, while M. intestinalis lacks population genetic structure and has an overall low genetic diversity. Therefore, the native origin of M. intestinalis remains unclear. With this study, we demonstrate that even highly related and biologically similar invasive species can differ in their invasion genetics. From this, we conclude that extrapolating invasion genetics dynamics from related invasive taxa may not always be possible.

How invasive oysters can affect parasite infection patterns in native mussels on a large spatial scale.

There are surprisingly few field studies on the role of invasive species on parasite infection patterns in native hosts. We investigated the role of invasive Pacific oysters (Magallana gigas) in determining parasite infection levels in native blue mussels (Mytilus edulis) in relation to other environmental and biotic factors. Using hierarchical field sampling covering three spatial scales along a large intertidal ecosystem (European Wadden Sea), we found strong spatial differences in infection levels of five parasite species associated with mussels and oysters. We applied mixed models to analyse the associations between parasite prevalence and abundance in mussels and oysters, and 12 biological and environmental factors. For each parasite-host relationship, an optimal model (either a null, one-factor or two-factor model) was selected based on AIC scores. We found that the density of invasive oysters contributed to three of the 12 models. Other biological factors such as host size (six models), and the density of target or alternative host species (five models) contributed more frequently to the best models. Furthermore, for parasite species infecting both mussels and oysters, parasite population densities were higher in native mussels, attributed to the higher densities of mussels. Our results indicate that invasive species can affect parasite infection patterns in native species in the field, but that their relative contribution may be further mediated by other biological and environmental parameters. These results stress the usefulness of large-scale field studies for detailed assessments of the mechanisms underlying the impacts of invasive species on native host communities.

Strong population structure but no equilibrium yet: Genetic connectivity and phylogeography in the kelp Saccharina latissima (Laminariales, Phaeophyta).

Kelp aquaculture is globally developing steadily as human food source, along with other applications. One of the newer crop species is Saccharina latissima, a northern hemisphere kelp inhabiting temperate to arctic rocky shores. To protect and document its natural genetic variation at the onset of this novel aquaculture, as well as increase knowledge on its taxonomy and phylogeography, we collected new genetic data, both nuclear and mitochondrial, and combined it with previous knowledge to estimate genetic connectivity and infer colonization history. Isolation-with-migration coalescent analyses demonstrate that gene flow among the sampled locations is virtually nonexistent. An updated scenario for the origin and colonization history of S. latissima is developed as follows: We propose that the species (or species complex) originated in the northwest Pacific, crossed to the northeast Pacific in the Miocene, and then crossed the Bering Strait after its opening ~5.5 Ma into the Arctic and northeast Atlantic. It subsequently crossed the Atlantic from east to west. During the Pleistocene, it was compressed in the south with evidence for northern refugia in Europe. Postglacial recolonization led to secondary contact in the Canadian Arctic. Saccharina cichorioides is shown to probably belong to the S. latissima species complex and to derive from ancestral populations in the Asian North Pacific. Our novel approach of comparing inferred gene flow based on coalescent analysis versus Wright's island model suggests that equilibrium levels of differentiation have not yet been reached in Europe and, hence, that genetic differentiation is expected to increase further if populations are left undisturbed.

Cryptic invasion of a parasitic copepod: Compromised identification when morphologically similar invaders co-occur in invaded ecosystems.

Despite their frequent occurrence and strong impacts on native biota, biological invasions can long remain undetected. One reason for this is that an invasive species can be morphologically similar to either native species or introduced species previously established in the same region, and thus be subject to mistaken identification. One recent case involves congeneric invasive parasites, copepods that now infect bivalve hosts along European Atlantic coasts, after having been introduced independently first from the Mediterranean Sea (Mytilicola intestinalis Steuer, 1902) and later from Japan (Mytilicola orientalis Mori, 1935). At least one report on M. intestinalis may have actually concerned M. orientalis, and M. orientalis thus qualifies as a "cryptic invader". Because these two parasitic copepods are morphologically similar, knowledge about their distribution, impact and interactions depends crucially on reliable species identification. In this study, we evaluated the reliability of morphological identification of these two species in parts of their invasive range in Europe (Dutch Delta and Wadden Sea) in comparison with molecular methods of well-established accuracy based on COI gene sequences and ITS1 restriction fragment length polymorphism. Based on seven easily measured or scored macro-morphological variables that were recorded for 182 individual copepods isolated from blue mussels (Mytilus edulis Linnaeus, 1758), principal component analysis showed two relatively distinct but overlapping morphological species groups for females, but no clear separation in males. Discriminant function analysis showed that the females can be discriminated reasonably well based on some of the morphological characteristics (identification error rate of 7%) while males cannot (error rate of 25%). The direction of the dorsolateral thoracic protuberances was identified as the most important trait for species discrimination, but among the morphological features checked, none could flawlessly discriminate between both species. We recommend the use of molecular techniques in future studies of invasive Mytilicola to reliably discriminate between the species. The morphological similarity of these two invaders suggests a more general problem of cryptic invasions and compromised identification of parasites in invaded ecosystems. This problem should be borne in mind whenever invasive parasites are investigated.

Dual transcriptomics reveals co-evolutionary mechanisms of intestinal parasite infections in blue mussels Mytilus edulis.

On theoretical grounds, antagonistic co-evolution between hosts and their parasites should be a widespread phenomenon but only received little empirical support so far. Consequently, the underlying molecular mechanisms and evolutionary steps remain elusive, especially in nonmodel systems. Here, we utilized the natural history of invasive parasites to document the molecular underpinnings of co-evolutionary trajectories. We applied a dual-species transcriptomics approach to experimental cross-infections of blue mussel Mytilus edulis hosts and their invasive parasitic copepods Mytilicola intestinalis from two invasion fronts in the Wadden Sea. We identified differentially regulated genes from an experimental infection contrast for hosts (infected vs. control) and a sympatry contrast (sympatric vs. allopatric combinations) for both hosts and parasites. The damage incurred by Mytilicola infection and the following immune response of the host were mainly reflected in cell division processes, wound healing, apoptosis and the production of reactive oxygen species (ROS). Furthermore, the functional coupling of host and parasite sympatry contrasts revealed the concerted regulation of chitin digestion by a Chitotriosidase 1 homolog in hosts with several cuticle proteins in the parasite. Together with the coupled regulation of ROS producers and antagonists, these genes represent candidates that mediate the different evolutionary trajectories within the parasite's invasion. The host-parasite combination-specific coupling of these effector mechanisms suggests that underlying recognition mechanisms create specificity and local adaptation. In this way, our study demonstrates the use of invasive species' natural history to elucidate molecular mechanisms of host-parasite co-evolution in the wild.

Trophic relationship between the invasive parasitic copepod Mytilicola orientalis and its native blue mussel (Mytilus edulis) host.

Invasive parasites can spill over to new hosts in invaded ecosystems with often unpredictable trophic relationships in the newly arising parasite-host interactions. In European seas, the intestinal copepod Mytilicola orientalis was co-introduced with Pacific oysters (Magallana gigas) and spilled over to native blue mussels (Mytilus edulis), with negative impacts on the condition of infected mussels. However, whether the parasite feeds on host tissue and/or stomach contents is yet unknown. To answer this question, we performed a stable isotope analysis in which we included mussel host tissue and the primary food sources of the mussels, microphytobenthos (MPB) and particulate organic matter (POM). The copepods were slightly enriched in δ15N (mean Δ15N ± s.d.; 1·22 ± 0·58‰) and δ13C (Δ13C 0·25 ± 0·32‰) with respect to their host. Stable isotope mixing models using a range of trophic fractionation factors indicated that host tissue was the main food resource with consistent additional contributions of MPB and POM. These results suggest that the trophic relationship of the invasive copepod with its mussel host is parasitic as well as commensalistic. Stable isotope studies such as this one may be a useful tool to unravel trophic relationships in new parasite-host associations in the course of invasions.

Where are the polyps? Molecular identification, distribution and population differentiation of Aurelia aurita jellyfish polyps in the southern North Sea area.

For many species of metagenic jellyfish the location of the benthic polyps is unknown. To gain insight in the distribution, species composition and population structure of scyphozoan jellyfish polyps in the southern North Sea area, polyp samples were collected from natural and artificial substrates (settling plates, marina floats and wrecks) at ten inshore locations in the Netherlands, seven offshore locations in the North Sea and in the Gullmar Fjord in Sweden. Polyps were identified to species level by sequencing both a fragment of 18S rDNA and a fragment of mitochondrial COI, and comparing these sequences to reference sequences available in GenBank and to newly obtained sequences from medusae collected in the area. All polyps sequenced did belong to Aurelia aurita. For this species, molecular diversity in mitochondrial COI was high, with 50 haplotypes among 183 polyps. Population differentiation was detected between the Dogger Bank and other-more coastal-locations, indicating extremely low connectivity. No significant differences were found between coastal samples. The location of polyps of Cyanea capillata, Cyanea lamarckii, Chrysaora hysoscella and Rhizostoma octopus in the study area remains unresolved.

Population Genetic Structure, Abundance, and Health Status of Two Dominant Benthic Species in the Saba Bank National Park, Caribbean Netherlands: Montastraea cavernosa and Xestospongia muta.

Saba Bank, a submerged atoll in the Caribbean Sea with an area of 2,200 km2, has attained international conservation status due to the rich diversity of species that reside on the bank. In order to assess the role of Saba Bank as a potential reservoir of diversity for the surrounding reefs, we examined the population genetic structure, abundance and health status of two prominent benthic species, the coral Montastraea cavernosa and the sponge Xestospongia muta. Sequence data were collected from 34 colonies of M. cavernosa (nDNA ITS1-5.8S-ITS2; 892 bp) and 68 X. muta sponges (mtDNA I3-M11 partition of COI; 544 bp) on Saba Bank and around Saba Island, and compared with published data across the wider Caribbean. Our data indicate that there is genetic connectivity between populations on Saba Bank and the nearby Saba Island as well as multiple locations in the wider Caribbean, ranging in distance from 100s-1000s km. The genetic diversity of Saba Bank populations of M. cavernosa (π = 0.055) and X. muta (π = 0.0010) was comparable to those in other regions in the western Atlantic. Densities and health status were determined along 11 transects of 50 m2 along the south-eastern rim of Saba Bank. The densities of M. cavernosa (0.27 ind. m-2, 95% CI: 0.12-0.52) were average, while the densities of X. muta (0.09 ind. m-2, 95% CI: 0.02-0.32) were generally higher with respect to other Caribbean locations. No disease or bleaching was present in any of the specimens of the coral M. cavernosa, however, we did observe partial tissue loss (77.9% of samples) as well as overgrowth (48.1%), predominantly by cyanobacteria. In contrast, the majority of observed X. muta (83.5%) showed signs of presumed bleaching. The combined results of apparent gene flow among populations on Saba Bank and surrounding reefs, the high abundance and unique genetic diversity, indicate that Saba Bank could function as an important buffer for the region. Either as a natural source of larvae to replenish genetic diversity or as a storehouse of diversity that can be utilized if needed for restoration practices.

An efficient method to find potentially universal population genetic markers, applied to metazoans.

BACKGROUND: Despite the impressive growth of sequence databases, the limited availability of nuclear markers that are sufficiently polymorphic for population genetics and phylogeography and applicable across various phyla restricts many potential studies, particularly in non-model organisms. Numerous introns have invariant positions among kingdoms, providing a potential source for such markers. Unfortunately, most of the few known EPIC (Exon Primed Intron Crossing) loci are restricted to vertebrates or belong to multigenic families. RESULTS: In order to develop markers with broad applicability, we designed a bioinformatic approach aimed at avoiding multigenic families while identifying intron positions conserved across metazoan phyla. We developed a program facilitating the identification of EPIC loci which allowed slight variation in intron position. From the Homolens databases we selected 29 gene families which contained 52 promising introns for which we designed 93 primer pairs. PCR tests were performed on several ascidians, echinoderms, bivalves and cnidarians. On average, 24 different introns per genus were amplified in bilaterians. Remarkably, five of the introns successfully amplified in all of the metazoan genera tested (a dozen genera, including cnidarians). The influence of several factors on amplification success was investigated. Success rate was not related to the phylogenetic relatedness of a taxon to the groups that most influenced primer design, showing that these EPIC markers are extremely conserved in animals. CONCLUSIONS: Our new method now makes it possible to (i) rapidly isolate a set of EPIC markers for any phylum, even outside the animal kingdom, and thus, (ii) compare genetic diversity at potentially homologous polymorphic loci between divergent taxa.

Different flooding responses in Rorippa amphibia and Rorippa sylvestris, and their modes of expression in F1 hybrids.

The river floodplain species Rorippa amphibia, Rorippa sylvestris, and their hybrid Rorippa x anceps were studied here, with the aim of identifying potential species differences with respect to flooding tolerance, and of assessing their expression in F1 hybrids. Parents and their F1 hybrids were subjected to three flooding treatments mimicking natural conditions, and growth-related and leaf morphological traits were compared. In contrast to R. sylvestris, R. amphibia responded to waterlogging by forming specialized roots, and its growth was not reduced. These traits were dominantly expressed in hybrids. Both species and the hybrids established shoot growth over 2 wk of complete submergence. Only in R. sylvestris was this not at the expense of root biomass, suggesting that R. sylvestris can photosynthesize underwater. Rorippa sylvestris also showed a hyponastic response. Hybrids were intermediate to the parents in this respect. This study shows that phenotypic expression of parental traits in F1 hybrids is mostly additive, but can also be dominant. This suggests that a large overlap in habitat use of parents and hybrids is likely. If such an overlap occurs, the main evolutionary consequences of hybridization in Rorippa will be the introgression of genes, as the hybrids are fully fertile.

Phylogeography of the common shrimp, Crangon crangon (L.) across its distribution range.

The common or brown shrimp Crangon crangon (L.) is a highly abundant and important taxon, both ecologically and commercially, yet knowledge on its population structure and historical biogeography is limited. We studied population genetic structure across the distribution range of this species by sequencing a 388 bp fragment of the cytochrome-c-oxidase I gene for 140 individuals from 25 locations. Strong population structuring and high levels of genetic diversity were observed. Four main phylogroups were uncovered: northeastern Atlantic, western Mediterranean, Adriatic Sea and Black Sea. Gene flow of these shrimp across known oceanographical barriers (e.g., the Strait of Gibraltar and/or Oran-Almeria front, Sicilian Straits, and Turkish Straits) is severely restricted. The oldest and most variable populations currently inhabit the western Mediterranean. The observed absence of structure across the entire northeastern Atlantic shelf is proposed not to be due to gene flow, but to relatively recent colonization following the glacial cycles of the late Pleistocene. Black Sea shrimp are currently disconnected from Mediterranean populations, and colonization is inferred, on the basis of coalescent analysis, to have happened relatively recently, but possibly earlier than 7000 years ago. We postulate the hypothesis that C. crangon survived the last brackish-water (<7 per thousand) period inside the Black Sea and/or one of the adjacent inland seas. We conclude that (1) common shrimp populations from different basins are strongly differentiated, (2) gene flow across basins is probably very limited, and (3) the biogeographic history of the taxon is largely in accordance with the geographic history of its distribution range. This study provides further evidence that high population connectivity of marine species (e.g., by policy makers) should not be assumed.

Genetic diversity in diploid vs. tetraploid Rorippa amphibia (Brassicaceae).

The frequency of polyploidy increases with latitude in the Northern Hemisphere, especially in deglaciated, recently colonized areas. The cause or causes of this pattern are largely unknown, but a greater genetic diversity of individual polyploid plants due to a doubled genome and/or a hybrid origin is seen as a likely factor underlying selective advantages related to life in extreme climates and/or colonization ability. A history of colonization in itself, as well as a recent origin, and possibly a limited number of polyploidization events would all predict less genetic diversity in polyploids than in diploids. The null hypothesis of higher gene diversity in polyploids has to date hardly been quantified and is here tested in self-incompatible Rorippa amphibia (Brassicaceae). The species occurs in diploid and tetraploid forms and displays clear geographical polyploidy in Europe. On the basis of eight microsatellite loci it can be concluded that the level of gene diversity is higher in tetraploids than in diploids, to an extent that is expected under neutral evolution when taking into account the larger effective population size in the doubled cytotype. There is thus no evidence for reduced genetic diversity in the tetraploids. The evidence presented here may mean that the tetraploids' origin is not recent, has not been affected by bottlenecks and/or that tetraploids were formed multiple times while an effect of introgression may also play a role.

Mosaic haploid-diploid embryos and polyspermy in the tellinid bivalve Macoma balthica.

We investigated meiosis, fertilization, and early development in eggs of the tellinid bivalve Macoma balthica (L.), which has external fertilization. Meiosis is standard but polyspermy is found to be very common. In all eight crosses examined, mosaic embryos consisting of a mixture of diploid (2n = 38) and haploid cells occur at a frequency ranging from 2.7 to 29.1%. The earliest mosaic found is in the two-cell stage. We propose that an androgenic haploid cell lineage can originate from one supernumerary sperm that decondenses into a functional haploid nucleus, starts mitosis, and is incorporated in the developing embryo.

Facilitation on an intertidal mudflat: the effect of siphon nipping by flatfish on burying depth of the bivalve Macoma balthica.

During deposit feeding on benthic micro-algae, the siphon of buried tellinid bivalves like Macoma balthica is vulnerable to nipping by plaice Pleuronectes platessa and other flatfish. We experimentally tested the hypothesis that siphon nippers facilitate predation on the entire bivalve (by shorebirds, for example) by inducing a decrease in burying depth. Three experiments in May and in September, during which likely siphon nippers (juvenile plaice P. platessa) were allowed to feed on M. balthica siphons, demonstrated that the bivalves indeed lost siphon mass and came closer to the surface. However, the strength of the burying response and the speed of recovery of siphons after nipping varied greatly between experiments, partly as a consequence of body condition-related differences in initial burying depth. Our experimental study confirms that siphon nippers can enhance the availability of bivalve prey to probing predators, and will thus facilitate avian predation on intertidal flats.