Comparison of the performance of an amplicon sequencing assay based on Oxford Nanopore technology to real-time PCR assays for detecting bacterial biodefense pathogens.

BACKGROUND: The state-of-the-art in nucleic acid based biodetection continues to be polymerase chain reaction (PCR), and many real-time PCR assays targeting biodefense pathogens for biosurveillance are in widespread use. These assays are predominantly singleplex; i.e. one assay tests for the presence of one target, found in a single organism, one sample at a time. Due to the intrinsic limitations of such tests, there exists a critical need for high-throughput multiplex assays to reduce the time and cost incurred when screening multiple targets, in multiple pathogens, and in multiple samples. Such assays allow users to make an actionable call while maximizing the utility of the small volumes of test samples. Unfortunately, current multiplex real-time PCR assays are limited in the number of targets that can be probed simultaneously due to the availability of fluorescence channels in real-time PCR instruments. RESULTS: To address this gap, we developed a pipeline in which the amplicons produced by a 14-plex end-point PCR assay using spiked samples were subsequently sequenced using Nanopore technology. We used bar codes to sequence multiple samples simultaneously, leading to the generation and subsequent analysis of sequence data resulting from a short sequencing run time (< 10 min). We compared the limits of detection (LoD) of real-time PCR assays to Oxford Nanopore Technologies (ONT)-based amplicon sequencing and estimated the sample-to-answer time needed for this approach. Overall, LoDs determined from the first 10 min of sequencing data were at least one to two orders of magnitude lower than real-time PCR. Given enough time, the amplicon sequencing approach is approximately 100 times more sensitive than real-time PCR, with detection of amplicon specific reads even at the lowest tested spiking concentration (around 2.5-50 Colony Forming Units (CFU)/ml). CONCLUSIONS: Based on these results, we propose amplicon sequencing assay as a viable alternative to replace the current real-time PCR based singleplex assays for higher throughput biodefense applications. We note, however, that targeted amplicon specific reads were not detectable even at the highest tested spike concentrations (2.5 X 104-5.0 X105 CFU/ml) without an initial amplification step, indicating that PCR is still necessary when utilizing this protocol.

A novel plasmid, pSAA0430-08, from Streptococcus anginosus subsp. anginosus strain 0430-08.

Mobile genetic elements (MGEs) are the genetic material often involved in the interspecies and intraspecies genetic transduction in bacteria. However, little is known about MGEs in the Anginosus group of streptococci (AGS), one of the streptococcal groups found in the oral cavity of humans. We looked for the presence of MGEs in Streptococcus anginosus subsp. anginosus (SAA), a representative species belonging to AGS, and found a novel plasmid from SAA strain 0430-08. This plasmid was 7038bp and ~31% G/C content which we named pSAA0430-08, and examined its genetic structure and characteristics. Open reading frame (ORF) prediction revealed that pSAA0430-08 was composed of 10 ORFs including a putative plasmid replication protein (ORF1) and a putative toxin-antitoxin system (ORF9 and ORF10). Between ORF10 and ORF 1, four tandem repeats of 22bp each, generally termed as iteron, were also observed. Using variant plasmids of pSAA0430-08, we confirmed that both ORF1 and iteron were necessary for replication in host cells. Interestingly, the region from ORF4 to ORF7 showed homology with a genomic DNA segment of S. gordonii strains. Thus, this plasmid may travel between the different species in Streptococci, i.e., S. gordonii and S. anginosus.

Draft Genome Assemblies of Phage AP50c-Resistant Derivatives of Bacillus anthracis Sterne Strain 7702 Lacking Plasmid pXO2.

Mutants of the attenuated Bacillus anthracis (Sterne) strain 7702 that are resistant to phage AP50c have been previously described. Here, we report the draft genome assemblies of the parent strain, several phage-resistant derivatives, and mutants of genes in the pathways for synthesis and assembly of the S-layer.

Genetic evidence for the interaction between Bacillus anthracis-encoded phage receptors and their cognate phage-encoded receptor binding proteins.

