Detecting genomic islands using bioinformatics approaches

Evolution of the pheV-tRNA integrated genomic island in Escherichia coli

Escherichia coli exhibit extensive genetic diversity at the genome level, particularly within their accessory genome. The tRNA integrated genomic islands (GIs), a part of the E. coli accessory genome, play an important role in pathogenicity. However, studies examining the evolution of GIs have been challenging due to their large size, considerable gene content variation and fragmented assembly in draft genomes. Here we examined the evolution of the GI integrated at pheV-tRNA (GI-pheV), with a primary focus on uropathogenic E. coli (UPEC) and the globally disseminated multidrug resistant ST131 clone. We show the gene content of GI-pheV is highly diverse and arranged in a modular configuration, with the P4 integrase encoding gene intP4 the only conserved gene. Despite this diversity, the GI-pheV gene content displayed conserved features among strains from the same pathotype. In ST131, GI-pheV corresponding to the reference strain EC958 (EC958_GI-pheV) was found in ~90% of strains. Phylogenetic analyses suggested that GI-pheV in ST131 has evolved together with the core genome, with the loss/gain of specific modules (or the entire GI) linked to strain specific events. Overall, we show GI-pheV exhibits a dynamic evolutionary pathway, in which modules and genes have evolved through multiple events including insertions, deletions and recombination.

Regional variation and adaptive evolution in Bifidobacterium pseudocatenulatum: Insights into genomic and functional diversity in human gut

Bifidobacterium pseudocatenulatum is a prevalent gut microbe in humans of all ages and plays a crucial role in host health. However, its adaptive evolutionary characteristics remain poorly understood. This study analyzed the genome of 247 B. pseudocatenulatum isolates from Chinese, Vietnamese, Japanese and other region populations using population genomics and functional genomics. Our findings revealed high genetic heterogeneity and regional clustering within B. pseudocatenulatum isolates. Significant differences were observed in genome characteristics, phylogeny, and functional genes. Specifically, Chinese and Vietnamese isolates exhibited a higher abundance of genes involved in the metabolism of plant-derived carbohydrates (GH13, GH43, and GH5 enzyme families), aligning with the predominantly vegetable-, wheat- and fruit-based diets of these populations. Additionally, we found widespread transmission of antibiotic resistance genes (tetO and tetW) through mobile genetic elements, such as genomic islands (GIs), resulting in substantial intra-regional differences. Our findings highlight distinct adaptive evolution in B. pseudocatenulatum driven by gene specialization, possibly in response to regional variations in diet and lifestyle. This study sheds light on bifidobacteria colonization mechanisms in the host gut. IMPORTANCE: Gut microbiota, as a key link in the gut-brain axis, helps to maintain the health of the organism, among which, Bifidobacterium pseudocatenulatum (B. pseudocatenulatum) is an important constituent member of the gut microbiota, which plays an important role in maintaining the balance of gut microbiota. The probiotic properties of B. pseudocatenulatum have been widely elaborated, and in order to excavate its evolutionary features at the genomic level, here we focused on the genetic background and evolutionary mechanism of the B. pseudocatenulatum genomes isolated from the intestinal tracts of different populations. Ultimately, based on the phylogenetic tree, we found that B. pseudocatenulatum has high genetic diversity and regional clustering phenomenon, in which plant-derived carbohydrate metabolism genes (GH13, GH43, GH5) showed significant regional differences, and this genetic differentiation drove the adaptive evolution, which likely shaped by diet and lifestyle.

Isolation of a Virulent Clostridium perfringens Strain from Elaphurus davidianus and Characterization by Whole-Genome Sequence Analysis

Clostridium perfringens (C. perfringens) is an important veterinary pathogen and a noteworthy threat to human and animal health. Recently, there has been a significant rise in the number of moose fatalities caused by this rare, endemic species in China. Currently, there is an increasing trend in conducting whole-genome analysis of C. perfringens strains originating from pigs and chickens, whereas fewer studies have been undertaken on Elaphurus davidianus-originating strains at the whole-genome level. Our laboratory has identified and isolated five C. perfringens type A from affected Elaphurus davidianus. The current study identified the most potent strain of C. perfringens, which originated from Elaphurus davidianus, and sequenced its genome to reveal virulence genes and pathogenicity. Our findings show that strain CX1-4 exhibits the highest levels of phospholipase activity, hemolytic activity, and mouse toxicity compared to the other four isolated C. perfringens type A strains. The chromosome sequence length of the CX1-4 strain was found to be 3,355,389 bp by complete genome sequencing. The current study unveils the genomic characteristics of C. perfringens type A originating from Elaphurus davidianus. It provides a core foundation for further investigation regarding the prevention and treatment of such infectious diseases in Elaphurus davidianus.

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.

Wolbachia symbionts control sex in a parasitoid wasp using a horizontally acquired gene

Host reproduction can be manipulated by bacterial symbionts in various ways. Parthenogenesis induction is the most effective type of reproduction manipulation by symbionts for their transmission. Insect sex is determined by regulation of doublesex (dsx) splicing through transformer2 (tra2) and transformer (tra) interaction. Although parthenogenesis induction by symbionts has been studied since the 1970s, its underlying molecular mechanism is unknown. Here we identify a Wolbachia parthenogenesis-induction feminization factor gene (piff) that targets sex-determining genes and causes female-producing parthenogenesis in the haplodiploid parasitoid Encarsia formosa. We found that Wolbachia elimination repressed expression of female-specific dsx and enhanced expression of male-specific dsx, which led to the production of wasp haploid male offspring. Furthermore, we found that E. formosa tra is truncated and non-functional, and Wolbachia has a functional tra homolog, termed piff, with an insect origin. Wolbachia PIFF can colocalize and interact with wasp TRA2. Moreover, Wolbachia piff has coordinated expression with tra2 and dsx of E. formosa. Our results demonstrate the bacterial symbiont Wolbachia has acquired an insect gene to manipulate the host sex determination cascade and induce parthenogenesis in wasps. This study reveals insect-to-bacteria horizontal gene transfer drives the evolution of animal sex determination systems, elucidating a striking mechanism of insect-microbe symbiosis.

Investigating the nature of prokaryotic genomic island locations within a genome

Horizontal gene transfer (HGT) is a powerful evolutionary force that considerably shapes the structure of prokaryotic genomes and is associated with genomic islands (GIs). A GI is a DNA segment composed of transferred genes that can be found within a prokaryotic genome, obtained through HGT. Much research has focused on detecting GIs in genomes, but here we pursue a new course, which is identifying possible preferred locations of GIs in the prokaryotic genome. Here, we identify the locations of the GIs within prokaryotic genomes to examine patterns in those locations. Prokaryotic GIs were analyzed according to the genome structure that they are located in, whether it be a circular or a linear genome. The analytical investigations employed are: (1) studying the GI locations in relation to the origin of replication (oriC); (2) exploring the distances between GIs; and (3) determining the distribution of GIs across the genomes. For each of the investigations, the analysis was performed on all of the GIs in the data set. Moreover, to void bias caused by the distribution of the genomes represented, the GIs in one genome from each species and the GIs of the most frequent species are also analyzed. Overall, the results showed that there are preferred sites for the GIs in the genome. In the linear genomes, these sites are usually located in the oriC region and terminus region, while in the circular genomes, they are located solely in the terminus region. These results also showed that the distance distribution between the GIs is almost exponential, which proves that GIs have preferred sites within genomes. The oriC and termniuns are preferred sites for the GIs and a possible natural explanation for this could be connected to the content of the oriC region. Moreover, the content of the GIs in terms of its protein families was studied and the results demonstrated that the majority of frequent protein families are close to identical in each section.

Genomic and environmental controls on Castellaniella biogeography in an anthropogenically disturbed subsurface

Castellaniella species have been isolated from a variety of mixed-waste environments including the nitrate and multiple metal-contaminated subsurface at the Oak Ridge Reservation (ORR). Previous studies examining microbial community composition and nitrate removal at ORR during biostimulation efforts reported increased abundances of members of the Castellaniella genus concurrent with increased denitrification rates. Thus, we asked how genomic and abiotic factors control the Castellaniella biogeography at the site to understand how these factors may influence nitrate transformation in an anthropogenically impacted setting. We report the isolation and characterization of several Castellaniella strains from the ORR subsurface. Five of these isolates match at 100% identity (at the 16S rRNA gene V4 region) to two Castellaniella amplicon sequence variants (ASVs), ASV1 and ASV2, that have persisted in the ORR subsurface for at least 2 decades. However, ASV2 has consistently higher relative abundance in samples taken from the site and was also the dominant blooming denitrifier population during a prior biostimulation effort. We found that the ASV2 representative strain has greater resistance to mixed metal stress than the ASV1 representative strains. We attribute this resistance, in part, to the large number of unique heavy metal resistance genes identified on a genomic island in the ASV2 representative genome. Additionally, we suggest that the relatively lower fitness of ASV1 may be connected to the loss of the nitrous oxide reductase (nos) operon (and associated nitrous oxide reductase activity) due to the insertion at this genomic locus of a mobile genetic element carrying copper resistance genes. This study demonstrates the value of integrating genomic, environmental, and phenotypic data to characterize the biogeography of key microorganisms in contaminated sites.

