Petabase-scale sequence alignment catalyses viral discovery.

Public databases contain a planetary collection of nucleic acid sequences, but their systematic exploration has been inhibited by a lack of efficient methods for searching this corpus, which (at the time of writing) exceeds 20 petabases and is growing exponentially1. Here we developed a cloud computing infrastructure, Serratus, to enable ultra-high-throughput sequence alignment at the petabase scale. We searched 5.7 million biologically diverse samples (10.2 petabases) for the hallmark gene RNA-dependent RNA polymerase and identified well over 105 novel RNA viruses, thereby expanding the number of known species by roughly an order of magnitude. We characterized novel viruses related to coronaviruses, hepatitis delta virus and huge phages, respectively, and analysed their environmental reservoirs. To catalyse the ongoing revolution of viral discovery, we established a free and comprehensive database of these data and tools. Expanding the known sequence diversity of viruses can reveal the evolutionary origins of emerging pathogens and improve pathogen surveillance for the anticipation and mitigation of future pandemics.

Photoacoustic discrimination of antibiotic-resistant and sensitive Staphylococcus aureus isolates.

OBJECTIVES: Bacteremia is a serious and potentially lethal condition. Staphylococcus aureus is a leading cause of bacteremia and methicillin-resistant S. aureus (MRSA) accounts for more than a third of the cases. Compared to methicillin-sensitive S. aureus, MRSA is more than twice as likely to be fatal. Furthermore, subpopulations of seemingly isogenic bacteria may exhibit a range of susceptibilities, often called heterogenous resistance. These heterogeneous antibiotic-resistant infections are often misdiagnosed as hospital-acquired secondary infections because there are no clinically used tests that can differentiate between homogeneous and heterogeneous antibiotic resistance. We describe the development and proof of concept of rapid bacterial identification using photoacoustic flow cytometry and labeled bacteriophages with the characterization and differentiation of heterogeneous antibiotic-resistant bacterial infections. METHODS: In photoacoustic flow cytometry, pulsed laser light is delivered to a sample flowing past a focused transducer and particles that absorb laser light create an acoustic response. Optically labeled bacteriophage are added to a bacterial mixture that flows through the photoacoustic chamber. The presence of target bacteria is determined by bound labeled phage which are detected photoacoustically. Incubation of bacterial samples in the presence and absence of the antibiotic daptomycin creates a difference in bacterial cell numbers that is quantified using photoacoustic flow cytometry. RESULTS: Four clinical isolates were tested in the presence and absence of daptomycin. Photoacoustic events for each isolate were recorded and compared to growth curves. Samples treated with daptomycin fell into three categories: resistant, susceptible, and heterogeneous resistant. CONCLUSIONS: Here we show a method to determine the presence of bacteria as a marker for bloodstream infection level and antibiotic sensitivity in less than 4 hours. Additionally, these results show an ability to identify heterogeneous resistant strains that are often misidentified.

Predicting Metastasis in Melanoma by Enumerating Circulating Tumor Cells Using Photoacoustic Flow Cytometry.

BACKGROUND AND OBJECTIVES: Enumerating circulating tumor cells has been used as a method of monitoring progression of various cancers. Various methods for detecting circulating melanoma cells (CMCs) have been reported, but none has had sufficient sensitivity to determine if the presence of rare CMCs in the blood of Stage I-III melanoma patients predicts if those patients eventually develop metastatic disease. STUDY DESIGN: We quantified CMCs in serial blood samples from 38 early stage melanoma patients to determine if CMC numbers predict development of metastatic melanoma. CMCs were enumerated using a photoacoustic flow cytometric detection system that uses a laser to induce high frequency acoustic signals in pigmented CMCs. RESULTS: We observed that detection of greater than 2 CMCs/ml of blood from patients with Stage I-III melanoma predicts metastatic disease. Of the 11 patients we studied who had two or fewer CMCs detected at all time points tested, none progressed to metastatic disease over a mean follow-up of 1288 days. In contrast, 18 of the 27 patients (67%) having more than 2 CMCs/ml at one or more time points progressed to metastatic disease over a mean follow-up of 850 days. CONCLUSIONS: Photoacoustic flow cytometry can detect rare CMCs in the blood of Stage I-III melanoma patients and detectionof these cells is predictive of subsequent development of metastatic disease. Lasers Surg. Med.

