Accurate multiplexing and filtering for high-throughput amplicon-sequencing

Is it worthy to use environmental DNA instead of scientific trawling or video survey to monitor taxa in soft-bottom habitats?

Non-extractive techniques such as video analysis are increasingly used by scientists to study marine communities instead of extractive methods such as trawling. Currently, environmental DNA (eDNA) analysis is seen as a revolutionary tool to study taxonomic diversity. We aimed to determine which method is the most appropriate to describe fish and commercial invertebrate diversity comparing bottom trawl hauls, video transects and seawater eDNA. Our results reveal that video detected the lowest number of taxa and trawling the highest. eDNA analysis is powerful to describe marine bony fish communities, but some taxa of importance for the ecosystem such as elasmobranchs, crustaceans or molluscs are poorly detected. This may be due to several factors such as marker specificity, incomplete reference gene databases or low DNA release in the environment. For now, the various methods provide different information and none is exhaustive enough to be used alone for biodiversity characterisation.

Quantitative assessment of reef foraminifera community from metabarcoding data

Describing living community compositions is essential to monitor ecosystems in a rapidly changing world, but it is challenging to produce fast and accurate depiction of ecosystems due to methodological limitations. Morphological methods provide absolute abundances with limited throughput, whereas metabarcoding provides relative abundances of genes that may not correctly represent living communities from environmental DNA assessed with morphological methods. However, it has the potential to deliver fast descriptions of living communities provided that it is interpreted with validated species-specific calibrations and reference databases. Here, we developed a quantitative approach to retrieve from metabarcoding data the assemblages of living large benthic foraminifera (LBF), photosymbiotic calcifying protists, from Indonesian coral reefs that are under increasing anthropogenic pressure. To depict the diversity, we calculated taxon-specific correction factors to reduce biological biases by comparing surface area, biovolume and calcite volume, and the number of mitochondrial gene copies in seven common LBF species. To validate the approach, we compared calibrated datasets of morphological communities from mock samples with bulk reef sediment; both sample types were metabarcoded. The calibration of the data significantly improved the estimations of genus relative abundance, with a difference of ±5% on average, allowing for comparison of past morphological datasets with future molecular ones. Our results also highlight the application of our quantitative approach to support reef monitoring operations by capturing fine-scale processes, such as seasonal and pollution-driven dynamics, that require high-throughput sampling treatment.

Are plant traits drivers of endophytic communities in seasonally flooded tropical forests?

PREMISE

In the Amazon basin, seasonally flooded (SF) forests offer varying water constraints, providing an excellent way to investigate the role of habitat selection on microbial communities within plants. However, variations in the microbial community among host plants cannot solely be attributed to environmental factors, and how plant traits contribute to microbial assemblages remains an open question.

METHODS

We described leaf- and root-associated microbial communities using ITS2 and 16 S high-throughput sequencing and investigated the stochastic-deterministic balance shaping these community assemblies using two null models. Plant ecophysiological functioning was evaluated by focusing on 10 leaf and root traits in 72 seedlings, belonging to seven tropical SF tree species in French Guiana. We then analyzed how root and leaf traits drove the assembly of endophytic communities.

RESULTS

While both stochastic and deterministic processes governed the endophyte assembly in the leaves and roots, stochasticity prevailed. Discrepancies were found between fungi and bacteria, highlighting that these microorganisms have distinct ecological strategies within plants. Traits, especially leaf traits, host species and spatial predictors better explained diversity than composition, but they were modest predictors overall.

CONCLUSIONS

This study widens our knowledge about tree species in SF forests, a habitat sensitive to climate change, through the combined analyses of their associated microbial communities with functional traits. We emphasize the need to investigate other plant traits to better disentangle the drivers of the relationship between seedlings and their associated microbiomes, ultimately enhancing their adaptive capacities to climate change.

Assessing the potential of nematode metabarcoding for benthic monitoring of offshore oil platforms

Environmental DNA metabarcoding is gaining momentum as a time and cost-effective tool for biomonitoring and environmental impact assessment. Yet, its use as a replacement for the conventional marine benthic monitoring based on morphological analysis of macrofauna is still challenging. Here we propose to study the meiofauna, which is much better represented in sediment DNA samples. We focus on nematodes, which are the most numerous and diverse group of meiofauna. Our aim is to assess the potential of nematode metabarcoding to monitor impacts associated with offshore oil platform activities. To achieve this goal, we used nematode-optimized marker (18S V1V2-Nema) and universal eukaryotic marker (18S V9) region to analyse 252 sediment DNA samples collected near three offshore oil platforms in the North Sea. For both markers, we analysed changes in alpha and beta diversity in relation to distance from the platforms and environmental variables. We also defined three impact classes based on selected environmental variables that are associated with oil extraction activities and used random forest classifiers to compare the predictive performance of both datasets. Our results show that alpha- and beta-diversity of nematodes varies with the increasing distance from the platforms. The variables directly related to platform activity, such as Ba and THC, strongly influence the nematode community. The nematode metabarcoding data provide more robust predictive models than eukaryotic data. Furthermore, the nematode community appears more stable in time and space, as illustrated by the overlap of nematode datasets obtained from the same platform three years apart. A significative negative correlation between distance and Shannon diversity also advocates for higher performance of the V1V2-Nema over the V9. Overall, these results suggest that the sensitivity of nematodes is higher compared to the eukaryotic community. Hence, nematode metabarcoding has the potential to become an effective tool for benthic monitoring in marine environment.

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.

Template-specific optimization of NGS genotyping pipelines reveals allele-specific variation in MHC gene expression

Using high-throughput sequencing for precise genotyping of multi-locus gene families, such as the major histocompatibility complex (MHC), remains challenging, due to the complexity of the data and difficulties in distinguishing genuine from erroneous variants. Several dedicated genotyping pipelines for data from high-throughput sequencing, such as next-generation sequencing (NGS), have been developed to tackle the ensuing risk of artificially inflated diversity. Here, we thoroughly assess three such multi-locus genotyping pipelines for NGS data, the DOC method, AmpliSAS and ACACIA, using MHC class IIβ data sets of three-spined stickleback gDNA, cDNA and "artificial" plasmid samples with known allelic diversity. We show that genotyping of gDNA and plasmid samples at optimal pipeline parameters was highly accurate and reproducible across methods. However, for cDNA data, the gDNA-optimal parameter configuration yielded decreased overall genotyping precision and consistency between pipelines. Further adjustments of key clustering parameters were required tο account for higher error rates and larger variation in sequencing depth per allele, highlighting the importance of template-specific pipeline optimization for reliable genotyping of multi-locus gene families. Through accurate paired gDNA-cDNA typing and MHC-II haplotype inference, we show that MHC-II allele-specific expression levels correlate negatively with allele number across haplotypes. Lastly, sibship-assisted cDNA-typing of MHC-I revealed novel variants linked in haplotype blocks, and a higher-than-previously-reported individual MHC-I allelic diversity. In conclusion, we provide novel genotyping protocols for the three-spined stickleback MHC-I and -II genes, and evaluate the performance of popular NGS-genotyping pipelines. We also show that fine-tuned genotyping of paired gDNA-cDNA samples facilitates amplification bias-corrected MHC allele expression analysis.

Effects of bottom trawling and environmental factors on benthic bacteria, meiofauna and macrofauna communities and benthic ecosystem processes

