Local stability of DNA and RNA secondary structure and its relation to biological functions

Biomarkers of food intake and their relevance to metabolic syndrome

Metabolic syndrome (MetS) constitutes a prevalent risk factor associated with non communicable diseases such as cardiovascular disease and type 2 diabetes. A major factor impacting the etiology of MetS is diet. Dietary patterns and several individual food constituents have been related to the risk of developing MetS or have been proposed as adjuvant treatment. However, traditional methods of dietary assessment such as 24 h recalls rely greatly on intensive user-interaction and are subject to bias. Hence, more objective methods are required for unbiased dietary assessment and efficient prevention. While it is accepted that some dietary-derived constituents in blood plasma are indicators for certain dietary patterns, these may be too unstable (such as vitamin C as a marker for fruits/vegetables) or too broad (e.g. polyphenols for plant-based diets) or reflect too short-term intake only to allow for strong associations with prolonged intake of individual food groups. In the present manuscript, commonly employed biomarkers of intake including those related to specific food items (e.g. genistein for soybean or astaxanthin and EPA for fish intake) and novel emerging ones (e.g. stable isotopes for meat intake or microRNA for plant foods) are emphasized and their suitability as biomarker for food intake discussed. Promising alternatives to plasma measures (e.g. ethyl glucuronide in hair for ethanol intake) are also emphasized. As many biomarkers (i.e. secondary plant metabolites) are not limited to dietary assessment but are also capable of regulating e.g. anti-inflammatory and antioxidant pathways, special attention will be given to biomarkers presenting a double function to assess both dietary patterns and MetS risk.

Variation in genomic traits of microbial communities among ecosystems

Free-living bacteria in nutrient limited environments often exhibit traits which may reduce the cost of reproduction, such as smaller genome size, low GC content and fewer sigma (σ) factor and 16S rRNA gene copies. Despite the potential utility of these traits to detect relationships between microbial communities and ecosystem-scale properties, few studies have assessed these traits on a community-scale. Here, we analysed these traits from publicly available metagenomes derived from marine, soil, host-associated and thermophilic communities. In marine and thermophilic communities, genome size and GC content declined in parallel, consistent with genomic streamlining, with GC content in thermophilic communities generally higher than in marine systems. In contrast, soil communities averaging smaller genomes featured higher GC content and were often from low-carbon environments, suggesting unique selection pressures in soil bacteria. The abundance of specific σ-factors varied with average genome size and ecosystem type. In oceans, abundance of fliA, a σ-factor controlling flagella biosynthesis, was positively correlated with community average genome size—reflecting known trade-offs between nutrient conservation and chemotaxis. In soils, a high abundance of the stress response σ-factor gene rpoS was associated with smaller average genome size and often located in harsh and/or carbon-limited environments—a result which tracks features observed in culture and indicates an increased capacity for stress response in nutrient-poor soils. This work shows how ecosystem-specific constraints are associated with trade-offs which are embedded in the genomic features of bacteria in microbial communities, and which can be detected at the community level, highlighting the importance of genomic features in microbial community analysis.

On the epigenetic role of guanosine oxidation

Chemical modifications of DNA and RNA regulate genome functions or trigger mutagenesis resulting in aging or cancer. Oxidations of macromolecules, including DNA, are common reactions in biological systems and often part of regulatory circuits rather than accidental events. DNA alterations are particularly relevant since the unique role of nuclear and mitochondrial genome is coding enduring and inheritable information. Therefore, an alteration in DNA may represent a relevant problem given its transmission to daughter cells. At the same time, the regulation of gene expression allows cells to continuously adapt to the environmental changes that occur throughout the life of the organism to ultimately maintain cellular homeostasis. Here we review the multiple ways that lead to DNA oxidation and the regulation of mechanisms activated by cells to repair this damage. Moreover, we present the recent evidence suggesting that DNA damage caused by physiological metabolism acts as epigenetic signal for regulation of gene expression. In particular, the predisposition of guanine to oxidation might reflect an adaptation to improve the genome plasticity to redox changes.

High temperature-induced proteomic and metabolomic profiles of a thermophilic Bacillus manusensis isolated from the deep-sea hydrothermal field of Manus Basin

Thermophiles are organisms that grow optimally at 50 °C-80 °C and studies on the survival mechanisms of thermophiles have drawn great attention. Bacillus manusensis S50-6 is the type strain of a new thermophilic species isolated from hydrothermal vent in Manus Basin. In this study, we examined the growth and global responses of S50-6 to high temperature on molecular level using multi-omics method (genomics, proteomics, and metabolomics). S50-6 grew optimally at 50 °C (Favorable, F) and poorly at 65 °C (Non-Favorable, NF); it formed spores at F but not at NF condition. At NF condition, S50-6 formed long filaments containing undivided cells. A total of 1621 proteins were identified at F and NF conditions, and 613 proteins were differentially expressed between F and NF. At NF condition, proteins of glycolysis, rRNA mature and modification, and DNA/protein repair were up-regulated, whereas proteins of sporulation and amino acid/nucleotide metabolism were down-regulated. Consistently, many metabolites associated with amino acid and nucleotide metabolic processes were down-regulated at NF condition. Our results revealed molecular strategies of deep-sea B. manusensis to survive at unfavorable high temperature and provided new insights into the thermotolerant mechanisms of thermophiles. SIGNIFICANCE: In this study, we systematically characterized the genomic, proteomic and metabolomic profiles of a thermophilic deep-sea Bacillus manusensis under different temperatures. Based on these analysis, we propose a model delineating the global responses of B. manusensis to unfavorable high temperature. Under unfavorable high temperature, glycolysis is a more important energy supply pathway; protein synthesis is subjected to more stringent regulation by increased tRNA modification; protein and DNA repair associated proteins are enhanced in production to promote heat survival. In contrast, energy-costing pathways, such as sporulation, are repressed, and basic metabolic pathways, such as amino acid and nucleotide metabolisms, are slowed down. Our results provide new insights into the thermotolerant mechanisms of thermophilic Bacillus.

GC bias lead to increased small amino acids and random coils of proteins in cold-water fishes

BACKGROUND

Temperature adaptation of biological molecules is fundamental in evolutionary studies but remains unsolved. Fishes living in cold water are adapted to low temperatures through adaptive modification of their biological molecules, which enables their functioning in extreme cold. To study nucleotide and amino acid preference in cold-water fishes, we investigated the substitution asymmetry of codons and amino acids in protein-coding DNA sequences between cold-water fishes and tropical fishes., The former includes two Antarctic fishes, Dissostichus mawsoni (Antarctic toothfish), Gymnodraco acuticeps (Antarctic dragonfish), and two temperate fishes, Gadus morhua (Atlantic cod) and Gasterosteus aculeatus (stickleback), and the latter includes three tropical fishes, including Danio rerio (zebrafish), Oreochromis niloticus (Nile tilapia) and Xiphophorus maculatus (Platyfish).

