Evidence of a large novel gene pool associated with prokaryotic genomic islands

DICEP: An integrative approach to augmenting genomic island detection

Mobilization of clusters of genes called genomic islands (GIs) across bacterial lineages facilitates dissemination of traits, such as, resistance against antibiotics, virulence or hypervirulence, and versatile metabolic capabilities. Robust delineation of GIs is critical to understanding bacterial evolution that has a vast impact on different life forms. Methods for identification of GIs exploit different evolutionary features or signals encoded within the genomes of bacteria, however, the current state-of-the-art in GI detection still leaves much to be desired. Here, we have taken a combinatorial approach that accounted for GI specific features such as compositional bias, aberrant phyletic pattern, and marker gene enrichment within an integrative framework to delineate GIs in bacterial genomes. Our GI prediction tool, DICEP, was assessed on simulated genomes and well-characterized bacterial genomes. DICEP compared favorably with current GI detection tools on real and synthetic datasets.

Escherichia coli possessing the dihydroxyacetone phosphate shunt utilize 5'-deoxynucleosides for growth

All organisms utilize S-adenosyl-l-methionine (SAM) as a key co-substrate for the methylation of biological molecules, the synthesis of polyamines, and radical SAM reactions. When these processes occur, 5'-deoxy-nucleosides are formed as byproducts such as S-adenosyl-l-homocysteine, 5'-methylthioadenosine (MTA), and 5'-deoxyadenosine (5dAdo). A prevalent pathway found in bacteria for the metabolism of MTA and 5dAdo is the dihydroxyacetone phosphate (DHAP) shunt, which converts these compounds into dihydroxyacetone phosphate and 2-methylthioacetaldehyde or acetaldehyde, respectively. Previous work in other organisms has shown that the DHAP shunt can enable methionine synthesis from MTA or serve as an MTA and 5dAdo detoxification pathway. Rather, the DHAP shunt in Escherichia coli ATCC 25922, when introduced into E. coli K-12, enables the use of 5dAdo and MTA as a carbon source for growth. When MTA is the substrate, the sulfur component is not significantly recycled back to methionine but rather accumulates as 2-methylthioethanol, which is slowly oxidized non-enzymatically under aerobic conditions. The DHAP shunt in ATCC 25922 is active under oxic and anoxic conditions. Growth using 5-deoxy-d-ribose was observed during aerobic respiration and anaerobic respiration with Trimethylamine N-oxide (TMAO), but not during fermentation or respiration with nitrate. This suggests the DHAP shunt may only be relevant for extraintestinal pathogenic E. coli lineages with the DHAP shunt that inhabit oxic or TMAO-rich extraintestinal environments. This reveals a heretofore overlooked role of the DHAP shunt in carbon and energy metabolism from ubiquitous SAM utilization byproducts and suggests a similar role may occur in other pathogenic and non-pathogenic bacteria with the DHAP shunt.

IMPORTANCE

The acquisition and utilization of organic compounds that serve as growth substrates are essential for Escherichia coli to grow and multiply. Ubiquitous enzymatic reactions involving S-adenosyl-l-methionine as a co-substrate by all organisms result in the formation of the 5'-deoxy-nucleoside byproducts, 5'-methylthioadenosine and 5'-deoxyadenosine. All E. coli possess a conserved nucleosidase that cleaves these 5'-deoxy-nucleosides into 5-deoxy-pentose sugars for adenine salvage. The DHAP shunt pathway is found in some extraintestinal pathogenic E. coli, but its function in E. coli possessing it has remained unknown. This study reveals that the DHAP shunt enables the utilization of 5'-deoxy-nucleosides and 5-deoxy-pentose sugars as growth substrates in E. coli strains with the pathway during aerobic respiration and anaerobic respiration with TMAO, but not fermentative growth. This provides an insight into the diversity of sugar compounds accessible by E. coli with the DHAP shunt and suggests that the DHAP shunt is primarily relevant in oxic or TMAO-rich extraintestinal environments.

DMF mineralization and substrate specificity mechanism of Aminobacter ciceronei DMFA1

N,N-dimethylformamide (DMF) is widely used in various industries, but its direct release into water poses high risks to human beings. Although a lot of DMF-degrading bacteria has been isolated, limited studies focus on the degradation preference among DMF and its analogues. In this study, an efficient DMF mineralization bacterium designated Aminobacter ciceronei DMFA1 was isolated from marine sediment. When exposed to a 0.2% DMF (∼1900 mg/L), strain DMFA1 exhibited a degradation efficiency of 100% within 4 days. The observed growth using formamide as the sole carbon source implied the possible DMF degradation pathway of strain DMFA1. Meanwhile,the strain DMFA1 possesses a broad-spectrum substrate degradation, which could effectively degraded 0.2% N,N-dimethylacetamide (DMAC) and N-methylformamide (NMF). Genomic analysis further confirmed the supposed pathway through annotating the genes encoding N, N-dimethylformamidase (DMFase), formamidase, and formate dehydrogenase. The existence of sole DMFase indicating its substrate specificity controlled the preference of DMAc of strain DMFA1. By integrating multiple sequence alignment, homology modeling and molecular docking, the preference of the DMFase in strain DMFA1 towards DMAc are related to: 1) Mutations in key active site residues; 2) the absence of small subunit; and 3) no energy barrier for substrates entering the active site.

Systematic comparison of genome information processing and boundary recognition tools used for genomic island detection

Genomic islands are fragments of foreign DNA that are found in bacterial and archaeal genomes, and are typically associated with symbiosis or pathogenesis. While numerous genomic island detection methods have been proposed, there has been limited evaluation of the efficiency of the genome information processing and boundary recognition tools. In this study, we conducted a review of the statistical methods involved in genomic signatures, host signature extraction, informative signature selection, divergence measures, and boundary detection steps in genomic island prediction. We compared the performances of these methods on simulated experiments using alien fragments obtained from both artificial and real genomes. Our results indicate that among the nine genomic signatures evaluated, genomic signature frequency and full probability performed the best. However, their performance declined when normalized to their expectations and variances, such as Z-score and composition vector. Based on our experiments of the E. coli genome, we found that the confidence intervals of the window variances achieved the best performance in the signature extraction of the host, with the best confidence interval being 1.5-2 times the standard error. Ordered kurtosis was most effective in selecting informative signatures from a single genome, without requiring prior knowledge from other datasets. Among the three divergence measures evaluated, the two-sample t-test was the most successful, and a non-overlapping window with a small eye window (size 2) was best suited for identifying compositionally distinct regions. Finally, the maximum of the Markovian Jensen-Shannon divergence score, in terms of GC-content bias, was found to make boundary detection faster while maintaining a similar error rate.

Utilization of 5'-deoxy-nucleosides as Growth Substrates by Extraintestinal Pathogenic E. coli via the Dihydroxyacetone Phosphate Shunt

All organisms utilize S-adenosyl-l-methionine (SAM) as a key co-substrate for methylation of biological molecules, synthesis of polyamines, and radical SAM reactions. When these processes occur, 5'-deoxy-nucleosides are formed as byproducts such as S-adenosyl-l-homocysteine (SAH), 5'-methylthioadenosine (MTA), and 5'-deoxyadenosine (5dAdo). One of the most prevalent pathways found in bacteria for the metabolism of MTA and 5dAdo is the DHAP shunt, which converts these compounds into dihydroxyacetone phosphate (DHAP) and 2-methylthioacetaldehyde or acetaldehyde, respectively. Previous work has shown that the DHAP shunt can enable methionine synthesis from MTA or serve as an MTA and 5dAdo detoxification pathway. Here we show that in Extraintestinal Pathogenic E. coil (ExPEC), the DHAP shunt serves none of these roles in any significant capacity, but rather physiologically functions as an assimilation pathway for use of MTA and 5dAdo as growth substrates. This is further supported by the observation that when MTA is the substrate for the ExPEC DHAP shunt, the sulfur components is not significantly recycled back to methionine, but rather accumulates as 2-methylthioethanol, which is slowly oxidized non-enzymatically under aerobic conditions. While the pathway is active both aerobically and anaerobically, it only supports aerobic ExPEC growth, suggesting that it primarily functions in oxygenic extraintestinal environments like blood and urine versus the predominantly anoxic gut. This reveals a heretofore overlooked role of the DHAP shunt in carbon assimilation and energy metabolism from ubiquitous SAM utilization byproducts and suggests a similar role may occur in other pathogenic and non-pathogenic bacteria with the DHAP shunt.

