Environmental genes and genomes: understanding the differences and challenges in the approaches and software for their analyses

Diverse patterns of correspondence between protist metabarcodes and protist metagenome-assembled genomes

Two common approaches to study the composition of environmental protist communities are metabarcoding and metagenomics. Raw metabarcoding data are usually processed into Operational Taxonomic Units (OTUs) or amplicon sequence variants (ASVs) through clustering or denoising approaches, respectively. Analogous approaches are used to assemble metagenomic reads into metagenome-assembled genomes (MAGs). Understanding the correspondence between the data produced by these two approaches can help to integrate information between the datasets and to explain how metabarcoding OTUs and MAGs are related with the underlying biological entities they are hypothesised to represent. MAGs do not contain the commonly used barcoding loci, therefore sequence homology approaches cannot be used to match OTUs and MAGs. We made an attempt to match V9 metabarcoding OTUs from the 18S rRNA gene (V9 OTUs) and MAGs from the Tara Oceans expedition based on the correspondence of their relative abundances across the same set of samples. We evaluated several metrics for detecting correspondence between features in these two datasets and developed controls to filter artefacts of data structure and processing. After selecting the best-performing metrics, ranking the V9 OTU/MAG matches by their proportionality/correlation coefficients and applying a set of selection criteria, we identified candidate matches between V9 OTUs and MAGs. In some cases, V9 OTUs and MAGs could be matched with a one-to-one correspondence, implying that they likely represent the same underlying biological entity. More generally, matches we observed could be classified into 4 scenarios: one V9 OTU matches many MAGs; many V9 OTUs match many MAGs; many V9 OTUs match one MAG; one V9 OTU matches one MAG. Notably, we found some instances in which different OTU-MAG matches from the same taxonomic group were not classified in the same scenario, with all four scenarios possible even within the same taxonomic group, illustrating that factors beyond taxonomic lineage influence the relationship between OTUs and MAGs. Overall, each scenario produces a different interpretation of V9 OTUs, MAGs and how they compare in terms of the genomic and ecological diversity they represent.

In situ metagenomics: A platform for rapid sequencing and analysis of metagenomes in less than one day

We present here a complete system for metagenomic analysis that allows performing the sequencing and analysis of a medium-size metagenome in less than one day. This unprecedented development was possible due to the conjunction of state-of-the-art experimental and computational advances: a portable laboratory suitable for DNA extraction and sequencing with nanopore technology; the powerful metagenomic analysis pipeline SqueezeMeta, capable to provide a complete analysis in a few hours and using scarce computational resources; and tools for the automatic inspection of the results via a graphical user interface, that can be coupled to a web server to allow remote visualization of data (SQMtools and SQMxplore). We have tested the feasibility of our approach in the sequencing of the microbiota associated to volcanic rocks in La Palma, Canary Islands. Also, we did a two-day sampling campaign of marine waters in which the results obtained on the first day guided the experimental design of the second day. We demonstrate that it is possible to generate metagenomic information in less than one day, making it feasible to obtain taxonomic and functional profiles fast and efficiently, even in field conditions. This capacity can be used in the further to perform real-time functional and taxonomic monitoring of microbial communities in remote areas.

Metabarcoding using nanopore long-read sequencing for the unbiased characterization of apicomplexan haemoparasites

Apicomplexan haemoparasites generate significant morbidity and mortality in humans and other animals, particularly in many low-to-middle income countries. Malaria caused by Plasmodium remains responsible for some of the highest numbers of annual deaths of any human pathogen, whilst piroplasmids, such as Babesia and Theileria can have immense negative economic effects through livestock loss. Diagnosing haemoparasites via traditional methods like microscopy is challenging due to low-level and transient parasitaemia. PCR-based diagnostics overcome these limitations by being both highly sensitive and specific, but they may be unable to accurately detect coinfections or identify novel species. In contrast, next-generation sequencing (NGS)-based methods can characterize all pathogens from a group of interest concurrently, although, the short-read platforms previously used have been limited in the taxonomic resolution achievable. Here, we used Oxford Nanopore Technologies' (ONT) long-read MinION™ sequencer to conduct apicomplexan haemoparasite metabarcoding via sequencing the near full-length 18S ribosomal RNA gene, demonstrating its ability to detect Babesia, Hepatozoon, Neospora, Plasmodium, Theileria and Toxoplasma species. This method was tested on blood-extracted DNA from 100 dogs and the results benchmarked against qPCR and Illumina-based metabarcoding. For two common haemoparasites, nanopore sequencing performed as well as qPCR (kappa agreement statistics > 0.98), whilst also detecting one pathogen, Hepatozoon felis, missed by the other techniques. The long-reads obtained by nanopore sequencing provide an improved species-level taxonomic resolution whilst the method's broad applicability mean it can be used to explore apicomplexan communities from diverse mammalian hosts, on a portable sequencer that easily permits adaptation to field use.

Development and validation of a long-read metabarcoding platform for the detection of filarial worm pathogens of animals and humans

BACKGROUND

Filarial worms are important vector-borne pathogens of a large range of animal hosts, including humans, and are responsible for numerous debilitating neglected tropical diseases such as, lymphatic filariasis caused by Wuchereria bancrofti and Brugia spp., as well as loiasis caused by Loa loa. Moreover, some emerging or difficult-to-eliminate filarioid pathogens are zoonotic using animals like canines as reservoir hosts, for example Dirofilaria sp. 'hongkongensis'. Diagnosis of filariasis through commonly available methods, like microscopy, can be challenging as microfilaremia may wane below the limit of detection. In contrast, conventional PCR methods are more sensitive and specific but may show limited ability to detect coinfections as well as emerging and/or novel pathogens. Use of deep-sequencing technologies obviate these challenges, providing sensitive detection of entire parasite communities, whilst also being better suited for the characterisation of rare or novel pathogens. Therefore, we developed a novel long-read metabarcoding assay for deep-sequencing the filarial nematode cytochrome c oxidase subunit I gene on Oxford Nanopore Technologies' (ONT) MinION™ sequencer. We assessed the overall performance of our assay using kappa statistics to compare it to commonly used diagnostic methods for filarial worm detection, such as conventional PCR (cPCR) with Sanger sequencing and the microscopy-based modified Knott's test (MKT).

RESULTS

We confirmed our metabarcoding assay can characterise filarial parasites from a diverse range of genera, including, Breinlia, Brugia, Cercopithifilaria, Dipetalonema, Dirofilaria, Onchocerca, Setaria, Stephanofilaria and Wuchereria. We demonstrated proof-of-concept for this assay by using blood samples from Sri Lankan dogs, whereby we identified infections with the filarioids Acanthocheilonema reconditum, Brugia sp. Sri Lanka genotype and zoonotic Dirofilaria sp. 'hongkongensis'. When compared to traditionally used diagnostics, such as the MKT and cPCR with Sanger sequencing, we identified an additional filarioid species and over 15% more mono- and coinfections.

CONCLUSIONS

Our developed metabarcoding assay may show broad applicability for the metabarcoding and diagnosis of the full spectrum of filarioids from a wide range of animal hosts, including mammals and vectors, whilst the utilisation of ONT' small and portable MinION™ means that such methods could be deployed for field use.

The Oxford Nanopore MinION as a Versatile Technology for the Diagnosis and Characterization of Emerging Plant Viruses

The emergence of novel viral epidemics that could affect major crops represents a serious threat to global food security. The early and accurate identification of the causative viral agent is the most important step for a rapid and effective response to disease outbreaks. Over the last years, the Oxford Nanopore Technologies (ONT) MinION sequencer has been proposed as an effective diagnostic tool for the early detection and identification of emerging viruses in plants, providing many advantages compared with different high-throughput sequencing (HTS) technologies. Here, we provide a step-by-step protocol that we optimized to obtain the virome of "Lamon bean" plants (Phaseolus vulgaris L.), an agricultural product with Protected Geographical Indication (PGI) in North-East of Italy, which is frequently subjected to multiple infections caused by different RNA viruses. The conversion of viral RNA in ds-cDNA enabled the use of Genomic DNA Ligation Sequencing Kit and Native Barcoding DNA Kit, which have been originally developed for DNA sequencing. This allowed the simultaneous diagnosis of both DNA- and RNA-based pathogens, providing a more versatile alternative to the use of direct RNA and/or direct cDNA sequencing kits.

