Heterozygosity and mutation rate: evidence for an interaction and its implications

If natural selection chose where new mutations occur it might well favour placing them near existing polymorphisms, thereby avoiding disruption of areas that work while adding novelty to regions where variation is tolerated or even beneficial. Such a system could operate if heterozygous sites are recognised and ‘repaired’ during the initial stages of crossing over. Such repairs involve an extra round of DNA replication, providing an opportunity for further mutations, thereby raising the local mutation rate. If so, the changes in heterozygosity that occur when populations grow or shrink could feed back to modulate both the rate and the distribution of mutations. Here, I review evidence from isozymes, microsatellites and single nucleotide polymorphisms that this potential is realised in real populations. I then consider the likely implications, focusing particularly on how these processes might affect microsatellites, concluding that heterozygosity does impact on the rate and distribution of mutations.     

The structure of interrupted human AC microsatellites

Microsatellite lengths change over evolutionary time through a process of replication slippage. A recently proposed model of this process holds that the expansionary tendencies of slippage mutation are balanced by point mutations breaking longer microsatellites into smaller units and that this process gives rise to the observed frequency distributions of uninterrupted microsatellite lengths. We refer to this as the slippage/point-mutation theory. Here we derive the theory's predictions for interrupted microsatellites comprising regions of perfect repeats, labeled segments, separated by dinucleotide interruptions containing point mutations. These predictions are tested by reference to the frequency distributions of segments of AC microsatellite in the human genome, and several predictions are shown not to be supported by the data, as follows. The estimated slippage rates are relatively low for the first four repeats, and then rise initially linearly with length, in accordance with previous work. However, contrary to expectation and the experimental evidence, the inferred slippage rates decline in segments above 10 repeats. Point mutation rates are also found to be higher within microsatellites than elsewhere. The theory provides an excellent fit to the frequency distribution of peripheral segment lengths but fails to explain why internal segments are shorter. Furthermore, there are fewer microsatellites with many segments than predicted. The frequencies of interrupted microsatellites decline geometrically with microsatellite size measured in number of segments, so that for each additional segment, the number of microsatellites is 33.6% less. Overall we conclude that the detailed structure of interrupted microsatellites cannot be reconciled with the existing slippage/point-mutation theory of microsatellite evolution, and we suggest that microsatellites are stabilized by processes acting on interior rather than on peripheral segments.     

The aggregation properties of Escherichia coli proteins associated with their cellular abundance

Proteins are key players in most cellular processes. Therefore, their abundances are thought to be tightly regulated at the gene-expression level. Recent studies indicate, however, that steady-state cellular-protein concentrations correlate better across species than the levels of the corresponding mRNAs; this supports the existence of selective forces to maintain precise cellular-protein concentrations and homeostasis, even if gene-expression levels diverge. One of these forces might be the avoidance of protein aggregation because, in the cell, the folding of proteins into functional conformations might be in competition with anomalous aggregation into non-functional and usually toxic structures in a concentration-dependent manner. The data in the present work provide support for this hypothesis because, in E. coli, the experimental solubility of proteins correlates better with the cellular abundance than with the gene-expression levels. We found that the divergence between protein and mRNAs levels is low for high-abundance proteins. This suggests that because abundant proteins are at higher risk of aggregation, cellular concentrations need to be stringently regulated by gene expression.     

Upper-limit mutation rate estimation for a plant RNA virus

It is generally accepted that mutation rates of RNA viruses are inherently high due to the lack of proofreading mechanisms. However, direct estimates of mutation rate are surprisingly scarce, in particular for plant viruses. Here, based on the analysis of in vivo mutation frequencies in tobacco etch virus, we calculate an upper-bound mutation rate estimation of 3×10⁻⁵ per site and per round of replication; a value which turns out to be undistinguishable from the methodological error. Nonetheless, the value is barely on the lower side of the range accepted for RNA viruses, although in good agreement with the only direct estimate obtained for other plant viruses. These observations suggest that, perhaps, differences in the selective pressures operating during plant virus evolution may have driven their mutation rates towards values lower than those characteristic of other RNA viruses infecting bacteria or animals.     

Heterogeneous evolution of microsatellites revealed by reconstruction of recent mutation history in an invasive apomictic snail, Potamopyrgus antipodarum

