Compositional Properties of Homologous Coding Sequences from Plants

In this work, we investigated (1) the compositional distributions of all available nuclear coding sequences (and of their three codon positions) of six dicots and four Gramineae; this considerably expanded our knowledge about the differences previously seen between these two groups of plants; (2) the compositional correlations of homologous genes from dicots and from Gramineae, as well as from both groups; all correlations were characterized by very good coefficients, with slopes close to unity in the former two cases and very high in the last; (3) the compositional transition that accompanied the emergence of Gramineae from an ancestral monocot; (4) the compositional correlations between exons and introns, which were very good in Gramineae, but only poor to good in dicots; and (5) the compositional profiles of homologous genes from angiosperms, which were characterized by a series of peaks (exons) and valleys (introns) separated by 15–20% GC. The conservative and transitional modes of compositional evolution in plant genes and their general implications are discussed. Received: 24 June 1997 / Accepted: 20 August 1997     

Compositional compartmentalization and compositional patterns in the nuclear genomes of plants.

We report here results which indicate (i) that the nuclear genomes of angiosperms is characterized by a compositional compartmentalization and an isochore structure; and (ii) that the nuclear genomes of some Gramineae exhibit strikingly different compositional patterns compared to those of many dicots. Indeed, the compositional distribution of nuclear DNA molecules (in the 50-100 Kb size range) from three dicots (pea, sunflower and tobacco) and three monocots (maize, rice and wheat) were found to be centered around lower (41%) and higher (45% for rice, 48% for maize and wheat) GC levels, respectively (and to trail towards even higher GC values in maize and wheat). Experiments on gene localization in density gradient fractions showed a remarkable compositional homogeneity in vast (greater than 100-200 Kb) regions surrounding the genes. On the other hand, the compositional distribution of coding sequences (GenBank and literature data) from dicots (several orders) was found to be narrow, symmetrical and centered around 46% GC, that from monocots (essentially barley, maize and wheat) to be broad, asymmetrical and characterized by an upward trend towards high GC values, with the majority of sequences between 60 and 70% GC. Introns exhibited a similar compositional distribution, but lower GC levels, compared to exons from the same genes.     

Gene distribution and isochore organization in the nuclear genome of plants.

The genomic distribution of 23 nuclear genes from three dicotyledons (pea, sunflower, tobacco) and five monocotyledons of the Gramineae family (barley, maize, rice, oat, wheat) was studied by localizing these genes in DNA fractions obtained by preparative centrifugation in Cs2SO4/BAMD density gradients. Each one of these genes (and of many other related genes and pseudogenes) was found to be located in DNA fragments (50-100 Kb in size) that were less than 1-2% GC apart from each other. This definitively demonstrates the existence of isochores in plant genomes, namely of compositionally homogeneous DNA regions at least 100-200 Kb in size. Moreover, the GC levels of the 23 coding sequences studied, of their first, second and third codon positions, and of the corresponding introns were found to be linearly correlated with the GC levels of the isochores harboring those genes. Compositional correlations displayed increasing slopes when going from second to first to third codon position with obvious effects on codon usage. Coding sequences for seed storage proteins and phytochrome of Gramineae deviate from the compositional correlations just described. Finally, CpG doublets of coding sequences were characterized by a shortage that decreased and vanished with increasing GC levels of the sequences. A number of these findings bear a striking similarity with results previously obtained for vertebrate genes.     

The compositional patterns of the avian genomes and their evolutionary implications

The compositional distributions of large (main-band) DNA fragments from eight birds belonging to eight different orders (including both paleognathous and neognathous species) are very broad and extremely close to each other. These findings, which are paralleled by the compositional similarity of homologous coding sequences and their codon positions, support the idea that birds are a monophyletic group.The compositional distribution of third-codon positions of genes from chicken, the only avian species for which a relatively large number of coding sequences is known, is very broad and bimodal, the minor GC-richer peak reaching 100% GC. The very high compositional heterogeneity of avian genomes is accompanied (as in the case of mammalian genomes) by a very high speciation rate compared to cold-blooded vertebrates which are characterized by genomes that are much less heterogeneous. The higher GC levels attained by avian compared to mammalian genomes might be correlated with the higher body temperature (41–43°C) of birds compared to mammals (37°C).A comparison of GC levels of coding sequences and codon positions from man and chicken revealed very close average GC levels and standard deviations. Homologous coding sequences and codon positions from man and chicken showed a surprisingly high degree of compositional similarity which was, however, higher for GC-poor than for GC-rich sequences. This indicates that GC-poor isochores of warm-blooded vertebrates reflect the composition of the isochores of the genome of the common reptilian ancestor of mammals and birds, which underwent only a small compositional change at the transition from cold- to warm-blooded vertebrates. In contrast, the GC-rich isochores of birds and mammals are the result of large compositional changes at the same evolutionary transition, where were in part different in the two classes of warm-blooded vertebrates.Correspondence to: G. Bernaadi     

