Environmental barcoding: a next-generation sequencing approach for biomonitoring applications using river benthos

Timely and accurate biodiversity analysis poses an ongoing challenge for the success of biomonitoring programs. Morphology-based identification of bioindicator taxa is time consuming, and rarely supports species-level resolution especially for immature life stages. Much work has been done in the past decade to develop alternative approaches for biodiversity analysis using DNA sequence-based approaches such as molecular phylogenetics and DNA barcoding. On-going assembly of DNA barcode reference libraries will provide the basis for a DNA-based identification system. The use of recently introduced next-generation sequencing (NGS) approaches in biodiversity science has the potential to further extend the application of DNA information for routine biomonitoring applications to an unprecedented scale. Here we demonstrate the feasibility of using 454 massively parallel pyrosequencing for species-level analysis of freshwater benthic macroinvertebrate taxa commonly used for biomonitoring. We designed our experiments in order to directly compare morphology-based, Sanger sequencing DNA barcoding, and next-generation environmental barcoding approaches. Our results show the ability of 454 pyrosequencing of mini-barcodes to accurately identify all species with more than 1% abundance in the pooled mixture. Although the approach failed to identify 6 rare species in the mixture, the presence of sequences from 9 species that were not represented by individuals in the mixture provides evidence that DNA based analysis may yet provide a valuable approach in finding rare species in bulk environmental samples. We further demonstrate the application of the environmental barcoding approach by comparing benthic macroinvertebrates from an urban region to those obtained from a conservation area. Although considerable effort will be required to robustly optimize NGS tools to identify species from bulk environmental samples, our results indicate the potential of an environmental barcoding approach for biomonitoring programs.     

Method for quick DNA barcode reference library construction

DNA barcoding has become one of the most important techniques in plant species identification. Successful application of this technology is dependent on the availability of reference database of high species coverage. Unfortunately, there are experimental and data processing challenges to construct such a library within a short time. Here, we present our solutions to these challenges. We sequenced six conventional DNA barcode fragments (ITS1, ITS2, matK1, matK2, rbcL1, and rbcL2) of 380 flowering plants on next‐generation sequencing (NGS) platforms (Illumina Hiseq 2500 and Ion Torrent S5) and the Sanger sequencing platform. After comparing the sequencing depths, read lengths, base qualities, and base accuracies, we conclude that Illumina Hiseq2500 PE250 run is suitable for conventional DNA barcoding. We developed a new “Cotu” method to create consensus sequences from NGS reads for longer output sequences and more reliable bases than the other three methods. Step‐by‐step instructions to our method are provided. By using high‐throughput machines (PCR and NGS), labeling PCR, and the Cotu method, it is possible to significantly reduce the cost and labor investments for DNA barcoding. A regional or even global DNA barcoding reference library with high species coverage is likely to be constructed in a few years.     

When to use next generation sequencing or diagnostic PCR in diet analyses

Next‐generation sequencing (NGS) is increasingly used for diet analyses; however, it may not always describe diet samples well. A reason for this is that diet samples contain mixtures of food DNA in different amounts as well as consumer DNA which can reduce the food DNA characterized. Because of this, detections will depend on the relative amount and identity of each type of DNA. For such samples, diagnostic PCR will most likely give more reliable results, as detection probability is only marginally dependent on other copresent DNA. We investigated the reliability of each method to test (a) whether predatory beetle regurgitates, supposed to be low in consumer DNA, allow to retrieve prey sequences using general barcoding primers that co‐amplify the consumer DNA, and (b) to assess the sequencing depth or replication needed for NGS and diagnostic PCR to give stable results. When consumer DNA is co‐amplified, NGS is better suited to discover the range of possible prey, than for comparing co‐occurrences of diet species between samples, as retested samples were repeatedly different in prey detections with this approach. This shows that samples were incompletely described, as prey detected by diagnostic PCR frequently were missed by NGS. As the sequencing depth needed to reliably describe the diet in such samples becomes very high, the cost‐efficiency and reliability of diagnostic PCR make diagnostic PCR better suited for testing large sample‐sets. Especially if the targeted prey taxa are thought to be of ecological importance, as diagnostic PCR gave more nested and consistent results in repeated testing of the same sample.     

