Search

The dominant organisms in modern oxic ecosystems rely on respiratory quinones with high redox potential (HPQs) for electron transport in aerobic respiration and photosynthesis. The diversification of quinones, from low redox potential (LPQ) in anaerobes to HPQs in aerobes, is assumed to have followed Earth's surface oxygenation ~2.3 billion years ago. However, the evolutionary origins of HPQs remain unresolved. Here, we characterize the structure and biosynthetic pathway of an ancestral HPQ, methyl-plastoquinone (mPQ), that is unique to bacteria of the phylum Nitrospirota . mPQ is structurally related to the two previously known HPQs, plastoquinone from Cyanobacteriota /chloroplasts and ubiquinone from Pseudomonadota /mitochondria, respectively. We demonstrate a common origin of the three HPQ biosynthetic pathways that predates the emergence of Nitrospirota , Cyanobacteriota , and Pseudomonadota . An ancestral HPQ biosynthetic pathway evolved ≥ 3.4 billion years ago in an extinct lineage and was laterally transferred to these three phyla ~2.5 to 3.2 billion years ago. We show that Cyanobacteriota and Pseudomonadota were ancestrally aerobic and thus propose that aerobic metabolism using HPQs significantly predates Earth's surface oxygenation. Two of the three HPQ pathways were later obtained by eukaryotes through endosymbiosis forming chloroplasts and mitochondria, enabling their rise to dominance in modern oxic ecosystems., Competing Interests: Competing interests statement:The authors declare no competing interest.

It is currently unclear how Pseudomonadota, a phylum that originated around the time of the Great Oxidation Event, became one of the most abundant and diverse bacterial phyla on Earth, with metabolically versatile members colonizing a wide range of environments with different O2 concentrations. Here, we address this question by studying isoprenoid quinones, which are central components of energy metabolism covering a wide range of redox potentials. We demonstrate that a dynamic repertoire of quinone biosynthetic pathways accompanied the diversification of Pseudomonadota. The low potential menaquinone (MK) was lost in an ancestor of Pseudomonadota while the high potential ubiquinone (UQ) emerged. We show that the O2-dependent and O2-independent UQ pathways were both present in the last common ancestor of Pseudomonadota, and transmitted vertically. The O2-independent pathway has a conserved genetic organization and displays signs of positive regulation by the master regulator "fumarate and nitrate reductase" (FNR), suggesting a conserved role for UQ in anaerobiosis across Pseudomonadota. The O2-independent pathway was lost in some lineages but maintained in others, where it favoured a secondary reacquisition of low potential quinones (MK or rhodoquinone), which promoted diversification towards aerobic facultative and anaerobic metabolisms. Our results support that the ecological success of Pseudomonadota is linked to the acquisition of the largest known repertoire of quinones, which allowed adaptation to oxic niches as O2 levels increased on Earth, and subsequent diversification into anoxic or O2-fluctuating environments., (© The Author(s) 2024. Published by Oxford University Press on behalf of the International Society for Microbial Ecology.)

Ubiquinone (UQ) is a redox polyisoprenoid lipid found in the membranes of bacteria and eukaryotes that has important roles, notably one in respiratory metabolism, which sustains cellular bioenergetics. In Escherichia coli , several steps of the UQ biosynthesis take place in the cytosol. To perform these reactions, a supramolecular assembly called Ubi metabolon is involved. This latter is composed of seven proteins (UbiE, UbiG, UbiF, UbiH, UbiI, UbiJ, and UbiK), and its structural organization is unknown as well as its protein stoichiometry. In this study, a computational framework has been designed to predict the structure of this macromolecular assembly. In several successive steps, we explored the possible protein interactions as well as the protein stoichiometry, to finally obtain a structural organization of the complex. The use of AlphaFold2-based methods combined with evolutionary information enabled us to predict several models whose quality and confidence were further analyzed using different metrics and scores. Our work led to the identification of a "core assembly" that will guide functional and structural characterization of the Ubi metabolon.

Protein secretion systems are complex molecular machineries that translocate proteins through the outer membrane and sometimes through multiple other barriers. They have evolved by co-option of components from other envelope-associated cellular machineries, making them sometimes difficult to identify and discriminate. Here, we describe how to identify protein secretion systems in bacterial genomes using the MacSyFinder program. This flexible computational tool uses the knowledge gathered from experimental studies to identify homologous systems in genome data. It can be used with a set of predefined MacSyFinder models, "TXSScan," to identify all major secretion systems of diderm bacteria (i.e., with inner and LPS-containing outer membranes) as well as evolutionarily related cell appendages (pili and flagella). For this, it identifies and clusters co-localized genes encoding proteins of secretion systems using sequence similarity search with Hidden Markov Model (HMM) protein profiles. Finally, it checks if the clusters' genetic content and genomic organization satisfy the constraints of the model. TXSScan models can be altered in the command line or customized to search for variants of known secretion systems. Models can also be built from scratch to identify novel systems. In this chapter, we describe a complete pipeline of analysis, starting from (i) the integration of information from a reference set of experimentally studied systems, (ii) the identification of conserved proteins and the construction of their HMM protein profiles, (iii) the definition and optimization of "macsy-models," and (iv) their use and online distribution as tools to search genomic data for secretion systems of interest. MacSyFinder is available here: https://github.com/gem-pasteur/macsyfinder, and MacSyFinder models here: https://github.com/macsy-models ., (© 2024. The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature.)

