Ultra-deep sequencing of foraminiferal microbarcodes unveils hidden richness of early monothalamous lineages in deep-sea sediments

Taxonomic and abundance biases affect the record of marine eukaryotic plankton communities in sediment DNA archives

Environmental DNA (eDNA) preserved in marine sediments is increasingly being used to study past ecosystems. However, little is known about how accurately marine biodiversity is recorded in sediment eDNA archives, especially planktonic taxa. Here, we address this question by comparing eukaryotic diversity in 273 eDNA samples from three water depths and the surface sediments of 24 stations in the Nordic Seas. Analysis of 18S-V9 metabarcoding data reveals distinct eukaryotic assemblages between water and sediment eDNA. Only 40% of Amplicon Sequence Variants (ASVs) detected in water were also found in sediment eDNA. Remarkably, the ASVs shared between water and sediment accounted for 80% of total sequence reads suggesting that a large amount of plankton DNA is transported to the seafloor, predominantly from abundant phytoplankton taxa. However, not all plankton taxa were equally archived on the seafloor. The plankton DNA deposited in the sediments was dominated by diatoms and showed an underrepresentation of certain nano- and picoplankton taxa (Picozoa or Prymnesiophyceae). Our study offers the first insights into the patterns of plankton diversity recorded in sediment in relation to seasonality and spatial variability of environmental conditions in the Nordic Seas. Our results suggest that the genetic composition and structure of the plankton community vary considerably throughout the water column and differ from what accumulates in the sediment. Hence, the interpretation of sedimentary eDNA archives should take into account potential taxonomic and abundance biases when reconstructing past changes in marine biodiversity.

Environmental DNA sheds new insight on molecular adaptation of foraminifera to temperature from laboratory-controlled culture experiment

Foraminifera is the most important temperature proxy of the ocean on long time scales. However, the absence of temperature-controlled experiments at different water depths hinders the advancement of paleotemperature reconstruction with foraminifera from the continental shelf. For the first time, this study investigated the response of benthic foraminifera to temperature change using microcosm culture and metabarcoding. Foraminiferal communities from three continental stations at varying water depths (6.0, 9.2, and 26.0 m) were cultured under five temperature gradients (6, 12, 18, 24, and 30°C), with each treatment performed in triplicate. The foraminifera were fed with microalgae every 4 days, and the filtered seawater (through 0.22 μm pores), acting as a medium, was changed accordingly. The experiment lasted for 80 days, and 47 DNA samples were obtained and analyzed, including three in situ samples. The results showed that foraminifera adjusted its growth rate within the low-temperature range and adopted an r-strategy to cope with high-temperature stress. In addition, the foraminifera from deeper water stations exhibited a pronounced vulnerability to diminishing read counts. The read counts, operational taxonomic units (OTU) counts and Margalef index of foraminifera and the read counts of Rotaliida exhibited a remarkably positive correlation with temperature. The recommended relationships were described as read counts = 1314.75*T + 44754.51; OTU counts = 1.13*T + 44.26; Margalef index =1.13*T + 44.26. This study established the first quantitative relationship between temperature and foraminifera molecular parameters that holds significant implications for long-time paleotemperature calibration in climate change.

The global genetic diversity of planktonic foraminifera reveals the structure of cryptic speciation in plankton

The nature and extent of diversity in the plankton has fascinated scientists for over a century. Initially, the discovery of many new species in the remarkably uniform and unstructured pelagic environment appeared to challenge the concept of ecological niches. Later, it became obvious that only a fraction of plankton diversity had been formally described, because plankton assemblages are dominated by understudied eukaryotic lineages with small size that lack clearly distinguishable morphological features. The high diversity of the plankton has been confirmed by comprehensive metabarcoding surveys, but interpretation of the underlying molecular taxonomies is hindered by insufficient integration of genetic diversity with morphological taxonomy and ecological observations. Here we use planktonic foraminifera as a study model and reveal the full extent of their genetic diversity and investigate geographical and ecological patterns in their distribution. To this end, we assembled a global data set of ~7600 ribosomal DNA sequences obtained from morphologically characterised individual foraminifera, established a robust molecular taxonomic framework for the observed diversity, and used it to query a global metabarcoding data set covering ~1700 samples with ~2.48 billion reads. This allowed us to extract and assign 1 million reads, enabling characterisation of the structure of the genetic diversity of the group across ~1100 oceanic stations worldwide. Our sampling revealed the existence of, at most, 94 distinct molecular operational taxonomic units (MOTUs) at a level of divergence indicative of biological species. The genetic diversity only doubles the number of formally described species identified by morphological features. Furthermore, we observed that the allocation of genetic diversity to morphospecies is uneven. Only 16 morphospecies disguise evolutionarily significant genetic diversity, and the proportion of morphospecies that show genetic diversity increases poleward. Finally, we observe that MOTUs have a narrower geographic distribution than morphospecies and that in some cases the MOTUs belonging to the same morphospecies (cryptic species) have different environmental preferences. Overall, our analysis reveals that even in the light of global genetic sampling, planktonic foraminifera diversity is modest and finite. However, the extent and structure of the cryptic diversity reveals that genetic diversification is decoupled from morphological diversification, hinting at different mechanisms acting at different levels of divergence.

Sediment DNA metabarcoding and morphology provide complementary insight into macrofauna and meiobenthos response to environmental gradients in an Arctic glacial fjord

Arctic fjords ecosystems are highly dynamic, with organisms exposed to various natural stressors along with productivity clines driven by advection of water masses from shelves. The benthic response to these environmental clines has been extensively studied using traditional, morphology-based approaches mostly focusing on macroinvertebrates. In this study we analyse the effects of glacially mediated disturbance on the biodiversity of benthic macrofauna and meiobenthos (meiofauna and Foraminifera) in a Svalbard fjord by comparing morphology and eDNA metabarcoding. Three genetic markers targeting metazoans (COI), meiofauna (18S V1V2) and Foraminifera (18S 37f) were analyzed. Univariate measures of alpha diversity and multivariate compositional dissimilarities were calculated and tested for similarities in response to environmental gradients using correlation analysis. Our study showed different taxonomic composition of morphological and molecular datasets for both macrofauna and meiobenthos. Some taxonomic groups while abundant in metabarcoding data were almost absent in morphology-based inventory and vice versa. In general, species richness and diversity measures in macrofauna morphological data were higher than in metabarcoding, and similar for the meiofauna. Both methodological approaches showed different patterns of response to the glacially mediated disturbance for the macrofauna and the meiobenthos. Macrofauna showed an evident distinction in taxonomic composition and a dramatic cline in alpha diversity indices between the outer and inner parts of fjord, while the meiobenthos showed a gradual change and more subtle responses to environmental changes along the fjord axis. The two methods can be seen as complementing rather than replacing each other. Morphological approach provides more accurate inventory of larger size species and more reliable quantitative data, while metabarcoding allows identification of inconspicuous taxa that are overlooked in morphology-based studies. As different taxa may show different sensitivities to environmental changes, both methods shall be used to monitor marine biodiversity in Arctic ecosystems and its response to dramatically changing environmental conditions.

Assigning the unassigned: A signature-based classification of rDNA metabarcodes reveals new deep-sea diversity

