Search

Bacterial symbionts, with their shorter generation times and capacity for horizontal gene transfer (HGT), play a critical role in allowing marine organisms to cope with environmental change. The closure of the Isthmus of Panama created distinct environmental conditions in the Tropical Eastern Pacific (TEP) and Caribbean, offering a "natural experiment" for studying how closely related animals evolve and adapt under environmental change. However, the role of bacterial symbionts in this process is often overlooked. We sequenced the genomes of endosymbiotic bacteria in two sets of sister species of chemosymbiotic bivalves from the genera Codakia and Ctena (family Lucinidae) collected on either side of the Isthmus, to investigate how differing environmental conditions have influenced the selection of symbionts and their metabolic capabilities. The lucinid sister species hosted different Candidatus Thiodiazotropha symbionts and only those from the Caribbean had the genetic potential for nitrogen fixation, while those from the TEP did not. Interestingly, this nitrogen-fixing ability did not correspond to symbiont phylogeny, suggesting convergent evolution of nitrogen fixation potential under nutrient-poor conditions. Reconstructing the evolutionary history of the nifHDKT operon by including other lucinid symbiont genomes from around the world further revealed that the last common ancestor (LCA) of Ca. Thiodiazotropha lacked nif genes, and populations in oligotrophic habitats later re-acquired the nif operon through HGT from the Sedimenticola symbiont lineage. Our study suggests that HGT of the nif operon has facilitated niche diversification of the globally distributed Ca. Thiodiazotropha endolucinida species clade. It highlights the importance of nitrogen availability in driving the ecological diversification of chemosynthetic symbiont species and the role that bacterial symbionts may play in the adaptation of marine organisms to changing environmental conditions., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Morel-Letelier et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Advances in the analysis of amplicon sequence datasets have introduced a methodological shift in how research teams investigate microbial biodiversity, away from sequence identity-based clustering (producing Operational Taxonomic Units, OTUs) to denoising methods (producing amplicon sequence variants, ASVs). While denoising methods have several inherent properties that make them desirable compared to clustering-based methods, questions remain as to the influence that these pipelines have on the ecological patterns being assessed, especially when compared to other methodological choices made when processing data (e.g. rarefaction) and computing diversity indices. We compared the respective influences of two widely used methods, namely DADA2 (a denoising method) vs. Mothur (a clustering method) on 16S rRNA gene amplicon datasets (hypervariable region v4), and compared such effects to the rarefaction of the community table and OTU identity threshold (97% vs. 99%) on the ecological signals detected. We used a dataset comprising freshwater invertebrate (three Unionidae species) gut and environmental (sediment, seston) communities sampled in six rivers in the southeastern USA. We ranked the respective effects of each methodological choice on alpha and beta diversity, and taxonomic composition. The choice of the pipeline significantly influenced alpha and beta diversities and changed the ecological signal detected, especially on presence/absence indices such as the richness index and unweighted Unifrac. Interestingly, the discrepancy between OTU and ASV-based diversity metrics could be attenuated by the use of rarefaction. The identification of major classes and genera also revealed significant discrepancies across pipelines. Compared to the pipeline's effect, OTU threshold and rarefaction had a minimal impact on all measurements., Competing Interests: The authors have declared that no competing interests exist.

Bacterial symbionts bring a wealth of functions to the associations they participate in, but by doing so, they endanger the genes and genomes underlying these abilities. When bacterial symbionts become obligately associated with their hosts, their genomes are thought to decay towards an organelle-like fate due to decreased homologous recombination and inefficient selection. However, numerous associations exist that counter these expectations, especially in marine environments, possibly due to ongoing horizontal gene flow. Despite extensive theoretical treatment, no empirical study thus far has connected these underlying population genetic processes with long-term evolutionary outcomes. By sampling marine chemosynthetic bacterial-bivalve endosymbioses that range from primarily vertical to strictly horizontal transmission, we tested this canonical theory. We found that transmission mode strongly predicts homologous recombination rates, and that exceedingly low recombination rates are associated with moderate genome degradation in the marine symbionts with nearly strict vertical transmission. Nonetheless, even the most degraded marine endosymbiont genomes are occasionally horizontally transmitted and are much larger than their terrestrial insect symbiont counterparts. Therefore, horizontal transmission and recombination enable efficient natural selection to maintain intermediate symbiont genome sizes and substantial functional genetic variation., Competing Interests: The authors have declared that no competing interests exist.

