SEARCH · Search PubMed

Results for “protein function annotation”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,675 records · Page 93Linked to original sources

PubMed2006

The consensus coding sequences of human breast and colorectal cancers.

The elucidation of the human genome sequence has made it possible to identify genetic alterations in cancers in unprecedented detail. To begin a systematic analysis of such alterations, we determined the sequence of well-annotated human protein-coding genes in two common tumor types. Analysis of 13,023 genes in 11 breast and 11 colorectal cancers revealed that individual tumors accumulate an average of approximately 90 mutant genes but that only a subset of these contribute to the neoplastic process. Using stringent criteria to delineate this subset, we identified 189 genes (average of 11 per tumor) that were mutated at significant frequency. The vast majority of these genes were not known to be genetically altered in tumors and are predicted to affect a wide range of cellular functions, including transcription, adhesion, and invasion. These data define the genetic landscape of two human cancer types, provide new targets for diagnostic and therapeutic intervention, and open fertile avenues for basic research in tumor biology.

Amino Acid Substitution↗

PubMed2004

The single-stranded DNA-binding protein of Deinococcus radiodurans.

BACKGROUND: Deinococcus radiodurans R1 is one of the most radiation-resistant organisms known and is able to repair an unusually large amount of DNA damage without induced mutation. Single-stranded DNA-binding (SSB) protein is an essential protein in all organisms and is involved in DNA replication, recombination and repair. The published genomic sequence from Deinococcus radiodurans includes a putative single-stranded DNA-binding protein gene (ssb; DR0100) requiring a translational frameshift for synthesis of a complete SSB protein. The apparently tripartite gene has inspired considerable speculation in the literature about potentially novel frameshifting or RNA editing mechanisms. Immediately upstream of the ssb gene is another gene (DR0099) given an ssb-like annotation, but left unexplored. RESULTS: A segment of the Deinococcus radiodurans strain R1 genome encompassing the ssb gene has been re-sequenced, and two errors involving omitted guanine nucleotides have been documented. The corrected sequence incorporates both of the open reading frames designated DR0099 and DR0100 into one contiguous ssb open reading frame (ORF). The corrected gene requires no translational frameshifts and contains two predicted oligonucleotide/oligosaccharide-binding (OB) folds. The protein has been purified and its sequence is closely related to the Thermus thermophilus and Thermus aquaticus SSB proteins. Like the Thermus SSB proteins, the SSBDr functions as a homodimer. The Deinococcus radiodurans SSB homodimer stimulates Deinococcus radiodurans RecA protein and Escherichia coli RecA protein-promoted DNA three-strand exchange reactions with at least the same efficiency as the Escherichia coli SSB homotetramer. CONCLUSIONS: The correct Deinococcus radiodurans ssb gene is a contiguous open reading frame that codes for the largest bacterial SSB monomer identified to date. The Deinococcus radiodurans SSB protein includes two OB folds per monomer and functions as a homodimer. The Deinococcus radiodurans SSB protein efficiently stimulates Deinococcus radiodurans RecA and also Escherichia coli RecA protein-promoted DNA strand exchange reactions. The identification and purification of Deinococcus radiodurans SSB protein not only allows for greater understanding of the SSB protein family but provides an essential yet previously missing player in the current efforts to understand the extraordinary DNA repair capacity of Deinococcus radiodurans.

DNA Repair↗

PubMed2003

The Comparative Toxicogenomics Database (CTD).

The Mount Desert Island Biological Laboratory in Salsbury Cove, Maine, USA, is developing the Comparative Toxicogenomics Database (CTD), a community-supported genomic resource devoted to genes and proteins of human toxicologic significance. CTD will be the first publicly available database to a) provide annotated associations among genes, proteins, references, and toxic agents, with a focus on annotating data from aquatic and mammalian organisms; b) include nucleotide and protein sequences from diverse species; c) offer a range of analysis tools for customized comparative studies; and d) provide information to investigators on available molecular reagents. This combination of features will facilitate cross-species comparisons of toxicologically significant genes and proteins. These comparisons will promote understanding of molecular evolution, the significance of conserved sequences, the genetic basis of variable sensitivity to environmental agents, and the complex interactions between the environment and human health. CTD is currently under development, and the planned scope and functions of the database are described herein. The intent of this report is to invite community participation in the development of CTD to ensure that it will be a valuable resource for environmental health, molecular biology, and toxicology research.

