MacSyFinder: a program to mine genomes for molecular systems with an application to CRISPR-Cas systems

Repeated horizontal gene transfers triggered parallel evolution of magnetotaxis in two evolutionary divergent lineages of magnetotactic bacteria

Under the same selection pressures, two genetically divergent populations may evolve in parallel toward the same adaptive solutions. Here, we hypothesized that magnetotaxis (i.e., magnetically guided chemotaxis) represents a key adaptation to micro-oxic habitats in aquatic sediments and that its parallel evolution homogenized the phenotypes of two evolutionary divergent clusters of freshwater spirilla. All magnetotactic bacteria affiliated to the Magnetospirillum genus (Alphaproteobacteria class) biomineralize the same magnetic particle chains and share highly similar physiological and ultrastructural features. We looked for the processes that could have contributed at shaping such an evolutionary pattern by reconciling species and gene trees using newly sequenced genomes of Magnetospirillum related bacteria. We showed that repeated horizontal gene transfers and homologous recombination of entire operons contributed to the parallel evolution of magnetotaxis. We propose that such processes could represent a more parsimonious and rapid solution for adaptation compared with independent and repeated de novo mutations, especially in the case of traits as complex as magnetotaxis involving tens of interacting proteins. Besides strengthening the idea about the importance of such a function in micro-oxic habitats, these results reinforce previous observations in experimental evolution suggesting that gene flow could alleviate clonal interference and speed up adaptation under some circumstances.

Phylogenetic background and habitat drive the genetic diversification of Escherichia coli

Escherichia coli is mostly a commensal of birds and mammals, including humans, where it can act as an opportunistic pathogen. It is also found in water and sediments. We investigated the phylogeny, genetic diversification, and habitat-association of 1,294 isolates representative of the phylogenetic diversity of more than 5,000 isolates from the Australian continent. Since many previous studies focused on clinical isolates, we investigated mostly other isolates originating from humans, poultry, wild animals and water. These strains represent the species genetic diversity and reveal widespread associations between phylogroups and isolation sources. The analysis of strains from the same sequence types revealed very rapid change of gene repertoires in the very early stages of divergence, driven by the acquisition of many different types of mobile genetic elements. These elements also lead to rapid variations in genome size, even if few of their genes rise to high frequency in the species. Variations in genome size are associated with phylogroup and isolation sources, but the latter determine the number of MGEs, a marker of recent transfer, suggesting that gene flow reinforces the association of certain genetic backgrounds with specific habitats. After a while, the divergence of gene repertoires becomes linear with phylogenetic distance, presumably reflecting the continuous turnover of mobile element and the occasional acquisition of adaptive genes. Surprisingly, the phylogroups with smallest genomes have the highest rates of gene repertoire diversification and fewer but more diverse mobile genetic elements. This suggests that smaller genomes are associated with higher, not lower, turnover of genetic information. Many of these genomes are from freshwater isolates and have peculiar traits, including a specific capsule, suggesting adaptation to this environment. Altogether, these data contribute to explain why epidemiological clones tend to emerge from specific phylogenetic groups in the presence of pervasive horizontal gene transfer across the species.

Previous large scale studies on the evolution of E. coli focused on clinical isolates emphasizing virulence and antibiotic resistance in medically important lineages. Yet, most E. coli strains are either human commensals or not associated with humans at all. Here, we analyzed a large collection of non-clinical isolates of the species to assess the mechanisms of gene repertoire diversification in the light of isolation sources and phylogeny. We show that gene repertoires evolve so rapidly by the high turnover of mobile genetic elements that epidemiologically indistinguishable strains can be phenotypically extremely heterogeneous, illustrating the velocity of bacterial adaptation and the importance of accounting for the information on the whole genome at the epidemiological scale. Phylogeny and habitat shape the genetic diversification of E. coli to similar extents. Surprisingly, freshwater strains seem specifically adapted to this environment, breaking the paradigm that E. coli environmental isolates are systematically fecal contaminations. As a consequence, the evolution of this species is also shaped by environmental habitats, and it may diversify by acquiring genes and mobile elements from environmental bacteria (and not just from gut bacteria). This may facilitate the acquisition of virulence factors and antibiotic resistance in the strains that become pathogenic.

CRISPRcasIdentifier: Machine learning for accurate identification and classification of CRISPR-Cas systems

BACKGROUND

CRISPR-Cas genes are extraordinarily diverse and evolve rapidly when compared to other prokaryotic genes. With the rapid increase in newly sequenced archaeal and bacterial genomes, manual identification of CRISPR-Cas systems is no longer viable. Thus, an automated approach is required for advancing our understanding of the evolution and diversity of these systems and for finding new candidates for genome engineering in eukaryotic models.

RESULTS

We introduce CRISPRcasIdentifier, a new machine learning-based tool that combines regression and classification models for the prediction of potentially missing proteins in instances of CRISPR-Cas systems and the prediction of their respective subtypes. In contrast to other available tools, CRISPRcasIdentifier can both detect cas genes and extract potential association rules that reveal functional modules for CRISPR-Cas systems. In our experimental benchmark on the most recently published and comprehensive CRISPR-Cas system dataset, CRISPRcasIdentifier was compared with recent and state-of-the-art tools. According to the experimental results, CRISPRcasIdentifier presented the best Cas protein identification and subtype classification performance.

CONCLUSIONS

Overall, our tool greatly extends the classification of CRISPR cassettes and, for the first time, predicts missing Cas proteins and association rules between Cas proteins. Additionally, we investigated the properties of CRISPR subtypes. The proposed tool relies not only on the knowledge of manual CRISPR annotation but also on models trained using machine learning.

Virulence Potential of a Multidrug-Resistant Escherichia coli Strain Belonging to the Emerging Clonal Group ST101-B1 Isolated from Bloodstream Infection

Escherichia coli EC121 is a multidrug-resistant (MDR) strain isolated from a bloodstream infection of an inpatient with persistent gastroenteritis and T-zone lymphoma that died due to septic shock. Despite causing an extraintestinal infection, previous studies showed that it did not have the usual characteristics of an extraintestinal pathogenic E. coli. Instead, it belonged to phylogenetic group B1 and harbored few known virulence genes. To evaluate the pathogenic potential of strain EC121, an extensive genome sequencing and in vitro characterization of various pathogenicity-associated properties were performed. The genomic analysis showed that strain EC121 harbors more than 50 complete virulence genetic clusters. It also displays the capacity to adhere to a variety of epithelial cell lineages and invade T24 bladder cells, as well as the ability to form biofilms on abiotic surfaces, and survive the bactericidal serum complement activity. Additionally, EC121 was shown to be virulent in the Galleria mellonella model. Furthermore, EC121 is an MDR strain harboring 14 antimicrobial resistance genes, including bla_CTX-M-2. Completing the scenario, it belongs to serotype O154:H25 and to sequence type 101-B1, which has been epidemiologically linked to extraintestinal infections as well as to antimicrobial resistance spread. This study with E. coli strain EC121 shows that clinical isolates considered opportunistic might be true pathogens that go underestimated.

CRISPRCasdb a successor of CRISPRdb containing CRISPR arrays and cas genes from complete genome sequences, and tools to download and query lists of repeats and spacers

In Archaea and Bacteria, the arrays called CRISPRs for 'clustered regularly interspaced short palindromic repeats' and the CRISPR associated genes or cas provide adaptive immunity against viruses, plasmids and transposable elements. Short sequences called spacers, corresponding to fragments of invading DNA, are stored in-between repeated sequences. The CRISPR-Cas systems target sequences homologous to spacers leading to their degradation. To facilitate investigations of CRISPRs, we developed 12 years ago a website holding the CRISPRdb. We now propose CRISPRCasdb, a completely new version giving access to both CRISPRs and cas genes. We used CRISPRCasFinder, a program that identifies CRISPR arrays and cas genes and determine the system's type and subtype, to process public whole genome assemblies. Strains are displayed either in an alphabetic list or in taxonomic order. The database is part of the CRISPR-Cas++ website which also offers the possibility to analyse submitted sequences and to download programs. A BLAST search against lists of repeats and spacers extracted from the database is proposed. To date, 16 990 complete prokaryote genomes (16 650 bacteria from 2973 species and 340 archaea from 300 species) are included. CRISPR-Cas systems were found in 36% of Bacteria and 75% of Archaea strains. CRISPRCasdb is freely accessible at https://crisprcas.i2bc.paris-saclay.fr/.

