Synonymous genes explore different evolutionary landscapes

L’impact des mutations neutres sur l’évolvabilité et l’évolution des génomes

Les mutations bénéfiques à forts effets sont rares et les mutations délétères sont éliminées par la sélection naturelle. La majorité des mutations qui s’accumulent dans les génomes ont donc des effets sélectifs très faibles, voire nuls ; elles sont alors appelées mutations neutres. Au cours des deux dernières décennies, il a été montré que les mutations, même en l’absence d’effet sur la valeur sélective des organismes, affectent leur évolvabilité, en donnant accès à de nouveaux phénotypes par le biais de mutations apparaissant ultérieurement, et qui n’auraient pas été disponibles autrement. En plus de cet effet, de nombreuses mutations neutres – indépendamment de leurs effets sélectifs – peuvent affecter la mutabilité de séquences d’ADN voisines, et moduler l’efficacité de la recombinaison homologue. De telles mutations ne modifient pas le spectre des phénotypes accessibles, mais plutôt la vitesse à laquelle de nouveaux phénotypes seront produits, un processus qui a des conséquences à long terme mais aussi potentiellement à court terme, en lien avec l’émergence de cancers.

Increased RNA virus population diversity improves adaptability

The replication machinery of most RNA viruses lacks proofreading mechanisms. As a result, RNA virus populations harbor a large amount of genetic diversity that confers them the ability to rapidly adapt to changes in their environment. In this work, we investigate whether further increasing the initial population diversity of a model RNA virus can improve adaptation to a single selection pressure, thermal inactivation. For this, we experimentally increased the diversity of coxsackievirus B3 (CVB3) populations across the capsid region. We then compared the ability of these high diversity CVB3 populations to achieve resistance to thermal inactivation relative to standard CVB3 populations in an experimental evolution setting. We find that viral populations with high diversity are better able to achieve resistance to thermal inactivation at both the temperature employed during experimental evolution as well as at a more extreme temperature. Moreover, we identify mutations in the CVB3 capsid that confer resistance to thermal inactivation, finding significant mutational epistasis. Our results indicate that even naturally diverse RNA virus populations can benefit from experimental augmentation of population diversity for optimal adaptation and support the use of such viral populations in directed evolution efforts that aim to select viruses with desired characteristics.

Primary and promiscuous functions coexist during evolutionary innovation through whole protein domain acquisitions

Molecular examples of evolutionary innovation are scarce and generally involve point mutations. Innovation can occur through larger rearrangements, but here experimental data is extremely limited. Integron integrases innovated from double-strand- toward single-strand-DNA recombination through the acquisition of the I2 α-helix. To investigate how this transition was possible, we have evolved integrase IntI1 to what should correspond to an early innovation state by selecting for its ancestral activity. Using synonymous alleles to enlarge sequence space exploration, we have retrieved 13 mutations affecting both I2 and the multimerization domains of IntI1. We circumvented epistasis constraints among them using a combinatorial library that revealed their individual and collective fitness effects. We obtained up to 10⁴-fold increases in ancestral activity with various asymmetrical trade-offs in single-strand-DNA recombination. We show that high levels of primary and promiscuous functions could have initially coexisted following I2 acquisition, paving the way for a gradual evolution toward innovation.

Large-scale in silico mutagenesis experiments reveal optimization of genetic code and codon usage for protein mutational robustness

Background

How, and the extent to which, evolution acts on DNA and protein sequences to ensure mutational robustness and evolvability is a long-standing open question in the field of molecular evolution. We addressed this issue through the first structurome-scale computational investigation, in which we estimated the change in folding free energy upon all possible single-site mutations introduced in more than 20,000 protein structures, as well as through available experimental stability and fitness data.

Results

At the amino acid level, we found the protein surface to be more robust against random mutations than the core, this difference being stronger for small proteins. The destabilizing and neutral mutations are more numerous in the core and on the surface, respectively, whereas the stabilizing mutations are about 4% in both regions. At the genetic code level, we observed smallest destabilization for mutations that are due to substitutions of base III in the codon, followed by base I, bases I+III, base II, and other multiple base substitutions. This ranking highly anticorrelates with the codon-anticodon mispairing frequency in the translation process. This suggests that the standard genetic code is optimized to limit the impact of random mutations, but even more so to limit translation errors. At the codon level, both the codon usage and the usage bias appear to optimize mutational robustness and translation accuracy, especially for surface residues.

