Relationship between operon preference and functional properties of persistent genes in bacterial genomes

The Fate of Duplicated Enzymes in Prokaryotes: The Case of Isomerases

The isomerases are a unique enzymatic class of enzymes that carry out a great diversity of chemical reactions at the intramolecular level. This class comprises about 300 members, most of which are involved in carbohydrate and terpenoid/polyketide metabolism. Along with oxidoreductases and translocases, isomerases are one of the classes with the highest ratio of paralogous enzymes. Due to its relatively small number of members, it is plausible to explore it in greater detail to identify specific cases of gene duplication. Here, we present an analysis at the level of individual isomerases and identify different members that seem to be involved in duplication events in prokaryotes. As was suggested in a previous study, there is no homogeneous distribution of paralogs, but rather they accumulate into a few subcategories, some of which differ between Archaea and Bacteria. As expected, the metabolic processes with more paralogous isomerases have to do with carbohydrate metabolism but also with RNA modification (a particular case involving an rRNA-modifying isomerase is thoroughly discussed and analyzed in detail). Overall, our findings suggest that the most common fate for paralogous enzymes is the retention of the original enzymatic function, either associated with a dosage effect or with differential expression in response to changing environments, followed by subfunctionalization and, to a much lesser degree, neofunctionalization, which is consistent with what has been reported elsewhere.

Unexpected absence of ribosomal protein genes from metagenome-assembled genomes

Metagenome-assembled genomes (MAGs) have revealed the hidden diversity and functions of uncultivated microbes, but their reconstruction from metagenomes remains a computationally difficult task. Repetitive or exogenous sequences, such as ribosomal RNA and horizontally transferred genes, are frequently absent from MAGs because of misassembly and binning errors. Here, we report that ribosomal protein genes are also often absent from MAGs, although they are neither repetitive nor exogenous. Comprehensive analyses of more than 190,000 MAGs revealed that these genes could be missing in more than 20-40% of near-complete (i.e., with completeness of 90% or higher) MAGs. While some uncultivated environmental microbes intrinsically lack some ribosomal protein genes, we found that this unexpected absence is largely due to special evolutionary patterns of codon usage bias in ribosomal protein genes and algorithmic characteristics of metagenomic binning, which is dependent on tetranucleotide frequencies of contigs. This problem reflects the microbial life-history strategy. Fast-growing microbes tend to have this difficulty, likely because of strong evolutionary pressures on ribosomal protein genes toward the efficient assembly of ribosomes. Our observations caution those who study genomics and phylogeny of uncultivated microbes, the diversity and evolution of microbial genes in the central dogma, and bioinformatics in metagenomics.

Plant metabolic gene clusters in the multi-omics era

Secondary metabolism in plants gives rise to a vast array of small-molecule natural products. The discovery of operon-like gene clusters in plants has provided a new perspective on the evolution of specialized metabolism and the opportunity to rapidly advance the metabolic engineering of natural product production. Here, we review historical aspects of the study of plant metabolic gene clusters as well as general strategies for identifying plant metabolic gene clusters in the multi-omics era. We also emphasize the exploration of their natural variation and evolution, as well as new strategies for the prospecting of plant metabolic gene clusters and a deeper understanding of how their structure influences their function.

Horizontal Gene Transfer of Genes Encoding Copper-Containing Membrane-Bound Monooxygenase (CuMMO) and Soluble Di-iron Monooxygenase (SDIMO) in Ethane- and Propane-Oxidizing Rhodococcus Bacteria