Phylogenetic diversity of functional genes in deep-sea cold seeps: a novel perspective on metagenomics

BACKGROUND

Leakages of cold, methane-rich fluids from subsurface reservoirs to the sea floor are termed cold seeps. Recent exploration of the deep sea has shed new light on the microbial communities in cold seeps. However, conventional metagenomic methods largely rely on reference databases and neglect the phylogeny of functional genes.

RESULTS

In this study, we developed the REMIRGE program to retrieve the full-length functional genes from shotgun metagenomic reads and fully explored the phylogenetic diversity in cold seep sediments. The abundance and diversity of functional genes involved in the methane, sulfur, and nitrogen cycles differed in the non-seep site and five cold seep sites. In one Haima cold seep site, the divergence of functional groups was observed at the centimeter scale of sediment depths, with the surface layer potentially acting as a reservoir of microbial species and functions. Additionally, positive correlations were found between specific gene sequence clusters of relevant genes, indicating coupling occurred within specific functional groups.

CONCLUSION

REMIRGE revealed divergent phylogenetic diversity of functional groups and functional pathway preferences in a deep-sea cold seep at finer scales, which could not be detected by conventional methods. Our work highlights that phylogenetic information is conducive to more comprehensive functional profiles, and REMIRGE has the potential to uncover more new insights from shotgun metagenomic data. Video Abstract.

Diversity, Taxonomic Novelty, and Encoded Functions of Salar de Ascotán Microbiota, as Revealed by Metagenome-Assembled Genomes

Salar de Ascotán is a high-altitude arsenic-rich salt flat exposed to high ultraviolet radiation in the Atacama Desert, Chile. It hosts unique endemic flora and fauna and is an essential habitat for migratory birds, making it an important site for conservation and protection. However, there is limited information on the resident microbiota's diversity, genomic features, metabolic potential, and molecular mechanisms that enable it to thrive in this extreme environment. We used long- and short-read metagenomics to investigate the microbial communities in Ascotán's water, sediment, and soil. Bacteria predominated, mainly Pseudomonadota, Acidobacteriota, and Bacteroidota, with a remarkable diversity of archaea in the soil. Following hybrid assembly, we recovered high-quality bacterial (101) and archaeal (6) metagenome-assembled genomes (MAGs), including representatives of two putative novel families of Patescibacteria and Pseudomonadota and two novel orders from the archaeal classes Halobacteriota and Thermoplasmata. We found different metabolic capabilities across distinct lineages and a widespread presence of genes related to stress response, DNA repair, and resistance to arsenic and other metals. These results highlight the remarkable diversity and taxonomic novelty of the Salar de Ascotán microbiota and its rich functional repertoire, making it able to resist different harsh conditions. The highly complete MAGs described here could serve future studies and bioprospection efforts focused on salt flat extremophiles, and contribute to enriching databases with microbial genome data from underrepresented regions of our planet.

Bioinformatic approaches for studying the microbiome of fermented food

High-throughput DNA sequencing-based approaches continue to revolutionise our understanding of microbial ecosystems, including those associated with fermented foods. Metagenomic and metatranscriptomic approaches are state-of-the-art biological profiling methods and are employed to investigate a wide variety of characteristics of microbial communities, such as taxonomic membership, gene content and the range and level at which these genes are expressed. Individual groups and consortia of researchers are utilising these approaches to produce increasingly large and complex datasets, representing vast populations of microorganisms. There is a corresponding requirement for the development and application of appropriate bioinformatic tools and pipelines to interpret this data. This review critically analyses the tools and pipelines that have been used or that could be applied to the analysis of metagenomic and metatranscriptomic data from fermented foods. In addition, we critically analyse a number of studies of fermented foods in which these tools have previously been applied, to highlight the insights that these approaches can provide.

Field and laboratory guidelines for reliable bioinformatic and statistical analysis of bacterial shotgun metagenomic data

Shotgun metagenomics is an increasingly cost-effective approach for profiling environmental and host-associated microbial communities. However, due to the complexity of both microbiomes and the molecular techniques required to analyze them, the reliability and representativeness of the results are contingent upon the field, laboratory, and bioinformatic procedures employed. Here, we consider 15 field and laboratory issues that critically impact downstream bioinformatic and statistical data processing, as well as result interpretation, in bacterial shotgun metagenomic studies. The issues we consider encompass intrinsic properties of samples, study design, and laboratory-processing strategies. We identify the links of field and laboratory steps with downstream analytical procedures, explain the means for detecting potential pitfalls, and propose mitigation measures to overcome or minimize their impact in metagenomic studies. We anticipate that our guidelines will assist data scientists in appropriately processing and interpreting their data, while aiding field and laboratory researchers to implement strategies for improving the quality of the generated results.

Archaeal lipids

The major archaeal membrane glycerolipids are distinguished from those of bacteria and eukaryotes by the contrasting stereochemistry of their glycerol backbones, and by the use of ether-linked isoprenoid-based alkyl chains rather than ester-linked fatty acyl chains for their hydrophobic moieties. These fascinating compounds play important roles in the extremophile lifestyles of many species, but are also present in the growing numbers of recently discovered mesophilic archaea. The past decade has witnessed significant advances in our understanding of archaea in general and their lipids in particular. Much of the new information has come from the ability to screen large microbial populations via environmental metagenomics, which has revolutionised our understanding of the extent of archaeal biodiversity that is coupled with a strict conservation of their membrane lipid compositions. Significant additional progress has come from new culturing and analytical techniques that are gradually enabling archaeal physiology and biochemistry to be studied in real time. These studies are beginning to shed light on the much-discussed and still-controversial process of eukaryogenesis, which probably involved both bacterial and archaeal progenitors. Puzzlingly, although eukaryotes retain many attributes of their putative archaeal ancestors, their lipid compositions only reflect their bacterial progenitors. Finally, elucidation of archaeal lipids and their metabolic pathways have revealed potentially interesting applications that have opened up new frontiers for biotechnological exploitation of these organisms. This review is concerned with the analysis, structure, function, evolution and biotechnology of archaeal lipids and their associated metabolic pathways.

Metabarcoding of fecal pellets in wild muskox populations reveals negative relationships between microbiome and diet alpha diversity

Microbiome diversity and diet composition concomitantly influence species health, fitness, immunity, and digestion. In environments where diet varies spatially and temporally, microbiome plasticity may promote rapid host adaptation to available resources. For northern ungulates in particular, metabarcoding of noninvasively collected fecal pellets presents unprecedented insights into their diverse ecological requirements and niches by clarifying the interrelationships of microbiomes, key to deriving nutrients, in context of altered forage availability in changing climates. Muskoxen (Ovibos moschatus) are Arctic-adapted species that experience fluctuating qualities and quantities of vegetation. Geography and seasonality have been noted to influence microbiome composition and diversity in muskoxen, yet it is unclear how their microbiomes intersect with diet. Following observations from other species, we hypothesized increasing diet diversity would result in higher microbiome diversity in muskoxen. We assessed diet composition in muskoxen using three common plant metabarcoding markers and explored correlations with microbiome data. Patterns of dietary diversity and composition were not fully concordant among the markers used, yet all reflected the primary consumption of willows and sedges. Individuals with similar diets had more similar microbiomes, yet in contrast to most literature, yielded negative relationships between microbiome and diet alpha diversity. This negative correlation may reflect the unique capacities of muskoxen to survive solely on high-fiber Arctic forage and provide insight into their resiliency to exploit changing dietary resources in a rapidly warming Arctic altering vegetation diversity.

Shepherd R, Slos D, O. Powers T, L. Porazinska D. 18S-NemaBase: Curated 18S rRNA Database of Nematode Sequences

