Protein evolution in deep sea bacteria: an analysis of amino acids substitution rates

The systematic codon usage bias has an important effect on genetic adaption in native species

Random mutations increase genetic variety and natural selection enhances adaption over generations. Codon usage biases (CUB) provide clues about the genome adaptation mechanisms of native species and extremophile species. Significant numbers of gene (CDS) of nine classes of endangered, native species, including extremophiles and mesophiles were utilised to compute CUB. Codon usage patterns differ among the lineages of endangered and extremophiles with native species. Polymorphic usage of nucleotides with codon burial suggests parallelism of native species within relatively confined taxonomic groups. Utilizing the deviation pattern of CUB of endangered and native species, I present a calculation parameter to estimate the extinction risk of endangered species. Species diversity and extinction risk are both positively associated with the propensity of random mutation in CDS (Coding DNA sequence). Codon bias tenet profoundly selected and it governs to adaptive evolution of native species.

Sulfate-dependant microbially induced corrosion of mild steel in the deep sea: a 10-year microbiome study

BACKGROUND

Metal corrosion in seawater has been extensively studied in surface and shallow waters. However, infrastructure is increasingly being installed in deep-sea environments, where extremes of temperature, salinity, and high hydrostatic pressure increase the costs and logistical challenges associated with monitoring corrosion. Moreover, there is currently only a rudimentary understanding of the role of microbially induced corrosion, which has rarely been studied in the deep-sea. We report here an integrative study of the biofilms growing on the surface of corroding mooring chain links that had been deployed for 10 years at ~2 km depth and developed a model of microbially induced corrosion based on flux-balance analysis.

METHODS

We used optical emission spectrometry to analyze the chemical composition of the mooring chain and energy-dispersive X-ray spectrometry coupled with scanning electron microscopy to identify corrosion products and ultrastructural features. The taxonomic structure of the microbiome was determined using shotgun metagenomics and was confirmed by 16S amplicon analysis and quantitative PCR of the dsrB gene. The functional capacity was further analyzed by generating binned, genomic assemblies and performing flux-balance analysis on the metabolism of the dominant taxa.

RESULTS

The surface of the chain links showed intensive and localized corrosion with structural features typical of microbially induced corrosion. The microbiome on the links differed considerably from that of the surrounding sediment, suggesting selection for specific metal-corroding biofilms dominated by sulfur-cycling bacteria. The core metabolism of the microbiome was reconstructed to generate a mechanistic model that combines biotic and abiotic corrosion. Based on this metabolic model, we propose that sulfate reduction and sulfur disproportionation might play key roles in deep-sea corrosion.

CONCLUSIONS

The corrosion rate observed was higher than what could be expected from abiotic corrosion mechanisms under these environmental conditions. High corrosion rate and the form of corrosion (deep pitting) suggest that the corrosion of the chain links was driven by both abiotic and biotic processes. We posit that the corrosion is driven by deep-sea sulfur-cycling microorganisms which may gain energy by accelerating the reaction between metallic iron and elemental sulfur. The results of this field study provide important new insights on the ecophysiology of the corrosion process in the deep sea.

Biotechnological and Ecological Potential of Micromonospora provocatoris sp. nov., a Gifted Strain Isolated from the Challenger Deep of the Mariana Trench

A Micromonospora strain, isolate MT25^T, was recovered from a sediment collected from the Challenger Deep of the Mariana Trench using a selective isolation procedure. The isolate produced two major metabolites, n-acetylglutaminyl glutamine amide and desferrioxamine B, the chemical structures of which were determined using 1D and 2D-NMR, including ¹H-¹⁵N HSQC and ¹H-¹⁵N HMBC 2D-NMR, as well as high resolution MS. A whole genome sequence of the strain showed the presence of ten natural product-biosynthetic gene clusters, including one responsible for the biosynthesis of desferrioxamine B. Whilst 16S rRNA gene sequence analyses showed that the isolate was most closely related to the type strain of Micromonospora chalcea, a whole genome sequence analysis revealed it to be most closely related to Micromonospora tulbaghiae 45142^T. The two strains were distinguished using a combination of genomic and phenotypic features. Based on these data, it is proposed that strain MT25^T (NCIMB 15245^T, TISTR 2834^T) be classified as Micromonospora provocatoris sp. nov. Analysis of the genome sequence of strain MT25^T (genome size 6.1 Mbp) revealed genes predicted to responsible for its adaptation to extreme environmental conditions that prevail in deep-sea sediments.

