A relationship between GC content and coding-sequence length

Since base composition of translational stop codons (TAG, TAA, and TGA) is biased toward a low G+C content, a differential density for these termination signals is expected in random DNA sequences of different base compositions. The expected length of reading frames (DNA segments of sense codons flanked by in-phase stop codons) in random sequences is thus a function of GC content. The analysis of DNA sequences from several genome databases stratified according to GC content reveals that the longest coding sequences—exons in vertebrates and genes in prokaryotes—are GC-rich, while the shortest ones are GC-poor. Exon lengthening in GC-rich vertebrate regions does not result, however, in longer vertebrate proteins, perhaps because of the lower number of exons in the genes located in these regions. The effects on coding-sequence lengths constitute a new evolutionary meaning for compositional variations in DNA GC content. Correspondence to: J. L. Oliver     

Rational genomics I: antisense open reading frames and codon bias in short-chain oxido reductase enzymes and the evolution of the genetic code

The short-chain oxidoreductase (SCOR) family of enzymes includes over 6000 members, extending from bacteria and archaea to humans. Nucleic acid sequence analysis reveals that significant numbers of these genes are remarkably free of stopcodons in reading frames other than the coding frame, including those on the antisense strand. The genes from this subset also use almost entirely the GC-rich half of the 64 codons. Analysis of a million hypothetical genes having random nucleotide composition shows that the percentage of SCOR genes having multiple open reading frames exceeds random by a factor of as much as 1 x 10(6). Nevertheless, screening the content of the SWISS-PROT TrEMBL database reveals that 15% of all genes contain multiple open reading frames. The SCOR genes having multiple open reading frames and a GC-rich coding bias exhibit a similar GC bias in the nucleotide triple composition of their DNA. This bias is not correlated with the GC content of the species in which the SCOR genes are found. One possible explanation for the conservation of multiple open reading frames and extreme bias in nucleic acid composition in the family of Rossman folds is that the primordial member of this family was encoded early using only very stable GC-rich DNA and that evolution proceeded with extremely limited introduction of any codons having two or more adenine or thymine nucleotides. These and other data suggest that the SCOR family of enzymes may even have diverged from a common ancestor before most of the AT-rich half of the genetic code was fully defined.     

Phenotypic expression of PCR-generated random mutations in a Pseudomonas putida gene after its introduction into an Acinetobacter chromosome by natural transformation

Localized sets of random point mutations generated by PCR amplification can be transferred efficiently to the chromosome of Acinetobacter ADP1 (also known as strain BD413) by natural transformation. The technique does not require cloning of PCR fragments in plasmids: PCR-amplified DNA fragments are internalized by cells and directly incorporated into their genomes by homologous recombination. Previously such procedures for random mutagenesis could be applied only to Acinetobacter genes affording the selection of mutant phenotypes. Here we describe the construction of a vector and recipient that allow for mutagenesis, recovery, and expression of heterologous genes that may lack a positive selection. The plasmid carries an Acinetobacter chromosomal segment interrupted by a multiple cloning site next to a kanamycin resistance marker. The insertion of heterologous DNA into the multiple cloning site prepares the insert as a target for PCR mutagenesis. PCR amplifies the kanamycin resistance marker and a flanking region of Acinetobacter DNA along with the insert of heterologous DNA. Nucleotide sequence identity between the flanking regions and corresponding chromosomal segments in an engineered Acinetobacter recipient allows homologous recombination of the PCR-amplified DNA fragments into a specific chromosomal docking site from which they can be expressed. The recipient strain contains only a portion of the kanamycin resistance gene, so donor DNA containing both this gene and the mutagenized insert can be selected by demanding growth of recombinants in the presence of kanamycin. The effectiveness of the technique was demonstrated with the relatively GC-rich Pseudomonas putida xylE gene. After only one round of PCR amplification (35 cycles), donor DNA produced transformants of which up to 30% carried a defective xylE gene after growth at 37 degrees C. Of recombinant clones that failed to express xylE at 37 degrees C, about 10% expressed the gene when grown at 22 degrees C. The techniques described here could be adapted to prepare colonies with an altered function in any gene for which either a selection or a suitable phenotypic screen exists.     

The maize gene space is compositionally compartimentalized

Previous investigations by Southern hybridization of cDNA with compositional DNA fractions showed that the majority of maize genes are located in a narrow GC range of DNA fragments and that the corresponding gene space was GC-richer than the region of the genome where zein genes are found. Here, we revisited the maize gene space using new data from the maize genome sequencing initiative. We found that the maize gene space itself is formed of two compositional compartments, i.e., a GC-poor and a GC-rich, characterized by a different distribution of Opie and Huck retrotransposons. The GC-rich compartment tends to be richer in GC-rich genes than the GC-poor compartment. However, the gene space compartimentalization of maize is much simpler than that of human.     