Apicomplexan parasites cause severe disease in both humans and their domesticated animals. Since these parasites readily develop drug resistance, development of new, effective drugs to treat infection caused by these parasites is an ongoing challenge for the medical and veterinary communities. We hypothesized that invertebrate-bacterial symbioses might be a rich source of anti-apicomplexan compounds because invertebrates are susceptible to infections with gregarines, parasites that are ancestral to all apicomplexans. We chose to explore the therapeutic potential of shipworm symbiotic bacteria as they are bona fide symbionts, are easily grown in axenic culture and have genomes rich in secondary metabolite loci [1,2]. Two strains of the shipworm symbiotic bacterium, Teredinibacter turnerae, were screened for activity against Toxoplasma gondii and one strain, T7901, exhibited activity against intracellular stages of the parasite. Bioassay-guided fractionation identified tartrolon E (trtE) as the source of the activity. TrtE has an EC50 of 3 nM against T. gondii, acts directly on the parasite itself and kills the parasites after two hours of treatment. TrtE exhibits nanomolar to picomolar level activity against Cryptosporidium, Plasmodium, Babesia, Theileria, and Sarcocystis; parasites representing all branches of the apicomplexan phylogenetic tree. The compound also proved effective against Cryptosporidium parvum infection in neonatal mice, indicating that trtE may be a potential lead compound for preclinical development. Identification of a promising new compound after such limited screening strongly encourages further mining of invertebrate symbionts for new anti-parasitic therapeutics., Competing Interests: The authors have declared that no competing interests exist.

Vibrio coralliilyticus is a pathogen of coral and mollusk, contributing to dramatic losses worldwide. In our study, we found that V. coralliilyticus challenge could directly affect adult Tridacna crocea survival; there were dead individuals appearing at 6 h post infection, and there were 45.56% and 56.78% mortality rates in challenged groups after 36 h of infection. The apoptosis rate of hemocytes was significantly increased by 1.8-fold at 6 h after V. coralliilyticus injection. To shed light on the mechanistic molecular responses of T. crocea to V. coralliilyticus infection, we used transcriptome sequencing analysis and other relevant techniques to analyze T. crocea hemocytes at 0 h, 6 h, 12 h and 24 h after V. coralliilyticus challenge. Our results revealed that the total numbers of unigenes and DEGs were 195651 and 3446, respectively. Additional details were found by KEGG pathway enrichment analysis, where DEGs were significantly enriched in immune-related signaling pathways, such as the TLR signaling pathway, and some were associated with signaling related to apoptosis. Quantitative validation results illustrated that with exposure to V. coralliilyticus, the expression of TLR pathway members, TLR, MyD88, IRAK4, TRAF6, and IкB-α, were significantly upregulated (by 22.9-, 9.6-, 4.0-, 3.6-, and 3.9-fold, respectively) at 6 h. The cytokine-related gene IL-17 exhibited an increase of 6.3-fold and 10.5-fold at 3 h and 6 h, respectively. The apoptosis-related gene IAP1 was dramatically increased by 2.99-fold at 6 h. These results indicate that adult T. crocea could initiate the TLR pathway to resist V. coralliilyticus, which promotes the release of inflammatory factors such as IL-17 and leads to the activation of a series of outcomes, such as apoptosis. The response mechanism is related to the T. crocea immunoreaction stimulated by V. coralliilyticus, providing a theoretical basis for understanding T. crocea immune response mechanisms., Competing Interests: The authors have declared that no competing interests exist.

