Analysis of invariant sequences in 266 complete genomes

The systematic analysis of ultraconserved genomic regions in the budding yeast

Motivation

In the evolution of species, a kind of special sequences, termed ultraconserved sequences (UCSs), have been inherited without any change, which strongly suggests those sequences should be crucial for the species to survive or adapt to the environment. However, the UCSs are still regarded as mysterious genetic sequences so far. Here, we present a systematic study of ultraconserved genomic regions in the budding yeast based on the publicly available genome sequences, in order to reveal their relationship with the adaptability or fitness advantages of the budding yeast.

Results

Our results indicate that, in addition to some fundamental biological functions, the UCSs play an important role in the adaptation of Saccharomyces cerevisiae to the acidic environment, which is backed up by the previous observation. Besides that, we also find the highly unchanged genes are enriched in some other pathways, such as the nutrient-sensitive signaling pathway. To facilitate the investigation of unique UCSs, the UCSC Genome Browser was utilized to visualize the chromosomal position and related annotations of UCSs in S.cerevisiae genome.

Availability and implementation

For more details on UCSs, please refer to the Supplementary information online, and the custom code is available on request.

Contact

fgao@tju.edu.cn.

Supplementary information

Supplementary data are available at Bioinformatics online.

Uniquely localized intra-molecular amino acid concentrations at the glycolytic enzyme catalytic/active centers of Archaea, Bacteria and Eukaryota are associated with their proposed temporal appearances on earth

The distributions of amino acids at most-conserved sites nearest catalytic/active centers (C/AC) in 4,645 sequences of ten enzymes of the glycolytic Embden-Meyerhof-Parnas pathway in Archaea, Bacteria and Eukaryota are similar to the proposed temporal order of their appearance on Earth. Glycine, isoleucine, leucine, valine, glutamic acid and possibly lysine often described as prebiotic, i.e., existing or occurring before the emergence of life, were localized in positional and conservational defined aggregations in all enzymes of all Domains. The distributions of all 20 biologic amino acids in most-conserved sites nearest their C/ACs were quite different either from distributions in sites less-conserved and further from their C/ACs or from all amino acids regardless of their position or conservation. The major concentrations of glycine, e.g., perhaps the earliest prebiotic amino acid, occupies ≈ 16 % of all the most-conserved sites within a volume of ≈ 7-8 Å radius from their C/ACs and decreases linearly towards the molecule's peripheries. Spatially localized major concentrations of isoleucine, leucine and valine are in the mid-conserved and mid-distant sites from their C/ACs in protein interiors. Lysine and glutamic acid comprise ≈ 25-30 % of all amino acids within an irregular volume bounded by ≈ 24-28 Å radii from their C/ACs at the most-distant least-conserved sites. The unreported characteristics of these amino acids: their spatially and conservationally identified concentrations in Archaea, Bacteria and Eukaryota, suggest some common structural organization of glycolytic enzymes that may be relevant to their evolution and that of other proteins. We discuss our data in relation to enzyme evolution, their reported prebiotic putative temporal appearances on Earth, abundances, biological "cost", neighbor-sequence preferences or "ordering" and some thermodynamic parameters.

Towards functional repertoire of the earliest proteins

The conserved protein sequence motifs present in all prokaryotic proteomes, "omnipresent motifs," presumably, correspond to the earliest proteins of the Last Universal Cellular Ancestor, from which all the proteomes have descended. Fifteen proteomes, each representing one of the total 15 diverse phyla of 131 Eubacteria and Archea, from which the omnipresent elements have been originally derived, are exhaustively screened. All those proteins which harbor the omnipresent motifs are identified. Six "omnipresent" protein types are revealed which are located in all 15 proteomes: ABC cassettes, FtsH proteases, translation initiation factors, translation elongation factors, isoleucyl-tRNA synthases, and RNA polymerases β'. In addition to the omnipresent motifs, these proteins also contain other highly conserved motifs, standing for additional modules of the proteins. Remarkably, the identified tentative earliest proteins are responsible for only three basic functions: supply of monomers (ABC transporters and proteases), protein synthesis (initiation and elongation factors, aminoacyl-tRNA synthases), and RNA synthesis (polymerases). No enzymes involved in metabolic activities are present in the list of the earliest proteins derived by this approach. Some of the omnipresent sequence motifs are found, indeed, in the metabolic enzymes (e.g. NTP binding motifs), but these enzymes do not make a sequence matching collection of 15 sequences, i.e. they are not omnipresent. Future analysis of less conserved sequence motifs may reveal at what degree of conservation (stage of evolution) the metabolic enzymes could have entered the scene.

Robustness assessment of whole bacterial genome segmentations

Comparison of closely related bacterial genomes has revealed the presence of highly conserved sequences forming a "backbone" that is interrupted by numerous, less conserved, DNA fragments. Segmentation of bacterial genomes into backbone and variable regions is particularly useful to investigate, among other things, bacterial genome evolution. Several software tools have been designed to compare complete bacterial chromosomes and a few online databases store pre-computed genome comparisons. However, very few statistical methods are available to evaluate the reliability of these software tools and to compare the results obtained with them. To fill this gap, we have developed two local scores to measure the robustness of bacterial genome segmentations. Our method uses a simulation procedure based on random perturbations of the compared genomes. The two scores described in this article provide useful information and are easy to implement, and their interpretation is intuitive. We show that they are suited to discriminate between robust and non-robust segmentations when genome aligners such as MAUVE and MGA are used.