Codon usage in Homo sapiens: evidence for a coding pattern on the non-coding strand and evolutionary implications of dinucleotide discrimination

This study reports the analysis of codon usage in 35 complete Homo sapiens genes. Both codon frequency and inter-codon interference exhibit patterns of evolutionary interest. There is a significant positive correlation between the frequency with which a given codon is used and the frequency with which its complement is used. Since the frequency of appearance of the complementary codon on the coding strand is equal to the frequency of appearance of the original codon on the non-coding strand, in the same phase, the non-coding strand is found to resemble the coding strand in triplet composition. The same effect has been observed in Escherichia coli. This preference for the use of certain complementary triplets as codons suggests that the evolution of the use of the genetic code depended to some extent upon the double-stranded nature of the coding material. In addition, the effect of discrimination against the use of two dinucleotides, CpG and UpA, is observed in codon usage and also in adjacent codon interference. Codons beginning with G, or A, are unlikely to be preceded by codons ending in C, or U, respectively. Consideration of codon assignment in the genetic code together with the observed CpG infrequency suggests that the evolution of the code may have been influenced by conditions in which the use of CpG dinucleotides was unfavorable. The infrequent use of UpA dinucleotides can be explained as the result of frameshift mutation during gene evolution.     

Synthetic neomycin-kanamycin phosphotransferase, type II coding sequence for gene targeting in mammalian cells

The bacterial neomycin-kanamycin phosphotransferase, type II enzyme is encoded by the neo gene and confers resistance to aminoglycoside drugs such as neomycin and kanamycin-bacterial selection and G418-eukaryotic cell selection. Although widely used in gene targeting in mouse embryonic stem cells, the neo coding sequence contains numerous cryptic splice sites and has a high CpG content. At least the former can cause unwanted effects in cis at the targeted locus. We describe a synthetic sequence, sneo, which encodes the same protein as that encoded by neo. This synthetic sequence has no predicted splice sites in either strand, low CpG content, and increased mammalian codon usage. In mouse embryonic stem cells sneo expressability is similar to neo. The use of sneo in gene targeting experiments should substantially reduce the probability of unwanted effects in cis due to splicing, and perhaps CpG methylation, within the coding sequence of the selectable marker.     

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.     

Multisequence comparisons in protein coding genes. Search for functional constraints

A very powerful method for detecting functional constraints operative in biological macromolecules is presented. This method entails performing a base permanence analysis of protein coding genes at each codon position simultaneously in different species. It calculates the degree of permanence of subregions of the gene by dividing it into segments, c codons long, counting how many sites remain unchanged in each segment among all species compared. By comparing the base permanence among several sequences with the expectations based on a stochastic evolutionary process, gene regions showing different degrees of conservation can be selected. This means that wherever the permanence deviates significantly from the expected value generated by the simulation, the corresponding regions are considered "constrained" or "hypervariable". The constrained regions are of two types: alpha and beta. The alpha regions result from constraints at the amino acid level, whereas the beta regions are those probably involved in "control" processing. The method has been applied to mitochondrial genes coding for subunit 6 of the ATPase and subunit 1 of the cytochrome oxidase in four mammalian species: human, rat, mouse, and cow. In the two mitochondrial genes a few regions that are highly conserved in all codon positions have been identified. Among these regions a sequence, common to both genes, that is complementary to a strongly conserved region of 12S rRNA has been found. This method can also be of great help in studying molecular evolution mechanisms.     

Molecular analysis of the gene of the alpha 1-antitrypsin deficiency variant, Mnichinan.

Mnichinan, a variant of alpha 1-antitrypsin (alpha 1-AT) was detected in a Japanese individual with serum alpha 1-AT deficiency (18 mg/dl), associated with aggregated alpha 1-AT molecules in the hepatocytes. Cloning and sequencing of the 10,627-bp-long region containing the Mnichinan gene and the normal M1(Val213) alpha 1-AT gene revealed all five exons of the Mnichinan gene to be identical with the M1(Val213) alpha 1-AT gene, except for two changes: a TTC trinucleotide deletion in the codon for amino acid Phe52 and a G-A substitution, by which the normal Gly148 (GGG) became Arg148 (AGG). Dot blot analysis of the polymerase chain-reaction-amplified DNA derived from the proband and other family members showed both mutations to be associated with an alpha 1-AT deficiency phenotype. Ninety-eight alpha 1-AT alleles were all negative for both changes. Comparison of the region, except for five exons between the Mnichinan and M1(Val213) genes, demonstrated one base difference in the 5' flanking region and 14 base changes in the introns. All exon-intron junctions were identical, and base changes in the 5' flanking region did not seem significant. The G-A substitution in codon 148 of the Mnichinan gene could not be responsible for the alpha 1-AT deficiency phenotype because Arg- and not Gly- was located at the corresponding position of the protein C inhibitor belonging to the serine protease inhibitor superfamily. The deletion of Phe52 may cause the newly synthesized alpha 1-AT protein to aggregate, resulting in alpha 1-AT deficiency. Comparison of the alpha 1-AT gene sequences available indicated that the C-T substitution at the CpG dinucleotide has an important role in generation of variants and nucleotide changes in the noncoding regions of the alpha 1-AT gene.     

