Comparative Genomic Study Reveals a Transition from TA Richness in Invertebrates to GC Richness in Vertebrates at CpG Flanking Sites： An Indication for Context-Dependent Mutagenicity of Methylated CpG Sites

Vertebrate genomes are characterized with CpG deficiency, particularly for GCpoor regions. The GC content-related CpG deficiency is probably caused by context-dependent deamination of methylated CpG sites. This hypothesis was examined in this study by comparing nucleotide frequencies at CpG flanking positions among invertebrate and vertebrate genomes. The finding is a transition of nucleotide preference of 5＇ T to 5＇ A at the invertebrate-vertebrate boundary, indicating that a large number of CpG sites with 5＇ Ts were depleted because of global DNA methylation developed in vertebrates. At genome level, we investigated CpG observed/expected （obs/exp） values in 500 bp fragments, and found that higher CpG obs/exp value is shown in GC-poor regions of invertebrate genomes （except sea urchin） but in GC-rich sequences of vertebrate genomes. We next compared GC content at CpG flanking positions with genomic average, showing that the GC content is lower than the average in invertebrate genomes, but higher than that in vertebrate genomes. These results indicate that although 5＇ T and 5＇ A are different in inducing deamination of methylated CpG sites, GC content is even more important in affecting the deamination rate. In all the tests, the results of sea urchin are similar to vertebrates perhaps due to its fractional DNA methylation. CpG deficiency is therefore suggested to be mainly a result of high mutation rates of methylated CpG sites in GC-poor regions.     

The distribution of the dinucleotide CpG and cytosine methylation in the vitellogenin gene family

Summary Sequence data from regions of five vertebrate vitellogenin genes were used to examine the frequency, distribution, and mutability of the dinucleotide CpG, the preferred modification site for eukaryotic DNA methyltransferases. The observed level of the CpG dinucleotide in all five genes was markedly lower than that expected from the known mononucleotide frequencies. CpG suppression was greater in introns than in exons. CpG-containing codons were found to be avoided in the vitellogenin genes, but not completely despite the redundancy of the genetic code. Frequency and distribution patterns of this dinucleotide varied dramatically among these otherwise closely related genes. Dense clusters of CpG dinucleotides tended to appear in regions of either functional or structural interest (e.g., in the transposon-like Vi-element ofXenopus) and these clusters contained 5-methylcytosine (5 mC). 5 mC is known to undergo deamination to form thymidine, but the extent to which this transition occurs in the heavily methylated genomes of vertebrates and its contribution to CpG suppression are still unclear. Sequence comparison of the methylated vitellogenin gene regions identified CT and GA substitutions that were found to occur at relatively high frequencies. The predicted products of CpG deamination, TpG and CpA, were elevated. These findings are consistent with the view that CpG distribution and methylation are interdependent and that deamination of 5 mC plays an important role in promoting evolutionary change at the nucleotide sequence level.     

mCCG methylation in angiosperms

Two different methods were used to investigate the abundance of cytosine methylation at the outer (5′) position in 5′-CCG-3′ trinucleotides in angiosperm genomes. Mspl is unable to cut its target site if the outer cytosine is methylated (5′-^mCCGG-3′). Using Mspl restriction analysis, it was shown that 5′-^mCCG-3′ is present in all angiosperm genomes examined, and that the amount of cytosine methylation at this site varies between species. Subsequently, direct measurements were made of the amount of methylation at both cytosines in a subset of 5′-CCG-3′ trinucleotides in the Arabidopsis thaliana genome. Based upon these analyses, it was estimated that approximately 20–30% of 5′-CCG-3′ trinucleotides in A. thaliana are methylated at the outer cytosine. Approximately 20% of the 5′-CCG-3′ trinucleotides contain 5-methyl-cytosine at the inner cytosine position, which corresponds to a previous determination of 5′-^mCG-3′ methylation in A. thaliana. The implications of 5′-^mCCG-3′ methylation are discussed.     

Asymmetrical distribution of CpG in an 'average' mammalian gene.

