Long-read powered viral metagenomics in the oligotrophic Sargasso Sea

Dominant microorganisms of the Sargasso Sea are key drivers of the global carbon cycle. However, associated viruses that shape microbial community structure and function are not well characterised. Here, we combined short and long read sequencing to survey Sargasso Sea phage communities in virus- and cellular fractions at viral maximum (80 m) and mesopelagic (200 m) depths. We identified 2,301 Sargasso Sea phage populations from 186 genera. Over half of the phage populations identified here lacked representation in global ocean viral metagenomes, whilst 177 of the 186 identified genera lacked representation in genomic databases of phage isolates. Viral fraction and cell-associated viral communities were decoupled, indicating viral turnover occurred across periods longer than the sampling period of three days. Inclusion of long-read data was critical for capturing the breadth of viral diversity. Phage isolates that infect the dominant bacterial taxa Prochlorococcus and Pelagibacter, usually regarded as cosmopolitan and abundant, were poorly represented.

Resistance to bacteriophage incurs a cost to virulence in drug-resistant Acinetobacter baumannii

Introduction . Acinetobacter baumannii is a critical priority pathogen for novel antimicrobials (World Health Organization) because of the rise in nosocomial infections and its ability to evolve resistance to last resort antibiotics. A. baumannii is thus a priority target for phage therapeutics. Two strains of a novel, virulent bacteriophage (LemonAid and Tonic) able to infect carbapenem-resistant A. baumannii (strain NCTC 13420), were isolated from environmental water samples collected through a citizen science programme.Gap statement. Phage-host coevolution can lead to emergence of host resistance, with a concomitant reduction in the virulence of host bacteria; a potential benefit to phage therapy applications.Methodology. In vitro and in vivo assays, genomics and microscopy techniques were used to characterize the phages; determine mechanisms and impact of phage resistance on host virulence, and the efficacy of the phages against A. baumannii.Results. A. baumannii developed resistance to both viruses, LemonAid and Tonic. Resistance came at a cost to virulence, with the resistant variants causing significantly reduced mortality in a Galleria mellonella larval in vivo model. A replicated 8 bp insertion increased in frequency (~40 % higher frequency than in the wild-type) within phage-resistant A. baumannii mutants, putatively resulting in early truncation of a protein of unknown function. Evidence from comparative genomics and an adsorption assay suggests this protein acts as a novel phage receptor site in A. baumannii. We find no evidence linking resistance to changes in capsule structure, a known virulence factor. LemonAid efficiently suppressed growth of A. baumanni in vitro across a wide range of titres. However, in vivo, while survival of A. baumannii infected larvae significantly increased with both remedial and prophylactic treatment with LemonAid (107 p.f.u. ml-1), the effect was weak and not sufficient to save larvae from morbidity and mortality.Conclusion. While LemonAid and Tonic did not prove effective as a treatment in a Galleria larvae model, there is potential to harness their ability to attenuate virulence in drug-resistant A. baumannii.

Influence of host phylogeny and water physicochemistry on microbial assemblages of the fish skin microbiome

The skin of fish contains a diverse microbiota that has symbiotic functions with the host, facilitating pathogen exclusion, immune system priming, and nutrient degradation. The composition of fish skin microbiomes varies across species and in response to a variety of stressors, however, there has been no systematic analysis across these studies to evaluate how these factors shape fish skin microbiomes. Here, we examined 1922 fish skin microbiomes from 36 studies that included 98 species and nine rearing conditions to investigate associations between fish skin microbiome, fish species, and water physiochemical factors. Proteobacteria, particularly the class Gammaproteobacteria, were present in all marine and freshwater fish skin microbiomes. Acinetobacter, Aeromonas, Ralstonia, Sphingomonas and Flavobacterium were the most abundant genera within freshwater fish skin microbiomes, and Alteromonas, Photobacterium, Pseudoalteromonas, Psychrobacter and Vibrio were the most abundant in saltwater fish. Our results show that different culturing (rearing) environments have a small but significant effect on the skin bacterial community compositions. Water temperature, pH, dissolved oxygen concentration, and salinity significantly correlated with differences in beta-diversity but not necessarily alpha-diversity. To improve study comparability on fish skin microbiomes, we provide recommendations for approaches to the analyses of sequencing data and improve study reproducibility.

The Citizen Phage Library: Rapid Isolation of Phages for the Treatment of Antibiotic Resistant Infections in the UK

Antimicrobial resistance poses one of the greatest threats to global health and there is an urgent need for new therapeutic options. Phages are viruses that infect and kill bacteria and phage therapy could provide a valuable tool for the treatment of multidrug-resistant infections. In this study, water samples collected by citizen scientists as part of the Citizen Phage Library (CPL) project, and wastewater samples from the Environment Agency yielded phages with activity against clinical strains Klebsiella pneumoniae BPRG1484 and Enterobacter cloacae BPRG1482. A total of 169 and 163 phages were found for K. pneumoniae and E. cloacae, respectively, within four days of receiving the strains. A third strain (Escherichia coli BPRG1486) demonstrated cross-reactivity with 42 E. coli phages already held in the CPL collection. Seed lots were prepared for four K. pneumoniae phages and a cocktail combining these phages was found to reduce melanisation in a Galleria mellonella infection model. The resources and protocols utilised by the Citizen Phage Library enabled the rapid isolation and characterisation of phages targeted against multiple strains. In the future, within a clearly defined regulatory framework, phage therapy could be made available on a named-patient basis within the UK.

Suboxic DOM is bioavailable to surface prokaryotes in a simulated overturn of an oxygen minimum zone, Devil's Hole, Bermuda

