Search

The emergence of the Coronavirus Disease 2019 (COVID-19) pandemic has fostered the use of high-throughput techniques to sequence the entire severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome and track its evolution. The present study proposes a rapid and relatively less expensive sequencing protocol for 384 samples by adapting the use of an Illumina NovaSeq library to an Illumina MiSeq flow cell instrument. The SARS-CoV-2 genome sequences obtained with Illumina NovaSeq and those obtained using MiSeq instruments were compared with the objective to validate the new, modified protocol. A total of 356 (94.6%) samples yielded interpretable sequences using the modified Illumina COVIDSeq protocol, with an average coverage of 91.6%. By comparison, 357 (94.9%) samples yielded interpretable sequences with the standard COVIDSeq protocol, with an average coverage of 95.6%. Our modified COVIDSeq protocol could save 14,155 euros per run and yield results from 384 samples in 53.5 h, compared to four times 55.5 h with the standard Illumina MiSeq protocol. The modified COVIDSeq protocol thus provides high quality results comparable to those obtained with the standard COVIDSeq protocol, four times faster, while saving money. [ABSTRACT FROM AUTHOR]

Traditional detection of aquatic invasive species via morphological identification is often time‐consuming and can require a high level of taxonomic expertise, leading to delayed mitigation responses. Environmental DNA (eDNA) detection approaches of multiple species using Illumina‐based sequencing technology have been used to overcome these hindrances, but sample processing is often lengthy. More recently, portable nanopore sequencing technology has become available, which has the potential to make molecular detection of invasive species more widely accessible and substantially decrease sample turnaround times. However, nanopore‐sequenced reads have a much higher error rate than those produced by Illumina platforms, which has so far hindered the adoption of this technology. We provide a detailed laboratory protocol and bioinformatic tools (msi package) to increase the reliability of nanopore sequencing to detect invasive species, and we test its application using invasive bivalves while comparing it with Illumina‐based sequencing. We sampled water from sites with pre‐existing bivalve occurrence and abundance data, and contrasting bivalve communities, in Italy and Portugal. Samples were extracted, amplified, and sequenced by the two platforms. The mean agreement between sequencing methods was 69% and the difference between methods was nonsignificant. The lack of detections of some species at some sites could be explained by their known low abundances. This is the first reported use of MinION to detect aquatic invasive species from eDNA samples. [ABSTRACT FROM AUTHOR]

Premise: Foliar fungal endophytes vary in their distributions across landscapes or plant host taxa, indicative of specialized ecologies and host specific adaptations. Accounts of specialization, however, depend on the taxonomic breadth and geographic range of the host plants included in each study. A broad region‐scale study or deep sampling of diverse potential host species still remains relatively rare but is becoming increasingly possible with high‐throughput sequencing. Methods: Amplicon sequencing was used to rapidly identify the fungal endophytic community among six pine (Pinus, Pinaceae) species co‐occurring across northeastern United States and to test for site and host specialization. We focused on the endophytic genus Lophodermium (Rhytismataceae), whose species members are thought to specialize on different pine species, to test if amplicon sequencing could rapidly verify previously implied or discover new patterns of host specificity. Results: While amplicon sequencing could analyze more samples at greater depths and recover greater numbers of unique Lophodermium taxa than when endophyte communities were surveyed with traditional culturing methods, patterns of specialization were not better supported. This may be because amplicon sequencing can indiscriminately capture non‐host specific organisms found incidentally from plant tissues or because we have overestimated host‐specificity in the past with biased culturing techniques. Conclusions: Amplicon sequencing can quickly identify patterns of host specificity by allowing large‐scale surveys but has limitations in quantifying the level of intimacy of these relationships. [ABSTRACT FROM AUTHOR]

Abstract Background Current library preparation protocols for Illumina HiSeq and MiSeq DNA sequencers require ≥2 nM initial library for subsequent loading of denatured cDNA onto flow cells. Such amounts are not always attainable from samples having a relatively low DNA or RNA input; or those for which a limited number of PCR amplification cycles is preferred (less PCR bias and/or more even coverage). A well-tested sub-nanomolar library preparation protocol for Illumina sequencers has however not been reported. The aim of this study is to provide a much needed working protocol for sub-nanomolar libraries to achieve outcomes as informative as those obtained with the higher library input (≥ 2 nM) recommended by Illumina’s protocols. Results Extensive studies were conducted to validate a robust sub-nanomolar (initial library of 100 pM) protocol using PhiX DNA (as a control), genomic DNA (Bordetella bronchiseptica and microbial mock community B for 16S rRNA gene sequencing), messenger RNA, microRNA, and other small noncoding RNA samples. The utility of our protocol was further explored for PhiX library concentrations as low as 25 pM, which generated only slightly fewer than 50% of the reads achieved under the standard Illumina protocol starting with > 2 nM. Conclusions A sub-nanomolar library preparation protocol (100 pM) could generate next generation sequencing (NGS) results as robust as the standard Illumina protocol. Following the sub-nanomolar protocol, libraries with initial concentrations as low as 25 pM could also be sequenced to yield satisfactory and reproducible sequencing results.

