Microeconomic principles explain an optimal genome size in bacteria

Exploring the enzymatic repertoires of Bacteria and Archaea and their associations with metabolic maps

The evolution, survival, and adaptation of microbes are consequences of gene duplication, acquisition, and divergence in response to environmental challenges. In this context, enzymes play a central role in the evolution of organisms, because they are fundamental in cell metabolism. Here, we analyzed the enzymatic repertoire in 6,467 microbial genomes, including their abundances, and their associations with metabolic maps. We found that the enzymes follow a power-law distribution, in relation to the genome sizes. Therefore, we evaluated the total proportion enzymatic classes in relation to the genomes, identifying a descending-order proportion: transferases (EC:2.-), hydrolases (EC:3.-), oxidoreductases (EC:1.-), ligases (EC:6.-), lyases (EC:4.-), isomerases (EC:5.-), and translocases (EC:7-.). In addition, we identified a preferential use of enzymatic classes in metabolism pathways for xenobiotics, cofactors and vitamins, carbohydrates, amino acids, glycans, and energy. Therefore, this analysis provides clues about the functional constraints associated with the enzymatic repertoire of functions in Bacteria and Archaea.

Genomic analyses reveal a low-temperature adapted clade in Halorubrum, a widespread haloarchaeon across global hypersaline environments

BACKGROUND

Cold-adapted archaea have diverse ecological roles in a wide range of low-temperature environments. Improving our knowledge of the genomic features that enable psychrophiles to grow in cold environments helps us to understand their adaptive responses. However, samples from typical cold regions such as the remote Arctic and Antarctic are rare, and the limited number of high-quality genomes available leaves us with little data on genomic traits that are statistically associated with cold environmental conditions.

RESULTS

In this study, we examined the haloarchaeal genus Halorubrum and defined a new clade that represents six isolates from polar and deep earth environments ('PD group' hereafter). The genomic G + C content and amino acid composition of this group distinguishes it from other Halorubrum and the trends are consistent with the established genomic optimization of psychrophiles. The cold adaptation of the PD group was further supported by observations of increased flexibility of proteins encoded across the genome and the findings of a growth test.

CONCLUSIONS

The PD group Halorubrum exhibited denser genome packing, which confers higher metabolic potential with constant genome size, relative to the reference group, resulting in significant differences in carbon, nitrogen and sulfur metabolic patterns. The most marked feature was the enrichment of genes involved in sulfur cycling, especially the production of sulfite from organic sulfur-containing compounds. Our study provides an updated view of the genomic traits and metabolic potential of Halorubrum and expands the range of sources of cold-adapted haloarchaea.

A Reduction of Transcriptional Regulation in Aquatic Oligotrophic Microorganisms Enhances Fitness in Nutrient-Poor Environments

In this review, we consider the regulatory strategies of aquatic oligotrophs, microbial cells that are adapted to thrive under low-nutrient concentrations in oceans, lakes, and other aquatic ecosystems. Many reports have concluded that oligotrophs use less transcriptional regulation than copiotrophic cells, which are adapted to high nutrient concentrations and are far more common subjects for laboratory investigations of regulation. It is theorized that oligotrophs have retained alternate mechanisms of regulation, such as riboswitches, that provide shorter response times and smaller amplitude responses and require fewer cellular resources. We examine the accumulated evidence for distinctive regulatory strategies in oligotrophs. We explore differences in the selective pressures copiotrophs and oligotrophs encounter and ask why, although evolutionary history gives copiotrophs and oligotrophs access to the same regulatory mechanisms, they might exhibit distinctly different patterns in how these mechanisms are used. We discuss the implications of these findings for understanding broad patterns in the evolution of microbial regulatory networks and their relationships to environmental niche and life history strategy. We ask whether these observations, which have emerged from a decade of increased investigation of the cell biology of oligotrophs, might be relevant to recent discoveries of many microbial cell lineages in nature that share with oligotrophs the property of reduced genome size.

A panoramic view of the genomic landscape of the genus Streptomyces

We delineate the evolutionary plasticity of the ecologically and biotechnologically important genus Streptomyces, by analysing the genomes of 213 species. Streptomycetes genomes demonstrate high levels of internal homology, whereas the genome of their last common ancestor was already complex. Importantly, we identify the species-specific fingerprint proteins that characterize each species. Even among closely related species, we observed high interspecies variability of chromosomal protein-coding genes, species-level core genes, accessory genes and fingerprints. Notably, secondary metabolite biosynthetic gene clusters (smBGCs), carbohydrate-active enzymes (CAZymes) and protein-coding genes bearing the rare TTA codon demonstrate high intraspecies and interspecies variability, which emphasizes the need for strain-specific genomic mining. Highly conserved genes, such as those specifying genus-level core proteins, tend to occur in the central region of the chromosome, whereas those encoding proteins with evolutionarily volatile species-level fingerprints, smBGCs, CAZymes and TTA-codon-bearing genes are often found towards the ends of the linear chromosome. Thus, the chromosomal arms emerge as the part of the genome that is mainly responsible for rapid adaptation at the species and strain level. Finally, we observed a moderate, but statistically significant, correlation between the total number of CAZymes and three categories of smBGCs (siderophores, e-Polylysin and type III lanthipeptides) that are related to competition among bacteria.

Abiotic selection of microbial genome size in the global ocean

Strong purifying selection is considered a major evolutionary force behind small microbial genomes in the resource-poor photic ocean. However, very little is currently known about how the size of prokaryotic genomes evolves in the global ocean and whether patterns reflect shifts in resource availability in the epipelagic and relatively stable deep-sea environmental conditions. Using 364 marine microbial metagenomes, we investigate how the average genome size of uncultured planktonic prokaryotes varies across the tropical and polar oceans to the hadal realm. We find that genome size is highest in the perennially cold polar ocean, reflecting elongation of coding genes and gene dosage effects due to duplications in the interior ocean microbiome. Moreover, the rate of change in genome size due to temperature is 16-fold higher than with depth up to 200 m. Our results demonstrate how environmental factors can influence marine microbial genome size selection and ecological strategies of the microbiome.

Ecological insights into soil health according to the genomic traits and environment-wide associations of bacteria in agricultural soils

