Algorithms and strategies in short-read shotgun metagenomic reconstruction of plant communities

New prospects of environmental RNA metabarcoding research in biological diversity, ecotoxicological monitoring, and detection of COVID-19: a critical review

Ecosystems are multifaceted and complex systems and understanding their composition is crucial for the implementation of efficient conservation and management. Conventional approaches to biodiversity surveys can have limitations in detecting the complete range of species present. In contrast, the study of environmental RNA (eRNA) offers a non-invasive and comprehensive method for monitoring and evaluating biodiversity across different ecosystems. Similar to eDNA, the examination of genetic material found in environmental samples can identify and measure many species, including ones that pose challenges to traditional methods. However, eRNA is degraded quickly and therefore shows promise in detection of living organisms closer to their actual location than eDNA methods. This method provides a comprehensive perspective on the well-being of ecosystems, facilitating the development of focused conservation approaches to save at-risk species and uphold ecological equilibrium. Furthermore, eRNA has been recognized as a valuable method for the identification of COVID-19 in the environment, besides its established uses in biodiversity protection. The SARS-CoV-2 virus, which is accountable for the worldwide epidemic, releases RNA particles into the surrounding environment via human waste, providing insights into the feasibility of detecting it in wastewater and other samples taken from the environment. In this article, we critically reviewed the recent research activities that use the eRNA method, including its utilization in biodiversity conservation, ecological surveillance, and ecotoxicological monitoring as well as its innovative potential in identifying COVID-19. Through this review, the reader can understand the recent developments, prospects, and challenges of eRNA research in ecosystem management and biodiversity conservation.

Characterization of Boxwood Shoot Bacterial Communities and Potential Impact from Fungicide Treatments

Phyllosphere bacterial communities play important roles in plant fitness and growth. The objective of this study was to characterize the epiphytic and endophytic bacterial communities of boxwood shoots and determine how they may respond to commonly used fungicides. In early summer and early fall, shoot samples were collected immediately before and 1, 7, and 14 days after three fungicides containing chlorothalonil and/or propiconazole were applied to the canopy. Total genomic DNA from shoot surface washings and surface-sterilized shoot tissues was used as the template for 16S rRNA metabarcoding, and the amplicons were sequenced on a Nanopore MinION sequencer to characterize the epiphytic and endophytic communities. The bacterial communities were phylogenetically more diverse on the boxwood shoot surface than in the internal tissue, although the two communities shared 12.7% of the total 1,649 identified genera. The most abundant epiphytes were Methylobacterium and Pantoea, while Stenotrophomonas and Brevundimonas were the dominant endophytes. Fungicide treatments had strong impacts on epiphytic bacterial community structure and composition. Analysis of compositions of microbiomes with bias correction (ANCOM-BC) and analysis of variance (ANOVA)-like differential expression (ALDEx2) together identified 312 and 1,362 epiphytes changed in abundance due to fungicide treatments in early summer and early fall, respectively, and over 50% of these epiphytes were negatively impacted by fungicide. The two chlorothalonil-based contact fungicides demonstrated more marked effects than the propiconazole-based systemic fungicide. These results are foundational for exploring and utilizing the full potential of the microbiome and fungicide applications and developing a systems approach to boxwood health and production. IMPORTANCE Agrochemicals are important tools for safeguarding plants from invasive pathogens, insects, mites, and weeds. How they may affect the plant microbiome, a critical component of crop health and production, was poorly understood. Here, we used boxwood, an iconic low-maintenance landscape plant, to characterize shoot epiphytic and endophytic bacterial communities and their responses to contact and systemic fungicides. This study expanded our understanding of the above-ground microbiome in ornamental plants and is foundational for utilizing the full benefits of the microbiome in concert with different fungicide chemistries to improve boxwood health. This study also sets an example for a more thorough evaluation of these and other agrochemicals for their effects on boxwood microbiomes during production and offers an expanded systems approach that could be used with other crops for enhanced integrated pest management.