The frequency and distribution of the rare dinucleotide CpG was examined in 15 mammalian genes. CpG is highly methylated at cytosine in mammalian DNA (1,2) and 5-methylcytosine (5mC) is thought to undergo a transition mutation via deamination to produce thymine (3). This would result in the accumulation of TpG and CpA and depletion of CpG during evolution (4). Consistent with this hypothesis, the gene sample of 26,541 dinucleotides contained CpG at 40% the frequency expected by base composition and the CpG transition products, TpG+CpA, were significantly elevated at 124% of expected random frequency. However, because CpG occurs at only 25% of expected random frequency in the genome, the sampled genes were considerably enriched in this dinucleotide. CpGs were asymmetrically distributed in sequences flanking the genes. 5'-flanking sequences were enriched in CpG at 135% of the frequency expected assuming a symmetrical distribution of all the CpGs in the sampled genes (p less than 0.01), while 3'-flanking regions were depleted in CpG at 40% of expected values (p less than 0.0001). This asymmetry may reflect the role of 5-methylcytosine in gene expression. In contrast the frequencies of GpC and GpT+ ApC did not differ significantly from that predicted by base composition and these dinucleotides were not asymmetrically distributed.     

Characterization of the level,target sites and inheritance of cytosine methylation in tomato nuclear DNA

The tomato nuclear genome was determined to have a G+C content of 37% which is among the lowest reported for any plant species. Non-coding regions have a G+C content even lower (32% average) whereas coding regions are considerably richer in G+C (46%).5-methyl cytosine was the only modified base detected and on average 23% of the cytosine residues are methylated. Immature tissues and protoplasts have significantly lower levels of cytosine methylation (average 20%) than mature tissues (average 25%). Mature pollen has an intermediate level of methylation (22%). Seeds gave the highest value (27%), suggesting de novo methylation after pollination and during seed development.Based on isoschizomer studies we estimate 55% of the CpG target sites (detected by Msp I/Hpa II) and 85% of the CpNpG target sites (detected by Bst NI/Eco RI)are methylated. Unmethylated target sites (both CpG and CpNpG) are not randomly distributed throughout the genome, but frequently occur in clusters. These clusters resemble CpG islands recently reported in maize and tobacco.The low G+C content and high levels of cytosine methylation in tomato may be due to previous transitions of 5mCT. This is supported by the fact that G+C levels are lowest in non-coding portions of the genome in which selection is relaxed and thus transitions are more likely to be tolerated. This hypothesis is also supported by the general deficiency of methylation target sites in the tomato genome, especially in non-coding regions.Using methylation isoschizomers and RFLP analysis we have also determined that polymorphism between plants, for cytosine methylation at allelic sites, is common in tomato. Comparing DNA from two tomato species, 20% of the polymorphisms detected by Bst NI/Eco RII could be attributed to differential methylation at the CpNpG target sites. With Msp I/Hpa II, 50% of the polymorphisms were attributable to methylation (CpG and CpNpG sites). Moreover, these polymorphisms were demonstrated to be inherited in a mendelian fashion and to co-segregate with the methylation target site and thus do not represent variation for transacting factors that might be involved in methylation of DNA. The potential role of heritable methylation polymorphism in evolution of gene regulation and in RFLP studies is discussed.     

Frequency of abnormal human haemoglobins caused by C----T transitions in CpG dinucleotides

A large part of human genetic disease apparently arises from deamination of cytosine residues in methylated CpG dinucleotides. Their mutation rate is known to be high when C is present as 5-methyl-cytosine, but is believed to be normal when it is unmethylated. The beta-globin gene contains five, the gamma-globin gene two, and each of the alpha-globin genes contains 35 CpG dinucleotides. The CpG dinucleotides in the beta and gamma-globin genes are methylated, while those in the alpha-globin genes are under-methylated. One would therefore have expected the CpG dinucleotides to be a frequent source of mutations in the beta and gamma-globin genes, but not in the alpha-globin genes. In fact, the evidence points to CpG dinucleotides being a frequent source of mutations in both the alpha and beta-globin genes. This suggests either that the mutation rates of both methylated and unmethylated CpG dinucleotides are abnormally high, which conflicts with published evidence, or that there is a finite chance of some of these in the alpha-globin genes of certain individuals being methylated and therefore subject to mutation.     