Comparative genomics and host range analysis of four Ralstonia pseudosolanacearum strains isolated from sunflower reveals genomic and phenotypic differences

BACKGROUND

Bacterial wilt caused by Ralstonia solanacearum species complex (RSSC) is one of the devastating diseases in crop production, seriously reducing the yield of crops. R. pseudosolanacearum, is known for its broad infrasubspecific diversity and comprises 36 sequevars that are currently known. Previous studies found that R. pseudosolanacearum contained four sequevars (13, 14, 17 and 54) isolated from sunflowers sown in the same field.

RESULTS

Here, we provided the complete genomes and the results of genome comparison of the four sequevars strains (RS639, RS642, RS647, and RS650). Four strains showed different pathogenicities to the same cultivars and different host ranges. Their genome sizes were about 5.84 ~ 5.94 Mb, encoding 5002 ~ 5079 genes and the average G + C content of 66.85% ~ 67%. Among the coding genes, 146 ~ 159 specific gene families (contained 150 ~ 160 genes) were found in the chromosomes and 34 ~ 77 specific gene families (contained 34 ~ 78 genes) in the megaplasmids from four strains. The average nucleotide identify (ANI) values between any two strains ranged from 99.05% ~ 99.71%, and the proportion of the total base length of collinear blocks accounts for the total gene length of corresponding genome was all more than 93.82%. Then, we performed a search for genomic islands, prophage sequences, the gene clusters macromolecular secretion systems, type III secreted effectors and other virulence factors in these strains, which provided detailed comparison results of their presence and distinctive features compared to the reference strain GMI1000. Among them, the number and types of T2SS gene clusters were different in the four strains, among which RS650 included all five types. T4SS gene cluster of RS639 and RS647 were missed. In the T6SS gene cluster, several genes were inserted in the RS639, RS647, and RS650, and gene deletion was also detected in the RS642. A total of 78 kinds of type III secreted effectors were found, which included 52 core and 9 specific effectors in four strains.

CONCLUSION

This study not only provided the complete genomes of multiple R. pseudosolanacearum strains isolated from a new host, but also revealed the differences in their genomic levels through comparative genomics. Furthermore, these findings expand human knowledge about the range of hosts that Ralstonia can infect, and potentially contribute to exploring rules and factors of the genetic evolution and analyzing its pathogenic mechanism.

Role of vertical and horizontal microbial transmission of antimicrobial resistance genes in early life: insights from maternal-infant dyads

Early life represents a critical window for metabolic, cognitive and immune system development, which is influenced by the maternal microbiome as well as the infant gut microbiome. Antibiotic exposure, mode of delivery and breastfeeding practices modulate the gut microbiome and the reservoir of antibiotic resistance genes (ARGs). Vertical and horizontal microbial gene transfer during early life and the mechanisms behind these transfers are being uncovered. In this review, we aim to provide an overview of the current knowledge on the transfer of antibiotic resistance in the mother-infant dyad through vertical and horizontal transmission and to highlight the main gaps and challenges in this area.

Biosynthetic gene cluster profiling from North Java Sea Virgibacillus salarius reveals hidden potential metabolites

Virgibacillus salarius 19.PP.SC1.6 is a coral symbiont isolated from Indonesia's North Java Sea; it has the ability to produce secondary metabolites that provide survival advantages and biological functions, such as ectoine, which is synthesized by an ectoine gene cluster. Apart from being an osmoprotectant for bacteria, ectoine is also known as a chemical chaperone with numerous biological activities such as maintaining protein stability, which makes ectoine in high demand in the market industry and makes it beneficial to investigate V. salarius ectoine. However, there has been no research on genome-based secondary metabolite and ectoine gene cluster characterization from Indonesian marine V. salarius. In this study, we performed a genomic analysis and ectoine identification of V. salarius. A high-quality draft genome with total size of 4.45 Mb and 4426 coding sequence (CDS) was characterized and then mapped into the Cluster of Orthologous Groups (COG) category. The genus Virgibacillus has an "open" pangenome type with total of 18 genomic islands inside the V. salarius 19.PP.SC1.6 genome. There were seven clusters of secondary metabolite-producing genes found, with a total of 80 genes classified as NRPS, PKS (type III), terpenes, and ectoine biosynthetic related genes. The ectoine gene cluster forms one operon consists of ectABC gene with 2190 bp gene cluster length, and is successfully characterized. The presence of ectoine in V. salarius was confirmed using UPLC-MS/MS operated in Multiple Reaction Monitoring (MRM) mode, which indicates that V. salarius has an intact ectoine gene clusters and is capable of producing ectoine as compatible solutes.

Novel approach toward the understanding of genetic diversity based on the two types of amino acid repeats in Erwinia amylovora

Erwinia amylovora is a notorious plant pathogenic bacterium of global concern that has devastated the apple and pear production industry worldwide. Nevertheless, the approaches available currently to understand the genetic diversity of E. amylovora remain unsatisfactory because of the lack of a trustworthy index and data covering the globally occurring E. amylovora strains; thus, their origin and distribution pattern remains ambiguous. Therefore, there is a growing need for robust approaches for obtaining this information via the comparison of the genomic structure of Amygdaloideae-infecting strains to understand their genetic diversity and distribution. Here, the whole-genome sequences of 245 E. amylovora strains available from the NCBI database were compared to identify intraspecific genes for use as an improved index for the simple classification of E. amylovora strains regarding their distribution. Finally, we discovered two kinds of strain-typing protein-encoding genes, i.e., the SAM-dependent methyltransferase and electron transport complex subunit RsxC. Interestingly, both of these proteins carried an amino acid repeat in these strains: SAM-dependent methyltransferase comprised a single-amino-acid repeat (asparagine), whereas RsxC carried a 40-amino-acid repeat, which was differentially distributed among the strains. These noteworthy findings and approaches may enable the exploration of the genetic diversity of E. amylovora from a global perspective.

A systematic approach to classify and characterize genomic islands driven by conjugative mobility using protein signatures

Genomic islands (GIs) play a crucial role in the spread of antibiotic resistance, virulence factors and antiviral defense systems in a broad range of bacterial species. However, the characterization and classification of GIs are challenging due to their relatively small size and considerable genetic diversity. Predicting their intercellular mobility is of utmost importance in the context of the emerging crisis of multidrug resistance. Here, we propose a large-scale classification method to categorize GIs according to their mobility profile and, subsequently, analyze their gene cargo. We based our classification decision scheme on a collection of mobility protein motif definitions available in publicly accessible databases. Our results show that the size distribution of GI classes correlates with their respective structure and complexity. Self-transmissible GIs are usually the largest, except in Bacillota and Actinomycetota, accumulate antibiotic and phage resistance genes, and favour the use of a tyrosine recombinase to insert into a host's replicon. Non-mobilizable GIs tend to use a DDE transposase instead. Finally, although tRNA genes are more frequently targeted as insertion sites by GIs encoding a tyrosine recombinase, most GIs insert in a protein-encoding gene. This study is a stepping stone toward a better characterization of mobile GIs in bacterial genomes and their mechanism of mobility.

Sequencing and characterization of Helcococcus ovis: a comprehensive comparative genomic analysis of virulence

BACKGROUND

Helcococcus ovis (H. ovis) is an emerging bacterial pathogen that commonly causes opportunistic respiratory, mammary, and uterine infections across mammalian hosts. This study applied long- and short-read whole genome sequencing technologies to identify virulence factors in five H. ovis isolates with low, medium, and high virulence phenotypes.

RESULTS

