Genomic landscapes of divergence among island bird populations: Evidence of parallel adaptation but at different loci?

When populations colonise new environments, they may be exposed to novel selection pressures but also suffer from extensive genetic drift due to founder effects, small population sizes and limited interpopulation gene flow. Genomic approaches enable us to study how these factors drive divergence, and disentangle neutral effects from differentiation at specific loci due to selection. Here, we investigate patterns of genetic diversity and divergence using whole-genome resequencing (>22× coverage) in Berthelot's pipit (Anthus berthelotii), a passerine endemic to the islands of three north Atlantic archipelagos. Strong environmental gradients, including in pathogen pressure, across populations in the species range, make it an excellent system in which to explore traits important in adaptation and/or incipient speciation. First, we quantify how genomic divergence accumulates across the speciation continuum, that is, among Berthelot's pipit populations, between sub species across archipelagos, and between Berthelot's pipit and its mainland ancestor, the tawny pipit (Anthus campestris). Across these colonisation timeframes (2.1 million-ca. 8000 years ago), we identify highly differentiated loci within genomic islands of divergence and conclude that the observed distributions align with expectations for non-neutral divergence. Characteristic signatures of selection are identified in loci associated with craniofacial/bone and eye development, metabolism and immune response between population comparisons. Interestingly, we find limited evidence for repeated divergence of the same loci across the colonisation range but do identify different loci putatively associated with the same biological traits in different populations, likely due to parallel adaptation. Incipient speciation across these island populations, in which founder effects and selective pressures are strong, may therefore be repeatedly associated with morphology, metabolism and immune defence.

Divergent Selection in Low Recombination Regions Shapes the Genomic Islands in Two Incipient Shorebird Species

Speciation in the face of gene flow is usually associated with a heterogeneous genomic landscape of divergence in nascent species pairs. However, multiple factors, such as divergent selection and local recombination rate variation, can influence the formation of these genomic islands. Examination of the genomic landscapes of species pairs that are still in the early stages of speciation provides an insight into this conundrum. In this study, population genomic analyses were undertaken using a wide range of sampling and whole-genome resequencing data from 96 unrelated individuals of Kentish plover (Charadrius alexandrinus) and white-faced plover (Charadrius dealbatus). We suggest that the two species exhibit varying levels of population admixture along the Chinese coast and on the Taiwan Island. Genome-wide analyses for introgression indicate that ancient introgression had occurred in Taiwan population, and gene flow is still ongoing in mainland coastal populations. Furthermore, we identified a few genomic regions with significant levels of interspecific differentiation and local recombination suppression, which contain several genes potentially associated with disease resistance, coloration, and regulation of plumage molting and thus may be relevant to the phenotypic and ecological divergence of the two nascent species. Overall, our findings suggest that divergent selection in low recombination regions may be a main force in shaping the genomic islands in two incipient shorebird species.

Genomic islands and their role in fitness traits of two key sepsis-causing bacterial pathogens

To survive and establish a niche for themselves, bacteria constantly evolve. Toward that, they not only insert point mutations and promote illegitimate recombinations within their genomes but also insert pieces of 'foreign' deoxyribonucleic acid, which are commonly referred to as 'genomic islands' (GEIs). The GEIs come in several forms, structures and types, often providing a fitness advantage to the harboring bacterium. In pathogenic bacteria, some GEIs may enhance virulence, thus altering disease burden, morbidity and mortality. Hence, delineating (i) the GEIs framework, (ii) their encoded functions, (iii) the triggers that help them move, (iv) the mechanisms they exploit to move among bacteria and (v) identification of their natural reservoirs will aid in superior tackling of several bacterial diseases, including sepsis. Given the vast array of comparative genomics data, in this short review, we provide an overview of the GEIs, their types and the compositions therein, especially highlighting GEIs harbored by two important pathogens, viz. Acinetobacter baumannii and Klebsiella pneumoniae, which prominently trigger sepsis in low- and middle-income countries. Our efforts help shed some light on the challenges these pathogens pose when equipped with GEIs. We hope that this review will provoke intense research into understanding GEIs, the cues that drive their mobility across bacteria and the ways and means to prevent their transfer, especially across pathogenic bacteria.

Flexible genomic island conservation across freshwater and marine Methylophilaceae

