On the Organization of Higher Chromosomes

OHTA and Kimura¹ have argued that only about 6% of the sequences in mammalian DNA can be under the intense selection that has characterized the evolutionary history of the cytochromes c, the globin chains and the histones. From the calculated mutation rate of fibrinopeptides A and B they show that if all genes are subjected to the same mutation rate 8.3 mutations would accumulate per genome per generation. Because 0,5 deleterious mutations per genome per generation is the maximum allowable in an equilibrium population², they conclude that the amount of DNA that codes for informational sequences such as the cytochromes, globins and histones must be no more than 0.5/8.3, or 6%. We are therefore left with the interesting observation that 94% of mammalian nuclear DNA serves a function not under strong selection. These authors make several assumptions, one of which is that the spontaneous mutation rate characteristic of a species is constant over all nucleotide sequences. I suggest here that this assumption is incorrect, for a variety of reasons and that by assuming that spontaneous mutation rates vary sequence by sequence, one can arrive at a plausible organizing principle for the structure of higher chromosomes.     

Patterns of Somatic Mutations in Immunoglobulin Variable Genes

The mechanism responsible for somatic mutation in the variable genes of antibodies is unknown and may differ from previously described mechanisms that produce mutation in DNA. We have analyzed 421 somatic mutations from the rearranged immunoglobulin variable genes of mice to determine if the nucleotide substitutions differ from those generated during meiosis and if the presence of nearby direct and inverted repeated sequences could template mutations around the variable gene. The results reveal a difference in the pattern of substitutions obtained from somatic mutations vs. meiotic mutations. An increased frequency of T:A to C:G transitions and a decreased frequency of mutations involving a G in the somatic mutants compared to the meiotic mutants is indicated. This suggests that the mutational processes responsible for somatic mutations in antibody genes differs from that responsible for mutation during meiosis. An analysis of the local DNA sequences revealed many direct repeats and palindromic sequences that were capable of templating some of the known mutations. Although additional factors may be involved in targeting mutations to the variable gene, mistemplating by nearby repeats may provide a mechanism for the enhancement of somatic mutation.     

Microsatellite mutations in the germline: implications for evolutionary inference

Microsatellite DNA sequences mutate at rates several orders of magnitude higher than that of the bulk of DNA. Such high rates mean that spontaneous mutations that form new-length variants can realistically be seen in pedigree analysis. Data on observed mutation events from various organisms are now accumulating, allowing inferences on DNA sequence evolution to be made through an unusually direct approach. Here I discuss and integrate microsatellite mutation data in an evolutionary context. A striking feature of the mutation process is that it seems highly heterogeneous, with distinct differences between species, repeat types, loci and alleles. Age and sex also affect the mutation rate. Within genomes at equilibrium, the microsatellite-length distribution is a delicate balance between biased mutation processes and point mutations acting towards the decay of repetitive DNA. Indeed, simple repeats do not evolve simply.     

The impact of population expansion and mutation rate heterogeneity on DNA sequence polymorphism

In order to study the effect of mutation rate heterogeneity on patterns of DNA polymorphism, we simulated samples of DNA sequences with gamma- distributed nucleotide substitution rates in stationary and expanding populations. We find that recent population expansions and mutation rate heterogeneity have similar effects on several polymorphism indicators, like the shape and the mean of the observed pairwise difference distribution, or the number of segregating sites. The inferred size of population expansion thus appears overestimated if nucleotides have dissimilar substitution rates. Interestingly, population expansion and uneven mutation rates have contrasting effects on Tajima's D statistic when acting separately, and the consequence on the associated test of selective neutrality is investigated. The patterns of polymorphism of several human populations analyzed for the mitochondrial control region are examined, mainly showing the difficulty in quantifying the respective contribution of past demographic history and uneven mutation rates from a single sampled evolutionary process. However, substitution rates appear more heterogeneous in the second hypervariable segment of the control region than in the first segment.      