Ruditapes philippinarum conglutination mud (RPM) is a byproduct from the aquiculture of an important commercially bivalve mollusk R. philippinarum and has been recently reported as a promising natural bioflocculant resource. However the origin of bioflocculation components within RPM is still a pending doubt and impedes its effective exploitation. This study investigated the probability that RPM bioflocculation components originate from its associated microbes. RPM samples from an aquaculture farm in Zhoushan of China were applied to characterize its microbial community structure, screen associated bioflocculant-producing strains, and explore the homology between extracellular polysaccharides (EPS) from bioflocculant-producing isolates and RPM flocculation components. Results showed that RPM exhibited high bacterial biodiversity, with Proteobacteria, Bacteroidetes and Actinobacteria as the most abundant phyla; hgcI_clade, CL500_29_marine_group, Fusibacter, MWH_UniP1_aquatic_group and Arcobacter as the dominant genera. Fourteen highly efficient bioflocculant-producing strains were screened and phylogenetically identified as Pseudoalteromonas sp. (5), Psychrobacter sp. (3), Halomonas sp. (2), Albirhodobacter sp. (1), Celeribacter sp. (1), Kocuria sp. (1) and Bacillus sp. (1), all of which except Bacillus sp. were reported for the first time for their excellent flocculation capability. Furthermore, EPS from the bioflocculant-producing strains exhibited highly similar monosaccharide composition to the reported flocculation-effective RPM polysaccharides. On the other hand, the existence of fungi in RPM was rare and showed no flocculation functionality. Findings from Zhoushan RPM strongly supported that RPM flocculation components were of bacterial origin and make RPM reproduction possible by fermentation approach., Competing Interests: The authors have declared that no competing interests exist.

Teredinidae are a family of highly adapted wood-feeding and wood-boring bivalves, commonly known as shipworms, whose evolution is linked to the acquisition of cellulolytic gammaproteobacterial symbionts harbored in bacteriocytes within the gills. In the present work we applied metagenomics to characterize microbiomes of the gills and digestive tract of Neoteredo reynei, a mangrove-adapted shipworm species found over a large range of the Brazilian coast. Comparative metagenomics grouped the gill symbiont community of different N. reynei specimens, indicating closely related bacterial types are shared. Similarly, the intestine and digestive gland communities were related, yet were more diverse than and showed no overlap with the gill community. Annotation of assembled metagenomic contigs revealed that the gill symbiotic community of N. reynei encodes a plethora of plant cell wall polysaccharides degrading glycoside hydrolase encoding genes, and Biosynthetic Gene Clusters (BGCs). In contrast, the digestive tract microbiomes seem to play little role in wood digestion and secondary metabolites biosynthesis. Metagenome binning recovered the nearly complete genome sequences of two symbiotic Teredinibacter strains from the gills, a representative of Teredinibacter turnerae "clade I" strain, and a yet to be cultivated Teredinibacter sp. type. These Teredinibacter genomes, as well as un-binned gill-derived gammaproteobacteria contigs, also include an endo-β-1,4-xylanase/acetylxylan esterase multi-catalytic carbohydrate-active enzyme, and a trans-acyltransferase polyketide synthase (trans-AT PKS) gene cluster with the gene cassette for generating β-branching on complex polyketides. Finally, we use multivariate analyses to show that the secondary metabolome from the genomes of Teredinibacter representatives, including genomes binned from N. reynei gills' metagenomes presented herein, stands out within the Cellvibrionaceae family by size, and enrichments for polyketide, nonribosomal peptide and hybrid BGCs. Results presented here add to the growing characterization of shipworm symbiotic microbiomes and indicate that the N. reynei gill gammaproteobacterial community is a prolific source of biotechnologically relevant enzymes for wood-digestion and bioactive compounds production., Competing Interests: The authors have declared that no competing interests exist.

Bacterial infectious diseases produced by Vibrio are the main cause of economic losses in aquaculture. During recent years it has been shown that the expression of virulence genes in some Vibrio species is controlled by a population-density dependent gene-expression mechanism known as quorum sensing (QS), which is mediated by the diffusion of signal molecules such as N-acylhomoserine lactones (AHLs). QS disruption, especially the enzymatic degradation of signalling molecules, known as quorum quenching (QQ), is one of the novel therapeutic strategies for the treatment of bacterial infections. In this study, we present the detection of AHLs in 34 marine Vibrionaceae strains. Three aquaculture-related pathogenic Vibrio strains, V. mediterranei VibC-Oc-097, V. owensii VibC-Oc-106 and V. coralliilyticus VibC-Oc-193 were selected for further studies based on their virulence and high production of AHLs. This is the first report where the signal molecules have been characterized in these emerging marine pathogens and correlated to the expression of virulence factors. Moreover, the results of AHL inactivation in the three selected strains have been confirmed in vivo against brine shrimps (Artemia salina) and Manila clams (Venerupis philippinarum). This research contributes to the development of future therapies based on AHL disruption, the most promising alternatives for fighting infectious diseases in aquaculture.

