Potential virus-mediated nitrogen cycling in oxygen-depleted oceanic waters

Viruses play an important role in the ecology and biogeochemistry of marine ecosystems. Beyond mortality and gene transfer, viruses can reprogram microbial metabolism during infection by expressing auxiliary metabolic genes (AMGs) involved in photosynthesis, central carbon metabolism, and nutrient cycling. While previous studies have focused on AMG diversity in the sunlit and dark ocean, less is known about the role of viruses in shaping metabolic networks along redox gradients associated with marine oxygen minimum zones (OMZs). Here, we analyzed relatively quantitative viral metagenomic datasets that profiled the oxygen gradient across Eastern Tropical South Pacific (ETSP) OMZ waters, assessing whether OMZ viruses might impact nitrogen (N) cycling via AMGs. Identified viral genomes encoded six N-cycle AMGs associated with denitrification, nitrification, assimilatory nitrate reduction, and nitrite transport. The majority of these AMGs (80%) were identified in T4-like Myoviridae phages, predicted to infect Cyanobacteria and Proteobacteria, or in unclassified archaeal viruses predicted to infect Thaumarchaeota. Four AMGs were exclusive to anoxic waters and had distributions that paralleled homologous microbial genes. Together, these findings suggest viruses modulate N-cycling processes within the ETSP OMZ and may contribute to nitrogen loss throughout the global oceans thus providing a baseline for their inclusion in the ecosystem and geochemical models.Subject terms: Virus-host interactions, Biogeochemistry, Microbial biooceanography, Microbial ecology     

Contrasting genomic patterns and infection strategies of two co‐existing Bacteroidetes podovirus genera

Bacterial viruses (phages) are abundant, ecologically important biological entities. However, our understanding of their impact is limited by model systems that are primarily not well represented in nature, e.g. Enterophages and their hosts. Here, we investigate genomic characteristics and infection strategies among six aquatic Bacteroidetes phages that represent two genera of exceptionally large (～70–75 kb genome) podoviruses, which were isolated from the same seawater sample using Cellulophaga baltica as host. Quantitative host range studies reveal that these genera have contrasting narrow (specialist) and broad (generalist) host ranges, with one‐step growth curves revealing reduced burst sizes for the generalist phages. Genomic comparisons suggest candidate genes in each genus that might explain this host range variation, as well as provide hypotheses about receptors in the hosts. One generalist phage, φ38:1, was more deeply characterized, as its infection strategy switched from lytic on its original host to either inefficient lytic or lysogenic on an alternative host. If lysogenic, this phage was maintained extrachromosomally in the alternative host and could not be induced by mitomycin C. This work provides fundamental knowledge regarding phage‐host ranges and their genomic drivers while also exploring the ‘host environment’ as a driver for switching phage replication mode.     

Diversity,evolution and life strategies of CbK-like phages

Caulobacter phage CbK has been extensively studied as a model system in virology and bacteriology. Lysogeny-related genes have been found in each CbK-like isolate, suggesting a life strategy of both lytic and lysogenic cycles. However, whether CbK-related phages can enter lysogeny is still undetermined. This study identified new CbK-like sequences and expanded the collection of CbK-related phages. A common ancestry with a temperate lifestyle was predicted for the group, however, which subsequently evolved into two clades of different genome sizes and host associations. Through the examination of phage recombinase genes, alignment of attachment sites on the phage and bacterial genomes (attP–attB pairing), and the experimental validation, different lifestyles were found among the different members. A majority of clade II members retain a lysogenic lifestyle, whereas all clade I members have evolved into an obligate lytic lifestyle via a loss of the gene encoding Cre-like recombinase and the coupled attP fragment. We postulated that the loss of lysogeny may be a by-product of the increase in phage genome size, and vice versa. Clade I is likely to overcome the costs through maintaining more auxiliary metabolic genes (AMGs), particularly for those involved in protein metabolism, to strengthen host takeover and further benefit virion production.     

Comparative whole genome analysis of six diagnostic brucellaphages

