Third Codon Position Research Articles

Relaxed Cleavage Specificity within the RelE Toxin Family

Bacterial type II toxin-antitoxin systems are widespread in bacteria. Among them, the RelE toxin family is one of the most abundant. The RelE(K-12) toxin of Escherichia coli K-12 represents the paradigm for this family and has been extensively studied, both in vivo and in vitro. RelE(K-12) is an endoribonuclease that cleaves mRNAs that are translated by the ribosome machinery as these transcripts enter the A site. Earlier in vivo reports showed that RelE(K-12) cleaves preferentially in the 5'-end coding region of the transcripts in a codon-independent manner. To investigate whether the molecular activity as well as the cleavage pattern are conserved within the members of this toxin family, RelE-like sequences were selected in Proteobacteria, Cyanobacteria, Actinobacteria, and Spirochaetes and tested in E. coli. Our results show that these RelE-like sequences are part of toxin-antitoxin gene pairs, and that they inhibit translation in E. coli by cleaving transcripts that are being translated. Primer extension analyses show that these toxins exhibit specific cleavage patterns in vivo, both in terms of frequency and location of cleavage sites. We did not observe codon-dependent cleavage but rather a trend to cleave upstream purines and between the second and third positions of codons, except for the actinobacterial toxin. Our results suggest that RelE-like toxins have evolved to rapidly and efficiently shut down translation in a large spectrum of bacterial species, which correlates with the observation that toxin-antitoxin systems are spreading by horizontal gene transfer.

An Extra Peptide within the Catalytic Module of a β-Agarase Affects the Agarose Degradation Pattern

Agarase hydrolyzes agarose into a series of oligosaccharides with repeating disaccharide units. The glycoside hydrolase (GH) module of agarase is known to be responsible for its catalytic activity. However, variations in the composition of the GH module and its effects on enzymatic functions have been minimally elucidated. The agaG4 gene, cloned from the genome of the agarolytic Flammeovirga strain MY04, encodes a 503-amino acid protein, AgaG4. Compared with elucidated agarases, AgaG4 contains an extra peptide (Asn(246)-Gly(302)) within its GH module. Heterologously expressed AgaG4 (recombinant AgaG4; rAgaG4) was determined to be an endo-type β-agarase. The protein degraded agarose into neoagarotetraose and neoagarohexaose at a final molar ratio of 1.5:1. Neoagarooctaose was the smallest substrate for rAgaG4, whereas neoagarotetraose was the minimal degradation product. Removing the extra fragment from the GH module led to the inability of the mutant (rAgaG4-T57) to degrade neoagarooctaose, and the final degradation products of agarose by the truncated protein were neoagarotetraose, neoagarohexaose, and neoagarooctaose at a final molar ratio of 2.7:2.8:1. The optimal temperature for agarose degradation also decreased to 40 °C for this mutant. Bioinformatic analysis suggested that tyrosine 276 within the extra fragment was a candidate active site residue for the enzymatic activity. Site-swapping experiments of Tyr(276) to 19 various other amino acids demonstrated that the characteristics of this residue were crucial for the AgaG4 degradation of agarose and the cleavage pattern of substrate.

Genome-wide survey of interindividual differences of RNA stability in human lymphoblastoid cell lines

The extent to which RNA stability differs between individuals and its contribution to the interindividual expression variation remain unknown. We conducted a genome-wide analysis of RNA stability in seven human HapMap lymphoblastoid cell lines (LCLs) and analyzed the effect of DNA sequence variation on RNA half-life differences. Twenty-six percent of the expressed genes exhibited RNA half-life differences between LCLs at a false discovery rate (FDR) < 0.05, which accounted for ~ 37% of the gene expression differences between individuals. Nonsense polymorphisms were associated with reduced RNA half-lives. In genes presenting interindividual RNA half-life differences, higher coding GC3 contents (G and C percentages at the third-codon positions) were correlated with increased RNA half-life. Consistently, G and C alleles of single nucleotide polymorphisms (SNPs) in protein coding sequences were associated with enhanced RNA stability. These results suggest widespread interindividual differences in RNA stability related to DNA sequence and composition variation.