Soft sediment marine benthic ecosystems comprise a diverse community of bacteria, meiofauna and macrofauna, which together support a range of ecosystem processes such as biogeochemical cycling. These ecosystems are also fishing grounds for demersal species that are often caught using bottom trawling. This fishing method can have deleterious effects on benthic communities by causing injury or mortality, and through alteration of sediment properties that in turn influence community structure. Although the impacts of bottom trawling on macrofauna are relatively well studied, less is known about the responses of meiofauna and bacteria to such disturbances, or how bottom trawling impacts benthic ecosystem processes. Quantifying trawling impacts against a background of natural environmental variability is also a challenge. To address these questions, we examined effects of bottom trawling and a range of environmental variables (e.g. water chemistry and physical and biochemical surface sediment properties) on a) bacterial, meiofaunal and macrofaunal community structure and b) benthic ecosystem processes (nutrient fluxes, extracellular enzyme activities and carbon turnover and degradation rates). We also investigated the link between the benthic macrofauna community and the same ecosystem processes. While there was a significant effect of bottom trawling intensity on macrofaunal community structure, the same was not seen for bacterial or meiofaunal community composition, which were more affected by environmental factors, such as surface sediment properties. The labile component of the surface sediment carbon pool was higher at highly trawled sites. Carbon degradation rates, extracellular enzyme activities, oxygen fluxes and some nutrient fluxes were significantly affected by trawling, but ecosystem processes were also strongly linked to the abundance of key bioturbators (Macoma balthica, Halicryptus spinulosus, Scoloplos armiger and Pontoporeia femorata). Although benthic ecosystems were affected by a combination of trawling and natural variability, disentangling these showed that the anthropogenic effects were clearest on the larger component of the community, i.e. macrofauna composition, and on ecosystem processes related to sedimentary carbon.

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.

DNA metabarcoding reveals spatial and temporal variation of fish eye fluke communities in lake ecosystems

Eye flukes (Diplostomidae) are diverse and abundant trematode parasites that form multi-species communities in fish with negative effects on host fitness and survival. However, the environmental factors and host-related characteristics that determine species diversity, composition, and coexistence in such communities remain poorly understood. Here, we developed a cost-effective cox1 region-specific DNA metabarcoding approach to characterize parasitic diplostomid communities in two common fish species (Eurasian perch and common roach) collected from seven temperate lakes in Estonia. We found considerable inter- and intra-lake, as well as inter-host species, variation in diplostomid communities. Sympatric host species characterization revealed that parasite communities were typically more diverse in roach than perch. Additionally, we detected five positive and two negative diplostomid species associations in roach, whereas only a single negative association was observed in perch. These results indicate that diplostomid communities in temperate lakes are complex and dynamic systems exhibiting both spatial and temporal heterogeneity. They are influenced by various environmental factors and by host-parasite and inter-parasite interactions. We expect that the described methodology facilitates ecological and biodiversity research of diplostomid parasites. It is also adaptable to other parasite groups where it could serve to improve current understanding of diversity, distribution, and interspecies interactions of other understudied taxa.

Molecular ecology of microbiomes in the wild: Common pitfalls, methodological advances and future directions

The study of microbiomes across organisms and environments has become a prominent focus in molecular ecology. This perspective article explores common challenges, methodological advancements, and future directions in the field. Key research areas include understanding the drivers of microbiome community assembly, linking microbiome composition to host genetics, exploring microbial functions, transience and spatial partitioning, and disentangling non-bacterial components of the microbiome. Methodological advancements, such as quantifying absolute abundances, sequencing complete genomes, and utilizing novel statistical approaches, are also useful tools for understanding complex microbial diversity patterns. Our aims are to encourage robust practices in microbiome studies and inspire researchers to explore the next frontier of this rapidly changing field.

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.

Plants, pollinators and their interactions under global ecological change: The role of pollen DNA metabarcoding

Anthropogenic activities are triggering global changes in the environment, causing entire communities of plants, pollinators and their interactions to restructure, and ultimately leading to species declines. To understand the mechanisms behind community shifts and declines, as well as monitoring and managing impacts, a global effort must be made to characterize plant-pollinator communities in detail, across different habitat types, latitudes, elevations, and levels and types of disturbances. Generating data of this scale will only be feasible with rapid, high-throughput methods. Pollen DNA metabarcoding provides advantages in throughput, efficiency and taxonomic resolution over traditional methods, such as microscopic pollen identification and visual observation of plant-pollinator interactions. This makes it ideal for understanding complex ecological networks and their responses to change. Pollen DNA metabarcoding is currently being applied to assess plant-pollinator interactions, survey ecosystem change and model the spatiotemporal distribution of allergenic pollen. Where samples are available from past collections, pollen DNA metabarcoding has been used to compare contemporary and past ecosystems. New avenues of research are possible with the expansion of pollen DNA metabarcoding to intraspecific identification, analysis of DNA in ancient pollen samples, and increased use of museum and herbarium specimens. Ongoing developments in sequencing technologies can accelerate progress towards these goals. Global ecological change is happening rapidly, and we anticipate that high-throughput methods such as pollen DNA metabarcoding are critical for understanding the evolutionary and ecological processes that support biodiversity, and predicting and responding to the impacts of change.

eDNA metabarcoding of avocado flowers: 'Hass' it got potential to survey arthropods in food production systems?

In the face of global biodiversity declines, surveys of beneficial and antagonistic arthropod diversity as well as the ecological services that they provide are increasingly important in both natural and agro-ecosystems. Conventional survey methods used to monitor these communities often require extensive taxonomic expertise and are time-intensive, potentially limiting their application in industries such as agriculture, where arthropods often play a critical role in productivity (e.g. pollinators, pests and predators). Environmental DNA (eDNA) metabarcoding of a novel substrate, crop flowers, may offer an accurate and high throughput alternative to aid in the detection of these managed and unmanaged taxa. Here, we compared the arthropod communities detected with eDNA metabarcoding of flowers, from an agricultural species (Persea americana-'Hass' avocado), with two conventional survey techniques: digital video recording (DVR) devices and pan traps. In total, 80 eDNA flower samples, 96 h of DVRs and 48 pan trap samples were collected. Across the three methods, 49 arthropod families were identified, of which 12 were unique to the eDNA dataset. Environmental DNA metabarcoding from flowers revealed potential arthropod pollinators, as well as plant pests and parasites. Alpha diversity levels did not differ across the three survey methods although taxonomic composition varied significantly, with only 12% of arthropod families found to be common across all three methods. eDNA metabarcoding of flowers has the potential to revolutionize the way arthropod communities are monitored in natural and agro-ecosystems, potentially detecting the response of pollinators and pests to climate change, diseases, habitat loss and other disturbances.

Amplicon Sequencing Minimal Information (ASqMI): Quality and Reporting Guidelines for Actionable Calls in Biodefense Applications

BACKGROUND

Accurate, high-confidence data is critical for assessing potential biothreat incidents. In a biothreat event, false-negative and -positive results have serious consequences. Worst case scenarios can result in unnecessary shutdowns or fatalities at an exorbitant monetary and psychological cost, respectively. Quantitative PCR assays for agents of interest have been successfully used for routine biosurveillance. Recently, there has been increased impetus for adoption of amplicon sequencing (AS) for biosurveillance because it enables discrimination of true positives from near-neighbor false positives, as well as broad, simultaneous detection of many targets in many pathogens in a high-throughput scheme. However, the high sensitivity of AS can lead to false positives. Appropriate controls and workflow reporting can help address these challenges.

OBJECTIVES

Data reporting standards are critical to data trustworthiness. The standards presented herein aim to provide a framework for method quality assessment in biodetection.

METHODS

We present a set of standards, Amplicon Sequencing Minimal Information (ASqMI), developed under the auspices of the AOAC INTERNATIONAL Stakeholder Program on Agent Detection Assays for making actionable calls in biosurveillance applications. In addition to the first minimum information guidelines for AS, we provide a controls checklist and scoring scheme to assure AS run quality and assess potential sample contamination.