Whole genome sequencing of six diagnostic brucellaphages, Tbilisi (Tb), Firenze (Fz), Weybridge (Wb), S708, Berkeley (Bk) and R/C, was followed with genomic comparisons including recently described genomes of the Tb phage from Mexico (Tb_M) and Pr phage to elucidate genomic diversity and candidate host range determinants. Comparative whole genome analysis revealed high sequence homogeneity among these brucellaphage genomes and resolved three genetic groups consistent with defined host range phenotypes. Group I was composed of Tb and Fz phages that are predominantly lytic for Brucella abortus and Brucella neotomae; Group II included Bk, R/C, and Pr phages that are lytic mainly for B. abortus, Brucella melitensis and Brucella suis; Group III was composed of Wb and S708 phages that are lytic for B. suis, B. abortus and B. neotomae. We found that the putative phage collar protein is a variable locus with features that may be contributing to the host specificities exhibited by different brucellaphage groups. The presence of several candidate host range determinants is illustrated herein for future dissection of the differential host specificity observed among these phages.     

Mobilization of Genomic Islands of Staphylococcus aureus by Temperate Bacteriophage

The virulence of Staphylococcus aureus, in both human and animal hosts, is largely influenced by the acquisition of mobile genetic elements (MGEs). Most S. aureus strains carry a variety of MGEs, including three genomic islands (νSaα, νSaβ, νSaγ) that are diverse in virulence gene content but conserved within strain lineages. Although the mobilization of pathogenicity islands, phages and plasmids has been well studied, the mobilization of genomic islands is poorly understood. We previously demonstrated the mobilization of νSaβ by the adjacent temperate bacteriophage ϕSaBov from strain RF122. In this study, we demonstrate that ϕSaBov mediates the mobilization of νSaα and νSaγ, which are located remotely from ϕSaBov, mostly to recipient strains belonging to ST151. Phage DNA sequence analysis revealed that chromosomal DNA excision events from RF122 were highly specific to MGEs, suggesting sequence-specific DNA excision and packaging events rather than generalized transduction by a temperate phage. Disruption of the int gene in ϕSaBov did not affect phage DNA excision, packaging, and integration events. However, disruption of the terL gene completely abolished phage DNA packing events, suggesting that the primary function of temperate phage in the transfer of genomic islands is to allow for phage DNA packaging by TerL and that transducing phage particles are the actual vehicle for transfer. These results extend our understanding of the important role of bacteriophage in the horizontal transfer and evolution of genomic islands in S. aureus.     

一株铅黄肠球菌噬菌体的生物学特性和全基因组分析

【目的】铅黄肠球菌是医源感染的机会致病菌,可引起危及生命的败血症、脑膜炎等,但针对其噬菌体的研究尚属空白。噬菌体作为细菌病毒,具有宿主特异性。本研究首次分离到可培养的铅黄肠球菌烈性噬菌体,对其基因组序列的分析和其他特征研究为进一步探讨噬菌体与宿主的作用机制及治疗应用提供参考。【方法】噬菌体Ecf_virus_SZ01以健康人粪便中分离的铅黄肠球菌(DO55)作为宿主菌,分离自深圳市南山区未经处理的生活污水样本,利用透射电镜观察噬菌体形态并对其生物学特征和基因组特点进行研究。【结果】透射电镜显示,噬菌体Ecf_virus_SZ01头部直径约为106 nm,尾部直径约为150 nm,尾长且无伸缩性尾鞘,属长尾噬菌体科;该噬菌体的最佳感染复数为0.01;一步生长曲线显示,潜伏期约为30 min,每个受感染细胞产生子代的平均数量为50 PFU/cell;抑菌曲线显示MOI=0.01时对宿主菌具有很好的抑制效果;宿主特异性强,不能实现跨属侵染;测序结果显示其基因组为dsDNA,长度为59 409 bp,GC含量为43.2%;该噬菌体共有102个开放阅读框,BLASTn比对显示该噬菌体与NCBI数据库中其他噬菌体相似性极低。【结论】首次分离到宿主为铅黄肠球菌的噬菌体,具有潜伏期短、裂解能力强、宿主专一的特征,基因组与数据库中现有噬菌体相比十分新颖,并对其生物学特性和基因组进行了分析。     