Bacteriophages such as γ and AP50c have been shown to infect strains of Bacillus anthracis with high specificity, and this feature has been exploited in the development of bacterial detection assays. To better understand the emergence of phage resistance, and thus the potential failure of such assays, it is important to identify the host and phage receptors necessary for attachment and entry. Using genetic approaches, the bacterial receptors of AP50c and γ have been identified as sap and GamR, respectively. A second AP50c-like phage, Wip1, also appears to use sap as a receptor. In parallel with this work, the cognate phage-encoded receptor binding proteins (RBPs) have also been identified (Gp14 for γ, P28 for AP50c, and P23 for Wip1); however, the strength of evidence supporting these protein-protein interactions varies, necessitating additional investigation. Here, we present genetic evidence further supporting the interaction between sap and the RBPs of AP50c and Wip1 using fluorescently tagged proteins and a panel of B. anthracis mutants. These results showed that the deletion of the sap gene, as well as the deletion of csaB, whose encoded protein anchors sap to the bacterial S-layer, resulted in the loss of RBP binding. Binding could then be rescued by expressing these genes in trans. We also found that the RBP of the γ-like prophage λBa03 relied on csaB activity for binding, possibly by a different mechanism. RBPλBa03 binding to B. anthracis cells was also unique in that it was not ablated by heat inactivation of vegetative cells, suggesting that its receptor is still functional following incubation at 98°C. These results extend our understanding of the diverse attachment and entry strategies used by B. anthracis phages, enabling future assay development.

Optimization of Oxford Nanopore Technology Sequencing Workflow for Detection of Amplicons in Real Time Using ONT-DART Tool.

An optimized, well-tested and validated targeted genomic sequencing-based high-throughput assay is currently not available ready for routine biodefense and biosurveillance applications. Earlier, we addressed this gap by developing and establishing baseline comparisons of a multiplex end-point Polymerase Chain Reaction (PCR) assay followed by Oxford Nanopore Technology (ONT) based amplicon sequencing to real time PCR and customized data processing. Here, we expand upon this effort by identifying the optimal ONT library preparation method for integration into a novel software platform ONT-DART (ONT-Detection of Amplicons in Real-Time). ONT-DART is a dockerized, real-time, amplicon-sequence analysis workflow that is used to reproducibly process and filter read data to support actionable amplicon detection calls based on alignment metrics, within sample statistics, and no-template control data. This analysis pipeline was used to compare four ONT library preparation protocols using R9 and Flongle (FL) flow cells. The two 4-Primer methods tested required the shortest preparation times (5.5 and 6.5 h) for 48 libraries but provided lower fidelity data. The Native Barcoding and Ligation methods required longer preparation times of 8 and 12 h, respectively, and resulted in higher overall data quality. On average, data derived from R9 flow cells produced true positive calls for target organisms more than twice as fast as the lower throughput FL flow cells. These results suggest that utilizing the R9 flowcell with an ONT Native Barcoding amplicon library method in combination with ONT-DART platform analytics provides the best sequencing-based alternative to current PCR-based biodetection methods.

The META tool optimizes metagenomic analyses across sequencing platforms and classifiers.

A major challenge in the field of metagenomics is the selection of the correct combination of sequencing platform and downstream metagenomic analysis algorithm, or "classifier". Here, we present the Metagenomic Evaluation Tool Analyzer (META), which produces simulated data and facilitates platform and algorithm selection for any given metagenomic use case. META-generated in silico read data are modular, scalable, and reflect user-defined community profiles, while the downstream analysis is done using a variety of metagenomic classifiers. Reported results include information on resource utilization, time-to-answer, and performance. Real-world data can also be analyzed using selected classifiers and results benchmarked against simulations. To test the utility of the META software, simulated data was compared to real-world viral and bacterial metagenomic samples run on four different sequencers and analyzed using 12 metagenomic classifiers. Lastly, we introduce "META Score": a unified, quantitative value which rates an analytic classifier's ability to both identify and count taxa in a representative sample.

A novel canis lupus familiaris reference genome improves variant resolution for use in breed-specific GWAS.

