Gene conversion as a focusing mechanism for correlated mutations: a hypothesis

Variation in heterozygosity predicts variation in human substitution rates between populations, individuals and genomic regions

The "heterozygote instability" (HI) hypothesis suggests that gene conversion events focused on heterozygous sites during meiosis locally increase the mutation rate, but this hypothesis remains largely untested. As humans left Africa they lost variability, which, if HI operates, should have reduced the mutation rate in non-Africans. Relative substitution rates were quantified in diverse humans using aligned whole genome sequences from the 1,000 genomes project. Substitution rate is consistently greater in Africans than in non-Africans, but only in diploid regions of the genome, consistent with a role for heterozygosity. Analysing the same data partitioned into a series of non-overlapping 2 Mb windows reveals a strong, non-linear correlation between the amount of heterozygosity lost "out of Africa" and the difference in substitution rate between Africans and non-Africans. Putative recent mutations, derived variants that occur only once among the 80 human chromosomes sampled, occur preferentially at the centre of 2 Kb windows that have elevated heterozygosity compared both with the same region in a closely related population and with an immediately adjacent region in the same population. More than half of all substitutions appear attributable to variation in heterozygosity. This observation provides strong support for HI with implications for many branches of evolutionary biology.

Two mechanisms produce mutation hotspots at DNA breaks in Escherichia coli

Mutation hotspots and showers occur across phylogeny and profoundly influence genome evolution, yet the mechanisms that produce hotspots remain obscure. We report that DNA double-strand breaks (DSBs) provoke mutation hotspots via stress-induced mutation in Escherichia coli. With tet reporters placed 2 kb to 2 Mb (half the genome) away from an I-SceI site, RpoS/DinB-dependent mutations occur maximally within the first 2 kb and decrease logarithmically to ∼60 kb. A weak mutation tail extends to 1 Mb. Hotspotting occurs independently of I-site/tet-reporter-pair position in the genome, upstream and downstream in the replication path. RecD, which allows RecBCD DSB-exonuclease activity, is required for strong local but not long-distance hotspotting, indicating that double-strand resection and gap-filling synthesis underlie local hotspotting, and newly illuminating DSB resection in vivo. Hotspotting near DSBs opens the possibility that specific genomic regions could be targeted for mutagenesis, and could also promote concerted evolution (coincident mutations) within genes/gene clusters, an important issue in the evolution of protein functions.

Population-specific links between heterozygosity and the rate human microsatellite evolution

Microsatellites form an abundant class of DNA sequences used widely as genetic markers. Surprisingly, the length of human microsatellites varies highly predictably with distance from Africa, apparently following the linear decline in variability that arose as we colonised the world. Such patterns have been used to argue that heterozygosity modulates the rate of microsatellite evolution. Here I test the ensuing prediction that variation in demographic history will cause individual populations predictably either to lead or to lag any given trend in length. I find that they do: larger populations with locally higher heterozygosity have microsatellites that are longer when a locus is expanding and shorter when a locus is contracting. These patterns remain even after controlling for the stepwise way in which heterozygosity and allele lengths decline across the world. This analysis provides support for a strongly discontinuous model for how human genetic variability is distributed and shows how individual populations differ in the average rate their microsatellites are evolving. Such patterns have the potential to provide a new window onto historical demography.

Using human demographic history to infer natural selection reveals contrasting patterns on different families of immune genes

Detecting regions of the human genome that are, or have been, influenced by natural selection remains an important goal for geneticists. Many methods are used to infer selection, but there is a general reliance on an accurate understanding of how mutation and recombination events are distributed, and the well-known link between these processes and their evolutionary transience introduces uncertainty into inferences. Here, we present and apply two new, independent approaches; one based on single nucleotide polymorphisms (SNPs) that exploits geographical patterns in how humans lost variability as we colonized the world, the other based on the relationship between microsatellite repeat number and heterozygosity. We show that the two methods give concordant results. Of these, the SNP-based method is both widely applicable and detects selection over a well-defined time interval, the last 50 000 years. Analysis of all human genes by their Gene Ontology codes reveals how accelerated and decelerated loss of variability are both preferentially associated with immune genes. Applied to 168 immune genes used as the focus of a previous study, we show that members of the same gene family tend to yield similar indices of selection, even when located on different chromosomes. We hope our approach will provide a useful tool with which to infer where selection has acted to shape the human genome.

Mutation biases and mutation rate variation around very short human microsatellites revealed by human-chimpanzee-orangutan genomic sequence alignments

