Rates of molecular evolution in bacteria are relatively constant despite spore dormancy

Rates of molecular evolution are known to vary considerably among lineages, partially due to differences in life-history traits such as generation time. The generation-time effect has been well documented in some eukaryotes, but its prevalence in prokaryotes is unknown. "Because many species of Firmicute bacteria spend long periods of time as metabolically dormant spores, which could result in fewer DNA substitutions per unit time, they present an excellent system for testing predictions of the molecular clock hypothesis." To test whether spore-forming bacteria evolve more slowly than their non-spore-forming relatives, I used phylogenetic methods to determine if there were differences in rates of amino acid substitution between spore-forming and non-spore-forming lineages of Firmicute bacteria. Although rates of evolution do vary among lineages, I find no evidence for an effect of spore-formation on evolutionary rate and, furthermore, evolutionary rates are similar to those calculated for enteric bacteria. These results support the notion that variation in generation time does not affect evolutionary rates in bacterial lineages.     

Male-driven biased gene conversion governs the evolution of base composition in human alu repeats

Regional biases in substitution pattern are likely to be responsible for the large-scale variation in base composition observed in vertebrate genomes. However, the evolutionary forces responsible for these biases are still not clearly defined. In order to study the processes of mutation and fixation across the entire human genome, we analyzed patterns of substitution in Alu repeats since their insertion. We also studied patterns of human polymorphism within the repeats. There is a highly significant effect of recombination rate on the pattern of substitution, whereas no such effect is seen on the pattern of polymorphism. These results suggest that regional biases in substitution are caused by biased gene conversion, a process that increases the probability of fixation of mutations that increase GC content. Furthermore, the strongest correlate of substitution patterns is found to be male recombination rates rather than female or sex-averaged recombination rates. This indicates that in addition to sexual dimorphism in recombination rates, the sexes also differ in the relative rates of crossover and gene conversion.     

Gene conversion drives the evolution of HINTW, an ampliconic gene on the female-specific avian W chromosome

The HINTW gene on the female-specific W chromosome of chicken and other birds is amplified and present in numerous copies. Moreover, as HINTW is distinctly different from its homolog on the Z chromosome (HINTZ), is a candidate gene in avian sex determination, and evolves rapidly under positive selection, it shows several common features to ampliconic and testis-specific genes on the mammalian Y chromosome. A phylogenetic analysis within galliform birds (chicken, turkey, quail, and pheasant) shows that individual HINTW copies within each species are more similar to each other than to gene copies of related species. Such convergent evolution is most easily explained by recurrent events of gene conversion, the rate of which we estimated at 10(-6)-10(-5) per site and generation. A significantly higher GC content of HINTW than of other W-linked genes is consistent with biased gene conversion increasing the fixation probability of mutations involving G and C nucleotides. Furthermore, and as a likely consequence, the neutral substitution rate is almost twice as high in HINTW as in other W-linked genes. The region on W encompassing the HINTW gene cluster is not covered in the initial assembly of the chicken genome, but analysis of raw sequence reads indicates that gene copy number is significantly higher than a previous estimate of 40. While sexual selection is one of several factors that potentially affect the evolution of ampliconic, male-specific genes on the mammalian Y chromosome, data from HINTW provide evidence that gene amplification followed by gene conversion can evolve in female-specific chromosomes in the absence of sexual selection. The presence of multiple and highly similar copies of HINTW may be related to protein function, but, more generally, amplification and conversion offers a means to the avoidance of accumulation of deleterious mutations in nonrecombining chromosomes.     

Mammalian gene evolution: Nucleotide sequence divergence between mouse and rat