Nematodes are the most abundant and diverse animals on the planet but lack representation in biodiversity research. This presents a problem for studying nematode diversity, particularly when molecular tools (i.e., barcoding and metabarcoding) rely on well-populated and curated reference databases, which are absent for nematodes. To improve molecular identification and the assessment of nematode diversity, we created and curated an 18S rRNA database specific to nematodes (18S-NemaBase) using sequences sourced from the most recent publicly available 18S rRNA SILVA v138 database. As part of the curation process, taxonomic strings were standardized to contain a fixed number of taxonomic ranks relevant to nematology and updated for the most recent accepted nematode classifications. In addition, apparent erroneous sequences were removed. To test the efficacy and accuracy of 18S-NemaBase, we compared it to an older but also curated SILVA v111 and the newest SILVA v138 by assigning taxonomies and analyzing the diversity of a nematode dataset from the Western Nebraska Sandhills. We showed that 18S-NemaBase provided more accurate taxonomic assignments and diversity assessments than either version of SILVA, with a much easier workflow and no need for manual corrections. Additionally, observed diversity further improved when 18S-NemaBase was supplemented with reference sequences from nematodes present in the study site. Although the 18S-NemaBase is a step in the right direction, a concerted effort to increase the number of high-quality, accessible, full-length nematode reference sequences is more important now than ever.

Advanced approaches for the diagnosis and chemoprevention of canine vector-borne pathogens and parasites-Implications for the Asia-Pacific region and beyond

Vector-borne pathogens (VBPs) of canines are a diverse range of infectious agents, including viruses, bacteria, protozoa and multicellular parasites, that are pernicious and potentially lethal to their hosts. Dogs across the globe are afflicted by canine VBPs, but the range of different ectoparasites and the VBPs that they transmit predominate in tropical regions. Countries within the Asia-Pacific have had limited prior research dedicated to exploring the epidemiology of canine VBPs, whilst the few studies that have been conducted show VBP prevalence to be high, with significant impacts on dog health. Moreover, such impacts are not restricted to dogs, as some canine VBPs are zoonotic. We reviewed the status of canine VBPs in the Asia-Pacific, with particular focus on nations in the tropics, whilst also investigating the history of VBP diagnosis and examining recent progress in the field, including advanced molecular methods, such as next-generation sequencing (NGS). These tools are rapidly changing the way parasites are detected and discovered, demonstrating a sensitivity equal to, or exceeding that of, conventional molecular diagnostics. We also provide a background to the armoury of chemopreventive products available for protecting dogs from VBP. Here, field-based research within high VBP pressure environments has underscored the importance of ectoparasiticide mode of action on their overall efficacy. The future of canine VBP diagnosis and prevention at a global level is also explored, highlighting how evolving portable sequencing technologies may permit diagnosis at point-of-care, whilst further research into chemopreventives will be essential if VBP transmission is to be effectively controlled.

Ecological and Oceanographic Perspectives in Future Marine Fungal Taxonomy

Marine fungi are an ecological rather than a taxonomic group that has been widely researched. Significant progress has been made in documenting their phylogeny, biodiversity, ultrastructure, ecology, physiology, and capacity for degradation of lignocellulosic compounds. This review (concept paper) summarizes the current knowledge of marine fungal diversity and provides an integrated and comprehensive view of their ecological roles in the world's oceans. Novel terms for 'semi marine fungi' and 'marine fungi' are proposed based on the existence of fungi in various oceanic environments. The major maritime currents and upwelling that affect species diversity are discussed. This paper also forecasts under-explored regions with a greater diversity of marine taxa based on oceanic currents. The prospects for marine and semi-marine mycology are highlighted, notably, technological developments in culture-independent sequencing approaches for strengthening our present understanding of marine fungi's ecological roles.

Moving towards improved surveillance and earlier diagnosis of aquatic pathogens: From traditional methods to emerging technologies

Early and accurate diagnosis is key to mitigating the impact of infectious diseases, along with efficient surveillance. This however is particularly challenging in aquatic environments due to hidden biodiversity and physical constraints. Traditional diagnostics, such as visual diagnosis and histopathology, are still widely used, but increasingly technological advances such as portable next generation sequencing (NGS) and artificial intelligence (AI) are being tested for early diagnosis. The most straightforward methodologies, based on visual diagnosis, rely on specialist knowledge and experience but provide a foundation for surveillance. Future computational remote sensing methods, such as AI image diagnosis and drone surveillance, will ultimately reduce labour costs whilst not compromising on sensitivity, but they require capital and infrastructural investment. Molecular techniques have advanced rapidly in the last 30 years, from standard PCR through loop-mediated isothermal amplification (LAMP) to NGS approaches, providing a range of technologies that support the currently popular eDNA diagnosis. There is now vast potential for transformative change driven by developments in human diagnostics. Here we compare current surveillance and diagnostic technologies with those that could be used or developed for use in the aquatic environment, against three gold standard ideals of high sensitivity, specificity, rapid diagnosis, and cost-effectiveness.

Trypanosomatid Richness Among Rats, Opossums, and Dogs in the Caatinga Biome, Northeast Brazil, a Former Endemic Area of Chagas Disease

Parasites are important components of the immense n-dimensional trophic network that connects all living beings because they, among others, forge biodiversity and deeply influence ecological evolution and host behavior. In this sense, the influence of Trypanosomatidae remains unknown. The aim of this study was to determine trypanosomatid infection and richness in rats, opossums, and dogs in the semiarid Caatinga biome. We submitted DNA samples from trypanosomatids obtained through axenic cultures of the blood of these mammals to mini exon multiplex-PCR, Sanger, and next-generation sequencing targeting the 18S rDNA gene. Phylogenetic analyses were performed to identify genetic diversity in the Trypanosomatidae family. Shannon, Simpson, equability, and beta-diversity indices were calculated per location and per mammalian host. Dogs were surveyed for trypanosomatid infection through hemocultures and serological assays. The examined mammal species of this area of the Caatinga biome exhibited an enormous trypanosomatid species/genotypes richness. Ten denoised Operational Taxonomic Units (ZOTUs), including three species (Trypanosoma cruzi, Trypanosoma rangeli and Crithidia mellificae) and one Trypanosoma sp. five genotypes/lineages (T. cruzi DTU TcI, TcII, and TcIV; T. rangeli A and B) and four DTU TcI haplotypes (ZOTU1, ZOTU2, ZOTU5, and ZOTU10 merged), as well as 13 Amplicon Sequence Variants (ASVs), including five species (T. cruzi, T. rangeli, C. mellificae, Trypanosoma dionisii, and Trypanosoma lainsoni), five genotypes/lineages (same as the ZOTUs) and six DTU TcI haplotypes (ASV, ASV1, ASV2, ASV3, ASV5 and ASV13), were identified in single and mixed infections. We observed that trypanosomatids present a broad host spectrum given that species related to a single host are found in other mammals from different taxa. Concomitant infections between trypanosomatids and new host-parasite relationships have been reported, and this immense diversity in mammals raised questions, such as how this can influence the course of the infection in these animals and its transmissibility. Dogs demonstrated a high infection rate by T. cruzi as observed by positive serological results (92% in 2005 and 76% in 2007). The absence of positive parasitological tests confirmed their poor infectivity potential but their importance as sentinel hosts of T. cruzi transmission.

DNA Metabarcoding and Isolation by Baiting Complement Each Other in Revealing Phytophthora Diversity in Anthropized and Natural Ecosystems

Isolation techniques supplemented by sequencing of DNA from axenic cultures have provided a robust methodology for the study of Phytophthora communities in agricultural and natural ecosystems. Recently, metabarcoding approaches have emerged as new paradigms for the detection of Phytophthora species in environmental samples. In this study, Illumina DNA metabarcoding and a conventional leaf baiting isolation technique were compared to unravel the variability of Phytophthora communities in different environments. Overall, 39 rhizosphere soil samples from a natural, a semi-natural and a horticultural small-scale ecosystem, respectively, were processed by both baiting and metabarcoding. Using both detection techniques, 28 out of 39 samples tested positive for Phytophthora. Overall, 1,406,613 Phytophthora internal transcribed spacer 1 (ITS1) sequences and 155 Phytophthora isolates were obtained, which grouped into 21 taxa, five retrieved exclusively by baiting (P. bilorbang; P. cryptogea; P. gonapodyides; P. parvispora and P. pseudocryptogea), 12 exclusively by metabarcoding (P. asparagi; P. occultans; P. psycrophila; P. syringae; P. aleatoria/P. cactorum; P. castanetorum/P. quercina; P. iranica-like; P. unknown sp. 1; P. unknown sp. 2; P. unknown sp. 3; P. unknown sp. 4; P. unknown sp. 5) and four with both techniques (P. citrophthora, P. multivora, P. nicotianae and P. plurivora). Both techniques complemented each other in describing the variability of Phytophthora communities from natural and managed ecosystems and revealing the presence of rare or undescribed Phytophthora taxa.