Compositional properties of nuclear genes from cold-blooded vertebrates

Summary We have investigated the compositional properties of coding sequences from cold-blooded vertebrates and we have compared them with those from warm-blooded vertebrates. Moreover, we have studied the compositional correlations of coding sequences with the genomes in which they are contained, as well as the compositional correlations among the codon positions of the genes analyzed.The distribution of GC levels of the third codon positions of genes from cold-blooded vertebrates are distinctly different from those of warm-blooded vertebrates in that they do not reach the high values attained by the latter. Moreover, coding sequences from cold-blooded vertebrates are either equal, or, in most cases, lower in GC (not only in third, but also in first and second codon positions) than homologous coding sequences from warm-blooded vertebrates; higher values are exceptional. These results at the gene level are in agreement with the compositional differences between cold-blooded and warm-blooded vertebrates previously found at the whole genome (DNA) level (Bernardi and Bernardi 1990a,b).Two linear correlations were found: one between the GC levels of coding sequences (or of their third codon positions) and the GC levels of the genomes of cold-blooded vertebrates containing them; and another between the GC levels of third and first+ second codon positions of genes from cold-blooded vertebrates. The first correlation applies to the genomes (or genome compartments) of all vertebrates and the second to the genes of all living organisms. These correlations are tantamount to a genomic code.     

The compositional distribution of coding sequences and DNA molecules in humans and murids

Summary The compositional distributions of coding sequences and DNA molecules (in the 50-100-kb range) are remarkably narrower in murids (rat and mouse) compared to humans (as well as to all other mammals explored so far). In murids, both distributions begin at higher and end at lower GC values. A comparison of homologous coding sequences from murids and humans revealed that their different compositional distributions are due to differences in GC levels in all three codon positions, particularly of genes located at both ends of the distribution. In turn, these differences are responsible for differences in both codon usage and amino acids. When GC levels at first+second codon positions and third codon positions, respectively, of murid genes are plotted against corresponding GC levels of homologous human genes, linear relationships (with very high correlation coefficients and slopes of about 0.78 and 0.60, respectively) are found. This indicates a conservation of the order of GC levels in homologous genes from humans and murids. (The same comparison for mouse and rat genes indicates a conservation of GC levels of homologous genes.) A similar linear relationship was observed when plotting GC levels of corresponding DNA fractions (as obtained by density gradient centrifugation in the presence of a sequence-specific ligand) from mouse and human. These findings indicate that orderly compositional changes affecting not only coding sequences but also noncoding sequences took place since the divergence of murids. Such directional fixations of mutations point to the existence of selective pressures affecting the genome as a whole.     

Increasing expression of P450 and P450-reductase proteins from monocots in heterologous systems

Monocotyledonous crop plants are usually more resistant to herbicides than grass weeds and most dicots. Their resistance to herbicides is mediated in many cases by P450 oxygenases. Monocots thus constitute an appealing source of P450 enzymes for manipulating herbicide resistance and recombinant forms of the major xenobiotic metabolizing mooxygenases are potential tools for the optimization of new active molecules. We report here the isolation and functional characterization of the first P450 and P450 reductase coding sequences from wheat. The first attempts at expressing these cDNAs in yeast and tobacco led to levels of protein, which were extremely low, often not even detectable. The wheat P450 cDNAs were efficiently transcribed, but no protein or activity was found. Wheat coding sequences, like those of other monocots, are characterized by a high GC content and by a related strong bias of codon usage, different from that observed in yeast or dicots. Complete recoding of genes being costly, the reengineering their 5'-end using a single PCR megaprimer designed to comply with codon usage of the host was attempted. It was sufficient to relieve translation inhibition and to obtain good levels of protein expression. The same strategy also resulted in a dramatic increase in protein expression in tobacco. A basis for the success of such a partial recoding strategy, much easier and cheaper than complete recoding of the cDNA, is proposed.     

