Adaptive role of increased frequency of polypurine tracts in mRNA sequences of thermophilic prokaryotes

High-quality genome of a novel Thermosynechococcaceae species from Namibia and characterization of its protein expression patterns at elevated temperatures

Thermophilic cyanobacteria thrive in extreme environments, making their thermoresistant enzymes valuable for industrial applications. Common habitats include hot springs, which act as evolutionary accelerators for speciation due to geographical isolation. The family Thermosynechococcaceae comprises thermophilic cyanobacteria known for their ability to thrive in high-temperature environments. These bacteria are notable for their photosynthetic capabilities, significantly contributing to primary production in extreme habitats. Members of Thermosynechococcaceae exhibit unique adaptations that allow them to perform photosynthesis efficiently at elevated temperatures, making them subjects of interest for studies on microbial ecology, evolution, and potential biotechnological applications. In this study, the genome of a thermophilic cyanobacterium, isolated from a hot spring near Okahandja in Namibia, was sequenced using a PacBio Sequel IIe long-read platform. Cultivations were performed at elevated temperatures of 40, 50, and 55°C, followed by proteome analyses based on the annotated genome. Phylogenetic investigations, informed by the 16S rRNA gene and aligned nucleotide identity (ANI), suggest that the novel cyanobacterium is a member of the family Thermosynechococcaceae. Furthermore, the new species was assigned to a separate branch, potentially representing a novel genus. Whole-genome alignments supported this finding, revealing few conserved regions and multiple genetic rearrangement events. Additionally, 129 proteins were identified as differentially expressed in a temperature-dependent manner. The results of this study broaden our understanding of cyanobacterial adaptation to extreme environments, providing a novel high-quality genome of Thermosynechococcaceae cyanobacterium sp. Okahandja and several promising candidate proteins for expression and characterization studies.

Extreme makeover: the incredible cell membrane adaptations of extremophiles to harsh environments

The existence of life beyond Earth has long captivated humanity, and the study of extremophiles-organisms surviving and thriving in extreme environments-provides crucial insights into this possibility. Extremophiles overcome severe challenges such as enzyme inactivity, protein denaturation, and damage of the cell membrane by adopting several strategies. This feature article focuses on the molecular strategies extremophiles use to maintain the cell membrane's structure and fluidity under external stress. Key strategies include homeoviscous adaptation (HVA), involving the regulation of lipid composition, and osmolyte-mediated adaptation (OMA), where small organic molecules protect the lipid membrane under stress. Proteins also have direct and indirect roles in protecting the lipid membrane. Examining the survival strategies of extremophiles provides scientists with crucial insights into how life can adapt and persist in harsh conditions, shedding light on the origins of life. This article examines HVA and OMA and their mechanisms in maintaining membrane stability, emphasizing our contributions to this field. It also provides a brief overview of the roles of proteins and concludes with recommendations for future research directions.

Idiosyncratic genome evolution of the thermophilic cyanobacterium Synechococcus at the limits of phototrophy

Thermophilic microorganisms are expected to have smaller cells and genomes compared with mesophiles, a higher proportion of horizontally acquired genes, and distinct nucleotide and amino acid composition signatures. Here, we took an integrative approach to investigate these apparent correlates of thermophily for Synechococcus A/B cyanobacteria, which include the most heat-tolerant phototrophs on the planet. Phylogenomics confirmed a unique origin of different thermotolerance ecotypes, with low levels of continued gene flow between ecologically divergent but overlapping populations, which has shaped the distribution of phenotypic traits along these geothermal gradients. More thermotolerant strains do have smaller genomes, but genome reduction is associated with a decrease in community richness and metabolic diversity, rather than with cell size. Horizontal gene transfer played only a limited role during Synechococcus evolution, but, the most thermotolerant strains have acquired a Thermus tRNA modification enzyme that may stabilize translation at high temperatures. Although nucleotide base composition was not associated with thermotolerance, we found a general replacement of aspartate with glutamate, as well as a dramatic remodeling of amino acid composition at the highest temperatures that substantially differed from previous predictions. We conclude that Synechococcus A/B genome diversification largely does not conform to the standard view of temperature adaptation. In addition, carbon fixation was more thermolabile than photosynthetic oxygen evolution for the most thermotolerant strains compared with less tolerant lineages. This suggests that increased flow of reducing power generated during the light reactions to an electron sink(s) beyond carbon dioxide has emerged during temperature adaptation of these bacteria.

Metagenomic Thermometer

Various microorganisms exist in environments, and each of them has its optimal growth temperature (OGT). The relationship between genomic information and OGT of each species has long been studied, and one such study revealed that OGT of prokaryotes can be accurately predicted based on the fraction of seven amino acids (IVYWREL) among all encoded amino-acid sequences in its genome. Extending this discovery, we developed a 'Metagenomic Thermometer' as a means of predicting environmental temperature based on metagenomic sequences. Temperature prediction of diverse environments using publicly available metagenomic data revealed that the Metagenomic Thermometer can predict environmental temperatures with small temperature changes and little influx of microorganisms from other environments. The accuracy of the Metagenomic Thermometer was also confirmed by a demonstration experiment using an artificial hot water canal. The Metagenomic Thermometer was also applied to human gut metagenomic samples, yielding a reasonably accurate value for human body temperature. The result further suggests that deep body temperature determines the dominant lineage of the gut community. Metagenomic Thermometer provides a new insight into temperature-driven community assembly based on amino-acid composition rather than microbial taxa.

Main Factors Shaping Amino Acid Usage Across Evolution

