Vector soup: high‐throughput identification of Neotropical phlebotomine sand flies using metabarcoding

Phlebotomine sand flies are haematophagous dipterans of primary medical importance. They represent the only proven vectors of leishmaniasis worldwide and are involved in the transmission of various other pathogens. Studying the ecology of sand flies is crucial to understand the epidemiology of leishmaniasis and further control this disease. A major limitation in this regard is that traditional morphological‐based methods for sand fly species identifications are time‐consuming and require taxonomic expertise. DNA metabarcoding holds great promise in overcoming this issue by allowing the identification of multiple species from a single bulk sample. Here, we assessed the reliability of a short insect metabarcode located in the mitochondrial 16S rRNA for the identification of Neotropical sand flies, and constructed a reference database for 40 species found in French Guiana. Then, we conducted a metabarcoding experiment on sand flies mixtures of known content and showed that the method allows an accurate identification of specimens in pools. Finally, we applied metabarcoding to field samples caught in a 1‐ha forest plot in French Guiana. Besides providing reliable molecular data for species‐level assignations of phlebotomine sand flies, our study proves the efficiency of metabarcoding based on the mitochondrial 16S rRNA for studying sand fly diversity from bulk samples. The application of this high‐throughput identification procedure to field samples can provide great opportunities for vector monitoring and eco‐epidemiological studies.     

Comparing the performance of three ancient DNA extraction methods for high‐throughput sequencing

The DNA molecules that can be extracted from archaeological and palaeontological remains are often degraded and massively contaminated with environmental microbial material. This reduces the efficacy of shotgun approaches for sequencing ancient genomes, despite the decreasing sequencing costs of high‐throughput sequencing (HTS). Improving the recovery of endogenous molecules from the DNA extraction and purification steps could, thus, help advance the characterization of ancient genomes. Here, we apply the three most commonly used DNA extraction methods to five ancient bone samples spanning a ~30 thousand year temporal range and originating from a diversity of environments, from South America to Alaska. We show that methods based on the purification of DNA fragments using silica columns are more advantageous than in solution methods and increase not only the total amount of DNA molecules retrieved but also the relative importance of endogenous DNA fragments and their molecular diversity. Therefore, these methods provide a cost‐effective solution for downstream applications, including DNA sequencing on HTS platforms.     

Not all are free‐living: high‐throughput DNA metabarcoding reveals a diverse community of protists parasitizing soil metazoa

Protists, the most diverse eukaryotes, are largely considered to be free‐living bacterivores, but vast numbers of taxa are known to parasitize plants or animals. High‐throughput sequencing (HTS) approaches now commonly replace cultivation‐based approaches in studying soil protists, but insights into common biases associated with this method are limited to aquatic taxa and samples. We created a mock community of common free‐living soil protists (amoebae, flagellates, ciliates), extracted DNA and amplified it in the presence of metazoan DNA using 454 HTS. We aimed at evaluating whether HTS quantitatively reveals true relative abundances of soil protists and at investigating whether the expected protist community structure is altered by the co‐amplification of metazoan‐associated protist taxa. Indeed, HTS revealed fundamentally different protist communities from those expected. Ciliate sequences were highly over‐represented, while those of most amoebae and flagellates were under‐represented or totally absent. These results underpin the biases introduced by HTS that prevent reliable quantitative estimations of free‐living protist communities. Furthermore, we detected a wide range of nonadded protist taxa probably introduced along with metazoan DNA, which altered the protist community structure. Among those, 20 taxa most closely resembled parasitic, often pathogenic taxa. Therewith, we provide the first HTS data in support of classical observational studies that showed that potential protist parasites are hosted by soil metazoa. Taken together, profound differences in amplification success between protist taxa and an inevitable co‐extraction of protist taxa parasitizing soil metazoa obscure the true diversity of free‐living soil protist communities.     

Universal and blocking primer mismatches limit the use of high‐throughput DNA sequencing for the quantitative metabarcoding of arthropods