Different hydrophobicities of orthologous proteins from Xenopus and human

A compositional transition was previously detected by comparing orthologous coding sequences from cold- and warm-blooded vertebrates (see Bernardi, G., Hughes, S., Mouchiroud, D., 1997. The major compositional transitions in the vertebrate genome. J. Mol. Evol. 44, S44-S51 for a review). The transition is characterized by higher GC levels (GC is the molar ratio of guanine+cytosine in DNA) and, especially, by higher GC3 levels (GC3 is the GC level of third codon positions) in coding sequences from warm-blooded vertebrates. This transition essentially affects GC-rich genes, although the nucleotide substitution rate is of the same order of magnitude in both GC-poor and GC-rich genes. In order to understand the evolutionary basis of the changes, we have compared the hydrophobicity of orthologous proteins from Xenopus and human. Although the differences are small in proteins encoded by coding sequences ranging from 0 to 65% in GC3, they are large in the proteins encoded by sequences characterized by GC3 values higher than 65%. The latter proteins are more hydrophobic in human than in Xenopus.     

The compositional properties of human genes

Summary The present work represents the first attempt to study in greater detail previously proposed compositional correlations in genomes, based on a body of additional data relating to gene localizations as well as to extended flanking sequences extracted from gene banks. We have investigated the correlations that exist between (1) the GC levels of exons of human genes, and (2) the GC levels of either intergenic sequences or introns associated with the genes under consideration. In both cases, linear relationships with slopes close to unity were found. The similarity of the linear relationships indicates similar GC levels in intergenic sequences and introns located in the same isochores. Moreover, both intergenic sequences and introns showed GC levels 5–10% lower than the corresponding exons. The above findings considerably strengthen the previously drawn conclusion that coding and noncoding sequences (both inter- and intragenic) from the same isochores of the human genome are compositionally correlated. In addition, we find linear correlations between the GC levels of codon positions and of the intergenic sequences or introns associated with the corresponding genes, as well as among the GC levels of codon positions of genes.     

Compositional Correlations in the Chicken Genome

This paper analyses the compositional correlations that hold in the chicken genome. Significant linear correlations were found among the regions studied—coding sequences (and their first, second, and third codon positions), flanking regions (5′ and 3′), and introns—as is the case in the human genome. We found that these compositional correlations are not limited to global GC levels but even extend to individual bases. Furthermore, an analysis of 1037 coding sequences has confirmed a correlation among GC₃, GC₂, and GC₁. The implications of these results are discussed. Received: 9 December 1998 / Accepted: 18 April 1999     

The compositional transition of vertebrate genomes: an analysis of the secondary structure of the proteins encoded by human genes

Fluctuations and increments of both C(3) and G(3) levels along the human coding sequences were investigated comparing two sets of Xenopus/human orthologous genes. The first set of genes shows minor differences of the GC(3) levels, the second shows considerable increments of the GC(3) levels in the human genes. In both data sets, the fluctuations of C(3) and G(3) levels along the coding sequences correlated with the secondary structures of the encoded proteins. The human genes that underwent the compositional transition showed a different increment of the C(3) and G(3) levels within and among the structural units of the proteins. The relative synonymous codon usage (RSCU) of several amino acids were also affected during the compositional transition, showing that there exists a correlation between RSCU and protein secondary structures in human genes. The importance of natural selection for the formation of isochore organization of the human genome has been discussed on the basis of these results.     

A universal compositional correlation among codon positions.

We have investigated the compositional distributions of third codon positions of genes from the 16 prokaryotes and seven eukaryotes for which the largest numbers of coding sequences are available in data banks. In prokaryotes, both narrow and broad distributions were found. In eukaryotes, distributions were very broad (except for Saccharomyces cerevisiae) and remarkably different for different genomes. In low-GC genomes, third codon positions were lower in GC than first + second codon positions and trailed towards high GC; the opposite situation was found for high-GC genomes. In all genomes, first codon positions were higher in GC than second codon positions. We then investigated the compositional correlations between third and first + second codon positions in prokaryotic genomes (the 16 mentioned above plus 87 additional ones) and in genome compartments of eukaryotes. A general, common relationship was found, which also holds within the same (heterogeneous) genomes. This universal correlation is due to the fact that the relative effects of compositional constraints on different codon positions are the same, on the average, whatever the genome under consideration.     

