Comparative genome analyses of Mycobacteroides immunogenum reveals two potential novel subspecies

Mycobacteroides immunogenum is an emerging opportunistic pathogen implicated in nosocomial infections. Comparative genome analyses may provide better insights into its genomic structure, functions and evolution. The present analysis showed that M. immunogenum has an open pan-genome. Approximately 36.8% of putative virulence genes were identified in the accessory regions of M. immunogenum. Phylogenetic analyses revealed two potential novel subspecies of M. immunogenum, supported by evidence from ANIb (average nucleotide identity using blast) and GGDC (Genome to Genome Distance Calculator) analyses. We identified 74 genomic islands (GIs) in Subspecies 1 and 23 GIs in Subspecies 2. All Subspecies 2-harboured GIs were not found in Subspecies 1, indicating that they might have been acquired by Subspecies 2 after their divergence. Subspecies 2 has more defence genes than Subspecies 1, suggesting that it might be more resistant to the insertion of foreign DNA and probably explaining why Subspecies 2 has fewer GIs. Positive selection analysis suggest that M. immunogenum has a lower selection pressure compared to non-pathogenic mycobacteria. Thirteen genes were positively selected and many were involved in virulence.

Effect of arsenite and growth in biofilm conditions on the evolution of Thiomonas sp. CB2

Thiomonas bacteria are ubiquitous at acid mine drainage sites and play key roles in the remediation of water at these locations by oxidizing arsenite to arsenate, favouring the sorption of arsenic by iron oxides and their coprecipitation. Understanding the adaptive capacities of these bacteria is crucial to revealing how they persist and remain active in such extreme conditions. Interestingly, it was previously observed that after exposure to arsenite, when grown in a biofilm, some strains of Thiomonas bacteria develop variants that are more resistant to arsenic. Here, we identified the mechanisms involved in the emergence of such variants in biofilms. We found that the percentage of variants generated increased in the presence of high concentrations of arsenite (5.33 mM), especially in the detached cells after growth under biofilm-forming conditions. Analysis of gene expression in the parent strain CB2 revealed that genes involved in DNA repair were upregulated in the conditions where variants were observed. Finally, we assessed the phenotypes and genomes of the subsequent variants generated to evaluate the number of mutations compared to the parent strain. We determined that multiple point mutations accumulated after exposure to arsenite when cells were grown under biofilm conditions. Some of these mutations were found in what is referred to as ICE19, a genomic island (GI) carrying arsenic-resistance genes, also harbouring characteristics of an integrative and conjugative element (ICE). The mutations likely favoured the excision and duplication of this GI. This research aids in understanding how Thiomonas bacteria adapt to highly toxic environments, and, more generally, provides a window to bacterial genome evolution in extreme environments.

Metagenome-assembled genome binning methods with short reads disproportionately fail for plasmids and genomic Islands

Metagenomic methods enable the simultaneous characterization of microbial communities without time-consuming and bias-inducing culturing. Metagenome-assembled genome (MAG) binning methods aim to reassemble individual genomes from this data. However, the recovery of mobile genetic elements (MGEs), such as plasmids and genomic islands (GIs), by binning has not been well characterized. Given the association of antimicrobial resistance (AMR) genes and virulence factor (VF) genes with MGEs, studying their transmission is a public-health priority. The variable copy number and sequence composition of MGEs makes them potentially problematic for MAG binning methods. To systematically investigate this issue, we simulated a low-complexity metagenome comprising 30 GI-rich and plasmid-containing bacterial genomes. MAGs were then recovered using 12 current prediction pipelines and evaluated. While 82-94 % of chromosomes could be correctly recovered and binned, only 38-44 % of GIs and 1-29 % of plasmid sequences were found. Strikingly, no plasmid-borne VF nor AMR genes were recovered, and only 0-45 % of AMR or VF genes within GIs. We conclude that short-read MAG approaches, without further optimization, are largely ineffective for the analysis of mobile genes, including those of public-health importance, such as AMR and VF genes. We propose that researchers should explore developing methods that optimize for this issue and consider also using unassembled short reads and/or long-read approaches to more fully characterize metagenomic data.

Evolutionary Mechanisms of Long-Term Genome Diversification Associated With Niche Partitioning in Marine Picocyanobacteria