RESULTS

Adoption of the ASqMI guidelines will improve data quality, help track workflow performance, and ultimately provide decision makers confidence to trust the results of this new and powerful technology.

CONCLUSION

AS workflows can provide robust, confident calls for biodetection; however, due diligence in reporting and controls are needed. The ASqMI guideline is the first AS minimum reporting guidance document that also provides the means for end users to evaluate their workflows to improve confidence.

HIGHLIGHTS

Standardized reporting guidance for actionable calls is critical to ensuring trustworthy data.

Ecological drivers of fine-scale distribution of arbuscular mycorrhizal fungi in a semiarid Mediterranean scrubland

BACKGROUND AND AIMS

Arbuscular mycorrhizal (AM) fungi enhance the uptake of water and minerals by the plant hosts, alleviating plant stress. Therefore, AM fungal-plant interactions are particularly important in drylands and other stressful ecosystems. We aimed to determine the combined and independent effects of above- and below-ground plant community attributes (i.e. diversity and composition), soil heterogeneity and spatial covariates on the spatial structure of the AM fungal communities in a semiarid Mediterranean scrubland. Furthermore, we evaluated how the phylogenetic relatedness of both plants and AM fungi shapes these symbiotic relationships.

METHODS

We characterized the composition and diversity of AM fungal and plant communities in a dry Mediterranean scrubland taxonomically and phylogenetically, using DNA metabarcoding and a spatially explicit sampling design at the plant neighbourhood scale.

KEY RESULTS

The above- and below-ground plant community attributes, soil physicochemical properties and spatial variables explained unique fractions of AM fungal diversity and composition. Mainly, variations in plant composition affected the AM fungal composition and diversity. Our results also showed that particular AM fungal taxa tended to be associated with closely related plant species, suggesting the existence of a phylogenetic signal. Although soil texture, fertility and pH affected AM fungal community assembly, spatial factors had a greater influence on AM fungal community composition and diversity than soil physicochemical properties.

CONCLUSIONS

Our results highlight that the more easily accessible above-ground vegetation is a reliable indicator of the linkages between plant roots and AM fungi. We also emphasize the importance of soil physicochemical properties in addition to below-ground plant information, while accounting for the phylogenetic relationships of both plants and fungi, because these factors improve our ability to predict the relationships between AM fungal and plant communities.

Monitoring mosquito richness in an understudied area: can environmental DNA metabarcoding be a complementary approach to adult trapping?

Mosquito surveillance programmes are essential to assess the risks of local vector-borne disease outbreaks as well as for early detection of mosquito invasion events. Surveys are usually performed with traditional sampling tools (i.e., ovitraps and dipping method for immature stages or light or decoy traps for adults). Over the past decade, numerous studies have highlighted that environmental DNA (eDNA) sampling can enhance invertebrate species detection and provide community composition metrics. However, the usefulness of eDNA for detection of mosquito species has, to date, been largely neglected. Here, we sampled water from potential larval breeding sites along a gradient of anthropogenic perturbations, from the core of an oil palm plantation to the rainforest on São Tomé Island (Gulf of Guinea, Africa). We showed that (i) species of mosquitoes could be detected via metabarcoding mostly when larvae were visible, (ii) larvae species richness was greater using eDNA than visual identification and (iii) new mosquito species were also detected by the eDNA approach. We provide a critical discussion of the pros and cons of eDNA metabarcoding for monitoring mosquito species diversity and recommendations for future research directions that could facilitate the adoption of eDNA as a tool for assessing insect vector communities.

Haplotype-level metabarcoding of freshwater macroinvertebrate species: A prospective tool for population genetic analysis

Metabarcoding is a molecular-based tool capable of large quantity high-throughput species identification from bulk samples that is a faster and more cost-effective alternative to conventional DNA-sequencing approaches. Still, further exploration and assessment of the laboratory and bioinformatics strategies are required to unlock the potential of metabarcoding-based inference of haplotype information. In this study, we assessed the inference of freshwater macroinvertebrate haplotypes from metabarcoding data in a mock sample. We also examined the influence of DNA template concentration and PCR cycle on detecting true and spurious haplotypes. We tested this strategy on a mock sample containing twenty individuals from four species with known haplotypes based on the 658-bp Folmer region of the mitochondrial cytochrome c oxidase gene. We recovered fourteen zero-radius operational taxonomic units (zOTUs) of 421-bp length, with twelve zOTUs having a 100% match with the Sanger haplotype sequences. High-quality reads relatively increased with increasing PCR cycles, and the relative abundance of each zOTU was consistent for each cycle. This suggests that increasing the PCR cycles from 24 to 64 did not affect the relative abundance of each zOTU. As metabarcoding becomes more established and laboratory protocols and bioinformatic pipelines are continuously being developed, our study demonstrated the method's ability to infer intraspecific variability while highlighting the challenges that must be addressed before its eventual application for population genetic studies.

Nanopore sequencing of PCR products enables multicopy gene family reconstruction

The importance of gene amplifications in evolution is more and more recognized. Yet, tools to study multi-copy gene families are still scarce, and many such families are overlooked using common sequencing methods. Haplotype reconstruction is even harder for polymorphic multi-copy gene families. Here, we show that all variants (or haplotypes) of a multi-copy gene family present in a single genome, can be obtained using Oxford Nanopore Technologies sequencing of PCR products, followed by steps of mapping, SNP calling and haplotyping. As a proof of concept, we acquired the sequences of highly similar variants of the cidA and cidB genes present in the genome of the Wolbachia wPip, a bacterium infecting Culex pipiens mosquitoes. Our method relies on a wide database of cid genes, previously acquired by cloning and Sanger sequencing. We addressed problems commonly faced when using mapping approaches for multi-copy gene families with highly similar variants. In addition, we confirmed that PCR amplification causes frequent chimeras which have to be carefully considered when working on families of recombinant genes. We tested the robustness of the method using a combination of bioinformatics (read simulations) and molecular biology approaches (sequence acquisitions through cloning and Sanger sequencing, specific PCRs and digital droplet PCR). When different haplotypes present within a single genome cannot be reconstructed from short reads sequencing, this pipeline confers a high throughput acquisition, gives reliable results as well as insights of the relative copy numbers of the different variants.

Changes in social behavior with MAPK2 and KCTD13/CUL3 pathways alterations in two new outbred rat models for the 16p11.2 syndromes with autism spectrum disorders

Copy number variations (CNVs) of the human 16p11.2 locus are associated with several developmental/neurocognitive syndromes. Particularly, deletion and duplication of this genetic interval are found in patients with autism spectrum disorders, intellectual disability and other psychiatric traits. The high gene density associated with the region and the strong phenotypic variability of incomplete penetrance, make the study of the 16p11.2 syndromes extremely complex. To systematically study the effect of 16p11.2 CNVs and identify candidate genes and molecular mechanisms involved in the pathophysiology, mouse models were generated previously and showed learning and memory, and to some extent social deficits. To go further in understanding the social deficits caused by 16p11.2 syndromes, we engineered deletion and duplication of the homologous region to the human 16p11.2 genetic interval in two rat outbred strains, Sprague Dawley (SD) and Long Evans (LE). The 16p11.2 rat models displayed convergent defects in social behavior and in the novel object test in male carriers from both genetic backgrounds. Interestingly major pathways affecting MAPK1 and CUL3 were found altered in the rat 16p11.2 models with additional changes in males compared to females. Altogether, the consequences of the 16p11.2 genetic region dosage on social behavior are now found in three different species: humans, mice and rats. In addition, the rat models pointed to sexual dimorphism with lower severity of phenotypes in rat females compared to male mutants. This phenomenon is also observed in humans. We are convinced that the two rat models will be key to further investigating social behavior and understanding the brain mechanisms and specific brain regions that are key to controlling social behavior.

