Conditions for the evolution of gene clusters in bacterial genomes

Noncontiguous operon atlas for the Staphylococcus aureus genome

Bacteria synchronize the expression of genes with related functions by organizing genes into operons so that they are cotranscribed together in a single polycistronic messenger RNA. However, some cellular processes may benefit if the simultaneous production of the operon proteins coincides with the inhibition of the expression of an antagonist gene. To coordinate such situations, bacteria have evolved noncontiguous operons (NcOs), a subtype of operons that contain one or more genes that are transcribed in the opposite direction to the other operon genes. This structure results in overlapping transcripts whose expression is mutually repressed. The presence of NcOs cannot be predicted computationally and their identification requires a detailed knowledge of the bacterial transcriptome. In this study, we used direct RNA sequencing methodology to determine the NcOs map in the Staphylococcus aureus genome. We detected the presence of 18 NcOs in the genome of S. aureus and four in the genome of the lysogenic prophage 80α. The identified NcOs comprise genes involved in energy metabolism, metal acquisition and transport, toxin-antitoxin systems, and control of the phage life cycle. Using the menaquinone operon as a proof of concept, we show that disarrangement of the NcO architecture results in a reduction of bacterial fitness due to an increase in menaquinone levels and a decrease in the rate of oxygen consumption. Our study demonstrates the significance of NcO structures in bacterial physiology and emphasizes the importance of combining operon maps with transcriptomic data to uncover previously unnoticed functional relationships between neighbouring genes.

Genome-guided prediction of acid resistance mechanisms in acidophilic methanotrophs of phylogenetically deep-rooted Verrucomicrobia isolated from geothermal environments

Verrucomicrobia are a group of microorganisms that have been proposed to be deeply rooted in the Tree of Life. Some are methanotrophs that oxidize the potent greenhouse gas methane and are thus important in decreasing atmospheric concentrations of the gas, potentially ameliorating climate change. They are widespread in various environments including soil and fresh or marine waters. Recently, a clade of extremely acidophilic Verrucomicrobia, flourishing at pH < 3, were described from high-temperature geothermal ecosystems. This novel group could be of interest for studies about the emergence of life on Earth and to astrobiologists as homologs for possible extraterrestrial life. In this paper, we describe predicted mechanisms for survival of this clade at low pH and suggest its possible evolutionary trajectory from an inferred neutrophilic ancestor. Extreme acidophiles are defined as organisms that thrive in extremely low pH environments (≤ pH 3). Many are polyextremophiles facing high temperatures and high salt as well as low pH. They are important to study for both providing fundamental insights into biological mechanisms of survival and evolution in such extreme environments and for understanding their roles in biotechnological applications such as industrial mineral recovery (bioleaching) and mitigation of acid mine drainage. They are also, potentially, a rich source of novel genes and pathways for the genetic engineering of microbial strains. Acidophiles of the Verrucomicrobia phylum are unique as they are the only known aerobic methanotrophs that can grow optimally under acidic (pH 2–3) and moderately thermophilic conditions (50–60°C). Three moderately thermophilic genera, namely Methylacidiphilum, Methylacidimicrobium, and Ca. Methylacidithermus, have been described in geothermal environments. Most of the investigations of these organisms have focused on their methane oxidizing capabilities (methanotrophy) and use of lanthanides as a protein cofactor, with no extensive study that sheds light on the mechanisms that they use to flourish at extremely low pH. In this paper, we extend the phylogenetic description of this group of acidophiles using whole genome information and we identify several mechanisms, potentially involved in acid resistance, including “first line of defense” mechanisms that impede the entry of protons into the cell. These include the presence of membrane-associated hopanoids, multiple copies of the outer membrane protein (Slp), and inner membrane potassium channels (kup, kdp) that generate a reversed membrane potential repelling the intrusion of protons. Acidophilic Verrucomicrobia also display a wide array of proteins potentially involved in the “second line of defense” where protons that evaded the first line of defense and entered the cell are expelled or neutralized, such as the glutamate decarboxylation (gadAB) and phosphate-uptake systems. An exclusive N-type ATPase F₀-F₁ was identified only in acidophiles of Verrucomicrobia and is predicted to be a specific adaptation in these organisms. Phylogenetic analyses suggest that many predicted mechanisms are evolutionarily conserved and most likely entered the acidophilic lineage of Verrucomicrobia by vertical descent from a common ancestor. However, it is likely that some defense mechanisms such as gadA and kup entered the acidophilic Verrucomicrobia lineage by horizontal gene transfer.

Gene cluster from plant to microbes: Their role in genome architecture, organism's development, specialized metabolism and drug discovery

Plants and microbes fulfil our daily requirements through different high-value chemicals, e.g., nutraceuticals, pharmaceuticals, cosmetics, and through varieties of fruits, crops, vegetables, and many more. Utmost care would therefore be taken for growth, development and sustainability of these important crops and medicinal plants and microbes. Homeobox genes and HOX clusters and their recently characterized expanded family members, including newly discovered homeobox, WOX gene from medicinal herb, Panax ginseng, significantly contributes in the growth and development of these organisms. On the other hand, secondary metabolites produced through secondary metabolism of plants and microbes are used as organisms defense as well as drugs/drug-like molecules for humans. Both the developmental HOX cluster and the biosynthetic gene-cluster (BGC) for secondary metabolites are organised in organisms genome. Genome mining and genomewide analysis of these clusters will definitely identify and characterize many more important molecules from unexplored plants and microbes and underexplored human microbiota and the evolution studies of these clusters will indicate their source of origin. Although genomics revolution now continues at a pace, till date only few hundred plant genome sequences are available. However, next-generation sequencing (NGS) technology now in market and may be applied even for plants with recalcitrant genomes, eventually may discover genomic potential towards production of secondary metabolites of diverse plants and micro-organisms present in the environment and microbiota. Additionally, the development of tools for genome mining e.g., antiSMASH, plantiSMASH, and more and more computational approaches that predicts hundreds of secondary metabolite BGCs will be discussed.

Phylogenetic analysis of the bacterial Pro-Pro-endopeptidase domain reveals a diverse family including secreted and membrane anchored proteins