The standard genetic code determines that in most species, including viruses, there are 20 amino acids that are coded by 61 codons, while the other three codons are stop triplets. Considering the whole proteome each species features its own amino acid frequencies, given the slow rate of change, closely related species display similar GC content and amino acids usage. In contrast, distantly related species display different amino acid frequencies. Furthermore, within certain multicellular species, as mammals, intragenomic differences in the usage of amino acids are evident. In this communication, we shall summarize some of the most prominent and well-established factors that determine the differences found in the amino acid usage, both across evolution and intragenomically.

The genomic basis of host and vector specificity in non-pathogenic trypanosomatids

Trypanosoma theileri, a non-pathogenic parasite of bovines, has a predicted surface protein architecture that likely aids survival in its mammalian host. Their surface proteins are encoded by genes which account for ∼10% of their genome. A non-pathogenic parasite of sheep, Trypanosoma melophagium, is transmitted by the sheep ked and is closely related to T. theileri. To explore host and vector specificity between these species, we sequenced the T. melophagium genome and transcriptome and an annotated draft genome was assembled. T. melophagium was compared to 43 kinetoplastid genomes, including T. theileri. T. melophagium and T. theileri have an AT biased genome, the greatest bias of publicly available trypanosomatids. This trend may result from selection acting to decrease the genomic nucleotide cost. The T. melophagium genome is 6.3Mb smaller than T. theileri and large families of proteins, characteristic of the predicted surface of T. theileri, were found to be absent or greatly reduced in T. melophagium. Instead, T. melophagium has modestly expanded protein families associated with the avoidance of complement-mediated lysis. We propose that the contrasting genomic features of these species is linked to their mode of transmission from their insect vector to their mammalian host. This article has an associated First Person interview with the first author of the paper.

Diversity and distribution of thermophilic microorganisms and their applications in biotechnology

Hot springs ecosystem is the most ancient continuously inhabited ecosystem on earth which harbors diverse thermophilic bacteria and archaea distributed worldwide. Life in extreme environments is very challenging so there is a great potential biological dark matter and their adaptation to harsh environments eventually producing thermostable enzymes which are very vital for the welfare of mankind. There is an enormous need for a new generation of stable enzymes that can endure harsh conditions in industrial processes and can either substitute or complement conventional chemical processes. Here, we review at the variety and distribution of thermophilic microbes, as well as the different thermostable enzymes that help them survive at high temperatures, such as proteases, amylases, lipases, cellulases, pullulanase, xylanases, and DNA polymerases, as well as their special properties, such as high-temperature stability. We have documented the novel isolated thermophilic and hyperthermophilic microorganisms, as well as the discovery of their enzymes, demonstrating their immense potential in the scientific community and in industry.

Crystal structure of a Cu,Zn superoxide dismutase from the thermophilic fungus Chaetomium thermophilum

BACKGROUND

Thermophilic fungi have recently emerged as a promising source of thermostable enzymes. Superoxide dismutases are key antioxidant metalloenzymes with promising therapeutic effects in various diseases, both acute and chronic. However, structural heterogeneity and low thermostability limit their therapeutic efficacy.

OBJECTIVE

Although several studies from hypethermophilic superoxide dismutases (SODs) have been reported, information about Cu,Zn-SODs from thermophilic fungi is scarce. Chaetomium thermophilum is a thermophilic fungus that could provide proteins with thermophilic properties.

METHOD

The enzyme was expressed in Pichia pastoris cells and crystallized using the vapor-diffusion method. X-ray data were collected, and the structure was determined and refined to 1.56 Å resolution. Structural analysis and comparisons were carried out.

RESULTS

The presence of 8 molecules (A through H) in the asymmetric unit resulted in four different interfaces. Molecules A and F form the typical homodimer which is also found in other Cu,Zn-SODs. Zinc was present in all subunits of the structure while copper was found in only four subunits with reduced occupancy (C, D, E and F).

CONCLUSION

The ability of the enzyme to form oligomers and the elevated Thr:Ser ratio may be contributing factors to its thermal stability. Two hydrophobic residues that participate in interface formation and are not present in other CuZn-SODs may play a role in the formation of new interfaces and the oligomerization process. The CtSOD crystal structure reported here is the first Cu,Zn-SOD structure from a thermophilic fungus.

PATH - Prediction of Amyloidogenicity by Threading and Machine Learning

Amyloids are protein aggregates observed in several diseases, for example in Alzheimer's and Parkinson's diseases. An aggregate has a very regular beta structure with a tightly packed core, which spontaneously assumes a steric zipper form. Experimental methods enable studying such peptides, however they are tedious and costly, therefore inappropriate for genomewide studies. Several bioinformatic methods have been proposed to evaluate protein propensity to form an amyloid. However, the knowledge of aggregate structures is usually not taken into account. We propose PATH (Prediction of Amyloidogenicity by THreading) - a novel structure-based method for predicting amyloidogenicity and show that involving available structures of amyloidogenic fragments enhances classification performance. Experimental aggregate structures were used in templatebased modeling to recognize the most stable representative structural class of a query peptide. Several machine learning methods were then applied on the structural models, using their energy terms. Finally, we identified the most important terms in classification of amyloidogenic peptides. The proposed method outperforms most of the currently available methods for predicting amyloidogenicity, with its area under ROC curve equal to 0.876. Furthermore, the method gave insight into significance of selected structural features and the potentially most stable structural class of a peptide fragment if subjected to crystallization.

Vesicles From Vibrio cholerae Contain AT-Rich DNA and Shorter mRNAs That Do Not Correlate With Their Protein Products

