Tracing the evolution of gene loss in obligate bacterial symbionts

Arsenophonus symbiosis with louse flies: multiple origins, coevolutionary dynamics, and metabolic significance

IMPORTANCE

Insects that live exclusively on vertebrate blood utilize symbiotic bacteria as a source of essential compounds, e.g., B vitamins. In louse flies, the most frequent symbiont originated in genus Arsenophonus, known from a wide range of insects. Here, we analyze genomic traits, phylogenetic origins, and metabolic capacities of 11 Arsenophonus strains associated with louse flies. We show that in louse flies, Arsenophonus established symbiosis in at least four independent events, reaching different stages of symbiogenesis. This allowed for comparative genomic analysis, including convergence of metabolic capacities. The significance of the results is twofold. First, based on a comparison of independently originated Arsenophonus symbioses, it determines the importance of individual B vitamins for the insect host. This expands our theoretical insight into insect-bacteria symbiosis. The second outcome is of methodological significance. We show that the comparative approach reveals artifacts that would be difficult to identify based on a single-genome analysis.

Genomic distribution of the insertion sequence IS711 reveal a potential role in Brucella genome plasticity and host preference

The Insertion Sequence 711 (IS711) is linked to the Brucella genus. Mapping the genomic distribution of IS711 can help understand this insertion element's biological and evolutionary role. This work aimed to delineate the genomic distribution of the IS711 element and to study its association with Brucella evolution. A total of 124 genomes representing 9 Brucella species were searched using BLASTn sequence alignment tool to identify complete and truncated copies of IS711. Based on the genomic context, each IS711 locus was assigned a code using the initial letters of its neighboring genes. Various tools were used to annotate the neighboring genes and determine the shared synteny around orthologous IS711 loci. The tool Islandviewer 4 was used to scan for genomic islands. The Codon Tree method was used to build phylogenetic trees of B. melitensis, B. abortus, and B. suis genomes. The phylogenetic trees of the three species were analyzed, taking into account the genomic distribution patterns of IS711. The result of IS711 frequency analysis showed a relatively conserved number of copies/genome for the different species and for some biovars. The analysis showed that Brucella species with a relatively low IS711 copy number (4-8 copies/genome) are linked to domestic animals as primary hosts and have potential for zoonotic transmission. However, species with a relatively higher copy number (12-30 copies/genome) are less zoonotic and tend to be linked with wild animals as primary hosts. Analyzing the genomic distribution map of IS711 loci showed several unique patterns of IS711 distribution that are correlated with the evolution of Brucella species and biovars. The results also showed that 46.2% of the conserved IS711 elements are located within genomic islands. Based on our results and previous data, we postulate a model explaining the IS711 role in Brucella evolution. We assume that during the transition from a free-living to an intracellular lifestyle, a descendant of the Brucella genus had acquired a progenitor sequence of the IS711. Subsequently, a burst in IS711 transposition occurred. This parasitic expansion can be deleterious and has to be counteracted by evolutionary forces to prevent lineage extension and to promote adaptation to host. Similar to other plasmid-free pathogenic α-Proteobacteria bacteria, the balance of expansion and reduction of insertion elements could be one of the mechanisms to control genome reduction and streamlining. We hypothesize that the IS711-mediated genomic changes and other small sequence nucleotide changes in specific orthologous genes could significantly contribute to Brucella's evolution and adaptation to different animal hosts.

Supply-driven evolution: Mutation bias and trait-fitness distributions can drive macro-evolutionary dynamics

Many well-documented macro-evolutionary phenomena still challenge current evolutionary theory. Examples include long-term evolutionary trends, major transitions in evolution, conservation of certain biological features such as hox genes, and the episodic creation of new taxa. Here, we present a framework that may explain these phenomena. We do so by introducing a probabilistic relationship between trait value and reproductive fitness. This integration allows mutation bias to become a robust driver of long-term evolutionary trends against environmental bias, in a way that is consistent with all current evolutionary theories. In cases where mutation bias is strong, such as when detrimental mutations are more common than beneficial mutations, a regime called “supply-driven” evolution can arise. This regime can explain the irreversible persistence of higher structural hierarchies, which happens in the major transitions in evolution. We further generalize this result in the long-term dynamics of phenotype spaces. We show how mutations that open new phenotype spaces can become frozen in time. At the same time, new possibilities may be observed as a burst in the creation of new taxa.

Evolutionarily recent dual obligatory symbiosis among adelgids indicates a transition between fungus- and insect-associated lifestyles

Adelgids (Insecta: Hemiptera: Adelgidae) form a small group of insects but harbor a surprisingly diverse set of bacteriocyte-associated endosymbionts, which suggest multiple replacement and acquisition of symbionts over evolutionary time. Specific pairs of symbionts have been associated with adelgid lineages specialized on different secondary host conifers. Using a metagenomic approach, we investigated the symbiosis of the Adelges laricis/Adelges tardus species complex containing betaproteobacterial ("Candidatus Vallotia tarda") and gammaproteobacterial ("Candidatus Profftia tarda") symbionts. Genomic characteristics and metabolic pathway reconstructions revealed that Vallotia and Profftia are evolutionary young endosymbionts, which complement each other's role in essential amino acid production. Phylogenomic analyses and a high level of genomic synteny indicate an origin of the betaproteobacterial symbiont from endosymbionts of Rhizopus fungi. This evolutionary transition was accompanied with substantial loss of functions related to transcription regulation, secondary metabolite production, bacterial defense mechanisms, host infection, and manipulation. The transition from fungus to insect endosymbionts extends our current framework about evolutionary trajectories of host-associated microbes.

Rapid molecular evolution of Spiroplasma symbionts of Drosophila

Spiroplasma is a genus of Mollicutes whose members include plant pathogens, insect pathogens and endosymbionts of animals. Spiroplasma phenotypes have been repeatedly observed to be spontaneously lost in Drosophila cultures, and several studies have documented a high genomic turnover in Spiroplasma symbionts and plant pathogens. These observations suggest that Spiroplasma evolves quickly in comparison to other insect symbionts. Here, we systematically assess evolutionary rates and patterns of Spiroplasma poulsonii, a natural symbiont of Drosophila. We analysed genomic evolution of sHy within flies, and sMel within in vitro culture over several years. We observed that S. poulsonii substitution rates are among the highest reported for any bacteria, and around two orders of magnitude higher compared with other inherited arthropod endosymbionts. The absence of mismatch repair loci mutS and mutL is conserved across Spiroplasma, and likely contributes to elevated substitution rates. Further, the closely related strains sMel and sHy (>99.5 % sequence identity in shared loci) show extensive structural genomic differences, which potentially indicates a higher degree of host adaptation in sHy, a protective symbiont of Drosophila hydei. Finally, comparison across diverse Spiroplasma lineages confirms previous reports of dynamic evolution of toxins, and identifies loci similar to the male-killing toxin Spaid in several Spiroplasma lineages and other endosymbionts. Overall, our results highlight the peculiar nature of Spiroplasma genome evolution, which may explain unusual features of its evolutionary ecology.

Proteome size reduction in Apicomplexans is linked with loss of DNA repair and host redundant pathways

