Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Codon optimality”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 505 records · Page 28Linked to original sources
PubMed1998

A decision tree system for finding genes in DNA.

MORGAN is an integrated system for finding genes in vertebrate DNA sequences. MORGAN uses a variety of techniques to accomplish this task, the most distinctive of which is a decision tree classifier. The decision tree system is combined with new methods for identifying start codons, donor sites, and acceptor sites, and these are brought together in a frame-sensitive dynamic programming algorithm that finds the optimal segmentation of a DNA sequence into coding and noncoding regions (exons and introns). The optimal segmentation is dependent on a separate scoring function that takes a subsequence and assigns to it a score reflecting the probability that the sequence is an exon. The scoring functions in MORGAN are sets of decision trees that are combined to give a probability estimate. Experimental results on a database of 570 vertebrate DNA sequences show that MORGAN has excellent performance by many different measures. On a separate test set, it achieves an overall accuracy of 95 %, with a correlation coefficient of 0.78, and a sensitivity and specificity for coding bases of 83 % and 79%. In addition, MORGAN identifies 58% of coding exons exactly; i.e., both the beginning and end of the coding regions are predicted correctly. This paper describes the MORGAN system, including its decision tree routines and the algorithms for site recognition, and its performance on a benchmark database of vertebrate DNA.

Algorithms↗
PubMedDate not supplied

Conversion of nucleotides sequences into genomic signals.

An original tetrahedral representation of the Genetic Code (GC) that better describes its structure, degeneration and evolution trends is defined. The possibility to reduce the dimension of the representation by projecting the GC tetrahedron on an adequately oriented plane is also analyzed, leading to some equivalent complex representations of the GC. On these bases, optimal symbolic-to-digital mappings of the linear, nucleic acid strands into real or complex genomic signals are derived at nucleotide, codon and amino acid levels. By converting the sequences of nucleotides and polypeptides into digital genomic signals, this approach offers the possibility to use a large variety of signal processing methods for their handling and analysis. It is also shown that some essential features of the nucleotide sequences can be better extracted using this representation. Specifically, the paper reports for the first time the existence of a global helicoidal wrapping of the complex representations of the bases along DNA sequences, a large scale trend of genomic signals. New tools for genomic signal analysis, including the use of phase, aggregated phase, unwrapped phase, sequence path, stem representation of components' relative frequencies, as well as analysis of the transitions are introduced at the nucleotide, codon and amino acid levels, and in a multiresolution approach.

Amino Acid Sequence↗
PubMed1988

Mutational analysis of the HIS4 translational initiator region in Saccharomyces cerevisiae.

We have mutated various features of the 5' noncoding region of the HIS4 mRNA in light of established Saccharomyces cerevisiae and mammalian consensus translational initiator regions. Our analysis indicates that insertion mutations that introduce G + C-rich sequences in the leader, particularly those that result in stable stem-loop structures in the 5' noncoding region of the HIS4 message, severely affect translation initiation. Mutations that alter the length of the HIS4 leader from 115 to 39 nucleotides had no effect on expression, and sequence context changes both 5' and 3' to the HIS4 AUG start codon resulted in no more than a twofold decrease of expression. Changing the normal context at HIS4 5'-AAUAAUGG-3' to the optimal sequence context proposed for mammalian initiator regions 5'-CACCAUGG-3' did not result in stimulation of HIS4 expression. These studies, in conjunction with comparative and genetic studies in S. cerevisiae, support a general mechanism of initiation of protein synthesis as proposed by the ribosomal scanning model.

Alleles↗
PubMed2003

Development of a dual-luciferase fusion gene as a sensitive marker for site-directed DNA repair strategies.