The Lost and Found: Unraveling the Functions of Orphan Genes

Orphan Genes (OGs) are a mysterious class of genes that have recently gained significant attention. Despite lacking a clear evolutionary history, they are found in nearly all living organisms, from bacteria to humans, and they play important roles in diverse biological processes. The discovery of OGs was first made through comparative genomics followed by the identification of unique genes across different species. OGs tend to be more prevalent in species with larger genomes, such as plants and animals, and their evolutionary origins remain unclear but potentially arise from gene duplication, horizontal gene transfer (HGT), or de novo origination. Although their precise function is not well understood, OGs have been implicated in crucial biological processes such as development, metabolism, and stress responses. To better understand their significance, researchers are using a variety of approaches, including transcriptomics, functional genomics, and molecular biology. This review offers a comprehensive overview of the current knowledge of OGs in all domains of life, highlighting the possible role of dark transcriptomics in their evolution. More research is needed to fully comprehend the role of OGs in biology and their impact on various biological processes.

Isolation and Genome-Based Characterization of Biocontrol Potential of Bacillus siamensis YB-1631 against Wheat Crown Rot Caused by Fusarium pseudograminearum

Fusarium crown rot (FCR) caused by Fusarium pseudograminearum is one of the most serious soil-borne diseases of wheat. Among 58 bacterial isolates from the rhizosphere soil of winter wheat seedlings, strain YB-1631 was found to have the highest in vitro antagonism to F. pseudograminearum growth. LB cell-free culture filtrates inhibited mycelial growth and conidia germination of F. pseudograminearum by 84.14% and 92.23%, respectively. The culture filtrate caused distortion and disruption of the cells. Using a face-to-face plate assay, volatile substances produced by YB-1631 inhibited F. pseudograminearum growth by 68.16%. In the greenhouse, YB-1631 reduced the incidence of FCR on wheat seedlings by 84.02% and increased root and shoot fresh weights by 20.94% and 9.63%, respectively. YB-1631 was identified as Bacillus siamensis based on the gyrB sequence and average nucleotide identity of the complete genome. The complete genome was 4,090,312 bp with 4357 genes and 45.92% GC content. In the genome, genes were identified for root colonization, including those for chemotaxis and biofilm production, genes for plant growth promotion, including those for phytohormones and nutrient assimilation, and genes for biocontrol activity, including those for siderophores, extracellular hydrolase, volatiles, nonribosomal peptides, polyketide antibiotics, and elicitors of induced systemic resistance. In vitro production of siderophore, β-1, 3-glucanase, amylase, protease, cellulase, phosphorus solubilization, and indole acetic acid were detected. Bacillus siamensis YB-1631 appears to have significant potential in promoting wheat growth and controlling wheat FCR caused by F. pseudograminearum.

Characterization and Genomic Analysis of the Naphthalene-Degrading Delftia tsuruhatensis ULwDis3 Isolated from Seawater

Strains of the genus Delftia are poorly studied microorganisms. In this work, the complete genome of the naphthalene-degrading Delftia tsuruhatensis strain ULwDis3 isolated from seawater of the Gulf of Finland of the Baltic Sea was assembled. For the first time, genes encoding naphthalene cleavage pathways via salicylate and gentisate were identified in a strain of the genus Delftia. The genes are part of one operon (nag genes). Three open reading frames (ORFs) were found in the genome of D. tsuruhatensis strain ULwDis3 that encode gentisate 1.2-dioxygenase. One of the ORFs is part of the nag operon. The physiological and biochemical characteristics of the strain ULwDis3 when cultured in mineral medium with naphthalene as the sole source of carbon and energy were also studied. It was found that after 22 h of growth, the strain stopped consuming naphthalene, and at the same time, naphthalene 1.2-dioxygenase and salicylate 5-hydroxylase activities were not detected. Later, a decrease in the number of living cells and the death of the culture were observed. Gentisate 1.2-dioxygenase activity was detected from the time of gentisate formation until culture death.

Pathogen-associated gene discovery workflows for novel antivirulence therapeutic development

Novel therapeutics to manage bacterial infections are urgently needed as the impact and prevalence of antimicrobial resistance (AMR) grows. Antivirulence therapeutics are an alternative approach to antibiotics that aim to attenuate virulence rather than target bacterial essential functions, while minimizing microbiota perturbation and the risk of AMR development. Beyond known virulence factors, pathogen-associated genes (PAGs; genes found only in pathogens to date) may play an important role in virulence or host association. Many identified PAGs encode uncharacterized hypothetical proteins and represent an untapped wealth of novel drug targets. Here, we review current advances in antivirulence drug research and development, including PAG identification, and provide a comprehensive workflow from the discovery of antivirulence drug targets to drug discovery. We highlight the importance of integrating bioinformatic/genomic-based methods for novel virulence factor discovery, coupled with experimental characterization, into existing drug screening platforms to develop novel and effective antivirulence drugs.

Connecting genomic islands across prokaryotic and phage genomes via protein families

Prokaryotic genomes evolve via horizontal gene transfer (HGT), mutations, and rearrangements. A noteworthy part of the HGT process is facilitated by genomic islands (GIs). While previous computational biology research has focused on developing tools to detect GIs in prokaryotic genomes, there has been little research investigating GI patterns and biological connections across species. We have pursued the novel idea of connecting GIs across prokaryotic and phage genomes via patterns of protein families. Such patterns are sequences of protein families frequently present in the genomes of multiple species. We combined the large data set from the IslandViewer4 database with protein families from Pfam while implementing a comprehensive strategy to identify patterns making use of HMMER, BLAST, and MUSCLE. we also implemented Python programs that link the analysis into a single pipeline. Research results demonstrated that related GIs often exist in species that are evolutionarily unrelated and in multiple bacterial phyla. Analysis of the discovered patterns led to the identification of biological connections among prokaryotes and phages. These connections suggest broad HGT connections across the bacterial kingdom and its associated phages. The discovered patterns and connections could provide the basis for additional analysis on HGT breadth and the patterns in pathogenic GIs.

Leveraging comparative genomics to uncover alien genes in bacterial genomes

A significant challenge in bacterial genomics is to catalogue genes acquired through the evolutionary process of horizontal gene transfer (HGT). Both comparative genomics and sequence composition-based methods have often been invoked to quantify horizontally acquired genes in bacterial genomes. Comparative genomics methods rely on completely sequenced genomes and therefore the confidence in their predictions increases as the databases become more enriched in completely sequenced genomes. Recent developments including in microbial genome sequencing call for reassessment of alien genes based on information-rich resources currently available. We revisited the comparative genomics approach and developed a new algorithm for alien gene detection. Our algorithm compared favourably with the existing comparative genomics-based methods and is capable of detecting both recent and ancient transfers. It can be used as a standalone tool or in concert with other complementary algorithms for comprehensively cataloguing alien genes in bacterial genomes.

Antibiotic and metal resistance of Stenotrophomonas maltophilia isolates from Eboling permafrost of the Tibetan Plateau

Whole-genome sequencing of pathogenic bacteria Stenotrophomonas maltophilia from a less polluted environment of permafrost can help understand the intrinsic resistome of both antibiotics and metals. This study aimed to examine the maximum minimum inhibitory concentration (MIC) of both antibiotics and metals, as well as antibiotic resistance genes and metal resistance genes annotated from whole-genome sequences. The permafrost S. maltophilia was sensitive to ciprofloxacin, tetracycline, streptomycin, and bacitracin, and resistant to chloramphenicol, trimethoprim-sulfamethoxazole, erythromycin, Zn²⁺, Ni²⁺, Cu²⁺, and Cr⁶⁺, with a lower maximum MIC, compared with clinical S. maltophilia. The former strain belonged to the lower antibiotic resistance gene (ARG) and metal resistance gene (MRG) clusters compared with the latter ones. The permafrost strain contained no or only one kind of ARG or MRG on a single genomic island, which explained the aforementioned lower maximum MIC and less diversity of ARGs or MRGs. The result indicated that the co-occurrence of antibiotic and metal resistance was due to a certain innate ability of S. maltophilia. The continuous human use of antibiotics or metals induced selective pressure, resulting in higher MIC and more diverse ARGs and MRGs in human-impacted environments.