Apicomplexans are alveolate parasites which include Plasmodium falciparum, the main cause of malaria, one of the world's biggest killers from infectious disease. Apicomplexans are characterized by a reduction in proteome size, which appears to result from metabolic and functional simplification, commensurate with their parasitic lifestyle. However, other factors may also help to explain gene loss such as population bottlenecks experienced during transmission, and the effect of reducing the overall genomic information content. The latter constitutes an 'informational constraint', which is proposed to exert a selective pressure to evolve and maintain genes involved in informational fidelity and error correction, proportional to the quantity of information in the genome (which approximates to proteome size). The dynamics of gene loss was examined in 41 Apicomplexan genomes using orthogroup analysis. We show that loss of genes involved in amino acid metabolism and steroid biosynthesis can be explained by metabolic redundancy with the host. We also show that there is a marked tendency to lose DNA repair genes as proteome size is reduced. This may be explained by a reduction in size of the informational constraint and can help to explain elevated mutation rates in pathogens with reduced genome size. Multiple Sequentially Markovian Coalescent (MSMC) analysis indicates a recent bottleneck, consistent with predictions generated using allele-based population genetics approaches, implying that relaxed selection pressure due to reduced population size might have contributed to gene loss. However, the non-randomness of pathways that are lost challenges this scenario. Lastly, we identify unique orthogroups in malaria-causing Plasmodium species that infect humans, with a high proportion of membrane associated proteins. Thus, orthogroup analysis appears useful for identifying novel candidate pathogenic factors in parasites, when there is a wide sample of genomes available.

Whitefly Endosymbionts: Biology, Evolution, and Plant Virus Interactions

Whiteflies (Hemiptera: Aleyrodidae) are sap-feeding global agricultural pests. These piercing-sucking insects have coevolved with intracellular endosymbiotic bacteria that help to supplement their nutrient-poor plant sap diets with essential amino acids and carotenoids. These obligate, primary endosymbionts have been incorporated into specialized organs called bacteriomes where they sometimes coexist with facultative, secondary endosymbionts. All whitefly species harbor the primary endosymbiont Candidatus Portiera aleyrodidarum and have a variable number of secondary endosymbionts. The secondary endosymbiont complement harbored by the cryptic whitefly species Bemisia tabaci is particularly complex with various assemblages of seven different genera identified to date. In this review, we discuss whitefly associated primary and secondary endosymbionts. We focus on those associated with the notorious B. tabaci species complex with emphasis on their biological characteristics and diversity. We also discuss their interactions with phytopathogenic begomoviruses (family Geminiviridae), which are transmitted exclusively by B. tabaci in a persistent-circulative manner. Unraveling the complex interactions of these endosymbionts with their insect hosts and plant viruses could lead to advancements in whitefly and whitefly transmitted virus management.

The Role of Symbioses in the Adaptation and Stress Responses of Marine Organisms

Ocean ecosystems are experiencing unprecedented rates of climate and anthropogenic change, which can often initiate stress in marine organisms. Symbioses, or associations between different organisms, are plentiful in the ocean and could play a significant role in facilitating organismal adaptations to stressful ocean conditions. This article reviews current knowledge about the role of symbiosis in marine organismal acclimation and adaptation. It discusses stress and adaptations in symbioses from coral reef ecosystems, which are among the most affected environments in the ocean, including the relationships between corals and microalgae, corals and bacteria, anemones and clownfish, and cleaner fish and client fish. Despite the importance of this subject, knowledge of how marine organisms adapt to stress is still limited, and there are vast opportunities for research and technological development in this area. Attention to this subject will enhance our understanding of the capacity of symbioses to alleviate organismal stress in the oceans.

Comparative genomics reveals a novel genetic organization of the sad cluster in the sulfonamide-degrader 'Candidatus Leucobacter sulfamidivorax' strain GP

Background

Microbial communities recurrently establish metabolic associations resulting in increased fitness and ability to perform complex tasks, such as xenobiotic degradation. In a previous study, we have described a sulfonamide-degrading consortium consisting of a novel low-abundant actinobacterium, named strain GP, and Achromobacter denitrificans PR1. However, we found that strain GP was unable to grow independently and could not be further purified.

Results

Previous studies suggested that strain GP might represent a new putative species within the Leucobacter genus (16S rRNA gene similarity < 97%). In this study, we found that average nucleotide identity (ANI) with other Leucobacter spp. ranged between 76.8 and 82.1%, further corroborating the affiliation of strain GP to a new provisional species. The average amino acid identity (AAI) and percentage of conserved genes (POCP) values were near the lower edge of the genus delimitation thresholds (65 and 55%, respectively). Phylogenetic analysis of core genes between strain GP and Leucobacter spp. corroborated these findings. Comparative genomic analysis indicates that strain GP may have lost genes related to tetrapyrrole biosynthesis and thiol transporters, both crucial for the correct assembly of cytochromes and aerobic growth. However, supplying exogenous heme and catalase was insufficient to abolish the dependent phenotype. The actinobacterium harbors at least two copies of a novel genetic element containing a sulfonamide monooxygenase (sadA) flanked by a single IS1380 family transposase. Additionally, two homologs of sadB (4-aminophenol monooxygenase) were identified in the metagenome-assembled draft genome of strain GP, but these were not located in the vicinity of sadA nor of mobile or integrative elements.

Conclusions

Comparative genomics of the genus Leucobacter suggested the absence of some genes encoding for important metabolic traits in strain GP. Nevertheless, although media and culture conditions were tailored to supply its potential metabolic needs, these conditions were insufficient to isolate the PR1-dependent actinobacterium further. This study gives important insights regarding strain GP metabolism; however, gene expression and functional studies are necessary to characterize and further isolate strain GP. Based on our data, we propose to classify strain GP in a provisional new species within the genus Leucobacter, ‘Candidatus Leucobacter sulfamidivorax‘.

Genomic diversity, lifestyles and evolutionary origins of DPANN archaea

Archaea-a primary domain of life besides Bacteria-have for a long time been regarded as peculiar organisms that play marginal roles in biogeochemical cycles. However, this picture changed with the discovery of a large diversity of archaea in non-extreme environments enabled by the use of cultivation-independent methods. These approaches have allowed the reconstruction of genomes of uncultivated microorganisms and revealed that archaea are diverse and broadly distributed in the biosphere and seemingly include a large diversity of putative symbiotic organisms, most of which belong to the tentative archaeal superphylum referred to as DPANN. This archaeal group encompasses at least 10 different lineages and includes organisms with extremely small cell and genome sizes and limited metabolic capabilities. Therefore, many members of DPANN may be obligately dependent on symbiotic interactions with other organisms and may even include novel parasites. In this contribution, we review the current knowledge of the gene repertoires and lifestyles of members of this group and discuss their placement in the tree of life, which is the basis for our understanding of the deep microbial roots and the role of symbiosis in the evolution of life on Earth.

Biological Control of Mosquito-Borne Diseases: The Potential of Wolbachia-Based Interventions in an IVM Framework

People living in the tropical and subtropical regions of the world face an enormous health burden due to mosquito-borne diseases such as malaria, dengue fever, and filariasis. Historically and today, targeting mosquito vectors with, primarily, insecticide-based control strategies have been a key control strategy against major mosquito-borne diseases. However, the success to date of such approaches is under threat from multiple insecticide resistance mechanisms while vector control (VC) options are still limited. The situation therefore requires the development of innovative control measures against major mosquito-borne diseases. Transinfecting mosquitos with symbiotic bacteria that can compete with targeted pathogens or manipulate host biology to reduce their vectorial capacity are a promising and innovative biological control approach. In this review, we discuss the current state of knowledge about the association between mosquitoes and Wolbachia, emphasizing the limitations of different mosquito control strategies and the use of mosquitoes' commensal microbiota as innovative approaches to control mosquito-borne diseases.