MicroScope: an integrated platform for the annotation and exploration of microbial gene functions through genomic, pangenomic and metabolic comparative analysis

Large-scale genome sequencing and the increasingly massive use of high-throughput approaches produce a vast amount of new information that completely transforms our understanding of thousands of microbial species. However, despite the development of powerful bioinformatics approaches, full interpretation of the content of these genomes remains a difficult task. Launched in 2005, the MicroScope platform (https://www.genoscope.cns.fr/agc/microscope) has been under continuous development and provides analysis for prokaryotic genome projects together with metabolic network reconstruction and post-genomic experiments allowing users to improve the understanding of gene functions. Here we present new improvements of the MicroScope user interface for genome selection, navigation and expert gene annotation. Automatic functional annotation procedures of the platform have also been updated and we added several new tools for the functional annotation of genes and genomic regions. We finally focus on new tools and pipeline developed to perform comparative analyses on hundreds of genomes based on pangenome graphs. To date, MicroScope contains data for >11 800 microbial genomes, part of which are manually curated and maintained by microbiologists (>4500 personal accounts in September 2019). The platform enables collaborative work in a rich comparative genomic context and improves community-based curation efforts.

The major subunit of widespread competence pili exhibits a novel and conserved type IV pilin fold

Type IV filaments (T4F), which are helical assemblies of type IV pilins, constitute a superfamily of filamentous nanomachines virtually ubiquitous in prokaryotes that mediate a wide variety of functions. The competence (Com) pilus is a widespread T4F, mediating DNA uptake (the first step in natural transformation) in bacteria with one membrane (monoderms), an important mechanism of horizontal gene transfer. Here, we report the results of genomic, phylogenetic, and structural analyses of ComGC, the major pilin subunit of Com pili. By performing a global comparative analysis, we show that Com pili genes are virtually ubiquitous in Bacilli, a major monoderm class of Firmicutes. This also revealed that ComGC displays extensive sequence conservation, defining a monophyletic group among type IV pilins. We further report ComGC solution structures from two naturally competent human pathogens, Streptococcus sanguinis (ComGCSS) and Streptococcus pneumoniae (ComGCSP), revealing that this pilin displays extensive structural conservation. Strikingly, ComGCSS and ComGCSP exhibit a novel type IV pilin fold that is purely helical. Results from homology modeling analyses suggest that the unusual structure of ComGC is compatible with helical filament assembly. Because ComGC displays such a widespread distribution, these results have implications for hundreds of monoderm species.

Identification of integrative and conjugative elements in pathogenic and commensal Neisseriaceae species via genomic distributions of DNA uptake sequence dialects

Mobile genetic elements (MGEs) are key factors responsible for dissemination of virulence determinants and antimicrobial-resistance genes amongst pathogenic bacteria. Conjugative MGEs are notable for their high gene loads donated per transfer event, broad host ranges and phylogenetic ubiquity amongst prokaryotes, with the subclass of chromosomally inserted integrative and conjugative elements (ICEs) being particularly abundant. The focus on a small number of model systems has biased the study of ICEs towards those conferring readily selectable phenotypes to host cells, whereas the identification and characterization of integrated cryptic elements remains challenging. Even though antimicrobial resistance and horizontally acquired virulence genes are major factors aggravating neisserial infection, conjugative MGEs of Neisseria gonorrhoeae and Neisseria meningitidis remain poorly characterized. Using a phenotype-independent approach based on atypical distributions of DNA uptake sequences (DUSs) in MGEs relative to the chromosomal background, we have identified two groups of chromosomally integrated conjugative elements in Neisseria: one found almost exclusively in pathogenic species possibly deriving from the genus Kingella, the other belonging to a group of Neisseria mucosa-like commensals. The former element appears to enable transfer of traditionally gonococcal-specific loci such as the virulence-associated toxin-antitoxin system fitAB to N. meningitidis chromosomes, whilst the circular form of the latter possesses a unique attachment site (attP) sequence seemingly adapted to exploit DUS motifs as chromosomal integration sites. In addition to validating the use of DUS distributions in Neisseriaceae MGE identification, the >170 identified ICE sequences provide a valuable resource for future studies of ICE evolution and host adaptation.

Comparative Genomic Analysis Provides Insights into the Phylogeny, Resistome, Virulome, and Host Adaptation in the Genus Ewingella

Ewingella americana is a cosmopolitan bacterial pathogen that has been isolated from many hosts. Here, we sequenced a high-quality genome of E. americana B6-1 isolated from Flammulina filiformis, an important cultivated mushroom, performed a comparative genomic analysis with four other E. americana strains from various origins, and tested the susceptibility of B6-1 to antibiotics. The genome size, predicted genes, and GC (guanine-cytosine) content of B6-1 was 4.67 Mb, 4301, and 53.80%, respectively. The origin of the strains did not significantly affect the phylogeny, but mobile genetic elements shaped the evolution of the genus Ewingella. The strains encoded a set of common genes for type secretion, virulence effectors, CAZymes, and toxins required for pathogenicity in all hosts. They also had antibiotic resistance, pigments to suppress or evade host defense responses, as well as genes for adaptation to different environmental conditions, including temperature, oxidation, and nutrients. These findings provide a better understanding of the virulence, antibiotic resistance, and host adaptation strategies of Ewingella, and they also contribute to the development of effective control strategies.

The Potato Yam Phyllosphere Ectosymbiont Paraburkholderia sp. Msb3 Is a Potent Growth Promotor in Tomato

The genus Paraburkholderia includes a variety of species with promising features for sustainable biotechnological solutions in agriculture through increasing crop productivity. Here, we present a novel Paraburkholderia isolate, a permanent and predominant member of the Dioscoreae bulbifera (yam family, Dioscoreaceae) phyllosphere, making up to 25% of the microbial community on leaf acumens. The 8.5 Mbp genome of isolate Msb3 encodes an unprecedented combination of features mediating a beneficial plant-associated lifestyle, including biological nitrogen fixation (BNF), plant hormone regulation, detoxification of various xenobiotics, degradation of aromatic compounds and multiple protein secretion systems including both T3SS and T6SS. The isolate exhibits significant growth promotion when applied to agriculturally important plants such as tomato, by increasing the total dry biomass by up to 40%. The open question about the “beneficial” nature of this strain led us to investigate ecological and generic boundaries in Burkholderia sensu lato. In a refined phylogeny including 279 Burkholderia sensu lato isolates strain Msb3 clusters within Clade I Paraburkholderia, which also includes few opportunistic strains that can potentially act as pathogens, as revealed by our ecological meta-data analysis. In fact, we demonstrate that all genera originating from the “plant beneficial and environmental” (PBE) Burkholderia species cluster include opportunists. This indicates that further functional examinations are needed before safe application of these strains in sustainable agricultural settings can be assured.

Phylogenetic profiling, an untapped resource for the prediction of secreted proteins and its complementation with sequence-based classifiers in bacterial type III, IV and VI secretion systems

In the establishment and maintenance of the interaction between pathogenic or symbiotic bacteria with a eukaryotic organism, protein substrates of specialized bacterial secretion systems called effectors play a critical role once translocated into the host cell. Proteins are also secreted to the extracellular medium by free-living bacteria or directly injected into other competing organisms to hinder or kill. In this work, we explore an approach based on the evolutionary dependence that most of the effectors maintain with their specific secretion system that analyzes the co-occurrence of any orthologous protein group and their corresponding secretion system across multiple genomes. We compared and complemented our methodology with sequence-based machine learning prediction tools for the type III, IV and VI secretion systems. Finally, we provide the predictive results for the three secretion systems in 1606 complete genomes at http://www.iib.unsam.edu.ar/orgsissec/.

Draft Genome Sequence of Pseudarthrobacter phenanthrenivorans Strain MHSD1, a Bacterial Endophyte Isolated From the Medicinal Plant Pellaea calomelanos

Pseudarthrobacter phenanthrenivorans strain MHSD1 is a bacterial endophyte isolated from sterilized leaves of Pellaea calomelanos, a medicinal plant capable of growing in arid environments. Here, we report the draft genome sequence and annotation of this bacterial endophyte. The draft genome sequence of P. phenanthrenivorans strain MHSD1 has 4 450 468 bp with a G + C content of 65.30%. The National Center for Biotechnology Information Prokaryotic Genome Annotation Pipeline identified a total of 4004 protein-coding genes, 56 genes coding for RNAs, and 82 pseudogenes. Biosynthesis pathways for various phytohormones such as auxin, salicylic acid, ethylene, cytokinin, jasmonic acid, abscisic acid, and gibberellins were identified. Putative genes involved in various characteristics of bacterial endophyte lifestyle such as transport, motility, adhesion, membrane proteins, secretion and delivery systems, plant cell wall modification, and detoxification were identified. Phylogenomic analysis showed P. phenanthrenivorans strain MHSD1 to be a subspecies of P. phenanthrenivorans Sphe3.

A matter of background: DNA repair pathways as a possible cause for the sparse distribution of CRISPR-Cas systems in bacteria

The absence of CRISPR-Cas systems in more than half of the sequenced bacterial genomes is intriguing, because their role in adaptive immunity and their frequent transfer between species should have made them almost ubiquitous, as is the case in Archaea. Here, we investigate the possibility that the success of CRISPR-Cas acquisition by horizontal gene transfer is affected by the interactions of these systems with the host genetic background and especially with components of double-strand break repair systems (DSB-RS). We first described the distribution of systems specialized in the repair of double-strand breaks in Bacteria: homologous recombination and non-homologous end joining. This allowed us to show that such systems are more often positively or negatively correlated with the frequency of CRISPR-Cas systems than random genes of similar frequency. The detailed analysis of these co-occurrence patterns shows that our method identifies previously known cases of mechanistic interactions between these systems. It also reveals other positive and negative patterns of co-occurrence between DSB-RS and CRISPR-Cas systems. Notably, it shows that the patterns of distribution of CRISPR-Cas systems in Proteobacteria are strongly dependent on the epistatic groups including RecBCD and AddAB. Our results suggest that the genetic background plays an important role in the success of adaptive immunity in different bacterial clades and provide insights to guide further experimental research on the interactions between CRISPR-Cas and DSB-RS. This article is part of a discussion meeting issue 'The ecology and evolution of prokaryotic CRISPR-Cas adaptive immune systems'.

The effect of bacterial mutation rate on the evolution of CRISPR-Cas adaptive immunity

CRISPR-Cas immune systems are present in around half of bacterial genomes. Given the specificity and adaptability of this immune mechanism, it is perhaps surprising that they are not more widespread. Recent insights into the requirement for specific host factors for the function of some CRISPR-Cas subtypes, as well as the negative epistasis between CRISPR-Cas and other host genes, have shed light on potential reasons for the partial distribution of this immune strategy in bacteria. In this study, we examined how mutations in the bacterial mismatch repair system, which are frequently observed in natural and clinical isolates and cause elevated host mutation rates, influence the evolution of CRISPR-Cas–mediated immunity. We found that hosts with a high mutation rate very rarely evolved CRISPR-based immunity to phage compared to wild-type hosts. We explored the reason for this effect and found that the higher frequency at which surface mutants pre-exist in the mutator host background causes them to rapidly become the dominant phenotype under phage infection. These findings suggest that natural variation in bacterial mutation rates may, therefore, influence the distribution of CRISPR-Cas adaptive immune systems.

This article is part of a discussion meeting issue ‘The ecology and evolution of prokaryotic CRISPR-Cas adaptive immune systems’.

Marine Sediments Illuminate Chlamydiae Diversity and Evolution

The bacterial phylum Chlamydiae is so far composed of obligate symbionts of eukaryotic hosts. Well known for Chlamydiaceae, pathogens of humans and other animals, Chlamydiae also include so-called environmental lineages that primarily infect microbial eukaryotes. Environmental surveys indicate that Chlamydiae are found in a wider range of environments than anticipated previously. However, the vast majority of this chlamydial diversity has been underexplored, biasing our current understanding of their biology, ecological importance, and evolution. Here, we report that previously undetected and active chlamydial lineages dominate microbial communities in deep anoxic marine sediments taken from the Arctic Mid-Ocean Ridge. Reaching relative abundances of up to 43% of the bacterial community, and a maximum diversity of 163 different species-level taxonomic units, these Chlamydiae represent important community members. Using genome-resolved metagenomics, we reconstructed 24 draft chlamydial genomes, expanding by over a third the known genomic diversity in this phylum. Phylogenomic analyses revealed several novel clades across the phylum, including a previously unknown sister lineage of the Chlamydiaceae, providing new insights into the origin of pathogenicity in this family. We were unable to identify putative eukaryotic hosts for these marine sediment chlamydiae, despite identifying genomic features that may be indicative of host-association. The high abundance and genomic diversity of Chlamydiae in these anoxic marine sediments indicate that some members could play an important, and thus far overlooked, ecological role in such environments and may indicate alternate lifestyle strategies.

Comparative Genomic Analysis Confirms Five Genetic Populations of the Select Agent, Rathayibacter toxicus

Rathayibacter toxicus is a Gram-positive, nematode-vectored bacterium that infects several grass species in the family Poaceae. Unique in its genus, R. toxicus has the smallest genome, possesses a complete CRISPR-Cas system, a vancomycin-resistance cassette, produces tunicamycin, a corynetoxin responsible for livestock deaths in Australia, and is designated a Select Agent in the United States. In-depth, genome-wide analyses performed in this study support the previously designated five genetic populations, with a core genome comprising approximately 80% of the genome for all populations. Results varied as a function of the type of analysis and when using different bioinformatics tools for the same analysis; e.g., some programs failed to identify specific genomic regions that were actually present. The software variance highlights the need to verify bioinformatics results by additional methods; e.g., PCR, mapping genes to genomes, use of multiple algorithms). These analyses suggest the following relationships among populations: RT-IV ↔ RT-I ↔ RT-II ↔ RT-III ↔ RT-V, with RT-IV and RT-V being the most unrelated. This is the most comprehensive analysis of R. toxicus that included populations RT-I and RT-V. Future studies require underrepresented populations and more recent isolates from varied hosts and geographic locations.