Type IV Pili Are a Critical Virulence Factor in Clinical Isolates of Paenibacillus thiaminolyticus

Hydrocephalus, the leading indication for childhood neurosurgery worldwide, is particularly prevalent in low- and middle-income countries. Hydrocephalus preceded by an infection, or postinfectious hydrocephalus, accounts for up to 60% of hydrocephalus in these areas. Since many children with hydrocephalus suffer poor long-term outcomes despite surgical intervention, prevention of hydrocephalus remains paramount. Our previous studies implicated a novel bacterial pathogen, Paenibacillus thiaminolyticus, as a causal agent of neonatal sepsis and postinfectious hydrocephalus in Uganda. Here, we report the isolation of three P. thiaminolyticus strains, Mbale, Mbale2, and Mbale3, from patients with postinfectious hydrocephalus. We constructed complete genome assemblies of the clinical isolates as well as the nonpathogenic P. thiaminolyticus reference strain and performed comparative genomic and proteomic analyses to identify potential virulence factors. All three isolates carry a unique beta-lactamase gene, and two of the three isolates exhibit resistance in culture to the beta-lactam antibiotics penicillin and ampicillin. In addition, a cluster of genes carried on a mobile genetic element that encodes a putative type IV pilus operon is present in all three clinical isolates but absent in the reference strain. CRISPR-mediated deletion of the gene cluster substantially reduced the virulence of the Mbale strain in mice. Comparative proteogenomic analysis identified various additional potential virulence factors likely acquired on mobile genetic elements in the virulent strains. These results provide insight into the emergence of virulence in P. thiaminolyticus and suggest avenues for the diagnosis and treatment of this novel bacterial pathogen. IMPORTANCE Postinfectious hydrocephalus, a devastating sequela of neonatal infection, is associated with increased childhood mortality and morbidity. A novel bacterial pathogen, Paenibacillus thiaminolyticus, is highly associated with postinfectious hydrocephalus in an African cohort. Whole-genome sequencing, RNA sequencing, and proteomics of clinical isolates and a reference strain in combination with CRISPR editing identified type IV pili as a critical virulence factor for P. thiaminolyticus infection. Acquisition of a type IV pilus-encoding mobile genetic element critically contributed to converting a nonpathogenic strain of P. thiaminolyticus into a pathogen capable of causing devastating diseases. Given the widespread presence of type IV pilus in pathogens, the presence of the type IV pilus operon could serve as a diagnostic and therapeutic target in P. thiaminolyticus and related bacteria.

Reconstructing horizontal gene flow network to understand prokaryotic evolution

Horizontal gene transfer (HGT) is a major source of phenotypic innovation and a mechanism of niche adaptation in prokaryotes. Quantification of HGT is critical to decipher its myriad roles in microbial evolution and adaptation. Advances in genome sequencing and bioinformatics have augmented our ability to understand the microbial world, particularly the direct or indirect influence of HGT on diverse life forms. Methods for detecting HGT can be classified into phylogenetic-based and parametric or composition-based approaches. Here, we exploited the complementary strengths of both the approaches to construct a high confidence horizontal gene flow network. Our network is unique in its ability to detect the transfer of native genes of a genome to genomes from other taxa, thus establishing donor and recipient organisms (taxa), rather than through a post hoc analysis as is the practice with several other approaches. The scale-free horizontal gene flow network presented here provides new insights into modes of transfer for the exchange of genetic information and also illuminates differential gene flow across phyla.

Reconstructing horizontal gene flow network to understand prokaryotic evolution

Horizontal gene transfer (HGT) is a major source of phenotypic innovation and a mechanism of niche adaptation in prokaryotes. Quantification of HGT is critical to decipher its myriad roles in microbial evolution and adaptation. Advances in genome sequencing and bioinformatics have augmented our ability to understand the microbial world, particularly the direct or indirect influence of HGT on diverse life forms. Methods for detecting HGT can be classified into phylogenetic-based and parametric or composition-based approaches. Here, we exploited the complementary strengths of both the approaches to construct a high confidence horizontal gene flow network. Our network is unique in its ability to detect the transfer of native genes of a genome to genomes from other taxa, thus establishing donor and recipient organisms (taxa), rather than through a post hoc analysis as is the practice with several other approaches. The scale-free horizontal gene flow network presented here provides new insights into modes of transfer for the exchange of genetic information and also illuminates differential gene flow across phyla.

High-quality pan-genome of Escherichia coli generated by excluding confounding and highly similar strains reveals an association between unique gene clusters and genomic islands

The pan-genome analysis of bacteria provides detailed insight into the diversity and evolution of a bacterial population. However, the genomes involved in the pan-genome analysis should be checked carefully, as the inclusion of confounding strains would have unfavorable effects on the identification of core genes, and the highly similar strains could bias the results of the pan-genome state (open versus closed). In this study, we found that the inclusion of highly similar strains also affects the results of unique genes in pan-genome analysis, which leads to a significant underestimation of the number of unique genes in the pan-genome. Therefore, these strains should be excluded from pan-genome analysis at the early stage of data processing. Currently, tens of thousands of genomes have been sequenced for Escherichia coli, which provides an unprecedented opportunity as well as a challenge for pan-genome analysis of this classical model organism. Using the proposed strategies, a high-quality E. coli pan-genome was obtained, and the unique genes was extracted and analyzed, revealing an association between the unique gene clusters and genomic islands from a pan-genome perspective, which may facilitate the identification of genomic islands.

Integrative Assessments on Molecular Taxonomy of Acidiferrobacter thiooxydans ZJ and Its Environmental Adaptation Based on Mobile Genetic Elements

Acidiferrobacter spp. are facultatively anaerobic acidophiles that belong to a distinctive Acidiferrobacteraceae family, which are similar to Ectothiorhodospiraceae phylogenetically, and are closely related to Acidithiobacillia class/subdivision physiologically. The limited genome information has kept them from being studied on molecular taxonomy and environmental adaptation in depth. Herein, Af. thiooxydans ZJ was isolated from acid mine drainage (AMD), and the complete genome sequence was reported to scan its genetic constitution for taxonomic and adaptative feature exploration. The genome has a single chromosome of 3,302,271 base pairs (bp), with a GC content of 63.61%. The phylogenetic tree based on OrthoANI highlighted the unique position of Af. thiooxydans ZJ, which harbored more unique genes among the strains from Ectothiorhodospiraceae and Acidithiobacillaceae by pan-genome analysis. The diverse mobile genetic elements (MGEs), such as insertion sequence (IS), clustered regularly interspaced short palindromic repeat (CRISPR), prophage, and genomic island (GI), have been identified and characterized in Af. thiooxydans ZJ. The results showed that Af. thiooxydans ZJ may effectively resist the infection of foreign viruses and gain functional gene fragments or clusters to shape its own genome advantageously. This study will offer more evidence of the genomic plasticity and improve our understanding of evolutionary adaptation mechanisms to extreme AMD environment, which could expand the potential utilization of Af. thiooxydans ZJ as an iron and sulfur oxidizer in industrial bioleaching.

Complete genome sequence and identification of polyunsaturated fatty acid biosynthesis genes of the myxobacterium Minicystis rosea DSM 24000T

Background

Myxobacteria harbor numerous biosynthetic gene clusters that can produce a diverse range of secondary metabolites. Minicystis rosea DSM 24000^T is a soil-dwelling myxobacterium belonging to the suborderSorangiineae and family Polyangiaceae and is known to produce various secondary metabolites as well as polyunsaturated fatty acids (PUFAs). Here, we use whole-genome sequencing to explore the diversity of biosynthetic gene clusters in M. rosea.

Results

