A gold standard dataset and evaluation of methods for lineage abundance estimation from wastewater.

During the SARS-CoV-2 pandemic, genome-based wastewater surveillance sequencing has been a powerful tool for public health to monitor circulating and emerging viral variants. As a medium, wastewater is very complex because of its mixed matrix nature, which makes the deconvolution of wastewater samples more difficult. Here we introduce a gold standard dataset constructed from synthetic viral control mixtures of known composition, spiked into a wastewater RNA matrix and sequenced on the Oxford Nanopore Technologies platform. We compare the performance of eight of the most commonly used deconvolution tools in identifying SARS-CoV-2 variants present in these mixtures. The software evaluated was primarily chosen for its relevance to the CDC wastewater surveillance reporting protocol, which until recently employed a pipeline that incorporates results from four deconvolution methods: Freyja, kallisto, Kraken 2/Bracken, and LCS. We also tested Lollipop, a deconvolution method used by the Swiss SARS-CoV-2 Sequencing Consortium, and three additional methods not used in the C-WAP pipeline: lineagespot, Alcov, and VaQuERo. We found that the commonly used software Freyja outperformed the other CDC pipeline tools in correct identification of lineages present in the control mixtures, and that the VaQuERo method was similarly accurate, with minor differences in the ability of the two methods to avoid false negatives and suppress false positives. Our results also provide insight into the effect of the tiling primer scheme and wastewater RNA extract matrix on viral sequencing and data deconvolution outcomes.

Performance evaluation of virus concentration methods for implementing SARS-CoV-2 wastewater based epidemiology emphasizing quick data turnaround.

Wastewater based epidemiology (WBE) has drawn significant attention as an early warning tool to detect and predict the trajectory of COVID-19 cases in a community, in conjunction with public health data. This means of monitoring for outbreaks has been used at municipal wastewater treatment centers to analyze COVID-19 trends in entire communities, as well as by universities and other community living environments to monitor COVID-19 spread in buildings. Sample concentration is crucial, especially when viral abundance in raw wastewater is below the threshold of detection by RT-qPCR analysis. We evaluated the performance of a rapid ultrafiltration-based virus concentration method using InnovaPrep Concentrating Pipette (CP) Select and compared this to the established electronegative membrane filtration (EMF) method. We evaluated sensitivity of SARS-CoV-2 quantification, surrogate virus recovery rate, and sample processing time. Results suggest that the CP Select concentrator is more efficient at concentrating SARS-CoV-2 from wastewater compared to the EMF method. About 25% of samples that tested negative when concentrated with the EMF method produced a positive signal with the CP Select protocol. Increased recovery of the surrogate virus control using the CP Select confirms this observation. We optimized the CP Select protocol by adding AVL lysis buffer and sonication, to increase the recovery of virus. Sonication increased Bovine Coronavirus (BCoV) recovery by 19%, which seems to compensate for viral loss during centrifugation. Filtration time decreases by approximately 30% when using the CP Select protocol, making this an optimal choice for building surveillance applications where quick turnaround time is necessary.

Characterization of Environmental and Cultivable Antibiotic-Resistant Microbial Communities Associated with Wastewater Treatment.

Bacterial resistance to antibiotics is a growing global concern, threatening human and environmental health, particularly among urban populations. Wastewater treatment plants (WWTPs) are thought to be "hotspots" for antibiotic resistance dissemination. The conditions of WWTPs, in conjunction with the persistence of commonly used antibiotics, may favor the selection and transfer of resistance genes among bacterial populations. WWTPs provide an important ecological niche to examine the spread of antibiotic resistance. We used heterotrophic plate count methods to identify phenotypically resistant cultivable portions of these bacterial communities and characterized the composition of the culturable subset of these populations. Resistant taxa were more abundant in raw sewage and wastewater before the biological aeration treatment stage. While some antibiotic-resistant bacteria (ARB) were detectable downstream of treated wastewater release, these organisms are not enriched relative to effluent-free upstream water, indicating efficient removal during treatment. Combined culture-dependent and -independent analyses revealed a stark difference in community composition between culturable fractions and the environmental source material, irrespective of culturing conditions. Higher proportions of the environmental populations were recovered than predicted by the widely accepted 1% culturability paradigm. These results represent baseline abundance and compositional data for ARB communities for reference in future studies addressing the dissemination of antibiotic resistance associated with urban wastewater treatment ecosystems.