Soil microbiomes are sensitive to current and previous soil conditions, and bacterial 'bioindicators' of biological, physical, and chemical soil properties have considerable potential for soil health assessment. However, the lack of ecological or physiological information for most soil microorganisms limits our ability to interpret the associations of bioindicators and, thus, their utility for guiding management. We identified bioindicators of tillage intensity and twelve soil properties used to rate soil health using a 16S rRNA gene-based survey of farmland across North America. We then inferred the genomic traits of bioindicators and evaluated their environment-wide associations (EWAS) with respect to agricultural management practice, disturbance, and plant associations with 89 studies from agroecosystems. Most bioindicators were either positively correlated with biological properties (e.g., organic matter) or negatively correlated with physical and chemical properties. Higher soil health ratings corresponded with smaller genome size and higher coding density, while lower ratings corresponded with larger genomes and higher rrn copy number. Community-weighted genome size explained most variation in health ratings. EWAS linked prominent bioindicators with the impacts of environmental disturbances. Our findings provide ecological insights into bioindicators of soil properties relevant to soil health management, illustrating the tight coupling of microbiome and soil function.

Ecological selection of bacterial taxa with larger genome sizes in response to polycyclic aromatic hydrocarbons stress

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous priority pollutants that cause great damage to the natural environment and health. Average genome size in a community is critical for shedding light on microbiome's functional response to pollution stress within an environment. Here, microcosms under different concentrations were performed to evaluate the selection of PAHs stress on the average genome size in a community. We found the distinct communities of significantly larger genome size with the increase of PAHs concentration gradients in soils, and consistent trends were discovered in soils at different latitudes. The abundance of Proteobacteria and Deinococcus-Thermus with relatively larger genomes increased along with PAHs stress and well adapted to polluted environments. In contrast, the abundance of Patescibacteria with a highly streamlined and smaller genome decreased, implying complex interactions between environmental selection and functional fitness resulted in bacteria with larger genomes becoming more abundant. Moreover, we confirmed the increased capacity for horizontal transfer of degrading genes between communities by showing an increased connection number per node positively related to the nidA gene along the concentration gradients in the co-occurrence network. Our findings suggest PAHs tend to select bacterial taxa with larger genome sizes, with significant consequences for community stability and potential biodegradation strategies.

Description and functional testing of four species of the novel phototrophic genus Chioneia gen. nov., isolated from different East Antarctic environments

Seven Gram-negative, aerobic, non-sporulating, motile strains were isolated from terrestrial (R-67880^T, R-67883, R-36501 and R-36677^T) and aquatic (R-39604, R-39161^T and R-39594^T) East Antarctic environments (i.e. soil and aquatic microbial mats), between 2007 and 2014. Analysis of near-complete 16S rRNA gene sequences revealed that the strains potentially form a novel genus in the family Sphingomonadaceae (Alphaproteobacteria). DNA-DNA reassociation and average nucleotide identity values indicated distinction from close neighbors in the family Sphingomonadaceae and showed that the seven isolates form four different species. The main central pathways present in the strains are the glycolysis, tricarboxylic acid cycle and pentose phosphate pathway. The strains can use only a limited number of carbon sources and mainly depend on ammonia and sulfate as a nitrogen and sulfur source, respectively. The novel strains showed the potential of aerobic anoxygenic phototrophy, based on the presence of bacteriochlorophyll a pigments, which was corroborated by the presence of genes for all building blocks for a type 2 photosynthetic reaction center in the annotated genomes. Based on the results of phenotypic, genomic, phylogenetic and chemotaxonomic analyses, the strains could be assigned four new species in the novel genus Chioneia gen. nov. in the family Sphingomonadaceae, for which the names C. frigida sp. nov. (R-67880^T, R-67883 and R-36501), C. hiemis sp. nov. (R-36677^T), C. brumae sp. nov. (R-39161^T and R-39604) and C. algoris sp. nov. (R-39594^T) are proposed.

The Ascidia Ciona robusta Provides Novel Insights on the Evolution of the AP-1 Transcriptional Complex

The Activator Protein-1 transcription factor family (AP-1) transcriptional complex is historically defined as an early response group of transcription factors formed by dimeric complexes of the Jun, Fos, Atf, and Maf bZIP proteins that control cell proliferation and differentiation by regulating gene expression. It has been greatly investigated in many model organisms across metazoan evolution. Nevertheless, its complexity and variability of action made its multiple functions difficult to be defined. Here, we place the foundations for understanding the complexity of AP-1 transcriptional members in tunicates. We investigated the gene members of this family in the ascidian Ciona robusta and identified single copies of Jun, Fos, Atf3, Atf2/7, and Maf bZIP-related factors that could have a role in the formation of the AP-1 complex. We highlight that mesenchyme is a common cellular population where all these factors are expressed during embryonic development, and that, moreover, Fos shows a wider pattern of expression including also notochord and neural cells. By ectopic expression in transgenic embryos of Jun and Fos genes alone or in combination, we investigated the phenotypic alterations induced by these factors and highlighted a degree of functional conservation of the AP-1 complex between Ciona and vertebrates. The lack of gene redundancy and the first pieces of evidence of conserved functions in the control of cell movements and structural organization exerted by these factors open the way for using Ciona as a helpful model system to uncover the multiple potentialities of this highly complex family of bZIP transcription factors.

Picking the right metaphors for addressing microbial systems: economic theory helps understanding biological complexity

Any descriptive language is necessarily metaphoric and interpretative. Two somewhat overlapping-but not identical-languages have been thoroughly employed in the last decade to address the issue of regulatory complexity in biological systems: the terminology of network theory and the jargon of electric circuitry. These approaches have found many formal equivalences between the layout of extant genetic circuits and the architecture of man-made counterparts. However, these languages still fail to describe accurately key features of biological objects, in particular the diversity of signal-transfer molecules and the diffusion that is inherent to any biochemical system. Furthermore, current formalisms associated with networks and circuits can hardly face the problem of multi-scale regulatory complexity-from single molecules to entire ecosystems. We argue that the language of economic theory might be instrumental not only to portray accurately many features of regulatory networks, but also to unveil aspects of the biological complexity problem that remain opaque to other types of analyses. The main perspective opened by the economic metaphor when applied to control of microbiological activities is a focus on metabolism, not gene selfishness, as the necessary background to make sense of regulatory phenomena. As an example, we analyse and reinterpret the widespread phenomenon of catabolite repression with the formal frame of the consumer's choice theory.

Ecological selection for small microbial genomes along a temperate-to-thermal soil gradient