•

Bacterial Pro-Pro-endopeptidase (PPEP) is the latest member of the metalloendopeptidase class (E.C. 3.4.24.89).

•

PPEP homologs are found in two firmicutes orders, clostridiales and bacillales spread over 9 genera and more than 130 species.

•

Some PPEP homologs have acquired additional anchor domains that bind noncovalently to various elements of the bacterial peptidoglycan cell wall.

•

Prototype family members, PPEP-1 and PPEP-2, target bacterial surface adhesion proteins, but homologs could target other extracellular proteins.

Pro-Pro-endopeptidases (PPEP, EC 3.4.24.89) are secreted, zinc metalloproteases that have the unusual capacity to cleave a peptide bond between two prolines, a bond that is generally less sensitive to proteolytic cleavage. Two well studied members of the family are PPEP-1 and PPEP-2, produced by Clostridioides difficile, a human pathogen, and Paenibacillus alvei, a bee secondary invader, respectively. Both proteases seem to be involved in mediating bacterial adhesion by cleaving cell surface anchor proteins on the bacterium itself.

By using basic alignment and phylogenetic profiling analysis, this work shows that the complete family of proteins that contain a PPEP domain includes proteins from more than 130 species spread over 9 genera. These analyses also suggest that the PPEP domain spread through horizontal gene transfer events between species within the Firmicutes’ classes Bacilli and Clostridia. Bacterial species containing PPEP homologs are found in diverse habitats, varying from human pathogens and gut microbiota to free-living bacteria, which were isolated from various environments, including extreme conditions such as hot springs, desert soil and salt lakes.

The phylogenetic tree reveals the relationships between family members and suggests that smaller subgroups could share cleavage specificity, substrates and functional similarity. Except for PPEP-1 and PPEP-2, no cleavage specificity, specific physiological target, or function has been assigned for any of the other PPEP-family members. Some PPEP proteins have acquired additional domains that recognize and bind noncovalently to various elements of the bacterial peptidoglycan cell-wall, anchoring these PPEPs. Secreted or anchored to the cell-wall surface PPEP proteins seem to perform various functions.

Bacterial cooperation through horizontal gene transfer

Cooperation exists across all scales of biological organization, from genetic elements to complex human societies. Bacteria cooperate by secreting molecules that benefit all individuals in the population (i.e., public goods). Genes associated with cooperation can spread among strains through horizontal gene transfer (HGT). We discuss recent findings on how HGT mediated by mobile genetic elements promotes bacterial cooperation, how cooperation in turn can facilitate more frequent HGT, and how the act of HGT itself may be considered as a form of cooperation. We propose that HGT is an important enforcement mechanism in bacterial populations, thus creating a positive feedback loop that further maintains cooperation. To enforce cooperation, HGT serves as a homogenizing force by transferring the cooperative trait, effectively eliminating cheaters.

Mining genomes to illuminate the specialized chemistry of life

All organisms produce specialized organic molecules, ranging from small volatile chemicals to large gene-encoded peptides, that have evolved to provide them with diverse cellular and ecological functions. As natural products, they are broadly applied in medicine, agriculture and nutrition. The rapid accumulation of genomic information has revealed that the metabolic capacity of virtually all organisms is vastly underappreciated. Pioneered mainly in bacteria and fungi, genome mining technologies are accelerating metabolite discovery. Recent efforts are now being expanded to all life forms, including protists, plants and animals, and new integrative omics technologies are enabling the increasingly effective mining of this molecular diversity.

Plant Metabolic Gene Clusters: Evolution, Organization, and Their Applications in Synthetic Biology

Plants are a remarkable source of high-value specialized metabolites having significant physiological and ecological functions. Genes responsible for synthesizing specialized metabolites are often clustered together for a coordinated expression, which is commonly observed in bacteria and filamentous fungi. Similar to prokaryotic gene clustering, plants do have gene clusters encoding enzymes involved in the biosynthesis of specialized metabolites. More than 20 gene clusters involved in the biosynthesis of diverse metabolites have been identified across the plant kingdom. Recent studies demonstrate that gene clusters are evolved through gene duplications and neofunctionalization of primary metabolic pathway genes. Often, these clusters are tightly regulated at nucleosome level. The prevalence of gene clusters related to specialized metabolites offers an attractive possibility of an untapped source of highly useful biomolecules. Accordingly, the identification and functional characterization of novel biosynthetic pathways in plants need to be worked out. In this review, we summarize insights into the evolution of gene clusters and discuss the organization and importance of specific gene clusters in the biosynthesis of specialized metabolites. Regulatory mechanisms which operate in some of the important gene clusters have also been briefly described. Finally, we highlight the importance of gene clusters to develop future metabolic engineering or synthetic biology strategies for the heterologous production of novel metabolites.

Treat Me Well or Will Resist: Uptake of Mobile Genetic Elements Determine the Resistome of Corynebacterium striatum

Corynebacterium striatum, a bacterium that is part of the normal skin microbiota, is also an opportunistic pathogen. In recent years, reports of infections and in-hospital and nosocomial outbreaks caused by antimicrobial multidrug-resistant C. striatum strains have been increasing worldwide. However, there are no studies about the genomic determinants related to antimicrobial resistance in C. striatum. This review updates global information related to antimicrobial resistance found in C. striatum and highlights the essential genomic aspects in its persistence and dissemination. The resistome of C. striatum comprises chromosomal and acquired elements. Resistance to fluoroquinolones and daptomycin are due to mutations in chromosomal genes. Conversely, resistance to macrolides, tetracyclines, phenicols, beta-lactams, and aminoglycosides are associated with mobile genomic elements such as plasmids and transposons. The presence and diversity of insertion sequences suggest an essential role in the expression of antimicrobial resistance genes (ARGs) in genomic rearrangements and their potential to transfer these elements to other pathogens. The present study underlines that the resistome of C. striatum is dynamic; it is in evident expansion and could be acting as a reservoir for ARGs.

The Association between Insertion Sequences and Antibiotic Resistance Genes

