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.     

Patterns of transitional mutation biases within and among mammalian genomes

Significant transition/transversion mutation bias is a well-appreciated aspect of mammalian nuclear genomes; however, patterns of bias among genes within a genome and among species remain largely uncharacterized. Understanding these patterns is important for understanding similarities and differences in mutational patterns among genomes and genomic regions. Therefore, we have conducted an analysis of 7,587 pairs of sequences of 4,347 mammalian protein-coding genes from seven species (human, mouse, rat, cow, sheep, pig, and macaque) and from the introns of 51 gene pairs and multiple intergenic regions (37 kbp, 52 kbp and 65 kbp) from the human, chimpanzee, and baboon genomes. Our analyses show that genes and regions with widely varying base composition exhibit uniformity of transition mutation rate both within and among mammalian lineages, as long as the transitional mutations caused by CpG hypermutability are excluded. The estimates show no relationship to potential intrachromosomal or interchromosomal effects. This uniformity points to similarity in point mutation processes in genomic regions with substantially different GC-content biases.     

Noncoding DNA, isochores and gene expression: nucleosome formation potential

The nucleosome formation potential of introns, intergenic spacers and exons of human genes is shown here to negatively correlate with among-tissues breadth of gene expression. The nucleosome formation potential is also found to negatively correlate with the GC content of genomic sequences; the slope of regression line is steeper in exons compared with noncoding DNA (introns and intergenic spacers). The correlation with GC content is independent of sequence length; in turn, the nucleosome formation potential of introns and intergenic spacers positively (albeit weakly) correlates with sequence length independently of GC content. These findings help explain the functional significance of the isochores (regions differing in GC content) in the human genome as a result of optimization of genomic structure for epigenetic complexity and support the notion that noncoding DNA is important for orderly chromatin condensation and chromatin-mediated suppression of tissue-specific genes.     

Modeling DNA Base Substitution in Large Genomic Regions from Two Organisms

We studied the substitution patterns in 7661 well-conserved human–mouse alignments corresponding to the intergenic regions of human chromosome 22. Alignments with a high average GC content tend to have a higher human GC content than mouse GC content, indicating a lack of stationarity. Segmenting the alignments into four groups of GC content and fitting the general reversible substitution model (REV) separately gave significantly better fits than the overall fit and the levels of fit are close to that expected under an REV model. In addition, most of the fitted rate matrices are not of the HKY type but are remarkably strand-symmetric, and we constructed a number of substitution matrices that should be useful for genomic DNA sequence alignment. We did not find obvious signs of temporal inhomogeneity in the substitution rates and concluded that the conserved intergenic regions in human chromosome 22 and mouse appear to have evolved from their common ancestors via a process that is approximately reversible and strand-symmetric, assuming site homogeneity and independence.     

CpG mutation rates in the human genome are highly dependent on local GC content

CpG dinucleotides mutate at a high rate because cytosine is vulnerable to deamination, cytosines in CpG dinucleotides are often methylated, and deamination of 5-methylcytosine (5mC) produces thymidine. Previous experiments have shown that DNA melting is the rate-limiting step in cytosine deamination. Here we show, through the analysis of human single-nucleotide polymorphisms (SNPs), that the mutation rate produced by 5mC deamination is highly dependent on local GC content. In fact, linear regression analysis showed that the log(10) of the 5mC mutation rates (inferred from SNP frequencies) had slopes of -3 when graphed with respect to the GC content of neighboring sequences. This is the ideal slope that would be expected if the correlation between CpG underrepresentation and GC content had been solely caused by DNA melting. Moreover, this same result was obtained regardless of the SNP locations (all SNPs versus only SNPs in noncoding intergenic regions, excluding CpG islands) and regardless of the lengths over which GC content was calculated (SNP sequences with a modal length of 564 bp versus genomic contigs with a modal length of 163 kb). Several alternative interpretations are discussed.     

CpG islands, genes and isochores in the genomes of vertebrates