The families of copper-containing membrane-bound monooxygenases (CuMMOs) and soluble di-iron monooxygenases (SDIMOs) are involved not only in methane oxidation but also in short-chain alkane oxidation. Here, we describe Rhodococcus sp. strain ZPP, a bacterium able to grow with ethane or propane as the sole carbon and energy source, and report on the horizontal gene transfer (HGT) of actinobacterial hydrocarbon monooxygenases (HMOs) of the CuMMO family and the sMMO (soluble methane monooxygenase)-like SDIMO in the genus Rhodococcus. The key function of HMO in strain ZPP for propane oxidation was verified by allylthiourea inhibition. The HMO genes (designated hmoCAB) and those encoding sMMO-like SDIMO (designated smoXYB1C1Z) are located on a linear megaplasmid (pRZP1) of strain ZPP. Comparative genomic analysis of similar plasmids indicated the mobility of these plasmids within the genus Rhodococcus. The plasmid pRZP1 in strain ZPP could be conjugatively transferred to a recipient Rhodococcus erythropolis strain in a mating experiment and showed similar ethane- and propane-consuming activities. Finally, our findings demonstrate that the horizontal transfer of plasmid-based CuMMO and SDIMO genes confers the ability to use ethane and propane on the recipient. IMPORTANCE CuMMOs and SDIMOs initiate the aerobic oxidation of alkanes in bacteria. Here, the supposition that horizontally transferred plasmid-based CuMMO and SDIMO genes confer on the recipient similar abilities to use ethane and propane was proposed and confirmed in Rhodococcus. This study is a living example of HGT of CuMMOs and SDIMOs and outlines the plasmid-borne properties responsible for gaseous alkane degradation. Our results indicate that plasmids can support the rapid evolution of enzyme-mediated biogeochemical processes.

Challenges and Perspective in Integrated Multi-Omics in Gut Microbiota Studies

The advent of omic technology has made it possible to identify viable but unculturable micro-organisms in the gut. Therefore, application of multi-omic technologies in gut microbiome studies has become invaluable for unveiling a comprehensive interaction between these commensals in health and disease. Meanwhile, despite the successful identification of many microbial and host-microbial cometabolites that have been reported so far, it remains difficult to clearly identify the origin and function of some proteins and metabolites that are detected in gut samples. However, the application of single omic techniques for studying the gut microbiome comes with its own challenges which may be overcome if a number of different omics techniques are combined. In this review, we discuss our current knowledge about multi-omic techniques, their challenges and future perspective in this field of gut microbiome studies.

Prediction and analysis of metagenomic operons via MetaRon: a pipeline for prediction of Metagenome and whole-genome opeRons

Background

Efficient regulation of bacterial genes in response to the environmental stimulus results in unique gene clusters known as operons. Lack of complete operonic reference and functional information makes the prediction of metagenomic operons a challenging task; thus, opening new perspectives on the interpretation of the host-microbe interactions.

Results

In this work, we identified whole-genome and metagenomic operons via MetaRon (Metagenome and whole-genome opeRon prediction pipeline). MetaRon identifies operons without any experimental or functional information. MetaRon was implemented on datasets with different levels of complexity and information. Starting from its application on whole-genome to simulated mixture of three whole-genomes (E. coli MG1655, Mycobacterium tuberculosis H37Rv and Bacillus subtilis str. 16), E. coli c20 draft genome extracted from chicken gut and finally on 145 whole-metagenome data samples from human gut. MetaRon consistently achieved high operon prediction sensitivity, specificity and accuracy across E. coli whole-genome (97.8, 94.1 and 92.4%), simulated genome (93.7, 75.5 and 88.1%) and E. coli c20 (87, 91 and 88%,), respectively. Finally, we identified 1,232,407 unique operons from 145 paired-end human gut metagenome samples. We also report strong association of type 2 diabetes with Maltose phosphorylase (K00691), 3-deoxy-D-glycero-D-galacto-nononate 9-phosphate synthase (K21279) and an uncharacterized protein (K07101).

Conclusion

With MetaRon, we were able to remove two notable limitations of existing whole-genome operon prediction methods: (1) generalizability (ability to predict operons in unrelated bacterial genomes), and (2) whole-genome and metagenomic data management. We also demonstrate the use of operons as a subset to represent the trends of secondary metabolites in whole-metagenome data and the role of secondary metabolites in the occurrence of disease condition. Using operonic data from metagenome to study secondary metabolic trends will significantly reduce the data volume to more precise data. Furthermore, the identification of metabolic pathways associated with the occurrence of type 2 diabetes (T2D) also presents another dimension of analyzing the human gut metagenome. Presumably, this study is the first organized effort to predict metagenomic operons and perform a detailed analysis in association with a disease, in this case type 2 diabetes. The application of MetaRon to metagenomic data at diverse scale will be beneficial to understand the gene regulation and therapeutic metagenomics.

The default cyanobacterial linked genome: an interactive platform based on cyanobacterial linkage networks to assist functional genomics

