The evolution of the heat-shock protein GroEL from Buchnera, the primary endosymbiont of aphids, is governed by positive selection

A new lysine biosynthetic enzyme from a bacterial endosymbiont shaped by genetic drift and genome reduction

The effect of population bottlenecks and genome reduction on enzyme function is poorly understood. Candidatus Liberibacter solanacearum is a bacterium with a reduced genome that is transmitted vertically to the egg of an infected psyllid-a population bottleneck that imposes genetic drift and is predicted to affect protein structure and function. Here, we define the function of Ca. L. solanacearum dihydrodipicolinate synthase (CLsoDHDPS), which catalyzes the committed branchpoint reaction in diaminopimelate and lysine biosynthesis. We demonstrate that CLsoDHDPS is expressed in Ca. L. solanacearum and expression is increased ~2-fold in the insect host compared to in planta. CLsoDHDPS has decreased thermal stability and increased aggregation propensity, implying mutations have destabilized the enzyme but are compensated for through elevated chaperone expression and a stabilized oligomeric state. CLsoDHDPS uses a ternary-complex kinetic mechanism, which is to date unique among DHDPS enzymes, has unusually low catalytic ability, but an unusually high substrate affinity. Structural studies demonstrate that the active site is more open, and the structure of CLsoDHDPS with both pyruvate and the substrate analogue succinic-semialdehyde reveals that the product is both structurally and energetically different and therefore evolution has in this case fashioned a new enzyme. Our study suggests the effects of genome reduction and genetic drift on the function of essential enzymes and provides insights on bacteria-host co-evolutionary associations. We propose that bacteria with endosymbiotic lifestyles present a rich vein of interesting enzymes useful for understanding enzyme function and/or informing protein engineering efforts.

In vivo interference of pea aphid endosymbiont Buchnera groEL gene by synthetic peptide nucleic acids

The unculturable nature of intracellular obligate symbionts presents a significant challenge for elucidating gene functionality, necessitating the development of gene manipulation techniques. One of the best-studied obligate symbioses is that between aphids and the bacterial endosymbiont Buchnera aphidicola. Given the extensive genome reduction observed in Buchnera, the remaining genes are crucial for understanding the host-symbiont relationship, but a lack of tools for manipulating gene function in the endosymbiont has significantly impeded the exploration of the molecular mechanisms underlying this mutualism. In this study, we introduced a novel gene manipulation technique employing synthetic single-stranded peptide nucleic acids (PNAs). We targeted the critical Buchnera groEL using specially designed antisense PNAs conjugated to an arginine-rich cell-penetrating peptide (CPP). Within 24 h of PNA administration via microinjection, we observed a significant reduction in groEL expression and Buchnera cell count. Notably, the interference of groEL led to profound morphological malformations in Buchnera, indicative of impaired cellular integrity. The gene knockdown technique developed in this study, involving the microinjection of CPP-conjugated antisense PNAs, provides a potent approach for in vivo gene manipulation of unculturable intracellular symbionts, offering valuable insights into their biology and interactions with hosts.

A complex interplay of evolutionary forces continues to shape ancient co-occurring symbiont genomes

Many insects depend on ancient associations with intracellular bacteria for essential nutrition. The genomes of these bacteria are often highly reduced. Although drift is a major driver of symbiont evolution, other evolutionary forces continue to influence them. To understand how ongoing molecular evolution and gene loss shape symbiont genomes, we sequenced two of the most ancient symbionts known, Sulcia and Nasuia, from 20 Hawaiian Nesophrosyne leafhoppers. We leveraged the parallel divergence of Nesophrosyne lineages throughout Hawaii as a natural experimental framework. Sulcia and Nasuia experience ongoing-but divergent-gene loss, often in a convergent fashion. Although some genes are under relaxed selection, purifying and positive selection are also important drivers of genome evolution, particularly in maintaining certain nutritional and cellular functions. Our results further demonstrate that symbionts experience dramatically different evolutionary environments, as evidenced by the finding that Sulcia and Nasuia have one of the slowest and fastest rates of molecular evolution known.

Potato leafroll virus reduces Buchnera aphidocola titer and alters vector transcriptome responses

Viruses in the Luteoviridae family, such as Potato leafroll virus (PLRV), are transmitted by aphids in a circulative and nonpropagative mode. This means the virions enter the aphid body through the gut when they feed from infected plants and then the virions circulate through the hemolymph to enter the salivary glands before being released into the saliva. Although these viruses do not replicate in their insect vectors, previous studies have demonstrated viruliferous aphid behavior is altered and the obligate symbiont of aphids, Buchnera aphidocola, may be involved in transmission. Here we provide the transcriptome of green peach aphids (Myzus persicae) carrying PLRV and virus-free control aphids using Illumina sequencing. Over 150 million paired-end reads were obtained through Illumina sequencing, with an average of 19 million reads per library. The comparative analysis identified 134 differentially expressed genes (DEGs) between the M. persicae transcriptomes, including 64 and 70 genes that were up- and down-regulated in aphids carrying PLRV, respectively. Using functional classification in the GO databases, 80 of the DEGs were assigned to 391 functional subcategories at category level 2. The most highly up-regulated genes in aphids carrying PLRV were cytochrome p450s, genes related to cuticle production, and genes related to development, while genes related to heat shock proteins, histones, and histone modification were the most down-regulated. PLRV aphids had reduced Buchnera titer and lower abundance of several Buchnera transcripts related to stress responses and metabolism. These results suggest carrying PLRV may reduce both aphid and Buchnera genes in response to stress. This work provides valuable basis for further investigation into the complicated mechanisms of circulative and nonpropagative transmission.

The holobiont transcriptome of teneral tsetse fly species of varying vector competence

Background

Tsetse flies are the obligate vectors of African trypanosomes, which cause Human and Animal African Trypanosomiasis. Teneral flies (newly eclosed adults) are especially susceptible to parasite establishment and development, yet our understanding of why remains fragmentary. The tsetse gut microbiome is dominated by two Gammaproteobacteria, an essential and ancient mutualist Wigglesworthia glossinidia and a commensal Sodalis glossinidius. Here, we characterize and compare the metatranscriptome of teneral Glossina morsitans to that of G. brevipalpis and describe unique immunological, physiological, and metabolic landscapes that may impact vector competence differences between these two species.

Results

An active expression profile was observed for Wigglesworthia immediately following host adult metamorphosis. Specifically, ‘translation, ribosomal structure and biogenesis’ followed by ‘coenzyme transport and metabolism’ were the most enriched clusters of orthologous genes (COGs), highlighting the importance of nutrient transport and metabolism even following host species diversification. Despite the significantly smaller Wigglesworthia genome more differentially expressed genes (DEGs) were identified between interspecific isolates (n = 326, ~ 55% of protein coding genes) than between the corresponding Sodalis isolates (n = 235, ~ 5% of protein coding genes) likely reflecting distinctions in host co-evolution and adaptation. DEGs between Sodalis isolates included genes involved in chitin degradation that may contribute towards trypanosome susceptibility by compromising the immunological protection provided by the peritrophic matrix. Lastly, G. brevipalpis tenerals demonstrate a more immunologically robust background with significant upregulation of IMD and melanization pathways.

Conclusions

These transcriptomic differences may collectively contribute to vector competence differences between tsetse species and offers translational relevance towards the design of novel vector control strategies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07729-5.

Genetic Constraints, Transcriptome Plasticity, and the Evolutionary Response to Climate Change