Abstract Background Knowledge on the capacity of Australian ticks to carry Borrelia species is currently limited or missing. To evaluate the potential of ticks to carry bacterial pathogens and their DNA, it is imperative to have a robust workflow that maximises recovery of bacterial DNA within ticks in order to enable accurate identification. By exploiting the bilateral anatomical symmetry of ticks, we were able to directly compare two DNA extraction methods for 16S rRNA gene diversity profiling and pathogen detection. We aimed to assess which combination of DNA extraction and 16S rRNA hypervariable region enables identification of the greatest bacterial diversity, whilst minimising bias, and providing the greatest capacity for the identification of Borrelia spp. Results We collected Australian endemic ticks (Bothriocroton undatum), isolated DNA from equal tick halves using two commercial DNA extraction methods and sequenced samples using V1-V3 and V3-V4 16S rRNA gene diversity profiling assays. Two distinct Borrelia spp. operational taxonomic units (OTUs) were detected using the V1-V3 16S rRNA hypervariable region and matching Borrelia spp. sequences were obtained using a conventional nested-PCR. The tick 16S rRNA gene diversity profile was dominated by Rickettsia spp. (98–99%), while the remaining OTUs belonged to Proteobacteria (51–81%), Actinobacteria (6–30%) and Firmicutes (2–7%). Multiple comparisons tests demonstrated biases in each of the DNA extraction kits towards different bacterial taxa. Conclusions Two distinct Borrelia species belonging to the reptile-associated Borrelia group were identified. Our results show that the method of DNA extraction can promote bias in the final microbiota identified. We determined an optimal DNA extraction method and 16S rRNA gene diversity profile assay that maximises detection of Borrelia species.

High-throughput amplicon sequencing that primarily targets the 16S ribosomal DNA (rDNA) (for bacteria and archaea) and the Internal Transcribed Spacer rDNA (for fungi) have facilitated microbial community discovery across diverse environments. A three-step PCR that utilizes flexible primer choices to construct the library for Illumina amplicon sequencing has been applied to several studies in forest and agricultural systems. The three-step PCR protocol, while producing high-quality reads, often yields a large number (up to 46%) of reads that are unable to be assigned to a specific sample according to its barcode. Here, we improve this technique through an optimized two-step PCR protocol. We tested and compared the improved two-step PCR meta-barcoding protocol against the three-step PCR protocol using four different primer pairs (fungal ITS: ITS1F-ITS2 and ITS1F-ITS4, and bacterial 16S: 515F-806R and 341F-806R). We demonstrate that the sequence quantity and recovery rate were significantly improved with the two-step PCR approach (fourfold more read counts per sample; determined reads ≈90% per run) while retaining high read quality (Q30 > 80%). Given that synthetic barcodes are incorporated independently from any specific primers, this two-step PCR protocol can be broadly adapted to different genomic regions and organisms of scientific interest.

Abstract: Despite the ever increase in rigorous control and monitoring measures to assure safe food along the entire farm‐to‐fork chain, the past decade has also witnessed an increase in microbial food alerts. Hence, research on food safety and quality remain of utmost importance. Complementary, and at least as important, is the necessity to be able to assess the potential microbial risks along the food chain. Risk assessment relies on sound scientific data. Unfortunately, often, quality data are limited if not lacking. High‐throughput tools such as next‐generation sequencing (NGS) could fill this gap. NGS approaches can be used to generate ample qualitative and quantitative data to be used in the risk assessment process. NGS applications are not new in food microbiology with applications ranging from pathogen detection along the food chain, food epidemiology studies, whole genome analysis of food‐associated microorganisms up to describing complete food microbiomes. Yet, its application in the area of microbial risk assessment is still at an early stage and faces important challenges. The possibilities of NGS for risk assessment are ample, but so are the questions on the subject. One of the major strengths of NGS lies in its capacity to generate a lot of data, but to what extend can this wealth be of use in hazard identification, hazard characterisation and exposure assessment to perform a sound risk characterisation, which in turn will make it possible to take substantiated risk management decisions. [ABSTRACT FROM AUTHOR]