The resulting assemblies contained one circular chromosome ranging from 1,744,566 to 1,850,083 bp in length and had a mean GC content of 27.6%. Phylogenetic and nucleotide identity analyses found low virulence strain KG38 to be part of a clade that forms an outgroup apart from the rest of the H. ovis taxon. Assembling the first complete genomes of the species revealed major genomic rearrangements in KG38. One to six prophage regions were identified in each genome. A novel pathogenicity island was found exclusively in the two high virulence strains (KG37 and KG104), along with two hypothetical transmembrane proteins designated as putative VFs. Finally, three zinc ABC transporters and three Type-II/IV secretion systems were identified as possible virulence determinants in this species. The low virulence strain KG38 has fewer intact paralogs of these operons in its genome compared to the higher virulence isolates, which strongly suggests a role in virulence. This strain is also missing four putative virulence factors (VFs) found in other isolates associated with adherence (collagen adhesin precursor), immune evasion (choline-binding protein A and a PspA-like hypothetical protein) and cell wall synthesis (glycerol-3-phosphate cytidylyltransferase).

CONCLUSIONS

In this study, we assembled reference-quality complete genomes for five H. ovis strains to identify putative virulence factors. Phylogenetic analyses of H. ovis isolates revealed the presence of a clade representing a potentially novel species within the genus Helcococcus. A novel pathogenicity island and two hypothetical transmembrane proteins were found exclusively in high-virulence strains. The identification of Zinc ABC transporters and Type-II/IV secretion systems as possible virulence determinants, along with the differences in operon content between the low and high virulence isolates, strongly suggests they also play a role in the bacterium's pathogenicity. Taken together, these findings are a valuable first step toward deciphering the pathogenesis of H. ovis infections.

Distinctive microbial community and genome structure in coastal seawater from a human-made port and nearby offshore island in northern Taiwan facing the Northwestern Pacific Ocean

Pollution in human-made fishing ports caused by petroleum from boats, dead fish, toxic chemicals, and effluent poses a challenge to the organisms in seawater. To decipher the impact of pollution on the microbiome, we collected surface water from a fishing port and a nearby offshore island in northern Taiwan facing the Northwestern Pacific Ocean. By employing 16S rRNA gene amplicon sequencing and whole-genome shotgun sequencing, we discovered that Rhodobacteraceae, Vibrionaceae, and Oceanospirillaceae emerged as the dominant species in the fishing port, where we found many genes harboring the functions of antibiotic resistance (ansamycin, nitroimidazole, and aminocoumarin), metal tolerance (copper, chromium, iron and multimetal), virulence factors (chemotaxis, flagella, T3SS1), carbohydrate metabolism (biofilm formation and remodeling of bacterial cell walls), nitrogen metabolism (denitrification, N2 fixation, and ammonium assimilation), and ABC transporters (phosphate, lipopolysaccharide, and branched-chain amino acids). The dominant bacteria at the nearby offshore island (Alteromonadaceae, Cryomorphaceae, Flavobacteriaceae, Litoricolaceae, and Rhodobacteraceae) were partly similar to those in the South China Sea and the East China Sea. Furthermore, we inferred that the microbial community network of the cooccurrence of dominant bacteria on the offshore island was connected to dominant bacteria in the fishing port by mutual exclusion. By examining the assembled microbial genomes collected from the coastal seawater of the fishing port, we revealed four genomic islands containing large gene-containing sequences, including phage integrase, DNA invertase, restriction enzyme, DNA gyrase inhibitor, and antitoxin HigA-1. In this study, we provided clues for the possibility of genomic islands as the units of horizontal transfer and as the tools of microbes for facilitating adaptation in a human-made port environment.

A Novel CovS Variant Harbored by a Colonization Strain Reduces Streptococcus pyogenes Virulence

Streptococcus pyogenes, also known as group A Streptococcus, causes a wide variety of diseases ranging from mild noninvasive to severe invasive infections. To identify possible causes of colonization-to-invasive switches, we determined the genomic sequences of 10 isolates from five pairs each composed of an invasive strain and a carriage strain originating from five infectious clusters. Among them, one pair displayed a single-nucleotide difference in covS, encoding the sensor histidine kinase of the two-component CovRS system that controls the expression of 15% of the genome. In contrast to previously described cases where the invasive strains harbor nonfunctional CovS proteins, the carriage strain possessed the mutation covST115C, leading to the replacement of the tyrosine at position 39 by a histidine. The CovSY39H mutation affected the expression of the genes from the CovR regulon in a unique fashion. Genes usually overexpressed in covS mutant strains were underexpressed and vice versa. Furthermore, the covS mutant strain barely responded to the addition of the CovS-signaling compounds Mg2+ and LL-37. The variations in the accumulation of two virulence factors paralleled the transcription modifications. In addition, the covST115C mutant strain showed less survival than its wild-type counterpart in murine macrophages. Finally, in two murine models of infection, the covS mutant strain was less virulent than the wild-type strain. Our study suggests that the CovSY39H protein compromises CovS phosphatase activity and that this yields a noninvasive strain. IMPORTANCE Streptococcus pyogenes, also known as group A Streptococcus, causes a wide variety of diseases, leading to 517,000 deaths yearly. The two-component CovRS system, which responds to MgCl2 and the antimicrobial peptide LL-37, controls the expression of 15% of the genome. Invasive strains may harbor nonfunctional CovS sensor proteins that lead to the derepression of most virulence genes. We isolated a colonization strain that harbors a novel covS mutation. This mutant strain harbored a transcriptome profile opposite that of other covS mutant strains, barely responded to environmental signals, and was less virulent than the wild-type strain. This supports the importance of the derepression of the expression of most virulence genes, via mutations that impact the phosphorylation of the regulator CovR, for favoring S. pyogenes invasive infections.

Genomic Islands Identified in Highly Resistant Serratia sp. HRI: A Pathway to Discover New Disinfectant Resistance Elements

Molecular insights into the mechanisms of resistance to disinfectants are severely limited, together with the roles of various mobile genetic elements. Genomic islands are a well-characterised molecular resistance element in antibiotic resistance, but it is unknown whether genomic islands play a role in disinfectant resistance. Through whole-genome sequencing and the bioinformatic analysis of Serratia sp. HRI, an isolate with high disinfectant resistance capabilities, nine resistance islands were predicted and annotated within the genome. Resistance genes active against several antimicrobials were annotated in these islands, most of which are multidrug efflux pumps belonging to the MFS, ABC and DMT efflux families. Antibiotic resistance islands containing genes encoding for multidrug resistance proteins ErmB (macrolide and erythromycin resistance) and biclomycin were also found. A metal fitness island harbouring 13 resistance and response genes to copper, silver, lead, cadmium, zinc, and mercury was identified. In the search for disinfectant resistance islands, two genomic islands were identified to harbour smr genes, notorious for conferring disinfectant resistance. This suggests that genomic islands are capable of conferring disinfectant resistance, a phenomenon that has not yet been observed in the study of biocide resistance and tolerance.

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.

The complete and closed genome of the facultative generalist Candidatus Endoriftia persephone from deep-sea hydrothermal vents

The mutualistic interactions between Riftia pachyptila and its endosymbiont Candidatus Endoriftia persephone (short Endoriftia) have been extensively researched. However, the closed Endoriftia genome is still lacking. Here, by employing single-molecule real-time sequencing we present the closed chromosomal sequence of Endoriftia. In contrast to theoretical predictions of enlarged and mobile genetic element-rich genomes related to facultative endosymbionts, the closed Endoriftia genome is streamlined with fewer than expected coding sequence regions, insertion-, prophage-sequences and transposase-coding sequences. Automated and manually curated functional analyses indicated that Endoriftia is more versatile regarding sulphur metabolism than previously reported. We identified the presence of two identical rRNA operons and two long CRISPR regions in the closed genome. Additionally, pangenome analyses revealed the presence of three types of secretion systems (II, IV and VI) in the different Endoriftia populations indicating lineage-specific adaptations. The in depth mobilome characterization identified the presence of shared genomic islands in the different Endoriftia drafts and in the closed genome, suggesting that the acquisition of foreign DNA predates the geographical dispersal of the different endosymbiont populations. Finally, we found no evidence of epigenetic regulation in Endoriftia, as revealed by gene screenings and absence of methylated modified base motifs in the genome. As a matter of fact, the restriction-modification system seems to be dysfunctional in Endoriftia, pointing to a higher importance of molecular memory-based immunity against phages via spacer incorporation into CRISPR system. The Endoriftia genome is the first closed tubeworm endosymbiont to date and will be valuable for future gene oriented and evolutionary comparative studies.

Reduced and Nonreduced Genomes in Paraburkholderia Symbionts of Social Amoebas