In situ adaptation to climate change will be critical for the persistence of many ectotherm species due to their relative lack of dispersal capacity. Climate change is causing increases in both the mean and the variance of environmental temperature, each of which may act as agents of selection on different traits. Importantly, these traits may not be heritable or have the capacity to evolve independently from one another. When genetic constraints prevent the "baseline" values of thermal performance traits from evolving rapidly, phenotypic plasticity driven by gene expression might become critical. We review the literature for evidence that thermal performance traits in ectotherms are heritable and have genetic architectures that permit their unconstrained evolution. Next, we examine the relationship between gene expression and both the magnitude and duration of thermal stress. Finally, we identify genes that are likely to be important for adaptation to a changing climate and determine whether they show patterns consistent with thermal adaptation. Although few studies have measured narrow-sense heritabilities of thermal performance traits, current evidence suggests that the end points of thermal reaction norms (tolerance limits) are moderately heritable and have the potential to evolve rapidly. By contrast, performance at intermediate temperatures has substantially lower evolutionary potential. Moreover, evolution in many species appears to be constrained by genetic correlations such that populations can adapt to either increases in mean temperature or temperature variability, but not both. Finally, many species have the capacity for plastic expression of the transcriptome in response to temperature shifts, with the number of differentially expressed genes increasing with the magnitude, but not the duration, of thermal stress. We use these observations to develop a conceptual model that describes the likely trajectory of genome evolution in response to changes in environmental temperature. Our results indicate that extreme weather events, rather than gradual increases in mean temperature, are more likely to drive genetic and phenotypic change in wild ectotherms.

Uniparental Inheritance Promotes Adaptive Evolution in Cytoplasmic Genomes

Eukaryotes carry numerous asexual cytoplasmic genomes (mitochondria and plastids). Lacking recombination, asexual genomes should theoretically suffer from impaired adaptive evolution. Yet, empirical evidence indicates that cytoplasmic genomes experience higher levels of adaptive evolution than predicted by theory. In this study, we use a computational model to show that the unique biology of cytoplasmic genomes-specifically their organization into host cells and their uniparental (maternal) inheritance-enable them to undergo effective adaptive evolution. Uniparental inheritance of cytoplasmic genomes decreases competition between different beneficial substitutions (clonal interference), promoting the accumulation of beneficial substitutions. Uniparental inheritance also facilitates selection against deleterious cytoplasmic substitutions, slowing Muller's ratchet. In addition, uniparental inheritance generally reduces genetic hitchhiking of deleterious substitutions during selective sweeps. Overall, uniparental inheritance promotes adaptive evolution by increasing the level of beneficial substitutions relative to deleterious substitutions. When we assume that cytoplasmic genome inheritance is biparental, decreasing the number of genomes transmitted during gametogenesis (bottleneck) aids adaptive evolution. Nevertheless, adaptive evolution is always more efficient when inheritance is uniparental. Our findings explain empirical observations that cytoplasmic genomes-despite their asexual mode of reproduction-can readily undergo adaptive evolution.

Effects of different temperature regimes on survival of Diaphorina citri and its endosymbiotic bacterial communities

The Asian citrus psyllid, Diaphorina citri, is a major pest of citrus and vector of citrus greening (huanglongbing) in Asian. In our field-collected psyllid samples, we discovered that Fuzhou (China) and Faisalabad (Pakistan), populations harbored an obligate primary endosymbiont Candidatus Carsonella (gen. nov.) with a single species, Candidatus Carsonella ruddii (sp. nov.) and a secondary endosymbiont, Wolbachia surface proteins (WSP) which are intracellular endosymbionts residing in the bacteriomes. Responses of these symbionts to different temperatures were examined and their host survival assessed. Diagnostic PCR assays showed that the endosymbionts infection rates were not significantly reduced in both D. citri populations after 24 h exposure to cold or heat treatments. Although quantitative PCR assays showed significant reduction of WSP relative densities at 40°C for 24 h, a substantial decrease occurred as the exposure duration increased beyond 3 days. Under the same temperature regimes, Ca. C. ruddii density was initially less affected during the first exposure day, but rapidly reduced at 3-5 days compared to WSP. However, the mortality of the psyllids increased rapidly as exposure time to heat treatment increased. The responses of the two symbionts to unfavorable temperature regimes highlight the complex host-symbionts interactions between D. citri and its associated endosymbionts.

Chance and necessity in the genome evolution of endosymbiotic bacteria of insects

An open question in evolutionary biology is how does the selection–drift balance determine the fates of biological interactions. We searched for signatures of selection and drift in genomes of five endosymbiotic bacterial groups known to evolve under strong genetic drift. Although most genes in endosymbiotic bacteria showed evidence of relaxed purifying selection, many genes in these bacteria exhibited stronger selective constraints than their orthologs in free-living bacterial relatives. Remarkably, most of these highly constrained genes had no role in the host–symbiont interactions but were involved in either buffering the deleterious consequences of drift or other host-unrelated functions, suggesting that they have either acquired new roles or their role became more central in endosymbiotic bacteria. Experimental evolution of Escherichia coli under strong genetic drift revealed remarkable similarities in the mutational spectrum, genome reduction patterns and gene losses to endosymbiotic bacteria of insects. Interestingly, the transcriptome of the experimentally evolved lines showed a generalized deregulation of the genome that affected genes encoding proteins involved in mutational buffering, regulation and amino acid biosynthesis, patterns identical to those found in endosymbiotic bacteria. Our results indicate that drift has shaped endosymbiotic associations through a change in the functional landscape of bacterial genes and that the host had only a small role in such a shift.

Bridging the physical scales in evolutionary biology: from protein sequence space to fitness of organisms and populations

Bridging the gap between the molecular properties of proteins and organismal/population fitness is essential for understanding evolutionary processes. This task requires the integration of the several physical scales of biological organization, each defined by a distinct set of mechanisms and constraints, into a single unifying model. The molecular scale is dominated by the constraints imposed by the physico-chemical properties of proteins and their substrates, which give rise to trade-offs and epistatic (non-additive) effects of mutations. At the systems scale, biological networks modulate protein expression and can either buffer or enhance the fitness effects of mutations. The population scale is influenced by the mutational input, selection regimes, and stochastic changes affecting the size and structure of populations, which eventually determine the evolutionary fate of mutations. Here, we summarize the recent advances in theory, computer simulations, and experiments that advance our understanding of the links between various physical scales in biology.

Repeated replacement of an intrabacterial symbiont in the tripartite nested mealybug symbiosis

Stable endosymbiosis of a bacterium into a host cell promotes cellular and genomic complexity. The mealybug Planococcus citri has two bacterial endosymbionts with an unusual nested arrangement: the γ-proteobacterium Moranella endobia lives in the cytoplasm of the β-proteobacterium Tremblaya princeps These two bacteria, along with genes horizontally transferred from other bacteria to the P. citri genome, encode gene sets that form an interdependent metabolic patchwork. Here, we test the stability of this three-way symbiosis by sequencing host and symbiont genomes for five diverse mealybug species and find marked fluidity over evolutionary time. Although Tremblaya is the result of a single infection in the ancestor of mealybugs, the γ-proteobacterial symbionts result from multiple replacements of inferred different ages from related but distinct bacterial lineages. Our data show that symbiont replacement can happen even in the most intricate symbiotic arrangements and that preexisting horizontally transferred genes can remain stable on genomes in the face of extensive symbiont turnover.