Using PacBio sequencing technology, we assembled the 16.04 Mbp complete genome of M. rosea DSM 24000^T, the largest bacterial genome sequenced to date. About 44% of its coding potential represents paralogous genes predominantly associated with signal transduction, transcriptional regulation, and protein folding. These genes are involved in various essential functions such as cellular organization, diverse niche adaptation, and bacterial cooperation, and enable social behavior like gliding motility, sporulation, and predation, typical of myxobacteria. A profusion of eukaryotic-like kinases (353) and an elevated ratio of phosphatases (8.2/1) in M. rosea as compared to other myxobacteria suggest gene duplication as one of the primary modes of genome expansion. About 7.7% of the genes are involved in the biosynthesis of a diverse array of secondary metabolites such as polyketides, terpenes, and bacteriocins. Phylogeny of the genes involved in PUFA biosynthesis (pfa) together with the conserved synteny of the complete pfa gene cluster suggests acquisition via horizontal gene transfer from Actinobacteria.

Conclusion

Overall, this study describes the complete genome sequence of M. rosea, comparative genomic analysis to explore the putative reasons for its large genome size, and explores the secondary metabolite potential, including the biosynthesis of polyunsaturated fatty acids.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07955-x.

Carbapenem Resistance Determinants Acquired through Novel Chromosomal Integrations in Extensively Drug-Resistant Pseudomonas aeruginosa

Two novel bla_DIM-1- or bla_IMP-1-containing genomic islands (GIs) were discovered by whole-genome sequence analyses in four extensively drug-resistant (XDR) Pseudomonas aeruginosa isolates from inpatients at a tertiary hospital in Ghana. The strains were of sequence type 234 (ST234) and formed a phylogenetic clade together with ST111, which is recognized as a global high-risk clone. Their carbapenem resistance was encoded by two Tn402-type integrons, In1592 (bla_DIM-1) and In1595 (bla_IMP-1), both carrying complete tni mobilization modules. In1595 was bound by conserved 25-bp inverted repeats (IRs) flanked by 5-bp direct repeats (DRs) associated with target site duplication. The integrons were embedded in two GIs that contained cognate integrases and were distinguished by a lower GC content than the chromosomal average. PAGI-97A (52.659 bp; In1592), which encoded a P4-type site-specific integrase of the tyrosine recombinase family in its 3′ border, was integrated into tRNA-Pro(ggg) and bracketed by a 49-bp perfect DR created by 3′-end target duplication. GIs with the same structural features, but diverse genetic content, were identified in 41/226 completed P. aeruginosa genomes. PAGI-97B (22,636 bp; In1595), which encoded an XerC/D superfamily integrase in its 5′ border, was inserted into the small RNA (sRNA) PrrF1/PrrF2 locus. Specific insertions into this highly conserved locus involved in iron-dependent regulation, all leaving PrrF1 intact, were identified in an additional six phylogenetically unrelated P. aeruginosa genomes. Our molecular analyses unveiled a hospital-associated clonal dissemination of carbapenem-resistant ST234 P. aeruginosa in which the XDR phenotype resulted from novel insertions of two GIs into specific chromosomal sites.

SSG-LUGIA: Single Sequence based Genome Level Unsupervised Genomic Island Prediction Algorithm

BACKGROUND

Genomic Islands (GIs) are clusters of genes that are mobilized through horizontal gene transfer. GIs play a pivotal role in bacterial evolution as a mechanism of diversification and adaptation to different niches. Therefore, identification and characterization of GIs in bacterial genomes is important for understanding bacterial evolution. However, quantifying GIs is inherently difficult, and the existing methods suffer from low prediction accuracy and precision-recall trade-off. Moreover, several of them are supervised in nature, and thus, their applications to newly sequenced genomes are riddled with their dependency on the functional annotation of existing genomes.

RESULTS

