Insights into the genomic traits of Yersinia frederiksenii, Yersinia intermedia and Yersinia kristensenii isolated from diverse sources in Brazil.

Yersinia is an important genus comprising foodborne, zoonotic and pathogenic bacteria. On the other hand, species of the so-called group Yersinia enterocolitica-like are understudied and mostly characterized as non-pathogenic, despite of some reports of human infections. The present study aimed to provide genomic insights of Yersinia frederiksenii (YF), Yersinia intermedia (YI) and Yersinia kristensenii (YK) isolated worldwide. A total of 22 YF, 20 YI and 14 YK genomes were searched for antimicrobial resistance genes, plasmids, prophages, and virulence factors. Their phylogenomic relatedness was analyzed by Gegenees and core-genome multi-locus sequence typing. Beta-lactam resistance gene blaTEM-116 and five plasmids replicons (pYE854, ColRNAI, ColE10, Col(pHAD28) and IncN3) were detected in less than five genomes. A total of 59 prophages, 106 virulence markers of the Yersinia genus, associated to adherence, antiphagocytosis, exoenzymes, invasion, iron uptake, proteases, secretion systems and the O-antigen, and virulence factors associated to other 20 bacterial genera were detected. Phylogenomic analysis revealed high inter-species distinction and four highly diverse YF clusters. In conclusion, the results obtained through the analyses of YF, YI and YK genomes suggest the virulence potential of these strains due to the broad diversity and high frequency of prophages and virulence factors found. Phylogenetic analyses were able to correctly distinguish these closely related species and show the presence of different genetic subgroups. These data contributed for a better understanding of YF, YI and YK virulence-associated features and global genetic diversity, and reinforced the need for better characterization of these Y. enterocolitica-like species considered non-pathogenic.

Genomic and taxonomic evaluation of 38 Treponema prophage sequences.

BACKGROUND: Despite Spirochetales being a ubiquitous and medically important order of bacteria infecting both humans and animals, there is extremely limited information regarding their bacteriophages. Of the genus Treponema, there is just a single reported characterised prophage. RESULTS: We applied a bioinformatic approach on 24 previously published Treponema genomes to identify and characterise putative treponemal prophages. Thirteen of the genomes did not contain any detectable prophage regions. The remaining eleven contained 38 prophage sequences, with between one and eight putative prophages in each bacterial genome. The prophage regions ranged from 12.4 to 75.1 kb, with between 27 and 171 protein coding sequences. Phylogenetic analysis revealed that 24 of the prophages formed three distinct sequence clusters, identifying putative myoviral and siphoviral morphology. ViPTree analysis demonstrated that the identified sequences were novel when compared to known double stranded DNA bacteriophage genomes. CONCLUSIONS: In this study, we have started to address the knowledge gap on treponeme bacteriophages by characterising 38 prophage sequences in 24 treponeme genomes. Using bioinformatic approaches, we have been able to identify and compare the prophage-like elements with respect to other bacteriophages, their gene content, and their potential to be a functional and inducible bacteriophage, which in turn can help focus our attention on specific prophages to investigate further.

The complete genome sequence of unculturable Mycoplasma faucium obtained through clinical metagenomic next-generation sequencing.

Introduction: Diagnosing Mycoplasma faucium poses challenges, and it's unclear if its rare isolation is due to infrequent occurrence or its fastidious nutritional requirements. Methods: This study analyzes the complete genome sequence of M. faucium, obtained directly from the pus of a sternum infection in a lung transplant patient using metagenomic sequencing. Results: Genome analysis revealed limited therapeutic options for the M. faucium infection, primarily susceptibility to tetracyclines. Three classes of mobile genetic elements were identified: two new insertion sequences, a new prophage (phiUMCG-1), and a species-specific variant of a mycoplasma integrative and conjugative element (MICE). Additionally, a Type I Restriction-Modification system was identified, featuring 5'-terminally truncated hsdS pseudogenes with overlapping repeats, indicating the potential for forming alternative hsdS variants through recombination. Conclusion: This study represents the first-ever acquisition of a complete circularized bacterial genome directly from a patient sample obtained from invasive infection of a primary sterile site using culture-independent, PCR-free clinical metagenomics.

Slings and arrows: sRNAs mediate intragenomic competition.

Bacterial genomes are littered with exogenous: competing DNA elements. Here, Sprenger et al. demonstrate that the Vibrio cholerae prophage VP882 modulates host functions via production of regulatory sRNAs to promote phage development. Alternatively, host sRNAs inhibit the VP882 lytic phase by specifically regulating phage genes.