A novel virulent bacteriophage, vB_VspP_pVa5, infecting a strain of Vibrio splendidus was isolated from a sea-cage aquaculture farm in Greece, and characterized using microbiological methods and genomic analysis. Bacteriophage vB_VspP_pVa5 is a N4-like podovirus with an icosahedral head measuring 85 nm in length and a short non-contractile tail. The phage had a narrow host range infecting only the bacterial host, a latent period of 30 min and a burst size of 24 virions per infected bacterium. Its genome size was 78,145 bp and genomic analysis identified 107 densely-packed genes, 40 of which could be annotated. In addition to the very large virion encapsulated DNA-dependent RNA polymerase which is the signature of the N4-like genus, an interesting feature of the novel phage is the presence of a self-splicing group I intron in the thymidylate synthase gene. A tRNAStop interrupted by a ~2.5kb open reading frame-containing area was also identified. The absence of genes related to lysogeny along with the high efficacy observed during in vitro cell lysis trials, indicate that the vB_VspP_pVa5 is a potential candidate component in a bacteriophage cocktail suitable for the biological control of V. splendidus in aquaculture.

The aim of this study was first to identify lysozymes paralogs in the deep sea mussel Bathymodiolus azoricus then to measure their relative expression or activity in different tissue or conditions. B. azoricus is a bivalve that lives close to hydrothermal chimney in the Mid-Atlantic Ridge (MAR). They harbour in specialized gill cells two types of endosymbiont (gram-bacteria): sulphide oxidizing bacteria (SOX) and methanotrophic bacteria (MOX). This association is thought to be ruled by specific mechanism or actors of regulation to deal with the presence of symbiont but these mechanisms are still poorly understood. Here, we focused on the implication of lysozyme, a bactericidal enzyme, in this endosymbiosis. The relative expression of Ba-lysozymes paralogs and the global anti-microbial activity, were measured in natural population (Lucky Strike--1700 m, Mid-Atlantic Ridge), and in in situ experimental conditions. B. azoricus individuals were moved away from the hydrothermal fluid to induce a loss of symbiont. Then after 6 days some mussels were brought back to the mussel bed to induce a re-acquisition of symbiotic bacteria. Results show the presence of 6 paralogs in B. azoricus. In absence of symbionts, 3 paralogs are up-regulated while others are not differentially expressed. Moreover the global activity of lysozyme is increasing with the loss of symbiont. All together these results suggest that lysozyme may play a crucial role in symbiont regulation.

The cornerstones of sunken wood ecosystems are microorganisms involved in cellulose degradation. These can either be free-living microorganisms in the wood matrix or symbiotic bacteria associated with wood-boring bivalves such as emblematic species of Xylophaga, the most common deep-sea woodborer. Here we use experimentally submerged pine wood, placed in and outside the Mediterranean submarine Blanes Canyon, to compare the microbial communities on the wood, in fecal pellets of Xylophaga spp. and associated with the gills of these animals. Analyses based on tag pyrosequencing of the 16S rRNA bacterial gene showed that sunken wood contained three distinct microbial communities. Wood and pellet communities were different from each other suggesting that Xylophaga spp. create new microbial niches by excreting fecal pellets into their burrows. In turn, gills of Xylophaga spp. contain potential bacterial symbionts, as illustrated by the presence of sequences closely related to symbiotic bacteria found in other wood eating marine invertebrates. Finally, we found that sunken wood communities inside the canyon were different and more diverse than the ones outside the canyon. This finding extends to the microbial world the view that submarine canyons are sites of diverse marine life.

