Nucleotide Sequence of a 28-kb Mouse Genomic Region Comprising the Imprinted Igf2 Gene

The mouse insulin-like growth factor II gene (Igf2) is physicallylinked to the insulin II gene (Ins2) and both are subject totissue-specific genomic imprinting. The paternal-specific expressionof Igf2 has been associated with hypermethylation of some CpGsites in the 5' flanking region and in the body of the gene.As a first step in analyzing the structural features of thisimprinted locus, we here report the complete nucleotide sequenceof Igf2, including all introns and the intergenic region adjacentto Ins2. This 28-kb segment of mouse chromosome 7 exhibits 80%overall identity with the corresponding rat sequence and hasa high GC content of 52%. In addition to the known CpG islandwithin the second Igf2 promoter, another island was identifiedapproximately 2 kb 5' to the first exon. Other features of thislocus include a 35-fold tandem repeat of an 11-bp sequence thatoverlaps Igf2 pseudo-exon 2, and a B2 repeat element in theintergenic region between Ins2 and Igf2. The GC-richness andthe presence of CpG islands associated with tandem repeats arecommon features of imprinted genes and thus may play a rolein the imprinting mechanism.     

Recombination drives the evolution of GC-content in the human genome

Unraveling the evolutionary forces responsible for variations of neutral substitution patterns among taxa or along genomes is a major issue in the identification of functional sequence features. Mammalian genomes show large-scale regional variations of GC-content (the isochores), but the substitution processes at the origin of this structure are poorly understood. We have analyzed the pattern of neutral substitutions in 14.3 Mb of primate noncoding regions. We show that the GC-content toward which sequences are evolving is strongly correlated (r(2) = 0.61, P isochore structure of our genome. This non-equilibrium situation suggests that changes of recombination rates occur relatively frequently during evolution, possibly as a consequence of karyotype rearrangements. These results have important implications for understanding the spatial and temporal variations of substitution processes in a broad range of sexual organisms, and for detecting the hallmarks of natural selection in DNA sequences.     

Exponential decay of GC content detected by strand-symmetric substitution rates influences the evolution of isochore structure

The distribution of guanine and cytosine nucleotides throughout a genome, or the GC content, is associated with numerous features in mammals; understanding the pattern and evolutionary history of GC content is crucial to our efforts to annotate the genome. The local GC content is decaying toward an equilibrium point, but the causes and rates of this decay, as well as the value of the equilibrium point, remain topics of debate. By comparing the results of 2 methods for estimating local substitution rates, we identify 620 Mb of the human genome in which the rates of the various types of nucleotide substitutions are the same on both strands. These strand-symmetric regions show an exponential decay of local GC content at a pace determined by local substitution rates. DNA segments subjected to higher rates experience disproportionately accelerated decay and are AT rich, whereas segments subjected to lower rates decay more slowly and are GC rich. Although we are unable to draw any conclusions about causal factors, the results support the hypothesis proposed by Khelifi A, Meunier J, Duret L, and Mouchiroud D (2006. GC content evolution of the human and mouse genomes: insights from the study of processed pseudogenes in regions of different recombination rates. J Mol Evol. 62:745-752.) that the isochore structure has been reshaped over time. If rate variation were a determining factor, then the current isochore structure of mammalian genomes could result from the local differences in substitution rates. We predict that under current conditions strand-symmetric portions of the human genome will stabilize at an average GC content of 30% (considerably less than the current 42%), thus confirming that the human genome has not yet reached equilibrium.     

Sequence analysis and compositional properties of untranslated regions of human mRNAs

A detailed computer analysis of the untranslated regions, 5′-UTR and 3′- UTR, of human mRNA sequences is reported. The compositional properties of these regions, compared with those of the corresponding coding regions, indicate that 5′-UTR and 3′-UTR are less affected by the isochore compartmentalization than the corresponding third codon positions of mRNAs. The presence of higher functional constraints in 5′-UTR is also reported. Dinucleotide analysis shows a depletion of CpG and TpA in both sequences. A search for significant sequence motifs using the WORDUP algorithm reveals the patterns already known to have a functional role in the mRNA UTR, and several other motifs whose functional roles remain to be demonstrated. This type of analysis may be particularly useful for guiding site-directed mutagenesis experiments. In addition, it can be used for assessing the nature of anonymous sequences now produced in large amounts in megabase sequencing projects.     

Interpretation of the temperature dependence of equilibrium and rate constants

The objective of this review is to draw attention to potential pitfalls in attempts to glean mechanistic information from the magnitudes of standard enthalpies and entropies derived from the temperature dependence of equilibrium and rate constants for protein interactions. Problems arise because the minimalist model that suffices to describe the energy differences between initial and final states usually comprises a set of linked equilibria, each of which is characterized by its own energetics. For example, because the overall standard enthalpy is a composite of those individual values, a positive magnitude for DeltaH(o) can still arise despite all reactions within the subset being characterized by negative enthalpy changes: designation of the reaction as being entropy driven is thus equivocal. An experimenter must always bear in mind the fact that any mechanistic interpretation of the magnitudes of thermodynamic parameters refers to the reaction model rather than the experimental system. For the same reason there is little point in subjecting the temperature dependence of rate constants for protein interactions to transition-state analysis. If comparisons with reported values of standard enthalpy and entropy of activation are needed, they are readily calculated from the empirical Arrhenius parameters.     

Evolution of the isochore structure in the scale of chromosome: insight from the mutation bias and fixation bias

Following the development of reliable methods for inferring the direction of mutations of the single nucleotide polymorphism (SNP), and the revealing of the human isochore map, it has become possible to investigate the evolution of the isochore structure in a continuous region. In this study, the recent evolution of the isochore structure on human chromosome 18, as inferred from the SNP, was examined. A remarkable mutation bias was found, which was destroying the present isochore structure. However, a fixation bias contributed by the biased gene conversion (BGC) effect and a rising fixation probability of derived alleles with increasing GC content was extending the present isochore structure. Combining the two opposing processes, the old isochore structure was declining and a more homogenous isochore structure with higher GC content was being formed on the chromosome. During this process, both the CpG and genic sites, which were present in the isochore but were paid little attention to before, played an important role. In addition, the recombination was confirmed to promote the GC alleles fixed in the genome because of the BGC effect. For the first time, it was observed that with the occurrence of little recombination, AT alleles had the identical fixation probability with GC alleles in the recombination cold spots.     

Compositional evolution of noncoding DNA in the human and chimpanzee genomes

We have examined the compositional evolution of noncoding DNA in the primate genome by comparison of lineage-specific substitutions observed in 1.8 Mb of genomic alignments of human, chimpanzee, and baboon with 6542 human single-nucleotide polymorphisms (SNPs) rooted using chimpanzee sequence. The pattern of compositional evolution, measured in terms of the numbers of GC-->AT and AT-->GC changes, differs significantly between fixed and polymorphic sites, and indicates that there is a bias toward fixation of AT-->GC mutations, which could result from weak directional selection or biased gene conversion in favor of high GC content. Comparison of the frequency distributions of a subset of the SNPs revealed no significant difference between GC-->AT and AT-->GC polymorphisms, although AT-->GC polymorphisms in regions of high GC segregate at slightly higher frequencies on average than GC-->AT polymorphisms, which is consistent with a fixation bias favoring high GC in these regions. However, the substitution data suggest that this fixation bias is relatively weak, because the compositional structure of the human and chimpanzee genomes is becoming homogenized, with regions of high GC decreasing in GC content and regions of low GC increasing in GC content. The rate and pattern of nucleotide substitution in 333 Alu repeats within the human-chimpanzee-baboon alignments are not significantly affected by the GC content of the region in which they are inserted, providing further evidence that, since the time of the human-chimpanzee ancestor, there has been little or no regional variation in mutation bias.     

Strong regional biases in nucleotide substitution in the chicken genome

Interspersed repeats have emerged as a valuable tool for studying neutral patterns of molecular evolution. Here we analyze variation in the rate and pattern of nucleotide substitution across all autosomes in the chicken genome by comparing the present-day CR1 repeat sequences with their ancestral copies and reconstructing nucleotide substitutions with a maximum likelihood model. The results shed light on the origin and evolution of large-scale heterogeneity in GC content found in the genomes of birds and mammals--the isochore structure. In contrast to mammals, where GC content is becoming homogenized, heterogeneity in GC content is being reinforced in the chicken genome. This is also supported by patterns of substitution inferred from alignments of introns in chicken, turkey, and quail. Analysis of individual substitution frequencies is consistent with the biased gene conversion (BGC) model of isochore evolution, and it is likely that patterns of evolution in the chicken genome closely resemble those in the ancestral amniote genome, when it is inferred that isochores originated. Microchromosomes and distal regions of macrochromosomes are found to have elevated substitution rates and a more GC-biased pattern of nucleotide substitution. This can largely be accounted for by a strong correlation between GC content and the rate and pattern of substitution. The results suggest that an interaction between increased mutability at CpG motifs and fixation biases due to BGC could explain increased levels of divergence in GC-rich regions.     

Heterogeneity in regional GC content and differential usage of codons and amino acids in GC-poor and GC-rich regions of the genome of Apis mellifera

The honeybee (Apis mellifera) has a genome with a wide variation in GC content showing 2 clear modal GC values, in some ways reminiscent of an isochore-like structure. To gain insight into causes and consequences of this pattern, we used a comparative approach to study the genome-wide alignment of primarily coding sequence of A. mellifera with Drosophila melanogaster and Anopheles gambiae. The latter 2 species show a higher average GC content than A. mellifera and no indications of bimodality, suggesting that the GC-poor mode is a derived condition in honeybee. In A. mellifera, synonymous sites of genes generally adopt the GC content of the region in which they reside. A large proportion of genes in GC-poor regions have not been assigned to the honeybee assembly because of the low sequence complexity of their genome neighborhood. The synonymous substitution rate between A. mellifera and the other species is very close to saturation, but analyses of nonsynonymous substitutions as well as amino acid substitutions indicate that the GC-poor regions are not evolving faster than the GC-rich regions. We describe the codon usage and amino acid usage and show that they are remarkably heterogeneous within the honeybee genome between the 2 different GC regions. Specifically, the genes located in GC-poor regions show a much larger deviation in both codon usage bias and amino acid usage from the Dipterans than the genes located in the GC-rich regions.