The Molecular Basis for Life in Extreme Environments

Sampling and genomic efforts over the past decade have revealed an enormous quantity and diversity of life in Earth's extreme environments. This new knowledge of life on Earth poses the challenge of understandingits molecular basis in such inhospitable conditions, given that such conditions lead to loss of structure and of function in biomolecules from mesophiles. In this review, we discuss the physicochemical properties of extreme environments. We present the state of recent progress in extreme environmental genomics. We then present an overview of our current understanding of the biomolecular adaptation to extreme conditions. As our current and future understanding of biomolecular structure-function relationships in extremophiles requires methodologies adapted to extremes of pressure, temperature, and chemical composition, advances in instrumentation for probing biophysical properties under extreme conditions are presented. Finally, we briefly discuss possible future directions in extreme biophysics.

Whole Genome Sequence of Dermacoccus abyssi MT1.1 Isolated from the Challenger Deep of the Mariana Trench Reveals Phenazine Biosynthesis Locus and Environmental Adaptation Factors

Dermacoccus abyssi strain MT1.1^T is a piezotolerant actinobacterium that was isolated from Mariana Trench sediment collected at a depth of 10898 m. The organism was found to produce ten dermacozines (A‒J) that belonged to a new phenazine family and which displayed various biological activities such as radical scavenging and cytotoxicity. Here, we report on the isolation and identification of a new dermacozine compound, dermacozine M, the chemical structure of which was determined using 1D and 2D-NMR, and high resolution MS. A whole genome sequence of the strain contained six secondary metabolite-biosynthetic gene clusters (BGCs), including one responsible for the biosynthesis of a family of phenazine compounds. A pathway leading to the biosynthesis of dermacozines is proposed. Bioinformatic analyses of key stress-related genes provide an insight into how the organism adapted to the environmental conditions that prevail in the deep-sea.

Protein Stability in Titan's Subsurface Water Ocean

Models of Titan predict that there is a subsurface ocean of water and ammonia under a layer of ice. Such an ocean would be important in the search for extraterrestrial life since it provides a potentially habitable environment. To evaluate how Earth-based proteins would behave in Titan's subsurface ocean environment, we used molecular dynamics simulations to calculate the properties of proteins with the most common secondary structure types (alpha helix and beta sheet) in both Earth and Titan-like conditions. The Titan environment was simulated by using a temperature of 300 K, a pressure of 1000 bar, and a eutectic mixture of water and ammonia. We analyzed protein compactness, flexibility, and backbone dihedral distributions to identify differences between the two environments. Secondary structures in the Titan environment were found to be less long-lasting, less flexible, and had small differences in backbone dihedral preferences (e.g., in one instance a pi helix formed). These environment-driven differences could lead to changes in how these proteins interact with other biomolecules and therefore changes in how evolution would potentially shape proteins to function in subsurface ocean environments.

Deciphering the rationale behind specific codon usage pattern in extremophiles

Protein stability is affected at different hierarchies – gene, RNA, amino acid sequence and structure. Gene is the first level which contributes via varying codon compositions. Codon selectivity of an organism differs with normal and extremophilic milieu. The present work attempts at detailing the codon usage pattern of six extremophilic classes and their harmony. Homologous gene datasets of thermophile-mesophile, psychrophile-mesophile, thermophile-psychrophile, acidophile-alkaliphile, halophile-nonhalophile and barophile-nonbarophile were analysed for filtering statistically significant attributes. Relative abundance analysis, 1–9 scale ranking, nucleotide compositions, attribute weighting and machine learning algorithms were employed to arrive at findings. AGG in thermophiles and barophiles, CAA in mesophiles and psychrophiles, TGG in acidophiles, GAG in alkaliphiles and GAC in halophiles had highest preference. Preference of GC-rich and G/C-ending codons were observed in halophiles and barophiles whereas, a decreasing trend was reflected in psychrophiles and alkaliphiles. GC-rich codons were found to decrease and G/C-ending codons increased in thermophiles whereas, acidophiles showed equal contents of GC-rich and G/C-ending codons. Codon usage patterns exhibited harmony among different extremophiles and has been detailed. However, the codon attribute preferences and their selectivity of extremophiles varied in comparison to non-extremophiles. The finding can be instrumental in codon optimization application for heterologous expression of extremophilic proteins.