RESULTS

Cold-water fishes showed preference for Guanines and cytosines (GCs) in both synonymous and nonsynonymous codon substitution when compared with tropical fishes. Amino acids coded by GC-rich codons are favored in the temperate fishes, while those coded by AT-rich codons are disfavored. Similar trends were discovered in Antarctic fishes but were statistically weaker. The preference of GC rich codons in nonsynonymous substitution tends to increase ratio of small amino acid in proteins, which was demonstrated by biased small amino acid substitutions in the cold-water species when compared with the tropical species, especially in the temperate species. Prediction and comparison of secondary structure of the proteomes showed that frequency of random coils are significantly larger in the cold-water fish proteomes than those of the tropical fishes.

CONCLUSIONS

Our results suggested that natural selection in cold temperature might favor biased GC content in the coding DNA sequences, which lead to increased frequency of small amino acids and consequently increased random coils in the proteomes of cold-water fishes.

Surprising Conformers of the Biologically Important A·T DNA Base Pairs: QM/QTAIM Proofs

For the first time novel high-energy conformers-A·T(wWC) (5.36), A·T(wrWC) (5.97), A·T(wH) (5.78), and A·T(wrH) (ΔG = 5.82 kcal·mol-1) (See Graphical Abstract) were revealed for each of the four biologically important A·T DNA base pairs - Watson-Crick A·T(WC), reverse Watson-Crick A·T(rWC), Hoogsteen A·T(H) and reverse Hoogsteen A·T(rH) at the MP2/aug-cc-pVDZ//B3LYP/6-311++G(d,p) level of quantum-mechanical theory in the continuum with ε = 4 under normal conditions. Each of these conformers possesses substantially non-planar wobble (w) structure and is stabilized by the participation of the two anti-parallel N6H/N6H'…O4/O2 and N3H…N6 H-bonds, involving the pyramidalized amino group of the A DNA base as an acceptor and a donor of the H-bonding. The transition states - TSA·T(WC)↔A·T(wWC), TSA·T(rWC)↔A·T(wrWC), TSA·T(H)↔A·T(wH), and TSA·T(rH)↔A·T(wrH), controlling the dipole-active transformations of the conformers from the main plane-symmetric state into the high-energy, significantly non-planar state and vice versa, were localized. They also possess wobble structures similarly to the high-energy conformers and are stabilized by the participation of the N6H/N6H'…O4/O2 and N3H…N6 H-bonds. Discovered conformers of the A·T DNA base pairs are dynamically stable short-lived structures [lifetime τ = (1.4-3.9) ps]. Their possible biological significance and future perspectives have been briefly discussed.

Improved biodegradation of textile dye effluent by coculture

The present study demonstrates the de-colorization and degradation of textile effluent by coculture consisting of three bacterial species isolated from textile effluent contaminated environment with an aim to reduce the treatment time. The isolates were identified as Ochrobactrum sp., Pseudomonas aeruginosa and Providencia vermicola by 16S rRNA analysis. Their secondary structure was predicted and GC content of the sequence was found to be 54.39, 52.10, and 52.53%. The co-culture showed a prominent increase in the degradation activity due to the action of oxidoreductase enzymatic mechanism of laccase, NADH-DCIP reductase and azoreductase activity. The biodegradability index of 0.75 was achieved with 95% chemical oxygen demand (COD) reduction in 16 h and 78 and 85% reduction in total organic carbon (TOC) and total solids was observed. Bioaccumulation of metals was identified by X-ray diffraction (XRD) analysis. The effective decolorization was confirmed from the results of UV-vis spectroscopy, high performance liquid chromatography and Fourier transformed infrared spectrometer analyzes. The possible degradation pathway was obtained from the analysis of liquid chromatography-mass spectroscopy analysis and the metabolites such as 2-amino naphthalene and N-phenyl-1.3,5 triazine were observed. The toxic nature of the effluent was analyzed using phyto-toxicity, cell-death assay and geno-toxicity tests.

A statistical thermodynamic model for investigating the stability of DNA sequences from oligonucleotides to genomes

We describe the development and testing of a simple statistical mechanics methodology for duplex DNA applicable to sequences of any composition and extensible to genomes. The microstates of a DNA sequence are modeled in terms of blocks of basepairs that are assumed to be fully closed (paired) or open. This approach generates an ensemble of bubblelike microstates that are used to calculate the corresponding partition function. The energies of the microstates are calculated as additive contributions from hydrogen bonding, basepair stacking, and solvation terms parameterized from a comprehensive series of molecular dynamics simulations including solvent and ions. Thermodynamic properties and nucleotide stability constants for DNA sequences follow directly from the partition function. The methodology was tested by comparing computed free energies per basepair with the experimental melting temperatures of 60 oligonucleotides, yielding a correlation coefficient of -0.96. The thermodynamic stability of genic/nongenic regions was tested in terms of nucleotide stability constants versus sequence for the Escherichia coli K-12 genome. It showed clear differentiation of the genes from promoters and captures genic regions with a sensitivity of 0.94. The statistical thermodynamic model presented here provides a seemingly new handle on the challenging problem of interpreting genomic sequences.

The separation between the 5'-3' ends in long RNA molecules is short and nearly constant

RNA molecules play different roles in coding, decoding and gene expression regulation. Such roles are often associated to the RNA secondary or tertiary structures. The folding dynamics lead to multiple secondary structures of long RNA molecules, since an RNA molecule might fold into multiple distinct native states. Despite an ensemble of different structures, it has been theoretically proposed that the separation between the 5′ and 3′ ends of long single-stranded RNA molecules (ssRNA) remains constant, independent of their base content and length. Here, we present the first experimental measurements of the end-to-end separation in long ssRNA molecules. To determine this separation, we use single molecule Fluorescence Resonance Energy Transfer of fluorescently end-labeled ssRNA molecules ranging from 500 to 5500 nucleotides in length, obtained from two viruses and a fungus. We found that the end-to-end separation is indeed short, within 5–9 nm. It is remarkable that the separation of the ends of all RNA molecules studied remains small and similar, despite the origin, length and differences in their secondary structure. This implies that the ssRNA molecules are ‘effectively circularized’ something that might be a general feature of RNAs, and could result in fine-tuning for translation and gene expression regulation.