A database of cyanobacterial linked genomes that can be accessed through an interactive platform (https://dfgm.ua.es/genetica/investigacion/cyanobacterial_genetics/Resources.html) was generated on the bases of conservation of gene neighborhood across 124 cyanobacterial species. It allows flexible generation of gene networks at different threshold values. The default cyanobacterial linked genome, whose global properties are analyzed here, connects most of the cyanobacterial core genes. The potential of the web tool is discussed in relation to other bioinformatics approaches based on guilty-by-association principles, with selected examples of networks illustrating its usefulness for genes found exclusively in cyanobacteria or in cyanobacteria and chloroplasts. We believe that this tool will provide useful predictions that are readily testable in Synechococcus elongatus PCC7942 and other model organisms performing oxygenic photosynthesis.

Prediction and large-scale analysis of primary operons in plastids reveals unique genetic features in the evolution of chloroplasts

While bacterial operons have been thoroughly studied, few analyses of chloroplast operons exist, limiting the ability to study fundamental elements of these structures and utilize them for synthetic biology. Here, we describe the creation of a plastome-specific operon database (link provided below) achieved by combining experimental tools and predictive modeling. Using a Reverse-Transcription-PCR based method and published data, we determined the transcription-state of 213 gene pairs from four plastomes of evolutionary distinct organisms. By analyzing sequence-based features computed for our dataset, we were able to highlight fundamental characteristics differentiating between operon pairs and non-operon pairs. These include an interesting tendency toward maintaining similar messenger RNA-folding profiles in operon gene pairs, a feature that failed to yield any informative separation in cyanobacteria, suggesting that it catches unique traits of operon gene expression, which have evolved post-endosymbiosis. Subsequently, we used this feature set to train a random-forest classifier for operon prediction. As our results demonstrate the ability of our predictor to obtain accurate (84%) and robust predictions on unlabeled datasets, we proceeded to building operon maps for 2018 sequenced plastids. Our database may now present new opportunities for promoting metabolic engineering and synthetic biology in chloroplasts.

Cold plasma effect on the proteome of Pseudomonas aeruginosa - Role for bacterioferritin

Cold atmospheric-pressure plasma (CAP) is a relatively new method used for bacterial inactivation. CAP is ionized gas that can be generated by applying an electric current to air or a feeding gas. It contains reactive species and emits UV radiation, which have antibacterial activity. Previous data suggests that CAP is effective in microbial inactivation and can decontaminate and sterilize surfaces, but its exact mode of action is still under debate. This study demonstrates the effect of CAP on the whole proteome of Pseudomonas aeruginosa PAO1 biofilms, which is a dominant pathogen in cystic fibrosis and medical device-related infections. Liquid chromatography-mass spectrometry (LC-MS) was used to identify differentially regulated proteins of whole cell P. aeruginosa extracts. A total of 16 proteins were identified to be affected by plasma treatment compared to the control. Eight of the identified proteins have functions in transcription and translation and their expression changes are likely to be part of a general physiological response instead of a CAP-specific adaptation. However, CAP also affected bacterioferritin (Bfr), Isocitrate dehydrogenase (Idh), Trigger factor (Tig) and a chemotaxis protein, which may be involved in P. aeruginosa’s specific response to CAP. We confirm that bacterioferritin B plays a role in the bacterial response to CAP because ΔbfrB mutants of both PAO1 and PA14 are more susceptible to plasma-induced cell-death than their corresponding wild-type strains. To our knowledge, this is the first study showing the effect of plasma on the whole proteome of a pathogenic microorganism. It will help our understanding of the mode of action of CAP-mediated bacterial inactivation and thus support a safe and effective routine use of CAP in clinical and industrial settings.

The urea carboxylase and allophanate hydrolase activities of urea amidolyase are functionally independent

Urea amidolyase (UAL) is a multifunctional biotin-dependent enzyme that contributes to both bacterial and fungal pathogenicity by catalyzing the ATP-dependent cleavage of urea into ammonia and CO2 . UAL is comprised of two enzymatic components: urea carboxylase (UC) and allophanate hydrolase (AH). These enzyme activities are encoded on separate but proximally related genes in prokaryotes while, in most fungi, they are encoded by a single gene that produces a fusion enzyme on a single polypeptide chain. It is unclear whether the UC and AH activities are connected through substrate channeling or other forms of direct communication. Here, we use multiple biochemical approaches to demonstrate that there is no substrate channeling or interdomain/intersubunit communication between UC and AH. Neither stable nor transient interactions can be detected between prokaryotic UC and AH and the catalytic efficiencies of UC and AH are independent of one another. Furthermore, an artificial fusion of UC and AH does not significantly alter the AH enzyme activity or catalytic efficiency. These results support the surprising functional independence of AH from UC in both the prokaryotic and fungal UAL enzymes and serve as an important reminder that the evolution of multifunctional enzymes through gene fusion events does not always correlate with enhanced catalytic function.

The rnc Gene Promotes Exopolysaccharide Synthesis and Represses the vicRKX Gene Expressions via MicroRNA-Size Small RNAs in Streptococcus mutans

Dental caries is a biofilm-dependent disease that largely relies on the ability of Streptococcus mutans to synthesize exopolysaccharides. Although the rnc gene is suggested to be involved in virulence mechanisms in many other bacteria, the information regarding it in S. mutans is very limited. Here, using deletion or overexpression mutant assay, we demonstrated that rnc in S. mutans significantly positively regulated exopolysaccharide synthesis and further altered biofilm formation. Meanwhile, the cariogenecity of S. mutans was decreased by deletion of rnc in a specific pathogen-free (SPF) rat model. Interestingly, analyzing the expression at mRNA level, we found the downstream vic locus was repressed by rnc in S. mutans. Using deep sequencing and bioinformatics analysis, for the first time, three putative microRNA-size small RNAs (msRNAs) targeting vicRKX were predicted in S. mutans. The expression levels of these msRNAs were negatively correlated with vicRKX but positively correlated with rnc, indicating rnc probably repressed vicRKX expression through msRNAs at the post-transcriptional level. In all, the results present that rnc has a potential role in the regulation of exopolysaccharide synthesis and can affect vicRKX expressions via post-transcriptional repression in S. mutans. This study provides an alternative avenue for further research aimed at preventing caries.

A Transcriptome Map of Actinobacillus pleuropneumoniae at Single-Nucleotide Resolution Using Deep RNA-Seq

Actinobacillus pleuropneumoniae is the pathogen of porcine contagious pleuropneumoniae, a highly contagious respiratory disease of swine. Although the genome of A. pleuropneumoniae was sequenced several years ago, limited information is available on the genome-wide transcriptional analysis to accurately annotate the gene structures and regulatory elements. High-throughput RNA sequencing (RNA-seq) has been applied to study the transcriptional landscape of bacteria, which can efficiently and accurately identify gene expression regions and unknown transcriptional units, especially small non-coding RNAs (sRNAs), UTRs and regulatory regions. The aim of this study is to comprehensively analyze the transcriptome of A. pleuropneumoniae by RNA-seq in order to improve the existing genome annotation and promote our understanding of A. pleuropneumoniae gene structures and RNA-based regulation. In this study, we utilized RNA-seq to construct a single nucleotide resolution transcriptome map of A. pleuropneumoniae. More than 3.8 million high-quality reads (average length ~90 bp) from a cDNA library were generated and aligned to the reference genome. We identified 32 open reading frames encoding novel proteins that were mis-annotated in the previous genome annotations. The start sites for 35 genes based on the current genome annotation were corrected. Furthermore, 51 sRNAs in the A. pleuropneumoniae genome were discovered, of which 40 sRNAs were never reported in previous studies. The transcriptome map also enabled visualization of 5'- and 3'-UTR regions, in which contained 11 sRNAs. In addition, 351 operons covering 1230 genes throughout the whole genome were identified. The RNA-Seq based transcriptome map validated annotated genes and corrected annotations of open reading frames in the genome, and led to the identification of many functional elements (e.g. regions encoding novel proteins, non-coding sRNAs and operon structures). The transcriptional units described in this study provide a foundation for future studies concerning the gene functions and the transcriptional regulatory architectures of this pathogen.

Modulation of Biofilm Exopolysaccharides by the Streptococcus mutans vicX Gene

The cariogenic pathogen Streptococcus mutans effectively utilizes dietary sucrose for the synthesis of exopolysaccharide, which act as a scaffold for its biofilm, thus contributing to its pathogenicity, environmental stress tolerance, and antimicrobial resistance. The two-component system VicRK of S. mutans regulates a group of virulence genes that are associated with biofilm matrix synthesis. Knockout of vicX affects biofilm formation, oxidative stress tolerance, and transformation of S. mutans. However, little is known regarding the vicX-modulated structural characteristics of the exopolysaccharides underlying the biofilm formation and the phenotypes of the vicX mutants. Here, we identified the role of vicX in the structural characteristics of the exopolysaccharide matrix and biofilm physiology. The vicX mutant (SmuvicX) biofilms seemingly exhibited "desertification" with architecturally impaired exopolysaccharide-enmeshed cell clusters, compared with the UA159 strain (S. mutans wild type strain). Concomitantly, SmuvicX showed a decrease in water-insoluble glucan (WIG) synthesis and in WIG/water-soluble glucan (WSG) ratio. Gel permeation chromatography (GPC) showed that the WIG isolated from the SmuvicX biofilms had a much lower molecular weight compared with the UA159 strain indicating differences in polysaccharide chain lengths. A monosaccharide composition analysis demonstrated the importance of the vicX gene in the glucose metabolism. We performed metabolite profiling via (1)H nuclear magnetic resonance spectroscopy, which showed that several chemical shifts were absent in both WSG and WIG of SmuvicX biofilms compared with the UA159 strain. Thus, the modulation of structural characteristics of exopolysaccharide by vicX provides new insights into the interaction between the exopolysaccharide structure, gene functions, and cariogenicity. Our results suggest that vicX gene modulates the structural characteristics of exopolysaccharide associated with cariogenicity, which may be explored as a potential target that contributes to dental caries management. Furthermore, the methods used to purify the EPS of S. mutans biofilms and to analyze multiple aspects of its structure (GPC, gas chromatography-mass spectrometry, and (1)H nuclear magnetic resonance spectroscopy) may be useful approaches to determine the roles of other virulence genes for dental caries prevention.

Functional Operons in Secondary Metabolic Gene Clusters in Glarea lozoyensis (Fungi, Ascomycota, Leotiomycetes)

Operons are multigene transcriptional units which occur mostly in prokaryotes but rarely in eukaryotes. Protein-coding operons have not been reported in the Fungi even though they represent a very diverse kingdom of organisms. Here, we report a functional operon involved in the secondary metabolism of the fungus Glarea lozoyensis belonging to Leotiomycetes (Ascomycota). Two contiguous genes, glpks3 and glnrps7, encoding polyketide synthase and nonribosomal peptide synthetase, respectively, are cotranscribed into one dicistronic mRNA under the control of the same promoter, and the mRNA is then translated into two individual proteins, GLPKS3 and GLNRPS7. Heterologous expression in Aspergillus nidulans shows that the GLPKS3-GLNRPS7 enzyme complex catalyzes the biosynthesis of a novel pyrrolidinedione-containing compound, xenolozoyenone (compound 1), which indicates the operon is functional. Although it is structurally similar to prokaryotic operons, the glpks3-glnrps7 operon locus has a monophylogenic origin from fungi rather than having been horizontally transferred from prokaryotes. Moreover, two additional operons, glpks28-glnrps8 and glpks29-glnrps9, were verified at the transcriptional level in the same fungus. This is the first report of protein-coding operons in a member of the Fungi.

Operons are multigene transcriptional units which occur mostly in prokaryotes but rarely in eukaryotes. Three operon-like gene structures for secondary metabolism that were discovered in the filamentous fungus Glarea lozoyensis are the first examples of protein-coding operons identified in a member of the Fungi. Among them, the glpks3-glnrps7 operon is responsible for the biosynthesis of xenolozoyenone, which is a novel tetramic acid-containing compound. Although structurally similar to prokaryotic operons, the glpks3-glnrps7 operon locus did not result from horizontal gene transfer from prokaryotes. In addition, operonlike structures have been predicted in silico to be common in other fungi. The common occurrence and operonlike structure in fungi provide evolutionary insight and essential data for eukaryotic gene transcription.

Missing genes, multiple ORFs, and C-to-U type RNA editing in Acrasis kona (Heterolobosea, Excavata) mitochondrial DNA

Discoba (Excavata) is an ancient group of eukaryotes with great morphological and ecological diversity. Unlike the other major divisions of Discoba (Jakobida and Euglenozoa), little is known about the mitochondrial DNAs (mtDNAs) of Heterolobosea. We have assembled a complete mtDNA genome from the aggregating heterolobosean amoeba, Acrasis kona, which consists of a single circular highly AT-rich (83.3%) molecule of 51.5 kb. Unexpectedly, A. kona mtDNA is missing roughly 40% of the protein-coding genes and nearly half of the transfer RNAs found in the only other sequenced heterolobosean mtDNAs, those of Naegleria spp. Instead, over a quarter of A. kona mtDNA consists of novel open reading frames. Eleven of the 16 protein-coding genes missing from A. kona mtDNA were identified in its nuclear DNA and polyA RNA, and phylogenetic analyses indicate that at least 10 of these 11 putative nuclear-encoded mitochondrial (NcMt) proteins arose by direct transfer from the mitochondrion. Acrasis kona mtDNA also employs C-to-U type RNA editing, and 12 homologs of DYW-type pentatricopeptide repeat (PPR) proteins implicated in plant organellar RNA editing are found in A. kona nuclear DNA. A mapping of mitochondrial gene content onto a consensus phylogeny reveals a sporadic pattern of relative stasis and rampant gene loss in Discoba. Rampant loss occurred independently in the unique common lineage leading to Heterolobosea + Tsukubamonadida and later in the unique lineage leading to Acrasis. Meanwhile, mtDNA gene content appears to be remarkably stable in the Acrasis sister lineage leading to Naegleria and in their distant relatives Jakobida.

[Current status of theoretical studies on essential genes in microbes]

Essential genes are indispensable for the survival of an organism in optimal conditions. Recently, study on essential gene is becoming a hot topic of microbiology, genomics, and bioinformatics. This paper described the experiments that determined essential genes in some microbes and the theoretical researches on essential genes were reviewed. The major content contained comparison of essential genes and non-essential genes based on information on evolutionary conservation and sequence composition, and in silico prediction of essential genes, and analysis of the chromosomal distributions of essential genes. Finally, related progresses were concluded and the open problems were pointed out.

Genomics of the proteorhodopsin-containing marine flavobacterium Dokdonia sp. strain MED134

Proteorhodopsin phototrophy is expected to have considerable impact on the ecology and biogeochemical roles of marine bacteria. However, the genetic features contributing to the success of proteorhodopsin-containing bacteria remain largely unknown. We investigated the genome of Dokdonia sp. strain MED134 (Bacteroidetes) for features potentially explaining its ability to grow better in light than darkness. MED134 has a relatively high number of peptidases, suggesting that amino acids are the main carbon and nitrogen sources. In addition, MED134 shares with other environmental genomes a reduction in gene copies at the expense of important ones, like membrane transporters, which might be compensated by the presence of the proteorhodopsin gene. The genome analyses suggest Dokdonia sp. MED134 is able to respond to light at least partly due to the presence of a strong flavobacterial consensus promoter sequence for the proteorhodopsin gene. Moreover, Dokdonia sp. MED134 has a complete set of anaplerotic enzymes likely to play a role in the adaptation of the carbon anabolism to the different sources of energy it can use, including light or various organic matter compounds. In addition to promoting growth, proteorhodopsin phototrophy could provide energy for the degradation of complex or recalcitrant organic matter, survival during periods of low nutrients, or uptake of amino acids and peptides at low concentrations. Our analysis suggests that the ability to harness light potentially makes MED134 less dependent on the amount and quality of organic matter or other nutrients. The genomic features reported here may well be among the keys to a successful photoheterotrophic lifestyle.

Parallel evolution of transcriptome architecture during genome reorganization

Assembly of genes into operons is generally viewed as an important process during the continual adaptation of microbes to changing environmental challenges. However, the genome reorganization events that drive this process are also the roots of instability for existing operons. We have determined that there exists a statistically significant trend that correlates the proportion of genes encoded in operons in archaea to their phylogenetic lineage. We have further characterized how microbes deal with operon instability by mapping and comparing transcriptome architectures of four phylogenetically diverse extremophiles that span the range of operon stabilities observed across archaeal lineages: a photoheterotrophic halophile (Halobacterium salinarum NRC-1), a hydrogenotrophic methanogen (Methanococcus maripaludis S2), an acidophilic and aerobic thermophile (Sulfolobus solfataricus P2), and an anaerobic hyperthermophile (Pyrococcus furiosus DSM 3638). We demonstrate how the evolution of transcriptional elements (promoters and terminators) generates new operons, restores the coordinated regulation of translocated, inverted, and newly acquired genes, and introduces completely novel regulation for even some of the most conserved operonic genes such as those encoding subunits of the ribosome. The inverse correlation (r=-0.92) between the proportion of operons with such internally located transcriptional elements and the fraction of conserved operons in each of the four archaea reveals an unprecedented view into varying stages of operon evolution. Importantly, our integrated analysis has revealed that organisms adapted to higher growth temperatures have lower tolerance for genome reorganization events that disrupt operon structures.

Modeling and comparing the organization of circular genomes

Motivation: Most prokaryotic genomes are circular with a single chromosome (called circular genomes), which consist of bacteria and archaea. Orthologous genes (abbreviated as orthologs) are genes directly evolved from an ancestor gene, and can be traced through different species in evolution. Shared orthologs between bacterial genomes have been used to measure their genome evolution. Here, organization of circular genomes is analyzed via distributions of shared orthologs between genomes. However, these distributions are often asymmetric and bimodal; to date, there is no joint distribution to model such data. This motivated us to develop a family of bivariate distributions with generalized von Mises marginals (BGVM) and its statistical inference.

Results: A new measure based on circular grade correlation and the fraction of shared orthologs is proposed for association between circular genomes, and a visualization tool developed to depict genome structure similarity. The proposed procedures are applied to eight pairs of prokaryotes separated from domain down to species, and 13 mycoplasma bacteria that are mammalian pathogens belonging to the same genus. We close with remarks on further applications to many features of genomic organization, e.g. shared transcription factor binding sites, between any pair of circular genomes. Thus, the proposed procedures may be applied to identifying conserved chromosome backbones, among others, for genome construction in synthetic biology.

Availability: All codes of the BGVM procedures and 1000+ prokaryotic genomes are available at http://www.stat.sinica.edu.tw/∼gshieh/bgvm.htm.

Contact:gshieh@stat.sinica.edu.tw

Supplementary information:Supplementary data are available at Bioinformatics online.

Gene duplications in prokaryotes can be associated with environmental adaptation

Background

Gene duplication is a normal evolutionary process. If there is no selective advantage in keeping the duplicated gene, it is usually reduced to a pseudogene and disappears from the genome. However, some paralogs are retained. These gene products are likely to be beneficial to the organism, e.g. in adaptation to new environmental conditions. The aim of our analysis is to investigate the properties of paralog-forming genes in prokaryotes, and to analyse the role of these retained paralogs by relating gene properties to life style of the corresponding prokaryotes.

Results

Paralogs were identified in a number of prokaryotes, and these paralogs were compared to singletons of persistent orthologs based on functional classification. This showed that the paralogs were associated with for example energy production, cell motility, ion transport, and defence mechanisms. A statistical overrepresentation analysis of gene and protein annotations was based on paralogs of the 200 prokaryotes with the highest fraction of paralog-forming genes. Biclustering of overrepresented gene ontology terms versus species was used to identify clusters of properties associated with clusters of species. The clusters were classified using similarity scores on properties and species to identify interesting clusters, and a subset of clusters were analysed by comparison to literature data. This analysis showed that paralogs often are associated with properties that are important for survival and proliferation of the specific organisms. This includes processes like ion transport, locomotion, chemotaxis and photosynthesis. However, the analysis also showed that the gene ontology terms sometimes were too general, imprecise or even misleading for automatic analysis.

Conclusions

Properties described by gene ontology terms identified in the overrepresentation analysis are often consistent with individual prokaryote lifestyles and are likely to give a competitive advantage to the organism. Paralogs and singletons dominate different categories of functional classification, where paralogs in particular seem to be associated with processes involving interaction with the environment.