Marine picocyanobacteria of the genera Prochlorococcus and Synechococcus are the most abundant photosynthetic organisms on Earth, an ecological success thought to be linked to the differential partitioning of distinct ecotypes into specific ecological niches. However, the underlying processes that governed the diversification of these microorganisms and the appearance of niche-related phenotypic traits are just starting to be elucidated. Here, by comparing 81 genomes, including 34 new Synechococcus, we explored the evolutionary processes that shaped the genomic diversity of picocyanobacteria. Time-calibration of a core-protein tree showed that gene gain/loss occurred at an unexpectedly low rate between the different lineages, with for instance 5.6 genes gained per million years (My) for the major Synechococcus lineage (sub-cluster 5.1), among which only 0.71/My have been fixed in the long term. Gene content comparisons revealed a number of candidates involved in nutrient adaptation, a large proportion of which are located in genomic islands shared between either closely or more distantly related strains, as identified using an original network construction approach. Interestingly, strains representative of the different ecotypes co-occurring in phosphorus-depleted waters (Synechococcus clades III, WPC1, and sub-cluster 5.3) were shown to display different adaptation strategies to this limitation. In contrast, we found few genes potentially involved in adaptation to temperature when comparing cold and warm thermotypes. Indeed, comparison of core protein sequences highlighted variants specific to cold thermotypes, notably involved in carotenoid biosynthesis and the oxidative stress response, revealing that long-term adaptation to thermal niches relies on amino acid substitutions rather than on gene content variation. Altogether, this study not only deciphers the respective roles of gene gains/losses and sequence variation but also uncovers numerous gene candidates likely involved in niche partitioning of two key members of the marine phytoplankton.

Molecular insights into the genome dynamics and interactions between core and acquired genomes of Vibrio cholerae

Bacterial species are hosts to horizontally acquired mobile genetic elements (MGEs), which encode virulence, toxin, antimicrobial resistance, and other metabolic functions. The bipartite genome of Vibrio cholerae harbors sporadic and conserved MGEs that contribute in the disease development and survival of the pathogens. For a comprehensive understanding of dynamics of MGEs in the bacterial genome, we engineered the genome of V. cholerae and examined in vitro and in vivo stability of genomic islands (GIs), integrative conjugative elements (ICEs), and prophages. Recombinant vectors carrying the integration module of these GIs, ICE and CTXΦ, helped us to understand the efficiency of integrations of MGEs in the V. cholerae chromosome. We have deleted more than 250 acquired genes from 6 different loci in the V. cholerae chromosome and showed contribution of CTX prophage in the essentiality of SOS response master regulator LexA, which is otherwise not essential for viability in other bacteria, including Escherichia coli In addition, we observed that the core genome-encoded RecA helps CTXΦ to bypass V. cholerae immunity and allow it to replicate in the host bacterium in the presence of similar prophage in the chromosome. Finally, our proteomics analysis reveals the importance of MGEs in modulating the levels of cellular proteome. This study engineered the genome of V. cholerae to remove all of the GIs, ICEs, and prophages and revealed important interactions between core and acquired genomes.

Antibiotic Resistance in Vibrio cholerae: Mechanistic Insights from IncC Plasmid-Mediated Dissemination of a Novel Family of Genomic Islands Inserted at trmE

The increasing association of the etiological agent of cholera, Vibrio cholerae serogroup O1 and O139, with multiple antibiotic resistance threatens to deprive health practitioners of this effective tool. Drug resistance in cholera results mainly from acquisition of mobile genetic elements. Genomic islands conferring multidrug resistance and mobilizable by IncC conjugative plasmids were reported to circulate in non-O1/non-O139 V. cholerae clinical strains isolated from the 2010 Haitian cholera outbreak. As these genomic islands can be transmitted to pandemic V. cholerae serogroups, their mechanism of transmission needed to be investigated. Our research revealed plasmid- and genomic island-encoded factors required for the resistance island excision, mobilization, and integration, as well as regulation of these functions. The discovery of related genomic islands carrying diverse phage resistance genes but lacking antibiotic resistance-conferring genes in a wide range of marine dwelling bacteria suggests that these elements are ancient and recently acquired drug resistance genes.

Cholera remains a formidable disease, and reports of multidrug-resistant strains of the causative agent Vibrio cholerae have become common during the last 3 decades. The pervasiveness of resistance determinants has largely been ascribed to mobile genetic elements, including SXT/R391 integrative conjugative elements, IncC plasmids, and genomic islands (GIs). Conjugative transfer of IncC plasmids is activated by the master activator AcaCD whose regulatory network extends to chromosomally integrated GIs. MGIVchHai6 is a multidrug resistance GI integrated at the 3′ end of trmE (mnmE or thdF) in chromosome 1 of non-O1/non-O139 V. cholerae clinical isolates from the 2010 Haitian cholera outbreak. In the presence of an IncC plasmid expressing AcaCD, MGIVchHai6 excises from the chromosome and transfers at high frequency. Herein, the mechanism of mobilization of MGIVchHai6 GIs by IncC plasmids was dissected. Our results show that AcaCD drives expression of GI-borne genes, including xis and mobI_M, involved in excision and mobilization. A 49-bp fragment upstream of mobI_M was found to serve as the minimal origin of transfer (oriT) of MGIVchHai6. The direction of transfer initiated at oriT was determined using IncC plasmid-driven mobilization of chromosomal markers via MGIVchHai6. In addition, IncC plasmid-encoded factors, including the relaxase TraI, were found to be required for GI transfer. Finally, in silico exploration of Gammaproteobacteria genomes identified 47 novel related and potentially AcaCD-responsive GIs in 13 different genera. Despite sharing conserved features, these GIs integrate at trmE, yicC, or dusA and carry a diverse cargo of genes involved in phage resistance.