Fluctuations in the DNA double helix: A critical review

A critical overview of the extensive literature on fluctuations in the DNA double helix is presented. Both theory and experiment are comprehensively reviewed and analyzed. Fluctuations, which open up the DNA double helix making bases accessible for hydrogen exchange and chemical modification, are the main focus of the review. Theoretical descriptions of the DNA fluctuations are discussed with special emphasis on most popular among them: the nonlinear-dynamic Peyrard-Bishop-Dauxois (PBD) model and the empirical two-state (or helix-coil) model. The experimental data on the issue are comprehensibly overviewed in the historical retrospective with main emphasis on the hydrogen exchange data and formaldehyde kinetics. The theoretical descriptions are critically evaluated from the viewpoint of their applicability to describe DNA in water environment and from the viewpoint of agreement of their predictions with the reliable experimental data. The presented analysis makes it possible to conclude that, while the two-state model is most adequate from theoretical viewpoint and its predictions, based on an empirical parametrization, agree with experimental data very well, the PBD model is inapplicable to DNA in water from theoretical viewpoint on one hand and it makes predictions totally incompatible with reliable experimental data on the other. In particular, it is argued that any oscillation movements of nucleotides, assumed by the PBD model, are severely damped in water, that no "bubbles", which the PBD model predicts, exist in reality in linear DNA well below the melting range and the lifetime of an open state in DNA is actually 5 orders of magnitude longer than the value predicted by the PBD model.

Bayesian prediction of bacterial growth temperature range based on genome sequences

Background

The preferred habitat of a given bacterium can provide a hint of which types of enzymes of potential industrial interest it might produce. These might include enzymes that are stable and active at very high or very low temperatures. Being able to accurately predict this based on a genomic sequence, would thus allow for an efficient and targeted search for production organisms, reducing the need for culturing experiments.

Results

This study found a total of 40 protein families useful for distinction between three thermophilicity classes (thermophiles, mesophiles and psychrophiles). The predictive performance of these protein families were compared to those of 87 basic sequence features (relative use of amino acids and codons, genomic and 16S rDNA AT content and genome size). When using naïve Bayesian inference, it was possible to correctly predict the optimal temperature range with a Matthews correlation coefficient of up to 0.68. The best predictive performance was always achieved by including protein families as well as structural features, compared to either of these alone. A dedicated computer program was created to perform these predictions.

Conclusions

This study shows that protein families associated with specific thermophilicity classes can provide effective input data for thermophilicity prediction, and that the naïve Bayesian approach is effective for such a task. The program created for this study is able to efficiently distinguish between thermophilic, mesophilic and psychrophilic adapted bacterial genomes.

Detection of putatively thermophilic anaerobic methanotrophs in diffuse hydrothermal vent fluids

The anaerobic oxidation of methane (AOM) is carried out by a globally distributed group of uncultivated Euryarchaeota, the anaerobic methanotrophic arachaea (ANME). In this work, we used G+C analysis of 16S rRNA genes to identify a putatively thermophilic ANME group and applied newly designed primers to study its distribution in low-temperature diffuse vent fluids from deep-sea hydrothermal vents. We found that the G+C content of the 16S rRNA genes (P(GC)) is significantly higher in the ANME-1GBa group than in other ANME groups. Based on the positive correlation between the P(GC) and optimal growth temperatures (T(opt)) of archaea, we hypothesize that the ANME-1GBa group is adapted to thrive at high temperatures. We designed specific 16S rRNA gene-targeted primers for the ANME-1 cluster to detect all phylogenetic groups within this cluster, including the deeply branching ANME-1GBa group. The primers were successfully tested both in silico and in experiments with sediment samples where ANME-1 phylotypes had previously been detected. The primers were further used to screen for the ANME-1 microorganisms in diffuse vent fluid samples from deep-sea hydrothermal vents in the Pacific Ocean, and sequences belonging to the ANME-1 cluster were detected in four individual vents. Phylotypes belonging to the ANME-1GBa group dominated in clone libraries from three of these vents. Our findings provide evidence of existence of a putatively extremely thermophilic group of methanotrophic archaea that occur in geographically and geologically distinct marine hydrothermal habitats.

DNA-energetics-based analyses suggest additional genes in prokaryotes

We present here a novel methodology for predicting new genes in prokaryotic genomes on the basis of inherent energetics of DNA. Regions of higher thermodynamic stability were identified, which were filtered based on already known annotations to yield a set of potentially new genes. These were then processed for their compatibility with the stereo-chemical properties of proteins and tripeptide frequencies of proteins in Swissprot data, which results in a reliable set of new genes in a genome. Quite surprisingly, the methodology identifies new genes even in well-annotated genomes. Also, the methodology can handle genomes of any GC-content, size and number of annotated genes.

Characterization of synonymous codon usage bias in the pseudorabies virus US1 gene

In the present study, we examined the codon usage bias between pseudorabies virus (PRV) US1 gene and the US1-like genes of 20 reference alphaherpesviruses. Comparative analysis showed noticeable disparities of the synonymous codon usage bias in the 21 alphaherpesviruses, indicated by codon adaptation index, effective number of codons (ENc) and GC3s value. The codon usage pattern of PRV US1 gene was phylogenetically conserved and similar to that of the US1-like genes of the genus Varicellovirus of alphaherpesvirus, with a strong bias towards the codons with C and G at the third codon position. Cluster analysis of codon usage pattern of PRV US1 gene with its reference alphaherpesviruses demonstrated that the codon usage bias of US1-like genes of 21 alphaherpesviruses had a very close relation with their gene functions. ENc-plot revealed that the genetic heterogeneity in PRV US1 gene and the 20 reference alphaherpesviruses was constrained by G+C content, as well as the gene length. In addition, comparison of codon preferences in the US1 gene of PRV with those of E. coli, yeast and human revealed that there were 50 codons showing distinct usage differences between PRV and yeast, 49 between PRV and human, but 48 between PRV and E. coli. Although there were slightly fewer differences in codon usages between E.coli and PRV, the difference is unlikely to be statistically significant, and experimental studies are necessary to establish the most suitable expression system for PRV US1. In conclusion, these results may improve our understanding of the evolution, pathogenesis and functional studies of PRV, as well as contributing to the area of herpesvirus research or even studies with other viruses.

Evolutionary force in confamiliar marine vertebrates of different temperature realms: adaptive trends in zoarcid fish transcriptomes

BACKGROUND

