Environmental DNA recovers fish composition turnover of the coral reefs of West Indian Ocean islands.

Islands have been used as model systems to study ecological and evolutionary processes, and they provide an ideal set-up for validating new biodiversity monitoring methods. The application of environmental DNA metabarcoding for monitoring marine biodiversity requires an understanding of the spatial scale of the eDNA signal, which is best tested in island systems. Here, we investigated the variation in Actinopterygii and Elasmobranchii species composition recovered from eDNA metabarcoding along a gradient of distance-to-reef in four of the five French Scattered Islands in the Western Indian Ocean. We collected surface water samples at an increasing distance from reefs (0 m, 250 m, 500 m, 750 m). We used a metabarcoding protocol based on the 'teleo' primers to target marine reef fishes and classified taxa according to their habitat types (benthic or pelagic). We investigated the effect of distance-to-reef on ß diversity variation using generalised linear mixed models and estimated species-specific distance-to-reef effects using a model-based approach for community data. Environmental DNA metabarcoding analyses recovered distinct fish species compositions across the four inventoried islands and variations along the distance-to-reef gradient. The analysis of ß-diversity variation showed significant taxa turnover between the eDNA samples on and away from the reefs. In agreement with a spatially localised signal from eDNA, benthic species were distributed closer to the reef than pelagic ones. Our findings demonstrate that the combination of eDNA inventories and spatial modelling can provide insights into species habitat preferences related to distance-to-reef gradients at a small scale. As such, eDNA can not only recover large compositional differences among islands but also help understand habitat selection and distribution of marine species at a finer spatial scale.

The Tree of Life eDNA metabarcoding reveals a similar taxonomic richness but dissimilar evolutionary lineages between seaports and marine reserves.

Coastal areas host a major part of marine biodiversity but are seriously threatened by ever-increasing human pressures. Transforming natural coastlines into urban seascapes through habitat artificialization may result in loss of biodiversity and key ecosystem functions. Yet, the extent to which seaports differ from nearby natural habitats and marine reserves across the whole Tree of Life is still unknown. This study aimed to assess the level of α and ß-diversity between seaports and reserves, and whether these biodiversity patterns are conserved across taxa and evolutionary lineages. For that, we used environmental DNA (eDNA) metabarcoding to survey six seaports on the French Mediterranean coast and four strictly no-take marine reserves nearby. By targeting four different groups-prokaryotes, eukaryotes, metazoans and fish-with appropriate markers, we provide a holistic view of biodiversity on contrasted habitats. In the absence of comprehensive reference databases, we used bioinformatic pipelines to gather similar sequences into molecular operational taxonomic units (MOTUs). In contrast to our expectations, we obtained no difference in MOTU richness (α-diversity) between habitats except for prokaryotes and threatened fishes with higher diversity in reserves than in seaports. However, we observed a marked dissimilarity (ß-diversity) between seaports and reserves for all taxa. Surprisingly, this biodiversity signature of seaports was preserved across the Tree of Life, up to the order. This result reveals that seaports and nearby marine reserves share few taxa and evolutionary lineages along urbanized coasts and suggests major differences in terms of ecosystem functioning between both habitats.

Rarity mediates species-specific responses of tropical reef fishes to protection.

Marine protected areas (MPAs) are the most widely applied tool for marine biodiversity conservation, yet many gaps remain in our understanding of their species-specific effects, partly because the socio-environmental context and spatial autocorrelation may blur and bias perceived conservation outcomes. Based on a large data set of nearly 3000 marine fish surveys spanning all tropical regions of the world, we build spatially explicit models for 658 fish species to estimate species-specific responses to protection while controlling for the environmental, habitat and socio-economic contexts experienced across their geographic ranges. We show that the species responses are highly variable, with ~40% of fishes not benefitting from protection. When investigating how traits influence species' responses, we find that rare top-predators and small herbivores benefit the most from MPAs while mid-trophic level species benefit to a lesser extent, and rare large herbivores experience adverse effects, indicating potential trophic cascades.

Optimizing detectability of the endangered fan mussel using eDNA and ddPCR.