Oxygen minimum zones (OMZs) are expanding due to increased sea surface temperatures, subsequent increased oxygen demand through respiration, reduced oxygen solubility, and thermal stratification driven in part by anthropogenic climate change. Devil's Hole, Bermuda is a model ecosystem to study OMZ microbial biogeochemistry because the formation and subsequent overturn of the suboxic zone occur annually. During thermally driven stratification, suboxic conditions develop, with organic matter and nutrients accumulating at depth. In this study, the bioavailability of the accumulated dissolved organic carbon (DOC) and the microbial community response to reoxygenation of suboxic waters was assessed using a simulated overturn experiment. The surface inoculated prokaryotic community responded to the deep (formerly suboxic) 0.2 μm filtrate with cell densities increasing 2.5-fold over 6 days while removing 5 μmol L-1 of DOC. After 12 days, the surface community began to shift, and DOC quality became less diagenetically altered along with an increase in SAR202, a Chloroflexi that can degrade recalcitrant dissolved organic matter (DOM). Labile DOC production after 12 days coincided with an increase of Nitrosopumilales, a chemoautotrophic ammonia oxidizing archaea (AOA) that converts ammonia to nitrite based on the ammonia monooxygenase (amoA) gene copy number and nutrient data. In comparison, the inoculation of the deep anaerobic prokaryotic community into surface 0.2 μm filtrate demonstrated a die-off of 25.5% of the initial inoculum community followed by a 1.5-fold increase in cell densities over 6 days. Within 2 days, the prokaryotic community shifted from a Chlorobiales dominated assemblage to a surface-like heterotrophic community devoid of Chlorobiales. The DOM quality changed to less diagenetically altered material and coincided with an increase in the ribulose-1,5-bisphosphate carboxylase/oxygenase form I (cbbL) gene number followed by an influx of labile DOM. Upon reoxygenation, the deep DOM that accumulated under suboxic conditions is bioavailable to surface prokaryotes that utilize the accumulated DOC initially before switching to a community that can both produce labile DOM via chemoautotrophy and degrade the more recalcitrant DOM.

Novel pelagiphage isolate Polarivirus skadi is a polar specialist that dominates SAR11-associated bacteriophage communities at high latitudes

The SAR11 clade are the most abundant members of surface marine bacterioplankton and a critical component of global biogeochemical cycles. Similarly, pelagiphages that infect SAR11 are ubiquitous and highly abundant in the oceans. Pelagiphages are predicted to shape SAR11 community structures and increase carbon turnover throughout the oceans. Yet, ecological drivers of host and niche specificity of pelagiphage populations are poorly understood. Here we report the global distribution of a novel pelagiphage called "Polarivirus skadi", which is the sole representative of a novel genus. P. skadi was isolated from the Western English Channel using a cold-water ecotype of SAR11 as bait. P. skadi is closely related to the globally dominant pelagiphage HTVC010P. Along with other HTVC010P-type viruses, P. skadi belongs to a distinct viral family within the order Caudovirales, for which we propose the name Ubiqueviridae. Metagenomic read recruitment identified P. skadi as one of the most abundant pelagiphages on Earth. P. skadi is a polar specialist, replacing HTVC010P at high latitudes. Experimental evaluation of P. skadi host range against cold- and warm-water SAR11 ecotypes supported cold-water specialism. Relative abundance of P. skadi in marine metagenomes correlated negatively with temperature, and positively with nutrients, available oxygen, and chlorophyll concentrations. In contrast, relative abundance of HTVC010P correlated negatively with oxygen and positively with salinity, with no significant correlation to temperature. The majority of other pelagiphages were scarce in most marine provinces, with a few representatives constrained to discrete ecological niches. Our results suggest that pelagiphage populations persist within a global viral seed bank, with environmental parameters and host availability selecting for a few ecotypes that dominate ocean viromes.

Lower viral evolutionary pressure under stable versus fluctuating conditions in subzero Arctic brines

BACKGROUND

Climate change threatens Earth's ice-based ecosystems which currently offer archives and eco-evolutionary experiments in the extreme. Arctic cryopeg brine (marine-derived, within permafrost) and sea ice brine, similar in subzero temperature and high salinity but different in temporal stability, are inhabited by microbes adapted to these extreme conditions. However, little is known about their viruses (community composition, diversity, interaction with hosts, or evolution) or how they might respond to geologically stable cryopeg versus fluctuating sea ice conditions.

RESULTS

We used long- and short-read viromics and metatranscriptomics to study viruses in Arctic cryopeg brine, sea ice brine, and underlying seawater, recovering 11,088 vOTUs (~species-level taxonomic unit), a 4.4-fold increase of known viruses in these brines. More specifically, the long-read-powered viromes doubled the number of longer (≥25 kb) vOTUs generated and recovered more hypervariable regions by >5-fold compared to short-read viromes. Distribution assessment, by comparing to known viruses in public databases, supported that cryopeg brine viruses were of marine origin yet distinct from either sea ice brine or seawater viruses, while 94% of sea ice brine viruses were also present in seawater. A virus-encoded, ecologically important exopolysaccharide biosynthesis gene was identified, and many viruses (~half of metatranscriptome-inferred "active" vOTUs) were predicted as actively infecting the dominant microbial genera Marinobacter and Polaribacter in cryopeg and sea ice brines, respectively. Evolutionarily, microdiversity (intra-species genetic variations) analyses suggested that viruses within the stable cryopeg brine were under significantly lower evolutionary pressures than those in the fluctuating sea ice environment, while many sea ice brine virus-tail genes were under positive selection, indicating virus-host co-evolutionary arms races.

CONCLUSIONS