Primate CpG islands are maintained by heterogeneous evolutionary regimes involving minimal selection

Mammalian CpG islands are key epigenomic elements that were first characterized experimentally as genomic fractions with low levels of DNA methylation. Currently, CpG islands are defined based on their genomic sequences alone. Here, we develop evolutionary models to show that several distinct evolutionary processes generate and maintain CpG islands. One central evolutionary regime resulting in enriched CpG content is driven by low levels of DNA methylation and consequentially low rates of CpG deamination. Another major force forming CpG islands is biased gene conversion that stabilizes constitutively methylated CpG islands by balancing rapid deamination with CpG fixation. Importantly, evolutionary analysis and population genetics data suggest that selection for high CpG content is not?a significant factor contributing to conservation of CpGs in differentially methylated regions. The heterogeneous, but not selective, origins of CpG islands have direct implications for the understanding of DNA methylation patterns in healthy and diseased cells.     

The rate heterogeneity of nonsynonymous substitutions in mammalian mitochondrial genes

Substitution rates at the three codon positions (r1, r2, and r3) of mammalian mitochondrial genes are in the order of r3 > r1 > r2, and the rate heterogeneity at the three positions, as measured by the shape parameter of the gamma distribution (alpha 1, alpha 2, and alpha 3), is in the order of alpha 3 > alpha 1 > alpha 2. The causes for the rate heterogeneity at the three codon positions remain unclear and, in particular, there has been no satisfactory explanation for the observation of alpha 1 > alpha 2. I attempted to dissect the causes of rate heterogeneity by studying the pattern of nonsynonymous substitutions with respect to codon positions in 10 mitochondrial genes from 19 mammalian species. Nonsynonymous substitutions involve more different amino acid replacements at the second than at the first codon position, which results in r1 > r2. The difference between r1 and r2 increases with the intensity of purifying selection, and so does the rate heterogeneity in nonsynonymous substitutions among sites at the same codon position. All mitochondrial genes appear to have functionally important and unimportant codons, with the latter having all three codon positions prone to nonsynonymous substitutions. Within the functionally important codons, the second codon position is much more conservative than the codon position. This explains why alpha 1 > alpha 2. The result suggests that overweighting of the second codon position in phylogenetic analysis may be a misguided practice.      

A third mutation at the CpG dinucleotide of codon 504 and a silent mutation at codon 506 of the HEX A gene.

Two CpG mutations at codon 504 of the gene encoding the alpha-subunit of beta-hexosaminidase (the HEX A gene) have been identified previously: (1) a C deletion resulting in premature termination of the alpha-subunit and (2) a G----A transition resulting in 504Arg----His substitution, in patients with infantile Tay-Sachs disease and juvenile GM2 gangliosidosis, respectively. This prompted a search for a C----T transition in the same dinucleotide, as would be expected from the mechanism of CpG mutagenesis. Such a mutation, which results in a substitution of cysteine for arginine, was found in a patient with chronic GM2 gangliosidosis, in compound heterozygosity with the known 269Gly----Ser allele. The biochemical phenotype of the 504Arg----Cys mutation was examined by site-directed mutagenesis of the alpha-subunit cDNA and transfection of Cos-1 cells. The expression of the mutagenized cDNA with the cysteine substitution gave rise to an alpha-subunit with the same defects as those resulting from expression of mutagenized cDNA with the histidine substitution--i.e., secretion primarily as the alpha-monomer rather than as the alpha alpha dimer, along with absence of enzymatic activity. The 504Arg----Cys/269Gly----Ser genotype of the chronic GM2 gangliosidosis patient is shared by her sibling, who has mild adult-onset GM2 gangliosidosis, implying that the clinical differences between them must be attributed to other factors. The family is unique in yet another respect--namely, that the normal allele of the mother and of a 504Arg----Cys heterozygous sibling has a silent mutation, a G----A transition in the wobble position of the glutamic acid codon at position 506.(ABSTRACT TRUNCATED AT 250 WORDS)     