Host Range and Genetic Plasticity Explain the Coexistence of Integrative and Extrachromosomal Mobile Genetic Elements

Self-transmissible mobile genetic elements drive horizontal gene transfer between prokaryotes. Some of these elements integrate in the chromosome, whereas others replicate autonomously as plasmids. Recent works showed the existence of few differences, and occasional interconversion, between the two types of elements. Here, we enquired on why evolutionary processes have maintained the two types of mobile genetic elements by comparing integrative and conjugative elements (ICE) with extrachromosomal ones (conjugative plasmids) of the highly abundant MPFT conjugative type. We observed that plasmids encode more replicases, partition systems, and antibiotic resistance genes, whereas ICEs encode more integrases and metabolism-associated genes. ICEs and plasmids have similar average sizes, but plasmids are much more variable, have more DNA repeats, and exchange genes more frequently. On the other hand, we found that ICEs are more frequently transferred between distant taxa. We propose a model where the different genetic plasticity and amplitude of host range between elements explain the co-occurrence of integrative and extrachromosomal elements in microbial populations. In particular, the conversion from ICE to plasmid allows ICE to be more plastic, while the conversion from plasmid to ICE allows the expansion of the element's host range.

Genomic diversity in flavobacterial pathogens of aquatic origin

Flavobacterium species are considered important fish pathogens in wild and cultured fish throughout the world. They can cause acute, subacute, and chronic infections, which are mainly characterized by gill damage, skin lesions, and deep necrotic ulcerations. Primarily, three Flavobacterium species, F. branchiophilum, F. columnare, and F. psychrophilum, have been reported to cause substantial losses to freshwater fish. In this study, we evaluated genomes of 86 Flavobacterium species isolated from aquatic hosts (mainly fish) to identify their unique and shared genome features. Our results showed that F. columnare genomes cluster into four different genetic groups. In silico secretion system analysis identified that all genomes carry type I (T1SS) and type IX (T9SS) secretion systems, but the number of type I secretion system genes shows diversity between species. F. branchiophilum, F. araucananum, F. chilense, F. spartansii, and F. tructae genomes have full type VI secretion system (T6SS). F. columnare, F. hydatis, and F. plurextorum carry partial T6SS with some of the T6SS genes missing. F. columnare, F. araucananum, F. chilense, F. spartansii, F. araucananum, F. tructae, Flavobacterium sp., F. crassostreae, F. succinicans, F. hydatis, and F. plurextorum carry most of the type IV secretion system genes (T4SS). F. columnare genetic groups 1 and 2, Flavobacterium sp., and F. crassostreae encode the least number of antibiotic resistance elements. F. hydatis, F. chilense, and F. plurextorum encode the greatest number of antibiotic resistance genes. Additionally, F. spartansii, F. araucananum, and chilense encode the greatest number of virulence genes while Flavobacterium sp. and F. crassostreae encode the least number of virulence genes. In conclusion, comparative genomics of Flavobacterium species of aquatic origin will help our understanding of Flavobacterium pathogenesis.