The Domestication of a Large DNA Virus by the Wasp Venturia canescens Involves Targeted Genome Reduction through Pseudogenization

Polydnaviruses (PDVs) are compelling examples of viral domestication, in which wasps express a large set of genes originating from a chromosomally integrated virus to produce particles necessary for their reproductive success. Parasitoid wasps generally use PDVs as a virulence gene delivery system allowing the protection of their progeny in the body of parasitized host. However, in the wasp Venturia canescens an independent viral domestication process led to an alternative strategy as the wasp incorporates virulence proteins in viral liposomes named virus-like particles (VLPs), instead of DNA molecules. Proteomic analysis of purified VLPs and transcriptome sequencing revealed the loss of some viral functions. In particular, the genes coding for capsid components are no longer expressed, which explains why VLPs do not incorporate DNA. Here a thorough examination of V. canescens genome revealed the presence of the pseudogenes corresponding to most of the genes involved in lost functions. This strongly suggests that an accumulation of mutations that leads to gene specific pseudogenization precedes the loss of viral genes observed during virus domestication. No evidence was found for block loss of collinear genes, although extensive gene order reshuffling of the viral genome was identified from comparisons between endogenous and exogenous viruses. These results provide the first insights on the early stages of large DNA virus domestication implicating massive genome reduction through gene-specific pseudogenization, a process which differs from the large deletions described for bacterial endosymbionts.

An Unbiased Genome-Wide View of the Mutation Rate and Spectrum of the Endosymbiotic Bacterium Teredinibacter turnerae

Mutations contribute to genetic variation in all living systems. Thus, precise estimates of mutation rates and spectra across a diversity of organisms are required for a full comprehension of evolution. Here, a mutation-accumulation (MA) assay was carried out on the endosymbiotic bacterium Teredinibacter turnerae. After ∼3,025 generations, base-pair substitutions (BPSs) and insertion–deletion (indel) events were characterized by whole-genome sequencing analysis of 47 independent MA lines, yielding a BPS rate of 1.14 × 10⁻⁹ per site per generation and indel rate of 1.55 × 10⁻¹⁰ events per site per generation, which are among the highest within free-living and facultative intracellular bacteria. As in other endosymbionts, a significant bias of BPSs toward A/T and an excess of deletion mutations over insertion mutations are observed for these MA lines. However, even with a deletion bias, the genome remains relatively large (∼5.2 Mb) for an endosymbiotic bacterium. The estimate of the effective population size (N_e) in T. turnerae is quite high and comparable to free-living bacteria (∼4.5 × 10⁷), suggesting that the heavy bottlenecking associated with many endosymbiotic relationships is not prevalent during the life of this endosymbiont. The efficiency of selection scales with increasing N_e and such strong selection may have been operating against the deletion bias, preventing genome erosion. The observed mutation rate in this endosymbiont is of the same order of magnitude of those with similar N_e, consistent with the idea that population size is a primary determinant of mutation-rate evolution within endosymbionts, and that not all endosymbionts have low N_e.

Chemoenzymatic Dissection of Polyketide β-Branching in the Bryostatin Pathway

β-Branching is an expansion upon canonical polyketide synthase extension that allows for the installation of diverse chemical moieties in several natural products. Several of these moieties are unique among natural products, including the two vinyl methylesters found in the core structure of bryostatins. This family of molecules is derived from an obligate bacterial symbiont of a sessile marine bryozoan, Bugula neritina. Within this family, bryostatin 1 has been investigated as an anticancer, neuroprotective, and immunomodulatory compound. We have turned to the biosynthetic gene cluster within the bacterial symbiont to investigate the biosynthesis of bryostatins. Recent sequencing efforts resulted in the annotation of two missing genes: bryT and bryU. Using novel chemoenzymatic techniques, we have validated these as the missing enoyl-CoA hydratase and donor acyl carrier protein, essential components of the β-branching cassette of the bryostatin pathway. Together, this cassette installs the vinyl methylester moieties essential to the activity of bryostatins.

Why Might Bacterial Pathogens Have Small Genomes?

Bacteria that cause serious disease often have smaller genomes, and fewer genes, than their nonpathogenic, or less pathogenic relatives. Here, we review evidence for the generality of this association, and summarise the various reasons why the association might hold. We focus on the population genetic processes that might lead to reductive genome evolution, and show how several of these could be connected to pathogenicity. We find some evidence for most of the processes having acted in bacterial pathogens, including several different modes of genome reduction acting in the same lineage. We argue that predictable processes of genome evolution might not reflect any common underlying process.

Assessment of the bimodality in the distribution of bacterial genome sizes

Bacterial genome sizes have previously been shown to exhibit a bimodal distribution. This phenomenon has prompted discussion regarding the evolutionary forces driving genome size in bacteria and its ecological significance. We investigated the level of inherent redundancy in the public database and the effect it has on the shape of the apparent bimodal distribution. Our study reveals that there is a significant bias in the genome sequencing efforts towards a certain group of species, and that correcting the bias using species nomenclature and clustering of the 16S rRNA gene, results in a unimodal rather than the previously published bimodal distribution. The true genome size distribution and its wider ecological implications will soon emerge as we are currently witnessing rapid growth in the number of sequenced genomes from diverse environmental niches across a range of habitats at an unprecedented rate.

Gene Turnover Contributes to the Evolutionary Adaptation of Acidithiobacillus caldus: Insights from Comparative Genomics

Acidithiobacillus caldus is an extremely acidophilic sulfur-oxidizer with specialized characteristics, such as tolerance to low pH and heavy metal resistance. To gain novel insights into its genetic complexity, we chosen six A. caldus strains for comparative survey. All strains analyzed in this study differ in geographic origins as well as in ecological preferences. Based on phylogenomic analysis, we clustered the six A. caldus strains isolated from various ecological niches into two groups: group 1 strains with smaller genomes and group 2 strains with larger genomes. We found no obvious intraspecific divergence with respect to predicted genes that are related to central metabolism and stress management strategies between these two groups. Although numerous highly homogeneous genes were observed, high genetic diversity was also detected. Preliminary inspection provided a first glimpse of the potential correlation between intraspecific diversity at the genome level and environmental variation, especially geochemical conditions. Evolutionary genetic analyses further showed evidence that the difference in environmental conditions might be a crucial factor to drive the divergent evolution of A. caldus species. We identified a diverse pool of mobile genetic elements including insertion sequences and genomic islands, which suggests a high frequency of genetic exchange in these harsh habitats. Comprehensive analysis revealed that gene gains and losses were both dominant evolutionary forces that directed the genomic diversification of A. caldus species. For instance, horizontal gene transfer and gene duplication events in group 2 strains might contribute to an increase in microbial DNA content and novel functions. Moreover, genomes undergo extensive changes in group 1 strains such as removal of potential non-functional DNA, which results in the formation of compact and streamlined genomes. Taken together, the findings presented herein show highly frequent gene turnover of A. caldus species that inhabit extremely acidic environments, and shed new light on the contribution of gene turnover to the evolutionary adaptation of acidophiles.

Host-Associated Genomic Features of the Novel Uncultured Intracellular Pathogen Ca. Ichthyocystis Revealed by Direct Sequencing of Epitheliocysts