I have studied mutation patterns around very short microsatellites, focusing mainly on sequences carrying only two repeat units. By using human-chimpanzee-orangutan alignments, inferences can be made about both the relative rates of mutations and which bases have mutated. I find remarkable non-randomness, with mutation rate depending on a base's position relative to the microsatellite, the identity of the base itself and the motif in the microsatellite. Comparing the patterns around AC2 with those around other four-base combinations reveals that AC2 does not stand out as being special in the sense that non-repetitive tetramers also generate strong mutation biases. However, comparing AC2 and AC3 with AC4 reveals a step change in both the rate and nature of mutations occurring, suggesting a transition state, AC4 exhibiting an alternating high-low mutation rate pattern consistent with the sequence patterning seen around longer microsatellites. Surprisingly, most changes in repeat number occur through base substitutions rather than slippage, and the relative probability of gaining versus losing a repeat in this way varies greatly with repeat number. Slippage mutations reveal rather similar patterns of mutability compared with point mutations, being rare at two repeats where most cause the loss of a repeat, with both mutation rate and the proportion of expansion mutations increasing up to 6-8 repeats. Inferences about longer repeat tracts are hampered by uncertainties about the proportion of multi-species alignments that fail due to multi-repeat mutations and other rearrangements.

Mutability and importance of a hypermutable cell subpopulation that produces stress-induced mutants in Escherichia coli

In bacterial, yeast, and human cells, stress-induced mutation mechanisms are induced in growth-limiting environments and produce non-adaptive and adaptive mutations. These mechanisms may accelerate evolution specifically when cells are maladapted to their environments, i.e., when they are are stressed. One mechanism of stress-induced mutagenesis in Escherichia coli occurs by error-prone DNA double-strand break (DSB) repair. This mechanism was linked previously to a differentiated subpopulation of cells with a transiently elevated mutation rate, a hypermutable cell subpopulation (HMS). The HMS could be important, producing essentially all stress-induced mutants. Alternatively, the HMS was proposed to produce only a minority of stress-induced mutants, i.e., it was proposed to be peripheral. We characterize three aspects of the HMS. First, using improved mutation-detection methods, we estimate the number of mutations per genome of HMS-derived cells and find that it is compatible with fitness after the HMS state. This implies that these mutants are not necessarily an evolutionary dead end, and could contribute to adaptive evolution. Second, we show that stress-induced Lac(+) mutants, with and without evidence of descent from the HMS, have similar Lac(+) mutation sequences. This provides evidence that HMS-descended and most stress-induced mutants form via a common mechanism. Third, mutation-stimulating DSBs introduced via I-SceI endonuclease in vivo do not promote Lac(+) mutation independently of the HMS. This and the previous finding support the hypothesis that the HMS underlies most stress-induced mutants, not just a minority of them, i.e., it is important. We consider a model in which HMS differentiation is controlled by stress responses. Differentiation of an HMS potentially limits the risks of mutagenesis in cell clones.

Too many mutants with multiple mutations

It has recently become clear that the classical notion of the random nature of mutation does not hold for the distribution of mutations among genes: most collections of mutants contain more isolates with two or more mutations than predicted by the mutant frequency on the assumption of a random distribution of mutations. Excesses of multiples are seen in a wide range of organisms, including riboviruses, DNA viruses, prokaryotes, yeasts, and higher eukaryotic cell lines and tissues. In addition, such excesses are produced by DNA polymerases in vitro. These "multiples" appear to be generated by transient, localized hypermutation rather than by heritable mutator mutations. The components of multiples are sometimes scattered at random and sometimes display an excess of smaller distances between mutations. As yet, almost nothing is known about the mechanisms that generate multiples, but such mutations have the capacity to accelerate those evolutionary pathways that require multiple mutations where the individual mutations are neutral or deleterious. Examples that impinge on human health may include carcinogenesis and the adaptation of microbial pathogens as they move between individual hosts.

Mutation as a stress response and the regulation of evolvability

Our concept of a stable genome is evolving to one in which genomes are plastic and responsive to environmental changes. Growing evidence shows that a variety of environmental stresses induce genomic instability in bacteria, yeast, and human cancer cells, generating occasional fitter mutants and potentially accelerating adaptive evolution. The emerging molecular mechanisms of stress-induced mutagenesis vary but share telling common components that underscore two common themes. The first is the regulation of mutagenesis in time by cellular stress responses, which promote random mutations specifically when cells are poorly adapted to their environments, i.e., when they are stressed. A second theme is the possible restriction of random mutagenesis in genomic space, achieved via coupling of mutation-generating machinery to local events such as DNA-break repair or transcription. Such localization may minimize accumulation of deleterious mutations in the genomes of rare fitter mutants, and promote local concerted evolution. Although mutagenesis induced by stresses other than direct damage to DNA was previously controversial, evidence for the existence of various stress-induced mutagenesis programs is now overwhelming and widespread. Such mechanisms probably fuel evolution of microbial pathogenesis and antibiotic-resistance, and tumor progression and chemotherapy resistance, all of which occur under stress, driven by mutations. The emerging commonalities in stress-induced-mutation mechanisms provide hope for new therapeutic interventions for all of these processes.