Updating the 97% identity threshold for 16S ribosomal RNA OTUs.

Motivation: The 16S ribosomal RNA (rRNA) gene is widely used to survey microbial communities. Sequences are often clustered into Operational Taxonomic Units (OTUs) as proxies for species. The canonical clustering threshold is 97% identity, which was proposed in 1994 when few 16S rRNA sequences were available, motivating a reassessment on current data. Results: Using a large set of high-quality 16S rRNA sequences from finished genomes, I assessed the correspondence of OTUs to species for five representative clustering algorithms using four accuracy metrics. All algorithms had comparable accuracy when tuned to a given metric. Optimal identity thresholds were ∼99% for full-length sequences and ∼100% for the V4 hypervariable region. Availability and implementation: Reference sequences and source code are provided in the Supplementary Material. Supplementary information: Supplementary data are available at Bioinformatics online.

Refugia under threat: Mass bleaching of coral assemblages in high-latitude eastern Australia.

Environmental anomalies that trigger adverse physiological responses and mortality are occurring with increasing frequency due to climate change. At species' range peripheries, environmental anomalies are particularly concerning because species often exist at their environmental tolerance limits and may not be able to migrate to escape unfavourable conditions. Here, we investigated the bleaching response and mortality of 14 coral genera across high-latitude eastern Australia during a global heat stress event in 2016. We evaluated whether the severity of assemblage-scale and genus-level bleaching responses was associated with cumulative heat stress and/or local environmental history, including long-term mean temperatures during the hottest month of each year (SSTLTMAX ), and annual fluctuations in water temperature (SSTVAR ) and solar irradiance (PARZVAR ). The most severely-bleached genera included species that were either endemic to the region (Pocillopora aliciae) or rare in the tropics (e.g. Porites heronensis). Pocillopora spp., in particular, showed high rates of immediate mortality. Bleaching severity of Pocillopora was high where SSTLTMAX was low or PARZVAR was high, whereas bleaching severity of Porites was directly associated with cumulative heat stress. While many tropical Acropora species are extremely vulnerable to bleaching, the Acropora species common at high latitudes, such as A. glauca and A. solitaryensis, showed little incidence of bleaching and immediate mortality. Two other regionally-abundant genera, Goniastrea and Turbinaria, were also largely unaffected by the thermal anomaly. The severity of assemblage-scale bleaching responses was poorly explained by the environmental parameters we examined. Instead, the severity of assemblage-scale bleaching was associated with local differences in species abundance and taxon-specific bleaching responses. The marked taxonomic disparity in bleaching severity, coupled with high mortality of high-latitude endemics, point to climate-driven simplification of assemblage structures and progressive homogenisation of reef functions at these high-latitude locations.

Defining the core Arabidopsis thaliana root microbiome.