Background: Current library preparation protocols for Illumina HiSeq and MiSeq DNA sequencers require ≥2 nM initial library for subsequent loading of denatured cDNA onto flow cells. Such amounts are not always attainable from samples having a relatively low DNA or RNA input; or those for which a limited number of PCR amplification cycles is preferred (less PCR bias and/or more even coverage). A well-tested sub-nanomolar library preparation protocol for Illumina sequencers has however not been reported. The aim of this study is to provide a much needed working protocol for sub-nanomolar libraries to achieve outcomes as informative as those obtained with the higher library input (≥ 2 nM) recommended by Illumina's protocols. Results: Extensive studies were conducted to validate a robust sub-nanomolar (initial library of 100 pM) protocol using PhiX DNA (as a control), genomic DNA (Bordetella bronchiseptica and microbial mock community B for 16S rRNA gene sequencing), messenger RNA, microRNA, and other small noncoding RNA samples. The utility of our protocol was further explored for PhiX library concentrations as low as 25 pM, which generated only slightly fewer than 50% of the reads achieved under the standard Illumina protocol starting with > 2 nM. Conclusions: A sub-nanomolar library preparation protocol (100 pM) could generate next generation sequencing (NGS) results as robust as the standard Illumina protocol. Following the sub-nanomolar protocol, libraries with initial concentrations as low as 25 pM could also be sequenced to yield satisfactory and reproducible sequencing results. [ABSTRACT FROM AUTHOR]

Background: Knowledge on the capacity of Australian ticks to carry Borrelia species is currently limited or missing. To evaluate the potential of ticks to carry bacterial pathogens and their DNA, it is imperative to have a robust workflow that maximises recovery of bacterial DNA within ticks in order to enable accurate identification. By exploiting the bilateral anatomical symmetry of ticks, we were able to directly compare two DNA extraction methods for 16S rRNA gene diversity profiling and pathogen detection. We aimed to assess which combination of DNA extraction and 16S rRNA hypervariable region enables identification of the greatest bacterial diversity, whilst minimising bias, and providing the greatest capacity for the identification of Borrelia spp. Results: We collected Australian endemic ticks (Bothriocroton undatum), isolated DNA from equal tick halves using two commercial DNA extraction methods and sequenced samples using V1-V3 and V3-V4 16S rRNA gene diversity profiling assays. Two distinct Borrelia spp. operational taxonomic units (OTUs) were detected using the V1-V3 16S rRNA hypervariable region and matching Borrelia spp. sequences were obtained using a conventional nested-PCR. The tick 16S rRNA gene diversity profile was dominated by Rickettsia spp. (98-99%), while the remaining OTUs belonged to Proteobacteria (51-81%), Actinobacteria (6-30%) and Firmicutes (2-7%). Multiple comparisons tests demonstrated biases in each of the DNA extraction kits towards different bacterial taxa. Conclusions: Two distinct Borrelia species belonging to the reptile-associated Borrelia group were identified. Our results show that the method of DNA extraction can promote bias in the final microbiota identified. We determined an optimal DNA extraction method and 16S rRNA gene diversity profile assay that maximises detection of Borrelia species. [ABSTRACT FROM AUTHOR]

Illumina library preparation methods for ultra-low input amounts were compared using genomic DNA from two foodborne parasites ( Angiostrongylus cantonensis and Cyclospora cayetanensis ) as examples. The Ovation Ultralow method resulted in libraries with the highest concentration and produced quality sequencing data, even when the input DNA was in the picogram range. [ABSTRACT FROM AUTHOR]

Classical Sanger sequencing is still frequently used to generate sequence data for phylogenetic and phylogeographic inference. In this contribution we present a novel approach to genotype whole mitogenomic haplotypes using Illumina MiSeq reads from indexed amplicons. Our new approach reduces preparation time by multiplexing loci within a single or few PCR reactions and by plate format library construction. The use of paired-end reads allows covering amplicons of about 0.5 kb and thus no nebulisation and assembly are necessary. We tested the power and effectiveness of this technique by analysing the mitogenomic diversity of East African white-eye bird species (Zosteropidae), a taxonomically highly diverse and complex species flock found in various ecosystems spread across major parts of Africa. We compare the newly generated mitogenomic data set with published data of three mitochondrial genes for a similar set of populations and taxa. The comparison demonstrates that our new procedure represents a cost effective use of NGS for medium throughput phylogenetic analyses. Using this method, we were able to increase the amount of phylogenetic information significantly, while reducing the costs and effort in the laboratory. The mitogenomic data show a higher resolution than previous studies providing higher support and new insights in the relationships of Zosterops species. Our data suggest to split Z. poliogaster into four distinct species, three of which had previously been proposed: Z. silvanus , Z. mbulensis , Z. kikyuensis and Z. kulalensis . Our approach allows the genotyping of whole mitogenomes for a large number of individuals and thus allows more reliable reconstruction of phylogenetic and phylogeographic relationships – also for non-model organisms. [ABSTRACT FROM AUTHOR]