IMPORTANCE The increasing association of the etiological agent of cholera, Vibrio cholerae serogroup O1 and O139, with multiple antibiotic resistance threatens to deprive health practitioners of this effective tool. Drug resistance in cholera results mainly from acquisition of mobile genetic elements. Genomic islands conferring multidrug resistance and mobilizable by IncC conjugative plasmids were reported to circulate in non-O1/non-O139 V. cholerae clinical strains isolated from the 2010 Haitian cholera outbreak. As these genomic islands can be transmitted to pandemic V. cholerae serogroups, their mechanism of transmission needed to be investigated. Our research revealed plasmid- and genomic island-encoded factors required for the resistance island excision, mobilization, and integration, as well as regulation of these functions. The discovery of related genomic islands carrying diverse phage resistance genes but lacking antibiotic resistance-conferring genes in a wide range of marine dwelling bacteria suggests that these elements are ancient and recently acquired drug resistance genes.

Multiple chromosomal inversions contribute to adaptive divergence of a dune sunflower ecotype

Both models and case studies suggest that chromosomal inversions can facilitate adaptation and speciation in the presence of gene flow by suppressing recombination between locally adapted alleles. Until recently, however, it has been laborious and time-consuming to identify and genotype inversions in natural populations. Here we apply RAD sequencing data and newly developed population genomic approaches to identify putative inversions that differentiate a sand dune ecotype of the prairie sunflower (Helianthus petiolaris) from populations found on the adjacent sand sheet. We detected seven large genomic regions that exhibit a different population structure than the rest of the genome and that vary in frequency between dune and nondune populations. These regions also show high linkage disequilibrium and high heterozygosity between, but not within, arrangements, consistent with the behaviour of large inversions, an inference subsequently validated in part by comparative genetic mapping. Genome-environment association analyses show that key environmental variables, including vegetation cover and soil nitrogen, are significantly associated with inversions. The inversions colocate with previously described "islands of differentiation," and appear to play an important role in adaptive divergence and incipient speciation within H. petiolaris.

Metabolic remodeling in Escherichia coli MG1655. A prophage e14-encoded small RNA, co293, post-transcriptionally regulates transcription factors HcaR and FadR

In previous studies, we have shown the existence of metabolic remodeling in glucose-grown Escherichia coli MG1655 cells expressing the esterase Orf306 from the opd island of Sphingobium fuliginis. We now show that Orf306-dependent metabolic remodeling is due to regulation of a novel small RNA (sRNA). Endogenous propionate, produced due to the esterase/lipase activity of Orf306, repressed expression of a novel E. coli sRNA, co293. This sRNA post-transcriptionally regulates expression of the transcription factors HcaR and FadR either by inhibiting translation or by destabilizing their transcripts. Hence, repression of co293 expression elevates the levels of HcaR and FadR with consequent activation of alternative carbon catabolic pathways. HcaR activates the hca and MHP operons leading to upregulation of the phenyl propionate and hydroxy phenyl propionate (HPP) degradation pathways. Similarly, FadR stimulates the expression of the transcription factor IclR which negatively regulates the glyoxylate bypass pathway genes, aceBAK.

Horizontally Acquired Homologs of Xenogeneic Silencers: Modulators of Gene Expression Encoded by Plasmids, Phages and Genomic Islands

Acquisition of mobile elements by horizontal gene transfer can play a major role in bacterial adaptation and genome evolution by providing traits that contribute to bacterial fitness. However, gaining foreign DNA can also impose significant fitness costs to the host bacteria and can even produce detrimental effects. The efficiency of horizontal acquisition of DNA is thought to be improved by the activity of xenogeneic silencers. These molecules are a functionally related group of proteins that possess affinity for the acquired DNA. Binding of xenogeneic silencers suppresses the otherwise uncontrolled expression of genes from the newly acquired nucleic acid, facilitating their integration to the bacterial regulatory networks. Even when the genes encoding for xenogeneic silencers are part of the core genome, homologs encoded by horizontally acquired elements have also been identified and studied. In this article, we discuss the current knowledge about horizontally acquired xenogeneic silencer homologs, focusing on those encoded by genomic islands, highlighting their distribution and the major traits that allow these proteins to become part of the host regulatory networks.

Comparative Genomics of Streptococcus thermophilus Support Important Traits Concerning the Evolution, Biology and Technological Properties of the Species

