MetaGeneAnnotator: detecting species-specific patterns of ribosomal binding site for precise gene prediction in anonymous prokaryotic and phage genomes

Complete Genome Sequence of Burkholderia cenocepacia Phage Paku

Burkholderia cenocepacia is able to cause infections in cystic fibrosis patients. B. cenocepacia phage Paku has a 42,727-bp genome sharing a phiKMV-like genome arrangement. T7-like tail components were identified in parallel with a tyrosine integrase, suggesting that Paku might exhibit a temperate lifestyle, an atypical feature for an Autographiviridae phage.

Complete Genome Sequence of Stenotrophomonas maltophilia Siphophage Sonora

Phage Sonora is a siphophage that was isolated against the opportunistic human pathogen Stenotrophomonas maltophilia. The genome of phage Sonora is 63,825 bp long and is not related to that of any phage at the nucleotide level. Sonora shares 46 of 97 total proteins with the Bordetella phages CN2, MW2, and FP1.

Complete Genome Sequence of Stenotrophomonas maltophilia Podophage Piffle

Stenotrophomonas maltophilia is an emerging multidrug-resistant opportunistic human pathogen causing various nosocomial infections. Here, we characterize the genome of S. maltophilia podophage Piffle. Its 76,332-bp genome is most closely related to the N4-like S. maltophilia podophage Pokken, with over 86% genome-wide nucleotide identity and 84 shared proteins.

Complete Genome Sequence of Stenotrophomonas maltophilia Siphophage Silvanus

Stenotrophomonas maltophilia is an opportunistic Gram-negative bacterium capable of causing respiratory infections. S. maltophilia siphophage Silvanus was isolated, and its 45,678-bp genome is not closely related to known phages based on whole-genome comparative genomics analysis. It is predicted to use cos-type packaging due to the similarity of its large terminase subunit to that of phage HK97.

Complete Genome Sequence of Stenotrophomonas maltophilia Myophage Marzo

Stenotrophomonas maltophilia is a Gram-negative opportunistic bacterium that is increasingly being associated with infections. Here, we report the complete genome of the S. maltophilia myophage Marzo, with a 159,384-bp genome encoding 268 proteins, 23 tRNAs, and 1 transfer-messenger RNA. Marzo is closely related to S. maltophilia phages IME-SM1 and Mendera.

Complete Genome Sequence of Stenotrophomonas maltophilia Siphophage Suso

Phage Suso is a temperate siphophage of Stenotrophomonas maltophilia with a 44,659-bp genome. This phage is closely related to Stenotrophomonas phage SM171, sharing 92% overall nucleotide identity as determined by BLASTn, and it shares 14 similar proteins (BLASTp, E value < 0.001) with some Pseudomonas phages from the genus Beetrevirus.

Complete Genome Sequence of Stenotrophomonas maltophilia Podophage Ptah

Stenotrophomonas maltophilia is an opportunistic pathogen demonstrating increasing drug resistance. Here, the genome of the T7-like S. maltophilia podophage Ptah is described. Its 42,593-bp genome is closely related to previously reported T7-like S. maltophilia podophages, including phage Ponderosa.

Complete Genome Sequence of Stenotrophomonas maltophilia Siphophage Suzuki

Stenotrophomonas maltophilia is a Gram-negative bacterium known to cause respiratory tract infections and other diseases in humans. Here, we describe the isolation and genome annotation of S. maltophilia siphophage Suzuki. Its 56,042-bp genome has 83 predicted protein-coding genes and demonstrates similarity with Xylella phages Sano and Salvo.

Complete Genome Sequence of Stenotrophomonas maltophilia Siphophage Summit

Stenotrophomonas maltophilia is an opportunistic pathogen exhibiting resistance to multiple antimicrobials. This study reports the complete genome of an S. maltophilia siphophage, Summit. Its genome of 95,728 bp has 148 protein-coding genes and 5 tRNAs, including 1 predicted suppressor tRNA. Possible target genes for the suppressor tRNA are not identified.

A Novel Phage Infecting the Marine Photoheterotrophic Bacterium Citromicrobium bathyomarinum

This study isolated and characterized a new phage infecting the marine photoheterotrophic bacterium Citromicrobium bathyomarinum, which fills the gap in research on phages targeting this ecologically important species. The phage vB_CbaS-RXM (RXM) has a dsDNA genome with a length of 104,206 bp and G+C content of 61.64%. The taxonomic analysis found a close evolutionary relationship between RXM, Erythrobacter phage vB_EliS-L02, and Sphingobium phage Lacusarx, and we propose that RXM represents a new species of the Lacusarxvirus genus. A one-step growth curve revealed a burst size of 75 plaque-forming units (PFUs) per cell in a 3-hour infection cycle. The lysis profile of RXM showed an intraspecific lethal rate of 26.3% against 38 citromicrobial strains. RXM contains 15 auxiliary metabolic genes (AMGs) related to diverse cellular processes, such as putative metabolic innovation and hijacking of host nucleotide metabolism to enhance its biosynthetic capacity. An in-depth analysis showed that phage functional genes strongly rely on the host for translation, while the translation of unique phage genes with less host dependency may be complemented by phage tRNA. Overall, our study investigated the infection kinetics, genetic traits, taxonomy, and predicted roles of AMGs and tRNA genes of this new phage, which contributes to a better understanding of phage diversity and phage–bacterium interactions.

Sulfur cycling and host-virus interactions in Aquificales-dominated biofilms from Yellowstone's hottest ecosystems

