TnpPred: A Web Service for the Robust Prediction of Prokaryotic Transposases

MOBHunter: a data integration platform for identification and classification of mobile genetic elements in microbial genomes

Horizontal gene transfer plays a critical role in microbial genome evolution and adaptation. Integrated foreign DNA fragments encompass various types of mobile genetic elements (MGEs), ranging from small transposons to conspicuous integrative and conjugative elements. These regions often confer advantageous traits, including antibiotic resistance or novel metabolic capabilities, and contain foreign sequence signatures and hallmark genes such as transposases, integrases, etc. While bioinformatic tools target specific MGE subsets using alignments, compositional signatures, or diagnostic gene mapping, no single platform offers a unified framework for comprehensive, evidence-based, MGE identification and classification. To address this challenge, we developed MOBHunter, an advanced bioinformatic pipeline that consolidates standalone tools and in-house algorithms. Unlike basic tool concatenation, MOBHunter yields consensus identifications, score-supported classifications, and enhanced web visualizations. The platform reduces end user analysis time by integrating data collection, processing, and interpretation into a unified workflow. It delivers robust classifications of MGEs into distinct families and provides a comprehensive overview of the flexible regions of any given input genome. URL: https://informatica.utem.cl/mobhunter/.

DRAMMA: a multifaceted machine learning approach for novel antimicrobial resistance gene detection in metagenomic data

BACKGROUND

Antibiotics are essential for medical procedures, food security, and public health. However, ill-advised usage leads to increased pathogen resistance to antimicrobial substances, posing a threat of fatal infections and limiting the benefits of antibiotics. Therefore, early detection of antimicrobial resistance genes (ARGs), especially in pathogens, is crucial for human health. Most computational methods for ARG detection rely on homology to a predefined gene database and therefore are limited in their ability to discover novel genes.

RESULTS

We introduce DRAMMA, a machine learning method for predicting new ARGs with no sequence similarity to known ARGs or any annotated gene. DRAMMA utilizes various features, including protein properties, genomic context, and evolutionary patterns. The model demonstrated robust predictive performance both in cross-validation and an external validation set annotated by an empirical ARG database. Analyses of the high-ranking model-generated candidates revealed a significant enrichment of candidates within the Bacteroidetes/Chlorobi and Betaproteobacteria taxonomic groups.

CONCLUSIONS

DRAMMA enables rapid ARG identification for global-scale genomic and metagenomic samples, thus holding promise for the discovery of novel ARGs that lack sequence similarity to any known resistance genes. Further, our model has the potential to facilitate early detection of specific ARGs, potentially influencing the selection of antibiotics administered to patients. Video Abstract.

Diverse anti-defence systems are encoded in the leading region of plasmids

Plasmids are major drivers of gene mobilization by means of horizontal gene transfer and play a key role in spreading antimicrobial resistance among pathogens1,2. Despite various bacterial defence mechanisms such as CRISPR-Cas, restriction-modification systems and SOS-response genes that prevent the invasion of mobile genetic elements3, plasmids robustly transfer within bacterial populations through conjugation4,5. Here we show that the leading region of plasmids, the first to enter recipient cells, is a hotspot for an extensive repertoire of anti-defence systems, encoding anti-CRISPR, anti-restriction, anti-SOS and other counter-defence proteins. We further identified in the leading region a prevalence of promoters known to allow expression from single-stranded DNA6, potentially facilitating rapid protection against bacterial immunity during the early stages of plasmid establishment. We demonstrated experimentally the importance of anti-defence gene localization in the leading region for efficient conjugation. These results indicate that focusing on the leading region of plasmids could lead to the discovery of diverse anti-defence genes. Combined, our findings show a new facet of plasmid dissemination and provide theoretical foundations for developing efficient conjugative delivery systems for natural microbial communities.

Co-localization of antibiotic resistance genes is widespread in the infant gut microbiome and associates with an immature gut microbial composition

BACKGROUND

In environmental bacteria, the selective advantage of antibiotic resistance genes (ARGs) can be increased through co-localization with genes such as other ARGs, biocide resistance genes, metal resistance genes, and virulence genes (VGs). The gut microbiome of infants has been shown to contain numerous ARGs, however, co-localization related to ARGs is unknown during early life despite frequent exposures to biocides and metals from an early age.

RESULTS

We conducted a comprehensive analysis of genetic co-localization of resistance genes in a cohort of 662 Danish children and examined the association between such co-localization and environmental factors as well as gut microbial maturation. Our study showed that co-localization of ARGs with other resistance and virulence genes is common in the early gut microbiome and is associated with gut bacteria that are indicative of low maturity. Statistical models showed that co-localization occurred mainly in the phylum Proteobacteria independent of high ARG content and contig length. We evaluated the stochasticity of co-localization occurrence using enrichment scores. The most common forms of co-localization involved tetracycline and fluoroquinolone resistance genes, and, on plasmids, co-localization predominantly occurred in the form of class 1 integrons. Antibiotic use caused a short-term increase in mobile ARGs, while non-mobile ARGs showed no significant change. Finally, we found that a high abundance of VGs was associated with low gut microbial maturity and that VGs showed even higher potential for mobility than ARGs.

CONCLUSIONS

We found that the phenomenon of co-localization between ARGs and other resistance and VGs was prevalent in the gut at the beginning of life. It reveals the diversity that sustains antibiotic resistance and therefore indirectly emphasizes the need to apply caution in the use of antimicrobial agents in clinical practice, animal husbandry, and daily life to mitigate the escalation of resistance. Video Abstract.

Differential responses of the gut microbiome and resistome to antibiotic exposures in infants and adults

Despite their crucial importance for human health, there is still relatively limited knowledge on how the gut resistome changes or responds to antibiotic treatment across ages, especially in the latter case. Here, we use fecal metagenomic data from 662 Danish infants and 217 young adults to fill this gap. The gut resistomes are characterized by a bimodal distribution driven by E. coli composition. The typical profile of the gut resistome differs significantly between adults and infants, with the latter distinguished by higher gene and plasmid abundances. However, the predominant antibiotic resistance genes (ARGs) are the same. Antibiotic treatment reduces bacterial diversity and increased ARG and plasmid abundances in both cohorts, especially core ARGs. The effects of antibiotic treatments on the gut microbiome last longer in adults than in infants, and different antibiotics are associated with distinct impacts. Overall, this study broadens our current understanding of gut resistome dynamics and the impact of antibiotic treatment across age groups.

Pangenome-level analysis of nucleoid-associated proteins in the Acidithiobacillia class: insights into their functional roles in mobile genetic elements biology