The availability of an ever-increasing diversity of prokaryotic genomes and metagenomes represents a major opportunity to understand and decipher the mechanisms behind the functional diversification of microbial biosynthetic pathways. However, it remains unclear to what extent a pathway producing a specific molecule from a specific precursor can diversify. In this study, we focus on the biosynthesis of ubiquinone (UQ), a crucial coenzyme that is central to the bioenergetics and to the functioning of a wide variety of enzymes in Eukarya and Pseudomonadota (a subgroup of the formerly named Proteobacteria). UQ biosynthesis involves three hydroxylation reactions on contiguous carbon atoms. We and others have previously shown that these reactions are catalyzed by different sets of UQ-hydroxylases that belong either to the iron-dependent Coq7 family or to the more widespread flavin monooxygenase (FMO) family. Here, we combine an experimental approach with comparative genomics and phylogenetics to reveal how UQ-hydroxylases evolved different selectivities within the constrained framework of the UQ pathway. It is shown that the UQ-FMOs diversified via at least three duplication events associated with two cases of neofunctionalization and one case of subfunctionalization, leading to six subfamilies with distinct hydroxylation selectivity. We also demonstrate multiple transfers of the UbiM enzyme and the convergent evolution of UQ-FMOs toward the same function, which resulted in two independent losses of the Coq7 ancestral enzyme. Diversification of this crucial biosynthetic pathway has therefore occurred via a combination of parallel evolution, gene duplications, transfers, and losses., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Plasmodium vivax malaria affects 14 million people each year. Its invasion requires interactions between the parasitic Duffy-binding protein ( Pv DBP) and the N-terminal extracellular domain (ECD1) of the host's Duffy antigen/receptor for chemokines (DARC). ECD1 is highly flexible and intrinsically disordered, therefore it can adopt different conformations. We computationally modeled the challenging ECD1 local structure. With T-REMD simulations, we sampled its dynamic behavior and collected its most representative conformations. Our results suggest that most of the DARC ECD1 domain remains in a disordered state during the simulated time. Globular local conformations are found in the analyzed local free-energy minima. These globular conformations share an α-helix spanning residues Ser18 to Ser29 and in many cases they comprise an antiparallel β-sheet, whose β-strands are formed around residues Leu10 and Ala49. The formation of a parallel β-sheet is almost negligible. So far, progress in understanding the mechanisms forming the basis of the P. vivax malaria infection of reticulocytes has been hampered by experimental difficulties, along with a lack of DARC structural information. Our collection of the most probable ECD1 structural conformations will help to advance modeling of the DARC structure and to explore DARC-ECD1 interactions with a range of physiological and pathological ligands.

Ubiquinone (UQ) is a polyisoprenoid lipid found in the membranes of bacteria and eukaryotes. UQ has important roles, notably in respiratory metabolisms which sustain cellular bioenergetics. Most steps of UQ biosynthesis take place in the cytosol of E. coli within a multiprotein complex called the Ubi metabolon, that contains five enzymes and two accessory proteins, UbiJ and UbiK. The SCP2 domain of UbiJ was proposed to bind the hydrophobic polyisoprenoid tail of UQ biosynthetic intermediates in the Ubi metabolon. How the newly synthesised UQ might be released in the membrane is currently unknown. In this paper, we focused on better understanding the role of the UbiJ-UbiK 2 heterotrimer forming part of the metabolon. Given the difficulties to gain functional insights using biophysical techniques, we applied a multiscale molecular modelling approach to study the UbiJ-UbiK 2 heterotrimer. Our data show that UbiJ-UbiK 2 interacts closely with the membrane and suggests possible pathways to enable the release of UQ into the membrane. This study highlights the UbiJ-UbiK 2 complex as the likely interface between the membrane and the enzymes of the Ubi metabolon and supports that the heterotrimer is key to the biosynthesis of UQ 8 and its release into the membrane of E. coli .

Marine sediments represent a vast habitat for complex microbiomes. Among these, ammonia oxidizing archaea (AOA) of the phylum Thaumarchaeota are one of the most common, yet little explored, inhabitants, which seem extraordinarily well adapted to the harsh conditions of the subsurface biosphere. We present 11 metagenome-assembled genomes of the most abundant AOA clades from sediment cores obtained from the Atlantic Mid-Ocean ridge flanks and Pacific abyssal plains. Their phylogenomic placement reveals three independently evolved clades within the order Nitrosopumilales, of which no cultured representative is known yet. In addition to the gene sets for ammonia oxidation and carbon fixation known from other AOA, all genomes encode an extended capacity for the conversion of fermentation products that can be channeled into the central carbon metabolism, as well as uptake of amino acids probably for protein maintenance or as an ammonia source. Two lineages encode an additional (V-type) ATPase and a large repertoire of DNA repair systems that may allow to overcome the challenges of high hydrostatic pressure. We suggest that the adaptive radiation of AOA into marine sediments occurred more than once in evolution and resulted in three distinct lineages with particular adaptations to this extremely energy-limiting and high-pressure environment., (© 2021. The Author(s).)