We have shown that human genes associated with CpG islands increase in number as they increase in % of guanine + cytosine (GC) levels, and that most genes associated with CpG islands are located in the GC-richest compartment of the human genome. This is an independent confirmation of the concentration gradient of CpG islands (detected as HpaII tiny fragments, or HTF) which was demonstrated in the genome of warm-blooded vertebrates [A?ssani and Bernardi, Gene 106 (1991) 173-183]. We then reassessed the location of CpG islands using the data currently available and confirmed that CpG islands are most frequently located in the 5'-flanking sequences of genes and that they overlap genes to variable extents. We have shown that such extents increase with the increasing GC levels of genes, the GC-richest genes being completely included in CpG islands. Under such circumstances, we have investigated the properties of the 'extragenic' CpG islands located in the 5'-flanking segments of homologous genes from both warm- and cold-blooded vertebrates. We have confirmed that, in cold-blooded vertebrates, CpG islands are often absent; when present, they have lower GC and CpG levels; the latter attain, however, statistically expected values. Finally, we have shown that CpG doublets increase with the increasing GC of exons, introns and intergenic sequences (including 'extragenic' CpG islands) in the genomes from both warm- and cold-blooded vertebrates. The correlations found are the same for both classes of vertebrates, and are similar for exons, introns and intergenic sequences (including 'extragenic' CpG islands). The findings just outlined indicate that the origin and evolution of CpG islands in the vertebrate genome are associated with compositional transitions (GC increases) in genes and isochores.     

Variation in mutation dynamics across the maize genome as a function of regional and flanking base composition

We examine variation in mutation dynamics across a single genome (Zea mays ssp. mays) in relation to regional and flanking base composition using a data set of 10,472 SNPs generated by resequencing 1776 transcribed regions. We report several relationships between flanking base composition and mutation pattern. The A + T content of the two sites immediately flanking the mutation site is correlated with rate, transition bias, and GC --> AT pressure. We also observe a significant CpG effect, or increase in transition rate at CpG sites. At the regional level we find that the strength of the CpG effect is correlated with regional A + T content, ranging from a 1.7-fold increase in transition rate in relatively G + C-rich regions to a 2.6-fold increase in A + T-rich regions. We also observe a relationship between locus A + T content and GC --> AT pressure. This regional effect is in opposition to the influence of the two immediate neighbors in that GC --> AT pressure increases with increasing locus A + T content but decreases with increasing flanking base A + T content and may represent a relationship between genome location and mutation bias. The data indicate multiple context effects on mutations, resulting in significant variation in mutation dynamics across the genome.     

Distinct changes of genomic biases in nucleotide substitution at the time of Mammalian radiation

Differences in the regional substitution patterns in the human genome created patterns of large-scale variation of base composition known as genomic isochores. To gain insight into the origin of the genomic isochores, we develop a maximum-likelihood approach to determine the history of substitution patterns in the human genome. This approach utilizes the vast amount of repetitive sequence deposited in the human genome over the past approximately 250 Myr. Using this approach, we estimate the frequencies of seven types of substitutions: the four transversions, two transitions, and the methyl-assisted transition of cytosine in CpG. Comparing substitutional patterns in repetitive elements of various ages, we reconstruct the history of the base-substitutional process in the different isochores for the past 250 Myr. At around 90 MYA (around the time of the mammalian radiation), we find an abrupt fourfold to eightfold increase of the cytosine transition rate in CpG pairs compared with that of the reptilian ancestor. Further analysis of nucleotide substitutions in regions with different GC content reveals concurrent changes in the substitutional patterns. Although the substitutional pattern was dependent on the regional GC content in such ways that it preserved the regional GC content before the mammalian radiation, it lost this dependence afterward. The substitutional pattern changed from an isochore-preserving to an isochore-degrading one. We conclude that isochores have been established before the radiation of the eutherian mammals and have been subject to the process of homogenization since then.     