A Review of Molecular Mechanisms in the Development of Hepatocellular Carcinoma by Aflatoxin and Hepatitis B and C Viruses

Aflatoxins are food-borne secondary fungal metabolites that are hepatotoxic, hepatocarcinogenic, and mutagenic. Urinary and serum biomarkers are more efficient in reflecting dietary exposure to aflatoxin B₁ (AFB₁) than other methods such as food sampling and dietary questionnaires. Chronic infection of the hepatitis B virus (HBV) and dietary exposure to AFB₁ are the major risk factors in a multifactorial etiology of hepatocellular carcinogenesis, raising the possibility of a synergistic interaction between 2 agents. These effects are due to the formation of DNA and protein adducts and lipid peroxidation. Most patients with hepatocellular carcinoma and HBV infection had prevalent GC → TA transversion mutation at the third position of codon 249 of the p53 gene. The HBx protein of HBV also promotes cell cycle progression, increases the expression of telomerase reverse transcriptase, inactivates negative growth regulators, and binds to and inhibits the expression of p53 (antiapoptotic activity) and other tumor suppressor genes and senescence-related factors. Some reports also evidence the role of hepatitis C virus in the pathogenesis of HCC. Inhibitors of AFB₁ adducts are found to be potent chemoprotective agents against AFB₁-induced HCC. This review focuses on the interaction of aflatoxin, HBV, and hepatitis C virus in the development of HCC.

Why has nature invented three stop codons of DNA and only one start codon?

We examine the standard genetic code with three stop codons. Assuming that the synchronization period of length 3 in DNA or RNA is violated during the transcription or translation processes, the probability of reading a frameshifted stop codon is higher than if the code would have only one stop codon. Consequently, the synthesis of RNA or proteins will soon terminate. In this way, cells do not produce undesirable proteins and essentially save energy. This hypothesis is tested on the AT-rich Drosophila genome, where the detection of frameshifted stop codons is even higher than the theoretical value. Using the binomial theorem, we establish the probability of reading a frameshifted stop codon within n steps. Since the genetic code is largely redundant, there is still space for some hidden secondary functions of this code. In particular, because stop codons do not contain cytosine, random C → U and C → T mutations in the third position of codons increase the number of hidden frameshifted stops and simultaneously the same amino acids are coded. This evolutionary advantage is demonstrated on the genomes of several simple species, e.g. Escherichia coli.

Complete genome sequence of Brachyspira intermedia reveals unique genomic features in Brachyspira species and phage-mediated horizontal gene transfer

BackgroundBrachyspira spp. colonize the intestines of some mammalian and avian species and show different degrees of enteropathogenicity. Brachyspira intermedia can cause production losses in chickens and strain PWS/AT now becomes the fourth genome to be completed in the genus Brachyspira.Results15 classes of unique and shared genes were analyzed in B. intermedia, B. murdochii, B. hyodysenteriae and B. pilosicoli. The largest number of unique genes was found in B. intermedia and B. murdochii. This indicates the presence of larger pan-genomes. In general, hypothetical protein annotations are overrepresented among the unique genes. A 3.2 kb plasmid was found in B. intermedia strain PWS/AT. The plasmid was also present in the B. murdochii strain but not in nine other Brachyspira isolates. Within the Brachyspira genomes, genes had been translocated and also frequently switched between leading and lagging strands, a process that can be followed by different AT-skews in the third positions of synonymous codons. We also found evidence that bacteriophages were being remodeled and genes incorporated into them.ConclusionsThe accessory gene pool shapes species-specific traits. It is also influenced by reductive genome evolution and horizontal gene transfer. Gene-transfer events can cross both species and genus boundaries and bacteriophages appear to play an important role in this process. A mechanism for horizontal gene transfer appears to be gene translocations leading to remodeling of bacteriophages in combination with broad tropism.

