Inference and analysis of the relative stability of bacterial chromosomes

Comparative analysis of mitochondrial genomes in Ceratocystis fimbriata complex across diverse hosts

The decline ofAcacia mangiumWilld. in Malaysia, especially in Sabah since 2010, is primarily due to Ceratocystiswilt and canker disease (CWCD) caused by theCeratocystis fimbriataEllis & Halst. complex. This study was aimed to investigate the mitochondrial genome architecture of two differentC. fimbriatacomplex isolates from Malaysia: one fromA. mangiumin Pahang (FRIM1162) and another fromEucalyptus pellitain Sarawak (FRIM1441). This research employed Next-Generation Sequencing (NGS) to contrast genomes from diverse hosts with nine additional mitochondrial sequences, identifying significant genetic diversity and mutational hotspots in the mitochondrial genome alignment. The mitochondrial genome-based phylogenetic analysis revealed a significant genetic relationship between the studied isolates and theC. fimbriatacomplex in the South American Subclade, indicating that theC. fimbriatacomplex discovered in Malaysia isC. manginecans. The comparative mitochondrial genome demonstrates the adaptability of the complex due to mobile genetic components and genomic rearrangements in the studiedfungal isolates. This research enhances our knowledge of the genetic diversity and evolutionary patterns within theC. fimbriatacomplex, aiding in a deeper understanding of fungal disease development and host adaption processes. The acquired insights are crucial for creating specific management strategies for CWCD, improving the overall understanding of fungal disease evolution and control.

A fast comparative genome browser for diverse bacteria and archaea

Genome sequencing has revealed an incredible diversity of bacteria and archaea, but there are no fast and convenient tools for browsing across these genomes. It is cumbersome to view the prevalence of homologs for a protein of interest, or the gene neighborhoods of those homologs, across the diversity of the prokaryotes. We developed a web-based tool, fast.genomics, that uses two strategies to support fast browsing across the diversity of prokaryotes. First, the database of genomes is split up. The main database contains one representative from each of the 6,377 genera that have a high-quality genome, and additional databases for each taxonomic order contain up to 10 representatives of each species. Second, homologs of proteins of interest are identified quickly by using accelerated searches, usually in a few seconds. Once homologs are identified, fast.genomics can quickly show their prevalence across taxa, view their neighboring genes, or compare the prevalence of two different proteins. Fast.genomics is available at https://fast.genomics.lbl.gov.

Relaxed purifying selection pressure drives accelerated and dynamic gene rearrangements in thrips (Insecta: Thysanoptera) mitochondrial genomes

Insect mitochondrial genomes (mitogenome) generally present a typical gene order, which is considered as the ancestral arrangement. All sequenced mitogenomes in the Thysanoptera display high levels of gene rearrangement. Due to limited number of thrips mitogenomes sequenced, how gene rearrangement may be shaped by evolution remain unclear. Here, we analyzed 33 thrips mitogenomes, including 14 newly sequenced. These mitogenomes were diverse in organization, nucleotides substitution and gene arrangements. We found 28 highly rearranged gene orders with the breakpoints of gene rearrangements from 25 to 33. Reconstruction of the ancestors mitochondrial gene arrangements states indicated that Tubulifera have more complex pathways than Terebrantia in the gene order evolution. Molecular calibration estimated that divergence of two suborders occurred in the middle Triassic while the radiation of thrips was associated with the arose and flourish of angiosperm. Our evolutionary hypothesis testing suggests that relaxation of selection pressure enabled the early phase of Thysanoptera evolution, followed by a stronger selective pressure fixed diversification. Our analyses found gene inversion increases the nonsynonymous substitution rates and provide an evolutionary hypothesis driving the diverse gene orders.

Positive selection during niche adaptation results in large-scale and irreversible rearrangement of chromosomal gene order in bacteria

Analysis of bacterial genomes shows that, whereas diverse species share many genes in common, their linear order on the chromosome is often not conserved. Whereas rearrangements in gene order could occur by genetic drift, an alternative hypothesis is rearrangement driven by positive selection during niche adaptation (SNAP). Here, we provide the first experimental support for the SNAP hypothesis. We evolved Salmonella to adapt to growth on malate as the sole carbon source and followed the evolutionary trajectories. The initial adaptation to growth in the new environment involved the duplication of 1.66 Mb, corresponding to one-third of the Salmonella chromosome. This duplication is selected to increase the copy number of a single gene, dctA, involved in the uptake of malate. Continuing selection led to the rapid loss or mutation of duplicate genes from either copy of the duplicated region. After 2000 generations, only 31% of the originally duplicated genes remained intact and the gene order within the Salmonella chromosome has been significantly and irreversibly altered. These results experientially validate predictions made by the SNAP hypothesis and show that SNAP can be a strong driving force for rearrangements in chromosomal gene order.

Deep conservation of histone variants in Thermococcales archaea

Histones are ubiquitous in eukaryotes where they assemble into nucleosomes, binding and wrapping DNA to form chromatin. One process to modify chromatin and regulate DNA accessibility is the replacement of histones in the nucleosome with paralogous variants. Histones are also present in archaea but whether and how histone variants contribute to the generation of different physiologically relevant chromatin states in these organisms remains largely unknown. Conservation of paralogs with distinct properties can provide prima facie evidence for defined functional roles. We recently revealed deep conservation of histone paralogs with different properties in the Methanobacteriales, but little is known experimentally about these histones. In contrast, the two histones of the model archaeon Thermococcus kodakarensis, HTkA and HTkB, have been examined in some depth, both in vitro and in vivo. HTkA and HTkB exhibit distinct DNA-binding behaviors and elicit unique transcriptional responses when deleted. Here, we consider the evolution of HTkA/B and their orthologs across the order Thermococcales. We find histones with signature HTkA- and HTkB-like properties to be present in almost all Thermococcales genomes. Phylogenetic analysis indicates the presence of one HTkA- and one HTkB-like histone in the ancestor of Thermococcales and long-term maintenance of these two paralogs throughout Thermococcales diversification. Our results support the notion that archaea and eukaryotes have convergently evolved histone variants that carry out distinct adaptive functions. Intriguingly, we also detect more highly diverged histone-fold proteins, related to those found in some bacteria, in several Thermococcales genomes. The functions of these bacteria-type histones remain unknown, but structural modeling suggests that they can form heterodimers with HTkA/B-like histones.

Reconstructing the Evolutionary History of a Highly Conserved Operon Cluster in Gammaproteobacteria and Bacilli

The evolution of gene order rearrangements within bacterial chromosomes is a fast process. Closely related species can have almost no conservation in long-range gene order. A prominent exception to this rule is a >40 kb long cluster of five core operons (secE-rpoBC-str-S10-spc-alpha) and three variable adjacent operons (cysS, tufB, and ecf) that together contain 57 genes of the transcriptional and translational machinery. Previous studies have indicated that at least part of this operon cluster might have been present in the last common ancestor of bacteria and archaea. Using 204 whole genome sequences, ∼2 Gy of evolution of the operon cluster were reconstructed back to the last common ancestors of the Gammaproteobacteria and of the Bacilli. A total of 163 independent evolutionary events were identified in which the operon cluster was altered. Further examination showed that the process of disconnecting two operons generally follows the same pattern. Initially, a small number of genes is inserted between the operons breaking the concatenation followed by a second event that fully disconnects the operons. While there is a general trend for loss of gene synteny over time, there are examples of increased alteration rates at specific branch points or within specific bacterial orders. This indicates the recurrence of relaxed selection on the gene order within bacterial chromosomes. The analysis of the alternation events indicates that segmental genome duplications and/or transposon-directed recombination play a crucial role in rearrangements of the operon cluster.

Herpesvirus-bacteria synergistic interaction in periodontitis

The etiopathogenesis of severe periodontitis includes herpesvirus-bacteria coinfection. This article evaluates the pathogenicity of herpesviruses (cytomegalovirus and Epstein-Barr virus) and periodontopathic bacteria (Aggregatibacter actinomycetemcomitans and Porphyromonas gingivalis) and coinfection of these infectious agents in the initiation and progression of periodontitis. Cytomegalovirus and A. actinomycetemcomitans/P. gingivalis exercise synergistic pathogenicity in the development of localized ("aggressive") juvenile periodontitis. Cytomegalovirus and Epstein-Barr virus are associated with P. gingivalis in adult types of periodontitis. Periodontal herpesviruses that enter the general circulation may also contribute to disease development in various organ systems. A 2-way interaction is likely to occur between periodontal herpesviruses and periodontopathic bacteria, with herpesviruses promoting bacterial upgrowth, and bacterial factors reactivating latent herpesviruses. Bacterial-induced gingivitis may facilitate herpesvirus colonization of the periodontium, and herpesvirus infections may impede the antibacterial host defense and alter periodontal cells to predispose for bacterial adherence and invasion. Herpesvirus-bacteria synergistic interactions, are likely to comprise an important pathogenic determinant of aggressive periodontitis. However, mechanistic investigations into the molecular and cellular interaction between periodontal herpesviruses and bacteria are still scarce. Herpesvirus-bacteria coinfection studies may yield significant new discoveries of pathogenic determinants, and drug and vaccine targets to minimize or prevent periodontitis and periodontitis-related systemic diseases.