Environmental DNA metabarcoding reveals a vast genetic diversity of marine eukaryotes. Yet, most of the metabarcoding data remain unassigned due to the paucity of reference databases. This is particularly true for the deep-sea meiofauna and eukaryotic microbiota, whose hidden diversity is largely unexplored. Here, we tackle this issue by using unique DNA signatures to classify unknown metabarcodes assigned to deep-sea foraminifera. We analyzed metabarcoding data obtained from 311 deep-sea sediment samples collected in the Clarion-Clipperton Fracture Zone, an area of potential polymetallic nodule exploitation in the Eastern Pacific Ocean. Using the signatures designed in the 37F hypervariable region of the 18S rRNA gene, we were able to classify 802 unassigned metabarcodes into 61 novel lineages, which have been placed in 27 phylogenetic clades. The comparison of new lineages with other foraminiferal datasets shows that most novel lineages are widely distributed in the deep sea. Five lineages are also present in the shallow-water datasets; however, phylogenetic analysis of these lineages separates deep-sea and shallow-water metabarcodes except in one case. While the signature-based classification does not solve the problem of gaps in reference databases, this taxonomy-free approach provides insight into the distribution and ecology of deep-sea species represented by unassigned metabarcodes, which could be useful in future applications of metabarcoding for environmental monitoring.

Unlocking foraminiferal genetic diversity on estuarine mudflats with eDNA metabarcoding

Environmental biomonitoring is a prerequisite for efficient evaluation and remediation of ecosystem degradation due to anthropogenic pressure or climate change. Estuaries are key habitats subject to multiple anthropogenic and natural stressors. Due to these multiple stressors, the detection of anthropogenic pressure is challenging. The fact that abundant natural stressors often lead to negative quality assessments has been coined the "estuarine quality paradox". To solve this issue, the application of molecular approaches with successful bioindicators like foraminifera is promising. However, sampling protocols, molecular procedures and data analyses need to be validated before such tools can be routinely applied. We conducted an environmental DNA survey of estuarine mudflats along the French Atlantic coast, using a metabarcoding approach targeting foraminifera. Our results demonstrate that estuarine environments have only a few active OTUs dominating the community composition and a large stock of dormant or propagule stages. This last genetic diversity components constitute an important reservoir, with different species which can potentially develop in response to the temporal variability of the multiple stressors. In fact, different OTUs were dominant in the studied estuaries. Our statistical model shows that the physical and chemical characteristics of the sediment and the climatic conditions explain only 43 % of the community composition variance. This suggests that other, less easily quantifiable factors, such as the history and use of the estuaries or the ecological drift could play an important role as well. Environmental DNA biomonitoring opens new perspectives to better characterize the genetic diversity in estuaries.

Novel contributions to the molecular and morphological diversity of freshwater monothalamid foraminifera: Description of six new species

Non-marine monothalamous foraminifera are common in freshwater and soil habitats. They comprise a poorly-known group lacking sufficient information about diversity, morphology, reproduction, distribution and ecology. Based on an integrative morphological and molecular approach we describe a new family, a new genus and six new species of freshwater monothalamids from different localities in the Netherlands and France. We establish Astroperula as a new genus of organic walled freshwater foraminifera that contains two species, Astroperula dumacki and Astroperula parvipila. Furthermore, two new agglutinated freshwater monothalamids are described and illustrated, Limnogromia leanneae and Lacogromia pawlowskii, one new organic walled, Velamentofex dujardini, and a new naked monothalamid, Haplomyxa retiforma. Additional information is provided about a special form of cell division in Claparedellus lachmannii and L. leanneae, and feeding behavior in cultures of V. dujardini. Morphological observations are added for an unidentified Limnogromia sampled from the Netherlands and two types of Lacogromia sampled from an alpine region in the French Pyrenees and a karst sinkhole in Bosnia and Herzegovina. The present study provides additional insight into the ecology, diversity and occurrence of freshwater monothalamids and emphasizes the fact that a combination of morphological and molecular methods is necessary to clearly distinguish species in this group.

Metabarcoding reveals high diversity of benthic foraminifera linked to water masses circulation at coastal Svalbard

Arctic marine biodiversity is undergoing rapid changes due to global warming and modifications of oceanic water masses circulation. These changes have been demonstrated in the case of mega- and macrofauna, but much less is known about their impact on the biodiversity of smaller size organisms, such as foraminifera that represent a main component of meiofauna in the Arctic. Several studies analyzed the distribution and diversity of Arctic foraminifera. However, all these studies are based exclusively on the morphological identification of specimens sorted from sediment samples. Here, we present the first assessment of Arctic foraminifera diversity based on metabarcoding of sediment DNA samples collected in fjords and open sea areas in the Svalbard Archipelago. We obtained a total of 5,968,786 reads that represented 1384 amplicon sequence variants (ASVs). More than half of the ASVs (51.7%) could not be assigned to any group in the reference database suggesting a high genetic novelty of Svalbard foraminifera. The sieved and unsieved samples resolved comparable communities, sharing 1023 ASVs, comprising over 97% of reads. Our analyses show that the foraminiferal assemblage differs between the localities, with communities distinctly separated between fjord and open sea stations. Each locality was characterized by a specific assemblage, with only a small overlap in the case of open sea areas. Our study demonstrates a clear pattern of the influence of water masses on the structure of foraminiferal communities. The stations situated on the western coast of Svalbard that are strongly influenced by warm and salty Atlantic water (AW) are characterized by much higher diversity than stations in the northern and eastern part, where the impact of AW is less pronounced. This high diversity and specificity of Svalbard foraminifera associated with water mass distribution indicate that the foraminiferal metabarcoding data can be very useful for inferring present and past environmental conditions in the Arctic.

Renewal of planktonic foraminifera diversity after the Cretaceous Paleogene mass extinction by benthic colonizers

The biotic crisis following the end-Cretaceous asteroid impact resulted in a dramatic renewal of pelagic biodiversity. Considering the severe and immediate effect of the asteroid impact on the pelagic environment, it is remarkable that some of the most affected pelagic groups, like the planktonic foraminifera, survived at all. Here we queried a surface ocean metabarcoding dataset to show that calcareous benthic foraminifera of the clade Globothalamea are able to disperse actively in the plankton, and we show using molecular clock phylogeny that the modern planktonic clades originated from different benthic ancestors that colonized the plankton after the end-Cretaceous crisis. We conclude that the diversity of planktonic foraminifera has been the result of a constant leakage of benthic foraminifera diversity into the plankton, continuously refueling the planktonic niche, and challenge the classical interpretation of the fossil record that suggests that Mesozoic planktonic foraminifera gave rise to the modern communities.

Newly designed foraminifera primers identify habitat-specific lineages through metabarcoding analyses

Foraminifera include diverse shell-building lineages found in a wide array of aquatic habitats from the deep-sea to intertidal zones to brackish and freshwater ecosystems. Recent estimates of morphological and molecular foraminifera diversity have increased the knowledge of foraminiferal diversity, which is critical as these lineages are used as bioindicators of past and present environmental perturbation. However, a comparative analysis of foraminiferal biodiversity between their major habitats (freshwater, brackish, intertidal, and marine) is underexplored, particularly using molecular tools. Here, we present metabarcoding survey of foraminiferal diversity across different ecosystems using newly designed foraminifera-specific primers that target the hypervariable regions of the foraminifera SSU-rRNA gene (~250-300bp long). We tested these primer sets on four foraminifera species and then across several environments: the intertidal zone, coastal ecosystems, and freshwater vernal pools. We retrieved 655 operational taxonomic units (OTUs); the majority are undetermined taxa that have no closely-matching sequences in the database. Furthermore, we identified 163 OTUs with distinct habitat preferences. Most of the observed OTUs belonged to lineages of single-chambered foraminifera, including poorly explored freshwater foraminifera which encompass a clade of Reticulomyxa-like forms. Our pilot study provides the community with an additional set of newly designed and taxon-specific primers to elucidate foraminiferal diversity across different habitats.

Patterns of eukaryotic diversity from the surface to the deep-ocean sediment