Shipworms are marine bivalve mollusks (Family Teredinidae) that use wood for shelter and food. They harbor a group of closely related, yet phylogenetically distinct, bacterial endosymbionts in bacteriocytes located in the gills. This endosymbiotic community is believed to support the host's nutrition in multiple ways, through the production of cellulolytic enzymes and the fixation of nitrogen. The genome of the shipworm endosymbiont Teredinibacter turnerae T7901 was recently sequenced and in addition to the potential for cellulolytic enzymes and diazotrophy, the genome also revealed a rich potential for secondary metabolites. With nine distinct biosynthetic gene clusters, nearly 7% of the genome is dedicated to secondary metabolites. Bioinformatic analyses predict that one of the gene clusters is responsible for the production of a catecholate siderophore. Here we describe this gene cluster in detail and present the siderophore product from this cluster. Genes similar to the entCEBA genes of enterobactin biosynthesis involved in the production and activation of dihydroxybenzoic acid (DHB) are present in this cluster, as well as a two-module non-ribosomal peptide synthetase (NRPS). A novel triscatecholate siderophore, turnerbactin, was isolated from the supernatant of iron-limited T. turnerae T7901 cultures. Turnerbactin is a trimer of N-(2,3-DHB)-L-Orn-L-Ser with the three monomeric units linked by Ser ester linkages. A monomer, dimer, dehydrated dimer, and dehydrated trimer of 2,3-DHB-L-Orn-L-Ser were also found in the supernatant. A link between the gene cluster and siderophore product was made by constructing a NRPS mutant, TtAH03. Siderophores could not be detected in cultures of TtAH03 by HPLC analysis and Fe-binding activity of culture supernatant was significantly reduced. Regulation of the pathway by iron is supported by identification of putative Fur box sequences and observation of increased Fe-binding activity under iron restriction. Evidence of a turnerbactin fragment was found in shipworm extracts, suggesting the production of turnerbactin in the symbiosis.

Vesicomyid clams harbor intracellular sulfur-oxidizing bacteria that are predominantly maternally inherited and co-speciate with their hosts. Genome recombination and the occurrence of non-parental strains were recently demonstrated in symbionts. However, mechanisms favoring such events remain to be identified. In this study, we investigated symbionts in two phylogenetically distant vesicomyid species, Christineconcha regab and Laubiericoncha chuni, which sometimes co-occur at a cold-seep site in the Gulf of Guinea. We showed that each of the two species harbored a single dominant bacterial symbiont strain. However, for both vesicomyid species, the symbiont from the other species was occasionally detected in the gills using fluorescence in situ hybridization and gene sequences analyses based on six symbiont marker genes. Symbiont strains co-occurred within a single host only at sites where both host species were found; whereas one single symbiont strain was detected in C. regab specimens from a site where no L. chuni individuals had been observed. These results suggest that physical proximity favored the acquisition of non-parental symbiont strains in Vesicomyidae. Over evolutionary time, this could potentially lead to genetic exchanges among symbiont species and eventually symbiont displacement. Symbiont densities estimated using 3D fluorescence in situ hybridization varied among host species and sites, suggesting flexibility in the association despite the fact that a similar type of metabolism is expected in all symbionts.

Transform faults are geological structures that interrupt the continuity of mid-ocean ridges and can act as dispersal barriers for hydrothermal vent organisms. In the equatorial Atlantic Ocean, it has been hypothesized that long transform faults impede gene flow between the northern and the southern Mid-Atlantic Ridge (MAR) and disconnect a northern from a southern biogeographic province. To test if there is a barrier effect in the equatorial Atlantic, we examined phylogenetic relationships of chemosynthetic bivalves and their bacterial symbionts from the recently discovered southern MAR hydrothermal vents at 5°S and 9°S. We examined Bathymodiolus spp. mussels and Abyssogena southwardae clams using the mitochondrial cytochrome c oxidase subunit I (COI) gene as a phylogenetic marker for the hosts and the bacterial 16S rRNA gene as a marker for the symbionts. Bathymodiolus spp. from the two southern sites were genetically divergent from the northern MAR species B. azoricus and B. puteoserpentis but all four host lineages form a monophyletic group indicating that they radiated after divergence from their northern Atlantic sister group, the B. boomerang species complex. This suggests dispersal of Bathymodiolus species from north to south across the equatorial belt. 16S rRNA genealogies of chemoautotrophic and methanotrophic symbionts of Bathymodiolus spp. were inconsistent and did not match the host COI genealogy indicating disconnected biogeography patterns. The vesicomyid clam Abyssogena southwardae from 5°S shared an identical COI haplotype with A. southwardae from the Logatchev vent field on the northern MAR and their symbionts shared identical 16S phylotypes, suggesting gene flow across the Equator. Our results indicate genetic connectivity between the northern and southern MAR and suggest that a strict dispersal barrier does not exist.