GC3 of Genes Can Be Used as a Proxy for Isochore Base Composition: A Reply to Elhaik et al.

In an article published in these pages, Elhaik et al. (Elhaik E, Landan G, Graur D. 2009. Can GC content at third-codon positions be used as a proxy for isochore composition? Mol Biol Evol. 26:1829-1833) asked if GC3, the GC level of the third-codon positions in protein-coding genes, can be used as a "proxy" to estimate the GC level of the surrounding isochore. We use available data to directly answer this simple question in the affirmative and show how the use of indirect methods can lead to apparently conflicting conclusions. The answer reasserts that in human and other vertebrates, genes have a strong tendency to reside in compositionally corresponding isochores, which has far-reaching implications for genome structure and evolution.

CpG island clusters and pro-epigenetic selection for CpGs in protein-coding exons of HOX and other transcription factors

CpG dinucleotides contribute to epigenetic mechanisms by being the only site for DNA methylation in mammalian somatic cells. They are also mutation hotspots and approximately 5-fold depleted genome-wide. We report here a study focused on CpG sites in the coding regions of Hox and other transcription factor genes, comparing methylated genomes of Homo sapiens, Mus musculus, and Danio rerio with nonmethylated genomes of Drosophila melanogaster and Caenorhabditis elegans. We analyzed 4-fold degenerate, synonymous codons with the potential for CpG. That is, we studied "silent" changes that do not affect protein products but could damage epigenetic marking. We find that DNA-binding transcription factors and other developmentally relevant genes show, only in methylated genomes, a bimodal distribution of CpG usage. Several genetic code-based tests indicate, again for methylated genomes only, that the frequency of silent CpGs in Hox genes is much greater than expectation. Also informative are NCG-GNN and NCC-GNN codon doublets, for which an unusually high rate of G to C and C to G transversions was observed at the third (silent) position of the first codon. Together these results are interpreted as evidence for strong "pro-epigenetic" selection acting to preserve CpG sites in coding regions of many genes controlling development. We also report that DNA-binding transcription factors and developmentally important genes are dramatically overrepresented in or near clusters of three or more CpG islands, suggesting a possible relationship between evolutionary preservation of CpG dinucleotides in both coding regions and CpG islands.

GENOME ANALYSIS OF GENBANK KNOWN RABBIT (Oryctolagus cuniculus) GENES

After downloading all known rabbit genes, key words specific to complete coding sequences were used to filter out partial coding sequences. As a result, 160 full-length, nuclear, protein-coding and functionally annotated genes were extracted from GenBank database. These genes were subjected to synonymous codon and amino acid usage analysis. The results showed a clear base composition bias in the genes analyzed. The effective number of codons used (Nc values) ranged between 30.07 and 59.98 with a mean of 51.31 with a standard deviation (SD) of 5.71. The frequency of G + C at the synonymous third position of codons (GC3s) varied between 0.3 and 0.96 with a mean of 0.55 and a SD of 0.14, clearly indicating marked variation and heterogeneity in codon usage patterns among the different genes. The distribution of relative synonymous codon usage (RSCU) values calculated for all the genes indicated that codons with a G or C in the third position are widely employed, and although RCSU values for A/T or G/C ending codons are almost equal, G/C appears to play a dominant role. This pattern of dominance was confirmed by the distribution of amino acid occurrence. Leucine (CUN) and serine (UCN) were the most frequent amino acids, followed by arginine (CGN), proline (CCN), glycine (GGN) and alanine (GCN), in relatively equal proportions. The heterogeneity observed in the analyzed genes was then further probed by multivariate statistical analyses. A major trend in codon usage that correlated with gene expression values was revealed. These findings suggest that translational efficiency exerts a stronger influence on codon usage preference than compositional bias in the sampled rabbit genes.