Spatial and temporal monitoring of species threatened with extinction is of critical importance for conservation and ecosystem management. In the Mediterranean coast, the fan mussel (Pinna nobilis) is listed as critically endangered after suffering from a mass mortality event since 2016, leading to 100% mortality in most marine populations. Conventional monitoring for this macroinvertebrate is done using scuba, which is challenging in dense meadows or with low visibility. Here we developed an environmental DNA assay targeting the fan mussel and assessed the influence of several environmental parameters on the species detectability in situ. We developed and tested an eDNA molecular marker and collected 48 water samples in two sites at the Thau lagoon (France) with distinct fan mussel density, depths and during two seasons (summer and autumn). Our marker can amplify fan mussel DNA but lacks specificity since it also amplifies a conspecific species (Pinna rudis). We successfully amplified fan mussel DNA from in situ samples with 46 positive samples (out of 48) using ddPCR, although the DNA concentrations measured were low over almost all samples. Deeper sampling depth slightly increased DNA concentrations, but no seasonal effect was found. We highlight a putative spawning event on a single summer day with much higher DNA concentration compared to all other samples. We present an eDNA molecular assay able to detect the endangered fan mussel and provide guidelines to optimize the sampling protocol to maximize detectability. Effective and non-invasive monitoring tools for endangered species are promising to monitor remaining populations and have the potential of ecological restoration or habitat recolonization following a mass mortality event.

Diversity across organisational scale emerges through dispersal ability and speciation dynamics in tropical fish.

BACKGROUND: Biodiversity exists at different levels of organisation: e.g. genetic, individual, population, species, and community. These levels of organisation all exist within the same system, with diversity patterns emerging across organisational scales through several key processes. Despite this inherent interconnectivity, observational studies reveal that diversity patterns across levels are not consistent and the underlying mechanisms for variable continuity in diversity across levels remain elusive. To investigate these mechanisms, we apply a spatially explicit simulation model to simulate the global diversification of tropical reef fishes at both the population and species levels through emergent population-level processes. RESULTS: We find significant relationships between the population and species levels of diversity which vary depending on both the measure of diversity and the spatial partitioning considered. In turn, these population-species relationships are driven by modelled biological trait parameters, especially the divergence threshold at which populations speciate. CONCLUSIONS: To explain variation in multi-level diversity patterns, we propose a simple, yet novel, population-to-species diversity partitioning mechanism through speciation which disrupts continuous diversity patterns across organisational levels. We expect that in real-world systems this mechanism is driven by the molecular dynamics that determine genetic incompatibility, and therefore reproductive isolation between individuals. We put forward a framework in which the mechanisms underlying patterns of diversity across organisational levels are universal, and through this show how variable patterns of diversity can emerge through organisational scale.

New deep learning-based methods for visualizing ecosystem properties using environmental DNA metabarcoding data.

Environmental DNA (eDNA) metabarcoding provides an efficient approach for documenting biodiversity patterns in marine and terrestrial ecosystems. The complexity of these data prevents current methods from extracting and analyzing all the relevant ecological information they contain, and new methods may provide better dimensionality reduction and clustering. Here we present two new deep learning-based methods that combine different types of neural networks (NNs) to ordinate eDNA samples and visualize ecosystem properties in a two-dimensional space: the first is based on variational autoencoders and the second on deep metric learning. The strength of our new methods lies in the combination of two inputs: the number of sequences found for each molecular operational taxonomic unit (MOTU) detected and their corresponding nucleotide sequence. Using three different datasets, we show that our methods accurately represent several biodiversity indicators in a two-dimensional latent space: MOTU richness per sample, sequence α-diversity per sample, Jaccard's and sequence ß-diversity between samples. We show that our nonlinear methods are better at extracting features from eDNA datasets while avoiding the major biases associated with eDNA. Our methods outperform traditional dimension reduction methods such as Principal Component Analysis, t-distributed Stochastic Neighbour Embedding, Nonmetric Multidimensional Scaling and Uniform Manifold Approximation and Projection for dimension reduction. Our results suggest that NNs provide a more efficient way of extracting structure from eDNA metabarcoding data, thereby improving their ecological interpretation and thus biodiversity monitoring.

The global depth range of marine fishes and their genetic coverage for environmental DNA metabarcoding.