Remote deep-ocean sediment (DOS) ecosystems are among the least explored biomes on Earth. Genomic assessments of their biodiversity have failed to separate indigenous benthic organisms from sinking plankton. Here, we compare global-scale eukaryotic DNA metabarcoding datasets (18S-V9) from abyssal and lower bathyal surficial sediments and euphotic and aphotic ocean pelagic layers to distinguish plankton from benthic diversity in sediment material. Based on 1685 samples collected throughout the world ocean, we show that DOS diversity is at least threefold that in pelagic realms, with nearly two-thirds represented by abundant yet unknown eukaryotes. These benthic communities are spatially structured by ocean basins and particulate organic carbon (POC) flux from the upper ocean. Plankton DNA reaching the DOS originates from abundant species, with maximal deposition at high latitudes. Its seafloor DNA signature predicts variations in POC export from the surface and reveals previously overlooked taxa that may drive the biological carbon pump.

Comparison of Deep-Sea Picoeukaryotic Composition Estimated from the V4 and V9 Regions of 18S rRNA Gene with a Focus on the Hadal Zone of the Mariana Trench

Diversity of microbial eukaryotes is estimated largely based on sequencing analysis of the hypervariable regions of 18S rRNA genes. But the use of different regions of 18S rRNA genes as molecular markers may generate bias in diversity estimation. Here, we compared the differences between the two most widely used markers, V4 and V9 regions of the 18S rRNA gene, in describing the diversity of epipelagic, bathypelagic, and hadal picoeukaryotes in the Challenger Deep of the Mariana Trench, which is a unique and little explored environment. Generally, the V9 region identified more OTUs in deeper waters than V4, while the V4 region provided greater Shannon diversity than V9. In the epipelagic zone, where Alveolata was the dominant group, picoeukaryotic community compositions identified by V4 and V9 markers are similar at different taxonomic levels. However, in the deep waters, the results of the two datasets show clear differences. These differences were mainly contributed by Retaria, Fungi, and Bicosoecida. The primer targeting the V9 region has an advantage in amplifying Bicosoecids in the bathypelagic and hadal zone of the Mariana Trench, and its high abundance in V9 dataset pointed out the possibility of Bicosoecids as a dominant group in this environment. Chrysophyceae, Fungi, MALV-I, and Retaria were identified as the dominant picoeukaryotes in the bathypelagic and hadal zone and potentially play important roles in deep-sea microbial food webs and biogeochemical cycling by their phagotrophic, saprotrophic, and parasitic life styles. Overall, the use of different markers of 18S rRNA gene allows a better assessment and understanding of the picoeukaryotic diversity in deep-sea environments.

High variability in SSU rDNA gene copy number among planktonic foraminifera revealed by single-cell qPCR

Metabarcoding has become the workhorse of community ecology. Sequencing a taxonomically informative DNA fragment from environmental samples gives fast access to community composition across taxonomic groups, but it relies on the assumption that the number of sequences for each taxon correlates with its abundance in the sampled community. However, gene copy number varies among and within taxa, and the extent of this variability must therefore be considered when interpreting community composition data derived from environmental sequencing. Here we measured with single-cell qPCR the SSU rDNA gene copy number of 139 specimens of five species of planktonic foraminifera. We found that the average gene copy number varied between of ~4000 to ~50,000 gene copies between species, and individuals of the same species can carry between ~300 to more than 350,000 gene copies. This variability cannot be explained by differences in cell size and considering all plausible sources of bias, we conclude that this variability likely reflects dynamic genomic processes acting during the life cycle. We used the observed variability to model its impact on metabarcoding and found that the application of a correcting factor at species level may correct the derived relative abundances, provided sufficiently large populations have been sampled.

Microbial Involvement in Carbon Transformation via CH4 and CO2 in Saline Sedimentary Pool

Simple Summary

Methane and carbon dioxide are commonly found in the environment and are considered the most important greenhouse gases. Transformation of these gases is in large carried by microorganisms, which occur even in extreme environments. This study presents methane-related biological processes in saline sediments of the Miocene Wieliczka Formation, Poland. Biological activity (carbon dioxide and methane production or methane oxidation), confirmed by stable isotope indices, occurred in all of the studied Wieliczka rocks. CH₄-utilizing microbes constituted 0.7–3.6% while methanogens (represented by Methanobacterium) only 0.01–0.5% of taxa present in the Wieliczka Salt Mine rocks. Water activity was the key factor regulating microbial activity in saline subsurface sediments. Generally, CO₂ respiration was higher in anaerobic conditions while methanogenic and methanotrophic activities were dependent on the type of rock.

Abstract

Methane and carbon dioxide are one of the most important greenhouse gases and significant components of the carbon cycle. Biogeochemical methane transformation may occur even in the extreme conditions of deep subsurface ecosystems. This study presents methane-related biological processes in saline sediments of the Miocene Wieliczka Formation, Poland. Rock samples (W2, W3, and W4) differed in lithology (clayey salt with veins of fibrous salt and lenses of gypsum and anhydrite; siltstone and sandstone; siltstone with veins of fibrous salt and lenses of anhydrite) and the accompanying salt type (spiza salts or green salt). Microbial communities present in the Miocene strata were studied using activity measurements and high throughput sequencing. Biological activity (i.e., carbon dioxide and methane production or methane oxidation) occurred in all of the studied clayey salt and siltstone samples but mainly under water-saturated conditions. Microcosm studies performed at elevated moisture created more convenient conditions for the activity of both methanogenic and methanotrophic microorganisms than the intact sediments. This points to the fact that water activity is an important factor regulating microbial activity in saline subsurface sediments. Generally, respiration was higher in anaerobic conditions and ranged from 36 ± 2 (W2_200%t.w.c) to 48 ± 4 (W3_200%t.w.c) nmol CO₂ gdw⁻¹ day⁻¹. Methanogenic activity was the highest in siltstone and sandstone (W3, 0.025 ± 0.018 nmol CH₄ gdw⁻¹ day⁻¹), while aerobic methanotrophic activity was the highest in siltstone with salt and anhydrite (W4, 220 ± 66 nmol CH₄ gdw⁻¹ day⁻¹). The relative abundance of CH₄-utilizing microorganisms (Methylomicrobium, Methylomonas, Methylocystis) constituted 0.7–3.6% of all taxa. Methanogens were represented by Methanobacterium (0.01–0.5%). The methane-related microbes were accompanied by a significant number of unclassified microorganisms (3–64%) and those of the Bacillus genus (4.5–91%). The stable isotope composition of the CO₂ and CH₄ trapped in the sediments suggests that methane oxidation could have influenced δ¹³C_CH4, especially in W3 and W4.

High and specific diversity of protists in the deep-sea basins dominated by diplonemids, kinetoplastids, ciliates and foraminiferans

Heterotrophic protists (unicellular eukaryotes) form a major link from bacteria and algae to higher trophic levels in the sunlit ocean. Their role on the deep seafloor, however, is only fragmentarily understood, despite their potential key function for global carbon cycling. Using the approach of combined DNA metabarcoding and cultivation-based surveys of 11 deep-sea regions, we show that protist communities, mostly overlooked in current deep-sea foodweb models, are highly specific, locally diverse and have little overlap to pelagic communities. Besides traditionally considered foraminiferans, tiny protists including diplonemids, kinetoplastids and ciliates were genetically highly diverse considerably exceeding the diversity of metazoans. Deep-sea protists, including many parasitic species, represent thus one of the most diverse biodiversity compartments of the Earth system, forming an essential link to metazoans.

Metadata standards and practical guidelines for specimen and DNA curation when building barcode reference libraries for aquatic life