Correlated mutations: advances and limitations. A study on fusion proteins and on the Cohesin-Dockerin families

Correlated mutations have been repeatedly exploited for intramolecular contact map prediction. Over the last decade these efforts yielded several methods for measuring correlated mutations. Nevertheless, the application of correlated mutations for the prediction of intermolecular interactions has not yet been explored. This gap is due to several obstacles, such as 3D complexes availability, paralog discrimination, and the availability of sequence pairs that are required for inter- but not intramolecular analyses. Here we selected for analysis fusion protein families that bypass some of these obstacles. We find that several correlated mutation measurements yield reasonable accuracy for intramolecular contact map prediction on the fusion dataset. However, the accuracy level drops sharply in intermolecular contacts prediction. This drop in accuracy does not occur always. In the Cohesin-Dockerin family, reasonable accuracy is achieved in the prediction of both intra- and intermolecular contacts. The Cohesin-Dockerin family is well suited for correlated mutation analysis. Because, however, this family constitutes a special case (it has radical mutations, has domain repeats, within each species each Dockerin domain interacts with each Cohesin domain, see below), the successful prediction in this family does not point to a general potential in using correlated mutations for predicting intermolecular contacts. Overall, the results of our study indicate that current methodologies of correlated mutations analysis are not suitable for large-scale intermolecular contact prediction, and thus cannot assist in docking. With current measurements, sequence availability, sequence annotations, and underdeveloped sequence pairing methods, correlated mutations can yield reasonable accuracy only for a handful of families.

Mechanism of action of the oxazolidinone antibacterial agents

Oxazolidinones represent a new class of synthetic antibacterial agents active against multiply-resistant Gram-positive pathogens, including methicillin-resistant Staphylococcus aureus (MRSA), penicillin-resistant streptococci and vancomycin-resistant enterococci. Eperezolid and linezolid are two novel analogues, which have respectively completed Phase I and Phase II clinical testing. The lack of cross-resistance between oxazolidinones and other antibiotics supports a novel mechanism of action. Oxazolidinones are protein synthesis inhibitors which target an early step involving the binding of N-formylmethionyl-tRNA to the ribosome. Binding studies demonstrate that these agents interact with the 50S subunit, but not the 30S subunit of the ribosome. Crosslinking experiments provide evidence for an interaction with both the 16S rRNA of the small subunit and the 23S rRNA of the large subunit. Development of resistance in the laboratory is slow, resulting in two independently isolated point mutations at G2447U and G2576U of the 23S rRNA. This review discusses the results of published studies involving oxazolidinone mechanism of action.

Dobzhansky-Muller incompatibilities in protein evolution

We study fitness landscape in the space of protein sequences by relating sets of human pathogenic missense mutations in 32 proteins to amino acid substitutions that occurred in the course of evolution of these proteins. On average, approximately 10% of deviations of a nonhuman protein from its human ortholog are compensated pathogenic deviations (CPDs), i.e., are caused by an amino acid substitution that, at this site, would be pathogenic to humans. Normal functioning of a CPD-containing protein must be caused by other, compensatory deviations of the nonhuman species from humans. Together, a CPD and the corresponding compensatory deviation form a Dobzhansky-Muller incompatibility that can be visualized as the corner on a fitness ridge. Thus, proteins evolve along fitness ridges which contain only approximately 10 steps between successive corners. The fraction of CPDs among all deviations of a protein from its human ortholog does not increase with the evolutionary distance between the proteins, indicating that substitutions that carry evolving proteins around these corners occur in rapid succession, driven by positive selection. Data on fitness of interspecies hybrids suggest that the compensatory change that makes a CPD fit usually occurs within the same protein. Data on protein structures and on cooccurrence of amino acids at different sites of multiple orthologous proteins often make it possible to provisionally identify the substitution that compensates a particular CPD.

Illusory defects and mismatches: why must DNA repair always be (slightly) error prone?