Cloning and characterization of the rat cytochrome P450 4F5 (CYP4F5) gene

The analysis of a non-redundant set of human proteins, for which both the crystallographic structures and the corresponding gene sequences are available, show that bases at third codon position are non-uniformly distributed along the coding sequences. Significant compositional differences are found by comparing the gene regions corresponding to the different secondary structures of the proteins. Inter-and intra-structure differences were most pronounced in the GC-richest genes. These results are not compatible with any proposed hypotheses based on a neutral process of formation/maintenance of the high GC₃ levels of the genes localized in the GC-richest isochores of the human genome.     

Cloning, sequencing, and enhanced expression of the dihydropteroate synthase gene of Escherichia coli MC4100.

The Escherichia coli gene coding for dihydropteroate synthase (DHPS) has been cloned and sequenced. The protein has 282 amino acids and a compositional molecular mass of 30,314 daltons. Increased expression of the enzyme was realized by using a T7 expression system. The enzyme was purified and crystallized. A temperature-sensitive mutant was isolated and found to express a DHPS with a lower specific activity and lower affinities for para-aminobenzoic acid and sulfathiazole. The allele had a point mutation that changed a phenylalanine codon to a leucine codon, and the mutation was in a codon that is conserved among published DHPS sequences.     

New insights into the interplay between codon bias determinants in plants

Codon bias is the non-random use of synonymous codons, a phenomenon that has been observed in species as diverse as bacteria, plants and mammals. The preferential use of particular synonymous codons may reflect neutral mechanisms (e.g. mutational bias, G|C-biased gene conversion, genetic drift) and/or selection for mRNA stability, translational efficiency and accuracy. The extent to which these different factors influence codon usage is unknown, so we dissected the contribution of mutational bias and selection towards codon bias in genes from 15 eudicots, 4 monocots and 2 mosses. We analysed the frequency of mononucleotides, dinucleotides and trinucleotides and investigated whether the compositional genomic background could account for the observed codon usage profiles. Neutral forces such as mutational pressure and G|C-biased gene conversion appeared to underlie most of the observed codon bias, although there was also evidence for the selection of optimal translational efficiency and mRNA folding. Our data confirmed the compositional differences between monocots and dicots, with the former featuring in general a lower background compositional bias but a higher overall codon bias.     

Compositional compartmentalization and gene composition in the genome of vertebrates

Summary The compositional distribution of coding sequences from five vertebrates (Xenopus, chicken, mouse, rat, and human) is shifted toward higher GC values compared to that of the DNA molecules (in the 35–85-kb size range) isolated from the corresponding genomes. This shift is due to the lower GC levels of intergenic sequences compared to coding sequences. In the cold-blooded vertebrate, the two distributions are similar in that GC-poor genes and GC-poor DNA molecules are largely predominant. In contrast, in the warm-blooded vertebrates, GC-rich genes are largely predominant over GC-poor genes, whereas GC-poor DNA molecules are largely predominant over GC-rich DNA molecules. As a consequence, the genomes of warm-blooded vertebrates show a compositional gradient of gene concentration. The compositional distributions of coding sequences (as well as of DNA molecules) showed remarkable differences between chicken and mammals, and between mouse (or rat) and human. Differences were also detected in the compositional distribution of housekeeping and tissue-specific genes, the former being more abundant among GC-rich genes.     

Directional fixation of mutations in vertebrate evolution