BACKGROUND: Several novel techniques have been developed recently for the site-specific repair of DNA as an approach to gene therapy. Correction efficiencies as high as 40% have been reported, well within the range of therapeutic impact for a number of genetic diseases. Unfortunately, many of the model systems in which these methods have been employed typically target genes that are not well suited for analyzing the various techniques. METHODS: To address this, we have constructed and characterized a dual-luciferase fusion gene as a sensitive marker for optimizing repair strategies. The genes encoding two distinct luciferase proteins were fused so that expression of one luciferase necessitated expression of the other. Engineering a stop codon in the downstream luciferase gene created an ideal tool to study the efficiency of various site-directed DNA repair techniques as one luciferase can act as an internal control while the other is targeted for correction. RESULTS: Fusing two luciferase genes resulted in a single protein that produces two bioluminescent activities in a constant ratio. The utility of this system as a target for site-directed DNA repair research was demonstrated using two of the recently developed gene repair techniques, small fragment homologous replacement and oligonucleotide-mediated repair, to mediate correction and by the ability to detect repair efficiencies of less than 5 x 10(-6) (<1 event in 200000). CONCLUSIONS: The ability to rapidly and accurately quantify the amount of correction using the dual-luciferase fusion system will allow the comparison and evaluation of the many factors involved in successful gene repair and lead to the optimization of these techniques, both in cell culture and in whole animals.

Amino Acid Sequence↗
PubMed1997

Two genes encoding an endoglucanase and a cellulose-binding protein are clustered and co-regulated by a TTA codon in Streptomyces halstedii JM8.

Streptomyces halstedii JM8 Cel2 is an endoglucanase of 28 kDa that is first produced as a protein of 42 kDa (p42) and is later processed at its C-terminus. Cel2 displays optimal activity towards CM-cellulose at pH6 and 50 degrees C and shows no activity against crystalline cellulose or xylan. The N-terminus of p42 shares similarity with cellulases included in family 12 of the beta-glycanases and the C-terminus shares similarity with bacterial cellulose-binding domains included in family II. This latter domain enables the precursor to bind so tightly to Avicel that it can only be eluted by boiling in 10% (w/v) SDS. Another open reading frame (ORF) situated 216 bp downstream from the p42 ORF encodes a protein of 40 kDa (p40) that does not have any clear hydrolytic activity against cellulosic or xylanosic compounds, but shows high affinity for Avicel (crystalline cellulose). The p40 protein is processed in old cultures to give a protein of 35 kDa that does not bind to Avicel. Translation of both ORFs is impaired in Streptomyces coelicolor bldA mutants, suggesting that a TTA codon situated at the fourth position of the first ORF is responsible for this regulation. S1 nuclease protection experiments demonstrate that both ORFs are co-transcribed.

Amino Acid Sequence↗
PubMed1987

Endosymbiotic origin and codon bias of the nuclear gene for chloroplast glyceraldehyde-3-phosphate dehydrogenase from maize.

The nuclei of plant cells harbor genes for two types of glyceraldehyde-3-phosphate dehydrogenases (GAPDH) displaying a sequence divergence corresponding to the prokaryote/eukaryote separation. This strongly supports the endosymbiotic theory of chloroplast evolution and in particular the gene transfer hypothesis suggesting that the gene for the chloroplast enzyme, initially located in the genome of the endosymbiotic chloroplast progenitor, was transferred during the course of evolution into the nuclear genome of the endosymbiotic host. Codon usage in the gene for chloroplast GAPDH of maize is radically different from that employed by present-day chloroplasts and from that of the cytosolic (glycolytic) enzyme from the same cell. This reveals the presence of subcellular selective pressures which appear to be involved in the optimization of gene expression in the economically important graminaceous monocots.

Amino Acid Sequence↗
PubMed2003

Block structure and stability of the genetic code.

It is known that different codons may be unified into larger groups related to the hierarchical structure, approximate hidden symmetries, and evolutionary origin of the universal genetic code. Using a simplified evolutionary motivated two-letter version of genetic code, the general principles of the most stable coding are discussed. By the complete enumeration in such a reduced code it is strictly proved that the maximum stability with respect to point mutations and shifts in the reading frame needs the fixation of the middle letters within codons in groups with different physico-chemical properties, thus, explaining a key feature of the universal genetic code. The translational stability of the genetic code is studied by the mapping of code onto de Bruijn graph providing both the compact visual representation of mutual relationships between different codons as well as between codons and protein coding DNA sequence and a powerful tool for the investigation of stability of protein coding. Then, the results are extended to four-letter codes. As is shown, the universal genetic code obeys mainly the principles of optimal coding. These results demonstrate the hierarchical character of optimization of universal genetic code with strictly optimal coding being evolved at the earliest stages of molecular evolution. Finally, the universal genetic code is compared with the other natural variants of genetic codes.

Amino Acids↗
PubMed1984

Genetic code and optimal resistance to the effects of mutations.