Microbial Community Structure-Function Relationships in Yaquina Bay Estuary Reveal Spatially Distinct Carbon and Nitrogen Cycling Capacities

Linking microbial community structure to ecological processes requires understanding of the functional roles among individual populations and the factors that influence their distributions. These structure–function relationships are particularly difficult to disentangle in estuaries, due to highly variable physico-chemical conditions. Yet, examining microbe-mediated turnover of resources in these “bioreactor” ecosystems is critical for understanding estuarine ecology. In this study, a combined metagenomics and metaproteomics approach was used to show that the unequal distribution of microbial populations across the Yaquina Bay estuary led to a habitat-specific taxonomic and functional structure and a clear spatial distribution in microbe-mediated capacities for cycling of carbon and nitrogen. For example, size-fractionation revealed that communities inhabiting suspended particulate material encoded more diverse types of metabolisms (e.g., fermentation and denitrification) than those with a planktonic lifestyle, suggesting that the metabolic reactions can differ between size fractions of the same parcel of an estuarine water column. Similarly, communities inhabiting oligotrophic conditions in the lower estuary were enriched in genes involved in central carbon metabolism (e.g., TCA cycle), while communities in the upper estuary were enriched in genes typical of copiotrophic populations (e.g., cell growth, cell division). Integrating gene and protein data revealed that abundant populations of Flavobacteriales and Rhodobacterales encoded similar genomic functions, yet differed significantly in protein expression, dedicating a large proportion of their respective proteomes to rapid growth and division versus metabolic versatility and resource acquisition. This suggested potentially distinct life-strategies between these two co-occurring lineages and was concomitant with differing patterns of positive evolutionary selection on their encoded genes. Microbial communities and their functions across Yaquina Bay appear to be structured by population-level habitat preferences, resulting in spatially distinct elemental cycling, while within each community, forces such as competitive exclusion and evolutionary selection influence species life-strategies and may help maintain microbial diversity.

Insights into the molecular basis of piezophilic adaptation: Extraction of piezophilic signatures

Piezophiles are the organisms which can successfully survive at extreme pressure conditions. However, the molecular basis of piezophilic adaptation is still poorly understood. Analysis of the protein sequence adjustments that had taken place during evolution can help to reveal the sequence adaptation parameters responsible for protein functional and structural adaptation at such high pressure conditions. In this current work we have used SVM classifier for filtering strong instances and generated human interpretable rules from these strong instances by using the PART algorithm. These generated rules were analyzed for getting insights into the molecular signature patterns present in the piezophilic proteins. The experiments were performed on three different temperature ranges piezophilic groups, namely psychrophilic-piezophilic, mesophilic-piezophilic, and thermophilic-piezophilic for the detailed comparative study. The best classification results were obtained as we move up the temperature range from psychrophilic-piezophilic to thermophilic-piezophilic. Based on the physicochemical classification of amino acids and using feature ranking algorithms, hydrophilic and polar amino acid groups have higher discriminative ability for psychrophilic-piezophilic and mesophilic-piezophilic groups along with hydrophobic and nonpolar amino acids for the thermophilic-piezophilic groups. We also observed an overrepresentation of polar, hydrophilic and small amino acid groups in the discriminatory rules of all the three temperature range piezophiles along with aliphatic, nonpolar and hydrophobic groups in the mesophilic-piezophilic and thermophilic-piezophilic groups.

Molecular dynamics simulations of the Nip7 proteins from the marine deep- and shallow-water Pyrococcus species

Background

The identification of the mechanisms of adaptation of protein structures to extreme environmental conditions is a challenging task of structural biology. We performed molecular dynamics (MD) simulations of the Nip7 protein involved in RNA processing from the shallow-water (P. furiosus) and the deep-water (P. abyssi) marine hyperthermophylic archaea at different temperatures (300 and 373 K) and pressures (0.1, 50 and 100 MPa). The aim was to disclose similarities and differences between the deep- and shallow-sea protein models at different temperatures and pressures.