The bathymetric and geographical distribution of marine species represent a key information in biodiversity conservation. Yet, deep-sea ecosystems are among the least explored on Earth and are increasingly impacted by human activities. Environmental DNA (eDNA) metabarcoding has emerged as a promising method to study fish biodiversity but applications to the deep-sea are still scarce. A major limitation in the application of eDNA metabarcoding is the incompleteness of species sequences available in public genetic databases which reduces the extent of detected species. This incompleteness by depth is still unknown. Here, we built the global bathymetric and geographical distribution of 10,826 actinopterygian and 960 chondrichthyan fish species. We assessed their genetic coverage by depth and by ocean for three main metabarcoding markers used in the literature: teleo and MiFish-U/E. We also estimated the number of primer mismatches per species amplified by in silico polymerase chain reaction which influence the probability of species detection. Actinopterygians show a stronger decrease in species richness with depth than Chondrichthyans. These richness gradients are accompanied by a continuous species turnover between depths. Fish species coverage with the MiFish-U/E markers is higher than with teleo while threatened species are more sequenced than the others. "Deep-endemic" species, those not ascending to the shallow depth layer, are less sequenced than not threatened species. The number of primer mismatches is not higher for deep-sea species than for shallower ones. eDNA metabarcoding is promising for species detection in the deep-sea to better account for the 3-dimensional structure of the ocean in marine biodiversity monitoring and conservation. However, we argue that sequencing efforts on "deep-endemic" species are needed.

Re-thinking the environment in landscape genomics.

Detecting the extrinsic selective pressures shaping genomic variation is critical for a better understanding of adaptation and for forecasting evolutionary responses of natural populations to changing environmental conditions. With increasing availability of geo-referenced environmental data, landscape genomics provides unprecedented insights into how genomic variation and underlying gene functions affect traits potentially under selection. Yet, the robustness of genotype-environment associations used in landscape genomics remains tempered due to various limitations, including the characteristics of environmental data used, sampling designs employed, and statistical frameworks applied. Here, we argue that using complementary or new environmental data sources and well-informed sampling designs may help improve the detection of selective pressures underlying patterns of local adaptation in various organisms and environments.

Applying convolutional neural networks to speed up environmental DNA annotation in a highly diverse ecosystem.

High-throughput DNA sequencing is becoming an increasingly important tool to monitor and better understand biodiversity responses to environmental changes in a standardized and reproducible way. Environmental DNA (eDNA) from organisms can be captured in ecosystem samples and sequenced using metabarcoding, but processing large volumes of eDNA data and annotating sequences to recognized taxa remains computationally expensive. Speed and accuracy are two major bottlenecks in this critical step. Here, we evaluated the ability of convolutional neural networks (CNNs) to process short eDNA sequences and associate them with taxonomic labels. Using a unique eDNA data set collected in highly diverse Tropical South America, we compared the speed and accuracy of CNNs with that of a well-known bioinformatic pipeline (OBITools) in processing a small region (60 bp) of the 12S ribosomal DNA targeting freshwater fishes. We found that the taxonomic labels from the CNNs were comparable to those from OBITools, with high correlation levels for the composition of the regional fish fauna. The CNNs enabled the processing of raw fastq files at a rate of approximately 1 million sequences per minute, which was about 150 times faster than with OBITools. Given the good performance of CNNs in the highly diverse ecosystem considered here, the development of more elaborate CNNs promises fast deployment for future biodiversity inventories using eDNA.

Ecological and genomic vulnerability to climate change across native populations of Robusta coffee (Coffea canephora).

The assessment of population vulnerability under climate change is crucial for planning conservation as well as for ensuring food security. Coffea canephora is, in its native habitat, an understorey tree that is mainly distributed in the lowland rainforests of tropical Africa. Also known as Robusta, its commercial value constitutes a significant revenue for many human populations in tropical countries. Comparing ecological and genomic vulnerabilities within the species' native range can provide valuable insights about habitat loss and the species' adaptive potential, allowing to identify genotypes that may act as a resource for varietal improvement. By applying species distribution models, we assessed ecological vulnerability as the decrease in climatic suitability under future climatic conditions from 492 occurrences. We then quantified genomic vulnerability (or risk of maladaptation) as the allelic composition change required to keep pace with predicted climate change. Genomic vulnerability was estimated from genomic environmental correlations throughout the native range. Suitable habitat was predicted to diminish to half its size by 2050, with populations near coastlines and around the Congo River being the most vulnerable. Whole-genome sequencing revealed 165 candidate SNPs associated with climatic adaptation in C. canephora, which were located in genes involved in plant response to biotic and abiotic stressors. Genomic vulnerability was higher for populations in West Africa and in the region at the border between DRC and Uganda. Despite an overall low correlation between genomic and ecological vulnerability at broad scale, these two components of vulnerability overlap spatially in ways that may become damaging. Genomic vulnerability was estimated to be 23% higher in populations where habitat will be lost in 2050 compared to regions where habitat will remain suitable. These results highlight how ecological and genomic vulnerabilities are relevant when planning on how to cope with climate change regarding an economically important species.

