Uncovering differences in the composition and function of phage communities and phage-bacterium interactions in raw soy sauce.

INTRODUCTION

Although the composition and succession of microbial communities in soy sauce fermentation have been well-characterized, the understanding of phage communities in soy sauce remains limited.

METHODS

This study determined the diversity, taxonomic composition, and predicted function of phage communities and the phage-host interactions in two types of raw soy sauce (Cantonese-type fermentation, NJ; Japanese-type fermentation, PJ) using shotgun metagenomics.

RESULTS AND DISCUSSION

These two raw soy sauces showed differences in phage composition (121 viral operational taxonomic units (vOTUs) in NJ and 387 vOTUs in PJ), with a higher abundance of the family (58.50%) in the NJ phage community and a higher abundance of (33.01%) in PJ. Auxiliary metabolic functional annotation analyses showed that phages in the raw soy sauces mostly encoded genes with unknown functions (accounting for 66.33% of COG profiles), but the NJ sample contained genes mostly annotated to conventional functions related to carbohydrate metabolism (0.74%) and lipid metabolism (0.84%), while the PJ sample presented a higher level of amino acid metabolism functions (0.12%). Thirty auxiliary metabolism genes (AMGs) were identified in phage genomes, which were associated with carbohydrate utilization, cysteine and methionine metabolism, and aspartic acid biosynthesis for the host. To identify phage-host interactions, 30 host genomes (affiliated with 22 genera) were also recruited from the metagenomic dataset. The phage-host interaction analysis revealed a wide range of phage hosts, for which a total of 57 phage contigs were associated with 17 host genomes, with and infected by the most phages. This study provides a comprehensive understanding of the phage community composition, auxiliary metabolic functions, and interactions with hosts in two different types of raw soy sauce.