Results

The current results demonstrate that the 3D models of the two proteins at all the examined values of pressures and temperatures are compact, stable and similar to the known crystal structure of the P. abyssi Nip7. The structural deviations and fluctuations in the polypeptide chain during the MD simulations were the most pronounced in the loop regions, their magnitude being larger for the C-terminal domain in both proteins. A number of highly mobile segments the protein globule presumably involved in protein-protein interactions were identified. Regions of the polypeptide chain with significant difference in conformational dynamics between the deep- and shallow-water proteins were identified.

Conclusions

The results of our analysis demonstrated that in the examined ranges of temperatures and pressures, increase in temperature has a stronger effect on change in the dynamic properties of the protein globule than the increase in pressure. The conformational changes of both the deep- and shallow-sea protein models under increasing temperature and pressure are non-uniform. Our current results indicate that amino acid substitutions between shallow- and deep-water proteins only slightly affect overall stability of two proteins. Rather, they may affect the interactions of the Nip7 protein with its protein or RNA partners.

Transcriptomics reveal several gene expression patterns in the piezophile Desulfovibrio hydrothermalis in response to hydrostatic pressure

RNA-seq was used to study the response of Desulfovibrio hydrothermalis, isolated from a deep-sea hydrothermal chimney on the East-Pacific Rise at a depth of 2,600 m, to various hydrostatic pressure growth conditions. The transcriptomic datasets obtained after growth at 26, 10 and 0.1 MPa identified only 65 differentially expressed genes that were distributed among four main categories: aromatic amino acid and glutamate metabolisms, energy metabolism, signal transduction, and unknown function. The gene expression patterns suggest that D. hydrothermalis uses at least three different adaptation mechanisms, according to a hydrostatic pressure threshold (HP_t) that was estimated to be above 10 MPa. Both glutamate and energy metabolism were found to play crucial roles in these mechanisms. Quantitation of the glutamate levels in cells revealed its accumulation at high hydrostatic pressure, suggesting its role as a piezolyte. ATP measurements showed that the energy metabolism of this bacterium is optimized for deep-sea life conditions. This study provides new insights into the molecular mechanisms linked to hydrostatic pressure adaptation in sulfate-reducing bacteria.

Comparative analysis of deep-sea bacterioplankton OMICS revealed the occurrence of habitat-specific genomic attributes

Bathyal aphotic ocean represents the largest biotope on our planet, which sustains highly diverse but low-density microbial communities, with yet untapped genomic attributes, potentially useful for discovery of new biomolecules, industrial enzymes and pathways. In the last two decades, culture-independent approaches of high-throughput sequencing have provided new insights into structure and function of marine bacterioplankton, leading to unprecedented opportunities to accurately characterize microbial communities and their interactions with the environments. In the present review we focused on the analysis of relatively few deep-sea OMICS studies, completed thus far, to find the specific genomic patterns determining the lifeway and adaptation mechanisms of prokaryotes thriving in the dark deep ocean below the depth of 1000m. Phylogenomic and omic studies provided clear evidence that the bathyal microbial communities are distinct from the epipelagic counterparts and, along with generally larger genomes, possess their own habitat-specific genomic attributes. The high abundance in the deep ocean OMICS of the systems for environmental sensing, signal transduction and metabolic versatility as compared to the epipelagic counterparts is thought to enable the deep-sea bacterioplankton to rapidly adapt to changing environmental conditions associated with resource scarcity and high diversity of energy and carbon substrates in the bathyal biotopes. Together with a versatile heterotrophy, mixotrophy and anaplerosis are thought to enable the deep-sea bacterioplankton to cope with these environmental conditions.

Exposure of Bacillus subtilis to low pressure (5 kilopascals) induces several global regulons, including those involved in the SigB-mediated general stress response