Animals↗

PubMed2026

How Not to be Seen: Predicting Unseen Enzyme Functions using Contrastive Learning.

MOTIVATION: Predicting enzyme function from its sequence is still an unsolved problem in the life sciences. Moreover, with the explosion of annotated genome data, we are inundated with potential enzymatic sequences that have not yet been biochemically characterized. While it is not possible to assign a not-yet-existing label to such a sequence, there is high value in placing the sequence as accurately as possible in known function space. Doing so can help provide more accurate falsifiable hypotheses for experimentalists wishing to characterize enzymes from specific functional families. RESULTS: Here we present a contrastive learning algorithm for predicting enzyme function from sequence. Our method, EnzPlacer, predicts the third, second, and first EC numbers for a protein whose fourth EC number is not in the training corpus. This novel prediction mechanism accurately places a protein sequence within a narrowed-down functional context, even if the precise function remains unknown. AVAILABILITY: EnzPlacer is available from https://github.com/drxiangma/EnzPlacer under a GPL3 license.

Contrastive learning↗

PubMed2003

Phylogenomic analysis of the Giardia intestinalis transcarboxylase reveals multiple instances of domain fusion and fission in the evolution of biotin-dependent enzymes.

Sequencing of the gene encoding a pyruvate carboxylase-like protein from the amitochondrial eukaryote Giardia intestinalis revealed a 1,338 aa protein composed of acetyl-CoA carboxyltransferase (ACCT), pyruvate carboxyltransferase (PycB), and biotin carboxyl carrier protein (BCCP) domains, linked in a single polypeptide chain. This particular domain combination has been previously seen only in the methylmalonyl-CoA:pyruvate transcarboxylase from Propionibacterium freudenreichii, where each of these domains is encoded by an individual gene and forms a separate subunit. To get an insight into the evolutionary origin and biochemical function of the G. intestinalis enzyme, we compared its domain composition to those of other biotin-dependent enzymes and performed a phylogenetic analysis of each of its domains. The results obtained indicate that: (1) evolution of the BCCP domain included several domain fusion events, leading to the ACCT-BCCP and PycB-BCCP domain combinations; (2) fusions of the PycB and BCCP domains in pyruvate carboxylases and oxaloacetate decarboxylases occurred on several independent occasions in different prokaryotic lineages, probably due to selective pressure towards co-expression of these genes, and (3) because newly sequenced biotin-dependent enzymes are often misannotated in sequence databases, their annotation as either carboxylases, decarboxylases, or transcarboxylases has to rely on detailed analysis of their domain composition, operon organization of the corresponding genes, gene content in the particular genome, and phylogenetic analysis.

Acetyl-CoA Carboxylase↗

PubMed2001

From first base: the sequence of the tip of the X chromosome of Drosophila melanogaster, a comparison of two sequencing strategies.

We present the sequence of a contiguous 2.63 Mb of DNA extending from the tip of the X chromosome of Drosophila melanogaster. Within this sequence, we predict 277 protein coding genes, of which 94 had been sequenced already in the course of studying the biology of their gene products, and examples of 12 different transposable elements. We show that an interval between bands 3A2 and 3C2, believed in the 1970s to show a correlation between the number of bands on the polytene chromosomes and the 20 genes identified by conventional genetics, is predicted to contain 45 genes from its DNA sequence. We have determined the insertion sites of P-elements from 111 mutant lines, about half of which are in a position likely to affect the expression of novel predicted genes, thus representing a resource for subsequent functional genomic analysis. We compare the European Drosophila Genome Project sequence with the corresponding part of the independently assembled and annotated Joint Sequence determined through "shotgun" sequencing. Discounting differences in the distribution of known transposable elements between the strains sequenced in the two projects, we detected three major sequence differences, two of which are probably explained by errors in assembly; the origin of the third major difference is unclear. In addition there are eight sequence gaps within the Joint Sequence. At least six of these eight gaps are likely to be sites of transposable elements; the other two are complex. Of the 275 genes in common to both projects, 60% are identical within 1% of their predicted amino-acid sequence and 31% show minor differences such as in choice of translation initiation or termination codons; the remaining 9% show major differences in interpretation.

Animals↗

PubMed2026

Whole-exome characterization of host genetic variation in HIV-associated genes across the high-prevalence Mizo population, Northeast India.