Mobile genetic elements (MGEs) are relevant agents in bacterial adaptation and evolutionary diversification. Stable appropriation of these DNA elements depends on host factors, among which are the nucleoid-associated proteins (NAPs). NAPs are highly abundant proteins that bind and bend DNA, altering its topology and folding, thus affecting all known cellular DNA processes from replication to expression. Even though NAP coding genes are found in most prokaryotic genomes, their functions in host chromosome biology and xenogeneic silencing are only known for a few NAP families. Less is known about the occurrence, abundance, and roles of MGE-encoded NAPs in foreign elements establishment and mobility. In this study, we used a combination of comparative genomics and phylogenetic strategies to gain insights into the diversity, distribution, and functional roles of NAPs within the class Acidithiobacillia with a special focus on their role in MGE biology. Acidithiobacillia class members are aerobic, chemolithoautotrophic, acidophilic sulfur-oxidizers, encompassing substantial genotypic diversity attributable to MGEs. Our search for NAP protein families (PFs) in more than 90 genomes of the different species that conform the class, revealed the presence of 1,197 proteins pertaining to 12 different NAP families, with differential occurrence and conservation across species. Pangenome-level analysis revealed 6 core NAP PFs that were highly conserved across the class, some of which also existed as variant forms of scattered occurrence, in addition to NAPs of taxa-restricted distribution. Core NAPs identified are reckoned as essential based on the conservation of genomic context and phylogenetic signals. In turn, various highly diversified NAPs pertaining to the flexible gene complement of the class, were found to be encoded in known plasmids or, larger integrated MGEs or, present in genomic loci associated with MGE-hallmark genes, pointing to their role in the stabilization/maintenance of these elements in strains and species with larger genomes. Both core and flexible NAPs identified proved valuable as markers, the former accurately recapitulating the phylogeny of the class, and the later, as seed in the bioinformatic identification of novel episomal and integrated mobile elements.

Influence of mobile genetic elements and insertion sequences in long- and short-term adaptive processes of Acidithiobacillus ferrooxidans strains

The recent revision of the Acidithiobacillia class using genomic taxonomy methods has shown that, in addition to the existence of previously unrecognized genera and species, some species of the class harbor levels of divergence that are congruent with ongoing differentiation processes. In this study, we have performed a subspecies-level analysis of sequenced strains of Acidithiobacillus ferrooxidans to prove the existence of distinct sublineages and identify the discriminant genomic/genetic characteristics linked to these sublineages, and to shed light on the processes driving such differentiation. Differences in the genomic relatedness metrics, levels of synteny, gene content, and both integrated and episomal mobile genetic elements (MGE) repertoires support the existence of two subspecies-level taxa within A. ferrooxidans. While sublineage 2A harbors a small plasmid related to pTF5, this episomal MGE is absent in sublineage 2B strains. Likewise, clear differences in the occurrence, coverage and conservation of integrated MGEs are apparent between sublineages. Differential MGE-associated gene cargo pertained to the functional categories of energy metabolism, ion transport, cell surface modification, and defense mechanisms. Inferred functional differences have the potential to impact long-term adaptive processes and may underpin the basis of the subspecies-level differentiation uncovered within A. ferrooxidans. Genome resequencing of iron- and sulfur-adapted cultures of a selected 2A sublineage strain (CCM 4253) showed that both episomal and large integrated MGEs are conserved over twenty generations in either growth condition. In turn, active insertion sequences profoundly impact short-term adaptive processes. The ISAfe1 element was found to be highly active in sublineage 2A strain CCM 4253. Phenotypic mutations caused by the transposition of ISAfe1 into the pstC2 encoding phosphate-transport system permease protein were detected in sulfur-adapted cultures and shown to impair growth on ferrous iron upon the switch of electron donor. The phenotypic manifestation of the △pstC2 mutation, such as a loss of the ability to oxidize ferrous iron, is likely related to the inability of the mutant to secure the phosphorous availability for electron transport-linked phosphorylation coupled to iron oxidation. Depletion of the transpositional △pstC2 mutation occurred concomitantly with a shortening of the iron-oxidation lag phase at later transfers on a ferrous iron-containing medium. Therefore, the pstII operon appears to play an essential role in A. ferrooxidans when cells oxidize ferrous iron. Results highlight the influence of insertion sequences and both integrated and episomal mobile genetic elements in the short- and long-term adaptive processes of A. ferrooxidans strains under changing growth conditions.

Metagenomic-based surveillance systems for antibiotic resistance in non-clinical settings

The success of antibiotics as a therapeutic agent has led to their ineffectiveness. The continuous use and misuse in clinical and non-clinical areas have led to the emergence and spread of antibiotic-resistant bacteria and its genetic determinants. This is a multi-dimensional problem that has now become a global health crisis. Antibiotic resistance research has primarily focused on the clinical healthcare sectors while overlooking the non-clinical sectors. The increasing antibiotic usage in the environment - including animals, plants, soil, and water - are drivers of antibiotic resistance and function as a transmission route for antibiotic resistant pathogens and is a source for resistance genes. These natural compartments are interconnected with each other and humans, allowing the spread of antibiotic resistance via horizontal gene transfer between commensal and pathogenic bacteria. Identifying and understanding genetic exchange within and between natural compartments can provide insight into the transmission, dissemination, and emergence mechanisms. The development of high-throughput DNA sequencing technologies has made antibiotic resistance research more accessible and feasible. In particular, the combination of metagenomics and powerful bioinformatic tools and platforms have facilitated the identification of microbial communities and has allowed access to genomic data by bypassing the need for isolating and culturing microorganisms. This review aimed to reflect on the different sequencing techniques, metagenomic approaches, and bioinformatics tools and pipelines with their respective advantages and limitations for antibiotic resistance research. These approaches can provide insight into resistance mechanisms, the microbial population, emerging pathogens, resistance genes, and their dissemination. This information can influence policies, develop preventative measures and alleviate the burden caused by antibiotic resistance.

Prediction and Inferred Evolution of Acid Tolerance Genes in the Biotechnologically Important Acidihalobacter Genus

Acidihalobacter is a genus of acidophilic, gram-negative bacteria known for its ability to oxidize pyrite minerals in the presence of elevated chloride ions, a capability rare in other iron-sulfur oxidizing acidophiles. Previous research involving Acidihalobacter spp. has focused on their applicability in saline biomining operations and their genetic arsenal that allows them to cope with chloride, metal and oxidative stress. However, an understanding of the molecular adaptations that enable Acidihalobacter spp. to thrive under both acid and chloride stress is needed to provide a more comprehensive understanding of how this genus can thrive in such extreme biomining conditions. Currently, four genomes of the Acidihalobacter genus have been sequenced: Acidihalobacter prosperus DSM 5130^T, Acidihalobacter yilgarnensis DSM 105917^T, Acidihalobacter aeolianus DSM 14174^T, and Acidihalobacter ferrooxydans DSM 14175^T. Phylogenetic analysis shows that the Acidihalobacter genus roots to the Chromatiales class consisting of mostly halophilic microorganisms. In this study, we aim to advance our knowledge of the genetic repertoire of the Acidihalobacter genus that has enabled it to cope with acidic stress. We provide evidence of gene gain events that are hypothesized to help the Acidihalobacter genus cope with acid stress. Potential acid tolerance mechanisms that were found in the Acidihalobacter genomes include multiple potassium transporters, chloride/proton antiporters, glutamate decarboxylase system, arginine decarboxylase system, urease system, slp genes, squalene synthesis, and hopanoid synthesis. Some of these genes are hypothesized to have entered the Acidihalobacter via vertical decent from an inferred non-acidophilic ancestor, however, horizontal gene transfer (HGT) from other acidophilic lineages is probably responsible for the introduction of many acid resistance genes.