DNA barcoding and metabarcoding is increasingly used to effectively and precisely assess and monitor biodiversity in aquatic ecosystems. As these methods rely on data availability and quality of barcode reference libraries, it is important to develop and follow best practices to ensure optimal quality and traceability of the metadata associated with the reference barcodes used for identification. Sufficient metadata, as well as vouchers, corresponding to each reference barcode must be available to ensure reliable barcode library curation and, thereby, provide trustworthy baselines for downstream molecular species identification. This document (1) specifies the data and metadata required to ensure the relevance, the accessibility and traceability of DNA barcodes and (2) specifies the recommendations for DNA harvesting and for the storage of both voucher specimens/samples and barcode data.

Looking for Hidden Enemies of Metabarcoding: Species Composition, Habitat and Management Can Strongly Influence DNA Extraction while Examining Grassland Communities

Despite the raising preoccupation, the critical question of how the plant community is composed belowground still remains unresolved, particularly for the conservation priority types of vegetation. The usefulness of metabarcoding analysis of the belowground parts of the plant community is subjected to a considerable bias, that often impedes detection of all species in a sample due to insufficient DNA quality or quantity. In the presented study we have attempted to find environmental factors that determine the amount and quality of DNA extracted from total plant tissue from above- and belowground samples (1000 and 10,000 cm2). We analyzed the influence of land use intensity, soil properties, species composition, and season on DNA extraction. The most important factors for DNA quality were vegetation type, soil conductometry (EC), and soil pH for the belowground samples. The species that significantly decreased the DNA quality were Calamagrostis epigejos, Coronilla varia, and Holcus lanatus. For the aboveground part of the vegetation, the season, management intensity, and certain species-with the most prominent being Centaurea rhenana and Cirsium canum-have the highest influence. Additionally, we found that sample size, soil granulation, MgO, organic C, K2O, and total soil N content are important for DNA extraction effectiveness. Both low EC and pH reduce significantly the yield and quality of DNA. Identifying the potential inhibitors of DNA isolation and predicting difficulties of sampling the vegetation plots for metabarcoding analysis will help to optimize the universal, low-cost multi-stage DNA extraction procedure in molecular ecology studies.

Integrating morphology and metagenomics to understand taxonomic variability of Amphisorus (Foraminifera, Miliolida) from Western Australia and Indonesia

Foraminifera are a group of mostly marine protists with high taxonomic diversity. Species identification is often complex, as both morphological and molecular approaches can be challenging due to a lack of unique characters and reference sequences. An integrative approach combining state of the art morphological and molecular tools is therefore promising. In this study, we analysed large benthic Foraminifera of the genus Amphisorus from Western Australia and Indonesia. Based on previous findings on high morphological variability observed in the Soritidae and the discontinuous distribution of Amphisorus along the coast of western Australia, we expected to find multiple morphologically and genetically unique Amphisorus types. In order to gain detailed insights into the diversity of Amphisorus, we applied micro CT scanning and shotgun metagenomic sequencing. We identified four distinct morphotypes of Amphisorus, two each in Australia and Indonesia, and showed that each morphotype is a distinct genotype. Furthermore, metagenomics revealed the presence of three dinoflagellate symbiont clades. The most common symbiont was Fugacium Fr5, and we could show that its genotypes were mostly specific to Amphisorus morphotypes. Finally, we assembled the microbial taxa associated with the two Western Australian morphotypes, and analysed their microbial community composition. Even though each Amphisorus morphotype harboured distinct bacterial communities, sampling location had a stronger influence on bacterial community composition, and we infer that the prokaryotic community is primarily shaped by the microhabitat rather than host identity. The integrated approach combining analyses of host morphology and genetics, dinoflagellate symbionts, and associated microbes leads to the conclusion that we identified distinct, yet undescribed taxa of Amphisorus. We argue that the combination of morphological and molecular methods provides unprecedented insights into the diversity of foraminifera, which paves the way for a deeper understanding of their biodiversity, and facilitates future taxonomic and ecological work.

Planktonic foraminifera eDNA signature deposited on the seafloor remains preserved after burial in marine sediments

Environmental DNA (eDNA) metabarcoding of marine sediments has revealed large amounts of sequences assigned to planktonic taxa. How this planktonic eDNA is delivered on the seafloor and preserved in the sediment is not well understood. We address these questions by comparing metabarcoding and microfossil foraminifera assemblages in sediment cores taken off Newfoundland across a strong ecological gradient. We detected planktonic foraminifera eDNA down to 30 cm and observed that the planktonic/benthic amplicon ratio changed with depth. The relative proportion of planktonic foraminiferal amplicons remained low from the surface down to 10 cm, likely due to the presence of DNA from living benthic foraminifera. Below 10 cm, the relative proportion of planktonic foraminifera amplicons rocketed, likely reflecting the higher proportion of planktonic eDNA in the DNA burial flux. In addition, the microfossil and metabarcoding assemblages showed a congruent pattern indicating that planktonic foraminifera eDNA is deposited without substantial lateral advection and preserves regional biogeographical patterns, indicating deposition by a similar mechanism as the foraminiferal shells. Our study shows that the planktonic eDNA preserved in marine sediments has the potential to record climatic and biotic changes in the pelagic community with the same spatial and temporal resolution as microfossils.

Benthic foraminiferal metabarcoding and morphology-based assessment around three offshore gas platforms: Congruence and complementarity

Since the 1960 s, there has been a rapid expansion of drilling activities in the central and northern Adriatic Sea to meet the increasing global energy demand. The discharges of organic and inorganic pollutants, as well as the alteration of the sediment substrate, are among the main impacts associated with these activities. In the present study, we evaluate the response of benthic foraminifera to the activities of three gas platforms in the northwestern Adriatic Sea, with a special focus on the Armida A platform for which extensive geochemical data (organic matter, trace elements, polycyclic aromatic hydrocarbons, other hydrocarbons, and volatile organic compounds) are available. The response to disturbance is assessed by analyzing the foraminiferal diversity using the traditional morphology-based approach and by 18S rDNA-based metabarcoding. The two methods give congruent results, showing relatively lower foraminiferal diversity and higher dominance values at stations closer to the platforms (<50 m). The taxonomic compositions of the morphological and metabarcoding datasets are very different, the latter being dominated by monothalamous, mainly soft-walled species. However, compositional changes consistently occur at 50 m from the platform and can be related to variations in sediment grain-size variation and higher concentrations of Ni, Zn, Ba, hydrocarbons and total organic carbon. Additionally, several morphospecies and Molecular Operational Taxonomic Units (MOTUs) show strong correlations with distance from the platform and with environmental parameters extracted from BIOENV analysis. Some of these MOTUs have the potential to become new bioindicators, complementing the assemblage of hard-shelled foraminiferal species detected through microscopic analyses. The congruence and complementarity between metabarcoding and morphological approaches support the application of foraminiferal metabarcoding in routine biomonitoring surveys as a reliable, time- and cost-effective methodology to assess the environmental impacts of marine industries.

Amoeboid protist systematics: A report on the "Systematics of amoeboid protists" symposium at the VIIIth ECOP/ISOP meeting in Rome, 2019

Amoeboid protists are extremely abundant and diverse in natural systems where they often play outstanding ecological roles. They can be found in almost all major eukaryotic divisions, and genomic approaches are bringing major changes in our perception of their deep evolutionary relationships. At fine taxonomic levels, the generalization of barcoding is revealing a considerable and unsuspected specific diversity that can be appreciated with careful morphometric analyses based on light and electron microscopic observations. We provide examples on the difficulties and advances in amoeboid protists systematics in a selection of groups that were presented at the VIIIth ECOP/ISOP meeting in Rome, 2019. We conclude that, in all studied groups, important taxonomical rearrangements will certainly take place in the next few years, and systematics must be adapted to incorporate these changes. Notably, nomenclature should be flexible enough to integrate many new high level taxa, and a unified policy must be adopted to species description and to the establishment of types.