Abiotic Stress and Belowground Microbiome: The Potential of Omics Approaches

Nowadays, the worldwide agriculture is experiencing a transition process toward more sustainable production, which requires the reduction of chemical inputs and the preservation of microbiomes’ richness and biodiversity. Plants are no longer considered as standalone entities, and the future of agriculture should be grounded on the study of plant-associated microorganisms and all their potentiality. Moreover, due to the climate change scenario and the resulting rising incidence of abiotic stresses, an innovative and environmentally friendly technique in agroecosystem management is required to support plants in facing hostile environments. Plant-associated microorganisms have shown a great attitude as a promising tool to improve agriculture sustainability and to deal with harsh environments. Several studies were carried out in recent years looking for some beneficial plant-associated microbes and, on the basis of them, it is evident that Actinomycetes and arbuscular mycorrhizal fungi (AMF) have shown a considerable number of positive effects on plants’ fitness and health. Given the potential of these microorganisms and the effects of climate change, this review will be focused on their ability to support the plant during the interaction with abiotic stresses and on multi-omics techniques which can support researchers in unearthing the hidden world of plant–microbiome interactions. These associated microorganisms can increase plants’ endurance of abiotic stresses through several mechanisms, such as growth-promoting traits or priming-mediated stress tolerance. Using a multi-omics approach, it will be possible to deepen these mechanisms and the dynamic of belowground microbiomes, gaining fundamental information to exploit them as staunch allies and innovative weapons against crop abiotic enemies threatening crops in the ongoing global climate change context.

PHYTO-THREATS: Addressing Threats to UK Forests and Woodlands from Phytophthora; Identifying Risks of Spread in Trade and Methods for Mitigation

The multidisciplinary ‘Phyto-threats’ project was initiated in 2016 to address the increasing risks to UK forest and woodland ecosystems from trade-disseminated Phytophthora. A major component of this project was to examine the risk of Phytophthora spread through nursery and trade practices. Close to 4000 water and root samples were collected from plant nurseries located across the UK over a three-year period. Approximately half of the samples tested positive for Phytophthora DNA using a metabarcoding approach with 63 Phytophthora species identified across nurseries, including quarantine-regulated pathogens and species not previously reported in the UK. Phytophthora diversity within nurseries was linked to high-risk management practices such as use of open rather than closed water sources. Analyses of global Phytophthora risks identified biological traits and trade pathways that explained global spread and host range, and which may be of value for horizon-scanning. Phytophthoras having a higher oospore wall index and faster growth rates had wider host ranges, whereas cold-tolerant species had broader geographic and latitudinal ranges. Annual workshops revealed how stakeholder and sector ‘appetite’ for nursery accreditation increased over three years, although an exploratory cost-benefit analysis indicated that the predicted benefits of introducing best practice expected by nurseries outweigh their costs only when a wider range of pests and diseases (for example, Xylella) is considered. However, scenario analyses demonstrated the significant potential carbon costs to society from the introduction and spread of a new tree-infecting Phytophthora: Thus, the overall net benefit to society from nurseries adopting best practice could be substantial.

Capture by hybridization for full-length barcode-based eukaryotic and prokaryotic biodiversity inventories of deep sea ecosystems

Biodiversity inventory of marine systems remains limited due to unbalanced access to the three ocean dimensions. The use of environmental DNA (eDNA) for metabarcoding allows fast and effective biodiversity inventory and is forecast as a future biodiversity research and biomonitoring tool. However, in poorly understood ecosystems, eDNA results remain difficult to interpret due to large gaps in reference databases and PCR bias limiting the detection of some major phyla. Here, we aimed to circumvent these limitations by avoiding PCR and recollecting larger DNA fragments to improve assignment of detected taxa through phylogenetic reconstruction. We applied capture by hybridization (CBH) to enrich DNA from deep-sea sediment samples and compared the results with those obtained through an up-to-date metabarcoding PCR-based approach (MTB). Originally developed for bacterial communities and targeting 16S rDNA, the CBH approach was applied to 18S rDNA to improve the detection of species forming benthic communities of eukaryotes, with a particular focus on metazoans. The results confirmed the possibility of extending CBH to metazoans with two major advantages: (i) CBH revealed a broader spectrum of prokaryotic, eukaryotic, and particularly metazoan diversity, and (ii) CBH allowed much more robust phylogenetic reconstructions of full-length barcodes with up to 1900 base pairs. This is particularly important for taxa whose assignment is hampered by gaps in reference databases. This study provides a database and probes to apply 18S CBH to diverse marine systems, confirming this promising new tool to improve biodiversity assessments in data-poor ecosystems such as those in the deep sea.

Integrating eDNA metabarcoding and simultaneous underwater visual surveys to describe complex fish communities in a marine biodiversity hotspot

Marine biodiversity can be surveyed using underwater visual censuses and recently with eDNA metabarcoding. Although a promising tool, eDNA studies have shown contrasting results related to its detection scale and the number of species identified compared to other survey methods. Also, its accuracy relies on complete reference databases used for taxonomic assignment and, as other survey methods, species detection may show false-negative and false-positive errors. Here, we compared results from underwater visual censuses and simultaneous eDNA metabarcoding fish surveys in terms of observed species and community composition. We also assess the effect of a custom reference database in the taxonomic assignment, and evaluate occupancy, capture and detection probabilities, as well as error rates of eDNA survey data. We amplified a 12S rRNA fish barcode from 24 sampling sites in the gulf of California. More species were detected with eDNA metabarcoding than with UVC. Because each survey method largely detected different sets of species, the combined approach doubled the number of species registered. Both survey methods recovered a known biodiversity gradient and a biogeographic break, but eDNA captured diversity over a broader geographic and bathymetric scale. Furthermore, the use of a modest-sized custom reference database significantly increased taxonomic assignment. In a subset of species, occupancy models revealed eDNA surveys provided similar or higher detection probabilities compared to UVC. The occupancy value of each species had a large influence on eDNA detectability, and in the false positive and negative error. Overall, these results highlight the potential of eDNA metabarcoding in complementing other established ecological methods for studies of marine fishes.

An Overview of DNA-Based Applications for the Assessment of Benthic Macroinvertebrates Biodiversity in Mediterranean Aquatic Ecosystems

The loss of aquatic biodiversity is increasing at a rapid rate globally. There is a worldwide effort to protect, preserve and restore aquatic ecosystems. For efficient biodiversity monitoring and reliable management tools, comprehensive biodiversity data are required. The abundance and species diversity of benthic macroinvertebrates are commonly used as indicators of the aquatic ecosystem condition. Currently, macroinvertebrate species biodiversity assessment is based on morpho-taxonomy, which could be enhanced by recent advances in DNA-based tools for species identification. In particular, DNA metabarcoding has the potential to identify simultaneously many different taxa in a pool of species and to improve aquatic biomonitoring significantly, especially for indicator species. This review is focused on the current state of DNA-based aquatic biomonitoring using benthic macroinvertebrates in the Mediterranean region.

Daring to be differential: metabarcoding analysis of soil and plant-related microbial communities using amplicon sequence variants and operational taxonomical units

Background

Microorganisms are not only indispensable to ecosystem functioning, they are also keystones for emerging technologies. In the last 15 years, the number of studies on environmental microbial communities has increased exponentially due to advances in sequencing technologies, but the large amount of data generated remains difficult to analyze and interpret. Recently, metabarcoding analysis has shifted from clustering reads using Operational Taxonomical Units (OTUs) to Amplicon Sequence Variants (ASVs). Differences between these methods can seriously affect the biological interpretation of metabarcoding data, especially in ecosystems with high microbial diversity, as the methods are benchmarked based on low diversity datasets.

Results