BACKGROUND: The Mizoram state of Northeast India has one of the highest HIV prevalence rates in Asia, yet the host genetic factors influencing HIV susceptibility in this Tibeto-Burman population remain uncharacterised. METHODS: We performed whole-exome sequencing using Illumina NovaSeq 6000, mean coverage 100X on 76 HIV-negative Mizo individuals. Variants were called using GATK HaplotypeCaller v4.3 against GRCh38p14, annotated with ANNOVAR, and filtered using hard-quality thresholds (QD ≥ 2, SOR ≤ 3, MQ ≥ 40, DP ≥ 10, GQ ≥ 20). The allele frequencies were compared against gnomAD v2.1.1 population databases. Hardy-Weinberg equilibrium was assessed using the Wigginton exact test with Bonferroni correction. RESULTS: Post-quality filtering resulted in 12,011 sample-variants across 2,821 unique positions from 36 HIV-associated loci (33 protein-coding genes, 2 chemokine ligands, and 3 lncRNA targets). Of these, 784 observations (51 unique positions) were high-impact nonsynonymous or loss-of-function variants. ADAR rs2229857 (p.K384R, NM_015840) was the most frequently observed variant (Mizo carrier frequency = 0.895; 95% CI: 0.806-0.946). CXCR1 rs16858808 (p.R335C) showed the greatest population enrichment (Mizo carrier frequency = 0.197; 95% CI: 0.123-0.300; 7.65-fold carrier-frequency enrichment versus gnomAD South Asian; CADD = 15.60). Sixteen of 20 tested variants deviated from Hardy-Weinberg equilibrium after Bonferroni correction (p < 0.0025), predominantly showing excess homozygosity consistent with the endogamous Mizo population. The protective variant CCR5-Δ32 was absent in all the 76 individuals tested. CONCLUSION: This first whole-exome characterization of HIV host genes in the Mizo population identifies CXCR1 rs16858808 as the most population-enriched functional variant and reveals a pervasive endogamy signature. These findings provide a population-specific genetic framework for future HIV susceptibility studies and ART pharmacogenomics research.

Humans↗

PubMed2004

Large-scale analysis of the barley transcriptome based on expressed sequence tags.

To provide resources for barley genomics, 110,981 expressed sequence tags (ESTs) were generated from 22 cDNA libraries representing tissues at various developmental stages. This EST collection corresponds to approximately one-third of the 380,000 publicly available barley ESTs. Clustering and assembly resulted in 14,151 tentative consensi (TCs) and 11 073 singletons, altogether representing 25 224 putatively unique sequences. Of these, 17.5% showed no significant similarity to other barley ESTs present in dbEST. More than 41% of all barley genes are supposed to belong to multigene families and approximately 4% of the barley genes undergo alternative splicing. Based on the functional annotation of the set of unique sequences, the functional category 'Energy' was further analysed to reveal tissue- and stage-specific differences in gene expression. Hierarchical clustering of 362 differentially expressed TCs resulted in the identification of seven major clusters. The clusters reflect biochemical pathways predominantly activated in specific tissues and at various developmental stages. During seed germination glycolysis could be identified as the most predominant biochemical pathway. Germination-specific glycolysis is characterized by the coordinated expression of phosphoenolpyruvate carboxylase and phosphoenolpyruvate carboxykinase, whose antagonistic actions possibly regulate the flux of amino acids into protein biosynthesis and gluconeogenesis respectively. The expression of defence-related and antioxidant genes during germination might be controlled by the ethylene-signalling pathway as concluded from the coordinated expression of those genes and the transcription factors (TF) EIN3 and EREBPG. Moreover, because of their predominant expression in germinating seeds, TF of the AP2 and MYB type are presumably major regulators of germination.

Expressed Sequence Tags↗

PubMed2004

Transcriptional analysis of biofilm formation processes in the anaerobic, hyperthermophilic bacterium Thermotoga maritima.

