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.

Cascabel: A Scalable and Versatile Amplicon Sequence Data Analysis Pipeline Delivering Reproducible and Documented Results.

Marker gene sequencing of the rRNA operon (16S, 18S, ITS) or cytochrome c oxidase I (CO1) is a popular means to assess microbial communities of the environment, microbiomes associated with plants and animals, as well as communities of multicellular organisms via environmental DNA sequencing. Since this technique is based on sequencing a single gene, or even only parts of a single gene rather than the entire genome, the number of reads needed per sample to assess the microbial community structure is lower than that required for metagenome sequencing. This makes marker gene sequencing affordable to nearly any laboratory. Despite the relative ease and cost-efficiency of data generation, analyzing the resulting sequence data requires computational skills that may go beyond the standard repertoire of a current molecular biologist/ecologist. We have developed Cascabel, a scalable, flexible, and easy-to-use amplicon sequence data analysis pipeline, which uses Snakemake and a combination of existing and newly developed solutions for its computational steps. Cascabel takes the raw data as input and delivers a table of operational taxonomic units (OTUs) or Amplicon Sequence Variants (ASVs) in BIOM and text format and representative sequences. Cascabel is a highly versatile software that allows users to customize several steps of the pipeline, such as selecting from a set of OTU clustering methods or performing ASV analysis. In addition, we designed Cascabel to run in any linux/unix computing environment from desktop computers to computing servers making use of parallel processing if possible. The analyses and results are fully reproducible and documented in an HTML and optional pdf report. Cascabel is freely available at Github: https://github.com/AlejandroAb/CASCABEL.

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.

Characterization and Temperature Dependence of Arctic Micromonas polaris Viruses.

Global climate change-induced warming of the Artic seas is predicted to shift the phytoplankton community towards dominance of smaller-sized species due to global warming. Yet, little is known about their viral mortality agents despite the ecological importance of viruses regulating phytoplankton host dynamics and diversity. Here we report the isolation and basic characterization of four prasinoviruses infectious to the common Arctic picophytoplankter Micromonas. We furthermore assessed how temperature influenced viral infectivity and production. Phylogenetic analysis indicated that the putative double-stranded DNA (dsDNA) Micromonas polaris viruses (MpoVs) are prasinoviruses (Phycodnaviridae) of approximately 120 nm in particle size. One MpoV showed intrinsic differences to the other three viruses, i.e., larger genome size (205 ± 2 vs. 191 ± 3 Kb), broader host range, and longer latent period (39 vs. 18 h). Temperature increase shortened the latent periods (up to 50%), increased the burst size (up to 40%), and affected viral infectivity. However, the variability in response to temperature was high for the different viruses and host strains assessed, likely affecting the Arctic picoeukaryote community structure both in the short term (seasonal cycles) and long term (global warming).

Virus production in phosphorus-limited Micromonas pusilla stimulated by a supply of naturally low concentrations of different phosphorus sources, far into the lytic cycle.

Earlier studies show that the proliferation of phytoplankton viruses can be inhibited by depletion of soluble reactive phosphorus (SRP; orthophosphate). In natural marine waters, phytoplankton phosphorus (P) availability is, however, largely determined by the supply rate of SRP (e.g. through remineralization) and potentially by the source of P as well (i.e. the utilization of soluble non-reactive P; SNP). Here we show how a steady low supply of P (mimicking natural P recycling) to virally infected P-limited Micromonas pusilla stimulates virus proliferation. Independent of the degree of P limitation prior to infection (0.32 and 0.97µmax chemostat cultures), SRP supply resulted in 2-fold higher viral burst sizes (viruses lysed per host cell) as compared with no addition (P starvation). Delaying these spikes during the infection cycle showed that the added SRP was utilized for extra M. pusilla virus (MpV) production far into the lytic cycle (18 h post-infection). Moreover, P-limited M. pusilla utilized several SNP compounds with high efficiency and with the same extent of burst size stimulation as for SRP. Finally, addition of virus-free MpV lysate (representing a complex SNP mixture) to newly infected cells enhanced MpV production, implicating host-associated alkaline phosphatase activity, and highlighting its important role in oligotrophic environments.

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.

Exposing the grey seal as a major predator of harbour porpoises.

Harbour porpoises (Phocoena phocoena) stranding in large numbers around the southern North Sea with fatal, sharp-edged mutilations have spurred controversy among scientists, the fishing industry and conservationists, whose views about the likely cause differ. The recent detection of grey seal (Halichoerus grypus) DNA in bite marks on three mutilated harbour porpoises, as well as direct observations of grey seal attacks on porpoises, have identified this seal species as a probable cause. Bite mark characteristics were assessed in a retrospective analysis of photographs of dead harbour porpoises that stranded between 2003 and 2013 (n = 1081) on the Dutch coastline. There were 271 animals that were sufficiently fresh to allow macroscopic assessment of grey seal-associated wounds with certainty. In 25% of these, bite and claw marks were identified that were consistent with the marks found on animals that had tested positive for grey seal DNA. Affected animals were mostly healthy juveniles that had a thick blubber layer and had recently fed. We conclude that the majority of the mutilated harbour porpoises were victims of grey seal attacks and that predation by this species is one of the main causes of death in harbour porpoises in The Netherlands. We provide a decision tree that will help in the identification of future cases of grey seal predation on porpoises.