The emergence and spread of antibiotic resistance genes (ARGs) among pathogens threaten the prevention and treatment of bacterial infections as well as our food production chains. Early knowledge about mobile ARGs that are present in pathogens or that have the potential to become clinically relevant could help mitigate potential negative consequences. Recently, exploring integron gene cassettes was shown to be successful for identifying novel mobilized ARGs, some of which were already circulating in pathogens. Still, only a subset of ARGs is mobilized by integrons, and the contexts of other mobile genetic elements associated with ARGs remain unexplored. This includes insertion sequences (ISs) responsible for the mobilization of many ARGs. Our analyses identified ISs, species, and environments where ARG-IS relationships are particularly strong. This could be a first step to guide the discovery of novel ARGs, while also providing insights into mechanisms involved in the mobilization and transfer of ARGs.

Insertion sequences (ISs) are abundant mobile genetic elements on bacterial genomes, responsible for mobilization of many genes, including antibiotic resistance genes (ARGs). As ARGs often occur in similar genetic contexts, understanding which ISs tend to be associated with known ARGs could be a first step toward discovering novel ARGs through predictive or experimental strategies. This could be valuable, as early identification of ARGs in pathogens could facilitate surveillance, confinement actions, molecular diagnostics, and drug development. Here, we present a comprehensive analysis of the association of specific ISs with known ARGs. A large collection of bacterial genomes was used to characterize the immediate context of 2,437 known ARGs and 3,768 ISs. While many ARGs were consistently found close to specific ISs, the contexts around all ISs were more variable. Nevertheless, a subset of individual ISs, as well as tentative composite transposons, showed significant associations with ARGs. These included, e.g., insertion sequences classified as IS6, Tn3, IS4, and IS1 that were not only strongly associated with diverse ARGs but also highly abundant in pathogens. Therefore, we conclude that the context of this subset of ISs and tentative composite transposons would be particularly valuable to discover novel ARGs through modeling or empirical approaches. A set of 1,891 metagenomes were analyzed to identify environments where those ISs commonly associated with ARGs were particularly abundant. The associations found in metagenomes were similar to those found in genomes.

IMPORTANCE The emergence and spread of antibiotic resistance genes (ARGs) among pathogens threaten the prevention and treatment of bacterial infections as well as our food production chains. Early knowledge about mobile ARGs that are present in pathogens or that have the potential to become clinically relevant could help mitigate potential negative consequences. Recently, exploring integron gene cassettes was shown to be successful for identifying novel mobilized ARGs, some of which were already circulating in pathogens. Still, only a subset of ARGs is mobilized by integrons, and the contexts of other mobile genetic elements associated with ARGs remain unexplored. This includes insertion sequences (ISs) responsible for the mobilization of many ARGs. Our analyses identified ISs, species, and environments where ARG-IS relationships are particularly strong. This could be a first step to guide the discovery of novel ARGs, while also providing insights into mechanisms involved in the mobilization and transfer of ARGs.

Host-glycan metabolism is regulated by a species-conserved two-component system in Streptococcus pneumoniae

Pathogens of the Streptococcus genus inhabit many different environmental niches during the course of an infection in a human host and the bacteria must adjust their metabolism according to available nutrients. Despite their lack of the citric-acid cycle, some streptococci proliferate in niches devoid of a readily available carbohydrate source. Instead they rely on carbohydrate scavenging for energy acquisition, which are obtained from the host. Here we discover a two-component system (TCS07) of Streptococcus pneumoniae that responds to glycoconjugated structures on proteins present on the host cells. Using next-generation RNA sequencing we find that the uncharacterized TCS07 regulon encodes proteins important for host-glycan processing and transporters of the released glycans, as well as intracellular carbohydrate catabolizing enzymes. We find that a functional TCS07 allele is required for growth on the glycoconjugated model protein fetuin. Consistently, we see a TCS07-dependent activation of the glycan degradation pathway. Thus, we pinpoint the molecular constituents responsible for sensing host derived glycans and link this to the induction of the proteins necessary for glycan degradation. Furthermore, we connect the TCS07 regulon to virulence in a mouse model, thereby establishing that host-derived glycan-metabolism is important for infection in vivo. Finally, a comparative phylogenomic analysis of strains from the Streptococcus genus reveal that TCS07 and most of its regulon is specifically conserved in species that utilize host-glycans for growth.

Worldwide, Streptococcus pneumoniae is the most common cause of community acquired pneumonia with high mortality rates. Interestingly, S. pneumoniae strictly relies on carbohydrate scavenging for energy acquisition, which are obtained from the host. This is a critical step in pathogenesis and a common mechanism among Streptococcal species. In this study, we discover an uncharacterized two-component system that responds to the carbohydrate structures present on the host cells. These are important findings as we describe the molecular mechanism responsible for sensing these host derived glycans, and how this mechanism is linked to virulence, thus highlighting that glycan metabolism is important for infection in vivo, thereby posing a novel target for intervention. Our phylogenetic analysis reveals that the two-component system and the genetic regulon co-occur and are specifically conserved among Streptococcal species capable of degrading host-glycans.

Transcriptome analysis of Xanthomonas oryzae pv. oryzicola exposed to H2O2 reveals horizontal gene transfer contributes to its oxidative stress response

Xanthomonas oryzae pv. oryzicola (Xoc), the causal agent of bacterial leaf streak, is one of the most severe seed-borne bacterial diseases of rice. However, the molecular mechanisms underlying Xoc in response to oxidative stress are still unknown. In this study, we performed a time-course RNA-seq analysis on the Xoc in response to H2O2, aiming to reveal its oxidative response network. Overall, our RNA sequence analysis of Xoc revealed a significant global gene expression profile when it was exposed to H2O2. There were 7, 177, and 246 genes that were differentially regulated at the early, middle, and late stages after exposure, respectively. Three genes (xoc_1643, xoc_1946, xoc_3249) showing significantly different expression levels had proven relationships with oxidative stress response and pathogenesis. Moreover, a hypothetical protein (XOC_2868) showed significantly differential expression, and the xoc_2868 mutants clearly displayed a greater H2O2 sensitivity and decreased pathogenicity than those of the wild-type. Gene localization and phylogeny analysis strongly suggests that this gene may have been horizontally transferred from a Burkholderiaceae ancestor. Our study not only provides a first glance of Xoc's global response against oxidative stress, but also reveals the impact of horizontal gene transfer in the evolutionary history of Xoc.

