Abstract
Reef-building corals depend on an intricate community of microorganisms for functioning and resilience. The infection of coral-associated bacteria by bacteriophages can modify bacterial ecological interactions, yet very little is known about phage functions in the holobiont. This gap stems from methodological limitations that have prevented the recovery of high-quality viral genomes and bacterial host assignment from coral samples. Here, we introduce a size fractionation approach that increased bacterial and viral recovery in coral metagenomes by 9-fold and 2-fold, respectively, and enabled the assembly and binning of bacterial and viral genomes at relatively low sequencing coverage. We combined these viral genomes with those derived from 677 publicly available metagenomes, viromes, and bacterial isolates from stony corals to build a global coral virus database of over 20,000 viral genomic sequences spanning four viral realms. The tailed bacteriophage families Kyanoviridae and Autographiviridae were the most abundant, replacing groups formerly referred to as Myoviridae and Podoviridae, respectively. Prophage and CRISPR spacer linkages between these viruses and 626 bacterial metagenome-assembled genomes and bacterial isolates showed that most viruses infected Alphaproteobacteria, the most abundant class, and less abundant taxa like Halanaerobiia and Bacteroidia. A host-phage-gene network identified keystone viruses with the genomic capacity to modulate bacterial metabolic pathways and direct molecular interactions with eukaryotic cells. This study reveals the genomic basis of nested symbioses between bacteriophage, bacteria, and the coral host and its endosymbiotic algae.