Studies of how microorganisms respond to pressure have been limited mostly to the extreme high pressures of the deep sea (i.e., the piezosphere). In contrast, despite the fact that the growth of most bacteria is inhibited at pressures below ∼2.5 kPa, little is known of microbial responses to low pressure (LP). To study the global LP response, we performed transcription microarrays on Bacillus subtilis cells grown under normal atmospheric pressure (∼101 kPa) and a nearly inhibitory LP (5 kPa), equivalent to the pressure found at an altitude of ∼20 km. Microarray analysis revealed altered levels of 363 transcripts belonging to several global regulons (AbrB, CcpA, CodY, Fur, IolR, ResD, Rok, SigH, Spo0A). Notably, the highest number of upregulated genes, 86, belonged to the SigB-mediated general stress response (GSR) regulon. Upregulation of the GSR by LP was confirmed by monitoring the expression of the SigB-dependent ctc-lacZ reporter fusion. Measuring transcriptome changes resulting from exposure of bacterial cells to LP reveals insights into cellular processes that may respond to LP exposure.

Explaining bathymetric diversity patterns in marine benthic invertebrates and demersal fishes: physiological contributions to adaptation of life at depth

Bathymetric biodiversity patterns of marine benthic invertebrates and demersal fishes have been identified in the extant fauna of the deep continental margins. Depth zonation is widespread and evident through a transition between shelf and slope fauna from the shelf break to 1000 m, and a transition between slope and abyssal fauna from 2000 to 3000 m; these transitions are characterised by high species turnover. A unimodal pattern of diversity with depth peaks between 1000 and 3000 m, despite the relatively low area represented by these depths. Zonation is thought to result from the colonisation of the deep sea by shallow-water organisms following multiple mass extinction events throughout the Phanerozoic. The effects of low temperature and high pressure act across hierarchical levels of biological organisation and appear sufficient to limit the distributions of such shallow-water species. Hydrostatic pressures of bathyal depths have consistently been identified experimentally as the maximum tolerated by shallow-water and upper bathyal benthic invertebrates at in situ temperatures, and adaptation appears required for passage to deeper water in both benthic invertebrates and demersal fishes. Together, this suggests that a hyperbaric and thermal physiological bottleneck at bathyal depths contributes to bathymetric zonation. The peak of the unimodal diversity-depth pattern typically occurs at these depths even though the area represented by these depths is relatively low. Although it is recognised that, over long evolutionary time scales, shallow-water diversity patterns are driven by speciation, little consideration has been given to the potential implications for species distribution patterns with depth. Molecular and morphological evidence indicates that cool bathyal waters are the primary site of adaptive radiation in the deep sea, and we hypothesise that bathymetric variation in speciation rates could drive the unimodal diversity-depth pattern over time. Thermal effects on metabolic-rate-dependent mutation and on generation times have been proposed to drive differences in speciation rates, which result in modern latitudinal biodiversity patterns over time. Clearly, this thermal mechanism alone cannot explain bathymetric patterns since temperature generally decreases with depth. We hypothesise that demonstrated physiological effects of high hydrostatic pressure and low temperature at bathyal depths, acting on shallow-water taxa invading the deep sea, may invoke a stress-evolution mechanism by increasing mutagenic activity in germ cells, by inactivating canalisation during embryonic or larval development, by releasing hidden variation or mutagenic activity, or by activating or releasing transposable elements in larvae or adults. In this scenario, increased variation at a physiological bottleneck at bathyal depths results in elevated speciation rate. Adaptation that increases tolerance to high hydrostatic pressure and low temperature allows colonisation of abyssal depths and reduces the stress-evolution response, consequently returning speciation of deeper taxa to the background rate. Over time this mechanism could contribute to the unimodal diversity-depth pattern.

Single-cell enabled comparative genomics of a deep ocean SAR11 bathytype