Next‐generation sequencing to inventory taxonomic diversity in eukaryotic communities: a test for freshwater diatoms

The recent emergence of barcoding approaches coupled to those of next‐generation sequencing (NGS) has raised new perspectives for studying environmental communities. In this framework, we tested the possibility to derive accurate inventories of diatom communities from pyrosequencing outputs with an available DNA reference library. We used three molecular markers targeting the nuclear, chloroplast and mitochondrial genomes (SSU rDNA, rbcL and cox1) and three samples of a mock community composed of 30 known diatom strains belonging to 21 species. In the goal to detect methodological biases, one sample was constituted directly from pooled cultures, whereas the others consisted of pooled PCR products. The NGS reads obtained by pyrosequencing (Roche 454) were compared first to a DNA reference library including the sequences of all the species used to constitute the mock community, and second to a complete DNA reference library with a larger taxonomic coverage. A stringent taxonomic assignation gave inventories that were compared to the real one. We detected biases due to DNA extraction and PCR amplification that resulted in false‐negative detection. Conversely, pyrosequencing errors appeared to generate false positives, especially in case of closely allied species. The taxonomic coverage of DNA reference libraries appears to be the most crucial factor, together with marker polymorphism which is essential to identify taxa at the species level. RbcL offers a high resolving power together with a large DNA reference library. Although needing further optimization, pyrosequencing is suitable for identifying diatom assemblages and may find applications in the field of freshwater biomonitoring.     

DNA barcodes from century‐old type specimens using next‐generation sequencing

Type specimens have high scientific importance because they provide the only certain connection between the application of a Linnean name and a physical specimen. Many other individuals may have been identified as a particular species, but their linkage to the taxon concept is inferential. Because type specimens are often more than a century old and have experienced conditions unfavourable for DNA preservation, success in sequence recovery has been uncertain. This study addresses this challenge by employing next‐generation sequencing (NGS) to recover sequences for the barcode region of the cytochrome c oxidase 1 gene from small amounts of template DNA. DNA quality was first screened in more than 1800 century‐old type specimens of Lepidoptera by attempting to recover 164‐bp and 94‐bp reads via Sanger sequencing. This analysis permitted the assignment of each specimen to one of three DNA quality categories – high (164‐bp sequence), medium (94‐bp sequence) or low (no sequence). Ten specimens from each category were subsequently analysed via a PCR‐based NGS protocol requiring very little template DNA. It recovered sequence information from all specimens with average read lengths ranging from 458 bp to 610 bp for the three DNA categories. By sequencing ten specimens in each NGS run, costs were similar to Sanger analysis. Future increases in the number of specimens processed in each run promise substantial reductions in cost, making it possible to anticipate a future where barcode sequences are available from most type specimens.     

木霉属补充DNA条形码筛选

木霉属真菌是一类重要的生物资源,在工农业、环境保护等方面具有较高经济价值,对其进行快速、准确的物种鉴定兼具理论意义和应用前景。以木霉属35个概念清晰的种为材料,选择ITS、rpb2和 tef1作为候选基因序列,利用TaxonGap对231个序列片段进行分析,将种内与种间序列差异以及序列获取难易程度作为评价指标,筛选该属的补充条形码片段。结果表明,rpb2具有适宜的种内与种间序列差异,其最小的种间差异（2.48%）,大于最大种内差异（1.8%）,种内、种间遗传距离存在明显的间隔区,并且该基因序列具有较高的PCR扩增与测序成功率（94.4%）;ITS和tef1基因序列的种内与种间序列差异之间存在交叉重叠。因此建议rpb2作为木霉属的补充DNA条形码序列。     

High-throughput multilocus sequence typing: bringing molecular typing to the next level

Multilocus sequence typing (MLST) is a widely used system for typing microorganisms by sequence analysis of housekeeping genes. The main advantage of MLST in comparison to other typing techniques is the unambiguity and transferability of sequence data. However, a main disadvantage is the high cost of DNA sequencing. Here we introduce a high-throughput MLST (HiMLST) method that employs next-generation sequencing (NGS) technology (Roche 454), to generate large quantities of high-quality MLST data at low costs. The HiMLST protocol consists of two steps. In the first step MLST target genes are amplified by PCR in multi-well plates. During this PCR the amplicons of each bacterial isolate are provided with a unique DNA barcode, the multiplex identifier (MID). In the second step all amplicons are pooled and sequenced in a single NGS-run. The MLST profile of each individual isolate can be retrieved easily using its unique MID. With HiMLST we have profiled 575 isolates of Legionella pneumophila, Staphylococcus aureus, Pseudomonas aeruginosa and Streptococcus pneumoniae in mixed species HiMLST experiments. In conclusion, the introduction of HiMLST paves the way for a broad employment of the MLST as a high-quality and cost-effective method for typing microbial species.     