The social amoeba Dictyostelium discoideum is a predatory soil protist frequently used for studying host-pathogen interactions. A subset of D. discoideum strains isolated from soil persistently carry symbiotic Paraburkholderia, recently formally described as P. agricolaris, P. bonniea, and P. hayleyella. The three facultative symbiont species of D. discoideum present a unique opportunity to study a naturally occurring symbiosis in a laboratory model protist. There is a large difference in genome size between P. agricolaris (8.7 million base pairs [Mbp]) versus P. hayleyella and P. bonniea (4.1 Mbp). We took a comparative genomics approach and compared the three genomes of D. discoideum symbionts to 12 additional Paraburkholderia genomes to test for genome evolution patterns that frequently accompany host adaptation. Overall, P. agricolaris is difficult to distinguish from other Paraburkholderia based on its genome size and content, but the reduced genomes of P. bonniea and P. hayleyella display characteristics indicative of genome streamlining rather than deterioration during adaptation to their protist hosts. In addition, D. discoideum-symbiont genomes have increased secretion system and motility genes that may mediate interactions with their host. Specifically, adjacent BurBor-like type 3 and T6SS-5-like type 6 secretion system operons shared among all three D. discoideum-symbiont genomes may be important for host interaction. Horizontal transfer of these secretion system operons within the amoeba host environment may have contributed to the unique ability of these symbionts to establish and maintain a symbiotic relationship with D. discoideum. IMPORTANCE Protists are a diverse group of typically single cell eukaryotes. Bacteria and archaea that form long-term symbiotic relationships with protists may evolve in additional ways than those in relationships with multicellular eukaryotes such as plants, animals, or fungi. Social amoebas are a predatory soil protist sometimes found with symbiotic bacteria living inside their cells. They present a unique opportunity to explore a naturally occurring symbiosis in a protist frequently used for studying host-pathogen interactions. We show that one amoeba-symbiont species is similar to other related bacteria in genome size and content, while the two reduced-genome-symbiont species show characteristics of genome streamlining rather than deterioration during adaptation to their host. We also identify sets of genes present in all three amoeba-symbiont genomes that are potentially used for host-symbiont interactions. Because the amoeba symbionts are distantly related, the amoeba host environment may be where these genes were shared among symbionts.

Sequence Analysis of Novel Staphylococcus aureus Lineages from Wild and Captive Macaques

Staphylococcus aureus is a widespread and common opportunistic bacterium that can colonise or infect humans as well as a wide range of animals. There are a few studies of both methicillin-susceptible S. aureus (MSSA) and methicillin-resistant S. aureus (MRSA) isolated from monkeys, apes, and lemurs, indicating a presence of a number of poorly or unknown lineages of the pathogen. In order to obtain insight into staphylococcal diversity, we sequenced strains from wild and captive individuals of three macaque species (Macaca mulatta, M. assamensis, and M. sylvanus) using Nanopore and Illumina technologies. These strains were previously identified by microarray as poorly or unknown strains. Isolates of novel lineages ST4168, ST7687, ST7688, ST7689, ST7690, ST7691, ST7692, ST7693, ST7694, ST7695, ST7745, ST7746, ST7747, ST7748, ST7749, ST7750, ST7751, ST7752, ST7753, and ST7754 were sequenced and characterised for the first time. In addition, isolates belonging to ST2990, a lineage also observed in humans, and ST3268, a MRSA strain already known from macaques, were also included into the study. Mobile genetic elements, genomic islands, and carriage of prophages were analysed. There was no evidence for novel host-specific virulence factors. However, a conspicuously high rate of carriage of a pathogenicity island harbouring edinB and etD2/etE as well as a higher number of repeat units within the gene sasG (encoding an adhesion factor) than in human isolates were observed. None of the strains harboured the genes encoding Panton-Valentine leukocidin. In conclusion, wildlife including macaques may harbour an unappreciated diversity of S. aureus lineages that may be of clinical relevance for humans, livestock, or for wildlife conservation, given the declining state of many wildlife populations.

Genome-wide identification of type III effectors and other virulence factors in Ralstonia pseudosolanacearum causing bacterial wilt in ginger (Zingiber officinale)

Ralstonia pseudosolanacearum causes bacterial wilt in ginger, reducing ginger production worldwide. We sequenced the whole genome of a highly virulent phylotype I, race 4, biovar 3 Ralstonia pseudosolanacearum strain GRsMep isolated from a severely infected ginger field in India. R. pseudosolanacearum GRsMep genome is organised into two replicons: chromosome and megaplasmid with a total genome size of 5,810,605 bp. This strain encodes approximately 72 effectors which include a combination of core effectors as well as highly variable, diverse repertoire of type III effectors. Comparative genome analysis with GMI1000 identified conservation in the genes involved in the general virulence mechanism. Our analysis identified type III effectors, RipBJ and RipBO as present in GRsMep but absent in the reported genomes of other strains infecting Zingiberaceae family. GRsMep contains 126 unique genes when compared to the pangenome of the Ralstonia strains that infect the Zingiberaceae family. The whole-genome data of R. pseudosolanacearum strain will serve as a resource for exploring the evolutionary processes that structure and regulate the virulence determinants of the strain. Pathogenicity testing of the transposon insertional mutant library of GRsMep through virulence assay on ginger plants identified a few candidate virulence determinants specific to bacterial wilt in ginger.

De Novo Arginine Synthesis Is Required for Full Virulence of Xanthomonas arboricola pv. juglandis During Walnut Bacterial Blight Disease

Walnut blight (WB) disease caused by Xanthomonas arboricola pv. juglandis (Xaj) threatens orchards worldwide. Nitrogen metabolism in this bacterial pathogen is dependent on arginine, a nitrogen-enriched amino acid that can either be synthesized or provided by the plant host. The arginine biosynthetic pathway uses argininosuccinate synthase (argG), associated with increased bacterial virulence. We examined the effects of bacterial arginine and nitrogen metabolism on the plant response during WB by proteomic analysis of the mutant strain Xaj argG^-. Phenotypically, the mutant strain produced 42% fewer symptoms and survived in the plant tissue with 2.5-fold reduced growth compared with wild type, while showing itself to be auxotrophic for arginine in vitro. Proteomic analysis of infected tissue enabled the profiling of 676 Xaj proteins and 3,296 walnut proteins using isobaric labeling in a data-dependent acquisition approach. Comparative analysis of differentially expressed proteins revealed distinct plant responses. Xaj wild type (WT) triggered processes of catabolism and oxidative stress in the host under observed disease symptoms, while most of the host biosynthetic processes triggered by Xaj WT were inhibited during Xaj argG^- infection. Overall, the Xaj proteins revealed a drastic shift in carbon and energy management induced by disruption of nitrogen metabolism while the top differentially expressed proteins included a Fis transcriptional regulator and a peptidyl-prolyl isomerase. Our results show the critical role of de novo arginine biosynthesis to sustain virulence and minimal growth during WB. This study is timely and critical as copper-based control methods are losing their effectiveness, and new sustainable methods are urgently needed in orchard environments.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Comparative genomic analysis of Canadian Mycoplasma bovis strains isolated from Bison and Cattle

Mycoplasma bovis (M. bovis) in cattle causes pneumonia, arthritis, otitis media, and mastitis. In addition, multiple outbreaks have been recorded in North American bison. The genomic data on Canadian M. bovis in bison and cattle to date is limited. Whole-genome sequencing (WGS) was used to assess the degree of genome conservation across four Canadian M. bovis strains recovered from bison and cattle. Whole-genome sequences of four M. bovis isolates (Mb1, Mb160, Mb300, Mb304) and the PG45 reference genome were utilized to identify the M. bovis genomic similarity, whole-genome single nucleotide polymorphism (WGS-SNP), virulence determinants, and genomic islands. The pan-genome analysis showed that M. bovis encodes a minimum of 971 genes, while the core genome contained 637 genes. Comparative genomics revealed limited diversity in gene content between bison and cattle isolates. Whole-genome SNP analysis showed that the four M. bovis isolates differed from each other and to PG45. A total of 40 putative virulence genes associated with adhesion, colonization, and destruction of tissues were found in the bison and cattle isolates using the virulence factors database (VFDB). These putative virulence factors were equally distributed among isolates. Genomic Islands (GIs) ranging from 4 to 9 and associated with transposases, restriction-modification, ribosomal hypothetical proteins, variable surface lipoproteins, and unknowns were also identified. Overall, the genomic characterization of these isolates may provide new insights into the mechanisms of pathogenicity in M. bovis.

Enabling genomic island prediction and comparison in multiple genomes to investigate bacterial evolution and outbreaks