Heterogeneous patterns of microsatellite evolution present a major challenge for the development of mutation models, and an improved understanding of the determinants of variation in mutation rates and patterns among loci, alleles and taxa is required. A 19th Century bottleneck associated with the introduction of clones of the snail Potamopyrgus antipodarum to Britain presented an opportunity to reconstruct recent microsatellite evolution within the most common apomictic lineage. There was significant variation in both the number and step size of mutations among the seven loci studied. Patterns of mutability were consistent with higher mutation rates for di- than trinucleotides and for longer alleles at a locus. Mutation size was influenced in a more complex way, decreasing with relative allele length much more strongly for tri-, than dinucleotides. We found support for this latter, highly novel result in the literature via reanalysis of data in a recent genome-scan study of human microsatellites, which showed a similarly disparate pattern of length-dependence between di- and trinucleotides. In spite of the apomictic form of reproduction and an unusually strong excess of microsatellite contractions in P. antipodarum, there were notable similarities with mutation processes of human microsatellites, supporting the wider taxonomic generality of such evolutionary mechanisms.     

Difference between evolutionarily effective and germ line mutation rate due to stochastically varying haplogroup size

Within a Y-chromosome haplogroup defined by unique event mutations, variation in microsatellites can accumulate due to their rapid mutation. Estimates based on pedigrees for the Y-chromosome microsatellite mutation rate are 3 or more times greater than the same estimates from evolutionary considerations. We show by simulation that the haplogroups that survive the stochastic processes of drift and extinction accumulate microsatellite variation at a lower rate than predicted from corresponding pedigree estimates; in particular, under constant total population size, the accumulated variance is on average 3-4 times smaller.     

Characterizing the time dependency of human mitochondrial DNA mutation rate estimates

Previous research has established a discrepancy of nearly anorder of magnitude between pedigree-based and phylogeny-based(human vs. chimpanzee) estimates of the mitochondrial DNA (mtDNA)control region mutation rate. We characterize the time dependencyof the human mitochondrial hypervariable region one mutationrate by generating 14 new phylogeny-based mutation rate estimatesusing within-human comparisons and archaeological dates. Rateestimates based on population events between 15,000 and 50,000years ago are at least 2-fold lower than pedigree-based estimates.These within-human estimates are also higher than estimatesgenerated from phylogeny-based human–chimpanzee comparisons.Our new estimates establish a rapid decay in evolutionary mutationrate between approximately 2,500 and 50,000 years ago and aslow decay from 50,000 to 6 Ma. We then extend this analysisto the mtDNA-coding region. Our within-human coding region mutationrate estimates display a similar, though less rapid, time-dependentdecay. We explore the possibility that multiple hits explainthe discrepancy between pedigree-based and phylogeny-based mutationrates. We conclude that whereas nucleotide substitution modelsincorporating multiple hits do provide a possible explanationfor the discrepancy between pedigree-based and human–chimpanzeemutation rate estimates, they do not explain the rapid declineof within-human rate estimates. We propose that demographicprocesses such as serial bottlenecks prior to the Holocene couldexplain the difference between rates estimated before and after15,000 years ago. Our findings suggest that human mtDNA estimatesof dates of population and phylogenetic events should be adjustedin light of this time dependency of the mutation rate estimates.     

An unusually low microsatellite mutation rate in Dictyostelium discoideum, an organism with unusually abundant microsatellites

The genome of the social amoeba Dictyostelium discoideum is known to have a very high density of microsatellite repeats, including thousands of triplet microsatellite repeats in coding regions that apparently code for long runs of single amino acids. We used a mutation accumulation study to see if unusually high microsatellite mutation rates contribute to this pattern. There was a modest bias toward mutations that increase repeat number, but because upward mutations were smaller than downward ones, this did not lead to a net average increase in size. Longer microsatellites had higher mutation rates than shorter ones, but did not show greater directional bias. The most striking finding is that the overall mutation rate is the lowest reported for microsatellites: approximately 1 x 10(-6) for 10 dinucleotide loci and 6 x 10(-6) for 52 trinucleotide loci (which were longer). High microsatellite mutation rates therefore do not explain the high incidence of microsatellites. The causal relation may in fact be reversed, with low mutation rates evolving to protect against deleterious fitness effects of mutation at the numerous microsatellites.     

Low base‐substitution mutation rate and predominance of insertion‐deletion events in the acidophilic bacterium Acidobacterium capsulatum