Small bacterial and archaeal genomes provide insights into the minimal requirements for life¹ and are phylogenetically widespread². However, the precise environmental pressures that constrain genome size in free-living microorganisms are unknown. A study including isolates has shown that thermophiles and other bacteria with high optimum growth temperatures often have small genomes³. It is unclear whether this relationship extends generally to microorganisms in nature^4,5 and more specifically to microorganisms that inhabit complex and highly variable environments, such as soils^3,6,7. To understand the genomic traits of thermally adapted microorganisms, here we investigated metagenomes from a 45 °C gradient of temperate-to-thermal soils that lie over the ongoing Centralia, Pennsylvania (USA) coal-seam fire. We found that hot soils harboured distinct communities with small genomes and small cell sizes relative to those in ambient soils. Hot soils notably lacked genes that encode known two-component regulatory systems, and antimicrobial production and resistance genes. Our results provide field evidence for the inverse relationship between microbial genome size and temperature in a diverse, free-living community over a wide range of temperatures that support microbial life.

Abundance, diversity and domain architecture variability in prokaryotic DNA-binding transcription factors

Gene regulation at the transcriptional level is a central process in all organisms, and DNA-binding transcription factors, known as TFs, play a fundamental role. This class of proteins usually binds at specific DNA sequences, activating or repressing gene expression. In general, TFs are composed of two domains: the DNA-binding domain (DBD) and an extra domain, which in this work we have named “companion domain” (CD). This latter could be involved in one or more functions such as ligand binding, protein-protein interactions or even with enzymatic activity. In contrast to DBDs, which have been widely characterized both experimentally and bioinformatically, information on the abundance, distribution, variability and possible role of the CDs is scarce. Here, we investigated these issues associated with the domain architectures of TFs in prokaryotic genomes. To this end, 19 families of TFs in 761 non-redundant bacterial and archaeal genomes were evaluated. In this regard we found four main groups based on the abundance and distribution in the analyzed genomes: i) LysR and TetR/AcrR; ii) AraC/XylS, SinR, and others; iii) Lrp, Fis, ArsR, and others; and iv) a group that included only two families, ArgR and BirA. Based on a classification of the organisms according to the life-styles, a major abundance of regulatory families in free-living organisms, in contrast with pathogenic, extremophilic or intracellular organisms, was identified. Finally, the protein architecture diversity associated to the 19 families considering a weight score for domain promiscuity evidenced which regulatory families were characterized by either a large diversity of CDs, here named as “promiscuous” families given the elevated number of variable domains found in those TFs, or a low diversity of CDs. Altogether this information helped us to understand the diversity and distribution of the 19 Prokaryotes TF families. Moreover, initial steps were taken to comprehend the variability of the extra domain in those TFs, which eventually might assist in evolutionary and functional studies.

Family-specific scaling laws in bacterial genomes

Among several quantitative invariants found in evolutionary genomics, one of the most striking is the scaling of the overall abundance of proteins, or protein domains, sharing a specific functional annotation across genomes of given size. The size of these functional categories change, on average, as power-laws in the total number of protein-coding genes. Here, we show that such regularities are not restricted to the overall behavior of high-level functional categories, but also exist systematically at the level of single evolutionary families of protein domains. Specifically, the number of proteins within each family follows family-specific scaling laws with genome size. Functionally similar sets of families tend to follow similar scaling laws, but this is not always the case. To understand this systematically, we provide a comprehensive classification of families based on their scaling properties. Additionally, we develop a quantitative score for the heterogeneity of the scaling of families belonging to a given category or predefined group. Under the common reasonable assumption that selection is driven solely or mainly by biological function, these findings point to fine-tuned and interdependent functional roles of specific protein domains, beyond our current functional annotations. This analysis provides a deeper view on the links between evolutionary expansion of protein families and the functional constraints shaping the gene repertoire of bacterial genomes.

Following Rapoport's Rule: the geographic range and genome size of bacterial taxa decline at warmer latitudes

We sought to test whether stream bacterial communities conform to Rapoport's Rule, a pattern commonly observed for plants and animals whereby taxa exhibit decreased latitudinal range sizes closer to the equator. Using a DNA sequencing approach, we explored the biogeography of biofilm bacterial communities in 204 streams across a ∼1000 km latitudinal gradient. The range sizes of bacterial taxa were strongly correlated with latitude, decreasing closer to the equator, which coincided with a greater than fivefold increase in bacterial taxonomic richness. The relative richness and range size of bacteria were associated with spatially correlated variation in temperature and rainfall. These patterns were observed despite enormous variability in catchment environmental characteristics. Similar results were obtained when restricting the same analyses to native forest catchments, thereby controlling for spatial biases in land use. We analysed genomic data from ∼500 taxa detected in this study, for which data were available and found that bacterial communities at cooler latitudes also tended to possess greater potential metabolic potential. Collectively, these data provide the first evidence of latitudinal variation in the range size distributions of freshwater bacteria, a trend which may be determined, in part, by a trade-off between bacterial genome size and local variation in climatic conditions.

Letting go: bacterial genome reduction solves the dilemma of adapting to predation mortality in a substrate-restricted environment

Resource limitation and predation mortality are major determinants of microbial population dynamics, and optimization for either aspect is considered to imply a trade-off with respect to the other. Adaptation to these selective factors may, moreover, lead to disadvantages at rich growth conditions. We present an example of a concomitant evolutionary optimization to both, substrate limitation and predation in an aggregate-forming freshwater bacterial isolate, and we elucidate an underlying genomic mechanism. Bacteria were propagated in serial batch culture in a nutrient-restricted environment either with or without a bacterivorous flagellate. Strains isolated after 26 growth cycles of the predator-prey co-cultures formed as much total biomass as the ancestor at ancestral growth conditions, albeit largely reallocated to cell aggregates. A ~273 kbp genome fragment was lost in three strains that had independently evolved with predators. These strains had significantly higher growth yield on substrate-restricted media than others that were isolated from the same treatment before the excision event. Under predation pressure, the isolates with the deletion outcompeted both, the ancestor and the strains evolved without predators even at rich growth conditions. At the same time, genome reduction led to a growth disadvantage in the presence of benzoate due to the loss of the respective degradation pathway, suggesting that niche constriction might be the price for the bidirectional optimization.

Genomic and phenotypic insights into the ecology of Arthrobacter from Antarctic soils

BACKGROUND