Genome Evolution in the Obligate but Environmentally Active Luminous Symbionts of Flashlight Fish

The luminous bacterial symbionts of anomalopid flashlight fish are thought to be obligately dependent on their hosts for growth and share several aspects of genome evolution with unrelated obligate symbionts, including genome reduction. However, in contrast to most obligate bacteria, anomalopid symbionts have an active environmental phase that may be important for symbiont transmission. Here we investigated patterns of evolution between anomalopid symbionts compared with patterns in free-living relatives and unrelated obligate symbionts to determine if trends common to obligate symbionts are also found in anomalopid symbionts. Two symbionts, "Candidatus Photodesmus katoptron" and "Candidatus Photodesmus blepharus," have genomes that are highly similar in gene content and order, suggesting genome stasis similar to ancient obligate symbionts present in insect lineages. This genome stasis exists in spite of the symbiont's inferred ability to recombine, which is frequently lacking in obligate symbionts with stable genomes. Additionally, we used genome comparisons and tests of selection to infer which genes may be particularly important for the symbiont's ecology compared with relatives. In keeping with obligate dependence, substitution patterns suggest that most symbiont genes are experiencing relaxed purifying selection compared with relatives. However, genes involved in motility and carbon storage, which are likely to be used outside the host, appear to be under increased purifying selection. Two chemoreceptor chemotaxis genes are retained by both species and show high conservation with amino acid sensing genes, suggesting that the bacteria may actively seek out hosts using chemotaxis toward amino acids, which the symbionts are not able to synthesize.

A Tale of Transmission: Aeromonas veronii Activity within Leech-Exuded Mucus

Transmission, critical to the establishment and persistence of host-associated microbiotas, also exposes symbionts to new environmental conditions. With horizontal transmission, these different conditions represent major lifestyle shifts. Yet genome-wide analyses of how microbes adjust their transcriptomes toward these dramatic shifts remain understudied. Here, we provide a comprehensive and comparative analysis of the global transcriptional profiles of a symbiont as it shifts between lifestyles during transmission. The gammaproteobacterium Aeromonas veronii is transmitted from the gut of the medicinal leech to other hosts via host mucosal castings, yet A. veronii can also transition from mucosal habitancy to a free-living lifestyle. These three lifestyles are characterized by distinct physiological constraints and consequently lifestyle-specific changes in the expression of stress-response genes. Mucus-bound A. veronii had the greatest expression in terms of both the number of loci and levels of transcription of stress-response mechanisms. However, these bacteria are still capable of proliferating within the mucus, suggesting the availability of nutrients within this environment. We found that A. veronii alters transcription of loci in a synthetic pathway that obtains and incorporates N-acetylglucosamine (NAG; a major component of mucus) into the bacterial cell wall, enabling proliferation. Our results demonstrate that symbionts undergo dramatic local adaptation, demonstrated by widespread transcriptional changes, throughout the process of transmission that allows them to thrive while they encounter new environments which further shape their ecology and evolution.

Hsp90: A Global Regulator of the Genotype-to-Phenotype Map in Cancers

Cancer cells have the unusual capacity to limit the cost of the mutation load that they harbor and simultaneously harness its evolutionary potential. This property fuels drug resistance, a key failure mode in oncogene-directed therapy. However, the factors that regulate this capacity might also provide an Achilles' heel that could be exploited therapeutically. Recently, insight has come from a seemingly distant field: protein folding. It is now clear that protein homeostasis broadly supports malignancy and fuels the rapid evolution of drug resistance. Among protein homeostatic mechanisms that influence cancer biology, the essential ATP-driven molecular chaperone heat-shock protein 90 (Hsp90) is especially important. Hsp90 catalyzes folding of many proteins that regulate growth and development. These "client" kinases, transcription factors, and ubiquitin ligases often play critical roles in human disease, especially cancer. Studies in a wide range of systems-from single-celled organisms to human tumor samples-suggest that Hsp90 can broadly reshape the map between genotype and phenotype, acting as a "capacitor" and "potentiator" of genetic variation. Indeed, it has likely done so to such a degree that it has left an impress on diverse genome sequences. Hsp90 can constitute as much as 5% of total protein in transformed cells and increased levels of heat-shock activation correlate with poor prognosis in breast cancer. These findings and others have motivated a flurry of interest in Hsp90 inhibitors as cancer therapeutics, which have met with rather limited success as single agents, but may eventually prove invaluable in limiting the emergence of resistance to other chemotherapeutics, both genotoxic and molecularly targeted. Here, we provide an overview of Hsp90 function, review its relationship to genetic variation and the evolution of new traits, and discuss the importance of these findings for cancer biology and future efforts to drug this pathway.

Trends in substitution models of molecular evolution

Substitution models of evolution describe the process of genetic variation through fixed mutations and constitute the basis of the evolutionary analysis at the molecular level. Almost 40 years after the development of first substitution models, highly sophisticated, and data-specific substitution models continue emerging with the aim of better mimicking real evolutionary processes. Here I describe current trends in substitution models of DNA, codon and amino acid sequence evolution, including advantages and pitfalls of the most popular models. The perspective concludes that despite the large number of currently available substitution models, further research is required for more realistic modeling, especially for DNA coding and amino acid data. Additionally, the development of more accurate complex models should be coupled with new implementations and improvements of methods and frameworks for substitution model selection and downstream evolutionary analysis.

Dynamics of Wolbachia pipientis Gene Expression Across the Drosophila melanogaster Life Cycle

Symbiotic interactions between microbes and their multicellular hosts have manifold biological consequences. To better understand how bacteria maintain symbiotic associations with animal hosts, we analyzed genome-wide gene expression for the endosymbiotic α-proteobacteria Wolbachia pipientis across the entire life cycle of Drosophila melanogaster. We found that the majority of Wolbachia genes are expressed stably across the D. melanogaster life cycle, but that 7.8% of Wolbachia genes exhibit robust stage- or sex-specific expression differences when studied in the whole-organism context. Differentially-expressed Wolbachia genes are typically up-regulated after Drosophila embryogenesis and include many bacterial membrane, secretion system, and ankyrin repeat-containing proteins. Sex-biased genes are often organized as small operons of uncharacterized genes and are mainly up-regulated in adult Drosophila males in an age-dependent manner. We also systematically investigated expression levels of previously-reported candidate genes thought to be involved in host-microbe interaction, including those in the WO-A and WO-B prophages and in the Octomom region, which has been implicated in regulating bacterial titer and pathogenicity. Our work provides comprehensive insight into the developmental dynamics of gene expression for a widespread endosymbiont in its natural host context, and shows that public gene expression data harbor rich resources to probe the functional basis of the Wolbachia-Drosophila symbiosis and annotate the transcriptional outputs of the Wolbachia genome.

Metabolic Coevolution in the Bacterial Symbiosis of Whiteflies and Related Plant Sap-Feeding Insects