Assessing without harvesting: Pros and cons of environmental DNA sampling and image analysis for marine biodiversity evaluation

Marine stock assessments or biodiversity monitoring studies, which historically relied on extractive techniques (e.g., trawl or grab surveys), are being progressively replaced by non-extractive approaches. For instance, species abundance indices can be calculated using data obtained from high-definition underwater cameras that enable to identify taxa at low taxonomical level. In biodiversity studies, environmental DNA (eDNA) has proven to be a useful tool for characterising fish species richness. However, several marine phyla remain poorly represented in reference gene databases or release limited amounts of DNA, restricting their detection. The absence of amplification of some invertebrate taxa might also reflect primer bias. We here explore and compare the performance of eDNA and image data in describing the marine communities of several sites in the Bay of Biscay. This was achieved by deploying a remotely operated vehicle to both record images and collect seawater samples. A total of 88 taxa were identified from the eDNA samples and 121 taxa from the images. For both methods, the best characterised phylum was Chordata, with 29 and 27 Actinopterygii species detected using image versus eDNA, respectively. Neither Bryozoa nor Cnidaria was detected in the eDNA samples while the phyla were easily identifiable by imagery. Similarly, Asteroidea (Echinodermata) and Cephalopoda (Mollusca) were scarcely detected in the eDNA samples but present on the images, while Annelida were mostly identified by eDNA (18 taxa vs 7 taxa from imagery). The complementary community descriptions we highlight from these two methods therefore advocate for using both eDNA and imagery in tandem in order to capture the macroscopic biodiversity of bentho-demersal marine communities.

The topological nature of tag jumping in environmental DNA metabarcoding studies

Metabarcoding of environmental DNA constitutes a state-of-the-art tool for environmental studies. One fundamental principle implicit in most metabarcoding studies is that individual sample amplicons can still be identified after being pooled with others-based on their unique combinations of tags-during the so-called demultiplexing step that follows sequencing. Nevertheless, it has been recognized that tags can sometimes be changed (i.e., tag jumping), which ultimately leads to sample crosstalk. Here, using four DNA metabarcoding data sets derived from the analysis of soils and sediments, we show that tag jumping follows very specific and systematic patterns. Specifically, we find a strong correlation between the number of reads in blank samples and their topological position in the tag matrix (described by vertical and horizontal vectors). This observed spatial pattern of artefactual sequences could be explained by polymerase activity, which leads to the exchange of the 3' tag of single stranded tagged sequences through the formation of heteroduplexes with mixed barcodes. Importantly, tag jumping substantially distorted our data sets-despite our use of methods suggested to minimize this error. We developed a topological model to estimate the noise based on the counts in our blanks, which suggested that 40%-80% of the taxa in our soil and sedimentary samples were likely false positives introduced through tag jumping. We highlight that the amount of false positive detections caused by tag jumping strongly biased our community analyses.

Genomic analyses of hair from Ludwig van Beethoven

Ludwig van Beethoven (1770-1827) remains among the most influential and popular classical music composers. Health problems significantly impacted his career as a composer and pianist, including progressive hearing loss, recurring gastrointestinal complaints, and liver disease. In 1802, Beethoven requested that following his death, his disease be described and made public. Medical biographers have since proposed numerous hypotheses, including many substantially heritable conditions. Here we attempt a genomic analysis of Beethoven in order to elucidate potential underlying genetic and infectious causes of his illnesses. We incorporated improvements in ancient DNA methods into existing protocols for ancient hair samples, enabling the sequencing of high-coverage genomes from small quantities of historical hair. We analyzed eight independently sourced locks of hair attributed to Beethoven, five of which originated from a single European male. We deemed these matching samples to be almost certainly authentic and sequenced Beethoven's genome to 24-fold genomic coverage. Although we could not identify a genetic explanation for Beethoven's hearing disorder or gastrointestinal problems, we found that Beethoven had a genetic predisposition for liver disease. Metagenomic analyses revealed furthermore that Beethoven had a hepatitis B infection during at least the months prior to his death. Together with the genetic predisposition and his broadly accepted alcohol consumption, these present plausible explanations for Beethoven's severe liver disease, which culminated in his death. Unexpectedly, an analysis of Y chromosomes sequenced from five living members of the Van Beethoven patrilineage revealed the occurrence of an extra-pair paternity event in Ludwig van Beethoven's patrilineal ancestry.

Encapsulated in sediments: eDNA deciphers the ecosystem history of one of the most polluted European marine sites

The Anthropocene is characterized by dramatic ecosystem changes driven by human activities. The impact of these activities can be assessed by different geochemical and paleontological proxies. However, each of these proxies provides only a fragmentary insight into the effects of anthropogenic impacts. It is highly challenging to reconstruct, with a holistic view, the state of the ecosystems from the preindustrial period to the present day, covering all biological components, from prokaryotes to multicellular eukaryotes. Here, we used sedimentary ancient DNA (sedaDNA) archives encompassing all trophic levels of biodiversity to reconstruct the two century-natural history in Bagnoli-Coroglio (Gulf of Pozzuoli, Tyrrhenian Sea), one of the most polluted marine-coastal sites in Europe. The site was characterized by seagrass meadows and high eukaryotic diversity until the beginning of the 20th century. Then, the ecosystem completely changed, with seagrasses and associated fauna as well as diverse groups of planktonic and benthic protists being replaced by low diversity biota dominated by dinophyceans and infaunal metazoan species. The sedaDNA analysis revealed a five-phase evolution of the area, where changes appear as the result of a multi-level cascade effect of impacts associated with industrial activities, urbanization, water circulation and land-use changes. The sedaDNA allowed to infer reference conditions that must be considered when restoration actions are to be implemented.

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.

Ts66Yah, a mouse model of Down syndrome with improved construct and face validity

Down syndrome (DS) is caused by trisomy of human chromosome 21 (Hsa21). The understanding of genotype-phenotype relationships, the identification of driver genes and various proofs of concept for therapeutics have benefited from mouse models. The premier model, named Ts(1716)65Dn/J (Ts65Dn), displayed phenotypes related to human DS features. It carries an additional minichromosome with the Mir155 to Zbtb21 region of mouse chromosome 16, homologous to Hsa21, encompassing around 90 genes, fused to the centromeric part of mouse chromosome 17 from Pisd-ps2/Scaf8 to Pde10a, containing 46 genes not related to Hsa21. Here, we report the investigation of a new model, Ts66Yah, generated by CRISPR/Cas9 without the genomic region unrelated to Hsa21 on the minichromosome. As expected, Ts66Yah replicated DS cognitive features. However, certain phenotypes related to increased activity, spatial learning and molecular signatures were changed, suggesting genetic interactions between the Mir155-Zbtb21 and Scaf8-Pde10a intervals. Thus, Ts66Yah mice have stronger construct and face validity than Ts65Dn mice for mimicking consequences of DS genetic overdosage. Furthermore, this study is the first to demonstrate genetic interactions between triplicated regions homologous to Hsa21 and others unrelated to Hsa21. This article has an associated First Person interview with the first author of the paper.

A next generation targeted amplicon sequencing method to screen for insecticide resistance mutations in Aedes aegypti populations reveals a rdl mutation in mosquitoes from Cabo Verde