Marine bivalves of the family Teredinidae (shipworms) are voracious consumers of wood in marine environments. In several shipworm species, dense communities of intracellular bacterial endosymbionts have been observed within specialized cells (bacteriocytes) of the gills (ctenidia). These bacteria are proposed to contribute to digestion of wood by the host. While the microbes of shipworm gills have been studied extensively in several species, the abundance and distribution of microbes in the digestive system have not been adequately addressed. Here we use Fluorescence In-Situ Hybridization (FISH) and laser scanning confocal microscopy with 16S rRNA directed oligonucleotide probes targeting all domains, domains Bacteria and Archaea, and other taxonomic groups to examine the digestive microbiota of 17 specimens from 5 shipworm species (Bankia setacea, Lyrodus pedicellatus, Lyrodus massa, Lyrodus sp. and Teredo aff. triangularis). These data reveal that the caecum, a large sac-like appendage of the stomach that typically contains large quantities of wood particles and is considered the primary site of wood digestion, harbors only very sparse microbial populations. However, a significant number of bacterial cells were observed in fecal pellets within the intestines. These results suggest that due to low abundance, bacteria in the caecum may contribute little to lignocellulose degradation. In contrast, the comparatively high population density of bacteria in the intestine suggests a possible role for intestinal bacteria in the degradation of lignocellulose.

The marine bacterium Vibrio tapetis is the causative agent of the brown ring disease, which affects the clam Ruditapes philippinarum and causes heavy economic losses in North of Europe and in Eastern Asia. Further characterization of V. tapetis isolates showed that all the investigated strains harbored at least one large plasmid. We determined the sequence of the 82,266 bp plasmid pVT1 from the CECT4600(T) reference strain and analyzed its genetic content. pVT1 is a mosaic plasmid closely related to several conjugative plasmids isolated from Vibrio vulnificus strains and was shown to be itself conjugative in Vibrios. In addition, it contains DNA regions that have similarity with several other plasmids from marine bacteria (Vibrio sp., Shewanella sp., Listonella anguillarum and Photobacterium profundum). pVT1 contains a number of mobile elements, including twelve Insertion Sequences or inactivated IS genes and an RS1 phage element related to the CTXphi phage of V. cholerae. The genetic organization of pVT1 underscores an important role of horizontal gene transfer through conjugative plasmid shuffling and transposition events in the acquisition of new genetic resources and in generating the pVT1 modular organization. In addition, pVT1 presents a copy number of 9, relatively high for a conjugative plasmid, and appears to belong to a new type of replicon, which may be specific to Vibrionaceae and Shewanelleacae.

Background: Deep-sea mussels harboring chemoautotrophic symbionts from hydrothermal vents and seeps are assumed to have evolved from shallow-water asymbiotic relatives by way of biogenic reducing environments such as sunken wood and whale falls. Such symbiotic associations have been well characterized in mussels collected from vents, seeps and sunken wood but in only a few from whale falls., Methodology/principal Finding: Here we report symbioses in the gill tissues of two mussels, Adipicola crypta and Adipicola pacifica, collected from whale-falls on the continental shelf in the northwestern Pacific. The molecular, morphological and stable isotopic characteristics of bacterial symbionts were analyzed. A single phylotype of thioautotrophic bacteria was found in A. crypta gill tissue and two distinct phylotypes of bacteria (referred to as Symbiont A and Symbiont C) in A. pacifica. Symbiont A and the A. crypta symbiont were affiliated with thioautotrophic symbionts of bathymodiolin mussels from deep-sea reducing environments, while Symbiont C was closely related to free-living heterotrophic bacteria. The symbionts in A. crypta were intracellular within epithelial cells of the apical region of the gills and were extracellular in A. pacifica. No spatial partitioning was observed between the two phylotypes in A. pacifica in fluorescence in situ hybridization experiments. Stable isotopic analyses of carbon and sulfur indicated the chemoautotrophic nature of A. crypta and mixotrophic nature of A. pacifica. Molecular phylogenetic analyses of the host mussels showed that A. crypta constituted a monophyletic clade with other intracellular symbiotic (endosymbiotic) mussels and that A. pacifica was the sister group of all endosymbiotic mussels., Conclusions/significance: These results strongly suggest that the symbiosis in A. pacifica is at an earlier stage in evolution than other endosymbiotic mussels. Whale falls and other modern biogenic reducing environments may act as refugia for primal chemoautotrophic symbioses between eukaryotes and prokaryotes since the extinction of ancient large marine vertebrates.