Genomic decay is a common feature of intracellular bacteria that have entered into symbiosis with plant sap-feeding insects. This study of the whitefly Bemisia tabaci and two bacteria (Portiera aleyrodidarum and Hamiltonella defensa) cohoused in each host cell investigated whether the decay of Portiera metabolism genes is complemented by host and Hamiltonella genes, and compared the metabolic traits of the whitefly symbiosis with other sap-feeding insects (aphids, psyllids, and mealybugs). Parallel genomic and transcriptomic analysis revealed that the host genome contributes multiple metabolic reactions that complement or duplicate Portiera function, and that Hamiltonella may contribute multiple cofactors and one essential amino acid, lysine. Homologs of the Bemisia metabolism genes of insect origin have also been implicated in essential amino acid synthesis in other sap-feeding insect hosts, indicative of parallel coevolution of shared metabolic pathways across multiple symbioses. Further metabolism genes coded in the Bemisia genome are of bacterial origin, but phylogenetically distinct from Portiera, Hamiltonella and horizontally transferred genes identified in other sap-feeding insects. Overall, 75% of the metabolism genes of bacterial origin are functionally unique to one symbiosis, indicating that the evolutionary history of metabolic integration in these symbioses is strongly contingent on the pattern of horizontally acquired genes. Our analysis, further, shows that bacteria with genomic decay enable host acquisition of complex metabolic pathways by multiple independent horizontal gene transfers from exogenous bacteria. Specifically, each horizontally acquired gene can function with other genes in the pathway coded by the symbiont, while facilitating the decay of the symbiont gene coding the same reaction.

The evolution of protein moonlighting: adaptive traps and promiscuity in the chaperonins

Moonlighting proteins exhibit functions that are alternative to their main role in the cell. Heat-shock proteins, also known as molecular chaperones, are now recognized for their wide range of activities in and/or outside the cell, being prominent examples of moonlighting proteins. Chaperonins are highly conserved molecular chaperones that fold other proteins into their native conformation allowing them to carry out essential functions in the cell. Activities alternative to folding have been reported for the chaperonin (Cpn) 60 protein. Preservation of various alternative functions in one protein conflicts with the optimization of each of the functions. What evolutionary mechanisms have allowed the persistence of moonlighting proteins, and in particular the chaperonins, remains a mystery. In the present article, I argue that mechanisms that increase the resistance of phenotypes to genetic and environmental perturbations enable the persistence of a reservoir of genetic variants, each potentially codifying for a distinct function. Gene duplication is one such mechanism that has characterized the expansion and has been concomitant with the emergence of novel functions in these protein families. Indeed, Cpn60 performs a large list of folding-independent functions, including roles in the transmission of viruses from insects to plants and stimulation of the immune system, among others. In addition to the innovation promoted by gene duplication, I discuss that the Cpn60 protein comprises a hidden amino acid combinatorial code that may well be responsible for its ability to develop novel functions while maintaining an optimized folding ability. The present review points to a complex model of evolution of protein moonlighting.

Fitness Trade-Offs Determine the Role of the Molecular Chaperonin GroEL in Buffering Mutations

Molecular chaperones fold many proteins and their mutated versions in a cell and can sometimes buffer the phenotypic effect of mutations that affect protein folding. Unanswered questions about this buffering include the nature of its mechanism, its influence on the genetic variation of a population, the fitness trade-offs constraining this mechanism, and its role in expediting evolution. Answering these questions is fundamental to understand the contribution of buffering to increase genetic variation and ecological diversification. Here, we performed experimental evolution, genome resequencing, and computational analyses to determine the trade-offs and evolutionary trajectories of Escherichia coli expressing high levels of the essential chaperonin GroEL. GroEL is abundantly present in bacteria, particularly in bacteria with large loads of deleterious mutations, suggesting its role in mutational buffering. We show that groEL overexpression is costly to large populations evolving in the laboratory, leading to groE expression decline within 66 generations. In contrast, populations evolving under the strong genetic drift characteristic of endosymbiotic bacteria avoid extinction or can be rescued in the presence of abundant GroEL. Genomes resequenced from cells evolved under strong genetic drift exhibited significantly higher tolerance to deleterious mutations at high GroEL levels than at native levels, revealing that GroEL is buffering mutations in these cells. GroEL buffered mutations in a highly diverse set of proteins that interact with the environment, including substrate and ion membrane transporters, hinting at its role in ecological diversification. Our results reveal the fitness trade-offs of mutational buffering and how genetic variation is maintained in populations.

Survival and innovation: The role of mutational robustness in evolution

Biological systems are resistant to perturbations caused by the environment and by the intrinsic noise of the system. Robustness to mutations is a particular aspect of robustness in which the phenotype is resistant to genotypic variation. Mutational robustness has been linked to the ability of the system to generate heritable genetic variation (a property known as evolvability). It is known that greater robustness leads to increased evolvability. Therefore, mechanisms that increase mutational robustness fuel evolvability. Two such mechanisms, molecular chaperones and gene duplication, have been credited with enormous importance in generating functional diversity through the increase of system's robustness to mutational insults. However, the way in which such mechanisms regulate robustness remains largely uncharacterized. In this review, I provide evidence in support of the role of molecular chaperones and gene duplication in innovation. Specifically, I present evidence that these mechanisms regulate robustness allowing unstable systems to survive long periods of time, and thus they provide opportunity for other mutations to compensate the destabilizing effects of functionally innovative mutations. The findings reported in this study set new questions with regards to the synergy between robustness mechanisms and how this synergy can alter the adaptive landscape of proteins. The ideas proposed in this article set the ground for future research in the understanding of the role of robustness in evolution.

The cost of antibiotic resistance depends on evolutionary history in Escherichia coli

BACKGROUND

The persistence of antibiotic resistance depends on the fitness effects of resistance elements in the absence of antibiotics. Recent work shows that the fitness effect of a given resistance mutation is influenced by other resistance mutations on the same genome. However, resistant bacteria acquire additional beneficial mutations during evolution in the absence of antibiotics that do not alter resistance directly but may modify the fitness effects of new resistance mutations.

RESULTS

We experimentally evolved rifampicin-resistant and sensitive Escherichia coli in a drug-free environment, before measuring the effects of new resistance elements on fitness in antibiotic-free conditions. Streptomycin-resistance mutations had small fitness effects in rifampicin-resistant genotypes that had adapted to antibiotic-free growth medium, compared to the same genotypes without adaptation. We observed a similar effect when resistance was encoded by a different mechanism and carried on a plasmid. Antibiotic-sensitive bacteria that adapted to the same conditions showed the same pattern for some resistance elements but not others.

CONCLUSIONS

Epistatic variation of costs of resistance can result from evolution in the absence of antibiotics, as well as the presence of other resistance mutations.

Coevolution analyses illuminate the dependencies between amino acid sites in the chaperonin system GroES-L

BACKGROUND

GroESL is a heat-shock protein ubiquitous in bacteria and eukaryotic organelles. This evolutionarily conserved protein is involved in the folding of a wide variety of other proteins in the cytosol, being essential to the cell. The folding activity proceeds through strong conformational changes mediated by the co-chaperonin GroES and ATP. Functions alternative to folding have been previously described for GroEL in different bacterial groups, supporting enormous functional and structural plasticity for this molecule and the existence of a hidden combinatorial code in the protein sequence enabling such functions. Describing this plasticity can shed light on the functional diversity of GroEL. We hypothesize that different overlapping sets of amino acids coevolve within GroEL, GroES and between both these proteins. Shifts in these coevolutionary relationships may inevitably lead to evolution of alternative functions.

RESULTS