MetaDamage tool: Examining post-mortem damage in sedaDNA on a metagenomic scale

The use of metagenomic datasets to support ancient sedimentary DNA (sedaDNA) for paleoecological reconstruction has been demonstrated to be a powerful tool to understand multi-organism responses to climatic shifts and events. Authentication remains integral to the ancient DNA discipline, and this extends to sedaDNA analysis. Furthermore, distinguishing authentic sedaDNA from contamination or modern material also allows for a better understanding of broader questions in sedaDNA research, such as formation processes, source and catchment, and post-depositional processes. Existing tools for the detection of damage signals are designed for single-taxon input, require a priori organism specification, and require a significant number of input sequences to establish a signal. It is therefore often difficult to identify an established cytosine deamination rate consistent with ancient DNA across a sediment sample. In this study, we present MetaDamage, a tool that examines cytosine deamination on a metagenomic (all organisms) scale for multiple previously undetermined taxa and can produce a damage profile based on a few hundred reads. We outline the development and testing of the MetaDamage tool using both authentic sedaDNA sequences and simulated data to demonstrate the resolution in which MetaDamage can identify deamination levels consistent with the presence of ancient DNA. The MetaDamage tool offers a method for the initial assessment of the presence of sedaDNA and a better understanding of key questions of preservation for paleoecological reconstruction.

Drastic reduction of false positive species in samples of insects by intersecting the default output of two popular metagenomic classifiers

The use of high-throughput sequencing to recover short DNA reads of many species has been widely applied on biodiversity studies, either as amplicon metabarcoding or shotgun metagenomics. These reads are assigned to taxa using classifiers. However, for different reasons, the results often contain many false positives. Here we focus on the reduction of false positive species attributable to the classifiers. We benchmarked two popular classifiers, BLASTn followed by MEGAN6 (BM) and Kraken2 (K2), to analyse shotgun sequenced artificial single-species samples of insects. To reduce the number of misclassified reads, we combined the output of the two classifiers in two different ways: (1) by keeping only the reads that were attributed to the same species by both classifiers (intersection approach); and (2) by keeping the reads assigned to some species by any classifier (union approach). In addition, we applied an analytical detection limit to further reduce the number of false positives species. As expected, both metagenomic classifiers used with default parameters generated an unacceptably high number of misidentified species (tens with BM, hundreds with K2). The false positive species were not necessarily phylogenetically close, as some of them belonged to different orders of insects. The union approach failed to reduce the number of false positives, but the intersection approach got rid of most of them. The addition of an analytic detection limit of 0.001 further reduced the number to ca. 0.5 false positive species per sample. The misidentification of species by most classifiers hampers the confidence of the DNA-based methods for assessing the biodiversity of biological samples. Our approach to alleviate the problem is straightforward and significantly reduced the number of reported false positive species.

Shotgun metagenomics of soil invertebrate communities reflects taxonomy, biomass, and reference genome properties