Integrative Genomics Sheds Light on Evolutionary Forces Shaping the Acidithiobacillia Class Acidophilic Lifestyle

Extreme acidophiles thrive in environments rich in protons (pH values <3) and often high levels of dissolved heavy metals. They are distributed across the three domains of the Tree of Life including members of the Proteobacteria. The Acidithiobacillia class is formed by the neutrophilic genus Thermithiobacillus along with the extremely acidophilic genera Fervidacidithiobacillus, Igneacidithiobacillus, Ambacidithiobacillus, and Acidithiobacillus. Phylogenomic reconstruction revealed a division in the Acidithiobacillia class correlating with the different pH optima that suggested that the acidophilic genera evolved from an ancestral neutrophile within the Acidithiobacillia. Genes and mechanisms denominated as "first line of defense" were key to explaining the Acidithiobacillia acidophilic lifestyle including preventing proton influx that allows the cell to maintain a near-neutral cytoplasmic pH and differ from the neutrophilic Acidithiobacillia ancestors that lacked these systems. Additional differences between the neutrophilic and acidophilic Acidithiobacillia included the higher number of gene copies in the acidophilic genera coding for "second line of defense" systems that neutralize and/or expel protons from cell. Gain of genes such as hopanoid biosynthesis involved in membrane stabilization at low pH and the functional redundancy for generating an internal positive membrane potential revealed the transition from neutrophilic properties to a new acidophilic lifestyle by shaping the Acidithiobacillaceae genomic structure. The presence of a pool of accessory genes with functional redundancy provides the opportunity to "hedge bet" in rapidly changing acidic environments. Although a core of mechanisms for acid resistance was inherited vertically from an inferred neutrophilic ancestor, the majority of mechanisms, especially those potentially involved in resistance to extremely low pH, were obtained from other extreme acidophiles by horizontal gene transfer (HGT) events.

Metagenomic evidence for co-occurrence of antibiotic, biocide and metal resistance genes in pigs

Antibiotic-resistant pathogens constitute an escalating public health concern. Hence a better understanding of the underlying processes responsible for this expansion is urgently needed. Co-selection of heavy metal/biocide and antibiotic resistance genes (ARGs) has been suggested as one potential mechanism promoting the proliferation of antimicrobial resistance (AMR). This paper aims to elucidate this interplay and exploit differences in antibiotic usage to infer patterns of co-selection by the non-antibiotic factors metals and biocides in the context of pig farming. We examined 278 gut metagenomes from pigs with continuous antibiotic exposure, only at weaning and at no exposure. Metals as growth promoters and biocides as disinfectants are currently used with little restrictions in stock farming. The pigs under continuous antibiotic exposure displayed the highest co-occurrence of ARGs and other genetic elements while the pigs under limited use of antibiotics still showed abundant co-occurrences. Pathogens belonging to Enterobacteriaceae displayed increased co-occurrence phenomena, suggesting that this maintenance is not a random selection process from a mobilized pool but pertains to specific phylogenetic clades. These results suggest that metals and biocides displayed strong selective pressures on ARGs exerted by intensive farming, regardless of the current use of antibiotics.

Assessment of the plasmidome of an extremophilic microbial community from the Diamante Lake, Argentina

Diamante Lake located at 4589 m.a.s.l. in the Andean Puna constitutes an extreme environment. It is exposed to multiple extreme conditions such as an unusually high concentration of arsenic (over 300 mg L^-1) and low oxygen pressure. Microorganisms thriving in the lake display specific genotypes that facilitate survival, which include at least a multitude of plasmid-encoded resistance traits. Hence, the genetic information provided by the plasmids essentially contributes to understand adaptation to different stressors. Though plasmids from cultivable organisms have already been analyzed to the sequence level, the impact of the entire plasmid-borne genetic information on such microbial ecosystem is not known. This study aims at assessing the plasmidome from Diamante Lake, which facilitates the identification of potential hosts and prediction of gene functions as well as the ecological impact of mobile genetic elements. The deep-sequencing analysis revealed a large fraction of previously unknown DNA sequences of which the majority encoded putative proteins of unknown function. Remarkably, functions related to the oxidative stress response, DNA repair, as well as arsenic- and antibiotic resistances were annotated. Additionally, all necessary capacities related to plasmid replication, mobilization and maintenance were detected. Sequences characteristic for megaplasmids and other already known plasmid-associated genes were identified as well. The study highlights the potential of the deep-sequencing approach specifically targeting plasmid populations as it allows to evaluate the ecological impact of plasmids from (cultivable and non-cultivable) microorganisms, thereby contributing to the understanding of the distribution of resistance factors within an extremophilic microbial community.

Prediction of prokaryotic transposases from protein features with machine learning approaches

Identification of prokaryotic transposases (Tnps) not only gives insight into the spread of antibiotic resistance and virulence but the process of DNA movement. This study aimed to develop a classifier for predicting Tnps in bacteria and archaea using machine learning (ML) approaches. We extracted a total of 2751 protein features from the training dataset including 14852 Tnps and 14852 controls, and selected 75 features as predictive signatures using the combined mutual information and least absolute shrinkage and selection operator algorithms. By aggregating these signatures, an ensemble classifier that integrated a collection of individual ML-based classifiers, was developed to identify Tnps. Further validation revealed that this classifier achieved good performance with an average AUC of 0.955, and met or exceeded other common methods. Based on this ensemble classifier, a stand-alone command-line tool designated TnpDiscovery was established to maximize the convenience for bioinformaticians and experimental researchers toward Tnp prediction. This study demonstrates the effectiveness of ML approaches in identifying Tnps, facilitating the discovery of novel Tnps in the future.

digIS: towards detecting distant and putative novel insertion sequence elements in prokaryotic genomes

Background

The insertion sequence elements (IS elements) represent the smallest and the most abundant mobile elements in prokaryotic genomes. It has been shown that they play a significant role in genome organization and evolution. To better understand their function in the host genome, it is desirable to have an effective detection and annotation tool. This need becomes even more crucial when considering rapid-growing genomic and metagenomic data. The existing tools for IS elements detection and annotation are usually based on comparing sequence similarity with a database of known IS families. Thus, they have limited ability to discover distant and putative novel IS elements.

Results

In this paper, we present digIS, a software tool based on profile hidden Markov models assembled from catalytic domains of transposases. It shows a very good performance in detecting known IS elements when tested on datasets with manually curated annotation. The main contribution of digIS is in its ability to detect distant and putative novel IS elements while maintaining a moderate level of false positives. In this category it outperforms existing tools, especially when tested on large datasets of archaeal and bacterial genomes.

Conclusion

We provide digIS, a software tool using a novel approach based on manually curated profile hidden Markov models, which is able to detect distant and putative novel IS elements. Although digIS can find known IS elements as well, we expect it to be used primarily by scientists interested in finding novel IS elements. The tool is available at https://github.com/janka2012/digIS.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04177-6.

A glimpse of antimicrobial resistance gene diversity in kefir and yoghurt