Genetic recombination-mediated evolutionary interactions between phages of potential industrial importance and prophages of their hosts within or across the domains of Escherichia, Listeria, Salmonella, Campylobacter, and Staphylococcus.

BACKGROUND: The in-depth understanding of the role of lateral genetic transfer (LGT) in phage-prophage interactions is essential to rationalizing phage applications for human and animal therapy, as well as for food and environmental safety. This in silico study aimed to detect LGT between phages of potential industrial importance and their hosts. METHODS: A large array of genetic recombination detection algorithms, implemented in SplitsTree and RDP4, was applied to detect LGT between various Escherichia, Listeria, Salmonella, Campylobacter, Staphylococcus, Pseudomonas, and Vibrio phages and their hosts. PHASTER and RAST were employed respectively to identify prophages across the host genome and to annotate LGT-affected genes with unknown functions. PhageAI was used to gain deeper insights into the life cycle history of recombined phages. RESULTS: The split decomposition inferences (bootstrap values: 91.3-100; fit: 91.433-100), coupled with the Phi (0.0-2.836E-12) and RDP4 (P being well below 0.05) statistics, provided strong evidence for LGT between certain Escherichia, Listeria, Salmonella, and Campylobacter virulent phages and prophages of their hosts. The LGT events entailed mainly the phage genes encoding for hypothetical proteins, while some of these genetic loci appeared to have been affected even by intergeneric recombination in specific E. coli and S. enterica virulent phages when interacting with their host prophages. Moreover, it is shown that certain L. monocytogenes virulent phages could serve at least as the donors of the gene loci, involved in encoding for the basal promoter specificity factor, for L. monocytogenes. In contrast, the large genetic clusters were determined to have been simultaneously exchanged by many S. aureus prophages and some Staphylococcus temperate phages proposed earlier as potential therapeutic candidates (in their native or modified state). The above genetic clusters were found to encompass multiple genes encoding for various proteins, such as e.g., phage tail proteins, the capsid and scaffold proteins, holins, and transcriptional terminator proteins. CONCLUSIONS: It is suggested that phage-prophage interactions, mediated by LGT (including intergeneric recombination), can have a far-reaching impact on the co-evolutionary trajectories of industrial phages and their hosts especially when excessively present across microbially rich environments.

Exploring the diversity and evolutionary strategies of prophages in Hyphomicrobiales, comparing animal-associated with non-animal-associated bacteria.

The Hyphomicrobiales bacterial order (previously Rhizobiales) exhibits a wide range of lifestyle characteristics, including free-living, plant-association, nitrogen-fixing, and association with animals (Bartonella and Brucella). This study explores the diversity and evolutionary strategies of bacteriophages within the Hyphomicrobiales order, comparing animal-associated (AAB) with non-animal-associated bacteria (NAAB). We curated 560 high-quality complete genomes of 58 genera from this order and used the PHASTER server for prophage annotation and classification. For 19 genera with representative genomes, we curated 96 genomes and used the Defense-Finder server to summarize the type of anti-phage systems (APS) found in this order. We analyzed the genetic repertoire and length distributions of prophages, estimating evolutionary rates and comparing intact, questionable, and incomplete prophages in both groups. Analyses of best-fit parameters and bootstrap sensitivity were used to understand the evolutionary processes driving prophage gene content. A total of 1860 prophages distributed in Hyphomicrobiales were found, 695 in AAB and 1165 in the NAAB genera. The results revealed a similar number of prophages per genome in AAB and NAAB and a similar length distribution, suggesting shared mechanisms of genetic acquisition of prophage genes. Changes in the frequency of specific gene classes were observed between incomplete and intact prophages, indicating preferential loss or enrichment in both groups. The analysis of best-fit parameters and bootstrap sensitivity tests indicated a higher selection coefficient, induction rate, and turnover in NAAB genomes. We found 68 types of APS in Hyphomicrobiales; restriction modification (RM) and abortive infection (Abi) were the most frequent APS found for all Hyphomicrobiales, and within the AAB group. This classification of APS showed that NAAB genomes have a greater diversity of defense systems compared to AAB, which could be related to the higher rates of prophage induction and turnover in the latter group. Our study provides insights into the distributions of both prophages and APS in Hyphomicrobiales genomes, demonstrating that NAAB carry more defense systems against phages, while AAB show increased prophage stability and an increased number of incomplete prophages. These results suggest a greater role for domesticated prophages within animal-associated bacteria in Hyphomicrobiales.