Implementing building-level SARS-CoV-2 wastewater surveillance on a university campus.

The COVID-19 pandemic has been a source of ongoing challenges and presents an increased risk of illness in group environments, including jails, long-term care facilities, schools, and residential college campuses. Early reports that the SARS-CoV-2 virus was detectable in wastewater in advance of confirmed cases sparked widespread interest in wastewater-based epidemiology (WBE) as a tool for mitigation of COVID-19 outbreaks. One hypothesis was that wastewater surveillance might provide a cost-effective alternative to other more expensive approaches such as pooled and random testing of groups. In this paper, we report the outcomes of a wastewater surveillance pilot program at the University of North Carolina at Charlotte, a large urban university with a substantial population of students living in on-campus dormitories. Surveillance was conducted at the building level on a thrice-weekly schedule throughout the university's fall residential semester. In multiple cases, wastewater surveillance enabled the identification of asymptomatic COVID-19 cases that were not detected by other components of the campus monitoring program, which also included in-house contact tracing, symptomatic testing, scheduled testing of student athletes, and daily symptom reporting. In the context of all cluster events reported to the University community during the fall semester, wastewater-based testing events resulted in the identification of smaller clusters than were reported in other types of cluster events. Wastewater surveillance was able to detect single asymptomatic individuals in dorms with resident populations of 150-200. While the strategy described was developed for COVID-19, it is likely to be applicable to mitigation of future pandemics in universities and other group-living environments.

Corrigendum: Phylogeny of Vibrio vulnificus From the Analysis of the Core-Genome: Implications for Intra-Species Taxonomy.

[This corrects the article DOI: 10.3389/fmicb.2017.02613.].

The quantitative impact of read mapping to non-native reference genomes in comparative RNA-Seq studies.

Sequence read alignment to a reference genome is a fundamental step in many genomics studies. Accuracy in this fundamental step is crucial for correct interpretation of biological data. In cases where two or more closely related bacterial strains are being studied, a common approach is to simply map reads from all strains to a common reference genome, whether because there is no closed reference for some strains or for ease of comparison. The assumption is that the differences between bacterial strains are insignificant enough that the results of differential expression analysis will not be influenced by choice of reference. Genes that are common among the strains under study are used for differential expression analysis, while the remaining genes, which may fail to express in one sample or the other because they are simply absent, are analyzed separately. In this study, we investigate the practice of using a common reference in transcriptomic analysis. We analyze two multi-strain transcriptomic data sets that were initially presented in the literature as comparisons based on a common reference, but which have available closed genomic sequence for all strains, allowing a detailed examination of the impact of reference choice. We provide a method for identifying regions that are most affected by non-native alignments, leading to false positives in differential expression analysis, and perform an in depth analysis identifying the extent of expression loss. We also simulate several data sets to identify best practices for non-native reference use.

The impact of RNA secondary structure on read start locations on the Illumina sequencing platform.

High-throughput sequencing is subject to sequence dependent bias, which must be accounted for if researchers are to make precise measurements and draw accurate conclusions from their data. A widely studied source of bias in sequencing is the GC content bias, in which levels of GC content in a genomic region effect the number of reads produced during sequencing. Although some research has been performed on methods to correct for GC bias, there has been little effort to understand the underlying mechanism. The availability of sequencing protocols that target the specific location of structure in nucleic acid molecules enables us to investigate the underlying molecular origin of observed GC bias in sequencing. By applying a parallel analysis of RNA structure (PARS) protocol to bacterial genomes of varying GC content, we are able to observe the relationship between local RNA secondary structure and sequencing outcome, and to establish RNA secondary structure as the significant contributing factor to observed GC bias.

Comparative Genomic Analysis of Two Vibrio toranzoniae Strains with Different Virulence Capacity Reveals Clues on Its Pathogenicity for Fish.

Vibrio toranzoniae is a Gram-negative bacterium of the Splendidus clade within the Vibrio genus. V. toranzoniae was first isolated from healthy clams in Galicia (Spain) but recently was also identified associated to disease outbreaks of red conger eel in Chile. Experimental challenges showed that the Chilean isolates were able to produce fish mortalities but not the strains isolated from clams. The aim of the present study was to determine the differences at the genomic level between the type strain of the species (CECT 7225T) and the strain R17, isolated from red conger eel in Chile, which could explain their different virulent capacity. The genome-based comparison showed high homology between both strains but differences were observed in certain gene clusters that include some virulence factors. Among these, we found that iron acquisition systems and capsule synthesis genes were the main differential features between both genomes that could explain the differences in the pathogenicity of the strains. Besides, the studied genomes presented genomic islands and toxins, and the R17 strain presented CRISPR sequences that are absent on the type strain. Taken together, this analysis provided important insights into virulence factors of V. toranzoniae that will lead to a better understanding of the pathogenic process.