Our results confirmed the benefits of long-read-powered viromics in understanding the environmental virosphere through significantly improved genomic recovery, expanding viral discovery and the potential for biological inference. Evidence of viruses actively infecting the dominant microbes in subzero brines and modulating host metabolism underscored the potential impact of viruses on these remote and underexplored extreme ecosystems. Microdiversity results shed light on different strategies viruses use to evolve and adapt when extreme conditions are stable versus fluctuating. Together, these findings verify the value of long-read-powered viromics and provide foundational data on viral evolution and virus-microbe interactions in Earth's destabilized and rapidly disappearing cryosphere. Video Abstract.

Correction: Alathari et al. A Multiplexed, Tiled PCR Method for Rapid Whole-Genome Sequencing of Infectious Spleen and Kidney Necrosis Virus (ISKNV) in Tilapia. Viruses 2023, 15, 965

In the original publication [...].

A Reduction of Transcriptional Regulation in Aquatic Oligotrophic Microorganisms Enhances Fitness in Nutrient-Poor Environments

In this review, we consider the regulatory strategies of aquatic oligotrophs, microbial cells that are adapted to thrive under low-nutrient concentrations in oceans, lakes, and other aquatic ecosystems. Many reports have concluded that oligotrophs use less transcriptional regulation than copiotrophic cells, which are adapted to high nutrient concentrations and are far more common subjects for laboratory investigations of regulation. It is theorized that oligotrophs have retained alternate mechanisms of regulation, such as riboswitches, that provide shorter response times and smaller amplitude responses and require fewer cellular resources. We examine the accumulated evidence for distinctive regulatory strategies in oligotrophs. We explore differences in the selective pressures copiotrophs and oligotrophs encounter and ask why, although evolutionary history gives copiotrophs and oligotrophs access to the same regulatory mechanisms, they might exhibit distinctly different patterns in how these mechanisms are used. We discuss the implications of these findings for understanding broad patterns in the evolution of microbial regulatory networks and their relationships to environmental niche and life history strategy. We ask whether these observations, which have emerged from a decade of increased investigation of the cell biology of oligotrophs, might be relevant to recent discoveries of many microbial cell lineages in nature that share with oligotrophs the property of reduced genome size.

A Multiplexed, Tiled PCR Method for Rapid Whole-Genome Sequencing of Infectious Spleen and Kidney Necrosis Virus (ISKNV) in Tilapia

Tilapia farming is one of the most important sectors in aquaculture worldwide and of major importance to global food security. Infectious spleen and kidney necrosis virus (ISKNV) has been identified as an agent of high morbidity and mortality, threatening tilapia aquaculture. ISKNV was detected in Lake Volta, Ghana, in September 2018 and spread rapidly, with mortality rates between 60 and 90% and losses of more than 10 tonnes of fish per day. Understanding the spread and evolution of viral pathogens is important for control strategies. Here, we developed a tiled-PCR sequencing approach for the whole-genome sequencing of ISKNV, using long read sequencing to enable field-based, real-time genomic surveillance. This work represents the first use of tiled-PCR for whole genome recovery of viruses in aquaculture, with the longest genome target (>110 kb dsDNA) to date. Our protocol was applied to field samples collected from the ISKNV outbreaks from four intensive tilapia cage culture systems across Lake Volta, between October 2018 and May 2022. Despite the low mutation rate of dsDNA viruses, 20 single nucleotide polymorphisms accumulated during the sampling period. Droplet digital PCR identified a minimum requirement of template in a sample to recover 50% of an ISKNV genome at 275 femtograms (2410 viral templates per 5 µL sequencing reaction). Overall, tiled-PCR sequencing of ISKNV provides an informative tool to assist in disease control in aquaculture.

Phage Therapy Administration Route, Regimen, and Need for Supplementary Antibiotics in Patients with Chronic Suppurative Lung Disease

Antimicrobial resistance is leading to increased mortality, posing risk to those with chronic suppurative lung disease (CSLD). One therapeutic option may be to target treatment-resistant bacteria using viruses (bacteriophages [phages]). Currently, patients receiving phage therapy on compassionate grounds may not be receiving optimal treatment as there is no defined approach for phage use. This review aims to explore administration route, regimen, and need for supplementary antibiotics in phage therapy to treat bacterial infection in CSLD. Twelve articles totaling 18 participants included details of numerous phage administration routes with varying regimens. All articles reported an initial reduction of bacterial load or an improvement in patient symptoms, highlighting the potential of phage therapy in CSLD. Fifteen out of 18 used supplementary antibiotics. Standardized protocols informed by high-quality research are necessary to ensure safe and effective phage therapy. In the interim, systematic recording of information within case reports may be useful.

Influence of short and long term processes on SAR11 communities in open ocean and coastal systems

SAR11 bacteria dominate the surface ocean and are major players in converting fixed carbon back to atmospheric carbon dioxide. The SAR11 clade is comprised of niche-specialized ecotypes that display distinctive spatiotemporal transitions. We analyzed SAR11 ecotype seasonality in two long-term 16S rRNA amplicon time series representing different North Atlantic regimes: the Sargasso Sea (subtropical ocean-gyre; BATS) and the temperate coastal Western English Channel (WEC). Using phylogenetically resolved amplicon sequence variants (ASVs), we evaluated seasonal environmental constraints on SAR11 ecotype periodicity. Despite large differences in temperature and nutrient availability between the two sites, at both SAR11 succession was defined by summer and winter clusters of ASVs. The summer cluster was dominated by ecotype Ia.3 in both sites. Winter clusters were dominated by ecotypes Ib and IIa.A at BATS and Ia.1 and IIa.B at WEC. A 2-year weekly analysis within the WEC time series showed that the response of SAR11 communities to short-term environmental fluctuations was variable. In 2016, community shifts were abrupt and synchronized to environmental shifts. However, in 2015, changes were gradual and decoupled from environmental fluctuations, likely due to increased mixing from strong winds. We demonstrate that interannual weather variability disturb the pace of SAR11 seasonal progression.

MetaPop: a pipeline for macro- and microdiversity analyses and visualization of microbial and viral metagenome-derived populations