QuikChange shuffling: a convenient and robust method for site-directed mutagenesis and random recombination of homologous genes

Here we describe a robust method, termed QuikChange shuffling, for efficient site-directed mutagenesis and random recombination of homologous genes. The homologous genes are fragmented, and the random fragments are reassembled in a self-priming polymerase reaction to obtain chimeric genes. The product is then mixed with linearized vector and two pairs of complementary mutagenic primers, followed by assembly of the chimeric genes and linearized vector through QuikChange-like amplification to introduce recombinant plasmids with a site-directed mutation. The method, which can yield 100% chimeric genes after library construction, is more convenient and efficient than current DNA shuffling methods.     

Turtle isochore structure is intermediate between amphibians and other amniotes

Vertebrate genomes are comprised of isochores that are relatively long (>100 kb) regions with a relatively homogenous (either GC-rich or AT-rich) base composition and with rather sharp boundaries with neighboring isochores. Mammals and living archosaurs (birds and crocodilians) have heterogeneous genomes that include very GC-rich isochores. In sharp contrast, the genomes of amphibians and fishes are more homogeneous and they have a lower overall GC content. Because DNA with higher GC content is more thermostable, the elevated GC content of mammalian and archosaurian DNA has been hypothesized to be an adaptation to higher body temperatures. This hypothesis can be tested by examining structure of isochores across the reptilian clade, which includes the archosaurs, testudines (turtles), and lepidosaurs (lizards and snakes), because reptiles exhibit diverse body sizes, metabolic rates, and patterns of thermoregulation. This study focuses on a comparative analysis of a new set of expressed genes of the red-eared slider turtle and orthologs of the turtle genes in mammalian (human, mouse, dog, and opossum), archosaurian (chicken and alligator), and amphibian (western clawed frog) genomes. EST (expressed sequence tag) data from a turtle cDNA library enriched for genes that have specialized functions (developmental genes) revealed using the GC content of the third-codon-position to examine isochore structure requires careful consideration of the types of genes examined. The more highly expressed genes (e.g., housekeeping genes) are more likely to be GC-rich than are genes with specialized functions. However, the set of highly expressed turtle genes demonstrated that the turtle genome has a GC content that is intermediate between the GC-poor amphibians and the GC-rich mammals and archosaurs. There was a strong correlation between the GC content of all turtle genes and the GC content of other vertebrate genes, with the slope of the line describing this relationship also indicating that the isochore structure of turtles is intermediate between that of amphibians and other amniotes. These data are consistent with some thermal hypotheses of isochore evolution, but we believe that the credible set of models for isochore evolution still includes a variety of models. These data expand the amount of genomic data available from reptiles upon which future studies of reptilian genomics can build.     

Alteration of ultraviolet-induced mutagenesis in yeast through molecular modulation of the REV3 and REV7 gene expression

DNA damage can lead to mutations during replication. The damage-induced mutagenesis pathway is an important mechanism that fixes DNA lesions into mutations. DNA polymerase zeta (Pol zeta), formed by Rev3 and Rev7 protein complex, and Rev1 are components of the damage-induced mutagenesis pathway. Since mutagenesis is an important factor during the initiation and progression of human cancer, we postulate that this mutagenesis pathway may provide an inhibiting target for cancer prevention and therapy. In this study, we tested if UV-induced mutagenesis can be altered by molecular modulation of Rev3 enzyme levels using the yeast Saccharomyces cerevisiae as a eukaryotic model system. Reducing the REV3 expression in yeast cells through molecular techniques was employed to mimic Pol zeta inhibition. Lower levels of Pol zeta significantly decreased UV-induced mutation frequency, thus achieving inhibition of mutagenesis. In contrast, elevating the Pol zeta level by enhanced expression of both REV3 and REV7 genes led to a approximately 3-fold increase in UV-induced mutagenesis as determined by the arg4-17 mutation reversion assays. In vivo, UV lesion bypass by Pol zeta requires the Rev1 protein. Even overexpression of Pol zeta could not alleviate the defective UV mutagenesis in the rev1 mutant cells. These observations provide evidence that the mutagenesis pathway could be used as a target for inhibiting damage-induced mutations.     

Detection of nucleolar organizer and mitochondrial DNA insertion regions based on the isochore map of Arabidopsis thaliana