Advances in microbial culturing conditions to activate silent biosynthetic gene clusters for novel metabolite production

Natural products (NPs) produced by bacteria and fungi are often used as therapeutic agents due to their complex structures and wide range of bioactivities. Enzymes that build NPs are encoded by co-localized biosynthetic gene clusters (BGCs), and genome sequencing has recently revealed that many BGCs are “silent” under standard laboratory conditions. There are numerous methods used to activate “silent” BGCs that rely either upon altering culture conditions or genetic modification. In this review, we discuss several recent microbial cultivation methods that have been used to expand the scope of NPs accessible in the laboratory. These approaches are divided into three categories: addition of a physical scaffold, addition of small molecule elicitors, and co-cultivation with another microbe.

The predominance of nucleotidyl activation in bacterial phosphonate biosynthesis

Phosphonates are rare and unusually bioactive natural products. However, most bacterial phosphonate biosynthetic capacity is dedicated to tailoring cell surfaces with molecules like 2-aminoethylphosphonate (AEP). Although phosphoenolpyruvate mutase (Ppm)-catalyzed installation of C-P bonds is known, subsequent phosphonyl tailoring (Pnt) pathway steps remain enigmatic. Here we identify nucleotidyltransferases in over two-thirds of phosphonate biosynthetic gene clusters, including direct fusions to ~60% of Ppm enzymes. We characterize two putative phosphonyl tailoring cytidylyltransferases (PntCs) that prefer AEP over phosphocholine (P-Cho) – a similar substrate used by the related enzyme LicC, which is a virulence factor in Streptococcus pneumoniae. PntC structural analyses reveal steric discrimination against phosphocholine. These findings highlight nucleotidyl activation as a predominant chemical logic in phosphonate biosynthesis and set the stage for probing diverse phosphonyl tailoring pathways.

Phosphonate modifications can be present on microbial cell surfaces. Here the authors perform bioinformatics analyses and observe a widespread occurrence of nucleotidyltransferase-encoding genes in bacterial phosphonate biosynthesis and functionally characterize two of the identified phosphonate specific cytidylyltransferases (PntCs) and determine the crystal structure of T. denticola PntC.

Motion, fixation probability and the choice of an evolutionary process

Seemingly minor details of mathematical and computational models of evolution are known to change the effect of population structure on the outcome of evolutionary processes. For example, birth-death dynamics often result in amplification of selection, while death-birth processes have been associated with suppression. In many biological populations the interaction structure is not static. Instead, members of the population are in motion and can interact with different individuals at different times. In this work we study populations embedded in a flowing medium; the interaction network is then time dependent. We use computer simulations to investigate how this dynamic structure affects the success of invading mutants, and compare these effects for different coupled birth and death processes. Specifically, we show how the speed of the motion impacts the fixation probability of an invading mutant. Flows of different speeds interpolate between evolutionary dynamics on fixed heterogeneous graphs and well-stirred populations; this allows us to systematically compare against known results for static structured populations. We find that motion has an active role in amplifying or suppressing selection by fragmenting and reconnecting the interaction graph. While increasing flow speeds suppress selection for most evolutionary models, we identify characteristic responses to flow for the different update rules we test. In particular we find that selection can be maximally enhanced or suppressed at intermediate flow speeds.

Whether a mutation spreads in a population or not is one of the most important questions in biology. The evolution of cancer and antibiotic resistance, for example, are mediated by invading mutants. Recent work has shown that population structure can have important consequences for the outcome of evolution. For instance, a mutant can have a higher or a lower chance of invasion than in unstructured populations. These effects can depend on seemingly minor details of the evolutionary model, such as the order of birth and death events. Many biological populations are in motion, for example due to external stirring. Experimentally this is known to be important; the performance of mutants in E. coli populations, for example, depends on the rate of mixing. Here, we focus on simulations of populations in a flowing medium, and compare the success of a mutant for different flow speeds. We contrast different evolutionary models, and identify what features of the evolutionary model affect mutant success for different speeds of the flow. We find that the chance of mutant invasion can be at its highest (or lowest) at intermediate flow speeds, depending on the order in which birth and death events occur in the evolutionary process.

Gene Gangs of the Chloroviruses: Conserved Clusters of Collinear Monocistronic Genes

Chloroviruses (family Phycodnaviridae) are dsDNA viruses found throughout the world's inland waters. The open reading frames in the genomes of 41 sequenced chloroviruses (330 ± 40 kbp each) representing three virus types were analyzed for evidence of evolutionarily conserved local genomic "contexts", the organization of biological information into units of a scale larger than a gene. Despite a general loss of synteny between virus types, we informatically detected a highly conserved genomic context defined by groups of three or more genes that we have termed "gene gangs". Unlike previously described local genomic contexts, the definition of gene gangs requires only that member genes be consistently co-localized and are not constrained by strand, regulatory sites, or intervening sequences (and therefore represent a new type of conserved structural genomic element). An analysis of functional annotations and transcriptomic data suggests that some of the gene gangs may organize genes involved in specific biochemical processes, but that this organization does not involve their coordinated expression.

FLIM-MAP: Gene Context Based Identification of Functional Modules in Bacterial Metabolic Pathways

Prediction of functional potential of bacteria can only be ascertained by the accurate annotation of its metabolic pathways. Homology based methods decipher metabolic gene content but ignore the fact that homologs of same protein can function in different pathways. Therefore, mere presence of all constituent genes in an organism is not sufficient to indicate a pathway. Contextual occurrence of genes belonging to a pathway on the bacterial genome can hence be exploited for an accurate estimation of functional potential of a bacterium. In this communication, we present a novel annotation resource to accurately identify pathway presence by using gene context. Our tool FLIM-MAP (Functionally Important Modules in bacterial Metabolic Pathways) predicts biologically relevant functional units called ‘GCMs’ (Gene Context based Modules) from a given metabolic reaction network. We benchmark the accuracy of our tool on amino acids and carbohydrate metabolism pathways.