Members of the bacterial genus Arthrobacter are both readily cultured and commonly identified in Antarctic soil communities. Currently, relatively little is known about the physiological traits that allow these bacteria to survive in the harsh Antarctic soil environment. The aim of this study is to investigate if Antarctic strains of Arthrobacter owe their resilience to substantial genomic changes compared to Arthrobacter spp. isolated from temperate soil environments.

RESULTS

Quantitative PCR-based analysis revealed that up to 4% of the soil bacterial communities were comprised of Arthrobacter spp. at four locations in the Ross Sea Region. Genome analysis of the seven Antarctic Arthrobacter isolates revealed several features that are commonly observed in psychrophilic/psychrotolerant bacteria. These include genes primarily associated with sigma factors, signal transduction pathways, the carotenoid biosynthesis pathway and genes induced by cold-shock, oxidative and osmotic stresses. However, these genes were also identified in genomes of seven temperate Arthrobacter spp., suggesting that these mechanisms are beneficial for growth and survival in a range of soil environments. Phenotypic characterisation revealed that Antarctic Arthrobacter isolates demonstrate significantly lower metabolic versatility and a narrower salinity tolerance range compared to temperate Arthrobacter species. Comparative analyses also revealed fewer protein-coding sequences and a significant decrease in genes associated with transcription and carbohydrate transport and metabolism in four of the seven Antarctic Arthrobacter isolates. Notwithstanding genome incompleteness, these differences together with the decreased metabolic versatility are indicative of genome content scaling.

CONCLUSIONS

The genomes of the seven Antarctic Arthrobacter isolates contained several features that may be beneficial for growth and survival in the Antarctic soil environment, although these features were not unique to the Antarctic isolates. These genome sequences allow further investigations into the expression of physiological traits that enable survival under extreme conditions and, more importantly, into the ability of these bacteria to respond to future perturbations including climate change and human impacts.

Genome-based comparative analyses of Antarctic and temperate species of Paenibacillus

Antarctic soils represent a unique environment characterised by extremes of temperature, salinity, elevated UV radiation, low nutrient and low water content. Despite the harshness of this environment, members of 15 bacterial phyla have been identified in soils of the Ross Sea Region (RSR). However, the survival mechanisms and ecological roles of these phyla are largely unknown. The aim of this study was to investigate whether strains of Paenibacillus darwinianus owe their resilience to substantial genomic changes. For this, genome-based comparative analyses were performed on three P. darwinianus strains, isolated from gamma-irradiated RSR soils, together with nine temperate, soil-dwelling Paenibacillus spp. The genome of each strain was sequenced to over 1,000-fold coverage, then assembled into contigs totalling approximately 3 Mbp per genome. Based on the occurrence of essential, single-copy genes, genome completeness was estimated at approximately 88%. Genome analysis revealed between 3,043-3,091 protein-coding sequences (CDSs), primarily associated with two-component systems, sigma factors, transporters, sporulation and genes induced by cold-shock, oxidative and osmotic stresses. These comparative analyses provide an insight into the metabolic potential of P. darwinianus, revealing potential adaptive mechanisms for survival in Antarctic soils. However, a large proportion of these mechanisms were also identified in temperate Paenibacillus spp., suggesting that these mechanisms are beneficial for growth and survival in a range of soil environments. These analyses have also revealed that the P. darwinianus genomes contain significantly fewer CDSs and have a lower paralogous content. Notwithstanding the incompleteness of the assemblies, the large differences in genome sizes, determined by the number of genes in paralogous clusters and the CDS content, are indicative of genome content scaling. Finally, these sequences are a resource for further investigations into the expression of physiological attributes that enable survival under extreme conditions and selection processes that affect prokaryotic genome evolution.

Reduced genome size of Helicobacter pylori originating from East Asia

Helicobacter pylori (H. pylori), a major pathogen colonizing the human stomach, shows great genetic variation. Comparative analysis of strains from different H. pylori populations revealed that the genome size of strains from East Asia decreased to 1.60 Mbp, which is significantly smaller than that from Europe or Africa. In parallel with the genome reduction, the number of protein coding genes was decreased, and the guanine-cytosine content was lowered to 38.9%. Elimination of non-essential genes by mutations is likely to be a major cause of the genome reduction. Bacteria with a small genome cost less energy. Thus, H. pylori strains from East Asia may have proliferation and growth advantages over those from Western countries. This could result in enhanced capacity of bacterial spreading. Therefore, the reduced genome size potentially contributes to the high prevalence of H. pylori in East Asia.

Implications of streamlining theory for microbial ecology

Whether a small cell, a small genome or a minimal set of chemical reactions with self-replicating properties, simplicity is beguiling. As Leonardo da Vinci reportedly said, 'simplicity is the ultimate sophistication'. Two diverging views of simplicity have emerged in accounts of symbiotic and commensal bacteria and cosmopolitan free-living bacteria with small genomes. The small genomes of obligate insect endosymbionts have been attributed to genetic drift caused by small effective population sizes (Ne). In contrast, streamlining theory attributes small cells and genomes to selection for efficient use of nutrients in populations where Ne is large and nutrients limit growth. Regardless of the cause of genome reduction, lost coding potential eventually dictates loss of function. Consequences of reductive evolution in streamlined organisms include atypical patterns of prototrophy and the absence of common regulatory systems, which have been linked to difficulty in culturing these cells. Recent evidence from metagenomics suggests that streamlining is commonplace, may broadly explain the phenomenon of the uncultured microbial majority, and might also explain the highly interdependent (connected) behavior of many microbial ecosystems. Streamlining theory is belied by the observation that many successful bacteria are large cells with complex genomes. To fully appreciate streamlining, we must look to the life histories and adaptive strategies of cells, which impose minimum requirements for complexity that vary with niche.

Growth temperature and genome size in bacteria are negatively correlated, suggesting genomic streamlining during thermal adaptation

Prokaryotic genomes are small and compact. Either this feature is caused by neutral evolution or by natural selection favoring small genomes—genome streamlining. Three separate prior lines of evidence argue against streamlining for most prokaryotes. We find that the same three lines of evidence argue for streamlining in the genomes of thermophile bacteria. Specifically, with increasing habitat temperature and decreasing genome size, the proportion of genomic DNA in intergenic regions decreases. Furthermore, with increasing habitat temperature, generation time decreases. Genome-wide selective constraints do not decrease as in the reduced genomes of host-associated species. Reduced habitat variability is not a likely explanation for the smaller genomes of thermophiles. Genome size may be an indirect target of selection due to its association with cell volume. We use metabolic modeling to demonstrate that known changes in cell structure and physiology at high temperature can provide a selective advantage to reduce cell volume at high temperatures.