Studies of temperature-induced adaptation on the basis of genomic sequence data were mainly done in extremophiles. Although the general hypothesis of an increased molecular flexibility in the cold is widely accepted, the results of thermal adaptation are still difficult to detect at proteomic down to the genomic sequence level. Approaches towards a more detailed picture emerge with the advent of new sequencing technologies. Only small changes in primary protein structure have been shown to modify kinetic and thermal properties of enzymes, but likewise for interspecies comparisons a high genetic identity is still essential to specify common principles. The present study uses comprehensive transcriptomic sequence information to uncover general patterns of thermal adaptation on the RNA as well as protein primary structure.

RESULTS

By comparing orthologous sequences of two closely related zoarcid fish inhabiting different latitudinal zones (Antarctica: Pachycara brachycephalum, temperate zone: Zoarces viviparus) we were able to detect significant differences in the codon usage. In the cold-adapted species a lower GC content in the wobble position prevailed for preserved amino acids. We were able to estimate 40-60% coverage of the functions represented within the two compared zoarcid cDNA-libraries on the basis of a reference genome of the phylogenetically closely related fish Gasterosteus aculeatus. A distinct pattern of amino acid substitutions could be identified for the non-synonymous codon exchanges, with a remarkable surplus of serine and reduction of glutamic acid and asparagine for the Antarctic species.

CONCLUSION

Based on the differences between orthologous sequences from confamiliar species, distinguished mainly by the temperature regimes of their habitats, we hypothesize that temperature leaves a signature on the composition of biological macromolecules (RNA, proteins) with implications for the transcription and translation level. As the observed pattern of amino acid substitutions only partly support the flexibility hypothesis further evolutionary forces may be effective at the global transcriptome level.

Does the growth temperature of a prokaryote influence the purine content of its mRNAs?

The formation and breaking of hydrogen bonds between nucleic acid bases are dependent on temperature. The high G+C content of organisms was surmised to be an adaptation for high temperature survival because of the thermal stability of G:C pairs. However, a survey of genomic GC% and optimum growth temperature (OGT) of several prokaryotes revoked any direct relation between them. Significantly high purine (R=A or G) content in mRNAs is also seen as a selective response for survival among thermophiles. Nevertheless, the biological relevance of thermophiles loading their unstable mRNAs with excess purines (purine-loading or R-loading) is not persuasive. Here, we analysed the mRNA sequences from the genomes of 168 prokaryotes (as obtained from NCBI Genome database) with their OGTs ranging from -5 °C to 100 °C to verify the relation between R-loading and OGT. Our analysis fails to demonstrate any correlation between R-loading of the mRNA pool and OGT of a prokaryote. The percentage of purine-loaded mRNAs in prokaryotes is found to be in a rough negative correlation with the genomic GC% (r(2)=0.655, slope=-1.478, P<000.1). We conclude that genomic GC% and bias against certain combinations of nucleotides drive the mRNA-synonymous (sense) strands of DNA towards variations in R-loading.

REGIONS OF RELATIVE GC% UNIFORMITY ARE RECOMBINATIONAL ISOLATORS

Among species within a phylogenetic group, genomic GC% values can cover a wide range that is particularly evident at third codon positions. However, among genes within a genome, genic GC% values can also cover a wide range that is, again, particularly evident at third codon positions. Individual genes and genomes each have a "homostabilizing propensity" to adopt a relatively uniform GC%. Each gene (a "microisochore") occupies a discrete GC% niche of relatively uniform base composition amongst its fellow genes, which can collectively span a wide GC% range. Homostabilization serves to recombinationally isolate both genome sectors (facilitating gene duplication and differentiation) and genomes (facilitating genome duplication and differentiation; e.g., speciation). Although they may sometimes be in conflict, the individualities of genomes, and of genes within those genomes, are separately sustained by a common mechanism, uniformity of GC%. The protection against inadvertent recombination afforded by GC% differentiation is, in the general case, a prerequisite for phenotypic differentiation.

Amino acid composition in endothermic vertebrates is biased in the same direction as in thermophilic prokaryotes

Background

Among bacteria and archaea, amino acid usage is correlated with habitat temperatures. In particular, protein surfaces in species thriving at higher temperatures appear to be enriched in amino acids that stabilize protein structure and depleted in amino acids that decrease thermostability. Does this observation reflect a causal relationship, or could the apparent trend be caused by phylogenetic relatedness among sampled organisms living at different temperatures? And do proteins from endothermic and exothermic vertebrates show similar differences?

Results

We find that the observed correlations between the frequencies of individual amino acids and prokaryotic habitat temperature are strongly influenced by evolutionary relatedness between the species analysed; however, a proteome-wide bias towards increased thermostability remains after controlling for phylogeny. Do eukaryotes show similar effects of thermal adaptation? A small shift of amino acid usage in the expected direction is observed in endothermic ('warm-blooded') mammals and chicken compared to ectothermic ('cold-blooded') vertebrates with lower body temperatures; this shift is not simply explained by nucleotide usage biases.

Conclusion

Protein homologs operating at different temperatures have different amino acid composition, both in prokaryotes and in vertebrates. Thus, during the transition from ectothermic to endothermic life styles, the ancestors of mammals and of birds may have experienced weak genome-wide positive selection to increase the thermostability of their proteins.

A phenomenological model for predicting melting temperatures of DNA sequences

We report here a novel method for predicting melting temperatures of DNA sequences based on a molecular-level hypothesis on the phenomena underlying the thermal denaturation of DNA. The model presented here attempts to quantify the energetic components stabilizing the structure of DNA such as base pairing, stacking, and ionic environment which are partially disrupted during the process of thermal denaturation. The model gives a Pearson product-moment correlation coefficient (r) of approximately 0.98 between experimental and predicted melting temperatures for over 300 sequences of varying lengths ranging from 15-mers to genomic level and at different salt concentrations. The approach is implemented as a web tool (www.scfbio-iitd.res.in/chemgenome/Tm_predictor.jsp) for the prediction of melting temperatures of DNA sequences.

8-Oxo-7,8-dihydroguanine: links to gene expression, aging, and defense against oxidative stress

The one-electron oxidation product of guanine, 8-oxo-7,8-dihydroguanine (8-oxoG), is an abundant lesion in genomic, mitochondrial, and telomeric DNA and RNA. It is considered to be a marker of oxidative stress that preferentially accumulates at the 5' end of guanine strings in the DNA helix, in guanine quadruplexes, and in RNA molecules. 8-OxoG has a lower oxidation potential compared to guanine; thus it is susceptible to oxidation/reduction and, along with its redox products, is traditionally considered to be a major mutagenic DNA base lesion. It does not change the architecture of the DNA double helix and it is specifically recognized and excised by 8-oxoguanine DNA glycosylase (OGG1) during the DNA base excision repair pathway. OGG1 null animals accumulate excess levels of 8-oxoG in their genome, yet they do not have shorter life span nor do they exhibit severe pathological symptoms including tumor formation. In fact they are increasingly resistant to inflammation. Here we address the rarely considered significance of 8-oxoG, such as its optimal levels in DNA and RNA under a given condition, essentiality for normal cellular physiology, evolutionary role, and ability to soften the effects of oxidative stress in DNA, and the harmful consequences of its repair, as well as its importance in transcriptional initiation and chromatin relaxation.