This paper deals with the notion of resistance of the genetic code to the effects of mutations. We measure the resistance of a group of t codons as the number of pairs of those which differ from each other in only one of their three bases. We find for each value of t the maximum possible value of the resistance and we describe some groups of codons giving this value. Important examples of such configurations are found in the genetic code, among these are the groups of synonymous codons, as observed elsewhere, and the cluster of codons which have an hydrophobic amino acid for translation.

Base Sequence↗
PubMed1994

Optimally parsing a sequence into different classes based on multiple types of evidence.

We consider the problem of parsing a sequence into different classes of subsequences. Two common examples are finding the exons and introns in genomic sequences and identifying the secondary structure domains of protein sequences. In each case there are various types of evidence that are relevant to the classification, but none are completely reliable, so we expect some weighted average of all the evidence to provide improved classifications. For example, in the problem of identifying coding regions in genomic DNA, the combined use of evidence such as codon bias and splice junction patterns can give more reliable predictions than either type of evidence alone. We show three main results: 1. For a given weighting of the evidence a dynamic programming algorithm returns the optimal parse and any number of sub-optimal parses. 2. For a given weighting of the evidence a dynamic programming algorithm determines the probability of the optimal parse and any number of sub-optimal parses under a natural Boltzmann-Gibbs distribution over the set of possible parses. 3. Given a set of sequences with known correct parses, a dynamic programming algorithm allows one to apply gradient descent to obtain the weights that maximize the probability of the correct parses of these sequences.

Algorithms↗
PubMed1995

Many combinations of amino acid sequences in a conserved region of the D1 protein satisfy photosystem II function.

The putative de helix of the D1 protein is located at the acceptor side of photosystem II (PS II) and serves as an indispensable part of a niche that binds the secondary plastoquinone QB. Combinatorial mutagenesis was applied to a stretch of four residues in a highly conserved region of this putative helix in order to reveal amino acid combinations that are able to support PS II function. An obligate photoheterotrophic mutant of the cyanobacterium Synechocystis sp. PCC 6803, missing four residues (delta YFGR254-7) in the de helix, was transformed with a D1-coding sequence carrying fully degenerate combinations of codons at the site of the deletion. Upon selection for photoautotrophy, 25 mutants with functional PS II were isolated. All mutants showed different codon combinations at positions 254 to 257; none was identical to the wild-type sequence, and none of the conserved residues was found to be mandatory for PS II function. However, 24 of the mutants contained Tyr of Phe at position 254 while at the other three positions many different amino acid combinations could be functionally accommodated. Most sequences maintained an amphiphilic arrangement of the helix that may align Tyr254 facing the QB binding pocket. This residue is proposed to be functionally analogous to Phe216 of the L subunit in purple bacteria which contributes to binding of QB. Most of the PS II properties were similar in the mutants compared to wild-type. Noticeable modifications in the mutants concerned the semiquinone equilibrium of electron transfer between QA and QB, and the affinity of PS II inhibitors. Differential effects on the semiquinone equilibrium were observed between two distinct quinones occupying the QB site (plastoquinone versus 2,5-dichloro-p-benzo- quinone), implying that residues in this domain are involved, directly or indirectly, with different binding determinants of the quinones. Even though many different combinations of amino acids in positions 254 to 257 of the D1 protein may satisfy the primary function of PS II, complex requirements need to be combined for optimized performance of the QB binding niche.

Amino Acid Sequence↗
PubMed1999

Local activation and inactivation of thyroid hormones: the deiodinase family.