Thermotoga maritima, a fermentative, anaerobic, hyperthermophilic bacterium, was found to attach to bioreactor glass walls, nylon mesh, and polycarbonate filters during chemostat cultivation on maltose-based media at 80 degrees C. A whole-genome cDNA microarray was used to examine differential expression patterns between biofilm and planktonic populations. Mixed-model statistical analysis revealed differential expression (twofold or more) of 114 open reading frames in sessile cells (6% of the genome), over a third of which were initially annotated as hypothetical proteins in the T. maritima genome. Among the previously annotated genes in the T. maritima genome, which showed expression changes during biofilm growth, were several that corresponded to biofilm formation genes identified in mesophilic bacteria (i.e., Pseudomonas species, Escherichia coli, and Staphylococcus epidermidis). Most notably, T. maritima biofilm-bound cells exhibited increased transcription of genes involved in iron and sulfur transport, as well as in biosynthesis of cysteine, thiamine, NAD, and isoprenoid side chains of quinones. These findings were all consistent with the up-regulation of iron-sulfur cluster assembly and repair functions in biofilm cells. Significant up-regulation of several beta-specific glycosidases was also noted in biofilm cells, despite the fact that maltose was the primary carbon source fed to the chemostat. The reasons for increased beta-glycosidase levels are unclear but are likely related to the processing of biofilm-based polysaccharides. In addition to revealing insights into the phenotype of sessile T. maritima communities, the methodology developed here can be extended to study other anaerobic biofilm formation processes as well as to examine aspects of microbial ecology in hydrothermal environments.

ATP-Binding Cassette Transporters↗

PubMed2005

Gene function prediction by a combined analysis of gene expression data and protein-protein interaction data.

Prediction of biological functions of genes is an important issue in basic biology research and has applications in drug discoveries and gene therapies. Previous studies have shown either gene expression data or protein-protein interaction data alone can be used for predicting gene functions. In particular, clustering gene expression profiles has been widely used for gene function prediction. In this paper, we first propose a new method for gene function prediction using protein-protein interaction data, which will facilitate combining prediction results based on clustering gene expression profiles. We then propose a new method to combine the prediction results based on either source of data by weighting on the evidence provided by each. Using protein-protein interaction data downloaded from the GRID database, published gene expression profiles from 300 microarray experiments for the yeast S. cerevisiae, we show that this new combined analysis provides improved predictive performance over that of using either data source alone in a cross-validated analysis of the MIPS gene annotations. Finally, we propose a logistic regression method that is flexible enough to combine information from any number of data sources while maintaining computational feasibility.

Algorithms↗

PubMed2005

A support vector machine approach to the identification of phosphorylation sites.

We describe a bioinformatics tool that can be used to predict the position of phosphorylation sites in proteins based only on sequence information. The method uses the support vector machine (SVM) statistical learning theory. The statistical models for phosphorylation by various types of kinases are built using a dataset of short (9-amino acid long) sequence fragments. The sequence segments are dissected around post-translationally modified sites of proteins that are on the current release of the Swiss-Prot database, and that were experimentally confirmed to be phosphorylated by any kinase. We represent them as vectors in a multidimensional abstract space of short sequence fragments. The prediction method is as follows. First, a given query protein sequence is dissected into overlapping short segments. All the fragments are then projected into the multidimensional space of sequence fragments via a collection of different representations. Those points are classified with pre-built statistical models (the SVM method with linear, polynomial and radial kernel functions) either as phosphorylated or inactive ones. The resulting list of plausible sites for phosphorylation by various types of kinases in the query protein is returned to the user. The efficiency of the method for each type of phosphorylation is estimated using leave-one-out tests and presented here. The sensitivities of the models can reach over 70%, depending on the type of kinase. The additional information from profile representations of short sequence fragments helps in gaining a higher degree of accuracy in some phosphorylation types. The further development of an automatic phosphorylation site annotation predictor based on our algorithm should yield a significant improvement when using statistical algorithms in order to quantify the results.

Algorithms↗

PubMedDate not supplied

Complement factor H: sequence analysis of 221 kb of human genomic DNA containing the entire fH, fHR-1 and fHR-3 genes.

Complement factor H (fH) is a member of a family of proteins involved in the regulation of complement activation (RCA). These proteins share a common structural motif, the Short Consensus Repeat (SCR), which is structurally conserved among related genes and between phylogenetically divergent species. fH is composed of 20 such SCRs and a variety of biological functions have been localised to specific SCR domains. The majority of individual SCRs identified are encoded by single exons, and processes such as gene conversion, duplication and exon shuffling have been implicated in the evolution and genomic radiation of SCR-encoding genes. We have analysed two GenBank sequence entries relating to two overlapping PAC clones sequenced at the Sanger Centre which contain the entire human fH gene and two adjacent fH-related (fHR) genes, fHR-1 and fHR-3. Here, we report the detailed analysis of the assembled 221 kb of contiguous, ungapped genomic sequence from human chromosome 1q32, in part employing the RUMMAGE-DP automated annotation tool. Genomic duplications involving fH and fHR exons were identified and Alu/L1 repeat dating established that the duplications occurred after the separation of rodent and primate lineages. The analysis indicates that retrotransposition as well as single and multiple exon duplication events are likely to have been involved in SCR radiation and RCA gene evolution, facilitated by conservation of splice-phasing and the single-exon, single-SCR nature of the encoded domains.