The quantification of the biological diversity in environmental samples using high‐throughput DNA sequencing is hindered by the PCR bias caused by variable primer–template mismatches of the individual species. In some dietary studies, there is the added problem that samples are enriched with predator DNA, so often a predator‐specific blocking oligonucleotide is used to alleviate the problem. However, specific blocking oligonucleotides could coblock nontarget species to some degree. Here, we accurately estimate the extent of the PCR biases induced by universal and blocking primers on a mock community prepared with DNA of twelve species of terrestrial arthropods. We also compare universal and blocking primer biases with those induced by variable annealing temperature and number of PCR cycles. The results show that reads of all species were recovered after PCR enrichment at our control conditions (no blocking oligonucleotide, 45 °C annealing temperature and 40 cycles) and high‐throughput sequencing. They also show that the four factors considered biased the final proportions of the species to some degree. Among these factors, the number of primer–template mismatches of each species had a disproportionate effect (up to five orders of magnitude) on the amplification efficiency. In particular, the number of primer–template mismatches explained most of the variation (~3/4) in the amplification efficiency of the species. The effect of blocking oligonucleotide concentration on nontarget species relative abundance was also significant, but less important (below one order of magnitude). Considering the results reported here, the quantitative potential of the technique is limited, and only qualitative results (the species list) are reliable, at least when targeting the barcoding COI region.     

Advances in clone selection using high‐throughput bioreactors

Effective clone selection is a crucial step toward developing a robust mammalian cell culture production platform. Currently, clone selection is done by culturing cells in well plates and picking the highest producers. Ideally, clone selection should be done in a stirred tank bioreactor as this would best replicate the eventual production environment. The actual number of clones selected for future evaluation in bioreactors at bench‐scale is limited by the scale‐up and operational costs involved. This study describes the application of miniaturized stirred high‐throughput bioreactors (35 mL working volume; HTBRs) with noninvasive optical sensors for clone screening and selection. We investigated a method for testing several subclones simultaneously in a stirred environment using our high throughput bioreactors (up to 12 clones per HTBR run) and compared it with a traditional well plate selection approach. Importantly, it was found that selecting clones solely based on results from stationary well plate cultures could result in the chance of missing higher producing clones. Our approach suggests that choosing a clone after analyzing its performance in a stirred bioreactor environment is an improved method for clone selection. © 2010 American Institute of Chemical Engineers Biotechnol. Prog., 2010     

A comparison of DNA extraction methods for high‐throughput DNA analyses

The inclusion of next‐generation sequencing technologies in population genetic and phylogenetic studies has elevated the need to balance time and cost of DNA extraction without compromising DNA quality. We tested eight extraction methods – ranging from low‐ to high‐throughput techniques – and eight phyla: Annelida, Arthropoda, Cnidaria, Chordata, Echinodermata, Mollusca, Ochrophyta and Porifera. We assessed DNA yield, purity, efficacy and cost of each method. Extraction efficacy was quantified using the proportion of successful polymerase chain reaction (PCR) amplification of two molecular markers for metazoans (mitochondrial COI and nuclear histone 3) and one for Ochrophyta (mitochondrial nad6) at four time points – 0.5, 1, 2 and 3 years following extraction. DNA yield and purity were quantified using NanoDrop absorbance ratios. Cost was estimated in terms of time and material expense. Results show differences in DNA yield, purity and PCR success between extraction methods and that performance also varied by taxon. The traditional time‐intensive, low‐throughput CTAB phenol–chloroform extraction performed well across taxa, but other methods also performed well and provide the opportunity to reduce time spent at the bench and increase throughput.     

A cost‐effective high‐throughput metabarcoding approach powerful enough to genotype ~44 000 year‐old rodent remains from Northern Africa