Bacterioplankton of the SAR11 clade are the most abundant microorganisms in marine systems, usually representing 25% or more of the total bacterial cells in seawater worldwide. SAR11 is divided into subclades with distinct spatiotemporal distributions (ecotypes), some of which appear to be specific to deep water. Here we examine the genomic basis for deep ocean distribution of one SAR11 bathytype (depth-specific ecotype), subclade Ic. Four single-cell Ic genomes, with estimated completeness of 55%-86%, were isolated from 770 m at station ALOHA and compared with eight SAR11 surface genomes and metagenomic datasets. Subclade Ic genomes dominated metagenomic fragment recruitment below the euphotic zone. They had similar COG distributions, high local synteny and shared a large number (69%) of orthologous clusters with SAR11 surface genomes, yet were distinct at the 16S rRNA gene and amino-acid level, and formed a separate, monophyletic group in phylogenetic trees. Subclade Ic genomes were enriched in genes associated with membrane/cell wall/envelope biosynthesis and showed evidence of unique phage defenses. The majority of subclade Ic-specfic genes were hypothetical, and some were highly abundant in deep ocean metagenomic data, potentially masking mechanisms for niche differentiation. However, the evidence suggests these organisms have a similar metabolism to their surface counterparts, and that subclade Ic adaptations to the deep ocean do not involve large variations in gene content, but rather more subtle differences previously observed deep ocean genomic data, like preferential amino-acid substitutions, larger coding regions among SAR11 clade orthologs, larger intergenic regions and larger estimated average genome size.

The impact of genomic variability on gene expression in environmental Saccharomyces cerevisiae strains

Environmental Saccharomyces cerevisiae strains are crucially important, as they represent the large pool from which domesticated industrial yeasts have been selected, and vineyard strains can be considered the genetic reservoir from which industrial wine strains with strong fermentative behaviour are selected. Four vineyard strains with different fermentation performances were chosen from a large collection of strains isolated from Italian vineyards. Their genomes were sequenced to identify how genetic variations influence gene expression during fermentation and to clarify the evolutionary relationship between vineyard isolates and industrial wine strains. RNA sequencing was performed on the four vineyard strains, as well as on the industrial wine yeast strain EC1118 and on the laboratory strain S288c, at two stages of fermentation. We showed that there was a large gene cluster with variable promoter regions modifying gene expression in the strains. Our results indicate that it is the evolvability of the yeast promoter regions, rather than structural variations or strain-specific genes, that is the main cause of the differences in gene expression. This promoter variability, determined by variable tandem repeats and a high number of single-nucleotide polymorphisms together with 49 differentially expressed transcription factors, explained the strong phenotypic differences in the strains.

Rapid evolution of the sequences and gene repertoires of secreted proteins in bacteria

Proteins secreted to the extracellular environment or to the periphery of the cell envelope, the secretome, play essential roles in foraging, antagonistic and mutualistic interactions. We hypothesize that arms races, genetic conflicts and varying selective pressures should lead to the rapid change of sequences and gene repertoires of the secretome. The analysis of 42 bacterial pan-genomes shows that secreted, and especially extracellular proteins, are predominantly encoded in the accessory genome, i.e. among genes not ubiquitous within the clade. Genes encoding outer membrane proteins might engage more frequently in intra-chromosomal gene conversion because they are more often in multi-genic families. The gene sequences encoding the secretome evolve faster than the rest of the genome and in particular at non-synonymous positions. Cell wall proteins in Firmicutes evolve particularly fast when compared with outer membrane proteins of Proteobacteria. Virulence factors are over-represented in the secretome, notably in outer membrane proteins, but cell localization explains more of the variance in substitution rates and gene repertoires than sequence homology to known virulence factors. Accordingly, the repertoires and sequences of the genes encoding the secretome change fast in the clades of obligatory and facultative pathogens and also in the clades of mutualists and free-living bacteria. Our study shows that cell localization shapes genome evolution. In agreement with our hypothesis, the repertoires and the sequences of genes encoding secreted proteins evolve fast. The particularly rapid change of extracellular proteins suggests that these public goods are key players in bacterial adaptation.

Duplication and evolution of devA-like genes in Streptomyces has resulted in distinct developmental roles

