Genome-Wide Identification of the Sulfate Transporters Gene Family in Blueberry (Vaccinium spp.) and Its Response to Ericoid Mycorrhizal Fungi.

Sulfur metabolism plays a major role in plant growth and development, environmental adaptation, and material synthesis, and the sulfate transporters are the beginning of sulfur metabolism. We identified 37 potential VcSULTR genes in the blueberry genome, encoding peptides with 534 to 766 amino acids. The genes were grouped into four subfamilies in an evolutionary analysis. The 37 putative VcSULTR proteins ranged in size from 60.03 to 83.87 kDa. These proteins were predicted to be hydrophobic and mostly localize to the plasma membrane. The VcSULTR genes were distributed on 30 chromosomes; VcSULTR3;5b and VcSULTR3;5c were the only tandemly repeated genes. The VcSULTR promoters contained cis-acting elements related to the fungal symbiosis and stress responses. The transcript levels of the VcSULTRs differed among blueberry organs and changed in response to ericoid mycorrhizal fungi and sulfate treatments. A subcellular localization analysis showed that VcSULTR2;1c localized to, and functioned in, the plasma membrane and chloroplast. The virus-induced gene knock-down of VcSULTR2;1c resulted in a significantly decreased endogenous sulfate content, and an up-regulation of genes encoding key enzymes in sulfur metabolism (VcATPS2 and VcSiR1). These findings enhance our understanding of mycorrhizal-fungi-mediated sulfate transport in blueberry, and lay the foundation for further research on blueberry-mycorrhizal symbiosis.

Microbial diversity and ecological complexity emerging from environmental variation and horizontal gene transfer in a simple mathematical model.

BACKGROUND: Microbiomes are generally characterized by high diversity of coexisting microbial species and strains, and microbiome composition typically remains stable across a broad range of conditions. However, under fixed conditions, microbial ecology conforms with the exclusion principle under which two populations competing for the same resource within the same niche cannot coexist because the less fit population inevitably goes extinct. Therefore, the long-term persistence of microbiome diversity calls for an explanation. RESULTS: To explore the conditions for stabilization of microbial diversity, we developed a simple mathematical model consisting of two competing populations that could exchange a single gene allele via horizontal gene transfer (HGT). We found that, although in a fixed environment, with unbiased HGT, the system obeyed the exclusion principle, in an oscillating environment, within large regions of the phase space bounded by the rates of reproduction and HGT, the two populations coexist. Moreover, depending on the parameter combination, all three major types of symbiosis were obtained, namely, pure competition, host-parasite relationship, and mutualism. In each of these regimes, certain parameter combinations provided for synergy, that is, a greater total abundance of both populations compared to the abundance of the winning population in the fixed environment. CONCLUSIONS: The results of this modeling study show that basic phenomena that are universal in microbial communities, namely, environmental variation and HGT, provide for stabilization and persistence of microbial diversity, and emergence of ecological complexity.

Transcriptomic responses of Mediterranean sponges upon encounter with symbiont microbial consortia.