Continental-scale distribution and diversity of Ceratobasidium orchid mycorrhizal fungi in Australia

Background and AimsMycorrhizal fungi are a critical component of the ecological niche of most plants and can potentially constrain their geographical range. Unlike other types of mycorrhizal fungi, the distributions of orchid mycorrhizal fungi (OMF) at large spatial scales are not well understood. Here, we investigate the distribution and diversity of Ceratobasidium OMF in orchids and soils across the Australian continent.MethodsWe sampled 217 Ceratobasidium isolates from 111 orchid species across southern Australia and combined these with 311 Ceratobasidium sequences from GenBank. To estimate the taxonomic diversity of Ceratobasidium associating with orchids, phylogenetic analysis of the ITS sequence locus was undertaken. Sequence data from the continent-wide Australian Microbiome Initiative were used to determine the geographical range of operational taxonomic units (OTUs) detected in orchids, with the distribution and climatic correlates of the two most frequently detected OTUs modelled using MaxEnt.Key ResultsWe identified 23 Ceratobasidium OTUs associating with Australian orchids, primarily from the orchid genera Pterostylis, Prasophyllum, Rhizanthella and Sarcochilus. OTUs isolated from orchids were closely related to, but distinct from, known pathogenic fungi. Data from soils and orchids revealed that ten of these OTUs occur on both east and west sides of the continent, while 13 OTUs were recorded at three locations or fewer. MaxEnt models suggested that the distributions of two widespread OTUs are correlated with temperature and soil moisture of the wettest quarter and far exceeded the distributions of their host orchid species.ConclusionsCeratobasidium OMF with cross-continental distributions are common in Australian soils and frequently have geographical ranges that exceed that of their host orchid species, suggesting these fungi are not limiting the distributions of their host orchids at large spatial scales. Most OTUs were distributed within southern Australia, although several OTUs had distributions extending into central and northern parts of the continent, illustrating their tolerance of an extraordinarily wide range of environmental conditions.     

Characterization of a temperate actinophage,MPphiWR-1, capable of infectingMicromonospora purpurea ATCC 15835

Summary A temperate actinophage was isolated from soil using the gentamicin-producing microorganism,Micromonospora purpurea ATCC 15835 as host. The characterization of the phage represents the initial step in its development as a cloning vector. The phage isolated, MPphiWR-1, formed red- to purple-pigmented turbid plaques. Cells isolated from these plaques were resistant to superinfection with lytic mutants of MPphiWR-1. Southern blots of genomic DNA from a resistant culture showed that MPphiWR-1 integrated into the host genome. The phage was UV- or Mitomycin C-inducible. The integration, resistance to superinfection and inducibility indicated a lysogenic relationship with the host.Using MPphiE-RCPM, a lytic derivative, the phage host range was demonstrated to include members of three genera: one species each ofAmpullariella andCatellatospora, and 12 species ofMicromonospora. The phage belonged to Ackerman's B1 morphotype having an isometric head and a flexible noncontractile tail. The density of the phage was 1.525 g/cc. Restriction site mapping demonstrated that the phage DNA was 57.9 kb long and had cohesive ends. Using EDTA enrichment, viable mutants with deletions of at least 3.5 kb were isolated and mapped. Phage adsorption, sensitivities and plating efficiency were investigated. Non-liposome PEG-mediated transfection was demonstrated.     

Geographical range and host breadth of Sebacina orchid mycorrhizal fungi associating with Caladenia in south‐western Australia