BACKGROUND

Microbes and their viruses are hidden engines driving Earth's ecosystems from the oceans and soils to humans and bioreactors. Though gene marker approaches can now be complemented by genome-resolved studies of inter-(macrodiversity) and intra-(microdiversity) population variation, analytical tools to do so remain scattered or under-developed.

RESULTS

Here, we introduce MetaPop, an open-source bioinformatic pipeline that provides a single interface to analyze and visualize microbial and viral community metagenomes at both the macro- and microdiversity levels. Macrodiversity estimates include population abundances and α- and β-diversity. Microdiversity calculations include identification of single nucleotide polymorphisms, novel codon-constrained linkage of SNPs, nucleotide diversity (π and θ), and selective pressures (pN/pS and Tajima's D) within and fixation indices (FST) between populations. MetaPop will also identify genes with distinct codon usage. Following rigorous validation, we applied MetaPop to the gut viromes of autistic children that underwent fecal microbiota transfers and their neurotypical peers. The macrodiversity results confirmed our prior findings for viral populations (microbial shotgun metagenomes were not available) that diversity did not significantly differ between autistic and neurotypical children. However, by also quantifying microdiversity, MetaPop revealed lower average viral nucleotide diversity (π) in autistic children. Analysis of the percentage of genomes detected under positive selection was also lower among autistic children, suggesting that higher viral π in neurotypical children may be beneficial because it allows populations to better "bet hedge" in changing environments. Further, comparisons of microdiversity pre- and post-FMT in autistic children revealed that the delivery FMT method (oral versus rectal) may influence viral activity and engraftment of microdiverse viral populations, with children who received their FMT rectally having higher microdiversity post-FMT. Overall, these results show that analyses at the macro level alone can miss important biological differences.

CONCLUSIONS

These findings suggest that standardized population and genetic variation analyses will be invaluable for maximizing biological inference, and MetaPop provides a convenient tool package to explore the dual impact of macro- and microdiversity across microbial communities. Video abstract.

Efficient dilution-to-extinction isolation of novel virus-host model systems for fastidious heterotrophic bacteria

Microbes and their associated viruses are key drivers of biogeochemical processes in marine and soil biomes. While viruses of phototrophic cyanobacteria are well-represented in model systems, challenges of isolating marine microbial heterotrophs and their viruses have hampered experimental approaches to quantify the importance of viruses in nutrient recycling. A resurgence in cultivation efforts has improved the availability of fastidious bacteria for hypothesis testing, but this has not been matched by similar efforts to cultivate their associated bacteriophages. Here, we describe a high-throughput method for isolating important virus-host systems for fastidious heterotrophic bacteria that couples advances in culturing of hosts with sequential enrichment and isolation of associated phages. Applied to six monthly samples from the Western English Channel, we first isolated one new member of the globally dominant bacterial SAR11 clade and three new members of the methylotrophic bacterial clade OM43. We used these as bait to isolate 117 new phages, including the first known siphophage-infecting SAR11, and the first isolated phage for OM43. Genomic analyses of 13 novel viruses revealed representatives of three new viral genera, and infection assays showed that the viruses infecting SAR11 have ecotype-specific host ranges. Similar to the abundant human-associated phage ɸCrAss001, infection dynamics within the majority of isolates suggested either prevalent lysogeny or chronic infection, despite a lack of associated genes, or host phenotypic bistability with lysis putatively maintained within a susceptible subpopulation. Broader representation of important virus-host systems in culture collections and genomic databases will improve both our understanding of virus-host interactions, and accuracy of computational approaches to evaluate ecological patterns from metagenomic data.

VirION2: a short- and long-read sequencing and informatics workflow to study the genomic diversity of viruses in nature

Microbes play fundamental roles in shaping natural ecosystem properties and functions, but do so under constraints imposed by their viral predators. However, studying viruses in nature can be challenging due to low biomass and the lack of universal gene markers. Though metagenomic short-read sequencing has greatly improved our virus ecology toolkit—and revealed many critical ecosystem roles for viruses—microdiverse populations and fine-scale genomic traits are missed. Some of these microdiverse populations are abundant and the missed regions may be of interest for identifying selection pressures that underpin evolutionary constraints associated with hosts and environments. Though long-read sequencing promises complete virus genomes on single reads, it currently suffers from high DNA requirements and sequencing errors that limit accurate gene prediction. Here we introduce VirION2, an integrated short- and long-read metagenomic wet-lab and informatics pipeline that updates our previous method (VirION) to further enhance the utility of long-read viral metagenomics. Using a viral mock community, we first optimized laboratory protocols (polymerase choice, DNA shearing size, PCR cycling) to enable 76% longer reads (now median length of 6,965 bp) from 100-fold less input DNA (now 1 nanogram). Using a virome from a natural seawater sample, we compared viromes generated with VirION2 against other library preparation options (unamplified, original VirION, and short-read), and optimized downstream informatics for improved long-read error correction and assembly. VirION2 assemblies combined with short-read based data (‘enhanced’ viromes), provided significant improvements over VirION libraries in the recovery of longer and more complete viral genomes, and our optimized error-correction strategy using long- and short-read data achieved 99.97% accuracy. In the seawater virome, VirION2 assemblies captured 5,161 viral populations (including all of the virus populations observed in the other assemblies), 30% of which were uniquely assembled through inclusion of long-reads, and 22% of the top 10% most abundant virus populations derived from assembly of long-reads. Viral populations unique to VirION2 assemblies had significantly higher microdiversity means, which may explain why short-read virome approaches failed to capture them. These findings suggest the VirION2 sample prep and workflow can help researchers better investigate the virosphere, even from challenging low-biomass samples. Our new protocols are available to the research community on protocols.io as a ‘living document’ to facilitate dissemination of updates to keep pace with the rapid evolution of long-read sequencing technology.

Contrasting impacts of a novel specialist vector on multihost viral pathogen epidemiology in wild and managed bees

Typically, pathogens infect multiple host species. Such multihost pathogens can show considerable variation in their degree of infection and transmission specificity, which has important implications for potential disease emergence. Transmission of multihost pathogens can be driven by key host species and changes in such transmission networks can lead to disease emergence. We study two viruses that show contrasting patterns of prevalence and specificity in managed honeybees and wild bumblebees, black queen cell virus (BQCV) and slow bee paralysis virus (SBPV), in the context of the novel transmission route provided by the virus‐vectoring Varroa destructor. Our key result is that viral communities and RNA virus genetic variation are structured by location, not host species or V. destructor presence. Interspecific transmission is pervasive with the same viral variants circulating between pollinator hosts in each location; yet, we found virus‐specific host differences in prevalence and viral load. Importantly, V. destructor presence increases the prevalence in honeybees and, indirectly, in wild bumblebees, but in contrast to its impact on deformed wing virus (DWV), BQCV and SBPV viral loads are not increased by Varroa presence, and do not show genetic evidence of recent emergence. Effective control of Varroa in managed honeybee colonies is necessary to mitigate further disease emergence, and alleviate disease pressure on our vital wild bee populations. More generally, our results highlight the over‐riding importance of geographical location to the epidemiological outcome despite the complexity of multihost‐parasite interactions.

Knock-on community impacts of a novel vector: spillover of emerging DWV-B from Varroa-infested honeybees to wild bumblebees

Novel transmission routes can directly impact the evolutionary ecology of infectious diseases, with potentially dramatic effect on host populations and knock‐on effects on the wider host community. The invasion of Varroa destructor, an ectoparasitic viral vector in Western honeybees, provides a unique opportunity to examine how a novel vector affects disease epidemiology in a host community. This specialist honeybee mite vectors deformed wing virus (DWV), an important re‐emerging honeybee pathogen that also infects wild bumblebees. Comparing island honeybee and wild bumblebee populations with and without V. destructor, we show that V. destructor drives DWV prevalence and titre in honeybees and sympatric bumblebees. Viral genotypes are shared across hosts, with the potentially more virulent DWV‐B overtaking DWV‐A in prevalence in a current epidemic. This demonstrates disease emergence across a host community driven by the acquisition of a specialist novel transmission route in one host, with dramatic community level knock‐on effects.

Marine DNA Viral Macro- and Microdiversity from Pole to Pole

Microbes drive most ecosystems and are modulated by viruses that impact their lifespan, gene flow, and metabolic outputs. However, ecosystem-level impacts of viral community diversity remain difficult to assess due to classification issues and few reference genomes. Here, we establish an ∼12-fold expanded global ocean DNA virome dataset of 195,728 viral populations, now including the Arctic Ocean, and validate that these populations form discrete genotypic clusters. Meta-community analyses revealed five ecological zones throughout the global ocean, including two distinct Arctic regions. Across the zones, local and global patterns and drivers in viral community diversity were established for both macrodiversity (inter-population diversity) and microdiversity (intra-population genetic variation). These patterns sometimes, but not always, paralleled those from macro-organisms and revealed temperate and tropical surface waters and the Arctic as biodiversity hotspots and mechanistic hypotheses to explain them. Such further understanding of ocean viruses is critical for broader inclusion in ecosystem models.

Long-read viral metagenomics captures abundant and microdiverse viral populations and their niche-defining genomic islands

Marine viruses impact global biogeochemical cycles via their influence on host community structure and function, yet our understanding of viral ecology is constrained by limitations in host culturing and a lack of reference genomes and 'universal' gene markers to facilitate community surveys. Short-read viral metagenomic studies have provided clues to viral function and first estimates of global viral gene abundance and distribution, but their assemblies are confounded by populations with high levels of strain evenness and nucleotide diversity (microdiversity), limiting assembly of some of the most abundant viruses on Earth. Such features also challenge assembly across genomic islands containing niche-defining genes that drive ecological speciation. These populations and features may be successfully captured by single-virus genomics and fosmid-based approaches, at least in abundant taxa, but at considerable cost and technical expertise. Here we established a low-cost, low-input, high throughput alternative sequencing and informatics workflow to improve viral metagenomic assemblies using short-read and long-read technology. The 'VirION' (Viral, long-read metagenomics via MinION sequencing) approach was first validated using mock communities where it was found to be as relatively quantitative as short-read methods and provided significant improvements in recovery of viral genomes. We then then applied VirION to the first metagenome from a natural viral community from the Western English Channel. In comparison to a short-read only approach, VirION: (i) increased number and completeness of assembled viral genomes; (ii) captured abundant, highly microdiverse virus populations, and (iii) captured more and longer genomic islands. Together, these findings suggest that VirION provides a high throughput and cost-effective alternative to fosmid and single-virus genomic approaches to more comprehensively explore viral communities in nature.

Host-hijacking and planktonic piracy: how phages command the microbial high seas

Microbial communities living in the oceans are major drivers of global biogeochemical cycles. With nutrients limited across vast swathes of the ocean, marine microbes eke out a living under constant assault from predatory viruses. Viral concentrations exceed those of their bacterial prey by an order of magnitude in surface water, making these obligate parasites the most abundant biological entities in the ocean. Like the pirates of the 17th and 18th centuries that hounded ships plying major trade and exploration routes, viruses have evolved mechanisms to hijack microbial cells and repurpose their cargo and indeed the vessels themselves to maximise viral propagation. Phenotypic reconfiguration of the host is often achieved through Auxiliary Metabolic Genes - genes originally derived from host genomes but maintained and adapted in viral genomes to redirect energy and substrates towards viral synthesis. In this review, we critically evaluate the literature describing the mechanisms used by bacteriophages to reconfigure host metabolism and to plunder intracellular resources to optimise viral production. We also highlight the mechanisms used when, in challenging environments, a 'batten down the hatches' strategy supersedes that of 'plunder and pillage'. Here, the infecting virus increases host fitness through phenotypic augmentation in order to ride out the metaphorical storm, with a concomitant impact on host substrate uptake and metabolism, and ultimately, their interactions with their wider microbial community. Thus, the traditional view of the virus-host relationship as predator and prey does not fully characterise the variety or significance of the interactions observed. Recent advances in viral metagenomics have provided a tantalising glimpse of novel mechanisms of viral metabolic reprogramming in global oceans. Incorporation of these new findings into global biogeochemical models requires experimental evidence from model systems and major improvements in our ability to accurately predict protein function from sequence data.

The abundant marine bacterium Pelagibacter simultaneously catabolizes dimethylsulfoniopropionate to the gases dimethyl sulfide and methanethiol

Marine phytoplankton produce ∼10(9) tonnes of dimethylsulfoniopropionate (DMSP) per year(1,2), an estimated 10% of which is catabolized by bacteria through the DMSP cleavage pathway to the climatically active gas dimethyl sulfide(3,4). SAR11 Alphaproteobacteria (order Pelagibacterales), the most abundant chemo-organotrophic bacteria in the oceans, have been shown to assimilate DMSP into biomass, thereby supplying this cell's unusual requirement for reduced sulfur(5,6). Here, we report that Pelagibacter HTCC1062 produces the gas methanethiol, and that a second DMSP catabolic pathway, mediated by a cupin-like DMSP lyase, DddK, simultaneously shunts as much as 59% of DMSP uptake to dimethyl sulfide production. We propose a model in which the allocation of DMSP between these pathways is kinetically controlled to release increasing amounts of dimethyl sulfide as the supply of DMSP exceeds cellular sulfur demands for biosynthesis.

Biogeochemical and Microbial Variation across 5500 km of Antarctic Surface Sediment Implicates Organic Matter as a Driver of Benthic Community Structure

Western Antarctica, one of the fastest warming locations on Earth, is a unique environment that is underexplored with regards to biodiversity. Although pelagic microbial communities in the Southern Ocean and coastal Antarctic waters have been well-studied, there are fewer investigations of benthic communities and most have a focused geographic range. We sampled surface sediment from 24 sites across a 5500 km region of Western Antarctica (covering the Ross Sea to the Weddell Sea) to examine relationships between microbial communities and sediment geochemistry. Sequencing of the 16S and 18S rRNA genes showed microbial communities in sediments from the Antarctic Peninsula (AP) and Western Antarctica (WA), including the Ross, Amundsen, and Bellingshausen Seas, could be distinguished by correlations with organic matter concentrations and stable isotope fractionation (total organic carbon; TOC, total nitrogen; TN, and δ¹³C). Overall, samples from the AP were higher in nutrient content (TOC, TN, and NH₄⁺) and communities in these samples had higher relative abundances of operational taxonomic units (OTUs) classified as the diatom, Chaetoceros, a marine cercozoan, and four OTUs classified as Flammeovirgaceae or Flavobacteria. As these OTUs were strongly correlated with TOC, the data suggests the diatoms could be a source of organic matter and the Bacteroidetes and cercozoan are grazers that consume the organic matter. Additionally, samples from WA have lower nutrients and were dominated by Thaumarchaeota, which could be related to their known ability to thrive as lithotrophs. This study documents the largest analysis of benthic microbial communities to date in the Southern Ocean, representing almost half the continental shoreline of Antarctica, and documents trophic interactions and coupling of pelagic and benthic communities. Our results indicate potential modifications in carbon sequestration processes related to change in community composition, identifying a prospective mechanism that links climate change to carbon availability.

Single-cell enabled comparative genomics of a deep ocean SAR11 bathytype

Bacterioplankton of the SAR11 clade are the most abundant microorganisms in marine systems, usually representing 25% or more of the total bacterial cells in seawater worldwide. SAR11 is divided into subclades with distinct spatiotemporal distributions (ecotypes), some of which appear to be specific to deep water. Here we examine the genomic basis for deep ocean distribution of one SAR11 bathytype (depth-specific ecotype), subclade Ic. Four single-cell Ic genomes, with estimated completeness of 55%-86%, were isolated from 770 m at station ALOHA and compared with eight SAR11 surface genomes and metagenomic datasets. Subclade Ic genomes dominated metagenomic fragment recruitment below the euphotic zone. They had similar COG distributions, high local synteny and shared a large number (69%) of orthologous clusters with SAR11 surface genomes, yet were distinct at the 16S rRNA gene and amino-acid level, and formed a separate, monophyletic group in phylogenetic trees. Subclade Ic genomes were enriched in genes associated with membrane/cell wall/envelope biosynthesis and showed evidence of unique phage defenses. The majority of subclade Ic-specfic genes were hypothetical, and some were highly abundant in deep ocean metagenomic data, potentially masking mechanisms for niche differentiation. However, the evidence suggests these organisms have a similar metabolism to their surface counterparts, and that subclade Ic adaptations to the deep ocean do not involve large variations in gene content, but rather more subtle differences previously observed deep ocean genomic data, like preferential amino-acid substitutions, larger coding regions among SAR11 clade orthologs, larger intergenic regions and larger estimated average genome size.

Giovannoni et al. reply

Replying to S. Våge, J. E. Storesund & T. F. Thingstad Nature 499, http://dx.doi.org/10.1038/nature12387 (2013). In the debate about top-down versus bottom-up control, most would argue that both are important. The most probable scenario is that in the absence of high rates of transport, the equilibrium between these competing processes determines bacterial abundance. In their Brief Communication Arising, Våge and co-workers respond to our Letter, arguing the case for defense specialism. ‘Defence specialism’ is not synonymous with ‘defence’. Because of the prevalence of viruses, it is probable that most bacteria have acquired defences against predation. ‘Defence specialism’ refers to success by defence, in the sense that a cell might allocate a large proportion of its resources to defence and thereby become successful.

High-resolution SAR11 ecotype dynamics at the Bermuda Atlantic Time-series Study site by phylogenetic placement of pyrosequences

Advances in next-generation sequencing technologies are providing longer nucleotide sequence reads that contain more information about phylogenetic relationships. We sought to use this information to understand the evolution and ecology of bacterioplankton at our long-term study site in the Western Sargasso Sea. A bioinformatics pipeline called PhyloAssigner was developed to align pyrosequencing reads to a reference multiple sequence alignment of 16S ribosomal RNA (rRNA) genes and assign them phylogenetic positions in a reference tree using a maximum likelihood algorithm. Here, we used this pipeline to investigate the ecologically important SAR11 clade of Alphaproteobacteria. A combined set of 2.7 million pyrosequencing reads from the 16S rRNA V1-V2 regions, representing 9 years at the Bermuda Atlantic Time-series Study (BATS) site, was quality checked and parsed into a comprehensive bacterial tree, yielding 929 036 Alphaproteobacteria reads. Phylogenetic structure within the SAR11 clade was linked to seasonally recurring spatiotemporal patterns. This analysis resolved four new SAR11 ecotypes in addition to five others that had been described previously at BATS. The data support a conclusion reached previously that the SAR11 clade diversified by subdivision of niche space in the ocean water column, but the new data reveal a more complex pattern in which deep branches of the clade diversified repeatedly across depth strata and seasonal regimes. The new data also revealed the presence of an unrecognized clade of Alphaproteobacteria, here named SMA-1 (Sargasso Mesopelagic Alphaproteobacteria, group 1), in the upper mesopelagic zone. The high-resolution phylogenetic analyses performed herein highlight significant, previously unknown, patterns of evolutionary diversification, within perhaps the most widely distributed heterotrophic marine bacterial clade, and strongly links to ecosystem regimes.

Metagenomics: microbial diversity through a scratched lens

Since nucleic acids were first extracted directly from the environment and sequenced, metagenomics has grown to one of the most data-rich and pervasive techniques for understanding the taxonomic and functional diversity of microbial communities. In the last decade, cheaper sequencing has democratized the application of metagenomics and generated billions of reads, revealing staggering microbial diversity and functional complexity. However, cheaper sequencing has come at the cost of reduced sequence length, resulting in poor gene annotation and overestimates of bacterial richness and abundance. Recent improvements in sequencing technology are beginning to provide reads of sufficient length for accurate annotation and assembly of whole operons and beyond, that will once again enable experimental testing of gene function and re-capture the early successes of metagenomic investigations.

Permanent draft genome sequence of Vibrio tubiashii strain NCIMB 1337 (ATCC19106)

Vibrio tubiashii NCIMB 1337 is a major and increasingly prevalent pathogen of bivalve mollusks, and shares a close phylogenetic relationship with both V. orientalis and V. coralliilyticus. It is a Gram-negative, curved rod-shaped bacterium, originally isolated from a moribund juvenile oyster, and is both oxidase and catalase positive. It is capable of growth under both aerobic and anaerobic conditions. Here we describe the features of this organism, together with the draft genome and annotation. The genome is 5,353,266 bp long, consisting of two chromosomes, and contains 4,864 protein-coding and 86 RNA genes.

Novel analysis of oceanic surface water metagenomes suggests importance of polyphosphate metabolism in oligotrophic environments

Polyphosphate is a ubiquitous linear homopolymer of phosphate residues linked by high-energy bonds similar to those found in ATP. It has been associated with many processes including pathogenicity, DNA uptake and multiple stress responses across all domains. Bacteria have also been shown to use polyphosphate as a way to store phosphate when transferred from phosphate-limited to phosphate-rich media--a process exploited in wastewater treatment and other environmental contaminant remediation. Despite this, there has, to date, been little research into the role of polyphosphate in the survival of marine bacterioplankton in oligotrophic environments. The three main proteins involved in polyphosphate metabolism, Ppk1, Ppk2 and Ppx are multi-domain and have differential inter-domain and inter-gene conservation, making unbiased analysis of relative abundance in metagenomic datasets difficult. This paper describes the development of a novel Isofunctional Homolog Annotation Tool (IHAT) to detect homologs of genes with a broad range of conservation without bias of traditional expect-value cutoffs. IHAT analysis of the Global Ocean Sampling (GOS) dataset revealed that genes associated with polyphosphate metabolism are more abundant in environments where available phosphate is limited, suggesting an important role for polyphosphate metabolism in marine oligotrophs.

Metagenomics

Metagenomics has evolved over the last 3 decades from the analysis of single genes and their apparent diversity in an ecosystem to the provision of complex genetic information relating to whole ecosystems. Metagenomics is a vast subject area in terms of methodology, which encompasses a suite of molecular technologies employed to investigate genomic information from all members of a microbial community. However, the relatively recent developments in high-throughput sequencing platforms have meant that metagenomic can be performed simply by extracting DNA and sequencing it. Here, we outline explicit methodologies for the extraction of metagenomic DNA from marine and sediments/soil environmental samples, the subsequent production and sequencing of large-insert metagenomic libraries, and also shotgun pyrosequencing considerations. We also provide relevant advice on bioinformatic analyses of the complex metagenomic datasets. We hope that the information provided here will be useful to establish the techniques in most reasonably equipped molecular biology laboratories.

The taxonomic and functional diversity of microbes at a temperate coastal site: a 'multi-omic' study of seasonal and diel temporal variation

How microbial communities change over time in response to the environment is poorly understood. Previously a six-year time series of 16S rRNA V6 data from the Western English Channel demonstrated robust seasonal structure within the bacterial community, with diversity negatively correlated with day-length. Here we determine whether metagenomes and metatranscriptomes follow similar patterns. We generated 16S rRNA datasets, metagenomes (1.2 GB) and metatranscriptomes (157 MB) for eight additional time points sampled in 2008, representing three seasons (Winter, Spring, Summer) and including day and night samples. This is the first microbial ‘multi-omic’ study to combine 16S rRNA amplicon sequencing with metagenomic and metatranscriptomic profiling. Five main conclusions can be drawn from analysis of these data: 1) Archaea follow the same seasonal patterns as Bacteria, but show lower relative diversity; 2) Higher 16S rRNA diversity also reflects a higher diversity of transcripts; 3) Diversity is highest in winter and at night; 4) Community-level changes in 16S-based diversity and metagenomic profiles are better explained by seasonal patterns (with samples closest in time being most similar), while metatranscriptomic profiles are better explained by diel patterns and shifts in particular categories (i.e., functional groups) of genes; 5) Changes in key genes occur among seasons and between day and night (i.e., photosynthesis); but these samples contain large numbers of orphan genes without known homologues and it is these unknown gene sets that appear to contribute most towards defining the differences observed between times. Despite the huge diversity of these microbial communities, there are clear signs of predictable patterns and detectable stability over time. Renewed and intensified efforts are required to reveal fundamental deterministic patterns in the most complex microbial communities. Further, the presence of a substantial proportion of orphan sequences underscores the need to determine the gene products of sequences with currently unknown function.