BACKGROUND: Sponges (phylum Porifera) constantly interact with microbes. They graze on microbes from the water column by filter-feeding and they harbor symbiotic partners within their bodies. In experimental setups, sponges take up symbionts at lower rates compared with seawater microbes. This suggests that sponges have the capacity to differentiate between microbes and preferentially graze in non-symbiotic microbes, although the underlying mechanisms of discrimination are still poorly understood. Genomic studies showed that, compared to other animal groups, sponges present an extended repertoire of immune receptors, in particular NLRs, SRCRs, and GPCRs, and a handful of experiments showed that sponges regulate the expression of these receptors upon encounter with microbial elicitors. We hypothesize that sponges may rely on differential expression of their diverse repertoire of poriferan immune receptors to sense different microbial consortia while filter-feeding. To test this, we characterized the transcriptomic response of two sponge species, Aplysina aerophoba and Dysidea avara, upon incubation with microbial consortia extracted from A. aerophoba in comparison with incubation with seawater microbes. The sponges were sampled after 1 h, 3 h, and 5 h for RNA-Seq differential gene expression analysis. RESULTS: D. avara incubated with A. aerophoba-symbionts regulated the expression of genes related to immunity, ubiquitination, and signaling. Within the set of differentially-expressed immune genes we identified different families of Nucleotide Oligomerization Domain (NOD)-Like Receptors (NLRs). These results represent the first experimental evidence that different types of NLRs are involved in microbial discrimination in a sponge. In contrast, the transcriptomic response of A. aerophoba to its own symbionts involved comparatively fewer genes and lacked genes encoding for immune receptors. CONCLUSION: Our work suggests that: (i) the transcriptomic response of sponges upon microbial exposure may imply "fine-tuning" of baseline gene expression as a result of their interaction with microbes, (ii) the differential response of sponges to microbial encounters varied between the species, probably due to species-specific characteristics or related to host's traits, and (iii) immune receptors belonging to different families of NLR-like genes played a role in the differential response to microbes, whether symbionts or food bacteria. The regulation of these receptors in sponges provides further evidence of the potential role of NLRs in invertebrate host-microbe interactions. The study of sponge responses to microbes exemplifies how investigating different animal groups broadens our knowledge of the evolution of immune specificity and symbiosis.

The receptor-like kinase ARK controls symbiotic balance across land plants.

The mutualistic arbuscular mycorrhizal (AM) symbiosis arose in land plants more than 450 million years ago and is still widely found in all major land plant lineages. Despite its broad taxonomic distribution, little is known about the molecular components underpinning symbiosis outside of flowering plants. The ARBUSCULAR RECEPTOR-LIKE KINASE (ARK) is required for sustaining AM symbiosis in distantly related angiosperms. Here, we demonstrate that ARK has an equivalent role in symbiosis maintenance in the bryophyte Marchantia paleacea and is part of a broad AM genetic program conserved among land plants. In addition, our comparative transcriptome analysis identified evolutionarily conserved expression patterns for several genes in the core symbiotic program required for presymbiotic signaling, intracellular colonization, and nutrient exchange. This study provides insights into the molecular pathways that consistently associate with AM symbiosis across land plants and identifies an ancestral role for ARK in governing symbiotic balance.

Domestication and evolutionary histories of specialized gut symbionts across cephalotine ants.

The evolution of animals and their gut symbionts is a complex phenomenon, obscured by lability and diversity. In social organisms, transmission of symbionts among relatives may yield systems with more stable associations. Here, we study the history of a social insect symbiosis involving cephalotine ants and their extracellular gut bacteria, which come predominantly from host-specialized lineages. We perform multi-locus phylogenetics for symbionts from nine bacterial orders, and map prior amplicon sequence data to lineage-assigned symbiont genomes, studying distributions of rigorously defined symbionts across 20 host species. Based on monophyly and additional hypothesis testing, we estimate that these specialized gut bacteria belong to 18 distinct lineages, of which 15 have been successfully isolated and cultured. Several symbiont lineages showed evidence for domestication events that occurred later in cephalotine evolutionary history, and only one lineage was ubiquitously detected in all 20 host species and 48 colonies sampled with amplicon 16S rRNA sequencing. We found evidence for phylogenetically constrained distributions in four symbionts, suggesting historical or genetic impacts on community composition. Two lineages showed evidence for frequent intra-lineage co-infections, highlighting the potential for niche divergence after initial domestication. Nearly all symbionts showed evidence for occasional host switching, but four may, more often, co-diversify with their hosts. Through our further assessment of symbiont localization and genomic functional profiles, we demonstrate distinct niches for symbionts with shared evolutionary histories, prompting further questions on the forces underlying the evolution of hosts and their gut microbiomes.

Whole-genome duplication in an algal symbiont bolsters coral heat tolerance.

The algal endosymbiont Durusdinium trenchii enhances the resilience of coral reefs under thermal stress. D. trenchii can live freely or in endosymbiosis, and the analysis of genetic markers suggests that this species has undergone whole-genome duplication (WGD). However, the evolutionary mechanisms that underpin the thermotolerance of this species are largely unknown. Here, we present genome assemblies for two D. trenchii isolates, confirm WGD in these taxa, and examine how selection has shaped the duplicated genome regions using gene expression data. We assess how the free-living versus endosymbiotic lifestyles have contributed to the retention and divergence of duplicated genes, and how these processes have enhanced the thermotolerance of D. trenchii. Our combined results suggest that lifestyle is the driver of post-WGD evolution in D. trenchii, with the free-living phase being the most important, followed by endosymbiosis. Adaptations to both lifestyles likely enabled D. trenchii to provide enhanced thermal stress protection to the host coral.