Aedes mosquito vectors transmit many viruses of global health concern, including dengue, chikungunya and Zika. These vector-borne viral diseases have a limited number of treatment options, and vaccines vary in their effectiveness. Consequently, integrated vector management is a primary strategy for disease control. However, the increasing emergence and spread of insecticide resistance is threatening the efficacy of vector control methods. Identifying mutations associated with resistance in vector populations is important to monitor the occurrence and evolution of insecticide resistance and inform control strategies. Rapid and cost-effective genome sequencing approaches are urgently needed. Here we present an adaptable targeted amplicon approach for cost-effective implementation within next generation sequencing platforms. This approach can identify single nucleotide polymorphisms (SNPs) and small insertions and deletions (indels) in genes involved in insecticide resistance in Aedes aegypti mosquitoes. We designed and tested eleven amplicons, which included segments of the ace-1 (carbamate target), the Voltage-Gated Sodium Channel (vgsc; pyrethroids, DDT and organochlorines), and rdl (dieldrin) genes; thereby covering established knockdown resistance (kdr) mutations (e.g., S989P, I1011M/V, V1016G/I and F1534C), with the potential to identify novel ones. The amplicon assays were designed with internal barcodes, to facilitate multiplexing of large numbers of mosquitoes at low cost, and were sequenced using an Illumina platform. Our approach was evaluated on 152 Ae. aegypti mosquitoes collected in Cabo Verde, an archipelago with a history of arbovirus outbreaks. The amplicon sequence data revealed 146 SNPs, including four non-synonymous polymorphisms in the vgsc gene, one in ace-1 and the 296S rdl mutation previously associated with resistance to organochlorines. The 296S rdl mutation was identified in 98% of mosquitoes screened, consistent with the past use of an organochlorine compound (e.g., DDT). Overall, our work shows that targeted amplicon sequencing is a rapid, robust, and cost-effective tool that can be used to perform high throughput monitoring of insecticide resistance.

Lack of ITS sequence homogenization in Erysimum species (Brassicaceae) with different ploidy levels

The internal transcribed spacers (ITS) exhibit concerted evolution by the fast homogenization of these sequences at the intragenomic level. However, the rate and extension of this process are unclear and might be conditioned by the number and divergence of the different ITS copies. In some cases, such as hybrid species and polyploids, ITS sequence homogenization appears incomplete, resulting in multiple haplotypes within the same organism. Here, we studied the dynamics of concerted evolution in 85 individuals of seven plant species of the genus Erysimum (Brassicaceae) with multiple ploidy levels. We estimated the rate of concerted evolution and the degree of sequence homogenization separately for ITS1 and ITS2 and whether these varied with ploidy. Our results showed incomplete sequence homogenization, especially for polyploid samples, indicating a lack of concerted evolution in these taxa. Homogenization was usually higher in ITS2 than in ITS1, suggesting that concerted evolution operates more efficiently on the former. Furthermore, the hybrid origin of several species appears to contribute to the maintenance of high haplotype diversity, regardless of the level of ploidy. These findings indicate that sequence homogenization of ITS is a dynamic and complex process that might result in varying intra- and inter-genomic diversity levels.

Novel High-Throughput DNA Part Characterization Technique for Synthetic Biology

This study presents a novel DNA part characterization technique that increases throughput by combinatorial DNA part assembly, solid plate-based quantitative fluorescence assay for phenotyping, and barcode tagging-based long-read sequencing for genotyping. We confirmed that the fluorescence intensities of colonies on plates were comparable to fluorescence at the single-cell level from a high-end, flow-cytometry device and developed a high-throughput image analysis pipeline. The barcode tagging-based long-read sequencing technique enabled rapid identification of all DNA parts and their combinations with a single sequencing experiment. Using our techniques, forty-four DNA parts (21 promoters and 23 RBSs) were successfully characterized in 72 h without any automated equipment. We anticipate that this high-throughput and easy-to-use part characterization technique will contribute to increasing part diversity and be useful for building genetic circuits and metabolic pathways in synthetic biology.

Sequencing introduced false positive rare taxa lead to biased microbial community diversity, assembly, and interaction interpretation in amplicon studies

BACKGROUND

Increasing studies have demonstrated potential disproportionate functional and ecological contributions of rare taxa in a microbial community. However, the study of the microbial rare biosphere is hampered by their inherent scarcity and the deficiency of currently available techniques. Sample-wise cross contaminations might be introduced by sample index misassignment in the most widely used metabarcoding amplicon sequencing approach. Although downstream bioinformatic quality control and clustering or denoising algorithms could remove sequencing errors and non-biological artifact reads, no algorithm could eliminate high quality reads from sample-wise cross contaminations introduced by index misassignment, making it difficult to distinguish between bona fide rare taxa and potential false positives in metabarcoding studies.

RESULTS

We thoroughly evaluated the rate of index misassignment of the widely used NovaSeq 6000 and DNBSEQ-G400 sequencing platforms using both commercial and customized mock communities, and observed significant lower (0.08% vs. 5.68%) fraction of potential false positive reads for DNBSEQ-G400 as compared to NovaSeq 6000. Significant batch effects could be caused by stochastically introduced false positive or false negative rare taxa. These false detections could also lead to inflated alpha diversity of relatively simple microbial communities and underestimated that of complex ones. Further test using a set of cow rumen samples reported differential rare taxa by different sequencing platforms. Correlation analysis of the rare taxa detected by each sequencing platform demonstrated that the rare taxa identified by DNBSEQ-G400 platform had a much higher possibility to be correlated with the physiochemical properties of rumen fluid as compared to NovaSeq 6000 platform. Community assembly mechanism and microbial network correlation analysis indicated that false positive or negative rare taxa detection could lead to biased community assembly mechanism and identification of fake keystone species of the community.

CONCLUSIONS

We highly suggest proper positive/negative/blank controls, technical replicate settings, and proper sequencing platform selection in future amplicon studies, especially when the microbial rare biosphere would be focused.

High-throughput barcoding method for the genetic surveillance of insecticide resistance and species identification in Anopheles gambiae complex malaria vectors

Surveillance of malaria vector species and the monitoring of insecticide resistance are essential to inform malaria control strategies and support the reduction of infections and disease. Genetic barcoding of mosquitoes is a useful tool to assist the high-throughput surveillance of insecticide resistance, discriminate between sibling species and to detect the presence of Plasmodium infections. In this study, we combined multiplex PCR, custom designed dual indexing, and Illumina next generation sequencing for high throughput single nucleotide polymorphism (SNP)-profiling of four species from the Anopheles (An.) gambiae complex (An. gambiae sensu stricto, An. coluzzii, An. arabiensis and An. melas). By amplifying and sequencing only 14 genetic fragments (500 bp each), we were able to simultaneously detect Plasmodium infection; insecticide resistance-conferring SNPs in ace1, gste2, vgsc and rdl genes; the partial sequences of nuclear ribosomal internal transcribed spacers (ITS1 and ITS2) and intergenic spacers (IGS), Short INterspersed Elements (SINE), as well as mitochondrial genes (cox1 and nd4) for species identification and genetic diversity. Using this amplicon sequencing approach with the four selected An. gambiae complex species, we identified a total of 15 non-synonymous mutations in the insecticide target genes, including previously described mutations associated with resistance and two new mutations (F1525L in vgsc and D148E in gste2). Overall, we present a reliable and cost-effective high-throughput panel for surveillance of An. gambiae complex mosquitoes in malaria endemic regions.

Benthic foraminifera as proxies for the environmental quality assessment of the Kuwait Bay (Kuwait, Arabian Gulf): Morphological and metabarcoding approaches