Plant DNA barcoding: from gene to genome

DNA barcoding is currently a widely used and effective tool that enables rapid and accurate identification of plant species; however, none of the available loci work across all species. Because single‐locus DNA barcodes lack adequate variations in closely related taxa, recent barcoding studies have placed high emphasis on the use of whole‐chloroplast genome sequences which are now more readily available as a consequence of improving sequencing technologies. While chloroplast genome sequencing can already deliver a reliable barcode for accurate plant identification it is not yet resource‐effective and does not yet offer the speed of analysis provided by single‐locus barcodes to unspecialized laboratory facilities. Here, we review the development of candidate barcodes and discuss the feasibility of using the chloroplast genome as a super‐barcode. We advocate a new approach for DNA barcoding that, for selected groups of taxa, combines the best use of single‐locus barcodes and super‐barcodes for efficient plant identification. Specific barcodes might enhance our ability to distinguish closely related plants at the species and population levels.     

Taxon-specific PCR for DNA barcoding arthropod prey in bat faeces

The application of DNA barcoding to dietary studies allows prey taxa to be identified in the absence of morphological evidence and permits a greater resolution of prey identity than is possible through direct examination of faecal material. For insectivorous bats, which typically eat a great diversity of prey and which chew and digest their prey thoroughly, DNA-based approaches to diet analysis may provide the only means of assessing the range and diversity of prey within faeces. Here, we investigated the effectiveness of DNA barcoding in determining the diets of bat species that specialize in eating different taxa of arthropod prey. We designed and tested a novel taxon-specific primer set and examined the performance of short barcode sequences in resolving prey species. We recovered prey DNA from all faecal samples and subsequent cloning and sequencing of PCR products, followed by a comparison of sequences to a reference database, provided species-level identifications for 149/207 (72%) clones. We detected a phylogenetically broad range of prey while completely avoiding detection of nontarget groups. In total, 37 unique prey taxa were identified from 15 faecal samples. A comparison of DNA data with parallel morphological analyses revealed a close correlation between the two methods. However, the sensitivity and taxonomic resolution of the DNA method were far superior. The methodology developed here provides new opportunities for the study of bat diets and will be of great benefit to the conservation of these ecologically important predators.     

Effectiveness of Annealing Blocking Primers versus Restriction Enzymes for Characterization of Generalist Diets: Unexpected Prey Revealed in the Gut Contents of Two Coral Reef Fish Species

Characterization of predator-prey interactions is challenging as researchers have to rely on indirect methods that can be costly, biased and too imprecise to elucidate the complexity of food webs. DNA amplification and sequencing techniques of gut and fecal contents are promising approaches, but their success largely depends on the ability to amplify the taxonomic array of prey consumed and then match prey amplicons with reference sequences. When little a priori information on diet is available or a generalist predator is targeted, versatile primer sets (also referred to as universal or general primers) as opposed to group- or species-specific primer sets are the most powerful to unveil the full range of prey consumed. However, versatile primers are likely to preferentially amplify the predominant, less degraded predator DNA if no manipulation is performed to exclude this confounding DNA template. In this study we compare two approaches that eliminate the confounding predator template: restriction digestion and the use of annealing blocking primers. First, we use a preliminary DNA barcode library provided by the Moorea BIOCODE project to 1) evaluate the cutting frequency of commercially available restriction enzymes and 2) design predator specific annealing blocking primers. We then compare the performance of the two predator removal strategies for the detection of prey templates using two versatile primer sets from the gut contents of two generalist coral reef fish species sampled in Moorea. Our study demonstrates that blocking primers should be preferentially used over restriction digestion for predator DNA removal as they recover greater prey diversity. We also emphasize that a combination of versatile primers may be required to best represent the breadth of a generalist''s diet.     