Specialized mycorrhizal interactions have the potential to limit the geographical range of plant species and contribute to reproductive isolation. We investigated these predictions in Caladenia (Orchidaceae) from south‐western Australia, a group known to have specialized mycorrhizal associations with the genus Sebacina s.l. Sequencing of fungal isolates from 47 of the 136 species of Western Australian Caladenia was undertaken to resolve the geographical range and habitat preferences of mycorrhizal fungal operational taxonomic units (OTUs) and their host breadth in Caladenia. Eight different fungal OTUs were used by Caladenia, with the more frequently detected OTUs occurring in a wide range of habitats and geographical regions. Given the comparatively narrow geographical ranges of most Western Australian Caladenia taxa, this suggests that the geographical ranges of fungal OTUs are unlikely to limit the geographical range of Caladenia spp. Extensive sharing of fungal OTUs between closely related orchid species was detected, suggesting that in the main there is little potential for mycorrhizal fungi to contribute to reproductive isolation between Caladenia spp. Our data mostly support previous work suggesting high mycorrhizal specificity in Caladenia, but this may not be the case in all subgenera, highlighting that Caladenia may offer powerful opportunities for investigating the evolution of specialized mycorrhizal associations.     

Genome-informed approach to identify genetic determinants of Flavobacterium psychrophilum phage susceptibility

The fish pathogen Flavobacterium psychrophilum infects farmed salmonids worldwide, and application of bacteriophages has been suggested for controlling disease outbreaks in aquaculture. Successful application of phages requires detailed knowledge about the variability in phage susceptibility of the host communities. In this study, we analysed the genetic diversity of F. psychrophilum hosts and phages from the Baltic Sea area to identify genetic determinants of phage-host interaction patterns. A host range analysis of 103 phages tested against 177 F. psychrophilum strains (18 231 phage–host interactions) identified nine phage clusters, infecting from 10% to 91% of the strain collection. The core genome-based comparison of 35 F. psychrophilum isolates revealed an extremely low overall genomic diversity (>99.5% similarity). However, a small subset of 16 ORFs, including genes involved in the type IX secretion system (T9SS), gliding motility and hypothetical cell-surface related proteins, exhibited a highly elevated genetic diversity. These specific genetic variations were linked to variability in phage infection patterns obtained from experimental studies, indicating that these genes are key determinants of phage susceptibility. These findings provide novel insights on the molecular mechanisms determining phage susceptibility in F. psychrophilum and emphasizes the importance of phages as drivers of core genomic diversity in this pathogen.     

The genera of bacteriophages and their receptors are the major determinants of host range

The host range of phages is a key to understand their impact on bacterial ecology and evolution. Because of the complexity of phage–host interactions, the variables that determine the breadth of a phage host range remain poorly understood. Here, we propose a novel holistic approach to identify the host range determinants of a new collection of phages infecting Salmonella, isolated from animal, environmental and wastewater samples that were able to infect 58 of the 71 Salmonella strains in our collection. By using a set of statistic approaches (non-metric dimensional scaling, Bray–Curtis distance, PERMANOVA), we analysed phenotypic (host range on wild-type and receptor mutants) and genetic data (taxonomic assignment and receptor binding proteins) to evaluate the impact of isolation strain and niche, phage receptor and genus on the host range. Statistical analysis revealed that two phage characteristics influence the host range by explaining the most variance: the receptor by 45% and the genus by 51%. Interestingly, phage genus and receptor in combination explained 79% of the variance, establishing these characteristics as the major determinants of the host range. This study demonstrates the power and the novelty of applying statistical approaches to phenotypic and genetic data to investigate the ecology of phage–host interactions.     

Isolation of lytic bacteriophage against Vibrio harveyi