The SNAP hypothesis: Chromosomal rearrangements could emerge from positive Selection during Niche Adaptation

The relative linear order of most genes on bacterial chromosomes is not conserved over evolutionary timescales. One explanation is that selection is weak, allowing recombination to randomize gene order by genetic drift. However, most chromosomal rearrangements are deleterious to fitness. In contrast, we propose the hypothesis that rearrangements in gene order are more likely the result of selection during niche adaptation (SNAP). Partial chromosomal duplications occur very frequently by recombination between direct repeat sequences. Duplicated regions may contain tens to hundreds of genes and segregate quickly unless maintained by selection. Bacteria exposed to non-lethal selections (for example, a requirement to grow on a poor nutrient) can adapt by maintaining a duplication that includes a gene that improves relative fitness. Further improvements in fitness result from the loss or inactivation of non-selected genes within each copy of the duplication. When genes that are essential in single copy are lost from different copies of the duplication, segregation is prevented even if the original selection is lifted. Functional gene loss continues until a new genetic equilibrium is reached. The outcome is a rearranged gene order. Mathematical modelling shows that this process of positive selection to adapt to a new niche can rapidly drive rearrangements in gene order to fixation. Signature features (duplication formation and divergence) of the SNAP model were identified in natural isolates from multiple species showing that the initial two steps in the SNAP process can occur with a remarkably high frequency. Further bioinformatic and experimental analyses are required to test if and to which extend the SNAP process acts on bacterial genomes.

All life on earth has evolved from a universal common ancestor with a specific order of genes on the chromosome. This order is not maintained in modern species and the standard hypothesis is that changes reflect a lack of strong selection on gene order. Here, we propose an alternative hypothesis, SNAP. The occupation of a novel environment by bacteria is generally a trade-off situation. For example, while the bacteria may not be adapted to grow well under the new conditions, they may benefit by not having to share available resources with other microorganisms. Bacterial populations frequently acquire duplications of chromosomal segments containing genes that can help them adapt to a new environment. Other genes that are also duplicated are not required in two copies so that over time a superfluous copy can be lost. Eventually, the process of duplication and gene loss can lead to the rearrangement of the gene order in the chromosomal segment. The major benefit of this model over the standard hypothesis is that the process is driven by positive selection and can reach fixation rapidly.

Subtelomeres are fast-evolving regions of the Streptomyces linear chromosome

Streptomyces possess a large linear chromosome (6-12 Mb) consisting of a conserved central region flanked by variable arms covering several megabases. In order to study the evolution of the chromosome across evolutionary times, a representative panel of Streptomyces strains and species (125) whose chromosomes are completely sequenced and assembled was selected. The pan-genome of the genus was modelled and shown to be open with a core-genome reaching 1018 genes. The evolution of Streptomyces chromosome was analysed by carrying out pairwise comparisons, and by monitoring indexes measuring the conservation of genes (presence/absence) and their synteny along the chromosome. Using the phylogenetic depth offered by the chosen panel, it was possible to infer that within the central region of the chromosome, the core-genes form a highly conserved organization, which can reveal the existence of an ancestral chromosomal skeleton. Conversely, the chromosomal arms, enriched in variable genes evolved faster than the central region under the combined effect of rearrangements and addition of new information from horizontal gene transfer. The genes hosted in these regions may be localized there because of the adaptive advantage that their rapid evolution may confer. We speculate that (i) within a bacterial population, the variability of these genes may contribute to the establishment of social characters by the production of 'public goods' (ii) at the evolutionary scale, this variability contributes to the diversification of the genetic pool of the bacteria.

Parallels between experimental and natural evolution of legume symbionts

The emergence of symbiotic interactions has been studied using population genomics in nature and experimental evolution in the laboratory, but the parallels between these processes remain unknown. Here we compare the emergence of rhizobia after the horizontal transfer of a symbiotic plasmid in natural populations of Cupriavidus taiwanensis, over 10 MY ago, with the experimental evolution of symbiotic Ralstonia solanacearum for a few hundred generations. In spite of major differences in terms of time span, environment, genetic background, and phenotypic achievement, both processes resulted in rapid genetic diversification dominated by purifying selection. We observe no adaptation in the plasmid carrying the genes responsible for the ecological transition. Instead, adaptation was associated with positive selection in a set of genes that led to the co-option of the same quorum-sensing system in both processes. Our results provide evidence for similarities in experimental and natural evolutionary transitions and highlight the potential of comparisons between both processes to understand symbiogenesis.

It is unclear if experimental evolution is a good model for natural processes. Here, Clerissi et al. find parallels between the evolution of symbiosis in rhizobia after horizontal transfer of a plasmid over 10 million years ago and experimentally evolved symbionts.

Inferring the Minimal Genome of Mesoplasma florum by Comparative Genomics and Transposon Mutagenesis

The creation and comparison of minimal genomes will help better define the most fundamental mechanisms supporting life. Mesoplasma florum is a near-minimal, fast-growing, nonpathogenic bacterium potentially amenable to genome reduction efforts. In a comparative genomic study of 13 M. florum strains, including 11 newly sequenced genomes, we have identified the core genome and open pangenome of this species. Our results show that all of the strains have approximately 80% of their gene content in common. Of the remaining 20%, 17% of the genes were found in multiple strains and 3% were unique to any given strain. On the basis of random transposon mutagenesis, we also estimated that ~290 out of 720 genes are essential for M. florum L1 in rich medium. We next evaluated different genome reduction scenarios for M. florum L1 by using gene conservation and essentiality data, as well as comparisons with the first working approximation of a minimal organism, Mycoplasma mycoides JCVI-syn3.0. Our results suggest that 409 of the 473 M. mycoides JCVI-syn3.0 genes have orthologs in M. florum L1. Conversely, 57 putatively essential M. florum L1 genes have no homolog in M. mycoides JCVI-syn3.0. This suggests differences in minimal genome compositions, even for these evolutionarily closely related bacteria. IMPORTANCE The last years have witnessed the development of whole-genome cloning and transplantation methods and the complete synthesis of entire chromosomes. Recently, the first minimal cell, Mycoplasma mycoides JCVI-syn3.0, was created. Despite these milestone achievements, several questions remain to be answered. For example, is the composition of minimal genomes virtually identical in phylogenetically related species? On the basis of comparative genomics and transposon mutagenesis, we investigated this question by using an alternative model, Mesoplasma florum, that is also amenable to genome reduction efforts. Our results suggest that the creation of additional minimal genomes could help reveal different gene compositions and strategies that can support life, even within closely related species.

Insights into origin and evolution of α-proteobacterial gene transfer agents

Several bacterial and archaeal lineages produce nanostructures that morphologically resemble small tailed viruses, but, unlike most viruses, contain apparently random pieces of the host genome. Since these elements can deliver the packaged DNA to other cells, they were dubbed gene transfer agents (GTAs). Because many genes involved in GTA production have viral homologs, it has been hypothesized that the GTA ancestor was a virus. Whether GTAs represent an atypical virus, a defective virus, or a virus co-opted by the prokaryotes for some function, remains to be elucidated. To evaluate these possibilities, we examined the distribution and evolutionary histories of genes that encode a GTA in the α-proteobacterium Rhodobacter capsulatus (RcGTA). We report that although homologs of many individual RcGTA genes are abundant across bacteria and their viruses, RcGTA-like genomes are mainly found in one subclade of α-proteobacteria. When compared with the viral homologs, genes of the RcGTA-like genomes evolve significantly slower, and do not have higher %A+T nucleotides than their host chromosomes. Moreover, they appear to reside in stable regions of the bacterial chromosomes that are generally conserved across taxonomic orders. These findings argue against RcGTA being an atypical or a defective virus. Our phylogenetic analyses suggest that RcGTA ancestor likely originated in the lineage that gave rise to contemporary α-proteobacterial orders Rhizobiales, Rhodobacterales, Caulobacterales, Parvularculales, and Sphingomonadales, and since that time the RcGTA-like element has co-evolved with its host chromosomes. Such evolutionary history is compatible with maintenance of these elements by bacteria due to some selective advantage. As for many other prokaryotic traits, horizontal gene transfer played a substantial role in the evolution of RcGTA-like elements, not only in shaping its genome components within the orders, but also in occasional dissemination of RcGTA-like regions across the orders and even to different bacterial phyla.