Eukaryotic genomes are composed of isochores, i.e. long sequences relatively homogeneous in GC content. In this paper, the isochore structure of Arabidopsis thaliana genome has been studied using a windowless technique based on the Z curve method and intuitive curves are drawn for all the five chromosomes. Using these curves, we can calculate the GC content at any resolution, even at the base level. It is observed that all the five chromosomes are composed of several GC-rich and AT-rich regions alternatively. Usually, these regions, named 'isochore-like regions', have large fluctuations in the GC content. Five isochores with little fluctuations are also observed. Detailed analyses have been performed for these isochores. A GC-rich 'isochore-like region' and a GC-isochore in chromosome II and IV, respectively, are the nucleolar organizer regions (NORs), and genes located in the two regions prefer to use GC-ending codons. Another GC-isochore located in chromosome II is a mitochondrial DNA insertion region, the position and size of this region is precisely predicted by the current method. The amino acid usage and codon preference of genes in this organellar-to-nuclear transfer region show significant difference from other regions. Moreover, the centromeres are located in GC-rich 'isochore-like regions' in all the five chromosomes. The current method can provide a useful tool for analyzing whole genomic sequences of eukaryotes.     

Theoretical analysis of mutation hotspots and their DNA sequence context specificity

Mutation frequencies vary significantly along nucleotide sequences such that mutations often concentrate at certain positions called hotspots. Mutation hotspots in DNA reflect intrinsic properties of the mutation process, such as sequence specificity, that manifests itself at the level of interaction between mutagens, DNA, and the action of the repair and replication machineries. The hotspots might also reflect structural and functional features of the respective DNA sequences. When mutations in a gene are identified using a particular experimental system, resulting hotspots could reflect the properties of the gene product and the mutant selection scheme. Analysis of the nucleotide sequence context of hotspots can provide information on the molecular mechanisms of mutagenesis. However, the determinants of mutation frequency and specificity are complex, and there are many analytical methods for their study. Here we review computational approaches for analyzing mutation spectra (distribution of mutations along the target genes) that include many mutable (detectable) positions. The following methods are reviewed: derivation of a consensus sequence, application of regression approaches to correlate nucleotide sequence features with mutation frequency, mutation hotspot prediction, analysis of oligonucleotide composition of regions containing mutations, pairwise comparison of mutation spectra, analysis of multiple spectra, and analysis of "context-free" characteristics. The advantages and pitfalls of these methods are discussed and illustrated by examples from the literature. The most reliable analyses were obtained when several methods were combined and information from theoretical analysis and experimental observations was considered simultaneously. Simple, robust approaches should be used with small samples of mutations, whereas combinations of simple and complex approaches may be required for large samples. We discuss several well-documented studies where analysis of mutation spectra has substantially contributed to the current understanding of molecular mechanisms of mutagenesis. The nucleotide sequence context of mutational hotspots is a fingerprint of interactions between DNA and DNA repair, replication, and modification enzymes, and the analysis of hotspot context provides evidence of such interactions.     

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.     

A Simple and Reproducible Method for Directed Evolution: Combination of Random Mutation with dITP and DNA Fragmentation with Endonuclease V

An alternative method to combine mutagenesis PCR with dITP and fragmentation by endonuclease V for directed evolution was developed. In comparison to the routine protocol for directed evolution, dITP was used as mutation reagent in the mutagenesis PCR. Subsequently, the incorporated dITP in the PCR products could represent as being the target of endonuclease V. Finally, the mutated dsDNA was fragmented by endonuclease V and then shuffled via assembly and reamplification as is usually done. In this study, the gene encoding kanamycin resistance has been used as reporter to verify the novel method for directed evolution. However, the mutation frequency could be easily adjusted by the amount of dITP used in the mutagenesis PCR reaction. Besides, this protocol yielded the mutation types with an obvious bias to transition substitutions as the normal error-prone PCR did. Conclusively, this novel method for directed evolution has been demonstrated to be efficient, reproducible, and easy to handle in actual practice. Using this protocol, we have successfully constructed a random mutation library for the gene encoding a serine alkaline protease.     

Evolution of chromosome bands: Molecular ecology of noncoding DNA

Summary Giemsa dark bands, G-bands, are a derived chromatin character that evolved along the chromosomes of early chordates. They are facultative heterochromatin reflecting acquisition of a late replication mechanism to repress tissue-specific genes. Subsequently, R-bands, the primitive chromatin state, became directionally GC rich as evidenced by Q-banding of mammalian and avian chromosomes. Contrary to predictions from the neutral mutation theory, noncoding DNA is positionally constrained along the banding pattern with short interspersed repeats in R-bands and long interspersed repeats in G-bands. Chromosomes seem dynamically stable: the banding pattern and gene arrangement along several human and murine autosomes has remained constant for 100 million years, whereas much of the noncoding DNA, especially retroposons, has changed. Several coding sequence attributes and probably mutation rates are determined more by where a gene lives than by what it does. R-band exons in homeotherms but not G-band exons have directionally acquired GC-rich wobble bases and the corresponding codon usage: CpG islands in mammals are specific to R-band exons, exons not facultatively heterochromatinized, and are independent of the tissue expression pattern of the gene. The dynamic organization of noncoding DNA suggests a feedback loop that could influence codon usage and stabilize the chromosome’s chromatin pattern: DNA sequences determine affinities of → proteins that together form → a chromatin that modulates → rate constants for DNA modification that determine → DNA sequences. Theories of hierarchical selection and molecular ecology show how selection can act on Darwinian units of noncoding DNA at the genome level thus creating positionally constrained DNA and contributing minimal genetic load at the individual level. Presented in part at the FEBS Symposium on Genome Organization and Evolution, held in Crete, Greece, September 1–5, 1986     