Phylogeny of Vibrio vulnificus from the Analysis of the Core-Genome: Implications for Intra-Species Taxonomy.

Vibrio vulnificus (Vv) is a multi-host pathogenic species currently subdivided into three biotypes (Bts). The three Bts are human-pathogens, but only Bt2 is also a fish-pathogen, an ability that is conferred by a transferable virulence-plasmid (pVvbt2). Here we present a phylogenomic analysis from the core genome of 80 Vv strains belonging to the three Bts recovered from a wide range of geographical and ecological sources. We have identified five well-supported phylogenetic groups or lineages (L). L1 comprises a mixture of clinical and environmental Bt1 strains, most of them involved in human clinical cases related to raw seafood ingestion. L2 is formed by a mixture of Bt1 and Bt2 strains from various sources, including diseased fish, and is related to the aquaculture industry. L3 is also linked to the aquaculture industry and includes Bt3 strains exclusively, mostly related to wound infections or secondary septicemia after farmed-fish handling. Lastly, L4 and L5 include a few strains of Bt1 associated with specific geographical areas. The phylogenetic trees for ChrI and II are not congruent to one another, which suggests that inter- and/or intra-chromosomal rearrangements have been produced along Vv evolution. Further, the phylogenetic trees for each chromosome and the virulence plasmid were also not congruent, which also suggests that pVvbt2 has been acquired independently by different clones, probably in fish farms. From all these clones, the one with zoonotic capabilities (Bt2-Serovar E) has successfully spread worldwide. Based on these results, we propose a new updated classification of the species based on phylogenetic lineages rather than on Bts, as well as the inclusion of all Bt2 strains in a pathovar with the particular ability to cause fish vibriosis, for which we suggest the name "piscis."

Draft Genome Sequence of Vibrio toranzoniae Strain CECT 7225T.

Vibrio toranzoniae(CECT 7225(T)) was isolated from healthy reared carpet shell clams in Galicia (Northwest Spain). In addition, this species has been recently identified as a potential pathogen of red conger eel in Chile. The draft genome sequence has 4.5 Mbp, a G+C content of 43.9%, and >3,800 protein-coding genes.

Impact of analytic provenance in genome analysis.

BACKGROUND: Many computational methods are available for assembly and annotation of newly sequenced microbial genomes. However, when new genomes are reported in the literature, there is frequently very little critical analysis of choices made during the sequence assembly and gene annotation stages. These choices have a direct impact on the biologically relevant products of a genomic analysis--for instance identification of common and differentiating regions among genomes in a comparison, or identification of enriched gene functional categories in a specific strain. Here, we examine the outcomes of different assembly and analysis steps in typical workflows in a comparison among strains of Vibrio vulnificus. RESULTS: Using six recently sequenced strains of V. vulnificus, we demonstrate the "alternate realities" of comparative genomics, and how they depend on the choice of a robust assembly method and accurate ab initio annotation. We apply several popular assemblers for paired-end Illumina data, and three well-regarded ab initio genefinders. We demonstrate significant differences in detected gene overlap among comparative genomics workflows that depend on these two steps. The divergence between workflows, even those using widely adopted methods, is obvious both at the single genome level and when a comparison is performed. In a typical example where multiple workflows are applied to the strain V. vulnificus CECT 4606, a workflow that uses the Velvet assembler and Glimmer gene finder identifies 3275 gene features, while a workflow that uses the Velvet assembler and the RAST annotation system identifies 5011 gene features. Only 3171 genes are identical between both workflows. When we examine 9 assembly/annotation workflow scenarios as input to a three-way genome comparison, differentiating genes and even differentially represented functional categories change significantly from scenario to scenario. CONCLUSIONS: Inconsistencies in genomic analysis can arise depending on the choices that are made during the assembly and annotation stages. These inconsistencies can have a significant impact on the interpretation of an individual genome's content. The impact is multiplied when comparison of content and function among multiple genomes is the goal. Tracking the analysis history of the data--its analytic provenance--is critical for reproducible analysis of genome data.