Modern linkages among magmatic, geochemical, and geobiological processes provide clues about the importance of thermophiles in the origin of biogeochemical cycles. The aim of this study was to identify the primary chemoautotrophs and host-virus interactions involved in microbial colonization and biogeochemical cycling at sublacustrine, vapor-dominated vents that represent the hottest measured ecosystems in Yellowstone National Park (~140 °C). Filamentous microbial communities exposed to extreme thermal and geochemical gradients were sampled using a remotely operated vehicle and subjected to random metagenome sequencing and microscopic analyses. Sulfurihydrogenibium (phylum Aquificae) was the predominant lineage (up to 84% relative abundance) detected at vents that discharged high levels of dissolved H₂, H₂S, and CO₂. Metabolic analyses indicated carbon fixation by Sulfurihydrogenibium spp. was powered by the oxidation of reduced sulfur and H₂, which provides organic carbon for heterotrophic community members. Highly variable Sulfurihydrogenibium genomes suggested the importance of intra-population diversity under extreme environmental and viral pressures. Numerous lytic viruses (primarily unclassified taxa) were associated with diverse archaea and bacteria in the vent community. Five circular dsDNA uncultivated virus genomes (UViGs) of ~40 kbp length were linked to the Sulfurihydrogenibium metagenome-assembled genome (MAG) by CRISPR spacer matches. Four UViGs contained consistent genome architecture and formed a monophyletic cluster with the recently proposed Pyrovirus genus within the Caudovirales. Sulfurihydrogenibium spp. also contained CRISPR arrays linked to plasmid DNA with genes for a novel type IV filament system and a highly expressed β-barrel porin. A diverse suite of transcribed secretion systems was consistent with direct microscopic analyses, which revealed an extensive extracellular matrix likely critical to community structure and function. We hypothesize these attributes are fundamental to the establishment and survival of microbial communities in highly turbulent, extreme-gradient environments.

Complete Genome Sequence of Luteitalea sp. Strain TBR-22

We report a complete genome sequence of a novel bacterial isolate, strain TBR-22, belonging to the class Vicinamibacteria of the phylum Acidobacteria, which was isolated from duckweed fronds. The genome expands our knowledge of the lifestyle of this abundant but rarely characterized phylum.

Assembly and Annotation of Escherichia coli Bacteriophage U115

Here, we present the annotated genome sequence of Escherichia coli bacteriophage U115, a T4-like bacteriophage. Phage U115 has a genome length of 166,986 bp and has 286 predicted genes.

Complete Genome Sequence of a Jumbo Bacteriophage, Escherichia Phage vB_EcoM_EC001

Here, we report the genome sequence of a jumbo Escherichia phage vB_EcoM_EC001, a myovirus isolated from primary sludge using enterohemorrhagic Escherichia coli O157:H7. The genome is 240,200 bp long and has 270 predicted coding sequences, including a tryptophanyl tRNA gene. It belongs to genus Seoulvirus.

Characterization of antimicrobial resistance genes and virulence factor genes in an Arctic permafrost region revealed by metagenomics

Antimicrobial resistance genes (ARGs) and virulence factor genes (VFGs) constitute a serious threat to public health, and climate change has been predicted to affect the increase in bacterial pathogens harboring ARGs and VFGs. However, studies on bacterial pathogens and their ARGs and VFGs in permafrost region have received limited attention. In this study, a metagenomic approach was applied to a comprehensive survey to detect potential ARGs, VFGs, and pathogenic antibiotic resistant bacteria (PARB) carrying both ARGs and VFGs in the active layer and permafrost. Overall, 70 unique ARGs against 18 antimicrobial drug classes and 599 VFGs classified as 38 virulence factors were detected in the Arctic permafrost region. Eight genes with mobile genetic elements (MGEs) carrying ARGs were identified; most MGEs were classified as phages. In the metagenome-assembled genomes, the presence of 15 PARB was confirmed. The soil profile showed that the transcripts per million (TPM) values of ARGs and VFGs in the sub-soil horizon were significantly lower than those in the top soil horizon. Based on the TPM value of each gene, major ARGs, VFGs, and these genes in PARB from the Arctic permafrost region were identified and their distribution was confirmed. The major host bacteria for ARGs and VFGs and PARB were identified. A comparison of the percentage identity distribution of ARGs and VFGs to reference databases indicated that ARGs and VFGs in the Arctic soils differ from previously identified genes. Our results may help understand the characteristics and distribution of ARGs, VFGs, and these genes in PARB in the Arctic permafrost region. This findings suggest that the Arctic permafrost region may serve as potential reservoirs for ARGs, VFGs, and PARB. These genes could pose a new threat to human health if they are released by permafrost thawing owing to global warming and propagate to other regions.

Draft Genome Sequence of Bryobacteraceae Strain F-183

Here, we report a draft genome sequence of a bacterial strain, F-183, isolated from a duckweed frond. Strain F-183 belongs to the family Bryobacteraceae of the phylum Acidobacteria, and its genomic information would contribute to understanding the ecophysiology of this abundant but rarely characterized phylum.

Mariniplasma anaerobium gen. nov., sp. nov., a novel anaerobic marine mollicute, and proposal of three novel genera to reclassify members of Acholeplasma clusters II-IV

A novel strictly anaerobic chemoorganotrophic bacterium, designated Mahy22^T, was isolated from sulfidic bottom water of a shallow brackish meromictic lake in Japan. Cells of the strain were Gram-stain-negative, non-motile and coccoid in shape with diameters of about 600-800 nm. The temperature range for growth was 15-37 °C, with optimum growth at 30-32 °C. The pH range for growth was pH 6.2-8.9, with optimum growth at pH 7.2-7.4. The strain grew with NaCl concentrations of 5% or below (optimum, 2-3%). Growth of the strain was enhanced by the addition of thiosulfate. The major cellular fatty acids were C_16:0 and anteiso-C_15:0. Respiratory quinones were not detected. The complete genome sequence of strain Mahy22^T possessed a 1 885 846 bp circular chromosome and a 12 782 bp circular genetic element. The G+C content of the genome sequence was 30.1 mol%. Phylogenetic analysis based on the 16S rRNA gene revealed that the novel strain belonged to the family Acholeplasmataceae, class Mollicutes. The closest relative of strain Mahy22^T with a validly published name was Acholeplasma palmae J233^T with a 16S rRNA gene sequence similarity of 90.5%. Based on the results of polyphasic analysis, the name Mariniplasma anaerobium gen. nov., sp. nov. is proposed to accommodate strain Mahy22^T, along with reclassification of some Acholeplasma species into Alteracholeplasma gen. nov., Haploplasma gen. nov. and Paracholeplasma gen. nov.

No one tool to rule them all: prokaryotic gene prediction tool annotations are highly dependent on the organism of study

MOTIVATION