Atypical organizations and epistatic interactions of CRISPRs and cas clusters in genomes and their mobile genetic elements

Prokaryotes use CRISPR-Cas systems for adaptive immunity, but the reasons for the frequent existence of multiple CRISPRs and cas clusters remain poorly understood. Here, we analysed the joint distribution of CRISPR and cas genes in a large set of fully sequenced bacterial genomes and their mobile genetic elements. Our analysis suggests few negative and many positive epistatic interactions between Cas subtypes. The latter often result in complex genetic organizations, where a locus has a single adaptation module and diverse interference mechanisms that might provide more effective immunity. We typed CRISPRs that could not be unambiguously associated with a cas cluster and found that such complex loci tend to have unique type I repeats in multiple CRISPRs. Many chromosomal CRISPRs lack a neighboring Cas system and they often have repeats compatible with the Cas systems encoded in trans. Phages and 25 000 prophages were almost devoid of CRISPR-Cas systems, whereas 3% of plasmids had CRISPR-Cas systems or isolated CRISPRs. The latter were often compatible with the chromosomal cas clusters, suggesting that plasmids can co-opt the latter. These results highlight the importance of interactions between CRISPRs and cas present in multiple copies and in distinct genomic locations in the function and evolution of bacterial immunity.

Highly Reduced Genomes of Protist Endosymbionts Show Evolutionary Convergence

Genome evolution in bacterial endosymbionts is notoriously extreme: the combined effects of strong genetic drift and unique selective pressures result in highly reduced genomes with distinctive adaptations to hosts [1-4]. These processes are mostly known from animal endosymbionts, where nutritional endosymbioses represent the best-studied systems. However, eukaryotic microbes, or protists, also harbor diverse bacterial endosymbionts, but their genome reduction and functional relationships with their hosts are largely unexplored [5-7]. We sequenced the genomes of four bacterial endosymbionts from three species of diplonemids, poorly studied but abundant and diverse heterotrophic protists [8-12]. The endosymbionts come from two bacterial families, Rickettsiaceae and Holosporaceae, that have invaded two families of diplonemids, and their genomes have converged on an extremely small size (605-632 kilobase pairs [kbp]), similar gene content (e.g., metabolite transporters and secretion systems), and reduced metabolic potential (e.g., loss of energy metabolism). These characteristics are generally found in both families, but the diplonemid endosymbionts have evolved greater extremes in parallel. They possess modified type VI secretion systems that could function in manipulating host metabolism or other intracellular interactions. Finally, modified cellular machinery like the ATP synthase without oxidative phosphorylation, and the reduced flagellar apparatus present in some diplonemid endosymbionts and nutritional animal endosymbionts, indicates that intracellular mechanisms have converged in bacterial endosymbionts with various functions and from different eukaryotic hosts across the tree of life.

CRISPRCasdb a successor of CRISPRdb containing CRISPR arrays and cas genes from complete genome sequences, and tools to download and query lists of repeats and spacers

In Archaea and Bacteria, the arrays called CRISPRs for 'clustered regularly interspaced short palindromic repeats' and the CRISPR associated genes or cas provide adaptive immunity against viruses, plasmids and transposable elements. Short sequences called spacers, corresponding to fragments of invading DNA, are stored in-between repeated sequences. The CRISPR-Cas systems target sequences homologous to spacers leading to their degradation. To facilitate investigations of CRISPRs, we developed 12 years ago a website holding the CRISPRdb. We now propose CRISPRCasdb, a completely new version giving access to both CRISPRs and cas genes. We used CRISPRCasFinder, a program that identifies CRISPR arrays and cas genes and determine the system's type and subtype, to process public whole genome assemblies. Strains are displayed either in an alphabetic list or in taxonomic order. The database is part of the CRISPR-Cas++ website which also offers the possibility to analyse submitted sequences and to download programs. A BLAST search against lists of repeats and spacers extracted from the database is proposed. To date, 16 990 complete prokaryote genomes (16 650 bacteria from 2973 species and 340 archaea from 300 species) are included. CRISPR-Cas systems were found in 36% of Bacteria and 75% of Archaea strains. CRISPRCasdb is freely accessible at https://crisprcas.i2bc.paris-saclay.fr/.

Identifying Conjugative Plasmids and Integrative Conjugative Elements with CONJscan

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.

Comparative Genomic Analysis of Citrobacter and Key Genes Essential for the Pathogenicity of Citrobacter koseri

Citrobacter species are opportunistic bacterial pathogens that have been implicated in both nosocomial and community-acquired infections. Among the genus Citrobacter, Citrobacter koseri is often isolated from clinical material, and has been known to cause meningitis and brain abscess in neonates and immunocompromised individuals. The virulence determinants of Citrobacter, however, remain largely unknown. Based on traditional methods, the genus Citrobacter has been divided into 11 species, but this has been problematic. Here, we determined an improved, detailed, and more accurate phylogeny of the genus Citrobacter based on whole genome sequence (WGS) data from 129 Citrobacter genomes, 31 of which were sequenced in this study. A maximum likelihood (ML) phylogeny constructed with core genome single-nucleotide polymorphisms (SNPs) classified all Citrobacter isolates into 11 distinct groups, with all C. koseri strains clustering into a single group. For comprehensive and systematic comparative genomic analyses, we investigated the distribution of virulence factors, resistance genes, and macromolecular secretion systems among the Citrobacter genus. Moreover, combined with group-specific genes analysis, we identified a key gene cluster for iron transport, which is present in the C. koseri group, but absent in other the groups, suggesting that the high-pathogenicity island (HPI) cluster may be important for the pathogenicity of C. koseri. Animal experiments showed that loss of the HPI cluster significantly decreased C. koseri virulence in mice and rat. Further, we provide evidence to explain why Citrobacter freundii is less susceptible than C. koseri to several antibiotics in silico. Overall, our data reveal novel virulence clusters specific to the predominantly pathogenic C. koseri strains, which form the basis for elucidating the virulence mechanisms underlying these important pathogens.

Implications of Mobile Genetic Elements for Salmonella enterica Single-Nucleotide Polymorphism Subtyping and Source Tracking Investigations