Summary We have made pairwise comparisons between the coding sequences of 21 genes from coldblooded vertebrates and 41 homologous sequences from warm-blooded vertebrates. In the case of 12 genes, GC levels were higher, especially in third codon positions, in warm-blooded vertebrates compared to cold-blooded vertebrates. Six genes showed no remarkable difference in GC level and three showed a lower level. In the first case, higher GC levels appear to be due to a directional fixation of mutations, presumably under the influence of body temperature (see Bernardi and Bernardi 1986b). These GC-richer genes of warm-blooded vertebrates were located, in all cases studied, in isochores higher in GC than those comprising the homologous genes of cold-blooded vertebrates. In the third case, increases appear to be due to a limited formation of GC-rich isochores which took place in some cold-blooded vertebrates after the divergence of warm-blooded vertebrates. The directional changes in the GC content of coding sequences and the evolutionary conservation of both increased and unchanged GC levels are in keeping with the existence of compositional constraints on the genome.     

DNA序列信息的一种新的测度

根据信息理论给出了测度ＤＮＡ序列信息的一种新的方法,获得ＤＮＡ序列４个层次的信息量测度：Ｉｂ,Ｉｆ（１）,Ｉｆ（２）ａｎｄＩｆ（３）,这４种信息测度可分别用来测度ＤＮＡ的碱基序列、密码子序列、编码蛋白质序列和功能蛋白质序列的信息量。从Ｍ．ｅｄｕｌｉｓ的线粒体基因组中两个较短的编码蛋白质的ＤＮＡ序列和使用具有不同倍性的间并密码子组组成的模拟ＤＮＡ序列中所获得计算结果表明,这些信息测度确实能用来揭示所     

Characterization and phylogenetic utility of the mammalian protamine p1 gene

We sequenced the protamine P1 gene (ca. 450 bp) from 20 bats (order Chiroptera) and the flying lemur (order Dermoptera). We compared these sequences with published sequences from 19 other mammals representing seven orders (Artiodactyla, Carnivora, Cetacea, Perissodactyla, Primates, Proboscidea, and Rodentia) to assess structure, base compositional bias, and phylogenetic utility. Approximately 80% of second codon positions were guanine, resulting in protamine proteins containing a high frequency of arginine residues. Our data indicate that codon usage for arginine differs among higher mammalian taxa. Parsimony analysis of 40 species representing nine orders produced a well-resolved tree in which most nodes were supported strongly, except at the lowest taxonomic levels (e.g., within Artiodactyla and Vespertilionidae). These data support monophyly of several taxa proposed by morphologic and molecular studies (all nine orders: Laurasiatheria, Cetartiodactytla, Yangochiroptera, Noctilionoidea, Rhinolophoidea, Vespertilionoidea, Phyllostomidae, Natalidae, and Vespertilionidae) and, in agreement with recent molecular studies, reject monophyly of Archonta, Volitantia, and Microchiroptera. Bats were sister to a clade containing Perissodactyla, Carnivora, and Cetartiodactyla, and, although not unequivocally, rhinolophoid bats (traditional microchiropterans) were sister to megachiropterans. Sequences of the protamine P1 gene are useful for resolving relationships at and above the familial level in bats, and generally within and among mammalian orders, but with some drawbacks. The coding and intervening sequences are small, producing few phylogenetically informative characters, and aligning the intron is difficult, even among closely related families. Given these caveats, the protamine P1 gene may be important to future systematic studies because its functional and evolutionary constraints differ from other genes currently used in systematic studies.     

Codon usage <Emphasis Type="BoldItalic">vis-a-vis</Emphasis> start and stop codon context analysis of three dicot species

To understand the variation in genomic composition and its effect on codon usage, we performed the comparative analysis of codon usage and nucleotide usage in the genes of three dicots, Glycine max, Arabidopsis thaliana and Medicago truncatula. The dicot genes were found to be A/T rich and have predominantly A-ending and/or T-ending codons. GC3s directly mimic the usage pattern of global GC content. Relative synonymous codon usage analysis suggests that the high usage frequency of A/T over G/C mononucleotide containing codons in AT-rich dicot genome is due to compositional constraint as a factor of codon usage bias. Odds ratio analysis identified the dinucleotides TpG, TpC, GpA, CpA and CpT as over-represented, where, CpG and TpA as under-represented dinucleotides. The results of (NcExp?NcObs)/NcExp plot suggests that selection pressure other than mutation played a significant role in influencing the pattern of codon usage in these dicots. PR2 analysis revealed the significant role of selection pressure on codon usage. Analysis of varience on codon usage at start and stop site showed variation in codon selection in these sites. This study provides evidence that the dicot genes were subjected to compositional selection pressure.