Host-symbiont stability and fast evolutionary rates in an ant-bacterium association: cospeciation of camponotus species and their endosymbionts, candidatus blochmannia

Bacterial endosymbionts are widespread across several insect orders and are involved in interactions ranging from obligate mutualism to reproductive parasitism. Candidatus Blochmannia gen. nov. (Blochmannia) is an obligate bacterial associate of Camponotus and related ant genera (Hymenoptera: Formicidae). The occurrence of Blochmannia in all Camponotus species sampled from field populations and its maternal transmission to host offspring suggest that this bacterium is engaged in a long-term, stable association with its ant hosts. However, evidence for cospeciation in this system is equivocal because previous phylogenetic studies were based on limited gene sampling, lacked statistical analysis of congruence, and have even suggested host switching. We compared phylogenies of host genes (the nuclear EF-1alphaF2 and mitochondrial COI/II) and Blochmannia genes (16S ribosomal DNA [rDNA], groEL, gidA, and rpsB), totaling more than 7 kilobases for each of 16 Camponotus species. Each data set was analyzed using maximum likelihood and Bayesian phylogenetic reconstruction methods. We found minimal conflict among host and symbiont phylogenies, and the few areas of discordance occurred at deep nodes that were poorly supported by individual data sets. Concatenated protein-coding genes produced a very well-resolved tree that, based on the Shimodaira-Hasegawa test, did not conflict with any host or symbiont data set. Correlated rates of synonymous substitution (d(S)) along corresponding branches of host and symbiont phylogenies further supported the hypothesis of cospeciation. These findings indicate that Blochmannia-Camponotus symbiosis has been evolutionarily stable throughout tens of millions of years. Based on inferred divergence times among the ant hosts, we estimated rates of sequence evolution of Blochmannia to be approximately 0.0024 substitutions per site per million years (s/s/MY) for the 16S rDNA gene and approximately 0.1094 s/s/MY at synonymous positions of the genes sampled. These rates are several-fold higher than those for related bacteria Buchnera aphidicola and Escherichia coli. Phylogenetic congruence among Blochmannia genes indicates genome stability that typifies primary endosymbionts of insects.     

A conservative test of genetic drift in the endosymbiotic bacterium Buchnera: slightly deleterious mutations in the chaperonin groEL

The obligate endosymbiotic bacterium Buchnera aphidicola shows elevated rates of sequence evolution compared to free-living relatives, particularly at nonsynonymous sites. Because Buchnera experiences population bottlenecks during transmission to the offspring of its aphid host, it is hypothesized that genetic drift and the accumulation of slightly deleterious mutations can explain this rate increase. Recent studies of intraspecific variation in Buchnera reveal patterns consistent with this hypothesis. In this study, we examine inter- and intraspecific nucleotide variation in groEL, a highly conserved chaperonin gene that is constitutively overexpressed in Buchnera. Maximum-likelihood estimates of nonsynonymous substitution rates across Buchnera species are strikingly low at groEL compared to other loci. Despite this evidence for strong purifying selection on groEL, our intraspecific analysis of this gene documents reduced synonymous polymorphism, elevated nonsynonymous polymorphism, and an excess of rare alleles relative to the neutral expectation, as found in recent studies of other Buchnera loci. Comparisons with Escherichia coli generally show patterns predicted by their differences in N(e). The sum of these observations is not expected under relaxed or balancing selection, selective sweeps, or increased mutation rate. Rather, they further support the hypothesis that drift is an important force driving accelerated protein evolution in this obligate mutualist.     

Genome evolution in bacterial endosymbionts of insects

Many insect species rely on intracellular bacterial symbionts for their viability and fecundity. Large-scale DNA-sequence analyses are revealing the forces that shape the evolution of these bacterial associates and the genetic basis of their specialization to an intracellular lifestyle. The full genome sequences of two obligate mutualists, Buchnera aphidicola of aphids and Wigglesworthia glossinidia of tsetse flies, reveal substantial gene loss and an integration of host and symbiont metabolic functions. Further genomic comparisons should reveal the generality of these features among bacterial mutualists and the extent to which they are shared with other intracellular bacteria, including obligate pathogens.     

Nonhomogeneous model of sequence evolution indicates independent origins of primary endosymbionts within the enterobacteriales (gamma-Proteobacteria)

Standard methods of phylogenetic reconstruction are based on models that assume homogeneity of nucleotide composition among taxa. However, this assumption is often violated in biological data sets. In this study, we examine possible effects of nucleotide heterogeneity among lineages on the phylogenetic reconstruction of a bacterial group that spans a wide range of genomic nucleotide contents: obligately endosymbiotic bacteria and free-living or commensal species in the gamma-Proteobacteria. We focus on AT-rich primary endosymbionts to better understand the origins of obligately intracellular lifestyles. Previous phylogenetic analyses of this bacterial group point to the importance of accounting for base compositional variation in estimating relationships, particularly between endosymbiotic and free-living taxa. Here, we develop an approach to compare susceptibility of various phylogenetic reconstruction methods to the effects of nucleotide heterogeneity. First, we identify candidate trees of gamma-Proteobacteria groEL and 16S rRNA using approaches that assume homogeneous and stationary base composition, including Bayesian, maximum likelihood, parsimony, and distance methods. We then create permutations of the resulting candidate trees by varying the placement of the AT-rich endosymbiont Buchnera. These permutations are evaluated under the nonhomogeneous and nonstationary maximum likelihood model of Galtier and Gouy, which allows equilibrium base content to vary among examined lineages. Our results show that commonly used phylogenetic methods produce incongruent trees of the Enterobacteriales, and that the placement of Buchnera is especially unstable. However, under a nonhomogeneous model, various groEL and 16S rRNA phylogenies that separate Buchnera from other AT-rich endosymbionts (Blochmannia and Wigglesworthia) have consistently and significantly higher likelihood scores. Blochmannia and Wigglesworthia appear to have evolved from secondary endosymbionts, and represent an origin of primary endosymbiosis that is independent from Buchnera. This application of a nonhomogeneous model offers a computationally feasible way to test specific phylogenetic hypotheses for taxa with heterogeneous and nonstationary base composition.     