Analyses of spontaneous mutation have shown that total genome‐wide mutation rates are quantitatively similar for most prokaryotic organisms. However, this view is mainly based on organisms that grow best around neutral pH values (6.0–8.0). In particular, the whole‐genome mutation rate has not been determined for an acidophilic organism. Here, we have determined the genome‐wide rate of spontaneous mutation in the acidophilic Acidobacterium capsulatum using a direct and unbiased method: a mutation‐accumulation experiment followed by whole‐genome sequencing. Evaluation of 69 mutation accumulation lines of A. capsulatum after an average of ~2900 cell divisions yielded a base‐substitution mutation rate of 1.22 × 10⁻¹⁰ per site per generation or 4 × 10⁻⁴ per genome per generation, which is significantly lower than the consensus value (2.5−4.6 × 10⁻³) of mesothermophilic (~15–40°C) and neutrophilic (pH 6–8) prokaryotic organisms. However, the insertion‐deletion rate (0.43 × 10⁻¹⁰ per site per generation) is high relative to the base‐substitution mutation rate. Organisms with a similar effective population size and a similar expected effect of genetic drift should have similar mutation rates. Because selection operates on the total mutation rate, it is suggested that the relatively high insertion‐deletion rate may be balanced by a low base‐substitution rate in A. capsulatum, with selection operating on the total mutation rate.     

Male-biased mutation, sex linkage, and the rate of adaptive evolution

An interaction between sex-linked inheritance and sex-biased mutation rates may affect the rate of adaptive evolution. Males have much higher mutation rates than females in several vertebrate and plant taxa. When evolutionary rates are limited by the supply of favorable new mutations, then genes will evolve faster when located on sex chromosomes that spend more time in males. For mutations with additive effects, Y-linked genes evolve fastest, followed by Z-linked genes, autosomal genes, X-linked genes, and finally W-linked and cytoplasmic genes. This ordering can change when mutations show dominance. The predicted differences in substitution rates may be detectable at the molecular level. Male-biased mutation could cause adaptive changes to accumulate more readily on certain kinds of chromosomes and favor animals with Z-W sex determination to have rapidly evolving male sexual displays.     

A genomic region evolving toward different GC contents in humans and chimpanzees indicates a recent and regionally limited shift in the mutation pattern

DNA sequences evolving differently in the human and chimpanzee genomes signal recent and regionally limited changes in the process of DNA sequence evolution. Here we present the comparison of 90 kb from the nonrecombining part of the human Y chromosome to the corresponding part of the chimpanzee genome using gorilla as out-group. Our results reveal a significant difference in the region-specific substitution process among the human and chimpanzee lineages. As a consequence, this region experiences a change in its GC content on the human lineage while it resides in compositional equilibrium on the chimpanzee lineage. Based on our analysis, we suggest a recent and species-specific shift in the region's mutation pattern as the cause of its differing evolution in humans and chimpanzees.     

High mutation rate of TPE repeats: a microsatellite in the putative transposase of the hobo element in Drosophila melanogaster

The hobo transposable element contains a polymorphic microsatellite sequence located in its coding region, the TPE repeats. Previous surveys of natural populations of Drosophila melanogaster have detected at least seven different hobo transposons. These natural populations are geographically structured with regard to TPE polymorphism, and a scenario has been proposed for the invasion process. Natural populations have recently been completely invaded by hobo elements with three TPE repeats. New elements then appeared by mutation, triggering a new stage of invasion by other elements. Since TPE polymorphism appeared over a short period of time, we focused on estimating the mutation rate of these TPE repeats. We used transgenic lines harboring three TPE and/or five TPE hobo elements that had been evolving for at least 16 generations to search for a new TPE repeat polymorphism. We detected three mutants, with four, seven, and eight TPE repeats, respectively. The estimated mutation rate of the TPE repeats is therefore higher than that of neutral microsatellites in D. melanogaster (4.2 x 10-4 versus 6.5 x 10-6). The role of the transposition mechanism and the particular structure of the TPE repeats of the hobo element in this increase in the mutation rate are discussed.     

fcsA29 mutation is an allele of polA gene of Escherichia coli

The cold-sensitive fcsA29 mutation of Escherichia coli was found to be a new type of cold-sensitive allele of the polA gene encoding DNA polymerase I, caused by an Asp(116)-->Asn change in the 5'-->3' exonuclease domain. The fcsA29 mutant showed typical polA mutant phenotypes such as UV sensitivity and unacceptability of recA mutation. Cold-sensitive growth of the mutant was suppressed by introduction of a sulA mutation, indicating that cell filamentation was due to the SOS response.     