Advances in single-cell and mini-metagenome sequencing have enabled important investigations into uncultured bacteria. In this study, we applied the mini-metagenome sequencing method to assemble genome drafts of the uncultured causative agents of epitheliocystis, an emerging infectious disease in the Mediterranean aquaculture species gilthead seabream. We sequenced multiple cyst samples and constructed 11 genome drafts from a novel beta-proteobacterial lineage, Candidatus Ichthyocystis. The draft genomes demonstrate features typical of pathogenic bacteria with an obligate intracellular lifestyle: a reduced genome of up to 2.6 Mb, reduced G + C content, and reduced metabolic capacity. Reconstruction of metabolic pathways reveals that Ca Ichthyocystis genomes lack all amino acid synthesis pathways, compelling them to scavenge from the fish host. All genomes encode type II, III, and IV secretion systems, a large repertoire of predicted effectors, and a type IV pilus. These are all considered to be virulence factors, required for adherence, invasion, and host manipulation. However, no evidence of lipopolysaccharide synthesis could be found. Beyond the core functions shared within the genus, alignments showed distinction into different species, characterized by alternative large gene families. These comprise up to a third of each genome, appear to have arisen through duplication and diversification, encode many effector proteins, and are seemingly critical for virulence. Thus, Ca Ichthyocystis represents a novel obligatory intracellular pathogenic beta-proteobacterial lineage. The methods used: mini-metagenome analysis and manual annotation, have generated important insights into the lifestyle and evolution of the novel, uncultured pathogens, elucidating many putative virulence factors including an unprecedented array of novel gene families.

Conservation and Role of Electrostatics in Thymidylate Synthase

Conservation of function across families of orthologous enzymes is generally accompanied by conservation of their active site electrostatic potentials. To study the electrostatic conservation in the highly conserved essential enzyme, thymidylate synthase (TS), we conducted a systematic species-based comparison of the electrostatic potential in the vicinity of its active site. Whereas the electrostatics of the active site of TS are generally well conserved, the TSs from minimal organisms do not conform to the overall trend. Since the genomes of minimal organisms have a high thymidine content compared to other organisms, the observation of non-conserved electrostatics was surprising. Analysis of the symbiotic relationship between minimal organisms and their hosts, and the genetic completeness of the thymidine synthesis pathway suggested that TS from the minimal organism Wigglesworthia glossinidia (W.g.b.) must be active. Four residues in the vicinity of the active site of Escherichia coli TS were mutated individually and simultaneously to mimic the electrostatics of W.g.b TS. The measured activities of the E. coli TS mutants imply that conservation of electrostatics in the region of the active site is important for the activity of TS, and suggest that the W.g.b. TS has the minimal activity necessary to support replication of its reduced genome.

The Phylogeny of Rickettsia Using Different Evolutionary Signatures: How Tree-Like is Bacterial Evolution?

Rickettsia is a genus of intracellular bacteria whose hosts and transmission strategies are both impressively diverse, and this is reflected in a highly dynamic genome. Some previous studies have described the evolutionary history of Rickettsia as non-tree-like, due to incongruity between phylogenetic reconstructions using different portions of the genome. Here, we reconstruct the Rickettsia phylogeny using whole-genome data, including two new genomes from previously unsampled host groups. We find that a single topology, which is supported by multiple sources of phylogenetic signal, well describes the evolutionary history of the core genome. We do observe extensive incongruence between individual gene trees, but analyses of simulations over a single topology and interspersed partitions of sites show that this is more plausibly attributed to systematic error than to horizontal gene transfer. Some conflicting placements also result from phylogenetic analyses of accessory genome content (i.e., gene presence/absence), but we argue that these are also due to systematic error, stemming from convergent genome reduction, which cannot be accommodated by existing phylogenetic methods. Our results show that, even within a single genus, tests for gene exchange based on phylogenetic incongruence may be susceptible to false positives.

Genomic Analysis of an Ascomycete Fungus from the Rice Planthopper Reveals How It Adapts to an Endosymbiotic Lifestyle

A number of sap-sucking insects harbor endosymbionts, which are thought to play an important role in the development of their hosts. One of the most important rice pests, the brown planthopper (BPH), Nilaparvata lugens (Stål), harbors an obligatory yeast-like symbiont (YLS) that cannot be cultured in vitro. Genomic information on this YLS would be useful to better understand its evolution. In this study, we performed genome sequencing of the YLS using both 454 and Illumina approaches, generating a draft genome that shows a slightly smaller genome size and relatively higher GC content than most ascomycete fungi. A phylogenomic analysis of the YLS supported its close relationship with insect pathogens. We analyzed YLS-specific genes and the categories of genes that are likely to have changed in the YLS during its evolution. The loss of mating type locus demonstrated in the YLS sheds light on the evolution of eukaryotic symbionts. This information about the YLS genome provides a helpful guide for further understanding endosymbiotic associations in hemiptera and the symbiotic replacement of ancient bacteria with a multifunctional YLS seems to have been a successful change.

DNA Repair Is Associated with Information Content in Bacteria, Archaea, and DNA Viruses

The concept of a "proteomic constraint" proposes that DNA repair capacity is positively correlated with the information content of a genome, which can be approximated to the size of the proteome (P). This in turn implies that DNA repair genes are more likely to be present in genomes with larger values of P. This stands in contrast to the common assumption that informational genes have a core function and so are evenly distributed across organisms. We examined the presence/absence of 18 DNA repair genes in bacterial genomes. A positive relationship between gene presence and P was observed for 17 genes in the total dataset, and 16 genes when only nonintracellular bacteria were examined. A marked reduction of DNA repair genes was observed in intracellular bacteria, consistent with their reduced value of P. We also examined archaeal and DNA virus genomes, and show that the presence of DNA repair genes is likewise related to a larger value of P. In addition, the products of the bacterial genes mutY, vsr, and ndk, involved in the correction of GC/AT mutations, are strongly associated with reduced genome GC content. We therefore propose that a reduction in information content leads to a loss of DNA repair genes and indirectly to a reduction in genome GC content in bacteria by exposure to the underlying AT mutation bias. The reduction in P may also indirectly lead to the increase in substitution rates observed in intracellular bacteria via loss of DNA repair genes.

Candidatus Frankia Datiscae Dg1, the Actinobacterial Microsymbiont of Datisca glomerata, Expresses the Canonical nod Genes nodABC in Symbiosis with Its Host Plant

Frankia strains are nitrogen-fixing soil actinobacteria that can form root symbioses with actinorhizal plants. Phylogenetically, symbiotic frankiae can be divided into three clusters, and this division also corresponds to host specificity groups. The strains of cluster II which form symbioses with actinorhizal Rosales and Cucurbitales, thus displaying a broad host range, show suprisingly low genetic diversity and to date can not be cultured. The genome of the first representative of this cluster, Candidatus Frankia datiscae Dg1 (Dg1), a microsymbiont of Datisca glomerata, was recently sequenced. A phylogenetic analysis of 50 different housekeeping genes of Dg1 and three published Frankia genomes showed that cluster II is basal among the symbiotic Frankia clusters. Detailed analysis showed that nodules of D. glomerata, independent of the origin of the inoculum, contain several closely related cluster II Frankia operational taxonomic units. Actinorhizal plants and legumes both belong to the nitrogen-fixing plant clade, and bacterial signaling in both groups involves the common symbiotic pathway also used by arbuscular mycorrhizal fungi. However, so far, no molecules resembling rhizobial Nod factors could be isolated from Frankia cultures. Alone among Frankia genomes available to date, the genome of Dg1 contains the canonical nod genes nodA, nodB and nodC known from rhizobia, and these genes are arranged in two operons which are expressed in D. glomerata nodules. Furthermore, Frankia Dg1 nodC was able to partially complement a Rhizobium leguminosarum A34 nodC::Tn5 mutant. Phylogenetic analysis showed that Dg1 Nod proteins are positioned at the root of both α- and β-rhizobial NodABC proteins. NodA-like acyl transferases were found across the phylum Actinobacteria, but among Proteobacteria only in nodulators. Taken together, our evidence indicates an Actinobacterial origin of rhizobial Nod factors.