Parallel acceleration of evolutionary rates in symbiont genes underlying host nutrition

The overproduction of essential amino acids by Buchnera aphidicola, the primary bacterial mutualist of aphids, is considered an adaptation for increased production of nutrients that are lacking in aphids' diet of plant sap. Given their shared role in host nutrition, amino acid biosynthetic genes of Buchnera are expected to experience parallel changes in selection that depend on host diet quality, growth rate, and population structure. This study evaluates the hypothesis of parallel selection across biosynthetic pathways by testing for correlated changes in evolutionary rates at biosynthetic genes of Buchnera. Previous studies show fast evolutionary rates at tryptophan biosynthetic genes among Buchnera associated with the aphid genus Uroleucon and suggest reduced purifying selection on symbiont nutritional functions in this aphid group. Here, we test for parallel rate acceleration at other amino acid biosynthetic genes of Buchnera-Uroleucon, including those for leucine (leuABC) and isoleucine/valine biosynthesis (ilvC). Ratios of nonsynonymous to synonymous substitutions (d(N)/d(S)) were estimated using codon-based maximum-likelihood methods that account for the extreme AT compositional bias of Buchnera sequences. A significant elevation in d(N)/d(S) at biosynthetic loci but not at two housekeeping genes sampled (dnaN and tuf) suggests reduced host-level selection on biosynthetic capabilities of Buchnera-Uroleucon. In addition, the discovery of trpEG pseudogenes in Buchnera-U. obscurum further supports reduced selection on amino acid biosynthesis.     

A strong effect of AT mutational bias on amino acid usage in Buchnera is mitigated at high-expression genes

The advent of full genome sequences provides exceptionally rich data sets to explore molecular and evolutionary mechanisms that shape divergence among and within genomes. In this study, we use multivariate analysis to determine the processes driving genome-wide patterns of amino usage in the obligate endosymbiont Buchnera and its close free-living relative Escherichia coli. In the AT-rich Buchnera genome, the primary source of variation in amino acid usage differentiates high- and low-expression genes. Amino acids of high-expression Buchnera genes are generally less aromatic and use relatively GC-rich codons, suggesting that selection against aromatic amino acids and against amino acids with AT-rich codons is stronger in high-expression genes. Selection to maintain hydrophobic amino acids in integral membrane proteins is a primary factor driving protein evolution in E. coli but is a secondary factor in Buchnera. In E. coli, gene expression is a secondary force driving amino acid usage, and a correlation with tRNA abundance suggests that translational selection contributes to this effect. Although this and previous studies demonstrate that AT mutational bias and genetic drift influence amino acid usage in Buchnera, this genome-wide analysis argues that selection is sufficient to affect the amino acid content of proteins with different expression and hydropathy levels.     

Mutualism meltdown in insects: bacteria constrain thermal adaptation

Predicting whether and how organisms will successfully cope with climate change presents critical questions for biologists and environmental scientists. Models require knowing how organisms interact with their abiotic environment, as well understanding biotic interactions that include a network of symbioses in which all species are embedded. Bacterial symbionts of insects offer valuable models to examine how microbes can facilitate and constrain adaptation to a changing environment. While some symbionts confer plasticity that accelerates adaptation, long-term bacterial mutualists of insects are characterized by tight lifestyle constraints, genome deterioration, and vulnerability to thermal stress. These essential bacterial partners are eliminated at high temperatures, analogous to the loss of zooanthellae during coral bleaching. Recent field-based studies suggest that thermal sensitivity of bacterial mutualists constrains insect responses. In this sense, highly dependent mutualisms may be the Achilles' heel of thermal responses in insects.     

Toward a Wolbachia Multilocus Sequence Typing System: Discrimination of Wolbachia Strains Present in Drosophila Species

Among the diverse maternally inherited symbionts in arthropods, Wolbachia are the most common and infect over 20% of all species. In a departure from traditional genotyping or phylogenetic methods relying on single Wolbachia genes, the present study represents an initial Multilocus Sequence Typing (MLST) analysis to discriminate closely related Wolbachia pipientis strains, and additional data on sequence diversity in Wolbachia. We report a new phylogenetic characterization of four genes (aspC, atpD, sucB, and pdhB), and provide an expanded analysis of markers described in previous studies (16S rDNA, ftsZ, groEL, dnaA, and gltA). MLST analysis of the bacterial strains present in 16 different Drosophila–Wolbachia associations detected four distinct clonal complexes that also corresponded to maximum-likelihood identified phylogenetic clades. Among the 16 associations analyzed, six could not be assigned to MLST clonal complexes and were also shown to be in conflict with relationships predicted by maximum-likelihood phylogenetic inferences. The results demonstrate the discriminatory power of MLST for identifying strains and clonal lineages of Wolbachia and provide a robust foundation for studying the ecology and evolution of this widespread endosymbiont.     

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.