Mechanisms of spontaneous mutagenesis: an analysis of the spectrum of spontaneous mutation in the Escherichia coli lacI gene

We have obtained via DNA sequence analysis a spectrum of 174 spontaneous mutations occurring in the lac I gene of Escherichia coli. The spectrum comprised base substitution, frameshift, deletion, duplication and insertion mutations, of which the relative contributions to spontaneous mutation could be estimated. Two thirds of all lacI mutations occurred in the frameshift hotspot site. An analysis of the local DNA sequence suggested that the intensity of this hotspot may depend on structural features of the DNA that extend beyond those permitted by the repeated tetramer at this site. Deletions comprised the largest non-hotspot class (37%). They could be divided into two subclasses, depending on whether they included the lac operator sequence; the latter was found to be a preferred site for deletion endpoints. Most of the deletions internal to the lacI gene were associated with the presence of directly or invertedly repeated sequences capable of accounting for their endpoints. Base substitutions comprised 34% of the non-hotspot events. Unlike the base substitution spectrum obtained via nonsense mutations, G . C----A . T transitions do not predominate. A new base substitution hotspot was discovered at position +6 in the lac operator; its intensity may reflect specific features of the operator DNA. IS1 insertion mutations contributed 12% of the non-hotspot mutations and occurred dispersed throughout the gene in both orientations. Since the lacI gene is not A + T-rich, the contribution of IS1 insertion to spontaneous mutation in general might be underestimated. Single-base frameshift mutations were found only infrequently. In general, they did not occur in runs of a common base. Instead, their occurrence seemed based on the "perfection" of direct or inverted repeats in the local DNA sequence. Three (tandem) duplication events were recovered. No repeated sequences were found that might have determined their endpoints.     

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.     

Somatic immunoglobulin sequence divergence and its implications for studies of evolutionary divergence

The divergence of immunoglobulin genes due to somatic mutation provides a natural example of DNA sequence divergence. This divergence was examined to gain insight into the processes of evolution and the determinants of the variance-to-mean ratio of sequence divergence. Normally, this ratio is found to be larger than expected (1.0 under Poisson assumptions) for the evolutionary divergence or most genes. Although not significantly less than one, all seven groups of immunoglobulin amino acid sequences have ratios smaller than expected, contrary to the evolutionary pattern generally observed. The substitutions in the immunoglobulin genes appear to be highly nonrandom and an excess of parallel changes (the major nonrandom feature of these mutations) is shown to cause smaller ratios. Because convergent or parallel mutations are often observed in the evolutionary divergence of genes, this suggests that forces causing the large observed ratios may actually have to be more powerful than previously expected. Further, since selection is one of the likely causes of parallel mutations, it should be noted that selection could significantly decrease the variance-to-mean ratio. The high frequency of parallel mutations and their resulting effects, as observed in the immunoglobulin genes, suggest that only poor inferences of sequence divergence can be made without actual knowledge of the ancestral sequence.     

Short-range template switching in great ape genomes explored using pair hidden Markov models

Many complex genomic rearrangements arise through template switch errors, which occur in DNA replication when there is a transient polymerase switch to an alternate template nearby in three-dimensional space. While typically investigated at kilobase-to-megabase scales, the genomic and evolutionary consequences of this mutational process are not well characterised at smaller scales, where they are often interpreted as clusters of independent substitutions, insertions and deletions. Here we present an improved statistical approach using pair hidden Markov models, and use it to detect and describe short-range template switches underlying clusters of mutations in the multi-way alignment of hominid genomes. Using robust statistics derived from evolutionary genomic simulations, we show that template switch events have been widespread in the evolution of the great apes’ genomes and provide a parsimonious explanation for the presence of many complex mutation clusters in their phylogenetic context. Larger-scale mechanisms of genome rearrangement are typically associated with structural features around breakpoints, and accordingly we show that atypical patterns of secondary structure formation and DNA bending are present at the initial template switch loci. Our methods improve on previous non-probabilistic approaches for computational detection of template switch mutations, allowing the statistical significance of events to be assessed. By specifying realistic evolutionary parameters based on the genomes and taxa involved, our methods can be readily adapted to other intra- or inter-species comparisons.     