Virus world as an evolutionary network of viruses and capsidless selfish elements

Viruses were defined as one of the two principal types of organisms in the biosphere, namely, as capsid-encoding organisms in contrast to ribosome-encoding organisms, i.e., all cellular life forms. Structurally similar, apparently homologous capsids are present in a huge variety of icosahedral viruses that infect bacteria, archaea, and eukaryotes. These findings prompted the concept of the capsid as the virus "self" that defines the identity of deep, ancient viral lineages. However, several other widespread viral "hallmark genes" encode key components of the viral replication apparatus (such as polymerases and helicases) and combine with different capsid proteins, given the inherently modular character of viral evolution. Furthermore, diverse, widespread, capsidless selfish genetic elements, such as plasmids and various types of transposons, share hallmark genes with viruses. Viruses appear to have evolved from capsidless selfish elements, and vice versa, on multiple occasions during evolution. At the earliest, precellular stage of life's evolution, capsidless genetic parasites most likely emerged first and subsequently gave rise to different classes of viruses. In this review, we develop the concept of a greater virus world which forms an evolutionary network that is held together by shared conserved genes and includes both bona fide capsid-encoding viruses and different classes of capsidless replicons. Theoretical studies indicate that selfish replicons (genetic parasites) inevitably emerge in any sufficiently complex evolving ensemble of replicators. Therefore, the key signature of the greater virus world is not the presence of a capsid but rather genetic, informational parasitism itself, i.e., various degrees of reliance on the information processing systems of the host.

Symbiotic factors in Burkholderia essential for establishing an association with the bean bug, Riptortus pedestris

Symbiotic bacteria are common in insects and intimately affect the various aspects of insect host biology. In a number of insect symbiosis models, it has been possible to elucidate the effects of the symbiont on host biology, whereas there is a limited understanding of the impact of the association on the bacterial symbiont, mainly due to the difficulty of cultivating insect symbionts in vitro. Furthermore, the molecular features that determine the establishment and persistence of the symbionts in their host (i.e., symbiotic factors) have remained elusive. However, the recently established model, the bean bug Riptortus pedestris, provides a good opportunity to study bacterial symbiotic factors at a molecular level through their cultivable symbionts. Bean bugs acquire genus Burkholderia cells from the environment and harbor them as gut symbionts in the specialized posterior midgut. The genome of the Burkholderia symbiont was sequenced, and the genomic information was used to generate genetically manipulated Burkholderia symbiont strains. Using mutant symbionts, we identified several novel symbiotic factors necessary for establishing a successful association with the host gut. In this review, these symbiotic factors are classified into three categories based on the colonization dynamics of the mutant symbiont strains: initiation, accommodation, and persistence factors. In addition, the molecular characteristics of the symbiotic factors are described. These newly identified symbiotic factors and on-going studies of the Riptortus-Burkholderia symbiosis are expected to contribute to the understanding of the molecular cross-talk between insects and bacterial symbionts that are of ecological and evolutionary importance.

Erosion of functional independence early in the evolution of a microbial mutualism

Many species have evolved to function as specialized mutualists, often to the detriment of their ability to survive independently. However, there are few, if any, well-controlled observations of the evolutionary processes underlying the genesis of new mutualisms. Here, we show that within the first 1,000 generations of initiating independent syntrophic interactions between a sulfate reducer (Desulfovibrio vulgaris) and a hydrogenotrophic methanogen (Methanococcus maripaludis), D. vulgaris frequently lost the capacity to grow by sulfate respiration, thus losing the primary physiological attribute of the genus. The loss of sulfate respiration was a consequence of mutations in one or more of three key genes in the pathway for sulfate respiration, required for sulfate activation (sat) and sulfate reduction to sulfite (apsA or apsB). Because loss-of-function mutations arose rapidly and independently in replicated experiments, and because these mutations were correlated with enhanced growth rate and productivity, gene loss could be attributed to natural selection, even though these mutations should significantly restrict the independence of the evolved D. vulgaris. Together, these data present an empirical demonstration that specialization for a mutualistic interaction can evolve by natural selection shortly after its origin. They also demonstrate that a sulfate-reducing bacterium can readily evolve to become a specialized syntroph, a situation that may have often occurred in nature.

Natural selection for operons depends on genome size

In prokaryotes, genome size is associated with metabolic versatility, regulatory complexity, effective population size, and horizontal transfer rates. We therefore analyzed the covariation of genome size and operon conservation to assess the evolutionary models of operon formation and maintenance. In agreement with previous results, intraoperonic pairs of essential and of highly expressed genes are more conserved. Interestingly, intraoperonic pairs of genes are also more conserved when they encode proteins at similar cell concentrations, suggesting a role of cotranscription in diminishing the cost of waste and shortfall in gene expression. Larger genomes have fewer and smaller operons that are also less conserved. Importantly, lower conservation in larger genomes was observed for all classes of operons in terms of gene expression, essentiality, and balanced protein concentration. We reached very similar conclusions in independent analyses of three major bacterial clades (α- and β-Proteobacteria and Firmicutes). Operon conservation is inversely correlated to the abundance of transcription factors in the genome when controlled for genome size. This suggests a negative association between the complexity of genetic networks and operon conservation. These results show that genome size and/or its proxies are key determinants of the intensity of natural selection for operon organization. Our data fit better the evolutionary models based on the advantage of coregulation than those based on genetic linkage or stochastic gene expression. We suggest that larger genomes with highly complex genetic networks and many transcription factors endure weaker selection for operons than smaller genomes with fewer alternative tools for genetic regulation.

Comparative genomics reveals surprising divergence of two closely related strains of uncultivated UCYN-A cyanobacteria

Marine planktonic cyanobacteria capable of fixing molecular nitrogen (termed 'diazotrophs') are key in biogeochemical cycling, and the nitrogen fixed is one of the major external sources of nitrogen to the open ocean. Candidatus Atelocyanobacterium thalassa (UCYN-A) is a diazotrophic cyanobacterium known for its widespread geographic distribution in tropical and subtropical oligotrophic oceans, unusually reduced genome and symbiosis with a single-celled prymnesiophyte alga. Recently a novel strain of this organism was also detected in coastal waters sampled from the Scripps Institute of Oceanography pier. We analyzed the metagenome of this UCYN-A2 population by concentrating cells by flow cytometry. Phylogenomic analysis provided strong bootstrap support for the monophyly of UCYN-A (here called UCYN-A1) and UCYN-A2 within the marine Crocosphaera sp. and Cyanothece sp. clade. UCYN-A2 shares 1159 of the 1200 UCYN-A1 protein-coding genes (96.6%) with high synteny, yet the average amino-acid sequence identity between these orthologs is only 86%. UCYN-A2 lacks the same major pathways and proteins that are absent in UCYN-A1, suggesting that both strains can be grouped at the same functional and ecological level. Our results suggest that UCYN-A1 and UCYN-A2 had a common ancestor and diverged after genome reduction. These two variants may reflect adaptation of the host to different niches, which could be coastal and open ocean habitats.