Reference genome fidelity is critically important for genome wide association studies, yet most vary widely from the study population. A typical whole genome sequencing approach implies short-read technologies resulting in fragmented assemblies with regions of ambiguity. Further information is lost by economic necessity when genotyping populations, as lower resolution technologies such as genotyping arrays are commonly used. Here, we present a phased reference genome for Canis lupus familiaris using high molecular weight DNA-sequencing technologies. We tested wet laboratory and bioinformatic approaches to demonstrate a minimum workflow to generate the 2.4 gigabase genome for a Labrador Retriever. The de novo assembly required eight Oxford Nanopore R9.4 flowcells (∼23X depth) and running a 10X Genomics library on the equivalent of one lane of an Illumina NovaSeq S1 flowcell (∼88X depth), bringing the cost of generating a nearly complete reference genome to less than $10K (USD). Mapping of short-read data from 10 Labrador Retrievers against this reference resulted in 1% more aligned reads versus the current reference (CanFam3.1, P < 0.001), and a 15% reduction of variant calls, increasing the chance of identifying true, low-effect size variants in a genome-wide association studies. We believe that by incorporating the cost to produce a full genome assembly into any large-scale genotyping project, an investigator can improve study power, decrease costs, and optimize the overall scientific value of their study.

Complete Genome Sequence of Francisella tularensis Live Vaccine Strain NR-28537 (BEI Master Cell Bank).

The genome of Francisella tularensis live vaccine strain NR-28537 was sequenced by a hybrid approach utilizing an Oxford Nanopore Technologies R9 flow cell and an Illumina MiSeq platform. De novo assembly of the resulting long and short reads produced a single-contig whole-genome sequence.

A bioinformatic filter for improved base-call accuracy and polymorphism detection using the Affymetrix GeneChip whole-genome resequencing platform.

DNA resequencing arrays enable rapid acquisition of high-quality sequence data. This technology represents a promising platform for rapid high-resolution genotyping of microorganisms. Traditional array-based resequencing methods have relied on the use of specific PCR-amplified fragments from the query samples as hybridization targets. While this specificity in the target DNA population reduces the potential for artifacts caused by cross-hybridization, the subsampling of the query genome limits the sequence coverage that can be obtained and therefore reduces the technique's resolution as a genotyping method. We have developed and validated an Affymetrix Inc. GeneChip(R) array-based, whole-genome resequencing platform for Francisella tularensis, the causative agent of tularemia. A set of bioinformatic filters that targeted systematic base-calling errors caused by cross-hybridization between the whole-genome sample and the array probes and by deletions in the sample DNA relative to the chip reference sequence were developed. Our approach eliminated 91% of the false-positive single-nucleotide polymorphism calls identified in the SCHU S4 query sample, at the cost of 10.7% of the true positives, yielding a total base-calling accuracy of 99.992%.

Extra-Chromosomal DNA Sequencing Reveals Episomal Prophages Capable of Impacting Virulence Factor Expression in Staphylococcus aureus.

Staphylococcus aureus is a major human pathogen with well-characterized bacteriophage contributions to its virulence potential. Recently, we identified plasmidial and episomal prophages in S. aureus strains using an extra-chromosomal DNA (exDNA) isolation and sequencing approach, uncovering the plasmidial phage ϕBU01, which was found to encode important virulence determinants. Here, we expanded our extra-chromosomal sequencing of S. aureus, selecting 15 diverse clinical isolates with known chromosomal sequences for exDNA isolation and next-generation sequencing. We uncovered the presence of additional episomal prophages in 5 of 15 samples, but did not identify any plasmidial prophages. exDNA isolation was found to enrich for circular prophage elements, and qPCR characterization of the strains revealed that such prophage enrichment is detectable only in exDNA samples and would likely be missed in whole-genome DNA preparations (e.g., detection of episomal prophages did not correlate with higher prophage excision rates nor higher excised prophage copy numbers in qPCR experiments using whole-genome DNA). In S. aureus MSSA476, we found that enrichment and excision of the prophage ϕSa4ms into the cytoplasm was temporal and that episomal prophage localization did not appear to be a precursor to lytic cycle replication, suggesting ϕSa4ms excision into the cytoplasm may be part of a novel lysogenic switch. For example, we show that ϕSa4ms excision alters the promoter and transcription of htrA2 , encoding a stress-response serine protease, and that alternative promotion of htrA2 confers increased heat-stress survival in S. aureus COL. Overall, exDNA isolation and focused sequencing may offer a more complete genomic picture for bacterial pathogens, offering insights into important chromosomal dynamics likely missed with whole-genome DNA-based approaches.

Iterative microarray and RNA interference-based interrogation of the SRC-induced invasive phenotype.