CpG Distribution and Methylation Pattern in Porcine Parvovirus

Based on GC content and the observed/expected CpG ratio (oCpGr), we found three major groups among the members of subfamily Parvovirinae: Group I parvoviruses with low GC content and low oCpGr values, Group II with low GC content and high oCpGr values and Group III with high GC content and high oCpGr values. Porcine parvovirus belongs to Group I and it features an ascendant CpG distribution by position in its coding regions similarly to the majority of the parvoviruses. The entire PPV genome remains hypomethylated during the viral lifecycle independently from the tissue of origin. In vitro CpG methylation of the genome has a modest inhibitory effect on PPV replication. The in vitro hypermethylation disappears from the replicating PPV genome suggesting that beside the maintenance DNMT1 the de novo DNMT3a and DNMT3b DNA methyltransferases can’t methylate replicating PPV DNA effectively either, despite that the PPV infection does not seem to influence the expression, translation or localization of the DNA methylases. SNP analysis revealed high mutability of the CpG sites in the PPV genome, while introduction of 29 extra CpG sites into the genome has no significant biological effects on PPV replication in vitro. These experiments raise the possibility that beyond natural selection mutational pressure may also significantly contribute to the low level of the CpG sites in the PPV genome.     

Evolutionary Basis of Codon Usage and Nucleotide Composition Bias in Vertebrate DNA Viruses

Understanding the extent and causes of biases in codon usage and nucleotide composition is essential to the study of viral evolution, particularly the interplay between viruses and host cells or immune responses. To understand the common features and differences among viruses we analyzed the genomic characteristics of a representative collection of all sequenced vertebrate-infecting DNA viruses. This revealed that patterns of codon usage bias are strongly correlated with overall genomic GC content, suggesting that genome-wide mutational pressure, rather than natural selection for specific coding triplets, is the main determinant of codon usage. Further, we observed a striking difference in CpG content between DNA viruses with large and small genomes. While the majority of large genome viruses show the expected frequency of CpG, most small genome viruses had CpG contents far below expected values. The exceptions to this generalization, the large gammaherpesviruses and iridoviruses and the small dependoviruses, have sufficiently different life-cycle characteristics that they may help reveal some of the factors shaping the evolution of CpG usage in viruses. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Nicolas Galtier]     

Mutation spectrum of 4-nitroquinoline N-oxide in the lacI transgenic Big Blue Rat2 cell line.

This paper describes the spectrum of mutations induced by 4-nitroquinoline N-oxide (4-NQO) in the lacI target gene of the transgenic Big Blue Rat2 cell line. There are only a few report for the mutational spectrum of 4-NQO in a mammalian system although its biological and genetic effects have been well studied. Big Blue Rat2 cells were treated with 0.03125, 0.0625 or 0.125 microg/ml of 4-NQO, the highest concentration giving 85% survival. Our results indicated that the mutant frequency (MF) induced by 4-NQO was dose-dependent with increases from three- to seven-fold. The DNA sequence analysis of lacI mutants from the control and 4-NQO treatment groups revealed an obvious difference in the spectra of mutations. In spontaneous mutants, transition (60%) mutations, especially G:C-->A:T transition (45%), were most frequent. However, the major type of base substitution after treatment of 4-NQO was transversions (68.8%), especially G:C-->T:A (43.8%), while only 25% of mutants were transitions. These results are consistent with those produced by 4-NQO in other systems and the transgenic assay system will be a powerful tool to postulate more accurately the mechanism of chemical carcinogenesis involved.     

Evolution of noncoding and silent coding sites in the Plasmodium falciparum and Plasmodium reichenowi genomes