Metagenomics - shotgun sequencing of all DNA fragments from a community DNA extract - is routinely used to describe the composition, structure, and function of microorganism communities. Advances in DNA sequencing and the availability of genome databases increasingly allow the use of shotgun metagenomics on eukaryotic communities. Metagenomics offers major advances in the recovery of biomass relationships in a sample, in comparison to taxonomic marker gene-based approaches (metabarcoding). However, little is known about the factors which influence metagenomics data from eukaryotic communities, such as differences among organism groups, the properties of reference genomes, and genome assemblies.We evaluated how shotgun metagenomics records composition and biomass in artificial soil invertebrate communities at different sequencing efforts. We generated mock communities of controlled biomass ratios from 28 species from all major soil mesofauna groups: mites, springtails, nematodes, tardigrades, and potworms. We shotgun sequenced these communities and taxonomically assigned them with a database of over 270 soil invertebrate genomes.We recovered over 95% of the species, and observed relatively high false-positive detection rates. We found strong differences in reads assigned to different taxa, with some groups (e.g., springtails) consistently attracting more hits than others (e.g., enchytraeids). Original biomass could be predicted from read counts after considering these taxon-specific differences. Species with larger genomes, and with more complete assemblies, consistently attracted more reads than species with smaller genomes. The GC content of the genome assemblies had no effect on the biomass-read relationships. Results were similar among different sequencing efforts.The results show considerable differences in taxon recovery and taxon specificity of biomass recovery from metagenomic sequence data. The properties of reference genomes and genome assemblies also influence biomass recovery, and they should be considered in metagenomic studies of eukaryotes. We show that low- and high-sequencing efforts yield similar results, suggesting high cost-efficiency of metagenomics for eukaryotic communities. We provide a brief roadmap for investigating factors which influence metagenomics-based eukaryotic community reconstructions. Understanding these factors is timely as accessibility of DNA sequencing and momentum for reference genomes projects show a future where the taxonomic assignment of DNA from any community sample becomes a reality.

Strategies for molecular authentication of herbal products: from experimental design to data analysis

Molecular herbal authentication has gained worldwide popularity in the past decade. DNA-based methods, including DNA barcoding and species-specific amplification, have been adopted for herbal identification by various pharmacopoeias. Development of next-generating sequencing (NGS) drastically increased the throughput of sequencing process and has sped up sequence collection and assembly of organelle genomes, making more and more reference sequences/genomes available. NGS allows simultaneous sequencing of multiple reads, opening up the opportunity of identifying multiple species from one sample in one go. Two major experimental approaches have been applied in recent publications of identification of herbal products by NGS, the PCR-dependent DNA metabarcoding and PCR-free genome skimming/shotgun metagenomics. This review provides a brief introduction of the use of DNA metabarcoding and genome skimming/shotgun metagenomics in authentication of herbal products and discusses some important considerations in experimental design for botanical identification by NGS, with a specific focus on quality control, reference sequence database and different taxon assignment programs. The potential of quantification or abundance estimation by NGS is discussed and new scientific findings that could potentially interfere with accurate taxon assignment and/or quantification is presented.

Hybridization capture of larch (Larix Mill.) chloroplast genomes from sedimentary ancient DNA reveals past changes of Siberian forest

Siberian larch (Larix Mill.) forests dominate vast areas of northern Russia and contribute important ecosystem services to the world. It is important to understand the past dynamics of larches in order to predict their likely response to a changing climate in the future. Sedimentary ancient DNA extracted from lake sediment cores can serve as archives to study past vegetation. However, the traditional method of studying sedimentary ancient DNA-metabarcoding-focuses on small fragments, which cannot resolve Larix to species level nor allow a detailed study of population dynamics. Here, we use shotgun sequencing and hybridization capture with long-range PCR-generated baits covering the complete Larix chloroplast genome to study Larix populations from a sediment core reaching back to 6700 years from the Taymyr region in northern Siberia. In comparison with shotgun sequencing, hybridization capture results in an increase in taxonomically classified reads by several orders of magnitude and the recovery of complete chloroplast genomes of Larix. Variation in the chloroplast reads corroborates an invasion of Larix gmelinii into the range of Larix sibirica before 6700 years ago. Since then, both species have been present at the site, although larch populations have decreased with only a few trees remaining in what was once a forested area. This study demonstrates for the first time that hybridization capture applied directly to ancient DNA of plants extracted from lake sediments can provide genome-scale information and is a viable tool for studying past genomic changes in populations of single species, irrespective of a preservation as macrofossil.

Shotgun sequencing of honey DNA can describe honey bee derived environmental signatures and the honey bee hologenome complexity