Outbreaks of virulent and/or drug-resistant bacteria have a significant impact on human health and major economic consequences. Genomic islands (GIs; defined as clusters of genes of probable horizontal origin) are of high interest because they disproportionately encode virulence factors, some antimicrobial-resistance (AMR) genes, and other adaptations of medical or environmental interest. While microbial genome sequencing has become rapid and inexpensive, current computational methods for GI analysis are not amenable for rapid, accurate, user-friendly and scalable comparative analysis of sets of related genomes. To help fill this gap, we have developed IslandCompare, an open-source computational pipeline for GI prediction and comparison across several to hundreds of bacterial genomes. A dynamic and interactive visualization strategy displays a bacterial core-genome phylogeny, with bacterial genomes linearly displayed at the phylogenetic tree leaves. Genomes are overlaid with GI predictions and AMR determinants from the Comprehensive Antibiotic Resistance Database (CARD), and regions of similarity between the genomes are also displayed. GI predictions are performed using Sigi-HMM and IslandPath-DIMOB, the two most precise GI prediction tools based on nucleotide composition biases, as well as a novel blast-based consistency step to improve cross-genome prediction consistency. GIs across genomes sharing sequence similarity are grouped into clusters, further aiding comparative analysis and visualization of acquisition and loss of mobile GIs in specific sub-clades. IslandCompare is an open-source software that is containerized for local use, plus available via a user-friendly, web-based interface to allow direct use by bioinformaticians, biologists and clinicians (at https://islandcompare.ca).

Enhanced production of polyhydroxyalkanoates in Pseudomonas putida KT2440 by a combination of genome streamlining and promoter engineering

Polyhydroxyalkanoates (PHAs), a class of bioplastics produced by a variety of microorganisms, have become the ideal alternatives for oil-derived plastics due to their superior physicochemical and material characteristics. Pseudomonas putida KT2440 can produce medium-chain-length PHA (mcl-PHA) from various substrates. In this study, a novel strategy of the large-scale deletion of genomic islands (GIs) coupling with promoter engineering was developed in P. putida KT2440 for constructing the minimal genome cell factories (MGF) capable of efficiently producing mcl-PHA. Firstly, P. putida KTU-U13, a 13 GIs- and upp-deleted mutant derived from the parental strain P. putida KT2440, was used as a starting strain for further deletion of GIs to generate a series of genome-reduced strains. Subsequently, the two minimal genome strains KTU-U24 and KTU-U27, which had a 7.19% and 8.35% reduction relative to the genome size of KT2440 and were advantageous over the strain KTU (KT2440∆upp) and KTU-U13 in several physiological traits such as the maximum specific growth rate, plasmid transformation efficiency, heterologous protein expression capacity and PHA production capacity, were selected as the chassis cells for PHA metabolic engineering. To prevent the formation of the by-product gluconic acid, the glucose dehydrogenase gene was deleted in KTU-U24 and KTU-U27, resulting in KTU-U24∆gcd and KTU-U27∆gcd. To enhance the transcriptional level of PHA synthase genes (phaC) and the supply of the precursor acetyl-CoA, a strong endogenous promoter P46 was inserted into upstream of the phaC operon and pyruvate dehydrogenase gene in the genome of KTU-U24∆gcd and KTU-U27∆gcd, to generate KTU-U24∆gcd-P46CA and KTU-U27∆gcd-P46CA, with the PHA yield of 50.5 wt% and 53.8 wt% (weight percent of PHA in cell dry weight). Finally, KTU-U27∆gcd-P46CA, the most minimal KT2440 chassis currently available, was able to accumulate the PHA to 55.82 wt% in a 5-l fermentor, which is the highest PHA yield obtained with P. putida KT2440 so far. This study suggests that genome streamlining in combination with promoter engineering may be a feasible strategy for the development of the MGF for the efficient production of high value products. Moreover, further streamlining of the P. putida KT2440 genome has great potential to create the optimal chassis for synthetic biology applications.

Complete genome sequences of the antibiotic sulfamethoxazole-mineralizing bacteria Paenarthrobacter sp. P27 and Norcardiodes sp. N27

Sulfonamide antibiotics (SAs) have been produced and consumed on a large scale over the last few decades. SAs are a typical class of refractory contaminants that are omnipresent in various environments. Although several [phenyl]-SA-degrading bacteria and their corresponding genomes have been documented, limited genetic information is available for the degraders of heterocyclic products (e.g., 3-amino-5-methylisoxazole [3A5MI] produced via sulfamethoxazole [SMX] catabolism). In this study, the previously isolated SMX-mineralizing bacterial partners, Paenarthrobacter sp. P27 (responsible for the initial cleavage of the -C-S-N- bond of SMX and further degradation of [phenyl]-SMX) and Norcardiodes sp. N27 (responsible for 3A5MI catabolism), were further studied and their complete genomes were sequenced. Complete degradation and bacterial growth were verified by pure-culture experiments with SMX or 3A5MI as the sole carbon, nitrogen, and energy source. By cross-feeding strains P27 and N27, complete catabolism of SMX could be achieved over a wide range of initial SMX concentrations. Moreover, strain P27 was capable of transforming the additional nine SA representatives into their corresponding nitrogen-containing heterocyclic products, strongly indicating the broad substrate spectrum and marked bioremediation potential of strain P27. The genome of strain P27 contained the highly homologous monooxygenase gene cluster, sadABC, which initially attacked the sulfonamide molecules. The complete genome sequences of the two important degraders will benefit future research centering on the molecular mechanism underlying advanced SMX mineralization and will aid in further understanding the interspecific interactions and metabolite exchanges for the optimization of artificially constructed synthetic functional microbiomes.

Distribution and Evolutionary History of Sialic Acid Catabolism in the Phylum Actinobacteria

Sialic acids are present in humans and other metazoans, playing essential roles in physiological and pathological processes. Commensal and pathogenic bacteria have evolved the capacity to utilize sialic acids as nutrient and energy sources. However, in some actinobacteria, sialic acid catabolism (SAC) is associated with free-living populations. To unravel the distribution and evolutionary history of SAC in the phylum Actinobacteria, we analyzed the presence and diversity of the putative SAC gene cluster (nan) in 7,180 high-quality, nonredundant actinobacterial genomes that covered 1,969 species. The results showed that ∼13% of actinobacterial species had the potential to utilize sialic acids, with 45 species capable of anhydro-SAC, all except two of them through the canonical pathway. These species belonged to 20 orders and 81 genera, with ∼36% of them from four genera, Actinomyces, Bifidobacterium, Corynebacterium, and Streptomyces. Moreover, ∼40% of the nan-positive species are free living. Phylogenetic analysis of the key nan genes, nanA, nanK, and nanE, revealed a strong signal of horizontal gene transfer (HGT), accompanied with vertical inheritance and gene loss. This evolutionary pattern led to high diversity and differential distribution of nan among actinobacterial taxa and might cause the cluster to spread to some free-living species while losing in some host-associated species. The evolution of SAC in actinobacteria probably represents the evolution of certain kinds of noncore bacterial functions for environmental adaptation and lifestyle switch, in which HGT plays a dominant role. IMPORTANCE Sialic acids play essential roles in the physiology of humans and other metazoan animals, and microbial sialic acid catabolism (SAC) is one of the processes critical for pathogenesis. To date, microbial SAC is studied mainly in commensals and pathogens, while its distribution in free-living microbes and evolutionary pathway remain largely unexplored. Here, by examining all actinobacterial genomes available, we demonstrate that putative SAC is present in a small proportion of actinobacterial species, of which, however, ∼40% are free-living species. We also reveal remarkable difference in the distribution of SAC among actinobacterial taxa and high diversity of the putative SAC gene clusters. HGT plays a significant role in the evolution of SAC, accompanied with vertical inheritance and gene loss. Our results provide a comprehensive and systematic picture of the distribution and evolutionary history of SAC in actinobacteria, expanding the current knowledge on bacterial adaptation and diversification.

Genome sequence of a multidrug-resistant Campylobacter coli strain isolated from a newborn with severe diarrhea in Lebanon

A multidrug-resistant (MDR) Campylobacter coli (C. coli) strain was isolated from a 2-month-old newborn who suffered from severe diarrhea in Lebanon. Here, whole-genome sequencing (WGS) analysis was deployed to determine the genetic basis of antimicrobial resistance and virulence in the C. coli isolate and to identify its epidemiological background (sequence type). The identity of the isolate was confirmed using API® Campy, MALDI-TOF, and 16S rRNA gene sequencing analysis. The antimicrobial susceptibility phenotype was determined using the disk diffusion assay. Our analysis showed that resistance to macrolide and quinolone was potentially associated with the presence of multiple point mutations in antibiotic targets on the chromosomal DNA. Furthermore, tetracycline and aminoglycoside resistance were encoded by genes on a pTet plasmid. The bla_OXA-61, which is associated with beta-lactam resistance, was also detected in the C. coli genome. A set of 30 genes associated with the virulence in C. coli was detected using WGS analysis. MLST analysis classified the isolate as belonging to a new sequence type (ST-9588), a member of ST-828 complex which is mainly associated with humans and chickens. Taking together, this study provides the first WGS analysis of Campylobacter isolated from Lebanon. The detection of a variety of AMR and virulence determinants strongly emphasizes the need for studying the burden of Campylobacter in Lebanon and the Middle East and North Africa (MENA) region, where information on campylobacteriosis is scant.

Genomic Features and Construction of Streamlined Genome Chassis of Nisin Z Producer Lactococcus lactis N8

Lactococcus lactis is a commonly used fermenting bacteria in cheese, beverages and meat products. Due to the lack of simplified chassis strains, it has not been widely used in the fields of synthetic biology. Thus, the construction of lactic acid bacteria chassis strains becomes more and more important. In this study, we performed whole genome sequencing, annotation and analysis of L. lactis N8. Based on the genome analysis, we found that L. lactis N8 contains two large plasmids, and the function prediction of the plasmids shows that some regions are related to carbohydrate transport/metabolism, multi-stress resistance and amino acid uptake. L. lactis N8 contains a total of seven prophage-related fragments and twelve genomic islands. A gene cluster encoding a hybrid NRPS–PKS system that was found in L. lactis N8 reveals that the strain has the potential to synthesize novel secondary metabolites. Furthermore, we have constructed a simplified genome chassis of L. lactis N8 and achieved the largest amount of deletion of L. lactis so far. Taken together, the present study offers further insights into the function and potential role of L. lactis N8 as a model strain of lactic acid bacteria and lays the foundation for its application in the field of synthetic biology.

Genome-based characterization of a plasmid-associated micrococcin P1 biosynthetic gene cluster and virulence factors in Mammaliicoccus sciuri IMDO-S72

Analysis of the de novo assembled genome of Mammaliicoccus sciuri IMDO-S72 revealed the genetically encoded machinery behind its earlier reported antibacterial phenotype and gave further insight into the repertoire of putative virulence factors of this recently reclassified species. A plasmid-encoded biosynthetic gene cluster was held responsible for the antimicrobial activity of M. sciuri IMDO-S72, comprising genes involved in thiopeptide production. The compound encoded by this gene cluster was structurally identified as micrococcin P1. Further examination of its genome highlighted the ubiquitous presence of innate virulence factors mainly involved in surface colonization. Determinants contributing to aggressive virulence were generally absent, with exception of a plasmid-associated ica cluster. The native antibiotic resistance genes sal(A) and mecA were detected within the genome, amongst others, but were not consistently linked with a resistant phenotype. While mobile genetic elements were identified within the genome, such as an untypeable SCC element, they proved to be generally free of virulence- and antibiotic-related genes. These results further suggest a commensal lifestyle of M. sciuri and indicate the association of antibiotic resistance determinants with mobile genetic elements, as an important factor in conferring antibiotic resistance, in addition to their unilateral annotation. Importance Mammaliicoccus sciuri has been put forward as an important carrier of virulence and antibiotic resistance genes, which can be transmitted to clinically important staphylococcal species such as Staphylococcus aureus. As a common inhabitant of mammal skin, this species is believed to have a predominant commensal lifestyle although it has been reported as an opportunistic pathogen in some cases. This study provides an extensive genome-wide description of its putative virulence potential taking into consideration the genomic context in which these genes appear, an aspect that is often overlooked during virulence analysis. Additional genome and biochemical analysis linked M. sciuri with the production of micrococcin P1, gaining further insight to which extent these biosynthetic gene cluster are distributed amongst different related species. The frequent plasmid-associated character hints that these traits can be horizontally transferred and might confer a competitive advantage to its recipient within its ecological niche.

Evolutionary origin and ecological implication of a unique nif island in free-living Bradyrhizobium lineages

The alphaproteobacterial genus Bradyrhizobium has been best known as N₂-fixing members that nodulate legumes, supported by the nif and nod gene clusters. Recent environmental surveys show that Bradyrhizobium represents one of the most abundant free-living bacterial lineages in the world's soils. However, our understanding of Bradyrhizobium comes largely from symbiotic members, biasing the current knowledge of their ecology and evolution. Here, we report the genomes of 88 Bradyrhizobium strains derived from diverse soil samples, including both nif-carrying and non-nif-carrying free-living (nod free) members. Phylogenomic analyses of these and 252 publicly available Bradyrhizobium genomes indicate that nif-carrying free-living members independently evolved from symbiotic ancestors (carrying both nif and nod) multiple times. Intriguingly, the nif phylogeny shows that the vast majority of nif-carrying free-living members comprise an independent cluster, indicating that horizontal gene transfer promotes nif expansion among the free-living Bradyrhizobium. Comparative genomics analysis identifies that the nif genes found in free-living Bradyrhizobium are located on a unique genomic island of ~50 kb equipped with genes potentially involved in coping with oxygen tension. We further analyze amplicon sequencing data to show that Bradyrhizobium members presumably carrying this nif island are widespread in a variety of environments. Given the dominance of Bradyrhizobium in world's soils, our findings have implications for global nitrogen cycles and agricultural research.

Diversity in the composition of the accessory genome of Mexican Pseudomonas aeruginosa strains

BACKGROUND

Pseudomonas aeruginosa is an important opportunistic pathogen especially in nosocomial infections due to its easy adaptation to different environments; this characteristic is due to the great genetic diversity that presents its genome. In addition, it is considered a pathogen of critical priority due to the high antimicrobial resistance.

OBJECTIVES

The aim of this study was to characterize the mobile genetic elements present in the chromosome of six Mexican P. aeruginosa strains isolated from adults with pneumonia and children with bacteremia.

METHODS

The genomic DNA of six P. aeruginosa strains were isolated and sequenced using PacBio RS-II platform. They were annotated using Prokaryotic Genome Annotation Pipeline and manually curated and analyzed for the presence of mobile genetic elements, antibiotic resistances genes, efflux pumps and virulence factors using several bioinformatics programs and databases.

RESULTS

The global analysis of the strains chromosomes showed a novel chromosomal rearrangement in two strains, possibly mediated by subsequent recombination and inversion events. They have a high content of mobile genetic elements: 21 genomic islands, four new islets, four different integrative conjugative elements, 28 different prophages, one CRISPR-Cas arrangements, and one class 1 integron. The acquisition of antimicrobials resistance genes into these elements are in concordance with their phenotype of multi-drug resistance.

CONCLUSION

The accessory genome increased the ability of the strains to adapt or survive to the hospital environment, promote genomic plasticity and chromosomal rearrangements, which may affect the expression or functionality of the gene and might influence the clinical outcome, having an impact on the treatment.

Systematic Detection of Large-Scale Multigene Horizontal Transfer in Prokaryotes

Horizontal gene transfer (HGT) is central to prokaryotic evolution. However, little is known about the “scale” of individual HGT events. In this work, we introduce the first computational framework to help answer the following fundamental question: How often does more than one gene get horizontally transferred in a single HGT event? Our method, called HoMer, uses phylogenetic reconciliation to infer single-gene HGT events across a given set of species/strains, employs several techniques to account for inference error and uncertainty, combines that information with gene order information from extant genomes, and uses statistical analysis to identify candidate horizontal multigene transfers (HMGTs) in both extant and ancestral species/strains. HoMer is highly scalable and can be easily used to infer HMGTs across hundreds of genomes. We apply HoMer to a genome-scale data set of over 22,000 gene families from 103 Aeromonas genomes and identify a large number of plausible HMGTs of various scales at both small and large phylogenetic distances. Analysis of these HMGTs reveals interesting relationships between gene function, phylogenetic distance, and frequency of multigene transfer. Among other insights, we find that 1) the observed relative frequency of HMGT increases as divergence between genomes increases, 2) HMGTs often have conserved gene functions, and 3) rare genes are frequently acquired through HMGT. We also analyze in detail HMGTs involving the zonula occludens toxin and type III secretion systems. By enabling the systematic inference of HMGTs on a large scale, HoMer will facilitate a more accurate and more complete understanding of HGT and microbial evolution.