Recruitment of a Lineage-Specific Virulence Regulatory Pathway Promotes Intracellular Infection by a Plant Pathogen Experimentally Evolved into a Legume Symbiont

Ecological transitions between different lifestyles, such as pathogenicity, mutualism and saprophytism, have been very frequent in the course of microbial evolution, and often driven by horizontal gene transfer. Yet, how genomes achieve the ecological transition initiated by the transfer of complex biological traits remains poorly known. Here, we used experimental evolution, genomics, transcriptomics and high-resolution phenotyping to analyze the evolution of the plant pathogen Ralstonia solanacearum into legume symbionts, following the transfer of a natural plasmid encoding the essential mutualistic genes. We show that a regulatory pathway of the recipient R. solanacearum genome involved in extracellular infection of natural hosts was reused to improve intracellular symbiosis with the Mimosa pudica legume. Optimization of intracellular infection capacity was gained through mutations affecting two components of a new regulatory pathway, the transcriptional regulator efpR and a region upstream from the RSc0965-0967 genes of unknown functions. Adaptive mutations caused the downregulation of efpR and the over-expression of a downstream regulatory module, the three unknown genes RSc3146-3148, two of which encoding proteins likely associated to the membrane. This over-expression led to important metabolic and transcriptomic changes and a drastic qualitative and quantitative improvement of nodule intracellular infection. In addition, these adaptive mutations decreased the virulence of the original pathogen. The complete efpR/RSc3146-3148 pathway could only be identified in the genomes of the pathogenic R. solanacearum species complex. Our findings illustrate how the rewiring of a genetic network regulating virulence allows a radically different type of symbiotic interaction and contributes to ecological transitions and trade-offs.

Mitochondrial genome evolution in the Saccharomyces sensu stricto complex

Exploring the evolutionary patterns of mitochondrial genomes is important for our understanding of the Saccharomyces sensu stricto (SSS) group, which is a model system for genomic evolution and ecological analysis. In this study, we first obtained the complete mitochondrial sequences of two important species, Saccharomyces mikatae and Saccharomyces kudriavzevii. We then compared the mitochondrial genomes in the SSS group with those of close relatives, and found that the non-coding regions evolved rapidly, including dramatic expansion of intergenic regions, fast evolution of introns and almost 20-fold higher rearrangement rates than those of the nuclear genomes. However, the coding regions, and especially the protein-coding genes, are more conserved than those in the nuclear genomes of the SSS group. The different evolutionary patterns of coding and non-coding regions in the mitochondrial and nuclear genomes may be related to the origin of the aerobic fermentation lifestyle in this group. Our analysis thus provides novel insights into the evolution of mitochondrial genomes.

Abundance and co-occurrence of extracellular capsules increase environmental breadth: Implications for the emergence of pathogens

Extracellular capsules constitute the outermost layer of many bacteria, are major virulence factors, and affect antimicrobial therapies. They have been used as epidemiological markers and recently became vaccination targets. Despite the efforts to biochemically serotype capsules in a few model pathogens, little is known of their taxonomic and environmental distribution. We developed, validated, and made available a computational tool, CapsuleFinder, to identify capsules in genomes. The analysis of over 2500 prokaryotic genomes, accessible in a database, revealed that ca. 50% of them—including Archaea—encode a capsule. The Wzx/Wzy-dependent capsular group was by far the most abundant. Surprisingly, a fifth of the genomes encode more than one capsule system—often from different groups—and their non-random co-occurrence suggests the existence of negative and positive epistatic interactions. To understand the role of multiple capsules, we queried more than 6700 metagenomes for the presence of species encoding capsules and showed that their distribution varied between environmental categories and, within the human microbiome, between body locations. Species encoding capsules, and especially those encoding multiple capsules, had larger environmental breadths than the other species. Accordingly, capsules were more frequent in environmental bacteria than in pathogens and, within the latter, they were more frequent among facultative pathogens. Nevertheless, capsules were frequent in clinical samples, and were usually associated with fast-growing bacteria with high infectious doses. Our results suggest that capsules increase the environmental range of bacteria and make them more resilient to environmental perturbations. Capsules might allow opportunistic pathogens to profit from empty ecological niches or environmental perturbations, such as those resulting from antibiotic therapy, to colonize the host. Capsule-associated virulence might thus be a by-product of environmental adaptation. Understanding the role of capsules in natural environments might enlighten their function in pathogenesis.

Extracellular capsules protect bacterial cells from external aggressions such as antibiotics or desiccation, but can also be targeted by vaccines. Since little was known about their frequency across Prokaryotes, we created and made freely available a computational tool, CapsuleFinder, to identify them from genomic data. Surprisingly, its use showed that many bacterial strains, especially those with the largest genomes, encode several capsules. The frequencies of the different combinations of capsule groups depended strongly on the phyla and the groups themselves, suggesting the existence of epistatic interactions between capsules. Bacteria encoding capsule systems were found in many natural environments, and were frequent in the human microbiome. In contrast to their frequent association with virulence, we found many more capsules in non-pathogens or facultative pathogens than among obligatory pathogens. We suggest that capsules increase the environmental breadth of bacteria thereby facilitating host colonization by opportunistic pathogens.

Contrasting evolutionary genome dynamics between domesticated and wild yeasts

Structural rearrangements have long been recognized as an important source of genetic variation, with implications in phenotypic diversity and disease, yet their detailed evolutionary dynamics remain elusive. Here we use long-read sequencing to generate end-to-end genome assemblies for 12 strains representing major subpopulations of the partially domesticated yeast Saccharomyces cerevisiae and its wild relative Saccharomyces paradoxus. These population-level high-quality genomes with comprehensive annotation enable precise definition of chromosomal boundaries between cores and subtelomeres and a high-resolution view of evolutionary genome dynamics. In chromosomal cores, S. paradoxus shows faster accumulation of balanced rearrangements (inversions, reciprocal translocations and transpositions), whereas S. cerevisiae accumulates unbalanced rearrangements (novel insertions, deletions and duplications) more rapidly. In subtelomeres, both species show extensive interchromosomal reshuffling, with a higher tempo in S. cerevisiae. Such striking contrasts between wild and domesticated yeasts are likely to reflect the influence of human activities on structural genome evolution.

Elevated Rate of Genome Rearrangements in Radiation-Resistant Bacteria

A number of bacterial, archaeal, and eukaryotic species are known for their resistance to ionizing radiation. One of the challenges these species face is a potent environmental source of DNA double-strand breaks, potential drivers of genome structure evolution. Efficient and accurate DNA double-strand break repair systems have been demonstrated in several unrelated radiation-resistant species and are putative adaptations to the DNA damaging environment. Such adaptations are expected to compensate for the genome-destabilizing effect of environmental DNA damage and may be expected to result in a more conserved gene order in radiation-resistant species. However, here we show that rates of genome rearrangements, measured as loss of gene order conservation with time, are higher in radiation-resistant species in multiple, phylogenetically independent groups of bacteria. Comparison of indicators of selection for genome organization between radiation-resistant and phylogenetically matched, nonresistant species argues against tolerance to disruption of genome structure as a strategy for radiation resistance. Interestingly, an important mechanism affecting genome rearrangements in prokaryotes, the symmetrical inversions around the origin of DNA replication, shapes genome structure of both radiation-resistant and nonresistant species. In conclusion, the opposing effects of environmental DNA damage and DNA repair result in elevated rates of genome rearrangements in radiation-resistant bacteria.

Genomic epidemiology and global diversity of the emerging bacterial pathogen Elizabethkingia anophelis