Clonality and Evolutionary History of Rhabdomyosarcoma

To infer the subclonality of rhabdomyosarcoma (RMS) and predict the temporal order of genetic events for the tumorigenic process, and to identify novel drivers, we applied a systematic method that takes into account germline and somatic alterations in 44 tumor-normal RMS pairs using deep whole-genome sequencing. Intriguingly, we find that loss of heterozygosity of 11p15.5 and mutations in RAS pathway genes occur early in the evolutionary history of the PAX-fusion-negative-RMS (PFN-RMS) subtype. We discover several early mutations in non-RAS mutated samples and predict them to be drivers in PFN-RMS including recurrent mutation of PKN1. In contrast, we find that PAX-fusion-positive (PFP) subtype tumors have undergone whole-genome duplication in the late stage of cancer evolutionary history and have acquired fewer mutations and subclones than PFN-RMS. Moreover we predict that the PAX3-FOXO1 fusion event occurs earlier than the whole genome duplication. Our findings provide information critical to the understanding of tumorigenesis of RMS.     

A Major Role for Bacteriophage T4 DNA Polymerase in Frameshift Mutagenesis

T4 DNA polymerase strongly influences the frequency and specificity of frameshift mutagenesis. Fifteen of 19 temperature-sensitive alleles of the DNA polymerase gene substantially influenced the reversion frequencies of frameshift mutations measured in the T4 rII genes. Most polymerase mutants increased frameshift frequencies, but a few alleles (previously noted as antimutators for base substitution mutations) decreased the frequencies of certain frameshifts while increasing the frequencies of others. The various patterns of enhanced or decreased frameshift mutation frequencies suggest that T4 DNA polymerase is likely to play a variety of roles in the metabolic events leading to frameshift mutation. A detailed genetic study of the specificity of the mutator properties of three DNA polymerase alleles (tsL56, tsL98 and tsL88) demonstrated that each produces a distinctive frameshift spectrum. Differences in frameshift frequencies at similar DNA sequences within the rII genes, the influence of mutant polymerase alleles on these frequencies, and the presence or absence of the dinucleotide sequence associated with initiation of Okazaki pieces at the frameshift site has led us to suggest that the discontinuities associated with discontinuous DNA replication may contribute to spontaneous frameshift mutation frequencies in T4.     

假基因的组成、分布及其分子进化

假基因(pseudogene)是指基因组中与正常基因序列相似, 但是缺乏功能的DNA 序列。通过序列同源性搜索, 可以收集基因组中假基因的群体特性、染色体分布和同源家族等特性。假基因很好地保留了数百万年前基因组中祖先基因的分子记录, 被视为“基因化石”, 因此假基因在进化和比较基因组学中是重要的资源。应用假基因和基因比较体系, 可以探究生物基因的进化史和基因组稳定性。如: 用Ka/Ks比值确定假基因的自然选择压、物种亲缘关系和进化距离, 分析假基因自身的进化趋势, 探讨DNA 突变的成因等。     

Mechanisms of spontaneous mutation in DNA repair-proficient Escherichia coli

This paper describes the DNA sequence analysis of 729 independent spontaneous lacI⁻ mutation This total is comprised of 478 novel mutations and 251 previously described events, and therefore should allow a more comprehensive view of spontaneous mutation in Escherichia coli. The spectrum is dominated by a hotspot (71% of all events). Mutations at this site consist of related addition and deletion events involving a number of repetitive sequences. Here we discuss how the frequency and proportion of these events vary in different DNA repair-deficient genetic backgrounds. The distribution of non-hotspot events includes base substitutions (38%), deletions (35%), frameshifts (14%), duplications (4%) and insertion elements (4%). G:C → A:T events dominate among base substitutions, while G:C → C:G events are the least common; the remaining types of base substitution are equally represented. Among deletions, a significant number do not display repeated sequences at their endpoints (26/72). However, almost all multiply recovered events (15/17) possess repeated sequences capable of accounting for the deletion endpoints. Similarily, over of all duplications recovered (5/7) display repeated endpoints. Single-base frameshifts are equally divided between A:T and G:C sites, in each case (−) 1 events occur 3-fold more frequently that (+)1 events. A comparative analysis of each mutational class recovered to lacI⁻ spectra available in a variety of DNA repair/metabolism-deficient strains is presented here in an attempt to assess possible contributions from chemical, physical and enzymic sources of damage.     