Evolution of Sphingomonad Gene Clusters Related to Pesticide Catabolism Revealed by Genome Sequence and Mobilomics of Sphingobium herbicidovorans MH

Bacterial degraders of chlorophenoxy herbicides have been isolated from various ecosystems, including pristine environments. Among these degraders, the sphingomonads constitute a prominent group that displays versatile xenobiotic-degradation capabilities. Four separate sequencing strategies were required to provide the complete sequence of the complex and plastic genome of the canonical chlorophenoxy herbicide-degrading Sphingobium herbicidovorans MH. The genome has an intricate organization of the chlorophenoxy-herbicide catabolic genes sdpA, rdpA, and cadABCD that encode the (R)- and (S)-enantiomer-specific 2,4-dichlorophenoxypropionate dioxygenases and four subunits of a Rieske non-heme iron oxygenase involved in 2-methyl-chlorophenoxyacetic acid degradation, respectively. Several major genomic rearrangements are proposed to help understand the evolution and mobility of these important genes and their genetic context. Single-strain mobilomic sequence analysis uncovered plasmids and insertion sequence-associated circular intermediates in this environmentally important bacterium and enabled the description of evolutionary models for pesticide degradation in strain MH and related organisms. The mobilome presented a complex mosaic of mobile genetic elements including four plasmids and several circular intermediate DNA molecules of insertion-sequence elements and transposons that are central to the evolution of xenobiotics degradation. Furthermore, two individual chromosomally integrated prophages were shown to excise and form free circular DNA molecules. This approach holds great potential for improving the understanding of genome plasticity, evolution, and microbial ecology.

The evolutionary life cycle of the polysaccharide biosynthetic gene cluster based on the Sphingomonadaceae

Although clustering of genes from the same metabolic pathway is a widespread phenomenon, the evolution of the polysaccharide biosynthetic gene cluster remains poorly understood. To determine the evolution of this pathway, we identified a scattered production pathway of the polysaccharide sanxan by Sphingomonas sanxanigenens NX02, and compared the distribution of genes between sphingan-producing and other Sphingomonadaceae strains. This allowed us to determine how the scattered sanxan pathway developed, and how the polysaccharide gene cluster evolved. Our findings suggested that the evolution of microbial polysaccharide biosynthesis gene clusters is a lengthy cyclic process comprising cluster 1 → scatter → cluster 2. The sanxan biosynthetic pathway proved the existence of a dispersive process. We also report the complete genome sequence of NX02, in which we identified many unstable genetic elements and powerful secretion systems. Furthermore, nine enzymes for the formation of activated precursors, four glycosyltransferases, four acyltransferases, and four polymerization and export proteins were identified. These genes were scattered in the NX02 genome, and the positive regulator SpnA of sphingans synthesis could not regulate sanxan production. Finally, we concluded that the evolution of the sanxan pathway was independent. NX02 evolved naturally as a polysaccharide producing strain over a long-time evolution involving gene acquisitions and adaptive mutations.

Conserved Units of Co-Expression in Bacterial Genomes: An Evolutionary Insight into Transcriptional Regulation

Genome-wide measurements of transcriptional activity in bacteria indicate that the transcription of successive genes is strongly correlated beyond the scale of operons. Here, we analyze hundreds of bacterial genomes to identify supra-operonic segments of genes that are proximal in a large number of genomes. We show that these synteny segments correspond to genomic units of strong transcriptional co-expression. Structurally, the segments contain operons with specific relative orientations (co-directional or divergent) and nucleoid-associated proteins are found to bind at their boundaries. Functionally, operons inside a same segment are highly co-expressed even in the apparent absence of regulatory factors at their promoter regions. Remote operons along DNA can also be co-expressed if their corresponding segments share a transcriptional or sigma factor, without requiring these factors to bind directly to the promoters of the operons. As evidence that these results apply across the bacterial kingdom, we demonstrate them both in the Gram-negative bacterium Escherichia coli and in the Gram-positive bacterium Bacillus subtilis. The underlying process that we propose involves only RNA-polymerases and DNA: it implies that the transcription of an operon mechanically enhances the transcription of adjacent operons. In support of a primary role of this regulation by facilitated co-transcription, we show that the transcription en bloc of successive operons as a result of transcriptional read-through is strongly and specifically enhanced in synteny segments. Finally, our analysis indicates that facilitated co-transcription may be evolutionary primitive and may apply beyond bacteria.

Genome rearrangements can make and break small RNA genes

Small RNAs (sRNAs) are short, transcribed regulatory elements that are typically encoded in the intergenic regions (IGRs) of bacterial genomes. Several sRNAs, first recognized in Escherichia coli, are conserved among enteric bacteria, but because of the regulatory roles of sRNAs, differences in sRNA repertoires might be responsible for features that differentiate closely related species. We scanned the E. coli MG1655 and Salmonella enterica Typhimurium genomes for nonsyntenic IGRs as a potential source of uncharacterized, species-specific sRNAs and found that genome rearrangements have reconfigured several IGRs causing the disruption and formation of sRNAs. Within an IGR that is present in E. coli but was disrupted in Salmonella by a translocation event is an sRNA that is associated with the FNR/CRP global regulators and influences E. coli biofilm formation. A Salmonella-specific sRNA evolved de novo through point mutations that generated a σ⁷⁰ promoter sequence in an IGR that arose through genome rearrangement events. The differences in the sRNA pools among bacterial species have previously been ascribed to duplication, deletion, or horizontal acquisition. Here, we show that genomic rearrangements also contribute to this process by either disrupting sRNA-containing IGRs or creating IGRs in which novel sRNAs may evolve.

Sioxanthin, a novel glycosylated carotenoid, reveals an unusual subclustered biosynthetic pathway