Single-nucleotide polymorphisms (SNPs) are widely used for whole-genome sequencing (WGS)-based subtyping of foodborne pathogens in outbreak and source tracking investigations. Mobile genetic elements (MGEs) are commonly present in bacterial genomes and may affect SNP subtyping results if their evolutionary history and dynamics differ from that of the bacterial chromosomes. Using Salmonella enterica as a model organism, we surveyed major categories of MGEs, including plasmids, phages, insertion sequences, integrons, and integrative and conjugative elements (ICEs), in 990 genomes representing 21 major serotypes of S. enterica We evaluated whether plasmids and chromosomal MGEs affect SNP subtyping with 9 outbreak clusters of different serotypes found in the United States in 2018. The median total length of chromosomal MGEs accounted for 2.5% of a typical S. enterica chromosome. Of the 990 analyzed S. enterica isolates, 68.9% contained at least one assembled plasmid sequence. The median total length of assembled plasmids in these isolates was 93,671 bp. Plasmids that carry high densities of SNPs were found to substantially affect both SNP phylogenies and SNP distances among closely related isolates if they were present in the reference genome for SNP subtyping. In comparison, chromosomal MGEs were found to have limited impact on SNP subtyping. We recommend the identification of plasmid sequences in the reference genome and the exclusion of plasmid-borne SNPs from SNP subtyping analysis.IMPORTANCE Despite increasingly routine use of WGS and SNP subtyping in outbreak and source tracking investigations, whether and how MGEs affect SNP subtyping has not been thoroughly investigated. Besides chromosomal MGEs, plasmids are frequently entangled in draft genome assemblies and yet to be assessed for their impact on SNP subtyping. This study provides evidence-based guidance on the treatment of MGEs in SNP analysis for Salmonella to infer phylogenetic relationship and SNP distance between isolates.

Draft Genome Sequence of Enterobacter hormaechei Strain MHSD6, a Plant Endophyte Isolated from Medicinal Plant Pellaea calomelanos

We describe here the draft genome sequence of Enterobacter hormaechei strain MHSD6, a bacterial endophyte isolated from the medicinal plant Pellaea calomelanos. The Enterobacter hormaechei strain MHSD6 draft genome is 4,817,102 bp in length, with a G+C content of 55.50%.

We describe here the draft genome sequence of Enterobacter hormaechei strain MHSD6, a bacterial endophyte isolated from the medicinal plant Pellaea calomelanos. The Enterobacter hormaechei strain MHSD6 draft genome is 4,817,102 bp in length, with a G+C content of 55.50%.

Ecological and functional capabilities of an uncultured Kordia sp

Cultivable bacteria represent only a fraction of the diversity in microbial communities. However, the official procedures for classification and characterization of a novel prokaryotic species still rely on isolates. Nevertheless, due to single cell genomics, it is possible to retrieve genomes from environmental samples by sequencing them individually, and to assign specific genes to a specific taxon, regardless of their ability to grow in culture. In this study, a complete description was performed for uncultured Kordia sp. TARA_039_SRF, a proposed novel species within the genus Kordia, using culture-independent techniques. The type material was a high-quality draft genome (94.97% complete, 4.65% gene redundancy) co-assembled using ten nearly identical single amplified genomes (SAGs) from surface seawater in the North Indian Ocean during the Tara Oceans Expedition. The assembly process was optimized to obtain the best possible assembly metrics and a less fragmented genome. The closest relative of the species was Kordia periserrulae, which shared 97.56% similarity of the 16S rRNA gene, 75% orthologs and 89.13% average nucleotide identity. The functional potential of the proposed novel species included proteorhodopsin, the ability to incorporate nitrate, cytochrome oxidases with high affinity for oxygen, and CAZymes that were unique features within the genus. Its abundance at different depths and size fractions was also evaluated together with its functional annotation, revealing that its putative ecological niche could be particles of phytoplanktonic origin. It could putatively attach to these particles and consume them while sinking to the deeper and oxygen depleted layers of the North Indian Ocean.

Comparative genomics reveals a novel genetic organization of the sad cluster in the sulfonamide-degrader 'Candidatus Leucobacter sulfamidivorax' strain GP

Background

Microbial communities recurrently establish metabolic associations resulting in increased fitness and ability to perform complex tasks, such as xenobiotic degradation. In a previous study, we have described a sulfonamide-degrading consortium consisting of a novel low-abundant actinobacterium, named strain GP, and Achromobacter denitrificans PR1. However, we found that strain GP was unable to grow independently and could not be further purified.

Results

Previous studies suggested that strain GP might represent a new putative species within the Leucobacter genus (16S rRNA gene similarity < 97%). In this study, we found that average nucleotide identity (ANI) with other Leucobacter spp. ranged between 76.8 and 82.1%, further corroborating the affiliation of strain GP to a new provisional species. The average amino acid identity (AAI) and percentage of conserved genes (POCP) values were near the lower edge of the genus delimitation thresholds (65 and 55%, respectively). Phylogenetic analysis of core genes between strain GP and Leucobacter spp. corroborated these findings. Comparative genomic analysis indicates that strain GP may have lost genes related to tetrapyrrole biosynthesis and thiol transporters, both crucial for the correct assembly of cytochromes and aerobic growth. However, supplying exogenous heme and catalase was insufficient to abolish the dependent phenotype. The actinobacterium harbors at least two copies of a novel genetic element containing a sulfonamide monooxygenase (sadA) flanked by a single IS1380 family transposase. Additionally, two homologs of sadB (4-aminophenol monooxygenase) were identified in the metagenome-assembled draft genome of strain GP, but these were not located in the vicinity of sadA nor of mobile or integrative elements.

Conclusions

Comparative genomics of the genus Leucobacter suggested the absence of some genes encoding for important metabolic traits in strain GP. Nevertheless, although media and culture conditions were tailored to supply its potential metabolic needs, these conditions were insufficient to isolate the PR1-dependent actinobacterium further. This study gives important insights regarding strain GP metabolism; however, gene expression and functional studies are necessary to characterize and further isolate strain GP. Based on our data, we propose to classify strain GP in a provisional new species within the genus Leucobacter, ‘Candidatus Leucobacter sulfamidivorax‘.

Population Gene Introgression and High Genome Plasticity for the Zoonotic Pathogen Streptococcus agalactiae

The influence that bacterial adaptation (or niche partitioning) within species has on gene spillover and transmission among bacterial populations occupying different niches is not well understood. Streptococcus agalactiae is an important bacterial pathogen that has a taxonomically diverse host range making it an excellent model system to study these processes. Here, we analyze a global set of 901 genome sequences from nine diverse host species to advance our understanding of these processes. Bayesian clustering analysis delineated 12 major populations that closely aligned with niches. Comparative genomics revealed extensive gene gain/loss among populations and a large pan genome of 9,527 genes, which remained open and was strongly partitioned among niches. As a result, the biochemical characteristics of 11 populations were highly distinctive (significantly enriched). Positive selection was detected and biochemical characteristics of the dispensable genes under selection were enriched in ten populations. Despite the strong gene partitioning, phylogenomics detected gene spillover. In particular, tetracycline resistance (which likely evolved in the human-associated population) from humans to bovine, canines, seals, and fish, demonstrating how a gene selected in one host can ultimately be transmitted into another, and biased transmission from humans to bovines was confirmed with a Bayesian migration analysis. Our findings show high bacterial genome plasticity acting in balance with selection pressure from distinct functional requirements of niches that is associated with an extensive and highly partitioned dispensable genome, likely facilitating continued and expansive adaptation.

Pan-genome analysis of Riemerella anatipestifer reveals its genomic diversity and acquired antibiotic resistance associated with genomic islands

Riemerella anatipestifer is a gram-negative bacterium that leads to severe contagious septicemia in ducks, turkeys, chickens, and wild waterfowl. Here, a pan-genome with 32 R. anatipestifer genomes is re-established, and the mathematical model is calculated to evaluate the expansion of R. anatipestifer genomes, which were determined to be open. Average nucleotide identity (ANI) and phylogenetic analysis preliminarily clarify intraspecies variation and distance. Comparative genomic analysis of R. anatipestifer found that horizontal gene transfer events, which provide an expressway for the recruitment of novel functionalities and facilitate genetic diversity in microbial genomes, play a key role in the process of acquiring and transmitting antibiotic-resistance genes in R. anatipestifer. Furthermore, a new antibiotic-resistance gene cluster was identified in the same loci in 14 genomes. The uneven distribution of virulence factors was also confirmed by our results. Our study suggests that the ability to acquire foreign genes (such as antibiotic-resistance genes) increases the adaptability of R. anatipestifer, and the virulence genes with little mobility are highly conserved in R. anatipestifer.

