Conservation of adjacency as evidence of paralogous operons

Abstract

Most of the analyses on the conservation of gene order are limited to orthologous genes. However, the organization of genes into operons might also result in the conservation of gene order of paralogous genes. Thus, we sought computational evidence that conservation of gene order of paralogous genes represents another level of conservation of genes in operons. We found that pairs of genes within experimentally characterized operons of Escherichia coli K12 and Bacillus subtilis tend to have more adjacently conserved paralogs than pairs of genes at transcription unit boundaries. The fraction of same strand gene pairs corresponding to conserved paralogs averages 0.07 with a maximum of 0.22 in Borrelia burgdorferi. The use of evidence from the conservation of adjacency of paralogous genes can improve the prediction of operons in E.coli K12 by ∼0.27 over predictions using conservation of adjacency of orthologous genes alone.