Nucleic Acids Research, 2003, Vol. 31, No. 7 1838-1844
© 2003 Oxford University Press
DNA helix: the importance of being GC-rich
Institute of Cytology, Russian Academy of Sciences, Tikhoretsky Avenue 4, St Petersburg 194064, Russia
*Tel: +7 812 247 5310; Fax: +7 812 247 0341; Email: aevin{at}mail.cytspb.rssi.ru
Received January 21, 2003; Revised and Accepted February 12, 2003
| ABSTRACT |
|---|
|
|
|---|
A new explanation for the emergence of heavy (GC-rich) isochores is proposed, based on the study of thermostability, bendability, ability to BZ transition and curvature of the DNA helix. The absolute values of thermostability, bendability and ability to BZ transition correlated positively with GC content, whereas curvature correlated negatively. The relative values of these parameters were determined as compared to randomized sequences. In genes and intergenic spacers of warm-blooded animals, both the relative bendability and ability to BZ transition increased with elevation of GC content, whereas the relative thermostability and curvature decreased. The usage of synonymous codons in GC-rich genes was also found to augment bendability and ability to BZ transition and to reduce thermostability of DNA (as compared to synonymous codons with the same GC content). The analysis of transposable elements (Alu and B2 repeats in the human and mouse) showed that the level of their divergence from the consensus sequence positively correlated with relative bendability and ability to BZ transition and negatively with relative thermostability. The bendability and ability to BZ transition are known to relate to open chromatin and active transcription, whereas curvature facilitates chromatin condensation. Because heavy isochores are known to be gene-rich and show a high level of transcription, it is suggested here that isochores arose not as an adaptation to elevated temperature but because of a certain grade of general organization and correspondingly advanced level of genomic organization, reflected in genome structuring, with physical properties of DNA in the gene-rich regions being optimized for active transcription and in the gene-poor regions for chromatin condensation (transcription/grade concept).
| INTRODUCTION |
|---|
|
|
|---|
Organisms which form the crowns of the animal and plant kingdoms, warm-blooded animals and monocot plants, show strong base compositional heterogeneity of genomes. They consist of relatively homogeneous regions (>300 kb long), called isochores, which differ in the GC content and gene concentration (15). It is the heavy (GC-rich and gene-rich) regions which have evolved lately and thus brought about the heterogeneity of the genome as a whole (24,6). There are different views on the origin of heavy isochores. According to the neutralist viewpoint, they were formed by variation of mutation bias (in a broad sense, including gene conversion) among different genomic regions (712). The selectionist hypothesis assumes that their emergence was caused by a need for a higher thermostability of the DNA helix (3,13). However, the case of prokaryotes where there is no correlation between GC content and habitat temperature (notwithstanding a very wide range of temperatures with the upper margin far exceeding 40°C) undermines the thermal hypothesis (14,15). This conclusion was supported by the latest studies on ectothermic (cold-blooded) vertebrates where no correlation between GC content and adaptation temperature was observed, with the temperature ranging from 2 to +45°C (16,17). Recently, it was shown that with elevation of GC content, the thermostability of genes in warm-blooded animals increased more slowly than in random sequences, whereas their bendability rose faster (18). This suggests that, if the emergence of heavy isochores was indeed due to selection, it was selection not for thermostability but for bendability of the DNA molecule, which may be related to active transcription in these gene-rich genomic regions.
Here this hypothesis is examined in several aspects. First, two other physical properties of the DNA helix (ability to BZ transition and curvature), which are known to relate to transcription and chromatin condensation (1922), were analyzed. Second, the intergenic spacers and the usage of synonymous codons in coding sequences were studied. Third, the evolution of transposable elements (TEs) was explored with regard to a change in the physical properties of DNA.
| MATERIALS AND METHODS |
|---|
|
|
|---|
The sequences of nuclear genes were extracted from GenBank. Only genes with complete coding sequences (CDSs) were taken. Genes were checked for duplicates on the basis of CDS similarity (>99%). All coding sequences were checked for the absence of internal stop codons to ensure they were not pseudogenes. The exon/intron boundaries were taken from annotations. Pseudogenes were not included in the gene level analysis but they were taken into account for determination of intergenic spacers. Only sequences where each nucleotide was known for at least 99% of positions were used. The intergenic spacers were studied only in the (better presented) human, mouse and rice genomes. The gene level analysis was also done for the chicken and maize.
The TEs in introns and intergenic spacers were determined using the RepeatMasker program (A.F.A.Smit and P.Green, http://ftp.genome.washington.edu/RM/RepeatMasker.html, with parameter -gccalc, without masking simple repeats and low complexity regions). The level of divergence (percent of nucleotide substitutions) with regard to consensus sequence for each repeat family, the libraries of which were based on the expanded versions of Repbase (23), were also calculated with this program. The transposable elements were studied in the human, mouse and rice genomes.
The bendability and curvature of nucleotide sequences were determined using the trinucleotide table of consensus values obtained from DNase I digestion and nucleosome positioning studies (24,25). The bendability is a local parameter representing the ability of DNA to bend (usually towards the major groove) as a result of thermal fluctuations or DNAprotein interactions (26). It was calculated in a sliding trinucleotide frame (with a 1 nt step), and averaged for each sequence. The curvature is a relatively macroscopic DNA bend, representing the intrinsic tendency of DNA to follow a non-linear pathway over an appreciable length, which is a result of variation of local bends in phase with the DNA helix (26). It is usually calculated for DNA segments of arbitrary length, which is a multiple of the number of helix turns (24,26). Different numbers of turns (one to three) were tried here; the results did not differ qualitatively. Therefore, calculations made with the shortest DNA segment (one turn, 10.5 bp) were chosen to minimize possible margin effects for short sequences. Thus, curvature was calculated in a sliding 10.5 nt frame (with a 1 nt step), and averaged for each sequence. The thermostability was determined using a unified dinucleotide table for free energy of melting (
G) (27), obtained from different physical studies, mainly the UV absorbance and calorimetric temperature profiles. The ability to BZ transition was determined using the dinucleotide table for free energy of BZ transition (the lower the free energy, the higher the ability to transition), obtained by an advanced modification of the JUMNA method (28). The thermostability and ability to BZ transition were calculated in a sliding dinucleotide frame (with a 1 nt step), and averaged for each sequence.
For each genomic sequence (exons, introns, intergenic spacers and TEs), 10 randomizations were made and the parameters under study, calculated for these randomized sequences, were averaged. For each protein coding sequence (CDS), all synonymous codon positions, which can be permuted with conservation of the same GC content of each synonymous position (i.e. only G
C and A
T replacements), were randomly permuted 10 times with preservation of the mean purine content of the total set of synonymous positions, and parameters calculated for these permuted sequences were averaged. Each CDS was taken as a whole sequence, therefore curvature was not determined because of its relatively long calculation frame, which could produce considerable errors on exon junctions. The relative value of each physical parameter under study was determined for each genomic sequence as the difference between the value for the genomic sequence itself and the average for its 10 randomized or permuted sequences. (The relative parameters determined as ratios between the values for genomic and randomized/permuted sequences were also analyzed and showed qualitatively the same results.)
The conventional statistical analyses were done with the Statgraphics (Statistical Graphics Co.) and Statistica (StatSoft Inc.) software.
| RESULTS |
|---|
|
|
|---|
Absolute values
For any type of DNA sequences (genomic or randomized) studied here, the bendability, ability to BZ transition and thermostability increased with elevation of GC content, whereas curvature decreased. Human intergenic spacers are presented as an example (Table 1). The bendability and ability to BZ transition always correlated positively and strongly between themselves, even independently of their link to GC content (Table 1). The thermostability and curvature usually correlated negatively (and more weakly) with bendability and ability to BZ transition if GC content was fixed (Table 1).
|
Relative values in introns, exons and intergenic spacers
In the human, mouse and chicken, the relative bendability and ability to BZ transition increased with elevation of GC percent in all types of genomic sequences (Table 2 and Fig. 1). The relative thermostability and curvature showed quite the opposite trend. All correlations were stronger in the introns as compared to exons and intergenic spacers (Table 2). In exons of rice and maize, the three parameters (relative bendability, ability to BZ transition and thermostability) increased with elevation of GC content, whereas relative curvature decreased (Table 2). For introns and intergenic spacers of monocots, the data were inconsistent.
|
|
Synonym-relative values in coding sequences
In the human, mouse and chicken, the synonym-relative values, obtained by permutation of synonymous codon positions with preservation of GC content of each position, behave in a similar way to the relative values obtained by complete randomization of their genomic sequences. The synonym-relative bendability and ability to BZ transition increased with elevation of GC percent, whereas the synonym-relative thermostability decreased (Table 3 and Fig. 2). In the monocots, the synonym-relative values of coding sequences behave similarly with the relative values for exons, i.e. bendability, ability to BZ transition and thermostability increased with elevation of GC content (Table 3 and Fig. 2).
|
|
Relative values in the transposable elements
The most numerous transposable elements in the human and mice genomes are Alu and B2 repeats. The level of their divergence from the consensus sequence positively correlated with the relative bendability and ability to BZ transition, and negatively with relative thermostability (Table 4). The most numerous transposable elements in rice (among those presented in Repbase) were the TNR1 and Wanderer repeats. With the growth of divergence from the consensus sequence, the TNR1 repeats showed an increase in three parameters, bendability, ability to BZ transition and thermostability (Table 4). The Wanderer repeats showed a picture similar to transposons of warm-blooded vertebrates, i.e. an increase in relative bendability and ability to BZ transition and a decrease in relative thermostability (Table 4).
|
| DISCUSSION |
|---|
|
|
|---|
The obtained results indicate that, if heavy isochores of warm-blooded vertebrates were indeed formed by selection for some physical properties of DNA, it was selection for traits associated with transcription (higher bendability and ability to BZ transition and lower curvature) rather than for thermostability. In the case of selection not for GC content itself but for a certain physical property related to it (which is quite reasonable because non-coding sequences of heavy isochores are also involved in the elevation of GC content), the enhancement of this physical property should overtake the increase in GC content. Therefore, the statistically significant deviation from random can serve as an indicator of this selective force. The conclusion about the increase in relative bendability and ability to BZ transition with elevation of GC content pertains not only to introns and intergenic spacers (and to the evolution of TEs located in them) but also to synonymous positions in protein-coding DNA as well. An unequal usage of synonymous codons is a long-standing problem (see for example 2931). In some organisms, it is associated with the level of gene expression and can be explained by relative abundance of corresponding isoacceptor tRNAs (29,32,33). However, in the warm-blooded vertebrates the causes of unequal codon usage remain unclear (29,30,34). In these organisms, codon usage is generally related to isochore structure (2,3,35). The data obtained here show that the usage of codons in the GC-rich genes of warm-blooded vertebrates augments bendability and ability to BZ transition of DNA (even as compared to synonymous codons with the same GC content).
With its non-neutralized electric charges, the DNA helix behaves as a rigid, steel-like rod, and even after its neutralization with histones, bending of such a wire into the orderly, compact chromatin structure remains a challenging task (36,37). In the multicellular eukaryotes, the main problem of the genetic machinery is to control repression and activation of genes, which is connected with bending and unbending of the DNA molecule. Although the GC-rich regions make up only a minor part of the mammalian genome (1015%), gene concentration in the heaviest isochore family (H3) is 20-fold higher than in the GC-poor regions (2,6). The heavy isochores constitute early replicated, actively transcribed chromatin (38,39). The expression level of genes located in the heavy isochores is on average higher (40). It is known that an increase in DNA bendability promotes the transition from nucleosome-wrapped to extended DNA (21,22), thereby it should facilitate transcription. In the warm-blooded vertebrates, all correlations found here are stronger in introns as compared to exons and intergenic spacers (Table 2). This is in agreement with the proposed concept, because exons are under strong selection for informational content, while intergenic spacers should be involved in transcription-related decondensation to a lesser extent as compared to genes.
The ability to BZ transition is also linked very tightly to transcription. It is known that antibodies to Z-DNA bind preferentially to actively transcribed genes (20). The association between the amount of Z-DNA and transcription activity has even been shown for individual genes (c-myc, corticotropin releasing-hormone gene). When the gene is up-regulated, the level of Z-form increases, when it is down-regulated, it decreases (20). At the same time, the level of Z-DNA is not affected by replication. A wave of Z-DNA is formed after the moving RNA polymerase, which may prevent the following RNA polymerase from transcribing that region of a gene (19,20). It should also facilitate untwisting of the DNA helix, absorbing supercoils and thereby buffering torsional stress generated by RNA polymerase until this stress is relieved by topoisomerases. In any case, the strong association between Z-DNA and transcription is well established.
In monocot plants, the picture is more complicated and selection for thermostability cannot be excluded. The coding DNA might have been selected both for bendability/ability to BZ transition and thermostability. Different TEs probably located predominantly in different genomic segments evolved either similarly to monocot CDSs (i.e. in the direction of both bendability/ability to BZ transition and thermostability) or similarly to the genomes of warm-blooded vertebrates (in the direction of bendability/ability to BZ transition only). It is noteworthy that although exons of GC-rich genes usually show a higher GC content as compared to corresponding introns in both groups of organisms, the exonintron contrast in GC content is much higher in monocots than in warm-blooded vertebrates (41,42). The same picture was observed for the coding DNAintergenic spacer contrast and for the silent codon positionnon-coding DNA contrast (Fig. 3). This can explain why non-coding DNA behaves more similarly to exons in warm-blooded vertebrates as compared to monocots (Table 2). The heavy isochores of monocots are not the same as in warm-blooded vertebrates: in the former they are more heterogeneous in GC content [although they have a similar base compositional smoothness on the intronexon margins (42)].
|
In general, the GC-rich DNA seems to be a good DNA. In the highest multicellular organisms, it constitutes gene-rich, actively transcribed genomic regions. It is thermostable, bendable and able to BZ transition. Then why are there the GC-poor parts of the genome and why are introns and intergenic spacers in the heavy isochores still GC-poorer than exons? The substitution pressure in the human genome is generally AT-biased (4,5,12,43). This helps explain why the most part of the genome is GC-poor. The recombination rate is much higher in the heavy isochores (4446), which should facilitate selection for a good DNA in these regions. One possible cause for a lower GC content in non-coding DNA is bombardment by transposable elements, since purifying selection against the integration of TEs should be much stronger in the coding DNA (47). The functional significance of a lower GC content can be related to a dual role of curvature and bendability in chromatin condensation. The curved DNA segments serving as primers for DNAhistone interaction are necessary for nucleosome formation (21,22). Since curvature is a long-range parameter (as compared to bendability), this property cannot be furnished in exons which are under strong pressure for informational content (and where silent positions do not constitute a continuous sequence). Therefore, it is introns which should provide the necessary segments of curvature. Several cases have been described where the experimental removal of an intron resulted in a failure of nucleosome formation (48,49). Bendability and curvature, which are shown here to correlate disparately with the GC content, affect chromatin condensation in the opposite manner. Bendability is more prone to destabilize nucleosomes, whereas curvature acts in the opposite direction (21,22). In the GC-poor isochores, which constitute the late replicated (i.e. rarely decondensed) chromatin, there should be a greater need in the curved DNA segments for a tight chromatin condensation, as compared to GC-rich regions.
The concept fostered here suggests that GC-rich isochores (and hence the compositional heterogeneity of the genome as a whole) arose in warm-blooded animals not as an adaptation to elevated temperature but because of the high level of their general organization and the corresponding level of genomic organization, with the physical properties of DNA in the gene-rich regions being optimized for active transcription and in the gene-poor regions for chromatin condensation. Thus, the structuring of the genome reflects its advanced organization.
The less pronounced isochores in rodents as compared to other mammals (5052) can be related to a lower organization of these r-selected, short-living animals [which is expressed, for instance, in a relatively smaller brain (53)]. It is known that rodents have a less efficient DNA repair (54,55), which suggests that precise genomic organization is not as important for them as for other mammals. Similarly, some lesser pronounced isochoric structure probably exists in the cold-blooded vertebrates (16,56), which may reflect the respective levels of their general organization.
For the isochoric structure to emerge, selection need not necessarily be on each nucleotide (especially in non-coding DNA) but could be selection of DNA-handling proteins, which in their turn provide an appropriate mutation bias. Recently, support was found for a biased gene conversion (BGC) explanation for elevation of GC content in highly recombining regions of the mammalian genome (11,5761). [At least some of these facts can be equally interpreted as a result of nucleotide selection (11,5760), albeit this would suggest a high genetic load (11).] More exactly, the elevation of GC content is assumed to be due to biased DNA repair after gene conversion (11,59), i.e. gene conversion simply plays a role of mutation. Recombination, which is the main cause of gene conversion, was shown to be generally mutagenic (62). One of the problems for a purely neutralist explanation of the origin of GC-rich isochores (e.g. due to BGC) is the increasing deviation of DNA physical properties from randomized sequences with the elevation of GC content, revealed here and in a previous work (18). This directional deviation is due to the sequence specificity of GC-rich regions, which suggests that even if it indeed arose due to some mutation/repair bias, a whole complex of mutation/repair pathways should be involved in this process. The proteins belonging to this complex could be under selection to provide an appropriate bias in different genomic regions. If so, such indirect selection should probably have an inertia, which could drive the process beyond the optimum. This may explain the recent proposal that in mammals (especially rodents) the isochores are lately deteriorating (35,59,63).
Whatever the actual mechanisms, it is the final result that matters. The independent origin of strongly pronounced isochores in mammals and birds (2), even not considering monocots, and their functional significance upheld here suggest that the isochoric structure of the genome is a grade trait, and even if actual mechanisms do involve mutation/repair bias, the emergence of isochores was not accidental.
| ACKNOWLEDGEMENTS |
|---|
I am grateful to Giorgio Bernardi and two anonymous reviewers for helpful comments. This work was supported by the Russian Foundation for Basic Research (RFBR).
| REFERENCES |
|---|
|
|
|---|
- Bernardi,G., Olofsson,B., Filipski,J., Zerial,M., Salinas,J., Cuny,G., Meunier-Rotival,M. and Rodier,F. (1985) The mosaic genome of warm-blooded vertebrates. Science, 228, 953958.
[Abstract/Free Full Text] - Bernardi,G. (2000) Isochores and the evolutionary genomics of vertebrates. Gene, 241, 317.[CrossRef][Web of Science][Medline]
- Bernardi,G. (2000) The compositional evolution of vertebrate genomes. Gene, 259, 3143.[CrossRef][Web of Science][Medline]
- Bernardi,G. (2001) Misunderstandings about isochores. Part 1. Gene, 276, 313.[CrossRef][Web of Science][Medline]
- Eyre-Walker,A. and Hurst,L.D. (2001) The evolution of isochores. Nature Rev. Genet., 2, 549555.
- DOnofrio,G., Jabbari,K., Musto,H., Alvarez-Valin,F., Cruveiller,S. and Bernardi,G. (1999) Evolutionary genomics of vertebrates and its implications. Ann. N. Y. Acad. Sci., 870, 8194.[CrossRef][Web of Science][Medline]
- Wolfe,K., Sharp,P.M. and Li,W.H. (1989) Mutation rates differ among regions of the mammalian genome. Nature, 337, 441456.
- Francino,M.P. and Ochman,H. (1999) Isochores result from mutation not selection. Nature, 400, 3031.[CrossRef][Medline]
- Fryxell,K.J. and Zuckerkandl,E. (2000) Cytosine deamination plays a primary role in the evolution of mammalian isochores. Mol. Biol. Evol., 17, 13711383.
[Abstract/Free Full Text] - Sueoka,N. and Kawanishi,Y. (2000) DNA G+C content of the third codon position and codon usage biases of human genes. Gene, 261, 5362.[CrossRef][Web of Science][Medline]
- Galtier,N., Piganeau,G., Mouchiroud,D. and Duret,L. (2001) GC content evolution in mammalian genomes, the biased gene conversion hypothesis. Genetics, 159, 907911.
[Free Full Text] - International Human Genome Sequencing Consortium (2001) Initial sequencing and analysis of the human genome. Nature, 409, 860921.[CrossRef][Medline]
- Bernardi,G. and Bernardi,G. (1986) Compositional constraints and genome evolution. J. Mol. Evol., 24, 111.[CrossRef][Web of Science][Medline]
- Galtier,N. and Lobry,J.R. (1997) Relationships between genomic G+C content, RNA secondary structures and optimal growth temperature in prokaryotes. J. Mol. Evol., 44, 632636.[CrossRef][Web of Science][Medline]
- Hurst,L.D. and Merchant,A.R. (2001) High guaninecytosine content is not an adaptation to high temperature, a comparative analysis amongst prokaryotes. Proc. R. Soc. Lond. B Biol. Sci., 268, 493497.[Medline]
- Belle,E.M., Smith,N. and Eyre-Walker,A. (2002) Analysis of the phylogenetic distribution of isochores in vertebrates and a test of the thermal stability hypothesis. J. Mol. Evol., 55, 356363.[CrossRef][Web of Science][Medline]
- Ream,R.A., Johns,G.C. and Somero,G.N. (2003) Base compositions of genes encoding alpha-actin and lactate dehydrogenase-A from differently adapted vertebrates show no temperature-adaptive variation in G + C content. Mol. Biol. Evol., 20, 105110.
[Abstract/Free Full Text] - Vinogradov,A.E. (2001) Bendable genes of warm-blooded vertebrates. Mol. Biol. Evol., 18, 21952200.
[Abstract/Free Full Text] - Herbert,A. and Rich,A. (1996) The biology of left-handed Z-DNA. J. Biol. Chem., 271, 1159511598.
[Free Full Text] - Herbert,A. and Rich,A. (1999) Left-handed Z-DNA, structure and function. Genetica, 106, 3747.[CrossRef][Web of Science][Medline]
- Anselmi,C., Bocchinfuso,G., De Santis,P., Savino,M. and Scipioni,A. (1999) Dual role of DNA intrinsic curvature and flexibility in determining nucleosome stability. J. Mol. Biol., 286, 12931301.[CrossRef][Web of Science][Medline]
- Anselmi,C., Bocchinfuso,G., De Santis,P., Savino,M. and Scipioni,A. (2000) A theoretical model for the prediction of sequence-dependent nucleosome thermodynamic stability. Biophys. J., 79, 601613.[Web of Science][Medline]
- Jurka,J. (2000) Repbase update, a database and an electronic journal of repetitive elements. Trends Genet., 16, 418420.[CrossRef][Web of Science][Medline]
- Gabrielian,A. and Pongor,S. (1996) Correlation of intrinsic DNA curvature with DNA property periodicity. FEBS Lett., 393, 6568.[CrossRef][Web of Science][Medline]
- Munteanu,M.G., Vlahovicek,K., Parthasarathy,S., Simon,I. and Pongor,S. (1998) Rod models of DNA, sequence-dependent anisotropic elastic modelling of local bending phenomena. Trends Biochem. Sci., 23, 341347.[CrossRef][Web of Science][Medline]
- Goodsell,D.S. and Dickerson,R.E. (1994) Bending and curvature calculations in B-DNA. Nucleic Acids Res., 22, 54975503.
[Free Full Text] - SantaLucia,J.,Jr (1998) A unified view of polymer, dumbbell and oligonucleotide DNA nearest-neighbor thermodynamics. Proc. Natl Acad. Sci. USA, 95, 14601465.
[Abstract/Free Full Text] - Lafontaine,I. and Lavery,R. (2000) Optimization of nucleic acid sequences. Biophys. J., 79, 680685.[Web of Science][Medline]
- Sharp,P.M., Averof,M., Lloyd,A.T., Matassi,G. and Peden,J.F. (1995) DNA sequence evolution, the sounds of silence. Phil. Trans. R. Soc. Lond. B Biol. Sci., 349, 241247.[Web of Science][Medline]
- Urrutia,A.O. and Hurst,L.D. (2001) Codon usage bias covaries with expression breadth and the rate of synonymous evolution in humans, but this is not evidence for selection. Genetics, 159, 11911199.
[Abstract/Free Full Text] - Vinogradov,A.E. (2001) Intron length and codon usage. J. Mol. Evol., 52, 25.[Web of Science][Medline]
- Xia,X. (1996) Maximizing transcription efficiency causes codon usage bias. Genetics, 144, 13091320.[Abstract]
- Chiapello,H., Lisacek,F., Caboche,M. and Henaut,A. (1998) Codon usage and gene function are related in sequences of Arabidopsis thaliana. Gene, 209, GC1GC38.[CrossRef][Web of Science][Medline]
- Hurst,L.D. and Williams,E.J. (2000) Covariation of GC content and the silent site substitution rate in rodents, implications for methodology and for the evolution of isochores. Gene, 261, 107114.[CrossRef][Web of Science][Medline]
- Galtier,N. and Mouchiroud,D. (1998) Isochore evolution in mammals, a human-like ancestral structure. Genetics, 150, 15771584.
[Abstract/Free Full Text] - Maher,L.J.,III (1998) Mechanisms of DNA bending. Curr. Opin. Chem. Biol., 2, 688694.[CrossRef][Web of Science][Medline]
- Williams,L.D. and Maher,L.J.,III (2000) Electrostatic mechanisms of DNA deformation. Annu. Rev. Biophys. Biomol. Struct., 29, 497521.[CrossRef][Web of Science][Medline]
- Federico,C., Saccone,S. and Bernardi,G. (1998) The gene-richest bands of human chromosomes replicate at the onset of the S-phase. Cytogenet. Cell Genet., 80, 8388.[CrossRef][Web of Science][Medline]
- Saccone,S., Federico,C., Solovei,I., Croquette,M.F., Della Valle,G. and Bernardi,G. (1999) Identification of the gene-richest bands in human prometaphase chromosomes. Chromosome Res., 7, 379386.[CrossRef][Web of Science][Medline]
- Konu,O.O. and Li,M.D. (2002) Correlations between mRNA expression levels and GC contents of coding and untranslated regions of genes in rodents. J. Mol. Evol., 54, 3541.[CrossRef][Web of Science][Medline]
- Carels,N., Hatey,P., Jabbari,K. and Bernardi,G. (1998) Compositional properties of homologous coding sequences from plants. J. Mol. Evol., 46, 4553.[CrossRef][Web of Science][Medline]
- Vinogradov,A.E. (2001) Within-intron correlation with base composition of adjacent exons in different genomes. Gene, 276, 143151.[CrossRef][Web of Science][Medline]
- Petrov,D.A. and Hartl,D.L. (1999) Patterns of nucleotide substitution in Drosophila and mammalian genomes. Proc. Natl Acad. Sci. USA, 96, 14751479.
[Abstract/Free Full Text] - Eisenbarth,I., Striebel,A.M., Moschgath,E., Vogel,W. and Assum,G. (2001) Long-range sequence composition mirrors linkage disequilibrium pattern in a 1.13 Mb region of human chromosome 22. Hum. Mol. Genet., 10, 28332839.
[Abstract/Free Full Text] - Hurst,L.D. (2000) Crossing over a boundary. Trends Genet., 16, 542.[CrossRef][Medline]
- Eisenbarth,I., Vogel,G., Krone,W., Vogel,W. and Assum,G. (2000) An isochore transition in the NF1 gene region coincides with a switch in the extent of linkage disequilibrium. Am. J. Hum. Genet., 67, 873880.[CrossRef][Web of Science][Medline]
- Duret,L. and Hurst,L.D. (2001) The elevated G and C content at exonic third sites is not evidence against neutralist models of isochore evolution. Mol. Biol. Evol., 18, 757762.
[Abstract/Free Full Text] - Lauderdale,J.D. and Stein,A. (1992) Introns of the chicken ovalbumin gene promote nucleosome alignment in vitro. Nucleic Acids Res., 20, 65896596.
[Abstract/Free Full Text] - Liu,K., Sandgren,E.P., Palmiter,R.D. and Stein,A. (1995) Rat growth hormone gene introns stimulate nucleosome alignment in vitro and in transgenic mice. Proc. Natl Acad. Sci. USA, 92, 77247728.
[Abstract/Free Full Text] - Sabeur,G., Macaya,G., Kadi,F. and Bernardi,G. (1993) The isochore patterns of mammalian genomes and their phylogenetic implications. J. Mol. Evol., 37, 93108.[CrossRef][Web of Science][Medline]
- Robinson,M., Gautier,C. and Mouchiroud,D. (1997) Evolution of isochores in rodents. Mol. Biol. Evol., 14, 823828.[Abstract]
- Douady,C., Carels,N., Clay,O., Catzeflis,F. and Bernardi,G. (2000) Diversity and phylogenetic implications of CsCl profiles from rodent DNAs. Mol. Phylogenet. Evol., 17, 219230.[CrossRef][Web of Science][Medline]
- Lewin,R. (1988) Living in the fast track makes for small brains. Science, 242, 513514.
[Free Full Text] - Adelman,R., Saul,R.L. and Ames,B.N. (1988) Oxidative damage to DNA: relation to species metabolic rate and life span. Proc. Natl Acad. Sci. USA, 85, 27062708.
[Abstract/Free Full Text] - Ames,B.N. (1989) Endogenous DNA damage as related to cancer and aging. Mutat. Res., 214, 4146.[CrossRef][Web of Science][Medline]
- Hughes,S., Zelus,D. and Mouchiroud,D. (1999) Warm-blooded isochore structure in Nile crocodile and turtle. Mol. Biol. Evol., 16, 15211527.[Abstract]
- Eyre-Walker,A. (1999) Evidence of selection on silent site base composition in mammals: potential implications for the evolution of isochores and junk DNA. Genetics, 152, 675683.
[Abstract/Free Full Text] - Smith,N.G. and Eyre-Walker,A. (2001) Synonymous codon bias is not caused by mutation bias in human. Mol. Biol. Evol., 18, 982986.
[Abstract/Free Full Text] - Duret,L., Semon,M., Piganeau,G., Mouchiroud,D. and Galtier,N. (2002) Vanishing GC-rich isochores in mammalian genomes. Genetics, 162, 18371847.
[Abstract/Free Full Text] - Lercher,M.J., Smith,N.G., Eyre-Walker,A. and Hurst,L.D. (2002) The evolution of isochores. Evidence from snp frequency distributions. Genetics, 162, 18051810.
[Abstract/Free Full Text] - Galtier,N. (2003) Gene conversion drives GC content evolution in mammalian histones. Trends Genet., 19, 6568.[CrossRef][Web of Science][Medline]
- Lercher,M.J. and Hurst,LD. (2002) Human SNP variability and mutation rate are higher in regions of high recombination. Trends Genet., 18, 337340.[CrossRef][Web of Science][Medline]
- Smith,N.G. and Eyre-Walker,A. (2002) The compositional evolution of the murid genome. J. Mol. Evol., 55, 197201.[CrossRef][Web of Science][Medline]
This article has been cited by other articles:
![]() |
T. Zhou, M. Weems, and C. O. Wilke Translationally Optimal Codons Associate with Structurally Sensitive Sites in Proteins Mol. Biol. Evol., July 1, 2009; 26(7): 1571 - 1580. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. V. Kharchenko, C. J. Woo, M. Y. Tolstorukov, R. E. Kingston, and P. J. Park Nucleosome positioning in human HOX gene clusters Genome Res., October 1, 2008; 18(10): 1554 - 1561. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Khuu, M. Sandor, J. DeYoung, and P. S. Ho Phylogenomic analysis of the emergence of GC-rich transcription elements PNAS, October 16, 2007; 104(42): 16528 - 16533. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. T. Webster, E. Axelsson, and H. Ellegren Strong Regional Biases in Nucleotide Substitution in the Chicken Genome Mol. Biol. Evol., June 1, 2006; 23(6): 1203 - 1216. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Semon, J. R. Lobry, and L. Duret No Evidence for Tissue-Specific Adaptation of Synonymous Codon Usage in Humans Mol. Biol. Evol., March 1, 2006; 23(3): 523 - 529. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. E Vinogradov and O. V Anatskaya Genome size and metabolic intensity in tetrapods: a tale of two lines Proc R Soc B, January 7, 2006; 273(1582): 27 - 32. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. E. Vinogradov Noncoding DNA, isochores and gene expression: nucleosome formation potential Nucleic Acids Res., January 26, 2005; 33(2): 559 - 563. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Zhang, S. Kasif, C. R. Cantor, and N. E. Broude GC/AT-content spikes as genomic punctuation marks PNAS, November 30, 2004; 101(48): 16855 - 16860. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. E. Vinogradov Isochores and tissue-specificity Nucleic Acids Res., September 1, 2003; 31(17): 5212 - 5220. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||