We present a cost‐effective metabarcoding approach, aMPlex Torrent, which relies on an improved multiplex PCR adapted to highly degraded DNA, combining barcoding and next‐generation sequencing to simultaneously analyse many heterogeneous samples. We demonstrate the strength of these improvements by generating a phylochronology through the genotyping of ancient rodent remains from a Moroccan cave whose stratigraphy covers the last 120 000 years. Rodents are important for epidemiology, agronomy and ecological investigations and can act as bioindicators for human‐ and/or climate‐induced environmental changes. Efficient and reliable genotyping of ancient rodent remains has the potential to deliver valuable phylogenetic and paleoecological information. The analysis of multiple ancient skeletal remains of very small size with poor DNA preservation, however, requires a sensitive high‐throughput method to generate sufficient data. We show this approach to be particularly adapted at accessing this otherwise difficult taxonomic and genetic resource. As a highly scalable, lower cost and less labour‐intensive alternative to targeted sequence capture approaches, we propose the aMPlex Torrent strategy to be a useful tool for the genetic analysis of multiple degraded samples in studies involving ecology, archaeology, conservation and evolutionary biology.     

Next‐generation monitoring of aquatic biodiversity using environmental DNA metabarcoding

Global biodiversity in freshwater and the oceans is declining at high rates. Reliable tools for assessing and monitoring aquatic biodiversity, especially for rare and secretive species, are important for efficient and timely management. Recent advances in DNA sequencing have provided a new tool for species detection from DNA present in the environment. In this study, we tested whether an environmental DNA (eDNA) metabarcoding approach, using water samples, can be used for addressing significant questions in ecology and conservation. Two key aquatic vertebrate groups were targeted: amphibians and bony fish. The reliability of this method was cautiously validated in silico, in vitro and in situ. When compared with traditional surveys or historical data, eDNA metabarcoding showed a much better detection probability overall. For amphibians, the detection probability with eDNA metabarcoding was 0.97 (CI = 0.90–0.99) vs. 0.58 (CI = 0.50–0.63) for traditional surveys. For fish, in 89% of the studied sites, the number of taxa detected using the eDNA metabarcoding approach was higher or identical to the number detected using traditional methods. We argue that the proposed DNA‐based approach has the potential to become the next‐generation tool for ecological studies and standardized biodiversity monitoring in a wide range of aquatic ecosystems.     

Promises and pitfalls of using high‐throughput sequencing for diet analysis

The application of high‐throughput sequencing‐based approaches to DNA extracted from environmental samples such as gut contents and faeces has become a popular tool for studying dietary habits of animals. Due to the high resolution and prey detection capacity they provide, both metabarcoding and shotgun sequencing are increasingly used to address ecological questions grounded in dietary relationships. Despite their great promise in this context, recent research has unveiled how a wealth of biological (related to the study system) and technical (related to the methodology) factors can distort the signal of taxonomic composition and diversity. Here, we review these studies in the light of high‐throughput sequencing‐based assessment of trophic interactions. We address how the study design can account for distortion factors, and how acknowledging limitations and biases inherent to sequencing‐based diet analyses are essential for obtaining reliable results, thus drawing appropriate conclusions. Furthermore, we suggest strategies to minimize the effect of distortion factors, measures to increase reproducibility, replicability and comparability of studies, and options to scale up DNA sequencing‐based diet analyses. In doing so, we aim to aid end‐users in designing reliable diet studies by informing them about the complexity and limitations of DNA sequencing‐based diet analyses, and encourage researchers to create and improve tools that will eventually drive this field to its maturity.     

Seasonal dynamics and co‐occurrence patterns of honey bee pathogens revealed by high‐throughput RT‐qPCR analysis

The health of the honey bee Apis mellifera is challenged by introduced parasites that interact with its inherent pathogens and cause elevated rates of colony losses. To elucidate co‐occurrence, population dynamics, and synergistic interactions of honey bee pathogens, we established an array of diagnostic assays for a high‐throughput qPCR platform. Assuming that interaction of pathogens requires co‐occurrence within the same individual, single worker bees were analyzed instead of collective samples. Eleven viruses, four parasites, and three pathogenic bacteria were quantified in more than one thousand single bees sampled from sixteen disease‐free apiaries in Southwest Germany. The most abundant viruses were black queen cell virus (84%), Lake Sinai virus 1 (42%), and deformed wing virus B (35%). Forager bees from asymptomatic colonies were infected with two different viruses in average, and simultaneous infection with four to six viruses was common (14%). Also, the intestinal parasites Nosema ceranae (96%) and Crithidia mellificae/Lotmaria passim (52%) occurred very frequently. These results indicate that low‐level infections in honey bees are more common than previously assumed. All viruses showed seasonal variation, while N. ceranae did not. The foulbrood bacteria Paenibacillus larvae and Melissococcus plutonius were regionally distributed. Spearman's correlations and multiple regression analysis indicated possible synergistic interactions between the common pathogens, particularly for black queen cell virus. Beyond its suitability for further studies on honeybees, this targeted approach may be, due to its precision, capacity, and flexibility, a viable alternative to more expensive, sequencing‐based approaches in nonmodel systems.     