Transcriptome sequencing reveals the virulence and environmental genetic programs of Vibrio vulnificus exposed to host and estuarine conditions.

Vibrio vulnificus is a natural inhabitant of estuarine waters worldwide and is of medical relevance due to its ability to cause grievous wound infections and/or fatal septicemia. Genetic polymorphisms within the virulence-correlated gene (vcg) serve as a primary feature to distinguish clinical (C-) genotypes from environmental (E-) genotypes. C-genotypes demonstrate superior survival in human serum relative to E-genotypes, and genome comparisons have allowed for the identification of several putative virulence factors that could potentially aid C-genotypes in disease progression. We used RNA sequencing to analyze the transcriptome of C-genotypes exposed to human serum relative to seawater, which revealed two divergent genetic programs under these two conditions. In human serum, cells displayed a distinct "virulence profile" in which a number of putative virulence factors were upregulated, including genes involved in intracellular signaling, substrate binding and transport, toxin and exoenzyme production, and the heat shock response. Conversely, the "environmental profile" exhibited by cells in seawater revealed upregulation of transcription factors such as rpoS, rpoN, and iscR, as well as genes involved in intracellular signaling, chemotaxis, adherence, and biofilm formation. This dichotomous genetic switch appears to be largely governed by cyclic-di-GMP signaling, and remarkably resembles the dual life-style of V. cholerae as it transitions from host to environment. Furthermore, we found a "general stress response" module, known as the stressosome, to be upregulated in seawater. This signaling system has been well characterized in Gram-positive bacteria, however its role in V. vulnificus is not clear. We examined temporal gene expression patterns of the stressosome and found it to be upregulated in natural estuarine waters indicating that this system plays a role in sensing and responding to the environment. This study advances our understanding of gene regulation in V. vulnificus, and brings to the forefront a number of previously overlooked genetic networks.

The illusion of specific capture: surface and solution studies of suboptimal oligonucleotide hybridization.

BACKGROUND: Hybridization based assays and capture systems depend on the specificity of hybridization between a probe and its intended target. A common guideline in the construction of DNA microarrays, for instance, is that avoiding complementary stretches of more than 15 nucleic acids in a 50 or 60-mer probe will eliminate sequence specific cross-hybridization reactions. Here we present a study of the behavior of partially matched oligonucleotide pairs with complementary stretches starting well below this threshold complementarity length - in silico, in solution, and at the microarray surface. The modeled behavior of pairs of oligonucleotide probes and their targets suggests that even a complementary stretch of sequence 12 nt in length would give rise to specific cross-hybridization. We designed a set of binding partners to a 50-mer oligonucleotide containing complementary stretches from 6 nt to 21 nt in length. RESULTS: Solution melting experiments demonstrate that stable partial duplexes can form when only 12 bp of complementary sequence are present; surface hybridization experiments confirm that a signal close in magnitude to full-strength signal can be obtained from hybridization of a 12 bp duplex within a 50mer oligonucleotide. CONCLUSIONS: Microarray and other molecular capture strategies that rely on a 15 nt lower complementarity bound for eliminating specific cross-hybridization may not be sufficiently conservative.

Physico-chemical foundations underpinning microarray and next-generation sequencing experiments.

Hybridization of nucleic acids on solid surfaces is a key process involved in high-throughput technologies such as microarrays and, in some cases, next-generation sequencing (NGS). A physical understanding of the hybridization process helps to determine the accuracy of these technologies. The goal of a widespread research program is to develop reliable transformations between the raw signals reported by the technologies and individual molecular concentrations from an ensemble of nucleic acids. This research has inputs from many areas, from bioinformatics and biostatistics, to theoretical and experimental biochemistry and biophysics, to computer simulations. A group of leading researchers met in Ploen Germany in 2011 to discuss present knowledge and limitations of our physico-chemical understanding of high-throughput nucleic acid technologies. This meeting inspired us to write this summary, which provides an overview of the state-of-the-art approaches based on physico-chemical foundation to modeling of the nucleic acids hybridization process on solid surfaces. In addition, practical application of current knowledge is emphasized.

GenoSets: visual analytic methods for comparative genomics.