We conducted the first coevolution analyses in an extensive bacterial phylogeny, revealing complex networks of evolutionary dependencies between residues in GroESL. These networks differed among bacterial groups and involved amino acid sites with functional importance and others with previously unsuspected functional potential. Coevolutionary networks formed statistically independent units among bacterial groups and map to structurally continuous regions in the protein, suggesting their functional link. Sites involved in coevolution fell within narrow structural regions, supporting dynamic combinatorial functional links involving similar protein domains. Moreover, coevolving sites within a bacterial group mapped to regions previously identified as involved in folding-unrelated functions, and thus, coevolution may mediate alternative functions.

CONCLUSIONS

Our results highlight the evolutionary plasticity of GroEL across the entire bacterial phylogeny. Evidence on the functional importance of coevolving sites illuminates the as yet unappreciated functional diversity of proteins.

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.

Genome fluctuations in cyanobacteria reflect evolutionary, developmental and adaptive traits

BACKGROUND

Cyanobacteria belong to an ancient group of photosynthetic prokaryotes with pronounced variations in their cellular differentiation strategies, physiological capacities and choice of habitat. Sequencing efforts have shown that genomes within this phylum are equally diverse in terms of size and protein-coding capacity. To increase our understanding of genomic changes in the lineage, the genomes of 58 contemporary cyanobacteria were analysed for shared and unique orthologs.

RESULTS

A total of 404 protein families, present in all cyanobacterial genomes, were identified. Two of these are unique to the phylum, corresponding to an AbrB family transcriptional regulator and a gene that escapes functional annotation although its genomic neighbourhood is conserved among the organisms examined. The evolution of cyanobacterial genome sizes involves a mix of gains and losses in the clade encompassing complex cyanobacteria, while a single event of reduction is evident in a clade dominated by unicellular cyanobacteria. Genome sizes and gene family copy numbers evolve at a higher rate in the former clade, and multi-copy genes were predominant in large genomes. Orthologs unique to cyanobacteria exhibiting specific characteristics, such as filament formation, heterocyst differentiation, diazotrophy and symbiotic competence, were also identified. An ancestral character reconstruction suggests that the most recent common ancestor of cyanobacteria had a genome size of approx. 4.5 Mbp and 1678 to 3291 protein-coding genes, 4%-6% of which are unique to cyanobacteria today.

CONCLUSIONS

The different rates of genome-size evolution and multi-copy gene abundance suggest two routes of genome development in the history of cyanobacteria. The expansion strategy is driven by gene-family enlargment and generates a broad adaptive potential; while the genome streamlining strategy imposes adaptations to highly specific niches, also reflected in their different functional capacities. A few genomes display extreme proliferation of non-coding nucleotides which is likely to be the result of initial expansion of genomes/gene copy number to gain adaptive potential, followed by a shift to a life-style in a highly specific niche (e.g. symbiosis). This transition results in redundancy of genes and gene families, leading to an increase in junk DNA and eventually to gene loss. A few orthologs can be correlated with specific phenotypes in cyanobacteria, such as filament formation and symbiotic competence; these constitute exciting exploratory targets.

Simultaneous analysis of foodborne pathogenic bacteria by an oligonucleotide microarray assay

A rapid and accurate method for simultaneous identification of foodborne infectious pathogens was developed based on oligonucleotide microarray technology. The proposed identification method is based on PCR amplification of the target region of the groEL genes with degenerate primers, followed by the PCR products hybridization with oligonucleotide probes specific for species. The groEL gene amplification products of seventeen species of pathogenic bacteria were hybridized to the oligonucleotide array. Hybridization results were analyzed with digoxigenin-linked enzyme reaction. Results indicated that fifteen species of pathogenic bacteria showed high sensitivity and specificity for the oligonucleotide array, while two other species gave cross-reaction with the E. coli. Our results suggested that microarray analysis of foodborne infectious pathogens might be very useful for simultaneous identification of bacterial pathogens. The oligonucleotide array can also be applied to samples collected in clinical settings of foodborne infections. The superiority of oligonucleotide array over other tests lies on its rapidity, accuracy and efficiency in the diagnosis, treatment and control of foodborne infections.

Protein homeostasis and the phenotypic manifestation of genetic diversity: principles and mechanisms

Changing a single nucleotide in a genome can have profound consequences under some conditions, but the same change can have no consequences under others. Indeed, organisms can be surprisingly robust to environmental and genetic perturbations. Yet, the mechanisms underlying such robustness are controversial. Moreover, how they might affect evolutionary change remains enigmatic. Here, we review the recently appreciated central role of protein homeostasis in buffering and potentiating genetic variation and discuss how these processes mediate the critical influence of the environment on the relationship between genotype and phenotype. Deciphering how robustness emerges from biological organization and the mechanisms by which it is overcome in changing environments will lead to a more complete understanding of both fundamental evolutionary processes and diverse human diseases.

Structural calibration of the rates of amino acid evolution in a search for Darwin in drifting biological systems

In the last two decades, many reports of proteins under positive selection have brought the neutral theory into question. However, the methods used to detect selection have ignored the evolvability of amino acids within proteins, which is fundamental to distinguishing positive selection from the relaxed constraints caused by genetic drift. Disentangling these two counterbalancing forces is essential to test the neutral theory. Here, we calibrate rates of amino acid divergence by using structural information from the full set of crystallized proteins in bacteria. In agreement with previous reports, we show that rates of amino acid evolution correlate negatively with the number of per-amino acid atomic interactions. Calibration of the rates of evolution allows identifying signatures of selection in biological systems that evolve under strong genetic drift, such as endosymbiotic bacteria. Application of this method identifies different rates and dynamics of evolution for highly connected amino acids in the structure compared with sparsely connected ones. We also unearth patterns of Darwinian selection in fundamental cellular proteins in endosymbiotic bacteria including the cochaperonin GroES, ribosomal proteins, proteins involved in cell cycle control, DNA-binding proteins, and proteins involved in DNA replication and repair. This is, to our knowledge, the first attempt to distinguish adaptive evolution from relaxed constraints in biological systems under genetic drift.

GroEL dependency affects codon usage--support for a critical role of misfolding in gene evolution

Integrating genome-scale sequence, expression, structural and protein interaction data from E. coli we establish an interaction between chaperone (GroEL) dependency and optimal codon usage.

Highly expressed sporadic substrates of GroEL employ more optimal codons than expected, show enrichment for optimal codons at structurally sensitive sites and greater conservation of codon optimality under conditions of relaxed purifying selection.

We suggest that highly expressed genes cannot routinely utilize GroEL for error control so that codon usage has evolved to provide complementary error limitation, whereas obligate GroEL substrates experience relaxed selection on codon usage.

Our results support a critical role of misfolding prevention in gene evolution.

Errors during gene expression are relatively commonplace, which has prompted speculations that many features of gene and genome anatomy and organization have evolved to reduce or mitigate such errors. One type of error that can be particularly costly occurs when the polypeptide chain that emerges from the ribosome fails to fold into its native structure. Some aberrantly folded proteins, exposing hydrophobic residues that would normally be buried, may begin to promiscuously interact with other proteins, become toxic to the cell and thus pose a substantial fitness concern (Gregersen et al, 2006).

In trans, molecular chaperones have long been recognized to play crucial roles in misfolding prevention and remedy. In cis, it has recently been suggested that the use of optimal codons limits mistranslation-induced protein misfolding (Drummond and Wilke, 2008). Evidence for the latter is centred on the argument that synonymous codons differ in their propensity to cause mistranslation. Translationally optimal codons, typically represented by more abundant cognate tRNAs (Duret, 2000), are thought less likely to cause ribosomal stalling and/or incorporation of the wrong amino acid.