Patterns of nucleotides that flank substitutions in human orthologous genes

BackgroundSequence context is an important aspect of base mutagenesis, and three-base periodicity is an intrinsic property of coding sequences. However, how three-base periodicity is influenced in the vicinity of substitutions is still unclear. The effect of context on mutagenesis should be revealed in the usage of nucleotides that flank substitutions. Relative entropy (also known as Kullback-Leibler divergence) is useful for finding unusual patterns in biological sequences.ResultsUsing relative entropy, we visualized the periodic patterns in the context of substitutions in human orthologous genes. Neighbouring patterns differed both among substitution categories and within a category that occurred at three codon positions. Transition tended to occur in periodic sequences relative to transversion. Periodic signals were stronger in a set of flanking sequences of substitutions that occurred at the third-codon positions than in those that occurred at the first- or second-codon positions. To determine how the three-base periodicity was affected near the substitution sites, we fitted a sine model to the values of the relative entropy. A sine of period equal to 3 is a good approximation for the three-base periodicity at sites not in close vicinity to some substitutions. These periods were interrupted near the substitution site and then reappeared away from substitutions. A comparative analysis between the native and codon-shuffled datasets suggested that the codon usage frequency was not the sole origin of the three-base periodicity, implying that the native order of codons also played an important role in this periodicity. Synonymous codon shuffling revealed that synonymous codon usage bias was one of the factors responsible for the observed three-base periodicity.ConclusionsOur results offer an efficient way to illustrate unusual periodic patterns in the context of substitutions and provide further insight into the origin of three-base periodicity. This periodicity is a result of the native codon order in the reading frame. The length of the period equal to 3 is caused by the usage bias of nucleotides in synonymous codons. The periodic features in nucleotides surrounding substitutions aid in further understanding genetic variation and nucleotide mutagenesis.

Construction and characterization of a normalized cDNA library of Nannochloropsis oculata (Eustigmatophyceae)

We constructed and characterized a normalized cDNA library of Nannochloropsis oculata CS-179, and obtained 905 nonredundant sequences (NRSs) ranging from 431–1 756 bp in length. Among them, 496 were very similar to nonredundant ones in the GenBank (E ≤1.0e-05), and 349 ESTs had significant hits with the clusters of eukaryotic orthologous groups (KOG). Bases G and/or C at the third position of codons of 14 amino acid residues suggested a strong bias in the conserved domain of 362 NRSs (>60%). We also identified the unigenes encoding phosphorus and nitrogen transporters, suggesting that N. oculata could efficiently transport and metabolize phosphorus and nitrogen, and recognized the unigenes that involved in biosynthesis and storage of both fatty acids and polyunsaturated fatty acids (PUFAs), which will facilitate the demonstration of eicosapentaenoic acid (EPA) biosynthesis pathway of N. oculata. In comparison with the original cDNA library, the normalized library significantly increased the efficiencies of random sequencing and rarely expressed genes discovering, and decreased the frequency of abundant gene sequences.

Contrasting GC-content dynamics across 33 mammalian genomes: Relationship with life-history traits and chromosome sizes

The origin, evolution, and functional relevance of genomic variations in GC content are a long-debated topic, especially in mammals. Most of the existing literature, however, has focused on a small number of model species and/or limited sequence data sets. We analyzed more than 1000 orthologous genes in 33 fully sequenced mammalian genomes, reconstructed their ancestral isochore organization in the maximum likelihood framework, and explored the evolution of third-codon position GC content in representatives of 16 orders and 27 families. We showed that the previously reported erosion of GC-rich isochores is not a general trend. Several species (e.g., shrew, microbat, tenrec, rabbit) have independently undergone a marked increase in GC content, with a widening gap between the GC-poorest and GC-richest classes of genes. The intensively studied apes and (especially) murids do not reflect the general placental pattern. We correlated GC-content evolution with species life-history traits and cytology. Significant effects of body mass and genome size were detected, with each being consistent with the GC-biased gene conversion model.