Conclusion

Our results highlight the non-universality of mutational robustness and its multiscale dependence on protein features, the structure of the genetic code, and the codon usage. Our analyses and approach are strongly supported by available experimental mutagenesis data.

The Impact of Neutral Mutations on Genome Evolvability

Beneficial mutations are rare and deleterious mutations are purged by natural selection. As a result, the vast majority of mutations that accumulate in genomes belong to the class of neutral mutations. Over the last two decades, neutral mutations, despite their null effect on fitness, have been shown to affect evolvability by providing access to new phenotypes through subsequent mutations that would not have been available otherwise. Here we propose that in addition, many mutations - independent of their selective effects - can affect the mutability of neighboring DNA sequences and modulate the efficacy of homologous recombination. Such mutations do not change the spectrum of accessible phenotypes, but rather the rate at which new phenotypes will be produced. Therefore, neutral mutations that accumulate in genomes have an important long-term impact on the evolutionary fate of genomes.

HIV-1 Protease Evolvability Is Affected by Synonymous Nucleotide Recoding

One unexplored aspect of HIV-1 genetic architecture is how codon choice influences population diversity and evolvability. Here we compared the levels of development of HIV-1 resistance to protease inhibitors (PIs) between wild-type (WT) virus and a synthetic virus (MAX) carrying a codon-pair-reengineered protease sequence including 38 (13%) synonymous mutations. The WT and MAX viruses showed indistinguishable replication in MT-4 cells or peripheral blood mononuclear cells (PBMCs). Both viruses were subjected to serial passages in MT-4 cells, with selective pressure from the PIs atazanavir (ATV) and darunavir (DRV). After 32 successive passages, both the WT and MAX viruses developed phenotypic resistance to PIs (50% inhibitory concentrations [IC₅₀s] of 14.6 ± 5.3 and 21.2 ± 9 nM, respectively, for ATV and 5.9 ± 1.0 and 9.3 ± 1.9, respectively, for DRV). Ultradeep sequence clonal analysis revealed that both viruses harbored previously described mutations conferring resistance to ATV and DRV. However, the WT and MAX virus proteases showed different resistance variant repertoires, with the G16E and V77I substitutions observed only in the WT and the L33F, S37P, G48L, Q58E/K, and L89I substitutions detected only in the MAX virus. Remarkably, the G48L and L89I substitutions are rarely found in vivo in PI-treated patients. The MAX virus showed significantly higher nucleotide and amino acid diversity of the propagated viruses with and without PIs (P < 0.0001), suggesting a higher selective pressure for change in this recoded virus. Our results indicate that the HIV-1 protease position in sequence space delineates the evolution of its mutant spectrum. Nevertheless, the investigated synonymously recoded variant showed mutational robustness and evolvability similar to those of the WT virus.IMPORTANCE Large-scale synonymous recoding of virus genomes is a new tool for exploring various aspects of virus biology. Synonymous virus genome recoding can be used to investigate how a virus's position in sequence space defines its mutant spectrum, evolutionary trajectory, and pathogenesis. In this study, we evaluated how synonymous recoding of the human immunodeficiency virus type 1 (HIV-1) protease affects the development of protease inhibitor (PI) resistance. HIV-1 protease is a main target of current antiretroviral therapies. Our present results demonstrate that the wild-type (WT) virus and a virus with recoded protease exhibited different patterns of resistance mutations after PI treatment. Nevertheless, the developed PI resistance phenotypes were indistinguishable between the recoded virus and the WT virus, suggesting that the HIV-1 strain with synonymously recoded protease and the WT virus are equally robust and evolvable.

Polymorphism of antibiotic-inactivating enzyme driven by ecology expands the environmental resistome

