Fig. 2

Genetic divergence of methane-oxidizing endosymbionts between vent and seep mussels. a ML phylogenetic reconstruction of the assembled endosymbiont genomes based on orthologous genes conserved across the pangenome. A total of five PacBio HiFi sequencing datasets were used for genome assembly. Every dataset was assembled separately, and the squares prepending each sequence label correspond to specific datasets derived from seep and vent mussel individuals. b GC content of the assembled genomes in the three clades. c Bin size of the assembled genomes in the three clades. d Density plot of the dN/dS values for each clade estimated from 1000 concatenated alignments constructed from twenty-five randomly chosen orthologs. Dash line showing the mean value of each clade. e Distribution of endosymbiont abundance across individual mussels collected from vent (Iheya North Knoll, IN, and Daiyon-Yonaguni Knoll, DY) and seep (Formosa Ridge, FR) environments