In this work we have thoroughly examined the differences in community diversity, structure, and complexity between the OTU and ASV methods. We have examined culture-based mock and simulated datasets as well as soil- and plant-associated bacterial and fungal environmental communities. Four key findings were revealed. First, analysis of microbial datasets at family level guaranteed both consistency and adequate coverage when using either method. Second, the performance of both methods used are related to community diversity and sample sequencing depth. Third, differences in the method used affected sample diversity and number of detected differentially abundant families upon treatment; this may lead researchers to draw different biological conclusions. Fourth, the observed differences can mostly be attributed to low abundant (relative abundance < 0.1%) families, thus extra care is recommended when studying rare species using metabarcoding. The ASV method used outperformed the adopted OTU method concerning community diversity, especially for fungus-related sequences, but only when the sequencing depth was sufficient to capture the community complexity.

Conclusions

Investigation of metabarcoding data should be done with care. Correct biological interpretation depends on several factors, including in-depth sequencing of the samples, choice of the most appropriate filtering strategy for the specific research goal, and use of family level for data clustering.

High throughput shotgun sequencing of eRNA reveals taxonomic and derived functional shifts across a benthic productivity gradient

Benthic macrofauna is regularly used in monitoring programmes, however the vast majority of benthic eukaryotic biodiversity lies mostly in microscopic organisms, such as meiofauna (invertebrates < 1 mm) and protists, that rapidly responds to environmental change. These communities have traditionally been hard to sample and handle in the laboratory, but DNA sequencing has made such work less time consuming. While DNA sequencing captures both alive and dead organisms, environmental RNA (eRNA) better targets living organisms or organisms of recent origin in the environment. Here, we assessed the biodiversity of three known bioindicator microeukaryote groups (nematodes, foraminifera, and ciliates) in sediment samples collected at seven coastal sites along an organic carbon (OC) gradient. We aimed to investigate if eRNA shotgun sequencing can be used to simultaneously detect differences in (i) biodiversity of multiple microeukaryotic communities; and (ii) functional feeding traits of nematodes. Results showed that biodiversity was lower for nematodes and foraminifera in high OC (6.2%-6.9%), when compared to low OC sediments (1.2%-2.8%). Dissimilarity in community composition increased for all three groups between Low OC and High OC, as well as the classified feeding type of nematode genera (with more nonselective deposit feeders in high OC sediment). High relative abundant genera included nematode Sabatieria and foraminifera Elphidium in high OC, and Cryptocaryon-like ciliates in low OC sediments. Considering that future sequencing technologies are likely to decrease in cost, the use of eRNA shotgun sequencing to assess biodiversity of benthic microeukaryotes could be a powerful tool in recurring monitoring programmes.

Identifying error and accurately interpreting environmental DNA metabarcoding results: A case study to detect vertebrates at arid zone waterholes

Environmental DNA (eDNA) metabarcoding surveys enable rapid, noninvasive identification of taxa from trace samples with wide-ranging applications from characterizing local biodiversity to identifying food-web interactions. However, the technique is prone to error from two major sources: (a) contamination through foreign DNA entering the workflow, and (b) misidentification of DNA within the workflow. Both types of error have the potential to obscure true taxon presence or to increase taxonomic richness by incorrectly identifying taxa as present at sample sites, but multiple error sources can remain unaccounted for in metabarcoding studies. Here, we use data from an eDNA metabarcoding study designed to detect vertebrate species at waterholes in Australia's arid zone to illustrate where and how in the workflow errors can arise, and how to mitigate those errors. We detected the DNA of 36 taxa spanning 34 families, 19 orders and five vertebrate classes in water samples from waterholes, demonstrating the potential for eDNA metabarcoding surveys to provide rapid, noninvasive detection in remote locations, and to widely sample taxonomic diversity from aquatic through to terrestrial taxa. However, we initially identified 152 taxa in the samples, meaning there were many false positive detections. We identified the sources of these errors, allowing us to design a stepwise process to detect and remove error, and provide a template to minimize similar errors that are likely to arise in other metabarcoding studies. Our findings suggest eDNA metabarcoding surveys need to be carefully conducted and screened for errors to ensure their accuracy.

A new method of metabarcoding Microsporidia and their hosts reveals high levels of microsporidian infections in mosquitoes (Culicidae)

DNA metabarcoding offers new perspectives, especially with regard to the high‐throughput identification and diagnostics of pathogens. Microsporidia are an example of widely distributed, opportunistic and pathogenic microorganisms in which molecular identification is important for both environmental research and clinical diagnostics. We have developed a method for parallel detection of both microsporidian infection and the host species. We designed new primer sets: one specific for the classical Microsporidia (targeting the hypervariable V5 region of small subunit [ssu] rDNA), and a second one targeting a shortened fragment of the COI gene (standard metazoan DNA‐barcode); both markers are well suited for next generation sequencing. Analysis of the ssu rDNA data set representing 607 microsporidian species (120 genera) indicated that the V5 region enables identification of >98% species in the data set (596/607). To test the method, we used microsporidians that infect mosquitoes in natural populations. Using mini‐COI data, all field‐collected mosquitoes were unambiguously assigned to seven species; among them almost 60% of specimens were positive for at least 11 different microsporidian species, including a new microsporidian ssu rDNA sequence (Microsporidium sp. PL01). Phylogenetic analysis showed that this species belongs to one of the two main clades in the Terresporidia. We found a high rate of microsporidian co‐infections (9.4%). The numbers of sequence reads for the operational taxonomic units suggest that the occurrence of Nosema spp. in co‐infections could benefit them; however, this observation should be retested using a more intensive host sampling. Our results show that DNA barcoding is a rapid and cost‐effective method for deciphering sample diversity in greater resolution, including the hidden biodiversity that may be overlooked using classical methodology.

Estimation of the relative abundance of species in artificial mixtures of insects using low-coverage shotgun metagenomics

Amplicon metabarcoding is an established technique to analyse the taxonomic composition of communities of organisms using high-throughput DNA sequencing, but there are doubts about its ability to quantify the relative proportions of the species, as opposed to the species list. Here, we bypass the enrichment step and avoid the PCR-bias, by directly sequencing the extracted DNA using shotgun metagenomics. This approach is common practice in prokaryotes, but not in eukaryotes, because of the low number of sequenced genomes of eukaryotic species. We tested the metagenomics approach using insect species whose genome is already sequenced and assembled to an advanced degree. We shotgun-sequenced, at low-coverage, 18 species of insects in 22 single-species and 6 mixed-species libraries and mapped the reads against 110 reference genomes of insects. We used the single-species libraries to calibrate the process of assignation of reads to species and the libraries created from species mixtures to evaluate the ability of the method to quantify the relative species abundance. Our results showed that the shotgun metagenomic method is easily able to set apart closely-related insect species, like four species of Drosophila included in the artificial libraries. However, to avoid the counting of rare misclassified reads in samples, it was necessary to use a rather stringent detection limit of 0.001, so species with a lower relative abundance are ignored. We also identified that approximately half the raw reads were informative for taxonomic purposes. Finally, using the mixed-species libraries, we showed that it was feasible to quantify with confidence the relative abundance of individual species in the mixtures.

A novel metabarcoding diagnostic tool to explore protozoan haemoparasite diversity in mammals: a proof-of-concept study using canines from the tropics

Haemoparasites are responsible for some of the most prevalent and debilitating canine illnesses across the globe, whilst also posing a significant zoonotic risk to humankind. Nowhere are the effects of such parasites more pronounced than in developing countries in the tropics where the abundance and diversity of ectoparasites that transmit these pathogens reaches its zenith. Here we describe the use of a novel next-generation sequencing (NGS) metabarcoding based approach to screen for a range of blood-borne apicomplexan and kinetoplastid parasites from populations of temple dogs in Bangkok, Thailand. Our methodology elucidated high rates of Hepatozoon canis and Babesia vogeli infection, whilst also being able to characterise co-infections. In addition, our approach was confirmed to be more sensitive than conventional endpoint PCR diagnostic methods. Two kinetoplastid infections were also detected, including one by Trypanosoma evansi, a pathogen that is rarely screened for in dogs and another by Parabodo caudatus, a poorly documented organism that has been previously reported inhabiting the urinary tract of a dog with haematuria. Such results demonstrate the power of NGS methodologies to unearth rare and unusual pathogens, especially in regions of the world where limited information on canine vector-borne haemoparasites exist.