Pyrosequencing of Mytilus galloprovincialis cDNAs: tissue-specific expression patterns

Background

Mytilus species are important in marine ecology and in environmental quality assessment, yet their molecular biology is poorly understood. Molecular aspects of their reproduction, hybridisation between species, mitochondrial inheritance, skewed sex ratios of offspring and adaptation to climatic and pollution factors are priority areas.

Methodology/Principal Findings

To start to address this situation, expressed genetic transcripts from M. galloprovincialis were pyrosequenced. Transcripts were isolated from the digestive gland, foot, gill and mantle of both male and female mussels. In total, 175,547 sequences were obtained and for foot and mantle, 90% of the sequences could be assembled into contiguous fragments but this reduced to 75% for the digestive gland and gill. Transcripts relating to protein metabolism and respiration dominated including ribosomal proteins, cytochrome oxidases and NADH dehydrogenase subunits. Tissue specific variation was identified in transcripts associated with mitochondrial energy metabolism, with the digestive gland and gill having the greatest transcript abundance. Using fragment recruitment it was also possible to identify sites of potential small RNAs involved in mitochondrial transcriptional regulation. Sex ratios based on Vitelline Envelop Receptor for Lysin and Vitelline Coat Lysin transcript abundances, indicated that an equal sex distribution was maintained. Taxonomic profiling of the M. galloprovincialis tissues highlighted an abundant microbial flora associated with the digestive gland. Profiling of the tissues for genes involved in intermediary metabolism demonstrated that the gill and digestive gland were more similar to each other than to the other two tissues, and specifically the foot transcriptome was most dissimilar.