The Escherichia coli dnaW mutation is an allele of the adk gene

Summary A dnaW mutant, isolated on the basis of inability to effect conjugal DNA transfer at high temperature, has been shown by complementation and enzyme assay to be defective in the adk (adenylate kinase; EC 2.7.4.3) locus. The adk mutant, known to have reduced ATP concentration at the nonpermissive temperature (Cousin and Belaich 1966), was used to demonstrate a donor energy requirement for stable aggregate formation and for chromosome transfer in conjugation.     

Robustness of coalescent estimators to between-lineage mutation rate variation

Data from HIV and from human neoplastic cells can show substantial between-lineage mutation rate variation even within a single population. Such variation may affect estimators of population quantities such as Theta = 4N(e)mu. Using simulated DNA data, I measured the effect of rate variation on recovery of Theta by the summary-statistic estimator of Watterson (Watterson GA. 1975. On the number of segregating sites in genetical systems without recombination. Theor Popul Biol. 7:256-276) and the coalescent maximum likelihood algorithm LAMARC (Kuhner MK. 2006. LAMARC 2.0: maximum likelihood and Bayesian estimation of population parameters. Bioinformatics. Advance Access doi: 10.1093/bioinformatics/btk051). Watterson's estimator showed a downward bias, as expected, with high values of Theta. LAMARC's mean estimate was accurate for all tested values of Theta and rate variation except for a downward bias when rate variation was maximal (i.e., the slow rate was zero). LAMARC had consistently narrower confidence intervals (CIs) than Watterson's estimator. Both methods tended to reject the truth too often when rate variation was 8x or greater and independent among branches, as well as when variation was 4x or greater and correlated among branches. In the case of Watterson's estimate, this excess rejection was fully attributable to variation among genealogies in the amount of total branch length associated with the fast and slow rates. However, in the case of LAMARC, some excess rejection was still observed even when between-genealogy variation was taken into account. Both estimators are robust to modest rate variation; however, their use should be coupled with a statistical test to rule out extreme rate variation as the resulting CIs may not be reliable.     

YqjD is an inner membrane protein associated with stationary-phase ribosomes in Escherichia coli

Here, we provide evidence that YqjD, a hypothetical protein of Escherichia coli, is an inner membrane and ribosome binding protein. This protein is expressed during the stationary growth phase, and expression is regulated by stress response sigma factor RpoS. YqjD possesses a transmembrane motif in the C-terminal region and associates with 70S and 100S ribosomes at the N-terminal region. Interestingly, E. coli possesses two paralogous proteins of YqjD, ElaB and YgaM, which are expressed and bind to ribosomes in a similar manner to YqjD. Overexpression of YqjD leads to inhibition of cell growth. It has been suggested that YqjD loses ribosomal activity and localizes ribosomes to the membrane during the stationary phase.     

含自杀基因回路的大肠杆菌的生长和变异规律

具有鲁棒性的基因回路构建是合成生物学的基础工作.基于群体感应的自杀基因回路转化大肠杆菌后赋予宿主菌在一定菌群密度时启动自杀的特性.为了使基因回路更具鲁棒性,在转化后以不同浓度的IPTG为诱导剂,观察宿主菌的表型特征以及在自杀过程中突变菌产生的规律,并通过基因测序确定变异位点.结果表明:IPTG浓度与细菌自杀率呈正相关,自杀强度愈大,突变菌株出现得越早,蔓延的速度也越快.基因测序结果表明,在基因回路质粒上,luxR基因中间有转座子插入,从而破坏了群体感应系统.实验也表明:仅需该突变,就足以使宿主菌逃避自杀.结果为下一步优化基因回路设计,实现细菌密度的连续振荡提供了思路.     

Directional evolution of size coupled with ascertainment bias for variation in Drosophila microsatellites

Species-specific differences in microsatellite locus length and ascertainment bias have both been proposed to explain differences in microsatellite variability and length usually observed when loci isolated in one species are used to survey variation in a related species. Here we provide a simple algebraic approach to independently estimate the contributions of true species-specific length differences and ascertainment bias. We apply this approach to a reciprocal-isolation microsatellite study and show contributions of both ascertainment bias and a true longer average microsatellite length in Drosophila melanogaster compared with D. simulans.