Comparative genomic analysis of the Hafnia genus reveals an explicit evolutionary relationship between the species alvei and paralvei and provides insights into pathogenicity

Background

The Hafnia genus is an opportunistic pathogen that has been implicated in both nosocomial and community-acquired infections. Although Hafnia is fairly often isolated from clinical material, its taxonomy has remained an unsolved riddle, and the involvement and importance of Hafnia in human disease is also uncertain. Here, we used comparative genomic analysis to define the taxonomy of Hafnia, identify species-specific genes that may be the result of ecological and pathogenic specialization, and reveal virulence-related genetic profiles that may contribute to pathogenesis.

Results

One complete genome sequence and 19 draft genome sequences for Hafnia strains were generated and combined with 27 publicly available genomes. We provided high-resolution typing methods by constructing phylogeny and population structure based on single-copy core genes in combination with whole genome average nucleotide identity to identify two distant Hafnia species (alvei and paralvei) and one mislabeled strain. The open pan-genome and the presence of numerous mobile genetic elements reveal that Hafnia has undergone massive gene rearrangements. Presence of species-specific core genomes associated with metabolism and transport suggests the putative niche differentiation between alvei and paralvei. We also identified possession of diverse virulence-related profiles in both Hafnia species., including the macromolecular secretion system, virulence, and antimicrobial resistance. In the macromolecular system, T1SS, Flagellum 1, Tad pilus and T6SS-1 were conserved in Hafnia, whereas T4SS, T5SS, and other T6SSs exhibited the evolution of diversity. The virulence factors in Hafnia are related to adherence, toxin, iron uptake, stress adaptation, and efflux pump. The identified resistance genes are associated with aminoglycoside, beta-lactam, bacitracin, cationic antimicrobial peptide, fluoroquinolone, and rifampin. These virulence-related profiles identified at the genomic level provide insights into Hafnia pathogenesis and the differentiation between alvei and paralvei.

Conclusions

Our research using core genome phylogeny and comparative genomics analysis of a larger collection of strains provides a comprehensive view of the taxonomy and species-specific traits between Hafnia species. Deciphering the genome of Hafnia strains possessing a reservoir of macromolecular secretion systems, virulence factors, and resistance genes related to pathogenicity may provide insights into addressing its numerous infections and devising strategies to combat the pathogen.

Comparative genomics of Aeromonas veronii: Identification of a pathotype impacting aquaculture globally

Aeromonas veronii is a gram-negative species abundant in aquatic environments that causes disease in humans as well as terrestrial and aquatic animals. In the current study, 41 publicly available A. veronii genomes were compared to investigate distribution of putative virulence genes, global dissemination of pathotypes, and potential mechanisms of virulence. The complete genome of A. veronii strain ML09-123 from an outbreak of motile aeromonas septicemia in farm-raised catfish in the southeastern United States was included. Dissemination of A. veronii strain types was discovered in dispersed geographical locations. Isolate ML09-123 is highly similar to Chinese isolate TH0426, suggesting the two strains have a common origin and may represent a pathotype impacting aquaculture in both countries. Virulence of strain ML09-123 in catfish in a dose-dependent manner was confirmed experimentally. Subsystem category disposition showed the majority of genomes exhibit similar distribution of genomic elements. The type I secretion system (T1SS), type II secretion system (T2SS), type 4 pilus (T4P), and flagellum core elements are conserved in all A. veronii genomes, whereas the type III secretion system (T3SS), type V secretion system (T5SS), type VI secretion system (T6SS), and tight adherence (TAD) system demonstrate variable dispersal. Distribution of mobile elements is dependent on host and geographic origin, suggesting this species has undergone considerable genetic exchange. The data presented here lends insight into the genomic variation of A. veronii and identifies a pathotype impacting aquaculture globally.

Evolution of a family of molecular Rube Goldberg contraptions

Case studies of the evolution of molecular machines remain scarce. One of the most diverse and widespread homologous families of machines is the type IV filament (TFF) superfamily, comprised of type IV pili, type II secretion systems (T2SSs), archaella, and other less-well-characterized families. These families have functions including twitching motility, effector export, rotary propulsion, nutrient uptake, DNA uptake, and even electrical conductance, but it is unclear how such diversity evolved from a common ancestor. In this issue, Denise and colleagues take a significant step toward understanding evolution of the TFF superfamily by determining a global phylogeny and using it to infer an evolutionary pathway. Results reveal that the superfamily predates the divergence of Bacteria and Archaea, and show how duplications, acquisitions, and losses coincide with changes in function. Surprises include that tight adherence (Tad) pili were horizontally acquired from Archaea and that T2SSs were relatively recently repurposed from type IV pili. Results also enable better understanding of the function of the ATPase family that powers the superfamily. The study highlights the role of tinkering by exaptation-the repurposing of pre-existing functions for new roles-in the diversification of molecular machines.

Comparative genome analysis of Lactobacillus mudanjiangensis, an understudied member of the Lactobacillus plantarum group

The genus Lactobacillus is known to be extremely diverse and consists of different phylogenetic groups that show a diversity that is roughly equal to the expected diversity of a typical bacterial genus. One of the most prominent phylogenetic groups within this genus is the Lactobacillus plantarum group, which contains the understudied Lactobacillus mudanjiangensis species. Before this study, only one L. mudanjiangensis strain, DSM 28402^T, had been described, but without whole-genome analysis. In this study, three strains classified as L. mudanjiangensis were isolated from three different carrot juice fermentations and their whole-genome sequence was determined, together with the genome sequence of the type strain. The genomes of all four strains were compared with publicly available L. plantarum group genome sequences. This analysis showed that L. mudanjiangensis harboured the second largest genome size and gene count of the whole L. plantarum group. In addition, all members of this species showed the presence of a gene coding for a cellulose-degrading enzyme. Finally, three of the four L. mudanjiangensis strains studied showed the presence of pili on scanning electron microscopy (SEM) images, which were linked to conjugative gene regions, coded on a plasmid in at least two of the strains studied.

Just the Two of Us? A Family of Pseudomonas Megaplasmids Offers a Rare Glimpse into the Evolution of Large Mobile Elements

Pseudomonads are ubiquitous group of environmental proteobacteria, well known for their roles in biogeochemical cycling, in the breakdown of xenobiotic materials, as plant growth promoters, and as pathogens of a variety of host organisms. We have previously identified a large megaplasmid present within one isolate of the plant pathogen Pseudomonas syringae, and here we report that a second member of this megaplasmid family is found within an environmental Pseudomonad isolate most closely related to Pseudomonas putida. Many of the shared genes are involved in critical cellular processes like replication, transcription, translation, and DNA repair. We argue that presence of these shared pathways sheds new light on discussions about the types of genes that undergo horizontal gene transfer (i.e., the complexity hypothesis) as well as the evolution of pangenomes. Furthermore, although both megaplasmids display a high level of synteny, genes that are shared differ by over 50% on average at the amino acid level. This combination of conservation in gene order despite divergence in gene sequence suggests that this Pseudomonad megaplasmid family is relatively old, that gene order is under strong selection within this family, and that there are likely many more members of this megaplasmid family waiting to be found in nature.

Ubiquinone Biosynthesis over the Entire O2 Range: Characterization of a Conserved O2-Independent Pathway