Land plants associate with a root microbiota distinct from the complex microbial community present in surrounding soil. The microbiota colonizing the rhizosphere (immediately surrounding the root) and the endophytic compartment (within the root) contribute to plant growth, productivity, carbon sequestration and phytoremediation. Colonization of the root occurs despite a sophisticated plant immune system, suggesting finely tuned discrimination of mutualists and commensals from pathogens. Genetic principles governing the derivation of host-specific endophyte communities from soil communities are poorly understood. Here we report the pyrosequencing of the bacterial 16S ribosomal RNA gene of more than 600 Arabidopsis thaliana plants to test the hypotheses that the root rhizosphere and endophytic compartment microbiota of plants grown under controlled conditions in natural soils are sufficiently dependent on the host to remain consistent across different soil types and developmental stages, and sufficiently dependent on host genotype to vary between inbred Arabidopsis accessions. We describe different bacterial communities in two geochemically distinct bulk soils and in rhizosphere and endophytic compartments prepared from roots grown in these soils. The communities in each compartment are strongly influenced by soil type. Endophytic compartments from both soils feature overlapping, low-complexity communities that are markedly enriched in Actinobacteria and specific families from other phyla, notably Proteobacteria. Some bacteria vary quantitatively between plants of different developmental stage and genotype. Our rigorous definition of an endophytic compartment microbiome should facilitate controlled dissection of plant-microbe interactions derived from complex soil communities.

UPARSE: highly accurate OTU sequences from microbial amplicon reads.

Amplified marker-gene sequences can be used to understand microbial community structure, but they suffer from a high level of sequencing and amplification artifacts. The UPARSE pipeline reports operational taxonomic unit (OTU) sequences with ≤1% incorrect bases in artificial microbial community tests, compared with >3% incorrect bases commonly reported by other methods. The improved accuracy results in far fewer OTUs, consistently closer to the expected number of species in a community.

Error filtering, pair assembly and error correction for next-generation sequencing reads.

MOTIVATION: Next-generation sequencing produces vast amounts of data with errors that are difficult to distinguish from true biological variation when coverage is low. RESULTS: We demonstrate large reductions in error frequencies, especially for high-error-rate reads, by three independent means: (i) filtering reads according to their expected number of errors, (ii) assembling overlapping read pairs and (iii) for amplicon reads, by exploiting unique sequence abundances to perform error correction. We also show that most published paired read assemblers calculate incorrect posterior quality scores. AVAILABILITY AND IMPLEMENTATION: These methods are implemented in the USEARCH package. Binaries are freely available at http://drive5.com/usearch. CONTACT: robert@drive5.com SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Correction: Ribovirus classification by a polymerase barcode sequence.

[This corrects the article DOI: 10.7717/peerj.14055.].

Viroid-like colonists of human microbiomes.

Here, we describe the "Obelisks," a previously unrecognised class of viroid-like elements that we first identified in human gut metatranscriptomic data. "Obelisks" share several properties: (i) apparently circular RNA ~1kb genome assemblies, (ii) predicted rod-like secondary structures encompassing the entire genome, and (iii) open reading frames coding for a novel protein superfamily, which we call the "Oblins". We find that Obelisks form their own distinct phylogenetic group with no detectable sequence or structural similarity to known biological agents. Further, Obelisks are prevalent in tested human microbiome metatranscriptomes with representatives detected in ~7% of analysed stool metatranscriptomes (29/440) and in ~50% of analysed oral metatranscriptomes (17/32). Obelisk compositions appear to differ between the anatomic sites and are capable of persisting in individuals, with continued presence over >300 days observed in one case. Large scale searches identified 29,959 Obelisks (clustered at 90% nucleotide identity), with examples from all seven continents and in diverse ecological niches. From this search, a subset of Obelisks are identified to code for Obelisk-specific variants of the hammerhead type-III self-cleaving ribozyme. Lastly, we identified one case of a bacterial species (Streptococcus sanguinis) in which a subset of defined laboratory strains harboured a specific Obelisk RNA population. As such, Obelisks comprise a class of diverse RNAs that have colonised, and gone unnoticed in, human, and global microbiomes.

Known phyla dominate the Tara Oceans RNA virome.