The evolutionary trajectory of Methylophilaceae includes habitat transitions from freshwater sediments to freshwater and marine pelagial that resulted in genome reduction (genome-streamlining) of the pelagic taxa. However, the extent of genetic similarities in the genomic structure and microdiversity of the two genome-streamlined pelagic lineages (freshwater "Ca. Methylopumilus" and the marine OM43 lineage) has so far never been compared. Here, we analyzed complete genomes of 91 "Ca. Methylopumilus" strains isolated from 14 lakes in Central Europe and 12 coastal marine OM43 strains. The two lineages showed a remarkable niche differentiation with clear species-specific differences in habitat preference and seasonal distribution. On the other hand, we observed a synteny preservation in their genomes by having similar locations and types of flexible genomic islands (fGIs). Three main fGIs were identified: a replacement fGI acting as phage defense, an additive fGI harboring metabolic and resistance-related functions, and a tycheposon containing nitrogen-, thiamine-, and heme-related functions. The fGIs differed in relative abundances in metagenomic datasets suggesting different levels of variability ranging from strain-specific to population-level adaptations. Moreover, variations in one gene seemed to be responsible for different growth at low substrate concentrations and a potential biogeographic separation within one species. Our study provides a first insight into genomic microdiversity of closely related taxa within the family Methylophilaceae and revealed remarkably similar dynamics involving mobile genetic elements and recombination between freshwater and marine family members.

Comprehensive Analysis of Virulence Determinants and Genomic Islands of blaNDM-1-Producing Enterobacter hormaechei Clinical Isolates from Greece

Nosocomial outbreaks of multidrug-resistant (MDR) Enterobacter cloacae complex (ECC) are often reported worldwide, mostly associated with a small number of multilocus-sequence types of E. hormaechei and E. cloacae strains. In Europe, the largest clonal outbreak of blaNDM-1-producing ECC has been recently reported, involving an ST182 E. hormaechei strain in a Greek teaching hospital. In the current study, we aimed to further investigate the genetic make-up of two representative outbreak isolates. Comparative genomics of whole genome sequences (WGS) was performed, including whole genome-based taxonomic analysis and in silico prediction of virulence determinants of the bacterial cell surface, plasmids, antibiotic resistance genes and virulence factors present on genomic islands. The enterobacterial common antigen and the colanic antigen of the cell surface were identified in both isolates, being similar to the gene clusters of the E. hormaechei ATCC 49162 and E. cloacae ATCC 13047 type strains, whereas the two strains possessed different gene clusters encoding lipopolysaccharide O-antigens. Other virulence factors of the bacterial cell surface, such as flagella, fimbriae and pili, were also predicted to be encoded by gene clusters similar to those found in Enterobacter spp. and other Enterobacterales. Secretion systems and toxin-antitoxin systems, which also contribute to pathogenicity, were identified. Both isolates harboured resistance genes to multiple antimicrobial classes, including β-lactams, aminoglycosides, quinolones, chloramphenicol, trimethoprim, sulfonamides and fosfomycin; they carried blaTEM-1, blaOXA-1, blaNDM-1, and one of them also carried blaCTXM-14, blaCTXM-15 and blaLAP-2 plasmidic alleles. Our comprehensive analysis of the WGS assemblies revealed that blaNDM-1-producing outbreak isolates possess components of the bacterial cell surface as well as genomic islands, harbouring resistance genes to several antimicrobial classes and various virulence factors. Differences in the plasmids carrying β-lactamase genes between the two strains have also shown diverse modes of acquisition and an ongoing evolution of these mobile elements.

Staphylococcus brunensis sp. nov. isolated from human clinical specimens with a staphylococcal cassette chromosome-related genomic island outside of the rlmH gene bearing the ccrDE recombinase gene complex