Antimicrobial resistance (AMR) is a global threat gaining more and more practical significance every year. The main determinants of AMR are the antimicrobial resistance genes (ARGs). Since bacteria can share genetic components via horizontal gene transfer, even non-pathogenic bacteria may provide ARG to any pathogens which they become physically close to (e.g. in the human gut). In addition, fermented food naturally contains bacteria in high amounts. In this study, we examined the diversity of ARG content in various kefir and yoghurt samples (products, grains, bacterial strains) using a unified metagenomic approach. We found numerous ARGs of commonly used fermenting bacteria. Even with the strictest filter restrictions, we identified ARGs undermining the efficacy of aminocoumarins, aminoglycosides, carbapenems, cephalosporins, cephamycins, diaminopyrimidines, elfamycins, fluoroquinolones, fosfomycins, glycylcyclines, lincosamides, macrolides, monobactams, nitrofurans, nitroimidazoles, penams, penems, peptides, phenicols, rifamycins, tetracyclines and triclosan. In the case of gene lmrD, we detected genetic environment providing mobility of this ARG. Our findings support the theory that during the fermentation process, the ARG content of foods can grow due to bacterial multiplication. The results presented suggest that the starting culture strains of fermented foods should be monitored and selected in order to decrease the intake of ARGs via foods.

Unlocking Survival Mechanisms for Metal and Oxidative Stress in the Extremely Acidophilic, Halotolerant Acidihalobacter Genus

Microorganisms used for the biohydrometallurgical extraction of metals from minerals must be able to survive high levels of metal and oxidative stress found in bioleaching environments. The Acidihalobacter genus consists of four species of halotolerant, iron–sulfur-oxidizing acidophiles that are unique in their ability to tolerate chloride and acid stress while simultaneously bioleaching minerals. This paper uses bioinformatic tools to predict the genes and mechanisms used by Acidihalobacter members in their defense against a wide range of metals and oxidative stress. Analysis revealed the presence of multiple conserved mechanisms of metal tolerance. Ac. yilgarnensis F5^T, the only member of this genus that oxidizes the mineral chalcopyrite, contained a 39.9 Kb gene cluster consisting of 40 genes encoding mobile elements and an array of proteins with direct functions in copper resistance. The analysis also revealed multiple strategies that the Acidihalobacter members can use to tolerate high levels of oxidative stress. Three of the Acidihalobacter genomes were found to contain genes encoding catalases, which are not common to acidophilic microorganisms. Of particular interest was a rubrerythrin genomic cluster containing genes that have a polyphyletic origin of stress-related functions.

The CRISPR-Cas systems were selectively inactivated during evolution of Bacillus cereus group for adaptation to diverse environments

CRISPR-Cas systems are considered as barriers to horizontal gene transfer (HGT). However, the influence of such systems on HGT within species is unclear. Also, little is known about the impact of CRISPR-Cas systems on bacterial evolution at the population level. Here, using Bacillus cereus sensu lato as model, we investigate the interplay between CRISPR-Cas systems and HGT at the population scale. We found that only a small fraction of the strains have CRISPR-Cas systems (13.9% of 1871), and most of such systems are defective based on their gene content analysis. Comparative genomic analysis revealed that the CRISPR-Cas systems are barriers to HGT within this group, since strains harboring active systems contain less mobile genetic elements (MGEs), have lower fraction of unique genes and also display limited environmental distributions than strains without active CRISPR-Cas systems. The introduction of a functional CRISPR-Cas system into a strain lacking the system resulted in reduced adaptability to various stresses and decreased pathogenicity of the transformant strain, indicating that B. cereus group strains could benefit from inactivating such systems. Our work provides a large-scale case to support that the CRISPR-Cas systems are barriers to HGT within species, and that in the B. cereus group the inactivation of CRISPR-Cas systems correlated with acquisition of MGEs that could result in better adaptation to diverse environments.

New candidates for regulated gene integrity revealed through precise mapping of integrative genetic elements

Integrative genetic elements (IGEs) are mobile multigene DNA units that integrate into and excise from host bacterial genomes. Each IGE usually targets a specific site within a conserved host gene, integrating in a manner that preserves target gene function. However, a small number of bacterial genes are known to be inactivated upon IGE integration and reactivated upon excision, regulating phenotypes of virulence, mutation rate, and terminal differentiation in multicellular bacteria. The list of regulated gene integrity (RGI) cases has been slow-growing because IGEs have been challenging to precisely and comprehensively locate in genomes. We present software (TIGER) that maps IGEs with unprecedented precision and without attB site bias. TIGER uses a comparative genomic, ping-pong BLAST approach, based on the principle that the IGE integration module (i.e. its int-attP region) is cohesive. The resultant IGEs from 2168 genomes, along with integrase phylogenetic analysis and gene inactivation tests, revealed 19 new cases of genes whose integrity is regulated by IGEs (including dut, eccCa1, gntT, hrpB, merA, ompN, prkA, tqsA, traG, yifB, yfaT and ynfE), as well as recovering previously known cases (in sigK, spsM, comK, mlrA and hlb genes). It also recovered known clades of site-promiscuous integrases and identified possible new ones.

Evolution of Predicted Acid Resistance Mechanisms in the Extremely Acidophilic Leptospirillum Genus

Organisms that thrive in extremely acidic environments (≤pH 3.5) are of widespread importance in industrial applications, environmental issues, and evolutionary studies. Leptospirillum spp. constitute the only extremely acidophilic microbes in the phylogenetically deep-rooted bacterial phylum Nitrospirae. Leptospirilli are Gram-negative, obligatory chemolithoautotrophic, aerobic, ferrous iron oxidizers. This paper predicts genes that Leptospirilli use to survive at low pH and infers their evolutionary trajectory. Phylogenetic and other bioinformatic approaches suggest that these genes can be classified into (i) "first line of defense", involved in the prevention of the entry of protons into the cell, and (ii) neutralization or expulsion of protons that enter the cell. The first line of defense includes potassium transporters, predicted to form an inside positive membrane potential, spermidines, hopanoids, and Slps (starvation-inducible outer membrane proteins). The "second line of defense" includes proton pumps and enzymes that consume protons. Maximum parsimony, clustering methods, and gene alignments are used to infer the evolutionary trajectory that potentially enabled the ancestral Leptospirillum to transition from a postulated circum-neutral pH environment to an extremely acidic one. The hypothesized trajectory includes gene gains/loss events driven extensively by horizontal gene transfer, gene duplications, gene mutations, and genomic rearrangements.

Comparative Genomic Analysis of Rhodococcus equi: An Insight into Genomic Diversity and Genome Evolution

Rhodococcus equi, a member of the Rhodococcus genus, is a gram-positive pathogenic bacterium. Rhodococcus possesses an open pan-genome that constitutes the basis of its high genomic diversity and allows for adaptation to specific niche conditions and the changing host environments. Our analysis further showed that the core genome of R. equi contributes to the pathogenicity and niche adaptation of R. equi. Comparative genomic analysis revealed that the genomes of R. equi shared identical collinearity relationship, and heterogeneity was mainly acquired by means of genomic islands and prophages. Moreover, genomic islands in R. equi were always involved in virulence, resistance, or niche adaptation and possibly working with prophages to cause the majority of genome expansion. These findings provide an insight into the genomic diversity, evolution, and structural variation of R. equi and a valuable resource for functional genomic studies.