Effect of target gene CpG content on spontaneous mutation in transgenic mice

Transgenic mutation assays utilizing bacterial target genes display a high frequency of spontaneous mutation at CpG sequences. This is believed to result from the fact that: (1) the prokaryotic genes currently being used as transgenic mutation targets have a high CpG content and (2) these sequences are methylated by mammalian cells to produce 5-methylcytosine (5MC), a known promutagenic base. To study the effect of CpG content on the frequency and type of spontaneous mutation, we have synthesized an analogue of the bacterial lacI target gene (mrkII) that contains a reduced number of CpG sequences. This gene was inserted into a lambda vector and used to construct trangenic mice that undergo vector rescue from genomic DNA upon in vitro packaging. Results on spontaneous mutation frequency and spectrum have been collected and compared to those observed at the lacI gene in Big Blue™ transgenic mice. Spontaneous mutations at the mrkII gene occurred at a frequency in the mid-10⁻⁵ range and were predominantly base pair substitutions, similar to results seen in Big Blue™. However, mrkII mutations were distributed toward the carboxyl end of the gene instead of the bias toward the amino terminus seen in lacI. Unexpectedly, 23% of the spontaneous mrkII mutations were GC → AT transitions at CpG sequences (compared to 32% in lacI), despite the reduction in CpG number from 95 in lacI to only 13 in mrkII. Nine of the CpG bases undergoing transition mutations in mrkII have not been recorded previously as spontaneous sites in Big Blue™. Therefore, substantial reduction of the number of CpG sequences in the lacI transgene did not significantly reduce the rate of spontaneous mutation or alter the contribution of CpG-related events. This suggests that other factors are also operating to establish frequency and composition of spontaneous mutations in transgenic targets.     

Estimates of DNA and protein sequence divergence: an examination of some assumptions

Some of the assumptions underlying estimates of DNA and protein sequence divergence are examined. A solution for the variance of these estimates that allows for different mutation rates and different population sizes in each species and for an arbitrary structure in the initial population is obtained. It is shown that these conditions do not strongly affect estimates of divergence. In general, they cause the variance of divergence to be smaller than a binomial variance. Thus, the binomial variance that is usually assumed for these estimates is safely conservative. It is shown that variability in the mutation rate among sites can have an effect as large as or larger than variability in the mutation rate among bases. Variability in the mutation rate among bases and among sites causes the number of substitutions between two sequences to be underestimated. Protein and DNA sequences from several species are collected to estimate the variability in mutation rates among sites. When many homologous sequences are known, standard methods to estimate this variability can be used. The estimates of this variability show that this factor is important when considering the spectrum of spontaneous mutations and is strongly reflected in the divergence of sequences. Smaller variability is found for the third position of codons than for the first and second codon positions. This may be because of less selective constraints on this position or because the third position has been saturated with mutations for the sequences examined.      

Evolution of mammalian chitinase(-like) members of family 18 glycosyl hydrolases