From barcodes to genomes: extending the concept of DNA barcoding

DNA barcoding has had a major impact on biodiversity science. The elegant simplicity of establishing massive scale databases for a few barcode loci is continuing to change our understanding of species diversity patterns, and continues to enhance human abilities to distinguish among species. Capitalizing on the developments of next generation sequencing technologies and decreasing costs of genome sequencing, there is now the opportunity for the DNA barcoding concept to be extended to new kinds of genomic data. We illustrate the benefits and capacity to do this, and also note the constraints and barriers to overcome before it is truly scalable. We advocate a twin track approach: (i) continuation and acceleration of global efforts to build the DNA barcode reference library of life on earth using standard DNA barcodes and (ii) active development and application of extended DNA barcodes using genome skimming to augment the standard barcoding approach.     

PCR-Free Enrichment of Mitochondrial DNA from Human Blood and Cell Lines for High Quality Next-Generation DNA Sequencing

Recent advances in sequencing technology allow for accurate detection of mitochondrial sequence variants, even those in low abundance at heteroplasmic sites. Considerable sequencing cost savings can be achieved by enriching samples for mitochondrial (relative to nuclear) DNA. Reduction in nuclear DNA (nDNA) content can also help to avoid false positive variants resulting from nuclear mitochondrial sequences (numts). We isolate intact mitochondrial organelles from both human cell lines and blood components using two separate methods: a magnetic bead binding protocol and differential centrifugation. DNA is extracted and further enriched for mitochondrial DNA (mtDNA) by an enzyme digest. Only 1 ng of the purified DNA is necessary for library preparation and next generation sequence (NGS) analysis. Enrichment methods are assessed and compared using mtDNA (versus nDNA) content as a metric, measured by using real-time quantitative PCR and NGS read analysis. Among the various strategies examined, the optimal is differential centrifugation isolation followed by exonuclease digest. This strategy yields >35% mtDNA reads in blood and cell lines, which corresponds to hundreds-fold enrichment over baseline. The strategy also avoids false variant calls that, as we show, can be induced by the long-range PCR approaches that are the current standard in enrichment procedures. This optimization procedure allows mtDNA enrichment for efficient and accurate massively parallel sequencing, enabling NGS from samples with small amounts of starting material. This will decrease costs by increasing the number of samples that may be multiplexed, ultimately facilitating efforts to better understand mitochondria-related diseases.     

Counting with DNA in metabarcoding studies: How should we convert sequence reads to dietary data?

Advances in DNA sequencing technology have revolutionized the field of molecular analysis of trophic interactions, and it is now possible to recover counts of food DNA sequences from a wide range of dietary samples. But what do these counts mean? To obtain an accurate estimate of a consumer's diet should we work strictly with data sets summarizing frequency of occurrence of different food taxa, or is it possible to use relative number of sequences? Both approaches are applied to obtain semi‐quantitative diet summaries, but occurrence data are often promoted as a more conservative and reliable option due to taxa‐specific biases in recovery of sequences. We explore representative dietary metabarcoding data sets and point out that diet summaries based on occurrence data often overestimate the importance of food consumed in small quantities (potentially including low‐level contaminants) and are sensitive to the count threshold used to define an occurrence. Our simulations indicate that using relative read abundance (RRA) information often provides a more accurate view of population‐level diet even with moderate recovery biases incorporated; however, RRA summaries are sensitive to recovery biases impacting common diet taxa. Both approaches are more accurate when the mean number of food taxa in samples is small. The ideas presented here highlight the need to consider all sources of bias and to justify the methods used to interpret count data in dietary metabarcoding studies. We encourage researchers to continue addressing methodological challenges and acknowledge unanswered questions to help spur future investigations in this rapidly developing area of research.     

Detecting invertebrate species in archived collections using next‐generation sequencing