A recent study proposed five new RNA virus phyla, two of which, 'Taraviricota' and 'Arctiviricota', were stated to be 'dominant in the oceans'. However, the study's assignments classify 28,353 putative RdRp-containing contigs to known phyla but only 886 (2.8%) to the five proposed new phyla combined. I re-mapped the reads to the contigs, finding that known phyla also account for a large majority (93.8%) of reads according to the study's classifications, and that contigs originally assigned to 'Arctiviricota' accounted for only a tiny fraction (0.01%) of reads from Arctic Ocean samples. Performing my own virus identification and classifications, I found that 99.95 per cent of reads could be assigned to known phyla. The most abundant species was Beihai picorna-like virus 34 (15% of reads), and the most abundant order-like cluster was classified as Picornavirales (45% of reads). Sequences in the claimed new phylum 'Pomiviricota' were placed inside a phylogenetic tree for established order Durnavirales with 100 per cent confidence. Moreover, two contigs assigned to the proposed phylum 'Taraviricota' were found to have high-identity alignments to dinoflagellate proteins, tentatively identifying this group of RdRp-like sequences as deriving from non-viral transcripts. Together, these results comprehensively contradict the claim that new phyla dominate the data.

Photoacoustic Flow Cytometry Using Functionalized Microspheres for Selective Detection of Bacteria.

Photoacoustic flow cytometry is a method to detect rare analytes in fluids. We developed photoacoustic flow cytometry to detect pathological cells in body fluids, such as circulating tumor cells or bacteria in blood. In order to induce specific optical absorption in bacteria, we use modified bacteriophage that precisely target bacterial species or subspecies for rapid identification. In order to reduce detection variability and to halt the lytic lifescycle that results in lysis of the bacteria, we attached dyed latex microspheres to the tail fibers of bacteriophage that retained the bacterial recognition binding sites. We tested these microsphere complexes using Salmonella enterica (Salmonella) and Escherichia coli (E. coli) bacteria and found robust and specific detection of targeted bacteria. In our work we used LT2, a strain of Salmonella, against K12, a strain of E. coli. Using Det7, a bacteriophage that binds to LT2 and not to K12, we detected an average of 109.3±9.0 of LT2 versus 2.0±1.7 of K12 using red microspheres and 86.7±13.2 of LT2 versus 0.3±0.6 of K12 using blue microspheres. These results confirmed our ability to selectively detect bacterial species using photoacoustic flow cytometry.

Differentiating methicillin resistant and susceptible Staphylococcus aureus from ocular infections using photoacoustic labeling.

Introduction: Antibiotic resistance in bacterial species constitutes a growing problem in the clinical management of infections. Not only does it limit therapeutic options, but application of ineffective antibiotics allows resistant species to progress prior to prescribing more effective treatment to patients. Methicillin resistance in Staphylococcus aureus is a major problem in clinical infections as it is the most common hospital acquired infection. Methods: We developed a photoacoustic flow cytometer using engineered bacteriophage as probes for rapid determination of methicillin resistance in Staphylococcus aureus with thirteen clinical samples obtained from keratitis patients. This method irradiates cells under flow with 532 nm laser light and selectively generates acoustic waves in labeled bacterial cells, thus enabling detection and enumeration of them. Staphylococcus aureus isolates were classified from culture isolation as either methicillin resistant or susceptible using cefoxitin disk diffusion testing. The photoacoustic method enumerates bacterial cells before and after treatment with antibiotics. Decreasing counts of bacteria after treatment indicate susceptible strains. We quantified the bacterial cells in the treated and untreated samples. Results: Using k-means clustering on the data, we achieved 100% concordance with the classification of Staphylococcus aureus resistance using culture. Discussion: Photoacoustics can be used to differentiate methicillin resistant and susceptible strains of bacteria from ocular infections. This method may be generalized to other bacterial species using appropriate bacteriophages and testing for resistance using other antibiotics.

Hybrids of RNA viruses and viroid-like elements replicate in fungi.