Streptococcus thermophilus is a major starter for the dairy industry with great economic importance. In this study we analyzed 23 fully sequenced genomes of S. thermophilus to highlight novel aspects of the evolution, biology and technological properties of this species. Pan/core genome analysis revealed that the species has an important number of conserved genes and that the pan genome is probably going to be closed soon. According to whole genome phylogeny and average nucleotide identity (ANI) analysis, most S. thermophilus strains were grouped in two major clusters (i.e., clusters A and B). More specifically, cluster A includes strains with chromosomes above 1.83 Mbp, while cluster B includes chromosomes below this threshold. This observation suggests that strains belonging to the two clusters may be differentiated by gene gain or gene loss events. Furthermore, certain strains of cluster A could be further subdivided in subgroups, i.e., subgroup I (ASCC 1275, DGCC 7710, KLDS SM, MN-BM-A02, and ND07), II (MN-BM-A01 and MN-ZLW-002), III (LMD-9 and SMQ-301), and IV (APC151 and ND03). In cluster B certain strains formed one distinct subgroup, i.e., subgroup I (CNRZ1066, CS8, EPS, and S9). Clusters and subgroups observed for S. thermophilus indicate the existence of lineages within the species, an observation which was further supported to a variable degree by the distribution and/or the architecture of several genomic traits. These would include exopolysaccharide (EPS) gene clusters, Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs)-CRISPR associated (Cas) systems, as well as restriction-modification (R-M) systems and genomic islands (GIs). Of note, the histidine biosynthetic cluster was found present in all cluster A strains (plus strain NCTC12958^T) but was absent from all strains in cluster B. Other loci related to lactose/galactose catabolism and urea metabolism, aminopeptidases, the majority of amino acid and peptide transporters, as well as amino acid biosynthetic pathways were found to be conserved in all strains suggesting their central role for the species. Our study highlights the necessity of sequencing and analyzing more S. thermophilus complete genomes to further elucidate important aspects of strain diversity within this starter culture that may be related to its application in the dairy industry.

Transcriptomic Analysis of Aggregatibacter actinomycetemcomitans Core and Accessory Genes in Different Growth Conditions

Aggregatibacter actinomycetemcomitans genome can be divided into an accessory gene pool (found in some but not all strains) and a core gene pool (found in all strains). The functions of the accessory genes (genomic islands and non-island accessory genes) are largely unknown. We hypothesize that accessory genes confer critical functions for A. actinomycetemcomitans in vivo. This study examined the expression patterns of accessory and core genes of A. actinomycetemcomitans in distinct growth conditions. We found similar expression patterns of island and non-island accessory genes, which were generally lower than the core genes in all growth conditions. The median expression levels of genomic islands were 29%–37% of the core genes in enriched medium but elevated to as high as 63% of the core genes in nutrient-limited media. Several putative virulence genes, including the cytolethal distending toxin operon, were found to be activated in nutrient-limited conditions. In conclusion, genomic islands and non-island accessory genes exhibited distinct patterns of expression from the core genes and may play a role in the survival of A. actinomycetemcomitans in nutrient-limited environments.

Targeted transposition with Tn7 elements: safe sites, mobile plasmids, CRISPR/Cas and beyond

Transposon Tn7 is notable for the control it exercises over where transposition events are directed. One Tn7 integration pathways recognizes a highly conserved attachment (att) site in the chromosome, while a second pathway specifically recognizes mobile plasmids that facilitate transfer of the element to new hosts. In this review, I discuss newly discovered families of Tn7-like elements with different targeting pathways. Perhaps the most exciting examples are multiple instances where Tn7-like elements have repurposed CRISPR/Cas systems. In these cases, the CRISPR/Cas systems have lost their canonical defensive function to destroy incoming mobile elements; instead, the systems have been naturally adapted to use guide RNAs to specifically direct transposition into these mobile elements. The new families of Tn7-like elements also include a variety of novel att sites in bacterial chromosomes where genome islands can form. Interesting families have also been revealed where proteins described in the prototypic Tn7 element are fused or otherwise repurposed for the new dual activities. This expanded understanding of Tn7-like elements broadens our view of how genetic systems are repurposed and provides potentially exciting new tools for genome modification and genomics. Future opportunities and challenges to understanding the impact of the new families of Tn7-like elements are discussed.

Characterization and Comparative Analysis of the Staphylococcus aureus Genomic Island vSaβ: an In Silico Approach