A method for high‐throughput production of sequence‐verified DNA libraries and strain collections

The low costs of array‐synthesized oligonucleotide libraries are empowering rapid advances in quantitative and synthetic biology. However, high synthesis error rates, uneven representation, and lack of access to individual oligonucleotides limit the true potential of these libraries. We have developed a cost‐effective method called Recombinase Directed Indexing (REDI), which involves integration of a complex library into yeast, site‐specific recombination to index library DNA, and next‐generation sequencing to identify desired clones. We used REDI to generate a library of ~3,300 DNA probes that exhibited > 96% purity and remarkable uniformity (> 95% of probes within twofold of the median abundance). Additionally, we created a collection of ~9,000 individually accessible CRISPR interference yeast strains for > 99% of genes required for either fermentative or respiratory growth, demonstrating the utility of REDI for rapid and cost‐effective creation of strain collections from oligonucleotide pools. Our approach is adaptable to any complex DNA library, and fundamentally changes how these libraries can be parsed, maintained, propagated, and characterized.     

Diet of the critically endangered blackchin guitarfish Glaucostegus cemiculus revealed using DNA metabarcoding

We present the first assessment of the diet of the blackchin guitarfish Glaucostegus cemiculus (Geoffroy Saint-Hilaire, 1817) for West Africa using DNA metabarcoding on stomach contents of individuals captured in the Bijagós Archipelago, Guinea-Bissau. The diet was dominated by crustaceans, particularly caramote prawn Penaeus kerathurus (frequency of occurrence [FO] = 74%, numerical frequency [NF] = 54%) and fiddler crab Afruca tangeri (FO = 74%, NF = 12%). Bony fishes were present in 30% of the stomachs. We highlight the importance of conservation action for intertidal habitats and their associated benthic invertebrates for the survival of the critically endangered blackchin guitarfish.     

Modular tagging of amplicons using a single PCR for high‐throughput sequencing

High‐throughput sequencing (HTS) of PCR amplicons is becoming the method of choice to sequence one or several targeted loci for phylogenetic and DNA barcoding studies. Although the development of HTS has allowed rapid generation of massive amounts of DNA sequence data, preparing amplicons for HTS remains a rate‐limiting step. For example, HTS platforms require platform‐specific adapter sequences to be present at the 5′ and 3′ end of the DNA fragment to be sequenced. In addition, short multiplex identifier (MID) tags are typically added to allow multiple samples to be pooled in a single HTS run. Existing methods to incorporate HTS adapters and MID tags into PCR amplicons are either inefficient, requiring multiple enzymatic reactions and clean‐up steps, or costly when applied to multiple samples or loci (fusion primers). We describe a method to amplify a target locus and add HTS adapters and MID tags via a linker sequence using a single PCR. We demonstrate our approach by generating reference sequence data for two mitochondrial loci (COI and 16S) for a diverse suite of insect taxa. Our approach provides a flexible, cost‐effective and efficient method to prepare amplicons for HTS.     

Responses of foraminifera communities to aquaculture‐derived organic enrichment as revealed by environmental DNA metabarcoding