Preferential carcinogen-DNA adduct formation at codons 12 and 14 in the human K-ras gene and their possible mechanisms

In the ras gene superfamily, codon 12 (-TGGTG-) of the K-ras gene is the most frequently mutated codon in human cancers. Recently, we have found that bulky chemical carcinogens preferentially form DNA adducts at codons 12 and 14 (-CGTAG-) in the K-ras gene in normal human bronchial epithelial (NHBE) cells. Furthermore, DNA adducts formed at codon 12 of the K-ras gene are poorly repaired compared with those at other codons including codon 14. These results suggest that targeted carcinogen-DNA adduct formation is a major reason for the observed high mutation frequency at codon 12 of the K-ras gene in human cancers. This preferential carcinogen-DNA adduct formation at codons 12 and 14 could result from effects of (1) primary sequences of these codons and their surrounding codons in the K-ras gene, (2) the chromatin structure, and/or (3) epigenetic factors such as C5 cytosine methylation or other DNA modifications at these codons and their surrounding codons. To distinguish these possibilities, we have introduced modifications with benzo[a]pyrene diol epoxide, N-hydroxy-2-aminofluorene, and aflatoxin B1 8,9-epoxide in (1) naked intact genomic DNA isolated from NHBE cells, (2) fragmented genomic DNA digested by restriction enzymes, and (3) in vitro synthesized DNA fragments containing the K-ras gene exon 1 sequence with or without methylation of the cytosines at CpG sites and the cytosines pairing with the guanines of codons 12 and 14. The distribution of carcinogen-DNA adducts in the K-ras gene was mapped at the nucleotide sequence level using the UvrABC nuclease incision method with or without the ligation-mediated polymerase chain reaction technique. We have found that carcinogens preferentially form adducts at codons 12 and 14 in the K-ras gene exon 1 in intact as well as in fragmented genomic DNA. In contrast, this preferential DNA adduct formation at codons 12 and 14 was not observed in PCR-amplified DNA fragments containing the K-ras gene exon 1 sequence. Methylation of the cytosine at the CpG site of codon 14, or the cytosine pairing with guanine of codon 14, greatly enhanced carcinogen-DNA adduct formation at codon 14 but did not affect carcinogen-DNA adduct formation at codon 12. Methylation of the cytosine pairing with the guanine of codon 12 also did not enhance carcinogen-DNA adduct formation at codon 12. Furthermore, we found that the cytosine at the CpG site of codon 14 is highly methylated in NHBE cells. These results suggest that cytosine methylation at the CpG site is the major reason for the preferential DNA damage at codon 14 and that epigenetic modification(s) other than cytosine methylation may contribute to the preferential DNA damage at codon 12 of the K-ras gene.     

Assignment of orthologous relationships among mammalian alpha-globin genes by examining flanking regions reveals a rapid rate of evolution

In order to study the relationships among mammalian alpha-globin genes, we have determined the sequence of the 3' flanking region of the human alpha 1 globin gene and have made pairwise comparisons between sequenced alpha-globin genes. The flanking regions were examined in detail because sequence matches in these regions could be interpreted with the least complication from the gene duplications and conversions that have occurred frequently in mammalian alpha-like globin gene clusters. We found good matches between the flanking regions of human alpha 1 and rabbit alpha 1, human psi alpha 1 and goat I alpha, human alpha 2 and goat II alpha, and horse alpha 1 and goat II alpha. These matches were used to align the alpha-globin genes in gene clusters from different mammals. This alignment shows that genes at equivalent positions in the gene clusters of different mammals can be functional or nonfunctional, depending on whether they corrected against a functional alpha-globin gene in recent evolutionary history. The number of alpha-globin genes (including pseudogenes) appears to differ among species, although highly divergent pseudogenes may not have been detected in all species examined. Although matching sequences could be found in interspecies comparisons of the flanking regions of alpha- globin genes, these matches are not as extensive as those found in the flanking regions of mammalian beta-like globin genes. This observation suggests that the noncoding sequences in the mammalian alpha-globin gene clusters are evolving at a faster rate than those in the beta-like globin gene clusters. The proposed faster rate of evolution fits with the poor conservation of the genetic linkage map around alpha-globin gene clusters when compared to that of the beta-like globin gene clusters. Analysis of the 3' flanking regions of alpha-globin genes has revealed a conserved sequence approximately 100-150 bp 3' to the polyadenylation site; this sequence may be involved in the expression or regulation of alpha-globin genes.      