Novel species of coagulase-negative staphylococci, which could serve as reservoirs of virulence and antimicrobial resistance factors for opportunistic pathogens from the genus Staphylococcus, are recognized in human and animal specimens due to advances in diagnostic techniques. Here, we used whole-genome sequencing, extensive biotyping, MALDI-TOF mass spectrometry, and chemotaxonomy to characterize five coagulase-negative strains from the Staphylococcus haemolyticus phylogenetic clade obtained from human ear swabs, wounds, and bile. Based on the results of polyphasic taxonomy, we propose the species Staphylococcus brunensis sp. nov. (type strain NRL/St 16/872T = CCM 9024T = LMG 31872T = DSM 111349T). The genomic analysis revealed numerous variable genomic elements, including staphylococcal cassette chromosome (SCC), prophages, plasmids, and a unique 18.8 kb-long genomic island SbCIccrDE integrated into the ribosomal protein L7 serine acetyltransferase gene rimL. SbCIccrDE has a cassette chromosome recombinase (ccr) gene complex with a typical structure found in SCCs. Based on nucleotide and amino acid identity to other known ccr genes and the distinct integration site that differs from the canonical methyltransferase gene rlmH exploited by SCCs, we classified the ccr genes as novel variants, ccrDE. The comparative genomic analysis of SbCIccrDE with related islands shows that they can accumulate virulence and antimicrobial resistance factors creating novel resistance elements, which reflects the evolution of SCC. The spread of these resistance islands into established pathogens such as Staphylococcus aureus would pose a great threat to the healthcare system. IMPORTANCE The coagulase-negative staphylococci are important opportunistic human pathogens, which cause bloodstream and foreign body infections, mainly in immunocompromised patients. The mobile elements, primarily the staphylococcal cassette chromosome mec, which confers resistance to methicillin, are the key to the successful dissemination of staphylococci into healthcare and community settings. Here, we present a novel species of the Staphylococcus genus isolated from human clinical material. The detailed analysis of its genome revealed a previously undescribed genomic island, which is closely related to the staphylococcal cassette chromosome and has the potential to accumulate and spread virulence and resistance determinants. The island harbors a set of conserved genes required for its mobilization, which we recognized as novel cassette chromosome recombinase genes ccrDE. Similar islands were revealed not only in the genomes of coagulase-negative staphylococci but also in S. aureus. The comparative genomic study contributes substantially to the understanding of the evolution and pathogenesis of staphylococci.

Antimicrobial resistance, pathogenic potential, and genomic features of carbapenem-resistant Klebsiella pneumoniae isolated in Chile: high-risk ST25 clones and novel mobile elements

Multidrug- and carbapenem-resistant Klebsiella pneumoniae (CR-Kp) are critical threats to global health and key traffickers of resistance genes to other pathogens. Despite the sustained increase in CR-Kp infections in Chile, few strains have been described at the genomic level, lacking details of their resistance and virulence determinants and the mobile elements mediating their dissemination. In this work, we studied the antimicrobial susceptibility and performed a comparative genomic analysis of 10 CR-Kp isolates from the Chilean surveillance of carbapenem-resistant Enterobacteriaceae. High resistance was observed among the isolates (five ST25, three ST11, one ST45, and one ST505), which harbored 44 plasmids, most carrying genes for conjugation and resistance to several antibiotics and biocides. Ten plasmids encoding carbapenemases were characterized, including novel plasmids or variants with additional resistance genes, a novel genetic environment for blaKPC-2, and plasmids widely disseminated in South America. ST25 K2 isolates belonging to CG10224, a clone traced back to 2012 in Chile, which recently acquired blaNDM-1, blaNDM-7, or blaKPC-2 plasmids stood out as high-risk clones. Moreover, this corresponds to the first report of ST25 and ST45 Kp producing NDM-7 in South America and ST505 CR-Kp producing both NDM-7 and KPC-2 worldwide. Also, we characterized a variety of genomic islands carrying virulence and fitness factors. These results provide baseline knowledge for a detailed understanding of molecular and genetic determinants behind antibiotic resistance and virulence of CR-Kp in Chile and South America. IMPORTANCE In the ongoing antimicrobial resistance crisis, carbapenem-resistant strains of Klebsiella pneumoniae are critical threats to public health. Besides globally disseminated clones, the burden of local problem clones remains substantial. Although genomic analysis is a powerful tool for improving pathogen and antimicrobial resistance surveillance, it is still restricted in low- to middle-income countries, including Chile, causing them to be underrepresented in genomic databases and epidemiology surveys. This study provided the first 10 complete genomes of the Chilean surveillance for carbapenem-resistant K. pneumoniae in healthcare settings, unveiling their resistance and virulence determinants and the mobile genetic elements mediating their dissemination, placed in the South American and global K. pneumoniae epidemiological context. We found ST25 with K2 capsule as an emerging high-risk clone, along with other lineages producing two carbapenemases and several other resistance and virulence genes encoded in novel plasmids and genomic islands.

Comparative Genomic Analysis of Cold-Water Coral-Derived Sulfitobacter faviae: Insights into Their Habitat Adaptation and Metabolism