Different dynamics of genome content shuffling among host-specificity groups of the symbiotic actinobacterium Frankia

BACKGROUND

Frankia is a genus of soil actinobacteria forming nitrogen-fixing root-nodule symbiotic relationships with non-leguminous woody plant species, collectively called actinorhizals, from eight dicotyledonous families. Frankia strains are classified into four host-specificity groups (HSGs), each of which exhibits a distinct host range. Genome sizes of representative strains of Alnus, Casuarina, and Elaeagnus HSGs are highly diverged and are positively correlated with the size of their host ranges.

RESULTS

The content and size of 12 Frankia genomes were investigated by in silico comparative genome hybridization and pulsed-field gel electrophoresis, respectively. Data were collected from four query strains of each HSG and compared with those of reference strains possessing completely sequenced genomes. The degree of difference in genome content between query and reference strains varied depending on HSG. Elaeagnus query strains were missing the greatest number (22-32%) of genes compared with the corresponding reference genome; Casuarina query strains lacked the fewest (0-4%), with Alnus query strains intermediate (14-18%). In spite of the remarkable gene loss, genome sizes of Alnus and Elaeagnus query strains were larger than would be expected based on total length of the absent genes. In contrast, Casuarina query strains had smaller genomes than expected.

CONCLUSIONS

The positive correlation between genome size and host range held true across all investigated strains, supporting the hypothesis that size and genome content differences are responsible for observed diversity in host plants and host plant biogeography among Frankia strains. In addition, our results suggest that different dynamics of shuffling of genome content have contributed to these symbiotic and biogeographic adaptations. Elaeagnus strains, and to a lesser extent Alnus strains, have gained and lost many genes to adapt to a wide range of environments and host plants. Conversely, rather than acquiring new genes, Casuarina strains have discarded genes to reduce genome size, suggesting an evolutionary orientation towards existence as specialist symbionts.

A single betaproteobacterium dominates the microbial community of the crambescidine-containing sponge Crambe crambe

Crambe crambe is a marine sponge that produces high concentrations of the pharmacologically significant pentacyclic guanidine alkaloids (PGAs), Crambescines and Crambescidines. Although bio-mimetic chemical synthesis of PGAs suggests involvement of microorganisms in their biosynthesis, there are conflicting reports on whether bacteria are associated with this sponge or not. Using 16S rRNA gene pyrosequencing we show that the associated bacterial community of C. crambe is dominated by a single bacterial species affiliated to the Betaproteobacteria. Microscopy analysis of sponge tissue sections using a specific probe and in situ hybridization confirmed its dominance in the sponge mesohyl and a single microbial morphology was observed by transmission electron microscopy. If confirmed the presence of a simple bacteria community in C. crambe makes this association a very pertinent model to study sponge-bacteria interactions and should allow further research into the possible implication of bacteria in PGA biosynthesis.

Amphibian skin may select for rare environmental microbes

Host-microbe symbioses rely on the successful transmission or acquisition of symbionts in each new generation. Amphibians host a diverse cutaneous microbiota, and many of these symbionts appear to be mutualistic and may limit infection by the chytrid fungus, Batrachochytrium dendrobatidis, which has caused global amphibian population declines and extinctions in recent decades. Using bar-coded 454 pyrosequencing of the 16S rRNA gene, we addressed the question of symbiont transmission by examining variation in amphibian skin microbiota across species and sites and in direct relation to environmental microbes. Although acquisition of environmental microbes occurs in some host-symbiont systems, this has not been extensively examined in free-living vertebrate-microbe symbioses. Juvenile bullfrogs (Rana catesbeiana), adult red-spotted newts (Notophthalmus viridescens), pond water and pond substrate were sampled at a single pond to examine host-specificity and potential environmental transmission of microbiota. To assess population level variation in skin microbiota, adult newts from two additional sites were also sampled. Cohabiting bullfrogs and newts had distinct microbial communities, as did newts across the three sites. The microbial communities of amphibians and the environment were distinct; there was very little overlap in the amphibians' core microbes and the most abundant environmental microbes, and the relative abundances of OTUs that were shared by amphibians and the environment were inversely related. These results suggest that, in a host species-specific manner, amphibian skin may select for microbes that are generally in low abundance in the environment.

Construction of an expressible BAC library of the unculturable insect microorganism, stink bug Plautia stali symbiont, for the search of biologically active and useful symbiont products

While gene products and metabolites of insect symbiotic bacteria may act as useful resources for insect-microbe studies and medicinal use, it is usually difficult to obtain the insect symbionts to some extent in quantity because most of them are unculturable. In this study, the possibility of using bacterial artificial chromosome (BAC) libraries as a heterologous gene expression tool for the discovery of novel symbiont metabolites was evaluated. A BAC library was constructed from the symbiont purified from the posterior midgut cecum of the stink bug Plautia stali. The BAC library, which consisted of 513 clones with an average insert size of 41 kb, represented greater than five-fold coverage of the genome. The ability of the BAC clones to express plural genes from large-sized insert DNA in Escherichia coli was examined by the growth of BAC-transformed leu operon-deficient DH10B cells on M9 minimal medium supplemented with glucose. Two BAC clones complemented leucine deficiency in DH10B cells; the clones contained the leu operon of the symbiont chromosome. The P. stali symbiont genes introduced into the BAC vector are functional in E. coli, and these genes are expressed in an operon unit. BAC libraries can be used to generate gene product- and metabolite-libraries, facilitating to characterize potential metabolites of the P. stali symbiont.

Evolution of small prokaryotic genomes

As revealed by genome sequencing, the biology of prokaryotes with reduced genomes is strikingly diverse. These include free-living prokaryotes with ∼800 genes as well as endosymbiotic bacteria with as few as ∼140 genes. Comparative genomics is revealing the evolutionary mechanisms that led to these small genomes. In the case of free-living prokaryotes, natural selection directly favored genome reduction, while in the case of endosymbiotic prokaryotes neutral processes played a more prominent role. However, new experimental data suggest that selective processes may be at operation as well for endosymbiotic prokaryotes at least during the first stages of genome reduction. Endosymbiotic prokaryotes have evolved diverse strategies for living with reduced gene sets inside a host-defined medium. These include utilization of host-encoded functions (some of them coded by genes acquired by gene transfer from the endosymbiont and/or other bacteria); metabolic complementation between co-symbionts; and forming consortiums with other bacteria within the host. Recent genome sequencing projects of intracellular mutualistic bacteria showed that previously believed universal evolutionary trends like reduced G+C content and conservation of genome synteny are not always present in highly reduced genomes. Finally, the simplified molecular machinery of some of these organisms with small genomes may be used to aid in the design of artificial minimal cells. Here we review recent genomic discoveries of the biology of prokaryotes endowed with small gene sets and discuss the evolutionary mechanisms that have been proposed to explain their peculiar nature.

Genomic deletions disrupt nitrogen metabolism pathways of a cyanobacterial diatom symbiont

Diatoms with symbiotic N₂-fixing cyanobacteria are often abundant in the oligotrophic open ocean gyres. The most abundant cyanobacterial symbionts form heterocysts (specialized cells for N₂ fixation) and provide nitrogen (N) to their hosts, but their morphology, cellular locations and abundances differ depending on the host. Here we show that the location of the symbiont and its dependency on the host are linked to the evolution of the symbiont genome. The genome of Richelia (found inside the siliceous frustule of Hemiaulus) is reduced and lacks ammonium transporters, nitrate/nitrite reductases and glutamine:2-oxoglutarate aminotransferase. In contrast, the genome of the closely related Calothrix (found outside the frustule of Chaetoceros) is more similar to those of free-living heterocyst-forming cyanobacteria. The genome of Richelia is an example of metabolic streamlining that has implications for the evolution of N₂-fixing symbiosis and potentially for manipulating plant–cyanobacterial interactions.

