The congruence between host and symbiont phylogenies reflects the evolutionary links among ecologically important interactions. As potential key symbionts, the members affiliated to the family Endozoicomonadaceae have previously been investigated for the cophylogenetic relationship with their hosts using their 16S rRNA gene sequences. However, this approach neglects the genomic features of symbionts that may influence the long-term associations between Endozoicomonadaceae members and their hosts. Here, we collected available high-quality genomes of Endozoicomonadaceae from diverse hosts and investigated their genomic features, including genome size, phages, insertion elements and the composition of functional genes. We also tested the host–Endozoicomonadaceae cophylogeny and examined the correlation between the cophylogenetic squared residuals and the genomic features of Endozoicomonadaceae members. Our results revealed a cophylogenetic pattern between members of the Endozoicomonadaceae family and their various hosts. Moreover, we found that the investigated genomes of Endozoicomonadaceae members were differentially eroded by phages and insertion elements. Additionally, Endozoicomonadaceae members with smaller, more eroded genomes tended to exhibit lower cophylogenetic residuals with their hosts. Gene function analysis further revealed that Endozoicomonadaceae members with closer associations with their hosts carried specific genes related to infection processes and host–symbiont interactions. This study suggests that the genomic features of Endozoicomonadaceae members may influence long-term host–Endozoicomonadaceae intimate associations.

Evolutionary history of the host may influence the skeletal morphology of scleractinian corals. However, its effects on the assembly and function of endolithic microbiomes remain unknown. We analyzed bacterial and archaeal microbiomes from the coral skeleton by using 16S rRNA gene sequencing. We collected the samples of seven coral genera distributed among the diverse “Complex” and “Robust” clades. In this study, bacterial α-diversity was significantly higher in the Complex clade relative to the Robust clade. Archaea, on the other hand, remained stable and showed no significant differences between the two host clades, and were most abundantly Nanoarchaeota and Thermoproteota. Analysis of the network topologies showed that network structures were different between the Complex group and the Robust clade. The Robust clade formed a dense and closely knit network among bacteria and archaea. The Com-plex group formed a more modular network structure. Functional predictions further highlighted lineage-specific metabolic strategies. Enrichment was apparent in both nitrification genes (amoB, amoC) and denitrification genes (nirK, nirS) in the Complex clade. This suggests that the coupling of these nitrogen cycles is possible. The opposite was observed for the Robust clade, which had low potential for both types of nitrogen cycling. This reflects the degree of diffusion limitation in the more massive skeleton of this host lineage. Overall, species evolutionary lineage is a pre-eminent driver for the selective filtering of endolithic assembly. It generates discrete skeletal micro-niches on which microbial strategies diverge. In particular, Complex corals favor fast metabolic flux, and Robust corals favor strong network connectivity.