Comparative genomics of the primary endosymbiont Buchnera aphidicola in aphid hosts and their coevolutionary relationships.

BACKGROUND: Coevolution between modern aphids and their primary obligate, bacterial endosymbiont, Buchnera aphidicola, has been previously reported at different classification levels based on molecular phylogenetic analyses. However, the Buchnera genome remains poorly understood within the Rhus gall aphids. RESULTS: We assembled the complete genome of the endosymbiont Buchnera in 16 aphid samples, representing 13 species in all six genera of Rhus gall aphids by shotgun genome skimming method. We compared the newly assembled genomes with those from GenBank to comprehensively investigate patterns of coevolution between the bacteria Buchnera and their aphid hosts. Buchnera genomes were mostly collinear, and the pan-genome contained 684 genes, in which the core genome contained 256 genes with some lineages having large numbers of tandem gene duplications. There has been substantial gene-loss in each Buchnera lineage. We also reconstructed the phylogeny for Buchnera and their host aphids, respectively, using 72 complete genomes of Buchnera, along with the complete mitochondrial genomes and three nuclear genes of 31 corresponding host aphid accessions. The cophylogenetic test demonstrated significant coevolution between these two partner groups at individual, species, generic, and tribal levels. CONCLUSIONS: Buchnera exhibits very high levels of genomic sequence divergence but relative stability in gene order. The relationship between the symbionts Buchnera and its aphid hosts shows a significant coevolutionary pattern and supports complexity of the obligate symbiotic relationship.

Convergent reductive evolution of cyanobacteria in symbiosis with Dinophysiales dinoflagellates.

The diversity of marine cyanobacteria has been extensively studied due to their vital roles in ocean primary production. However, little is understood about the diversity of cyanobacterial species involved in symbiotic relationships. In this study, we successfully sequenced the complete genome of a cyanobacterium in symbiosis with Citharistes regius, a dinoflagellate species thriving in the open ocean. A phylogenomic analysis revealed that the cyanobacterium (CregCyn) belongs to the marine picocyanobacterial lineage, akin to another cyanobacterial symbiont (OmCyn) of a different dinoflagellate closely related to Citharistes. Nevertheless, these two symbionts are representing distinct lineages, suggesting independent origins of their symbiotic lifestyles. Despite the distinct origins, the genome analyses of CregCyn revealed shared characteristics with OmCyn, including an obligate symbiotic relationship with the host dinoflagellates and a degree of genome reduction. In contrast, a detailed analysis of genome subregions unveiled that the CregCyn genome carries genomic islands that are not found in the OmCyn genome. The presence of the genomic islands implies that exogenous genes have been integrated into the CregCyn genome at some point in its evolution. This study contributes to our understanding of the complex history of the symbiosis between dinoflagellates and cyanobacteria, as well as the genomic diversity of marine picocyanobacteria.

The Nitroplast and Its Relatives Support a Universal Model of Features Predicting Gene Retention in Endosymbiont and Organelle Genomes.

Endosymbiotic relationships have shaped eukaryotic life. As endosymbionts coevolve with their host, toward full integration as organelles, their genomes tend to shrink, with genes being completely lost or transferred to the host nucleus. Modern endosymbionts and organelles show diverse patterns of gene retention, and why some genes and not others are retained in these genomes is not fully understood. Recent bioinformatic study has explored hypothesized influences on these evolutionary processes, finding that hydrophobicity and amino acid chemistry predict patterns of gene retention, both in organelles across eukaryotes and in less mature endosymbiotic relationships. The exciting ongoing elucidation of endosymbiotic relationships affords an independent set of instances to test this theory. Here, we compare the properties of retained genes in the nitroplast, recently reported to be an integrated organelle, two related cyanobacterial endosymbionts that form "spheroid bodies" in their host cells, and a range of other endosymbionts, with free-living relatives of each. We find that in each case, the symbiont's genome encodes proteins with higher hydrophobicity and lower amino pKa than their free-living relative, supporting the data-derived model predicting the retention propensity of genes across endosymbiont and organelle genomes.