The loss of CpC dinucleotides from DNA. II. Methylated and non-methylated genes of vertebrates

The frequencies of neighboring b.p. in more than 1100 genes of vertebrates in the EMBL bank (1000 kb) have been analysed. It has been found that the majority of such genes exhibit a lack of CpG duplexes and an excess of TpG+CpA. The loss of CpG may indicate that the major part of these sites in the genome is methylated and has been subjected to the pressure of CpG----TpG+CpA mutations. The methylated genes grouped into compartment M+ are represented by a fraction of repeated sequences and by genes of the most rapidly diverging families of proteins (globins, immunoglobulins, structural proteins, etc.). The genes of this compartment are characterized by a correlation between the G+C content and the value of CpG-suppression. A group of genes has been detected in which the CpG mutation process has gone so far that nearly all of these dinucleotides have disappeared from DNA. Judging by the value of CpG-suppression, these genes, grouped in the Mo+ compartment, used to be strongly methylated before. However, in the now extant vertebrates they have fully depleted their CpG reserve and for this reason lost the methylation capacity. Transitions in methylated CpG may be one of the sources of spontaneous mutagenesis resulting in the enhanced genetic instability of the cell. A gene compartment has been detected with an intermediate level of CpG deficiency; this compartment has been designated as M+. In these genes only a few of the available CpGs have been steadily methylated (and subjected to mutation). It has been found that the genome of vertebrates contains a specific CpG-rich fraction which exhibits no CpG-suppression, irrespective of the overall content of G+C. Probably, CpG sites have persisted unmethylated throughout the existence of these genes. We suggest them to constitute a M- compartment. This compartment comprises the genes of tRNA and rRNA (5S, 5.8S, 18S, 28S) and small nuclear RNAs U2-U6, as well as the genes of core histones, some enzymes, viruses and 5'-flanking sequences of certain protein-coding genes. In the genome of vertebrates, the genes of the evolutionary most conserved proteins and RNAs have not undergone methylation. A list of genes, belonging to different compartments of the vertebrate genome, is given. Compartment Mo+ constitutes 19%, M(+)--35%, M(+/-)--28% and M(-)--8% of all the vertebrate genes studied. Possible mechanisms, protecting the functionally most significant genes of vertebrates from methylation, and discussed.     

Methylation levels at selected CpG sites in the factor VIII and FGFR3 genes, in mature female and male germ cells: implications for male-driven evolution.

Transitional mutations at CpG dinucleotides account for approximately a third of all point mutations. These mutations probably arise through spontaneous deamination of 5-methylcytosine. Studies of CpG mutation rates in disease-linked genes, such as factor VIII and FGFR3, have indicated that they more frequently originate in male than in female germ cells. It has been speculated that these sex-biased mutation rates might be a consequence of sex-specific methylation differences between the female and the male germ lines. Using the bisulfite-based genomic-sequencing method, we investigated the methylation status of the human factor VIII and FGFR3 genes in mature male and female germ cells. With the exception of a single CpG, both genes were found to be equally and highly methylated in oocytes and spermatocytes. Whereas these observations strongly support the notion that DNA methylation is the major determining factor for recurrent CpG germ-line mutations in patients with hemophilia and achondroplasia, the higher mutation rate in the male germ line is apparently not a simple reflection of sex-specific methylation differences.     

Novel methylation at GpC dinucleotide in the fish Sparus aurata genome