Many important questions in biology are, fundamentally, comparative, and this extends to our analysis of a growing number of sequenced genomes. Existing genomic analysis tools are often organized around literal views of genomes as linear strings. Even when information is highly condensed, these views grow cumbersome as larger numbers of genomes are added. Data aggregation and summarization methods from the field of visual analytics can provide abstracted comparative views, suitable for sifting large multi-genome datasets to identify critical similarities and differences. We introduce a software system for visual analysis of comparative genomics data. The system automates the process of data integration, and provides the analysis platform to identify and explore features of interest within these large datasets. GenoSets borrows techniques from business intelligence and visual analytics to provide a rich interface of interactive visualizations supported by a multi-dimensional data warehouse. In GenoSets, visual analytic approaches are used to enable querying based on orthology, functional assignment, and taxonomic or user-defined groupings of genomes. GenoSets links this information together with coordinated, interactive visualizations for both detailed and high-level categorical analysis of summarized data. GenoSets has been designed to simplify the exploration of multiple genome datasets and to facilitate reasoning about genomic comparisons. Case examples are included showing the use of this system in the analysis of 12 Brucella genomes. GenoSets software and the case study dataset are freely available at http://genosets.uncc.edu. We demonstrate that the integration of genomic data using a coordinated multiple view approach can simplify the exploration of large comparative genomic data sets, and facilitate reasoning about comparisons and features of interest.

Pyrosequencing-based comparative genome analysis of Vibrio vulnificus environmental isolates.

Between 1996 and 2006, the US Centers for Disease Control reported that the only category of food-borne infections increasing in frequency were those caused by members of the genus Vibrio. The gram-negative bacterium Vibrio vulnificus is a ubiquitous inhabitant of estuarine waters, and is the number one cause of seafood-related deaths in the US. Many V. vulnificus isolates have been studied, and it has been shown that two genetically distinct subtypes, distinguished by 16S rDNA and other gene polymorphisms, are associated predominantly with either environmental or clinical isolation. While local genetic differences between the subtypes have been probed, only the genomes of clinical isolates have so far been completely sequenced. In order to better understand V. vulnificus as an agent of disease and to identify the molecular components of its virulence mechanisms, we have completed whole genome shotgun sequencing of three diverse environmental genotypes using a pyrosequencing approach. V. vulnificus strain JY1305 was sequenced to a depth of 33×, and strains E64MW and JY1701 were sequenced to lesser depth, covering approximately 99.9% of each genome. We have performed a comparative analysis of these sequences against the previously published sequences of three V. vulnificus clinical isolates. We find that the genome of V. vulnificus is dynamic, with 1.27% of genes in the C-genotype genomes not found in the E- genotype genomes. We identified key genes that differentiate between the genomes of the clinical and environmental genotypes. 167 genes were found to be specifically associated with environmental genotypes and 278 genes with clinical genotypes. Genes specific to the clinical strains include components of sialic acid catabolism, mannitol fermentation, and a component of a Type IV secretory pathway VirB4, as well as several other genes with potential significance for human virulence. Genes specific to environmental strains included several that may have implications for the balance between self-preservation under stress and nutritional competence.

Comparative transcriptomics among floral organs of the basal eudicot Eschscholzia californica as reference for floral evolutionary developmental studies.

BACKGROUND: Molecular genetic studies of floral development have concentrated on several core eudicots and grasses (monocots), which have canalized floral forms. Basal eudicots possess a wider range of floral morphologies than the core eudicots and grasses and can serve as an evolutionary link between core eudicots and monocots, and provide a reference for studies of other basal angiosperms. Recent advances in genomics have enabled researchers to profile gene activities during floral development, primarily in the eudicot Arabidopsis thaliana and the monocots rice and maize. However, our understanding of floral developmental processes among the basal eudicots remains limited. RESULTS: Using a recently generated expressed sequence tag (EST) set, we have designed an oligonucleotide microarray for the basal eudicot Eschscholzia californica (California poppy). We performed microarray experiments with an interwoven-loop design in order to characterize the E. californica floral transcriptome and to identify differentially expressed genes in flower buds with pre-meiotic and meiotic cells, four floral organs at preanthesis stages (sepals, petals, stamens and carpels), developing fruits, and leaves. CONCLUSIONS: Our results provide a foundation for comparative gene expression studies between eudicots and basal angiosperms. We identified whorl-specific gene expression patterns in E. californica and examined the floral expression of several gene families. Interestingly, most E. californica homologs of Arabidopsis genes important for flower development, except for genes encoding MADS-box transcription factors, show different expression patterns between the two species. Our comparative transcriptomics study highlights the unique evolutionary position of E. californica compared with basal angiosperms and core eudicots.