Evolving spatial conservation prioritization with intraspecific genetic data.

Spatial conservation prioritization (SCP) is a planning framework used to identify new conservation areas on the basis of the spatial distribution of species, ecosystems, and their services to human societies. The ongoing accumulation of intraspecific genetic data on a variety of species offers a way to gain knowledge of intraspecific genetic diversity and to estimate several population characteristics useful in conservation, such as dispersal and population size. Here, we review how intraspecific genetic data have been integrated into SCP and highlight their potential for identifying conservation area networks that represent intraspecific genetic diversity comprehensively and that ensure the long-term persistence of biodiversity in the face of global change.

Cross-ocean patterns and processes in fish biodiversity on coral reefs through the lens of eDNA metabarcoding.

Increasing speed and magnitude of global change threaten the world's biodiversity and particularly coral reef fishes. A better understanding of large-scale patterns and processes on coral reefs is essential to prevent fish biodiversity decline but it requires new monitoring approaches. Here, we use environmental DNA metabarcoding to reconstruct well-known patterns of fish biodiversity on coral reefs and uncover hidden patterns on these highly diverse and threatened ecosystems. We analysed 226 environmental DNA (eDNA) seawater samples from 100 stations in five tropical regions (Caribbean, Central and Southwest Pacific, Coral Triangle and Western Indian Ocean) and compared those to 2047 underwater visual censuses from the Reef Life Survey in 1224 stations. Environmental DNA reveals a higher (16%) fish biodiversity, with 2650 taxa, and 25% more families than underwater visual surveys. By identifying more pelagic, reef-associated and crypto-benthic species, eDNA offers a fresh view on assembly rules across spatial scales. Nevertheless, the reef life survey identified more species than eDNA in 47 shared families, which can be due to incomplete sequence assignment, possibly combined with incomplete detection in the environment, for some species. Combining eDNA metabarcoding and extensive visual census offers novel insights on the spatial organization of the richest marine ecosystems.

Genomic insights into the historical and contemporary demographics of the grey reef shark.

Analyses of genetic diversity can shed light on both the origins of biodiversity hotspots, as well as the conservation status of species that are impacted by human activities. With these objectives, we assembled a genomic dataset of 14,935 single nucleotide polymorphisms from 513 grey reef sharks (Carcharhinus amblyrhynchos) sampled across 17 locations in the tropical Indo-Pacific. We analysed geographic variation in genetic diversity, estimated ancient and contemporary effective population size (Ne) across sampling locations (using coalescent and linkage disequilibrium methods) and modelled the history of gene flow between the Coral Triangle and the Coral Sea. Genetic diversity decreased with distance away from the Coral Triangle and north-western Australia, implying that C. amblyrhynchos may have originated in this region. Increases in Ne were detected across almost all sampling locations 40,000-90,000 generations ago (approximately 0.6-1.5 mya, given an estimated generation time of 16.4 years), suggesting a range expansion around this time. More recent, secondary increases in Ne were inferred for the Misool and North Great Barrier Reef sampling locations, but joint modelling did not clarify whether these were due to population growth, migration, or both. Despite the greater genetic diversity and ancient Ne observed at sites around Australia and the Coral Triangle, remote reefs around north-western New Caledonia had the highest contemporary Ne, demonstrating the importance of using multiple population size assessment methods. This study provides insight into both the past and present demographics of C. amblyrhynchos and contributes to our understanding of evolution in marine biodiversity hotspots.

Adaptive potential of Coffea canephora from Uganda in response to climate change.