No other environment hosts as many microbial cells as the marine sedimentary biosphere. While the majority of these cells are expected to be alive, they are speculated to be persisting in a state of maintenance without net growth due to extreme starvation. Here, we report evidence for in situ growth of anaerobic ammonium-oxidizing (anammox) bacteria in ∼80,000-y-old subsurface sediments from the Arctic Mid-Ocean Ridge. The growth is confined to the nitrate-ammonium transition zone (NATZ), a widespread geochemical transition zone where most of the upward ammonium flux from deep anoxic sediments is being consumed. In this zone the anammox bacteria abundances, assessed by quantification of marker genes, consistently displayed a four order of magnitude increase relative to adjacent layers in four cores. This subsurface cell increase coincides with a markedly higher power supply driven mainly by intensified anammox reaction rates, thereby providing a quantitative link between microbial proliferation and energy availability. The reconstructed draft genome of the dominant anammox bacterium showed an index of replication (iRep) of 1.32, suggesting that 32% of this population was actively replicating. The genome belongs to a Scalindua species which we name Candidatus Scalindua sediminis , so far exclusively found in marine sediments. It has the capacity to utilize urea and cyanate and a mixotrophic lifestyle. Our results demonstrate that specific microbial groups are not only able to survive unfavorable conditions over geological timescales, but can proliferate in situ when encountering ideal conditions with significant consequences for biogeochemical nitrogen cycling., Competing Interests: The authors declare no competing interest., (Copyright © 2020 the Author(s). Published by PNAS.)

Ubiquinone is an important component of the electron transfer chains in proteobacteria and eukaryotes. The biosynthesis of ubiquinone requires multiple steps, most of which are common to bacteria and eukaryotes. Whereas the enzymes of the mitochondrial pathway that produces ubiquinone are highly similar across eukaryotes, recent results point to a rather high diversity of pathways in bacteria. This review focuses on ubiquinone in bacteria, highlighting newly discovered functions and detailing the proteins that are known to participate to its biosynthetic pathways. Novel results showing that ubiquinone can be produced by a pathway independent of dioxygen suggest that ubiquinone may participate to anaerobiosis, in addition to its well-established role for aerobiosis. We also discuss the supramolecular organization of ubiquinone biosynthesis proteins and we summarize the current understanding of the evolution of the ubiquinone pathways relative to those of other isoprenoid quinones like menaquinone and plastoquinone., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2020 Elsevier B.V. All rights reserved.)

Ammonia-oxidizing archaea (AOA) are widespread in nature and are involved in nitrification, an essential process in the global nitrogen cycle. The enzymes for ammonia oxidation and electron transport rely heavily on copper (Cu), which can be limited in nature. In this study the model soil archaeon Nitrososphaera viennensis was investigated via transcriptomic analysis to gain insight regarding possible Cu uptake mechanisms and compensation strategies when Cu becomes limiting. Upon Cu limitation, N. viennensis exhibited impaired nitrite production and thus growth, which was paralleled by downregulation of ammonia oxidation, electron transport, carbon fixation, nucleotide, and lipid biosynthesis pathway genes. Under Cu-limitation, 1547 out of 3180 detected genes were differentially expressed, with 784 genes upregulated and 763 downregulated. The most highly upregulated genes encoded proteins with a possible role in Cu binding and uptake, such as the Cu chelator and transporter CopC/D, disulfide bond oxidoreductase D (dsbD), and multicopper oxidases. While this response differs from the marine strain Nitrosopumilus maritimus, conserved sequence motifs in some of the Cu-responsive genes suggest conserved transcriptional regulation in terrestrial AOA. This study provides possible gene regulation and energy conservation mechanisms linked to Cu bioavailability and presents the first model for Cu uptake by a soil AOA.

Unlike all other archaeal lineages, ammonia-oxidizing archaea (AOA) of the phylum Thaumarchaeota are widespread and abundant in all moderate and oxic environments on Earth. The evolutionary adaptations that led to such unprecedented ecological success of a microbial clade characterized by highly conserved energy and carbon metabolisms have, however, remained underexplored. Here, we reconstructed the genomic content and growth temperature of the ancestor of all AOA, as well as the ancestors of the marine and soil lineages, based on 39 available complete or nearly complete genomes of AOA. Our evolutionary scenario depicts an extremely thermophilic, autotrophic, aerobic ancestor from which three independent lineages of a marine and two terrestrial groups radiated into moderate environments. Their emergence was paralleled by (i) a continuous acquisition of an extensive collection of stress tolerance genes mostly involved in redox maintenance and oxygen detoxification, (ii) an expansion of regulatory capacities in transcription and central metabolic functions, and (iii) an extended repertoire of cell appendages and modifications related to adherence and interactions with the environment. Our analysis provides insights into the evolutionary transitions and key processes that enabled the conquest of the diverse environments in which contemporary AOA are found., (Copyright © 2020 Abby et al.)