Cyanobacterial symbionts of marine diatoms can localize intracellularly or externally to their host partners. Here Hilton et al. describe the genomes of two diazotroph cyanobacterial symbionts of diatoms and show that the location of the symbiont affects expression of nitrogen assimilation genes.

Polynucleobacter necessarius, a model for genome reduction in both free-living and symbiotic bacteria

We present the complete genomic sequence of the essential symbiont Polynucleobacter necessarius (Betaproteobacteria), which is a valuable case study for several reasons. First, it is hosted by a ciliated protist, Euplotes; bacterial symbionts of ciliates are still poorly known because of a lack of extensive molecular data. Second, the single species P. necessarius contains both symbiotic and free-living strains, allowing for a comparison between closely related organisms with different ecologies. Third, free-living P. necessarius strains are exceptional by themselves because of their small genome size, reduced metabolic flexibility, and high worldwide abundance in freshwater systems. We provide a comparative analysis of P. necessarius metabolism and explore the peculiar features of a genome reduction that occurred on an already streamlined genome. We compare this unusual system with current hypotheses for genome erosion in symbionts and free-living bacteria, propose modifications to the presently accepted model, and discuss the potential consequences of translesion DNA polymerase loss.

Bacterial cell wall synthesis gene uppP is required for Burkholderia colonization of the Stinkbug Gut

To establish a host-bacterium symbiotic association, a number of factors involved in symbiosis must operate in a coordinated manner. In insects, bacterial factors for symbiosis have been poorly characterized at the molecular and biochemical levels, since many symbionts have not yet been cultured or are as yet genetically intractable. Recently, the symbiotic association between a stinkbug, Riptortus pedestris, and its beneficial gut bacterium, Burkholderia sp., has emerged as a promising experimental model system, providing opportunities to study insect symbiosis using genetically manipulated symbiotic bacteria. Here, in search of bacterial symbiotic factors, we targeted cell wall components of the Burkholderia symbiont by disruption of uppP gene, which encodes undecaprenyl pyrophosphate phosphatase involved in biosynthesis of various bacterial cell wall components. Under culture conditions, the ΔuppP mutant showed higher susceptibility to lysozyme than the wild-type strain, indicating impaired integrity of peptidoglycan of the mutant. When administered to the host insect, the ΔuppP mutant failed to establish normal symbiotic association: the bacterial cells reached to the symbiotic midgut but neither proliferated nor persisted there. Transformation of the ΔuppP mutant with uppP-encoding plasmid complemented these phenotypic defects: lysozyme susceptibility in vitro was restored, and normal infection and proliferation in the midgut symbiotic organ were observed in vivo. The ΔuppP mutant also exhibited susceptibility to hypotonic, hypertonic, and centrifugal stresses. These results suggest that peptidoglycan cell wall integrity is a stress resistance factor relevant to the successful colonization of the stinkbug midgut by Burkholderia symbiont.

Genomic analysis of the regulatory elements and links with intrinsic DNA structural properties in the shrunken genome of Buchnera

Background

Buchnera aphidicola is an obligate symbiotic bacterium, associated with most of the aphididae, whose genome has drastically shrunk during intracellular evolution. Gene regulation in Buchnera has been a matter of controversy in recent years as the combination of genomic information with the experimental results has been contradictory, refuting or arguing in favour of a functional and responsive transcription regulation in Buchnera.

The goal of this study was to describe the gene transcription regulation capabilities of Buchnera based on the inventory of cis- and trans-regulators encoded in the genomes of five strains from different aphids (Acyrthosiphon pisum, Schizaphis graminum, Baizongia pistacea, Cinara cedri and Cinara tujafilina), as well as on the characterisation of some intrinsic structural properties of the DNA molecule in these bacteria.

Results

Interaction graph analysis shows that gene neighbourhoods are conserved between E. coli and Buchnera in structures called transcriptons, interactons and metabolons, indicating that selective pressures have acted on the evolution of transcriptional, protein-protein interaction and metabolic networks in Buchnera. The transcriptional regulatory network in Buchnera is composed of a few general DNA-topological regulators (Nucleoid Associated Proteins and topoisomerases), with the quasi-absence of any specific ones (except for multifunctional enzymes with a known gene expression regulatory role in Escherichia coli, such as AlaS, PepA and BolA, and the uncharacterized hypothetical regulators YchA and YrbA). The relative positioning of regulatory genes along the chromosome of Buchnera seems to have conserved its ancestral state, despite the genome erosion. Sigma-70 promoters with canonical thermodynamic sequence profiles were detected upstream of about 94% of the CDS of Buchnera in the different aphids. Based on Stress-Induced Duplex Destabilization (SIDD) measurements, unstable σ⁷⁰ promoters were found specifically associated with the regulator and transporter genes.

Conclusions

This genomic analysis provides supporting evidence of a selection of functional regulatory structures and it has enabled us to propose hypotheses concerning possible links between these regulatory elements and the DNA-topology (i.e., supercoiling, curvature, flexibility and base-pair stability) in the regulation of gene expression in the shrunken genome of Buchnera.

Decay of genes encoding the oomycete flagellar proteome in the downy mildew Hyaloperonospora arabidopsidis

Zoospores are central to the life cycles of most of the eukaryotic microbes known as oomycetes, but some genera have lost the ability to form these flagellated cells. In the plant pathogen Phytophthora infestans, genes encoding 257 proteins associated with flagella were identified by comparative genomics. These included the main structural components of the axoneme and basal body, proteins involved in intraflagellar transport, regulatory proteins, enzymes for maintaining ATP levels, and others. Transcripts for over three-quarters of the genes were up-regulated during sporulation, and persisted to varying degrees in the pre-zoospore stage (sporangia) and motile zoospores. Nearly all of these genes had orthologs in other eukaryotes that form flagella or cilia, but not species that lack the organelle. Orthologs of 211 of the genes were also absent from a sister taxon to P. infestans that lost the ability to form flagella, the downy mildew Hyaloperonospora arabidopsidis. Many of the genes retained in H. arabidopsidis were also present in other non-flagellates, suggesting that they play roles both in flagella and other cellular processes. Remnants of the missing genes were often detected in the H. arabidopsidis genome. Degradation of the genes was associated with local compaction of the chromosome and a heightened propensity towards genome rearrangements, as such regions were less likely to share synteny with P. infestans.

The diversity of actinorhizal symbiosis

Filamentous aerobic soil actinobacteria of the genus Frankia can induce the formation of nitrogen-fixing nodules on the roots of a diverse group of plants from eight dicotyledonous families, collectively called actinorhizal plants. Within nodules, Frankia can fix nitrogen while being hosted inside plant cells. Like in legume/rhizobia symbioses, bacteria can enter the plant root either intracellularly through an infection thread formed in a curled root hair, or intercellularly without root hair involvement, and the entry mechanism is determined by the host plant species. Nodule primordium formation is induced in the root pericycle as for lateral root primordia. Mature actinorhizal nodules are coralloid structures consisting of multiple lobes, each of which represents a modified lateral root without a root cap, a superficial periderm and with infected cells in the expanded cortex. In this review, an overview of nodule induction mechanisms and nodule structure is presented including comparisons with the corresponding mechanisms in legume symbioses.