Family 18 of glycosyl hydrolases encompasses chitinases and so-called chi-lectins lacking enzymatic activity due to amino acid substitutions in their active site. Both types of proteins widely occur in mammals although these organisms lack endogenous chitin. Their physiological function(s) as well as evolutionary relationships are still largely enigmatic. An overview of all family members is presented and their relationships are described. Molecular phylogenetic analyses suggest that both active chitinases (chitotriosidase and AMCase) result from an early gene duplication event. Further duplication events, followed by mutations leading to loss of chitinase activity, allowed evolution of the chi-lectins. The homologous genes encoding chitinase(-like) proteins are clustered in two distinct loci that display a high degree of synteny among mammals. Despite the shared chromosomal location and high homology, individual genes have evolved independently. Orthologs are more closely related than paralogues, and calculated substitution rate ratios indicate that protein-coding sequences underwent purifying selection. Substantial gene specialization has occurred in time, allowing for tissue-specific expression of pH optimized chitinases and chi-lectins. Finally, several family 18 chitinase-like proteins are present only in certain lineages of mammals, exemplifying recent evolutionary events in the chitinase protein family.     

Likelihood models of somatic mutation and codon substitution in cancer genes

The role of somatic mutation in cancer is well established and several genes have been identified that are frequent targets. This has enabled large-scale screening studies of the spectrum of somatic mutations in cancers of particular organs. Cancer gene mutation databases compile the results of many studies and can provide insight into the importance of specific amino acid sequences and functional domains in cancer, as well as elucidate aspects of the mutation process. Past studies of the spectrum of cancer mutations (in particular genes) have examined overall frequencies of mutation (at specific nucleotides) and of missense, nonsense, and silent substitution (at specific codons) both in the sequence as a whole and in a specific functional domain. Existing methods ignore features of the genetic code that allow some codons to mutate to missense, or stop, codons more readily than others (i.e., by one nucleotide change, vs. two or three). A new codon-based method to estimate the relative rate of substitution (fixation of a somatic mutation in a cancer cell lineage) of nonsense vs. missense mutations in different functional domains and in different tumor tissues is presented. Models that account for several potential influences on rates of somatic mutation and substitution in cancer progenitor cells and allow biases of mutation rates for particular dinucleotide sequences (CGs and dipyrimidines), transition vs. transversion bias, and variable rates of silent substitution across functional domains (useful in detecting investigator sampling bias) are considered. Likelihood-ratio tests are used to choose among models, using cancer gene mutation data. The method is applied to analyze published data on the spectrum of p53 mutations in cancers. A novel finding is that the ratio of the probability of nonsense to missense substitution is much lower in the DNA-binding and transactivation domains (ratios near 1) than in structural domains such as the linker, tetramerization (oligomerization), and proline-rich domains (ratios exceeding 100 in some tissues), implying that the specific amino acid sequence may be less critical in structural domains (e.g., amino acid changes less often lead to cancer). The transition vs. transversion bias and effect of CpG dinucleotides on mutation rates in p53 varied greatly across cancers of different organs, likely reflecting effects of different endogenous and exogenous factors influencing mutation in specific organs.     

Estimating the frequency of events that cause multiple-nucleotide changes

Existing mathematical models of DNA sequence evolution assume that all substitutions derive from point mutations. There is, however, increasing evidence that larger-scale events, involving two or more consecutive sites, may also be important. We describe a model, denoted SDT, that allows for single-nucleotide, doublet, and triplet mutations. Applied to protein-coding DNA, the SDT model allows doublet and triplet mutations to overlap codon boundaries but still permits data to be analyzed using the simplifying assumption of independence of sites. We have implemented the SDT model for maximum-likelihood phylogenetic inference and have applied it to an alignment of mammalian globin sequences and to 258 other protein-coding sequence alignments from the Pandit database. We find the SDT model's inclusion of doublet and triplet mutations to be overwhelmingly successful in giving statistically significant improvements in fit of model to data, indicating that larger-scale mutation events do occur. Distributions of inferred parameter values over all alignments analyzed suggest that these events are far more prevalent than previously thought. Detailed consideration of our results and the absence of any known mechanism causing three adjacent nucleotides to be substituted simultaneously, however, leads us to suggest that the actual evolutionary events occurring may include still-larger-scale events, such as gene conversion, inversion, or recombination, or a series of rapid compensatory changes.     

DNA polymerase zeta introduces multiple mutations when bypassing spontaneous DNA damage in Saccharomyces cerevisiae