Microbial adaptation to spaceflight is correlated with bacteriophage-encoded functions.

Evidence from the International Space Station suggests microbial populations are rapidly adapting to the spacecraft environment; however, the mechanism of this adaptation is not understood. Bacteriophages are prolific mediators of bacterial adaptation on Earth. Here we survey 245 genomes sequenced from bacterial strains isolated on the International Space Station for dormant (lysogenic) bacteriophages. Our analysis indicates phage-associated genes are significantly different between spaceflight strains and their terrestrial counterparts. In addition, we identify 283 complete prophages, those that could initiate bacterial lysis and infect additional hosts, of which 21% are novel. These prophage regions encode functions that correlate with increased persistence in extreme environments, such as spaceflight, to include antimicrobial resistance and virulence, DNA damage repair, and dormancy. Our results correlate microbial adaptation in spaceflight to bacteriophage-encoded functions that may impact human health in spaceflight.

Prophages carrying Zot toxins on different Vibrio genomes: A comprehensive assessment using multilayer networks.

Vibrios, a group of bacteria that are among the most abundant in marine environments, include several species such as Vibrio cholerae and Vibrio parahaemolyticus, which can be pathogenic to humans. Some species of Vibrio contain prophages within their genomes. These prophages can carry genes that code for toxins, such as the zonula occludens toxin (Zot), which contribute to bacterial virulence. Understanding the association between different Vibrio species, prophages and Zot genes can provide insights into their ecological interactions. In this study, we evaluated 4619 Vibrio genomes from 127 species to detect the presence of prophages carrying the Zot toxin. We found 2030 potential prophages with zot-like genes in 43 Vibrio species, showing a non-random association within a primarily modular interaction network. Some prophages, such as CTX or Vf33, were associated with specific species. In contrast, prophages phiVCY and VfO3K6 were found in 28 and 20 Vibrio species, respectively. We also identified six clusters of Zot-like sequences in prophages, with the ZOT2 cluster being the most frequent, present in 34 Vibrio species. This analysis helps to understand the distribution patterns of zot-containing prophages across Vibrio genomes and the potential routes of Zot-like toxin dissemination.

Genomic analysis of an outbreak of Shiga toxin-producing Escherichia coli O183:H18 in the United Kingdom, 2023.

In June 2023, UKHSA surveillance systems detected an outbreak of severe gastrointestinal symptoms caused by a rare serotype of Shiga toxin-producing Escherichia coli, STEC O183:H18. There were 26 cases aged 6 months to 74 years (42 % cases were aged 0-9 years), distributed across the UK with onset dates range between 22 May 2023 and 4 July 2023. The epidemiological and food chain investigations were inconclusive, although meat products made from beef mince were implicated as a potential vehicle. The outbreak strain belonged to sequence type (ST) 657 and harboured a Shiga toxin (stx) subtype stx2a located on a prophage that was unique in the UKHSA stx-encoding bacteriophage database. Plasmid encoded, putative virulence genes subA, ehxA, saa, iha, lpfA and iss were detected, however, the established STEC virulence genes involved in attachment to the gut mucosa (eae and aggR) were absent. The acquisition of stx across the global population structure of ST657 appeared to correspond with the presence of subA, ehxA, saa, iha, lpfA and iss. During the outbreak investigation, we used long read sequencing to characterise the plasmid and prophage content of this atypical STEC, to look for evidence to explain its recent emergence. Although we were unable to determine source and transmission route of the outbreak strain, the genomic analysis revealed potential clues as to how novel strains for STEC evolve. With the implementation of PCR capable of detecting all STEC, and genome sequencing for typing and virulence profiling, we have the tools to enable us to monitor the changing landscape of STEC. Improvements in the standardised collection of epidemiological data and trace-back strategies within the food industry, will ensure we have a surveillance system capable of alerting us to emerging threats to public health.

Ligand cross-feeding resolves bacterial vitamin B12 auxotrophies.