Extracellular vesicles secreted by Gram-negative bacteria have proven to be important in bacterial defense, communication and host–pathogen relationships. They resemble smaller versions of the bacterial mother cell, with similar contents of proteins, LPS, DNA, and RNA. Vesicles can elicit a protective immune response in a range of hosts, and as vaccine candidates, it is of interest to properly characterize their cargo. Genetic sequencing data is already available for vesicles from several bacterial strains, but it is not yet clear how the genetic makeup of vesicles differ from that of their parent cells, and which properties may characterize enriched genetic material. The present study provides evidence for DNA inside vesicles from Vibrio cholerae O395, and key characteristics of their genetic and proteomic content are compared to that of whole cells. DNA analysis reveals enrichment of fragments containing ToxR binding sites, as well as a positive correlation between AT-content and enrichment. Some mRNAs were highly enriched in the vesicle fraction, such as membrane protein genes ompV, ompK, and ompU, DNA-binding protein genes hupA, hupB, ihfB, fis, and ssb, and a negative correlation was found between mRNA enrichment and transcript length, suggesting mRNA inclusion in vesicles may be a size-dependent process. Certain non-coding and functional RNAs were found to be enriched, such as VrrA, GcvB, tmRNA, RNase P, CsrB2, and CsrB3. Mass spectrometry revealed enrichment of outer membrane proteins, known virulence factors, phage components, flagella and extracellular proteins in the vesicle fraction, and a low, negative correlation was found between transcript-, and protein enrichment. This result opposes the hypothesis that a significant degree of protein translation occurs in vesicles after budding. The abundance of viral-, and flagellar proteins in the vesicle fraction underlines the importance of purification during vesicle isolation.

Selection-driven cost-efficiency optimization of transcripts modulates gene evolutionary rate in bacteria

BACKGROUND

Most amino acids are encoded by multiple synonymous codons. However, synonymous codons are not used equally, and this biased codon use varies between different organisms. It has previously been shown that both selection acting to increase codon translational efficiency and selection acting to decrease codon biosynthetic cost contribute to differences in codon bias. However, it is unknown how these two factors interact or how they affect molecular sequence evolution.

RESULTS

Through analysis of 1320 bacterial genomes, we show that bacterial genes are subject to multi-objective selection-driven optimization of codon use. Here, selection acts to simultaneously decrease transcript biosynthetic cost and increase transcript translational efficiency, with highly expressed genes under the greatest selection. This optimization is not simply a consequence of the more translationally efficient codons being less expensive to synthesize. Instead, we show that transfer RNA gene copy number alters the cost-efficiency trade-off of synonymous codons such that, for many species, selection acting on transcript biosynthetic cost and translational efficiency act in opposition. Finally, we show that genes highly optimized to reduce cost and increase efficiency show reduced rates of synonymous and non-synonymous mutation.

CONCLUSIONS

This analysis provides a simple mechanistic explanation for variation in evolutionary rate between genes that depends on selection-driven cost-efficiency optimization of the transcript. These findings reveal how optimization of resource allocation to messenger RNA synthesis is a critical factor that determines both the evolution and composition of genes.

The Amount of Nitrogen Used for Photosynthesis Modulates Molecular Evolution in Plants

Genome and transcript sequences are composed of long strings of nucleotide monomers (A, C, G, and T/U) that require different quantities of nitrogen atoms for biosynthesis. Here, it is shown that the strength of selection acting on transcript nitrogen content is influenced by the amount of nitrogen plants require to conduct photosynthesis. Specifically, plants that require more nitrogen to conduct photosynthesis experience stronger selection on transcript sequences to use synonymous codons that cost less nitrogen to biosynthesize. It is further shown that the strength of selection acting on transcript nitrogen cost constrains molecular sequence evolution such that genes experiencing stronger selection evolve at a slower rate. Together these findings reveal that the plant molecular clock is set by photosynthetic efficiency, and provide a mechanistic explanation for changes in plant speciation rates that occur concomitant with improvements in photosynthetic efficiency and changes in the environment such as light, temperature, and atmospheric CO2 concentration.

Proteome Evolution of Deep-Sea Hydrothermal Vent Alvinellid Polychaetes Supports the Ancestry of Thermophily and Subsequent Adaptation to Cold in Some Lineages

Temperature, perhaps more than any other environmental factor, is likely to influence the evolution of all organisms. It is also a very interesting factor to understand how genomes are shaped by selection over evolutionary timescales, as it potentially affects the whole genome. Among thermophilic prokaryotes, temperature affects both codon usage and protein composition to increase the stability of the transcriptional/translational machinery, and the resulting proteins need to be functional at high temperatures. Among eukaryotes less is known about genome evolution, and the tube-dwelling worms of the family Alvinellidae represent an excellent opportunity to test hypotheses about the emergence of thermophily in ectothermic metazoans. The Alvinellidae are a group of worms that experience varying thermal regimes, presumably having evolved into these niches over evolutionary times. Here we analyzed 423 putative orthologous loci derived from 6 alvinellid species including the thermophilic Alvinella pompejana and Paralvinella sulfincola. This comparative approach allowed us to assess amino acid composition, codon usage, divergence, direction of residue changes and the strength of selection along the alvinellid phylogeny, and to design a new eukaryotic thermophilic criterion based on significant differences in the residue composition of proteins. Contrary to expectations, the alvinellid ancestor of all present-day species seems to have been thermophilic, a trait subsequently maintained by purifying selection in lineages that still inhabit higher temperature environments. In contrast, lineages currently living in colder habitats likely evolved under selective relaxation, with some degree of positive selection for low-temperature adaptation at the protein level.

Structural signatures of thermal adaptation of bacterial ribosomal RNA, transfer RNA, and messenger RNA