Invertebrate biodiversity measured at mostly family level is widely used in biological monitoring programmes to assess anthropogenic impacts on ecosystems. However, next‐generation sequencing (NGS) could allow development of new more sensitive biomonitoring tools by allowing rapid species identification. This could be accelerated if archived invertebrate collections and environmental information from past programmes are used to understand species distributions and their environmental responses. In this study, we take archived macroinvertebrate samples from two sites collected on multiple occasions and test whether NGS can successfully detect species. Samples had been stored in 70% ethanol at room temperature for up to 12 years. Three amplicons ranging from 197 to 274 bps within the DNA barcode region were amplified from samples and compared to DNA barcoding libraries to identify species. We were able to amplify partial DNA barcodes from most samples, and species were often detected with multiple amplicons. However, some singletons and taxa poorly covered by DNA barcoding were missed. This suggests additional DNA barcodes will be required to fill ‘gaps’ in current DNA barcode libraries for aquatic macroinvertebrates and/or that it may not be possible to detect all taxa in a sample. Furthermore, older samples often detected fewer taxa and were less reliable for amplification, suggesting NGS is best used on samples within 8 years of collection. Nevertheless, many common taxa with existing DNA barcodes were reliably identified with NGS and were often present at sites across multiple years, showing the potential of NGS for detecting common and abundant species in archived material.     

Deep Sequencing to Identify the Causes of Viral Encephalitis

Deep sequencing allows for a rapid, accurate characterization of microbial DNA and RNA sequences in many types of samples. Deep sequencing (also called next generation sequencing or NGS) is being developed to assist with the diagnosis of a wide variety of infectious diseases. In this study, seven frozen brain samples from deceased subjects with recent encephalitis were investigated. RNA from each sample was extracted, randomly reverse transcribed and sequenced. The sequence analysis was performed in a blinded fashion and confirmed with pathogen-specific PCR. This analysis successfully identified measles virus sequences in two brain samples and herpes simplex virus type-1 sequences in three brain samples. No pathogen was identified in the other two brain specimens. These results were concordant with pathogen-specific PCR and partially concordant with prior neuropathological examinations, demonstrating that deep sequencing can accurately identify viral infections in frozen brain tissue.     

Messy eaters: Swabbing prey DNA from the exterior of inconspicuous predators when foraging cannot be observed

Complex coevolutionary relationships among competitors, predators, and prey have shaped taxa diversity, life history strategies, and even the avian migratory patterns we see today. Consequently, accurate documentation of prey selection is often critical for understanding these ecological and evolutionary processes. Conventional diet study methods lack the ability to document the diet of inconspicuous or difficult‐to‐study predators, such as those with large home ranges and those that move vast distances over short amounts of time, leaving gaps in our knowledge of trophic interactions in many systems. Migratory raptors represent one such group of predators where detailed diet studies have been logistically challenging. To address knowledge gaps in the foraging ecology of migrant raptors and provide a broadly applicable tool for the study of enigmatic predators, we developed a minimally invasive method to collect dietary information by swabbing beaks and talons of raptors to collect trace prey DNA. Using previously published COI primers, we were able to isolate and reference gene sequences in an open‐access barcode database to identify prey to species. This method creates a novel avenue to use trace molecular evidence to study prey selection of migrating raptors and will ultimately lead to a better understanding of raptor migration ecology. In addition, this technique has broad applicability and can be used with any wildlife species where even trace amounts of prey debris remain on the exterior of the predator after feeding.     

Pyrosequencing for mini-barcoding of fresh and old museum specimens

DNA barcoding is an effective approach for species identification and for discovery of new and/or cryptic species. Sanger sequencing technology is the method of choice for obtaining standard 650 bp cytochrome c oxidase subunit I (COI) barcodes. However, DNA degradation/fragmentation makes it difficult to obtain a full-length barcode from old specimens. Mini-barcodes of 130 bp from the standard barcode region have been shown to be effective for accurate identification in many animal groups and may be readily obtained from museum samples. Here we demonstrate the application of an alternative sequencing technology, the four-enzymes single-specimen pyrosequencing, in rapid, cost-effective mini-barcode analysis. We were able to generate sequences of up to 100 bp from mini-barcode fragments of COI in 135 fresh and 50 old Lepidoptera specimens (ranging from 53-97 year-old). The sequences obtained using pyrosequencing were of high quality and we were able to robustly match all the tested pyro-sequenced samples to their respective Sanger-sequenced standard barcode sequences, where available. Simplicity of the protocol and instrumentation coupled with higher speed and lower cost per sequence than Sanger sequencing makes this approach potentially useful in efforts to link standard barcode sequences from unidentified specimens to known museum specimens with only short DNA fragments.     