As a paradigm of mammalian gene evolution, the nature and extent of DNA sequence divergence between homologous protein-coding genes from mouse and rat have been investigated. The data set examined includes 363 genes totalling 411 kilobases, making this by far the largest comparison conducted between a single pair of species. Mouse and rat genes are on average 93.4% identical in nucleotide sequence and 93.9% identical in amino acid sequence. Individual genes vary substantially in the extent of nonsynonymous nucleotide substitution, as expected from protein evolution studies; here the variation is characterized. The extent of synonymous (or silent) substitution also varies considerably among genes, though the coefficient of variation is about four times smaller than for nonsynonymous substitutions. A small number of genes mapped to the X-chromosome have a slower rate of molecular evolution than average, as predicted if molecular evolution is male-driven. Base composition at silent sites varies from 33% to 95% G + C in different genes; mouse and rat homologues differ on average by only 1.7% in silent-site G + C, but it is shown that this is not necessarily due to any selective constraint on their base composition. Synonymous substitution rates and silent site base composition appear to be related (genes at intermediate G + C have on average higher rates), but the relationship is not as strong as in our earlier analyses. Rates of synonymous and nonsynonymous substitution are correlated, apparently because of an excess of substitutions involving adjacent pairs of nucleotides. Several factors suggest that synonymous codon usage in rodent genes is not subject to selection.     

Experimental evolution of RNA versus DNA viruses

Based on their extremely high mutation rates, RNA viruses have been traditionally considered as the fastest evolving entities in nature. However, recent work has revealed that, despite their greater replication fidelity, single-stranded (ss) DNA viruses can evolve fast in a similar way. To further investigate this issue, we have compared the rates of adaptation and molecular evolution of ssRNA and ssDNA viruses under highly controlled laboratory conditions using the bacteriophages ΦX174, G4, f1, Qβ, SP, and MS2 as model systems. Our results indicate that ssRNA phages evolve faster than ssDNA phages under strong selective pressure, and that their extremely high mutation rates appear to be optimal for this kind of scenario. However, their performance becomes similar to that of ssDNA phages over the longer term or when the population is moderately well-adapted. Interestingly, the roughly 100-fold difference between the mutation rates of ssRNA and ssDNA phages yields less than a fivefold difference in adaptation and nucleotide substitution rates. The results are therefore consistent with the observation that, despite their lower mutation rates, ssDNA viruses can sometimes match the evolvability of RNA viruses.     

Time-dependent rates of molecular evolution

For over half a century, it has been known that the rate of morphological evolution appears to vary with the time frame of measurement. Rates of microevolutionary change, measured between successive generations, were found to be far higher than rates of macroevolutionary change inferred from the fossil record. More recently, it has been suggested that rates of molecular evolution are also time dependent, with the estimated rate depending on the timescale of measurement. This followed surprising observations that estimates of mutation rates, obtained in studies of pedigrees and laboratory mutation-accumulation lines, exceeded long-term substitution rates by an order of magnitude or more. Although a range of studies have provided evidence for such a pattern, the hypothesis remains relatively contentious. Furthermore, there is ongoing discussion about the factors that can cause molecular rate estimates to be dependent on time. Here we present an overview of our current understanding of time-dependent rates. We provide a summary of the evidence for time-dependent rates in animals, bacteria and viruses. We review the various biological and methodological factors that can cause rates to be time dependent, including the effects of natural selection, calibration errors, model misspecification and other artefacts. We also describe the challenges in calibrating estimates of molecular rates, particularly on the intermediate timescales that are critical for an accurate characterization of time-dependent rates. This has important consequences for the use of molecular-clock methods to estimate timescales of recent evolutionary events.     

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.     

Temporal patterns of fruit fly (Drosophila) evolution revealed by mutation clocks