Members of the marine actinomycete genus Salinispora constitutively produce a characteristic orange pigment during vegetative growth. Contrary to the understanding of widespread carotenoid biosynthesis pathways in bacteria, Salinispora carotenoid biosynthesis genes are not confined to a single cluster. Instead, bioinformatic and genetic investigations confirm that four regions of the Salinispora tropica CNB-440 genome, consisting of two gene clusters and two independent genes, contribute to the in vivo production of a single carotenoid. This compound, namely (2'S)-1'-(β-D-glucopyranosyloxy)-3',4'-didehydro-1',2'-dihydro-φ,ψ-caroten-2'-ol, is novel and has been given the trivial name 'sioxanthin'. Sioxanthin is a C40 -carotenoid, glycosylated on one end of the molecule and containing an aryl moiety on the opposite end. Glycosylation is unusual among actinomycete carotenoids, and sioxanthin joins a rare group of carotenoids with polar and non-polar head groups. Gene sequence homology predicts that the sioxanthin biosynthetic pathway is present in all of the Salinispora as well as other members of the family Micromonosporaceae. Additionally, this study's investigations of clustering of carotenoid biosynthetic genes in heterotrophic bacteria show that a non-clustered genome arrangement is more common than previously suggested, with nearly half of the investigated genomes showing a non-clustered architecture.

Plastid-LCGbase: a collection of evolutionarily conserved plastid-associated gene pairs

Plastids carry their own genetic material that encodes a variable set of genes that are limited in number but functionally important. Aside from orthology, the lineage-specific order and orientation of these genes are also relevant. Here, we develop a database, Plastid-LCGbase (http://lcgbase.big.ac.cn/plastid-LCGbase/), which focuses on organizational variability of plastid genes and genomes from diverse taxonomic groups. The current Plastid-LCGbase contains information from 470 plastid genomes and exhibits several unique features. First, through a genome-overview page generated from OrganellarGenomeDRAW, it displays general arrangement of all plastid genes (circular or linear). Second, it shows patterns and modes of all paired plastid genes and their physical distances across user-defined lineages, which are facilitated by a step-wise stratification of taxonomic groups. Third, it divides the paired genes into three categories (co-directionally-paired genes or CDPGs, convergently-paired genes or CPGs and divergently-paired genes or DPGs) and three patterns (separation, overlap and inclusion) and provides basic statistics for each species. Fourth, the gene pairing scheme is expandable, where neighboring genes can also be included in species-/lineage-specific comparisons. We hope that Plastid-LCGbase facilitates gene variation (insertion-deletion, translocation and rearrangement) and transcription-level studies of plastid genomes.

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.

Integrative data-mining tools to link gene and function

Information derived from genomic and post-genomic data can be efficiently used to link gene and function. Several web-based platforms have been developed to mine these types of data by integrating different tools. This method paper is designed to allow the user to navigate these platforms in order to make functional predictions. The main focus is on phylogenetic distribution and physical clustering tools, but other tools such as pathway reconstruction, gene fusions, and analysis of high-throughput experimental data are also surveyed.

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.

Bacterial genome instability

Bacterial genomes are remarkably stable from one generation to the next but are plastic on an evolutionary time scale, substantially shaped by horizontal gene transfer, genome rearrangement, and the activities of mobile DNA elements. This implies the existence of a delicate balance between the maintenance of genome stability and the tolerance of genome instability. In this review, we describe the specialized genetic elements and the endogenous processes that contribute to genome instability. We then discuss the consequences of genome instability at the physiological level, where cells have harnessed instability to mediate phase and antigenic variation, and at the evolutionary level, where horizontal gene transfer has played an important role. Indeed, this ability to share DNA sequences has played a major part in the evolution of life on Earth. The evolutionary plasticity of bacterial genomes, coupled with the vast numbers of bacteria on the planet, substantially limits our ability to control disease.

Genetic Diversity of Tick-Borne Rickettsial Pathogens; Insights Gained from Distant Strains

The ability to capture genetic variation with unprecedented resolution improves our understanding of bacterial populations and their ability to cause disease. The goal of the pathogenomics era is to define genetic diversity that results in disease. Despite the economic losses caused by vector-borne bacteria in the Order Rickettsiales, little is known about the genetic variants responsible for observed phenotypes. The tick-transmitted rickettsial pathogen Anaplasma marginale infects cattle in tropical and subtropical regions worldwide, including Australia. Genomic analysis of North American A. marginale strains reveals a closed core genome defined by high levels of Single Nucleotide Polymorphisms (SNPs). Here we report the first genome sequences and comparative analysis for Australian strains that differ in virulence and transmissibility. A list of genetic differences that segregate with phenotype was evaluated for the ability to distinguish the attenuated strain from virulent field strains. Phylogenetic analyses of the Australian strains revealed a marked evolutionary distance from all previously sequenced strains. SNP analysis showed a strikingly reduced genetic diversity between these strains, with the smallest number of SNPs detected between any two A. marginale strains. The low diversity between these phenotypically distinct bacteria presents a unique opportunity to identify the genetic determinants of virulence and transmission.

Inference of phenotype-defining functional modules of protein families for microbial plant biomass degraders

BACKGROUND

Efficient industrial processes for converting plant lignocellulosic materials into biofuels are a key to global efforts to come up with alternative energy sources to fossil fuels. Novel cellulolytic enzymes have been discovered in microbial genomes and metagenomes of microbial communities. However, the identification of relevant genes without known homologs, and the elucidation of the lignocellulolytic pathways and protein complexes for different microorganisms remain challenging.

RESULTS

We describe a new computational method for the targeted discovery of functional modules of plant biomass-degrading protein families, based on their co-occurrence patterns across genomes and metagenome datasets, and the strength of association of these modules with the genomes of known degraders. From approximately 6.4 million family annotations for 2,884 microbial genomes, and 332 taxonomic bins from 18 metagenomes, we identified 5 functional modules that are distinctive for plant biomass degraders, which we term "plant biomass degradation modules" (PDMs). These modules incorporate protein families involved in the degradation of cellulose, hemicelluloses, and pectins, structural components of the cellulosome, and additional families with potential functions in plant biomass degradation. The PDMs were linked to 81 gene clusters in genomes of known lignocellulose degraders, including previously described clusters of lignocellulolytic genes. On average, 70% of the families of each PDM were found to map to gene clusters in known degraders, which served as an additional confirmation of their functional relationships. The presence of a PDM in a genome or taxonomic metagenome bin furthermore allowed us to accurately predict the ability of any particular organism to degrade plant biomass. For 15 draft genomes of a cow rumen metagenome, we used cross-referencing to confirmed cellulolytic enzymes to validate that the PDMs identified plant biomass degraders within a complex microbial community.

CONCLUSIONS

Functional modules of protein families that are involved in different aspects of plant cell wall degradation can be inferred from co-occurrence patterns across (meta-)genomes with a probabilistic topic model. PDMs represent a new resource of protein families and candidate genes implicated in microbial plant biomass degradation. They can also be used to predict the plant biomass degradation ability for a genome or taxonomic bin. The method is also suitable for characterizing other microbial phenotypes.

Metagenomic insights into the carbohydrate-active enzymes carried by the microorganisms adhering to solid digesta in the rumen of cows

The ruminal microbial community is a unique source of enzymes that underpin the conversion of cellulosic biomass. In this study, the microbial consortia adherent on solid digesta in the rumen of Jersey cattle were subjected to an activity-based metagenomic study to explore the genetic diversity of carbohydrolytic enzymes in Jersey cows, with a particular focus on cellulases and xylanases. Pyrosequencing and bioinformatic analyses of 120 carbohydrate-active fosmids identified genes encoding 575 putative Carbohydrate-Active Enzymes (CAZymes) and proteins putatively related to transcriptional regulation, transporters, and signal transduction coupled with polysaccharide degradation and metabolism. Most of these genes shared little similarity to sequences archived in databases. Genes that were predicted to encode glycoside hydrolases (GH) involved in xylan and cellulose hydrolysis (e.g., GH3, 5, 9, 10, 39 and 43) were well represented. A new subfamily (S-8) of GH5 was identified from contigs assigned to Firmicutes. These subfamilies of GH5 proteins also showed significant phylum-dependent distribution. A number of polysaccharide utilization loci (PULs) were found, and two of them contained genes encoding Sus-like proteins and cellulases that have not been reported in previous metagenomic studies of samples from the rumens of cows or other herbivores. Comparison with the large metagenomic datasets previously reported of other ruminant species (or cattle breeds) and wallabies showed that the rumen microbiome of Jersey cows might contain differing CAZymes. Future studies are needed to further explore how host genetics and diets affect the diversity and distribution of CAZymes and utilization of plant cell wall materials.

Evolutionary concepts in natural products discovery: what actinomycetes have taught us

Actinomycetes are a very important source of natural products for the pharmaceutical industry and other applications. Most of the strains belong to Streptomyces or related genera, partly because they are particularly amenable to growth in the laboratory and industrial fermenters. It is unlikely that chemical synthesis can fulfil the needs of the pharmaceutical industry for novel compounds so there is a continuing need to find novel natural products. An evolutionary perspective can help this process in several ways. Genome mining attempts to identify secondary metabolite biosynthetic clusters in DNA sequences, which are likely to produce interesting chemical entities. There are often technical problems in assembling the DNA sequences of large modular clusters in genome and metagenome projects, which can be overcome partially using information about the evolution of the domain sequences. Understanding the evolutionary mechanisms of modular clusters should allow simulation of evolutionary pathways in the laboratory to generate novel compounds.

Genome-wide association study identifies vitamin B5 biosynthesis as a host specificity factor in Campylobacter

Genome-wide association studies have the potential to identify causal genetic factors underlying important phenotypes but have rarely been performed in bacteria. We present an association mapping method that takes into account the clonal population structure of bacteria and is applicable to both core and accessory genome variation. Campylobacter is a common cause of human gastroenteritis as a consequence of its proliferation in multiple farm animal species and its transmission via contaminated meat and poultry. We applied our association mapping method to identify the factors responsible for adaptation to cattle and chickens among 192 Campylobacter isolates from these and other host sources. Phylogenetic analysis implied frequent host switching but also showed that some lineages were strongly associated with particular hosts. A seven-gene region with a host association signal was found. Genes in this region were almost universally present in cattle but were frequently absent in isolates from chickens and wild birds. Three of the seven genes encoded vitamin B5 biosynthesis. We found that isolates from cattle were better able to grow in vitamin B5-depleted media and propose that this difference may be an adaptation to host diet.

Diversification of the Salmonella fimbriae: a model of macro- and microevolution

Bacteria of the genus Salmonella comprise a large and evolutionary related population of zoonotic pathogens that can infect mammals, including humans and domestic animals, birds, reptiles and amphibians. Salmonella carries a plethora of virulence genes, including fimbrial adhesins, some of them known to participate in mammalian or avian host colonization. Each type of fimbria has its structural subunit and biogenesis genes encoded by one fimbrial gene cluster (FGC). The accumulation of new genomic information offered a timely opportunity to better evaluate the number and types of FGCs in the Salmonella pangenome, to test the use of current classifications based on phylogeny, and to infer potential correlations between FGC evolution in various Salmonella serovars and host niches. This study focused on the FGCs of the currently deciphered 90 genomes and 60 plasmids of Salmonella. The analysis highlighted a fimbriome consisting of 35 different FGCs, of which 16 were new, each strain carrying between 5 and 14 FGCs. The Salmonella fimbriome was extremely diverse with FGC representatives in 8 out of 9 previously categorized fimbrial clades and subclades. Phylogenetic analysis of Salmonella suggested macroevolutionary shifts detectable by extensive FGC deletion and acquisition. In addition, microevolutionary drifts were best depicted by the high level of allelic variation in predicted or known adhesins, such as the type 1 fimbrial adhesin FimH for which 67 different natural alleles were identified in S. enterica subsp. I. Together with strain-specific collections of FGCs, allelic variation among adhesins attested to the pathoadaptive evolution of Salmonella towards specific hosts and tissues, potentially modulating host range, strain virulence, disease progression, and transmission efficiency. Further understanding of how each Salmonella strain utilizes its panel of FGCs and specific adhesin alleles for survival and infection will support the development of new approaches for the control of Salmonellosis.

Evolution of variation in presence and absence of genes in bacterial pathways

Background

Bacterial genomes exhibit a remarkable degree of variation in the presence and absence of genes, which probably extends to the level of individual pathways. This variation may be a consequence of the significant evolutionary role played by horizontal gene transfer, but might also be explained by the loss of genes through mutation. A challenge is to understand why there would be variation in gene presence within pathways if they confer a benefit only when complete.

Results

Here, we develop a mathematical model to study how variation in pathway content is produced by horizontal transfer, gene loss and partial exposure of a population to a novel environment.

Conclusions

We discuss the possibility that variation in gene presence acts as cryptic genetic variation on which selection acts when the appropriate environment occurs. We find that a high level of variation in gene presence can be readily explained by decay of the pathway through mutation when there is no longer exposure to the selective environment, or when selection becomes too weak to maintain the genes. In the context of pathway variation the role of horizontal gene transfer is probably the initial introduction of a complete novel pathway rather than in building up the variation in a genome without the pathway.

Conservation of gene order and content in the circular chromosomes of 'Candidatus Liberibacter asiaticus' and other Rhizobiales

'Ca. Liberibacter asiaticus,' an insect-vectored, obligate intracellular bacterium associated with citrus-greening disease, also called "HLB," is a member of the Rhizobiales along with nitrogen-fixing microsymbionts Sinorhizobium meliloti and Bradyrhizobium japonicum, plant pathogen Agrobacterium tumefaciens and facultative intracellular mammalian pathogen Bartonella henselae. Comparative analyses of their circular chromosomes identified 514 orthologous genes shared among all five species. Shared among all five species are 50 identical blocks of microsyntenous orthologous genes (MOGs), containing a total of 283 genes. While retaining highly conserved genomic blocks of microsynteny, divergent evolution, horizontal gene transfer and niche specialization have disrupted macrosynteny among the five circular chromosomes compared. Highly conserved microsyntenous gene clusters help define the Rhizobiales, an order previously defined by 16S RNA gene similarity and herein represented by the three families: Bartonellaceae, Bradyrhizobiaceae and Rhizobiaceae. Genes without orthologs in the other four species help define individual species. The circular chromosomes of each of the five Rhizobiales species examined had genes lacking orthologs in the other four species. For example, 63 proteins are encoded by genes of 'Ca. Liberibacter asiaticus' not shared with other members of the Rhizobiales. Of these 63 proteins, 17 have predicted functions related to DNA replication or RNA transcription, and some of these may have roles related to low genomic GC content. An additional 17 proteins have predicted functions relevant to cellular processes, particularly modifications of the cell surface. Seventeen unshared proteins have specific metabolic functions including a pathway to synthesize cholesterol encoded by a seven-gene operon. The remaining 12 proteins encoded by 'Ca. Liberibacter asiaticus' genes not shared with other Rhizobiales are of bacteriophage origin. 'Ca. Liberibacter asiaticus' shares 11 genes with only Sinorhizobium meliloti and 12 genes are shared with only Bartonella henselae.

Evolution of catabolic pathways and their regulatory systems in synthetic nitroaromatic compounds degrading bacteria

Evolution of catabolic pathways for the degradation of synthetic nitroaromatic compounds is currently ongoing process because these compounds have been in nature only for a short time. Bacteria isolated from contaminated areas contain pathways for the degradation of nitroaromatic compounds at different stages of progression. Therefore, the emergence of pathways for the degradation of such chemicals provides a good opportunity to investigate evolutionary processes leading to the emergence of new metabolic routes and their regulatory systems. In Burkholderia sp. strain DNT the regulatory gene encoding the LysR-type transcriptional regulator DntR is placed divergently of the dinitrotoluene (DNT) dioxygenase genes. This regulator still recognizes salicylate, an effector of its NagR-like ancestor but not DNT. In this issue of Molecular Microbiology, de las Heras et al. demonstrate that the DntR does not respond to any metabolic intermediates of the DNT catabolic pathway. The results of this study suggest that the catabolic pathway for the degradation of DNT has reached to an early stage of evolution when novel specificities of the catabolic enzymes have already acquired but the cognate regulatory system is still missing. This research addresses some fundamental questions about bottlenecks to be solved during evolution of new catabolic operons.

Evolution and multiplicity of arginine decarboxylases in polyamine biosynthesis and essential role in Bacillus subtilis biofilm formation

Arginine decarboxylases (ADCs; EC 4.1.1.19) from four different protein fold families are important for polyamine biosynthesis in bacteria, archaea, and plants. Biosynthetic alanine racemase fold (AR-fold) ADC is widespread in bacteria and plants. We report the discovery and characterization of an ancestral form of the AR-fold ADC in the bacterial Chloroflexi and Bacteroidetes phyla. The ancestral AR-fold ADC lacks a large insertion found in Escherichia coli and plant AR-fold ADC and is more similar to the lysine biosynthetic enzyme meso-diaminopimelate decarboxylase, from which it has evolved. An E. coli acid-inducible ADC belonging to the aspartate aminotransferase fold (AAT-fold) is involved in acid resistance but not polyamine biosynthesis. We report here that the acid-inducible AAT-fold ADC has evolved from a shorter, ancestral biosynthetic AAT-fold ADC by fusion of a response regulator receiver domain protein to the N terminus. Ancestral biosynthetic AAT-fold ADC appears to be limited to firmicute bacteria. The phylogenetic distribution of different forms of ADC distinguishes bacteria from archaea, euryarchaeota from crenarchaeota, double-membraned from single-membraned bacteria, and firmicutes from actinobacteria. Our findings extend to eight the different enzyme forms carrying out the activity described by EC 4.1.1.19. ADC gene clustering reveals that polyamine biosynthesis employs diverse and exchangeable synthetic modules. We show that in Bacillus subtilis, ADC and polyamines are essential for biofilm formation, and this appears to be an ancient, evolutionarily conserved function of polyamines in bacteria. Also of relevance to human health, we found that arginine decarboxylation is the dominant pathway for polyamine biosynthesis in human gut microbiota.