In order to colonize environments with large O₂ gradients or fluctuating O₂ 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₂. In contrast, the long-described pathway for ubiquinone biosynthesis requires O₂ as a substrate. In fact, we find that many proteobacteria are equipped with the O₂-dependent and O₂-independent pathways, supporting that they are able to synthesize ubiquinone over the entire O₂ range. Overall, we propose that the novel O₂-independent pathway is part of the metabolic plasticity developed by proteobacteria to face various environmental O₂ levels.

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₂) 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₂-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₂-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₂. Together, the O₂-dependent and O₂-independent ubiquinone biosynthesis pathways contribute to optimizing bacterial metabolism over the entire O₂ range.

CasLocusAnno: a web-based server for annotating cas loci and their corresponding (sub)types

In prokaryotes, Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) and CRISPR-associated protein (Cas) systems constitute adaptive immune systems against mobile genetic elements (MGEs). Here, we introduce the Markov cluster algorithm (MCL) to Makarova et al.'s method in order to select a more reasonable profile. Additionally, our new Maximum Continuous Cas Subcluster (MCCS) method helps identification of tightly clustered loci. The comparison with two other commonly used programs shows that the method could identify Cas proteins with higher accuracy and lower Additional Prediction Rate (APR). Moreover, we developed a web-based server, CasLocusAnno (http://cefg.uestc.cn/CasLocusAnno), capable of annotating Cas proteins, cas loci and their (sub)types less than ~ 28 s following the whole proteome sequence submission. Its standalone version can be downloaded at https://github.com/RiversDong/CasLocusAnno.

Diversification of the type IV filament superfamily into machines for adhesion, protein secretion, DNA uptake, and motility

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.

First Report of Integrative Conjugative Elements in Riemerella anatipestifer Isolates From Ducks in China

We report for the first time the occurrence of integrative conjugative elements (ICEs) in Riemerella anatipestifer (R.anatipestifer) isolated from diseased ducks in China. For this purpose, a total of 48 genome sequences were investigated, which comprised 30 publicly available R. anatipestifer genome sequences, and 18 clinical isolates genomes sequences. Two ICEs, named ICERanRCAD0133-1 and ICERanRCAD0179-1 following the classic nomenclature system, were identified in R. anatipestifer through the use of bioinformatics tools. Comparative analysis revealed that three ICEs in Ornithobacterium rhinotracheale showed a high degree of conservation with the core genes of ICERanRCAD0133-1, while 13 ICEs with high similarity to ICERanRCAD0179-1 were found in Bacteroidetes. Based on the definition of ICE family, ICERanRCAD0179-1 was grouped in CTnDOT/ERL family; however, ICERanRCAD0133-1, which had no significant similarity with known ICEs, might be classified into a novel ICE family. The sequences of ICERanRCAD0133-1 and ICERanRCAD0179-1 were 70890 bp and 49166 bp in length, had 33.14 and 50.34% GC content, and contained 77 CDSs and 51 CDSs, respectively. Cargo genes carried by these two ICEs were predicted to encode: R-M systems, IS elements, a putative TonB-dependent receptor, a bacteriocin/lantibiotic efflux ABC transporter, a tetracycline resistance gene and more. In addition, phylogenetic analyses revealed that ICERanRCAD0179-1 and related ICEs were derived from a common ancestor, which may have undergone divergence prior to integartation into the host bacterial chromosome, and that the core genes co-evolved via a related evolutionary process or experienced only a low degree of recombination events during spread from a common CTnDOT/ERL family ancestor. Collectively, this study is the first identification and characterization of ICEs in R. anatipestifer; and provides new insights into the genetic diversity, evolution, adaptation, antimicrobial resistance, and virulence of R. anatipestifer.

Insight into the genome and brackish water adaptation strategies of toxic and bloom-forming Baltic Sea Dolichospermum sp. UHCC 0315

The Baltic Sea is a shallow basin of brackish water in which the spatial salinity gradient is one of the most important factors contributing to species distribution. The Baltic Sea is infamous for its annual cyanobacterial blooms comprised of Nodularia spumigena, Aphanizomenon spp., and Dolichospermum spp. that cause harm, especially for recreational users. To broaden our knowledge of the cyanobacterial adaptation strategies for brackish water environments, we sequenced the entire genome of Dolichospermum sp. UHCC 0315, a species occurring not only in freshwater environments but also in brackish water. Comparative genomics analyses revealed a close association with Dolichospermum sp. UHCC 0090 isolated from a lake in Finland. The genome closure of Dolichospermum sp. UHCC 0315 unraveled a mixture of two subtypes in the original culture, and subtypes exhibited distinct buoyancy phenotypes. Salinity less than 3 g L-1 NaCl enabled proper growth of Dolichospermum sp. UHCC 0315, whereas growth was arrested at moderate salinity (6 g L-1 NaCl). The concentrations of toxins, microcystins, increased at moderate salinity, whereas RNA sequencing data implied that Dolichospermum remodeled its primary metabolism in unfavorable high salinity. Based on our results, the predicted salinity decrease in the Baltic Sea may favor toxic blooms of Dolichospermum spp.

Wide diversity of methane and short-chain alkane metabolisms in uncultured archaea

Methanogenesis is an ancient metabolism of key ecological relevance, with direct impact on the evolution of Earth’s climate. Recent results suggest that the diversity of methane metabolisms and their derivations have probably been vastly underestimated. Here, by probing thousands of publicly available metagenomes for homologues of methyl-coenzyme M reductase complex (MCR), we have obtained ten metagenome-assembled genomes (MAGs) belonging to potential methanogenic, anaerobic methanotrophic and short-chain alkane oxidizing archaea. Five of these MAGs represent under-sampled (e.g., Verstraetearchaeota, Methanonatronarchaeia, ANME-1) or previously genomically undescribed (ANME-2c) archaeal lineages. The remaining five MAGs correspond to lineages that are only distantly related to previously known methanogens and span the entire archaeal phylogeny. Comprehensive comparative annotation significantly expands the metabolic diversity and energy conservation systems of MCR-bearing archaea. It also suggests the potential existence of a yet uncharacterized type of methanogenesis linked to short-chain alkane/fatty acid oxidation in a previously undescribed class of archaea (‘Ca. Methanoliparia’). We redefine a common core of marker genes specific to methanogenic, anaerobic methanotrophic and short-chain alkane-oxidizing archaea, and propose a possible scenario for the evolutionary and functional transitions that led to the emergence of such metabolic diversity.

Comparative genomics of Paraburkholderia kururiensis and its potential in bioremediation, biofertilization, and biocontrol of plant pathogens

Burkholderia harbors versatile Gram-negative species and is β-Proteobacteria. Recently, it was proposed to split the genus in two main branches: one of animal and plant pathogens and another, Paraburkholderia, harboring environmental and plant-beneficial species. Currently, Paraburkholderia comprises more than 70 species with ability to occupy very diverse environmental niches. Herein, we sequenced and analyzed the genome of Paraburkholderia kururiensis type strain KP23^T , and compared to P. kururiensis M130, isolated in Brazil, and P. kururiensis susbp. thiooxydans, from Korea. This study focused on the gene content of the three genomes with special emphasis on their potential of plant-association, biocontrol, and bioremediation. The comparative analyses revealed several genes related to plant benefits, including biosynthesis of IAA, ACC deaminase, multiple efflux pumps, dioxygenases, and degradation of aromatic compounds. Importantly, a range of genes for protein secretion systems (type III, IV, V, and VI) were characterized, potentially involved in P. kururiensis well documented ability to establish endophytic association with plants. These findings shed light onto bacteria-plant interaction mechanisms at molecular level, adding novel information that supports their potential application in bioremediation, biofertilization, and biocontrol of plant pathogens. P. kururiensis emerges as a promising model to investigate adaptation mechanisms in different ecological niches.

Comparative Genomics of Aeromonas hydrophila Secretion Systems and Mutational Analysis of hcp1 and vgrG1 Genes From T6SS