CpG dinucleotides are mutation hot spots in phenylketonuria

The coding region of the phenylalanine hydroxylase (PAH) gene contains 22 CpG dinucleotides, including five doublets in the seventh exon of the gene. We hypothesized that CpG doublets could represent mutation hot spots in PAH deficiencies and we carried out the systematic sequence analysis of exon 7 in 20 unrelated PAH-deficient kindreds of Mediterranean ancestry. This procedure resulted in the detection of two novel missense mutations whose location and nature (CG to CA and CG to TG) were consistent with the accidental deamination of a 5-methylcytosine in a CpG doublet (codon 261arg----gln and codon 252arg----trp). Moreover, the codon 261 mutation was found to be associated with mutant restriction fragment length polymorphism (RFLP) haplotype 1, the most frequent mutant RFLP haplotype at the PAH locus in the studies reported thus far. However, since the mutation was detected in only 36% of haplotype 1 mutant alleles, it appears that this haplotype at the PAH locus is genotypically heterogeneous in Mediterranean countries.     

Species-specific codon context rules unveil non-neutrality effects of synonymous mutations

BackgroundCodon pair usage (codon context) is a species specific gene primary structure feature whose evolutionary and functional roles are poorly understood. The data available show that codon-context has direct impact on both translation accuracy and efficiency, but one does not yet understand how it affects these two translation variables or whether context biases shape gene evolution.ConclusionsSince in vivo studies provide evidence for a role of codon context on decoding fidelity in E. coli and for decoding efficiency in mammalian cells, our data support the hypothesis that, like codon usage, codon context modulates the evolution of gene primary structure and fine tunes the structure of open reading frames for high genome translational fidelity and efficiency in the 3 domains of life.     

Codon-improved Cre recombinase (iCre) expression in the mouse

By applying the mammalian codon usage to Cre recombinase, we improved Cre expression, as determined by immunoblot and functional analysis, in three different mammalian cell lines. The improved Cre (iCre) gene was also designed to reduce the high CpG content of the prokaryotic coding sequence, thereby reducing the chances of epigenetic silencing in mammals. Transgenic iCre expressing mice were obtained with good frequency, and in these mice loxP-mediated DNA recombination was observed in all cells expressing iCre. Moreover, iCre fused to two estrogen receptor hormone binding domains for temporal control of Cre activity could also be expressed in transgenic mice. However, Cre induction after administration of tamoxifen yielded only low Cre activity. Thus, whereas efficient activation of Cre fusion proteins in the brain needs further improvements, our studies indicate that iCre should facilitate genetic experiments in the mouse.     

Extensive methylation of a part of the CpG island located 3.0-4.5 kbp upstream to the chicken alpha-globin gene cluster may contribute to silencing the globin genes in non-erythroid cells

Here, we show that in the chicken genome, the domain of alpha-globin genes is preceded by a CpG island of which the downstream part ( approximately 0.65 kbp) is heavily methylated in lymphoid cells; it is either non-methylated or undermethylated in erythroid cells. Recombinant plasmids were constructed with the corresponding DNA fragment (called "uCpG") placed upstream to a reporter CAT gene expressed from the promoter of the alpha(D) chicken globin gene. Selective methylation of CpG dinucleotides within the uCpG fragment suppressed fivefold the expression of the CAT gene, when neither this gene itself nor the alpha(D) promoter were methylated. Methylation of CpG dinucleotides within the alpha(D) gene promoter did not modify the suppression effect exerted by methylated uCpG. We interpret these results within the frame of the hypothesis postulating, that methylation of the upstream CpG island of the chicken alpha-globin gene domain may play an essential role in silencing the alpha-globin genes in non-erythroid cells.     

CpG and TpA frequencies in the plant system.

Higher plant nuclear sequences reveal avoidance of CpG and TpA doublets. Chloroplast sequences avoid the TpA doublet in all codon positions. The chloroplast genome is not methylated but codon positions II-III and untranslated regions avoid CpG. The mitochondrial genome, also unmethylated, avoids CpG in all codon positions. We therefore deduce that methylation is not sufficient to explain CpG avoidance in the higher plant systems. Other factors must be taken into account such as amino acid composition, codon choices and perhaps stability of the DNA helix.