Drosophila melanogaster has been a canonical model organism to study genetics, development, behavior, physiology, evolution, and population genetics for nearly a century. Despite this emphasis and the completion of its nuclear genome sequence, the timing of major speciation events leading to the origin of this fruit fly remain elusive because of the paucity of extensive fossil records and biogeographic data. Use of molecular clocks as an alternative has been fraught with non-clock-like accumulation of nucleotide and amino-acid substitutions. Here we present a novel methodology in which genomic mutation distances are used to overcome these limitations and to make use of all available gene sequence data for constructing a fruit fly molecular time scale. Our analysis of 2977 pairwise sequence comparisons from 176 nuclear genes reveals a long-term fruit fly mutation clock ticking at a rate of 11.1 mutations per kilobase pair per Myr. Genomic mutation clock-based timings of the landmark speciation events leading to the evolution of D. melanogaster show that it shared most recent common ancestry 5.4 MYA with D. simulans, 12.6 MYA with D. erecta+D. orena, 12.8 MYA with D. yakuba+D. teisseri, 35.6 MYA with the takahashii subgroup, 41.3 MYA with the montium subgroup, 44.2 MYA with the ananassae subgroup, 54.9 MYA with the obscura group, 62.2 MYA with the willistoni group, and 62.9 MYA with the subgenus Drosophila. These and other estimates are compatible with those known from limited biogeographic and fossil records. The inferred temporal pattern of fruit fly evolution shows correspondence with the cooling patterns of paleoclimate changes and habitat fragmentation in the Cenozoic.     

Effects of spontaneous mutation accumulation on sex ratio traits in a parasitoid wasp

Sex allocation theory has proved extremely successful at predicting when individuals should adjust the sex of their offspring in response to environmental conditions. However, we know rather little about the underlying genetics of sex ratio or how genetic architecture might constrain adaptive sex-ratio behavior. We examined how mutation influenced genetic variation in the sex ratios produced by the parasitoid wasp Nasonia vitripennis. In a mutation accumulation experiment, we determined the mutability of sex ratio, and compared this with the amount of genetic variation observed in natural populations. We found that the mutability (h(2)(m)) ranges from 0.001 to 0.002, similar to estimates for life-history traits in other organisms. These estimates suggest one mutation every 5-60 generations, which shift the sex ratio by approximately 0.01 (proportion males). In this and other studies, the genetic variation in N. vitripennis sex ratio ranged from 0.02 to 0.17 (broad-sense heritability, H(2)). If sex ratio is maintained by mutation-selection balance, a higher genetic variance would be expected given our mutational parameters. Instead, the observed genetic variance perhaps suggests additional selection against sex-ratio mutations with deleterious effects on other fitness traits as well as sex ratio (i.e., pleiotropy), as has been argued to be the case more generally.     

Methylation-dependent transition rates are dependent on local sequence lengths and genomic regions

Recently, Fryxell and Moon (2005) examined methylation-dependent transition rates (5mC deamination rates), which were calculated by the difference between the CpG transition and GpC transition rates, using 4,437 transition mutations in CpG or GpC dinucleotides. They concluded that 5mC deamination rates were highly dependent on local GC content but not on local sequence lengths over which GC content was calculated or the genomic regions where the mutations occurred. Here, we reexamined these statements by using 292,216 CpG-->TpG/CpA and GpC-->GpT/ApC mutations, an increase of 66 times as much data. Contrary to Fryxell and Moon's conclusions, our analysis indicated that 5mC deamination rates in the human genome were dependent on both the local sequence length and the genomic region. Some explanations for their conclusions were provided.     

Variable Spontaneous Mutation and Loss of Heterozygosity among Heterozygous Genomes in Yeast

Mutation and recombination are the primary sources of genetic variation. To better understand the evolution of genetic variation, it is crucial to comprehensively investigate the processes involving mutation accumulation and recombination. In this study, we performed mutation accumulation experiments on four heterozygous diploid yeast species in the Saccharomycodaceae family to determine spontaneous mutation rates, mutation spectra, and losses of heterozygosity (LOH). We observed substantial variation in mutation rates and mutation spectra. We also observed high LOH rates (1.65–11.07×10⁻⁶ events per heterozygous site per cell division). Biases in spontaneous mutation and LOH together with selection ultimately shape the variable genome-wide nucleotide landscape in yeast species.     

