Mutators and sex in bacteria: Conflict between adaptive strategies

Bacterial mutation rates can increase and produce genetic novelty, as shown by in vitro and in silico experiments. Despite the cost due to a heavy deleterious mutation load, mutator alleles, which increase the mutation rate, can spread in asexual populations during adaptation because they remain associated with the rare favorable mutations they generate. This indirect selection for a genetic system generating diversity (second-order selection) is expected to be highly sensitive to changes in the dynamics of adaptation. Here we show by a simulation approach that even rare genetic exchanges, such as bacterial conjugation or transformation, can dramatically reduce the selection of mutators. Moreover, drift or competition between the processes of mutation and recombination in the course of adaptation reveal how second-order selection is unable to optimize the rate of generation of novelty.

Mitochondrial mutational spectra in human cells and tissues

We have found that human organs such as colon, lung, and muscle, as well as their derived tumors, share nearly all mitochondrial hotspot point mutations. Seventeen hotspots, primarily G → A and A → G transitions, have been identified in the mitochondrial sequence of base pairs 10,030–10,130. Mutant fractions increase with the number of cell generations in a human B cell line, TK6, indicating that they are heritable changes. The mitochondrial point mutation rate appears to be more than two orders of magnitude higher than the nuclear point mutation rate in TK6 cells and in human tissues. The similarity of the hotspot sets in vivo and in vitro leads us to conclude that human mitochondrial point mutations in the sequence studied are primarily spontaneous in origin and arise either from DNA replication error or reactions of DNA with endogenous metabolites. The predominance of transition mutations and the high number of hotspots in this short sequence resembles spectra produced by DNA polymerases in vitro.

A thermostable endonuclease III homolog from the archaeon Pyrobaculum aerophilum

Pyrimidine adducts in cellular DNA arise from modification of the pyrimidine 5,6-double bond by oxidation, reduction or hydration. The biological outcome includes increased mutation rate and potential lethality. A major DNA N-glycosylase responsible for the excision of modified pyrimidine bases is the base excision repair (BER) glycosylase endonuclease III, for which functional homologs have been identified and characterized in Escherichia coli, yeast and humans. So far, little is known about how hyperthermophilic Archaea cope with such pyrimidine damage. Here we report characterization of an endonuclease III homolog, PaNth, from the hyperthermophilic archaeon Pyrobaculum aerophilum, whose optimal growth temperature is 100°C. The predicted product of 223 amino acids shares significant sequence homology with several [4Fe-4S]-containing DNA N-glycosylases including E.coli endonuclease III (EcNth). The histidine-tagged recombinant protein was expressed in E.coli and purified. Under optimal conditions of 80–160 mM NaCl and 70°C, PaNth displays DNA glycosylase/β-lyase activity with the modified pyrimidine base 5,6-dihydrothymine (DHT). This activity is enhanced when DHT is paired with G. Our data, showing the structural and functional similarity between PaNth and EcNth, suggests that BER of modified pyrimidines may be a conserved repair mechanism in Archaea. Conserved amino acid residues are identified for five subfamilies of endonuclease III/UV endonuclease homologs clustered by phylogenetic analysis.

Genes, peoples, and languages

The genetic history of a group of populations is usually analyzed by reconstructing a tree of their origins. Reliability of the reconstruction depends on the validity of the hypothesis that genetic differentiation of the populations is mostly due to population fissions followed by independent evolution. If necessary, adjustment for major population admixtures can be made. Dating the fissions requires comparisons with paleoanthropological and paleontological dates, which are few and uncertain. A method of absolute genetic dating recently introduced uses mutation rates as molecular clocks; it was applied to human evolution using microsatellites, which have a sufficiently high mutation rate. Results are comparable with those of other methods and agree with a recent expansion of modern humans from Africa. An alternative method of analysis, useful when there is adequate geographic coverage of regions, is the geographic study of frequencies of alleles or haplotypes. As in the case of trees, it is necessary to summarize data from many loci for conclusions to be acceptable. Results must be independent from the loci used. Multivariate analyses like principal components or multidimensional scaling reveal a number of hidden patterns and evaluate their relative importance. Most patterns found in the analysis of human living populations are likely to be consequences of demographic expansions, determined by technological developments affecting food availability, transportation, or military power. During such expansions, both genes and languages are spread to potentially vast areas. In principle, this tends to create a correlation between the respective evolutionary trees. The correlation is usually positive and often remarkably high. It can be decreased or hidden by phenomena of language replacement and also of gene replacement, usually partial, due to gene flow.