Extraordinary expansion of a Sorangium cellulosum genome from an alkaline milieu

Complex environmental conditions can significantly affect bacterial genome size by unknown mechanisms. The So0157-2 strain of Sorangium cellulosum is an alkaline-adaptive epothilone producer that grows across a wide pH range. Here, we show that the genome of this strain is 14,782,125 base pairs, 1.75-megabases larger than the largest bacterial genome from S. cellulosum reported previously. The total 11,599 coding sequences (CDSs) include massive duplications and horizontally transferred genes, regulated by lots of protein kinases, sigma factors and related transcriptional regulation co-factors, providing the So0157-2 strain abundant resources and flexibility for ecological adaptation. The comparative transcriptomics approach, which detected 90.7% of the total CDSs, not only demonstrates complex expression patterns under varying environmental conditions but also suggests an alkaline-improved pathway of the insertion and duplication, which has been genetically testified, in this strain. These results provide insights into and a paradigm for how environmental conditions can affect bacterial genome expansion.

Increments and duplication events of enzymes and transcription factors influence metabolic and regulatory diversity in prokaryotes

In this work, the content of enzymes and DNA-binding transcription factors (TFs) in 794 non-redundant prokaryotic genomes was evaluated. The identification of enzymes was based on annotations deposited in the KEGG database as well as in databases of functional domains (COG and PFAM) and structural domains (Superfamily). For identifications of the TFs, hidden Markov profiles were constructed based on well-known transcriptional regulatory families. From these analyses, we obtained diverse and interesting results, such as the negative rate of incremental changes in the number of detected enzymes with respect to the genome size. On the contrary, for TFs the rate incremented as the complexity of genome increased. This inverse related performance shapes the diversity of metabolic and regulatory networks and impacts the availability of enzymes and TFs. Furthermore, the intersection of the derivatives between enzymes and TFs was identified at 9,659 genes, after this point, the regulatory complexity grows faster than metabolic complexity. In addition, TFs have a low number of duplications, in contrast to the apparent high number of duplications associated with enzymes. Despite the greater number of duplicated enzymes versus TFs, the increment by which duplicates appear is higher in TFs. A lower proportion of enzymes among archaeal genomes (22%) than in the bacterial ones (27%) was also found. This low proportion might be compensated by the interconnection between the metabolic pathways in Archaea. A similar proportion was also found for the archaeal TFs, for which the formation of regulatory complexes has been proposed. Finally, an enrichment of multifunctional enzymes in Bacteria, as a mechanism of ecological adaptation, was detected.

The repertoire of DNA-binding transcription factors in prokaryotes: functional and evolutionary lessons

The capabilities of organisms to contend with environmental changes depend on their genes and their ability to regulate their expression. DNA-binding transcription factors (TFs) play a central role in this process, because they regulate gene expression positively and/or negatively, depending on the operator context and ligand-binding status. In this review, we summarise recent findings regarding the function and evolution of TFs in prokaryotes. We consider the abundance of TFs in bacteria and archaea, the role of DNA-binding domains and their partner domains, and the effects of duplication events in the evolution of regulatory networks. Finally, a comprehensive picture for how regulatory networks have evolved in prokaryotes is provided.

Comparative genomics suggests differential deployment of linear and branched signaling across bacteria

A major mode of signal transduction in bacteria is the two-component system, which involves phosphorylation of an output-generating receiver protein by a signal-sensing histidine kinase. This differs from the more common one-component system--where both signal sensing and output generation are performed by the same protein--in the spatial separation of the two activities and the obligate need for post-translational modification (phosphorylation). Many described two-component systems involve a linear structure where a single kinase phosphorylates a cognate receiver. However, inherently branched network structures are being increasingly discovered, though their prevalence is unknown. Though the simpler one-component systems are more common than two-component systems, some organisms encode a disproportionately high number of the latter; though these organisms are generally described as having 'complex' lifestyles, no systematic description of their signaling networks has been proposed. Finally, the relative contributions of the two modes of signal transduction towards achieving an optimal regulatory cost for growth and survival in an environment remain poorly understood. Here we present a comparative genomics survey of ~165,000 regulatory proteins from ~850 prokaryotic genomes and suggest that organisms with elevated occurrence of two-component systems--which generally belong to phylogenetic classes with relatively poor representation in genomic databases--also code for more complex and branched two-component networks. Such branched signaling might compensate for the apparent paucity in the total number of regulatory proteins these organisms encode. Finally, such interconnected signaling networks might be more common than anticipated, indicating the pressing need for genome-scale experimental studies of signaling networks in many understudied phylogenetic groups of organisms.

Measuring the evolutionary rewiring of biological networks

We have accumulated a large amount of biological network data and expect even more to come. Soon, we anticipate being able to compare many different biological networks as we commonly do for molecular sequences. It has long been believed that many of these networks change, or "rewire", at different rates. It is therefore important to develop a framework to quantify the differences between networks in a unified fashion. We developed such a formalism based on analogy to simple models of sequence evolution, and used it to conduct a systematic study of network rewiring on all the currently available biological networks. We found that, similar to sequences, biological networks show a decreased rate of change at large time divergences, because of saturation in potential substitutions. However, different types of biological networks consistently rewire at different rates. Using comparative genomics and proteomics data, we found a consistent ordering of the rewiring rates: transcription regulatory, phosphorylation regulatory, genetic interaction, miRNA regulatory, protein interaction, and metabolic pathway network, from fast to slow. This ordering was found in all comparisons we did of matched networks between organisms. To gain further intuition on network rewiring, we compared our observed rewirings with those obtained from simulation. We also investigated how readily our formalism could be mapped to other network contexts; in particular, we showed how it could be applied to analyze changes in a range of "commonplace" networks such as family trees, co-authorships and linux-kernel function dependencies.

An overview of prokaryotic transcription factors : a summary of function and occurrence in bacterial genomes