Protein secretion is important for many biotic and abiotic interactions. The evolution of protein secretion systems of bacteria, and related nanomachines, occurred by the co-option of machineries for motility, conjugation, injection, or adhesion. Some of these secretion systems emerged many times, whereas others are unique. In most cases, their evolution occurred by successive rounds of gene accretion, deletion, and horizontal transfer, resulting in machines that can be very different from the original ones. The frequency with which such co-option processes occurred seems to depend on the complexity of the systems, their differences to the ancestral machines, the availability of genetic material to tinker with, and possibly on the mechanisms of effector recognition. Understanding the evolution of secretion systems illuminates their functional diversification and could drive the discovery of novel systems., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

We present a computational method to identify conjugative systems in plasmids and chromosomes using the CONJscan module of MacSyFinder. The method relies on the identification of the protein components of the system using hidden Markov model profiles and then checking that the composition and genetic organization of the system is consistent with that expected from a conjugative system. The method can be assessed online using the Galaxy workflow or locally using a standalone software. The latter version allows to modify the models of the module (i.e., to change the expected components, their number, and their organization).CONJscan identifies conjugative systems, but when the mobile genetic element is integrative (ICE), one often also wants to delimit it from the chromosome. We present a method, with a script, to use the results of CONJscan and comparative genomics to delimit ICE in chromosomes. The method provides a visual representation of the ICE location. Together, these methods facilitate the identification of conjugative elements in bacterial genomes.

Processes of molecular innovation require tinkering and shifting in the function of existing genes. How this occurs in terms of molecular evolution at long evolutionary scales remains poorly understood. Here, we analyse the natural history of a vast group of membrane-associated molecular systems in Bacteria and Archaea-the type IV filament (TFF) superfamily-that diversified in systems involved in flagellar or twitching motility, adhesion, protein secretion, and DNA uptake. The phylogeny of the thousands of detected systems suggests they may have been present in the last universal common ancestor. From there, two lineages-a bacterial and an archaeal-diversified by multiple gene duplications, gene fissions and deletions, and accretion of novel components. Surprisingly, we find that the 'tight adherence' (Tad) systems originated from the interkingdom transfer from Archaea to Bacteria of a system resembling the 'EppA-dependent' (Epd) pilus and were associated with the acquisition of a secretin. The phylogeny and content of ancestral systems suggest that initial bacterial pili were engaged in cell motility and/or DNA uptake. In contrast, specialised protein secretion systems arose several times independently and much later in natural history. The functional diversification of the TFF superfamily was accompanied by genetic rearrangements with implications for genetic regulation and horizontal gene transfer: systems encoded in fewer loci were more frequently exchanged between taxa. This may have contributed to their rapid evolution and spread across Bacteria and Archaea. Hence, the evolutionary history of the superfamily reveals an impressive catalogue of molecular evolution mechanisms that resulted in remarkable functional innovation and specialisation from a relatively small set of components., Competing Interests: The authors have declared that no competing interests exist.

Climate change is causing arctic regions to warm disproportionally faster than those at lower latitudes, leading to alterations in carbon and nitrogen cycling, and potentially higher greenhouse gas emissions. It is thus increasingly important to better characterize the microorganisms driving arctic biogeochemical processes and their potential responses to changing conditions. Here, we describe a novel thaumarchaeon enriched from an arctic soil, Candidatus Nitrosocosmicus arcticus strain Kfb, which has been maintained for seven years in stable laboratory enrichment cultures as an aerobic ammonia oxidizer, with ammonium or urea as substrates. Genomic analyses show that this organism harbors all genes involved in ammonia oxidation and in carbon fixation via the 3-hydroxypropionate/4-hydroxybutyrate cycle, characteristic of all AOA, as well as the capability for urea utilization and potentially also for heterotrophic metabolism, similar to other AOA. Ca . N. arcticus oxidizes ammonia optimally between 20 and 28°C, well above average temperatures in its native high arctic environment (-13-4°C). Ammonia oxidation rates were nevertheless much lower than those of most cultivated mesophilic AOA (20-45°C). Intriguingly, we repeatedly observed apparent faster growth rates (based on marker gene counts) at lower temperatures (4-8°C) but without detectable nitrite production. Together with potential metabolisms predicted from its genome content, these observations indicate that Ca . N. arcticus is not a strict chemolithotrophic ammonia oxidizer and add to cumulating evidence for a greater metabolic and physiological versatility of AOA. The physiology of Ca . N. arcticus suggests that increasing temperatures might drastically affect nitrification in arctic soils by stimulating archaeal ammonia oxidation.

