- Research article
- Open Access
Gene socialization: gene order, GC content and gene silencing in Salmonella
© Papanikolaou et al; licensee BioMed Central Ltd. 2009
- Received: 10 June 2009
- Accepted: 11 December 2009
- Published: 11 December 2009
Genes of conserved order in bacterial genomes tend to evolve slower than genes whose order is not conserved. In addition, genes with a GC content lower than the GC content of the resident genome are known to be selectively silenced by the histone-like nucleoid structuring protein (H-NS) in Salmonella.
In this study, we use a comparative genomics approach to demonstrate that in Salmonella, genes whose order is not conserved (or genes without homologs) in closely related bacteria possess a significantly lower average GC content in comparison to genes that preserve their relative position in the genome. Moreover, these genes are more frequently targeted by H-NS than genes that have conserved their genomic neighborhood. We also observed that duplicated genes that do not preserve their genomic neighborhood are, on average, under less selective pressure.
We establish a strong association between gene order, GC content and gene silencing in a model bacterial species. This analysis suggests that genes that are not under strong selective pressure (evolve faster than others) in Salmonella tend to accumulate more AT-rich mutations and are eventually silenced by H-NS. Our findings may establish new approaches for a better understanding of bacterial genome evolution and function, using information from functional and comparative genomics.
- Duplication Event
- Gene Order
- Duplicate Gene
- Essential Gene
- Lateral Gene Transfer
The conservation of gene order has been found to play an essential role in genome evolution. More specifically, proteins encoded by genes of conserved order in bacteria tend to evolve more slowly when compared to proteins encoded by genes without a conserved order [1, 2] and genes with similar or related functions tend to occur in adjacent chromosomal positions in yeast . Moreover, genes with conserved order were found to evolve at similar rates  and, in prokaryotes, proteins encoded by genes with conserved order appear to interact physically . It has also been shown that in eukaryotes essential genes are clustered in regions with low recombination rates , whereas in bacteria essential genes are more conserved than non essential genes . In addition, it has been reported that the number of interactions involved in a protein network is directly correlated with the rate of evolution among these proteins  and that highly expressed genes evolve slowly [8, 9]. A case in point for the role of gene order in evolution can be illustrated by duplicated genes . Following a duplication event one of the two (paralog) genes might keep its original function, whereas the other one might be under less selective pressure. Yet, it is not always readily apparent which duplicated genes evolve faster. There have been reports that have marginally correlated sequence conservation with genome context , but there must be other, yet unknown, functional features that determine the fate of duplicate genes.
Recently, two research groups observed independently that the histone-like nucleoid structuring protein (H-NS) plays an important role as a general transcriptional repressor of a large number of genes in Salmonella enterica serovar Typhimurium LT2 (S. Typhimurium) [12, 13]. H-NS is a protein that is believed to play an essential role in the organization and compaction of bacterial chromatin as well as in transcriptional regulation for many bacterial genes [14–16]. H-NS binds to these genes and silences them transcriptionally. An apparent common feature of genes silenced by H-NS is that their GC content is significantly lower than the overall GC content of the Salmonella genome . Additionally, a large proportion of these H-NS repressed genes is predicted to have been acquired from a foreign source (horizontally transferred genes, HTGs) [13, 15], a fact in agreement with the observation that HTGs are relatively AT-rich [17, 18]. It has been suggested that this may be a defensive mechanism against foreign genetic material without loosing the benefit of future usage of this material if necessary .
Herein, we aimed to test the correlation between gene order conservation, gene duplication and H-NS dependent silencing in S. Typhimurium, using Escherichia coli K12 as a reference genome in order to identify the conservation or loss of gene order along the bacterial chromosome. We have also attempted to associate the above mentioned features with GC content and gene essentiality.
Evolutionary rate and GC content are related to genomic neighborhood conservation
In order to check the former hypothesis, we cross-checked the three groups for the existence of HTGs in the Salmonella genome using HGT-DB, a database of putative horizontally transferred genes in prokaryotes . As expected, the majority of predicted HTGs fall into the last two categories (nGCO, no homolog), but after removing HTGs from the complete gene set, the results were not affected (Additional file 2). However, this observation could be biased because the methodological approach used in HGT-DB to infer lateral gene transfer events, is largely dependent on GC content . The latter hypothesis could be explained by our finding that the last two categories of genes may be under weaker selective pressure than the first one and, consequently, accumulating a larger number of mutations, especially C-to-T mutations which are known to be the most common ones under neutral evolution . The observation that genomes of intracellular pathogens tend to be AT rich for energetic and resources availability reasons  could also explain the aforementioned result. After a speciation event occurs, some genes might remain GCO and continue to evolve slowly. This is probably the result of the strong evolutionary pressure to maintain their relative order, assisting the preservation of protein function or interaction with neighboring genes. On the other hand, other genes might lose their relative genomic position and start evolving faster. Genes in Salmonella with no homologs in E. coli K12 could either be HTGs or genes that have lost their genomic neighborhood and subsequently evolved to such an extent that no significant sequence similarities could be detected to support homology to any E. coli K12 gene.
The fate of duplicated genes is determined by gene order conservation
The most striking evidence that GCO genes evolve, on average, significantly slower than the nGCO genes, is inferred from an analysis of the subset of duplicated genes. We performed an all-against-all self-comparison of proteins encoded in the Salmonella genome and detected all groups of genes that could be assigned as homologs in the Salmonella genome that share the same best hit in the E. coli K12 genome (Additional file 3). We identified 687 such genes (15.5% of the complete Salmonella gene set) 343 of which are GCO (49.9%) and 344 are nGCO (51.1%). When compared to the same figures for the complete Salmonella gene set (68.3% and 12.6% respectively) the enrichment of the duplicate gene set in nGCO genes is evident. This enrichment probably takes place because, after a duplication event occurs, the selective pressure on one gene to retain its genomic neighborhood becomes less stringent. The proteins of GCO genes in the 343 duplicates share a 74.4% identity (standard deviation sd = 20.2), whereas for the 344 proteins of nGCO genes this figure drops to 37.0% (sd = 14.6; Wilcoxon rank-sum test: W = 106315.5, P-value = 0) (Fig. 1A). This result leads us to two major conclusions: First, as it may be deduced from a comparison of these identity percentages with the result for the complete genome set, the proteins encoded by this subset of duplicated genes evolve faster (in agreement with previous studies ). Second, whenever a gene duplication event occurs, the gene that retains its gene order tends to evolve more slowly than the one that travels away from the ancestral gene neighborhood and probably it is the former that keeps the original function. This finding complements an earlier study , although different experimental approaches have been employed here. Our results, even though based on the average conservation of putative homologs, corroborate the finding that genomic neighborhood conservation is an important additional criterion that should be taken into account in protein sequence similarity searches. We reached the same conclusions after calculating the Ka/Ks ratio for the two different duplicated gene groups - GCO and nGCO - with respective values 0.09 and 0.43 (Fig. 2). In most cases where all duplicated genes were nGCO, none of the genes exhibited particularly high identity to the common E. coli K12 homolog (Additional file 4). An analysis of the average GC content for the duplicated genes (GCO and nGCO) showed a similar pattern with the one of the overall genome (Fig. 1B).
Essential genes tend to conserve their order
H-NS preferentially regulates genes that do not conserve their order
Driven by recent results demonstrating that H-NS silences a significant amount of genes thought to have been acquired by lateral gene transfer, we calculated the number of GCO and nGCO genes in the dataset of 359 genes that were found to be silenced by the action of H-NS in the Salmonella chromosome . We ended up with the following numbers: 119 genes were GCO (33.1%), 99 (27.6%) were nGCO and 141 (39.3%) were without any detectable homolog in E. coli K12. The difference in the gene order status of genes compared to the overall genome is striking, with the representation being 68.3%, 12.6% and 19% respectively (see also fig. 3). This result demonstrates a preference of H-NS to repress genes that are either nGCO or that have no homologs in closely related species. The 119 GCO genes shared 74.7% identity and the nGCO 34.2% (P-value < 2.2e-16), a finding consistent with the results for both the overall genome set and the set of duplicated genes (Fig. 1A).
Gene order conservation of HNS repressed genes not predicted to be HTGs is associated to their GC content.
Gene order status
Number of Genes (%)
Average %GC content (sd)
No E. coli K12 homolog
In order to clarify whether the observed differences are due to pathogenicity related genes, we also compared Salmonella with an enteropathogenic strain of E. coli (O127:H6 str. E2348/69) (EPEC) ftp://ftp.ncbi.nih.gov/genomes/Bacteria/. We identified the complement of genes that had no homolog between Salmonella and E. coli K12 and had a homolog between Salmonella and EPEC (Additional file 5). Among the H-NS repressed genes within this dataset (14 genes out of 174 in total), we detected genes of Salmonella Pathogenicity Island-2 (SPI-2) such as ssaL, S , sseB and sscA . We also identified Salmonella genes whose order switches from nGCO to GCO when the reference genome is K12 and EPEC, respectively (Additional file 6). In this dataset, the H-NS repressed genes (11 genes out of 55 in total) include genes of SPI-2 such as ssaV, N, R, T, U .
Another interesting observation in the 359 H-NS repressed genes was the identification of 98 duplicated genes within this dataset (27.3% compared to 15.5% of duplicated genes in the whole genome) out of which only 29 were GCO. We subsequently analyzed the GC content for the duplicated genes of which only one paralog is targeted by H-NS. Of the 98 duplicated genes, a subset of 74 (out of which 22 were GCO) were identified to have 89 paralogs which are not H-NS repressed. After comparing the GC content of the above subsets of paralogs, we observed a striking difference. More specifically, the GC content of the 74 H-NS repressed genes was 46.5% (sd = 5.5) and for their 89 paralogs the GC content was 51.4% (sd = 4.2) (Wilcoxon rank-sum test: W = 1545, P-value = 5.663e-09). The above results provide a strong indication that, apart from repressing horizontally transferred genes, H-NS may also play a role in the transcriptional regulation of duplicated genes. Interestingly, genes that play an important role in the pathogenicity of Salmonella such as ssrB (essential for the expression of SPI-2 genes  genes) and ssaN, R, T, U and V (genes of SPI-2 that encode components of the second type III secretion apparatus ) were found to be nGCO, with low GC content (43.1 average GC%) whereas their corresponding duplicates were GCO with significantly higher GC content (54.79 average GC%). Moreover, the aforementioned virulence associated genes were identified to be repressed by H-NS , whereas their duplicates' expression was not significantly affected by H-NS.
Surprisingly, 20 genes out of those 359 were found to be essential for Salmonella. This paradox may be explained either by the growth conditions under which experiments where carried out  or by introducing H-NS as a part of a general, yet currently unknown, bacterial regulatory mechanism.
Taken together, this work shows that in Salmonella: (a) the proteins encoded by GCO genes tend to be more conserved (Fig. 1A and Fig. 2). Probably this is due to a correlation of conservation of gene order and the number of interacting proteins and/or expression levels of proteins encoded by genes belonging to the conserved clusters. Although it is known that protein expression levels are correlated to protein evolutionary rates [8, 9], and it is rather controversial whether there is a correlation between the connectivity of nodes in protein interaction networks and evolutionary rate [2, 7, 33–35], it still remains to be shown if any of these findings correlate to gene order conservation as well. (b) GCO genes have on average higher GC content when compared to genes that are nGCO and even higher than genes that do not share any homology to E. coli K12 genes (Fig. 1B). We speculate that this is a result of a mechanism that leads genes which evolve faster than others (and are not under strong selective pressure) to accumulate AT-rich mutations. (c) Genes that remain GCO have on average significantly lower Ka/Ks ratio than the nGCO ones, indicating that they are under higher selective pressure (Fig. 2). (d) Whenever a duplication event takes place, the duplicate that keeps its original order, is the one that keeps the highest sequence similarity (at least on average) and possibly the original function (in agreement to ). On the other hand, when the duplicated gene does not preserve its original order, it is subject to much faster evolution and less selective pressure as it is shown by comparing the respective Ka/Ks ratios (Fig. 2). Overrepresentation of nGCO genes in the duplicated group demonstrates that selective pressure on genes to keep their order is more relaxed after a duplication event. Interestingly, after a duplication event occurs and none of the paralogs retain their ancestral position in the genome, we failed to find any preference for sequence conservation (Additional file 4). (e) Finally, we demonstrated that the silencing mechanism via the H-NS protein in Salmonella significantly affects genes that are nGCO or genes that do not have any (detectable) homolog in E. coli K12. This could be another indication that H-NS might play a role, not only as a defensive instrument against foreign acquired DNA [12, 13, 36] but also as an internal control of native genes , as it could be of benefit to the organism to have them repressed for a variety of reasons associated with genome structure, for duplicated genes, non-essential genes or genes that did not retain their order.
In our analysis we tried to correlate gene order conservation with gene duplication, GC content and gene silencing in Salmonella enterica serovar Typhimurium LT2 based on two publications concerning H-NS silencing in the aforementioned species. We found that genes of conserved order not only evolve slower, as previously reported, but also tend to have higher GC content than the ones that do not conserve their order. We also found that essential genes tend to retain their genome position. Another finding was that after a duplication event takes place the gene that keeps its original position evolves slower than the one that looses it and probably keeps its primary function. Our study also shows that H-NS protein, a general repressor of a large number of genes in Salmonella, tends to avoid binding on genes with a conserved order. There is a clear indication that from the duplicated genes H-NS "prefers" to repress those ones that do not have a conserved gene order.
Definition of genes of conserved order
As a measure of conservation at the amino acid sequence level, we used the percentage of identical residues reported by BLASTP, while disregarding conservative substitutions and gaps. We follow the same criteria of gene order conservation as in . More specifically, given two genomes A and B encoding Na and Nb genes respectively, let us assume these genes are ordered a1, a2, ..., aNa and b1, b2, ..., bNb, respectively. Then, two genes ai and bj (1<i<Na, 1<j<Nb) are defined as GCO when the following two conditions are met: (1) ai, bj exhibit statistically significant sequence similarity; (2) at least one of ai-1, ai+1 exhibits statistically significant sequence similarity with one of bj-1, bj+1.
We applied a relatively strict E-value cut-off (E ≤ 10-6) and an additional requirement of 50% length coverage of both query and target sequences by the BLAST alignment. This approach was chosen in order to avoid spurious matches due to short conserved sequence motifs. The same criteria were applied both in the cross-genome comparison and in the self- comparison as presented in text. Differences from the results reported in  originate from the fact that we have used the latest versions of the E. coli K12 and Salmonella genomes, the difference in the E-value cut-off and the introduction of the alignment coverage criterion in this work.
We used the pairwise alignments reported by BLASTP to guide the construction of alignments of the corresponding coding regions, inserting in-frame gaps where necessary. For the Ka/Ks ratio calculation we used the yn00 program of the PAML package [25, 26] with the following parameters:
noisy = 1
icode = 10 * 0:universal code; 1:mammalian mt; 2-10:see below
weighting = 0 * weighting pathways between codons (0/1)?
commonf3 × 4 = 1 * use one set of codon freqs for all pairs (0/1)?
We have removed two genes (NP_462323 and NP_463455) for which the program couldn't provide us with a reliable Ka/Ks ratio, as well as all genes with Ka/Ks = 0 (Additional file 7). In order to avoid probable cases of saturated Ks values, where the estimated Ka/Ks ratio may be unreliable, we have further filtered our dataset to contain only those homologs with Ks<1.5 according to . Nevertheless, since sequences of nGCO genes on average seem to accumulate more non-synonymous substitutions than those of GCO genes, the difference in the Ka/Ks ratio between GCO and nGCO genes would probably be higher than the one reported here (Additional file 8, 9, 10.
All statistical calculations have been performed with R (Version 2.3.1) R Development Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria, 2004. URL http://www.R-project.org. 3-900051-07-0.
We thank Christos Ouzounis (Kings College, London), Pantelis Georgiades (University of Cyprus) and Ioannis Kyrmitzoglou (University of Cyprus) for critically reading of the manuscript, useful discussions and suggestions. Funding to pay the Open Access publication charges for this article was provided by the University of Cyprus. We thank Mpriki for great inspiration. This paper is dedicated to the memory of Prof. Costas Kastritsis.
- Dandekar T, Snel B, Huynen M, Bork P: Conservation of gene order: a fingerprint of proteins that physically interact. Trends Biochem Sci. 1998, 23 (9): 324-328. 10.1016/S0968-0004(98)01274-2.View ArticlePubMedGoogle Scholar
- Theodosiou T, Iliopoulos I: Protein sequences of linked genes are highly conserved in two bacterial species. J Evol Biol. 2006, 19 (4): 1343-1345. 10.1111/j.1420-9101.2006.01093.x.View ArticlePubMedGoogle Scholar
- Cohen BA, Mitra RD, Hughes JD, Church GM: A computational analysis of whole-genome expression data reveals chromosomal domains of gene expression. Nat Genet. 2000, 26 (2): 183-186. 10.1038/79896.View ArticlePubMedGoogle Scholar
- Williams EJ, Hurst LD: The proteins of linked genes evolve at similar rates. Nature. 2000, 407 (6806): 900-903. 10.1038/35038066.View ArticlePubMedGoogle Scholar
- Pal C, Hurst LD: Evidence for co-evolution of gene order and recombination rate. Nat Genet. 2003, 33 (3): 392-395. 10.1038/ng1111.View ArticlePubMedGoogle Scholar
- Jordan IK, Rogozin IB, Wolf YI, Koonin EV: Essential genes are more evolutionarily conserved than are nonessential genes in bacteria. Genome Res. 2002, 12 (6): 962-968.PubMed CentralView ArticlePubMedGoogle Scholar
- Fraser HB, Hirsh AE, Steinmetz LM, Scharfe C, Feldman MW: Evolutionary rate in the protein interaction network. Science. 2002, 296 (5568): 750-752. 10.1126/science.1068696.View ArticlePubMedGoogle Scholar
- Drummond DA, Bloom JD, Adami C, Wilke CO, Arnold FH: Why highly expressed proteins evolve slowly. Proc Natl Acad Sci USA. 2005, 102 (40): 14338-14343. 10.1073/pnas.0504070102.PubMed CentralView ArticlePubMedGoogle Scholar
- Pal C, Papp B, Hurst LD: Highly expressed genes in yeast evolve slowly. Genetics. 2001, 158 (2): 927-931.PubMed CentralPubMedGoogle Scholar
- Ohno S: Evolution by gene duplication. 1970, Berlin, New York: Springer-VerlagView ArticleGoogle Scholar
- Notebaart RA, Huynen MA, Teusink B, Siezen RJ, Snel B: Correlation between sequence conservation and the genomic context after gene duplication. Nucleic Acids Res. 2005, 33 (19): 6164-6171. 10.1093/nar/gki913.PubMed CentralView ArticlePubMedGoogle Scholar
- Lucchini S, Rowley G, Goldberg MD, Hurd D, Harrison M, Hinton JC: H-NS mediates the silencing of laterally acquired genes in bacteria. PLoS Pathog. 2006, 2 (8): e81-10.1371/journal.ppat.0020081.PubMed CentralView ArticlePubMedGoogle Scholar
- Navarre WW, Porwollik S, Wang Y, McClelland M, Rosen H, Libby SJ, Fang FC: Selective silencing of foreign DNA with low GC content by the H-NS protein in Salmonella. Science. 2006, 313 (5784): 236-238. 10.1126/science.1128794.View ArticlePubMedGoogle Scholar
- Dorman CJ: H-NS: a universal regulator for a dynamic genome. Nat Rev Microbiol. 2004, 2 (5): 391-400. 10.1038/nrmicro883.View ArticlePubMedGoogle Scholar
- Dorman CJ: H-NS, the genome sentinel. Nat Rev Microbiol. 2007, 5 (2): 157-161. 10.1038/nrmicro1598.View ArticlePubMedGoogle Scholar
- Pflum MK: H-NS gives invading DNA the silent treatment. Nat Chem Biol. 2006, 2 (8): 400-401. 10.1038/nchembio0806-400.View ArticlePubMedGoogle Scholar
- Lawrence JG, Ochman H: Amelioration of bacterial genomes: rates of change and exchange. J Mol Evol. 1997, 44 (4): 383-397. 10.1007/PL00006158.View ArticlePubMedGoogle Scholar
- Medigue C, Rouxel T, Vigier P, Henaut A, Danchin A: Evidence for horizontal gene transfer in Escherichia coli speciation. J Mol Biol. 1991, 222 (4): 851-856. 10.1016/0022-2836(91)90575-Q.View ArticlePubMedGoogle Scholar
- McClelland M, Sanderson KE, Spieth J, Clifton SW, Latreille P, Courtney L, Porwollik S, Ali J, Dante M, Du F, et al: Complete genome sequence of Salmonella enterica serovar Typhimurium LT2. Nature. 2001, 413 (6858): 852-856. 10.1038/35101614.View ArticlePubMedGoogle Scholar
- Blattner FR, Plunkett G, Bloch CA, Perna NT, Burland V, Riley M, Collado-Vides J, Glasner JD, Rode CK, Mayhew GF, et al: The complete genome sequence of Escherichia coli K-12. Science. 1997, 277 (5331): 1453-1474. 10.1126/science.277.5331.1453.View ArticlePubMedGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralView ArticlePubMedGoogle Scholar
- Garcia-Vallve S, Guzman E, Montero MA, Romeu A: HGT-DB: a database of putative horizontally transferred genes in prokaryotic complete genomes. Nucleic Acids Res. 2003, 31 (1): 187-189. 10.1093/nar/gkg004.PubMed CentralView ArticlePubMedGoogle Scholar
- Coulondre C, Miller JH, Farabaugh PJ, Gilbert W: Molecular basis of base substitution hotspots in Escherichia coli. Nature. 1978, 274 (5673): 775-780. 10.1038/274775a0.View ArticlePubMedGoogle Scholar
- Rocha EP, Danchin A: Base composition bias might result from competition for metabolic resources. Trends Genet. 2002, 18 (6): 291-294. 10.1016/S0168-9525(02)02690-2.View ArticlePubMedGoogle Scholar
- Yang Z, Nielsen R: Estimating synonymous and nonsynonymous substitution rates under realistic evolutionary models. Mol Biol Evol. 2000, 17 (1): 32-43.View ArticlePubMedGoogle Scholar
- Yang Z: PAML: a program package for phylogenetic analysis by maximum likelihood. Comput Appl Biosci. 1997, 13 (5): 555-556.PubMedGoogle Scholar
- Kondrashov FA, Rogozin IB, Wolf YI, Koonin EV: Selection in the evolution of gene duplications. Genome Biol. 2002, 3 (2): RESEARCH0008-10.1186/gb-2002-3-2-research0008.PubMed CentralView ArticlePubMedGoogle Scholar
- Knuth K, Niesalla H, Hueck CJ, Fuchs TM: Large-scale identification of essential Salmonella genes by trapping lethal insertions. Mol Microbiol. 2004, 51 (6): 1729-1744. 10.1046/j.1365-2958.2003.03944.x.View ArticlePubMedGoogle Scholar
- Iguchi A, Thomson NR, Ogura Y, Saunders D, Ooka T, Henderson IR, Harris D, Asadulghani M, Kurokawa K, Dean P, et al: Complete genome sequence and comparative genome analysis of enteropathogenic Escherichia coli O127:H6 strain E2348/69. J Bacteriol. 2009, 191 (1): 347-354. 10.1128/JB.01238-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Hensel M, Shea JE, Raupach B, Monack D, Falkow S, Gleeson C, Kubo T, Holden DW: Functional analysis of ssaJ and the ssaK/U operon, 13 genes encoding components of the type III secretion apparatus of Salmonella Pathogenicity Island 2. Mol Microbiol. 1997, 24 (1): 155-167. 10.1046/j.1365-2958.1997.3271699.x.View ArticlePubMedGoogle Scholar
- Ruiz-Albert J, Mundy R, Yu XJ, Beuzon CR, Holden DW: SseA is a chaperone for the SseB and SseD translocon components of the Salmonella pathogenicity-island-2-encoded type III secretion system. Microbiology. 2003, 149 (Pt 5): 1103-1111. 10.1099/mic.0.26190-0.View ArticlePubMedGoogle Scholar
- Fass E, Groisman EA: Control of Salmonella pathogenicity island-2 gene expression. Curr Opin Microbiol. 2009, 12 (2): 199-204. 10.1016/j.mib.2009.01.004.PubMed CentralView ArticlePubMedGoogle Scholar
- Batada NN, Hurst LD, Tyers M: Evolutionary and physiological importance of hub proteins. PLoS Comput Biol. 2006, 2 (7): e88-10.1371/journal.pcbi.0020088.PubMed CentralView ArticlePubMedGoogle Scholar
- Brown KR, Jurisica I: Unequal evolutionary conservation of human protein interactions in interologous networks. Genome Biol. 2007, 8 (5): R95-10.1186/gb-2007-8-5-r95.PubMed CentralView ArticlePubMedGoogle Scholar
- Jordan IK, Wolf YI, Koonin EV: No simple dependence between protein evolution rate and the number of protein-protein interactions: only the most prolific interactors tend to evolve slowly. BMC Evol Biol. 2003, 3: 1-10.1186/1471-2148-3-1.PubMed CentralView ArticlePubMedGoogle Scholar
- Navarre WW, McClelland M, Libby SJ, Fang FC: Silencing of xenogeneic DNA by H-NS-facilitation of lateral gene transfer in bacteria by a defense system that recognizes foreign DNA. Genes Dev. 2007, 21 (12): 1456-1471. 10.1101/gad.1543107.View ArticlePubMedGoogle Scholar
- Castillo-Davis CI, Hartl DL, Achaz G: cis-Regulatory and protein evolution in orthologous and duplicate genes. Genome Res. 2004, 14 (8): 1530-1536. 10.1101/gr.2662504.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.