Elizabethkingia anophelis is an emerging pathogen involved in human infections and outbreaks in distinct world regions. We investigated the phylogenetic relationships and pathogenesis-associated genomic features of two neonatal meningitis isolates isolated 5 years apart from one hospital in Central African Republic and compared them with Elizabethkingia from other regions and sources. Average nucleotide identity firmly confirmed that E. anophelis, E. meningoseptica and E. miricola represent demarcated genomic species. A core genome multilocus sequence typing scheme, broadly applicable to Elizabethkingia species, was developed and made publicly available (http://bigsdb.pasteur.fr/elizabethkingia). Phylogenetic analysis revealed distinct E. anophelis sublineages and demonstrated high genetic relatedness between the African isolates, compatible with persistence of the strain in the hospital environment. CRISPR spacer variation between the African isolates was mirrored by the presence of a large mobile genetic element. The pan-genome of E. anophelis comprised 6,880 gene families, underlining genomic heterogeneity of this species. African isolates carried unique resistance genes acquired by horizontal transfer. We demonstrated the presence of extensive variation of the capsular polysaccharide synthesis gene cluster in E. anophelis. Our results demonstrate the dynamic evolution of this emerging pathogen and the power of genomic approaches for Elizabethkingia identification, population biology and epidemiology.

Functional and Evolutionary Characterization of a Gene Transfer Agent's Multilocus "Genome"

Gene transfer agents (GTAs) are phage-like particles that can package and transfer a random piece of the producing cell’s genome, but are unable to transfer all the genes required for their own production. As such, GTAs represent an evolutionary conundrum: are they selfish genetic elements propagating through an unknown mechanism, defective viruses, or viral structures “repurposed” by cells for gene exchange, as their name implies? In Rhodobacter capsulatus, production of the R. capsulatus GTA (RcGTA) particles is associated with a cluster of genes resembling a small prophage. Utilizing transcriptomic, genetic and biochemical approaches, we report that the RcGTA “genome” consists of at least 24 genes distributed across five distinct loci. We demonstrate that, of these additional loci, two are involved in cell recognition and binding and one in the production and maturation of RcGTA particles. The five RcGTA “genome” loci are widespread within Rhodobacterales, but not all loci have the same evolutionary histories. Specifically, two of the loci have been subject to frequent, probably virus-mediated, gene transfer events. We argue that it is unlikely that RcGTA is a selfish genetic element. Instead, our findings are compatible with the scenario that RcGTA is a virus-derived element maintained by the producing organism due to a selective advantage of within-population gene exchange. The modularity of the RcGTA “genome” is presumably a result of selection on the host organism to retain GTA functionality.

Evolution by gene loss

The recent increase in genomic data is revealing an unexpected perspective of gene loss as a pervasive source of genetic variation that can cause adaptive phenotypic diversity. This novel perspective of gene loss is raising new fundamental questions. How relevant has gene loss been in the divergence of phyla? How do genes change from being essential to dispensable and finally to being lost? Is gene loss mostly neutral, or can it be an effective way of adaptation? These questions are addressed, and insights are discussed from genomic studies of gene loss in populations and their relevance in evolutionary biology and biomedicine.

Coevolution of the Organization and Structure of Prokaryotic Genomes

The cytoplasm of prokaryotes contains many molecular machines interacting directly with the chromosome. These vital interactions depend on the chromosome structure, as a molecule, and on the genome organization, as a unit of genetic information. Strong selection for the organization of the genetic elements implicated in these interactions drives replicon ploidy, gene distribution, operon conservation, and the formation of replication-associated traits. The genomes of prokaryotes are also very plastic with high rates of horizontal gene transfer and gene loss. The evolutionary conflicts between plasticity and organization lead to the formation of regions with high genetic diversity whose impact on chromosome structure is poorly understood. Prokaryotic genomes are remarkable documents of natural history because they carry the imprint of all of these selective and mutational forces. Their study allows a better understanding of molecular mechanisms, their impact on microbial evolution, and how they can be tinkered in synthetic biology.

High variability of mitochondrial gene order among fungi

From their origin as an early alpha proteobacterial endosymbiont to their current state as cellular organelles, large-scale genomic reorganization has taken place in the mitochondria of all main eukaryotic lineages. So far, most studies have focused on plant and animal mitochondrial (mt) genomes (mtDNA), but fungi provide new opportunities to study highly differentiated mtDNAs. Here, we analyzed 38 complete fungal mt genomes to investigate the evolution of mtDNA gene order among fungi. In particular, we looked for evidence of nonhomologous intrachromosomal recombination and investigated the dynamics of gene rearrangements. We investigated the effect that introns, intronic open reading frames (ORFs), and repeats may have on gene order. Additionally, we asked whether the distribution of transfer RNAs (tRNAs) evolves independently to that of mt protein-coding genes. We found that fungal mt genomes display remarkable variation between and within the major fungal phyla in terms of gene order, genome size, composition of intergenic regions, and presence of repeats, introns, and associated ORFs. Our results support previous evidence for the presence of mt recombination in all fungal phyla, a process conspicuously lacking in most Metazoa. Overall, the patterns of rearrangements may be explained by the combined influences of recombination (i.e., most likely nonhomologous and intrachromosomal), accumulated repeats, especially at intergenic regions, and to a lesser extent, mobile element dynamics.

Examination of prokaryotic multipartite genome evolution through experimental genome reduction

Many bacteria carry two or more chromosome-like replicons. This occurs in pathogens such as Vibrio cholerea and Brucella abortis as well as in many N₂-fixing plant symbionts including all isolates of the alfalfa root-nodule bacteria Sinorhizobium meliloti. Understanding the evolution and role of this multipartite genome organization will provide significant insight into these important organisms; yet this knowledge remains incomplete, in part, because technical challenges of large-scale genome manipulations have limited experimental analyses. The distinct evolutionary histories and characteristics of the three replicons that constitute the S. meliloti genome (the chromosome (3.65 Mb), pSymA megaplasmid (1.35 Mb), and pSymB chromid (1.68 Mb)) makes this a good model to examine this topic. We transferred essential genes from pSymB into the chromosome, and constructed strains that lack pSymB as well as both pSymA and pSymB. This is the largest reduction (45.4%, 3.04 megabases, 2866 genes) of a prokaryotic genome to date and the first removal of an essential chromid. Strikingly, strains lacking pSymA and pSymB (ΔpSymAB) lost the ability to utilize 55 of 74 carbon sources and various sources of nitrogen, phosphorous and sulfur, yet the ΔpSymAB strain grew well in minimal salts media and in sterile soil. This suggests that the core chromosome is sufficient for growth in a bulk soil environment and that the pSymA and pSymB replicons carry genes with more specialized functions such as growth in the rhizosphere and interaction with the plant. These experimental data support a generalized evolutionary model, in which non-chromosomal replicons primarily carry genes with more specialized functions. These large secondary replicons increase the organism's niche range, which offsets their metabolic burden on the cell (e.g. pSymA). Subsequent co-evolution with the chromosome then leads to the formation of a chromid through the acquisition of functions core to all niches (e.g. pSymB).

Rhizobia are free-living bacteria of agricultural and environmental importance that form root-nodules on leguminous plants and provide these plants with fixed nitrogen. Many of the rhizobia have a multipartite genome, as do several plant and animal pathogens. All isolates of the alfalfa symbiont, Sinorhizobium meliloti, carry three large replicons, the chromosome (∼3.7 Mb), pSymA megaplasmid (∼1.4 Mb), and pSymB chromid (∼1.7 Mb). To gain insight into the role and evolutionary history of these replicons, we have ‘reversed evolution’ by constructing a S. meliloti strain consisting solely of the chromosome and lacking the pSymB chromid and pSymA megaplasmid. As the resulting strain was viable, we could perform a detailed phenotypic analysis and these data provided significant insight into the biology and metabolism of S. meliloti. The data lend direct experimental evidence in understanding the evolution and role of the multipartite genome. Specifically the large secondary replicons increase the organism's niche range, and this advantage offsets the metabolic burden of these replicons on the cell. Additionally, the single-chromosome strain offers a useful platform to facilitate future forward genetic approaches to understanding and manipulating the symbiosis and plant-microbe interactions.

Natural selection for operons depends on genome size

In prokaryotes, genome size is associated with metabolic versatility, regulatory complexity, effective population size, and horizontal transfer rates. We therefore analyzed the covariation of genome size and operon conservation to assess the evolutionary models of operon formation and maintenance. In agreement with previous results, intraoperonic pairs of essential and of highly expressed genes are more conserved. Interestingly, intraoperonic pairs of genes are also more conserved when they encode proteins at similar cell concentrations, suggesting a role of cotranscription in diminishing the cost of waste and shortfall in gene expression. Larger genomes have fewer and smaller operons that are also less conserved. Importantly, lower conservation in larger genomes was observed for all classes of operons in terms of gene expression, essentiality, and balanced protein concentration. We reached very similar conclusions in independent analyses of three major bacterial clades (α- and β-Proteobacteria and Firmicutes). Operon conservation is inversely correlated to the abundance of transcription factors in the genome when controlled for genome size. This suggests a negative association between the complexity of genetic networks and operon conservation. These results show that genome size and/or its proxies are key determinants of the intensity of natural selection for operon organization. Our data fit better the evolutionary models based on the advantage of coregulation than those based on genetic linkage or stochastic gene expression. We suggest that larger genomes with highly complex genetic networks and many transcription factors endure weaker selection for operons than smaller genomes with fewer alternative tools for genetic regulation.

The genomic diversification of the whole Acinetobacter genus: origins, mechanisms, and consequences

Bacterial genomics has greatly expanded our understanding of microdiversification patterns within a species, but analyses at higher taxonomical levels are necessary to understand and predict the independent rise of pathogens in a genus. We have sampled, sequenced, and assessed the diversity of genomes of validly named and tentative species of the Acinetobacter genus, a clade including major nosocomial pathogens and biotechnologically important species. We inferred a robust global phylogeny and delimited several new putative species. The genus is very ancient and extremely diverse: Genomes of highly divergent species share more orthologs than certain strains within a species. We systematically characterized elements and mechanisms driving genome diversification, such as conjugative elements, insertion sequences, and natural transformation. We found many error-prone polymerases that may play a role in resistance to toxins, antibiotics, and in the generation of genetic variation. Surprisingly, temperate phages, poorly studied in Acinetobacter, were found to account for a significant fraction of most genomes. Accordingly, many genomes encode clustered regularly interspaced short palindromic repeats (CRISPR)-Cas systems with some of the largest CRISPR-arrays found so far in bacteria. Integrons are strongly overrepresented in Acinetobacter baumannii, which correlates with its frequent resistance to antibiotics. Our data suggest that A. baumannii arose from an ancient population bottleneck followed by population expansion under strong purifying selection. The outstanding diversification of the species occurred largely by horizontal transfer, including some allelic recombination, at specific hotspots preferentially located close to the replication terminus. Our work sets a quantitative basis to understand the diversification of Acinetobacter into emerging resistant and versatile pathogens.

Conserved patterns in bacterial genomes: a conundrum physically tailored by evolutionary tinkering

The proper functioning of bacteria is encoded in their genome at multiple levels or scales, each of which is constrained by specific physical forces. At the smallest spatial scales, interatomic forces dictate the folding and function of proteins and nucleic acids. On longer length scales, stochastic forces emerging from the thermal jiggling of proteins and RNAs impose strong constraints on the organization of genes along chromosomes, more particularly in the context of the building of nucleoprotein complexes and the operational mode of regulatory agents. At the cellular level, transcription, replication and cell division activities generate forces that act on both the internal structure and cellular location of chromosomes. The overall result is a complex multi-scale organization of genomes that reflects the evolutionary tinkering of bacteria. The goal of this review is to highlight avenues for deciphering this complexity by focusing on patterns that are conserved among evolutionarily distant bacteria. To this end, I discuss three different organizational scales: the protein structures, the chromosomal organization of genes and the global structure of chromosomes.

Elucidation of operon structures across closely related bacterial genomes

About half of the protein-coding genes in prokaryotic genomes are organized into operons to facilitate co-regulation during transcription. With the evolution of genomes, operon structures are undergoing changes which could coordinate diverse gene expression patterns in response to various stimuli during the life cycle of a bacterial cell. Here we developed a graph-based model to elucidate the diversity of operon structures across a set of closely related bacterial genomes. In the constructed graph, each node represents one orthologous gene group (OGG) and a pair of nodes will be connected if any two genes, from the corresponding two OGGs respectively, are located in the same operon as immediate neighbors in any of the considered genomes. Through identifying the connected components in the above graph, we found that genes in a connected component are likely to be functionally related and these identified components tend to form treelike topology, such as paths and stars, corresponding to different biological mechanisms in transcriptional regulation as follows. Specifically, (i) a path-structure component integrates genes encoding a protein complex, such as ribosome; and (ii) a star-structure component not only groups related genes together, but also reflects the key functional roles of the central node of this component, such as the ABC transporter with a transporter permease and substrate-binding proteins surrounding it. Most interestingly, the genes from organisms with highly diverse living environments, i.e., biomass degraders and animal pathogens of clostridia in our study, can be clearly classified into different topological groups on some connected components.

Temperature-driven differential gene expression by RNA thermosensors

Many prokaryotic genes are organized in operons. Genes organized in such transcription units are co-transcribed into a polycistronic mRNA. Despite being clustered in a single mRNA, individual genes can be subjected to differential regulation, which is mainly achieved at the level of translation depending on initiation and elongation. Efficiency of translation initiation is primarily determined by the structural accessibility of the ribosome binding site (RBS). Structured cis-regulatory elements like RNA thermometers (RNATs) can contribute to differential regulation of individual genes within a polycistronic mRNA. RNATs are riboregulators that mediate temperature-responsive regulation of a downstream gene by modulating the accessibility of its RBS. At low temperature, the RBS is trapped by intra-molecular base pairing prohibiting translation initiation. The secondary structure melts with increasing temperature thus liberating the RBS. Here, we present an overview of different RNAT types and specifically highlight recently discovered RNATs. The main focus of this review is on RNAT-based differential control of polycistronic operons. Finally, we discuss the influence of temperature on other riboregulators and the potential of RNATs in synthetic RNA biology. This article is part of a Special Issue entitled: Riboswitches.

Safety in numbers: multiple occurrences of highly similar homologs among Azotobacter vinelandii carbohydrate metabolism proteins probably confer adaptive benefits

BACKGROUND

Gene duplication and horizontal gene transfer are common processes in bacterial and archaeal genomes, and are generally assumed to result in either diversification or loss of the redundant gene copies. However, a recent analysis of the genome of the soil bacterium Azotobacter vinelandii DJ revealed an abundance of highly similar homologs among carbohydrate metabolism genes. In many cases these multiple genes did not appear to be the result of recent duplications, or to function only as a means of stimulating expression by increasing gene dosage, as the homologs were located in varying functional genetic contexts. Based on these initial findings we here report in-depth bioinformatic analyses focusing specifically on highly similar intra-genome homologs, or synologs, among carbohydrate metabolism genes, as well as an analysis of the general occurrence of very similar synologs in prokaryotes.

RESULTS

Approximately 900 bacterial and archaeal genomes were analysed for the occurrence of synologs, both in general and among carbohydrate metabolism genes specifically. This showed that large numbers of highly similar synologs among carbohydrate metabolism genes are very rare in bacterial and archaeal genomes, and that the A. vinelandii DJ genome contains an unusually large amount of such synologs. The majority of these synologs were found to be non-tandemly organized and localized in varying but metabolically relevant genomic contexts. The same observation was made for other genomes harbouring high levels of such synologs. It was also shown that highly similar synologs generally constitute a very small fraction of the protein-coding genes in prokaryotic genomes. The overall synolog fraction of the A. vinelandii DJ genome was well above the data set average, but not nearly as remarkable as the levels observed when only carbohydrate metabolism synologs were considered.

CONCLUSIONS

Large numbers of highly similar synologs are rare in bacterial and archaeal genomes, both in general and among carbohydrate metabolism genes. However, A. vinelandii and several other soil bacteria harbour large numbers of highly similar carbohydrate metabolism synologs which seem not to result from recent duplication or transfer events. These genes may confer adaptive benefits with respect to certain lifestyles and environmental factors, most likely due to increased regulatory flexibility and/or increased gene dosage.

Variation of the virus-related elements within syntenic genomes of the hyperthermophilic Archaeon Aeropyrum

The increasing number of genome sequences of archaea and bacteria show their adaptation to different environmental conditions at the genomic level. Aeropyrum spp. are aerobic and hyperthermophilic archaea. Aeropyrum camini was isolated from a deep-sea hydrothermal vent, and Aeropyrum pernix was isolated from a coastal solfataric vent. To investigate the adaptation strategy in each habitat, we compared the genomes of the two species. Shared genome features were a small genome size, a high GC content, and a large portion of orthologous genes (86 to 88%). The genomes also showed high synteny. These shared features may have been derived from the small number of mobile genetic elements and the lack of a RecBCD system, a recombinational enzyme complex. In addition, the specialized physiology (aerobic and hyperthermophilic) of Aeropyrum spp. may also contribute to the entire-genome similarity. Despite having stable genomes, interference of synteny occurred with two proviruses, A. pernix spindle-shaped virus 1 (APSV1) and A. pernix ovoid virus 1 (APOV1), and clustered regularly interspaced short palindromic repeat (CRISPR) elements. Spacer sequences derived from the A. camini CRISPR showed significant matches with protospacers of the two proviruses infecting A. pernix, indicating that A. camini interacted with viruses closely related to APSV1 and APOV1. Furthermore, a significant fraction of the nonorthologous genes (41 to 45%) were proviral genes or ORFans probably originating from viruses. Although the genomes of A. camini and A. pernix were conserved, we observed nonsynteny that was attributed primarily to virus-related elements. Our findings indicated that the genomic diversification of Aeropyrum spp. is substantially caused by viruses.

Short and long-term genome stability analysis of prokaryotic genomes

Background

Gene organization dynamics is actively studied because it provides useful evolutionary information, makes functional annotation easier and often enables to characterize pathogens. There is therefore a strong interest in understanding the variability of this trait and the possible correlations with life-style. Two kinds of events affect genome organization: on one hand translocations and recombinations change the relative position of genes shared by two genomes (i.e. the backbone gene order); on the other, insertions and deletions leave the backbone gene order unchanged but they alter the gene neighborhoods by breaking the syntenic regions. A complete picture about genome organization evolution therefore requires to account for both kinds of events.

Results

We developed an approach where we model chromosomes as graphs on which we compute different stability estimators; we consider genome rearrangements as well as the effect of gene insertions and deletions. In a first part of the paper, we fit a measure of backbone gene order conservation (hereinafter called backbone stability) against phylogenetic distance for over 3000 genome comparisons, improving existing models for the divergence in time of backbone stability. Intra- and inter-specific comparisons were treated separately to focus on different time-scales. The use of multiple genomes of a same species allowed to identify genomes with diverging gene order with respect to their conspecific. The inter-species analysis indicates that pathogens are more often unstable with respect to non-pathogens. In a second part of the text, we show that in pathogens, gene content dynamics (insertions and deletions) have a much more dramatic effect on genome organization stability than backbone rearrangements.

Conclusion

In this work, we studied genome organization divergence taking into account the contribution of both genome order rearrangements and genome content dynamics. By studying species with multiple sequenced genomes available, we were able to explore genome organization stability at different time-scales and to find significant differences for pathogen and non-pathogen species. The output of our framework also allows to identify the conserved gene clusters and/or partial occurrences thereof, making possible to explore how gene clusters assembled during evolution.

Streamlining and core genome conservation among highly divergent members of the SAR11 clade

SAR11 is an ancient and diverse clade of heterotrophic bacteria that are abundant throughout the world’s oceans, where they play a major role in the ocean carbon cycle. Correlations between the phylogenetic branching order and spatiotemporal patterns in cell distributions from planktonic ocean environments indicate that SAR11 has evolved into perhaps a dozen or more specialized ecotypes that span evolutionary distances equivalent to a bacterial order. We isolated and sequenced genomes from diverse SAR11 cultures that represent three major lineages and encompass the full breadth of the clade. The new data expand observations about genome evolution and gene content that previously had been restricted to the SAR11 Ia subclade, providing a much broader perspective on the clade’s origins, evolution, and ecology. We found small genomes throughout the clade and a very high proportion of core genome genes (48 to 56%), indicating that small genome size is probably an ancestral characteristic. In their level of core genome conservation, the members of SAR11 are outliers, the most conserved free-living bacteria known. Shared features of the clade include low GC content, high gene synteny, a large hypervariable region bounded by rRNA genes, and low numbers of paralogs. Variation among the genomes included genes for phosphorus metabolism, glycolysis, and C1 metabolism, suggesting that adaptive specialization in nutrient resource utilization is important to niche partitioning and ecotype divergence within the clade. These data provide support for the conclusion that streamlining selection for efficient cell replication in the planktonic habitat has occurred throughout the evolution and diversification of this clade.

The SAR11 clade is the most abundant group of marine microorganisms worldwide, making them key players in the global carbon cycle. Growing knowledge about their biochemistry and metabolism is leading to a more mechanistic understanding of organic carbon oxidation and sequestration in the oceans. The discovery of small genomes in SAR11 provided crucial support for the theory that streamlining selection can drive genome reduction in low-nutrient environments. Study of isolates in culture revealed atypical organic nutrient requirements that can be attributed to genome reduction, such as conditional auxotrophy for glycine and its precursors, a requirement for reduced sulfur compounds, and evidence for widespread cycling of C1 compounds in marine environments. However, understanding the genetic variation and distribution of such pathways and characteristics like streamlining throughout the group has required the isolation and genome sequencing of diverse SAR11 representatives, an analysis of which we provide here.

Evolution of conjugation and type IV secretion systems

Genetic exchange by conjugation is responsible for the spread of resistance, virulence, and social traits among prokaryotes. Recent works unraveled the functioning of the underlying type IV secretion systems (T4SS) and its distribution and recruitment for other biological processes (exaptation), notably pathogenesis. We analyzed the phylogeny of key conjugation proteins to infer the evolutionary history of conjugation and T4SS. We show that single-stranded DNA (ssDNA) and double-stranded DNA (dsDNA) conjugation, while both based on a key AAA⁺ ATPase, diverged before the last common ancestor of bacteria. The two key ATPases of ssDNA conjugation are monophyletic, having diverged at an early stage from dsDNA translocases. Our data suggest that ssDNA conjugation arose first in diderm bacteria, possibly Proteobacteria, and then spread to other bacterial phyla, including bacterial monoderms and Archaea. Identifiable T4SS fall within the eight monophyletic groups, determined by both taxonomy and structure of the cell envelope. Transfer to monoderms might have occurred only once, but followed diverse adaptive paths. Remarkably, some Firmicutes developed a new conjugation system based on an atypical relaxase and an ATPase derived from a dsDNA translocase. The observed evolutionary rates and patterns of presence/absence of specific T4SS proteins show that conjugation systems are often and independently exapted for other functions. This work brings a natural basis for the classification of all kinds of conjugative systems, thus tackling a problem that is growing as fast as genomic databases. Our analysis provides the first global picture of the evolution of conjugation and shows how a self-transferrable complex multiprotein system has adapted to different taxa and often been recruited by the host. As conjugation systems became specific to certain clades and cell envelopes, they may have biased the rate and direction of gene transfer by conjugation within prokaryotes.

Phylogenetic incongruence in E. coli O104: understanding the evolutionary relationships of emerging pathogens in the face of homologous recombination

Escherichia coli O104:H4 was identified as an emerging pathogen during the spring and summer of 2011 and was responsible for a widespread outbreak that resulted in the deaths of 50 people and sickened over 4075. Traditional phenotypic and genotypic assays, such as serotyping, pulsed field gel electrophoresis (PFGE), and multilocus sequence typing (MLST), permit identification and classification of bacterial pathogens, but cannot accurately resolve relationships among genotypically similar but pathotypically different isolates. To understand the evolutionary origins of E. coli O104:H4, we sequenced two strains isolated in Ontario, Canada. One was epidemiologically linked to the 2011 outbreak, and the second, unrelated isolate, was obtained in 2010. MLST analysis indicated that both isolates are of the same sequence type (ST678), but whole-genome sequencing revealed differences in chromosomal and plasmid content. Through comprehensive phylogenetic analysis of five O104:H4 ST678 genomes, we identified 167 genes in three gene clusters that have undergone homologous recombination with distantly related E. coli strains. These recombination events have resulted in unexpectedly high sequence diversity within the same sequence type. Failure to recognize or adjust for homologous recombination can result in phylogenetic incongruence. Understanding the extent of homologous recombination among different strains of the same sequence type may explain the pathotypic differences between the ON2010 and ON2011 strains and help shed new light on the emergence of this new pathogen.

A semi-quantitative, synteny-based method to improve functional predictions for hypothetical and poorly annotated bacterial and archaeal genes

During microbial evolution, genome rearrangement increases with increasing sequence divergence. If the relationship between synteny and sequence divergence can be modeled, gene clusters in genomes of distantly related organisms exhibiting anomalous synteny can be identified and used to infer functional conservation. We applied the phylogenetic pairwise comparison method to establish and model a strong correlation between synteny and sequence divergence in all 634 available Archaeal and Bacterial genomes from the NCBI database and four newly assembled genomes of uncultivated Archaea from an acid mine drainage (AMD) community. In parallel, we established and modeled the trend between synteny and functional relatedness in the 118 genomes available in the STRING database. By combining these models, we developed a gene functional annotation method that weights evolutionary distance to estimate the probability of functional associations of syntenous proteins between genome pairs. The method was applied to the hypothetical proteins and poorly annotated genes in newly assembled acid mine drainage Archaeal genomes to add or improve gene annotations. This is the first method to assign possible functions to poorly annotated genes through quantification of the probability of gene functional relationships based on synteny at a significant evolutionary distance, and has the potential for broad application.

Community ecology of hot spring cyanobacterial mats: predominant populations and their functional potential

Phototrophic microbial mat communities from 60°C and 65°C regions in the effluent channels of Mushroom and Octopus Springs (Yellowstone National Park, WY, USA) were investigated by shotgun metagenomic sequencing. Analyses of assembled metagenomic sequences resolved six dominant chlorophototrophic populations and permitted the discovery and characterization of undescribed but predominant community members and their physiological potential. Linkage of phylogenetic marker genes and functional genes showed novel chlorophototrophic bacteria belonging to uncharacterized lineages within the order Chlorobiales and within the Kingdom Chloroflexi. The latter is the first chlorophototrophic member of Kingdom Chloroflexi that lies outside the monophyletic group of chlorophototrophs of the Order Chloroflexales. Direct comparison of unassembled metagenomic sequences to genomes of representative isolates showed extensive genetic diversity, genomic rearrangements and novel physiological potential in native populations as compared with genomic references. Synechococcus spp. metagenomic sequences showed a high degree of synteny with the reference genomes of Synechococcus spp. strains A and B', but synteny declined with decreasing sequence relatedness to these references. There was evidence of horizontal gene transfer among native populations, but the frequency of these events was inversely proportional to phylogenetic relatedness.

Horizontal transfer, not duplication, drives the expansion of protein families in prokaryotes

Gene duplication followed by neo- or sub-functionalization deeply impacts the evolution of protein families and is regarded as the main source of adaptive functional novelty in eukaryotes. While there is ample evidence of adaptive gene duplication in prokaryotes, it is not clear whether duplication outweighs the contribution of horizontal gene transfer in the expansion of protein families. We analyzed closely related prokaryote strains or species with small genomes (Helicobacter, Neisseria, Streptococcus, Sulfolobus), average-sized genomes (Bacillus, Enterobacteriaceae), and large genomes (Pseudomonas, Bradyrhizobiaceae) to untangle the effects of duplication and horizontal transfer. After removing the effects of transposable elements and phages, we show that the vast majority of expansions of protein families are due to transfer, even among large genomes. Transferred genes--xenologs--persist longer in prokaryotic lineages possibly due to a higher/longer adaptive role. On the other hand, duplicated genes--paralogs--are expressed more, and, when persistent, they evolve slower. This suggests that gene transfer and gene duplication have very different roles in shaping the evolution of biological systems: transfer allows the acquisition of new functions and duplication leads to higher gene dosage. Accordingly, we show that paralogs share most protein-protein interactions and genetic regulators, whereas xenologs share very few of them. Prokaryotes invented most of life's biochemical diversity. Therefore, the study of the evolution of biology systems should explicitly account for the predominant role of horizontal gene transfer in the diversification of protein families.

Markedly different genome arrangements between serotype a strains and serotypes b or c strains of Aggregatibacter actinomycetemcomitans

BACKGROUND

Bacterial phenotype may be profoundly affected by the physical arrangement of their genes in the genome. The Gram-negative species Aggregatibacter actinomycetemcomitans is a major etiologic agent of human periodontitis. Individual clonal types of A. actinomycetemcomitans may exhibit variable virulence and different patterns of disease association. This study examined the genome arrangement of A. actinomycetemcomitans using the genome sequences of serotypes a-c strains. The genome alignment and rearrangement were analyzed by the MAUVE and the GRIMM algorithms. The distribution patterns of genes along the leading/lagging strands were investigated. The occurrence and the location of repeat sequences relative to the genome rearrangement breakpoints were also determined.

RESULTS

The genome arrangement of the serotype a strain D7S-1 is markedly different from the serotype b strain HK1651 or the serotype c strain D11S-1. Specific genome arrangements appear to be conserved among strains of the same serotypes. The reversal distance between D7S-1 and HK1651 by GRIMM analysis is also higher than the within-species comparisons of 7 randomly selected bacterial species. The locations of the orthologous genes are largely preserved between HK1651 and D11S-1 but not between D7S-1 and HK1651 (or D11S-1), irrespective of whether the genes are categorized as essential/nonessential or highly/nonhighly expressed. However, genome rearrangement did not disrupt the operons of the A. actinomycetemcomitans strains. A higher proportion of the genome in strain D7S-1 is occupied by repeat sequences than in strains HK1651 or D11S-1.

CONCLUSION

The results suggest a significant evolutionary divergence between serotype a strains and serotypes b/c strains of A. actinomycetemcomitans. The distinct patterns of genome arrangement may suggest phenotypic differences between serotype a and serotypes b/c strains.

Conditions for the evolution of gene clusters in bacterial genomes

Genes encoding proteins in a common pathway are often found near each other along bacterial chromosomes. Several explanations have been proposed to account for the evolution of these structures. For instance, natural selection may directly favour gene clusters through a variety of mechanisms, such as increased efficiency of coregulation. An alternative and controversial hypothesis is the selfish operon model, which asserts that clustered arrangements of genes are more easily transferred to other species, thus improving the prospects for survival of the cluster. According to another hypothesis (the persistence model), genes that are in close proximity are less likely to be disrupted by deletions. Here we develop computational models to study the conditions under which gene clusters can evolve and persist. First, we examine the selfish operon model by re-implementing the simulation and running it under a wide range of conditions. Second, we introduce and study a Moran process in which there is natural selection for gene clustering and rearrangement occurs by genome inversion events. Finally, we develop and study a model that includes selection and inversion, which tracks the occurrence and fixation of rearrangements. Surprisingly, gene clusters fail to evolve under a wide range of conditions. Factors that promote the evolution of gene clusters include a low number of genes in the pathway, a high population size, and in the case of the selfish operon model, a high horizontal transfer rate. The computational analysis here has shown that the evolution of gene clusters can occur under both direct and indirect selection as long as certain conditions hold. Under these conditions the selfish operon model is still viable as an explanation for the evolution of gene clusters.

Genes involved in a common pathway or function are frequently found near each other on bacterial chromosomes. A number of hypotheses have been previously presented to explain this observation. A particularly influential theory is the selfish operon model, which posits that horizontal transfer could promote gene clustering by favouring transfer of arrangements of genes that are close together. Subsequent theoretical development and analysis of genomic data have contributed to the debate about the plausibility of this model. Here, by re-examining the evolutionary dynamics of gene clusters, we provide and discuss conditions under which gene clusters can evolve. We find that first, some form of bias for clustering is required for clusters to evolve. This bias can be in the form of bias in horizontal transfer towards genes that are close together, or direct natural selection for gene proximity. Our computational work does not present a theoretical obstacle to the selfish operon model as a possible explanation for the evolution of gene clusters.

The birth of new genes by RNA- and DNA-mediated duplication during mammalian evolution

Gene duplication has long been recognized as a major force in genome evolution and has recently been recognized as an important source of individual variation. For many years, the origin of functional gene duplicates was assumed to be whole or partial genome duplication events, but recently retrotransposition has also been shown to contribute new functional protein coding genes and siRNA's. In this study, we utilize pseudogenes to recreate more complete gene family histories, and compare the rates of RNA and DNA-mediated duplication and new functional gene formation in five mammalian genomes. We find that RNA-mediated duplication occurs at a much higher and more variable rate than DNA-mediated duplication, and gives rise to many more duplicated sequences over time. We show that, while the chance of RNA-mediated duplicates becoming functional is much lower than that of their DNA-mediated counterparts, the higher rate of retrotransposition leads to nearly equal contributions of new genes by each mechanism. We also find that functional RNA-mediated duplicates are closer to neighboring genes than non-functional RNA-mediated copies, consistent with co-option of regulatory elements at the site of insertion. Overall, new genes derived from DNA and RNA-mediated duplication mechanisms are under similar levels of purifying selective pressure, but have broadly different functions. RNA-mediated duplication gives rise to a diversity of genes but is dominated by the highly expressed genes of RNA metabolic pathways. DNA-mediated duplication can copy regulatory material along with the protein coding region of the gene and often gives rise to classes of genes whose function are dependent on complex regulatory information. This mechanistic difference may in part explain why we find that mammalian protein families tend to evolve by either one mechanism or the other, but rarely by both. Supplementary Material has been provided (see online Supplementary Material at www.liebertonline.com ).

Genome flexibility in Neisseria meningitidis

Neisseria meningitidis usually lives as a commensal bacterium in the upper airways of humans. However, occasionally some strains can also cause life-threatening diseases such as sepsis and bacterial meningitis. Comparative genomics demonstrates that only very subtle genetic differences between carriage and disease strains might be responsible for the observed virulence differences and that N. meningitidis is, evolutionarily, a very recent species. Comparative genome sequencing also revealed a panoply of genetic mechanisms underlying its enormous genomic flexibility which also might affect the virulence of particular strains. From these studies, N. meningitidis emerges as a paradigm for organisms that use genome variability as an adaptation to changing and thus challenging environments.

Stochasticity in protein levels drives colinearity of gene order in metabolic operons of Escherichia coli

Gene order in some bacterial metabolic operons reflects ordering in the metabolic pathway. That this is true uniquely for operons expressed at low levels highlights the selective importance of fluctuations in protein levels.

In bacterial genomes, gene order is not random. This is most evident when looking at operons, these often encoding enzymes involved in the same metabolic pathway or proteins from the same complex. Is gene order within operons nonrandom, however, and if so why? We examine this issue using metabolic operons as a case study. Using the metabolic network of Escherichia coli, we define the temporal order of reactions. We find a pronounced trend for genes to appear in operons in the same order as they are needed in metabolism (colinearity). This is paradoxical as, at steady state, enzymes abundance should be independent of order within the operon. We consider three extensions of the steady-state model that could potentially account for colinearity: (1) increased productivity associated with higher expression levels of the most 5′ genes, (2) a faster metabolic processing immediately after up-regulation, and (3) metabolic stalling owing to stochastic protein loss. We establish the validity of these hypotheses by employing deterministic and stochastic models of enzyme kinetics. The stochastic stalling hypothesis correctly and uniquely predicts that colinearity is more pronounced both for lowly expressed operons and for genes that are not physically adjacent. The alternative models fail to find any support. These results support the view that stochasticity is a pervasive problem to a cell and that gene order evolution can be driven by the selective consequences of fluctuations in protein levels.

In bacteria, different enzymes from the same metabolic pathway are often encoded within one transcriptional unit, an operon. There is also, we show, a tendency for the enzymes that are needed earlier in the pathway to feature earlier in the operon, so-called colinearity. Why might this be? We test three ideas, one old and two new. The prior suggestion supposes that proteins of genes early in operons will be at a higher dose. Although some operons are like this, in general, we see no relationship of protein dose with colinearity. We also find no evidence that operons that frequently need up-regulation are any more likely to be colinear. A third model is, however, supported. If an operon is rarely expressed, then all the proteins for this part of metabolism can be lost by chance. Rebooting such metabolism is fastest if the operon is colinear. This model predicts, correctly, that colinearity should be more frequent in operons that are expressed at a low level. This result is important for at least two reasons. First, it supports the view that chance events (such as protein loss) within cells are important on a day-to-day basis. Second, it challenges the supposition that natural selection will be weakest on lowly expressed genes. Where chance events are concerned, natural selection can be strong on genes expressed at a low level.

Genesis, effects and fates of repeats in prokaryotic genomes

DNA repeats are causes and consequences of genome plasticity. Repeats are created by intrachromosomal recombination or horizontal transfer. They are targeted by recombination processes leading to amplifications, deletions and rearrangements of genetic material. The identification and analysis of repeats in nearly 700 genomes of bacteria and archaea is facilitated by the existence of sequence data and adequate bioinformatic tools. These have revealed the immense diversity of repeats in genomes, from those created by selfish elements to the ones used for protection against selfish elements, from those arising from transient gene amplifications to the ones leading to stable duplications. Experimental works have shown that some repeats do not carry any adaptive value, while others allow functional diversification and increased expression. All repeats carry some potential to disorganize and destabilize genomes. Because recombination and selection for repeats vary between genomes, the number and types of repeats are also quite diverse and in line with ecological variables, such as host-dependent associations or population sizes, and with genetic variables, such as the recombination machinery. From an evolutionary point of view, repeats represent both opportunities and problems. We describe how repeats are created and how they can be found in genomes. We then focus on the functional and genomic consequences of repeats that dictate their fate.

Genome sequences of three agrobacterium biovars help elucidate the evolution of multichromosome genomes in bacteria

The family Rhizobiaceae contains plant-associated bacteria with critical roles in ecology and agriculture. Within this family, many Rhizobium and Sinorhizobium strains are nitrogen-fixing plant mutualists, while many strains designated as Agrobacterium are plant pathogens. These contrasting lifestyles are primarily dependent on the transmissible plasmids each strain harbors. Members of the Rhizobiaceae also have diverse genome architectures that include single chromosomes, multiple chromosomes, and plasmids of various sizes. Agrobacterium strains have been divided into three biovars, based on physiological and biochemical properties. The genome of a biovar I strain, A. tumefaciens C58, has been previously sequenced. In this study, the genomes of the biovar II strain A. radiobacter K84, a commercially available biological control strain that inhibits certain pathogenic agrobacteria, and the biovar III strain A. vitis S4, a narrow-host-range strain that infects grapes and invokes a hypersensitive response on nonhost plants, were fully sequenced and annotated. Comparison with other sequenced members of the Alphaproteobacteria provides new data on the evolution of multipartite bacterial genomes. Primary chromosomes show extensive conservation of both gene content and order. In contrast, secondary chromosomes share smaller percentages of genes, and conserved gene order is restricted to short blocks. We propose that secondary chromosomes originated from an ancestral plasmid to which genes have been transferred from a progenitor primary chromosome. Similar patterns are observed in select Beta- and Gammaproteobacteria species. Together, these results define the evolution of chromosome architecture and gene content among the Rhizobiaceae and support a generalized mechanism for second-chromosome formation among bacteria.

The organization of the bacterial genome

Many bacterial cellular processes interact intimately with the chromosome. Such interplay is the major driving force of genome structure or organization. Interactions take place at different scales-local for gene expression, global for replication-and lead to the differentiation of the chromosome into organizational units such as operons, replichores, or macrodomains. These processes are intermingled in the cell and create complex higher-level organizational features that are adaptive because they favor the interplay between the processes. The surprising result of selection for genome organization is that gene repertoires change much more quickly than chromosomal structure. Comparative genomics and experimental genomic manipulations are untangling the different cellular and evolutionary mechanisms causing such resilience to change. Since organization results from cellular processes, a better understanding of chromosome organization will help unravel the underlying cellular processes and their diversity.

Systematic identification of stem-loop containing sequence families in bacterial genomes

BACKGROUND

Analysis of non-coding sequences in several bacterial genomes brought to the identification of families of repeated sequences, able to fold as secondary structures. These sequences have often been claimed to be transcribed and fulfill a functional role. A previous systematic analysis of a representative set of 40 bacterial genomes produced a large collection of sequences, potentially able to fold as stem-loop structures (SLS). Computational analysis of these sequences was carried out by searching for families of repetitive nucleic acid elements sharing a common secondary structure.

RESULTS

The initial clustering procedure identified clusters of similar sequences in 29 genomes, corresponding to about 1% of the whole population. Sequences selected in this way have a substantially higher aptitude to fold into a stable secondary structure than the initial set. Removal of redundancies and regrouping of the selected sequences resulted in a final set of 92 families, defined by HMM analysis. 25 of them include all well-known SLS containing repeats and others reported in literature, but not analyzed in detail. The remaining 67 families have not been previously described. Two thirds of the families share a common predicted secondary structure and are located within intergenic regions.

CONCLUSION

Systematic analysis of 40 bacterial genomes revealed a large number of repeated sequence families, including known and novel ones. Their predicted structure and genomic location suggest that, even in compact bacterial genomes, a relatively large fraction of the genome consists of non-protein-coding sequences, possibly functioning at the RNA level.

Persistence drives gene clustering in bacterial genomes

BACKGROUND

Gene clustering plays an important role in the organization of the bacterial chromosome and several mechanisms have been proposed to explain its extent. However, the controversies raised about the validity of each of these mechanisms remind us that the cause of this gene organization remains an open question. Models proposed to explain clustering did not take into account the function of the gene products nor the likely presence or absence of a given gene in a genome. However, genomes harbor two very different categories of genes: those genes present in a majority of organisms - persistent genes - and those present in very few organisms - rare genes.

RESULTS

We show that two classes of genes are significantly clustered in bacterial genomes: the highly persistent and the rare genes. The clustering of rare genes is readily explained by the selfish operon theory. Yet, genes persistently present in bacterial genomes are also clustered and we try to understand why. We propose a model accounting specifically for such clustering, and show that indispensability in a genome with frequent gene deletion and insertion leads to the transient clustering of these genes. The model describes how clusters are created via the gene flux that continuously introduces new genes while deleting others. We then test if known selective processes, such as co-transcription, physical interaction or functional neighborhood, account for the stabilization of these clusters.

CONCLUSION

We show that the strong selective pressure acting on the function of persistent genes, in a permanent state of flux of genes in bacterial genomes, maintaining their size fairly constant, that drives persistent genes clustering. A further selective stabilization process might contribute to maintaining the clustering.