We present SSG-LUGIA, a completely automated and unsupervised approach for identifying GIs and horizontally transferred genes. SSG-LUGIA is a novel method based on unsupervised anomaly detection technique, accompanied by further refinement using cues from signal processing literature. SSG-LUGIA leverages the atypical compositional biases of the alien genes to localize GIs in prokaryotic genomes. SSG-LUGIA was assessed on a large benchmark dataset `IslandPick' and on a set of 15 well-studied genomes in the literature and followed by a thorough analysis on the well-understood Salmonella typhi CT18 genome. Furthermore, the efficacy of SSG-LUGIA in identifying horizontally transferred genes was evaluated on two additional bacterial genomes, namely, those of Corynebacterium diphtheria NCTC13129 and Pseudomonas aeruginosa LESB58. SSG-LUGIA was examined on draft genomes and was demonstrated to be efficient as an ensemble method.

CONCLUSIONS

Our results indicate that SSG-LUGIA achieved superior performance in comparison to frequently used existing methods. Importantly, it yielded a better trade-off between precision and recall than the existing methods. Its nondependency on the functional annotation of genomes makes it suitable for analyzing newly sequenced, yet uncharacterized genomes. Thus, our study is a significant advance in identification of GIs and horizontally transferred genes. SSG-LUGIA is available as an open source software at https://nibtehaz.github.io/SSG-LUGIA/.

Genomic Island Prediction via Chi-Square Test and Random Forest Algorithm

Genomic islands are related to microbial adaptation and carry different genomic characteristics from the host. Therefore, many methods have been proposed to detect genomic islands from the rest of the genome by evaluating its sequence composition. Many sequence features have been proposed, but many of them have not been applied to the identification of genomic islands. In this paper, we present a scheme to predict genomic islands using the chi-square test and random forest algorithm. We extract seven kinds of sequence features and select the important features with the chi-square test. All the selected features are then input into the random forest to predict the genome islands. Three experiments and comparison show that the proposed method achieves the best performance. This understanding can be useful to design more powerful method for the genomic island prediction.

Finding orthologous gene blocks in bacteria: the computational hardness of the problem and novel methods to address it

Motivation

The evolution of complexity is one of the most fascinating and challenging problems in modern biology, and tracing the evolution of complex traits is an open problem. In bacteria, operons and gene blocks provide a model of tractable evolutionary complexity at the genomic level. Gene blocks are structures of co-located genes with related functions, and operons are gene blocks whose genes are co-transcribed on a single mRNA molecule. The genes in operons and gene blocks typically work together in the same system or molecular complex. Previously, we proposed a method that explains the evolution of orthologous gene blocks (orthoblocks) as a combination of a small set of events that take place in vertical evolution from common ancestors. A heuristic method was proposed to solve this problem. However, no study was done to identify the complexity of the problem.

Results

Here, we establish that finding the homologous gene block problem is NP-hard and APX-hard. We have developed a greedy algorithm that runs in polynomial time and guarantees an O(ln⁡n) approximation. In addition, we formalize our problem as an integer linear program problem and solve it using the PuLP package and the standard CPLEX algorithm. Our exploration of several candidate operons reveals that our new method provides more optimal results than the results from the heuristic approach, and is significantly faster.

Availability and implementation

The software and data accompanying this paper are available under the GPLv3 and CC0 license respectively on: https://github.com/nguyenngochuy91/Relevant-Operon.

A hydrogen-oxidizing bacterium enriched from the open ocean resembling a symbiont

A new autotrophic hydrogen-oxidizing Chromatiaceae bacterium, namely bacterium CTD079, was enriched from a water column sample at 1500 m water depth in the southern Pacific Ocean. Based on the phylogeny of 16S rRNA genes, it was closely related to a scaly snail endosymbiont (99.2% DNA sequence identity) whose host so far is only known to colonize hydrothermal vents along the Indian ridge. The average nucleotide identity between the genomes of CTD079 and the snail endosymbiont was 91%. The observed differences likely reflect adaptations to their specific habitats. For example, CTD079 encodes additional enzymes like the formate dehydrogenase increasing the organism's spectrum of energy generation pathways. Other additional physiological features of CTD079 included the increase of viral defence strategies, secretion systems and specific transporters for essential elements. These important genome characteristics suggest an adaptation to life in the open ocean.

Potential evolutionary impact of integrative and conjugative elements (ICEs) and genomic islands in the Ralstonia solanacearum species complex

Ralstonia solanacearum, a soil-borne plant pathogen, encompasses a large number of strains known as R. solanacearum species complex (RSSC). Although it has been suggested that mobile genetic elements (MGEs) may play an important role in the RSSC genome, the evolutionary impact of these elements remains unknown. Here, we identified and analysed Integrative and Conjugative Elements (ICEs) and Genomic Islands (GIs) in the 121 genomes published for Ralstonia spp., including RSSC strains and three other non-plant pathogenic Ralstonia spp. Our results provided a dataset of 12 ICEs and 31 GIs distributed throughout Ralstonia spp. Four novel ICEs in RSSC were found. Some of these elements cover 5% of the host genome and carry accessory genes with a potential impact on the fitness and pathogenicity of RSSC. In addition, phylogenetic analysis revealed that these MGEs clustered to the same species, but there is evidence of strains from different countries that host the same element. Our results provide novel insight into the RSSC adaptation, opening new paths to a better understanding of how these elements affect this soil-borne plant pathogen.

Microbial genomic island discovery, visualization and analysis

Horizontal gene transfer (also called lateral gene transfer) is a major mechanism for microbial genome evolution, enabling rapid adaptation and survival in specific niches. Genomic islands (GIs), commonly defined as clusters of bacterial or archaeal genes of probable horizontal origin, are of particular medical, environmental and/or industrial interest, as they disproportionately encode virulence factors and some antimicrobial resistance genes and may harbor entire metabolic pathways that confer a specific adaptation (solvent resistance, symbiosis properties, etc). As large-scale analyses of microbial genomes increases, such as for genomic epidemiology investigations of infectious disease outbreaks in public health, there is increased appreciation of the need to accurately predict and track GIs. Over the past decade, numerous computational tools have been developed to tackle the challenges inherent in accurate GI prediction. We review here the main types of GI prediction methods and discuss their advantages and limitations for a routine analysis of microbial genomes in this era of rapid whole-genome sequencing. An assessment is provided of 20 GI prediction software methods that use sequence-composition bias to identify the GIs, using a reference GI data set from 104 genomes obtained using an independent comparative genomics approach. Finally, we present guidelines to assist researchers in effectively identifying these key genomic regions.

Genome-Wide Analysis Reveals Genetic Potential for Aromatic Compounds Biodegradation of Sphingopyxis

Members of genus Sphingopyxis are frequently found in diverse eco-environments worldwide and have been traditionally considered to play vital roles in the degradation of aromatic compounds. Over recent decades, many aromatic-degrading Sphingopyxis strains have been isolated and recorded, but little is known about their genetic nature related to aromatic compounds biodegradation. In this study, bacterial genomes of 19 Sphingopyxis strains were used for comparative analyses. Phylogeny showed an ambiguous relatedness between bacterial strains and their habitat specificity, while clustering based on Cluster of Orthologous Groups suggested the potential link of functional profile with substrate-specific traits. Pan-genome analysis revealed that 19 individuals were predicted to share 1,066 orthologous genes, indicating a high genetic homogeneity among Sphingopyxis strains. Notably, KEGG Automatic Annotation Server results suggested that most genes pertaining aromatic compounds biodegradation were predicted to be involved in benzoate, phenylalanine, and aminobenzoate metabolism. Among them, β-ketoadipate biodegradation might be the main pathway in Sphingopyxis strains. Further inspection showed that a number of mobile genetic elements varied in Sphingopyxis genomes, and plasmid-mediated gene transfer coupled with prophage- and transposon-mediated rearrangements might play prominent roles in the evolution of bacterial genomes. Collectively, our findings presented that Sphingopyxis isolates might be the promising candidates for biodegradation of aromatic compounds in pollution sites.

Tracing the ancestry of operons in bacteria

MOTIVATION

Complexity is a fundamental attribute of life. Complex systems are made of parts that together perform functions that a single component, or subsets of components, cannot. Examples of complex molecular systems include protein structures such as the F1Fo-ATPase, the ribosome, or the flagellar motor: each one of these structures requires most or all of its components to function properly. Given the ubiquity of complex systems in the biosphere, understanding the evolution of complexity is central to biology. At the molecular level, operons are classic examples of a complex system. An operon's genes are co-transcribed under the control of a single promoter to a polycistronic mRNA molecule, and the operon's gene products often form molecular complexes or metabolic pathways. With the large number of complete bacterial genomes available, we now have the opportunity to explore the evolution of these complex entities, by identifying possible intermediate states of operons.

RESULTS

In this work, we developed a maximum parsimony algorithm to reconstruct ancestral operon states, and show a simple vertical evolution model of how operons may evolve from the individual component genes. We describe several ancestral states that are plausible functional intermediate forms leading to the full operon. We also offer Reconstruction of Ancestral Gene blocks Using Events or ROAGUE as a software tool for those interested in exploring gene block and operon evolution.

AVAILABILITY AND IMPLEMENTATION

The software accompanying this paper is available under GPLv3 license on: https://github.com/nguyenngochuy91/Ancestral-Blocks-Reconstruction.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

2SigFinder: the combined use of small-scale and large-scale statistical testing for genomic island detection from a single genome

BACKGROUND

Genomic islands are associated with microbial adaptations, carrying genomic signatures different from the host. Some methods perform an overall test to identify genomic islands based on their local features. However, regions of different scales will display different genomic features.

RESULTS

We proposed here a novel method "2SigFinder ", the first combined use of small-scale and large-scale statistical testing for genomic island detection. The proposed method was tested by genomic island boundary detection and identification of genomic islands or functional features of real biological data. We also compared the proposed method with the comparative genomics and composition-based approaches. The results indicate that the proposed 2SigFinder is more efficient in identifying genomic islands.

CONCLUSIONS

From real biological data, 2SigFinder identified genomic islands from a single genome and reported robust results across different experiments, without annotated information of genomes or prior knowledge from other datasets. 2SigHunter identified 25 Pathogenicity, 1 tRNA, 2 Virulence and 2 Repeats from 27 Pathogenicity, 1 tRNA, 2 Virulence and 2 Repeats, and detected 101 Phage and 28 HEG out of 130 Phage and 36 HEGs in S. enterica Typhi CT18, which shows that it is more efficient in detecting functional features associated with GIs.

Analysis of the biodegradative and adaptive potential of the novel polychlorinated biphenyl degrader Rhodococcus sp. WAY2 revealed by its complete genome sequence

The complete genome sequence of Rhodococcus sp. WAY2 (WAY2) consists of a circular chromosome, three linear replicons and a small circular plasmid. The linear replicons contain typical actinobacterial invertron-type telomeres with the central CGTXCGC motif. Comparative phylogenetic analysis of the 16S rRNA gene along with phylogenomic analysis based on the genome-to-genome blast distance phylogeny (GBDP) algorithm and digital DNA–DNA hybridization (dDDH) with other Rhodococcus type strains resulted in a clear differentiation of WAY2, which is likely a new species. The genome of WAY2 contains five distinct clusters of bph, etb and nah genes, putatively involved in the degradation of several aromatic compounds. These clusters are distributed throughout the linear plasmids. The high sequence homology of the ring-hydroxylating subunits of these systems with other known enzymes has allowed us to model the range of aromatic substrates they could degrade. Further functional characterization revealed that WAY2 was able to grow with biphenyl, naphthalene and xylene as sole carbon and energy sources, and could oxidize multiple aromatic compounds, including ethylbenzene, phenanthrene, dibenzofuran and toluene. In addition, WAY2 was able to co-metabolize 23 polychlorinated biphenyl congeners, consistent with the five different ring-hydroxylating systems encoded by its genome. WAY2 could also use n-alkanes of various chain-lengths as a sole carbon source, probably due to the presence of alkB and ladA gene copies, which are only found in its chromosome. These results show that WAY2 has a potential to be used for the biodegradation of multiple organic compounds.

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.

VCGIDB: A Database and Web Resource for the Genomic Islands from Vibrio Cholerae

Vibrio cholerae is the causative agent of cholera, which is a severe, life-threatening diarrheal disease. The current seventh pandemic has not been eradicated and the outbreak is still ongoing around the world. The evolution of the pandemic-causing strain has been greatly influenced by lateral gene transfer, and the mechanisms of acquisition of pathogenicity in V. cholerae are mainly involved with genomic islands (GIs). Thus, detecting GIs and their comprehensive information is necessary to understand the continuing resurgence and newly emerging pathogenic V. cholerae strains. In this study, 798 V. cholerae strains were tested using the GI-Scanner algorithm, which was developed to detect candidate GIs and identify them in a comparative genomics approach. The algorithm predicted 435 highly possible genomic islands, and we built a database, called Vibrio cholerae Genomic Island Database (VCGIDB). This database shows advanced results that were acquired from a large genome set using phylogeny-based predictions. Moreover, VCGIDB is a highly expendable database that does not require intensive computation, which enables us to update it with a greater number of genomes using a novel genomic island prediction method. The VCGIDB website allows the user to browse the data and presents the results in a visual manner.

IslandCafe: Compositional Anomaly and Feature Enrichment Assessment for Delineation of Genomic Islands

One of the evolutionary forces driving bacterial genome evolution is the acquisition of clusters of genes through horizontal gene transfer (HGT). These genomic islands may confer adaptive advantages to the recipient bacteria, such as, the ability to thwart antibiotics, become virulent or hypervirulent, or acquire novel metabolic traits. Methods for detecting genomic islands either search for markers or features typical of islands or examine anomaly in oligonucleotide composition against the genome background. The former tends to underestimate, missing islands that have the markers either lost or degraded, while the latter tends to overestimate, due to their inability to discriminate compositional atypicality arising because of HGT from those that are a consequence of other biological factors. We propose here a framework that exploits the strengths of both these approaches while bypassing the pitfalls of either. Genomic islands lacking markers are identified by their association with genomic islands with markers. This was made possible by performing marker enrichment and phyletic pattern analyses within an integrated framework of recursive segmentation and clustering. The proposed method, IslandCafe, compared favorably with frequently used methods for genomic island detection on synthetic test datasets and on a test-set of known islands from 15 well-characterized bacterial species. Furthermore, IslandCafe identified novel islands with imprints of likely horizontal acquisition.

Identification of an Integrase That Responsible for Precise Integration and Excision of Riemerella anatipestifer Genomic Island

Riemerella anatipestifer is a Gram-negative, pathogenic bacterium, which is harmful to poultry. However, the genomic islands (GIs) in R. anatipestifer are not well-studied. In this study, a 10K genomic island was predicted by the bioinformatics analysis of R. anatipestifer ATCC 11845, which is widely found in other R. anatipestifer genomes. We had first reported the genomic island integration and excision function in R. anatipestifer. We successfully constructed the integration plasmid by using the integrase and 53 bp attP elements. The 10K GI was found integrated at the 53 bp attB located in the Arg-tRNA of the R. anatipestifer RA-YM chromosome. We identified an integrase that helped in the precise integration and excision in R. anatipestifer and elucidated the molecular mechanism of the 10K genomic island integration and excision. Furthermore, we provided a new method for the gene expression and construction of complementary strain.

Improved genomic island predictions with IslandPath-DIMOB

Motivation

Genomic islands (GIs) are clusters of genes of probable horizontal origin that play a major role in bacterial and archaeal genome evolution and microbial adaptability. They are of high medical and industrial interest, due to their enrichment in virulence factors, some antimicrobial resistance genes and adaptive metabolic pathways. The development of more sensitive but precise prediction tools, using either sequence composition-based methods or comparative genomics, is needed as large-scale analyses of microbial genomes increase.

Results

IslandPath-DIMOB, a leading GI prediction tool in the IslandViewer webserver, has now been significantly improved by modifying both the decision algorithm to determine sequence composition biases, and the underlying database of HMM profiles for associated mobility genes. The accuracy of IslandPath-DIMOB and other major software has been assessed using a reference GI dataset predicted by comparative genomics, plus a manually curated dataset from literature review. Compared to the previous version (v0.2.0), this IslandPath-DIMOB v1.0.0 achieves 11.7% and 5.3% increase in recall and precision, respectively. IslandPath-DIMOB has the highest Matthews correlation coefficient among individual prediction methods tested, combining one of the highest recall measures (46.9%) at high precision (87.4%). The only method with higher recall had notably lower precision (55.1%). This new IslandPath-DIMOB v1.0.0 will facilitate more accurate studies of GIs, including their key roles in microbial adaptability of medical, environmental and industrial interest.

Availability and implementation

IslandPath-DIMOB v1.0.0 is freely available through the IslandViewer webserver {{http://www.pathogenomics.sfu.ca/islandviewer/}} and as standalone software {{https://github.com/brinkmanlab/islandpath/}} under the GNU-GPLv3.

Supplementary information

Supplementary data are available at Bioinformatics online.

IslandViewer 4: expanded prediction of genomic islands for larger-scale datasets

IslandViewer (http://www.pathogenomics.sfu.ca/islandviewer/) is a widely-used webserver for the prediction and interactive visualization of genomic islands (GIs, regions of probable horizontal origin) in bacterial and archaeal genomes. GIs disproportionately encode factors that enhance the adaptability and competitiveness of the microbe within a niche, including virulence factors and other medically or environmentally important adaptations. We report here the release of IslandViewer 4, with novel features to accommodate the needs of larger-scale microbial genomics analysis, while expanding GI predictions and improving its flexible visualization interface. A user management web interface as well as an HTTP API for batch analyses are now provided with a secured authentication to facilitate the submission of larger numbers of genomes and the retrieval of results. In addition, IslandViewer's integrated GI predictions from multiple methods have been improved and expanded by integrating the precise Islander method for pre-computed genomes, as well as an updated IslandPath-DIMOB for both pre-computed and user-supplied custom genome analysis. Finally, pre-computed predictions including virulence factors and antimicrobial resistance are now available for 6193 complete bacterial and archaeal strains publicly available in RefSeq. IslandViewer 4 provides key enhancements to facilitate the analysis of GIs and better understand their role in the evolution of successful environmental microbes and pathogens.

Horizontal acquisition of hydrogen conversion ability and other habitat adaptations in the Hydrogenovibrio strains SP-41 and XCL-2

BACKGROUND

Obligate sulfur oxidizing chemolithoauthotrophic strains of Hydrogenovibrio crunogenus have been isolated from multiple hydrothermal vent associated habitats. However, a hydrogenase gene cluster (encoding the hydrogen converting enzyme and its maturation/assembly machinery) detected on the first sequenced H. crunogenus strain (XCL-2) suggested that hydrogen conversion may also play a role in this organism. Yet, numerous experiments have underlined XCL-2's inability to consume hydrogen under the tested conditions. A recent study showed that the closely related strain SP-41 contains a homolog of the XCL-2 hydrogenase (a group 1b [NiFe]-hydrogenase), but that it can indeed use hydrogen. Hence, the question remained unresolved, why SP-41 is capable of using hydrogen, while XCL-2 is not.

RESULTS

Here, we present the genome sequence of the SP-41 strain and compare it to that of the XCL-2 strain. We show that the chromosome of SP-41 codes for a further hydrogenase gene cluster, including two additional hydrogenases: the first appears to be a group 1d periplasmic membrane-anchored hydrogenase, and the second a group 2b sensory hydrogenase. The region where these genes are located was likely acquired horizontally and exhibits similarity to other Hydrogenovibrio species (H. thermophilus MA2-6 and H. marinus MH-110 T) and other hydrogen oxidizing Proteobacteria (Cupriavidus necator H16 and Ghiorsea bivora TAG-1 T). The genomes of XCL-2 and SP-41 show a strong conservation in gene order. However, several short genomic regions are not contained in the genome of the other strain. These exclusive regions are often associated with signs of DNA mobility, such as genes coding for transposases. They code for transport systems and/or extend the metabolic potential of the strains.

CONCLUSIONS

Our results suggest that horizontal gene transfer plays an important role in shaping the genomes of these strains, as a likely mechanism for habitat adaptation, including, but not limited to the transfer of the hydrogen conversion ability.

MTGIpick allows robust identification of genomic islands from a single genome

Genomic islands (GIs) that are associated with microbial adaptations and carry sequence patterns different from that of the host are sporadically distributed among closely related species. This bias can dominate the signal of interest in GI detection. However, variations still exist among the segments of the host, although no uniform standard exists regarding the best methods of discriminating GIs from the rest of the genome in terms of compositional bias. In the present work, we proposed a robust software, MTGIpick, which used regions with pattern bias showing multiscale difference levels to identify GIs from the host. MTGIpick can identify GIs from a single genome without annotated information of genomes or prior knowledge from other data sets. When real biological data were used, MTGIpick demonstrated better performance than existing methods, as well as revealed potential GIs with accurate sizes missed by existing methods because of a uniform standard. Software and supplementary are freely available at http://bioinfo.zstu.edu.cn/MTGI or https://github.com/bioinfo0706/MTGIpick.

Nucleotide Second Messenger-Based Signaling in Extreme Acidophiles of the Acidithiobacillus Species Complex: Partition Between the Core and Variable Gene Complements

Cyclic and linear nucleotides are key elements of the signal transduction networks linking perception of the environment to specific cellular behavior of prokaryotes. These molecular mechanisms are particularly important in bacteria exposed to different, and frequently simultaneous, types of extreme conditions. This is the case in acidithiobacilli, a group of extremophilic bacteria thriving in highly acidic biotopes, that must also cope with significant variations in temperature, osmotic potentials and concentrations of various transition metals and metalloids. Environmental cues sensed by bacteria are transduced into differential levels of nucleotides acting as intracellular second messengers, promoting the activation or inhibition of target components and eliciting different output phenotypes. Cyclic (c) di-GMP, one of the most common bacterial second messengers, plays a key role in lifestyle changes in many bacteria, including acidithiobacilli. The presence of functional c-di-GMP-dependent signal transduction pathways in representative strains of the best-known linages of this species complex has been reported. However, a comprehensive panorama of the c-di-GMP modulated networks, the cognate input signals and output responses, are still missing for this group of extremophiles. Moreover, little fundamental understanding has been gathered for other nucleotides acting as second messengers. Taking advantage of the increasing number of sequenced genomes of the taxon, here we address the challenge of disentangling the nucleotide-driven signal transduction pathways in this group of polyextremophiles using comparative genomic tools and strategies. Results indicate that the acidithiobacilli possess all the genetic elements required to establish functional transduction pathways based in three different nucleotide-second messengers: (p)ppGpp, cyclic AMP (cAMP), and c-di-GMP. The elements related with the metabolism and transduction of (p)ppGpp and cAMP appear highly conserved, integrating signals related with nutrient starvation and polyphosphate metabolism, respectively. In contrast, c-di-GMP networks appear diverse and complex, differing both at the species and strain levels. Molecular elements of c-di-GMP metabolism and transduction were mostly found scattered along the flexible genome of the acidithiobacilli, allowing the identification of probable control modules that could be critical for substrate colonization, biofilm development and intercellular interactions. These may ultimately convey increased endurance to environmental stress and increased potential for gene sharing and adaptation to changing conditions.

An efficient genomic signature ranking method for genomic island prediction from a single genome

Genomic islands that are associated with microbial adaptations and carry genomic signatures different from that of the host, and thus many methods have been proposed to select the informative genomic signatures from a range of organisms and discriminate genomic islands from the rest of the genome in terms of these signature biases. However, they are of limited use when closely related genomes are unavailable. In the present work, we proposed a kurtosis-based ranking method to select the informative genomic signatures from a single genome. In simulations with alien fragments from artificial and real genomes, the proposed kurtosis-based ranking method efficiently selected the informative genomic signatures from a single genome, without annotated information of genomes or prior knowledge from other datasets. This understanding can be useful to design more powerful method for genomic island detection.

Orphan Genes Shared by Pathogenic Genomes Are More Associated with Bacterial Pathogenicity

Recent pangenome analyses of numerous bacterial species have suggested that each genome of a single species may have a significant fraction of its gene content unique or shared by a very few genomes (i.e., ORFans). We selected nine bacterial genera, each containing at least five pathogenic and five nonpathogenic genomes, to compare their ORFans in relation to pathogenicity-related genes. Pathogens in these genera are known to cause a number of common and devastating human diseases such as pneumonia, diphtheria, melioidosis, and tuberculosis. Thus, they are worthy of in-depth systems microbiology investigations, including the comparative study of ORFans between pathogens and nonpathogens. We provide direct evidence to suggest that ORFans shared by more pathogens are more associated with pathogenicity-related genes and thus are more important targets for development of new diagnostic markers or therapeutic drugs for bacterial infectious diseases.

Orphan genes (also known as ORFans [i.e., orphan open reading frames]) are new genes that enable an organism to adapt to its specific living environment. Our focus in this study is to compare ORFans between pathogens (P) and nonpathogens (NP) of the same genus. Using the pangenome idea, we have identified 130,169 ORFans in nine bacterial genera (505 genomes) and classified these ORFans into four groups: (i) SS-ORFans (P), which are only found in a single pathogenic genome; (ii) SS-ORFans (NP), which are only found in a single nonpathogenic genome; (iii) PS-ORFans (P), which are found in multiple pathogenic genomes; and (iv) NS-ORFans (NP), which are found in multiple nonpathogenic genomes. Within the same genus, pathogens do not always have more genes, more ORFans, or more pathogenicity-related genes (PRGs)—including prophages, pathogenicity islands (PAIs), virulence factors (VFs), and horizontal gene transfers (HGTs)—than nonpathogens. Interestingly, in pathogens of the nine genera, the percentages of PS-ORFans are consistently higher than those of SS-ORFans, which is not true in nonpathogens. Similarly, in pathogens of the nine genera, the percentages of PS-ORFans matching the four types of PRGs are also always higher than those of SS-ORFans, but this is not true in nonpathogens. All of these findings suggest the greater importance of PS-ORFans for bacterial pathogenicity.

IMPORTANCE Recent pangenome analyses of numerous bacterial species have suggested that each genome of a single species may have a significant fraction of its gene content unique or shared by a very few genomes (i.e., ORFans). We selected nine bacterial genera, each containing at least five pathogenic and five nonpathogenic genomes, to compare their ORFans in relation to pathogenicity-related genes. Pathogens in these genera are known to cause a number of common and devastating human diseases such as pneumonia, diphtheria, melioidosis, and tuberculosis. Thus, they are worthy of in-depth systems microbiology investigations, including the comparative study of ORFans between pathogens and nonpathogens. We provide direct evidence to suggest that ORFans shared by more pathogens are more associated with pathogenicity-related genes and thus are more important targets for development of new diagnostic markers or therapeutic drugs for bacterial infectious diseases.

Genome Sequence and Metabolic Analysis of a Fluoranthene-Degrading Strain Pseudomonas aeruginosa DN1

Pseudomonas aeruginosa DN1, isolated from petroleum-contaminated soil, showed excellent degradation ability toward diverse polycyclic aromatic hydrocarbons (PAHs). Many studies have been done to improve its degradation ability. However, the molecular mechanisms of PAHs degradation in DN1 strain are unclear. In this study, the whole genome of DN1 strain was sequenced and analyzed. Its genome contains 6,641,902 bp and encodes 6,684 putative open reading frames (ORFs), which has the largest genome in almost all the comparative Pseudomonas strains. Results of gene annotation showed that this strain harbored over 100 candidate genes involved in PAHs degradation, including those encoding 25 dioxygenases, four ring-hydroxylating dioxygenases, five ring-cleaving dioxygenases, and various catabolic enzymes, transcriptional regulators, and transporters in the degradation pathways. In addition, gene knockout experiments revealed that the disruption of some key PAHs degradation genes in DN1 strain, such as catA, pcaG, pcaH, and rhdA, did not completely inhibit fluoranthene degradation, even though their degradative rate reduced to some extent. Three intermediate metabolites, including 9-hydroxyfluorene, 1-acenaphthenone, and 1, 8-naphthalic anhydride, were identified as the dominating intermediates in presence of 50 μg/mL fluoranthene as the sole carbon source according to gas chromatography mass spectrometry analysis. Taken together, the genomic and metabolic analysis indicated that the fluoranthene degradation by DN1 strain was initiated by dioxygenation at the C-1, 2-, C-2, 3-, and C-7, 8- positions. These results provide new insights into the genomic plasticity and environmental adaptation of DN1 strain.

Complete Genomic Analysis of a Kingdom-Crossing Klebsiella variicola Isolate

Bacterial isolate X39 was isolated from a community-acquired pneumonia patient in Beijing, China. A phylogenetic tree based on rpoB genes and average nucleotide identity data confirmed that isolate X39 belonged to Klebsiella variicola. The genome of K. variicola X39 contained one circular chromosome and nine plasmids. Comparative genomic analyses with other K. variicola isolates revealed that K. variicola X39 contained the most unique genes. Of these unique genes, many were prophages and transposases. Many virulence factors were shared between K. variicola X39 and Klebsiella pneumoniae F1. The pathogenicity of K. variicola X39 was compared with that of K. pneumoniae F1 in an abdominal infection model. The results indicated that K. variicola X39 was less virulent than typical clinical K. pneumoniae F1. The genome of K. variicola X39 also contained some genes involved in plant colonization, nitrogen fixation, and defense against oxidative stress. GFP-labeled K. variicola X39 could colonize maize as an endophytic bacterium. We concluded that K. variicola X39 was a kingdom-crossing strain.

Integrase-Controlled Excision of Metal-Resistance Genomic Islands in Acinetobacter baumannii

Genomic islands (GIs) are discrete gene clusters encoding for a variety of functions including antibiotic and heavy metal resistance, some of which are tightly associated to lineages of the core genome phylogenetic tree. We have investigated the functions of two distinct integrase genes in the mobilization of two metal resistant GIs, G08 and G62, of Acinetobacter baumannii. Real-time PCR demonstrated integrase-dependent GI excision, utilizing isopropyl β-d-1-thiogalactopyranoside IPTG-inducible integrase genes in plasmid-based mini-GIs in Escherichia coli. In A. baumannii, integrase-dependent excision of the original chromosomal GIs could be observed after mitomycin C induction. In both E. coli plasmids and A. baumannii chromosome, the rate of excision and circularization was found to be dependent on the expression level of the integrases. Susceptibility testing in A. baumannii strain ATCC 17978, A424, and their respective ΔG62 and ΔG08 mutants confirmed the contribution of the GI-encoded efflux transporters to heavy metal decreased susceptibility. In summary, the data evidenced the functionality of two integrases in the excision and circularization of the two Acinetobacter heavy-metal resistance GIs, G08 and G62, in E. coli, as well as when chromosomally located in their natural host. These recombination events occur at different frequencies resulting in genome plasticity and may participate in the spread of resistance determinants in A. baumannii.

Genome Variation in the Model Halophilic Bacterium Salinibacter ruber

The halophilic bacterium Salinibacter ruber is an abundant and ecologically important member of halophilic communities worldwide. Given its broad distribution and high intraspecific genetic diversity, S. ruber is considered one of the main models for ecological and evolutionary studies of bacterial adaptation to hypersaline environments. However, current insights on the genomic diversity of this species is limited to the comparison of the genomes of two co-isolated strains. Here, we present a comparative genomic analysis of eight S. ruber strains isolated at two different time points in each of two different Mediterranean solar salterns. Our results show an open pangenome with contrasting evolutionary patterns in the core and accessory genomes. We found that the core genome is shaped by extensive homologous recombination (HR), which results in limited sequence variation within population clusters. In contrast, the accessory genome is modulated by horizontal gene transfer (HGT), with genomic islands and plasmids acting as gateways to the rest of the genome. In addition, both types of genetic exchange are modulated by restriction and modification (RM) or CRISPR-Cas systems. Finally, genes differentially impacted by such processes reveal functional processes potentially relevant for environmental interactions and adaptation to extremophilic conditions. Altogether, our results support scenarios that conciliate “Neutral” and “Constant Diversity” models of bacterial evolution.

Genus-wide comparison of Pseudovibrio bacterial genomes reveal diverse adaptations to different marine invertebrate hosts

Bacteria belonging to the genus Pseudovibrio have been frequently found in association with a wide variety of marine eukaryotic invertebrate hosts, indicative of their versatile and symbiotic lifestyle. A recent comparison of the sponge-associated Pseudovibrio genomes has shed light on the mechanisms influencing a successful symbiotic association with sponges. In contrast, the genomic architecture of Pseudovibrio bacteria associated with other marine hosts has received less attention. Here, we performed genus-wide comparative analyses of 18 Pseudovibrio isolated from sponges, coral, tunicates, flatworm, and seawater. The analyses revealed a certain degree of commonality among the majority of sponge- and coral-associated bacteria. Isolates from other marine invertebrate host, tunicates, exhibited a genetic repertoire for cold adaptation and specific metabolic abilities including mucin degradation in the Antarctic tunicate-associated bacterium Pseudovibrio sp. Tun.PHSC04_5.I4. Reductive genome evolution was simultaneously detected in the flatworm-associated bacteria and the sponge-associated bacterium P. axinellae AD2, through the loss of major secretion systems (type III/VI) and virulence/symbioses factors such as proteins involved in adhesion and attachment to the host. Our study also unraveled the presence of a CRISPR-Cas system in P. stylochi UST20140214-052 a flatworm-associated bacterium possibly suggesting the role of CRISPR-based adaptive immune system against the invading virus particles. Detection of mobile elements and genomic islands (GIs) in all bacterial members highlighted the role of horizontal gene transfer for the acquisition of novel genetic features, likely enhancing the bacterial ecological fitness. These findings are insightful to understand the role of genome diversity in Pseudovibrio as an evolutionary strategy to increase their colonizing success across a wide range of marine eukaryotic hosts.

xenoGI: reconstructing the history of genomic island insertions in clades of closely related bacteria

Background

Genomic islands play an important role in microbial genome evolution, providing a mechanism for strains to adapt to new ecological conditions. A variety of computational methods, both genome-composition based and comparative, have been developed to identify them. Some of these methods are explicitly designed to work in single strains, while others make use of multiple strains. In general, existing methods do not identify islands in the context of the phylogeny in which they evolved. Even multiple strain approaches are best suited to identifying genomic islands that are present in one strain but absent in others. They do not automatically recognize islands which are shared between some strains in the clade or determine the branch on which these islands inserted within the phylogenetic tree.

Results

We have developed a software package, xenoGI, that identifies genomic islands and maps their origin within a clade of closely related bacteria, determining which branch they inserted on. It takes as input a set of sequenced genomes and a tree specifying their phylogenetic relationships. Making heavy use of synteny information, the package builds gene families in a species-tree-aware way, and then attempts to combine into islands those families whose members are adjacent and whose most recent common ancestor is shared. The package provides a variety of text-based analysis functions, as well as the ability to export genomic islands into formats suitable for viewing in a genome browser. We demonstrate the capabilities of the package with several examples from enteric bacteria, including an examination of the evolution of the acid fitness island in the genus Escherichia. In addition we use output from simulations and a set of known genomic islands from the literature to show that xenoGI can accurately identify genomic islands and place them on a phylogenetic tree.

Conclusions

xenoGI is an effective tool for studying the history of genomic island insertions in a clade of microbes. It identifies genomic islands, and determines which branch they inserted on within the phylogenetic tree for the clade. Such information is valuable because it helps us understand the adaptive path that has produced living species.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2038-0) contains supplementary material, which is available to authorized users.

Genomic Analysis of Clavibacter michiganensis Reveals Insight Into Virulence Strategies and Genetic Diversity of a Gram-Positive Bacterial Pathogen

Clavibacter michiganensis subsp. michiganensis is a gram-positive bacterial pathogen that proliferates in the xylem vessels of tomato, causing bacterial canker disease. In this study, we sequenced and assembled genomes of 11 C. michiganensis subsp. michiganensis strains isolated from infected tomato fields in California as well as five Clavibacter strains that colonize tomato endophytically but are not pathogenic in this host. The analysis of the C. michiganensis subsp. michiganensis genomes supported the monophyletic nature of this pathogen but revealed genetic diversity among strains, consistent with multiple introduction events. Two tomato endophytes that clustered phylogenetically with C. michiganensis strains capable of infecting wheat and pepper and were also able to cause disease in these plants. Plasmid profiles of the California strains were variable and supported the essential role of the pCM1-like plasmid and the CelA cellulase in virulence, whereas the absence of the pCM2-like plasmid in some pathogenic C. michiganensis subsp. michiganensis strains revealed it is not essential. A large number of secreted C. michiganensis subsp. michiganensis proteins were carbohydrate-active enzymes (CAZymes). Glycome profiling revealed that C. michiganensis subsp. michiganensis but not endophytic Clavibacter strains is able to extensively alter tomato cell-wall composition. Two secreted CAZymes found in all C. michiganensis subsp. michiganensis strains, CelA and PelA1, enhanced pathogenicity on tomato. Collectively, these results provide a deeper understanding of C. michiganensis subsp. michiganensis diversity and virulence strategies.