Spontaneous DNA damage can be dealt with by multiple repair/bypass pathways that have overlapping specificities. We have used a frameshift reversion assay to examine spontaneous mutations that accumulate in yeast strains defective for the high-fidelity nucleotide excision repair or recombination pathways. In contrast to the simple frameshift mutations that occur in wild-type strains, the reversion events in mutant strains are often complex in nature, with the selected frameshift mutation being accompanied by one or more base substitutions. Genetic analyses demonstrate that the complex events are dependent on the Pol zeta translesion polymerase, thus implicating the DNA damage bypass activity of low-fidelity translesion polymerases in hypermutation phenomena.     

Direct estimation of the mitochondrial DNA mutation rate in Drosophila melanogaster

Mitochondrial DNA (mtDNA) variants are widely used in evolutionary genetics as markers for population history and to estimate divergence times among taxa. Inferences of species history are generally based on phylogenetic comparisons, which assume that molecular evolution is clock-like. Between-species comparisons have also been used to estimate the mutation rate, using sites that are thought to evolve neutrally. We directly estimated the mtDNA mutation rate by scanning the mitochondrial genome of Drosophila melanogaster lines that had undergone approximately 200 generations of spontaneous mutation accumulation (MA). We detected a total of 28 point mutations and eight insertion-deletion (indel) mutations, yielding an estimate for the single-nucleotide mutation rate of 6.2 × 10⁻⁸ per site per fly generation. Most mutations were heteroplasmic within a line, and their frequency distribution suggests that the effective number of mitochondrial genomes transmitted per female per generation is about 30. We observed repeated occurrences of some indel mutations, suggesting that indel mutational hotspots are common. Among the point mutations, there is a large excess of G→A mutations on the major strand (the sense strand for the majority of mitochondrial genes). These mutations tend to occur at nonsynonymous sites of protein-coding genes, and they are expected to be deleterious, so do not become fixed between species. The overall mtDNA mutation rate per base pair per fly generation in Drosophila is estimated to be about 10× higher than the nuclear mutation rate, but the mitochondrial major strand G→A mutation rate is about 70× higher than the nuclear rate. Silent sites are substantially more strongly biased towards A and T than nonsynonymous sites, consistent with the extreme mutation bias towards A+T. Strand-asymmetric mutation bias, coupled with selection to maintain specific nonsynonymous bases, therefore provides an explanation for the extreme base composition of the mitochondrial genome of Drosophila.     

Evolution of the major histocompatibility complex: independent origin of nonclassical class I genes in different groups of mammals

The class I major histocompatibility complex genes are composed of classical and nonclassical genes, the latter being largely nonfunctional. To understand the evolutionary relationships of the two groups of class I genes, a phylogenetic analysis of DNA sequences was conducted using 45 genes from six mammalian and one avian species. The results indicate that nonclassical genes in one species are more closely related to classical genes from the same species than to nonclassical genes from a species belonging to a different order or family. This indicates that the differentiation of classical and nonclassical genes occurs rather rapidly in the genome. Classical genes are apparently duplicated with a high frequency in the evolutionary process, and many of the duplicated genes seem to degenerate into nonclassical genes as a result of deleterious mutation. The nonclassical Qa genes in the mouse have sequences homologous to regulatory sequences involved in the universal expression of classical class I genes, but they have accumulated numerous nucleotide substitutions in these sequences. The pattern of nucleotide substitution in nonclassical genes is different from that in classical genes. In nonclassical genes, the rate of nonsynonymous substitution is higher in the antigen recognition site than in other gene regions, as is true of classical genes. However, unlike the case of classical genes, the nonsynonymous rate does not always exceed the synonymous rate in the antigen recognition site. Nonclassical proteins further differ from classical proteins in having amino acid replacements in conserved antigen recognition site positions. These observations are consistent with the hypothesis that nonclassical genes have originated from classical genes but have lost classical class I function because of deleterious mutation.