Conclusions

Pyrosequencing has provided extensive genomic information for M. galloprovincialis and generated novel observations on expression of different tissues, mitochondria and associated microorganisms. It will also facilitate the much needed production of an oligonucleotide microarray for the organism.

Evidence for phosphonate usage in the coral holobiont

Phosphonates are characterized by a stable carbon-phosphorus bond and commonly occur as lipid conjugates in invertebrate cell membranes. Phosphonoacetate hydrolase encoded by the phnA gene, catalyses the cleavage of phosphonoacetate to acetate and phosphate. In this study, we demonstrate the unusually high phnA diversity in coral-associated bacteria. The holobiont of eight coral species tested positive when screened for phnA using degenerate primers. In two soft coral species, Sinularia and Discosoma, sequencing of the phnA gene showed 13 distinct groups on the basis of 90% sequence identity across 100% of the sequence. A total of 16 bacterial taxa capable of using phosphonoacetate as the sole carbon and phosphorus source were isolated; 8 of which had a phnA+ genotype. This study enhances our understanding of the wide taxonomic and environmental distribution of phnA, and highlights the importance of phosphonates in marine ecosystems.

Bias in assessments of marine microbial biodiversity in fosmid libraries as evaluated by pyrosequencing

On the basis of 16S rRNA gene sequencing, the SAR11 clade of marine bacteria has an almost universal distribution, being detected as abundant sequences in all marine provinces. Yet, SAR11 sequences are rarely detected in fosmid libraries, suggesting that the widespread abundance may be an artefact of PCR cloning and that SAR11 has a relatively low abundance. Here the relative abundance of SAR11 is explored in both a fosmid library and a metagenomic sequence data set from the same biological community taken from fjord surface water from Bergen, Norway. Pyrosequenced data and 16S clone data confirmed an 11-15% relative abundance of SAR11 within the community. In contrast, not a single SAR11 fosmid was identified in a pooled shotgun sequence data set of 100 fosmid clones. This underrepresentation was evidenced by comparative abundances of SAR11 sequences assessed by taxonomic annotation and fragment recruitment. Analysis revealed a similar underrepresentation of low-GC Flavobacteriaceae. We speculate that a contributing factor towards the fosmid bias may be DNA fragmentation during preparation because of the low GC content of SAR11 sequences and other underrepresented taxa. This study suggests that, although fosmid libraries can be extremely useful, caution must be taken when directly inferring community composition from metagenomic fosmid libraries.