The rapid urbanization and industrialization of Kuwait and the consequent effluent discharges into marine environments have resulted in a degradation of water and sediment quality in the coastal marine ecosystems such as in the Kuwait Bay. This study investigates the ecological response of benthic foraminifera (protists) to environmental stress in the Kuwait Bay. The traditional morphological approach was compared to the innovative environmental DNA (eDNA) metabarcoding to evaluate the ecological quality status (EcoQS). Forty-six surface sediment samples were collected from selected stations in the Kuwait Bay. To detect the pollution gradient, environmental parameters from water (e.g., salinity, pH, dissolved oxygen) and sediment (e.g., grain-size, trace metals, total organic carbon, total petroleum hydrocarbons) were measured at each station. Although the foraminiferal assemblages were different in the morphological and molecular datasets, the species turnover was congruent and statistically significant. Diversity-based biotic indices derived from both morphological and metabarcoding approaches, reflect the environmental stress gradient (i.e., organic and metal contaminations) in the Kuwait Bay. The lowest values of EcoQS (i.e., bad to poor) are found in the innermost part (i.e., Sulaibikhat Bay and Ras Kazmah), while higher EcoQS values occur in the outer part of the bay. This study constitutes the first attempt to apply the foraminiferal metabarcoding to assess the EcoQS within the Arabian Gulf and presents its advantages compared to the conventional morphological approach.

The use of environmental DNA to monitor impacted coastal estuaries

Environmental DNA (eDNA) metabarcoding is increasingly being used to assess community composition in coastal ecosystems. In this study, we chose to examine temporal and spatial changes in the aquatic community of Manly Lagoon - one of the most heavily developed and polluted estuaries in eastern Australia. Based on metabarcoding of the 16S mitochondrial gene (for fish) and the 18S nuclear gene (for macroinvertebrates), we identified seasonal differences in fish and macroinvertebrate community composition as well as species richness, which correlated, in some cases, with the environmental parameters of sea surface temperature and freshwater input. Moreover, given the greater taxonomic resolution of fish versus macroinvertebrate assignments, we identified several known migratory fish species of management importance that contributed significantly to the overall patterns observed. Overall, our data support the use of eDNA metabarcoding to track fish assemblages shifting in response to environmental drivers in polluted estuaries with increased sampling and consultation with historical data.

Organelle 16S rRNA amplicon sequencing enables profiling of active gut microbiota in murine model

High-throughput sequencing of ribosomal RNA (rRNA) amplicons has served as a cornerstone in microbiome studies. Despite crucial implication of organelle 16S rRNA measurements to host gut microbial activities, genomic DNA (gDNA) was overwhelmingly targeted for amplicon sequencings. Although gDNA could be a reliable resource for gene existing validation, little information is revealed in regard to the activity of microorganisms owing to the limited changes gDNA undertaken in inactive, dormant, and dead bacteria. We applied both rRNA- and gDNA-derived sequencings on mouse cecal contents. Respective experimental designs were verified to be suitable for nucleic acid (NA) purification. Via benchmarking, mainstream 16S rRNA hypervariable region targets and reference databases were proven adequate for respective amplicon sequencing study. In phylogenetic studies, significant microbial composition differences were observed between two methods. Desulfovibrio spp. (an important group of anaerobic gut microorganisms that has caused analytical difficulties), Pediococcus spp., and Proteobacteria were drastically lower as represented by gDNA-derived compositions, while microbes like Firmicutes were higher as represented by gDNA-derived microbiome compositions. Also, using PICRUSt2 as an example, we illustrated that rRNA-derived sequencing might be more suitable for microbiome function predictions since pathways like sugar metabolism were lower as represented by rRNA-derived results. The findings of this study demonstrated that rRNA-derived amplicon sequencing could improve identification capability of specific gut microorganisms and might be more suitable for in silico microbiome function predictions. Therefore, rRNA-derived amplicon sequencings, preferably coupled with gDNA-derived ones, could be used as a capable tool to unveil active microbial components in host gut. KEY POINTS: • Conventional pipelines were adequate for the respective amplicon sequencing study • Groups, such as Desulfovibrio spp., were differently represented by two methods • Comparative amplicon sequencings could be useful in host active microbiota studies.

Toward global integration of biodiversity big data: a harmonized metabarcode data generation module for terrestrial arthropods

Metazoan metabarcoding is emerging as an essential strategy for inventorying biodiversity, with diverse projects currently generating massive quantities of community-level data. The potential for integrating across such data sets offers new opportunities to better understand biodiversity and how it might respond to global change. However, large-scale syntheses may be compromised if metabarcoding workflows differ from each other. There are ongoing efforts to improve standardization for the reporting of inventory data. However, harmonization at the stage of generating metabarcode data has yet to be addressed. A modular framework for harmonized data generation offers a pathway to navigate the complex structure of terrestrial metazoan biodiversity. Here, through our collective expertise as practitioners, method developers, and researchers leading metabarcoding initiatives to inventory terrestrial biodiversity, we seek to initiate a harmonized framework for metabarcode data generation, with a terrestrial arthropod module. We develop an initial set of submodules covering the 5 main steps of metabarcode data generation: (i) sample acquisition; (ii) sample processing; (iii) DNA extraction; (iv) polymerase chain reaction amplification, library preparation, and sequencing; and (v) DNA sequence and metadata deposition, providing a backbone for a terrestrial arthropod module. To achieve this, we (i) identified key points for harmonization, (ii) reviewed the current state of the art, and (iii) distilled existing knowledge within submodules, thus promoting best practice by providing guidelines and recommendations to reduce the universe of methodological options. We advocate the adoption and further development of the terrestrial arthropod module. We further encourage the development of modules for other biodiversity fractions as an essential step toward large-scale biodiversity synthesis through harmonization.

Drivers of Spatiotemporal Eukaryote Plankton Distribution in a Trans-Basin Water Transfer Canal in China

Planktonic eukaryotes are important components of aquatic ecosystems, and analyses of the whole eukaryotic planktonic community composition and function have far-reaching significance for water resource management. We aimed to understand the spatiotemporal variation and drivers of eukaryotic plankton distribution in the Middle Route Project of the South-to-North Water Diversion in Henan Province, China. Specifically, we examined planktonic assemblages and water quality at five stations along the canal and another one located before the dam in March, June, September, and December 2019. High-throughput sequencing revealed that the eukaryotic plankton community was primarily composed of 53 phyla, 200 genera, and 277 species, with Cryptophyta, Ciliophora, and norank_k_Cryptophyta being the dominant phyla. Redundancy analysis of the eukaryotic community and environmental factors showed that five vital factors affecting eukaryotic plankton distribution were oxidation-reduction potential, nitrate nitrogen, pH, total phosphorus, and water flow velocity. Furthermore, the geographical distribution of eukaryotic communities was consistent with the distance decay model. Importantly, environmental selection dominantly shaped the geographical distribution of the eukaryotic community. In summary, our study elucidates the ecological response of planktonic eukaryotes by identifying the diversity and ecological distribution of planktonic eukaryotes in trans-basin diversion channels.

Environmental RNA outperforms eDNA metabarcoding in assessing impact of marine pollution: A chromium-spiked mesocosm test