Here, we suggest that the role, if any, of error limitation in cis can be revealed by studying its interaction with well-established error management systems in trans (chaperones). If codon usage does indeed play a tangible role in misfolding prevention, we would expect selection on codon identity to vary with the degree to which a protein can rely on other error control mechanisms, namely chaperones. We use the E. coli chaperonin GroEL as a model system to explore whether there is any interaction between optimal codon usage and chaperone dependency.

Kerner et al (2005) had previously determined GroEL substrates on a genome-wide scale. Based on enrichment in GroEL complexes the authors assigned ∼250 proteins to three classes reflecting GroEL dependency: class-I proteins, only a small fraction of which (<1%) associates with GroEL and which spontaneously regain some activity; class-II proteins, which only exhibit spontaneous refolding at more permissive temperatures and class-III proteins, which are obligate substrates of GroEL and largely fail to refold even under more benign conditions. Notably, although on average less abundant than class-I/II proteins (‘sporadic clients'), class-III proteins (‘obligate clients') occupy ∼80% of GroEL's capacity in vivo. Consequently, a higher proportion (∼100% versus ∼20% for class-II and ∼1% for class-I) of these proteins is routinely processed by the GroEL system.

We demonstrate that sporadic but not obligate clients of GroEL exhibit enhanced codon adaptation, carefully controlling for possible confounding factors, notably expression level and protein length (Figure 1). We also point out that genes that recently entered the E. coli genome via horizontal gene transfer will distort equilibrium analyses of codon usage in bacteria and should thus be routinely eliminated from analysis.

Building on earlier work by Zhou et al (2009), we further show that sporadic substrates are conspicuously enriched for optimal codons at structurally sensitive sites, consistent with more severe fitness implications of codon choice for these proteins.

Lastly, we reveal that codon optimality in sporadic clients is more highly conserved in S. dysenteriae. S. dysenteriae is an E. coli clone that has diverged relatively recently from the E. coli K12 strain and has adopted an intracellular lifestyle (Balbi et al, 2009). Concomitant with that lifestyle, Shigella has experienced a lower effective population size and therefore reduced efficiency of purifying selection. This has generated conditions where, overall, codon optimality has started to decay. However, when we followed the fate of ancestrally optimal codons at buried sites in the S. dysenteriae and E. coli K12 genomes, we found that a lower fraction of buried sites has lost codon optimality in sporadic substrates (Figure 4), again consistent with greater structural importance of codon choice in these substrates.

Based on the these findings, we suggest the following explanation: As mentioned above, class-III substrates are defined not only by GroEL being critical for proper folding, but also by occupying most of GroEL's capacity (∼80%). With a high proportion of class-III protein passaged through the GroEL system, mistranslation errors in these proteins weigh less severely as GroEL can remedy at least some misfolding that ensues. In contrast, class-I and II genes are more highly expressed and cannot routinely rely on GroEL to rectify folding errors. Yet class-I/II proteins are clearly liable to misfold as testified by their sporadic association with GroEL. We argue that augmenting GroEL's capacity to address the misfolding propensity of these genes would be prohibitively costly to the organism and that, as an alternative strategy, these genes employ optimal codons to reduce the rate of misfolding error.

Our findings (a) reveal a cis–trans interaction between codon usage and chaperones in providing an integrated error management system, (b) provide independent evidence for a role of misfolding in shaping gene evolution and (c) suggest that the burden of deleterious mutations in long-term bottlenecking populations like that of the insect endosymbiont Buchnera not only comprises unfavourable amino-acid (Moran, 1996) but also synonymous substitutions.

It has recently been suggested that the use of optimal codons limits mistranslation-induced protein misfolding, yet evidence for this remains largely circumstantial. In contrast, molecular chaperones have long been recognized to play crucial roles in misfolding prevention and remedy. We propose that putative error limitation in cis can be elucidated by examining the interaction between codon usage and chaperoning processes. Using Escherichia coli as a model system, we find that codon optimality covaries with dependency on the chaperonin GroEL. Sporadic but not obligate substrates of GroEL exhibit higher average codon adaptation and are conspicuously enriched for optimal codons at structurally sensitive sites. Further, codon optimality of sporadic clients is more conserved in the E. coli clone Shigella dysenteriae. We suggest that highly expressed genes cannot routinely use GroEL for error control so that codon usage has evolved to provide complementary error limitation. These findings provide independent evidence for a role of misfolding in shaping gene evolution and highlight the need to co-characterize adaptations in cis and trans to unravel the workings of integrated molecular systems.

Systemic analysis of the symbiotic function of Buchnera aphidicola, the primary endosymbiont of the pea aphid Acyrthosiphon pisum

Buchnera aphidicola is the primary obligate intracellular symbiont of most aphid species. B. aphidicola and aphids have been evolving in parallel since their association started, about 150 Myr ago. Both partners have lost their autonomy, and aphid diversification has been confined to smaller ecological niches by this co-evolution. B. aphidicola has undergone major genomic and biochemical changes as a result of adapting to intracellular life. Several genomes of B. aphidicola from different aphid species have been sequenced in the last decade, making it possible to carry out analyses and comparative studies using system-level in silico methods. This review attempts to provide a systemic description of the symbiotic function of aphid endosymbionts, particularly of B. aphidicola from the pea aphid Acyrthosiphon pisum, by analyzing their structural genomic properties, as well as their genetic and metabolic networks.

Genome-wide functional divergence after the symbiosis of proteobacteria with insects unraveled through a novel computational approach

Symbiosis has been among the most important evolutionary steps to generate biological complexity. The establishment of symbiosis required an intimate metabolic link between biological systems with different complexity levels. The strict endo-cellular symbiotic bacteria of insects are beautiful examples of the metabolic coupling between organisms belonging to different kingdoms, a eukaryote and a prokaryote. The host (eukaryote) provides the endosymbiont (prokaryote) with a stable cellular environment while the endosymbiont supplements the host's diet with essential metabolites. For such communication to take place, endosymbionts' genomes have suffered dramatic modifications and reconfigurations of proteins' functions. Two of the main modifications, loss of genes redundant for endosymbiotic bacteria or the host and bacterial genome streamlining, have been extensively studied. However, no studies have accounted for possible functional shifts in the endosymbiotic proteomes. Here, we develop a simple method to screen genomes for evidence of functional divergence between two species clusters, and we apply it to identify functional shifts in the endosymbiotic proteomes. Despite the strong effects of genetic drift in the endosymbiotic systems, we unexpectedly identified genes to be under stronger selective constraints in endosymbionts of aphids and ants than in their free-living bacterial relatives. These genes are directly involved in supplementing the host's diet with essential metabolites. A test of functional divergence supports a strong relationship between the endosymbiosis and the functional shifts of proteins involved in the metabolic communication with the insect host. The correlation between functional divergence in the endosymbiotic bacterium and the ecological requirements of the host uncovers their intimate biochemical and metabolic communication and provides insights on the role of symbiosis in generating species diversity.

Biological complexity has emerged on earth by the combination of living forms. This combination, called symbiosis, had to overcome the problems caused by the uncoupled metabolisms of the organisms involved. One way to do so was through the loss of genes that were no longer needed for the endosymbiont in the protected cellular environment provided by the host. Another step necessary to adjust both metabolisms was through the change in the function of bacterial proteins to perform new roles in the symbiotic system. In this article, we test such events in symbiotic systems involving an insect and a bacterium by developing a new and simple method to identify proteome-wide functional shifts. Our results show that most of the functional changes occurred at genes involved in metabolic communication with the host and are correlated with the host's ecological traits.