To date, vertebrate DNA has been found methylated at the 5 position of cytosine exclusively in dinucleotide CpG or CpNpG stretches. On the the other hand, we determined that cytosine was methylated unusually in dinucleotide GpC at 5-GGCC-3 sequences in the teleost Sparus aurata EcoRI satellite DNA family. This finding is the first example of methylated GpC sequences in the eukaryotic genomes. At this regard, we have examined the relative methylation levels at this site of the highly repetitive EcoRI satellite DNA family from Sparus aurata different tissues. The EcoRI repeat was remarkably more methylated in male germ cells but hypomethylated in female germ cells at the Hae III restriction site ( GpC). The novel modification and the differential methylation pattern suggest that EcoRI satellite could have a structural and/or functional role at the centromeres of Sparus aurata.     

A human CpG island randomly inserted into a plant genome is protected from methylation

In vertebrate genomes the dinucleotide CpG is heavily methylated, except in CpG islands, which are normally unmethylated. It is not clear why the CpG islands are such poor substrates for DNA methyltransferase. Plant genomes display methylation, but otherwise the genomes of plants and animals represent two very divergent evolutionary lines. To gain a further understanding of the resistance of CpG islands to methylation, we introduced a human CpG island from the proteasome-like subunit I gene into the genome of the plant Arabidopsis thaliana. Our results show that prevention of methylation is an intrinsic property of CpG islands, recognized even if a human CpG island is transferred to a plant genome. Two different parts of the human CpG island – the promoter region/ first exon and exon2–4 – both displayed resistance against methylation, but the promoter/ exon1 construct seemed to be most resistant. In contrast, certain sites in a plant CpG-rich region used as a control transgene were always methylated. The frequency of silencing of the adjacent nptII (Km^R) gene in the human CpG constructs was lower than observed for the plant CpG-rich region. These results have implications for understanding DNA methylation, and for construction of vectors that will reduce transgene silencing.     

Methylation of CpG dinucleotides in the lacI gene of the Big Blue® transgenic mouse

Cytosine residues at CpG dinucleotides can be methylated by endogenous methyltransferases in mammalian cells. The resulting 5-methylcytosine base may undergo spontaneous deamination to form thymine causing G/C to A/T transition mutations. Methylated CpGs also can form preferential targets for environmental mutagens and carcinogens. The Big Blue® transgenic mouse has been used to investigate tissue and organ specificity of mutations and to deduce mutational mechanisms in a mammal in vivo. The transgenic mouse contains approximately 40 concatenated lambda-like shuttle vectors, each of which contains one copy of an Escherichia coli lacI gene as a mutational target. lacI mutations in lambda transgenic mice are characterized by a high frequency of spontaneous mutations targeted to CpG dinucleotides suggesting an important contribution from methylation-mediated events. To study the methylation status of CpGs in the lacI gene, we have mapped the distribution of 5-methylcytosines along the DNA-binding domain and flanking sequences of the lacI gene of transgenic mice. We analyzed genomic DNA from various tissues including thymus, liver, testis, and DNA derived from two thymic lymphomas. The mouse genomic DNAs and methylated and unmethylated control DNAs were chemically cleaved, then the positions of 5-methylcytosines were mapped by ligation-mediated PCR which can be used to distinguish methylated from unmethylated cytosines. Our data show that most CpG dinucleotides in the DNA binding domain of the lacI gene are methylated to a high extent (>98%) in all tissues tested; only a few sites are partially (70–90%) methylated. We conclude that tissue-specific methylation is unlikely to contribute significantly to tissue-specific mutational patterns, and that the occurrence of common mutation sites at specific CpGs in the lacI gene is not related to selective methylation of only these sequences. The data confirm previous suggestions that the high frequency of CpG mutations in lacI transgenes is related to the presence of 5-methylcytosine bases.     

Depletion of CpG Dinucleotides in Papillomaviruses and Polyomaviruses: A Role for Divergent Evolutionary Pressures

Background

Papillomaviruses and polyomaviruses are small ds-DNA viruses infecting a wide-range of vertebrate hosts. Evidence supporting co-evolution of the virus with the host does not fully explain the evolutionary path of papillomaviruses and polyomaviruses. Studies analyzing CpG dinucleotide frequencies in virus genomes have provided interesting insights on virus evolution. CpG dinucleotide depletion has not been extensively studied among papillomaviruses and polyomaviruses. We sought to analyze the relative abundance of dinucleotides and the relative roles of evolutionary pressures in papillomaviruses and polyomaviruses.