Tissue-specific activation and inactivation of ligands of nuclear receptors which belong to the steroid retinoid-thyroid hormone superfamily of transcription factors represents an important principle of development- and tissue-specific local modulation of hormone action. Recently, several enzyme families have been identified which act as 'guardians of the gate' of ligand-activated transcription modulation. Three monodeiodinase isoenzymes which are involved in activation the 'prohormone' L-thyroxine (T4), the main secretory product of the thyroid gland, have been identified, characterized, and cloned. Both, type I and type II 5'-deiodinase generate the thyromimetically active hormone 3,3',5-triiodothyronine (T3) by reductive deiodination of the phenolic ring of T4. Inactivation of T4 and its product T3 occurs by deiodination of iodothyronines at the tyrosyl ring. This reaction is catalyzed both the type III 5-deiodinase and also by the type I enzyme, which has a broader substrate specificity. The three deiodinases appear to constitute a newly discovered family of selenocysteine-containing proteins and the presence of selenocysteine in the protein is critical for enzyme activity. Whereas the selenoenzyme characteristics of the type I and type III deiodinases are definitively established some controversy still exists for the type II 5'-deiodinase in mammals. The mRNA probably encoding the type II 5'-deiodinase subunit is markedly longer than those of the two other deiodinases and its selenocysteine-insertion element is located more than 5 kB downstream of the UGA-codon in the 3'-untranslated region. The three deiodinase isoenzymes show a distinct development- and tissue-specific pattern of expression, operate at individual optimal substrate levels, are differently regulated and modulated by hormones, cytokines, signaling pathways, natural factors, and pharmaceuticals. Whereas circulating T3 mainly originates from hepatic production via the type I 5'-deiodinase, the local cellular thyroid hormone concentration in various tissues including the central nervous system is controlled by complex para-, auto-, and intracrine interactions of all three deiodinases. Local thyroid hormone availability is further modulated by conjugation reactions of the phenolic 4'-OH-group of iodothyronines, which also inactivate the thyroid hormones.

Animals↗
PubMedDate not supplied

Using proteomics to mine genome sequences.

We present a method for mining unannotated or annotated genome sequences with proteomic data to identify open reading frames. The region of a genome coding for a protein sequence is identified by using information from the analysis of proteins and peptides with MALDI-TOF mass spectrometry. The raw genome sequence or any unassembled contigs of an organism are theoretically cleaved into a number of equal sized but overlapping fragments, and these are then translated in all six frames into a series of virtual proteins. Each virtual protein is then subjected to a theoretical enzymatic digestion. Standard proteomic sample preparation methods are used to separate, array, and digest the proteins of interest to peptides. The masses of the resulting peptides are measured using mass spectrometry and compared to the theoretical peptide masses of the virtual proteins. The region of the genome responsible for coding for a particular protein can then be identified when there are a large number of hits between peptides from the protein and peptides from the virtual protein. The method makes no assumptions about the location of a protein in a particular gene sequence or the positions or types of start and stop codons. To illustrate this approach, all 773 proteins of Pseudomonas aeruginosa contained in SWISS-PROT were used to theoretically test the method and optimize parameters. Increasing the size of the virtual proteins results in an overall improvement in the ability to detect the coding region, at the cost of decreasing the sensitivity of the method for smaller proteins. Increasing the minimum number of matching peptides, lowering the mass error tolerance, or increasing the signal-to-noise ratio of the simulated mass spectrum, improves the ability to detect coding regions. The method is further demonstrated on experimental data from Mycobacterium tuberculosis and is also shown to work with eukaryotic organisms (e.g., Homo sapiens).

Amino Acid Sequence↗
PubMed2003

Expression of the Bacillus thuringiensis Cyt2Aa1 toxin in Pichia pastoris using a synthetic gene construct.

The nucleotide sequence data corresponding to the syncyt2Aa1 open reading frame was deposited with the EMBL/GenBank Nucleotide Sequence Databases under the accession number AF398463. Bacillus thuringiensis delta-endotoxins are membrane-active, pore-forming proteins with highly specific insecticidal activities. In addition to a well-established role in the biological control of a wide variety of crop pests and disease vectors, these toxins also have great potential for the development of anti-tumour agents called immunotoxins (ITs), chimaeric molecules consisting of a cell-binding ligand coupled to a toxin or its subunits. The ultimate goal of our study was the recombinant production of such ITs based on the Cyt2Aa1 toxin from B. thuringiensis subspecies kyushuensis. We explored the use of Pichia pastoris for recombinant IT production because earlier attempts in our laboratory using the Escherichia coli expression system or various chemical conjugation strategies yielded only low levels of functional product. However, our initial attempts were not successful because the A+T-rich bacterial cyt2Aa1 gene contained fortuitous polyadenylation sites, causing premature transcription termination in this yeast. Accordingly, we designed and constructed a synthetic cyt2Aa1 gene (syncyt2Aa1) optimized for heterologous expression in P. pastoris. This was achieved by increasing the overall G+C content of the bacterial cyt2Aa1 while changing its codon usage to that preferred by the methylotrophic yeast. Here we describe in detail the design, synthesis and requisite PCR repair of syncyt2Aa1, then present analyses of recombinant Cyt2Aa1 expression in P. pastoris using this synthetic gene. Following the results presented in this paper, the syncyt2Aa1 gene was also successfully used for the recombinant production of a Cyt2Aa1-based IT in the same expression host.