5' Untranslated Regions↗

PubMed2002

A study of the middle-scale nucleotide clustering in DNA sequences of various origin and functionality, by means of a method based on a modified standard deviation.

The deviation from randomness in the distribution of nucleotides in genomic sequences is quantified and studied, using a modified standard deviation (MSD). This method implies a "per block" computation of the standard deviation of the nucleotide frequencies of occurrence, using local means (means taken in a neighborhood of each block). This quantity may serve as a scale-dependent measure of the nucleotide clustering. In the present work, the meso-scale of tenths of nucleotides is principally explored, by means of suitably adjusted filter parameters. This length scale is of an order of magnitude not directly affected by the grammar and syntax rules of the protein-coding procedure, remaining shorter than the scale of appearance of large-scale characteristics of the genome. MSD has been found to distinguish systematically between the sequences of different origin and functionality. The most near-random are found to be coding sequences of prokaryotes, while in intronic and intergenic regions of eukaryotic genomes, extended clustering of similar nucleotides is observed. The distributions of MSD values of large collections of sequences are found to be in most cases characteristic of their biological role and origin. Protein- and non-coding, prokaryotic and eukaryotic DNA as well as promoter, rRNA, viral and organelle sequences have been examined. The presented results corroborate a recently proposed model for genome evolution. The method is also applied for an assessment of the annotation of ORFs taken from the complete genome of Saccharomyces cerevisiae.

Animals↗

PubMed2004

Prediction of protein subcellular locations using fuzzy k-NN method.

MOTIVATION: Protein localization data are a valuable information resource helpful in elucidating protein functions. It is highly desirable to predict a protein's subcellular locations automatically from its sequence. RESULTS: In this paper, fuzzy k-nearest neighbors (k-NN) algorithm has been introduced to predict proteins' subcellular locations from their dipeptide composition. The prediction is performed with a new data set derived from version 41.0 SWISS-PROT databank, the overall predictive accuracy about 80% has been achieved in a jackknife test. The result demonstrates the applicability of this relative simple method and possible improvement of prediction accuracy for the protein subcellular locations. We also applied this method to annotate six entirely sequenced proteomes, namely Saccharomyces cerevisiae, Caenorhabditis elegans, Drosophila melanogaster, Oryza sativa, Arabidopsis thaliana and a subset of all human proteins. AVAILABILITY: Supplementary information and subcellular location annotations for eukaryotes are available at http://166.111.30.65/hying/fuzzy_loc.htm

Algorithms↗

PubMed2004

VISTA: computational tools for comparative genomics.

Comparison of DNA sequences from different species is a fundamental method for identifying functional elements in genomes. Here, we describe the VISTA family of tools created to assist biologists in carrying out this task. Our first VISTA server at http://www-gsd.lbl.gov/vista/ was launched in the summer of 2000 and was designed to align long genomic sequences and visualize these alignments with associated functional annotations. Currently the VISTA site includes multiple comparative genomics tools and provides users with rich capabilities to browse pre-computed whole-genome alignments of large vertebrate genomes and other groups of organisms with VISTA Browser, to submit their own sequences of interest to several VISTA servers for various types of comparative analysis and to obtain detailed comparative analysis results for a set of cardiovascular genes. We illustrate capabilities of the VISTA site by the analysis of a 180 kb interval on human chromosome 5 that encodes for the kinesin family member 3A (KIF3A) protein.

Binding Sites↗

PubMed2006

SUBA: the Arabidopsis Subcellular Database.