Strong heterogeneity in nucleotidic composition and codon bias in the pea aphid (Acyrthosiphon pisum) shown by EST-based coding genome reconstruction

The aim of this study was to analyze patterns of nucleotidic composition and codon usage in the pea aphid genome (Acyrthosiphon pisum). A collection of 60,000 expressed sequence tags (ESTs) in the pea aphid has been used to automatically reconstruct 5809 coding sequences (CDSs), based on similarity with known proteins and on coding style recognition. Reconstructions were manually checked for ribosomal proteins, leading to tentatively reconstruct the nea-complete set of this category. Pea aphid coding sequences showed a shift toward AT (especially at the third codon position) compared to drosophila homologues. Genes with a putative high level of expression (ribosomal and other genes with high EST support) remained more GC3-rich and had a distinct codon usage from bulk sequences: they exhibited a preference for C-ending codons and CGT (for arginine), which thus appeared optimal for translation. However, the discrimination was not as strong as in drosophila, suggesting a reduced degree of translational selection. The space of variation in codon usage for A. pisum appeared to be larger than in drosophila, with a substantial fraction of genes that remained GC3-rich. Some of those (in particular some structural proteins) also showed high levels of codon bias and a very strong preference for C-ending codons, which could be explained either by strong translational selection or by other mechanisms. Finally, genomic traces were analyzed to build 206 fragments containing a full CDS, which allowed studying the correlations between GC contents of coding and those of noncoding (flanking and introns) sequences.

Stacking heterogeneity: a model for the sequence dependent melting cooperativity of duplex DNA

A microscopic Potts-like one-dimensional model with many particle interactions [referred as the generalized model of polypeptide chains (GMPCs)] is developed to investigate cooperativity of DNA sequence dependent melting. For modeling sequence, regular homogeneous sequences were arranged in heterogeneous blocks of various lengths. Within the framework of the GMPC the authors show that the inclusion of stacking interaction heterogeneity relative to homogeneous hydrogen bond interactions leads to an unexpected and quite remarkable increase in melting cooperativity for small blocks. In some cases this tendency persists for long blocks having sharp sequence heterogeneity.

Reinvestigating the codon and amino acid usage of S. cerevisiae genome: a new insight from protein secondary structure analysis

Biased usage of synonymous codons has been elucidated under the perspective of cellular tRNA abundance for quite a long time now. Taking advantage of publicly available gene expression data for Saccharomyces cerevisiae, a systematic analysis of the codon and amino acid usages in two different coding regions corresponding to the regular (helix and strand) as well as the irregular (coil) protein secondary structures, have been performed. Our analyses suggest that apart from tRNA abundance, mRNA folding stability is another major evolutionary force in shaping the codon and amino acid usage differences between the highly and lowly expressed genes in S. cerevisiae genome and surprisingly it depends on the coding regions corresponding to the secondary structures of the encoded proteins. This is obviously a new paradigm in understanding the codon usage in S. cerevisiae. Differential amino acid usage between highly and lowly expressed genes in the regions coding for the irregular protein secondary structure in S. cerevisiae is expounded by the stability of the mRNA folded structure. Irrespective of the protein secondary structural type, the highly expressed genes always tend to encode cheaper amino acids in order to reduce the overall biosynthetic cost of production of the corresponding protein. This study supports the hypothesis that the tRNA abundance is a consequence of and not a reason for the biased usage of amino acid between highly and lowly expressed genes.

Local mutagenic impact of insertions of LTR retrotransposons on the mouse genome

Solitary LTR loci are the predominant form of LTR retrotransposons in most eukaryotic genomes. They originate from recombination between the two LTRs of an ancestral retrovirus and are therefore incapable of transposition. Despite this inactivity, they appear to have a substantial impact on the host genome. Here we use the murine RMER10 LTR family as an example to describe how such elements can reshape regions of the genome through multiple mutations on an evolutionary time scale. Specifically, we use phylogenetic analysis of multiple copies of RMER10 in rodent species, as well as comparisons of orthologous pairs in mouse and rat, to argue that insertions of members of this family have locally induced the emergence of tandem repeat loci as well as many indels. Analysis of structural aspects of these sequences (secondary structures and transcription factors signals) may explain why RMER10 can become endogenous "mutagenic" factors through induction of replication fork blockages and/or error-prone repair of aberrant DNA structures. This hypothesis is also consistent with features of other interspersed repeated elements.

Analysis of Nanoarchaeum equitans genome and proteome composition: indications for hyperthermophilic and parasitic adaptation

Background

Nanoarchaeum equitans, the only known hyperthermophilic archaeon exhibiting parasitic life style, has raised some new questions about the evolution of the Archaea and provided a model of choice to study the genome landmarks correlated with thermo-parasitic adaptation. In this context, we have analyzed the genome and proteome composition of N. equitans and compared the same with those of other mesophiles, hyperthermophiles and obligatory host-associated organisms.

Results

Analysis of nucleotide, codon and amino acid usage patterns in N. equitans indicates the presence of distinct selective constraints, probably due to its adaptation to a thermo-parasitic life-style. Among the conspicuous characteristics featuring its hyperthermophilic adaptation are overrepresentation of purine bases in protein coding sequences, higher GC-content in tRNA/rRNA sequences, distinct synonymous codon usage, enhanced usage of aromatic and positively charged residues, and decreased frequencies of polar uncharged residues, as compared to those in mesophilic organisms. Positively charged amino acid residues are relatively abundant in the encoded gene-products of N. equitans and other hyperthermophiles, which is reflected in their isoelectric point distribution. Pairwise comparison of 105 orthologous protein sequences shows a strong bias towards replacement of uncharged polar residues of mesophilic proteins by Lys/Arg, Tyr and some hydrophobic residues in their Nanoarchaeal orthologs. The traits potentially attributable to the symbiotic/parasitic life-style of the organism include the presence of apparently weak translational selection in synonymous codon usage and a marked heterogeneity in membrane-associated proteins, which may be important for N. equitans to interact with the host and hence, may help the organism to adapt to the strictly host-associated life style. Despite being strictly host-dependent, N. equitans follows cost minimization hypothesis.