Assessment of a metabarcoding approach for the characterisation of vector-borne bacteria in canines from Bangkok, Thailand

BACKGROUND

Globally, bacterial vector-borne disease (VBD) exerts a large toll on dogs in terms of morbidity and mortality but nowhere is this more pronounced than in the tropics. Tropical environments permit a burgeoning diversity and abundance of ectoparasites some of which can transmit an extensive range of infectious agents, including bacteria, amongst others. Although some of these vector-borne bacteria are responsible for both animal and human diseases in the tropics, there is a scarcity of epidemiological investigation into these pathogens' prevalence. The situation is further exacerbated by frequent canine co-infection, complicating symptomatology that regular diagnostic techniques may miss or be unable to fully characterise. Such limitations draw attention to the need to develop screening tools capable of detecting a wide range of pathogens from a host simultaneously.

RESULTS

Here, we detail the employment of a next-generation sequencing (NGS) metabarcoding methodology to screen for the spectrum of bacterial VBD that are infecting semi-domesticated dogs across temple communities in Bangkok, Thailand. Our NGS detection protocol was able to find high levels of Ehrlichia canis, Mycoplasma haemocanis and Anaplasma platys infection rates as well as less common pathogens, such as "Candidatus Mycoplasma haematoparvum", Mycoplasma turicensis and Bartonella spp. We also compared our high-throughput approach to conventional endpoint PCR methods, demonstrating an improved detection ability for some bacterial infections, such as A. platys but a reduced ability to detect Rickettsia.

CONCLUSIONS

Our methodology demonstrated great strength at detecting coinfections of vector-borne bacteria and rare pathogens that are seldom screened for in canines in the tropics, highlighting its advantages over traditional diagnostics to better characterise bacterial pathogens in environments where there is a dearth of research.

Peat substrate amended with chitin modulates the N-cycle, siderophore and chitinase responses in the lettuce rhizobiome

Chitin is a valuable peat substrate amendment by increasing lettuce growth and reducing the survival of the zoonotic pathogen Salmonella enterica on lettuce leaves. The production of chitin-catabolic enzymes (chitinases) play a crucial role and are mediated through the microbial community. A higher abundance of plant-growth promoting microorganisms and genera involved in N and chitin metabolism are present in a chitin-enriched substrate. In this study, we hypothesize that chitin addition to peat substrate stimulates the microbial chitinase production. The degradation of chitin leads to nutrient release and the production of small chitin oligomers that are related to plant growth promotion and activation of the plant's defense response. First a shotgun metagenomics approach was used to decipher the potential rhizosphere microbial functions then the nutritional content of the peat substrate was measured. Our results show that chitin addition increases chitin-catabolic enzymes, bacterial ammonium oxidizing and siderophore genes. Lettuce growth promotion can be explained by a cascade degradation of chitin to N-acetylglucosamine and eventually ammonium. The occurrence of increased ammonium oxidizing bacteria, Nitrosospira, and amoA genes results in an elevated concentration of plant-available nitrate. In addition, the increase in chitinase and siderophore genes may have stimulated the plant's systemic resistance.

Metabarcoding reveals a high diversity of woody host-associated Phytophthora spp. in soils at public gardens and amenity woodlands in Britain

Forests and woodlands worldwide are being severely impacted by invasive Phytophthora species, with initial outbreaks in some cases occurring on host trees located in public parks and gardens. These highly disturbed sites with diverse planting practices may indeed act as harbours for invasive Phytophthora pathogens which are particularly well adapted to surviving in soil. High throughput Illumina sequencing was used to analyse Phytophthora species diversity in soil samples collected from 14 public garden/amenity woodland sites in northern Britain. Bioinformatic analyses revealed some limitations to using internal transcribed spacer as the barcode region; namely reporting of false positives and ambiguous species matches. Taking this into account, 35 distinct sequences were amplified across the sites, corresponding to 23 known Phytophthora species as well as twelve oomycete sequences with no match to any known Phytophthora species. Phytophthora pseudosyringae and P. austrocedri, both of which cause serious damage to trees and are regarded as fairly recent introductions to Britain, were the two most abundant Phytophthora species detected. There was no evidence that any of the detected Phytophthora species were more associated with any one type of host, healthy or otherwise. This study has demonstrated the ubiquity and diversity of Phytophthora species endemic in highly managed, extensively planted soil environments in Britain. Suggested improvements to the methodology and the practical implications of the findings in terms of mitigating Phytophthora spread and impact are discussed.

Sequence clustering in bioinformatics: an empirical study

Sequence clustering is a basic bioinformatics task that is attracting renewed attention with the development of metagenomics and microbiomics. The latest sequencing techniques have decreased costs and as a result, massive amounts of DNA/RNA sequences are being produced. The challenge is to cluster the sequence data using stable, quick and accurate methods. For microbiome sequencing data, 16S ribosomal RNA operational taxonomic units are typically used. However, there is often a gap between algorithm developers and bioinformatics users. Different software tools can produce diverse results and users can find them difficult to analyze. Understanding the different clustering mechanisms is crucial to understanding the results that they produce. In this review, we selected several popular clustering tools, briefly explained the key computing principles, analyzed their characters and compared them using two independent benchmark datasets. Our aim is to assist bioinformatics users in employing suitable clustering tools effectively to analyze big sequencing data. Related data, codes and software tools were accessible at the link http://lab.malab.cn/∼lg/clustering/.

MinION™ nanopore sequencing of environmental metagenomes: a synthetic approach

Environmental metagenomic analysis is typically accomplished by assigning taxonomy and/or function from whole genome sequencing or 16S amplicon sequences. Both of these approaches are limited, however, by read length, among other technical and biological factors. A nanopore-based sequencing platform, MinION™, produces reads that are ≥1 × 104 bp in length, potentially providing for more precise assignment, thereby alleviating some of the limitations inherent in determining metagenome composition from short reads. We tested the ability of sequence data produced by MinION (R7.3 flow cells) to correctly assign taxonomy in single bacterial species runs and in three types of low-complexity synthetic communities: a mixture of DNA using equal mass from four species, a community with one relatively rare (1%) and three abundant (33% each) components, and a mixture of genomic DNA from 20 bacterial strains of staggered representation. Taxonomic composition of the low-complexity communities was assessed by analyzing the MinION sequence data with three different bioinformatic approaches: Kraken, MG-RAST, and One Codex. Results: Long read sequences generated from libraries prepared from single strains using the version 5 kit and chemistry, run on the original MinION device, yielded as few as 224 to as many as 3497 bidirectional high-quality (2D) reads with an average overall study length of 6000 bp. For the single-strain analyses, assignment of reads to the correct genus by different methods ranged from 53.1% to 99.5%, assignment to the correct species ranged from 23.9% to 99.5%, and the majority of misassigned reads were to closely related organisms. A synthetic metagenome sequenced with the same setup yielded 714 high quality 2D reads of approximately 5500 bp that were up to 98% correctly assigned to the species level. Synthetic metagenome MinION libraries generated using version 6 kit and chemistry yielded from 899 to 3497 2D reads with lengths averaging 5700 bp with up to 98% assignment accuracy at the species level. The observed community proportions for “equal” and “rare” synthetic libraries were close to the known proportions, deviating from 0.1% to 10% across all tests. For a 20-species mock community with staggered contributions, a sequencing run detected all but 3 species (each included at <0.05% of DNA in the total mixture), 91% of reads were assigned to the correct species, 93% of reads were assigned to the correct genus, and >99% of reads were assigned to the correct family. Conclusions: At the current level of output and sequence quality (just under 4 × 103 2D reads for a synthetic metagenome), MinION sequencing followed by Kraken or One Codex analysis has the potential to provide rapid and accurate metagenomic analysis where the consortium is comprised of a limited number of taxa. Important considerations noted in this study included: high sensitivity of the MinION platform to the quality of input DNA, high variability of sequencing results across libraries and flow cells, and relatively small numbers of 2D reads per analysis limit. Together, these limited detection of very rare components of the microbial consortia, and would likely limit the utility of MinION for the sequencing of high-complexity metagenomic communities where thousands of taxa are expected. Furthermore, the limitations of the currently available data analysis tools suggest there is considerable room for improvement in the analytical approaches for the characterization of microbial communities using long reads. Nevertheless, the fact that the accurate taxonomic assignment of high-quality reads generated by MinION is approaching 99.5% and, in most cases, the inferred community structure mirrors the known proportions of a synthetic mixture warrants further exploration of practical application to environmental metagenomics as the platform continues to develop and improve. With further improvement in sequence throughput and error rate reduction, this platform shows great promise for precise real-time analysis of the composition and structure of more complex microbial communities.