Transcriptional initiation is arguably the most important control point for gene expression. It is regulated by a combination of factors, including DNA sequence and its three-dimensional topology, proteins and small molecules. In this chapter, we focus on the trans-acting factors of bacterial regulation. Initiation begins with the recruitment of the RNA polymerase holoenzyme to a specific locus upstream of the gene known as its promoter. The sigma factor, which is a component of the holoenzyme, provides the most fundamental mechanisms for orchestrating broad changes in gene expression state. It is responsible for promoter recognition as well as recruiting the holoenzyme to the promoter. Distinct sigma factors compete with for binding to a common pool of RNA polymerases, thus achieving condition-dependent differential expression. Another important class of bacterial regulators is transcription factors, which activate or repress transcription of target genes typically in response to an environmental or cellular trigger. These factors may be global or local depending on the number of genes and range of cellular functions that they target. The activities of both global and local transcription factors may be regulated either at a post-transcriptional level via signal-sensing protein domains or at the level of their own expression. In addition to modulating polymerase recruitment to promoters, several global factors are considered as "nucleoid-associated proteins" that impose structural constraints on the chromosome by altering the conformation of the bound DNA, thus influencing other processes involving DNA such as replication and recombination. This chapter concludes with a discussion of how regulatory interactions between transcription factors and their target genes can be represented as a network.

Genomic repertoires of DNA-binding transcription factors across the tree of life

Sequence-specific transcription factors (TFs) are important to genetic regulation in all organisms because they recognize and directly bind to regulatory regions on DNA. Here, we survey and summarize the TF resources available. We outline the organisms for which TF annotation is provided, and discuss the criteria and methods used to annotate TFs by different databases. By using genomic TF repertoires from ∼700 genomes across the tree of life, covering Bacteria, Archaea and Eukaryota, we review TF abundance with respect to the number of genes, as well as their structural complexity in diverse lineages. While typical eukaryotic TFs are longer than the average eukaryotic proteins, the inverse is true for prokaryotes. Only in eukaryotes does the same family of DNA-binding domain (DBD) occur multiple times within one polypeptide chain. This potentially increases the length and diversity of DNA-recognition sequence by reusing DBDs from the same family. We examined the increase in TF abundance with the number of genes in genomes, using the largest set of prokaryotic and eukaryotic genomes to date. As pointed out before, prokaryotic TFs increase faster than linearly. We further observe a similar relationship in eukaryotic genomes with a slower increase in TFs.

Translational selection is ubiquitous in prokaryotes

Codon usage bias in prokaryotic genomes is largely a consequence of background substitution patterns in DNA, but highly expressed genes may show a preference towards codons that enable more efficient and/or accurate translation. We introduce a novel approach based on supervised machine learning that detects effects of translational selection on genes, while controlling for local variation in nucleotide substitution patterns represented as sequence composition of intergenic DNA. A cornerstone of our method is a Random Forest classifier that outperformed previous distance measure-based approaches, such as the codon adaptation index, in the task of discerning the (highly expressed) ribosomal protein genes by their codon frequencies. Unlike previous reports, we show evidence that translational selection in prokaryotes is practically universal: in 460 of 461 examined microbial genomes, we find that a subset of genes shows a higher codon usage similarity to the ribosomal proteins than would be expected from the local sequence composition. These genes constitute a substantial part of the genome—between 5% and 33%, depending on genome size—while also exhibiting higher experimentally measured mRNA abundances and tending toward codons that match tRNA anticodons by canonical base pairing. Certain gene functional categories are generally enriched with, or depleted of codon-optimized genes, the trends of enrichment/depletion being conserved between Archaea and Bacteria. Prominent exceptions from these trends might indicate genes with alternative physiological roles; we speculate on specific examples related to detoxication of oxygen radicals and ammonia and to possible misannotations of asparaginyl–tRNA synthetases. Since the presence of codon optimizations on genes is a valid proxy for expression levels in fully sequenced genomes, we provide an example of an “adaptome” by highlighting gene functions with expression levels elevated specifically in thermophilic Bacteria and Archaea.

Synonymous codons are not equally common in genomes. The main causes of unequal codon usage are varying nucleotide substitution patterns, as manifested in the wide range of genomic nucleotide compositions. However, since the first E. coli and yeast genes were sequenced, it became evident that there was also a bias towards codons that can be translated to protein faster and more accurately. This bias was stronger in highly expressed genes, and its driving force was termed translational selection. Researchers sought for effects of translational selection in microbial genomes as they became available, employing a flurry of mathematical approaches which sometimes led to contradictory conclusions. We introduce a sensitive and accurate machine learning-based methodology and find that highly expressed genes have a recognizable codon usage pattern in almost every bacterial and archaeal genome analyzed, even after accounting for large differences in background nucleotide composition. We also show that the gene functional category has a great bearing on whether that gene is subject to translational selection. Since presence of codon optimizations can be used as a purely sequence-derived proxy for expression levels, we can delineate “adaptomes” by relating predicted gene activity to organisms' phenotypes, which we demonstrate on genomes of temperature-resistant Bacteria and Archaea.

On the need for widespread horizontal gene transfers under genome size constraint

Background

While eukaryotes primarily evolve by duplication-divergence expansion (and reduction) of their own gene repertoire with only rare horizontal gene transfers, prokaryotes appear to evolve under both gene duplications and widespread horizontal gene transfers over long evolutionary time scales. But, the evolutionary origin of this striking difference in the importance of horizontal gene transfers remains by and large a mystery.

Hypothesis

We propose that the abundance of horizontal gene transfers in free-living prokaryotes is a simple but necessary consequence of two opposite effects: i) their apparent genome size constraint compared to typical eukaryote genomes and ii) their underlying genome expansion dynamics through gene duplication-divergence evolution, as demonstrated by the presence of many tandem and block repeated genes. In principle, this combination of genome size constraint and underlying duplication expansion should lead to a coalescent-like process with extensive turnover of functional genes. This would, however, imply the unlikely, systematic reinvention of functions from discarded genes within independent phylogenetic lineages. Instead, we propose that the long-term evolutionary adaptation of free-living prokaryotes must have resulted in the emergence of efficient non-phylogenetic pathways to circumvent gene loss.

Implications

This need for widespread horizontal gene transfers due to genome size constraint implies, in particular, that prokaryotes must remain under strong selection pressure in order to maintain the long-term evolutionary adaptation of their "mutualized" gene pool, beyond the inevitable turnover of individual prokaryote species. By contrast, the absence of genome size constraint for typical eukaryotes has presumably relaxed their need for widespread horizontal gene transfers and strong selection pressure. Yet, the resulting loss of genetic functions, due to weak selection pressure and inefficient gene recovery mechanisms, must have ultimately favored the emergence of more complex life styles and ecological integration of many eukaryotes.

Reviewers