Multilocus sequence typing of Salmonella strains by high-throughput sequencing of selectively amplified target genes

Rapid development of next generation sequencing (NGS) technologies in recent years has made whole genome sequencing of bacterial genomes widely accessible. However, it is often unnecessary or not feasible to sequence the whole genome for most applications of genetic analyses in bacteria. Selectively capturing defined genomic regions followed by NGS analysis could be a promising approach for high-resolution molecular typing of a large set of strains. In this study, we describe a novel and straightforward PCR-based target-capturing method, hairpin-primed multiplex amplification (HPMA), which allows for simultaneous amplification of numerous target genes. To test the feasibility of NGS-based strain typing using HPMA, 20 target gene sequences were simultaneously amplified with barcode tagging in each of 41 Salmonella strains. The amplicons were then pooled and analyzed by 454 pyrosequencing. Analysis of the sequence data, as an extension of multilocus sequence typing (MLST), demonstrated the utility and potential of this novel typing method, MLST-seq, as a high-resolution strain typing method. With the rapidly increasing sequencing capacity of NGS, MLST-seq or its variations using different target enrichment methods can be expected to become a high-resolution typing method in the near future for high-throughput analysis of a large collection of bacterial strains.     

Integration of Wet and Dry Bench Processes Optimizes Targeted Next-generation Sequencing of Low-quality and Low-quantity Tumor Biopsies

All next-generation sequencing (NGS) procedures include assays performed at the laboratory bench ("wet bench") and data analyses conducted using bioinformatics pipelines ("dry bench"). Both elements are essential to produce accurate and reliable results, which are particularly critical for clinical laboratories. Targeted NGS technologies have increasingly found favor in oncology applications to help advance precision medicine objectives, yet the methods often involve disconnected and variable wet and dry bench workflows and uncoordinated reagent sets. In this report, we describe a method for sequencing challenging cancer specimens with a 21-gene panel as an example of a comprehensive targeted NGS system. The system integrates functional DNA quantification and qualification, single-tube multiplexed PCR enrichment, and library purification and normalization using analytically-verified, single-source reagents with a standalone bioinformatics suite. As a result, accurate variant calls from low-quality and low-quantity formalin-fixed, paraffin-embedded (FFPE) and fine-needle aspiration (FNA) tumor biopsies can be achieved. The method can routinely assess cancer-associated variants from an input of 400 amplifiable DNA copies, and is modular in design to accommodate new gene content. Two different types of analytically-defined controls provide quality assurance and help safeguard call accuracy with clinically-relevant samples. A flexible "tag" PCR step embeds platform-specific adaptors and index codes to allow sample barcoding and compatibility with common benchtop NGS instruments. Importantly, the protocol is streamlined and can produce 24 sequence-ready libraries in a single day. Finally, the approach links wet and dry bench processes by incorporating pre-analytical sample quality control results directly into the variant calling algorithms to improve mutation detection accuracy and differentiate false-negative and indeterminate calls. This targeted NGS method uses advances in both wetware and software to achieve high-depth, multiplexed sequencing and sensitive analysis of heterogeneous cancer samples for diagnostic applications.     

SinoDuplex: An Improved Duplex Sequencing Approach to Detect Low-frequency Variants in Plasma cfDNA Samples

Accurate detection of low frequency mutations from plasma cell-free DNA in blood using targeted next generation sequencing technology has shown promising benefits in clinical settings. Duplex sequencing technology is the most commonly used approach in liquid biopsies. Unique molecular identifiers are attached to each double-stranded DNA template, followed by production of low-error consensus sequences to detect low frequency variants. However, high sequencing costs have hindered application of this approach in clinical practice. Here, we have developed an improved duplex sequencing approach called Sino Duplex, which utilizes a pool of adapters containing pre-defined barcode sequences to generate far fewer barcode combinations than with random sequences, and implemented a novel computational analysis algorithm to generate duplex consensus sequences more precisely. Sino Duplex increased the output of duplex sequencing technology, making it more cost-effective. We evaluated our approach using reference standard samples and cell-free DNA samples from lung cancer patients. Our results showed that Sino Duplex has high sensitivity and specificity in detecting very low allele frequency mutations. The source code for Sino Duplex is freely available at https://github.com/Sin Oncology/sinoduplex.