In vivo requirement for RecJ, ExoVII, ExoI, and ExoX in methyl-directed mismatch repair

Biochemical studies with model DNA heteroduplexes have implicated RecJ exonuclease, exonuclease VII, exonuclease I, and exonuclease X in Escherichia coli methyl-directed mismatch correction. However, strains deficient in the four exonucleases display only a modest increase in mutation rate, raising questions concerning involvement of these activities in mismatch repair in vivo. The quadruple mutant deficient in the four exonucleases, as well as the triple mutant deficient in RecJ exonuclease, exonuclease VII, and exonuclease I, grow poorly in the presence of the base analogue 2-aminopurine, and exposure to the base analogue results in filament formation, indicative of induction of SOS DNA damage response. The growth defect and filamentation phenotypes associated with 2-aminopurine exposure are effectively suppressed by null mutations in mutH, mutL, mutS, or uvrD/mutU, which encode activities that act upstream of the four exonucleases in the mechanism for the methyl-directed reaction that has been proposed based on in vitro studies. The quadruple exonuclease mutant is also cold-sensitive, having a severe growth defect at 30°C. This phenotype is suppressed by a uvrD/mutU defect, and partially suppressed by mutH, mutL, or mutS mutations. These observations confirm involvement of the four exonucleases in methyl-directed mismatch repair in vivo and suggest that the low mutability of exonuclease-deficient strains is a consequence of under recovery of mutants due to a reduction in viability and/or chromosome loss associated with activation of the mismatch repair system in the absence of RecJ exonuclease, exonuclease VII, exonuclease I, and exonuclease X.

The notion of the Cambrian pananimalia genome.

The toil by photosynthesizing cyanobacteria and blue-green algae of nearly three billion years appeared to have finally resulted in the sufficient accumulation of molecular oxygen. So, the stage was set for the emergence, at the ocean bottom, of diverse animals that were consumers of molecular oxygen. It now appears that this Cambrian explosion, during which nearly all the extant animal phyla have emerged, was of an astonishingly short duration, lasting only 6-10 million years. Inasmuch as only a 1% DNA base sequence change is expected in 10 million years under the standard spontaneous mutation rate, I propose that all those diverse animals of the early Cambrian period, some 550 million years ago, were endowed with nearly identical genomes, with differential usage of the same set of genes accounting for the extreme diversities of body forms. Some of the more pertinent genes that are thought to be included in the Cambrian pananimalia genome are as follows. (i) A gene for lysyloxidase that, in the presence of molecular oxygen, crosslinked collagen triple helices to produce ligaments and tendons, thus contributing to the stout bodies of the Cambrian animals. (ii) Genes for hemoglobin; these internal transporters of molecular oxygen are today seen sporadically in members of diverse animal phyla. (iii) The Pax-6 gene for eye formation; the eyes of a ribbon worm to a human are organized by this gene. In animals without eyes, the same gene organizes other sensory systems and organs. (iv) A series of Hox genes for the anterior-posterior (cranio-caudal) body plans: these genes are also present in all phyla of the kingdom Animalia.

Replication infidelity during a single cycle of Ty1 retrotransposition.

Retroviruses undergo a high frequency of genetic alterations during the process of copying their RNA genomes. However, little is known about the replication fidelity of other elements that transpose via reverse transcription of an RNA intermediate. The complete sequence of 29 independently integrated copies of the yeast retrotransposon Ty1 (173,043 nt) was determined, and the mutation rate during a single cycle of replication was calculated. The observed base substitution rate of 2.5 x 10(-5) bp per replication cycle suggests that this intracellular element can mutate as rapidly as retroviruses. The pattern and distribution of errors in the Ty1 genome is nonrandom and provides clues to potential in vivo molecular mechanisms of reverse transcriptase-mediated error generation, including heterogeneous RNase H cleavage of Ty1 RNA, addition of terminal nontemplated bases, and transient dislocation and realignment of primer-templates. Overall, analysis of errors generated during Ty1 replication underscores the utility of a genetically tractable model system for the study of reverse transcriptase fidelity.

An immunoglobulin mutator that targets G.C base pairs.