The genome as a life-history character: why rate of molecular evolution varies between mammal species

DNA sequences evolve at different rates in different species. This rate variation has been most closely examined in mammals, revealing a large number of characteristics that can shape the rate of molecular evolution. Many of these traits are part of the mammalian life-history continuum: species with small body size, rapid generation turnover, high fecundity and short lifespans tend to have faster rates of molecular evolution. In addition, rate of molecular evolution in mammals might be influenced by behaviour (such as mating system), ecological factors (such as range restriction) and evolutionary history (such as diversification rate). I discuss the evidence for these patterns of rate variation, and the possible explanations of these correlations. I also consider the impact of these systematic patterns of rate variation on the reliability of the molecular date estimates that have been used to suggest a Cretaceous radiation of modern mammals, before the final extinction of the dinosaurs.     

Endosymbiont phylogenesis in the dryophthoridae weevils: evidence for bacterial replacement

Intracellular symbiosis is widespread in the insect world where it plays an important role in evolution and adaptation. The weevil family Dryophthoridae (Curculionoidea) is of particular interest in intracellular symbiosis evolution with regard to the great economical and ecological features of these invasive insects, and the potential for comparative studies across a wide range of host plants and environments. Here, we have analyzed the intracellular symbiotic bacteria of 19 Dryophthoridae species collected worldwide, representing a wide range of plant species and tissues. All except one (Sitophilus linearis) harbor symbiotic bacteria within specialized cells (the bacteriocytes) assembled as an organ, the bacteriome. Phylogenetic analysis of the 16S rDNA gene sequence of the Dryophthoridae endosymbionts revealed three endosymbiotic clades belonging to gamma3-Proteobacteria and characterized by different GC contents and evolutionary rate. The genus name Candidatus Nardonella was proposed for the ancestral clade infesting Dryophthoridae 100 MYA and represented by five of nine bacterial genera studied. For this clade showing low GC content (40.5% GC) and high evolutionary rate (0.128 substitutions/site per 100 Myr), a single infection and subsequent cospeciation of the host and the endosymbionts was observed. In the two other insect lineage endosymbionts, with relatively high GC content (53.4% and 53.8% GC), competition with ancestral pathogenic bacteria might have occurred, leading to endosymbiont replacement in present-day last insects.     

Genic mutation rates in mammals: local similarity,chromosomal heterogeneity,and X-versus-autosome disparity

The reduction of mutation rates on the mammalian X chromosome relative to autosomes is most often explained in the literature as evidence of male-driven evolution. This hypothesis attributes lowered mutation rates on the X chromosome to the fact that this chromosome spends less time in the germline of males than in the germline of females. In contrast to this majority view, two articles argued that the patterns of mutation rates across chromosomes are inconsistent with male-driven evolution. One article reported a 40% reduction in synonymous substitution rates (Ks) for X-linked genes relative to autosomes in the mouse-rat lineage. The authors argued that this reduction is too dramatic to be explained by male-driven evolution and concluded that selection has systematically reduced mutation rate on the X chromosome to a level optimal for this male-hemizygous chromosome. More recently, a second article found that chromosomal mutation rates in both the human-mouse and mouse-rat lineages were so heterogeneous that the X chromosome was not an outlier. Here again, the authors argued that this is at odds with male-driven evolution and suggested that selection has modulated chromosomal mutation rates to locally optimal levels, thus extending the argument of the first mentioned article to include autosomes. Here, we reexamine these conclusions using mouse-rat and human-mouse coding-region data. We find a more modest reduction of Ks on the X chromosome, but our results contradict the finding that the X chromosome is not distinct from autosomes. Multiple statistical tests show that Ks rates on the X chromosome differ systematically from the autosomes in both lineages. We conclude that the moderate reduction of mutation rate on the X chromosome of both lineages is consistent with male-driven evolution; however, the large variance in mutation rates across chromosomes suggests that mutation rates are affected by additional factors besides male-driven evolution. Investigation of mutation rates by synteny reveals that synteny blocks, rather than entire chromosomes, might represent the unit of mutation rate variation.     