Sequence analysis for the 3' end of genomes from enterovirus 71 isolated in Beijing during 2008-2009

Objective To investigate the sequences of the 3' end of genomes from enterovirus 71 isolated from pediatric patients with different symptoms in Beijing during 2008-2009.Methods Clinical specimens were collected from pediatric patients suffering from hand,foot and mouth disease(HFMD)and/or with neurological complications who visited the affiliated Children's Hospital during the epidemic seasons of 2008 and 2009 in Beijing.The samples were inoculated into the Vero cell line,and the virus isolates were further identified by nested reverse transcription-nested PCR(RT-nPCR)assay using universal enterovirus primers,type specific primers for EV71 and CA16.The 3' end of genomes(including 3D and 3' UTR regions)from 10 EV71 strains derived from various clinical presentations were amplified and sequenced.Results Analysis of the nucleotide sequences of the amplified fragments showed that the 3' end of genomes of 10 EV71 isolates include the 3D region of 1386 nucleotides(nt)encoding the 3D polymerase(3Dpol,462 amino acids),the terminator codon TGA and 3' UTR of 81 nt.Nucleotide sequence identities among 3D regions from these 10 EV71 isolates were in the range of 95.8%-99.6%,while the nucleotide sequence identities for 3' UTR were 96.3%-100%.The majority of nucleotide changes were located at the third codon positions which caused silent mutations,thus the amino acid sequence changes of 3Dpol among those 10 EV71isolates were scanty.The residues 140 and 263 which were R and 1 were substituted by K and V,respectively in 3 of 4 neurovirulent strains,whereas only 1 of 6 strains from mild cases had these 2 amino acid changes.The sequences of the 3D and 3' UTR regions of 10 EV71 isolates were compared to the representative strains of known genotypes from the GenBank.The nucleotide and amino acid sequences of 10 EV71 isolates in the 3D region exhibited highest homology to the subgenotype C4 of EV71(92.7%-94.2%and 96.8%-97.6%.respectively).However,3' UTR of 10 EV71 isolates shared the highest nucleotide identity with CA16/G10(88.9%-91.4%).The phylogenetic analysis based on the 3D regions demonstrated that 10 EV71 isolates had the closest genetic relationship with the representative isolate of CA sub-genotype of EV71 and shared more homology with CA16/G10 than other known genotypes of EV71.Conclusion Genetic analysis of the 3' end of genomes from 10 EV71 strains indicated that the 3' end of genome may play a role in the evolution of EV71. Key words: Hand,foot and mouth disease; Enterovirus 71; 3' End; Sequence analysis

Factors affecting synonymous codon usage bias in the gC gene of DPV CHv strain

The aim was to identify the factors affecting synonymous codon usage bias in the newly discovered gC gene of duck plague virus (DPV) CHv strain and a comparative analysis of the codon usage bias in the gC genes of 27 other reference herpesviruses was performed by using CAI, CHIPS and CUSP program of EMBOSS. The results showed that base composition, mutational bias and natural selection were the major determinants of the codon usage variation in the DPV gC gene. The primary codon usage trends in the DPV gC gene showed obvious difference with a strong bias towards the synonymous codons with A and T at the third codon position. Characterization of nucleotide composition in the DPV gC gene was related to dinucleotide usage bias characterized by the variation of CpG islands. The ENc-plot and GC12s-GC3s revealed that the genetic heterogeneity in the gC genes of the 28 herpesviruses were constrained by G + C content. Moreover, the phylogenetic analysis demonstrated that codon usage patterns of the DPV gC gene were phylogenetically conserved and similar to the gC genes of the avian alphaherpesvirinae. Furthermore, comparisons of the codon preferences in the DPV gC gene with those of Escherichia coli, yeast and human revealed that there was a statistically positive correlation between DPV and yeast (r = 0.646, P < 0.05). The above results could provide useful information for the synonymous codons usage bias of the DPV gC gene and promote the relevant mechanism for evolution, pathogenesis and functional studies in the area of DPV research and possibly studies with other herpesvirus viruses. Key words: Duck plague virus (DPV), gC gene, codon usage bias.