Aims: The isolation of lytic bacteriophage of Vibrio harveyi with potential for phage therapy of bacterial pathogens of phyllosoma larvae from the tropical rock lobster Panulirus ornatus. Methods and Results: Water samples from discharge channels and grow‐out ponds of a prawn farm in northeastern Australia were enriched for 24 h in a broth containing four V. harveyi strains. The bacteriophage‐enriched filtrates were spotted onto bacterial lawns demonstrating that the bacteriophage host range for the samples included strains of V. harveyi, Vibrio campbellii, Vibrio rotiferianus, Vibrio parahaemolyticus and Vibrio proteolyticus. Bacteriophage were isolated from eight enriched samples through triple plaque purification. The host range of purified phage included V. harveyi, V. campbellii, V. rotiferianus and V. parahaemolyticus. Transmission electron microscope examination revealed that six purified phage belonged to the family Siphoviridae, whilst two belonged to the family Myoviridae. The Myoviridae appeared to induce bacteriocin production in a limited number of host bacterial strains, suggesting that they were lysogenic rather than lytic. A purified Siphoviridae phage could delay the entry of a broth culture of V. harveyi strain 12 into exponential growth, but could not prevent the overall growth of the bacterial strain. Conclusions: Bacteriophage with lytic activity against V. harveyi were isolated from prawn farm samples. Purified phage of the family Siphoviridae had a clear lytic ability and no apparent transducing properties, indicating they are appropriate for phage therapy. Phage resistance is potentially a major constraint to the use of phage therapy in aquaculture as bacteria are not completely eliminated. Significance and Impact of the Study: Phage therapy is emerging as a potential antibacterial agent that can be used to control pathogenic bacteria in aquaculture systems. The development of phage therapy for aquaculture requires initial isolation and determination of the bacteriophage host range, with subsequent creation of suitable phage cocktails.     

基于高通量测序分析盐角草根部内生细菌多样性及动态规律

【目的】探讨盐角草根部内生细菌群落多样性特征,揭示内生细菌群落结构在宿主关键发育期动态变化规律。【方法】通过罗氏454高通量测序获得内生细菌16S r RNA片段,然后进行生物信息分析。【结果】共获得20363条16S r RNA基因序列。各样品中可操作分类单元(operational taxonomic units,OTUs)在552–941之间。根部内生细菌群落主要包括4个门,其中Proteobacteri门占主导地位,其余依次是Firmicutes,Actinobacteria,Bacteroidetes。在Proteobacteria门中,Gammaproteobacteria是第一大纲,其后是Betaproteobacteria纲。宿主5个发育时期共同拥有7个细菌属,包括Azomonas,Serratia,Pantoea,Serpens,Pseudomonas,Halomonas,Kushneria。整体上看,Gammaproteobacteria纲在宿主5个发育时期呈现增长趋势。优势菌属在5个发育期存在差异,分别为Delftia,Kushneria,Serratia,Pantoea,Erwinia。所有文库总共含2108个特异OTUs,共同拥有5个相同OTUs。花期OTUs数量最多,结种期内生细菌多样性降低。在宿主的5个发育时期中,土壤p H、月均温和土壤盐含量这3个环境因子组成的集合对其内生细菌群落变化具有显著影响。【结论】盐角草内生细菌群落多样性丰富,宿主发育期决定了内生细菌群落结构。     

Identification of Specialists and Abundance-Occupancy Relationships among Intestinal Bacteria of Aves,Mammalia, and Actinopterygii

The coalescence of next-generation DNA sequencing methods, ecological perspectives, and bioinformatics analysis tools is rapidly advancing our understanding of the evolution and function of vertebrate-associated bacterial communities. Delineation of host-microbe associations has applied benefits ranging from clinical treatments to protecting our natural waters. Microbial communities follow some broad-scale patterns observed for macroorganisms, but it remains unclear how the specialization of intestinal vertebrate-associated communities to a particular host environment influences broad-scale patterns in microbial abundance and distribution. We analyzed the V6 region of 16S rRNA genes amplified from 106 fecal samples spanning Aves, Mammalia, and Actinopterygii (ray-finned fish). We investigated the interspecific abundance-occupancy relationship, where widespread taxa tend to be more abundant than narrowly distributed taxa, among operational taxonomic units (OTUs) within and among host species. In a separate analysis, we identified specialist OTUs that were highly abundant in a single host and rare in all other hosts by using a multinomial model without excluding undersampled OTUs a priori. We show that intestinal microbes in humans and other vertebrates display abundance-occupancy relationships, but because intestinal host-associated communities have undergone intense specialization, this trend is violated by a disproportionately large number of specialist taxa. Although it is difficult to distinguish the effects of dispersal limitations, host selection, historical contingency, and stochastic processes on community assembly, results suggest that intestinal bacteria can be shared among diverse hosts in ways that resemble the distribution of “free-living” bacteria in the extraintestinal environment.     