Environmental (e)DNA metabarcoding holds great promise for biomonitoring and ecotoxicological applications. However, few studies have compared the performance of eDNA versus eRNA metabarcoding in assessing organismal response to marine pollution, in experimental conditions. Here, we performed a chromium (Cr)-spiked mesocosm experimental test on benthic foraminiferal community to investigate the effects on species diversity by analysing both eDNA and eRNA metabarcoding data across different Cr concentrations in the sediment. Foraminiferal diversity in the eRNA data showed a significant negative correlation with the Cr concentration in the sediment, while a positive response was observed in the eDNA data. The foraminiferal OTUs exhibited a higher turnover rate in eRNA than in the eDNA-derived community. Furthermore, in the eRNA samples, OTUs abundance was significantly affected by the Cr gradient in the sediment (Pseudo-R² = 0.28, p = 0.05), while no significant trend was observed in the eDNA samples. The correlation between Cr concentration and foraminiferal diversity in eRNA datasets was stronger when the less abundant OTUs (<100 reads) were removed and the analyses were conducted exclusively on OTUs shared between eRNA and eDNA datasets. This indicates the importance of metabarcoding data filtering to capture ecological impacts, in addition to using the putatively active organisms in the eRNA dataset. The comparative analyses on foraminiferal diversity revealed that eRNA-based metabarcoding can better assess the response to heavy metal exposure in presence of subtle concentrations of the pollutant. Furthermore, our results suggest that to unlock the full potential for ecosystem assessment, eDNA and eRNA should be studied in parallel to control for potential sequence artifacts in routine ecosystem surveys.

Strategies for sample labelling and library preparation in DNA metabarcoding studies

Metabarcoding of DNA extracted from environmental or bulk specimen samples is increasingly used to profile biota in basic and applied biodiversity research because of its targeted nature that allows sequencing of genetic markers from many samples in parallel. To achieve this, PCR amplification is carried out with primers designed to target a taxonomically informative marker within a taxonomic group, and sample-specific nucleotide identifiers are added to the amplicons prior to sequencing. The latter enables assignment of the sequences back to the samples they originated from. Nucleotide identifiers can be added during the metabarcoding PCR and during "library preparation", that is, when amplicons are prepared for sequencing. Different strategies to achieve this labelling exist. All have advantages, challenges and limitations, some of which can lead to misleading results, and in the worst case compromise the fidelity of the metabarcoding data. Given the range of questions addressed using metabarcoding, ensuring that data generation is robust and fit for the chosen purpose is critically important for practitioners seeking to employ metabarcoding for biodiversity assessments. Here, we present an overview of the three main workflows for sample-specific labelling and library preparation in metabarcoding studies on Illumina sequencing platforms; one-step PCR, two-step PCR, and tagged PCR. Further, we distill the key considerations for researchers seeking to select an appropriate metabarcoding strategy for their specific study. Ultimately, by gaining insights into the consequences of different metabarcoding workflows, we hope to further consolidate the power of metabarcoding as a tool to assess biodiversity across a range of applications.

Best practices in metabarcoding of fungi: From experimental design to results

The development of high-throughput sequencing (HTS) technologies has greatly improved our capacity to identify fungi and unveil their ecological roles across a variety of ecosystems. Here we provide an overview of current best practices in metabarcoding analysis of fungal communities, from experimental design through molecular and computational analyses. By reanalysing published data sets, we demonstrate that operational taxonomic units (OTUs) outperform amplified sequence variants (ASVs) in recovering fungal diversity, a finding that is particularly evident for long markers. Additionally, analysis of the full-length ITS region allows more accurate taxonomic placement of fungi and other eukaryotes compared to the ITS2 subregion. Finally, we show that specific methods for compositional data analyses provide more reliable estimates of shifts in community structure. We conclude that metabarcoding analyses of fungi are especially promising for integrating fungi into the full microbiome and broader ecosystem functioning context, recovery of novel fungal lineages and ancient organisms as well as barcoding of old specimens including type material.

Biodiversity and vector-borne diseases: host dilution and vector amplification occur simultaneously for Amazonian leishmaniases

Changes in biodiversity may impact infectious disease transmission through multiple mechanisms. We explored the impact of biodiversity changes on the transmission of Amazonian leishmaniases, a group of wild zoonoses transmitted by phlebotomine sand flies (Psychodidae), which represent an important health burden in a region where biodiversity is both rich and threatened. Using molecular analyses of sand fly pools and blood-fed dipterans, we characterized the disease system in forest sites in French Guiana undergoing different levels of human-induced disturbance. We show that the prevalence of Leishmania parasites in sand flies correlates positively with the relative abundance of mammal species known as Leishmania reservoirs. In addition, Leishmania reservoirs tend to dominate in less diverse mammal communities, in accordance with the dilution effect hypothesis. This results in a negative relationship between Leishmania prevalence and mammal diversity. On the other hand, higher mammal diversity is associated with higher sand fly density, possibly because more diverse mammal communities harbor higher biomass and more abundant feeding resources for sand flies, although more research is needed to identify the factors that shape sand fly communities. As a consequence of these antagonistic effects, decreased mammal diversity comes with an increase of parasite prevalence in sand flies, but has no detectable impact on the density of infected sand flies. These results represent additional evidence that biodiversity changes may simultaneously dilute and amplify vector-borne disease transmission through different mechanisms that need to be better understood before drawing generalities on the biodiversity-disease relationship.

Monitoring newt communities in urban area using eDNA metabarcoding

Newts are amphibians commonly present in small ponds or garden pools in urban areas. They are protected in many countries and their presence is monitored through visual observation and/or trapping. However, newts are not easy to spot as they are small, elusive and often hidden at the bottom of water bodies. In recent years, environmental DNA (eDNA) has become a popular tool for detecting newts, with a focus on individual species using qPCR assays. Here, we assess the effectiveness of eDNA metabarcoding compared to conventional visual surveys of newt diversity in 45 ponds within urban areas of Geneva canton, Switzerland. We designed newt-specific mitochondrial 16S rRNA primers, which assign the majority of amplicons to newts, and were able to detect four species known to be present in the region, including the invasive subspecies Lissotriton vulgaris meridionalis, native to the Italian peninsula, that has been introduced in the Geneva area recently. The obtained eDNA results were congruent overall with conventional surveys, confirming the morphological observations in the majority of cases (67%). In 25% of cases, a species was only detected genetically, while in 8% of cases, the observations were not supported by eDNA metabarcoding. Our study confirms the usefulness of eDNA metabarcoding as a tool for the effective and non-invasive monitoring of newt community and suggests its broader use for the survey of newt diversity in urban area at larger scales.

A microfluidic-enabled combinatorial formulation and integrated inkjet printing platform for evaluating functionally graded material blends

Sample library preparation is a central step in the process of evaluating materials with the general aim of efficient library formulation while minimizing resource consumption. We demonstrate here the first implementation of a microfluidic-enabled thin film sample library formulation platform with integrated inkjet printing capability for directly patterning these libraries with reduced material wastage. System development and general performance screening protocol for these patterned thin films are described. We study the combinatorial formulation capabilities of this system by focusing on some practical case studies for probing the electrical conductivity in organic, biocompatible and electroactive polymer/additive (PEDOT:PSS/DMSO and PEDOT:PSS/EG) blends. Functionally-graded thin film libraries are prepared by mixing ink components and directly dispensing the processed blends into programmed geometries using the integrated platform. Electrical and morphological characterization of these printed thin film libraries is conducted to validate the formulation efficacy of the platform. Interrogating these printed libraries, we were able to iteratively identify the location of conductivity maxima for the studied blends and corroborate the morphological basis of this enhancement with established theories.

Music of metagenomics-a review of its applications, analysis pipeline, and associated tools

This humble effort highlights the intricate details of metagenomics in a simple, poetic, and rhythmic way. The paper enforces the significance of the research area, provides details about major analytical methods, examines the taxonomy and assembly of genomes, emphasizes some tools, and concludes by celebrating the richness of the ecosystem populated by the "metagenome."

Capacity of soil bacteria to reach the phyllosphere and convergence of floral communities despite soil microbiota variation