Codon Usage Patterns inCorynebacterium glutamicum: Mutational Bias, Natural Selection and Amino Acid Conservation

The alternative synonymous codons in Corynebacterium glutamicum, a well-known bacterium used in industry for the production of amino acid, have been investigated by multivariate analysis. As C. glutamicum is a GC-rich organism, G and C are expected to predominate at the third position of codons. Indeed, overall codon usage analyses have indicated that C and/or G ending codons are predominant in this organism. Through multivariate statistical analysis, apart from mutational selection, we identified three other trends of codon usage variation among the genes. Firstly, the majority of highly expressed genes are scattered towards the positive end of the first axis, whereas the majority of lowly expressed genes are clustered towards the other end of the first axis. Furthermore, the distinct difference in the two sets of genes was that the C ending codons are predominate in putatively highly expressed genes, suggesting that the C ending codons are translationally optimal in this organism. Secondly, the majority of the putatively highly expressed genes have a tendency to locate on the leading strand, which indicates that replicational and transciptional selection might be invoked. Thirdly, highly expressed genes are more conserved than lowly expressed genes by synonymous and nonsynonymous substitutions among orthologous genes fromthe genomes of C. glutamicum and C. diphtheriae. We also analyzed other factors such as the length of genes and hydrophobicity that might influence codon usage and found their contributions to be weak.

Allele-specific SNP markers for the new low linolenic mutant genotype of winter oilseed rape

With 5 figures and 1 table Abstract One of the goals in oilseed rape breeding is the identification of genotypes with low linolenic acid content in seed oil. Here, we present new genetic markers for mutant alleles of a chemically induced low linolenic rapeseed line DH219/05. Genomic clones comprising fatty acids desaturase 3 (FAD3) genes from the mutant and wild-type rapeseed lines were sequenced. Two statistically important single nucleotide polymorphisms were detected: (1) a C to T substitution in the third position of the sixth codon of the seventh exon in the BnaA.FAD3 gene and (2) a G to A transition in the 5′ splice donor site of the sixth intron in the BnaC.FAD3 gene. Allele-specific SNP markers were designed involving detection of the wild-type and mutant alleles by SNaPshot analysis and locus-specific PCR primers. Strong negative correlation between the presence of mutant alleles in the A and C genomes and linolenic acid content was revealed by analysis of variance. Sequence analysis of transcript variants confirmed predictions on possible negative effects of mutations on FAD3 gene expression.

Clonal diversity in Giardia duodenalis isolates from Thailand: Evidences for intragenic recombination and purifying selection at the beta giardin locus

The beta giardin locus of Giardia duodenalis encodes a structural component of ventral disc and exhibits sequence variation among isolates rendering it a useful marker for genotypic analysis. To determine the distribution of genotypes of G. duodenalis in Thailand and to explore the extent of sequence variation in this locus, we deployed the PCR-RFLP method and sequence analysis of recombinant subclones from 30 clinical isolates. In total, assemblage B was more prevalent than assemblage A. Sequence analysis revealed that 13 isolates had clonal mixtures, comprising three to five distinct sequences per isolate. Nucleotide diversity of assemblage B was greater than that of assemblage A. A striking transitional bias was noted at the first and the third positions of codons in both assemblages; however, they differed in the patterns of nucleotide diversity at 0-fold and 4-fold-degenerate sites. Most amino acid exchanges were conservative in terms of polarity, charge and volume. Both assemblages have evolved under purifying selection as evidenced by a significantly greater number of mean synonymous substitutions per synonymous site ( d S) than that of nonsynonymous substitutions per nonsynonymous site ( d N) as well as significant negative Tajima's D values and its related statistics. The significant negative Tajima's D test at nonsynonymous sites further suggests that elimination of slightly deleterious mutations at these sites by purifying selection is ongoing as predicted in the nearly neutral theory. Furthermore, a minimum number of seven recombination sites was detected by the four gamete test in assemblage B, consistent with previous reports on meiotic recombination in G. duodenalis. Therefore, accurate subassemblage assignment of clinical isolates that has practical consequence for disease control could be complicated by the presence of intra-isolate clonal diversity and interallelic recombination.