Methods

We studied 127 full-length sequences from papillomaviruses and 56 full-length sequences from polyomaviruses. We analyzed the relative abundance of dinucleotides, effective codon number (ENC), differences in synonymous codon usage. We examined the association, if any, between the extent of CpG dinucleotide depletion and the evolutionary lineage of the infected host. We also investigated the contribution of mutational pressure and translational selection to the evolution of papillomaviruses and polyomaviruses.

Results

All papillomaviruses and polyomaviruses are CpG depleted. Interestingly, the evolutionary lineage of the infected host determines the extent of CpG depletion among papillomaviruses and polyomaviruses. CpG dinucleotide depletion was more pronounced among papillomaviruses and polyomaviruses infecting human and other mammals as compared to those infecting birds. Our findings demonstrate that CpG depletion among papillomaviruses is linked to mutational pressure; while CpG depletion among polyomaviruses is linked to translational selection. We also present evidence that suggests methylation of CpG dinucleotides may explain, at least in part, the depletion of CpG dinucleotides among papillomaviruses but not polyomaviruses.

Conclusions

The extent of CpG depletion among papillomaviruses and polyomaviruses is linked to the evolutionary lineage of the infected host. Our results highlight the existence of divergent evolutionary pressures leading to CpG dinucleotide depletion among small ds-DNA viruses infecting vertebrate hosts.     

The DeltauvrB mutations in the Ames strains of Salmonella span 15 to 119 genes

The lacI transgene used in the Big Blue (BB) mouse and rat mutation assays typically displays spontaneous mutation frequencies in the 5x10(-5) range. Recently, the bone marrow and bladder of the Big Blue rat were reported to have, by an order of magnitude, the lowest spontaneous mutation frequencies ever observed for lacI in a transgenic animal, approaching the value for endogenous targets such as hprt ( approximately 10(-6)). Since spontaneous mutations in transgenes have been attributed in part to deamination of 5-methylcytosine in CpG sequences, we have investigated the methylation status of the lacI transgene in bone marrow of BB rats and compared it to that present in other tissues including liver, spleen, and breast. The first 400 bases of the lacI gene were investigated using bisulfite genomic sequencing since this region contains the majority of both spontaneous and induced mutations. Surprisingly, all the CpG cytosines in the lacI sequence were fully methylated in all the tissues examined from both 2- and 14-week-old rats. Thus, there is no correlation between 5-methylcytosine content at CpG sites in lacI and the frequency of spontaneous mutation at this marker. We also investigated the methylation status of another widely used transgenic mutation target, the cII gene. The CpG sites in cII in BB rats were fully methylated while those in BB mice were partially methylated (each site approximately 50% methylated). Since spontaneous mutation frequency at cII is comparable in rat and mouse, the methylation status of CpG sequences in this gene also does not correlate with spontaneous frequency. We conclude that other mechanisms besides spontaneous deamination of 5-methylcytosine at CpG sites are driving spontaneous mutation at BB transgenic loci.     

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/CpNpG motifs in the coding region are preferred sites for mutagenesis in the breast cancer susceptibility genes

The range of BRCA1/BRCA2 gene mutations is diverse and the mechanism accounting for this heterogeneity is obscure. To gain insight into the endogenous mutational mechanisms involved, we evaluated the association of specific sequences (i.e. CpG/CpNpG motifs, homonucleotides, short repeats) and mutations within the genes. We classified 1337 published mutations in BRCA1 (1765 BRCA2 mutations) for each specific sequence, and employed computer simulation combined with mathematical calculations to estimate the true underlying tendency of mutation occurrence. Interestingly, we found no mutational bias to homonucleotides and repeats in deletions/insertions and substitutions but striking bias to CpG/CpNpG in substitutions in both genes. This suggests that methylation-dependent DNA alterations would be a major mechanism for mutagenesis.