Most bacteria can generate ATP by respiratory metabolism, in which electrons are shuttled from reduced substrates to terminal electron acceptors, via quinone molecules like ubiquinone. Dioxygen (O 2 ) is the terminal electron acceptor of aerobic respiration and serves as a co-substrate in the biosynthesis of ubiquinone. Here, we characterize a novel, O 2 -independent pathway for the biosynthesis of ubiquinone. This pathway relies on three proteins, UbiT (YhbT), UbiU (YhbU), and UbiV (YhbV). UbiT contains an SCP2 lipid-binding domain and is likely an accessory factor of the biosynthetic pathway, while UbiU and UbiV (UbiU-UbiV) are involved in hydroxylation reactions and represent a novel class of O 2 -independent hydroxylases. We demonstrate that UbiU-UbiV form a heterodimer, wherein each protein binds a 4Fe-4S cluster via conserved cysteines that are essential for activity. The UbiT, -U, and -V proteins are found in alpha-, beta-, and gammaproteobacterial clades, including several human pathogens, supporting the widespread distribution of a previously unrecognized capacity to synthesize ubiquinone in the absence of O 2 Together, the O 2 -dependent and O 2 -independent ubiquinone biosynthesis pathways contribute to optimizing bacterial metabolism over the entire O 2 range. IMPORTANCE In order to colonize environments with large O 2 gradients or fluctuating O 2 levels, bacteria have developed metabolic responses that remain incompletely understood. Such adaptations have been recently linked to antibiotic resistance, virulence, and the capacity to develop in complex ecosystems like the microbiota. Here, we identify a novel pathway for the biosynthesis of ubiquinone, a molecule with a key role in cellular bioenergetics. We link three uncharacterized genes of Escherichia coli to this pathway and show that the pathway functions independently from O 2 In contrast, the long-described pathway for ubiquinone biosynthesis requires O 2 as a substrate. In fact, we find that many proteobacteria are equipped with the O 2 -dependent and O 2 -independent pathways, supporting that they are able to synthesize ubiquinone over the entire O 2 range. Overall, we propose that the novel O 2 -independent pathway is part of the metabolic plasticity developed by proteobacteria to face various environmental O 2 levels., (Copyright © 2019 Pelosi et al.)

Diazotrophic Actinobacteria of the genus Frankia represent a challenge to classical bacterial taxonomy as they include many unculturable strains. As a consequence, we still have a poor understanding of their diversity, evolution and biogeography. In this study, a Multi-Locus Sequence Analysis (MLSA) using atpD, dnaA, ftsZ, pgk, and rpoB loci was done on a large set of cultured and uncultured strains, compared to 16S rRNA and correlated to Average Nucleotide Identity (ANI) from available Frankia genomes. MLSA provided a robust resolution of Frankia genus phylogeny and clarified the status of unresolved species and complex of species. The robustness of single-gene topologies and their congruence with the MLSA tree were tested. Lateral Gene Transfers (LGT) were few and scattered, suggesting they had no impact on the concatenate topology. The pgk marker - providing the longest sequence, highest mean genetic divergence and least occurrence of LGT - was used to survey an unequalled number of Alnus-infective Frankia - mainly uncultured strains from a broad range of host-species and geographic origins. This marker allowed reliable Single-Locus Strain Typing (SLST) below the species level, revealed an undiscovered taxonomical diversity, and highlighted the effect of cultivation, sporulation phenotype and host plant species on symbiont richness, diversity and phylogeny., (Copyright © 2018 Elsevier GmbH. All rights reserved.)

Ammonia oxidizing archaea (AOA) of the phylum Thaumarchaeota are widespread in moderate environments but their occurrence and activity has also been demonstrated in hot springs. Here we present the first enrichment of a thermophilic representative with a sequenced genome, which facilitates the search for adaptive strategies and for traits that shape the evolution of Thaumarchaeota. Candidatus Nitrosocaldus cavascurensis has been enriched from a hot spring in Ischia, Italy. It grows optimally at 68°C under chemolithoautotrophic conditions on ammonia or urea converting ammonia stoichiometrically into nitrite with a generation time of approximately 23 h. Phylogenetic analyses based on ribosomal proteins place the organism as a sister group to all known mesophilic AOA. The 1.58 Mb genome of Ca. N. cavascurensis harbors an amo AXCB gene cluster encoding ammonia monooxygenase and genes for a 3-hydroxypropionate/4-hydroxybutyrate pathway for autotrophic carbon fixation, but also genes that indicate potential alternative energy metabolisms. Although a bona fide gene for nitrite reductase is missing, the organism is sensitive to NO-scavenging, underlining the potential importance of this compound for AOA metabolism. Ca. N. cavascurensis is distinct from all other AOA in its gene repertoire for replication, cell division and repair. Its genome has an impressive array of mobile genetic elements and other recently acquired gene sets, including conjugative systems, a provirus, transposons and cell appendages. Some of these elements indicate recent exchange with the environment, whereas others seem to have been domesticated and might convey crucial metabolic traits.

Protein secretion systems are complex molecular machineries that translocate proteins through the outer membrane, and sometimes through multiple other barriers. They have evolved by co-option of components from other envelope-associated cellular machineries, making them sometimes difficult to identify and discriminate. Here, we describe how to identify protein secretion systems in bacterial genomes using MacSyFinder. This flexible computational tool uses the knowledge stemming from experimental studies to identify homologous systems in genome data. It can be used with a set of predefined models-"TXSScan"-to identify all major secretion systems of diderm bacteria (i.e., with inner and with LPS-containing outer membranes). For this, it identifies and clusters colocalized components of secretion systems using sequence similarity searches with hidden Markov model protein profiles. Finally, it checks whether the genetic content and organization of clusters satisfy the constraints of the model. TXSScan models can be customized to search for variants of known systems. The models can also be built from scratch to identify novel systems. In this chapter, we describe a complete pipeline of analysis, including the identification of a reference set of experimentally studied systems, the identification of components and the construction of their protein profiles, the definition of the models, their optimization, and, finally, their use as tools to search genomic data.