Diversity hotspot and unique community structure of foraminifera in the world's deepest marine blue hole - Sansha Yongle Blue Hole

Marine blue holes are precious geological heritages with high scientific research values. Their physical and chemical characteristics are unique because of the steep-walled structure and isolated water column which create isolated ecosystems in geographically restricted areas. The Sansha Yongle Blue Hole (SYBH) is the world's deepest marine blue hole. Here, we generated the first DNA metabarcoding dataset from SYBH sediment focusing on foraminifera, a group of protists that have colonized various marine environments. We collected sediment samples from SYBH along a depth gradient to characterize the foraminiferal diversity and compared them with the foraminiferal diversity of the costal Jiaozhou Bay (JZB) and the abyssal Northwest Pacific Ocean (NWP). We amplified the SSU rDNA of foraminifera and sequenced them with high-throughput sequencing. The results showed that the foraminiferal assemblages in SYBH were vertically structured in response to the abiotic gradients and diversity was higher than in JZB and NWP. This study illustrates the capacity of foraminifera to colonize hostile environments and shows that blue holes are natural laboratories to explore physiological innovation associated with anoxia.

Protist diversity and function in the dark ocean - Challenging the paradigms of deep-sea ecology with special emphasis on foraminiferans and naked protists

The dark ocean and the underlying deep seafloor together represent the largest environment on this planet, comprising about 80% of the oceanic volume and covering more than two-thirds of the Earth's surface, as well as hosting a major part of the total biosphere. Emerging evidence suggests that these vast pelagic and benthic habitats play a major role in ocean biogeochemistry and represent an "untapped reservoir" of high genetic and metabolic microbial diversity. Due to its huge volume, the water column of the dark ocean is the largest reservoir of organic carbon in the biosphere and likely plays a major role in the global carbon budget. The dark ocean and the seafloor beneath it are also home to a largely enigmatic food web comprising little-known and sometimes spectacular organisms, mainly prokaryotes and protists. This review considers the globally important role of pelagic and benthic protists across all protistan size classes in the deep-sea realm, with a focus on their taxonomy, diversity, and physiological properties, including their role in deep microbial food webs. We argue that, given the important contribution that protists must make to deep-sea biodiversity and ecosystem processes, they should not be overlooked in biological studies of the deep ocean.

Modern Benthic Foraminiferal Diversity: An Initial Insight into the Total Foraminiferal Diversity along the Kuwait Coastal Water

Kuwait territorial water hosts an important part of national biodiversity (i.e., zooplankton and phytoplankton), but very limited information exists on the overall diversity of benthic foraminifera. On the basis of the integration of publications, reports and theses with new available data from the Kuwait Bay and the northern islands, this study infers the total benthic foraminiferal diversity within Kuwait territorial water. This new literature survey documents the presence of 451 species belonging to 156 genera, 64 families, 31 superfamilies and 9 orders. These values are relatively high in consideration of the limited extension and the shallow depth of the Kuwait territorial water. Kuwait waters offer a variety of different environments and sub-environments (low salinity/muddy areas in the northern part, embayment, rocky tidal flats, coral reef systems, islands and shelf slope) that all together host largely diversified benthic foraminiferal communities. These figures are herein considered as underestimated because of the grouping of unassigned species due to the lack of reference collections and materials, as well as the neglection of the soft-shell monothalamids (‘allogromiids’).

Metabarcoding Insights Into the Trophic Behavior and Identity of Intertidal Benthic Foraminifera

Foraminifera are ubiquitous marine protists with an important role in the benthic carbon cycle. However, morphological observations often fail to resolve their exact taxonomic placement and there is a lack of field studies on their particular trophic preferences. Here, we propose the application of metabarcoding as a tool for the elucidation of the in situ feeding behavior of benthic foraminifera, while also allowing the correct taxonomic assignment of the feeder, using the V9 region of the 18S (small subunit; SSU) rRNA gene. Living foraminiferal specimens were collected from two intertidal mudflats of the Wadden Sea and DNA was extracted from foraminiferal individuals and from the surrounding sediments. Molecular analysis allowed us to confirm that our foraminiferal specimens belong to three genetic types: Ammonia sp. T6, Elphidium sp. S5 and Haynesina sp. S16. Foraminiferal intracellular eukaryote communities reflected to an extent those of the surrounding sediments but at different relative abundances. Unlike sediment eukaryote communities, which were largely determined by the sampling site, foraminiferal intracellular eukaryote communities were driven by foraminiferal species, followed by sediment depth. Our data suggests that Ammonia sp. T6 can predate on metazoan classes, whereas Elphidium sp. S5 and Haynesina sp. S16 are more likely to ingest diatoms. These observations, alongside the use of metabarcoding in similar ecological studies, significantly contribute to our overall understanding of the ecological roles of these protists in intertidal benthic environments and their position and function in the benthic food webs.

Sample size effects on the assessment of eukaryotic diversity and community structure in aquatic sediments using high-throughput sequencing

Understanding how biodiversity changes in time and space is vital to assess the effects of environmental change on benthic ecosystems. Due to the limitations of morphological methods, there has been a rapid expansion in the application of high-throughput sequencing methods to study benthic eukaryotic communities. However, the effect of sample size and small-scale spatial variation on the assessment of benthic eukaryotic diversity is still not well understood. Here, we investigate the effect of different sample volumes in the genetic assessment of benthic metazoan and non-metazoan eukaryotic community composition. Accordingly, DNA was extracted from five different cumulative sediment volumes comprising 100% of the top 2 cm of five benthic sampling cores, and used as template for Ilumina MiSeq sequencing of 18 S rRNA amplicons. Sample volumes strongly impacted diversity metrics for both metazoans and non-metazoan eukaryotes. Beta-diversity of treatments using smaller sample volumes was significantly different from the beta-diversity of the 100% sampled area. Overall our findings indicate that sample volumes of 0.2 g (1% of the sampled area) are insufficient to account for spatial heterogeneity at small spatial scales, and that relatively large percentages of sediment core samples are needed for obtaining robust diversity measurement of both metazoan and non-metazoan eukaryotes.

Eyes of the Deep-sea Floor: The Integrative Taxonomy of the Foraminiferal Genus Vanhoeffenella

Vanhoeffenella is a common deep-sea monothalamous foraminifer, some species of which have a unique eye-like test morphology. Owing to its world-wide distribution, it has been recorded numerous times since the "heroic age" of the deep-sea exploration in the early 20th century. So far, only 4 species have been described, and no attempts have been made to estimate the real diversity of this peculiar genus. Over the last fifteen years, we have collected specimens of Vanhoeffenella from various deep-sea areas, providing the basis for an integrative taxonomy and biogeography of this genus. Here, we clarify the phylogenetic position of Vanhoeffenella and give an account of its diversity in the Atlantic, Arctic and Southern Oceans (the Weddell Sea) as revealed by genetic marker (SSU rDNA) and morphology. Our study shows that Vanhoeffenella branches within Clade F of monothalamids and incorporates at least 10 putative species. Some could be distinguished by either morphological or molecular features, but only the integrative taxonomic approach provides a robust way to assess their diversity. We examine the new material of the type species (V. gaussi Rhumbler), redescribe the poorly-known V. oculus Earland and describe formally a fifth species, V. dilatata sp. nov.

Censusing marine eukaryotic diversity in the twenty-first century