Extraordinary diversity of viruses in deep-sea sediments as revealed by metagenomics without prior virion separation

Our current knowledge of the virosphere in deep-sea sediments remains rudimentary. Here we investigated viral diversity at both gene and genomic levels in deep-sea sediments of Southwest Indian Ocean. Analysis of 19 676 106 non-redundant genes from the metagenomic DNA sequences revealed a large number of unclassified viral groups in these samples. A total of 1106 high-confidence viral contigs were obtained after two runs of assemblies, and 217 of these contigs with sizes up to ~120 kb were shown to represent complete viral genomes. These contigs are clustered with no known viral genomes, and over 2/3 of the ORFs on the viral contigs encode no known functions. Furthermore, most of the complete viral contigs show limited similarity to known viral genomes in genome organization. Most of the classified viral contigs are derived from dsDNA viruses belonging to the order Caudovirales, including primarily members of the families Myoviridae, Podoviridae and Siphoviridae. Most of these viruses infect Proteobacteria and, less frequently, Planctomycetes, Firmicutes, Chloroflexi, etc. Auxiliary metabolic genes (AMGs), present in abundance on the viral contigs, appear to function in modulating the host ability to sense environmental gradients and community changes, and to uptake and metabolize nutrients.     

Metamicrobiomics in herbivore beetles of the genus Cryptocephalus (Chrysomelidae): toward the understanding of ecological determinants in insect symbiosis

The Cryptocephalus marginellus (Coleoptera: Chrysomelidae) complex is composed by six species that are supposed to have originated by events of allo‐ or parapatric speciation. In the present study we investigated the alternative hypotheses that the bacterial communities associated with six populations of this species complex are shaped by environmental factors, or reflect the proposed pattern of speciation. The microbiota associated with the six populations, from five species of the complex, have been characterized through 16S rRNA pyrotag sequencing. Based on a 97% sequence similarity threshold, data were clustered into 381 OTUs, which were analyzed using a variety of diversity indices. The microbiota of C. acquitanus and C. marginellus (Calanques) were the most diverse (over 100 OTUs), while that from C. zoiai yielded less bacterial diversity (45 OTUs). Taxonomic assignment revealed Proteobacteria, Tenericutes and Firmicutes as the dominant components of these beetles’ microbiota. The most abundant genera were Ralstonia, Sphingomonas, Rickettsia, and Pseudomonas. Different strains of Rickettsia were detected in C. eridani and C. renatae. The analysis of β‐diversity revealed high OTU turnover among the populations of C. marginellus complex, with only few shared species. Hierarchical clustering taking into account relative abundances of OTUs does not match the phylogeny of the beetles, therefore we hypothesize that factors other than phylogenetic constraints play a role in shaping the insects’ microbiota. Environmental factors that could potentially affect the composition of bacterial communities were tested by fitting them on the results of a multi‐dimensional scaling analysis. No significant correlations were observed towards the geographic distances or the host plants, while the composition of the microbiota appeared associated with altitude. The metabolic profiles of the microbiotas associated with each population were inferred from bacterial taxonomy, and interestingly, the obtained clustering pattern was consistent with the host phylogeny.     

Disentangling the assembly mechanisms of ant cuticular bacterial communities of two Amazonian ant species sharing a common arboreal nest

Bacteria living on the cuticle of ants are generally studied for their protective role against pathogens, especially in the clade of fungus‐growing ants. However, little is known regarding the diversity of cuticular bacteria in other ant host species, as well as the mechanisms leading to the composition of these communities. Here, we used 16S rRNA gene amplicon sequencing to study the influence of host species, species interactions and the pool of bacteria from the environment on the assembly of cuticular bacterial communities on two phylogenetically distant Amazonian ant species that frequently nest together inside the roots system of epiphytic plants, Camponotus femoratus and Crematogaster levior. Our results show that (a) the vast majority of the bacterial community on the cuticle is shared with the nest, suggesting that most bacteria on the cuticle are acquired through environmental acquisition, (b) 5.2% and 2.0% of operational taxonomic units (OTUs) are respectively specific to Ca. femoratus and Cr. levior, probably representing their respective core cuticular bacterial community, and (c) 3.6% of OTUs are shared between the two ant species. Additionally, mass spectrometry metabolomics analysis of metabolites on the cuticle of ants, which excludes the detection of cuticular hydrocarbons produced by the host, were conducted to evaluate correlations among bacterial OTUs and m/z ion mass. Although some positive and negative correlations are found, the cuticular chemical composition was weakly species‐specific, suggesting that cuticular bacterial communities are prominently environmentally acquired. Overall, our results suggest the environment is the dominant source of bacteria found on the cuticle of ants.     