Dropping Hints: Estimating the diets of livestock in rangelands using DNA metabarcoding of faeces

The introduction of domesticated animals into new environments can lead to considerable ecological disruption, and it can be difficult to predict their impact on the new ecosystem. In this study, we use faecal metabarcoding to characterize the diets of three ruminant taxa in the rangelands of south-western New South Wales, Australia. Our study organisms included goats (Capraaegagrushircus) and two breeds of sheep (Ovisaries): Merinos, which have been present in Australia for over two hundred years, and Dorpers, which were introduced in the 1990s. We used High-Throughput Sequencing methods to sequence therbcLand ITS2 genes of plants in the faecal samples, and identified the samples using the GenBank and BOLD online databases, as well as a reference collection of sequences from plants collected in the study area. We found that the diets of all three taxa were dominated by the family Malvaceae, and that the Dorper diet was more diverse than the Merino diet at both the family and the species level. We conclude that Dorpers, like Merinos, are potentially a threat to some vulnerable species in the rangelands of New South Wales.

Is there foul play in the leaf pocket? The metagenome of floating fern Azolla reveals endophytes that do not fix N2 but may denitrify

Dinitrogen fixation by Nostoc azollae residing in specialized leaf pockets supports prolific growth of the floating fern Azolla filiculoides. To evaluate contributions by further microorganisms, the A. filiculoides microbiome and nitrogen metabolism in bacteria persistently associated with Azolla ferns were characterized. A metagenomic approach was taken complemented by detection of N₂ O released and nitrogen isotope determinations of fern biomass. Ribosomal RNA genes in sequenced DNA of natural ferns, their enriched leaf pockets and water filtrate from the surrounding ditch established that bacteria of A. filiculoides differed entirely from surrounding water and revealed species of the order Rhizobiales. Analyses of seven cultivated Azolla species confirmed persistent association with Rhizobiales. Two distinct nearly full-length Rhizobiales genomes were identified in leaf-pocket-enriched samples from ditch grown A. filiculoides. Their annotation revealed genes for denitrification but not N₂ -fixation. ¹⁵ N₂ incorporation was active in ferns with N. azollae but not in ferns without. N₂ O was not detectably released from surface-sterilized ferns with the Rhizobiales. N₂ -fixing N. azollae, we conclude, dominated the microbiome of Azolla ferns. The persistent but less abundant heterotrophic Rhizobiales bacteria possibly contributed to lowering O₂ levels in leaf pockets but did not release detectable amounts of the strong greenhouse gas N₂ O.

Is there foul play in the leaf pocket? The metagenome of floating fern Azolla reveals endophytes that do not fix N2 but may denitrify

Dinitrogen fixation by Nostoc azollae residing in specialized leaf pockets supports prolific growth of the floating fern Azolla filiculoides. To evaluate contributions by further microorganisms, the A. filiculoides microbiome and nitrogen metabolism in bacteria persistently associated with Azolla ferns were characterized. A metagenomic approach was taken complemented by detection of N₂ O released and nitrogen isotope determinations of fern biomass. Ribosomal RNA genes in sequenced DNA of natural ferns, their enriched leaf pockets and water filtrate from the surrounding ditch established that bacteria of A. filiculoides differed entirely from surrounding water and revealed species of the order Rhizobiales. Analyses of seven cultivated Azolla species confirmed persistent association with Rhizobiales. Two distinct nearly full-length Rhizobiales genomes were identified in leaf-pocket-enriched samples from ditch grown A. filiculoides. Their annotation revealed genes for denitrification but not N₂ -fixation. ¹⁵ N₂ incorporation was active in ferns with N. azollae but not in ferns without. N₂ O was not detectably released from surface-sterilized ferns with the Rhizobiales. N₂ -fixing N. azollae, we conclude, dominated the microbiome of Azolla ferns. The persistent but less abundant heterotrophic Rhizobiales bacteria possibly contributed to lowering O₂ levels in leaf pockets but did not release detectable amounts of the strong greenhouse gas N₂ O.

MetaMeta: integrating metagenome analysis tools to improve taxonomic profiling

BACKGROUND

Many metagenome analysis tools are presently available to classify sequences and profile environmental samples. In particular, taxonomic profiling and binning methods are commonly used for such tasks. Tools available among these two categories make use of several techniques, e.g., read mapping, k-mer alignment, and composition analysis. Variations on the construction of the corresponding reference sequence databases are also common. In addition, different tools provide good results in different datasets and configurations. All this variation creates a complicated scenario to researchers to decide which methods to use. Installation, configuration and execution can also be difficult especially when dealing with multiple datasets and tools.

RESULTS

We propose MetaMeta: a pipeline to execute and integrate results from metagenome analysis tools. MetaMeta provides an easy workflow to run multiple tools with multiple samples, producing a single enhanced output profile for each sample. MetaMeta includes a database generation, pre-processing, execution, and integration steps, allowing easy execution and parallelization. The integration relies on the co-occurrence of organisms from different methods as the main feature to improve community profiling while accounting for differences in their databases.

CONCLUSIONS

In a controlled case with simulated and real data, we show that the integrated profiles of MetaMeta overcome the best single profile. Using the same input data, it provides more sensitive and reliable results with the presence of each organism being supported by several methods. MetaMeta uses Snakemake and has six pre-configured tools, all available at BioConda channel for easy installation (conda install -c bioconda metameta). The MetaMeta pipeline is open-source and can be downloaded at: https://gitlab.com/rki_bioinformatics .

Opening the treasure chest: A DNA-barcoding primer set for most higher taxa of Central European birds and mammals from museum collections

DNA-barcoding is a rapidly developing method for efficiently identifying samples to species level by means of short standard DNA sequences. However, reliable species assignment requires the availability of a comprehensive DNA barcode reference library, and hence numerous initiatives aim at generating such barcode databases for particular taxa or geographic regions. Historical museum collections represent a potentially invaluable source for the DNA-barcoding of many taxa. This is particularly true for birds and mammals, for which collecting fresh (voucher) material is often very difficult to (nearly) impossible due to the special animal welfare and conservation regulations that apply to vertebrates in general, and birds and mammals in particular. Moreover, even great efforts might not guarantee sufficiently complete sampling of fresh material in a short period of time. DNA extracted from historical samples is usually degraded, such that only short fragments can be amplified, rendering the recovery of the barcoding region as a single fragment impossible. Here, we present a new set of primers that allows the efficient amplification and sequencing of the entire barcoding region in most higher taxa of Central European birds and mammals in six overlapping fragments, thus greatly increasing the value of historical museum collections for generating DNA barcode reference libraries. Applying our new primer set in recently established NGS protocols promises to further increase the efficiency of barcoding old bird and mammal specimens.

Advances in the application of high-throughput sequencing in invertebrate virology

Over the last decade, advances in high-throughput sequencing technologies have revolutionised biological research, making it possible for DNA/RNA sequencing of any organism of interest to be undertaken. Sequencing approaches are now routinely used in the detection and characterisation of (novel) viruses, investigation of host-pathogen interactions, and effective development of disease treatment strategies. For the sequencing and identification of viruses of interest, metagenomics approaches using infected host tissue are frequently used, as it is not always possible to culture and isolate these pathogens. High-throughput sequencing can also be used to investigate host-pathogen interactions by investigating (temporal) transcriptomic responses of both the host and virus, potentially leading to the discovery of novel opportunities for treatment and drug targets. In addition, viruses in environmental samples (e.g. water or soil samples) can be identified using eDNA/metagenomics approaches. The promise that recent developments in sequencing brings to the field of invertebrate virology are not devoid of technical challenges, including the need for better laboratory and bioinformatics strategies to sequence and assemble virus genomes within complex tissue or environmental samples, and the difficulties associated with the annotation of the large number of novel viruses being discovered.