The ocean constitutes one of the vastest and richest biomes on our planet. Most recent estimations, all based on indirect approaches, suggest that there are millions of marine eukaryotic species. Moreover, a large majority of these are small (less than 1 mm), cryptic and still unknown to science. However, this knowledge gap, caused by the lack of diagnostic morphological features in small organisms and the limited sampling of the global ocean, is currently being filled, thanks to new DNA-based approaches. The molecular technique of PCR amplification of homologous gene regions combined with high-throughput sequencing, routinely used to census unculturable prokaryotes, is now also being used to characterize whole communities of marine eukaryotes. Here, we review how this methodological advancement has helped to better quantify the magnitude and patterns of marine eukaryotic diversity, with an emphasis on taxonomic groups previously largely overlooked. We then discuss obstacles remaining to achieve a global understanding of marine eukaryotic diversity. In particular, we argue that 18S variable regions do not provide sufficient taxonomic resolution to census marine life, and suggest combining broad eukaryotic surveys targeting the 18S rRNA region with more taxon-focused analyses of hypervariable regions to improve our understanding of the diversity of species, the functional units of marine ecosystems.

This article is part of the themed issue ‘From DNA barcodes to biomes’.

Assessing SSU rDNA Barcodes in Foraminifera: A Case Study using Bolivina quadrata

The Small Subunit Ribosomal RNA gene (SSU rDNA) is a widely used tool to reconstruct phylogenetic relationships among foraminiferal species. Recently, the highly variable regions of this gene have been proposed as DNA barcodes to identify foraminiferal species. However, the resolution of these barcodes has not been well established, yet. In this study, we evaluate four SSU rDNA hypervariable regions (37/f, 41/f, 43/e, and 45/e) as DNA barcodes to distinguish among species of the genus Bolivina, with particular emphasis on Bolivina quadrata for which ten new sequences (KY468817-KY468826) were obtained during this study. Our analyses show that a single SSU rDNA hypervariable sequence is insufficient to resolve all Bolivina species and that some regions (37/f and 41/f) are more useful than others (43/e and 45/e) to distinguish among closely related species. In addition, polymorphism analyses reveal a high degree of variability. In the context of barcoding studies, these results emphasize the need to assess the range of intraspecific variability of DNA barcodes prior to their application to identify foraminiferal species in environmental samples; our results also highlight the possibility that a longer SSU rDNA region might be required to distinguish among species belonging to the same taxonomic group (i.e. genus).

Preferential suppression of Anopheles gambiae host sequences allows detection of the mosquito eukaryotic microbiome

Anopheles mosquitoes are vectors of the human malaria parasite, Plasmodium falciparum. The vector microbiota is a likely factor influencing parasite transmission. The prokaryotic microbiota of mosquitoes is efficiently surveyed by sequencing of hypervariable regions of the 16s ribosomal RNA (rRNA) gene. However, identification of the eukaryotic microbiota by targeting the 18s rRNA gene is challenging due to simultaneous amplification of the abundant 18s rRNA gene target in the mosquito host. Consequently, the eukaryotic microbial diversity of mosquitoes is vastly underexplored. An efficient methodology is needed to identify this component of the microbiota, expected to include relatives of Plasmodium. Here, we use defined panels of Anopheles samples from West Africa to test two experimental PCR clamp approaches to maximize the specific amplification of 18s rRNA gene hypervariable regions from eukaryotic microbes: anneal-inhibiting blocking primers and peptide-nucleic acid (PNA) oligonucleotide blockers. Of the two, PNA blockers were the only efficient blocking strategy, allowing a reduction of mosquito 18s rRNA gene sequences by more than 80% for the V4 hypervariable region. These PNA blockers will facilitate taxonomic profiling of the eukaryotic microbiota of the A. gambiae species complex, and contribute to a better understanding of microbial influence upon immunity and pathogen infection.

Novel benthic foraminifera are abundant and diverse in an area of the abyssal equatorial Pacific licensed for polymetallic nodule exploration

The benthic biota of the Clarion–Clipperton Zone (CCZ, abyssal eastern equatorial Pacific) is the focus of a major research effort linked to possible future mining of polymetallic nodules. Within the framework of ABYSSLINE, a biological baseline study conducted on behalf of Seabed Resources Development Ltd. in the UK-1 exploration contract area (eastern CCZ, ~4,080 m water depth), we analysed foraminifera (testate protists), including ‘live’ (Rose Bengal stained) and dead tests, in 5 cores (0–1 cm layer, >150-μm fraction) recovered during separate megacorer deployments inside a 30 by 30 km seafloor area. In both categories (live and dead) we distinguished between complete and fragmented specimens. The outstanding feature of these assemblages is the overwhelming predominance of monothalamids, a group often ignored in foraminiferal studies. These single-chambered foraminifera, which include agglutinated tubes, spheres and komokiaceans, represented 79% of 3,607 complete tests, 98% of 1,798 fragments and 76% of the 416 morphospecies (live and dead combined) in our samples. Only 3.1% of monothalamid species and 9.8% of all species in the UK-1 assemblages are scientifically described and many are rare (29% singletons). Our results emphasise how little is known about foraminifera in abyssal areas that may experience major impacts from future mining activities.

Spatio-temporal monitoring of deep-sea communities using metabarcoding of sediment DNA and RNA

We assessed spatio-temporal patterns of diversity in deep-sea sediment communities using metabarcoding. We chose a recently developed eukaryotic marker based on the v7 region of the 18S rRNA gene. Our study was performed in a submarine canyon and its adjacent slope in the Northwestern Mediterranean Sea, sampled along a depth gradient at two different seasons. We found a total of 5,569 molecular operational taxonomic units (MOTUs), dominated by Metazoa, Alveolata and Rhizaria. Among metazoans, Nematoda, Arthropoda and Annelida were the most diverse. We found a marked heterogeneity at all scales, with important differences between layers of sediment and significant changes in community composition with zone (canyon vs slope), depth, and season. We compared the information obtained from metabarcoding DNA and RNA and found more total MOTUs and more MOTUs per sample with DNA (ca. 20% and 40% increase, respectively). Both datasets showed overall similar spatial trends, but most groups had higher MOTU richness with the DNA template, while others, such as nematodes, were more diverse in the RNA dataset. We provide metabarcoding protocols and guidelines for biomonitoring of these key communities in order to generate information applicable to management efforts.

First evaluation of foraminiferal metabarcoding for monitoring environmental impact from an offshore oil drilling site

At present, environmental impacts from offshore oil and gas activities are partly determined by measuring changes in macrofauna diversity. Morphological identification of macrofauna is time-consuming, expensive and dependent on taxonomic expertise. In this study, we evaluated the applicability of using foraminiferal-specific metabarcoding for routine monitoring. Sediment samples were collected along distance gradients from two oil platforms off Taranaki (New Zealand) and their physico-chemical properties, foraminiferal environmental DNA/RNA, and macrofaunal composition analyzed. Macrofaunal and foraminiferal assemblages showed similar shifts along impact gradients, but responded differently to environmental perturbations. Macrofauna were affected by hypoxia, whereas sediment grain size appeared to drive shifts in foraminifera. We identified eight foraminiferal molecular operational taxonomic units that have potential to be used as bioindicator taxa. Our results show that metabarcoding represents an effective tool for assessing foraminiferal communities near offshore oil and gas platforms, and that it can be used to complement current monitoring techniques.

Evaluation of next generation sequencing for the analysis of Eimeria communities in wildlife