Conclusion

The present study reveals that the genome and proteome composition of N. equitans are marked with the signatures of dual adaptation – one to high temperature and the other to obligatory parasitism. While the analysis of nucleotide/amino acid preferences in N. equitans offers an insight into the molecular strategies taken by the archaeon for thermo-parasitic adaptation, the comparative study of the compositional characteristics of mesophiles, hyperthermophiles and obligatory host-associated organisms demonstrates the generality of such strategies in the microbial world.

Sequence-dependent base pair opening in DNA double helix

Preservation of genetic information in DNA relies on shielding the nucleobases from damage within the double helix. Thermal fluctuations lead to infrequent events of the Watson-Crick basepair opening, or DNA "breathing", thus making normally buried groups available for modification and interaction with proteins. Fluctuational basepair opening implies the disruption of hydrogen bonds between the complementary bases and flipping of the base out of the helical stack. Prediction of sequence-dependent basepair opening probabilities in DNA is based on separation of the two major contributions to the stability of the double helix: lateral pairing between the complementary bases and stacking of the pairs along the helical axis. The partition function calculates the basepair opening probability at every position based on the loss of two stacking interactions and one base-pairing. Our model also includes a term accounting for the unfavorable positioning of the exposed base, which proceeds through a formation of a highly constrained small loop, or a ring. Quantitatively, the ring factor is found as an adjustable parameter from the comparison of the theoretical basepair opening probabilities and the experimental data on short DNA duplexes measured by NMR spectroscopy. We find that these thermodynamic parameters suggest nonobvious sequence dependent basepair opening probabilities.

Genomic conflict settled in favour of the species rather than the gene at extreme GC percentage values

Wada and colleagues have shown that, whether prokaryotic or eukaryotic, each gene has a "homostabilising propensity" to adopt a relatively uniform GC percentage (GC%). Accordingly, each gene can be viewed as a "microisochore" occupying a discrete GC% niche of relatively uniform base composition amongst its fellow genes. Although first, second and third codon positions usually differ in GC%, each position tends to maintain a uniform, gene-specific GC% value. Thus, within a genome, genic GC% values can cover a wide range. This is most evident at third codon positions, which are least constrained by amino acid encoding needs. In 1991, Wada and colleagues further noted that, within a phylogenetic group, genomic GC% values can also cover a wide range. This is again most evident at third codon positions. Thus, the dispersion of GC% values among genes within a genome matches the dispersion of GC% values among genomes within a phylogenetic group. Wada described the context-independence of plots of different codon position GC% values against total GC% as a "universal" characteristic. Several studies relate this to recombination. We have confirmed that third codon positions usually relate more to the genes that contain them than to the species. However, in genomes with extreme GC% values (low or high), third codon positions tend to maintain a constant GC%, thus relating more to the species than to the genes that contain them. Genes in an extreme-GC% genome collectively span a smaller GC% range, and mainly rely on first and second codon positions for differentiation as "microisochores". Our results are consistent with the view that differences in GC% serve to recombinationally isolate both genome sectors (facilitating gene duplication) and genomes (facilitating genome duplication, e.g. speciation). In intermediate-GC% genomes, conflict between the needs of the species and the needs of individual genes within that species is minimal. However, in extreme-GC% genomes there is a conflict, which is settled in favour of the species (i.e. group selection) rather than in favour of the gene (genic selection).

The helix-coil transition in heterogeneous double stranded DNA: Microcanonical method

A microscopic Potts-like one-dimensional model with many-particle interactions is developed to construct a statistical mechanical description of the melting of heterogeneous sequence duplex DNA. For this model, referred as the generalized model of polypeptide chains (GMPC), a closed-form expression for the free energy is derived. The characteristic equation of the model enables estimates on the melting temperature and transition interval, consistent with results obtained from more classical approaches. From the characteristic equation of the model, the temperature-dependent statistical weight parameter for helical states is evaluated. This parameter is shown to change throughout the transition from a harmonic form in early regions of the transition to an arithmetic form in later stages. The GMPC is extended to consider the influence of sequence heterogeneity in the melting of duplex DNA.

On the origin of genomic adaptation at high temperature for prokaryotic organisms

For a long time, the central issue of evolutionary genomics was to find out the adaptive strategy of nucleic acid molecules of various microorganisms having different optimal growth temperatures (Topt). Long-standing controversies exist regarding the correlations between genomic G+C content and Topt, and this debate has not been yet settled. We address this problem by considering the fact that adaptation to growth at high temperature requires a coordinated set of evolutionary changes affecting: (i) nucleic acid thermostability and (ii) stability of codon-anticodon interactions. In the present study, we analyzed 16 prokaryotic genomes having intermediate G+C content and widely varying optimal growth temperatures. Results show that elevated growth temperature imposes selective constraints not only on nucleic acid level but also affects the stability of codon-anticodon interaction. We observed a decrease in the frequency of SSC and SSG codons with the increase in Topt to avoid the formation of side-by-side GC base pairs in the codon-anticodon interaction, thereby making it impossible for a genome to increase GC composition uniformly through the whole coding sequence. Thus, we suggest that any attempt to obtain a generalized relation between genomic GC composition and optimal growth temperature would hardly evolve any satisfactory result.

Mutation exposed: a neutral explanation for extreme base composition of an endosymbiont genome

The influence of neutral mutation pressure versus selection on base composition evolution is a subject of considerable controversy. Yet the present study represents the first explicit population genetic analysis of this issue in prokaryotes, the group in which base composition variation is most dramatic. Here, we explore the impact of mutation and selection on the dynamics of synonymous changes in Buchnera aphidicola, the AT-rich bacterial endosymbiont of aphids. Specifically, we evaluated three forms of evidence. (i) We compared the frequencies of directional base changes (AT-->GC vs. GC-->AT) at synonymous sites within and between Buchnera species, to test for selective preference versus effective neutrality of these mutational categories. Reconstructed mutational changes across a robust intraspecific phylogeny showed a nearly 1:1 AT-->GC:GC-->AT ratio. Likewise, stationarity of base composition among Buchnera species indicated equal rates of AT-->GC and GC-->AT substitutions. The similarity of these patterns within and between species supported the neutral model. (ii) We observed an equivalence of relative per-site AT mutation rate and current AT content at synonymous sites, indicating that base composition is at mutational equilibrium. (iii) We demonstrated statistically greater equality in the frequency of mutational categories in Buchnera than in parallel mammalian studies that documented selection on synonymous sites. Our results indicate that effectively neutral mutational pressure, rather than selection, represents the major force driving base composition evolution in Buchnera. Thus they further corroborate recent evidence for the critical role of reduced N(e) in the molecular evolution of bacterial endosymbionts.