Understanding vulnerabilities of plant populations to climate change could help preserve their biodiversity and reveal new elite parents for future breeding programmes. To this end, landscape genomics is a useful approach for assessing putative adaptations to future climatic conditions, especially in long-lived species such as trees. We conducted a population genomics study of 207 Coffea canephora trees from seven forests along different climate gradients in Uganda. For this, we sequenced 323 candidate genes involved in key metabolic and defence pathways in coffee. Seventy-one single nucleotide polymorphisms (SNPs) were found to be significantly associated with bioclimatic variables, and were thereby considered as putatively adaptive loci. These SNPs were linked to key candidate genes, including transcription factors, like DREB-like and MYB family genes controlling plant responses to abiotic stresses, as well as other genes of organoleptic interest, such as the DXMT gene involved in caffeine biosynthesis and a putative pest repellent. These climate-associated genetic markers were used to compute genetic offsets, predicting population responses to future climatic conditions based on local climate change forecasts. Using these measures of maladaptation to future conditions, substantial levels of genetic differentiation between present and future diversity were estimated for all populations and scenarios considered. The populations from the forests Zoka and Budongo, in the northernmost zone of Uganda, appeared to have the lowest genetic offsets under all predicted climate change patterns, while populations from Kalangala and Mabira, in the Lake Victoria region, exhibited the highest genetic offsets. The potential of these findings in terms of ex situ conservation strategies are discussed.

Climate differently influences the genomic patterns of two sympatric marine fish species.

Climate influences population genetic variation in marine species. Capturing these impacts remains challenging for marine fishes which disperse over large geographical scales spanning steep environmental gradients. It requires the extensive spatial sampling of individuals or populations, representative of seascape heterogeneity, combined with a set of highly informative molecular markers capable of revealing climatic-associated genetic variations. We explored how space, dispersal and environment shape the genomic patterns of two sympatric fish species in the Mediterranean Sea, which ranks among the oceanic basins most affected by climate change and human pressure. We hypothesized that the population structure and climate-associated genomic signatures of selection would be stronger in the less mobile species, as restricted gene flow tends to facilitate the fixation of locally adapted alleles. To test our hypothesis, we genotyped two species with contrasting dispersal abilities: the white seabream Diplodus sargus and the striped red mullet Mullus surmuletus. We collected 823 individuals and used genotyping by sequencing (GBS) to detect 8,206 single nucleotide polymorphisms (SNPs) for the seabream and 2,794 for the mullet. For each species, we identified highly differentiated genomic regions (i.e. outliers) and disentangled the relative contribution of space, dispersal and environmental variables (climate, marine primary productivity) on the outliers' genetic structure to test the prevalence of gene flow and local adaptation. We observed contrasting patterns of gene flow and adaptive genetic variation between the two species. The seabream showed a distinct Alboran sea population and panmixia across the Mediterranean Sea. The mullet revealed additional differentiation within the Mediterranean Sea that was significantly correlated to summer and winter temperatures, as well as marine primary productivity. Functional annotation of the climate-associated outlier SNPs then identified candidate genes involved in heat tolerance that could be examined to further predict species' responses to climate change. Our results illustrate the key steps of a comparative seascape genomics study aiming to unravel the evolutionary processes at play in marine species, to better anticipate their response to climate change. Defining population adaptation capacities and environmental niches can then serve to incorporate evolutionary processes into species conservation planning.

How many replicates to accurately estimate fish biodiversity using environmental DNA on coral reefs?

Quantifying fish species diversity in rich tropical marine environments remains challenging. Environmental DNA (eDNA) metabarcoding is a promising tool to face this challenge through the filtering, amplification, and sequencing of DNA traces from water samples. However, because eDNA concentration is low in marine environments, the reliability of eDNA to detect species diversity can be limited. Using an eDNA metabarcoding approach to identify fish Molecular Taxonomic Units (MOTUs) with a single 12S marker, we aimed to assess how the number of sampling replicates and filtered water volume affect biodiversity estimates. We used a paired sampling design of 30 L per replicate on 68 reef transects from 8 sites in 3 tropical regions. We quantified local and regional sampling variability by comparing MOTU richness, compositional turnover, and compositional nestedness. We found strong turnover of MOTUs between replicated pairs of samples undertaken in the same location, time, and conditions. Paired samples contained non-overlapping assemblages rather than subsets of one another. As a result, non-saturated localized diversity accumulation curves suggest that even 6 replicates (180 L) in the same location can underestimate local diversity (for an area <1 km). However, sampling regional diversity using ~25 replicates in variable locations (often covering 10 s of km) often saturated biodiversity accumulation curves. Our results demonstrate variability of diversity estimates possibly arising from heterogeneous distribution of eDNA in seawater, highly skewed frequencies of eDNA traces per MOTU, in addition to variability in eDNA processing. This high compositional variability has consequences for using eDNA to monitor temporal and spatial biodiversity changes in local assemblages. Avoiding false-negative detections in future biomonitoring efforts requires increasing replicates or sampled water volume to better inform management of marine biodiversity using eDNA.