Functional Traits Co-Occurring with Mobile Genetic Elements in the Microbiome of the Atacama Desert

Mobile genetic elements (MGEs) play an essential role in bacterial adaptation and evolution. These elements are enriched within bacterial communities from extreme environments. However, very little is known if specific genes co-occur with MGEs in extreme environments and, if so, what their function is. We used shotgun-sequencing to analyse the metagenomes of 12 soil samples and characterized the composition of MGEs and the genes co-occurring with them. The samples ranged from less arid coastal sites to the inland hyperarid core of the Atacama Desert, as well as from sediments below boulders, protected from UV-irradiation. MGEs were enriched at the hyperarid sites compared with sediments from below boulders and less arid sites. MGEs were mostly co-occurring with genes belonging to the Cluster Orthologous Group (COG) categories “replication, recombination and repair,” “transcription” and “signal transduction mechanisms.” In general, genes coding for transcriptional regulators and histidine kinases were the most abundant genes proximal to MGEs. Genes involved in energy production were significantly enriched close to MGEs at the hyperarid sites. For example, dehydrogenases, reductases, hydrolases and chlorite dismutase and other enzymes linked to nitrogen metabolism such as nitrite- and nitro-reductase. Stress response genes, including genes involved in antimicrobial and heavy metal resistance genes, were rarely found near MGEs. The present study suggests that MGEs could play an essential role in the adaptation of the soil microbiome in hyperarid desert soils by the modulation of housekeeping genes such as those involved in energy production.

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.

Oral administration of antibiotics increased the potential mobility of bacterial resistance genes in the gut of the fish Piaractus mesopotamicus

BACKGROUND

Aquaculture is on the rise worldwide, and the use of antibiotics is fostering higher production intensity. However, recent findings suggest that the use of antibiotics comes at the price of increased antibiotic resistance. Yet, the effect of the oral administration of antibiotics on the mobility of microbial resistance genes in the fish gut is not well understood. In the present study, Piaractus mesopotamicus was used as a model to evaluate the effect of the antimicrobial florfenicol on the diversity of the gut microbiome as well as antibiotic resistance genes (ARGs) and mobile genetic elements (MGEs) using a metagenomic approach.

RESULTS

The total relative abundance of ARGs and MGEs significantly increased during the antibiotic exposure. Additionally, phage integrases, transposases, and transposons flanking ARGs accumulated in the gut microbiome of P. mesopotamicus because of the antibiotic exposure. MGEs co-occurring with ARGs showed a significant positive correlation with the total ARGs found. Furthermore, shifts in the gut microbiome towards well-known putative pathogens such as Salmonella, Plesiomonas, and Citrobacter were observed following florfenicol treatment. Mainly Plesiomonas and Citrobacter harbored genes that code for multidrug and phenicol efflux pumps. Moreover, several genes related to RNA processing and modification, cell motility, SOS response, and extracellular structure were enriched due to the antibiotic application. The observed effects were visible during the complete application phase and disappeared at the post-exposure phase.

CONCLUSIONS

Our findings suggest that the oral administration of antibiotics increases the potential for MGE-mediated exchange of ARGs in the gut of fish and could contribute to the enrichment and dispersion of ARGs in aquaculture systems. Importantly, this increase in the potential for ARGs exchange could be an effect of changes in community structure and/or ARG mobilization.

ISEScan: automated identification of insertion sequence elements in prokaryotic genomes

Motivation

The insertion sequence (IS) elements are the smallest but most abundant autonomous transposable elements in prokaryotic genomes, which play a key role in prokaryotic genome organization and evolution. With the fast growing genomic data, it is becoming increasingly critical for biology researchers to be able to accurately and automatically annotate ISs in prokaryotic genome sequences. The available automatic IS annotation systems are either providing only incomplete IS annotation or relying on the availability of existing genome annotations. Here, we present a new IS elements annotation pipeline to address these issues.

Results

ISEScan is a highly sensitive software pipeline based on profile hidden Markov models constructed from manually curated IS elements. ISEScan performs better than existing IS annotation systems when tested on prokaryotic genomes with curated annotations of IS elements. Applying it to 2784 prokaryotic genomes, we report the global distribution of IS families across taxonomic clades in Archaea and Bacteria.

Availability and implementation

ISEScan is implemented in Python and released as an open source software at https://github.com/xiezhq/ISEScan.

Contact

hatang@indiana.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

A closely-related clade of globally distributed bloom-forming cyanobacteria within the Nostocales

In order to better understand the relationships among current Nostocales cyanobacterial blooms, eight genomes were sequenced from cultured isolates or from environmental metagenomes of recent planktonic Nostocales blooms. Phylogenomic analysis of publicly available sequences placed the new genomes among a group of 15 genomes from four continents in a distinct ADA clade (Anabaena/Dolichospermum/Aphanizomenon) within the Nostocales. This clade contains four species-level groups, two of which include members with both Anabaena-like and Aphanizomenon flos-aquae-like morphology. The genomes contain many repetitive genetic elements and a sizable pangenome, in which ABC-type transporters are highly represented. Alongside common core genes for photosynthesis, the differentiation of N₂-fixing heterocysts, and the uptake and incorporation of the major nutrients P, N and S, we identified several gene pathways in the pangenome that may contribute to niche partitioning. Genes for problematic secondary metabolites-cyanotoxins and taste-and-odor compounds-were sporadically present, as were other polyketide synthase (PKS) and nonribosomal peptide synthetase (NRPS) gene clusters. By contrast, genes predicted to encode the ribosomally generated bacteriocin peptides were found in all genomes.

A Practical Guide for Comparative Genomics of Mobile Genetic Elements in Prokaryotic Genomes

Mobile genetic elements (MGEs) are an important feature of prokaryote genomes but are seldom well annotated and, consequently, are often underestimated. MGEs include transposons (Tn), insertion sequences (ISs), prophages, genomic islands (GEIs), integrons, and integrative and conjugative elements (ICEs). They are intimately involved in genome evolution and promote phenomena such as genomic expansion and rearrangement, emergence of virulence and pathogenicity, and symbiosis. In spite of the annotation bottleneck, there are so far at least 75 different programs and databases dedicated to prokaryotic MGE analysis and annotation, and this number is rapidly growing. Here, we present a practical guide to explore, compare, and visualize prokaryote MGEs using a combination of available software and databases tailored to small scale genome analyses. This protocol can be coupled with expert MGE annotation and exploited for evolutionary and comparative genomic analyses.

Phage Genome Annotation Using the RAST Pipeline

Phages are complex biomolecular machineries that have to survive in a bacterial world. Phage genomes show many adaptations to their lifestyle such as shorter genes, reduced capacity for redundant DNA sequences, and the inclusion of tRNAs in their genomes. In addition, phages are not free-living, they require a host for replication and survival. These unique adaptations provide challenges for the bioinformatics analysis of phage genomes. In particular, ORF calling, genome annotation, noncoding RNA (ncRNA) identification, and the identification of transposons and insertions are all complicated in phage genome analysis. We provide a road map through the phage genome annotation pipeline, and discuss the challenges and solutions for phage genome annotation as we have implemented in the rapid annotation using subsystems (RAST) pipeline.