Staphylococcus aureus is a widespread opportunistic pathogen to humans and animals. Of its genome, 20 to 25% varies between strains and consists of phages, pathogenicity islands, transposons, and genomic islands. S. aureus harbors up to three genomic islands, vSaα, vSaβ, and vSaγ. The vSaβ region of S. aureus can encode a number of virulence-associated factors, such as serine proteases, leukocidins, enterotoxins, bacteriocins, or a hyaluronate lyase. In this study, the vSaβ regions of 103 clinically relevant S. aureus strains were characterized in silico and compared to the three predefined vSaβ types. We here suggest a superordinate system of 15 different vSaβ types, of which 12 were newly defined. Each vSaβ type has a distinct structure with a distinct set of genes, which are both highly conserved. Between the different types, gene content and composition vary substantially. Based on our data, a strain's vSaβ type is strongly coupled with its clonal complex, suggesting that vSaβ was acquired in an ancestral S. aureus strain, arguably by phage mediation, before differentiation into clonal complexes. In addition, we addressed the issue of ambiguous nomenclature in the serine protease gene cluster and propose a novel, phylogeny-based nomenclature of the cluster contained in the vSaβ region.IMPORTANCE With the rapid increase of available sequencing data on clinically relevant bacterial species such as S. aureus, the genomic basis of clinical phenotypes can be investigated in much more detail, allowing a much deeper understanding of the mechanisms involved in disease. We characterized in detail the S. aureus genomic island vSaβ and defined a superordinate system to categorize S. aureus strains based on their vSaβ type, providing information about the strains' virulence-associated genes and clinical potential.

Comparative Analysis of Genomic Island Prediction Tools

Tools for genomic island prediction use strategies for genomic comparison analysis and sequence composition analysis. The goal of comparative analysis is to identify unique regions in the genomes of related organisms, whereas sequence composition analysis evaluates and relates the composition of specific regions with other regions in the genome. The goal of this study was to qualitatively and quantitatively evaluate extant genomic island predictors. We chose tools reported to produce significant results using sequence composition prediction, comparative genomics, and hybrid genomics methods. To maintain diversity, the tools were applied to eight complete genomes of organisms with distinct characteristics and belonging to different families. Escherichia coli CFT073 was used as a control and considered as the gold standard because its islands were previously curated in vitro. The results of predictions with the gold standard were manually curated, and the content and characteristics of each predicted island were analyzed. For other organisms, we created GenBank (GBK) files using Artemis software for each predicted island. We copied only the amino acid sequences from the coding sequence and constructed a multi-FASTA file for each predictor. We used BLASTp to compare all results and generate hits to evaluate similarities and differences among the predictions. Comparison of the results with the gold standard revealed that GIPSy produced the best results, covering ~91% of the composition and regions of the islands, followed by Alien Hunter (81%), IslandViewer (47.8%), Predict Bias (31%), GI Hunter (17%), and Zisland Explorer (16%). The tools with the best results in the analyzes of the set of organisms were the same ones that presented better performance in the tests with the gold standard.

Unintentional Genomic Changes Endow Cupriavidus metallidurans with an Augmented Heavy-Metal Resistance

For the past three decades, Cupriavidus metallidurans has been one of the major model organisms for bacterial tolerance to heavy metals. Its type strain CH34 contains at least 24 gene clusters distributed over four replicons, allowing for intricate and multilayered metal responses. To gain organic mercury resistance in CH34, broad-spectrum mer genes were introduced in a previous work via conjugation of the IncP-1β plasmid pTP6. However, we recently noted that this CH34-derived strain, MSR33, unexpectedly showed an increased resistance to other metals (i.e., Co2+, Ni2+, and Cd2+). To thoroughly investigate this phenomenon, we resequenced the entire genome of MSR33 and compared its DNA sequence and basal gene expression profile to those of its parental strain CH34. Genome comparison identified 11 insertions or deletions (INDELs) and nine single nucleotide polymorphisms (SNPs), whereas transcriptomic analysis displayed 107 differentially expressed genes. Sequence data implicated the transposition of IS1088 in higher Co2+ and Ni2+ resistances and altered gene expression, although the precise mechanisms of the augmented Cd2+ resistance in MSR33 remains elusive. Our work indicates that conjugation procedures involving large complex genomes and extensive mobilomes may pose a considerable risk toward the introduction of unwanted, undocumented genetic changes. Special efforts are needed for the applied use and further development of small nonconjugative broad-host plasmid vectors, ideally involving CRISPR-related and advanced biosynthetic technologies.

Genome Structure of the Opportunistic Pathogen Paracoccus yeei (Alphaproteobacteria) and Identification of Putative Virulence Factors