Src kinase has long been recognized as a factor in the progression of colorectal cancer and seems to play a specific role in the development of the metastatic phenotype. In spite of numerous studies conducted to elucidate the exact role of Src in cancer progression, downstream targets of Src remain poorly understood. Gene expression profiling has permitted the identification of large sets of genes that may be functionally interrelated but it is often unclear as to which molecular pathways they belong. Here we have developed an iterative approach to experimentally reconstruct a network of gene activity regulated by Src and contributing to the invasive phenotype. Our strategy uses a combination of phenotypic anchoring of gene expression profiles and loss-of-function screening by way of RNA-mediated interference. Using a panel of human colon cancer cell lines exhibiting differential Src-specific activity and invasivity, we identify the first two levels of gene transcription responsible for the invasive phenotype, where first-tier genes are controlled by Src activity and the second-tier genes are under the influence of the first tier. Specifically, perturbation of first-tier gene activity by either pharmacologic inhibition of Src or RNA-mediated interference-directed knockdown leads to a loss of invasivity and decline of second-tier gene activity. The targeting of first-tier genes may be bypassed altogether because knockdown of second-tier genes led to a similar loss of invasive potential. In this manner, numerous members of a "transcriptional cascade" pathway for metastatic activity have been identified and functionally validated.

Transcriptional profiling with a blood pressure QTL interval-specific oligonucleotide array.

Although the evidence for a genetic predisposition to human essential hypertension is compelling, the genetic control of blood pressure (BP) is poorly understood. The Dahl salt-sensitive (S) rat is a model for studying the genetic component of BP. Using this model, we previously reported the identification of 16 different genomic regions that contain one or more BP quantitative trait loci (QTLs). The proximal region of rat chromosome 1 contains multiple BP QTLs. Of these, we have localized the BP QTL1b region to a 13.5-cM (20.92 Mb) region. Interestingly, five additional independent studies in rats and four independent studies in humans have reported genetic linkage for BP control by regions homologous to QTL1b. To view the overall renal transcriptional topography of the positional candidate genes for this QTL, we sought a comparative gene expression profiling between a congenic strain containing QTL1b and control S rats by employing 1) a saturated QTL1b interval-specific oligonucleotide array and 2) a whole genome cDNA microarray representing 20,465 unique genes that are positioned outside the QTL. Results indicated that 17 of the 231 positional candidate genes for this QTL are differentially expressed between the two strains tested. Surprisingly, >1,500 genes outside of QTL1b were differentially expressed between the two rat strains. Integrating the results from the two approaches revealed at least one complex network of transcriptional control initiated by the positional candidate Nr2f2. This network appears to account for the majority of gene expression differences occurring outside of the QTL interval. Further substitution mapping is currently underway to test the validity of each of these differentially expressed positional candidate genes. These results demonstrate the importance of using a saturated oligonucleotide array for identifying and prioritizing differentially expressed positional candidate genes of a BP QTL.

Genome sequencing of 15 clinical Vibrio isolates, including 13 non-o1/non-o139 serogroup strains.

We present draft genome sequences of 15 clinical Vibrio isolates of various serogroups. These are valuable data for use in studying Vibrio cholerae genetic diversity, epidemic potential, and strain attribution.

Beyond the chromosome: the prevalence of unique extra-chromosomal bacteriophages with integrated virulence genes in pathogenic Staphylococcus aureus.

In Staphylococcus aureus, the disease impact of chromosomally integrated prophages on virulence is well described. However, the existence of extra-chromosomal prophages, both plasmidial and episomal, remains obscure. Despite the recent explosion in bacterial and bacteriophage genomic sequencing, studies have failed to specifically focus on extra-chromosomal elements. We selectively enriched and sequenced extra-chromosomal DNA from S. aureus isolates using Roche-454 technology and uncovered evidence for the widespread distribution of multiple extra-chromosomal prophages (ExPΦs) throughout both antibiotic-sensitive and -resistant strains. We completely sequenced one such element comprised of a 43.8 kbp, circular ExPΦ (designated ФBU01) from a vancomycin-intermediate S. aureus (VISA) strain. Assembly and annotation of ФBU01 revealed a number of putative virulence determinants encoded within a bacteriophage immune evasion cluster (IEC). Our identification of several potential ExPΦs and mobile genetic elements (MGEs) also revealed numerous putative virulence factors and antibiotic resistance genes. We describe here a previously unidentified level of genetic diversity of stealth extra-chromosomal elements in S. aureus, including phages with a larger presence outside the chromosome that likely play a prominent role in pathogenesis and strain diversity driven by horizontal gene transfer (HGT).