Cobalamin (vitamin B12, herein referred to as B12) is an essential cofactor for most marine prokaryotes and eukaryotes1,2. Synthesized by a limited number of prokaryotes, its scarcity affects microbial interactions and community dynamics2-4. Here we show that two bacterial B12 auxotrophs can salvage different B12 building blocks and cooperate to synthesize B12. A Colwellia sp. synthesizes and releases the activated lower ligand α-ribazole, which is used by another B12 auxotroph, a Roseovarius sp., to produce the corrin ring and synthesize B12. Release of B12 by Roseovarius sp. happens only in co-culture with Colwellia sp. and only coincidently with the induction of a prophage encoded in Roseovarius sp. Subsequent growth of Colwellia sp. in these conditions may be due to the provision of B12 by lysed cells of Roseovarius sp. Further evidence is required to support a causative role for prophage induction in the release of B12. These complex microbial interactions of ligand cross-feeding and joint B12 biosynthesis seem to be widespread in marine pelagic ecosystems. In the western and northern tropical Atlantic Ocean, bacteria predicted to be capable of salvaging cobinamide and synthesizing only the activated lower ligand outnumber B12 producers. These findings add new players to our understanding of B12 supply to auxotrophic microorganisms in the ocean and possibly in other ecosystems.

Active prophages in coral-associated Halomonas capable of lateral transduction.

Temperate phages can interact with bacterial hosts through lytic and lysogenic cycles via different mechanisms. Lysogeny has been identified as the major form of bacteria-phage interaction in the coral-associated microbiome. However, the lysogenic-to-lytic switch of temperate phages in ecologically important coral-associated bacteria and its ecological impact have not been extensively investigated. By studying the prophages in coral-associated Halomonas meridiana, we found that two prophages, Phm1 and Phm3, are inducible by the DNA-damaging agent mitomycin C and that Phm3 is spontaneously activated under normal cultivation conditions. Furthermore, Phm3 undergoes an atypical lytic pathway that can amplify and package adjacent host DNA, potentially resulting in lateral transduction. The induction of Phm3 triggered a process of cell lysis accompanied by the formation of outer membrane vesicles (OMVs) and Phm3 attached to OMVs. This unique cell-lysis process was controlled by a four-gene lytic module within Phm3. Further analysis of the Tara Ocean dataset revealed that Phm3 represents a new group of temperate phages that are widely distributed and transcriptionally active in the ocean. Therefore, the combination of lateral transduction mediated by temperate phages and OMV transmission offers a versatile strategy for host-phage coevolution in marine ecosystems.

Integrative and Conjugative Elements and Prophage DNA as Carriers of Resistance Genes in Erysipelothrix rhusiopathiae Strains from Domestic Geese in Poland.

Goose erysipelas is a serious problem in waterfowl breeding in Poland. However, knowledge of the characteristics of Erysipelothrix rhusiopathiae strains causing this disease is limited. In this study, the antimicrobial susceptibility and serotypes of four E. rhusiopathiae strains from domestic geese were determined, and their whole-genome sequences (WGSs) were analyzed to detect resistance genes, integrative and conjugative elements (ICEs), and prophage DNA. Sequence type and the presence of resistance genes and transposons were compared with 363 publicly available E. rhusiopathiae strains, as well as 13 strains of other Erysipelothrix species. Four strains tested represented serotypes 2 and 5 and the MLST groups ST 4, 32, 242, and 243. Their assembled circular genomes ranged from 1.8 to 1.9 kb with a GC content of 36-37%; a small plasmid was detected in strain 1023. Strains 1023 and 267 were multidrug-resistant. The resistance genes detected in the genome of strain 1023 were erm47, tetM, and lsaE-lnuB-ant(6)-Ia-spw cluster, while strain 267 contained the tetM and ermB genes. Mutations in the gyrA gene were detected in both strains. The tetM gene was embedded in a Tn916-like transposon, which in strain 1023, together with the other resistance genes, was located on a large integrative and conjugative-like element of 130 kb designated as ICEEr1023. A minor integrative element of 74 kb was identified in strain 1012 (ICEEr1012). This work contributes to knowledge about the characteristics of E. rhusiopathiae bacteria and, for the first time, reveals the occurrence of erm47 and ermB resistance genes in strains of this species. Phage infection appears to be responsible for the introduction of the ermB gene into the genome of strain 267, while ICEs most likely play a key role in the spread of the other resistance genes identified in E. rhusiopathiae.

Assessing phage-host population dynamics by reintroducing virulent viruses to synthetic microbiomes.

Microbiomes feature complex interactions between diverse bacteria and bacteriophages. Synthetic microbiomes offer a powerful way to study these interactions; however, a major challenge is obtaining a representative bacteriophage population during the bacterial isolation process. We demonstrate that colony isolation reliably excludes virulent viruses from sample sources with low virion-to-bacteria ratios such as feces, creating "virulent virus-free" controls. When the virulent dsDNA virome is reintroduced to a 73-strain synthetic gut microbiome in a bioreactor model of the human colon, virulent viruses target susceptible strains without significantly altering community structure or metabolism. In addition, we detected signals of prophage induction that associate with virulent predation. Overall, our findings indicate that dilution-based isolation methods generate synthetic gut microbiomes that are heavily depleted, if not devoid, of virulent viruses and that such viruses, if reintroduced, have a targeted effect on community assembly, metabolism, and prophage replication.

