Loss of genes related to Nucleotide Excision Repair (NER) and implications for reductive genome evolution in symbionts of deep-sea vesicomyid clams

Divergent paths in the evolutionary history of maternally transmitted clam symbionts

Vertical transmission of bacterial endosymbionts is accompanied by virtually irreversible gene loss that results in a progressive reduction in genome size. While the evolutionary processes of genome reduction have been well described in some terrestrial symbioses, they are less understood in marine systems where vertical transmission is rarely observed. The association between deep-sea vesicomyid clams and chemosynthetic Gammaproteobacteria is one example of maternally inherited symbioses in the ocean. Here, we assessed the contributions of drift, recombination and selection to genome evolution in two extant vesicomyid symbiont clades by comparing 15 representative symbiont genomes (1.017-1.586 Mb) to those of closely related bacteria and the hosts' mitochondria. Our analyses suggest that drift is a significant force driving genome evolution in vesicomyid symbionts, though selection and interspecific recombination appear to be critical for maintaining symbiont functional integrity and creating divergent patterns of gene conservation. Notably, the two symbiont clades possess putative functional differences in sulfide physiology, anaerobic respiration and dependency on environmental vitamin B12, which probably reflect adaptations to different ecological habitats available to each symbiont group. Overall, these results contribute to our understanding of the eco-evolutionary processes shaping reductive genome evolution in vertically transmitted symbioses.

Host hybridization as a potential mechanism of lateral symbiont transfer in deep-sea vesicomyid clams

Deep-sea vesicomyid clams live in mutualistic symbiosis with chemosynthetic bacteria that are inherited through the maternal germ line. On evolutionary timescales, strictly vertical transmission should lead to cospeciation of host mitochondrial and symbiont lineages; nonetheless, examples of incongruent phylogenies have been reported, suggesting that symbionts are occasionally horizontally transmitted between host species. The current paradigm for vesicomyid clams holds that direct transfers cause host shifts or mixtures of symbionts. An alternative hypothesis suggests that hybridization between host species might explain symbiont transfers. Two clam species, Archivesica gigas and Phreagena soyoae, frequently co-occur at deep-sea hydrocarbon seeps in the eastern Pacific Ocean. Although the two species typically host gammaproteobacterial symbiont lineages marked by divergent 16S rRNA phylotypes, we identified a number of clams with the A. gigas mitotype that hosted symbionts with the P. soyoae phylotype. Demographic inference models based on genome-wide SNP data and three Sanger sequenced gene markers provided evidence that A. gigas and P. soyoae hybridized in the past, supporting the hypothesis that hybridization might be a viable mechanism of interspecific symbiont transfer. These findings provide new perspectives on the evolution of vertically transmitted symbionts and their hosts in deep-sea chemosynthetic environments.

Comparative mitochondrial genomic analyses of three chemosynthetic vesicomyid clams from deep-sea habitats

Vesicomyid clams of the subfamily Pliocardinae are among the dominant chemosymbiotic bivalves found in sulfide-rich deep-sea habitats. Plastic morphologies and present molecular data could not resolve taxonomic uncertainties. The complete mitochondrial (mt) genomes will provide more data for comparative studies on molecular phylogeny and systematics of this taxonomically uncertain group, and help to clarify generic classifications. In this study, we analyze the features and evolutionary dynamics of mt genomes from three Archivesica species (Archivesica sp., Ar. gigas and Ar. pacifica) pertaining to subfamily Pliocardinae. Sequence coverage is nearly complete for the three newly sequenced mt genomes, with only the control region and some tRNA genes missing. Gene content, base composition, and codon usage are highly conserved in these pliocardiin species. Comparative analysis revealed the vesicomyid have a relatively lower ratio of Ka/Ks, and all 13 protein-coding genes (PGCs) are under strong purifying selection with a ratio of Ka/Ks far lower than one. Minimal changes in gene arrangement among vesicomyid species are due to the translocation trnaG in Isorropodon fossajaponicum. Additional tRNA genes were detected between trnaG and nad2 in Abyssogena mariana (trnaL3), Ab. phaseoliformis (trnaS3), and Phreagena okutanii (trnaM2), and display high similarity to other pliocardiin sequences at the same location. Single base insertion in multiple sites of this location could result in new tRNA genes, suggesting a possible tRNA arising from nongeneic sequence. Phylogenetic analysis based on 12 PCGs (excluding atp8) supports the monophyly of Pliocardiinae. These nearly complete mitogenomes provide relevant data for further comparative studies on molecular phylogeny and systematics of this taxonomically uncertain group of chemosymbiotic bivalves.

Ancient Occasional Host Switching of Maternally Transmitted Bacterial Symbionts of Chemosynthetic Vesicomyid Clams

Vesicomyid clams in deep-sea chemosynthetic ecosystems harbor sulfur-oxidizing bacteria in their gill epithelial cells. These symbionts, which are vertically transmitted, are species-specific and thought to have cospeciated with their hosts. However, recent studies indicate incongruent phylogenies between some vesicomyid clams and their symbionts, suggesting that symbionts are horizontally transmitted. To more precisely understand the evolution of vesicomyid clams and their symbionts, we compared the evolution of vesicomyid clams and their symbionts through phylogenetic analyses using multi-gene data sets. Many clades in the phylogenetic trees of 13 host species (Abyssogena mariana, Ab. phaseoliformis, Akebiconcha kawamurai, Calyptogena fausta, C. laubieri, C. magnifica, C. nautilei, C. pacifica, Isorropodon fossajaponicum, Phreagena kilmeri, Ph. okutanii, Ph. soyoae, and Pliocardia stearnsii) and their symbionts were well resolved. Six of the 13 host-symbiont pairs (C. fausta, C. magnifica, C. pacifica, Ph. kilmeri, Ph. okutanii, and Ph. soyoae, and their respective symbionts) showed topological congruence. However, the remaining seven pairs (Ak. kawamurai, Ab mariana, Ab. phaseoliformis, C. laubieri, C. nautilei, I. fossajaponicum, and Pl. stearnsii and their corresponding symbionts) showed incongruent topologies, which were supported by the approximately unbiased and Bayes factor tests. Coevolution analyses indicated that six pairs cospeciated, whereas host switching events occurred in the remaining seven pairs. Markedly, multiple host switching events may have occurred in the lineages from the common ancestral symbiont of C. pacifica and C. fausta. Our phylogenetic and coevolution analyses provide additional evidence for host switching during the evolution of vesicomyids.