A random mutation capture assay to detect genomic point mutations in mouse tissue

Herein, a detailed protocol for a random mutation capture (RMC) assay to measure nuclear point mutation frequency in mouse tissue is described. This protocol is a simplified version of the original method developed for human tissue that is easier to perform, yet retains a high sensitivity of detection. In contrast to assays relying on phenotypic selection of reporter genes in transgenic mice, the RMC assay allows direct detection of mutations in endogenous genes in any mouse strain. Measuring mutation frequency within an intron of a transcribed gene, we show this assay to be highly reproducible. We analyzed mutation frequencies from the liver tissue of animals with a mutation within the intrinsic exonuclease domains of the two major DNA polymerases, δ and ε. These mice exhibited significantly higher mutation frequencies than did wild-type animals. A comparison with a previous analysis of these genotypes in Big Blue mice revealed the RMC assay to be more sensitive than the Big Blue assay for this application. As RMC does not require analysis of a particular gene, simultaneous analysis of mutation frequency at multiple genetic loci is feasible. This assay provides a versatile alternative to transgenic mouse models for the study of mutagenesis in vivo.     

GC-biased segregation of noncoding polymorphisms in Drosophila

The study of base composition evolution in Drosophila has been achieved mostly through the analysis of coding sequences. Third codon position GC content, however, is influenced by both neutral forces (e.g., mutation bias) and natural selection for codon usage optimization. In this article, large data sets of noncoding DNA sequence polymorphism in D. melanogaster and D. simulans were gathered from public databases to try to disentangle these two factors-noncoding sequences are not affected by selection for codon usage. Allele frequency analyses revealed an asymmetric pattern of AT vs. GC noncoding polymorphisms: AT --> GC mutations are less numerous, and tend to segregate at a higher frequency, than GC --> AT ones, especially at GC-rich loci. This is indicative of nonstationary evolution of base composition and/or of GC-biased allele transmission. Fitting population genetics models to the allele frequency spectra confirmed this result and favored the hypothesis of a biased transmission. These results, together with previous reports, suggest that GC-biased gene conversion has influenced base composition evolution in Drosophila and explain the correlation between intron and exon GC content.     

Bovine serum albumin further enhances the effects of organic solvents on increased yield of polymerase chain reaction of GC-rich templates

ABSTRACT: BACKGROUND: While being a standard powerful molecular biology technique, applications of the PCR to the amplification of high GC-rich DNA samples still present challenges which include limited yield and poor specificity of the reaction. Organic solvents, including DMSO and formamide, have been often employed as additives to increase the efficiency of amplification of high GC content (GC > 60%) DNA sequences. Bovine serum albumin (BSA) has been used as an additive in several applications, including restriction enzyme digestions as well as in PCR amplification of templates from environmental samples that contain potential inhibitors such as phenolic compounds. FINDINGS: Significant increase in PCR amplification yields of GC-rich DNA targets ranging in sizes from 0.4 kb to 7.1 kb were achieved by using BSA as a co-additive along with DMSO and formamide. Notably, enhancing effects of BSA occurs in the initial PCR cycles with BSA additions having no detrimental impact on PCR yield or specificity. When a PCR was set up such that the cycling parameters paused after every ten cycles to allow for supplementation of BSA, combining BSA and organic solvent produced significantly higher yields relative to conditions using the solvent alone. The co-enhancing effects of BSA in presence of organic solvents were also obtained in other PCR applications, including site-directed mutagenesis and overlap extension PCR. CONCLUSIONS: BSA significantly enhances PCR amplification yield when used in combination with organic solvents, DMSO or formamide. BSA enhancing effects were obtained in several PCR applications, with DNA templates of high GC content and spanning a broad size range. When added to the reaction buffer, promoting effects of BSA were seen in the first cycles of the PCR, regardless of the size of the DNA to amplify. The strategy outlined here provides a cost-effective alternative for increasing the efficiency of PCR amplification of GC-rich DNA targets over a broad size range.