Temperature adaptation of bacterial RNAs is a subject of both fundamental and practical interest because it will allow a better understanding of molecular mechanism of RNA folding with potential industrial application of functional thermophilic or psychrophilic RNAs. Here, we performed a comprehensive study of rRNA, tRNA, and mRNA of more than 200 bacterial species with optimal growth temperatures (OGT) ranging from 4°C to 95°C. We investigated temperature adaptation at primary, secondary and tertiary structure levels. We showed that unlike mRNA, tRNA and rRNA were optimized for their structures at compositional levels with significant tertiary structural features even for their corresponding randomly permutated sequences. tRNA and rRNA are more exposed to solvent but remain structured for hyperthermophiles with nearly OGT-independent fluctuation of solvent accessible surface area within a single RNA chain. mRNA in hyperthermophiles is essentially the same as random sequences without tertiary structures although many mRNA in mesophiles and psychrophiles have well-defined tertiary structures based on their low overall solvent exposure with clear separation of deeply buried from partly exposed bases as in tRNA and rRNA. These results provide new insight into temperature adaptation of different RNAs.

Dietary nitrogen alters codon bias and genome composition in parasitic microorganisms

BACKGROUND

Genomes are composed of long strings of nucleotide monomers (A, C, G and T) that are either scavenged from the organism's environment or built from metabolic precursors. The biosynthesis of each nucleotide differs in atomic requirements with different nucleotides requiring different quantities of nitrogen atoms. However, the impact of the relative availability of dietary nitrogen on genome composition and codon bias is poorly understood.

RESULTS

Here we show that differential nitrogen availability, due to differences in environment and dietary inputs, is a major determinant of genome nucleotide composition and synonymous codon use in both bacterial and eukaryotic microorganisms. Specifically, low nitrogen availability species use nucleotides that require fewer nitrogen atoms to encode the same genes compared to high nitrogen availability species. Furthermore, we provide a novel selection-mutation framework for the evaluation of the impact of metabolism on gene sequence evolution and show that it is possible to predict the metabolic inputs of related organisms from an analysis of the raw nucleotide sequence of their genes.

CONCLUSIONS

Taken together, these results reveal a previously hidden relationship between cellular metabolism and genome evolution and provide new insight into how genome sequence evolution can be influenced by adaptation to different diets and environments.

Insight into thermophiles and their wide-spectrum applications

The deconstruction of biomass is a pivotal process for the manufacture of target products using microbial cells and their enzymes. But the enzymes that possess a significant role in the breakdown of biomass remain relatively unexplored. Thermophilic microorganisms are of special interest as a source of novel thermostable enzymes. Many thermophilic microorganisms possess properties suitable for biotechnological and commercial use. There is, indeed, a considerable demand for a new generation of stable enzymes that are able to withstand severe conditions in industrial processes by replacing or supplementing traditional chemical processes. This manuscript reviews the pertinent role of thermophilic microorganisms as a source for production of thermostable enzymes, factors afftecting them, recent patents on thermophiles and moreso their wide spectrum applications for commercial and biotechnological use.

The survival mechanisms of thermophiles at high temperatures: an angle of omics

Thermophiles are referred to as microorganisms with optimal growth temperatures of >60 °C. Over the past few years, a number of studies have been conducted regarding thermophiles, especially using the omics strategies. This review provides a systematic view of the survival physiology of thermophiles from an "omics" perspective, which suggests that the adaptive ability of thermophiles is based on a cooperative mode with multi-dimensional regulations integrating genomics, transcriptomics, and proteomics.

Sequence determinants of prokaryotic gene expression level under heat stress

Prokaryotic gene expression is environment-dependent and temperature plays an important role in shaping the gene expression profile. Revealing the regulation mechanisms of gene expression pertaining to temperature has attracted tremendous efforts in recent years particularly owning to the yielding of transcriptome and proteome data by high-throughput techniques. However, most of the previous works concentrated on the characterization of the gene expression profile of individual organism and little effort has been made to disclose the commonality among organisms, especially for the gene sequence features. In this report, we collected the transcriptome and proteome data measured under heat stress condition from recently published literature and studied the sequence determinants for the expression level of heat-responsive genes on multiple layers. Our results showed that there indeed exist commonness and consistent patterns of the sequence features among organisms for the differentially expressed genes under heat stress condition. Some features are attributed to the requirement of thermostability while some are dominated by gene function. The revealed sequence determinants of bacterial gene expression level under heat stress complement the knowledge about the regulation factors of prokaryotic gene expression responding to the change of environmental conditions. Furthermore, comparisons to thermophilic adaption have been performed to reveal the similarity and dissimilarity of the sequence determinants for the response to heat stress and for the adaption to high habitat temperature, which elucidates the complex landscape of gene expression related to the same physical factor of temperature.

Assessing diversity of DNA structure-related sequence features in prokaryotic genomes

Prokaryotic genomes are diverse in terms of their nucleotide and oligonucleotide composition as well as presence of various sequence features that can affect physical properties of the DNA molecule. We present a survey of local sequence patterns which have a potential to promote non-canonical DNA conformations (i.e. different from standard B-DNA double helix) and interpret the results in terms of relationships with organisms' habitats, phylogenetic classifications, and other characteristics. Our present work differs from earlier similar surveys not only by investigating a wider range of sequence patterns in a large number of genomes but also by using a more realistic null model to assess significant deviations. Our results show that simple sequence repeats and Z-DNA-promoting patterns are generally suppressed in prokaryotic genomes, whereas palindromes and inverted repeats are over-represented. Representation of patterns that promote Z-DNA and intrinsic DNA curvature increases with increasing optimal growth temperature (OGT), and decreases with increasing oxygen requirement. Additionally, representations of close direct repeats, palindromes and inverted repeats exhibit clear negative trends with increasing OGT. The observed relationships with environmental characteristics, particularly OGT, suggest possible evolutionary scenarios of structural adaptation of DNA to particular environmental niches.