Comparative genomic analysis of Bradyrhizobium strains with natural variability in the efficiency of nitrogen fixation, competitiveness, and adaptation to stressful edaphoclimatic conditions.

Bradyrhizobium is known for fixing atmospheric nitrogen in symbiosis with agronomically important crops. This study focused on two groups of strains, each containing eight natural variants of the parental strains, Bradyrhizobium japonicum SEMIA 586 (=CNPSo 17) or Bradyrhizobium diazoefficiens SEMIA 566 (=CNPSo 10). CNPSo 17 and CNPSo 10 were used as commercial inoculants for soybean crops in Brazil at the beginning of the crop expansion in the southern region in the 1960s-1970s. Variants derived from these parental strains were obtained in the late 1980s through a strain selection program aimed at identifying elite strains adapted to a new cropping frontier in the central-western Cerrado region, with a higher capacity of biological nitrogen fixation (BNF) and competitiveness. Here, we aimed to detect genetic variations possibly related to BNF, competitiveness for nodule occupancy, and adaptation to the stressful conditions of the Brazilian Cerrado soils. High-quality genome assemblies were produced for all strains. The core genome phylogeny revealed that strains of each group are closely related, as confirmed by high average nucleotide identity values. However, variants accumulated divergences resulting from horizontal gene transfer, genomic rearrangements, and nucleotide polymorphisms. The B. japonicum group presented a larger pangenome and a higher number of nucleotide polymorphisms than the B. diazoefficiens group, possibly due to its longer adaptation time to the Cerrado soil. Interestingly, five strains of the B. japonicum group carry two plasmids. The genetic variability found in both groups is discussed considering the observed differences in their BNF capacity, competitiveness for nodule occupancy, and environmental adaptation.IMPORTANCEToday, Brazil is a global leader in the study and use of biological nitrogen fixation with soybean crops. As Brazilian soils are naturally void of soybean-compatible bradyrhizobia, strain selection programs were established, starting with foreign isolates. Selection searched for adaptation to the local edaphoclimatic conditions, higher efficiency of nitrogen fixation, and strong competitiveness for nodule occupancy. We analyzed the genomes of two parental strains of Bradyrhizobium japonicum and Bradyrhizobium diazoefficiens and eight variant strains derived from each parental strain. We detected two plasmids in five strains and several genetic differences that might be related to adaptation to the stressful conditions of the soils of the Brazilian Cerrado biome. We also detected genetic variations in specific regions that may impact symbiotic nitrogen fixation. Our analysis contributes to new insights into the evolution of Bradyrhizobium, and some of the identified differences may be applied as genetic markers to assist strain selection programs.

The genomic and cellular basis of biosynthetic innovation in rove beetles.

How evolution at the cellular level potentiates macroevolutionary change is central to understanding biological diversification. The >66,000 rove beetle species (Staphylinidae) form the largest metazoan family. Combining genomic and cell type transcriptomic insights spanning the largest clade, Aleocharinae, we retrace evolution of two cell types comprising a defensive gland-a putative catalyst behind staphylinid megadiversity. We identify molecular evolutionary steps leading to benzoquinone production by one cell type via a mechanism convergent with plant toxin release systems, and synthesis by the second cell type of a solvent that weaponizes the total secretion. This cooperative system has been conserved since the Early Cretaceous as Aleocharinae radiated into tens of thousands of lineages. Reprogramming each cell type yielded biochemical novelties enabling ecological specialization-most dramatically in symbionts that infiltrate social insect colonies via host-manipulating secretions. Our findings uncover cell type evolutionary processes underlying the origin and evolvability of a beetle chemical innovation.

Ecological divergence despite common mating sites: Genotypes and symbiotypes shed light on cryptic diversity in the black bean aphid species complex.