Adjuvants, Immunologic↗
PubMed1992

Rapid and direct detection of the most frequent Mediterranean beta-thalassemic mutations by multiplex allele-specific enzymatic amplification.

A rapid nonradioactive method for the diagnosis of the most frequent Mediterranean beta-thalassemic mutations is described based on a multiplex allele-specific polymerase chain reaction (PCR). This method allows direct detection of normal or mutated alleles on genomic DNA. We have used this approach to detect the most frequent Mediterranean mutations: IVS-1 nt 110 (G----A) and 39 nonsense (C----T). For each mutation three allele-specific oligonucleotides were used: one common upstream primer and two downstream primers differing in their terminal 3' nucleotide (one specific for the normal allele and one for the mutant allele). For each sample two PCR reactions were performed in parallel using in one case IVS-1 nt 110 and codon 39 normal primers and in the second case using the corresponding mutated primers. In both cases the different PCR fragments were visualized. After optimization these primers directed only amplification of their complementary allele. A single blind study was performed on the DNA of 18 individuals who were homozygous or heterozygous for these mutations. In comparison with a parallel investigation, using oligonucleotide probes, all the results were unambiguous. This diagnosis method, which is rapid, easy, direct, and inexpensive, allows the screening of a population group, including heterozygotes, which is required from an epidemiological and anthropological point of view. It could be extended to the large series screening of haplotypes before targeted diagnosis of various genetic diseases.

Algeria↗
PubMedDate not supplied

[Cloning of the gene for thermostable Thermus aquaticus YT1 DNA polymerase and its expression in Escherichia coli].

Using the phasmid vector pSL5, the genomic DNA fragment of T. aquaticus YT1 which contained the thermostable DNA polymerase (Taq-polymerase) gene was cloned. The BglII fragment of this genome locus was subcloned in the BamHI site of the pUC19 plasmid. To optimize the Taq-polymerase gene expression in E. coli cells, the gene was cloned in the correct reading frame regarding the initiation ATG codon of the pPR-TGATG-1 expression vector. The gene expression in this vector was controlled by the phage lambda PR promoter and the temperature-sensitive phage lambda repressor. We used PCR to amplify the short 5'-end fragment of the Taq-polymerase gene coding for the part into which an artificial SacI site was introduced. This site has been used for cloning the PCR product into the pPR-TGATG-1 vector, and the missing gene part was cloned into the KpnI site of the PCR product from the natural cloned gene. The cells of the E. coli PVG-A1 strain, which was obtained in the end, expressed efficiently the Taq-polymerase gene at the nonpermissive temperature. The content of the recombinant Taq-polymerase in the cells was about 1-2% of total proteins. The purified nearly homogeneous Taq-polymerase amplified efficiently in the PCR DNA fragments up to 5.5 kb long and was useful in DNA sequencing the by Sanger method. The half-life of the purified Taq polymerase was about 60 min at 95 degrees C, it was active for at least 65 standard PCR circles. The specific activity of recombinant enzyme preparations was about 180-200,000 units per mg of protein. The E. coli PVG-A1 strain enables one to isolate up to 500,000 units of purified enzyme from 2 l of bacterial culture.

Bacteriophage lambda↗
PubMed2001

Efficient translation initiation is required for replication of bovine viral diarrhea virus subgenomic replicons.

An internal ribosome entry site (IRES) mediates translation initiation of bovine viral diarrhea virus (BVDV) RNA. Studies have suggested that a portion of the N(pro) open reading frame (ORF) is required, although its exact function has not been defined. Here we show that a subgenomic (sg) BVDV RNA in which the NS3 ORF is preceded only by the 5' nontranslated region did not replicate to detectable levels following transfection. However, RNA synthesis and cytopathic effects were observed following serial passage in the presence of a noncytopathic helper virus. Five sg clones derived from the passaged virus contained an identical, silent substitution near the beginning of the NS3 coding sequence (G400U), as well as additional mutations. Four of the reconstructed mutant RNAs replicated in transfected cells, and in vitro translation showed increased levels of NS3 for the mutant RNAs compared to that of wild-type (wt) MetNS3. To more precisely dissect the role of these mutations, we constructed two sg derivatives: ad3.10, which contains only the G400U mutation, and ad3.7, with silent substitutions designed to minimize RNA secondary structure downstream of the initiator AUG. Both RNAs replicated and were translated in vitro to similar levels. Moreover, ad3.7 and ad3.10, but not wt MetNS3, formed toeprints downstream of the initiator AUG codon in an assay for detecting the binding of 40S ribosomal subunits and 43S ribosomal complexes to the IRES. These results suggest that a lack of stable RNA secondary structure(s), rather than a specific RNA sequence, immediately downstream of the initiator AUG is important for optimal translation initiation of pestivirus RNAs.