Species ecology explains the spatial components of genetic diversity in tropical reef fishes.

Generating genomic data for 19 tropical reef fish species of the Western Indian Ocean, we investigate how species ecology influences genetic diversity patterns from local to regional scales. We distinguish between the α, ß and γ components of genetic diversity, which we subsequently link to six ecological traits. We find that the α and γ components of genetic diversity are strongly correlated so that species with a high total regional genetic diversity display systematically high local diversity. The α and γ diversity components are negatively associated with species abundance recorded using underwater visual surveys and positively associated with body size. Pelagic larval duration is found to be negatively related to genetic ß diversity supporting its role as a dispersal trait in marine fishes. Deviation from the neutral theory of molecular evolution motivates further effort to understand the processes shaping genetic diversity and ultimately the diversification of the exceptional diversity of tropical reef fishes.

Opportunities and challenges of macrogenetic studies.

The rapidly emerging field of macrogenetics focuses on analysing publicly accessible genetic datasets from thousands of species to explore large-scale patterns and predictors of intraspecific genetic variation. Facilitated by advances in evolutionary biology, technology, data infrastructure, statistics and open science, macrogenetics addresses core evolutionary hypotheses (such as disentangling environmental and life-history effects on genetic variation) with a global focus. Yet, there are important, often overlooked, limitations to this approach and best practices need to be considered and adopted if macrogenetics is to continue its exciting trajectory and reach its full potential in fields such as biodiversity monitoring and conservation. Here, we review the history of this rapidly growing field, highlight knowledge gaps and future directions, and provide guidelines for further research.

Maximizing regional biodiversity requires a mosaic of protection levels.

Protected areas are the flagship management tools to secure biodiversity from anthropogenic impacts. However, the extent to which adjacent areas with distinct protection levels host different species numbers and compositions remains uncertain. Here, using reef fishes, European alpine plants, and North American birds, we show that the composition of species in adjacent Strictly Protected, Restricted, and Non-Protected areas is highly dissimilar, whereas the number of species is similar, after controlling for environmental conditions, sample size, and rarity. We find that between 12% and 15% of species are only recorded in Non-Protected areas, suggesting that a non-negligible part of regional biodiversity occurs where human activities are less regulated. For imperiled species, the proportion only recorded in Strictly Protected areas reaches 58% for fishes, 11% for birds, and 7% for plants, highlighting the fundamental and unique role of protected areas and their environmental conditions in biodiversity conservation.

Benchmarking bioinformatic tools for fast and accurate eDNA metabarcoding species identification.

Bioinformatic analysis of eDNA metabarcoding data is a crucial step toward rigorously assessing biodiversity. Many programs are now available for each step of the required analyses, but their relative abilities at providing fast and accurate species lists have seldom been evaluated. We used simulated mock communities and real fish eDNA metabarcoding data to evaluate the performance of 13 bioinformatic programs and pipelines to retrieve fish occurrence and read abundance using the 12S mt rRNA gene marker. We used four indices to compare the outputs of each program with the simulated samples: sensitivity, F-measure, root-mean-square error (RMSE) on read relative abundances, and execution time. We found marked differences among programs only for the taxonomic assignment step, both in terms of sensitivity, F-measure and RMSE. Running time was highly different between programs for each step. The fastest programs with best indices for each step were assembled into a pipeline. We compared this pipeline to pipelines constructed from existing toolboxes (OBITools, Barque, and QIIME 2). Our pipeline and Barque obtained the best performance for all indices and appear to be better alternatives to highly used pipelines for analysing fish eDNA metabarcoding data when a complete reference database is available. Analysis on real eDNA metabarcoding data also indicated differences for taxonomic assignment and execution time only. This study reveals major differences between programs during the taxonomic assignment step. The choice of algorithm for the taxonomic assignment can have a significant impact on diversity estimates and should be made according to the objectives of the study.