The genomic content and context of auxiliary metabolic genes in roseophages

Marine bacteriophages frequently possess auxiliary metabolic genes (AMGs) that accelerate host metabolism during phage infection. The significance of AMGs in phage infecting the ecologically important Roseobacter clade, found predominantly in marine environments, remains to be determined. Here, we analysed the distribution and genomic context of 180 AMGs, annotated into 20 types, across 50 roseophage genomes. Roseophages share seven high-frequency AMGs (trx, grx, RNR, thyX, DCD, phoH, and mazG)， most of them involved in the nucleotide biosynthesis pathway that represent conserved intra and inter operational taxonomic units (OTUs), and share ≥97% full-length DNA sequence similarity. Sporadic AMGs (dUTPase, lexA, degS, Que, NAPRT, AHL, pcnB, ctrA, RTX, RNR-nrdA, RNR-nrdE, wclP, and flgJ), present in only one or two OTUs, show high functional diversity. The roseophage AMG repertoire weakly correlates with environmental factors, while host range partially explains the sporadic AMG distribution. Locally co-linear blocks distribution index (LDI) analysis indicated that high-frequency roseopodovirus AMGs are restricted to particular genomic islands, possibly originating from limited historical acquisition events. Low-frequency roseopodovirus AMGs and all roseosiphovirus AMGs have high LDI values, implying multiple historical acquisition events. In summary, roseophages have acquired a range of AMGs through horizontal gene transfer, and the forces shaping the evolution of roseophages are described.