Prophage terminase with tRNase activity sensitizes Salmonella enterica to oxidative stress.

Phage viruses shape the evolution and virulence of their bacterial hosts. The Salmonella enterica genome encodes several stress-inducible prophages. The Gifsy-1 prophage terminase protein, whose canonical function is to process phage DNA for packaging in the virus head, unexpectedly acts as a transfer ribonuclease (tRNase) under oxidative stress, cleaving the anticodon loop of tRNALeu. The ensuing RNA fragmentation compromises bacterial translation, intracellular survival, and recovery from oxidative stress in the vertebrate host. S. enterica adapts to this transfer RNA (tRNA) fragmentation by transcribing the RNA repair Rtc system. The counterintuitive translational arrest provided by tRNA cleavage may subvert prophage mobilization and give the host an opportunity for repair as a way of maintaining bacterial genome integrity and ultimately survival in animals.

Targeted deletion of Pf prophages from diverse Pseudomonas aeruginosa isolates has differential impacts on quorum sensing and virulence traits.

Pseudomonas aeruginosa is an opportunistic bacterial pathogen that commonly causes medical hardware, wound, and respiratory infections. Temperate filamentous Pf phages that infect P. aeruginosa impact numerous virulence phenotypes. Most work on Pf phages has focused on Pf4 and its host P. aeruginosa PAO1. Expanding from Pf4 and PAO1, this study explores diverse Pf phages infecting P. aeruginosa clinical isolates. We describe a simple technique targeting the Pf lysogeny maintenance gene, pflM (PA0718), that enables the effective elimination of Pf prophages from diverse P. aeruginosa hosts. The pflM gene shows diversity among different Pf phage isolates; however, all examined pflM alleles encode the DUF5447 domain. We demonstrate that pflM deletion results in prophage excision but not replication, leading to total prophage loss, indicating a role for lysis/lysogeny decisions for the DUF5447 domain. This study also assesses the effects different Pf phages have on host quorum sensing, biofilm formation, pigment production, and virulence against the bacterivorous nematode Caenorhabditis elegans. We find that Pf phages have strain-specific impacts on quorum sensing and biofilm formation, but nearly all suppress pigment production and increase C. elegans avoidance behavior. Collectively, this research not only introduces a valuable tool for Pf prophage elimination from diverse P. aeruginosa isolates but also advances our understanding of the complex relationship between P. aeruginosa and filamentous Pf phages.IMPORTANCEPseudomonas aeruginosa is an opportunistic bacterial pathogen that is frequently infected by filamentous Pf phages (viruses) that integrate into its chromosome, affecting behavior. Although prior work has focused on Pf4 and PAO1, this study investigates diverse Pf in clinical isolates. A simple method targeting the deletion of the Pf lysogeny maintenance gene pflM (PA0718) effectively eliminates Pf prophages from clinical isolates. The research evaluates the impact Pf prophages have on bacterial quorum sensing, biofilm formation, and virulence phenotypes. This work introduces a valuable tool to eliminate Pf prophages from clinical isolates and advances our understanding of P. aeruginosa and filamentous Pf phage interactions.

Prophage maintenance is determined by environment-dependent selective sweeps rather than mutational availability.

Prophages, viral sequences integrated into bacterial genomes, can be beneficial and costly. Despite the risk of prophage activation and subsequent bacterial death, active prophages are present in most bacterial genomes. However, our understanding of the selective forces that maintain prophages in bacterial populations is limited. Combining experimental evolution with stochastic modeling, we show that prophage maintenance and loss are primarily determined by environmental conditions that alter the net fitness effect of a prophage on its bacterial host. When prophages are too costly, they are rapidly lost through environment-specific sequences of selective sweeps. Conflicting selection pressures that select against the prophage but for a prophage-encoded accessory gene can maintain prophages. The dynamics of prophage maintenance additionally depend on the sociality of this accessory gene. Prophage-encoded genes that exclusively benefit the lysogen maintain prophages at higher frequencies compared with genes that benefit the entire population. That is because the latter can protect phage-free "cheaters," reducing the benefit of maintaining the prophage. Our simulations suggest that environmental variation plays a larger role than mutation rates in determining prophage maintenance. These findings highlight the complexity of selection pressures that act on mobile genetic elements and challenge our understanding of the role of environmental factors relative to random chance events in shaping the evolutionary trajectory of bacterial populations. By shedding light on the key factors that shape microbial populations in the face of environmental changes, our study significantly advances our understanding of the complex dynamics of microbial evolution and diversification.