Different host plants represent ecologically dissimilar environments for phytophagous insects. The resulting divergent selection can promote the evolution of specialized host races, provided that gene flow is reduced between populations feeding on different plants. In black bean aphids belonging to the Aphis fabae complex, several morphologically cryptic taxa have been described based on their distinct host plant preferences. However, host choice and mate choice are largely decoupled in these insects: they are host-alternating and migrate between specific summer host plants and shared winter hosts, with mating occurring on the shared hosts. This provides a yearly opportunity for gene flow among aphids using different summer hosts, and raises the question if and to what extent the ecologically defined taxa are reproductively isolated. Here, we analyzed a geographically and temporally structured dataset of microsatellite genotypes from A. fabae that were mostly collected from their main winter host Euonymus europaeus, and additionally from another winter host and fourteen summer hosts. The data reveals multiple, strongly differentiated genetic clusters, which differ in their association with different summer and winter hosts. The clusters also differ in the frequency of infection with two heritable, facultative endosymbionts, separately hinting at reproductive isolation and divergent ecological selection. Furthermore, we found evidence for occasional hybridization among genetic clusters, with putative hybrids collected more frequently in spring than in autumn. This suggests that similar to host races in other phytophagous insects, both prezygotic and postzygotic barriers including selection against hybrids maintain genetic differentiation among A. fabae taxa, despite a common mating habitat.

Of sea, rivers and symbiosis: Diversity, systematics, biogeography and evolution of the deeply diverging florideophycean order Hildenbrandiales (Rhodophyta).

The Hildenbrandiales, a typically saxicolous red algal order, is an early diverging florideophycean group with global significance in marine and freshwater ecosystems across diverse temperature zones. To comprehensively elucidate the diversity, phylogeny, biogeography, and evolution of this order, we conducted a thorough re-examination employing molecular data derived from nearly 700 specimens. Employing a species delimitation method, we identified Evolutionary Species Units (ESUs) within the Hildenbrandiales aiming to enhance our understanding of species diversity and generate the first time-calibrated tree and ancestral area reconstruction for this order. Mitochondrial cox1 and chloroplast rbcL markers were used to infer species boundaries, and subsequent phylogenetic reconstructions involved concatenated sequences of cox1, rbcL, and 18S rDNA. Time calibration of the resulting phylogenetic tree used a fossil record from a Triassic purportedly freshwater Hildenbrandia species and three secondary time points from the literature. Our species delimitation analysis revealed an astounding 97 distinct ESUs, quintupling the known diversity within this order. Our time-calibration analysis placed the origin of Hildenbrandiales (crown age) in the Ediacaran period, with freshwater species emerging as a monophyletic group during the later Permian to early Triassic. Phylogenetic reconstructions identified seven major clades, experiencing early diversification during the Silurian to Carboniferous period. Two major evolutionary events-colonization of freshwater habitats and obligate systemic symbiosis with a marine fungus-marked this order, leading to significant morphological alterations without a commensurate increase in species diversification. Despite the remarkable newly discovered diversity, the extant taxon diversity appears relatively constrained when viewed against an evolutionary timeline spanning over 800 million years. This limitation may stem from restricted geographic sampling or the prevalence of asexual reproduction. However, species richness estimation and rarefaction analyses suggest a substantially larger diversity yet to be uncovered-potentially four times greater. These findings drastically reshape our understanding of the deeply diverging florideophycean order Hildenbrandiales species diversity, and contribute valuable insights into this order's evolutionary history and ecological adaptations. Supported by phylogenetic, ecological and morphological evidence, we established the genus Riverina gen. nov. to accommodate freshwater species of Hildenbrandiales, which form a monophyletic clade in our analyses. This marks the first step toward refining the taxonomy of the Hildenbrandiales, an order demanding thorough revisions, notably with the creation of several genera to address the polyphyletic status of Hildenbrandia. However, the limited diagnostic features pose a challenge, necessitating a fresh approach to defining genera. A potential solution lies in embracing a molecular systematic perspective, which can offer precise delineations of taxonomic boundaries.

A molecular atlas of plastid and mitochondrial proteins reveals organellar remodeling during plant evolutionary transitions from algae to angiosperms.