Characterization of Aquifex aeolicus RNase E/G

The RNase E/G homologue from the thermophilic eubacterium Aquifex aeolicus has been overexpressed in Escherichia coli, purified, and characterized in vitro. We show that A. aeolicus RNase E/G has a temperature-dependent, endoribonucleolytic activity. The enzyme site-specifically cleaves oligonucleotides and structured RNAs at locations that are partly overlapping or completely different when compared to the positions of E. coli RNase E and RNase G cleavage sites. The efficiency of cleavage by A. aeolicus RNase E/G is dependent on the 5'-phosphorylation status of RNA suggesting differential susceptibility of primary transcripts and their degradative intermediates to the nuclease activity of this enzyme in vivo. Similar to E. coli RNase E, A. aeolicus RNase E/G is able to selectively cleave internucleotide bonds in the 3'-5' direction, and to cut in intercistronic regions of putative tRNA precursors, thus suggesting a common function for RNase E/G homologues in eubacteria.

Investigation on the causes of codon and amino acid usages variation between thermophilic Aquifex aeolicus and mesophilic Bacillus subtilis

Base composition, codon usages and amino acid usages have been analyzed by taking 529 orthologous sequences of Aquifex aeolicus and Bacillus subtilis, having different optimal growth temperatures. These two bacteria do not have significant difference in overall GC composition, but GC(1+2) and GC3 levels were found to vary significantly. Significant increments in purine content and GC3 composition have been observed in the coding sequences of Aquifex aeolicus than its Bacillus subtilis counterparts. Correspondence analyses on codon and amino acid usages reveal that variation in base composition actually influences their codon and amino acid usages. Two selection pressures acting on the nucleotide level (GC3 and purine enrichment), causes variation in the amino acid usage differently in different protein secondary structures. Our results suggest that adaptation of amino acid usages in coil structure of Aquifex aeolicus proteins is under the control of both purine increment and GC3 composition, whereas the adaptation of the amino acids in the helical region of thermophilic bacteria is strongly influenced by the purine content. Evolutionary perspectives concerning the temperature adaptation of DNA and protein molecules of these two bacteria have been discussed on the basis of these results.

Correlations between genomic GC levels and optimal growth temperatures: some comments

Regarding the existence of any specific correlation between optimal growth temperature and genomic GC levels, Musto et al. [FEBS Lett. 573 (2004) 73] have recently performed analysis on 20 prokaryotic families and showed that in most of the families there exists a positive correlation between these two parameters. On the basis of these results they claimed that optimal growth temperature is one of the factors that influence genomic GC composition in prokaryotes. In a subsequent article, Marashi and Ghalanbor [Biochem. Biophys. Res. Commun. 325 (2004) 381] have demonstrated that the correlation values change substantially when very few points in some of the families were excluded from the data set of Musto et al. [FEBS Lett. 573 (2004) 73]. But Marashi and Ghalanbor have not provided any reason behind this. The points excluded by Marashi and Ghalanbor are actually the outliers in the data set, which strongly affect the correlation coefficients. But the presence of outliers in large data set hardly had any effect on the correlation values. Marashi and Ghalanbor have excluded points from only those families that have small sample sizes and observed a substantial change in correlation coefficient values. Therefore, we argue that any conclusion drawn for a small sample size having outliers is always questionable. Although Musto's approach is a novel one, but to make any generalization one needs to be careful about the flawlessness in the data set.

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.

Chromosomal speciation: a reply

The "genic" and the "non-genic" (chromosomal) hypotheses for the predominant mechanism by which species diverge into two have long been in contention. In 1998 Coyne and Orr attacked certain formulations of the chromosomal hypothesis on the grounds that they required macromutations (structural changes in chromosomes). In 1999 I replied that numerous independent micromutations (single DNA base changes) should suffice (GC% hypothesis). Kliman et al., with the support of Coyne and Charlesworth, have presented various counterarguments, to which the present paper responds with evidence that GC% differences are primary to genic differences and would operate by changing the structure of stem-loops extruded from duplex DNAs. Chromosomes attempting to align by means of complementary loop-loop interactions would fail if GC% differences exceeded a critical threshold. This would disrupt meiosis (hybrid sterility) and the parents of organisms with failed meiosis would be reproductively isolated from each other. If they could find new mates with which they were GC-compatible, then new species could emerge. The model leads to predictions consistent with several lines of evidence. The GC% version of the chromosomal hypothesis has a sound basis and deserves at least as much attention as its genic rival.

Gene structure and transcriptional regulation of dnaK and dnaJ genes from a psychrophilic bacterium, Colwellia maris

The dnaK and dnaJ genes, encoding heat shock proteins, were cloned from a psychrophilic bacterium, Colwellia maris. Significant homology was evident comparing DnaK and DnaJ of the psychrophilile with the counterparts of mesophilic and thermophilic bacteria. In the DnaJ protein, three conserved regions of the Hsp40 family were observed. A putative promoter similar to the sigma32 consensus sequence was found upstream of the dnaK gene. The G+C content in the 5'-untranslated region of the dnaK gene was much lower than that in the corresponding region of mesophilic bacteria. Northern-blot analysis and primer-extension analysis showed that both genes were transcribed separately as monocistronic mRNAs. Following several temperature upshifts from 10 to 26 degrees C, maximum induction of the dnaK and dnaJ mRNAs was detected at 20 degrees C, suggesting that this temperature induces the heat shock response in this bacterium. In addition, the level of the induction of the dnaJ gene was much lower than that of the dnaK gene. These findings together revealed several specific features of the heat shock response at a relatively low temperature in psychrophiles.

Stacking and hydrogen bonding: DNA cooperativity at melting

By taking into account base-base stacking interactions we improve the Generalized Model of Polypeptide Chain (GMPC). Based on a one-dimensional Potts-like model with many-particle interactions, the GMPC describes the helix-coil transition in both polypeptides and polynucleotides. In the framework of the GMPC we show that correctly introduced nearest-neighbor stacking interactions against the background of hydrogen bonding lead to increased stability (melting temperature) and, unexpectedly, to decreased cooperativity (maximal correlation length). The increase in stability is explained as due to an additional stabilizing interaction (stacking) and the surprising decrease in cooperativity is seen as a result of mixing of contributions of hydrogen bonding and stacking.

The base composition of the genes is correlated with the secondary structures of the encoded proteins

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(3) levels of the genes localized in the GC-richest isochores of the human genome.