This article was reviewed by Pierre Pontarotti, Eugene V Koonin and Sergei Maslov.

Darwinian evolution in the light of genomics

Comparative genomics and systems biology offer unprecedented opportunities for testing central tenets of evolutionary biology formulated by Darwin in the Origin of Species in 1859 and expanded in the Modern Synthesis 100 years later. Evolutionary-genomic studies show that natural selection is only one of the forces that shape genome evolution and is not quantitatively dominant, whereas non-adaptive processes are much more prominent than previously suspected. Major contributions of horizontal gene transfer and diverse selfish genetic elements to genome evolution undermine the Tree of Life concept. An adequate depiction of evolution requires the more complex concept of a network or ‘forest’ of life. There is no consistent tendency of evolution towards increased genomic complexity, and when complexity increases, this appears to be a non-adaptive consequence of evolution under weak purifying selection rather than an adaptation. Several universals of genome evolution were discovered including the invariant distributions of evolutionary rates among orthologous genes from diverse genomes and of paralogous gene family sizes, and the negative correlation between gene expression level and sequence evolution rate. Simple, non-adaptive models of evolution explain some of these universals, suggesting that a new synthesis of evolutionary biology might become feasible in a not so remote future.

Genomics of bacteria and archaea: the emerging dynamic view of the prokaryotic world

The first bacterial genome was sequenced in 1995, and the first archaeal genome in 1996. Soon after these breakthroughs, an exponential rate of genome sequencing was established, with a doubling time of approximately 20 months for bacteria and approximately 34 months for archaea. Comparative analysis of the hundreds of sequenced bacterial and dozens of archaeal genomes leads to several generalizations on the principles of genome organization and evolution. A crucial finding that enables functional characterization of the sequenced genomes and evolutionary reconstruction is that the majority of archaeal and bacterial genes have conserved orthologs in other, often, distant organisms. However, comparative genomics also shows that horizontal gene transfer (HGT) is a dominant force of prokaryotic evolution, along with the loss of genetic material resulting in genome contraction. A crucial component of the prokaryotic world is the mobilome, the enormous collection of viruses, plasmids and other selfish elements, which are in constant exchange with more stable chromosomes and serve as HGT vehicles. Thus, the prokaryotic genome space is a tightly connected, although compartmentalized, network, a novel notion that undermines the ‘Tree of Life’ model of evolution and requires a new conceptual framework and tools for the study of prokaryotic evolution.

Principles of transcriptional regulation and evolution of the metabolic system in E. coli

Organisms must adapt to make optimal use of the metabolic system in response to environmental changes. In the long-term, this involves evolution of the genomic repertoire of enzymes; in the short-term, transcriptional control ensures that appropriate enzymes are expressed in response to transitory extracellular conditions. Unicellular organisms are particularly susceptible to environmental changes; however, genome-scale impact of these modulatory effects has not been explored so far in bacteria. Here, we integrate genome-scale data to investigate the evolutionary trends and transcriptional control of metabolism in Escherichia coli K12. Globally, the regulatory system is organized in a clear hierarchy of general and specific transcription factors (TFs) that control differing ranges of metabolic functions. Further, catabolic, anabolic, and central metabolic pathways are targeted by distinct combinations of these TFs. Locally, enzymes catalyzing sequential reactions in a metabolic pathway are co-regulated by the same TFs. Regulation is more complex at junctions: General TFs control the overall activity of all connecting reactions, whereas specific TFs control individual enzymes. Divergent junctions play a special role in delineating metabolic pathways and decouple the regulation of incoming and outgoing reactions. We find little evidence for differential usage of isozymes, which are generally co-expressed in similar conditions, and thus are likely to reinforce the metabolic system through redundancy. Finally, we show that enzymes controlled by the same TFs have a strong tendency to co-evolve, suggesting a significant constraint to maintain similar regulatory regimes during evolution. Catabolic, anabolic, and central energy pathways evolve differently, emphasizing the role of the environment in shaping the metabolic system. Many of the observations also occur in yeast, and our findings may apply across large evolutionary distances.

Do phosphorus requirements for RNA limit genome size in crustacean zooplankton?

As for most other organisms, genome size in zooplankton differs widely. This may have a range of consequences for growth rate, development, and life history strategies, yet the causes of this pronounced variability are not settled. Here we propose that small genome size may be an evolutionary consequence of phosphorus (P) allocation from DNA to RNA under P deficiency. To test this hypothesis we have compared the two major groups of zooplankton, copepods and cladocerans, that have overlapping niches and body size. Relative to the cladocerans, copepods have a more complex life history and a lower mass-specific P content, while cladocerans tend to have higher P and RNA contents and higher specific growth rates and frequently experience P-limited growth, likely due to a shortage of P for ribosome synthesis. Cladocerans also generally have smaller genomes than copepods (1C = 0.17–0.63 pg DNA·cell–1vs. 1C = 0.10–10 pg DNA·cell–1). Furthermore, cladocerans have a higher slope of the relationship of body size with DNA content (1.5 vs. 0.28 in copepods) and present almost 15-fold higher RNA:DNA ratios (24.8 in cladocerans vs. 1.6 in copepods). Hence, small genome size in cladocerans could reflect an evolutionary pressure towards “efficient” genomes to conserve a key element needed to maximize growth rate. We do not claim that this is a universal cause of genome size variability, but propose that streamlining of genomes could be related to P conservation rather than energy conservation. This could be relevant for a range of organisms that may suffer P-limited growth rates.

Large changes in regulome size herald the main prokaryotic lineages

Using a large-scale reconstruction of ancestral gene content, we show that radical changes in regulome size occur at the origins of major prokaryotic lineages. Subsequently, the duplication and deletion of regulators slows down in most lineages, except proteobacteria, significantly reducing the scaling of regulators and keeping their average proportion lineage-specific. Our results also suggest that major transitions in prokaryote evolution are related to changes in regulatory capacity rather than proteome innovations.

A protein interaction atlas for the nuclear receptors: properties and quality of a hub-based dimerisation network

BACKGROUND

The nuclear receptors are a large family of eukaryotic transcription factors that constitute major pharmacological targets. They exert their combinatorial control through homotypic heterodimerisation. Elucidation of this dimerisation network is vital in order to understand the complex dynamics and potential cross-talk involved.

RESULTS