Organochlorine contamination enriches virus-encoded metabolism and pesticide degradation associated auxiliary genes in soil microbiomes

Viruses significantly influence local and global biogeochemical cycles and help bacteria to survive in different environments by encoding various auxiliary metabolic genes (AMGs) associated with energy acquisition, stress tolerance and degradation of xenobiotics. Here we studied whether bacterial (dsDNA) virus encoded AMGs are enriched in organochlorine pesticide (OCP) contaminated soil in China and if viral AMGs include genes linked to OCP biodegradation. Using metagenomics, we found that OCP-contaminated soils displayed a lower bacterial, but higher diversity of viruses that harbored a higher relative abundance of AMGs linked to pesticide degradation and metabolism. Furthermore, the diversity and relative abundance of AMGs significantly increased along with the severity of pesticide contamination, and several biodegradation genes were identified bioinformatically in viral metagenomes. Functional assays were conducted to experimentally demonstrate that virus-encoded L-2-haloacid dehalogenase gene (L-DEX) is responsible for the degradation of L-2-haloacid pesticide precursors, improving bacterial growth at sub-inhibitory pesticide concentrations. Taken together, these results demonstrate that virus-encoded AMGs are linked to bacterial metabolism and biodegradation, being more abundant and diverse in soils contaminated with pesticides. Moreover, our findings highlight the importance of virus-encoded accessory genes for bacterial ecology in stressful environments, providing a novel avenue for using viruses in the bioremediation of contaminated soils.Subject terms: Metagenomics, Soil microbiology, Microbial ecology