Hypermutation can be defined as an enhancement of the spontaneous mutation rate which the organism uses in certain types of differentiated cells where a high mutation rate is advantageous. At the immunoglobulin loci this process increases the mutation rate > 10(5)-fold over the normal, spontaneous rate. Its proximate cause is called the immunoglobulin mutator system. The most important function of this system is to improve antibody affinity in an ongoing response; it is turned on and off during the differentiation of B lymphocytes. We have established an in vitro system to study hypermutation by transfecting a rearranged mu gene into a cell line in which an immunoglobulin mutator has been demonstrated. A construct containing the mu gene and the 3' kappa enhancer has all the cis-acting elements necessary for hypermutation of the endogenous gene segments encoding the variable region. The activity of the mutator does not seem to depend strongly on the position of the transfected gene in the genome. The mutator is not active in transformed cells of a later differentiation stage. It is also not active on a transfected lacZ gene. These results are consistent with the specificity of the mutator system being maintained and make it possible to delineate cis and trans mutator elements in vitro. Surprisingly, the mutator preferentially targets G-C base pairs. Two hypotheses are discussed: (i) the immunoglobulin mutator system in mammals consists of several mutators, of which the mutator described here is only one; or (ii) the primary specificity of the system is biased toward mutation of G-C base pairs, but this specificity is obscured by antigenic selection.

Muller's ratchet decreases fitness of a DNA-based microbe.

Muller proposed that an asexual organism will inevitably accumulate deleterious mutations, resulting in an increase of the mutational load and an inexorable, ratchet-like, loss of the least mutated class [Muller, H.J. (1964) Mutat. Res. 1, 2-9]. The operation of Muller's ratchet on real populations has been experimentally demonstrated only in RNA viruses. However, these cases are exceptional in that the mutation rates of the RNA viruses are extremely high. We have examined whether Muller's ratchet operates in Salmonella typhimurium, a DNA-based organism with a more typical genomic mutation rate. Cells were grown asexually under conditions expected to result in high genetic drift, and the increase in mutational load was determined. S. typhimurium accumulated mutations under these conditions such that after 1700 generations, 1% of the 444 lineages tested had suffered an obvious loss of fitness, as determined by decreased growth rate. These results suggest that in the absence of sex and with high genetic drift, genetic mechanisms, such as back or compensatory mutations, cannot compensate for the accumulation of deleterious mutations. In addition, we measured the appearance of auxotrophs, which allowed us to calculate an average spontaneous mutation rate of approximately 0.3-1.5 x 10(-9) mutations per base pair per generation. This rate is measured for the largest genetic target studied so far, a collection of about 200 genes.

Structural variation of the pseudoautosomal region between and within inbred mouse strains.

The pseudoautosomal region (PAR) is a segment of shared homology between the sex chromosomes. Here we report additional probes for this region of the mouse genome. Genetic and fluorescence in situ hybridization analyses indicate that one probe, PAR-4, hybridizes to the pseudoautosomal telomere and a minor locus at the telomere of chromosome 9 and that a PCR assay based on the PAR-4 sequence amplifies only the pseudoautosomal locus (DXYHgu1). The region detected by PAR-4 is structurally unstable; it shows polymorphism both between mouse strains and between animals of the same inbred strain, which implies an unusually high mutation rate. Variation occurs in the region adjacent to a (TTAGGG)n array. Two pseudoautosomal probes can also hybridize to the distal telomeres of chromosomes 9 and 13, and all three telomeres contain DXYMov15. The similarity between these telomeres may reflect ancestral telomere-telomere exchange.

Genetic heterogeneity in isogenic homozygous clonal zebrafish.

The C32 isogenic homozygous diploid (IHD) strain of the zebrafish (Danio rerio) was found to be polyallelic at a malate dehydrogenase locus (sMdh-A). A variant allele is thought to have arisen via mutation within the past 10 bisexual generations that have maintained the strain since its last gynogenetic cloning event; this unique allele now predominates at the sMdh-A locus. The estimated mutation rate in this species is sufficiently high that long-term genetic homogeneity of its IHD clones cannot be assumed. Researchers using such bisexually maintained clones should be aware that they are not necessarily using genetically uniform subjects. Genetic uniformity of cloned IHD zebrafish will be maximized if experimental subjects are obtained soon after a cloning event.

Microsatellite variability and genetic distances.