Phylogeny, protein-protein interactions, protein-DNA interactions and gene expression data have been integrated to provide a comprehensive and up-to-date description of the topology and properties of the nuclear receptor interaction network in humans. We discriminate between DNA-binding and non-DNA-binding dimers, and provide a comprehensive interaction map, that identifies potential cross-talk between the various pathways of nuclear receptors.

CONCLUSION

We infer that the topology of this network is hub-based, and much more connected than previously thought. The hub-based topology of the network and the wide tissue expression pattern of NRs create a highly competitive environment for the common heterodimerising partners. Furthermore, a significant number of negative feedback loops is present, with the hub protein SHP [NR0B2] playing a major role. We also compare the evolution, topology and properties of the nuclear receptor network with the hub-based dimerisation network of the bHLH transcription factors in order to identify both unique themes and ubiquitous properties in gene regulation. In terms of methodology, we conclude that such a comprehensive picture can only be assembled by semi-automated text-mining, manual curation and integration of data from various sources.

Streamlining and simplification of microbial genome architecture

The genomes of unicellular species, particularly prokaryotes, are greatly reduced in size and simplified in terms of gene structure relative to those of multicellular eukaryotes. Arguments proposed to explain this disparity include selection for metabolic efficiency and elevated rates of deletion in microbes, but the evidence in support of these hypotheses is at best equivocal. An alternative explanation based on fundamental population-genetic principles is proposed here. By increasing the mutational target sizes of associated genes, most forms of nonfunctional DNA are opposed by weak selection. Free-living microbial species have elevated effective population sizes, and the consequent reduction in the power of random genetic drift appears to be sufficient to enable natural selection to inhibit the accumulation of excess DNA. This hypothesis provides a potentially unifying explanation for the continuity in genomic scaling from prokaryotes to multicellular eukaryotes, the divergent patterns of mitochondrial evolution in animals and land plants, and various aspects of genomic modification in microbial endosymbionts.

Protein Superfamily Evolution and the Last Universal Common Ancestor (LUCA)

By exploiting three-dimensional structure comparison, which is more sensitive than conventional sequence-based methods for detecting remote homology, we have identified a set of 140 ancestral protein domains using very restrictive criteria to minimize the potential error introduced by horizontal gene transfer. These domains are highly likely to have been present in the Last Universal Common Ancestor (LUCA) based on their universality in almost all of 114 completed prokaryotic (Bacteria and Archaea) and eukaryotic genomes. Functional analysis of these ancestral domains reveals a genetically complex LUCA with practically all the essential functional systems present in extant organisms, supporting the theory that life achieved its modern cellular status much before the main kingdom separation (Doolittle 2000). In addition, we have calculated different estimations of the genetic and functional versatility of all the superfamilies and functional groups in the prokaryote subsample. These estimations reveal that some ancestral superfamilies have been more versatile than others during evolution allowing more genetic and functional variation. Furthermore, the differences in genetic versatility between protein families are more attributable to their functional nature rather than the time that they have been evolving. These differences in tolerance to mutation suggest that some protein families have eroded their phylogenetic signal faster than others, hiding in many cases, their ancestral origin and suggesting that the calculation of 140 ancestral domains is probably an underestimate.

Exploiting protein structure data to explore the evolution of protein function and biological complexity

New directions in biology are being driven by the complete sequencing of genomes, which has given us the protein repertoires of diverse organisms from all kingdoms of life. In tandem with this accumulation of sequence data, worldwide structural genomics initiatives, advanced by the development of improved technologies in X-ray crystallography and NMR, are expanding our knowledge of structural families and increasing our fold libraries. Methods for detecting remote sequence similarities have also been made more sensitive and this means that we can map domains from these structural families onto genome sequences to understand how these families are distributed throughout the genomes and reveal how they might influence the functional repertoires and biological complexities of the organisms. We have used robust protocols to assign sequences from completed genomes to domain structures in the CATH database, allowing up to 60% of domain sequences in these genomes, depending on the organism, to be assigned to a domain family of known structure. Analysis of the distribution of these families throughout bacterial genomes identified more than 300 universal families, some of which had expanded significantly in proportion to genome size. These highly expanded families are primarily involved in metabolism and regulation and appear to make major contributions to the functional repertoire and complexity of bacterial organisms. When comparisons are made across all kingdoms of life, we find a smaller set of universal domain families (approx. 140), of which families involved in protein biosynthesis are the largest conserved component. Analysis of the behaviour of other families reveals that some (e.g. those involved in metabolism, regulation) have remained highly innovative during evolution, making it harder to trace their evolutionary ancestry. Structural analyses of metabolic families provide some insights into the mechanisms of functional innovation, which include changes in domain partnerships and significant structural embellishments leading to modulation of active sites and protein interactions.

The complement of enzymatic sets in different species

We present here a comprehensive analysis of the complement of enzymes in a large variety of species. As enzymes are a relatively conserved group there are several classification systems available that are common to all species and link a protein sequence to an enzymatic function. Enzymes are therefore an ideal functional group to study the relationship between sequence expansion, functional divergence and phenotypic changes. By using information retrieved from the well annotated SWISS-PROT database together with sequence information from a variety of fully sequenced genomes and information from the EC functional scheme we have aimed here to estimate the fraction of enzymes in genomes, to determine the extent of their functional redundancy in different domains of life and to identify functional innovations and lineage specific expansions in the metazoa lineage. We found that prokaryote and eukaryote species differ both in the fraction of enzymes in their genomes and in the pattern of expansion of their enzymatic sets. We observe an increase in functional redundancy accompanying an increase in species complexity. A quantitative assessment was performed in order to determine the degree of functional redundancy in different species. Finally, we report a massive expansion in the number of mammalian enzymes involved in signalling and degradation.

PROTEIN FAMILIES AND THEIR EVOLUTION—A STRUCTURAL PERSPECTIVE

We can now assign about two thirds of the sequences from completed genomes to as few as 1400 domain families for which structures are known and thus more ancient evolutionary relationships established. About 200 of these domain families are common to all kingdoms of life and account for nearly 50% of domain structure annotations in the genomes. Some of these domain families have been very extensively duplicated within a genome and combined with different domain partners giving rise to different multidomain proteins. The ways in which these domain combinations evolve tend to be specific to the organism so that less than 15% of the protein families found within a genome appear to be common to all kingdoms of life. Recent analyses of completed genomes, exploiting the structural data, have revealed the extent to which duplication of these domains and modifications of their functions can expand the functional repertoire of the organism, contributing to increasing complexity.