Management of DNA reference libraries for barcoding and metabarcoding studies with the R package refdb

DNA barcoding and metabarcoding are revolutionizing the study and survey of biodiversity. In order to assign taxonomic labels to the DNA sequence data retrieved, these methods are strongly dependent on comprehensive and accurate reference databases. Producing reliable databases linking biological sequences and taxonomic data can be—and often has been—done using mainstream tools such as spreadsheet software. However, spreadsheets quickly become insufficient when the amount of data increases to thousands of taxa and sequences to be matched, and validation operations become more complex and are error prone if done in a manual way. Thus, there is a clear need for providing scientists with user-friendly, reliable and powerful tools to manipulate and manage DNA reference databases in tractable, sound and efficient ways. Here, we introduce the R package refdb as an environment for semi-automatic and assisted construction of DNA reference libraries. The refdb package is a reference database manager offering a set of powerful functions to import, organize, clean, filter, audit and export the data. It is broadly applicable in metabarcoding data generally obtained in biodiversity and biomonitoring studies. We present the main features of the package and outline how refdb can speed up reference database generation, management and handling, and thus contribute to standardization and repeatability in barcoding and metabarcoding studies.     

Discovering the role of Patagonian birds in the dispersal of truffles and other mycorrhizal fungi

1. Download : Download high-res image (172KB)
2. Download : Download full-size image

     

Efficient and sensitive identification and quantification of airborne pollen using next‐generation DNA sequencing

Pollen monitoring is an important and widely used tool in allergy research and creation of awareness in pollen‐allergic patients. Current pollen monitoring methods are microscope‐based, labour intensive and cannot identify pollen to the genus level in some relevant allergenic plant groups. Therefore, a more efficient, cost‐effective and sensitive method is needed. Here, we present a method for identification and quantification of airborne pollen using DNA sequencing. Pollen is collected from ambient air using standard techniques. DNA is extracted from the collected pollen, and a fragment of the chloroplast gene trnL is amplified using PCR. The PCR product is subsequently sequenced on a next‐generation sequencing platform (Ion Torrent). Amplicon molecules are sequenced individually, allowing identification of different sequences from a mixed sample. We show that this method provides an accurate qualitative and quantitative view of the species composition of samples of airborne pollen grains. We also show that it correctly identifies the individual grass genera present in a mixed sample of grass pollen, which cannot be achieved using microscopic pollen identification. We conclude that our method is more efficient and sensitive than current pollen monitoring techniques and therefore has the potential to increase the throughput of pollen monitoring.     

PhytoREF: a reference database of the plastidial 16S rRNA gene of photosynthetic eukaryotes with curated taxonomy