Maternal smoking during pregnancy and cord blood DNA methylation: new insight on sex differences and effect modification by maternal folate levels

Maternal smoking during pregnancy may affect newborn DNA methylation (DNAm). However, little is known about how these associations vary by a newborn’s sex and/or maternal nutrition. To fill in this research gap, we investigated epigenome-wide DNAm associations with maternal smoking during pregnancy in African American mother-newborn pairs. DNAm profiling in cord (n = 379) and maternal blood (n = 300) were performed using the Illumina HumanMethylation450 BeadChip array. We identified 12 CpG sites whose DNAm levels in cord blood were associated with maternal smoking, at a false discovery rate <5%. The identified associations in the GFI1 gene were more pronounced in male newborns than in females (P = 0.002 for maternal smoking × sex interaction at cg18146737). We further observed that maternal smoking and folate level may interactively affect cord blood DNAm level at cg05575921 in the AHRR gene (P = 5.0 × 10^?4 for interaction): compared to newborns unexposed to maternal smoking and with a high maternal folate level (>19.2 nmol/L), the DNAm level was about 0.03 lower (P = 3.6 × 10^?4) in exposed newborns with a high maternal folate level, but was 0.08 lower (P = 1.2 × 10^?8) in exposed newborns with a low maternal folate level. Our data suggest that adequate maternal folate levels may partly counteract the impact of maternal smoking on DNAm. These findings may open new avenues of inquiry regarding sex differences in response to environmental insults and novel strategies to mitigate their intergenerational health effects through optimization of maternal nutrition.     

蕨类植物叶绿体rps12基因的分子进化研究

以68种蕨类植物和2种石松类植物的rps12基因为对象,在系统发育背景下,结合最大似然法,使用HyPhy和PAML软件对该基因进行进化速率和适应性进化研究。结果显示:位于IR区的外显子2~3,其替换率明显降低,rps12基因编码序列的替换率也随之降低,且rps12基因密码子第3位的GC含量明显升高;在蕨类植物的进化过程中,3′-rps12更倾向定位于IR区,以保持较低的替换率;rps12基因编码的123个氨基酸位点中,共检测到4个正选择位点和116个负选择位点。研究结果表明基因序列进入到IR区后,显示出降低的替换率;强烈的负选择压力表明RPS12蛋白的高度保守性以及rps12基因的功能和结构已经趋于稳定。     

A model for the correlation of mutation rate with GC content and the origin of GC-rich isochores

Based on the biochemical kinetics of DNA replication and mutagenesis, including misincorporation and correction, a model has been developed for studying the relationships among the mutation rate (u), the G + C content of the sequence (f), and the G + C proportion in the nucleotide precursor pool (N). Also a measure for the next-nucleotide effect, called the maximum capacity of the next-nucleotide effect (MC), has been proposed. Under the normal physiological conditions of mammalian germ cells, our results indicate: (1) the equilibrium G + C content in a sequence is approximately equal to the G + C proportion in the nucleotide precursor pool, i.e., f ≈ N, which is independent of the next-nucleotide effect; (2) an inverted-V-shaped distribution of mutation rates with respect to G + C contents is predicted, when the next-nucleotide effect is week, i.e., MC ≈ 1; (3) the distribution becomes flatter (i.e., inverted-U-shaped) as MC increases, but the peak at 50% GC is still observed when MC < 2; and (4) the peak disappears when MC > 2.8, that is, when the next-nucleotide effect becomes strong. Our results suggest that changes in the relative concentrations of nucleotide precursors can cause variations among genes both in mutation rate and in G + C content and that compositional isochores (DNA segments with a homogeneous G + C content) can arise in a genome due to differences in replication times of DNA segments. Correspondence to: W.-H. Li