Next-generation sequencing (NGS) techniques are well-established for studying bacterial communities but not yet for microbial eukaryotes. Parasite communities remain poorly studied, due in part to the lack of reliable and accessible molecular methods to analyse eukaryotic communities. We aimed to develop and evaluate a methodology to analyse communities of the protozoan parasite Eimeria from populations of the Australian marsupial Petrogale penicillata (brush-tailed rock-wallaby) using NGS. An oocyst purification method for small sample sizes and polymerase chain reaction (PCR) protocol for the 18S rRNA locus targeting Eimeria was developed and optimised prior to sequencing on the Illumina MiSeq platform. A data analysis approach was developed by modifying methods from bacterial metagenomics and utilising existing Eimeria sequences in GenBank. Operational taxonomic unit (OTU) assignment at a high similarity threshold (97%) was more accurate at assigning Eimeria contigs into Eimeria OTUs but at a lower threshold (95%) there was greater resolution between OTU consensus sequences. The assessment of two amplification PCR methods prior to Illumina MiSeq, single and nested PCR, determined that single PCR was more sensitive to Eimeria as more Eimeria OTUs were detected in single amplicons. We have developed a simple and cost-effective approach to a data analysis pipeline for community analysis of eukaryotic organisms using Eimeria communities as a model. The pipeline provides a basis for evaluation using other eukaryotic organisms and potential for diverse community analysis studies.

Protistan community analysis: key findings of a large-scale molecular sampling

Protists are perhaps the most lineage-rich of microbial lifeforms, but remain largely unknown. High-throughput sequencing technologies provide opportunities to screen whole habitats in depth and enable detailed comparisons of different habitats to measure, compare and map protistan diversity. Such comparisons are often limited by low sample numbers within single studies and a lack of standardisation between studies. Here, we analysed 232 samples from 10 sampling campaigns using a standardised PCR protocol and bioinformatics pipeline. We show that protistan community patterns are highly consistent within habitat types and geographic regions, provided that sample processing is standardised. Community profiles are only weakly affected by fluctuations of the abundances of the most abundant taxa and, therefore, provide a sound basis for habitat comparison beyond random short-term fluctuations in the community composition. Further, we provide evidence that distribution patterns are not solely resulting from random processes. Distinct habitat types and distinct taxonomic groups are dominated by taxa with distinct distribution patterns that reflect their ecology with respect to dispersal and habitat colonisation. However, there is no systematic shift of the distribution pattern with taxon abundance.

Patterns and processes in microbial biogeography: do molecules and morphologies give the same answers?

Our knowledge on microbial biogeography depends on the way we define and study diversity. In contrast to most microbes, some protist lineages have conspicuous structures that allow comparisons of diversity concepts and measures-those based on molecules and those based on morphology. We analyzed a group of shell-bearing planktonic ciliates, the tintinnids, in a coast-to-ocean gradient using high-throughput sequencing and microscopy. First, we compared molecular operational taxonomic units (OTUs) and morphospecies in terms of assemblage composition, distribution and relationships with the environment. OTUs revealed potentially novel and rare taxa, while morphospecies showed clearer correlations with environmental factors, and both approaches coincided in supporting a coastal versus oceanic pattern. Second, we explored which processes influence assembly across the environmental gradient examined. Assemblage fluctuations were associated with significant distance-decay and changes in morphospecies size and prey proxies, thus suggesting niche partitioning as a key structuring mechanism. Our conclusion is that molecules and morphologies generally agreed, but they provided complementary data, the first revealing hidden diversity, and the latter making better connections between distribution patterns and ecological processes. This highlights the importance of linking genotypes and phenotypes (using multidisciplinary analyses and/or reliable databases of barcoded species), to understand the diversity, biogeography and ecological roles of microbes.

Microbial eukaryotic distributions and diversity patterns in a deep-sea methane seep ecosystem

Although chemosynthetic ecosystems are known to support diverse assemblages of microorganisms, the ecological and environmental factors that structure microbial eukaryotes (heterotrophic protists and fungi) are poorly characterized. In this study, we examined the geographic, geochemical and ecological factors that influence microbial eukaryotic composition and distribution patterns within Hydrate Ridge, a methane seep ecosystem off the coast of Oregon using a combination of high-throughput 18S rRNA tag sequencing, terminal restriction fragment length polymorphism fingerprinting, and cloning and sequencing of full-length 18S rRNA genes. Microbial eukaryotic composition and diversity varied as a function of substrate (carbonate versus sediment), activity (low activity versus active seep sites), sulfide concentration, and region (North versus South Hydrate Ridge). Sulfide concentration was correlated with changes in microbial eukaryotic composition and richness. This work also revealed the influence of oxygen content in the overlying water column and water depth on microbial eukaryotic composition and diversity, and identified distinct patterns from those previously observed for bacteria, archaea and macrofauna in methane seep ecosystems. Characterizing the structure of microbial eukaryotic communities in response to environmental variability is a key step towards understanding if and how microbial eukaryotes influence seep ecosystem structure and function.

Unveiling the Biodiversity of Deep-Sea Nematodes through Metabarcoding: Are We Ready to Bypass the Classical Taxonomy?

Nematodes inhabiting benthic deep-sea ecosystems account for >90% of the total metazoan abundances and they have been hypothesised to be hyper-diverse, but their biodiversity is still largely unknown. Metabarcoding could facilitate the census of biodiversity, especially for those tiny metazoans for which morphological identification is difficult. We compared, for the first time, different DNA extraction procedures based on the use of two commercial kits and a previously published laboratory protocol and tested their suitability for sequencing analyses of 18S rDNA of marine nematodes. We also investigated the reliability of Roche 454 sequencing analyses for assessing the biodiversity of deep-sea nematode assemblages previously morphologically identified. Finally, intra-genomic variation in 18S rRNA gene repeats was investigated by Illumina MiSeq in different deep-sea nematode morphospecies to assess the influence of polymorphisms on nematode biodiversity estimates. Our results indicate that the two commercial kits should be preferred for the molecular analysis of biodiversity of deep-sea nematodes since they consistently provide amplifiable DNA suitable for sequencing. We report that the morphological identification of deep-sea nematodes matches the results obtained by metabarcoding analysis only at the order-family level and that a large portion of Operational Clustered Taxonomic Units (OCTUs) was not assigned. We also show that independently from the cut-off criteria and bioinformatic pipelines used, the number of OCTUs largely exceeds the number of individuals and that 18S rRNA gene of different morpho-species of nematodes displayed intra-genomic polymorphisms. Our results indicate that metabarcoding is an important tool to explore the diversity of deep-sea nematodes, but still fails in identifying most of the species due to limited number of sequences deposited in the public databases, and in providing quantitative data on the species encountered. These aspects should be carefully taken into account before using metabarcoding in quantitative ecological research and monitoring programmes of marine biodiversity.

Accurate assessment of the impact of salmon farming on benthic sediment enrichment using foraminiferal metabarcoding

Assessing the environmental impact of salmon farms on benthic systems is traditionally undertaken using biotic indices derived from microscopic analyses of macrobenthic infaunal (MI) communities. In this study, we tested the applicability of using foraminiferal-specific high-throughput sequencing (HTS) metabarcoding for monitoring these habitats. Sediment samples and physico-chemical data were collected along an enrichment gradient radiating out from three Chinook salmon (Oncorhynchus tshawytscha) farms in New Zealand. HTS of environmental DNA and RNA (eDNA/eRNA) resulted in 1,875,300 sequences that clustered into 349 Operational Taxonomic Units. Strong correlations were observed among various biotic indices calculated from MI data and normalized fourth-root transformed HTS data. Correlations were stronger using eRNA compared to eDNA data. Quantile regression spline analyses identified 12 key foraminiferal taxa that have potential to be used as bioindicator species. This study demonstrates the huge potential for using this method for biomonitoring of fish-farming and other marine industrial activities.

Highly divergent ancient gene families in metagenomic samples are compatible with additional divisions of life