Bacteria of the genus Paracoccus are common components of the microbiomes of many naturally- and anthropogenically shaped environments. One species, Paracoccus yeei, is unique within the genus because it is associated with opportunistic human infections. Therefore, strains of P. yeei may serve as an interesting model to study the transition from a saprophytic to a pathogenic lifestyle in environmental bacteria. Unfortunately, knowledge concerning the biology, genetics and genomic content of P. yeei is fragmentary; also the mechanisms of pathogenicity of this bacterium remain unclear. In this study we provide the first insight into the genome composition and metabolic potential of a clinical isolate, P. yeei CCUG 32053. This strain has a multipartite genome (4,632,079 bp) composed of a circular chromosome plus eight extrachromosomal replicons pYEE1–8: 3 chromids and 5 plasmids, with a total size of 1,247,173 bp. The genome has been significantly shaped by the acquisition of genomic islands, prophages (Myoviridae and Siphoviridae phage families) and numerous insertion sequences (ISs) representing seven IS families. Detailed comparative analysis with other complete genomic sequences of Paracoccus spp. (including P. yeei FDAARGOS_252 and TT13, as well as non-pathogenic strains of other species in this genus) enabled us to identify P. yeei species-specific genes and to predict putative determinants of virulence. This is the first attempt to identify pathoadaptive genetic information of P. yeei and to estimate the role of the mobilome in the evolution of pathogenicity in this species.

Deciphering pathogenicity and antibiotic resistance islands in methicillin-resistant Staphylococcus aureus genomes

Staphylococcus aureus is a versatile pathogen that is capable of causing infections in both humans and animals. It can cause furuncles, septicaemia, pneumonia and endocarditis. Adaptation of S. aureus to the modern hospital environment has been facilitated, in part, by the horizontal acquisition of drug resistance genes, such as mecA gene that imparts resistance to methicillin. Horizontal acquisitions of islands of genes harbouring virulence and antibiotic resistance genes have made S. aureus resistant to commonly used antibiotics. To decipher genomic islands (GIs) in 22 hospital- and 9 community-associated methicillin-resistant S. aureus strains and classify a subset of GIs carrying virulence and resistance genes as pathogenicity and resistance islands respectively, we applied a host of methods for localizing genomic islands in prokaryotic genomes. Surprisingly, none of the frequently used GI prediction methods could perform well in delineating the resistance islands in the S. aureus genomes. Rather, a gene clustering procedure exploiting biases in codon usage for identifying horizontally transferred genes outperformed the current methods for GI detection, in particular in identifying the known islands in S. aureus including the SCCmec island that harbours the mecA resistance gene. The gene clustering approach also identified novel, as yet unreported islands, with many of these found to harbour virulence and/or resistance genes. These as yet unexplored islands may provide valuable information on the evolution of drug resistance in S. aureus.

Insights into the genome diversity and virulence of two clinical isolates of Burkholderia cenocepacia

PURPOSE

Burkholderia cenocepacia is among the most common members of the Burkholderia cepacia complex (Bcc) isolated from patients with cystic fibrosis (CF). The factors triggering the high rates of morbidity and mortality in CF patients are not well elucidated. In this study, we aim to highlight the genome diversity of two clinical isolates of B. cenocepacia through comparative genome analysis.

METHODOLOGY

The repertoire of virulence factors and resistance genes compared to reference strains J2315 and K56-2 was elucidated. The isolates were screened for the presence of phages and insertion sequences. Two methods were combined to obtain an accurate prediction of genomic islands (GIs): the cumulative GC profile and the IslandViewer web tool. To study evolutionary relatedness, whole genome-based single-nucleotide polymorphism (wgSNP) analysis was also performed with 43 publically available strains of the Bcc of various sequence types.Results/Key findings. Genome-based species identification of the two isolates BC-AUH and BC-BMEH confirmed the species as B. cenocepacia. Both belonged to ST-602, a double-locus variant of ST-32 (CC31), genomovar IIIA, and carried a large number of antibiotic resistance genes. Eighteen GIs were predicted in BC-AUH and BC-BMEH, occupying 9.3 and 6.1 % of the respective genomes. Comparison to J2315 revealed 89 and 85 genes unique to BC-BMEH and BC-AUH, respectively. Additionally, 1823 intergenic SNPs were detected between BC-BMEH and BC-AUH.

CONCLUSION

This study mapped existing genetic variations in B. cenocepacia associated with notorious outcomes in CF patients, and the data obtained provide comprehensive, genome-inferred insights and multifactorial examination of an important human pathogen.

Unexpected diversity in the mobilome of a Pseudomonas aeruginosa strain isolated from a dental unit waterline revealed by SMRT Sequencing

The Gram-negative bacterium Pseudomonas aeruginosa is found in several habitats, both natural and human-made, and is particularly known for its recurrent presence as a pathogen in the lungs of patients suffering from cystic fibrosis, a genetic disease. Given its clinical importance, several major studies have investigated the genomic adaptation of P. aeruginosa in lungs and its transition as acute infections become chronic. However, our knowledge about the diversity and adaptation of the P. aeruginosa genome to non-clinical environments is still fragmentary, in part due to the lack of accurate reference genomes of strains from the numerous environments colonized by the bacterium. Here, we used PacBio long-read technology to sequence the genome of PPF-1, a strain of P. aeruginosa isolated from a dental unit waterline. Generating this closed genome was an opportunity to investigate genomic features that are difficult to accurately study in a draft genome (contigs state). It was possible to shed light on putative genomic islands, some shared with other reference genomes, new prophages, and the complete content of insertion sequences. In addition, four different group II introns were also found, including two characterized here and not listed in the specialized group II intron database.