Comparative analysis to identify determinants of changing life style in Thermosynechococcus elongatus BP-1, a thermophilic cyanobacterium

A comparative genomics analysis among all forty whole genome sequences available for cyanobacteria (3 thermophiles- Thermosynechococcus elongatus BP-1, Synechococcus sp. JA-2-3B'a (2-13), Synechococcus sp. JA-3-3Ab and 37 mesophiles) was performed to identify genomic and proteomic factors responsible for the behaviour of T. elongatus BP-1, a thermophilic unicellular cyanobacterium with optimum growth temperature [OGT] of 55°C. Majority of genomic and proteomic characteristics for this cyanobacterium indicated contrasting features indicating its mesophilic behaviour while the role of mutational biasness and selection pressure is thought to be responsible for high OGT. Contradictory results were obtained for T. elongatus for synonymous codon usage, CvP-bias and amino acid composition with respect to thermophilic behaviour. Calculated J2 index is lowest among all cyanobacterial genomes. Except for proline and termination codons, T. elongatus showed synonymous codon usage pattern which is expected for mesophiles. Results indicated that among cyanobacterial genomes, majority of genomic and proteomic determinants put T. elongatus very close to mesophiles and the whole genome of this organism represents continuous gain of mesophilic rather than thermophilic behavior.

The temperature dependent proteomic analysis of Thermotoga maritima

Thermotoga maritima (T. maritima) is a typical thermophile, and its proteome response to environmental temperature changes has yet to be explored. This study aims to uncover the temperature-dependent proteins of T. maritima using comparative proteomic approach. T. maritima was cultured under four temperatures, 60°C, 70°C, 80°C and 90°C, and the bacterial proteins were extracted and electrophoresed in two-dimensional mode. After analysis of gel images, a total of 224 spots, either cytoplasm or membrane, were defined as temperature-dependent. Of these spots, 75 unique bacterial proteins were identified using MALDI TOF/TOF MS. As is well known, the chaperone proteins such as heat shock protein 60 and elongation factor Tu, were up-regulated in abundance due to increased temperature. However, several temperature-dependent proteins of T. maritima responded very differently when compared to responses of the thermophile T. tengcongensis. Intriguingly, a number of proteins involved in central carbohydrate metabolism were significantly up-regulated at higher temperature. Their corresponding mRNA levels were elevated accordingly. The increase in abundance of several key enzymes indicates that a number of central carbohydrate metabolism pathways of T. maritima are activated at higher temperatures.

Effect of codon message on xylanase thermal activity

Because the genetic codon is known for degeneracy, its effect on enzyme thermal property is seldom investigated. A dataset was constructed for GH10 xylanase coding sequences and optimal temperatures for activity (T(opt)). Codon contents and relative synonymous codon usages were calculated and respectively correlated with the enzyme T(opt) values, which were used to describe the xylanase thermophilic tendencies without dividing them into two thermophilic and mesophilic groups. After analyses of codon content and relative synonymous codon usages were checked by the Bonferroni correction, we found five codons, with three (AUA, AGA, and AGG) correlating positively and two (CGU and AGC) correlating negatively with the T(opt) value. The three positive codons are purine-rich codons, and the two negative codons have A-ends. The two negative codons are pyridine-rich codons, and one has a C-end. Comparable with the codon C- and A-ending features, C- and A-content within mRNA correlated negatively and positively with the T(opt) value, respectively. Thereby, codons have effects on enzyme thermal property. When the issue is analyzed at the residual level, the effect of codon message is lost. The codons relating to enzyme thermal property are selected by thermophilic force at nucleotide level.

Microbial lifestyle and genome signatures

Microbes are known for their unique ability to adapt to varying lifestyle and environment, even to the extreme or adverse ones. The genomic architecture of a microbe may bear the signatures not only of its phylogenetic position, but also of the kind of lifestyle to which it is adapted. The present review aims to provide an account of the specific genome signatures observed in microbes acclimatized to distinct lifestyles or ecological niches. Niche-specific signatures identified at different levels of microbial genome organization like base composition, GC-skew, purine-pyrimidine ratio, dinucleotide abundance, codon bias, oligonucleotide composition etc. have been discussed. Among the specific cases highlighted in the review are the phenomena of genome shrinkage in obligatory host-restricted microbes, genome expansion in strictly intra-amoebal pathogens, strand-specific codon usage in intracellular species, acquisition of genome islands in pathogenic or symbiotic organisms, discriminatory genomic traits of marine microbes with distinct trophic strategies, and conspicuous sequence features of certain extremophiles like those adapted to high temperature or high salinity.

Proteome adaptation to high temperatures in the ectothermic hydrothermal vent Pompeii worm

Taking advantage of the massive genome sequencing effort made on thermophilic prokaryotes, thermal adaptation has been extensively studied by analysing amino acid replacements and codon usage in these unicellular organisms. In most cases, adaptation to thermophily is associated with greater residue hydrophobicity and more charged residues. Both of these characteristics are positively correlated with the optimal growth temperature of prokaryotes. In contrast, little information has been collected on the molecular 'adaptive' strategy of thermophilic eukaryotes. The Pompeii worm A. pompejana, whose transcriptome has recently been sequenced, is currently considered as the most thermotolerant eukaryote on Earth, withstanding the greatest thermal and chemical ranges known. We investigated the amino-acid composition bias of ribosomal proteins in the Pompeii worm when compared to other lophotrochozoans and checked for putative adaptive changes during the course of evolution using codon-based Maximum likelihood analyses. We then provided a comparative analysis of codon usage and amino-acid replacements from a greater set of orthologous genes between the Pompeii worm and Paralvinella grasslei, one of its closest relatives living in a much cooler habitat. Analyses reveal that both species display the same high GC-biased codon usage and amino-acid patterns favoring both positively-charged residues and protein hydrophobicity. These patterns may be indicative of an ancestral adaptation to the deep sea and/or thermophily. In addition, the Pompeii worm displays a set of amino-acid change patterns that may explain its greater thermotolerance, with a significant increase in Tyr, Lys and Ala against Val, Met and Gly. Present results indicate that, together with a high content in charged residues, greater proportion of smaller aliphatic residues, and especially alanine, may be a different path for metazoans to face relatively 'high' temperatures and thus a novelty in thermophilic metazoans.

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.

POSITIVE DARWINIAN SELECTION: DOES THE COMPARATIVE METHOD RULE?

To detect positive Darwinian selection it is thought essential to compare two sequences. Despite its defects, "the comparative method rules." However, genes evolving rapidly under positive selection conflict more with internal forces (the genome phenotype) than genes evolving slowly under negative selection. In particular, there is conflict with stem-loop potential. The conflict between protein-encoding potential (primary information) and stem-loop potential (secondary information) permits detection of positive selection in a single sequence. The degree to which secondary information is compromised provides a measure of the speed of transmission of primary information. Thus, the sovereignty of the comparative method is challenged not only by its own defects, but also by the availability of a single-sequence method. However, while of limited utility for positive selection, the comparative method casts new light on Darwin's great question — the origin of species. Comparison of rates of synonymous and non-synonymous mutation suggests that branching into new species begins with synonymous mutations.

Differences in dinucleotide frequencies of thermophilic genes encoding water soluble and membrane proteins

The occurrence frequencies of the dinucleotides of genes of three thermophilic and three mesophilic species from both archaea and eubacteria were investigated in this study. The genes encoding water soluble proteins were rich in the dinucleotides of purine dimers, whereas the genes encoding membrane proteins were rich in pyrimidine dimers. The dinucleotides of purine dimers are the counterparts of pyrimidine dimers in a double-stranded DNA. The purine/pyrimidine dimers were favored in the thermophiles but not in the mesophiles, based on comparisons of observed and expected frequencies. This finding is in agreement with our previous study which showed that purine/pyrimidine dimers are positive factors that increase the thermal stability of DNA. The dinucleotides AA, AG, and GA are components of the codons of charged residues of Glu, Asp, Lys, and Arg, and the dinucleotides TT, CT, and TC are components of the codons of hydrophobic residues of Leu, Ile, and Phe. This is consistent with the suitabilities of the different amino acid residues for water soluble and membrane proteins. Our analysis provides a picture of how thermophilic species produce water soluble and membrane proteins with distinctive characters: the genes encoding water soluble proteins use DNA sequences rich in purine dimers, and the genes encoding membrane proteins use DNA sequences rich in pyrimidine dimers on the opposite strand.

alpha-Amylase: an ideal representative of thermostable enzymes

The conditions prevailing in the industrial applications in which enzymes are used are rather extreme, especially with respect to temperature and pH. Therefore, there is a continuing demand to improve the stability of enzymes and to meet the requirements set by specific applications. In this respect, thermostable enzymes have been proposed to be industrially relevant. In this review, alpha-amylase, a well-established representative of thermostable enzymes, providing an attractive model for the investigation of the structural basis of thermostability of proteins, has been discussed.

Stretches of alternating pyrimidine/purines and purines are respectively linked with pathogenicity and growth temperature in prokaryotes

Background

The genomic fractions of purine (RR) and alternating pyrimidine/purine (YR) stretches of 10 base pairs or more, have been linked to genomic AT content, the formation of different DNA helices, strand-biased gene distribution, DNA structure, and more. Although some of these factors are a consequence of the chemical properties of purines and pyrimidines, a thorough statistical examination of the distributions of YR/RR stretches in sequenced prokaryotic chromosomes has to the best of our knowledge, not been undertaken. The aim of this study is to expand upon previous research by using regression analysis to investigate how AT content, habitat, growth temperature, pathogenicity, phyla, oxygen requirement and halotolerance correlated with the distribution of RR and YR stretches in prokaryotes.

Results

Our results indicate that RR and YR-stretches are differently distributed in prokaryotic phyla. RR stretches are overrepresented in all phyla except for the Actinobacteria and β-Proteobacteria. In contrast, YR tracts are underrepresented in all phyla except for the β-Proteobacterial group. YR-stretches are associated with phylum, pathogenicity and habitat, whilst RR-tracts are associated with phylum, AT content, oxygen requirement, growth temperature and halotolerance. All associations described were statistically significant with p < 0.001.

Conclusion

Analysis of chromosomal distributions of RR/YR sequences in prokaryotes reveals a set of associations with environmental factors not observed with mono- and oligonucleotide frequencies. This implies that important information can be found in the distribution of RR/YR stretches that is more difficult to obtain from genomic mono- and oligonucleotide frequencies. The association between pathogenicity and fractions of YR stretches is assumed to be linked to recombination and horizontal transfer.

Darwinian selection for sites of Asn-linked glycosylation in phylogenetically disparate eukaryotes and viruses

Numerous protists and rare fungi have truncated Asn-linked glycan precursors and lack N-glycan-dependent quality control (QC) systems for glycoprotein folding in the endoplasmic reticulum. Here, we show that the abundance of sequons (NXT or NXS), which are sites for N-glycosylation of secreted and membrane proteins, varies by more than a factor of 4 among phylogenetically diverse eukaryotes, based on a few variables. There is positive correlation between the density of sequons and the AT content of coding regions, although no causality can be inferred. In contrast, there appears to be Darwinian selection for sequons containing Thr, but not Ser, in eukaryotes that have N-glycan-dependent QC systems. Selection for sequons with Thr, which nearly doubles the sequon density in human secreted and membrane proteins, occurs by an increased conditional probability that Asn and Thr are present in sequons rather than elsewhere. Increasing sequon densities of the hemagglutinin (HA) of influenza viruses A/H3N2 and A/H1N1 during the past few decades of human infection also result from an increased conditional probability that Asn, Thr, and Ser are present in sequons rather than elsewhere. In contrast, there is no selection on sequons by this mechanism in HA of A/H5N1 or 2009 A/H1N1 (Swine flu). Very strong selection for sequons with both Thr and Ser in glycoprotein of M(r) 120,000 (gp120) of HIV and related retroviruses results from this same mechanism, as well as amino acid composition bias and increases in AT content. We conclude that there is Darwinian selection for sequons in phylogenetically disparate eukaryotes and viruses.

Classification and regression tree (CART) analyses of genomic signatures reveal sets of tetramers that discriminate temperature optima of archaea and bacteria

Classification and regression tree (CART) analysis was applied to genome-wide tetranucleotide frequencies (genomic signatures) of 195 archaea and bacteria. Although genomic signatures have typically been used to classify evolutionary divergence, in this study, convergent evolution was the focus. Temperature optima for most of the organisms examined could be distinguished by CART analyses of tetranucleotide frequencies. This suggests that pervasive (nonlinear) qualities of genomes may reflect certain environmental conditions (such as temperature) in which those genomes evolved. The predominant use of GAGA and AGGA as the discriminating tetramers in CART models suggests that purine-loading and codon biases of thermophiles may explain some of the results.

Different clustering of genomes across life using the A-T-C-G and degenerate R-Y alphabets: early and late signaling on genome evolution?

In this study, we have calculated distances between genomes based on our previously developed compositional spectra (CS) analysis. The study was conducted using genomes of 39 species of Eukarya, Eubacteria, and Archaea. Based on CS distances, we produced two different consensus dendrograms for four- and two-letter (purine-pyrimidine) alphabets. A comparison of the obtained structure using purine-pyrimidine alphabet with the standard three-kingdom (3K) scheme reveals substantial similarity. Surprisingly, this is not the case when the same procedure is based on the four-letter alphabet. In this situation, we also found three main clusters-but different from those in the 3K scheme. In particular, one of the clusters includes Eukarya and thermophilic bacteria and a part of the considered Archaea species. We speculate that the key factor in the last classification (based on the A-T-G-C alphabet) is related to ecology: two ecological parameters, temperature and oxygen, distinctly explain the clustering revealed by compositional spectra in the four-letter alphabet. Therefore, we assume that this result reflects two interdependent processes: evolutionary divergence and superimposed ecological convergence of the genomes, albeit another process, horizontal transfer, cannot be excluded as an important contributing factor.

Mutations Stabilize Small Subunit Ribosomal RNA in Desiccation-Tolerant Cyanobacteria Nostoc

The ribosomal RNA molecule is an ideal model for evaluating the stability of a gene product under desiccation stress. We isolated 8 Nostoc strains that had the capacity to withstand desiccation in habitats and sequenced their 16S rRNA genes. The stabilities of 16S rRNAs secondary structures, indicated by free energy change of folding, were compared among Nostoc and other related species. The results suggested that 16S rRNA secondary structures of the desiccation-tolerant Nostoc strains were more stable than that of planktonic Nostocaceae species. The stabilizing mutations were divided into two categories: (1) those causing GC to replace other types of base pairs in stems and (2) those causing extension of stems. By mapping stabilizing mutations onto the Nostoc phylogenetic tree based on 16S rRNA gene, it was shown that most of stabilizing mutations had evolved during adaptive radiation among Nostoc spp. The evolution of 16S rRNA along the Nostoc lineage is suggested to be selectively advantageous under desiccation stress.

Looking for organization patterns of highly expressed genes: purine-pyrimidine composition of precursor mRNAs

We analyzed precursor messenger RNAs (pre-mRNAs) of 12 eukaryotic species. In each species, three groups of highly expressed genes, ribosomal proteins, heat shock proteins, and amino-acyl tRNA synthetases, were compared with a control group (randomly selected genes). The purine-pyrimidine (R-Y) composition of pre-mRNAs of the three targeted gene groups proved to differ significantly from the control. The exons of the three groups tested have higher purine contents and R-tract abundance and lower abundance of Y-tracts compared to the control (R-tract-tract of sequential purines with Rn>or=5; Y-tract-tract of sequential pyrimidines with Yn>or=5). In species widely employing "intron definition" in the splicing process, the Y content of introns of the three targeted groups appeared to be higher compared to the control group. Furthermore, in all examined species, the introns of the targeted genes have a lower abundance of R-tracts compared to the control. We hypothesized that the R-Y composition of the targeted gene groups contributes to high rate and efficiency of both splicing and translation, in addition to the mRNA coding role. This is presumably achieved by (1) reducing the possibility of the formation of secondary structures in the mRNA, (2) using the R-tracts and R-biased sequences as exonic splicing enhancers, (3) lowering the amount of targets for pyrimidine tract binding protein in the exons, and (4) reducing the amount of target sequences for binding of serine/arginine-rich (SR) proteins in the introns, thereby allowing SR proteins to bind to proper (exonic) targets.

Prokaryotes that grow optimally in acid have purine-poor codons in long open reading frames

In nucleic acids the N-glycosyl bonds between purines and their ribose sugar moities are broken under acid conditions. If one strand of a duplex DNA segment were more vulnerable to mutation than the other, then the archaeon Picrophilus torridus, with an optimum growth pH near zero, could have adapted by decreasing the purine content of that strand. Yet, P. torridus has an optimum growth temperature near 60 degrees C, and thermophiles prefer purine-rich codons. We found that, as in other thermophiles, high growth temperature correlates with the use of purine-rich codons. The extra purines are often in third, non-amino acid determining, codon positions. However, as in other acidophiles, as open reading frame lengths increase, there is increased use of purine-poor codons, particularly those without purines in second, amino acid-determining, codon positions. Thus, P. torridus can be seen as adapting (a) to temperature by increasing its purines in all open reading frames without greatly impacting protein amino acid compositions, and (b) to pH by decreasing purines in longer open reading frames, thereby potentially impacting protein amino acid compositions. It is proposed that longer open reading frames, being larger mutational targets, have become less vulnerable to depurination by virtue of pyrimidine for purine substitutions.

Synonymous codon usage and its potential link with optimal growth temperature in prokaryotes

The relationship between codon usage in prokaryotes and their ability to grow at extreme temperatures has been given much attention over the past years. Previous studies have suggested that the difference in synonymous codon usage between (hyper)thermophiles and mesophiles is a consequence of a selective pressure linked to growth temperature. Here, we performed an updated analysis of the variation in synonymous codon usage with growth temperature; our study includes a large number of species from a wide taxonomic and growth temperature range. The presence of psychrophilic species in our study allowed us to test whether the same selective pressure acts on synonymous codon usage at very low growth temperature. Our results show that the synonymous codon usage for Arg (through the AGG, AGA and CGT codons) is the most discriminating factor between (hyper)thermophilic and non-thermophilic species, thus confirming previous studies. We report the unusual clustering of an Archaeal psychrophile with the thermophilic and hyperthermophilic species on the synonymous codon usage factorial map; the other psychrophiles in our study cluster with the mesophilic species. Our conclusion is that the difference in synonymous codon usage between (hyper)thermophilic and non-thermophilic species cannot be clearly attributed to a selective pressure linked to growth at high temperatures.

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.

Codon Signature Extremes In Eukaryote genomes

Twenty-one complete eukaryotic genomes are compared for codon signature biases. The codon signature refers to the dinucleotide relative abundance values at codon sites {1, 2}, {2, 3}, and {3, 4} (4 = 1 of the next codon site). The genomes under study include human, mouse, chicken, three invertebrates, one plant species, eight fungi, and six protists. The dinucleotide CpG is significantly underrepresented at all contiguous codon sites and drastically suppressed in noncoding regions in mammalian species, in yeast-like genomes, in the dicotArabidopsis thaliana, but not in the filamentous fungiNeurospora crassaandAsperigillus fumigatus, and in the protistEntamoeba histolytica.The dinucleotide TpA, probably due to DNA structural weaknesses, is underrepresented genome-wide and significantly underrepresented in the codon signature for all contiguous codon sites in mammals, inverterbrates, plants, and fungi, but somewhat restricted to codon sites {1, 2} among protists helping in avoidance of stop codons. The amino acid Ser, not of abundance in bacterial genomes, generally ranks among the two most used amino acids among eukaryotes ostensibly resulting from greater activity in the nucleus. The observed differences are linked to specifics of methylation, context-dependent mutation, DNA repair, and replication. For example, the amino acid Leu is broadly abundant in all life domains generally resulting from extra occurrences of the codon TTR, R purine. The malarial protistPlasmodium falciparumshows many codon signature extremes.

On the correlation between genomic G+C content and optimal growth temperature in prokaryotes: data quality and confounding factors

The correlation between genomic G+C content and optimal growth temperature in prokaryotes has gained renewed interest after Musto et al. [H. Musto, H. Naya, A. Zavala, H. Romero, F. Alvarex-Valin, G. Bernardi, Correlations between genomic GC levels and optimal growth temperatures in prokaryotes, FEBS Lett. 573 (2004) 73-77], reported that positive correlations exist in 15 families studied. We have reanalyzed their data and found that when genome size and data quality were adjusted for, there was no significant evidence of relationship between optimal temperature and GC content for two of the families that had previously shown strongly significant correlations. Using updated temperature optima for Halobacteriaceae species we found the correlation is insignificant in this family. For the family Enterobacteriaceae when genome size and optimal temperature are included in a multiple linear regression, only genome size is significant as a predictor of GC content. We showed that more profound statistical methods than simple two factor correlation analysis should be used for analyzing complex intrinsic and extrinsic factors that affect genomic GC content. We further found that a positive correlation between temperature and genomic GC is only evident in free-living species of low optimal growth temperatures.

Environments shape the nucleotide composition of genomes

To test the impact of environments on genome evolution, we analysed the relative abundance of the nucleotides guanine and cytosine ('GC content') of large numbers of sequences from four distinct environmental samples (ocean surface water, farm soil, an acidophilic mine drainage biofilm and deep-sea whale carcasses). We show that the GC content of complex microbial communities seems to be globally and actively influenced by the environment. The observed nucleotide compositions cannot be easily explained by distinct phylogenetic origins of the species in the environments; the genomic GC content may change faster than was previously thought, and is also reflected in the amino-acid composition of the proteins in these habitats.

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).