Current forensic DNA analysis predominantly involves identification of human donors by analysis of short tandem repeats (STRs) using Capillary Electrophoresis (CE). Recent developments in Massively Parallel Sequencing (MPS) technologies offer new possibilities in analysis of STRs since they might overcome some of the limitations of CE analysis. In this study 17 STRs and Amelogenin were sequenced in high coverage using a prototype version of the Promega PowerSeq™ system for 297 population samples from the Netherlands, Nepal, Bhutan and Central African Pygmies. In addition, 45 two-person mixtures with different minor contributions down to 1% were analysed to investigate the performance of this system for mixed samples. Regarding fragment length, complete concordance between the MPS and CE-based data was found, marking the reliability of MPS PowerSeq™ system. As expected, MPS presented a broader allele range and higher power of discrimination and exclusion rate. The high coverage sequencing data were used to determine stutter characteristics for all loci and stutter ratios were compared to CE data. The separation of alleles with the same length but exhibiting different stutter ratios lowers the overall variation in stutter ratio and helps in differentiation of stutters from genuine alleles in mixed samples. All alleles of the minor contributors were detected in the sequence reads even for the 1% contributions, but analysis of mixtures below 5% without prior information of the mixture ratio is complicated by PCR and sequencing artefacts. [ABSTRACT FROM AUTHOR]

High-throughput sequencing platforms are continuing to increase resulting read lengths, which is allowing for a deeper and more accurate depiction of environmental microbial diversity. With the nascent Reagent Kit v3, Illumina MiSeq now has the ability to sequence the eukaryotic hyper-variable V4 region of the SSU-r DNA locus with paired-end reads. Using DNA collected from soils with analyses of strictly- and nearly identical amplicons, here we ask how the new Illumina MiSeq data compares with what we can obtain with Roche/454 GS FLX with regard to quantity and quality, presence and absence, and abundance perspectives. We show that there is an easy qualitative transition from the Roche/454 to the Illumina MiSeq platforms. The ease of this transition is more nuanced quantitatively for low-abundant amplicons, although estimates of abundances are known to also vary within platforms. [ABSTRACT FROM AUTHOR]

Routine use of massively parallel sequencing (MPS) for forensic genomics is on the horizon. The last few years, several algorithms and workflows have been developed to analyze forensic MPS data. However, none have yet been tailored to the needs of the forensic analyst who does not possess an extensive bioinformatics background. We developed our previously published forensic MPS data analysis framework MyFLq (My-Forensic-Loci-queries) into an open-source, user-friendly, web-based application. It can be installed as a standalone web application, or run directly from the Illumina BaseSpace environment. In the former, laboratories can keep their data on-site, while in the latter, data from forensic samples that are sequenced on an Illumina sequencer can be uploaded to Basespace during acquisition, and can subsequently be analyzed using the published MyFLq BaseSpace application. Additional features were implemented such as an interactive graphical report of the results, an interactive threshold selection bar, and an allele length-based analysis in addition to the sequenced-based analysis. Practical use of the application is demonstrated through the analysis of four 16-plex short tandem repeat (STR) samples, showing the complementarity between the sequence- and length-based analysis of the same MPS data. [ABSTRACT FROM AUTHOR]