Preparation of an ochre suppressor tRNA recognizing exclusively UAA codon by using the molecular surgery technique

In order to create an ochre suppressor tRNA which exclusively recognizes UAA codon, we replaced the G34 at the first position of yeast tRNA(Tyr)[GPsiA] anticodon with pseudouridine34 (Psi34) by using the molecular surgery technique. This tRNA(Tyr)[PsiPsiA] recognized only the UAA codon as expectedly, but tRNA(Tyr)[UPsiA] made as a control also behaved similarly. This result may suggest that U34 must be somehow modified to facilitate the wobble-pairing to G at the third position of codon.

Loss of matKRNA editing in seed plant chloroplasts

BackgroundRNA editing in chloroplasts of angiosperms proceeds by C-to-U conversions at specific sites. Nuclear-encoded factors are required for the recognition of cis-elements located immediately upstream of editing sites. The ensemble of editing sites in a chloroplast genome differs widely between species, and editing sites are thought to evolve rapidly. However, large-scale analyses of the evolution of individual editing sites have not yet been undertaken.ResultsHere, we analyzed the evolution of two chloroplast editing sites, matK-2 and matK-3, for which DNA sequences from thousands of angiosperm species are available. Both sites are found in most major taxa, including deep-branching families such as the nymphaeaceae. However, 36 isolated taxa scattered across the entire tree lack a C at one of the two matK editing sites. Tests of several exemplary species from this in silico analysis of matK processing unexpectedly revealed that one of the two sites remain unedited in almost half of all species examined. A comparison of sequences between editors and non-editors showed that specific nucleotides co-evolve with the C at the matK editing sites, suggesting that these nucleotides are critical for editing-site recognition.Conclusion(i) Both matK editing sites were present in the common ancestor of all angiosperms and have been independently lost multiple times during angiosperm evolution.(ii) The editing activities corresponding to matK-2 and matK-3 are unstable.(iii) A small number of third-codon positions in the vicinity of editing sites are selectively constrained independent of the presence of the editing site, most likely because of interacting RNA-binding proteins.

Can GC Content at Third-Codon Positions Be Used as a Proxy for Isochore Composition?

The isochore theory depicts the genomes of warm-blooded vertebrates as a mosaic of long genomic regions that are characterized by relatively homogeneous GC content. In the absence of genomic data, the GC content at third-codon positions of protein-coding genes (GC3) was commonly used as a proxy for the GC content of isochores. Oddly, in the postgenomic era, GC3 is still sometimes used as a proxy for the GC composition of isochores. Here, we use genic and genomic sequences from human, chimpanzee, cow, mouse, rat, chicken, and zebrafish to show that GC3 only explains a very small proportion of the variation in GC content of long genomic sequences flanking the genes (GCf), and what little correlation there is between GC3 and GCf was found to decay rapidly with distance from the gene. The coefficient of variation of GC3 was found to be much larger than that of GCf and, therefore, GC3 and GCf values are not comparable with each other. Comparisons of orthologous gene pairs from 1) human and chimpanzee and 2) mouse and rat show strong correlations between their GC3 values, but very weak correlations between their GCf values. We conclude that the GC content of third-codon position cannot be used as stand-in for isochoric composition.