CRISPR-Cas-Mediated Phage Resistance Enhances Horizontal Gene Transfer by Transduction

A powerful contributor to prokaryotic evolution is horizontal gene transfer (HGT) through transformation, conjugation, and transduction, which can be advantageous, neutral, or detrimental to fitness. Bacteria and archaea control HGT and phage infection through CRISPR-Cas (clustered regularly interspaced short palindromic repeats-CRISPR-associated proteins) adaptive immunity. Although the benefits of resisting phage infection are evident, this can come at a cost of inhibiting the acquisition of other beneficial genes through HGT. Despite the ability of CRISPR-Cas to limit HGT through conjugation and transformation, its role in transduction is largely overlooked. Transduction is the phage-mediated transfer of bacterial DNA between cells and arguably has the greatest impact on HGT. We demonstrate that in Pectobacterium atrosepticum, CRISPR-Cas can inhibit the transduction of plasmids and chromosomal loci. In addition, we detected phage-mediated transfer of a large plant pathogenicity genomic island and show that CRISPR-Cas can inhibit its transduction. Despite these inhibitory effects of CRISPR-Cas on transduction, its more common role in phage resistance promotes rather than diminishes HGT via transduction by protecting bacteria from phage infection. This protective effect can also increase transduction of phage-sensitive members of mixed populations. CRISPR-Cas systems themselves display evidence of HGT, but little is known about their lateral dissemination between bacteria and whether transduction can contribute. We show that, through transduction, bacteria can acquire an entire chromosomal CRISPR-Cas system, including cas genes and phage-targeting spacers. We propose that the positive effect of CRISPR-Cas phage immunity on enhancing transduction surpasses the rarer cases where gene flow by transduction is restricted.IMPORTANCE The generation of genetic diversity through acquisition of DNA is a powerful contributor to microbial evolution and occurs through transformation, conjugation, and transduction. Of these, transduction, the phage-mediated transfer of bacterial DNA, is arguably the major route for genetic exchange. CRISPR-Cas adaptive immune systems control gene transfer by conjugation and transformation, but transduction has been mostly overlooked. Our results indicate that CRISPR-Cas can impede, but typically enhances the transduction of plasmids, chromosomal genes, and pathogenicity islands. By limiting wild-type phage replication, CRISPR-Cas immunity increases transduction in both phage-resistant and -sensitive members of mixed populations. Furthermore, we demonstrate mobilization of a chromosomal CRISPR-Cas system containing phage-targeting spacers by generalized transduction, which might partly account for the uneven distribution of these systems in nature. Overall, the ability of CRISPR-Cas to promote transduction reveals an unexpected impact of adaptive immunity on horizontal gene transfer, with broader implications for microbial evolution.

Zisland Explorer: detect genomic islands by combining homogeneity and heterogeneity properties

Genomic islands are genomic fragments of alien origin in bacterial and archaeal genomes, usually involved in symbiosis or pathogenesis. In this work, we described Zisland Explorer, a novel tool to predict genomic islands based on the segmental cumulative GC profile. Zisland Explorer was designed with a novel strategy, as well as a combination of the homogeneity and heterogeneity of genomic sequences. While the sequence homogeneity reflects the composition consistence within each island, the heterogeneity measures the composition bias between an island and the core genome. The performance of Zisland Explorer was evaluated on the data sets of 11 different organisms. Our results suggested that the true-positive rate (TPR) of Zisland Explorer was at least 10.3% higher than that of four other widely used tools. On the other hand, the new tool did not lose overall accuracy with the improvement in the TPR and showed better equilibrium among various evaluation indexes. Also, Zisland Explorer showed better accuracy in the prediction of experimental island data. Overall, the tool provides an alternative solution over other tools, which expands the field of island prediction and offers a supplement to increase the performance of the distinct predicting strategy. We have provided a web service as well as a graphical user interface and open-source code across multiple platforms for Zisland Explorer, which is available at http://cefg.uestc.edu.cn/Zisland_Explorer/ or http://tubic.tju.edu.cn/Zisland_Explorer/.

The Importance of the Expendable: Toxin-Antitoxin Genes in Plasmids and Chromosomes

Toxin–antitoxin (TA) genes were first reported in plasmids and were considered expendable genetic cassettes involved in the stable maintenance of the plasmid replicon by interfering with growth and/or viability of bacteria in which the plasmid was lost. TAs were later found in bacterial chromosomes and also in integrated mobile genetic elements; they were proposed to be involved in the bacterial response to stressful situations. At present, 100s of TAs have been identified and classified in up to six families (I to VI), with those belonging to the type II (constituted by two protein components) being the most studied. Based on well-characterized examples of several type II TAs, we discuss in this review that irrespective of their locations in plasmids or chromosomes, TAs functionally overlap as indicated by: (i) in both locations they can mediate the maintenance of genetic elements to which they are physical linked, and (ii) they can induce persistence or virulence in response to stress situations. Examples of functional confluences in homologous TA systems with different locations are also given. We also consider whether the physiological role of TAs is due to their genetic organization as operons or to their inherent properties, like the short lifespan of the antitoxin component.