Selection for translational robustness in Buchnera aphidicola, endosymbiotic bacteria of aphids

Its strong intergenerational bottlenecks and effectively asexual reproduction have led Buchnera aphidicola, the endocellular symbiotic bacterium of aphids, to spectacular evolutionary and genomic changes in comparison with its free-living bacterial cousins. These changes summarize into high fixation rates of mildly deleterious destabilizing mutations. This predicts a sharp decline of its fitness and the consequent early demise of this endosymbiotic bacterium. Its survival for hundreds of millions of years casts doubt on genetic drift as the sole evolutionary force and seeks further explanation. We identify in Buchnera selection to increase the robustness of proteins to misfolding translation errors. Translational robustness varies between Buchnera lineages and protein functional categories. Metabolic proteins have been under selection for translational efficiency, whereas evolutionary rates of proteins involved in fundamental cellular processes have been largely determined by selection for translational robustness. We detect the strongest signal of translational robustness in B. aphidicola Cinara cedri with a very similar pattern to that inferred for the most common symbiotic ancestor of Buchnera lineages. This indicates that B. aphidicola Cinara cedri lineage may have probably reached the minimum evolutionary stable gene composition for endosymbiotic lifestyle. The evolutionary patterns from the comparative genomic analyses of these endosymbionts support a paradoxically complex dynamic for apparently simple genomes.

Theoretical approaches to the evolution of development and genetic architecture

Developmental evolutionary biology has, in the past decade, started to move beyond simply adapting traditional population and quantitative genetics models and has begun to develop mathematical approaches that are designed specifically to study the evolution of complex, nonadditive systems. This article first reviews some of these methods, discussing their strengths and shortcomings. The article then considers some of the principal questions to which these theoretical methods have been applied, including the evolution of canalization, modularity, and developmental associations between traits. I briefly discuss the kinds of data that could be used to test and apply the theories, as well as some consequences for other approaches to phenotypic evolution of discoveries from theoretical studies of developmental evolution.

An insecticidal GroEL protein with chitin binding activity from Xenorhabdus nematophila

Xenorhabdus nematophila secretes insecticidal proteins to kill its larval prey. We have isolated an approximately 58-kDa GroEL homolog, secreted in the culture medium through outer membrane vesicles. The protein was orally insecticidal to the major crop pest Helicoverpa armigera with an LC50 of approximately 3.6 microg/g diet. For optimal insecticidal activity all three domains of the protein, apical, intermediate, and equatorial, were necessary. The apical domain alone was able to bind to the larval gut membranes and manifest low level insecticidal activity. At equimolar concentrations, the apical domain contained approximately one-third and the apical-intermediate domain approximately one-half bioactivity of that of the full-length protein. Interaction of the protein with the larval gut membrane was specifically inhibited by N-acetylglucosamine and chito-oligosaccharides. Treatment of the larval gut membranes with chitinase abolished protein binding. Based on the three-dimensional structural model, mutational analysis demonstrated that surface-exposed residues Thr-347 and Ser-356 in the apical domain were crucial for both binding to the gut epithelium and insecticidal activity. Double mutant T347A,S356A was 80% less toxic (p < 0.001) than the wild type protein. The GroEL homolog showed alpha-chitin binding activity with Kd approximately 0.64 microm and Bmax approximately 4.68 micromol/g chitin. The variation in chitin binding activity of the mutant proteins was in good agreement with membrane binding characteristics and insecticidal activity. The less toxic double mutant XnGroEL showed an approximately 8-fold increase of Kd in chitin binding assay. Our results demonstrate that X. nematophila secretes an insecticidal GroEL protein with chitin binding activity.

Population structure, levels of selection, and the evolution of intracellular symbionts

Many obligately intracellular symbionts exhibit a characteristic set of genetic changes that include an increase in substitution rates, loss of many genes, and apparent destabilization of many proteins and structural RNAs. Authors have suggested that these changes are due to increased mutation rates, or, more commonly, decreased effective population size due to population bottlenecks at the symbiont or, perhaps, host level. I propose that the increase in substitution rates and accumulation of deleterious mutations is a consequence of the population structure imposed on the endosymbionts by strict host association, loss of horizontal transmission and potentially conflicting levels of selection. I analyze a population genetic model of endosymbiont evolution, and demonstrate that substitution rates will increase, and the effect of those substitutions on endosymbiont fitness will become more deleterious as horizontal transmission among hosts decreases. Additionally, I find that there is a critical level of horizontal transmission below which natural selection cannot effectively purge deleterious mutations, leading to an expected loss of fitness over time. This critical level varies across loci with the degree of correlation between host and endosymbiont fitness, and may help explain differential retention and loss of certain genes.

Proteomic profiling of aphid Macrosiphum euphorbiae responses to host-plant-mediated stress induced by defoliation and water deficit

Abiotic and biotic host-plant stress, such as desiccation and herbivory, may strongly affect sap-sucking insects such as aphids via changes in plant chemicals of insect nutritional or plant defensive value. Here, we examined (i) water deprivation and (ii) defoliation by the beetle Leptinotarsa decemlineata as stresses indirectly affecting the aphid Macrosiphum euphorbiae via its host plant Solanum tuberosum. For plant-induced stress, aphids were reared on healthy vs. continuously stressed potato for 14 days (no watering; defoliation maintained at approximately 40%). Aphid performance under stress was correlated with metabolic responses monitored by profiling of the aphid proteome. M. euphorbiae was strongly affected by water stress, as adult survival, total aphid number and biomass were reduced by 67%, 64%, and 79%, respectively. Aphids performed normally on defoliated potato, indicating that they were unaffected or able to compensate any stress induced by plant defoliation. Stressed aphid proteomes revealed 419-453 protein spots, including 27 that were modulated specifically or jointly under each kind of host-plant stress. Reduced aphid fitness on water-stressed plants mostly correlated with modulation of proteins involved in energy metabolism, apparently to conserve energy in order to prioritize survival. Despite normal performance, several aphid proteins that are known to be implicated in cell communication were modulated on defoliated plants, possibly suggesting modified aphid behaviour. The GroEL protein (or symbionin) of the endosymbiont Buchnera aphidicola was predominant under all conditions in M. euphorbiae. Its expression level was not significantly affected by aphid host-plant stresses, which is consistent with the high priority of symbiosis in stressed aphids.

Structural divergence and adaptive evolution in mammalian cytochromes P450 2C

Cytochromes P450 (CYPs) comprise a superfamily of enzymes involved in various physiological functions, including the metabolism of drugs and carcinogenic compounds present in food, making them of great importance for human health. The possibility that CYPs could be broadening or changing substrate specificity in accordance to the high diversity of xenobiotics compounds environmentally available suggests that their metabolic function could be under adaptive evolution. We evaluated the existence of functional divergence and signatures of selection on mammalian genes from the drug-metabolizing CYP2 family. Thirteen of the sites found to be functionally divergent and the eight found to be under strong positive selection occurred in important functional domains, namely on the substrate entrance channel and within the active site. Our results provide insight into CYPs evolution and the role of molecular adaptation in enzyme substrate-specificity diversification.

Genotypic diversity of Haemophilus parasuis field strains