Sulfitobacter is one of the major sulfite-oxidizing alphaproteobacterial groups and is often associated with marine algae and corals. Their association with the eukaryotic host cell may have important ecological contexts due to their complex lifestyle and metabolism. However, the role of Sulfitobacter in cold-water corals remains largely unexplored. In this study, we explored the metabolism and mobile genetic elements (MGEs) in two closely related Sulfitobacter faviae strains isolated from cold-water black corals at a depth of ~1000 m by comparative genomic analysis. The two strains shared high sequence similarity in chromosomes, including two megaplasmids and two prophages, while both contained several distinct MGEs, including prophages and megaplasmids. Additionally, several toxin-antitoxin systems and other types of antiphage elements were also identified in both strains, potentially helping Sulfitobacter faviae overcome the threat of diverse lytic phages. Furthermore, the two strains shared similar secondary metabolite biosynthetic gene clusters and genes involved in dimethylsulfoniopropionate (DMSP) degradation pathways. Our results provide insight into the adaptive strategy of Sulfitobacter strains to thrive in ecological niches such as cold-water corals at the genomic level.

Complete Genome Sequence and Analysis of a ST573 Multidrug-Resistant Methicillin-Resistant Staphylococcus aureus SauR3 Clinical Isolate from Terengganu, Malaysia