Phage biology: The ins and outs of prophages in bacterial populations.

Bacterial genomes often harbor integrated viruses (prophages), which provide novel functions but also lyse cells under stressful conditions. A new paper combines mathematical models with experimental evolution to determine how prophages are maintained in bacterial populations despite their fitness costs.

Bacteriophages avoid autoimmunity from cognate immune systems as an intrinsic part of their life cycles.

Dormant prophages protect lysogenic cells by expressing diverse immune systems, which must avoid targeting their cognate prophages upon activation. Here we report that multiple Staphylococcus aureus prophages encode Tha (tail-activated, HEPN (higher eukaryotes and prokaryotes nucleotide-binding) domain-containing anti-phage system), a defence system activated by structural tail proteins of incoming phages. We demonstrate the function of two Tha systems, Tha-1 and Tha-2, activated by distinct tail proteins. Interestingly, Tha systems can also block reproduction of the induced tha-positive prophages. To prevent autoimmunity after prophage induction, these systems are inhibited by the product of a small overlapping antisense gene previously believed to encode an excisionase. This genetic organization, conserved in S. aureus prophages, allows Tha systems to protect prophages and their bacterial hosts against phage predation and to be turned off during prophage induction, balancing immunity and autoimmunity. Our results show that the fine regulation of these processes is essential for the correct development of prophages' life cycle.

Genome-based characterization of a novel prophage of Vibrio parahaemolyticus, VPS05ph1, a novel member of Peduoviridae.

Vibrio parahaemolyticus is a globally important bacterium related to climate warming and health threat to human and marine animals. Yet, there is limited knowledge about its polylysogeny harboring multiple prophages and the genetic information. In this study, two prophages (VPS05ph1 and VPS05ph2) were identified in a V. parahaemolyticus isolate through genomic and transcriptional analyses. Both prophages were determined as HP1-like phages, located in a novel phylogenetic lineage of Peduoviridae. They shared a moderate genome-wide sequence similarity with each other and high synteny with the closest relatives, but showed low identities to the repressor counterparts of the representative phages within the family. In addition, no bacterial virulence genes, antibiotic resistance genes and known phage-encoded lytic proteins were identified on both prophage genomes. Moreover, the V. parahaemolyticus isolate was induced with mitomycin, which caused aberrant cellular morphology and nonviability of bacterial cells and excision of prophage VPS05ph1, accompanied by the respective inhibition and promotion of transcriptions of the cI-like and cox-like regulator genes for phage decision making. Results in this study provide the genetic context of polylysogeny in the V. parahaemolyticus isolate, support the diversity and prevalence of HP1-like phages in vibrios, and promote to explore interactions between the HP1-like prophage and its vibrio host.

Genome Analysis of Epsilon CrAss-like Phages.

CrAss-like phages play an important role in maintaining ecological balance in the human intestinal microbiome. However, their genetic diversity and lifestyle are still insufficiently studied. In this study, a novel CrAssE-Sib phage genome belonging to the epsilon crAss-like phage genomes was found. Comparative analysis indicated that epsilon crAss-like phages are divided into two putative genera, which were proposed to be named Epsilonunovirus and Epsilonduovirus; CrAssE-Sib belongs to the former. The crAssE-Sib genome contains a diversity-generating retroelement (DGR) cassette with all essential elements, including the reverse transcriptase (RT) and receptor binding protein (RBP) genes. However, this RT contains the GxxxSP motif in its fourth domain instead of the usual GxxxSQ motif found in all known phage and bacterial DGRs. RBP encoded by CrAssE-Sib and other Epsilonunoviruses has an unusual structure, and no similar phage proteins were found. In addition, crAssE-Sib and other Epsilonunoviruses encode conserved prophage repressor and anti-repressors that could be involved in lysogenic-to-lytic cycle switches. Notably, DNA primase sequences of epsilon crAss-like phages are not included in the monophyletic group formed by the DNA primases of all other crAss-like phages. Therefore, epsilon crAss-like phage substantially differ from other crAss-like phages, indicating the need to classify these phages into a separate family.