Diversity, Prevalence, and Longitudinal Occurrence of Type II Toxin-Antitoxin Systems of Pseudomonas aeruginosa Infecting Cystic Fibrosis Lungs

Type II toxin-antitoxin (TA) systems are most commonly composed of two genes encoding a stable toxin, which harms the cell, and an unstable antitoxin that can inactivate it. TA systems were initially characterized as selfish elements, but have recently gained attention for regulating general stress responses responsible for pathogen virulence, formation of drug-tolerant persister cells and biofilms—all implicated in causing recalcitrant chronic infections. We use a bioinformatics approach to explore the distribution and evolution of type II TA loci of the opportunistic pathogen, Pseudomonas aeruginosa, across longitudinally sampled isolates from cystic fibrosis lungs. We identify their location in the genome, mutations, and gain/loss during infection to elucidate their function(s) in stabilizing selfish elements and pathogenesis. We found (1) 26 distinct TA systems, where all isolates harbor four in their core genome and a variable number of the remaining 22 on genomic islands; (2) limited mutations in core genome TA loci, suggesting they are not under negative selection; (3) no evidence for horizontal transmission of elements with TA systems between clone types within patients, despite their ability to mobilize; (4) no gain and limited loss of TA-bearing genomic islands, and of those elements partially lost, the remnant regions carry the TA systems supporting their role in genomic stabilization; (5) no significant correlation between frequency of TA systems and strain ability to establish as chronic infection, but those with a particular TA, are more successful in establishing a chronic infection.

Identification of Novel Genomic Islands in Liverpool Epidemic Strain of Pseudomonas aeruginosa Using Segmentation and Clustering

Pseudomonas aeruginosa is an opportunistic pathogen implicated in a myriad of infections and a leading pathogen responsible for mortality in patients with cystic fibrosis (CF). Horizontal transfers of genes among the microorganisms living within CF patients have led to highly virulent and multi-drug resistant strains such as the Liverpool epidemic strain of P. aeruginosa, namely the LESB58 strain that has the propensity to acquire virulence and antibiotic resistance genes. Often these genes are acquired in large clusters, referred to as "genomic islands (GIs)." To decipher GIs and understand their contributions to the evolution of virulence and antibiotic resistance in P. aeruginosa LESB58, we utilized a recursive segmentation and clustering procedure, presented here as a genome-mining tool, "GEMINI." GEMINI was validated on experimentally verified islands in the LESB58 strain before examining its potential to decipher novel islands. Of the 6062 genes in P. aeruginosa LESB58, 596 genes were identified to be resident on 20 GIs of which 12 have not been previously reported. Comparative genomics provided evidence in support of our novel predictions. Furthermore, GEMINI unraveled the mosaic structure of islands that are composed of segments of likely different evolutionary origins, and demonstrated its ability to identify potential strain biomarkers. These newly found islands likely have contributed to the hyper-virulence and multidrug resistance of the Liverpool epidemic strain of P. aeruginosa.

A hypervariable genomic island identified in clinical and environmental Mycobacterium avium subsp. hominissuis isolates from Germany

Mycobacterium avium subsp. hominissuis (MAH) is an opportunistic human pathogen widespread in the environment. Genomic islands (GI)s represent a part of the accessory genome of bacteria and influence virulence, drug-resistance or fitness and trigger bacterial evolution. We previously identified a novel GI in four MAH genomes. Here, we further explored this GI in a larger collection of MAH isolates from Germany (n=41), including 20 clinical and 21 environmental isolates. Based on comparative whole genome analysis, we detected this GI in 39/41 (95.1%) isolates. Although all these GIs integrated in the same insertion hotspot, there is high variability in the genetic structure of this GI: eight different types of GI have been identified, designated A-H (sized 6.2-73.3kb). These GIs were arranged as single GI (23/41, 56.1%), combination of two different GIs (14/41, 34.1%) or combination of three different GIs (2/41, 4.9%) in the insertion hotspot. Moreover, two GI types shared more than 80% sequence identity with sequences of M. canettii, responsible for Tuberculosis. A total of 253 different genes were identified in all GIs, among which the previously documented virulence-related genes mmpL10 and mce. The diversity of the GI and the sequence similarity with other mycobacteria suggests cross-species transfer, involving also highly pathogenic species. Shuffling of potential virulence genes such as mmpL10 via this GI may create new pathogens that can cause future outbreaks.