The reproduction mode of uncultivable microorganisms deserves investigation as it can largely diverge from conventional transverse binary fission. Here, we show that the rod-shaped gammaproteobacterium thriving on the surface of the Robbea hypermnestra nematode divides by FtsZ-based, non-synchronous invagination of its poles-that is, the host-attached and fimbriae-rich pole invaginates earlier than the distal one. We conclude that, in a naturally occurring animal symbiont, binary fission is host-oriented and does not require native FtsZ to polymerize into a ring at any septation stage.

Evolutionary events co-occurring along phylogenetic trees usually point to complex adaptive phenomena, sometimes implicating epistasis. While a number of methods have been developed to account for co-occurrence of events on the same internal or external branch of an evolutionary tree, there is a need to account for the larger diversity of possible relative positions of events in a tree. Here we propose a method to quantify to what extent two or more evolutionary events are associated on a phylogenetic tree. The method is applicable to any discrete character, like substitutions within a coding sequence or gains/losses of a biological function. Our method uses a general approach to statistically test for significant associations between events along the tree, which encompasses both events inseparable on the same branch, and events genealogically ordered on different branches. It assumes that the phylogeny and themapping of branches is known without errors. We address this problem from the statistical viewpoint by a linear algebra representation of the localization of the evolutionary events on the tree.We compute the full probability distribution of the number of paired events occurring in the same branch or in different branches of the tree, under a null model of independence where each type of event occurs at a constant rate uniformly inthephylogenetic tree. The strengths andweaknesses of themethodare assessed via simulations;we then apply the method to explore the loss of cell motility in intracellular pathogens., (© The Author(s) 2016. Published by Oxford University Press, on behalf of the Society of Systematic Biologists. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Bacteria with two cell membranes (diderms) have evolved complex systems for protein secretion. These systems were extensively studied in some model bacteria, but the characterisation of their diversity has lagged behind due to lack of standard annotation tools. We built online and standalone computational tools to accurately predict protein secretion systems and related appendages in bacteria with LPS-containing outer membranes. They consist of models describing the systems' components and genetic organization to be used with MacSyFinder to search for T1SS-T6SS, T9SS, flagella, Type IV pili and Tad pili. We identified ~10,000 candidate systems in bacterial genomes, where T1SS and T5SS were by far the most abundant and widespread. All these data are made available in a public database. The recently described T6SS(iii) and T9SS were restricted to Bacteroidetes, and T6SS(ii) to Francisella. The T2SS, T3SS, and T4SS were frequently encoded in single-copy in one locus, whereas most T1SS were encoded in two loci. The secretion systems of diderm Firmicutes were similar to those found in other diderms. Novel systems may remain to be discovered, since some clades of environmental bacteria lacked all known protein secretion systems. Our models can be fully customized, which should facilitate the identification of novel systems.

Marine phytoplankton produce half of the oxygen we breathe and their astounding diversity is just starting to be unraveled. Many microbial phytoplankton are thought to be phototrophic, depending solely on inorganic sources of carbon and minerals for growth rather than preying on other planktonic cells. However, there is increasing evidence that symbiotic associations, to a large extent with bacteria, are required for vitamin or nutrient uptake for many eukaryotic microalgae. Here, we use in silico approaches to look for putative symbiotic interactions by analysing the gene content of microbial communities associated with 13 different Ostreococcus tauri (Chlorophyta, Mamilleophyceae) cultures sampled from the Mediterranean Sea. While we find evidence for bacteria in all cultures, there is no ubiquitous bacterial group, and the most prevalent group, Flavobacteria, is present in 10 out of 13 cultures. Among seven of the microbiomes, we detected genes predicted to encode type 3 secretion systems (T3SS, in 6/7 microbiomes) and/or putative type 6 secretion systems (T6SS, in 4/7 microbiomes). Phylogenetic analyses show that the corresponding genes are closely related to genes of systems identified in bacterial-plant interactions, suggesting that these T3SS might be involved in cell-to-cell interactions with O. tauri.

Motivation: Biologists often wish to use their knowledge on a few experimental models of a given molecular system to identify homologs in genomic data. We developed a generic tool for this purpose., Results: Macromolecular System Finder (MacSyFinder) provides a flexible framework to model the properties of molecular systems (cellular machinery or pathway) including their components, evolutionary associations with other systems and genetic architecture. Modelled features also include functional analogs, and the multiple uses of a same component by different systems. Models are used to search for molecular systems in complete genomes or in unstructured data like metagenomes. The components of the systems are searched by sequence similarity using Hidden Markov model (HMM) protein profiles. The assignment of hits to a given system is decided based on compliance with the content and organization of the system model. A graphical interface, MacSyView, facilitates the analysis of the results by showing overviews of component content and genomic context. To exemplify the use of MacSyFinder we built models to detect and class CRISPR-Cas systems following a previously established classification. We show that MacSyFinder allows to easily define an accurate "Cas-finder" using publicly available protein profiles., Availability and Implementation: MacSyFinder is a standalone application implemented in Python. It requires Python 2.7, Hmmer and makeblastdb (version 2.2.28 or higher). It is freely available with its source code under a GPLv3 license at https://github.com/gem-pasteur/macsyfinder. It is compatible with all platforms supporting Python and Hmmer/makeblastdb. The "Cas-finder" (models and HMM profiles) is distributed as a compressed tarball archive as Supporting Information.

Conjugation of DNA through a type IV secretion system (T4SS) drives horizontal gene transfer. Yet little is known on the diversity of these nanomachines. We previously found that T4SS can be divided in eight classes based on the phylogeny of the only ubiquitous protein of T4SS (VirB4). Here, we use an ab initio approach to identify protein families systematically and specifically associated with VirB4 in each class. We built profiles for these proteins and used them to scan 2262 genomes for the presence of T4SS. Our analysis led to the identification of thousands of occurrences of 116 protein families for a total of 1623 T4SS. Importantly, we could identify almost always in our profiles the essential genes of well-studied T4SS. This allowed us to build a database with the largest number of T4SS described to date. Using profile-profile alignments, we reveal many new cases of homology between components of distant classes of T4SS. We mapped these similarities on the T4SS phylogenetic tree and thus obtained the patterns of acquisition and loss of these protein families in the history of T4SS. The identification of the key VirB4-associated proteins paves the way toward experimental analysis of poorly characterized T4SS classes., (© The Author(s) 2014. Published by Oxford University Press.)

The timing of the evolution of microbial life has largely remained elusive due to the scarcity of prokaryotic fossil record and the confounding effects of the exchange of genes among possibly distant species. The history of gene transfer events, however, is not a series of individual oddities; it records which lineages were concurrent and thus provides information on the timing of species diversification. Here, we use a probabilistic model of genome evolution that accounts for differences between gene phylogenies and the species tree as series of duplication, transfer, and loss events to reconstruct chronologically ordered species phylogenies. Using simulations we show that we can robustly recover accurate chronologically ordered species phylogenies in the presence of gene tree reconstruction errors and realistic rates of duplication, transfer, and loss. Using genomic data we demonstrate that we can infer rooted species phylogenies using homologous gene families from complete genomes of 10 bacterial and archaeal groups. Focusing on cyanobacteria, distinguished among prokaryotes by a relative abundance of fossils, we infer the maximum likelihood chronologically ordered species phylogeny based on 36 genomes with 8,332 homologous gene families. We find the order of speciation events to be in full agreement with the fossil record and the inferred phylogeny of cyanobacteria to be consistent with the phylogeny recovered from established phylogenomics methods. Our results demonstrate that lateral gene transfers, detected by probabilistic models of genome evolution, can be used as a source of information on the timing of evolution, providing a valuable complement to the limited prokaryotic fossil record.

Type 3 secretion systems (T3SSs) are essential components of two complex bacterial machineries: the flagellum, which drives cell motility, and the non-flagellar T3SS (NF-T3SS), which delivers effectors into eukaryotic cells. Yet the origin, specialization, and diversification of these machineries remained unclear. We developed computational tools to identify homologous components of the two systems and to discriminate between them. Our analysis of >1,000 genomes identified 921 T3SSs, including 222 NF-T3SSs. Phylogenomic and comparative analyses of these systems argue that the NF-T3SS arose from an exaptation of the flagellum, i.e. the recruitment of part of the flagellum structure for the evolution of the new protein delivery function. This reconstructed chronology of the exaptation process proceeded in at least two steps. An intermediate ancestral form of NF-T3SS, whose descendants still exist in Myxococcales, lacked elements that are essential for motility and included a subset of NF-T3SS features. We argue that this ancestral version was involved in protein translocation. A second major step in the evolution of NF-T3SSs occurred via recruitment of secretins to the NF-T3SS, an event that occurred at least three times from different systems. In rhizobiales, a partial homologous gene replacement of the secretin resulted in two genes of complementary function. Acquisition of a secretin was followed by the rapid adaptation of the resulting NF-T3SSs to multiple, distinct eukaryotic cell envelopes where they became key in parasitic and mutualistic associations between prokaryotes and eukaryotes. Our work elucidates major steps of the evolutionary scenario leading to extant NF-T3SSs. It demonstrates how molecular evolution can convert one complex molecular machine into a second, equally complex machine by successive deletions, innovations, and recruitment from other molecular systems., Competing Interests: The authors have declared that no competing interests exist.

Secretins form large multimeric complexes in the outer membranes of many Gram-negative bacteria, where they function as dedicated gateways that allow proteins to access the extracellular environment. Despite their overall relatedness, different secretins use different specific and general mechanisms for their targeting, assembly, and membrane insertion. We report that all tested secretins from several type II secretion systems and from the filamentous bacteriophage f1 can spontaneously multimerize and insert into liposomes in an in vitro transcription-translation system. Phylogenetic analyses indicate that these secretins form a group distinct from the secretins of the type IV piliation and type III secretion systems, which do not autoassemble in vitro. A mutation causing a proline-to-leucine substitution allowed PilQ secretins from two different type IV piliation systems to assemble in vitro, albeit with very low efficiency, suggesting that autoassembly is an inherent property of all secretins.

Lateral gene transfer (LGT), the acquisition of genes from other species, is a major evolutionary force. However, its success as an adaptive process makes the reconstruction of the history of life an intricate puzzle: If no gene has remained unaffected during the course of life's evolution, how can one rely on molecular markers to reconstruct the relationships among species? Here, we take a completely different look at LGT and its impact for the reconstruction of the history of life. Rather than trying to remove the effect of LGT in phylogenies, and ignoring as a result most of the information of gene histories, we use an explicit phylogenetic model of gene transfer to reconcile gene histories with the tree of species. We studied 16 bacterial and archaeal phyla, representing a dataset of 12,000 gene families distributed in 336 genomes. Our results show that, in most phyla, LGT provides an abundant phylogenetic signal on the pattern of species diversification and that this signal is robust to the choice of gene families under study. We also find that LGT brings an abundant signal on the location of the root of species trees, which has been previously overlooked. Our results quantify the great variety of gene transfer rates among lineages of the tree of life and provide strong support for the "complexity hypothesis," which states that genes whose products participate to macromolecular protein complexes are relatively resistant to transfer.

Many studies have been devoted to understand the mechanisms used by pathogenic bacteria to exploit human hosts. These mechanisms are very diverse in the detail, but share commonalities whose quantification should enlighten the evolution of virulence from both a molecular and an ecological perspective. We mined the literature for experimental data on infectious dose of bacterial pathogens in humans (ID50) and also for traits with which ID50 might be associated. These compilations were checked and complemented with genome analyses. We observed that ID50 varies in a continuous way by over 10 orders of magnitude. Low ID50 values are very strongly associated with the capacity of the bacteria to kill professional phagocytes or to survive in the intracellular milieu of these cells. Inversely, high ID50 values are associated with motile and fast-growing bacteria that use quorum-sensing based regulation of virulence factors expression. Infectious dose is not associated with genome size and shows insignificant phylogenetic inertia, in line with frequent virulence shifts associated with the horizontal gene transfer of a small number of virulence factors. Contrary to previous proposals, infectious dose shows little dependence on contact-dependent secretion systems and on the natural route of exposure. When all variables are combined, immune subversion and quorum-sensing are sufficient to explain two thirds of the variance in infectious dose. Our results show the key role of immune subversion in effective human infection by small bacterial populations. They also suggest that cooperative processes might be important for successful infection by bacteria with high ID50. Our results suggest that trade-offs between selection for population growth-related traits and selection for the ability to subvert the immune system shape bacterial infectiousness. Understanding these trade-offs provides guidelines to study the evolution of virulence and in particular the micro-evolutionary paths of emerging pathogens.

Proteins evolve at very different rates and, most notably, at rates inversely proportional to the level at which they are produced. The relative frequency of highly expressed proteins in the proteome, and thus their impact on the cell budget, increases steeply with growth rate. The maximal growth rate is a key life-history trait reflecting trade-offs between rapid growth and other fitness components. We show that the maximal growth rate is weakly affected by genetic drift. The negative correlation between protein expression levels and evolutionary rate and the positive correlation between expression levels of highly expressed proteins and growth rates, suggest that investment in growth affects the evolutionary rate of proteins, especially the highly expressed ones. Accordingly, analysis of 61 families of orthologs in 74 proteobacteria shows that differences in evolutionary rates between lowly and highly expressed proteins depend on maximal growth rates. Analyses of complexes with key roles in bacterial growth and strikingly different expression levels, the ribosome and the replisome, confirm these patterns and suggest that the growth-related sequence conservation is associated with protein synthesis. Maximal growth rates also shape protein evolution in the other bacterial clades. Long-branch attractions associated with this effect might explain why clades with persistent history of slow growth are attracted to the root when the tree of prokaryotes is inferred using highly, but not lowly, expressed proteins. These results indicate that reconstruction of deep phylogenies can be strongly affected by maximal growth rates, and highlight the importance of life-history traits and their physiological consequences for protein evolution.

Background: To understand the evolutionary role of Lateral Gene Transfer (LGT), accurate methods are needed to identify transferred genes and infer their timing of acquisition. Phylogenetic methods are particularly promising for this purpose, but the reconciliation of a gene tree with a reference (species) tree is computationally hard. In addition, the application of these methods to real data raises the problem of sorting out real and artifactual phylogenetic conflict., Results: We present Prunier, a new method for phylogenetic detection of LGT based on the search for a maximum statistical agreement forest (MSAF) between a gene tree and a reference tree. The program is flexible as it can use any definition of "agreement" among trees. We evaluate the performance of Prunier and two other programs (EEEP and RIATA-HGT) for their ability to detect transferred genes in realistic simulations where gene trees are reconstructed from sequences. Prunier proposes a single scenario that compares to the other methods in terms of sensitivity, but shows higher specificity. We show that LGT scenarios carry a strong signal about the position of the root of the species tree and could be used to identify the direction of evolutionary time on the species tree. We use Prunier on a biological dataset of 23 universal proteins and discuss their suitability for inferring the tree of life., Conclusions: The ability of Prunier to take into account branch support in the process of reconciliation allows a gain in complexity, in comparison to EEEP, and in accuracy in comparison to RIATA-HGT. Prunier's greedy algorithm proposes a single scenario of LGT for a gene family, but its quality always compares to the best solutions provided by the other algorithms. When the root position is uncertain in the species tree, Prunier is able to infer a scenario per root at a limited additional computational cost and can easily run on large datasets.Prunier is implemented in C++, using the Bio++ library and the phylogeny program Treefinder. It is available at: http://pbil.univ-lyon1.fr/software/prunier.