As the most abundant biological entities on earth, viruses of bacteria (bacteriophages referred as viruses from here on) play a critical role in modulating the ecology of microbial communities through lytic infection and lysogenic conversion of their bacterial hosts [1, 2]. Viruses significantly influence the biogeochemical cycles via the release of organic carbon and nutrients through host cell lysis, and in addition to core viral genes (i.e., genes encoding viral structural proteins [3]), they also encode various auxiliary metabolic genes (AMGs [4, 5]), which contribute the metabolic capacity and survival of their bacterial hosts. The role of AMGs has been especially well demonstrated with marine viruses that encode a diversity of AMGs involved in photosynthesis [6], translation machinery [7], carbon metabolism [8], phosphate metabolism [9] and sulfur cycle [10, 11]. Furthermore, sequencing of whole marine viral communities has revealed a clear involvement of viral AMGs in central carbon metabolism of host bacteria [10, 12]. Compared with the study of viral communities in marine ecosystem, the diversity and functional role of viral AMGs in soils are less well understood.In soils, viruses reach abundances of up to ~10⁹ per gram of soil leading to frequent encounters with their host bacteria [13]. Similar to aquatic environments, viruses can regulate host bacterial densities, leading to indirect changes in the relative abundance of non-target bacterial taxa likely via release of niche space [14, 15]. Moreover, over longer time periods, viruses can coevolve with their host, following fluctuating selection dynamics [16] or patterns of local adaptation [17]. Viruses are also important mediators of horizontal gene transfer, promoting the transfer of antibiotic resistance genes, virulence factors and AMGs [18, 19]. However, these effects are less well understood at viral community level. Recent advances in viral purification have enabled a glimpse into soil viral communities of permafrost peatland [20, 21] and agricultural ecosystems [22, 23] based on metagenomics. These studies have demonstrated that viruses may alter the biogeochemical nutrient cycling [1, 2] and bacterial adaptation and evolution by carrying genes linked to carbon and nitrogen metabolism [20, 21]. Moreover, recent identification of atrazine chlorohydrolase trzN [24] and arsenic methyltransferase arsM [25] genes in soil-associated lysogenic viruses suggest that virus-encoded AMGs could shape bacterial metabolism under pollutant exposure. Therefore, we hypothesize, that contaminated soil microbiomes could contain a relatively higher abundance of viruses carrying AMGs linked to the degradation of pesticides and xenobiotics due to their potential benefit for the host bacteria.Pesticide contamination imposes a serious threat to natural ecosystems and public health globally. China is the leading producer of organochlorine pesticides (OCPs), which are synthetic pesticides with vast applications in chemical and agricultural industries. OCPs are especially notorious due to their high toxicity, slow degradation and bioaccumulation [26]. Following the implementation of the Stockholm Convention, hundreds of pesticide plants in China were closed or re-located, and contaminated soils around the plants left untreated. As microbial communities are often capable of degrading OCPs, there is growing biotechnological interest to identify important genes and microbial taxa behind pesticide biodegradation. Heavy OCP contaminations have previously been shown to adversely impact soil bacterial diversity, composition, and activity [27, 28]. Prolonged exposure to contaminants has resulted in selection for bacteria that have evolved their own degradation enzymes, such as dehalogenases, which protect from the toxic effects of pesticides [29]. Interestingly, if also viruses can carry and encode such genes, pesticide exposure could create a strong positive selection for virus-encoded AMGs associated with pesticide degradation, potentially shifting soil microbiome community composition [30] by favoring bacterial and virus taxa that carry these genes.To address this, we used a combination of metagenomics and direct experimentation to explore how pesticide exposure affects the abundance and type of bacterial and virus-encoded AMGs in the soil of former OCP production factory in Yangtze River Delta (China). We found that contaminated and clean control soils harbored very distinct bacterial and viral communities, and crucially, pesticide exposure was linked to higher diversity and abundance of virus-encoded metabolism and pesticide degradation AMGs. The functional activity of one candidate viral AMG, L-2-haloacid dehalogenase (L-DEX), was experimentally shown to improve bacterial growth at sub-inhibitory concentrations of haloacid, which is an important precursor of herbicides and insecticides. Together, our findings suggest that virus-encoded auxiliary genes could help bacteria to counteract pesticide stress, potentially explaining the benefits of virus carriage in stressful soil microbiomes.     

Silage Collected from Dairy Farms Harbors an Abundance of Listeriaphages with Considerable Host Range and Genome Size Diversity

Since the food-borne pathogen Listeria monocytogenes is common in dairy farm environments, it is likely that phages infecting this bacterium (“listeriaphages”) are abundant on dairy farms. To better understand the ecology and diversity of listeriaphages on dairy farms and to develop a diverse phage collection for further studies, silage samples collected on two dairy farms were screened for L. monocytogenes and listeriaphages. While only 4.5% of silage samples tested positive for L. monocytogenes, 47.8% of samples were positive for listeriaphages, containing up to >1.5 × 10⁴ PFU/g. Host range characterization of the 114 phage isolates obtained, with a reference set of 13 L. monocytogenes strains representing the nine major serotypes and four lineages, revealed considerable host range diversity; phage isolates were classified into nine lysis groups. While one serotype 3c strain was not lysed by any phage isolates, serotype 4 strains were highly susceptible to phages and were lysed by 63.2 to 88.6% of phages tested. Overall, 12.3% of phage isolates showed a narrow host range (lysing 1 to 5 strains), while 28.9% of phages represented broad host range (lysing ≥11 strains). Genome sizes of the phage isolates were estimated to range from approximately 26 to 140 kb. The extensive host range and genomic diversity of phages observed here suggest an important role of phages in the ecology of L. monocytogenes on dairy farms. In addition, the phage collection developed here has the potential to facilitate further development of phage-based biocontrol strategies (e.g., in silage) and other phage-based tools.