Genome Sequencing of Four Strains of Rickettsia prowazekii, the Causative Agent of Epidemic Typhus, Including One Flying Squirrel Isolate.

Rickettsia prowazekii is a notable intracellular pathogen, the agent of epidemic typhus, and a potential biothreat agent. We present here whole-genome sequence data for four strains of R. prowazekii, including one from a flying squirrel.

A Yersinia pestis-specific, lytic phage preparation significantly reduces viable Y. pestis on various hard surfaces experimentally contaminated with the bacterium.

Five Y. pestis bacteriophages obtained from various sources were characterized to determine their biological properties, including their taxonomic classification, host range and genomic diversity. Four of the phages (YpP-G, Y, R and YpsP-G) belong to the Podoviridae family, and the fifth phage (YpsP-PST) belongs to the Myoviridae family, of the order Caudovirales comprising of double-stranded DNA phages. The genomes of the four Podoviridae phages were fully sequenced and found to be almost identical to each other and to those of two previously characterized Y. pestis phages Yepe2 and φA1122. However, despite their genomic homogeneity, they varied in their ability to lyse Y. pestis and Y. pseudotuberculosis strains. The five phages were combined to yield a "phage cocktail" (tentatively designated "YPP-100") capable of lysing the 59 Y. pestis strains in our collection. YPP-100 was examined for its ability to decontaminate three different hard surfaces (glass, gypsum board and stainless steel) experimentally contaminated with a mixture of three genetically diverse Y. pestis strains CO92, KIM and 1670G. Five minutes of exposure to YPP-100 preparations containing phage concentrations of ca. 10(9), 10(8) and 10(7) PFU/mL completely eliminated all viable Y. pestis cells from all three surfaces, but a few viable cells were recovered from the stainless steel coupons treated with YPP-100 diluted to contain ca. 10(6) PFU/mL. However, even that highly diluted preparation significantly (p = < 0.05) reduced Y. pestis levels by ≥ 99.97%. Our data support the idea that Y. pestis phages may be useful for decontaminating various hard surfaces naturally- or intentionally-contaminated with Y. pestis.

Genomic comparison of Escherichia coli O104:H4 isolates from 2009 and 2011 reveals plasmid, and prophage heterogeneity, including shiga toxin encoding phage stx2.

In May of 2011, an enteroaggregative Escherichia coli O104:H4 strain that had acquired a Shiga toxin 2-converting phage caused a large outbreak of bloody diarrhea in Europe which was notable for its high prevalence of hemolytic uremic syndrome cases. Several studies have described the genomic inventory and phylogenies of strains associated with the outbreak and a collection of historical E. coli O104:H4 isolates using draft genome assemblies. We present the complete, closed genome sequences of an isolate from the 2011 outbreak (2011C-3493) and two isolates from cases of bloody diarrhea that occurred in the Republic of Georgia in 2009 (2009EL-2050 and 2009EL-2071). Comparative genome analysis indicates that, while the Georgian strains are the nearest neighbors to the 2011 outbreak isolates sequenced to date, structural and nucleotide-level differences are evident in the Stx2 phage genomes, the mer/tet antibiotic resistance island, and in the prophage and plasmid profiles of the strains, including a previously undescribed plasmid with homology to the pMT virulence plasmid of Yersinia pestis. In addition, multiphenotype analysis showed that 2009EL-2071 possessed higher resistance to polymyxin and membrane-disrupting agents. Finally, we show evidence by electron microscopy of the presence of a common phage morphotype among the European and Georgian strains and a second phage morphotype among the Georgian strains. The presence of at least two stx2 phage genotypes in host genetic backgrounds that may derive from a recent common ancestor of the 2011 outbreak isolates indicates that the emergence of stx2 phage-containing E. coli O104:H4 strains probably occurred more than once, or that the current outbreak isolates may be the result of a recent transfer of a new stx2 phage element into a pre-existing stx2-positive genetic background.