Amino Acid Sequence↗
PubMed2000

Use of mitogenomic information in teleostean molecular phylogenetics: a tree-based exploration under the maximum-parsimony optimality criterion.

We explored the phylogenetic utility and limits of the individual and concatenated mitochondrial genes for reconstructing the higher-level relationships of teleosts, using the complete (or nearly complete) mitochondrial DNA sequences of eight teleosts (including three newly determined sequences), whose relative phylogenetic positions were noncontroversial. Maximum-parsimony analyses of the nucleotide and amino acid sequences of 13 protein-coding genes from the above eight teleosts, plus two outgroups (bichir and shark), indicated that all of the individual protein-coding genes, with the exception of ND5, failed to recover the expected phylogeny, although unambiguously aligned sequences from 22 concatenated transfer RNA (tRNA) genes (stem regions only) recovered the expected phylogeny successfully with moderate statistical support. The phylogenetic performance of the 13 protein-coding genes in recovering the expected phylogeny was roughly classified into five groups, viz. very good (ND5, ND4, COIII, COI), good (COII, cyt b), medium (ND3, ND2), poor (ND1, ATPase 6), and very poor (ND4L, ND6, ATPase 8). Although the universality of this observation was unclear, analysis of successive concatenation of the 13 protein-coding genes in the same ranking order revealed that the combined data sets comprising nucleotide sequences from the several top-ranked protein-coding genes (no 3rd codon positions) plus the 22 concatenated tRNA genes (stem regions only) best recovered the expected phylogeny, with all internal branches being supported by bootstrap values >90%. We conclude that judicious choice of mitochondrial genes and appropriate data weighting, in conjunction with purposeful taxonomic sampling, are prerequisites for resolving higher-level relationships in teleosts under the maximum-parsimony optimality criterion.

Animals↗
PubMed2001

High guanine-cytosine content is not an adaptation to high temperature: a comparative analysis amongst prokaryotes.

The causes of the variation between genomes in their guanine (G) and cytosine (C) content is one of the central issues in evolutionary genomics. The thermal adaptation hypothesis conjectures that, as G:C pairs in DNA are more thermally stable than adenonine:thymine pairs, high GC content may he a selective response to high temperature. A compilation of data on genomic GC content and optimal growth temperature for numerous prokaryotes failed to demonstrate the predicted correlation. By contrast, the GC content of Structural RNAs is higher at high temperatures. The issue that we address here is whether more freely evolving sites in exons (i.e. codonic third positions) evolve in the same manner as genomic DNA as a whole, Showing no correlated response, or like structural RNAs showing a strong correlation. The latter pattern would provide strong support for the thermal adaptation hypothesis, as the variation in GC content between orthologous genes is typically most profoundly seen at codon third sites (GC3). Simple analysis of completely sequenced prokaryotic genomes shows that GC3, but not genomic GC, is higher on average in thermophilic species. This demonstrates, if nothing else, that the results from the two measures cannot be presumed to be the same. A proper analysis, however, requires phylogenetic control. Here, therefore, we report the results of a comparative analysis of GC composition and optimal growth temperature for over 100 prokaryotes. Comparative analysis fails to show, in either Archea or Eubacteria, any hint of connection between optimal growth temperature and GC content in the genome as a whole, in protein-coding regions or, more crucially at GC. Conversely, comparable analysis confirms that GC content of structural RNA is strongly correlated with optimal temperature. Against the expectations of the thermal adaptation hypothesis, within prokaryotes GC content in protein-coding genies, even at relatively freely evolving sites, cannot be considered an adaptation to the thermal environment.

Adaptation, Physiological↗