Background

Microbial genetic diversity is often investigated via the comparison of relatively similar 16S molecules through multiple alignments between reference sequences and novel environmental samples using phylogenetic trees, direct BLAST matches, or phylotypes counts. However, are we missing novel lineages in the microbial dark universe by relying on standard phylogenetic and BLAST methods? If so, how can we probe that universe using alternative approaches? We performed a novel type of multi-marker analysis of genetic diversity exploiting the topology of inclusive sequence similarity networks.

Results

Our protocol identified 86 ancient gene families, well distributed and rarely transferred across the 3 domains of life, and retrieved their environmental homologs among 10 million predicted ORFs from human gut samples and other metagenomic projects. Numerous highly divergent environmental homologs were observed in gut samples, although the most divergent genes were over-represented in non-gut environments. In our networks, most divergent environmental genes grouped exclusively with uncultured relatives, in maximal cliques. Sequences within these groups were under strong purifying selection and presented a range of genetic variation comparable to that of a prokaryotic domain.

Conclusions

Many genes families included environmental homologs that were highly divergent from cultured homologs: in 79 gene families (including 18 ribosomal proteins), Bacteria and Archaea were less divergent than some groups of environmental sequences were to any cultured or viral homologs. Moreover, some groups of environmental homologs branched very deeply in phylogenetic trees of life, when they were not too divergent to be aligned. These results underline how limited our understanding of the most diverse elements of the microbial world remains, and encourage a deeper exploration of natural communities and their genetic resources, hinting at the possibility that still unknown yet major divisions of life have yet to be discovered.

Reviewers

This article was reviewed by Eugene Koonin, William Martin and James McInerney.

Electronic supplementary material

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

High-Throughput Sequencing-The Key to Rapid Biodiversity Assessment of Marine Metazoa?

The applications of traditional morphological and molecular methods for species identification are greatly restricted by processing speed and on a regional or greater scale are generally considered unfeasible. In this context, high-throughput sequencing, or metagenetics, has been proposed as an efficient tool to document biodiversity. Here we evaluated the effectiveness of 454 pyrosequencing in marine metazoan community analysis using the 18S rDNA: V1-V2 region. Multiplex pyrosequencing of the V1-V2 region was used to analyze two pooled samples of DNA, one comprising 118 and the other 37 morphologically identified species, and one natural sample taken directly from a North Sea zooplankton community. A DNA reference library comprising all species represented in the pooled samples was created by Sanger sequencing, and this was then used to determine the optimal similarity threshold for species delineation. The optimal threshold was found at 99% species similarity, with 85% identification success. Pyrosequencing was able to identify between fewer species: 67% and 78% of the species in the two pooled samples. Also, a large number of sequences for three species that were not included in the pooled samples were amplified by pyrosequencing, suggesting preferential amplification of some genotypes and the sensitivity of this approach to even low levels of contamination. Conversely, metagenetic analysis of the natural zooplankton sample identified many more species (particularly gelatinous zooplankton and meroplankton) than morphological analysis of a formalin-fixed sample from the same sampling site, suggesting an increased level of taxonomic resolution with pyrosequencing. The study demonstrated that, based on the V1-V2 region, 454 sequencing does not provide accurate species differentiation and reliable taxonomic classification, as it is required in most biodiversity monitoring. The analysis of artificially prepared samples indicated that species detection in pyrosequencing datasets is complicated by potential PCR-based biases and that the V1-V2 marker is poorly resolved for some taxa.

Deep-Sea, Deep-Sequencing: Metabarcoding Extracellular DNA from Sediments of Marine Canyons

Marine sediments are home to one of the richest species pools on Earth, but logistics and a dearth of taxonomic work-force hinders the knowledge of their biodiversity. We characterized α- and β-diversity of deep-sea assemblages from submarine canyons in the western Mediterranean using an environmental DNA metabarcoding. We used a new primer set targeting a short eukaryotic 18S sequence (ca. 110 bp). We applied a protocol designed to obtain extractions enriched in extracellular DNA from replicated sediment corers. With this strategy we captured information from DNA (local or deposited from the water column) that persists adsorbed to inorganic particles and buffered short-term spatial and temporal heterogeneity. We analysed replicated samples from 20 localities including 2 deep-sea canyons, 1 shallower canal, and two open slopes (depth range 100–2,250 m). We identified 1,629 MOTUs, among which the dominant groups were Metazoa (with representatives of 19 phyla), Alveolata, Stramenopiles, and Rhizaria. There was a marked small-scale heterogeneity as shown by differences in replicates within corers and within localities. The spatial variability between canyons was significant, as was the depth component in one of the canyons where it was tested. Likewise, the composition of the first layer (1 cm) of sediment was significantly different from deeper layers. We found that qualitative (presence-absence) and quantitative (relative number of reads) data showed consistent trends of differentiation between samples and geographic areas. The subset of exclusively benthic MOTUs showed similar patterns of β-diversity and community structure as the whole dataset. Separate analyses of the main metazoan phyla (in number of MOTUs) showed some differences in distribution attributable to different lifestyles. Our results highlight the differentiation that can be found even between geographically close assemblages, and sets the ground for future monitoring and conservation efforts on these bottoms of ecological and economic importance.

Marine Nematode Taxonomy in Africa: Promising Prospects Against Scarcity of Information

From the late 19th century, Africa has faced heavy exploitation of its natural resources with increasing land/water pollution, and several described species have already become extinct or close to extinction. This could also be the case for marine nematodes, which are the most abundant and diverse benthic group in marine sediments, and play major roles in ecosystem functioning. Compared to Europe and North America, only a handful of investigations on marine nematodes have been conducted to date in Africa. This is due to the scarcity of experienced taxonomists, absence of identification guides, as well as local appropriate infrastructures. A pivotal project has started recently between nematologists from Africa (Tunisia), India, and Europe (Italy) to promote taxonomic study and biodiversity estimation of marine nematodes in the African continent. To do this, as a first step, collection of permanent slides of marine nematodes (235 nominal species and 14 new to science but not yet described) was recently established at the Faculty of Sciences of Bizerte (Tunisia). Capacity building of next generation of African taxonomists have been carried out at level of both traditional and molecular taxonomy (DNA barcoding and next-generation sequencing [NGS]), but they need to be implemented. Indeed, the integration of these two approaches appears crucial to overcome lack of information on the taxonomy, ecology, and biodiversity of marine nematodes from African coastal waters.

Insights into Diversity and Imputed Metabolic Potential of Bacterial Communities in the Continental Shelf of Agatti Island

Marine microbes play a key role and contribute largely to the global biogeochemical cycles. This study aims to explore microbial diversity from one such ecological hotspot, the continental shelf of Agatti Island. Sediment samples from various depths of the continental shelf were analyzed for bacterial diversity using deep sequencing technology along with the culturable approach. Additionally, imputed metagenomic approach was carried out to understand the functional aspects of microbial community especially for microbial genes important in nutrient uptake, survival and biogeochemical cycling in the marine environment. Using culturable approach, 28 bacterial strains representing 9 genera were isolated from various depths of continental shelf. The microbial community structure throughout the samples was dominated by phylum Proteobacteria and harbored various bacterioplanktons as well. Significant differences were observed in bacterial diversity within a short region of the continental shelf (1-40 meters) i.e. between upper continental shelf samples (UCS) with lesser depths (i.e. 1-20 meters) and lower continental shelf samples (LCS) with greater depths (i.e. 25-40 meters). By using imputed metagenomic approach, this study also discusses several adaptive mechanisms which enable microbes to survive in nutritionally deprived conditions, and also help to understand the influence of nutrition availability on bacterial diversity.