The biases in CoDing Sequence (CDS) prediction tools, which have been based on historic genomic annotations from model organisms, impact our understanding of novel genomes and metagenomes. This hinders the discovery of new genomic information as it results in predictions being biased towards existing knowledge. To date, users have lacked a systematic and replicable approach to identify the strengths and weaknesses of any CDS prediction tool and allow them to choose the right tool for their analysis.

RESULTS

We present an evaluation framework (ORForise) based on a comprehensive set of 12 primary and 60 secondary metrics that facilitate the assessment of the performance of CDS prediction tools. This makes it possible to identify which performs better for specific use-cases. We use this to assess 15 ab initio- and model-based tools representing those most widely used (historically and currently) to generate the knowledge in genomic databases. We find that the performance of any tool is dependent on the genome being analysed, and no individual tool ranked as the most accurate across all genomes or metrics analysed. Even the top-ranked tools produced conflicting gene collections, which could not be resolved by aggregation. The ORForise evaluation framework provides users with a replicable, data-led approach to make informed tool choices for novel genome annotations and for refining historical annotations.

AVAILABILITY AND IMPLEMENTATION

Code and datasets for reproduction and customisation are available at https://github.com/NickJD/ORForise.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Complete Genome Sequence of the Virulent Klebsiella pneumoniae Phage Geezett Infecting Multidrug-Resistant Clinical Strains

Geezett was isolated from hospital sewage in Hangzhou, China, and exhibits lytic activity against clinical isolates of the nosocomial pathogen Klebsiella pneumoniae. The bacteriophage is a myovirus and has a double-stranded DNA (dsDNA) genome 50,707 bp long, containing 79 open reading frames (ORFs).

Complete Genome Sequences of Bacteriophages Kaya, Guyu, Kopi, and TehO, Which Target Clinical Strains of Pseudomonas aeruginosa

Pseudomonas aeruginosa is a major public health concern, as drug-resistant strains increase mortality in hospital-acquired infections. We report the isolation and complete genome sequences of four lytic bacteriophages that target clinical multidrug-resistant P. aeruginosa strains.

Phage Genome Annotation: Where to Begin and End

With the renewed interest in phage research, coupled with the rising accessibility to affordable sequencing, ever increasing numbers of phage genomes are being sequenced. Therefore, there is an increased need to assemble and annotate phage genomes. There is a plethora of tools and platforms that allow phage genomes to be assembled and annotated. The choice of tools can often be bewildering for those new to phage genome assembly. Here we provide an overview of the assembly and annotation process from obtaining raw reads to genome submission, with worked examples, providing those new to genome assembly and annotation with a guided pathway to genome submission. We focus on the use of open access tools that can be incorporated into workflows to allow easy repetition of steps, highlighting multiple tools that can be used and common pitfalls that may occur.

A review of computational tools for generating metagenome-assembled genomes from metagenomic sequencing data

Metagenomic sequencing provides a culture-independent avenue to investigate the complex microbial communities by constructing metagenome-assembled genomes (MAGs). A MAG represents a microbial genome by a group of sequences from genome assembly with similar characteristics. It enables us to identify novel species and understand their potential functions in a dynamic ecosystem. Many computational tools have been developed to construct and annotate MAGs from metagenomic sequencing, however, there is a prominent gap to comprehensively introduce their background and practical performance. In this paper, we have thoroughly investigated the computational tools designed for both upstream and downstream analyses, including metagenome assembly, metagenome binning, gene prediction, functional annotation, taxonomic classification, and profiling. We have categorized the commonly used tools into unique groups based on their functional background and introduced the underlying core algorithms and associated information to demonstrate a comparative outlook. Furthermore, we have emphasized the computational requisition and offered guidance to the users to select the most efficient tools. Finally, we have indicated current limitations, potential solutions, and future perspectives for further improving the tools of MAG construction and annotation. We believe that our work provides a consolidated resource for the current stage of MAG studies and shed light on the future development of more effective MAG analysis tools on metagenomic sequencing.

Intercrop mulch affects soil biology and microbial diversity in rainfed transgenic Bt cotton hybrids