Current monitoring methods to assess benthic impacts of marine finfish aquaculture are based on complex biological indices and/or geochemistry data. The former requires benthic macrofauna morpho‐taxonomic characterization that is time‐ and cost‐intensive, while the latter provides rapid assessment of the organic enrichment status of sediments but does not directly measure biotic impacts. In this study, sediment samples were collected from seven stations at six salmon farms in British Columbia, Canada, and analyzed for geochemical parameters and by eDNA metabarcoding to investigate linkages between geochemistry and foraminifera. Sediment texture across farm sites ranged from sand to silty loam, while the maximum sediment pore‐water sulphide concentration at each site ranged from 1,000 to 13,000 μM. Foraminifera alpha diversity generally increased with distance from cage edge. Adonis analyses revealed that farm site explained the most variation in foraminifera community, followed by sediment type, enrichment status, and distance from cage edge. Farm‐specific responses were observed in diversity analyses, taxonomic difference analyses, and correlation analyses. Results demonstrated that species diversity and composition of foraminifera characterized by eDNA metabarcoding generated signals consistent with benthic biodiversity being impacted by finfish farming activities. This substantiates the validity of eDNA metabarcoding for augmenting current approaches to benthic impact assessments by providing more cost‐effective and practicable biotic measures than traditional morpho‐taxonomy. To capitalize on this potential, further work is needed to design a new nomogram that combines eDNA metabarcoding data and geochemistry data to enable accurate monitoring of benthic impacts of fish farming in a time‐ and cost‐efficient way.     

High‐throughput metabarcoding of eukaryotic diversity for environmental monitoring of offshore oil‐drilling activities

As global exploitation of available resources increases, operations extend towards sensitive and previously protected ecosystems. It is important to monitor such areas in order to detect, understand and remediate environmental responses to stressors. The natural heterogeneity and complexity of communities means that accurate monitoring requires high resolution, both temporally and spatially, as well as more complete assessments of taxa. Increased resolution and taxonomic coverage is economically challenging using current microscopy‐based monitoring practices. Alternatively, DNA sequencing‐based methods have been suggested for cost‐efficient monitoring, offering additional insights into ecosystem function and disturbance. Here, we applied DNA metabarcoding of eukaryotic communities in marine sediments, in areas of offshore drilling on the Norwegian continental shelf. Forty‐five samples, collected from seven drilling sites in the Troll/Oseberg region, were assessed, using the small subunit ribosomal RNA gene as a taxonomic marker. In agreement with results based on classical morphology‐based monitoring, we were able to identify changes in sediment communities surrounding oil platforms. In addition to overall changes in community structure, we identified several potential indicator taxa, responding to pollutants associated with drilling fluids. These included the metazoan orders Macrodasyida, Macrostomida and Ceriantharia, as well as several ciliates and other protist taxa, typically not targeted by environmental monitoring programmes. Analysis of a co‐occurrence network to study the distribution of taxa across samples provided a framework for better understanding the impact of anthropogenic activities on the benthic food web, generating novel, testable hypotheses of trophic interactions structuring benthic communities.     

Introgression browser: high‐throughput whole‐genome SNP visualization

Breeding by introgressive hybridization is a pivotal strategy to broaden the genetic basis of crops. Usually, the desired traits are monitored in consecutive crossing generations by marker‐assisted selection, but their analyses fail in chromosome regions where crossover recombinants are rare or not viable. Here, we present the Introgression Browser (iBrowser ), a bioinformatics tool aimed at visualizing introgressions at nucleotide or SNP (Single Nucleotide Polymorphisms) accuracy. The software selects homozygous SNPs from Variant Call Format (VCF) information and filters out heterozygous SNPs, multi‐nucleotide polymorphisms (MNPs) and insertion–deletions (InDels). For data analysis iBrowser makes use of sliding windows, but if needed it can generate any desired fragmentation pattern through General Feature Format (GFF) information. In an example of tomato (Solanum lycopersicum) accessions we visualize SNP patterns and elucidate both position and boundaries of the introgressions. We also show that our tool is capable of identifying alien DNA in a panel of the closely related S. pimpinellifolium by examining phylogenetic relationships of the introgressed segments in tomato. In a third example, we demonstrate the power of the iBrowser in a panel of 597 Arabidopsis accessions, detecting the boundaries of a SNP‐free region around a polymorphic 1.17 Mbp inverted segment on the short arm of chromosome 4. The architecture and functionality of iBrowser makes the software appropriate for a broad set of analyses including SNP mining, genome structure analysis, and pedigree analysis. Its functionality, together with the capability to process large data sets and efficient visualization of sequence variation, makes iBrowser a valuable breeding tool.