Understanding morphological transformations is essential to elucidating the evolution and developmental biology of many organisms. The Gram-positive soil bacterium, Streptomyces coelicolor has a complex lifecycle which lends itself well to such studies. We recently identified a transcriptional regulator, devA, which is required for correct sporulation in this organism, with mutants forming short, mis-septate aerial hyphae. devA is highly conserved within the Streptomyces genus along with a duplicate copy, devE. Disruption of devE indicates this gene also plays a role in sporulation; however the phenotype of a devE mutant differs from a devA mutant, forming long un-septate aerial hyphae. Transcriptional analysis of devA and devE indicates that they are expressed at different stages of the lifecycle. This suggests that following duplication they have diverged in regulation and function. Analysis of fully sequenced actinomycete genomes shows that devA is found in a single copy in morphologically simpler actinobacteria, suggesting that duplication has lead to increased morphological complexity. Complementation studies with devA from Salinispora, which sporulates but does not form aerial hyphae, indicates the ancestral gene cannot complement devA or devE, suggesting neo-functionalisation has occurred. Analysis of the synonymous and non-synonymous nucleotide changes within the devA paralogues suggest subfunctionalisation has occurred as both copies have diverged from the ancestral sequences. Divergence is also asymmetric with a higher level of functional constraint observed in the DNA binding domain compared with the effector binding/oligomerisation domain, suggesting diversification in the substrate specificity of these paralogues has contributed to their evolution.

Evolutionary conservation of essential and highly expressed genes in Pseudomonas aeruginosa

Background

The constant increase in development and spread of bacterial resistance to antibiotics poses a serious threat to human health. New sequencing technologies are now on the horizon that will yield massive increases in our capacity for DNA sequencing and will revolutionize the drug discovery process. Since essential genes are promising novel antibiotic targets, the prediction of gene essentiality based on genomic information has become a major focus.

Results

In this study we demonstrate that pooled sequencing is applicable for the analysis of sequence variations of strain collections with more than 10 individual isolates. Pooled sequencing of 36 clinical Pseudomonas aeruginosa isolates revealed that essential and highly expressed proteins evolve at lower rates, whereas extracellular proteins evolve at higher rates. We furthermore refined the list of experimentally essential P. aeruginosa genes, and identified 980 genes that show no sequence variation at all. Among the conserved nonessential genes we found several that are involved in regulation, motility and virulence, indicating that they represent factors of evolutionary importance for the lifestyle of a successful environmental bacterium and opportunistic pathogen.

Conclusion

The detailed analysis of a comprehensive set of P. aeruginosa genomes in this study clearly disclosed detailed information of the genomic makeup and revealed a large set of highly conserved genes that play an important role for the lifestyle of this microorganism. Sequencing strain collections enables for a detailed and extensive identification of sequence variations as potential bacterial adaptation processes, e.g., during the development of antibiotic resistance in the clinical setting and thus may be the basis to uncover putative targets for novel treatment strategies.

Molecular evolution of the hyperthermophilic archaea of the Pyrococcus genus: analysis of adaptation to different environmental conditions

Background

Prokaryotic microorganisms are able to survive and proliferate in severe environmental conditions. The increasing number of complete sequences of prokaryotic genomes has provided the basis for studying the molecular mechanisms of their adaptation at the genomic level. We apply here a computer-based approach to compare the genomes and proteomes from P. furiosus, P. horikoshii, and P. abyssi to identify features of their molecular evolution related to adaptation strategy to diverse environmental conditions.

Results

Phylogenetic analysis of rRNA genes from 26 Pyrococcus strains suggested that the divergence of P. furiosus, P. horikoshii and P. abyssi might have occurred from ancestral deep-sea organisms. It was demonstrated that the function of genes that have been subject to positive Darwinian selection is closely related to abiotic and biotic conditions to which archaea managed to become adapted. Divergence of the P. furiosus archaea might have been due to loss of some genes involved in cell motility or signal transduction, and/or to evolution under positive selection of the genes for translation machinery. In the course of P. horikoshii divergence, positive selection was found to operate mainly on the transcription machinery; divergence of P. abyssi was related with positive selection for the genes mainly involved in inorganic ion transport. Analysis of radical amino acid replacement rate in evolving P. furiosus, P. horikoshii and P. abyssi showed that the fixation rate was higher for radical substitutions relative to the volume of amino acid side-chain.

Conclusions

The current results give due credit to the important role of hydrostatic pressure as a cause of variability in the P. furiosus, P. horikoshii and P. abyssi genomes evolving in different habitats. Nevertheless, adaptation to pressure does not appear to be the sole factor ensuring adaptation to environment. For example, at the stage of the divergence of P. horikoshii and P. abyssi, an essential evolutionary role may be assigned to changes in the trophic chain, namely, acquisition of a consumer status at a high (P. horikoshii) or low level (P. abyssi).