An Agile Functional Analysis of Metagenomic Data Using SUPER-FOCUS

One of the main goals in metagenomics is to identify the functional profile of a microbial community from unannotated shotgun sequencing reads. Functional annotation is important in biological research because it enables researchers to identify the abundance of functional genes of the organisms present in the sample, answering the question, "What can the organisms in the sample do?" Most currently available approaches do not scale with increasing data volumes, which is important because both the number and lengths of the reads provided by sequencing platforms keep increasing. Here, we present SUPER-FOCUS, SUbsystems Profile by databasE Reduction using FOCUS, an agile homology-based approach using a reduced reference database to report the subsystems present in metagenomic datasets and profile their abundances. SUPER-FOCUS was tested with real metagenomes, and the results show that it accurately predicts the subsystems present in the profiled microbial communities, is computationally efficient, and up to 1000 times faster than other tools. SUPER-FOCUS is freely available at http://edwards.sdsu.edu/SUPERFOCUS .

Do Vertebrate Gut Metagenomes Confer Rapid Ecological Adaptation?

During times of rapid environmental change, survival of most vertebrate populations depends on their phenomic plasticity. Although differential gene-expression and post-transcriptional processes of the host genome receive focus as the main molecular mechanisms, growing evidence points to the gut microbiota as a key driver defining hosts' phenotypes. We propose that the plasticity of the gut microbiota might be an essential factor determining phenomic plasticity of vertebrates, and that it might play a pivotal role when vertebrates acclimate and adapt to fast environmental variation. We contemplate some key questions and suggest methodological approaches and experimental designs that can be used to evaluate whether gut microorganisms provide a boost of plasticity to vertebrates' phenomes, thereby increasing their acclimation and adaptation capacity.

Metagenomics and Bioinformatics in Microbial Ecology: Current Status and Beyond

Metagenomic approaches are now commonly used in microbial ecology to study microbial communities in more detail, including many strains that cannot be cultivated in the laboratory. Bioinformatic analyses make it possible to mine huge metagenomic datasets and discover general patterns that govern microbial ecosystems. However, the findings of typical metagenomic and bioinformatic analyses still do not completely describe the ecology and evolution of microbes in their environments. Most analyses still depend on straightforward sequence similarity searches against reference databases. We herein review the current state of metagenomics and bioinformatics in microbial ecology and discuss future directions for the field. New techniques will allow us to go beyond routine analyses and broaden our knowledge of microbial ecosystems. We need to enrich reference databases, promote platforms that enable meta- or comprehensive analyses of diverse metagenomic datasets, devise methods that utilize long-read sequence information, and develop more powerful bioinformatic methods to analyze data from diverse perspectives.

A metabarcoding framework for facilitated survey of endolithic phototrophs with tufA

Background

In spite of their ecological importance as primary producers and microbioeroders of marine calcium carbonate (CaCO₃) substrata, endolithic phototrophs spanning both prokaryotic (the cyanobacteria) and eukaryotic algae lack established molecular resources for their facilitated survey with high throughput sequencing. Here, the development of a metabarcoding framework for the elongation factor EF-Ttu (tufA) was tested on four Illumina-sequenced marine CaCO₃ microfloras for the characterization of their endolithic phototrophs, especially the abundant bioeroding Ostreobium spp. (Ulvophyceae). The framework consists of novel tufA degenerate primers and a comprehensive database enabling Operational Taxonomic Unit (OTU) identification at multiple taxonomic ranks with percent identity thresholds determined herein.

Results

The newly established tufA database comprises 4057 non-redundant sequences (from 1339 eukaryotic and prokaryotic phototrophs, and 2718 prokaryotic heterotrophs) including 27 classes in 10 phyla of phototrophic diversity summarized from data mining on GenBank^®, our barcoding of >150 clones produced from coral reef microfloras, and >300 eukaryotic phototrophs (>230 Ulvophyceae including >100 ‘Ostreobium’ spp., and >70 Florideophyceae, Phaeophyceae and miscellaneous taxa). Illumina metabarcoding with the newly designed primers resulted in 802 robust OTUs including 618 phototrophs and 184 heterotrophs (77 and 23 % of OTUs, respectively). Phototrophic OTUs belonged to 14 classes of phototrophs found in seven phyla, and represented ~98 % of all reads. The phylogenetic profiles of coral reef microfloras showed few OTUs in large abundance (proportion of reads) for the Chlorophyta (Ulvophyceae, i.e. Ostreobium and Phaeophila), the Rhodophyta (Florideophyceae) and Haptophyta (Coccolithophyceae), and a large diversity (richness) of OTUs in lower abundance for the Cyanophyta (Cyanophyceae) and the Ochrophyta (the diatoms, ‘Bacillariophyta’). The bioerosive ‘Ostreobium’ spp. represented four families in a large clade of subordinal divergence, i.e. the Ostreobidineae, and a fifth, phylogenetically remote family in the suborder Halimedineae (provisionally assigned as the ‘Pseudostreobiaceae’). Together they harbor 85–95 delimited cryptic species of endolithic microsiphons.

Conclusions

The novel degenerate primers allowed for amplification of endolithic phototrophs across a wide phylogenetic breadth as well as their recovery in very large proportions of reads (overall 98 %) and diversity (overall 77 % of OTUs). The established companion tufA database and determined identity thresholds allow for OTU identification at multiple taxonomic ranks to facilitate the monitoring of phototrophic assemblages via metabarcoding, especially endolithic communities rich in bioeroding Ulvophyceae, such as those harboring ‘Ostreobium’ spp., Phaeophila spp. and associated algal diversity.

Electronic supplementary material

The online version of this article (doi:10.1186/s12898-016-0068-x) contains supplementary material, which is available to authorized users.

Analysis of Illumina MiSeq Metabarcoding Data: Application to Benthic Indices for Environmental Monitoring

This protocol details the analysis of Illumina MiSeq amplicon libraries derived from marine benthic macroinvertebrate samples and based on two barcodes of the mitochondrial cytochrome oxidase 1 (CO1) gene: a "short region," covered by overlapping forward and reverse reads and a "long region" for which forward and reverse reads do not overlap. Aside from providing guidelines for analyzing both types of amplicons, we show how amplicon reads can be used for the calculation of benthic indices for environmental monitoring.

Exploring the environmental diversity of kinetoplastid flagellates in the high-throughput DNA sequencing era

The class Kinetoplastea encompasses both free-living and parasitic species from a wide range of hosts. Several representatives of this group are responsible for severe human diseases and for economic losses in agriculture and livestock. While this group encompasses over 30 genera, most of the available information has been derived from the vertebrate pathogenic genera Leishmaniaand Trypanosoma. Recent studies of the previously neglected groups of Kinetoplastea indicated that the actual diversity is much higher than previously thought. This article discusses the known segment of kinetoplastid diversity and how gene-directed Sanger sequencing and next-generation sequencing methods can help to deepen our knowledge of these interesting protists.

SUPER-FOCUS: a tool for agile functional analysis of shotgun metagenomic data

SUMMARY

Analyzing the functional profile of a microbial community from unannotated shotgun sequencing reads is one of the important goals in metagenomics. Functional profiling has valuable applications in biological research because it identifies the abundances of the functional genes of the organisms present in the original sample, answering the question what they can do. Currently, available tools do not scale well with increasing data volumes, which is important because both the number and lengths of the reads produced by sequencing platforms keep increasing. Here, we introduce SUPER-FOCUS, SUbsystems Profile by databasE Reduction using FOCUS, an agile homology-based approach using a reduced reference database to report the subsystems present in metagenomic datasets and profile their abundances. SUPER-FOCUS was tested with over 70 real metagenomes, the results showing that it accurately predicts the subsystems present in the profiled microbial communities, and is up to 1000 times faster than other tools.

AVAILABILITY AND IMPLEMENTATION

SUPER-FOCUS was implemented in Python, and its source code and the tool website are freely available at https://edwards.sdsu.edu/SUPERFOCUS.

CONTACT

redwards@mail.sdsu.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.