The role of flowers as environmental filters for bacterial communities and the provenance of bacteria in the phyllosphere are currently poorly understood. We experimentally tested the effect of induced variation in soil communities on the microbiota of plant organs. We identified soil-derived bacteria in the phyllosphere and show a strong convergence of floral communities with an enrichment of members of the Burkholderiaceae family. This finding highlights a potential role of the flower in shaping the interaction between plants and a bacterial family known to harbor both plant pathogens and growth-promoting strains. Because the flower involves host–symbiont feedback, the selection of specific bacteria by the reproductive organs of angiosperms could be relevant for the modulation of fruit and seed production.

Leaves and flowers are colonized by diverse bacteria that impact plant fitness and evolution. Although the structure of these microbial communities is becoming well-characterized, various aspects of their environmental origin and selection by plants remain uncertain, such as the relative proportion of soilborne bacteria in phyllosphere communities. Here, to address this issue and to provide experimental support for bacteria being filtered by flowers, we conducted common-garden experiments outside and under gnotobiotic conditions. We grew Arabidopsis thaliana in a soil substitute and added two microbial communities from natural soils. We estimated that at least 25% of the phyllosphere bacteria collected from the plants grown in the open environment were also detected in the controlled conditions, in which bacteria could reach leaves and flowers only from the soil. These taxa represented more than 40% of the communities based on amplicon sequencing. Unsupervised hierarchical clustering approaches supported the convergence of all floral microbiota, and 24 of the 28 bacteria responsible for this pattern belonged to the Burkholderiaceae family, which includes known plant pathogens and plant growth-promoting members. We anticipate that our study will foster future investigations regarding the routes used by soil microbes to reach leaves and flowers, the ubiquity of the environmental filtering of Burkholderiaceae across plant species and environments, and the potential functional effects of the accumulation of these bacteria in the reproductive organs of flowering plants.

Assessing the ecological quality status of the highly polluted Bagnoli area (Tyrrhenian Sea, Italy) using foraminiferal eDNA metabarcoding

Morphology-based benthic foraminifera indices are increasingly used worldwide for biomonitoring the ecological quality of marine sediments. The recent development of foraminiferal eDNA metabarcoding offers a reliable, time-, and cost-effective alternative to morphology-based foraminiferal biomonitoring. However, the practical applications of these new tools are still highly limited. In the present study, we evaluate the response of benthic foraminifera and define the ecological quality status (EcoQS) in the Bagnoli area (Tyrrhenian Sea, Italy) based on a traditional morphology-based approach and eDNA metabarcoding. The geochemical data show that several sites in front of the former industrial plant contain higher concentrations of potentially toxic elements than the effect range median and are characterized by the highest total organic carbon (TOC) content, whereas the distantly located sites can be considered relatively low- to unpolluted. Significant differences (i.e., diversity and assemblage composition) in both morphological and molecular datasets were found between the relatively low- to unpolluted and the most polluted areas. Similarly, the selected ecological indices of both morphological and molecular datasets strikingly and congruently resulted in a clear separation following the environmental stress gradient. The molecular indices (i.e., g-exp(H'_bc), g-Foram AMBI, and g-Foram AMBI-MOTUs) reliably identified poor-to-bad EcoQS in the polluted area in front of the former industrial plant. On the other hand, the Foram-AMBI based on morphology well identified an overall trend but seemed to overestimate the EcoQS if the traditional class boundaries were considered. The congruent and complementary trends between morphological and metabarcoding data observed in the case of the Bagnoli site further support the application of foraminiferal metabarcoding in routine biomonitoring to assess the environmental impacts of heavily polluted marine areas.

Ten millennia of hepatitis B virus evolution

[Figure: see text].

Highly comparable metabarcoding results from MGI-Tech and Illumina sequencing platforms

With the developments in DNA nanoball sequencing technologies and the emergence of new platforms, there is an increasing interest in their performance in comparison with the widely used sequencing-by-synthesis methods. Here, we test the consistency of metabarcoding results from DNBSEQ-G400RS (DNA nanoball sequencing platform by MGI-Tech) and NovaSeq 6000 (sequencing-by-synthesis platform by Illumina) platforms using technical replicates of DNA libraries that consist of COI gene amplicons from 120 soil DNA samples. By subjecting raw sequencing data from both platforms to a uniform bioinformatics processing, we found that the proportion of high-quality reads passing through the filtering steps was similar in both datasets. Per-sample operational taxonomic unit (OTU) and amplicon sequence variant (ASV) richness patterns were highly correlated, but sequencing data from DNBSEQ-G400RS harbored a higher number of OTUs. This may be related to the lower dominance of most common OTUs in DNBSEQ data set (thus revealing higher richness by detecting rare taxa) and/or to a lower effective read quality leading to generation of spurious OTUs. However, there was no statistical difference in the ASV and post-clustered ASV richness between platforms, suggesting that additional denoising step in the ASV workflow had effectively removed the 'noisy' reads. Both OTU-based and ASV-based composition were strongly correlated between the sequencing platforms, with essentially interchangeable results. Therefore, we conclude that DNBSEQ-G400RS and NovaSeq 6000 are both equally efficient high-throughput sequencing platforms to be utilized in studies aiming to apply the metabarcoding approach, but the main benefit of the former is related to lower sequencing cost.

Symbiotic germination and development of fully mycoheterotrophic plants convergently targeting similar Glomeraceae taxa

Plants producing dust seeds often meet their carbon demands by exploiting fungi at the seedling stage. This germination strategy (i.e. mycoheterotrophic germination) has been investigated among orchidaceous and ericaceous plants exploiting Ascomycota or Basidiomycota. Although several other angiosperm lineages have evolved fully mycoheterotrophic relationships with Glomeromycota, the fungal identities involved in mycoheterotrophic germination remain largely unknown. Here, we conducted in situ seed baiting and high-throughput DNA barcoding to identify mycobionts associated with seedlings of Burmannia championii (Burmanniaceae: Dioscoreales) and Sciaphila megastyla (Triuridaceae: Pandanales), which have independently evolved full mycoheterotrophy. Subsequently, we revealed that both seedlings and adults in B. championii and S. megastyla predominantly associate with Glomeraceae. However, mycorrhizal communities are somewhat distinct between seedling and adult stages, particularly in S. megastyla. Notably, the dissimilarity of mycorrhizal communities between S. megastyla adult samples and S. megastyla seedling samples is significantly higher than that between B. championi adult samples and S. megastyla adult samples, based on some indices. This pattern is possibly due to both mycorrhizal shifts during ontogenetic development and convergent recruitment of cheating-susceptible fungi. The extensive fungal overlap in two unrelated mycoheterotrophic plants indicates that both species convergently exploit specific AM fungal phylotypes.

An Economical and Flexible Dual Barcoding, Two-Step PCR Approach for Highly Multiplexed Amplicon Sequencing

In microbiome research, phylogenetic and functional marker gene amplicon sequencing is the most commonly-used community profiling approach. Consequently, a plethora of protocols for the preparation and multiplexing of samples for amplicon sequencing have been developed. Here, we present two economical high-throughput gene amplification and sequencing workflows that are implemented as standard operating procedures at the Joint Microbiome Facility of the Medical University of Vienna and the University of Vienna. These workflows are based on a previously-published two-step PCR approach, but have been updated to either increase the accuracy of results, or alternatively to achieve orders of magnitude higher numbers of samples to be multiplexed in a single sequencing run. The high-accuracy workflow relies on unique dual sample barcoding. It allows the same level of sample multiplexing as the previously-published two-step PCR approach, but effectively eliminates residual read missasignments between samples (crosstalk) which are inherent to single barcoding approaches. The high-multiplexing workflow is based on combinatorial dual sample barcoding, which theoretically allows for multiplexing up to 299,756 amplicon libraries of the same target gene in a single massively-parallelized amplicon sequencing run. Both workflows presented here are highly economical, easy to implement, and can, without significant modifications or cost, be applied to any target gene of interest.