Earth's life may have originated as self-replicating RNA, and it has been argued that RNA viruses and viroid-like elements are remnants of such pre-cellular RNA world. RNA viruses are defined by linear RNA genomes encoding an RNA-dependent RNA polymerase (RdRp), whereas viroid-like elements consist of small, single-stranded, circular RNA genomes that, in some cases, encode paired self-cleaving ribozymes. Here we show that the number of candidate viroid-like elements occurring in geographically and ecologically diverse niches is much higher than previously thought. We report that, amongst these circular genomes, fungal ambiviruses are viroid-like elements that undergo rolling circle replication and encode their own viral RdRp. Thus, ambiviruses are distinct infectious RNAs showing hybrid features of viroid-like RNAs and viruses. We also detected similar circular RNAs, containing active ribozymes and encoding RdRps, related to mitochondrial-like fungal viruses, highlighting fungi as an evolutionary hub for RNA viruses and viroid-like elements. Our findings point to a deep co-evolutionary history between RNA viruses and subviral elements and offer new perspectives in the origin and evolution of primordial infectious agents, and RNA life.

UCHIME improves sensitivity and speed of chimera detection.

MOTIVATION: Chimeric DNA sequences often form during polymerase chain reaction amplification, especially when sequencing single regions (e.g. 16S rRNA or fungal Internal Transcribed Spacer) to assess diversity or compare populations. Undetected chimeras may be misinterpreted as novel species, causing inflated estimates of diversity and spurious inferences of differences between populations. Detection and removal of chimeras is therefore of critical importance in such experiments. RESULTS: We describe UCHIME, a new program that detects chimeric sequences with two or more segments. UCHIME either uses a database of chimera-free sequences or detects chimeras de novo by exploiting abundance data. UCHIME has better sensitivity than ChimeraSlayer (previously the most sensitive database method), especially with short, noisy sequences. In testing on artificial bacterial communities with known composition, UCHIME de novo sensitivity is shown to be comparable to Perseus. UCHIME is >100× faster than Perseus and >1000× faster than ChimeraSlayer. CONTACT: robert@drive5.com AVAILABILITY: Source, binaries and data: http://drive5.com/uchime. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Quality measures for protein alignment benchmarks.

Multiple protein sequence alignment methods are central to many applications in molecular biology. These methods are typically assessed on benchmark datasets including BALIBASE, OXBENCH, PREFAB and SABMARK, which are important to biologists in making informed choices between programs. In this article, annotations of domain homology and secondary structure are used to define new measures of alignment quality and are used to make the first systematic, independent evaluation of these benchmarks. These measures indicate sensitivity and specificity while avoiding the ambiguous residue correspondences and arbitrary distance cutoffs inherent to structural superpositions. Alignments by selected methods that indicate high-confidence columns (ALIGN-M, DIALIGN-T, FSA and MUSCLE) are also assessed. Fold space coverage and effective benchmark database sizes are estimated by reference to domain annotations, and significant redundancy is found in all benchmarks except SABMARK. Questionable alignments are found in all benchmarks, especially in BALIBASE where 87% of sequences have unknown structure, 20% of columns contain different folds according to SUPERFAMILY and 30% of 'core block' columns have conflicting secondary structure according to DSSP. A careful analysis of current protein multiple alignment benchmarks calls into question their ability to determine reliable algorithm rankings.

Characterization and distribution of retrotransposons and simple sequence repeats in the bovine genome.

Interspersed repeat composition and distribution in mammals have been best characterized in the human and mouse genomes. The bovine genome contains typical eutherian mammal repeats, but also has a significant number of long interspersed nuclear element RTE (BovB) elements proposed to have been horizontally transferred from squamata. Our analysis of the BovB repeats has indicated that only a few of them are currently likely to retrotranspose in cattle. However, bovine L1 repeats (L1 BT) have many likely active copies. Comparison of substitution rates for BovB and L1 BT indicates that L1 BT is a younger repeat family than BovB. In contrast to mouse and human, L1 occurrence is not negatively correlated with G+C content. However, BovB, Bov A2, ART2A, and Bov-tA are negatively correlated with G+C, although Bov-tAs correlation is weaker. Also, by performing genome wide correlation analysis of interspersed and simple sequence repeats, we have identified genome territories by repeat content that appear to define ancestral vs. ruminant-specific genomic regions. These ancestral regions, enriched with L2 and MIR repeats, are largely conserved between bovine and human.