Whole genome sequence of two Rathayibacter toxicus strains reveals a tunicamycin biosynthetic cluster similar to Streptomyces chartreusis

Rathayibacter toxicus is a forage grass associated Gram-positive bacterium of major concern to food safety and agriculture. This species is listed by USDA-APHIS as a plant pathogen select agent because it produces a tunicamycin-like toxin that is lethal to livestock and may be vectored by nematode species native to the U.S. The complete genomes of two strains of R. toxicus, including the type strain FH-79, were sequenced and analyzed in comparison with all available, complete R. toxicus genomes. Genome sizes ranged from 2,343,780 to 2,394,755 nucleotides, with 2079 to 2137 predicted open reading frames; all four strains showed remarkable synteny over nearly the entire genome, with only a small transposed region. A cluster of genes with similarity to the tunicamycin biosynthetic cluster from Streptomyces chartreusis was identified. The tunicamycin gene cluster (TGC) in R. toxicus contained 14 genes in two transcriptional units, with all of the functional elements for tunicamycin biosynthesis present. The TGC had a significantly lower GC content (52%) than the rest of the genome (61.5%), suggesting that the TGC may have originated from a horizontal transfer event. Further analysis indicated numerous remnants of other potential horizontal transfer events are present in the genome. In addition to the TGC, genes potentially associated with carotenoid and exopolysaccharide production, bacteriocins and secondary metabolites were identified. A CRISPR array is evident. There were relatively few plant-associated cell-wall hydrolyzing enzymes, but there were numerous secreted serine proteases that share sequence homology to the pathogenicity-associated protein Pat-1 of Clavibacter michiganensis. Overall, the genome provides clear insight into the possible mechanisms for toxin production in R. toxicus, providing a basis for future genetic approaches.

Comparative Genomics of Bacillus thuringiensis Reveals a Path to Specialized Exploitation of Multiple Invertebrate Hosts

Understanding the genetic basis of host shifts is a key genomic question for pathogen and parasite biology. The Bacillus cereus group, which encompasses Bacillus thuringiensis and Bacillus anthracis, contains pathogens that can infect insects, nematodes, and vertebrates. Since the target range of the essential virulence factors (Cry toxins) and many isolates is well known, this group presents a powerful system for investigating how pathogens can diversify and adapt to phylogenetically distant hosts. Specialization to exploit insects occurs at the level of the major clade and is associated with substantial changes in the core genome, and host switching between insect orders has occurred repeatedly within subclades. The transfer of plasmids with linked cry genes may account for much of the adaptation to particular insect orders, and network analysis implies that host specialization has produced strong associations between key toxin genes with similar targets. Analysis of the distribution of plasmid minireplicons shows that plasmids with orf156 and orf157, which carry genes encoding toxins against Lepidoptera or Diptera, were contained only by B. thuringiensis in the specialized insect clade (clade 2), indicating that tight genome/plasmid associations have been important in adaptation to invertebrate hosts. Moreover, the accumulation of multiple virulence factors on transposable elements suggests that cotransfer of diverse virulence factors is advantageous in terms of expanding the insecticidal spectrum, overcoming insect resistance, or through gains in pathogenicity via synergistic interactions between toxins.IMPORTANCE Population genomics have provided many new insights into the formation, evolution, and dynamics of bacterial pathogens of humans and other higher animals, but these pathogens usually have very narrow host ranges. As a pathogen of insects and nematodes, Bacillus thuringiensis, which produces toxins showing toxicity to many orders of insects and other invertebrates, can be used as a model to study the evolution of pathogens with wide host ranges. Phylogenomic analysis revealed that host specialization and switching occur at the level of the major clade and subclade, respectively. A toxin gene co-occurrence network indicates that multiple toxins with similar targets were accumulated by the same cell in the whole species. This accumulation may be one of the strategies that B. thuringiensis has used to fight against host resistance. This kind of formation and evolution of pathogens represents a different path used against multiple invertebrate hosts from that used against higher animals.

Bioinformatic Analyses of Unique (Orphan) Core Genes of the Genus Acidithiobacillus: Functional Inferences and Use As Molecular Probes for Genomic and Metagenomic/Transcriptomic Interrogation

Using phylogenomic and gene compositional analyses, five highly conserved gene families have been detected in the core genome of the phylogenetically coherent genus Acidithiobacillus of the class Acidithiobacillia. These core gene families are absent in the closest extant genus Thermithiobacillus tepidarius that subtends the Acidithiobacillus genus and roots the deepest in this class. The predicted proteins encoded by these core gene families are not detected by a BLAST search in the NCBI non-redundant database of more than 90 million proteins using a relaxed cut-off of 1.0e⁻⁵. None of the five families has a clear functional prediction. However, bioinformatic scrutiny, using pI prediction, motif/domain searches, cellular location predictions, genomic context analyses, and chromosome topology studies together with previously published transcriptomic and proteomic data, suggests that some may have functions associated with membrane remodeling during cell division perhaps in response to pH stress. Despite the high level of amino acid sequence conservation within each family, there is sufficient nucleotide variation of the respective genes to permit the use of the DNA sequences to distinguish different species of Acidithiobacillus, making them useful additions to the armamentarium of tools for phylogenetic analysis. Since the protein families are unique to the Acidithiobacillus genus, they can also be leveraged as probes to detect the genus in environmental metagenomes and metatranscriptomes, including industrial biomining operations, and acid mine drainage (AMD).

Gene Loss and Horizontal Gene Transfer Contributed to the Genome Evolution of the Extreme Acidophile "Ferrovum"

Acid mine drainage (AMD), associated with active and abandoned mining sites, is a habitat for acidophilic microorganisms that gain energy from the oxidation of reduced sulfur compounds and ferrous iron and that thrive at pH below 4. Members of the recently proposed genus “Ferrovum” are the first acidophilic iron oxidizers to be described within the Betaproteobacteria. Although they have been detected as typical community members in AMD habitats worldwide, knowledge of their phylogenetic and metabolic diversity is scarce. Genomics approaches appear to be most promising in addressing this lacuna since isolation and cultivation of “Ferrovum” has proven to be extremely difficult and has so far only been successful for the designated type strain “Ferrovum myxofaciens” P3G. In this study, the genomes of two novel strains of “Ferrovum” (PN-J185 and Z-31) derived from water samples of a mine water treatment plant were sequenced. These genomes were compared with those of “Ferrovum” sp. JA12 that also originated from the mine water treatment plant, and of the type strain (P3G). Phylogenomic scrutiny suggests that the four strains represent three “Ferrovum” species that cluster in two groups (1 and 2). Comprehensive analysis of their predicted metabolic pathways revealed that these groups harbor characteristic metabolic profiles, notably with respect to motility, chemotaxis, nitrogen metabolism, biofilm formation and their potential strategies to cope with the acidic environment. For example, while the “F. myxofaciens” strains (group 1) appear to be motile and diazotrophic, the non-motile group 2 strains have the predicted potential to use a greater variety of fixed nitrogen sources. Furthermore, analysis of their genome synteny provides first insights into their genome evolution, suggesting that horizontal gene transfer and genome reduction in the group 2 strains by loss of genes encoding complete metabolic pathways or physiological features contributed to the observed diversification.