Knowledge of protein localisation contributes towards our understanding of protein function and of biological inter-relationships. A variety of experimental methods are currently being used to produce localisation data that need to be made accessible in an integrated manner. Chimeric fluorescent fusion proteins have been used to define subcellular localisations with at least 1100 related experiments completed in Arabidopsis. More recently, many studies have employed mass spectrometry to undertake proteomic surveys of subcellular components in Arabidopsis yielding localisation information for approximately 2600 proteins. Further protein localisation information may be obtained from other literature references to analysis of locations (AmiGO: approximately 900 proteins), location information from Swiss-Prot annotations (approximately 2000 proteins); and location inferred from gene descriptions (approximately 2700 proteins). Additionally, an increasing volume of available software provides location prediction information for proteins based on amino acid sequence. We have undertaken to bring these various data sources together to build SUBA, a SUBcellular location database for Arabidopsis proteins. The localisation data in SUBA encompasses 10 distinct subcellular locations, >6743 non-redundant proteins and represents the proteins encoded in the transcripts responsible for 51% of Arabidopsis expressed sequence tags. The SUBA database provides a powerful means by which to assess protein subcellular localisation in Arabidopsis (http://www.suba.bcs.uwa.edu.au).

Arabidopsis Proteins↗

PubMed2004

Insights into the genomic basis of niche specificity of Pseudomonas putida KT2440.

A major challenge in microbiology is the elucidation of the genetic and ecophysiological basis of habitat specificity of microbes. Pseudomonas putida is a paradigm of a ubiquitous metabolically versatile soil bacterium. Strain KT2440, a safety strain that has become a laboratory workhorse worldwide, has been recently sequenced and its genome annotated. By drawing on both published information and on original in silico analysis of its genome, we address here the question of what genomic features of KT2440 could explain or are consistent with its ubiquity, metabolic versatility and adaptability. The genome of KT2440 exhibits combinations of features characteristic of terrestrial, rhizosphere and aquatic bacteria, which thrive in either copiotrophic or oligotrophic habitats, and suggests that P. putida has evolved and acquired functions that equip it to thrive in diverse, often inhospitable environments, either free-living, or in close association with plants. The high diversity of protein families encoded by its genome, the large number and variety of small aralogous families, insertion elements, repetitive extragenic palindromic sequences, as well as the mosaic structure of the genome (with many regions of 'atypical' composition) and the multiplicity of mobile elements, reflect a high functional diversity in P. putida and are indicative of its evolutionary trajectory and adaptation to the diverse habitats in which it thrives. The unusual wealth of determinants for high affinity nutrient acquisition systems, mono- and di-oxygenases, oxido-reductases, ferredoxins and cytochromes, dehydrogenases, sulfur metabolism proteins, for efflux pumps and glutathione-S-transfereases, and for the extensive array of extracytoplasmatic function sigma factors, regulators, and stress response systems, constitute the genomic basis for the exceptional nutritional versatility and opportunism of P. putida , its ubiquity in diverse soil, rhizosphere and aquatic systems, and its renowned tolerance of natural and anthropogenic stresses. This metabolic diversity is also the basis of the impressive evolutionary potential of KT2440, and its utility for the experimental design of novel pathways for the catabolism of organic, particularly aromatic, pollutants, and its potential for bioremediation of soils contaminated with such compounds as well as for its application in the production of high-added value compounds.

Adaptation, Physiological↗

PubMed2006

Cyanobacterial two-component proteins: structure, diversity, distribution, and evolution.

A survey of the already characterized and potential two-component protein sequences that exist in the nine complete and seven partially annotated cyanobacterial genome sequences available (as of May 2005) showed that the cyanobacteria possess a much larger repertoire of such proteins than most other bacteria. By analysis of the domain structure of the 1,171 potential histidine kinases, response regulators, and hybrid kinases, many various arrangements of about thirty different modules could be distinguished. The number of two-component proteins is related in part to genome size but also to the variety of physiological properties and ecophysiologies of the different strains. Groups of orthologues were defined, only a few of which have representatives with known physiological functions. Based on comparisons with the proposed phylogenetic relationships between the strains, the orthology groups show that (i) a few genes, some of them clustered on the genome, have been conserved by all species, suggesting their very ancient origin and an essential role for the corresponding proteins, and (ii) duplications, fusions, gene losses, insertions, and deletions, as well as domain shuffling, occurred during evolution, leading to the extant repertoire. These mechanisms are put in perspective with the different genetic properties that cyanobacteria have to achieve genome plasticity. This review is designed to serve as a basis for orienting further research aimed at defining the most ancient regulatory mechanisms and understanding how evolution worked to select and keep the most appropriate systems for cyanobacteria to develop in the quite different environments that they have successfully colonized.

Bacterial Proteins↗