Teleost introns are characterized by a high A+T content

We previously observed that Antarctic fish genes contain intron sequences of high A+T content (60-70% average A+T) which are in stark contrast with adjacent protein coding-sequences. Here, we report that this disparity in intron/exon base composition is a common feature among teleosts. We analyzed 483 teleost genomic DNA sequences, containing 2583 introns, from 80 teleost genera that populate polar, temperate, or tropical habitats. Eighty-nine percent of teleost introns display an A+T content between 50-84% A+T with a mean of 60% A+T. In contrast, only 37% of teleost exons have an A+T content greater-than 50% with a mean of 48% A+T. A comparison to homologous mammalian genes showed a striking difference; in this case, introns and exons have similar base compositions, averaging 45-47% A+T. This indicates that most teleost genes exhibit a large difference in base composition between their introns and exons. There was no correlation of teleost intron A+T content to intron length or habitat temperature range. Thus, teleost intron sequences tend to show the common feature of being much higher in A+T content then neighboring exons.

Evidence for strong selective constraint acting on the nucleotide composition of 16S ribosomal RNA genes

Previous studies have shown that the guanine plus cytosine (G+C) content of ribosomal RNAs (rRNAs) is highly correlated with bacterial growth temperatures. This correlation is strongest in the double-stranded stem regions of the rRNA, a fact that can be explained by selection for increased structural stability at high growth temperatures. In this study, we examined the single-stranded regions of 16S rRNAs. We reasoned that, since these regions of the molecule are subject to less structural constraint than the stem regions, their nucleotide content might simply reflect the overall nucleotide content of the genome. Contrary to this expectation, however, we found that all of the single-stranded regions are characterized by very high adenine (A) and relatively low cytosine (C) contents. Moreover, the nucleotide content of these single-stranded regions is surprisingly constant between species, despite dramatic differences in optimal growth temperatures, and despite large differences in the overall genomic G+C content. This provides compelling evidence for strong stabilizing selection acting on 16S rRNA single-stranded regions. We found that selection favors purines (A+G), and especially adenine (A), in the single-stranded regions of these rRNAs.

Evolution of transcript structure and base composition of rDNA expansion segment D3 in ticks

Four to thirty-two copies of the rDNA 28S gene expansion segment D3 and flanking H14 stem were sequenced in six species of ticks (Ixodes: Ixodidae: Acari). Sequence match among species varied from 66% to 97%. Sequence length averaged 130 bases in I. persulcatus across eight Eurasian sites and averaged 186 bases in five other species across 19 Eurasian and North American sites. The difference in length represents one or more deletions totalling about 60 bases that correspond to stems S3 or S4 of the folded transcript. The typical transcript conformation was observed as one possible low energy structure in the five species of longer D3. The structure entails a basal loop with four stem/loop structures, S1-S4 (moving 5' to 3') atop stem H14. A secondary structure lacking S4 but possessing all other putative standard features of the D3 transcript is possible with the shorter I. persulcatus sequences. Interspecific sequence differences occur at higher frequency in loops and bulges vs. complementary pairing regions of stems. Insertion/deletion events (indels) and base substitutions accounted equally for sequence differences. Indels are flanked by similar sequences, suggesting that they occur by slippage during replication. The D3 of Ixodes species is composed of a degenerate set of subrepeats. Thus, unequal exchange among subrepeats may have caused the reduction in length of the I. persulcatus D3. Compensatory base substitution and compensatory insertion/deletion events are indicated by the failure of mutations to affect secondary structure. Transversions accounted for 64% of sequence differences and were biased toward the gain of G and U and the loss of A and C. This bias could re-establish intramolecular base pairing when disrupted by insertions or deletions that shift one side of a stem relative to the other. The distribution of sequence differences, biased substitution, and conservation of transcript conformation in D3 suggest selective constraint.

The relationship between third-codon position nucleotide content, codon bias, mRNA secondary structure and gene expression in the drosophilid alcohol dehydrogenase genes Adh and Adhr

To gain insights into the relationship between codon bias, mRNA secondary structure, third-codon position nucleotide distribution, and gene expression, we predicted secondary structures in two related drosophilid genes, Adh and Adhr, which differ in degree of codon bias and level of gene expression. Individual structural elements (helices) were inferred using the comparative method. For each gene, four types of randomization simulations were performed to maintain/remove codon bias and/or to maintain or alter third-codon position nucleotide composition (N3). In the weakly expressed, weakly biased gene Adhr, the potential for secondary structure formation was found to be much stronger than in the highly expressed, highly biased gene Adh. This is consistent with the observation of approximately equal G and C percentages in Adhr ( approximately 31% across species), whereas in Adh the N3 distribution is shifted toward C (42% across species). Perturbing the N3 distribution to approximately equal amounts of A, G, C, and T increases the potential for secondary structure formation in Adh, but decreases it in Adhr. On the other hand, simulations that reduce codon bias without changing N3 content indicate that codon bias per se has only a weak effect on the formation of secondary structures. These results suggest that, for these two drosophilid genes, secondary structure is a relatively independent, negative regulator of gene expression. Whereas the degree of codon bias is positively correlated with level of gene expression, strong individual secondary structural elements may be selected for to retard mRNA translation and to decrease gene expression.

Features of distamycin preferential binding sites on natural DNA predicted using differential scanning calorimetry

The interaction of distamycin with ColE1 DNA was examined by using differential scanning calorimetry (DSC) taking the helix-coil transition theory of DNA into consideration. Our results here strongly indicate that the affinity of distamycin to DNA, at a low distamycin concentration, depends highly on the DNA sequence, and preferential binding occurs to the sites of four to six successive A-T pairs having two or more successive G-C pairs on both their ends.

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.

The compositional evolution of vertebrate genomes

The compositional evolution of vertebrate genomes is characterized: (i) by one predominant conservative mode, in which nucleotide changes occur, but the base composition of DNA sequences in general, and of coding sequences in particular, does not change; and (ii) by three different shifting or transitional modes, in which nucleotide changes are accompanied by changes in the base composition of sequences. Investigations on these evolutionary modes have shed new light on a central problem in molecular evolution, namely the role played by natural selection in modulating the mutational input. This review will present first the intragenomic shifts, the 'major shifts' and the 'minor shift', and then the 'whole-genome', or 'horizontal', shift. In each case, the shifts were preceded and followed by a conservative mode of evolution. This review expands on a previous one [Bernardi, Gene 241 (2000) 3-17], and summarizes the evidence that the changes of the compositional patterns of the genome and their maintenance are controlled by Darwinian natural selection.