Algae and plants carry 2 organelles of endosymbiotic origin that have been co-evolving in their host cells for more than a billion years. The biology of plastids and mitochondria can differ significantly across major lineages and organelle changes likely accompanied the adaptation to new ecological niches such as the terrestrial habitat. Based on organelle proteome data and the genomes of 168 phototrophic (Archaeplastida) versus a broad range of 518 non-phototrophic eukaryotes, we screened for changes in plastid and mitochondrial biology across 1 billion years of evolution. Taking into account 331,571 protein families (or orthogroups), we identify 31,625 protein families that are unique to primary plastid-bearing eukaryotes. The 1,906 and 825 protein families are predicted to operate in plastids and mitochondria, respectively. Tracing the evolutionary history of these protein families through evolutionary time uncovers the significant remodeling the organelles experienced from algae to land plants. The analyses of gained orthogroups identifies molecular changes of organelle biology that connect to the diversification of major lineages and facilitated major transitions from chlorophytes en route to the global greening and origin of angiosperms.

Wolbachia symbionts control sex in a parasitoid wasp using a horizontally acquired gene.

Host reproduction can be manipulated by bacterial symbionts in various ways. Parthenogenesis induction is the most effective type of reproduction manipulation by symbionts for their transmission. Insect sex is determined by regulation of doublesex (dsx) splicing through transformer2 (tra2) and transformer (tra) interaction. Although parthenogenesis induction by symbionts has been studied since the 1970s, its underlying molecular mechanism is unknown. Here we identify a Wolbachia parthenogenesis-induction feminization factor gene (piff) that targets sex-determining genes and causes female-producing parthenogenesis in the haplodiploid parasitoid Encarsia formosa. We found that Wolbachia elimination repressed expression of female-specific dsx and enhanced expression of male-specific dsx, which led to the production of wasp haploid male offspring. Furthermore, we found that E. formosa tra is truncated and non-functional, and Wolbachia has a functional tra homolog, termed piff, with an insect origin. Wolbachia PIFF can colocalize and interact with wasp TRA2. Moreover, Wolbachia piff has coordinated expression with tra2 and dsx of E. formosa. Our results demonstrate the bacterial symbiont Wolbachia has acquired an insect gene to manipulate the host sex determination cascade and induce parthenogenesis in wasps. This study reveals insect-to-bacteria horizontal gene transfer drives the evolution of animal sex determination systems, elucidating a striking mechanism of insect-microbe symbiosis.

Experimental evidence that root-associated fungi improve plant growth at high altitude.

Unravelling how species communities change along environmental gradients requires a dual understanding: the direct responses of the species to their abiotic surroundings and the indirect variation of these responses through biotic interactions. Here, we focus on the interactive relationships between plants and their symbiotic root-associated fungi (RAF) along stressful abiotic gradients. We investigate whether variations in RAF community composition along altitudinal gradients influence plant growth at high altitudes, where both plants and fungi face harsher abiotic conditions. We established a translocation experiment between pairs of Bistorta vivipara populations across altitudinal gradients. To separate the impact of shifting fungal communities from the overall influence of changing abiotic conditions, we used a root barrier to prevent new colonization by RAF following translocation. To characterize the RAF communities, we applied DNA barcoding to the root samples. Through the utilization of joint species distribution modelling, we assessed the relationship between changes in plant functional traits resulting from experimental treatments and the corresponding changes in the RAF communities. Our findings indicate that RAF communities influence plant responses to stressful abiotic conditions. Plants translocated from low to high altitudes grew more when they were able to associate with the resident high-altitude RAF compared to those plants that were not allowed to associate with the resident RAF. We conclude that interactions with RAF impact how plants respond to stressful abiotic conditions. Our results provide experimental support that interactions with RAF improve plant stress tolerance to altitudinal stressors such as colder temperatures and less nutrient availability.

Rhizobial-induced phosphatase GmPP2C61A positively regulates soybean nodulation.