Growing live mulch between the wide-row spaced transgenic Bt cotton hybrids is a low-cost option to control weeds compared to the use of plastic mulch. However, nothing is known about their effects on soil biology. Therefore, soil samples were collected from a long-term field study (2014-15 to 2018-19) to investigate the soil biological activities as well as the microbial diversity (soil metagenomic analysis). In general, mulching enhanced soil biological activity and influenced the microbial diversity in Bt-cotton. Mulch of sunnhemp (Crotalaria juncea L.), desmodium (Desmodium triflorum L.), sorghum (Sorghum bicolor L.) and plastic sheet recorded significantly higher soil biological activities such as, basal respiration, microbial biomass carbon, and soil enzymes than the other mulch treatments. Aromatic crops (bitter cumin (Centratherum anthelminticum (L.) Kuntze), carom (Trachyspermum ammi (L.) Sprague ex Turrill), coriander (Coriandrum sativum L.), fennel (Foeniculum vulgare Mill.), and fenugreek (Trigonella foenum-graecum L.) had a significant adverse effect on soil biological activity compared to the farmers' practice (no mulch or intercrop). The rarefaction curves as a measure of alpha diversity indicated higher species richness in plastic and newspaper sheet mulch treatments compared to the intercrop mulch. The soil metagenome data indicated Proteobacteria (28%-36%), Actinobacteria (10%-35%), and Acidobacteria (10%-26%) were highly abundant phyla in the mulch treatments. The phyla, Chloroflexi (4%-5%), Gemmatimonadetes (2%-6%), Planctomycetes (2%-4%), and Bacteroidetes (2%-3%) were recorded at lower frequencies in all mulch treatments. The sunnhemp and newspaper mulch treatments recorded low frequency (0.06%-0.07%) of the fungal phyla, Ascomycota. Compared to the bare soil, mulching positively improves soil biological activity. Furthermore, our study identifies some crops that could be grown as an intercrop with a viewpoint to improve soil biology and provide an alternative to the expensive plastic mulch.

Bakta: rapid and standardized annotation of bacterial genomes via alignment-free sequence identification

Command-line annotation software tools have continuously gained popularity compared to centralized online services due to the worldwide increase of sequenced bacterial genomes. However, results of existing command-line software pipelines heavily depend on taxon-specific databases or sufficiently well annotated reference genomes. Here, we introduce Bakta, a new command-line software tool for the robust, taxon-independent, thorough and, nonetheless, fast annotation of bacterial genomes. Bakta conducts a comprehensive annotation workflow including the detection of small proteins taking into account replicon metadata. The annotation of coding sequences is accelerated via an alignment-free sequence identification approach that in addition facilitates the precise assignment of public database cross-references. Annotation results are exported in GFF3 and International Nucleotide Sequence Database Collaboration (INSDC)-compliant flat files, as well as comprehensive JSON files, facilitating automated downstream analysis. We compared Bakta to other rapid contemporary command-line annotation software tools in both targeted and taxonomically broad benchmarks including isolates and metagenomic-assembled genomes. We demonstrated that Bakta outperforms other tools in terms of functional annotations, the assignment of functional categories and database cross-references, whilst providing comparable wall-clock runtimes. Bakta is implemented in Python 3 and runs on MacOS and Linux systems. It is freely available under a GPLv3 license at https://github.com/oschwengers/bakta. An accompanying web version is available at https://bakta.computational.bio.

Viruses infecting a warm water picoeukaryote shed light on spatial co-occurrence dynamics of marine viruses and their hosts

The marine picoeukaryote Bathycoccus prasinos has been considered a cosmopolitan alga, although recent studies indicate two ecotypes exist, Clade BI (B. prasinos) and Clade BII. Viruses that infect Bathycoccus Clade BI are known (BpVs), but not that infect BII. We isolated three dsDNA prasinoviruses from the Sargasso Sea against Clade BII isolate RCC716. The BII-Vs do not infect BI, and two (BII-V2 and BII-V3) have larger genomes (~210 kb) than BI-Viruses and BII-V1. BII-Vs share ~90% of their proteins, and between 65% to 83% of their proteins with sequenced BpVs. Phylogenomic reconstructions and PolB analyses establish close-relatedness of BII-V2 and BII-V3, yet BII-V2 has 10-fold higher infectivity and induces greater mortality on host isolate RCC716. BII-V1 is more distant, has a shorter latent period, and infects both available BII isolates, RCC716 and RCC715, while BII-V2 and BII-V3 do not exhibit productive infection of the latter in our experiments. Global metagenome analyses show Clade BI and BII algal relative abundances correlate positively with their respective viruses. The distributions delineate BI/BpVs as occupying lower temperature mesotrophic and coastal systems, whereas BII/BII-Vs occupy warmer temperature, higher salinity ecosystems. Accordingly, with molecular diagnostic support, we name Clade BII Bathycoccus calidus sp. nov. and propose that molecular diversity within this new species likely connects to the differentiated host-virus dynamics observed in our time course experiments. Overall, the tightly linked biogeography of Bathycoccus host and virus clades observed herein supports species-level host specificity, with strain-level variations in infection parameters.

Nitrospina-like Bacteria Are Dominant Potential Mercury Methylators in Both the Oyashio and Kuroshio Regions of the Western North Pacific

Highly neurotoxic methylmercury (MeHg) accumulates in marine organisms, thereby negatively affecting human and environmental health. Recent studies have revealed that oceanic prokaryotes harboring the hgcAB gene pair are involved in Hg methylation. Presently, little is known about the distribution and phylogeny of these genes in distinct oceanic regions of the western North Pacific. In this study, we used metagenomics to survey the distribution of hgcAB genes in the seawater columns of the subarctic Oyashio region and the subtropical Kuroshio region. The hgcAB genes were detected in the MeHg-rich offshore mesopelagic layers of both the Oyashio region, which is a highly productive area in the western North Pacific, and the Kuroshio region, which has low productivity. Comparative analysis revealed that hgcAB genes belonging to the Nitrospina-like lineage were dominant in the MeHg-rich mesopelagic layers of both regions. These results indicate that Nitrospina-like bacteria are the dominant Hg methylators in the mesopelagic layers throughout the western North Pacific. IMPORTANCE MeHg is highly neurotoxic and accumulates in marine organisms. Thus, understanding MeHg production in seawater is critical for environmental and human health. Recent studies have shown that microorganisms harboring mercury-methylating genes (hgcA and hgcB) are involved in MeHg production in several marine environments. Knowing the distribution and phylogeny of hgcAB genes in seawater columns can facilitate assessment of microbial MeHg production in the ocean. We report that hgcAB genes affiliated with the microaerophilic Nitrospina lineage were detected in the MeHg-rich mesopelagic layers of two hydrologically distinct oceanic regions of the western North Pacific. This finding facilitates understanding of the microbial Hg methylation and accumulation in seawater columns of the western North Pacific.

Complete Genome Sequence of the Lytic Bacteriophage Phab24, Which Infects Clinical Strains of the Nosocomial Pathogen Acinetobacter baumannii

Phab24 was isolated from river water in Zhejiang Province, China, and exhibits lytic activity against clinical isolates of the nosocomial pathogen Acinetobacter baumannii (X. Wang, B. Loh, Y. Yu, X. Hua, et al., bioRxiv 2021.07.23.453473, 2021, https://doi.org/10.1101/2021.07.23.453473). The bacteriophage belongs to the Myoviridae family and has a double-stranded DNA (dsDNA) genome sequence that is 93,604 bp long, containing 172 open reading frames (ORFs).

Assembly and Annotation of the Complete Genome Sequence of T4-Like Bacteriophage 132

Here, we describe the complete sequence of bacteriophage 132, a T4-like Escherichia coli phage. Phage 132 has a genome of 166,922-bp length, with 286 predicted genes.

Perspective on the clone library method for infectious diseases

Recently, culture-independent molecular methods, such as DNA sequencing techniques targeting the 16S-ribosomal RNA (rRNA) gene and/or other housekeeping genes with Sanger method-based technologies, next generation sequencing (NGS), and metagenomic analysis, have been developed for detecting microorganisms in the human body; these can provide information on microbiomes of samples from individuals with or without infectious diseases. Determining the bacterial species is crucial in identifying causative bacteria of upper and lower respiratory tract infections, especially for Streptococcus species, but NGS analysis is often not precise enough to identify bacteria at the species level. This review briefly introduces previous observations of the microbiome of samples from various respiratory and other infections assessed using the clone library method with Sanger sequencing of the 16S-rRNA gene. On analysis of 16S-rRNA gene-sequence data of bronchoalveolar lavage fluid obtained from pneumonia lesions in patients with bacterial pneumonia and lung abscess, anaerobes are often detected in non-elderly patients with pneumonia, and the detection rate of Staphylococcus aureus in patients with hospital-acquired pneumonia is lower than that previously reported. Analysis of pleural effusion samples from patients with pleurisy indicated a more important role of anaerobes than previous believed. The other topics reviewed include microbiomes of nontuberculous mycobacteriosis and lower respiratory tract infections in children with permanent tracheostomy due to neuromuscular disorders, in nasal discharge, in bacterial vaginosis, in the intracystic fluid of postoperative maxillary cyst, and in bacterial conjunctivitis; urine microbiota in urethritis; fecal microbiota; and newly detected infectious organisms in the human respiratory tract.

Comparative Genomics of Three Novel Jumbo Bacteriophages Infecting Staphylococcus aureus

The majority of previously described Staphylococcus aureus bacteriophages belong to three major groups: P68-like podophages, Twort-like or K-like myophages, and a more diverse group of temperate siphophages. Here we present three novel S. aureus "jumbo" phages: MarsHill, Madawaska, and Machias. These phages were isolated from swine production environments in the United States and represent a novel clade of S. aureus myophage. The average genome size for these phages is ∼269 kb with each genome encoding ∼263 predicted protein-coding genes. Phage genome organization and content is similar to known jumbo phages of Bacillus, including AR9 and vB_BpuM-BpSp. All three phages possess genes encoding complete virion and non-virion RNA polymerases, multiple homing endonucleases, and a retron-like reverse transcriptase. Like AR9, all of these phages are presumed to have uracil-substituted DNA which interferes with DNA sequencing. These phages are also able to transduce host plasmids, which is significant as these phages were found circulating in swine production environments and can also infect human S. aureus isolates. Importance of work: This study describes the comparative genomics of three novel S. aureus jumbo phages: MarsHill, Madawaska, and Machias. These three S. aureus myophages represent an emerging class of S. aureus phage. These genomes contain abundant introns which show a pattern consistent with repeated acquisition rather than vertical inheritance, suggesting intron acquisition and loss is an active process in the evolution of these phages. These phages have presumably hypermodified DNA which inhibits sequencing by several different common platforms. Therefore, these phages also represent potential genomic diversity that has been missed due to the limitations of standard sequencing techniques. In particular, such hypermodified genomes may be missed by metagenomic studies due to their resistance to standard sequencing techniques. Phage MarsHill was found to be able to transduce host DNA at levels comparable to that found for other transducing S. aureus phages, making them a potential vector for horizontal gene transfer in the environment.

The physiological potential of anammox bacteria as revealed by their core genome structure

We present here the second complete genome of anaerobic ammonium oxidation (anammox) bacterium, Candidatus (Ca.) Brocadia pituitae, along with those of a nitrite oxidizer and two incomplete denitrifiers from the anammox bacterial community (ABC) metagenome. Although NO2− reduction to NO is considered to be the first step in anammox, Ca. B. pituitae lacks nitrite reductase genes (nirK and nirS) responsible for this reaction. Comparative genomics of Ca. B. pituitae with Ca. Kuenenia stuttgartiensis and six other anammox bacteria with nearly complete genomes revealed that their core genome structure contains 1,152 syntenic orthologues. But nitrite reductase genes were absent from the core, whereas two other Brocadia species possess nirK and these genes were horizontally acquired from multiple lineages. In contrast, at least five paralogous hydroxylamine oxidoreductase genes containing candidate ones (hao2 and hao3) encoding another nitrite reductase were observed in the core. Indeed, these two genes were also significantly expressed in Ca. B. pituitae as in other anammox bacteria. Because many nirS and nirK genes have been detected in the ABC metagenome, Ca. B. pituitae presumably utilises not only NO supplied by the ABC members but also NO and/or NH₂OH by self-production for anammox metabolism.

A Novel Phage Infecting Alteromonas Represents a Distinct Group of Siphophages Infecting Diverse Aquatic Copiotrophs

Bacteriophages play critical roles in impacting microbial community succession both ecologically and evolutionarily. Although the majority of phage genetic diversity has been increasingly unveiled, phages infecting members of the ecologically important genus Alteromonas remain poorly understood. Here, we present a comprehensive analysis of a newly isolated alterophage, vB_AcoS-R7M (R7M), to characterize its life cycle traits, genomic features, and putative evolutionary origin. R7M harbors abundant genes identified as host-like auxiliary metabolic genes facilitating viral propagation. Genomic analysis suggested that R7M is distinct from currently known alterophages. Interestingly, R7M was found to share a set of similar characteristics with a number of siphophages infecting diverse aquatic opportunistic copiotrophs. We therefore proposed the creation of one new subfamily (Queuovirinae) to group with these evolutionarily related phages. Notably, tail genes were less likely to be shared among them, and baseplate-related genes varied the most. In-depth analyses indicated that R7M has replaced its distal tail with a Rhodobacter capsulatus gene transfer agent (RcGTA)-like baseplate and further acquired a putative receptor interaction site targeting Alteromonas. These findings suggest that horizontal exchanges of viral tail adsorption apparatuses are widespread and vital for phages to hunt new hosts and to adapt to new niches. IMPORTANCE The evolution and ecology of phages infecting members of Alteromonas, a marine opportunistic genus that is widely distributed and of great ecological significance, remain poorly understood. The present study integrates physiological and genomic evidence to characterize the properties and putative phage-host interactions of a newly isolated Alteromonas phage, vB_AcoS-R7M (R7M). A taxonomic study reveals close evolutionary relationships among R7M and a number of siphophages infecting various aquatic copiotrophs. Their similar head morphology and overall genetic framework suggest their putative common ancestry and the grouping of a new viral subfamily. However, their major difference lies in the viral tail adsorption apparatuses and the horizontal exchanges of which possibly account for variations in host specificity. These findings outline an evolutionary scenario for the emergence of diverse viral lineages of a shared genetic pool and give insights into the genetics and ecology of viral host jumps.

Enhanced mutualistic symbiosis between soil phages and bacteria with elevated chromium-induced environmental stress

BACKGROUND

Microbe-virus interactions have broad implications on the composition, function, and evolution of microbiomes. Elucidating the effects of environmental stresses on these interactions is critical to identify the ecological function of viral communities and understand microbiome environmental adaptation. Heavy metal-contaminated soils represent a relevant ecosystem to study the interplay between microbes, viruses, and environmental stressors.

RESULTS

Metagenomic analysis revealed that Cr pollution adversely altered the abundance, diversity, and composition of viral and bacterial communities. Host-phage linkage based on CRISPR indicated that, in soils with high Cr contamination, the abundance of phages associated with heavy metal-tolerant hosts increased, as did the relative abundance of phages with broad host ranges (identified as host-phage linkages across genera), which would facilitate transfection and broader distribution of heavy metal resistance genes in the bacterial community. Examining variations along the pollutant gradient, enhanced mutualistic phage-bacterium interactions were observed in the face of greater environmental stresses. Specifically, the fractions of lysogens in bacterial communities (identified by integrase genes within bacterial genomes and prophage induction assay by mitomycin-C) were positively correlated with Cr contamination levels. Furthermore, viral genomic analysis demonstrated that lysogenic phages under higher Cr-induced stresses carried more auxiliary metabolic genes regulating microbial heavy metal detoxification.

CONCLUSION

With the intensification of Cr-induced environmental stresses, the composition, replication strategy, and ecological function of the phage community all evolve alongside the bacterial community to adapt to extreme habitats. These result in a transformation of the phage-bacterium interaction from parasitism to mutualism in extreme environments and underscore the influential role of phages in bacterial adaptation to pollution-related stress and in related biogeochemical processes. Video Abstract.

Differential Bacteriophage Efficacy in Controlling Salmonella in Cattle Hide and Soil Models

Asymptomatic Salmonella carriage in beef cattle is a food safety concern and the beef feedlot environment and cattle hides are reservoirs of this pathogen. Bacteriophages present an attractive non-antibiotic strategy for control of Salmonella in beef. In this study, four diverse and genetically unrelated Salmonella phages, Sergei, Season12, Sw2, and Munch, were characterized and tested alone and in combination for their ability to control Salmonella in cattle hide and soil systems, which are relevant models for Salmonella control in beef production. Phage Sergei is a member of the genus Sashavirus, phage Season12 was identified as a member of the Chivirus genus, Sw2 was identified as a member of the T5-like Epseptimavirus genus, and Munch was found to be a novel “jumbo” myovirus. Observed pathogen reductions in the model systems ranged from 0.50 to 1.75 log₁₀ CFU/cm² in hides and from 0.53 to 1.38 log₁₀ CFU/g in soil, with phages Sergei and Sw2 producing greater reductions (∼1 log₁₀ CFU/cm² or CFU/g) than Season12 and Munch. These findings are in accordance with previous observations of phage virulence, suggesting the simple ability of a phage to form plaques on a bacterial strain is not a strong indicator of antimicrobial activity, but performance in liquid culture assays provides a better predictor. The antimicrobial efficacies of phage treatments were found to be phage-specific across model systems, implying that a phage capable of achieving bacterial reduction in one model is more likely to perform well in another. Phage combinations did not produce significantly greater efficacy than single phages even after 24 h in the soil model, and phage-insensitive colonies were not isolated from treated samples, suggesting that the emergence of phage resistance was not a major factor limiting efficacy in this system.

GeneMark-HM: improving gene prediction in DNA sequences of human microbiome

Computational reconstruction of nearly complete genomes from metagenomic reads may identify thousands of new uncultured candidate bacterial species. We have shown that reconstructed prokaryotic genomes along with genomes of sequenced microbial isolates can be used to support more accurate gene prediction in novel metagenomic sequences. We have proposed an approach that used three types of gene prediction algorithms and found for all contigs in a metagenome nearly optimal models of protein-coding regions either in libraries of pre-computed models or constructed de novo. The model selection process and gene annotation were done by the new GeneMark-HM pipeline. We have created a database of the species level pan-genomes for the human microbiome. To create a library of models representing each pan-genome we used a self-training algorithm GeneMarkS-2. Genes initially predicted in each contig served as queries for a fast similarity search through the pan-genome database. The best matches led to selection of the model for gene prediction. Contigs not assigned to pan-genomes were analyzed by crude, but still accurate models designed for sequences with particular GC compositions. Tests of GeneMark-HM on simulated metagenomes demonstrated improvement in gene annotation of human metagenomic sequences in comparison with the current state-of-the-art gene prediction tools.

Complete Genome Sequence of Klebsiella pneumoniae Podophage Pone

Klebsiella pneumoniae is a Gram-negative pathogen that has become increasingly antibiotic resistant. Phage therapy is potentially a useful approach to controlling this pathogen. Here, we present the genome sequence of the phiKMV-like K. pneumoniae podophage Pone.

A human respiratory tract-associated bacterium with an extremely small genome

Recent advances in culture-independent microbiological analyses have greatly expanded our understanding of the diversity of unculturable microbes. However, human pathogenic bacteria differing significantly from known taxa have rarely been discovered. Here, we present the complete genome sequence of an uncultured bacterium detected in human respiratory tract named IOLA, which was determined by developing a protocol to selectively amplify extremely AT-rich genomes. The IOLA genome is 303,838 bp in size with a 20.7% GC content, making it the smallest and most AT-rich genome among known human-associated bacterial genomes to our best knowledge and comparable to those of insect endosymbionts. While IOLA belongs to order Rickettsiales (mostly intracellular parasites), the gene content suggests an epicellular parasitic lifestyle. Surveillance of clinical samples provides evidence that IOLA can be predominantly detected in patients with respiratory bacterial infections and can persist for at least 15 months in the respiratory tract, suggesting that IOLA is a human respiratory tract-associated bacterium.

Kazumasa Fukuda et al. complete a new genome sequence for an uncultured bacterium detected in human respiratory tract named IOLA. The IOLA genome is found to be among the smallest and most AT-rich of known human-associated bacterial genomes and surveillance of clinical samples indicates that IOLA is in fact a human respiratory tract-associated bacterium.

Complete Genome Sequence of Burkholderia gladioli Myophage Mana

Burkholderia gladioli is known to cause respiratory tract infections in cystic fibrosis patients. Here, we describe the annotation of the 38,038-bp genome sequence of Mana, a P2-like phage of B. gladioli Understanding the genomic characteristics of phages infecting pathogens like B. gladioli can lead to advancements in phage therapy.

A guide to human microbiome research: study design, sample collection, and bioinformatics analysis

The purpose of this review is to provide medical researchers, especially those without a bioinformatics background, with an easy-to-understand summary of the concepts and technologies used in microbiome research. First, we define primary concepts such as microbiota, microbiome, and metagenome. Then, we discuss study design schemes, the methods of sample size calculation, and the methods for improving the reliability of research. We emphasize the importance of negative and positive controls in this section. Next, we discuss statistical analysis methods used in microbiome research, focusing on problems with multiple comparisons and ways to compare β-diversity between groups. Finally, we provide step-by-step pipelines for bioinformatics analysis. In summary, the meticulous study design is a key step to obtaining meaningful results, and appropriate statistical methods are important for accurate interpretation of microbiome data. The step-by-step pipelines provide researchers with insights into newly developed bioinformatics analysis methods.

Infection with Bacteroides Phage BV01 Alters the Host Transcriptome and Bile Acid Metabolism in a Common Human Gut Microbe

Gut-associated phages are hypothesized to alter the abundance and activity of their bacterial hosts, contributing to human health and disease. Although temperate phages constitute a significant fraction of the gut virome, the effects of lysogenic infection are underexplored. We report that the temperate phage, Bacteroides phage BV01, broadly alters its host's transcriptome, the prominent human gut symbiont Bacteroides vulgatus. This alteration occurs through phage-induced repression of a tryptophan-rich sensory protein (TspO) and represses bile acid deconjugation. Because microbially modified bile acids are important signals for the mammalian host, this is a mechanism by which a phage may influence mammalian phenotypes. Furthermore, BV01 and its relatives in the proposed phage family Salyersviridae are ubiquitous in human gut metagenomes, infecting a broad range of Bacteroides hosts. These results demonstrate the complexity of phage-bacteria-mammal relationships and emphasize a need to better understand the role of temperate phages in the gut microbiome.

Seasonality of archaeal proteorhodopsin and associated Marine Group IIb ecotypes (Ca. Poseidoniales) in the North Western Mediterranean Sea

The Archaea Marine Group II (MGII) is widespread in the world's ocean where it plays an important role in the carbon cycle. Despite recent discoveries on the group's metabolisms, the ecology of this newly proposed order (Candidatus Poseidoniales) remains poorly understood. Here we used a combination of time-series metagenome-assembled genomes (MAGs) and high-frequency 16S rRNA data from the NW Mediterranean Sea to test if the taxonomic diversity within the MGIIb family (Candidatus Thalassarchaeaceae) reflects the presence of different ecotypes. The MAGs' seasonality revealed a MGIIb family composed of different subclades that have distinct lifestyles and physiologies. The vitamin metabolisms were notably different between ecotypes with, in some, a possible link to sunlight's energy. Diverse archaeal proteorhodopsin variants, with unusual signature in key amino acid residues, had distinct seasonal patterns corresponding to changing day length. In addition, we show that in summer, archaea, as opposed to bacteria, disappeared completely from surface waters. Our results shed light on the diversity and the distribution of the euryarchaeotal proteorhodopsin, and highlight that MGIIb is a diverse ecological group. The work shows that time-series based studies of the taxonomy, seasonality, and metabolisms of marine prokaryotes is critical to uncover their diverse role in the ocean.

Complete Genome Sequence of Escherichia coli Bacteriophage U136B

We report the genome sequence of bacteriophage U136B, which is reliant on the lipopolysaccharide and the antibiotic efflux protein TolC for infection of Escherichia coli and is a useful model for studying trade-offs and trade-ups that shape evolution. Phage U136B has a 49,233-bp genome with 87 predicted genes.

Diversity and functional landscapes in the microbiota of animals in the wild

Animals in the wild are able to subsist on pathogen-infected and poisonous food and show immunity to various diseases. These may be due to their microbiota, yet we have a poor understanding of animal microbial diversity and function. We used metagenomics to analyze the gut microbiota of more than 180 species in the wild, covering diverse classes, feeding behaviors, geographies, and traits. Using de novo metagenome assembly, we constructed and functionally annotated a database of more than 5000 genomes, comprising 1209 bacterial species of which 75% are unknown. The microbial composition, diversity, and functional content exhibit associations with animal taxonomy, diet, activity, social structure, and life span. We identify the gut microbiota of wild animals as a largely untapped resource for the discovery of therapeutics and biotechnology applications.

Complete Genome Sequence of Stenotrophomonas maltophilia Siphophage Salva

Stenotrophomonas maltophilia is a Gram-negative pathogen causing severe and often refractory illnesses such as pneumonia and bacteremia. We present the genome of phage Salva, a novel S. maltophilia phage that is not closely related to any phages currently deposited in GenBank. The genome is 60,789 bp, containing 102 putative protein-coding genes.

Stenotrophomonas maltophilia is a Gram-negative pathogen causing severe and often refractory illnesses such as pneumonia and bacteremia. We present the genome of phage Salva, a novel S. maltophilia phage that is not closely related to any phages currently deposited in GenBank. The genome is 60,789 bp, containing 102 putative protein-coding genes.

Genome Sequence and Characterization of Acinetobacter Phage DMU1

Background: There has been a recent resurgence of research on the characterization of Acinetobacter phage for therapeutic use due to the morbidity and mortality associated with treatment failures in cases of multidrug-resistant Acinetobacter baumannii infections. Materials and Methods: A bacteriophage isolated from activated sludge that targets A. baumannii ATCC19606 was characterized by electron microscopy, genome sequencing, comparative genomics, and a host range analysis. Results: The morphology of Acinetobacter phage DMU1 resembles phages in Siphoviridae. Comparative genomic and phylogenetic analyses reveal that DMU1 is a siphophage and is most closely related to Acinetobacter phage SH-Ab 15497. Out of the strains tested, DMU1 was found to only infect A. baumannii strains ATCC19606 and ATCC17978. Conclusion: Acinetobacter phage DMU1 belongs to the Siphoviridae family and is most closely related to Acinetobacter phage SH-Ab 15497. Small-scale host-range analysis of DMU1 indicates a host range that is likely limited to specific A. baumannii strains.

Characterization and complete genome sequence of Privateer, a highly prolate Proteus mirabilis podophage

The Gram-negative bacterium Proteus mirabilis causes a large proportion of catheter-associated urinary tract infections, which are among the world's most common nosocomial infections. Here, we characterize P. mirabilis bacteriophage Privateer, a prolate podophage of the C3 morphotype isolated from Texas wastewater treatment plant activated sludge. Basic characterization assays demonstrated Privateer has a latent period of ~40 min and average burst size around 140. In the 90.7 kb Privateer genome, 43 functions were assigned for the 144 predicted protein-coding genes. Genes encoding DNA replication proteins, DNA modification proteins, four tRNAs, lysis proteins, and structural proteins were identified. Cesium-gradient purified Privateer particles analyzed via LC-MS/MS verified the presence of several predicted structural proteins, including a longer, minor capsid protein apparently produced by translational frameshift. Comparative analysis demonstrated Privateer shares 83% nucleotide similarity with Cronobacter phage vB_CsaP_009, but low nucleotide similarity with other known phages. Predicted structural proteins in Privateer appear to have evolutionary relationships with other prolate podophages, in particular the Kuraviruses.

Emergence of nosocomial associated opportunistic pathogens in the gut microbiome after antibiotic treatment

Introduction

According to the Centers for Disease Control’s 2015 Hospital Acquired Infection Hospital Prevalence Survey, 1 in 31 hospital patients was infected with at least one nosocomial pathogen while being treated for unrelated issues. Many studies associate antibiotic administration with nosocomial infection occurrence. However, to our knowledge, there is little to no direct evidence of antibiotic administration selecting for nosocomial opportunistic pathogens.

Aim

This study aims to confirm gut microbiota shifts in an animal model of antibiotic treatment to determine whether antibiotic use favors pathogenic bacteria.

Methodology

We utilized next-generation sequencing and in-house metagenomic assembly and taxonomic assignment pipelines on the fecal microbiota of a urinary tract infection mouse model with and without antibiotic treatment.

Results

Antibiotic therapy decreased the number of detectable species of bacteria by at least 20-fold. Furthermore, the gut microbiota of antibiotic treated mice had a significant increase of opportunistic pathogens that have been implicated in nosocomial infections, like Acinetobacter calcoaceticus/baumannii complex, Chlamydia abortus, Bacteroides fragilis, and Bacteroides thetaiotaomicron. Moreover, antibiotic treatment selected for antibiotic resistant gene enriched subpopulations for many of these opportunistic pathogens.

Conclusions

Oral antibiotic therapy may select for common opportunistic pathogens responsible for nosocomial infections. In this study opportunistic pathogens present after antibiotic therapy harbored more antibiotic resistant genes than populations of opportunistic pathogens before treatment. Our results demonstrate the effects of antibiotic therapy on induced dysbiosis and expansion of opportunistic pathogen populations and antibiotic resistant subpopulations of those pathogens. Follow-up studies with larger samples sizes and potentially controlled clinical investigations should be performed to confirm our findings.

Agricultural Soil Management Practices Differentially Shape the Bacterial and Fungal Microbiome of Sorghum bicolor

Soils play important roles in biological productivity. While past work suggests that microbes affect soil health and respond to agricultural practices, it is not well known how soil management shapes crop host microbiomes. To elucidate the impact of management on microbial composition and function in the sorghum microbiome, we performed 16S rRNA gene and ITS2 amplicon sequencing and metatranscriptomics on soil and root samples collected from a site in California's San Joaquin Valley that is under long-term cultivation with 1) standard (ST) or no tilling (NT) and 2) cover-cropping (CC) or leaving the field fallow (NO). Our results revealed that microbial diversity, composition, and function change across tillage and cover type, with a heightened response in fungal communities, versus bacterial. Surprisingly, ST harbored greater microbial alpha diversity than NT, indicating that tillage may open niche spaces for broad colonization. Across management regimes, we observed class-level taxonomic level shifts. Additionally, we found significant functional restructuring across treatments, including enrichment for microbial lipid and carbohydrate transport and metabolism and cell motility with NT. Differences in carbon cycling were also observed, with increased prevalence of glycosyltransferase and glycoside hydrolase carbohydrate active enzyme families with CC. Lastly, treatment significantly influenced arbuscular mycorrhizal fungi, which had the greatest prevalence and activity under ST, suggesting that soil practices mediate known beneficial plant-microbe relationships. Collectively, our results demonstrate how agronomic practices impact critical interactions within the plant microbiome and inform future efforts to configure trait-associated microbiomes in crops.Importance While numerous studies show that farming practices can influence the soil microbiome, there are often conflicting results on how microbial diversity and activity respond to treatment. In addition, there is very little work published on how the corresponding crop plant microbiome is impacted. With bacteria and fungi known to critically affect soil health and plant growth, we concurrently compared how the practices of no and standard tillage, in combination with either cover-cropping or fallow fields, shape soil and plant-associated microbiomes between the two classifications. In determining not only the response to treatment in microbial diversity and composition, but for activity as well, this work demonstrates the significance of agronomic practice in modulating plant-microbe interactions, as well as encourages future work on the mechanisms involved in community assemblages supporting similar crop outcomes.

Complete Genome Sequence of Burkholderia cenocepacia Phage Mica

Burkholderia cenocepacia is a multidrug-resistant Gram-negative pathogen known to colonize patients with chronic granulomatous disease and cystic fibrosis. Here, we describe Burkholderia phage Mica, which is predicted to be a lysogenic myophage based on the similarity of its structural proteins to Enterobacteria phage P2 and Burkholderia phage KL3.