Abstract: Forensic scientists are currently investigating how to transition from capillary electrophoresis (CE) to massive parallel sequencing (MPS) for analysis of forensic DNA profiles. MPS offers several advantages over CE such as virtually unlimited multiplexy of loci, combining both short tandem repeat (STR) and single nucleotide polymorphism (SNP) loci, small amplicons without constraints of size separation, more discrimination power, deep mixture resolution and sample multiplexing. We present our bioinformatic framework My-Forensic-Loci-queries (MyFLq) for analysis of MPS forensic data. For allele calling, the framework uses a MySQL reference allele database with automatically determined regions of interest (ROIs) by a generic maximal flanking algorithm which makes it possible to use any STR or SNP forensic locus. Python scripts were designed to automatically make allele calls starting from raw MPS data. We also present a method to assess the usefulness and overall performance of a forensic locus with respect to MPS, as well as methods to estimate whether an unknown allele, which sequence is not present in the MySQL database, is in fact a new allele or a sequencing error. The MyFLq framework was applied to an Illumina MiSeq dataset of a forensic Illumina amplicon library, generated from multilocus STR polymerase chain reaction (PCR) on both single contributor samples and multiple person DNA mixtures. Although the multilocus PCR was not yet optimized for MPS in terms of amplicon length or locus selection, the results show excellent results for most loci. The results show a high signal-to-noise ratio, correct allele calls, and a low limit of detection for minor DNA contributors in mixed DNA samples. Technically, forensic MPS affords great promise for routine implementation in forensic genomics. The method is also applicable to adjacent disciplines such as molecular autopsy in legal medicine and in mitochondrial DNA research. [Copyright &y& Elsevier]

Background: Next generation sequencing (NGS) technology has revolutionized genomic and genetic research. The pace of change in this area is rapid with three major new sequencing platforms having been released in 2011: Ion Torrent's PGM, Pacific Biosciences' RS and the Illumina MiSeq. Here we compare the results obtained with those platforms to the performance of the Illumina HiSeq, the current market leader. In order to compare these platforms, and get sufficient coverage depth to allow meaningful analysis, we have sequenced a set of 4 microbial genomes with mean GC content ranging from 19.3 to 67.7%. Together, these represent a comprehensive range of genome content. Here we report our analysis of that sequence data in terms of coverage distribution, bias, GC distribution, variant detection and accuracy. Results: Sequence generated by Ion Torrent, MiSeq and Pacific Biosciences technologies displays near perfect coverage behaviour on GC-rich, neutral and moderately AT-rich genomes, but a profound bias was observed upon sequencing the extremely AT-rich genome of Plasmodium falciparum on the PGM, resulting in no coverage for approximately 30% of the genome. We analysed the ability to call variants from each platform and found that we could call slightly more variants from Ion Torrent data compared to MiSeq data, but at the expense of a higher false positive rate. Variant calling from Pacific Biosciences data was possible but higher coverage depth was required. Context specific errors were observed in both PGM and MiSeq data, but not in that from the Pacific Biosciences platform. Conclusions: All three fast turnaround sequencers evaluated here were able to generate usable sequence. However there are key differences between the quality of that data and the applications it will support. [ABSTRACT FROM AUTHOR]

Kudzu ( Pueraria montana var. lobata) is a perennial vine native to eastern Asia and an invasive in the United States, representing a significant threat to biodiversity. We characterized 15 polymorphic microsatellite loci in four populations from both the native and invasive ranges of kudzu. The number of alleles per locus ranged from 2 to 8 and observed heterozygosity ranged from 0.0 to 0.9. Consistent with vegetative reproduction, loci were initially not in Hardy-Weinberg equilibrium (HWE) and in linkage disequilibrium but when clones were removed, markers were in HWE and in linkage equilibrium. These loci provide valuable tools to study the diversity and invasion history of an important invasive species. [ABSTRACT FROM AUTHOR]

High-throughput amplicon sequencing that primarily targets the 16S ribosomal DNA (rDNA) (for bacteria and archaea) and the Internal Transcribed Spacer rDNA (for fungi) have facilitated microbial community discovery across diverse environments. A three-step PCR that utilizes flexible primer choices to construct the library for Illumina amplicon sequencing has been applied to several studies in forest and agricultural systems. The three-step PCR protocol, while producing high-quality reads, often yields a large number (up to 46%) of reads that are unable to be assigned to a specific sample according to its barcode. Here, we improve this technique through an optimized two-step PCR protocol. We tested and compared the improved two-step PCR meta-barcoding protocol against the three-step PCR protocol using four different primer pairs (fungal ITS: ITS1F-ITS2 and ITS1F-ITS4, and bacterial 16S: 515F-806R and 341F-806R). We demonstrate that the sequence quantity and recovery rate were significantly improved with the two-step PCR approach (fourfold more read counts per sample; determined reads ≈90% per run) while retaining high read quality (Q30 > 80%). Given that synthetic barcodes are incorporated independently from any specific primers, this two-step PCR protocol can be broadly adapted to different genomic regions and organisms of scientific interest. [ABSTRACT FROM AUTHOR]