Honey bees are large-scale monitoring tools due to their extensive environmental exploration. In their activities and from the hive ecosystem complex, they get in close contact with many organisms whose traces can be transferred into the honey, which can represent an interesting reservoir of environmental DNA (eDNA) signatures and information useful to analyse the honey bee hologenome complexity. In this study, we tested a deep shotgun sequencing approach of honey DNA coupled with a specifically adapted bioinformatic pipeline. This methodology was applied to a few honey samples pointing out DNA sequences from 191 organisms spanning different kingdoms or phyla (viruses, bacteria, plants, fungi, protozoans, arthropods, mammals). Bacteria included the largest number of species. These multi-kingdom signatures listed common hive and honey bee gut microorganisms, honey bee pathogens, parasites and pests, which resembled a complex interplay that might provide a general picture of the honey bee pathosphere. Based on the Apis mellifera filamentous virus genome diversity (the most abundant detected DNA source) we obtained information that could define the origin of the honey at the apiary level. Mining Apis mellifera sequences made it possible to identify the honey bee subspecies both at the mitochondrial and nuclear genome levels.

Paleo-metagenomics of North American fossil packrat middens: Past biodiversity revealed by ancient DNA

Fossil rodent middens are powerful tools in paleoecology. In arid parts of western North America, packrat (Neotoma spp.) middens preserve plant and animal remains for tens of thousands of years. Midden contents are so well preserved that fragments of endogenous ancient DNA (aDNA) can be extracted and analyzed across millennia. Here, we explore the use of shotgun metagenomics to study the aDNA obtained from packrat middens up to 32,000 C14 years old. Eleven Illumina HiSeq 2500 libraries were successfully sequenced, and between 0.11% and 6.7% of reads were classified using Centrifuge against the NCBI "nt" database. Eukaryotic taxa identified belonged primarily to vascular plants with smaller proportions mapping to ascomycete fungi, arthropods, chordates, and nematodes. Plant taxonomic diversity in the middens is shown to change through time and tracks changes in assemblages determined by morphological examination of the plant remains. Amplicon sequencing of ITS2 and rbcL provided minimal data for some middens, but failed at amplifying the highly fragmented DNA present in others. With repeated sampling and deep sequencing, analysis of packrat midden aDNA from well-preserved midden material can provide highly detailed characterizations of past communities of plants, animals, bacteria, and fungi present as trace DNA fossils. The prospects for gaining more paleoecological insights from aDNA for rodent middens will continue to improve with optimization of laboratory methods, decreasing sequencing costs, and increasing computational power.

Metagenomic Composition Analysis of an Ancient Sequenced Polar Bear Jawbone from Svalbard

The sequencing of ancient DNA samples provides a novel way to find, characterize, and distinguish exogenous genomes of endogenous targets. After sequencing, computational composition analysis enables filtering of undesired sources in the focal organism, with the purpose of improving the quality of assemblies and subsequent data analysis. More importantly, such analysis allows extinct and extant species to be identified without requiring a specific or new sequencing run. However, the identification of exogenous organisms is a complex task, given the nature and degradation of the samples, and the evident necessity of using efficient computational tools, which rely on algorithms that are both fast and highly sensitive. In this work, we relied on a fast and highly sensitive tool, FALCON-meta, which measures similarity against whole-genome reference databases, to analyse the metagenomic composition of an ancient polar bear (Ursus maritimus) jawbone fossil. The fossil was collected in Svalbard, Norway, and has an estimated age of 110,000 to 130,000 years. The FASTQ samples contained 349 GB of nonamplified shotgun sequencing data. We identified and localized, relative to the FASTQ samples, the genomes with significant similarities to reference microbial genomes, including those of viruses, bacteria, and archaea, and to fungal, mitochondrial, and plastidial sequences. Among other striking features, we found significant similarities between modern-human, some bacterial and viral sequences (contamination) and the organelle sequences of wild carrot and tomato relative to the whole samples. For each exogenous candidate, we ran a damage pattern analysis, which in addition to revealing shallow levels of damage in the plant candidates, identified the source as contamination.