Taxonomy classification using genomic footprint of mitochondrial sequences

BACKGROUND

Advancement in the sequencing technology yields a huge number of genomes of a multitude of organisms in our planet. One of the fundamental tasks for processing and analyzing these sequences is to organize them in the existing taxonomic orders. <P> Method: Recently we proposed a novel approach, GenFooT, of taxonomy classification using the concept of genomic footprint (GFP). The technique is further refined and enhanced in this work leading to improved accuracies in the task of taxonomic classification on various benchmark datasets. GenFooT maps a genome sequence in a 2D coordinate space and extracts features from that representation. It uses two hyper-parameters, namely block size and number of fragments of genomic sequence while computing the feature. In this work, we propose an analysis for choosing values of those parameters adaptively from the sequences. The enhanced version of GenFooT is named GenFooT2. <P> Results and Conclusion: We have experimented GenFooT2 on ten different biological datasets of genomic sequences of various organisms belonging to different taxonomy ranks. Our experimental results indicate more than 3% improved classification performance of the proposed features with Logistic regression classifier than the GenFooT. We also performed the statistical test to compare the performance of GenFooT2 with the state-of-the-art methods including our previous method GenFooT. The experimental results as well as the statistical test exhibit that the performance of the proposed GenFooT2 is significantly better.