We compared levels of sequence divergence between fourfold synonymous coding sites and noncoding sites from the intergenic and intronic regions of the Plasmodium falciparum and Plasmodium reichenowi genomes. We observed significant differences in the level of divergence between these classes of silent sites. Fourfold synonymous coding sites exhibited the highest level of sequence divergence, followed by introns, and then intergenic sequences. This pattern of relative divergence rates has been observed in primate genomes but was unexpected in Plasmodium due to a paucity of variation at silent sites in P. falciparum and the corollary hypothesis that silent sites in this genome may be subject to atypical selective constraints. Exclusion of hypermutable CpG dinucleotides reduces the divergence level of synonymous coding sites to that of intergenic sites but does not diminish the significantly higher divergence level of introns relative to intergenic sites. A greater than expected incidence of CpG dinucleotides in intergenic regions less than 500 bp from genes may indicate selective maintenance of regulatory motifs containing CpGs. Divergence rates of different classes of silent sites in these Plasmodium genomes are determined by a combination of mutational and selective pressures.     

The Influence of Local DNA Sequence and DNA Repair Background on the Mutational Specificity of 1-Nitroso-8-Nitropyrene in Escherichia Coli: Inferences for Mutagenic Mechanisms

We have examined the mutational specificity of 1-nitroso-8-nitropyrene (1,8-NONP), an activated metabolite of the carcinogen 1,8-dinitropyrene, in the lacI gene of Escherichia coli strains which differ with respect to nucleotide excision repair (+/- delta uvrB) and MucA/B-mediated error-prone translesion synthesis (+/- pKM101). Several different classes of mutation were recovered, of which frameshifts, base substitutions, and deletions were clearly induced by 1,8-NONP treatment. The high proportion of point mutations (> 92%) which occurred at G.C sites correlates with the percentage of 1,8-NONP-DNA adducts which occur at the C(8) position of guanine. The most prominent frameshift mutations were -(G.C) events, which were induced by 1,8-NONP treatment in all strains, occurred preferentially in runs of guanine residues, and whose frequency increased markedly with the length of the reiterated sequence. Of the base substitution mutations G.C-->T.A transversions were induced to the greatest extent by 1,8-NONP. The distribution of the G.C-->T.A transversions was not influenced by the nature of flanking bases, nor was there a strand preference for these events. The presence of plasmid pKM101 specifically increased the frequency of G.C-->T.A transversions by a factor of 30-60. In contrast, the -(G.C) frameshift mutation frequency was increased only 2-4-fold in strains harboring pKM101 as compared to strains lacking this plasmid. There was, however, a marked influence of pKM101 on the strand specificity of frameshift mutation; a preference was observed for -G events on the transcribed strand. The ability of the bacteria to carry out nucleotide excision repair had a strong effect on the frequency of all classes of mutation but did not significantly influence either the overall distribution of mutational classes or the strand specificity of G.C-->T.A transversions and -(G.C) frameshifts. Deletion mutations were induced in the delta uvr, pKM101 strain. The endpoints of the majority of the deletion mutations were G.C rich and contained regions of considerable homology. The specificity of 1,8-NONP-induced mutation suggests that DNA containing 1,8-NONP adducts can be processed through different mutational pathways depending on the DNA sequence context of the adduct and the DNA repair background of the cell.     

Regional base composition variation along yeast chromosome III: evolution of chromosome primary structure.

The recent determination of the complete sequence of chromosome III from the yeast Saccharomyces cerevisiae allows, for the first time, the investigation of the long range primary structure of a eukaryotic chromosome. We have found that, against a background G+C level of about 35%, there are two regions (one in each chromosome arm) in which G+C values rise to over 50%. This effect is seen in silent sites within genes, but not in noncoding intergenic sequences. The variation in G+C content is not related to differential selection of synonymous codons, and probably reflects mutational biases. That the intergenic regions do not exhibit the same phenomenon is particularly interesting, and suggests that they are under substantial constraint. The yeast chromosome may be a model of the structure of the human genome, since there is evidence that it is also a mosaic of long regions of different base compositions, reflected in wide variation of G+C content at silent sites among genes. Two possible causes of this regional effect, replication timing, and recombination frequency, are discussed.     