There is growing evidence that recombination is mu;tagenic and that some forms of DNA repair synthesis are error prone. DNA synthesis in mismatch repair might also be error prone. DNA-repair systems detect structural defects in DNA with high efficiency but they occasionally also strike at normal sections of DNA. Considering the diversity of local DNA structure, some DNA sections with complementary sequences are bound to act as preferential false targets for a repair system (i.e. as "illusory defects"). However, if the repair system never changes the sequence upon repair, it will be solicited again and again by the illusory defect, a potentially harmful situation. It is therefore advantageous for a repair system to be, to some extent, error prone. Strong illusory defects may arise at the decanucleotide level and could be the cause of local increases in mutation levels. They might be used to initiate somatic hypermutation pathways.

Experimental evolution of a dense cluster of residues in tyrosyl-tRNA synthetase: quantitative effects on activity, stability and dimerization

A dense cluster of eight residues was identified at the crossing of two alpha-helices in tyrosyl-tRNA synthetase (TyrRS) from the thermophile Bacillus stearothermophilus. Its mechanism of evolution was characterized. Four residues of this cluster are not conserved in TyrRS from the mesophile Escherichia coli. The corresponding mutations were constructed in TyrRS(Delta1), a derivative of TyrRS from B. stearothermophilus in which the anticodon binding domain is deleted. Mutations I52L (i.e. Ile52 into Leu), M55L and L105V did not affect the activity of TyrRS(Delta1) in the pyrophosphate exchange reaction whereas T51P increased it. The kinetic stabilities of TyrRS(Delta1) and its mutant derivatives at 68.5 degreesC were determined from experiments of irreversible thermal precipitation. They were in the order L105V<I52L<T51P<Wild Type</=M55L; mutation I52L partially compensated L105V in these experiments whereas M55L was coupled neither to I52L nor to L105V. Mutations I52L and L105V affected the stability of the dimeric TyrRS(Delta1) at different steps of its unfolding by urea, monitored under equilibrium conditions by spectrofluorometry or size exclusion chromatography. I52L destabilized the association between the subunits even though residue Ile52 is more than 20 A away from the subunit interface. L105V destabilized the monomeric intermediate of unfolding. The two mutational pathways, going from the wild-type TyrRS(Delta1) to the I52L-L105V double mutant through each of the single mutants were not equivalent for the stability of the monomeric intermediate and for the total stability of the dimer. One pathway contained two neutral steps whereas the other pathway contained a destabilizing step followed by a stabilizing step. Mutation I52L allowed L105V along the first pathway and compensated it along the second pathway. Thus, the effects of I52L and L105V on stability depended on the structural context. The gain in activity due to T51P was at the expense of a slight destabilization.

The evolutionary design of error-rates, and the fast fixation enigma

Genetic and non-genetic error-rates are analyzed in parallel for a lower and a higher organism (E. coli and man, respectively). From the comparison of mutation with fixation rates, contrasting proposals are made, concerning the arrangement of error-rates in the two organisms. In E. coli, reproduction is very conservative, but genetic variability is high within populations. Most mutations are discarded by selection, yet single mutational variants of a gene have, on average, little impact on fitness. In man, the mutation rate per generation is high, the variability generated in the population is comparatively low, and most mutations are fixed by drift rather than selection. The variants of a gene are in general more deleterious than in E. coli. There is a discrepancy in the published mutation rates: the rate of mutation fixations in human populations is twice or four times higher than the individual rate of mutation production, a feature which is not consistent with current population genetics models. Two, not mutually exclusive, hypotheses may explain this 'fast fixation enigma': (i) Mutation rates have substantially decreased in recent human evolution and (ii) A substantial fraction of the fixed mutations were generated in a process-such as gene conversion-that violates the principle of independence of mutation events.

Clinical resistance to erythromycin and clindamycin in cutaneous propionibacteria isolated from acne patients is associated with mutations in 23S rRNA

The genetic basis of erythromycin resistance in cutaneous propionibacteria was determined by comparing the nucleotide sequences of the peptidyl transferase region in the 23S rRNAs from 9 susceptible and 26 resistant clinical isolates as well as 4 laboratory-selected erythromycin-resistant mutants of a susceptible strain. In 13 isolates and the 4 laboratory mutants, cross-resistance to macrolides, lincosamides, and B-type streptogramins was associated with an A-->G transition at a position cognate with Escherichia coli 23S rRNA base 2058. These strains were resistant to > or = 512 microg of erythromycin per ml. Two other mutations were identified, an A-->G transition at base 2059 in seven strains, associated with high-level resistance to all macrolides, and a G-->A transition at base 2057 in six strains, associated with low-level resistance to erythromycin. These mutations correspond to three of four phenotypic classes previously identified by using MIC determinations.