Mussels are conspicuous and often abundant members of rocky shores and may constitute an important site for the nitrogen cycle due to their feeding and excretion activities. We used shotgun metagenomics of the microbial community associated with the surface of mussels (Mytilus californianus) on Tatoosh Island in Washington state to test whether there is a nitrogen-based microbial assemblage associated with mussels. Analyses of both tidepool mussels and those on emergent benches revealed a diverse community of Bacteria and Archaea with approximately 31 million bp from 6 mussels in each habitat. Using MG-RAST, between 22.5-25.6% were identifiable using the SEED non-redundant database for proteins. Of those fragments that were identifiable through MG-RAST, the composition was dominated by Cyanobacteria and Alpha- and Gamma-proteobacteria. Microbial composition was highly similar between the tidepool and emergent bench mussels, suggesting similar functions across these different microhabitats. One percent of the proteins identified in each sample were related to nitrogen cycling. When normalized to protein discovery rate, the high diversity and abundance of enzymes related to the nitrogen cycle in mussel-associated microbes is as great or greater than that described for other marine metagenomes. In some instances, the nitrogen-utilizing profile of this assemblage was more concordant with soil metagenomes in the Midwestern U.S. than for open ocean system. Carbon fixation and Calvin cycle enzymes further represented 0.65 and 1.26% of all proteins and their abundance was comparable to a number of open ocean marine metagenomes. In sum, the diversity and abundance of nitrogen and carbon cycle related enzymes in the microbes occupying the shells of Mytilus californianus suggest these mussels provide a node for microbial populations and thus biogeochemical processes.

Here we report the complete genome sequence of Teredinibacter turnerae T7901. T. turnerae is a marine gamma proteobacterium that occurs as an intracellular endosymbiont in the gills of wood-boring marine bivalves of the family Teredinidae (shipworms). This species is the sole cultivated member of an endosymbiotic consortium thought to provide the host with enzymes, including cellulases and nitrogenase, critical for digestion of wood and supplementation of the host's nitrogen-deficient diet. T. turnerae is closely related to the free-living marine polysaccharide degrading bacterium Saccharophagus degradans str. 2-40 and to as yet uncultivated endosymbionts with which it coexists in shipworm cells. Like S. degradans, the T. turnerae genome encodes a large number of enzymes predicted to be involved in complex polysaccharide degradation (>100). However, unlike S. degradans, which degrades a broad spectrum (>10 classes) of complex plant, fungal and algal polysaccharides, T. turnerae primarily encodes enzymes associated with deconstruction of terrestrial woody plant material. Also unlike S. degradans and many other eubacteria, T. turnerae dedicates a large proportion of its genome to genes predicted to function in secondary metabolism. Despite its intracellular niche, the T. turnerae genome lacks many features associated with obligate intracellular existence (e.g. reduced genome size, reduced %G+C, loss of genes of core metabolism) and displays evidence of adaptations common to free-living bacteria (e.g. defense against bacteriophage infection). These results suggest that T. turnerae is likely a facultative intracellular ensosymbiont whose niche presently includes, or recently included, free-living existence. As such, the T. turnerae genome provides insights into the range of genomic adaptations associated with intracellular endosymbiosis as well as enzymatic mechanisms relevant to the recycling of plant materials in marine environments and the production of cellulose-derived biofuels.