Ribovirus classification by a polymerase barcode sequence.

RNA viruses encoding a polymerase gene (riboviruses) dominate the known eukaryotic virome. High-throughput sequencing is revealing a wealth of new riboviruses known only from sequence, precluding classification by traditional taxonomic methods. Sequence classification is often based on polymerase sequences, but standardised methods to support this approach are currently lacking. To address this need, we describe the polymerase palmprint, a segment of the palm sub-domain robustly delineated by well-conserved catalytic motifs. We present an algorithm, Palmscan, which identifies palmprints in nucleotide and amino acid sequences; PALMdb, a collection of palmprints derived from public sequence databases; and palmID, a public website implementing palmprint identification, search, and annotation. Together, these methods demonstrate a proof-of-concept workflow for high-throughput characterisation of RNA viruses, paving the path for the continued rapid growth in RNA virus discovery anticipated in the coming decade.

Muscle5: High-accuracy alignment ensembles enable unbiased assessments of sequence homology and phylogeny.

Multiple sequence alignments are widely used to infer evolutionary relationships, enabling inferences of structure, function, and phylogeny. Standard practice is to construct one alignment by some preferred method and use it in further analysis; however, undetected alignment bias can be problematic. I describe Muscle5, a novel algorithm which constructs an ensemble of high-accuracy alignment with diverse biases by perturbing a hidden Markov model and permuting its guide tree. Confidence in an inference is assessed as the fraction of the ensemble which supports it. Applied to phylogenetic tree estimation, I show that ensembles can confidently resolve topologies with low bootstrap according to standard methods, and conversely that some topologies with high bootstraps are incorrect. Applied to the phylogeny of RNA viruses, ensemble analysis shows that recently adopted taxonomic phyla are probably polyphyletic. Ensemble analysis can improve confidence assessment in any inference from an alignment.

GMP Manufacturing and IND-Enabling Studies of a Recombinant Hyperimmune Globulin Targeting SARS-CoV-2.

Conventionally, hyperimmune globulin drugs manufactured from pooled immunoglobulins from vaccinated or convalescent donors have been used in treating infections where no treatment is available. This is especially important where multi-epitope neutralization is required to prevent the development of immune-evading viral mutants that can emerge upon treatment with monoclonal antibodies. Using microfluidics, flow sorting, and a targeted integration cell line, a first-in-class recombinant hyperimmune globulin therapeutic against SARS-CoV-2 (GIGA-2050) was generated. Using processes similar to conventional monoclonal antibody manufacturing, GIGA-2050, comprising 12,500 antibodies, was scaled-up for clinical manufacturing and multiple development/tox lots were assessed for consistency. Antibody sequence diversity, cell growth, productivity, and product quality were assessed across different manufacturing sites and production scales. GIGA-2050 was purified and tested for good laboratory procedures (GLP) toxicology, pharmacokinetics, and in vivo efficacy against natural SARS-CoV-2 infection in mice. The GIGA-2050 master cell bank was highly stable, producing material at consistent yield and product quality up to >70 generations. Good manufacturing practices (GMP) and development batches of GIGA-2050 showed consistent product quality, impurity clearance, potency, and protection in an in vivo efficacy model. Nonhuman primate toxicology and pharmacokinetics studies suggest that GIGA-2050 is safe and has a half-life similar to other recombinant human IgG1 antibodies. These results supported a successful investigational new drug application for GIGA-2050. This study demonstrates that a new class of drugs, recombinant hyperimmune globulins, can be manufactured consistently at the clinical scale and presents a new approach to treating infectious diseases that targets multiple epitopes of a virus.