The Identification of Novel Diagnostic Marker Genes for the Detection of Beer Spoiling Pediococcus damnosus Strains Using the BlAst Diagnostic Gene findEr

As the number of bacterial genomes increases dramatically, the demand for easy to use tools with transparent functionality and comprehensible output for applied comparative genomics grows as well. We present BlAst Diagnostic Gene findEr (BADGE), a tool for the rapid prediction of diagnostic marker genes (DMGs) for the differentiation of bacterial groups (e.g. pathogenic / nonpathogenic). DMG identification settings can be modified easily and installing and running BADGE does not require specific bioinformatics skills. During the BADGE run the user is informed step by step about the DMG finding process, thus making it easy to evaluate the impact of chosen settings and options. On the basis of an example with relevance for beer brewing, being one of the oldest biotechnological processes known, we show a straightforward procedure, from phenotyping, genome sequencing, assembly and annotation, up to a discriminant marker gene PCR assay, making comparative genomics a means to an end. The value and the functionality of BADGE were thoroughly examined, resulting in the successful identification and validation of an outstanding novel DMG (fabZ) for the discrimination of harmless and harmful contaminations of Pediococcus damnosus, which can be applied for spoilage risk determination in breweries. Concomitantly, we present and compare five complete P. damnosus genomes sequenced in this study, finding that the ability to produce the unwanted, spoilage associated off-flavor diacetyl is a plasmid encoded trait in this important beer spoiling species.

Towards the bridging of molecular genetics data across Xenopus species

BACKGROUND

The clawed African frog Xenopus laevis has been one of the main vertebrate models for studies in developmental biology. However, for genetic studies, Xenopus tropicalis has been the experimental model of choice because it shorter life cycle and due to a more tractable genome that does not result from genome duplication as in the case of X. laevis. Today, although still organized in a large number of scaffolds, nearly 85% of X. tropicalis and 89% of X. laevis genomes have been sequenced. There is expectation for a comparative physical map that can be used as a Rosetta Stone between X. laevis genetic studies and X. tropicalis genomic research.

RESULTS

In this work, we have mapped using coarse-grained alignment the 18 chromosomes of X. laevis, release 9.1, on the 10 reference scaffolds representing the haploid genome of X. tropicalis, release 9.0. After validating the mapping with theoretical data, and estimating reference averages of genome sequence identity, 37 to 44% between the two species, we have carried out a synteny analysis for 2,112 orthologous genes. We found that 99.6% of genes are in the same organization.

CONCLUSIONS

Taken together, our results make possible to establish the correspondence between 62 and 65.5% of both genomes, percentage of identity, synteny and automatic annotation of transcripts of both species, providing a new and more comprehensive tool for comparative analysis of these two species, by allowing to bridge molecular genetics data among them.

Genome Analysis of the Biotechnologically Relevant Acidophilic Iron Oxidising Strain JA12 Indicates Phylogenetic and Metabolic Diversity within the Novel Genus "Ferrovum"

Background

Members of the genus “Ferrovum” are ubiquitously distributed in acid mine drainage (AMD) waters which are characterised by their high metal and sulfate loads. So far isolation and microbiological characterisation have only been successful for the designated type strain “Ferrovum myxofaciens” P3G. Thus, knowledge about physiological characteristics and the phylogeny of the genus “Ferrovum” is extremely scarce.

Objective

In order to access the wider genetic pool of the genus “Ferrovum” we sequenced the genome of a “Ferrovum”-containing mixed culture and successfully assembled the almost complete genome sequence of the novel “Ferrovum” strain JA12.

Phylogeny and Lifestyle

The genome-based phylogenetic analysis indicates that strain JA12 and the type strain represent two distinct “Ferrovum” species. “Ferrovum” strain JA12 is characterised by an unusually small genome in comparison to the type strain and other iron oxidising bacteria. The prediction of nutrient assimilation pathways suggests that “Ferrovum” strain JA12 maintains a chemolithoautotrophic lifestyle utilising carbon dioxide and bicarbonate, ammonium and urea, sulfate, phosphate and ferrous iron as carbon, nitrogen, sulfur, phosphorous and energy sources, respectively.

Unique Metabolic Features

The potential utilisation of urea by “Ferrovum” strain JA12 is moreover remarkable since it may furthermore represent a strategy among extreme acidophiles to cope with the acidic environment. Unlike other acidophilic chemolithoautotrophs “Ferrovum” strain JA12 exhibits a complete tricarboxylic acid cycle, a metabolic feature shared with the closer related neutrophilic iron oxidisers among the Betaproteobacteria including Sideroxydans lithotrophicus and Thiobacillus denitrificans. Furthermore, the absence of characteristic redox proteins involved in iron oxidation in the well-studied acidophiles Acidithiobacillus ferrooxidans (rusticyanin) and Acidithiobacillus ferrivorans (iron oxidase) indicates the existence of a modified pathway in “Ferrovum” strain JA12. Therefore, the results of the present study extend our understanding of the genus “Ferrovum” and provide a comprehensive framework for future comparative genome and metagenome studies.

Degradation of Benzene by Pseudomonas veronii 1YdBTEX2 and 1YB2 Is Catalyzed by Enzymes Encoded in Distinct Catabolism Gene Clusters

Pseudomonas veronii 1YdBTEX2, a benzene and toluene degrader, and Pseudomonas veronii 1YB2, a benzene degrader, have previously been shown to be key players in a benzene-contaminated site. These strains harbor unique catabolic pathways for the degradation of benzene comprising a gene cluster encoding an isopropylbenzene dioxygenase where genes encoding downstream enzymes were interrupted by stop codons. Extradiol dioxygenases were recruited from gene clusters comprising genes encoding a 2-hydroxymuconic semialdehyde dehydrogenase necessary for benzene degradation but typically absent from isopropylbenzene dioxygenase-encoding gene clusters. The benzene dihydrodiol dehydrogenase-encoding gene was not clustered with any other aromatic degradation genes, and the encoded protein was only distantly related to dehydrogenases of aromatic degradation pathways. The involvement of the different gene clusters in the degradation pathways was suggested by real-time quantitative reverse transcription PCR.

ISQuest: finding insertion sequences in prokaryotic sequence fragment data

MOTIVATION

Insertion sequences (ISs) are transposable elements present in most bacterial and archaeal genomes that play an important role in genomic evolution. The increasing availability of sequenced prokaryotic genomes offers the opportunity to study ISs comprehensively, but development of efficient and accurate tools is required for discovery and annotation. Additionally, prokaryotic genomes are frequently deposited as incomplete, or draft stage because of the substantial cost and effort required to finish genome assembly projects. Development of methods to identify IS directly from raw sequence reads or draft genomes are therefore desirable. Software tools such as Optimized Annotation System for Insertion Sequences and IScan currently identify IS elements in completely assembled and annotated genomes; however, to our knowledge no methods have been developed to identify ISs from raw fragment data or partially assembled genomes. We have developed novel methods to solve this computationally challenging problem, and implemented these methods in the software package ISQuest. This software identifies bacterial ISs and their sequence elements-inverted and direct repeats-in raw read data or contigs using flexible search parameters. ISQuest is capable of finding ISs in hundreds of partially assembled genomes within hours, making it a valuable high-throughput tool for a global search of IS elements. We tested ISQuest on simulated read libraries of 3810 complete bacterial genomes and plasmids in GenBank and were capable of detecting 82% of the ISs and transposases annotated in GenBank with 80% sequence identity.