Symbiotic nitrogen fixation (SNF) is crucial for legumes, providing them with the nitrogen necessary for plant growth and development. Nodulation is the first step in the establishment of SNF. However, the determinant genes in soybean nodulation and the understanding of the underlying molecular mechanisms governing nodulation are still limited. Herein, we identified a phosphatase, GmPP2C61A, which was specifically induced by rhizobia inoculation. Using transgenic hairy roots harboring GmPP2C61A::GUS, we showed that GmPP2C61A was mainly induced in epidermal cells following rhizobia inoculation. Functional analysis revealed that knockdown or knock-out of GmPP2C61A significantly reduced the number of nodules, while overexpression of GmPP2C61A promoted nodule formation. Additionally, GmPP2C61A protein was mainly localized in the cytoplasm and exhibited conserved phosphatase activity in vitro. Our findings suggest that phosphatase GmPP2C61A serves as a critical regulator in soybean nodulation, highlighting its potential significance in enhancing symbiotic nitrogen fixation.

SoxY gene family expansion underpins adaptation to diverse hosts and environments in symbiotic sulfide oxidizers.

Sulfur-oxidizing bacteria (SOB) have developed distinct ecological strategies to obtain reduced sulfur compounds for growth. These range from specialists that can only use a limited range of reduced sulfur compounds to generalists that can use many different forms as electron donors. Forming intimate symbioses with animal hosts is another highly successful ecological strategy for SOB, as animals, through their behavior and physiology, can enable access to sulfur compounds. Symbioses have evolved multiple times in a range of animal hosts and from several lineages of SOB. They have successfully colonized a wide range of habitats, from seagrass beds to hydrothermal vents, with varying availability of symbiont energy sources. Our extensive analyses of sulfur transformation pathways in 234 genomes of symbiotic and free-living SOB revealed widespread conservation in metabolic pathways for sulfur oxidation in symbionts from different host species and environments, raising the question of how they have adapted to such a wide range of distinct habitats. We discovered a gene family expansion of soxY in these genomes, with up to five distinct copies per genome. Symbionts harboring only the "canonical" soxY were typically ecological "specialists" that are associated with specific host subfamilies or environments (e.g., hydrothermal vents, mangroves). Conversely, symbionts with multiple divergent soxY genes formed versatile associations across diverse hosts in various marine environments. We hypothesize that expansion and diversification of the soxY gene family could be one genomic mechanism supporting the metabolic flexibility of symbiotic SOB enabling them and their hosts to thrive in a range of different and dynamic environments.IMPORTANCESulfur metabolism is thought to be one of the most ancient mechanisms for energy generation in microorganisms. A diverse range of microorganisms today rely on sulfur oxidation for their metabolism. They can be free-living, or they can live in symbiosis with animal hosts, where they power entire ecosystems in the absence of light, such as in the deep sea. In the millions of years since they evolved, sulfur-oxidizing bacteria have adopted several highly successful strategies; some are ecological "specialists," and some are "generalists," but which genetic features underpin these ecological strategies are not well understood. We discovered a gene family that has become expanded in those species that also seem to be "generalists," revealing that duplication, repurposing, and reshuffling existing genes can be a powerful mechanism driving ecological lifestyle shifts.

Adaptations to nitrogen availability drive ecological divergence of chemosynthetic symbionts.

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.

Stochasticity, determinism, and contingency shape genome evolution of endosymbiotic bacteria.

Evolution results from the interaction of stochastic and deterministic processes that create a web of historical contingency, shaping gene content and organismal function. To understand the scope of this interaction, we examine the relative contributions of stochasticity, determinism, and contingency in shaping gene inactivation in 34 lineages of endosymbiotic bacteria, Sodalis, found in parasitic lice, Columbicola, that are independently undergoing genome degeneration. Here we show that the process of genome degeneration in this system is largely deterministic: genes involved in amino acid biosynthesis are lost while those involved in providing B-vitamins to the host are retained. In contrast, many genes encoding redundant functions, including components of the respiratory chain and DNA repair pathways, are subject to stochastic loss, yielding historical contingencies that constrain subsequent losses. Thus, while selection results in functional convergence between symbiont lineages, stochastic mutations initiate distinct evolutionary trajectories, generating diverse gene inventories that lack the functional redundancy typically found in free-living relatives.