Uncovering the co-evolutionary network among prokaryotic genes

MOTIVATION

Correlated events of gains and losses enable inference of co-evolution relations. The reconstruction of the co-evolutionary interactions network in prokaryotic species may elucidate functional associations among genes.

RESULTS

We developed a novel probabilistic methodology for the detection of co-evolutionary interactions between pairs of genes. Using this method we inferred the co-evolutionary network among 4593 Clusters of Orthologous Genes (COGs). The number of co-evolutionary interactions substantially differed among COGs. Over 40% were found to co-evolve with at least one partner. We partitioned the network of co-evolutionary relations into clusters and uncovered multiple modular assemblies of genes with clearly defined functions. Finally, we measured the extent to which co-evolutionary relations coincide with other cellular relations such as genomic proximity, gene fusion propensity, co-expression, protein-protein interactions and metabolic connections. Our results show that co-evolutionary relations only partially overlap with these other types of networks. Our results suggest that the inferred co-evolutionary network in prokaryotes is highly informative towards revealing functional relations among genes, often showing signals that cannot be extracted from other network types.

AVAILABILITY AND IMPLEMENTATION

Available under GPL license as open source.

CONTACT

talp@post.tau.ac.il.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The irreversible loss of a decomposition pathway marks the single origin of an ectomycorrhizal symbiosis

Microbial symbioses have evolved repeatedly across the tree of life, but the genetic changes underlying transitions to symbiosis are largely unknown, especially for eukaryotic microbial symbionts. We used the genus Amanita, an iconic group of mushroom-forming fungi engaged in ectomycorrhizal symbioses with plants, to identify both the origins and potential genetic changes maintaining the stability of this mutualism. A multi-gene phylogeny reveals one origin of the symbiosis within Amanita, with a single transition from saprotrophic decomposition of dead organic matter to biotrophic dependence on host plants for carbon. Associated with this transition are the losses of two cellulase genes, each of which plays a critical role in extracellular decomposition of organic matter. However a third gene, which acts at later stages in cellulose decomposition, is retained by many, but not all, ectomycorrhizal species. Experiments confirm that symbiotic Amanita species have lost the ability to grow on complex organic matter and have therefore lost the capacity to live in forest soils without carbon supplied by a host plant. Irreversible losses of decomposition pathways are likely to play key roles in the evolutionary stability of these ubiquitous mutualisms.

Selection-driven gene loss in bacteria

Gene loss by deletion is a common evolutionary process in bacteria, as exemplified by bacteria with small genomes that have evolved from bacteria with larger genomes by reductive processes. The driving force(s) for genome reduction remains unclear, and here we examined the hypothesis that gene loss is selected because carriage of superfluous genes confers a fitness cost to the bacterium. In the bacterium Salmonella enterica, we measured deletion rates at 11 chromosomal positions and the fitness effects of several spontaneous deletions. Deletion rates varied over 200-fold between different regions with the replication terminus region showing the highest rates. Approximately 25% of the examined deletions caused an increase in fitness under one or several growth conditions, and after serial passage of wild-type bacteria in rich medium for 1,000 generations we observed fixation of deletions that substantially increased bacterial fitness when reconstructed in a non-evolved bacterium. These results suggest that selection could be a significant driver of gene loss and reductive genome evolution.

On the evolutionary ecology of symbioses between chemosynthetic bacteria and bivalves

Mutualistic associations between bacteria and eukaryotes occur ubiquitously in nature, forming the basis for key ecological and evolutionary innovations. Some of the most prominent examples of these symbioses are chemosynthetic bacteria and marine invertebrates living in the absence of sunlight at deep-sea hydrothermal vents and in sediments rich in reduced sulfur compounds. Here, chemosynthetic bacteria living in close association with their hosts convert CO₂ or CH₄ into organic compounds and provide the host with necessary nutrients. The dominant macrofauna of hydrothermal vent and cold seep ecosystems all depend on the metabolic activity of chemosynthetic bacteria, which accounts for almost all primary production in these complex ecosystems. Many of these enigmatic mutualistic associations are found within the molluscan class Bivalvia. Currently, chemosynthetic symbioses have been reported from five distinct bivalve families (Lucinidae, Mytilidae, Solemyidae, Thyasiridae, and Vesicomyidae). This brief review aims to provide an overview of the diverse physiological and genetic adaptations of symbiotic chemosynthetic bacteria and their bivalve hosts.

The SAR11 group of alpha-proteobacteria is not related to the origin of mitochondria

Although free living, members of the successful SAR11 group of marine alpha-proteobacteria contain a very small and A+T rich genome, two features that are typical of mitochondria and related obligate intracellular parasites such as the Rickettsiales. Previous phylogenetic analyses have suggested that Candidatus Pelagibacter ubique, the first cultured member of this group, is related to the Rickettsiales+mitochondria clade whereas others disagree with this conclusion. In order to determine the evolutionary position of the SAR11 group and its relationship to the origin of mitochondria, we have performed phylogenetic analyses on the concatenation of 24 proteins from 5 mitochondria and 71 proteobacteria. Our results support that SAR11 group is not the sistergroup of the Rickettsiales+mitochondria clade and confirm that the position of this group in the alpha-proteobacterial tree is strongly affected by tree reconstruction artefacts due to compositional bias. As a consequence, genome reduction and bias toward a high A+T content may have evolved independently in the SAR11 species, which points to a different direction in the quest for the closest relatives to mitochondria and Rickettsiales. In addition, our analyses raise doubts about the monophyly of the newly proposed Pelagibacteraceae family.

Reduced selective constraint in endosymbionts: elevation in radical amino acid replacements occurs genome-wide

As predicted by the nearly neutral model of evolution, numerous studies have shown that reduced N(e) accelerates the accumulation of slightly deleterious changes under genetic drift. While such studies have mostly focused on eukaryotes, bacteria also offer excellent models to explore the effects of N(e). Most notably, the genomes of host-dependent bacteria with small N(e) show signatures of genetic drift, including elevated K(a)/K(s). Here, I explore the utility of an alternative measure of selective constraint: the per-site rate of radical and conservative amino acid substitutions (D(r)/D(c)). I test the hypothesis that purifying selection against radical amino acid changes is less effective in two insect endosymbiont groups (Blochmannia of ants and Buchnera of aphids), compared to related gamma-Proteobacteria. Genome comparisons demonstrate a significant elevation in D(r)/D(c) in endosymbionts that affects the majority (66-79%) of shared orthologs examined. The elevation of D(r)/D(c) in endosymbionts affects all functional categories examined. Simulations indicate that D(r)/D(c) estimates are sensitive to codon frequencies and mutational parameters; however, estimation biases occur in the opposite direction as the patterns observed in genome comparisons, thereby making the inference of elevated D(r)/D(c) more conservative. Increased D(r)/D(c) and other signatures of genome degradation in endosymbionts are consistent with strong effects of genetic drift in their small populations, as well as linkage to selected sites in these asexual bacteria. While relaxed selection against radical substitutions may contribute, genome-wide processes such as genetic drift and linkage best explain the pervasive elevation in D(r)/D(c) across diverse functional categories that include basic cellular processes. Although the current study focuses on a few bacterial lineages, it suggests D(r)/D(c) is a useful gauge of selective constraint and may provide a valuable alternative to K(a)/K(s) when high sequence divergences preclude estimates of K(s). Broader application of D(r)/D(c) will benefit from approaches less prone to estimation biases.