Tri-nucleotide repeats and their association with genes in rice genome

Tri-nucleotide repeats (TNRs) are extremely abundant in rice genome, of which CCG/CGG repeats have an advantage over other repeats, with approximate half of all the TNRs in the genome. Our results show that rice genome has relatively abundant TNRs with high GC content, and containing only purines or pyrimidines under the same GC content. The AAT/ATT repeats that occur predominantly in intergenic and intronic regions have a considerably higher average length than that of other repeats. The highest frequency of TNRs occurs in 5'-UTR regions, followed by in coding and 5'-flanking regions. Purines-rich TNRs prefer to the coding regions, but pyrimidines-rich TNRs exhibit a stronger bias to upstream regions, suggesting that they might be considered as the regulatory elements in gene expression. As if TNRs located predominantly near the start of coding regions do not significantly influence on the protein function.     

Phylogenomic analysis of the uracil-DNA glycosylase superfamily

The spontaneous deamination of cytosine produces uracil mispaired with guanine in DNA, which will produce a mutation, unless repaired. In all domains of life, uracil-DNA glycosylases (UDGs) are responsible for the elimination of uracil from DNA. Thus, UDGs contribute to the integrity of the genetic information and their loss results in mutator phenotypes. We are interested in understanding the role of UDG genes in the evolutionary variation of the rate and the spectrum of spontaneous mutations. To this end, we determined the presence or absence of the five main UDG families in more than 1,000 completely sequenced genomes and analyzed their patterns of gene loss and gain in eubacterial lineages. We observe nonindependent patterns of gene loss and gain between UDG families in Eubacteria, suggesting extensive functional overlap in an evolutionary timescale. Given that UDGs prevent transitions at G:C sites, we expected the loss of UDG genes to bias the mutational spectrum toward a lower equilibrium G + C content. To test this hypothesis, we used phylogenetically independent contrasts to compare the G + C content at intergenic and 4-fold redundant sites between lineages where UDG genes have been lost and their sister clades. None of the main UDG families present in Eubacteria was associated with a higher G + C content at intergenic or 4-fold redundant sites. We discuss the reasons of this negative result and report several features of the evolution of the UDG superfamily with implications for their functional study. uracil-DNA glycosylase, mutation rate evolution, mutational bias, GC content, DNA repair, mutator gene.     

Purifying selection maintains highly conserved noncoding sequences in Drosophila

The majority of metazoan genomes consist of nonprotein-coding regions, although the functional significance of most noncoding DNA sequences remains unknown. Highly conserved noncoding sequences (CNSs) have proven to be reliable indicators of functionally constrained sequences such as cis-regulatory elements and noncoding RNA genes. However, CNSs may arise from nonselective evolutionary processes such as genomic regions with extremely low mutation rates known as mutation "cold spots." Here we combine comparative genomic data from recently completed insect genome projects with population genetic data in Drosophila melanogaster to test predictions of the mutational cold spot model of CNS evolution in the genus Drosophila. We find that point mutations in intronic and intergenic CNSs exhibit a significant reduction in levels of divergence relative to levels of polymorphism, as well as a significant excess of rare derived alleles, compared with either the nonconserved spacer regions between CNSs or with 4-fold silent sites in coding regions. Controlling for the effects of purifying selection, we find no evidence of positive selection acting on Drosophila CNSs, although we do find evidence for the action of recurrent positive selection in the spacer regions between CNSs. We estimate that approximately 85% of sites in Drosophila CNSs are under constraint with selection coefficients (N(e)s) on the order of 10-100, and thus, the estimated strength and number of sites under purifying selection is greater for Drosophila CNSs relative to those in the human genome. These patterns of nonneutral molecular evolution are incompatible with the mutational cold spot hypothesis to explain the existence of CNSs in Drosophila and, coupled with similar findings in mammals, argue against the general likelihood that CNSs are generated by mutational cold spots in any metazoan genome.