Genome Mining and Comparative Genome Analysis Revealed Niche-Specific Genome Expansion in Antibacterial Bacillus pumilus Strain SF-4

The present study reports the isolation of antibacterial exhibiting Bacillus pumilus (B. pumilus) SF-4 from soil field. The genome of this strain SF-4 was sequenced and analyzed to acquire in-depth genomic level insight related to functional diversity, evolutionary history, and biosynthetic potential. The genome of the strain SF-4 harbor 12 Biosynthetic Gene Clusters (BGCs) including four Non-ribosomal peptide synthetases (NRPSs), two terpenes, and one each of Type III polyketide synthases (PKSs), hybrid (NRPS/PKS), lipopeptide, β-lactone, and bacteriocin clusters. Plant growth-promoting genes associated with de-nitrification, iron acquisition, phosphate solubilization, and nitrogen metabolism were also observed in the genome. Furthermore, all the available complete genomes of B. pumilus strains were used to highlight species boundaries and diverse niche adaptation strategies. Phylogenetic analyses revealed local diversification and indicate that strain SF-4 is a sister group to SAFR-032 and 150a. Pan-genome analyses of 12 targeted strains showed regions of genome plasticity which regulate function of these strains and proposed direct strain adaptations to specific habitats. The unique genome pool carries genes mostly associated with “biosynthesis of secondary metabolites, transport, and catabolism” (Q), “replication, recombination and repair” (L), and “unknown function” (S) clusters of orthologous groups (COG) categories. Moreover, a total of 952 unique genes and 168 exclusively absent genes were prioritized across the 12 genomes. While newly sequenced B. pumilus SF-4 genome consists of 520 accessory, 59 unique, and seven exclusively absent genes. The current study demonstrates genomic differences among 12 B. pumilus strains and offers comprehensive knowledge of the respective genome architecture which may assist in the agronomic application of this strain in future.

Identifying genomic islands with deep neural networks

Background

Horizontal gene transfer is the main source of adaptability for bacteria, through which genes are obtained from different sources including bacteria, archaea, viruses, and eukaryotes. This process promotes the rapid spread of genetic information across lineages, typically in the form of clusters of genes referred to as genomic islands (GIs). Different types of GIs exist, and are often classified by the content of their cargo genes or their means of integration and mobility. While various computational methods have been devised to detect different types of GIs, no single method is capable of detecting all types.

Results

We propose a method, which we call Shutter Island, that uses a deep learning model (Inception V3, widely used in computer vision) to detect genomic islands. The intrinsic value of deep learning methods lies in their ability to generalize. Via a technique called transfer learning, the model is pre-trained on a large generic dataset and then re-trained on images that we generate to represent genomic fragments. We demonstrate that this image-based approach generalizes better than the existing tools.

Conclusions

We used a deep neural network and an image-based approach to detect the most out of the correct GI predictions made by other tools, in addition to making novel GI predictions. The fact that the deep neural network was re-trained on only a limited number of GI datasets and then successfully generalized indicates that this approach could be applied to other problems in the field where data is still lacking or hard to curate.

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

Transcriptional Profiling of Three Pseudomonas syringae pv. actinidiae Biovars Reveals Different Responses to Apoplast-Like Conditions Related to Strain Virulence on the Host

Pseudomonas syringae pv. actinidiae is a phytopathogen that causes devastating bacterial canker in kiwifruit. Among five biovars defined by genetic, biochemical, and virulence traits, P. syringae pv. actinidiae biovar 3 (Psa3) is the most aggressive and is responsible for the most recent reported outbreaks; however, the molecular basis of its heightened virulence is unclear. Therefore, we designed the first P. syringae multistrain whole-genome microarray, encompassing biovars Psa1, Psa2, and Psa3 and the well-established model P. syringae pv. tomato, and analyzed early bacterial responses to an apoplast-like minimal medium. Transcriptomic profiling revealed i) the strong activation in Psa3 of all hypersensitive reaction and pathogenicity (hrp) and hrp conserved (hrc) cluster genes, encoding components of the type III secretion system required for bacterial pathogenicity and involved in responses to environmental signals; ii) potential repression of the hrp/hrc cluster in Psa2; and iii) activation of flagellum-dependent cell motility and chemotaxis genes in Psa1. The detailed investigation of three gene families encoding upstream regulatory proteins (histidine kinases, their cognate response regulators, and proteins with diguanylate cyclase or phosphodiesterase domains) indicated that cyclic di-GMP may be a key regulator of virulence in P. syringae pv. actinidiae biovars. The gene expression data were supported by the quantification of biofilm formation. Our findings suggest that diverse early responses to the host apoplast, even among bacteria belonging to the same pathovar, can lead to different virulence strategies and may explain the differing outcomes of infections. Based on our detailed structural analysis of hrp operons, we also propose a revision of hrp cluster organization and operon regulation in P. syringae.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Detection of horizontal gene transfer in the genome of the choanoflagellate Salpingoeca rosetta

Horizontal gene transfer (HGT), the movement of heritable materials between distantly related organisms, is crucial in eukaryotic evolution. However, the scale of HGT in choanoflagellates, the closest unicellular relatives of metazoans, and its possible roles in the evolution of animal multicellularity remains unexplored. We identified at least 175 candidate HGTs in the genome of the colonial choanoflagellate Salpingoeca rosetta using sequence-based tests. The majority of these were orthologous to genes in bacterial and microalgal lineages, yet displayed genomic features consistent with the rest of the S. rosetta genome-evidence of ancient acquisition events. Putative functions include enzymes involved in amino acid and carbohydrate metabolism, cell signaling, and the synthesis of extracellular matrix components. Functions of candidate HGTs may have contributed to the ability of choanoflagellates to assimilate novel metabolites, thereby supporting adaptation, survival in diverse ecological niches, and response to external cues that are possibly critical in the evolution of multicellularity in choanoflagellates.

Oxford nanopore sequencing in clinical microbiology and infection diagnostics

Extended turnaround times and large economic costs hinder the usage of currently applied screening methods for bacterial pathogen identification (ID) and antimicrobial susceptibility testing. This review provides an overview of current detection methods and their usage in a clinical setting. Issues of timeliness and cost could soon be circumvented, however, with the emergence of detection methods involving single molecule sequencing technology. In the context of bringing diagnostics closer to the point of care, we examine the current state of Oxford Nanopore Technologies (ONT) products and their interaction with third-party software/databases to assess their capabilities for ID and antimicrobial resistance (AMR) prediction. We outline and discuss a potential diagnostic workflow, enumerating (1) rapid sample prep kits, (2) ONT hardware/software and (3) third-party software and databases to improve the cost, accuracy and turnaround times for ID and AMR. Multiple studies across a range of infection types support that the speed and accuracy of ONT sequencing is now such that established ID and AMR prediction tools can be used on its outputs, and so it can be harnessed for near real time, close to the point-of-care diagnostics in common clinical circumstances.

Comparative Genomic Analysis of Mycobacteriaceae Reveals Horizontal Gene Transfer-Mediated Evolution of the CRISPR-Cas System in the Mycobacterium tuberculosis Complex

Clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated (Cas) genes are conserved genetic elements in many prokaryotes, including Mycobacterium tuberculosis, the causative agent of tuberculosis. Although knowledge of CRISPR locus variability has been utilized in M. tuberculosis strain genotyping, its evolutionary path in Mycobacteriaceae is not well understood. In this study, we have performed a comparative analysis of 141 mycobacterial genomes and identified the exclusive presence of the CRISPR-Cas type III-A system in M. tuberculosis complex (MTBC). Our global phylogenetic analysis of CRISPR repeats and Cas10 proteins offers evidence of horizontal gene transfer (HGT) of the CRISPR-Cas module in the last common ancestor of MTBC and Mycobacterium canettii from a Streptococcus-like environmental bacterium. Additionally, our results show that the variation of CRISPR-Cas organization in M. tuberculosis lineages, especially in the Beijing sublineage of lineage 2, is due to the transposition of insertion sequence IS6110 The direct repeat (DR) region of the CRISPR-Cas locus acts as a hot spot for IS6110 insertion. We show in M. tuberculosis H37Rv that the repeat at the 5' end of CRISPR1 of the forward strand is an atypical repeat made up partly of IS-terminal inverted repeat and partly CRISPR DR. By tracing an undetectable spacer sequence in the DR region, the two CRISPR loci could theoretically be joined to reconstruct the ancestral single CRISPR-Cas locus organization, as seen in M. canettii This study retracing the evolutionary events of HGT and IS6110-driven genomic deletions helps us to better understand the strain-specific variations in M. tuberculosis lineages.IMPORTANCE Comparative genomic analysis of prokaryotes has led to a better understanding of the biology of several pathogenic microorganisms. One such clinically important pathogen is M. tuberculosis, the leading cause of bacterial infection worldwide. Recent evidence on the functionality of the CRISPR-Cas system in M. tuberculosis has brought back focus on these conserved genetic elements, present in many prokaryotes. Our study advances understanding of mycobacterial CRISPR-Cas origin and its diversity among the different species. We provide phylogenetic evidence of acquisition of CRISPR-Cas type III-A in the last common ancestor shared between MTBC and M. canettii, by HGT-mediated events. The most likely source of HGT was an environmental Firmicutes bacterium. Genomic mapping of the CRISPR loci showed the IS6110 transposition-driven variations in M. tuberculosis strains. Thus, this study offers insights into events related to the evolution of CRISPR-Cas in M. tuberculosis lineages.

Complete Genome Sequence of Ovine Mycobacterium avium subsp. paratuberculosis Strain JIII-386 (MAP-S/type III) and Its Comparison to MAP-S/type I, MAP-C, and M. avium Complex Genomes

Mycobacterium avium (M. a.) subsp. paratuberculosis (MAP) is a worldwide-distributed obligate pathogen in ruminants causing Johne’s disease. Due to a lack of complete subtype III genome sequences, there is not yet conclusive information about genetic differences between strains of cattle (MAP-C, type II) and sheep (MAP-S) type, and especially between MAP-S subtypes I, and III. Here we present the complete, circular genome of MAP-S/type III strain JIII-386 (DE) closed by Nanopore-technology and its comparison with MAP-S/type I closed genome of strain Telford (AUS), MAP-S/type III draft genome of strain S397 (U.S.), twelve closed MAP-C strains, and eight closed M.-a.-complex-strains. Structural comparative alignments revealed clearly the mosaic nature of MAP, emphasized differences between the subtypes and the higher diversity of MAP-S genomes. The comparison of various genomic elements including transposases and genomic islands provide new insights in MAP genomics. MAP type specific phenotypic features may be attributed to genes of known large sequence polymorphisms (LSP^Ss) regions I–IV and deletions #1 and #2, confirmed here, but could also result from identified frameshifts or interruptions of various virulence-associated genes (e.g., mbtC in MAP-S). Comprehensive core and pan genome analysis uncovered unique genes (e.g., cytochromes) and genes probably acquired by horizontal gene transfer in different MAP-types and subtypes, but also emphasized the highly conserved and close relationship, and the complex evolution of M.-a.-strains.

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

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

Pathogenicity Factors of Genomic Islands in Intestinal and Extraintestinal Escherichia coli

Escherichia coli is a versatile bacterial species that includes both harmless commensal strains and pathogenic strains found in the gastrointestinal tract in humans and warm-blooded animals. The growing amount of DNA sequence information generated in the era of "genomics" has helped to increase our understanding of the factors and mechanisms involved in the diversification of this bacterial species. The pathogenic side of E. coli that is afforded through horizontal transfers of genes encoding virulence factors enables this bacterium to become a highly diverse and adapted pathogen that is responsible for intestinal or extraintestinal diseases in humans and animals. Many of the accessory genes acquired by horizontal transfers form syntenic blocks and are recognized as genomic islands (GIs). These genomic regions contribute to the rapid evolution, diversification and adaptation of E. coli variants because they are frequently subject to rearrangements, excision and transfer, as well as to further acquisition of additional DNA. Here, we review a subgroup of GIs from E. coli termed pathogenicity islands (PAIs), a concept defined in the late 1980s by Jörg Hacker and colleagues in Werner Goebel's group at the University of Würzburg, Würzburg, Germany. As with other GIs, the PAIs comprise large genomic regions that differ from the rest of the genome by their G + C content, by their typical insertion within transfer RNA genes, and by their harboring of direct repeats (at their ends), integrase determinants, or other mobility loci. The hallmark of PAIs is their contribution to the emergence of virulent bacteria and to the development of intestinal and extraintestinal diseases. This review summarizes the current knowledge on the structure and functional features of PAIs, on PAI-encoded E. coli pathogenicity factors and on the role of PAIs in host-pathogen interactions.

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

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

Computational Basis for On-Demand Production of Diversified Therapeutic Phage Cocktails

The antibiotic resistance crisis has led to renewed interest in phage therapy as an alternative means of treating infection. However, conventional methods for isolating pathogen-specific phage are slow, labor-intensive, and frequently unsuccessful. We have demonstrated that computationally identified prophages carried by near-neighbor bacteria can serve as starting material for production of engineered phages that kill the target pathogen. Our approach and technology platform offer new opportunity for rapid development of phage therapies against most, if not all, bacterial pathogens, a foundational advance for use of phage in treating infectious disease.

New therapies are necessary to combat increasingly antibiotic-resistant bacterial pathogens. We have developed a technology platform of computational, molecular biology, and microbiology tools which together enable on-demand production of phages that target virtually any given bacterial isolate. Two complementary computational tools that identify and precisely map prophages and other integrative genetic elements in bacterial genomes are used to identify prophage-laden bacteria that are close relatives of the target strain. Phage genomes are engineered to disable lysogeny, through use of long amplicon PCR and Gibson assembly. Finally, the engineered phage genomes are introduced into host bacteria for phage production. As an initial demonstration, we used this approach to produce a phage cocktail against the opportunistic pathogen Pseudomonas aeruginosa PAO1. Two prophage-laden P. aeruginosa strains closely related to PAO1 were identified, ATCC 39324 and ATCC 27853. Deep sequencing revealed that mitomycin C treatment of these strains induced seven phages that grow on P. aeruginosa PAO1. The most diverse five phages were engineered for nonlysogeny by deleting the integrase gene (int), which is readily identifiable and typically conveniently located at one end of the prophage. The Δint phages, individually and in cocktails, killed P. aeruginosa PAO1 in liquid culture as well as in a waxworm (Galleria mellonella) model of infection.

IMPORTANCE The antibiotic resistance crisis has led to renewed interest in phage therapy as an alternative means of treating infection. However, conventional methods for isolating pathogen-specific phage are slow, labor-intensive, and frequently unsuccessful. We have demonstrated that computationally identified prophages carried by near-neighbor bacteria can serve as starting material for production of engineered phages that kill the target pathogen. Our approach and technology platform offer new opportunity for rapid development of phage therapies against most, if not all, bacterial pathogens, a foundational advance for use of phage in treating infectious disease.

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.

Complete genome sequence of the aromatic-hydrocarbon-degrading bacterium Achromobacter xylosoxidans DN002

Achromobacter xylosoxidans DN002 is capable of utilizing numerous aromatic hydrocarbons as sole carbon and energy resource. In this study, the whole genome of strain DN002 was sequenced and analyzed, which consisted of one circular chromosome of 5,943,204 bp and a 278,917 bp plasmid with an average GC content of 65.46 mol%, 5694 protein-coding genes, 13 rRNA genes and 57 tRNA genes. Analysis of cluster of orthologous group (COG) demonstrated that strain DN002 had remarkable gene abundance foramino acid transport and metabolism, transcription, inorganic ion transport and metabolism, energy production and conversion, and carbohydrate transport and metabolism. Genes related to biodegradation of aromatic hydrocarbons, chemotaxis and flagella were identified from the genome, which will advance our fundamental understanding the molecular mechanism for degradation and metabolizing of aromatic hydrocarbons.