CONTACT

abiswas@cs.odu.edu.

The capsule biosynthesis locus of Haemophilus influenzae shows conspicuous similarity to the corresponding locus in Haemophilus sputorum and may have been recruited from this species by horizontal gene transfer

The newly described species Haemophilus sputorum has been cultured from the upper respiratory tract of humans and appears to have little pathogenic potential. The species encodes a capsular biosynthesis locus of approximately 12  kb composed of three distinct regions. Region I and III genes, involved in export and processing of the capsular material, show high similarity to the corresponding genes in capsulate lineages of the pathogenic species Haemophilus influenzae; indeed, standard bexA and bexB PCRs for detection of capsulated strains of H. influenzae give positive results with strains of H. sputorum. Three ORFs are present in region II of the sequenced strain of H. sputorum, of which a putative phosphotransferase showed homology with corresponding genes from H. influenzae serotype c and f. Phylogenetic analysis of housekeeping genes from 24 Pasteurellaceae species showed that H. sputorum was only distantly related to H. influenzae. In contrast to H. influenzae, the capsule locus in H. sputorum is not associated with transposases or other transposable elements. Our data suggest that the capsule locus of capsulate lineages of H. influenzae may have been recruited relatively recently from the commensal species H. sputorum by horizontal gene transfer.

Draft Genome Sequence of the Iron-Oxidizing, Acidophilic, and Halotolerant "Thiobacillus prosperus" Type Strain DSM 5130

“Thiobacillus prosperus” is a halotolerant mesophilic acidophile that gains energy through iron and sulfur oxidation. Its physiology is poorly understood. Here, we describe the principal genomic features of the type strain of T. prosperus, DSM 5130. This is the first public genome sequence of an acidophilic halotolerant bacterium.

A new genome of Acidithiobacillus thiooxidans provides insights into adaptation to a bioleaching environment

Acidithiobacillus thiooxidans is a sulfur oxidizing acidophilic bacterium found in many sulfur-rich environments. It is particularly interesting due to its role in bioleaching of sulphide minerals. In this work, we report the genome sequence of At. thiooxidans Licanantay, the first strain from a copper mine to be sequenced and currently used in bioleaching industrial processes. Through comparative genomic analysis with two other At. thiooxidans non-metal mining strains (ATCC 19377 and A01) we determined that these strains share a large core genome of 2109 coding sequences and a high average nucleotide identity over 98%. Nevertheless, the presence of 841 strain-specific genes (absent in other At. thiooxidans strains) suggests a particular adaptation of Licanantay to its specific biomining environment. Among this group, we highlight genes encoding for proteins involved in heavy metal tolerance, mineral cell attachment and cysteine biosynthesis. Several of these genes were located near genetic motility genes (e.g. transposases and integrases) in genomic regions of over 10 kbp absent in the other strains, suggesting the presence of genomic islands in the Licanantay genome probably produced by horizontal gene transfer in mining environments.

Comparative genome analysis of pathogenic and non-pathogenic Clavibacter strains reveals adaptations to their lifestyle

BACKGROUND

The genus Clavibacter harbors economically important plant pathogens infecting agricultural crops such as potato and tomato. Although the vast majority of Clavibacter strains are pathogenic, there is an increasing number of non-pathogenic isolates reported. Non-pathogenic Clavibacter strains isolated from tomato seeds are particularly problematic because they affect the current detection and identification tests for Clavibacter michiganensis subsp. michiganensis (Cmm), which is regulated with a zero tolerance in tomato seed. Their misidentification as pathogenic Cmm hampers a clear judgment on the seed quality and health.

RESULTS

To get more insight in the genetic features linked to the lifestyle of these bacteria, a whole-genome sequence of the tomato seed-borne non-pathogenic Clavibacter LMG 26808 was determined. To gain a better understanding of the molecular determinants of pathogenicity, the genome sequence of LMG 26808 was compared with that of the pathogenic Cmm strain (NCPPB 382). The comparative analysis revealed that LMG 26808 does not contain plasmids pCM1 and pCM2 and also lacks the majority of important virulence factors described so far for pathogenic Cmm. This explains its apparent non-pathogenic nature in tomato plants. Moreover, the genome analysis of LMG 26808 detected sequences from a plasmid originating from a member of Enterobacteriaceae/Klebsiella relative. Genes received that way and coding for antibiotic resistance may provide a competitive advantage for survival of LMG 26808 in its ecological niche. Genetically, LMG 26808 was the most similar to the pathogenic Cmm NCPPB 382 but contained more mobile genetic elements. The genome of this non-pathogenic Clavibacter strain contained also a high number of transporters and regulatory genes.

CONCLUSIONS

The genome sequence of the non-pathogenic Clavibacter strain LMG 26808 and the comparative analyses with other pathogenic Clavibacter strains provided a better understanding of the genetic bases of virulence and adaptation mechanisms present in the genus Clavibacter.

Architecture and gene repertoire of the flexible genome of the extreme acidophile Acidithiobacillus caldus

BACKGROUND

Acidithiobacillus caldus is a sulfur oxidizing extreme acidophile and the only known mesothermophile within the Acidithiobacillales. As such, it is one of the preferred microbes for mineral bioprocessing at moderately high temperatures. In this study, we explore the genomic diversity of A. caldus strains using a combination of bioinformatic and experimental techniques, thus contributing first insights into the elucidation of the species pangenome.

PRINCIPAL FINDINGS

Comparative sequence analysis of A. caldus ATCC 51756 and SM-1 indicate that, despite sharing a conserved and highly syntenic genomic core, both strains have unique gene complements encompassing nearly 20% of their respective genomes. The differential gene complement of each strain is distributed between the chromosomal compartment, one megaplasmid and a variable number of smaller plasmids, and is directly associated to a diverse pool of mobile genetic elements (MGE). These include integrative conjugative and mobilizable elements, genomic islands and insertion sequences. Some of the accessory functions associated to these MGEs have been linked previously to the flexible gene pool in microorganisms inhabiting completely different econiches. Yet, others had not been unambiguously mapped to the flexible gene pool prior to this report and clearly reflect strain-specific adaption to local environmental conditions.

SIGNIFICANCE

For many years, and because of DNA instability at low pH and recurrent failure to genetically transform acidophilic bacteria, gene transfer in acidic environments was considered negligible. Findings presented herein imply that a more or less conserved pool of actively excising MGEs occurs in the A. caldus population and point to a greater frequency of gene exchange in this econiche than previously recognized. Also, the data suggest that these elements endow the species with capacities to withstand the diverse abiotic and biotic stresses of natural environments, in particular those associated with its extreme econiche.