Virulent Aeromonas hydrophila causes severe motile Aeromonas septicemia in warmwater fishes. In recent years, channel catfish farming in the U.S.A. and carp farming in China have been affected by virulent A. hydrophila, and genome comparisons revealed that these virulent A. hydrophila strains belong to the same clonal group. Bacterial secretion systems are often important virulence factors; in the current study, we investigated whether secretion systems contribute to the virulent phenotype of these strains. Thus, we conducted comparative secretion system analysis using 55 A. hydrophila genomes, including virulent A. hydrophila strains from U.S.A. and China. Interestingly, tight adherence (TaD) system is consistently encoded in all the vAh strains. The majority of U.S.A. isolates do not possess a complete type VI secretion system, but three core elements [tssD (hcp), tssH, and tssI (vgrG)] are encoded. On the other hand, Chinese isolates have a complete type VI secretion system operon. None of the virulent A. hydrophila isolates have a type III secretion system. Deletion of two genes encoding type VI secretion system proteins (hcp1 and vgrG1) from virulent A. hydrophila isolate ML09-119 reduced virulence 2.24-fold in catfish fingerlings compared to the parent strain ML09-119. By determining the distribution of genes encoding secretion systems in A. hydrophila strains, our study clarifies which systems may contribute to core A. hydrophila functions and which may contribute to more specialized adaptations such as virulence. Our study also clarifies the role of type VI secretion system in A. hydrophila virulence.

Genetic exchanges are more frequent in bacteria encoding capsules

Capsules allow bacteria to colonize novel environments, to withstand numerous stresses, and to resist antibiotics. Yet, even though genetic exchanges with other cells should be adaptive under such circumstances, it has been suggested that capsules lower the rates of homologous recombination and horizontal gene transfer. We analysed over one hundred pan-genomes and thousands of bacterial genomes for the evidence of an association between genetic exchanges (or lack thereof) and the presence of a capsule system. We found that bacteria encoding capsules have larger pan-genomes, higher rates of horizontal gene transfer, and higher rates of homologous recombination in their core genomes. Accordingly, genomes encoding capsules have more plasmids, conjugative elements, transposases, prophages, and integrons. Furthermore, capsular loci are frequent in plasmids, and can be found in prophages. These results are valid for Bacteria, independently of their ability to be naturally transformable. Since we have shown previously that capsules are commonly present in nosocomial pathogens, we analysed their co-occurrence with antibiotic resistance genes. Genomes encoding capsules have more antibiotic resistance genes, especially those encoding efflux pumps, and they constitute the majority of the most worrisome nosocomial bacteria. We conclude that bacteria with capsule systems are more genetically diverse and have fast-evolving gene repertoires, which may further contribute to their success in colonizing novel niches such as humans under antibiotic therapy.

Previous works showed that bacteria encoding capsules are better colonizers and are dominant in most environments suggesting a positive role for capsules in the genetic diversification of bacteria. Yet, it has been repeatedly suggested, based almost exclusively studies in few model species, that such bacteria are less diverse and engage in fewer genetic exchanges. Here, we reverse the current paradigm and show that bacteria encoding capsules have larger and more diverse gene repertoires, which change faster by horizontal gene transfer and recombination. Our study alters the traditional view of the capsule as a barrier to gene flow and raises novel questions about the role of capsules in bacterial adaptation.

The Genome Sequence of "Candidatus Fokinia solitaria": Insights on Reductive Evolution in Rickettsiales

"Candidatus Fokinia solitaria" is an obligate intracellular endosymbiont of a unicellular eukaryote, a ciliate of the genus Paramecium. Here, we present the genome sequence of this bacterium and subsequent analysis. Phylogenomic analysis confirmed the previously reported positioning of the symbiont within the "Candidatus Midichloriaceae" family (order Rickettsiales), as well as its high sequence divergence from other members of the family, indicative of fast sequence evolution. Consistently with this high evolutionary rate, a comparative genomic analysis revealed that the genome of this symbiont is the smallest of the Rickettsiales to date. The reduced genome does not present flagellar genes, nor the pathway for the biosynthesis of lipopolysaccharides (present in all the other so far sequenced members of the family "Candidatus Midichloriaceae") or genes for the Krebs cycle (present, although not always complete, in Rickettsiales). These results indicate an evolutionary trend toward a stronger dependence on the host, in comparison with other members of the family. Two alternative scenarios are compatible with our results; "Candidatus Fokinia solitaria" could be either a recently evolved, vertically transmitted mutualist, or a parasite with a high host-specificity.

BIOINFORMATIC SEARCH OF CRISPR/CAS SYSTEM STRUCTURES IN GENOME OF PCT281 PLASMID OF BACILLUS THURINGIENSIS SUBSP. CHINENSIS STRAIN CT-43

Background. CRISPR/Cas systems loci are one of the functionally important patterns in bacterial genome which perform the role of “adaptive immune defense” from foreign nucleic acids. The study of CRISPR/Cas systems structure in genomes of plasmids and phages provide new information about the evolution of this systems in bacterial hosts.

Aims. A search of CRISPR/Cas systems structures in pCT281 plasmid from Bacillus thuringiensis subsp. chinensis strain CT-43 using bioinformatic methods.

Materials and methods. Search studies using bioinformatics methods were performed with the genome of pCT281 plasmid of B. thuringiensis subsp. chinensis strain CT-43 from the RefSeq database. To search for the CRISPR/Cas system structure MacSyFinder (ver. 1.0.5) and three combined algorithms were used: CRISPRFinder; PILER-CR; CRISPR Recognition Tool (CRT). The consensus repeat sequence was generated in WebLogo 3.

Results and discussion. In pCT281 plasmid we detected one locus of CRISPR/Cas system of the type I-C which contains 2 CRISPR-cassettes and 4 cas-genes located between them. The CRISPR-cassette 1 includes 10 spacers from 32 to 35 bp and 11 repeats 32bp in length. 5 spacers (33–35 bp) separated by 6 repeats 32 bp in length were detected in the CRISPR-cassette 2.

Conclusions. The bioinformatic methods used in this study enable to conduct a search of CRISPR/Cas systems structures in plasmid genomes. The presence of the CRISPR-Cas locus in pCT281 plasmid confirms a possible transfer of this system from the nucleoid to this plasmid. The detected spacers provide information about phages this bacteria was encountered.

CRISPRCasFinder, an update of CRISRFinder, includes a portable version, enhanced performance and integrates search for Cas proteins

CRISPR (clustered regularly interspaced short palindromic repeats) arrays and their associated (Cas) proteins confer bacteria and archaea adaptive immunity against exogenous mobile genetic elements, such as phages or plasmids. CRISPRCasFinder allows the identification of both CRISPR arrays and Cas proteins. The program includes: (i) an improved CRISPR array detection tool facilitating expert validation based on a rating system, (ii) prediction of CRISPR orientation and (iii) a Cas protein detection and typing tool updated to match the latest classification scheme of these systems. CRISPRCasFinder can either be used online or as a standalone tool compatible with Linux operating system. All third-party software packages employed by the program are freely available. CRISPRCasFinder is available at https://crisprcas.i2bc.paris-saclay.fr.

Candidatus Nitrosocaldus cavascurensis, an Ammonia Oxidizing, Extremely Thermophilic Archaeon with a Highly Mobile Genome

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 amoAXCB 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.

Inhibition of NHEJ repair by type II-A CRISPR-Cas systems in bacteria

Type II CRISPR-Cas systems introduce double-strand breaks into DNA of invading genetic material and use DNA fragments to acquire novel spacers during adaptation. These breaks can be the substrate of several DNA repair pathways, paving the way for interactions. We report that non-homologous end-joining (NHEJ) and type II-A CRISPR-Cas systems only co-occur once among 5563 fully sequenced prokaryotic genomes. We investigated experimentally the possible molecular interactions using the NHEJ pathway from Bacillus subtilis and the type II-A CRISPR-Cas systems from Streptococcus thermophilus and Streptococcus pyogenes. Our results suggest that the NHEJ system has no effect on CRISPR immunity. On the other hand, we provide evidence for the inhibition of NHEJ repair by the Csn2 protein. Our findings give insights on the complex interactions between CRISPR-Cas systems and repair mechanisms in bacteria, contributing to explain the scattered distribution of CRISPR-Cas systems in bacterial genome.