Photosynthetic eukaryotes have a critical role as the main producers in most ecosystems of the biosphere. The ongoing environmental metabarcoding revolution opens the perspective for holistic ecosystems biological studies of these organisms, in particular the unicellular microalgae that often lack distinctive morphological characters and have complex life cycles. To interpret environmental sequences, metabarcoding necessarily relies on taxonomically curated databases containing reference sequences of the targeted gene (or barcode) from identified organisms. To date, no such reference framework exists for photosynthetic eukaryotes. In this study, we built the PhytoREF database that contains 6490 plastidial 16S rDNA reference sequences that originate from a large diversity of eukaryotes representing all known major photosynthetic lineages. We compiled 3333 amplicon sequences available from public databases and 879 sequences extracted from plastidial genomes, and generated 411 novel sequences from cultured marine microalgal strains belonging to different eukaryotic lineages. A total of 1867 environmental Sanger 16S rDNA sequences were also included in the database. Stringent quality filtering and a phylogeny‐based taxonomic classification were applied for each 16S rDNA sequence. The database mainly focuses on marine microalgae, but sequences from land plants (representing half of the PhytoREF sequences) and freshwater taxa were also included to broaden the applicability of PhytoREF to different aquatic and terrestrial habitats. PhytoREF, accessible via a web interface ( http://phytoref.fr ), is a new resource in molecular ecology to foster the discovery, assessment and monitoring of the diversity of photosynthetic eukaryotes using high‐throughput sequencing.     

Tag jumps illuminated – reducing sequence‐to‐sample misidentifications in metabarcoding studies

Metabarcoding of environmental samples on second‐generation sequencing platforms has rapidly become a valuable tool for ecological studies. A fundamental assumption of this approach is the reliance on being able to track tagged amplicons back to the samples from which they originated. In this study, we address the problem of sequences in metabarcoding sequencing outputs with false combinations of used tags (tag jumps). Unless these sequences can be identified and excluded from downstream analyses, tag jumps creating sequences with false, but already used tag combinations, can cause incorrect assignment of sequences to samples and artificially inflate diversity. In this study, we document and investigate tag jumping in metabarcoding studies on Illumina sequencing platforms by amplifying mixed‐template extracts obtained from bat droppings and leech gut contents with tagged generic arthropod and mammal primers, respectively. We found that an average of 2.6% and 2.1% of sequences had tag combinations, which could be explained by tag jumping in the leech and bat diet study, respectively. We suggest that tag jumping can happen during blunt‐ending of pools of tagged amplicons during library build and as a consequence of chimera formation during bulk amplification of tagged amplicons during library index PCR. We argue that tag jumping and contamination between libraries represents a considerable challenge for Illumina‐based metabarcoding studies, and suggest measures to avoid false assignment of tag jumping‐derived sequences to samples.     

Toward accurate molecular identification of species in complex environmental samples: testing the performance of sequence filtering and clustering methods

Metabarcoding has the potential to become a rapid, sensitive, and effective approach for identifying species in complex environmental samples. Accurate molecular identification of species depends on the ability to generate operational taxonomic units (OTUs) that correspond to biological species. Due to the sometimes enormous estimates of biodiversity using this method, there is a great need to test the efficacy of data analysis methods used to derive OTUs. Here, we evaluate the performance of various methods for clustering length variable 18S amplicons from complex samples into OTUs using a mock community and a natural community of zooplankton species. We compare analytic procedures consisting of a combination of (1) stringent and relaxed data filtering, (2) singleton sequences included and removed, (3) three commonly used clustering algorithms (mothur, UCLUST, and UPARSE), and (4) three methods of treating alignment gaps when calculating sequence divergence. Depending on the combination of methods used, the number of OTUs varied by nearly two orders of magnitude for the mock community (60–5068 OTUs) and three orders of magnitude for the natural community (22–22191 OTUs). The use of relaxed filtering and the inclusion of singletons greatly inflated OTU numbers without increasing the ability to recover species. Our results also suggest that the method used to treat gaps when calculating sequence divergence can have a great impact on the number of OTUs. Our findings are particularly relevant to studies that cover taxonomically diverse species and employ markers such as rRNA genes in which length variation is extensive.     

Winogradsky columns as a strategy to study typically rare microbial eukaryotes

Winogradsky columns have been widely used to study soil microbial communities, but the vast majority of those investigations have focused on the ecology and diversity of bacteria. In contrast, microbial eukaryotes (ME) have been regularly overlooked in studies based on experimental soil columns. Despite the recognized ecological relevance of ME in soil communities, investigations focused on ME diversity and the abundance of certain groups of interest are still scarce. In the present study, we used DNA metabarcoding (high-throughput sequencing of the V4 region of the 18S rRNA locus) to survey the ME diversity and abundance in an experimental Winogradsky soil column. Consistent with previous surveys in natural soils, our survey identified members of Cercozoa (Rhizaria; 31.2%), Apicomplexa and Ciliophora (Alveolata; 12.5%) as the predominant ME groups, but at particular depths we also detected the abundant presence of ME lineages that are typically rare in natural environments, such as members of the Vampyrellida (Rhizaria) and Breviatea (Amorphea). Our survey demonstrates that experimental soil columns are an efficient enrichment-culture approach that can enhance investigations about the diversity and ecology of ME in soils.     

Assessment of littoral algal diversity from the northern Gulf of Mexico using environmental DNA metabarcoding

Traditional methods for algal biodiversity monitoring are costly and time inefficient because they rely on high‐level taxonomic expertise to address species identity problems involving phenotypic plasticity and morphological convergence. These problems are exacerbated in regions such as the Gulf of Mexico, that has a limited history of phycological exploration, but that are economically important or threatened by numerous anthropogenic stressors. Given the high pace of disturbance to natural systems, there is a critical need for expedient and cost‐effective tools for the study of benthic algal communities. Here we document the use of environmental DNA metabarcoding, using the partial LSU rDNA and 23S rDNA plastid molecular markers, to elucidate littoral algal diversity in the Northern Gulf of Mexico. We assigned 73.7% of algal OTUs to genus and 59.6% to species ranks. Our current study detected molecular signals for 35 algal/protist species with no previous reports in the Gulf of Mexico, thus providing an important, molecular‐validated, baseline of species richness for this region. We also make several bioinformatic recommendations for the efficient use of high‐throughput sequence data to assess biological communities.     

DNA metabarcoding provides insights into seasonal diet variations in Chinese mole shrew (Anourosorex squamipes) with potential implications for evaluating crop impacts

Diet analysis of potential small mammals pest species is important for understanding feeding ecology and evaluating their impact on crops and stored foods. Chinese mole shrew (Anourosorex squamipes), distributed in Southwest China, has previously been reported as a farmland pest. Effective population management of this species requires a better understanding of its diet, which can be difficult to determine with high taxonomic resolution using conventional microhistological methods. In this study, we used two DNA metabarcoding assays to identify 38 animal species and 65 plant genera from shrew stomach contents, which suggest that A. squamipes is an omnivorous generalist. Earthworms are the most prevalent (>90%) and abundant (>80%) food items in the diverse diet of A. squamipes. Species of the Fabaceae (frequency of occurrence [FO]: 88%; such as peanuts) and Poaceae (FO: 71%; such as rice) families were the most common plant foods identified in the diet of A. squamipes. Additionally, we found a seasonal decrease in the diversity and abundance of invertebrate foods from spring and summer to winter. Chinese mole shrew has a diverse and flexible diet throughout the year to adapt to seasonal variations in food availability, contributing to its survival even when food resources are limited. This study provides a higher resolution identification of the diet of A. squamipes than has been previously described and is valuable for understanding shrew feeding ecology as well as evaluating possible species impacts on crops.     

DNA metabarcoding of feces to infer summer diet of Pacific walruses

Environmental conditions in the Chukchi Sea are changing rapidly and may alter the abundance and distribution of marine species and their benthic prey. We used a metabarcoding approach to identify potentially important prey taxa from Pacific walrus (Odobenus rosmarus divergens) fecal samples (n = 87). Bivalvia was the most dominant class of prey (66% of all normalized counts) and occurred in 98% of the samples. Polychaeta and Gastropoda occurred in 70% and 62% of the samples, respectively. The remaining nine invertebrate classes comprised <21% of all normalized counts. The common occurrence of these three prey classes is consistent with examinations of walrus stomach contents. Despite these consistencies, biases in the metabarcoding approach to determine diet from feces have been highlighted in other studies and require further study, in addition to biases that may have arisen from our opportunistic sampling. However, this noninvasive approach provides accurate identification of prey taxa from degraded samples and could yield much-needed information on shifts in walrus diet in a rapidly changing Arctic.