The environmental resistome has been recognized as the origin and reservoir of antibiotic resistance genes and considered to be dynamic and ever expanding. In this study, a targeted gene sequencing approach revealed that the polymorphic diversity of the aminoglycoside-inactivating enzyme AAC(6')-Ib was ecological niche-specific. AAC(6')-Ib-cr, previously known as a clinical variant, was prevalent in various soils and the intestines of chickens and humans, suggesting that this variant might not have arisen from adaptive mutations in the clinic but instead originated from the environment. Furthermore, ecologically dominant polymorphic variants of AAC(6')-Ib were characterized and found to display different substrate specificities for quinolones and aminoglycosides, conferring the altered resistance spectra. Interestingly, a novel variant with the D179Y substitution showed an extended resistance spectrum to the recently developed fluoroquinolone gemifloxacin. Our results suggest that soil and animal microbiomes could be major reservoirs of antibiotic resistance; polymorphic diversity expands the antibiotic resistome in the environment, resulting in the potential emergence of novel resistance.

The Integron: Adaptation On Demand

The integron is a powerful system which, by capturing, stockpiling, and rearranging new functions carried by gene encoding cassettes, confers upon bacteria a rapid adaptation capability in changing environments. Chromosomally located integrons (CI) have been identified in a large number of environmental Gram-negative bacteria. Integron evolutionary history suggests that these sedentary CIs acquired mobility among bacterial species through their association with transposable elements and conjugative plasmids. As a result of massive antibiotic use, these so-called mobile integrons are now widespread in clinically relevant bacteria and are considered to be the principal agent in the emergence and rise of antibiotic multiresistance in Gram-negative bacteria. Cassette rearrangements are catalyzed by the integron integrase, a site-specific tyrosine recombinase. Central to these reactions is the single-stranded DNA nature of one of the recombination partners, the attC site. This makes the integron a unique recombination system. This review describes the current knowledge on this atypical recombination mechanism, its implications in the reactions involving the different types of sites, attC and attI, and focuses on the tight regulation exerted by the host on integron activity through the control of attC site folding. Furthermore, cassette and integrase expression are also highly controlled by host regulatory networks and the bacterial stress (SOS) response. These intimate connections to the host make the integron a genetically stable and efficient system, granting the bacteria a low cost, highly adaptive evolution potential "on demand".

Molecular characterization influencing metal resistance in the Cupriavidus/Ralstonia genomes

Our environment is stressed with a load of heavy and toxic metals. Microbes, abundant in our environment, are found to adapt well to this metal-stressed condition. A comparative study among five Cupriavidus/Ralstonia genomes can offer a better perception of their evolutionary mechanisms to adapt to these conditions. We have studied codon usage among 1051 genes common to all these organisms and identified 15 optimal codons frequently used in highly expressed genes present within 1051 genes. We found the core genes of Cupriavidus metallidurans CH34 have a different optimal codon choice for arginine, glycine and alanine in comparison with the other four bacteria. We also found that the synonymous codon usage bias within these 1051 core genes is highly correlated with their gene expression. This supports that translational selection drives synonymous codon usage in the core genes of these genomes. Synonymous codon usage is highly conserved in the core genes of these five genomes. The only exception among them is C. metallidurans CH34. This genomewide shift in synonymous codon choice in C. metallidurans CH34 may have taken place due to the insertion of new genes in its genomes facilitating them to survive in heavy metal containing environment and the co-evolution of the other genes in its genome to achieve a balance in gene expression. Structural studies indicated the presence of a longer N-terminal region containing a copper-binding domain in the cupC proteins of C. metallidurans CH3 that helps it to attain higher binding efficacy with copper in comparison with its orthologs.

Encoding asymmetry of the N-glycosylation motif facilitates glycoprotein evolution

Protein N-glycosylation is found in all domains of life and has a conserved role in glycoprotein folding and stability. In animals, glycoproteins transit through the Golgi where the N-glycans are trimmed and rebuilt with sequences that bind lectins, an innovation that greatly increases structural diversity and redundancy of glycoprotein-lectin interaction at the cell surface. Here we ask whether the natural tension between increasing diversity (glycan-protein interactions) and site multiplicity (backup and status quo) might be revealed by a phylogenic examination of glycoproteins and NXS/T(X≠P) N-glycosylation sites. Site loss is more likely by mutation at Asn encoded by two adenosine (A)-rich codons, while site gain is more probable by generating Ser or Thr downstream of an existing Asn. Thus mutations produce sites at novel positions more frequently than the reversal of recently lost sites, and therefore more paths though sequence space are made available to natural selection. An intra-species comparison of secretory and cytosolic proteins revealed a departure from equilibrium in sequences one-mutation-away from NXS/T and in (A) content, indicating strong selective pressures and exploration of N-glycosylation positions during vertebrate evolution. Furthermore, secretory proteins have evolved at rates proportional to N-glycosylation site number, indicating adaptive interactions between the N-glycans and underlying protein. Given the topology of the genetic code, mutation of (A) is more often nonsynonomous, and Lys, another target of many PTMs, is also encoded by two (A)-rich codons. An examination of acetyl-Lys sites in proteins indicated similar evolutionary dynamics, consistent with asymmetry of the target and recognition portions of modified sites. Our results suggest that encoding asymmetry is an ancient mechanism of evolvability that increases diversity and experimentation with PTM site positions. Strong selective pressures on PTMs may have contributed to the A+T→G+C shift in genome-wide nucleotide composition during metazoan radiation.

D-Tailor: automated analysis and design of DNA sequences

Motivation: Current advances in DNA synthesis, cloning and sequencing technologies afford high-throughput implementation of artificial sequences into living cells. However, flexible computational tools for multi-objective sequence design are lacking, limiting the potential of these technologies.

Results: We developed DNA-Tailor (D-Tailor), a fully extendable software framework, for property-based design of synthetic DNA sequences. D-Tailor permits the seamless integration of multiple sequence analysis tools into a generic Monte Carlo simulation that evolves sequences toward any combination of rationally defined properties. As proof of principle, we show that D-Tailor is capable of designing sequence libraries comprising all possible combinations among three different sequence properties influencing translation efficiency in Escherichia coli. The capacity to design artificial sequences that systematically sample any given parameter space should support the implementation of more rigorous experimental designs.

Availability: Source code is available for download at https://sourceforge.net/projects/dtailor/

Contact:aparkin@lbl.gov or cambray.guillaume@gmail.com

Supplementary information:Supplementary data are available at Bioinformatics online (D-Tailor Tutorial).

The genetic code constrains yet facilitates Darwinian evolution

An important goal of evolutionary biology is to understand the constraints that shape the dynamics and outcomes of evolution. Here, we address the extent to which the structure of the standard genetic code constrains evolution by analyzing adaptive mutations of the antibiotic resistance gene TEM-1 β-lactamase and the fitness distribution of codon substitutions in two influenza hemagglutinin inhibitor genes. We find that the architecture of the genetic code significantly constrains the adaptive exploration of sequence space. However, the constraints endow the code with two advantages: the ability to restrict access to amino acid mutations with a strong negative effect and, most remarkably, the ability to enrich for adaptive mutations. Our findings support the hypothesis that the standard genetic code was shaped by selective pressure to minimize the deleterious effects of mutation yet facilitate the evolution of proteins through imposing an adaptive mutation bias.

Rise and dissemination of aminoglycoside resistance: the aac(6')-Ib paradigm

Enzymatic modification is a prevalent mechanism by which bacteria defeat the action of antibiotics. Aminoglycosides are often inactivated by aminoglycoside modifying enzymes encoded by genes present in the chromosome, plasmids, and other genetic elements. The AAC(6′)-Ib (aminoglycoside 6′-N-acetyltransferase type Ib) is an enzyme of clinical importance found in a wide variety of gram-negative pathogens. The AAC(6′)-Ib enzyme is of interest not only because of his ubiquity but also because of other characteristics, it presents significant microheterogeneity at the N-termini and the aac(6′)-Ib gene is often present in integrons, transposons, plasmids, genomic islands, and other genetic structures. Excluding the highly heterogeneous N-termini, there are 45 non-identical AAC(6′)-Ib related entries in the NCBI database, 32 of which have identical name in spite of not having identical amino acid sequence. While some variants conserved similar properties, others show dramatic differences in specificity, including the case of AAC(6′)-Ib-cr that mediates acetylation of ciprofloxacin representing a rare case where a resistance enzyme acquires the ability to utilize an antibiotic of a different class as substrate. Efforts to utilize antisense technologies to turn off expression of the gene or to identify enzymatic inhibitors to induce phenotypic conversion to susceptibility are under way.

Codon usage determines the mutational robustness, evolutionary capacity, and virulence of an RNA virus

RNA viruses exist as dynamic and diverse populations shaped by constant mutation and selection. Yet little is known about how the mutant spectrum contributes to virus evolvability and pathogenesis. Because several codon choices are available for a given amino acid, a central question concerns whether viral sequences have evolved to optimize not only the protein coding consensus, but also the DNA/RNA sequences accessible through mutation. Here we directly test this hypothesis by comparing wild-type poliovirus to synthetic viruses carrying re-engineered capsid sequences with hundreds of synonymous mutations. Strikingly, such rewiring of the population's mutant network reduced its robustness and attenuated the virus in an animal model of infection. We conclude that the position of a virus in sequence space defines its mutant spectrum, evolutionary trajectory, and pathogenicity. This organizing principle for RNA virus populations confers tolerance to mutations and facilitates replication and spread within the dynamic host environment.

Directed evolution of a model primordial enzyme provides insights into the development of the genetic code

The contemporary proteinogenic repertoire contains 20 amino acids with diverse functional groups and side chain geometries. Primordial proteins, in contrast, were presumably constructed from a subset of these building blocks. Subsequent expansion of the proteinogenic alphabet would have enhanced their capabilities, fostering the metabolic prowess and organismal fitness of early living systems. While the addition of amino acids bearing innovative functional groups directly enhances the chemical repertoire of proteomes, the inclusion of chemically redundant monomers is difficult to rationalize. Here, we studied how a simplified chorismate mutase evolves upon expanding its amino acid alphabet from nine to potentially 20 letters. Continuous evolution provided an enhanced enzyme variant that has only two point mutations, both of which extend the alphabet and jointly improve protein stability by >4 kcal/mol and catalytic activity tenfold. The same, seemingly innocuous substitutions (Ile→Thr, Leu→Val) occurred in several independent evolutionary trajectories. The increase in fitness they confer indicates that building blocks with very similar side chain structures are highly beneficial for fine-tuning protein structure and function.

Proteins are linear polymers of a set of typically 20 different amino acid building blocks. The amino acid sequence—encoded by a genetic template—directs the folding of newly synthesized proteins into compact 3D structures and dictates the function of the protein product. Monomers containing distinct physico-chemical properties and geometries allow the formation of highly sophisticated architectures, and diverse functional groups enable enzymes to catalyze a plethora of chemical transformations. Nevertheless, the biochemical rationale for the exact composition (and particularly the redundancy) of the proteinogenic amino acid alphabet, which contains multiple building blocks that are chemically similar, remains enigmatic. By subjecting a simplified enzyme—constructed from only nine different amino acids—to directed evolution, we were able to investigate the impact of amino acid diversity on protein function. The most prolific variant selected in the course of the experiments expanded its amino acid alphabet, albeit through two surprisingly subtle mutations (isoleucine to threonine and leucine to valine). The mutations improve both stability and catalytic activity of the enzyme, thereby demonstrating that the presence of structurally similar amino acids specified by the genetic code is highly beneficial for protein fitness.

Directional evolution of Chlamydia trachomatis towards niche-specific adaptation

On behalf of the host-pathogen "arms race," a cutting-edge approach for elucidating genotype-phenotype relationships relies on the identification of positively selected loci involved in pathoadaptation. We studied the obligate intracellular bacterium Chlamydia trachomatis, for which same-species strains display a nearly identical core and pan genome, while presenting a wide range of tissue tropism and ecological success. We sought to evaluate the evolutionary patterns underlying species separation (divergence) and C. trachomatis serovar radiation (polymorphism) and to establish genotype-phenotype associations. By analyzing 60 Chlamydia strains, we detected traces of Muller's ratchet as a result of speciation and identified positively selected genes and codons hypothetically involved in the infection of different human cell types (e.g., columnar epithelial cells of ocular or genital mucosae and mononuclear phagocytes) and also events likely driving pathogenic and ecological success dissimilarities. In general, these genes code for proteins involved in immune response elicitation, proteolysis, and the subversion of host-cell functions, and also for proteins with unknown function(s). Several genes are potentially involved in more than one adaptive process, suggesting multiple functions or a distinct modus operandi for a specific function, and thus should be considered as crucial research targets. In addition, six of the nine genes encoding the putative antigen/adhesin polymorphic membrane proteins seem to be under positive selection along specific serovars, which sustains an essential biological role of this extra-large paralogue family in chlamydial pathobiology. This study provides insight into how evolutionary inferences illuminate ecological processes such as adaptation to different niches, pathogenicity, or ecological success driven by arms races.

Epistasis increases the rate of conditionally neutral substitution in an adapting population

Kimura observed that the rate of neutral substitution should equal the neutral mutation rate. This classic result is central to our understanding of molecular evolution, and it continues to influence phylogenetics, genomics, and the interpretation of evolution experiments. By demonstrating that neutral mutations substitute at a rate independent of population size and selection at linked sites, Kimura provided an influential justification for the idea of a molecular clock and emphasized the importance of genetic drift in shaping molecular evolution. But when epistasis among sites is common, as numerous empirical studies suggest, do neutral mutations substitute according to Kimura's expectation? Here we study simulated, asexual populations of RNA molecules, and we observe that conditionally neutral mutations--i.e., mutations that do not alter the fitness of the individual in which they arise, but that may alter the fitness effects of subsequent mutations--substitute much more often than expected while a population is adapting. We quantify these effects using a simple population-genetic model that elucidates how the substitution rate at conditionally neutral sites depends on the population size, mutation rate, strength of selection, and prevalence of epistasis. We discuss the implications of these results for our understanding of the molecular clock, and for the interpretation of molecular variation in laboratory and natural populations.

Synonymous but not the same: the causes and consequences of codon bias

Despite their name, synonymous mutations have significant consequences for cellular processes in all taxa. As a result, an understanding of codon bias is central to fields as diverse as molecular evolution and biotechnology. Although recent advances in sequencing and synthetic biology have helped to resolve longstanding questions about codon bias, they have also uncovered striking patterns that suggest new hypotheses about protein synthesis. Ongoing work to quantify the dynamics of initiation and elongation is as important for understanding natural synonymous variation as it is for designing transgenes in applied contexts.

The origins of cancer robustness and evolvability

Unless diagnosed early, many adult cancers remain incurable diseases. This is despite an intense global research effort to develop effective anticancer therapies, calling into question the use of rational drug design strategies in targeting complex disease states such as cancer. A fundamental challenge facing researchers and clinicians is that cancers are inherently robust biological systems, able to survive, adapt and proliferate despite the perturbations resulting from anticancer drugs. It is essential that the mechanisms underlying tumor robustness be formally studied and characterized, as without a thorough understanding of the principles of tumor robustness, strategies to overcome therapy resistance are unlikely to be found. Degeneracy describes the ability of structurally distinct system components (e.g. proteins, pathways, cells, organisms) to be conditionally interchangeable in their contribution to system traits and it has been broadly implicated in the robustness and evolvability of complex biological systems. Here we focus on one of the most important mechanisms underpinning tumor robustness and degeneracy, the cellular heterogeneity that is the hallmark of most solid tumors. Based on a combination of computational, experimental and clinical studies we argue that stochastic noise is an underlying cause of tumor heterogeneity and particularly degeneracy. Drawing from a number of recent data sets, we propose an integrative model for the evolution of therapy resistance, and discuss recent computational studies that propose new therapeutic strategies aimed at defeating the adaptable cancer phenotype.

Adaptive evolution of the Chlamydia trachomatis dominant antigen reveals distinct evolutionary scenarios for B- and T-cell epitopes: worldwide survey

Background

Chlamydia trachomatis is one of the most disseminated human pathogens, for which no vaccine is available yet. Understanding the impact of the host pressure on pathogen antigens is crucial, but so far it was only assessed for highly-restricted geographic areas. We aimed to evaluate the evolutionary picture of the chlamydial key antigen (MOMP), which is one of the leading multi-subunit vaccine candidates, in a worldwide basis.

Methodology/Principal Findings

Using genetics, molecular evolution methods and mathematical modelling, we analyzed all MOMP sequences reported worldwide, composed by 5026 strains from 33 geographic regions of five continents. Overall, 35.9% of variants were detected. The evolutionary pattern of MOMP amino acid gains/losses was found to differ from the remaining chromosome, reflecting the demanding constraints of this porin, adhesin and dominant antigen. Amino acid changes were 4.3-fold more frequent in host-interacting domains (P<10⁻¹²), specifically within B-cell epitopes (P<10⁻⁵), where 25% of them are at fixation (P<10⁻⁵). According to the typical pathogen-host arms race, this rampant B-cell antigenic variation likely represents neutralization escape mutants, as some mutations were previously shown to abrogate neutralization of chlamydial infectivity in vitro. In contrast, T-cell clusters of diverse HLA specificities are under purifying selection, suggesting a strategy that may lead to immune subversion. Moreover, several silent mutations are at fixation, generating preferential codons that may influence expression, and may also reflect recombination-derived ‘hitchhiking-effect’ from favourable nonsilent changes. Interestingly, the most prevalent C. trachomatis genotypes, E and F, showed a mutation rate 22.3-fold lower than that of the remainder (P<10⁻²⁰), suggesting more fitted antigenic profiles.

Conclusions/Significance

Globally, the adaptive evolution of the C. trachomatis dominant antigen is likely driven by its complex pathogenesis-related function and reflects distinct evolutionary antigenic scenarios that may benefit the pathogen, and thus should be taking into account in the development of a MOMP-based vaccine.

The synthetic integron: an in vivo genetic shuffling device

As the field of synthetic biology expands, strategies and tools for the rapid construction of new biochemical pathways will become increasingly valuable. Purely rational design of complex biological pathways is inherently limited by the current state of our knowledge. Selection of optimal arrangements of genetic elements from randomized libraries may well be a useful approach for successful engineering. Here, we propose the construction and optimization of metabolic pathways using the inherent gene shuffling activity of a natural bacterial site-specific recombination system, the integron. As a proof of principle, we constructed and optimized a functional tryptophan biosynthetic operon in Escherichia coli. The trpA-E genes along with ‘regulatory’ elements were delivered as individual recombination cassettes in a synthetic integron platform. Integrase-mediated recombination generated thousands of genetic combinations overnight. We were able to isolate a large number of arrangements displaying varying fitness and tryptophan production capacities. Several assemblages required as many as six recombination events and produced as much as 11-fold more tryptophan than the natural gene order in the same context.

The IGF1 small dog haplotype is derived from Middle Eastern grey wolves

Background

A selective sweep containing the insulin-like growth factor 1 (IGF1) gene is associated with size variation in domestic dogs. Intron 2 of IGF1 contains a SINE element and single nucleotide polymorphism (SNP) found in all small dog breeds that is almost entirely absent from large breeds. In this study, we surveyed a large sample of grey wolf populations to better understand the ancestral pattern of variation at IGF1 with a particular focus on the distribution of the small dog haplotype and its relationship to the origin of the dog.

Results

We present DNA sequence data that confirms the absence of the derived small SNP allele in the intron 2 region of IGF1 in a large sample of grey wolves and further establishes the absence of a small dog associated SINE element in all wild canids and most large dog breeds. Grey wolf haplotypes from the Middle East have higher nucleotide diversity suggesting an origin there. Additionally, PCA and phylogenetic analyses suggests a closer kinship of the small domestic dog IGF1 haplotype with those from Middle Eastern grey wolves.

Conclusions

The absence of both the SINE element and SNP allele in grey wolves suggests that the mutation for small body size post-dates the domestication of dogs. However, because all small dogs possess these diagnostic mutations, the mutations likely arose early in the history of domestic dogs. Our results show that the small dog haplotype is closely related to those in Middle Eastern wolves and is consistent with an ancient origin of the small dog haplotype there. Thus, in concordance with past archeological studies, our molecular analysis is consistent with the early evolution of small size in dogs from the Middle East.

See associated opinion by Driscoll and Macdonald: http://jbiol.com/content/9/2/10

Mutational robustness can facilitate adaptation

Robustness seems to be the opposite of evolvability. If phenotypes are robust against mutation, we might expect that a population will have difficulty adapting to an environmental change, as several studies have suggested. However, other studies contend that robust organisms are more adaptable. A quantitative understanding of the relationship between robustness and evolvability will help resolve these conflicting reports and will clarify outstanding problems in molecular and experimental evolution, evolutionary developmental biology and protein engineering. Here we demonstrate, using a general population genetics model, that mutational robustness can either impede or facilitate adaptation, depending on the population size, the mutation rate and the structure of the fitness landscape. In particular, neutral diversity in a robust population can accelerate adaptation as long as the number of phenotypes accessible to an individual by mutation is smaller than the total number of phenotypes in the fitness landscape. These results provide a quantitative resolution to a significant ambiguity in evolutionary theory.

Mutational neighbourhood and mutation supply rate constrain adaptation in Pseudomonas aeruginosa

Understanding adaptation by natural selection requires understanding the genetic factors that determine which beneficial mutations are available for selection. Here, using experimental evolution of rifampicin-resistant Pseudomonas aeruginosa, we show that different genotypes vary in their capacity for adaptation to the cost of antibiotic resistance. We then use sequence data to show that the beneficial mutations associated with fitness recovery were specific to particular genetic backgrounds, suggesting that genotypes had access to different sets of beneficial mutations. When we manipulated the supply rate of beneficial mutations, by altering effective population size during evolution, we found that it constrained adaptation in some selection lines by restricting access to rare beneficial mutations, but that the effect varied among the genotypes in our experiment. These results suggest that mutational neighbourhood varies even among genotypes that differ by a single amino acid change, and this determines their capacity for adaptation as well as the influence of population biology processes that alter mutation supply rate.

Evolutionary dynamics of ompA, the gene encoding the Chlamydia trachomatis key antigen

Chlamydia trachomatis is the trachoma agent and causes most bacterial sexually transmitted infections worldwide. Its major outer membrane protein (MOMP) is a well-known porin and adhesin and is the dominant antigen. So far, investigation of MOMP variability has been focused mainly on molecular epidemiological surveys. In contrast, we aimed to evaluate the impact of the host pressure on this key antigen by analyzing its evolutionary dynamics in 795 isolates from urogenital infections, taking into account the MOMP secondary structure and the sizes/positions of antigenic regions. One-third of the specimens showed a mutational drift from the corresponding genotype, where approximately 42% of the mutations had never been described. Amino acid alterations were sixfold more frequent within B-cell epitopes than in the remaining protein (P = 0.027), and some mutations were also found within or close to T-cell antigenic clusters. Interestingly, the two most ecologically successful genotypes, E and F, showed a mutation rate 60.3-fold lower than that of the other genotypes (P < 10(-8)), suggesting that their efficacy may be the result of a better fitness in dealing with the host immune system rather than of specific virulence factors. Furthermore, the variability exhibited by some genetic variants involved residues that are known to play a critical role during the membrane mechanical movements, contributing to a more stable and flexible porin conformation, which suggests some plasticity to deal with environmental pressure. Globally, these MOMP mutational trends yielded no mosaic structures or important phylogenetic changes, but instead yielded point mutations on specific protein domains, which may enhance pathogen's infectivity, persistence, and transmission.

The bifunctional enzymes of antibiotic resistance

The evolutionary union of two genes--each encoding proteins of complementary enzymatic activity--into a single gene so as to allow the coordinated expression of these activities as a fusion polypeptide, is an increasingly recognized biological occurrence. The result of this genetic union is the bifunctional enzyme. This fusion of separate catalytic activities into a single protein, whose gene is regulated by a single promoter, is seen especially where the coordinated expression of the separate activities is highly desirable. Increasingly, a circumstance driving the evolution of the bifunctional enzyme in bacteria is the resistance response of bacteria to antibiotic chemotherapy. We summarize the knowledge on bifunctional antibiotic-resistance enzymes, as possible harbingers of clinically significant resistance mechanisms of the future.