We analyze the within- and between-population dynamics of the distribution of the number of repeats at multiple microsatellite DNA loci subject to stepwise mutation. Analytical expressions for moments up to the fourth order within a locus and the variance of between-locus variance at mutation-drift equilibrium have been obtained. These statistics may be used to test the appropriateness of the one-step mutation model and to detect between-locus variation in the mutation rate. Published data are compatible with the one-step mutation model, although they do not reject the two-step model. Using both multinomial sampling and diffusion approximations for the analysis of the genetic distance introduced by Goldstein et al. [Goldstein, D. B., Linares, A. R., Cavalli-Sforza, L. L. & Feldman, M. W. (1995) Proc. Natl. Acad. Sci. USA 92, 6723-6727], we show that this distance follows a chi 2 distribution with degrees of freedom equal to the number of loci when there is no variation in mutation rates among the loci. In the presence of such variation, the variance of the distance is obtained. We conclude that the number of microsatellite loci required for the construction of phylogenetic trees with reliable branch lengths may be several hundred. Also, mutations that change repeat scores by several units, even though extremely rare, may dramatically influence estimates of population parameters.

Genetic absolute dating based on microsatellites and the origin of modern humans.

We introduce a new genetic distance for microsatellite loci, incorporating features of the stepwise mutation model, and test its performance on microsatellite polymorphisms in humans, chimpanzees, and gorillas. We find that it performs well in determining the relations among the primates, but less well than other distance measures (not based on the stepwise mutation model) in determining the relations among closely related human populations. However, the deepest split in the human phylogeny seems to be accurately reconstructed by the new distance and separates African and non-African populations. The new distance is independent of population size and therefore allows direct estimation of divergence times if the mutation rate is known. Based on 30 microsatellite polymorphisms and a recently reported average mutation rate of 5.6 x 10(-4) at 15 dinucleotide microsatellites, we estimate that the deepest split in the human phylogeny occurred about 156,000 years ago. Unlike most previous estimates, ours requires no external calibration of the rate of molecular evolution. We can use such calibrations, however, to test our estimate.

Adaptive mutations in Escherichia coli as a model for the multiple mutational origins of tumors.

The cells in most tumors are found to carry multiple mutations; however, based upon mutation rates determined by fluctuation tests, the frequency of such multiple mutations should be so low that tumors are never detected within human populations. Fluctuation tests, which determine the cell-division-dependent mutation rate per cell generation in growing cells, may not be appropriate for estimating mutation rates in nondividing or very slowly dividing cells. Recent studies of time-dependent, "adaptive" mutations in nondividing populations of microorganisms suggest that similar measurements may be more appropriate to understanding the mutation origins of tumors. Here I use the ebgR and ebgA genes of Escherichia coli to measure adaptive mutation rates where multiple mutations are required for rapid growth. Mutations in either ebgA or ebgR allow very slow growth on lactulose (4-O-beta-D-galactosyl-D-fructose), with doubling times of 3.2 and 17.3 days, respectively. However, when both mutations are present, cells can grow rapidly with doubling times of 2.7 hr. I show that during prolonged (28-day) selection for growth on lactulose, the number of lactulose-utilizing mutants that accumulate is 40,000 times greater than can be accounted for on the basis of mutation rates measured by fluctuation tests, but is entirely consistent with the time-dependent adaptive mutation rates measured under the same conditions of prolonged selection.

Mutation in bone morphogenetic protein receptor-IB is associated with increased ovulation rate in Booroola Mérino ewes

Ewes from the Booroola strain of Australian Mérino sheep are characterized by high ovulation rate and litter size. This phenotype is due to the action of the FecBB allele of a major gene named FecB, as determined by statistical analysis of phenotypic data. By genetic analysis of 31 informative half-sib families from heterozygous sires, we showed that the FecB locus is situated in the region of ovine chromosome 6 corresponding to the human chromosome 4q22–23 that contains the bone morphogenetic protein receptor IB (BMPR-IB) gene encoding a member of the transforming growth factor-β (TGF-β) receptor family. A nonconservative substitution (Q249R) in the BMPR-IB coding sequence was found to be associated fully with the hyperprolificacy phenotype of Booroola ewes. In vitro, ovarian granulosa cells from FecBB/FecBB ewes were less responsive than granulosa cells from FecB+/FecB+ ewes to the inhibitory effect on steroidogenesis of GDF-5 and BMP-4, natural ligands of BMPR-IB. It is suggested that in FecBB/FecBB ewes, BMPR-IB would be inactivated partially, leading to an advanced differentiation of granulosa cells and an advanced maturation of ovulatory follicles.