Methicillin-resistant Staphylococcus aureus (MRSA) is a World Health Organization-listed priority pathogen. Scarce genomic data are available for MRSA isolates from Malaysia. Here, we present the complete genome sequence of a multidrug-resistant MRSA strain SauR3, isolated from the blood of a 6-year-old patient hospitalized in Terengganu, Malaysia, in 2016. S. aureus SauR3 was resistant to five antimicrobial classes comprising nine antibiotics. The genome was sequenced on the Illumina and Oxford Nanopore platforms and hybrid assembly was performed to obtain its complete genome sequence. The SauR3 genome consists of a circular chromosome of 2,800,017 bp and three plasmids designated pSauR3-1 (42,928 bp), pSauR3-2 (3011 bp), and pSauR3-3 (2473 bp). SauR3 belongs to sequence type 573 (ST573), a rarely reported sequence type of the staphylococcal clonal complex 1 (CC1) lineage, and harbors a variant of the staphylococcal cassette chromosome mec (SCCmec) type V (5C2&5) element which also contains the aac(6')-aph(2″) aminoglycoside-resistance genes. pSauR3-1 harbors several antibiotic resistance genes in a 14,095 bp genomic island (GI), previously reported in the chromosome of other staphylococci. pSauR3-2 is cryptic, whereas pSauR3-3 encodes the ermC gene that mediates inducible resistance to macrolide-lincosamide-streptogramin B (iMLS_B). The SauR3 genome can potentially be used as a reference genome for other ST573 isolates.

Burkholderia thailandensis Isolated from the Environment, United States

Burkholderia thailandensis, an opportunistic pathogen found in the environment, is a bacterium closely related to B. pseudomallei, the cause of melioidosis. Human B. thailandensis infections are uncommon. We isolated B. thailandensis from water in Texas and Puerto Rico and soil in Mississippi in the United States, demonstrating a potential public health risk.

Novel Plant-Associated Brevibacillus and Lysinibacillus Genomospecies Harbor a Rich Biosynthetic Potential of Antimicrobial Compounds

We have previously reported the draft genome sequences of 59 endospore-forming Gram-positive bacterial strains isolated from Vietnamese crop plants due to their ability to suppress plant pathogens. Based on their draft genome sequence, eleven of them were assigned to the Brevibacillus and one to the Lysinibacillus genus. Further analysis including full genome sequencing revealed that several of these strains represent novel genomospecies. In vitro and in vivo assays demonstrated their ability to promote plant growth, as well as the strong biocontrol potential of Brevibacilli directed against phytopathogenic bacteria, fungi, and nematodes. Genome mining identified 157 natural product biosynthesis gene clusters (BGCs), including 36 novel BGCs not present in the MIBiG data bank. Our findings indicate that plant-associated Brevibacilli are a rich source of putative antimicrobial compounds and might serve as a valuable starting point for the development of novel biocontrol agents.

Novel integrative elements and genomic plasticity in ocean ecosystems

Horizontal gene transfer accelerates microbial evolution. The marine picocyanobacterium Prochlorococcus exhibits high genomic plasticity, yet the underlying mechanisms are elusive. Here, we report a novel family of DNA transposons-"tycheposons"-some of which are viral satellites while others carry cargo, such as nutrient-acquisition genes, which shape the genetic variability in this globally abundant genus. Tycheposons share distinctive mobile-lifecycle-linked hallmark genes, including a deep-branching site-specific tyrosine recombinase. Their excision and integration at tRNA genes appear to drive the remodeling of genomic islands-key reservoirs for flexible genes in bacteria. In a selection experiment, tycheposons harboring a nitrate assimilation cassette were dynamically gained and lost, thereby promoting chromosomal rearrangements and host adaptation. Vesicles and phage particles harvested from seawater are enriched in tycheposons, providing a means for their dispersal in the wild. Similar elements are found in microbes co-occurring with Prochlorococcus, suggesting a common mechanism for microbial diversification in the vast oligotrophic oceans.

Sympatric or micro-allopatric speciation in a glacial lake? Genomic islands support neither

Apparent cases of sympatric speciation may actually be due to micro-allopatric or micro-parapatric speciation. One way to distinguish between these models is to examine the existence and nature of genomic islands of divergence, wherein divergent DNA segments are interspersed with low-divergence segments. Such islands should be rare or absent under micro-allopatric speciation but common in cases of speciation with gene flow. Sympatric divergence of endemic fishes is known from isolated saline, crater, postglacial, and ancient lakes. Two morphologically distinct cyprinid fishes, Gymnocypris eckloni scoliostomus (GS) and G. eckloni eckloni (GE), in a small glacial lake on the Qinghai-Tibet Plateau, Lake Sunmcuo, match the biogeographic criteria of sympatric speciation. In this study, we examined genome-wide variation in 46 individuals from these two groups. The divergence time between the GS and GE lineages was estimated to be 20-60 Kya. We identified 54 large genomic islands (≥100 kb) of speciation, which accounted for 89.4% of the total length of all genomic islands. These islands harboured divergent genes related to olfactory receptors and olfaction signals that may play important roles in food selection and assortative mating in fishes. Although the genomic islands clearly indicated speciation with gene flow and rejected micro-allopatric speciation, they were too large to support the hypothesis of sympatric speciation. Theoretical and recent empirical studies suggested that continual gene flow in sympatry should give rise to many small genomic islands (as small as a few kilobases in size). Thus, the observed pattern is consistent with the extensive evidence on parapatric speciation, in which adjacent habitats facilitate divergent selection but also permit gene flow during speciation. We suggest that many, if not most, of the reported cases of sympatric speciation are likely to be micro-parapatric speciation.

Comparative genomics of Bordetella pertussis and prediction of new vaccines and drug targets

Pertussis is a highly contagious respiratory disease caused by Bordetella pertussis, a Gram-negative bacterium described over a century ago. Despite broad vaccine coverage and treatment options, the disease is remerging as a public health problem especially in infants and older children. Recent data indicate re-emergence of the disease is related to bacterial resistance to immune defences and decreased vaccine effectiveness, which obviously suggests the need of new effective vaccines and drugs. In an attempt to contribute with solutions to this great challenge, bioinformatics tools were used to genetically comprehend the species of these bacteria and predict new vaccines and drug targets. In fact, approaches were used to analysis genomic plasticity, gene synteny and species similarities between the 20 genomes of Bordetella pertussis already available. Furthermore, it was conducted reverse vaccinology and docking analysis to identify proteins with potential to become vaccine and drug targets, respectively. The analyses showed the 20 genomes belongs to a homogeneous group that has preserved most of the genes over time. Besides that, were found genomics islands and good proteins to be candidates for vaccine and drugs. Taken together, these results suggests new possibilities that may be useful to develop new vaccines and drugs that will help the prevention and treatment strategies of pertussis disease caused by these Bordetella strains. Communicated by Ramaswamy H. Sarma.

Comparative Genomic Analysis of a Multidrug-Resistant Staphylococcus hominis ShoR14 Clinical Isolate from Terengganu, Malaysia, Led to the Discovery of Novel Mobile Genetic Elements

Staphylococcus hominis is a coagulase-negative Staphylococcus (CoNS) commensal capable of causing serious systemic infections in humans. The emergence of multidrug-resistant S. hominis strains is of concern but little is known about the characteristics of this organism, particularly from Malaysia. Here, we present the comparative genome analysis of S. hominis ShoR14, a multidrug-resistant, methicillin-resistant blood isolate from Terengganu, Malaysia. Genomic DNA of S. hominis ShoR14 was sequenced on the Illumina platform and assembled using Unicycler v0.4.8. ShoR14 belonged to sequence type (ST) 1 which is the most prevalent ST of the S. hominis subsp. hominis. Comparative genomic analysis with closely related strains in the database with complete genome sequences, led to the discovery of a novel variant of the staphylococcal chromosome cassette mec (SCCmec) type VIII element harboring the mecA methicillin-resistance gene in ShoR14 and its possible carriage of a SCCfus element that encodes the fusidic acid resistance gene (fusC). Up to seven possible ShoR14 plasmid contigs were identified, three of which harbored resistance genes for tetracycline (tetK), chloramphenicol (catA7), macrolides, lincosamides, and streptogramin B (ermC). Additionally, we report the discovery of a novel mercury-resistant transposon, Tn7456, other genomic islands, and prophages which make up the S. hominis mobilome.

Linked selection, differential introgression and recombination rate variation promote heterogeneous divergence in a pair of yellow croakers

Understanding the mechanisms underlying heterogeneous genomic divergence is of particular interest in evolutionary biology. Highly differentiated genomic regions, known as genomic islands, often evolve between diverging lineages. These genomic islands may be related to selection promoting adaptation or reproductive isolation. Based on whole genome assembly and genome-wide RAD sequencing in a pair of yellow croakers (genus: Larimichthys), we investigated the evolutionary processes shaping genomic landscapes of divergence. Demographic modelling indicated that the two species diverged following a secondary contact scenario, where differential introgression and linked selection were suggested to be involved in heterogeneous genomic divergence. We identified reduced recombination rate in genomic islands and a relatively good conservation of both genetic diversity and recombination landscapes between species, which highlight the roles of linked selection and recombination rate variation in promoting heterogeneous divergence in the common ancestral lineage of the two species. In addition, we found a positive correlation between differentiation (F_ST ) and absolute sequence divergence (D_xy ), and elevated D_xy in genomic islands, indicating that the genomic landscape of divergence was not shaped by linked selection alone. Restricted gene flow in highly differentiated regions has probably remodelled the landscape of heterogeneous genomic divergence. This study highlights that highly differentiated genomic regions can also arise from a combination of linked selection and differential gene flow in interaction with varying recombination rates.

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

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

Genomic Analysis of Carbapenem-Resistant Comamonas in Water Matrices: Implications for Public Health and Wastewater Treatments

Comamonas spp. are Gram-negative bacteria that catabolize a wide range of organic and inorganic substrates. Comamonas spp. are abundant in aquatic and soil environments, including wastewater, and can cause opportunistic infections in humans. Because of their potential in wastewater bioaugmentation and bioremediation strategies, the identification of Comamonas species harboring genes encoding carbapenemases and other clinically important antibiotic resistance genes warrant further investigation. Here, we present an analysis of 39 whole-genome sequences comprising three Comamonas species from aquatic environments in South Australia that were recovered on media supplemented with carbapenems. The analysis includes a detailed description of 33 Comamonas denitrificans isolates, some of which carried chromosomally acquired bla_GES-5, bla_OXA, and aminoglycoside resistance (aadA) genes located on putative genomic islands (GIs). All bla_GES-5- and bla_OXA-containing GIs appear to be unique to this Australian collection of C. denitrificans. Notably, most open reading frames (ORFs) within the GIs, including all antimicrobial resistance (AMR) genes, had adjacent attC sites, indicating that these ORFs are mobile gene cassettes. One C. denitrificans isolate carried an IncP-1 plasmid with genes involved in xenobiotic degradation and response to oxidative stress. Our assessment of the sequences highlights the very distant nature of C. denitrificans to the other Comamonas species and its apparent disposition to acquire antimicrobial resistance genes on putative genomic islands. IMPORTANCE Antimicrobial resistance (AMR) poses a global public health threat, and the increase in resistance to "last-resort drugs," such as carbapenems, is alarming. Wastewater has been flagged as a hot spot for AMR evolution. Comamonas spp. are among the most common bacteria in wastewater and play a role in its bioaugmentation. While the ability of Comamonas species to catabolize a wide range of organic and inorganic substrates is well documented, some species are also opportunistic pathogens. However, data regarding AMR in Comamonas spp. are limited. Here, through the genomic analyses of 39 carbapenem-resistant Comamonas isolates, we make several key observations, including the identification of a subset of C. denitrificans isolates that harbored genomic islands encoding carbapenemase bla_GES-5 or extended-spectrum β-lactamase bla_OXA alleles. Given the importance of Comamonas species in potential wastewater bioaugmentation and bioremediation strategies, as well as their status as emerging pathogens, the acquisition of critically important antibiotic resistance genes on genomic islands warrants future monitoring.

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

Insight Into Whole Genome of Aeromonas veronii Isolated From Freshwater Fish by Resistome Analysis Reveal Extensively Antibiotic Resistant Traits

Aeromonas veronii outbreaks in tilapia farming caused relatively high mortalities, and the bacteria was resistant to many kinds of antimicrobials used in Thailand aquaculture. According to the CLSI standard, the determination of antimicrobials efficacy has been limited to phenotypic analyses, and a genomics study is required. This research aimed to analyze the resistome of A. veronii isolated from diseased tilapia in Chainat, Nong Khai, and Uttaradit provinces in Thailand. A total of 12 isolates of A. veronii were identified based on the gyrB sequencing and then, the MIC values to eight antimicrobials (AMP, AML, GEN, ENR, OXO, OTC, SXT, and FFC) were determined. According to the MIC patterns, whole genome sequencing (WGS) of five representatives and resistome analysis were performed, including 15 genomes of A. veronii isolated from freshwater fish available in the NCBI. All tilapia isolates were susceptible to FFC but resistant to AML and AMP while OTC resistance was the most dominant. In addition to the WGS analysis, 4.5 Mbp of A. veronii was characterized. A total of 20 ARGs were detected by resistome analysis and 16 genes were shared among the A. veronii population. In conclusion, A. veronii strains isolated from tilapia exhibited a resistance to several antimicrobials and multidrug resistance (MDR) which was related to the presence of multiple ARGs. Aeromonas veronii shared the ARGs in their population worldwide with a possibility of a plasmid-mediated acquisition due to the presence of resistance islands.

Genomic Insights into Drug Resistance Determinants in Cedecea neteri, A Rare Opportunistic Pathogen

Cedecea, a genus in the Enterobacteriaceae family, includes several opportunistic pathogens reported to cause an array of sporadic acute infections, most notably of the lung and bloodstream. One species, Cedecea neteri, is associated with cases of bacteremia in immunocompromised hosts and has documented resistance to different antibiotics, including β-lactams and colistin. Despite the potential to inflict serious infections, knowledge about drug resistance determinants in Cedecea is limited. In this study, we utilized whole-genome sequence data available for three environmental strains (SSMD04, M006, ND14a) of C. neteri and various bioinformatics tools to analyze drug resistance genes in this bacterium. All three genomes harbor multiple chromosome-encoded β-lactamase genes. A deeper analysis of β-lactamase genes in SSMD04 revealed four metallo-β-lactamases, a novel variant, and a CMY/ACT-type AmpC putatively regulated by a divergently transcribed AmpR. Homologs of known resistance-nodulation-cell division (RND)-type multidrug efflux pumps such as OqxB, AcrB, AcrD, and MdtBC were also identified. Genomic island prediction for SSMD04 indicated that tolC, involved in drug and toxin export across the outer membrane of Gram-negative bacteria, was acquired by a transposase-mediated genetic transfer mechanism. Our study provides new insights into drug resistance mechanisms of an environmental microorganism capable of behaving as a clinically relevant opportunistic pathogen.

Experimental evolution supports signatures of sexual selection in genomic divergence

Comparative genomics has contributed to the growing evidence that sexual selection is an important component of evolutionary divergence and speciation. Divergence by sexual selection is implicated in faster rates of divergence of the X chromosome and of genes thought to underlie sexually selected traits, including genes that are sex biased in expression. However, accurately inferring the relative importance of complex and interacting forms of natural selection, demography, and neutral processes that occurred in the evolutionary past is challenging. Experimental evolution provides an opportunity to apply controlled treatments for multiple generations and examine the consequent genomic divergence. Here, we altered sexual selection intensity, elevating sexual selection in polyandrous lines and eliminating it in monogamous lines, and examined patterns of allele frequency divergence in the genome of Drosophila pseudoobscura after more than 160 generations of experimental evolution. Divergence is not uniform across the genome but concentrated in "islands," many of which contain candidate genes implicated in mating behaviors and other sexually selected phenotypes. These are more often seen on the X chromosome, which also shows greater divergence in F ST than neutral expectations. There are characteristic signatures of selection seen in these regions, with lower diversity on the X chromosome than the autosomes, and differences in diversity on the autosomes between selection regimes. Reduced Tajima's D within some of the divergent regions may imply that selective sweeps have occurred, despite considerable recombination. These changes are associated with both differential gene expression between the lines and sex-biased gene expression within the lines. Our results are very similar to those thought to implicate sexual selection in divergence between species and natural populations, and hence provide experimental support for the likely role of sexual selection in driving such types of genetic divergence, but also illustrate how variable outcomes can be for different genomic regions.

Genomic islands and the evolution of livestock-associated Staphylococcus aureus genomes

BACKGROUND

Genomic Islands (GIs) are commonly believed to be relics of horizontal transfer and associated with specific metabolic capacities, including virulence of the strain. Horizontal gene transfer (HGT) plays a vital role in the acquisition of GIs and the evolution and adaptation of bacterial genomes.

OBJECTIVE

The present study was designed to predict the GIs and role of HGT in evolution of livestock-associated Staphylococcus aureus (LA-SA).

METHODS

GIs were predicted with two methods namely, Ensemble algorithm for Genomic Island Detection (EGID) tool, and Seq word Sniffer script. Functional characterization of GI elements was performed with clustering of orthologs. The putative donor predictions of GIs was done with the aid of the pre_GI database.

RESULTS

The present study predicted a pan of 46 GIs across the LA-SA genomes. Functional characterization of GI sequences revealed few unique results like the presence of metabolic operons like leuABCD and folPK genes in GIs and showed the importance of GIs in the adaptation to the host niche. The developed framework for GI donor prediction results revealed Rickettsia and Mycoplasma as the major donors of GI elements.

CONCLUSIONS

The role of GIs during the evolutionary race of LA-SA could be concluded from the present study. Niche adaptation of LA-SA enhanced presumably due to these GIs. Future studies could focus on the evolutionary relationships between Rickettsia and Mycoplasma sp. with S. aureus and also the evolution of Leucine/Isoleucine mosaic operon (leuABCD).

Mobile Carbapenemase Genes in Pseudomonas aeruginosa

Carbapenem-resistant Pseudomonas aeruginosa is one of the major concerns in clinical settings impelling a great challenge to antimicrobial therapy for patients with infections caused by the pathogen. While membrane permeability, together with derepression of the intrinsic beta-lactamase gene, is the global prevailing mechanism of carbapenem resistance in P. aeruginosa, the acquired genes for carbapenemases need special attention because horizontal gene transfer through mobile genetic elements, such as integrons, transposons, plasmids, and integrative and conjugative elements, could accelerate the dissemination of the carbapenem-resistant P. aeruginosa. This review aimed to illustrate epidemiologically the carbapenem resistance in P. aeruginosa, including the resistance rates worldwide and the carbapenemase-encoding genes along with the mobile genetic elements responsible for the horizontal dissemination of the drug resistance determinants. Moreover, the modular mobile elements including the carbapenemase-encoding gene, also known as the P. aeruginosa resistance islands, are scrutinized mostly for their structures.

Whole Genome Sequence-Based Prediction of Resistance Determinants in High-Level Multidrug-Resistant Campylobacter jejuni Isolates in Lithuania

Spread of antibiotic resistance via mobile genetic elements associates with transfer of genes providing resistance against multiple antibiotics. Use of various comparative genomics analysis techniques enables to find intrinsic and acquired genes associated with phenotypic antimicrobial resistance (AMR) in Campylobacter jejuni genome sequences with exceptionally high-level multidrug resistance. In this study, we used whole genome sequences of seven C. jejuni to identify isolate-specific genomic features associated with resistance and virulence determinants and their role in multidrug resistance (MDR). All isolates were phenotypically highly resistant to tetracycline, ciprofloxacin, and ceftriaxone (MIC range from 64 to ≥256 µg/mL). Besides, two C. jejuni isolates were resistant to gentamicin, and one was resistant to erythromycin. The extensive drug-resistance profiles were confirmed for the two C. jejuni isolates assigned to ST-4447 (CC179). The most occurring genetic antimicrobial-resistance determinants were tetO, beta-lactamase, and multidrug efflux pumps. In this study, mobile genetic elements (MGEs) were detected in genomic islands carrying genes that confer resistance to MDR, underline their importance for disseminating antibiotic resistance in C. jejuni. The genomic approach showed a diverse distribution of virulence markers, including both plasmids and phage sequences that serve as horizontal gene transfer tools. The study findings describe in silico prediction of AMR and virulence genetics determinants combined with phenotypic AMR detection in multidrug-resistant C. jejuni isolates from Lithuania.

Comparative genome analyses of Mycobacteroides immunogenum reveals two potential novel subspecies

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Comparative Analysis of Genomic Island Prediction Tools

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

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

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

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

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

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

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

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

PURPOSE

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

METHODOLOGY

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

CONCLUSION

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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