Haemophilus parasuis is the cause of Glässer's disease and other clinical disorders in pigs. It can also be isolated from the upper respiratory tracts of healthy pigs, and isolates can have significant differences in virulence. In this work, a partial sequence from the 60-kDa heat shock protein (Hsp60) gene was assessed as an epidemiological marker. We analyzed partial sequences of hsp60 and 16S rRNA genes from 103 strains of H. parasuis and other related species to obtain a better classification of the strains and examine the correlation with virulence. The results were compared with those obtained by enterobacterial repetitive intergenic consensus PCR. Our results showed that hsp60 is a reliable marker for epidemiological studies of H. parasuis and that the analysis of its sequence is a better approach than fingerprinting methods. Furthermore, the analysis of the hsp60 and 16S rRNA gene sequences revealed the presence of a separate lineage of virulent strains and indicated the occurrence of lateral gene transfer among H. parasuis and Actinobacillus strains.

GRAST: a new way of genome reduction analysis using comparative genomics

MOTIVATION

Establishment of intra-cellular life involved a profound re-configuration of the genetic characteristics of bacteria, including genome reduction and rearrangements. Understanding the mechanisms underlying these phenomena will shed light on the genome rearrangements essential for the development of an intra-cellular lifestyle. Comparison of genomes with differences in their sizes poses statistical as well as computational problems. Little efforts have been made to develop flexible computational tools with which to analyse genome reduction and rearrangements.

RESULTS

Investigation of genome reduction and rearrangements in endosymbionts using a novel computational tool (GRAST) identified gathering of genes with similar functions. Conserved clusters of functionally related genes (CGSCs) were detected. Heterogeneous gene and gene cluster non-functionalization/loss are identified between genome regions, functional gene categories and during evolution. Results show that gene non-functionalisation has accelerated during the last 50 MY of Buchnera's evolution while CGSCs have been static.

A novel method for detecting intramolecular coevolution: adding a further dimension to selective constraints analyses

Protein evolution depends on intramolecular coevolutionary networks whose complexity is proportional to the underlying functional and structural interactions among sites. Here we present a novel approach that vastly improves the sensitivity of previous methods for detecting coevolution through a weighted comparison of divergence between amino acid sites. The analysis of the HIV-1 Gag protein detected convergent adaptive coevolutionary events responsible for the selective variability emerging between subtypes. Coevolution analysis and functional data for heat-shock proteins, Hsp90 and GroEL, highlight that almost all detected coevolving sites are functionally or structurally important. The results support previous suggestions pinpointing the complex interdomain functional interactions within these proteins and we propose new amino acid sites as important for interdomain functional communication. Three-dimensional information sheds light on the functional and structural constraints governing the coevolution between sites. Our covariation analyses propose two types of coevolving sites in agreement with previous reports: pairs of sites spatially proximal, where compensatory mutations could maintain the local structure stability, and clusters of distant sites located in functional domains, suggesting a functional dependency between them. All sites detected under adaptive evolution in these proteins belong to coevolution groups, further underlining the importance of testing for coevolution in selective constraints analyses.

The roles of positive and negative selection in the molecular evolution of insect endosymbionts

The evolutionary rate acceleration observed in most endosymbiotic bacteria may be explained by higher mutation rates, changes in selective pressure, and increased fixation of deleterious mutations by genetic drift. Here, we explore the forces influencing molecular evolution in Blochmannia, an obligate endosymbiont of Camponotus and related ant genera. Our goals were to compare rates of sequence evolution in Blochmannia with related bacteria, to explore variation in the strength and efficacy of negative (purifying) selection, and to evaluate the effect of positive selection. For six Blochmannia pairs, plus Buchnera and related enterobacteria, estimates of sequence divergence at four genes confirm faster rates of synonymous evolution in the ant mutualist. This conclusion is based on higher dS between Blochmannia lineages despite their more recent divergence. Likewise, generally higher dN in Blochmannia indicates faster rates of nonsynonymous substitution in this group. One exception is the groEL gene, for which lower dN and dN/dS compared to Buchnera indicate exceptionally strong negative selection in Blochmannia. In addition, we explored evidence for positive selection in Blochmannia using both site-and lineage-based maximum likelihood models. These approaches confirmed heterogeneity of dN/dS among codon sites and revealed significant variation in dN/dS across Blochmannia lineages for three genes. Lineage variation affected genes independently, with no evidence of parallel changes in dN/dS across genes along a given branch. Our data also reveal instances of dN/dS greater than one; however, we do not interpret these large dN/dS ratios as evidence for positive selection. In sum, while drift may contribute to an overall rate acceleration at nonsynonymous sites in Blochmannia, variable selective pressures best explain the apparent gene-specific changes in dN/dS across lineages of this ant mutualist. In the course of this study, we reanalyzed variation at Buchnera groEL and found no evidence of positive selection that was previously reported.

Adaptive evolution in GroEL from distantly related endosymbiotic bacteria of insects

Many symbioses between bacteria and insects resulted from ancient infections followed by strict vertical transmission within host lineages. The strong bottlenecks under which this transmission occurs promote the neutral fixation of slightly deleterious mutations by genetic drift. As predicted by Muller's ratchet, this fixation will drive endosymbiotic bacteria through an irreversible dynamics of fitness loss. The chaperonin GroEL has been proposed as a compensatory mechanism whereby endosymbiotic bacteria of aphids persist. Here, we show that endosymbiotic bacteria of insects from two phylogenetically very distant bacterial phyla have fixed amino acid substitutions by positive selection in functionally important GroEL regions involved in either GroES/peptide binding or in the en bloc movement of the GroEL apical domain. These results, together with the high levels of constitutive expression of GroEL in these endosymbionts, provide valuable insights into the evolution of a molecular mechanism responsible for the maintenance of the symbiotic lifestyle.

Bacterial endosymbionts of insects: insights from comparative genomics

The development of molecular techniques for the study of uncultured bacteria allowed the extensive study of the widespread association between insects and intracellular symbiotic bacteria. Most of the bacterial endosymbionts involved in such associations are gamma-proteobacteria, closely related to Escherichia coli. In recent years, five genomes from insect endosymbionts have been sequenced, allowing the performance of extensive genome comparative analysis that, as a complement of phylogenetic studies, and analysis on individual genes, can help to understand the different traits of this particular association, including how the symbiotic process is established, the explanation of the special features of these microbial genomes, the bases of this intimate association and the possible future that awaits the endosymbionts with extremely reduced genomes.

Evolution of the leucine gene cluster in Buchnera aphidicola: insights from chromosomal versions of the cluster

In Buchnera aphidicola strains associated with the aphid subfamilies Thelaxinae, Lachninae, Pterocommatinae, and Aphidinae, the four leucine genes (leuA, -B, -C, and -D) are located on a plasmid. However, these genes are located on the main chromosome in B. aphidicola strains associated with the subfamilies Pemphiginae and Chaitophorinae. The sequence of the chromosomal fragment containing the leucine cluster and flanking genes has different positions in the chromosome in B. aphidicola strains associated with three tribes of the subfamily Pemphiginae and one tribe of the subfamily Chaitophorinae. Due to the extreme gene order conservation of the B. aphidicola genomes, the variability in the position of the leucine cluster in the chromosome may be interpreted as resulting from independent insertions from an ancestral plasmid-borne leucine gene. These findings do not support a chromosomal origin for the leucine genes in the ancestral B. aphidicola and do support a back transfer evolutionary scenario from a plasmid to the main chromosome.