The Home Environment Is a Reservoir for Methicillin-Resistant Coagulase-Negative Staphylococci and Mammaliicocci

Coagulase-negative staphylococci (CoNS) and mammaliicocci are opportunistic human and animal pathogens, often resistant to multiple antimicrobials, including methicillin. Methicillin-resistant CoNS (MRCoNS) have traditionally been linked to hospitals and healthcare facilities, where they are significant contributors to nosocomial infections. However, screenings of non-hospital environments have linked MRCoNS and methicillin-resistant mammaliicocci (MRM) to other ecological niches. The aim of this study was to explore the home environment as a reservoir for MRCoNS and MRM. A total of 33 households, including households with a dog with a methicillin-resistant staphylococcal infection, households with healthy dogs or cats and households without pets, were screened for MRCoNS and MRM by sampling one human, one pet (if present) and the environment. Samples were analyzed by a selective culture-based method, and bacterial species were identified by MALDI-TOF MS and tested for antibiotic susceptibility by the agar disk diffusion method. Following whole-genome sequencing, a large diversity of SCCmec elements and sequence types was revealed, which did not indicate any clonal dissemination of specific strains. Virulome and mobilome analyses indicated a high degree of species specificity. Altogether, this study documents that the home environment is a reservoir for a variety of MRCoNS and MRM regardless of the type of household.

Mechanisms Underlying the Overlooked Chiral Fungicide-Driven Enantioselective Proliferation of Antibiotic Resistance in Earthworm Intestinal Microbiome

From a "One Health" perspective, the global threat of antibiotic resistance genes (ARGs) is associated with modern agriculture practices including agrochemicals application. Chiral fungicides account for a considerable proportion of wildly used agrochemicals; however, whether and how their enantiomers lead to differential proliferation of antibiotic resistance in agricultural environments remain overlooked. Focused on the soil-earthworm ecosystem, we for the first time deciphered the mechanisms underlying the enantioselective proliferation of antibiotic resistance driven by the enantiomers of a typical chiral fungicide mandipropamid (i.e., R-MDP and S-MDP) utilizing a multiomic approach. Time-series metagenomic analysis revealed that R-MDP led to a significant enhancement of ARGs with potential mobility (particularly the plasmid-borne ARGs) in the earthworm intestinal microbiome. We further demonstrated that R-MDP induced a concentration-dependent facilitation of plasmid-mediated ARG transfer among microbes. In addition, transcriptomic analysis with verification identified the key aspects involved, where R-MDP enhanced cell membrane permeability, transfer ability, biofilm formation and quorum sensing, rebalanced energy production, and decreased cell mobility versus S-MDP. Overall, the findings provide novel insights into the enantioselective disruption of microbiome and resistome in earthworm gut by chiral fungicides and offer significant contributions to the comprehensive risk assessment of chiral agrochemicals in agroecosystems.

Detecting horizontal gene transfer among microbiota: an innovative pipeline for identifying co-shared genes within the mobilome through advanced comparative analysis

Horizontal gene transfer (HGT) is a key driver in the evolution of bacterial genomes. The acquisition of genes mediated by HGT may enable bacteria to adapt to ever-changing environmental conditions. Long-term application of antibiotics in intensive agriculture is associated with the dissemination of antibiotic resistance genes among bacteria with the consequences causing public health concern. Commensal farm-animal-associated gut microbiota are considered the reservoir of the resistance genes. Therefore, in this study, we identified known and not-yet characterized mobilized genes originating from chicken and porcine fecal samples using our innovative pipeline followed by network analysis to provide appropriate visualization to support proper interpretation.

Structure and dispersion of the conjugative mobilome in surface ocean bacterioplankton

Mobile genetic elements (MGEs), collectively referred to as the "mobilome", can have a significant impact on the fitness of microbial communities and therefore on ecological processes. Marine MGEs have mainly been associated with wide geographical and phylogenetic dispersal of adaptative traits. However, whether the structure of this mobilome exhibits deterministic patterns in the natural community is still an open question. The aim of this study was to characterize the structure of the conjugative mobilome in the ocean surface bacterioplankton by searching the publicly available marine metagenomes from the TARA Oceans survey, together with molecular markers, such as relaxases and type IV coupling proteins of the type IV secretion system (T4SS). The T4SS machinery was retrieved in more abundance than relaxases in the surface marine bacterioplankton. Moreover, among the identified MGEs, mobilizable elements were the most abundant, outnumbering self-conjugative sequences. Detection of a high number of incomplete T4SSs provides insight into possible strategies related to trans-acting activity between MGEs, and accessory functions of the T4SS (e.g. protein secretion), allowing the host to maintain a lower metabolic burden in the highly dynamic marine system. Additionally, the results demonstrate a wide geographical dispersion of MGEs throughout oceanic regions, while the Southern Ocean appears segregated from other regions. The marine mobilome also showed a high similarity of functions present in known plasmid databases. Moreover, cargo genes were mostly related to DNA processing, but scarcely associated with antibiotic resistance. Finally, within the MGEs, integrative and conjugative elements showed wider marine geographic dispersion than plasmids.

Effect of Antimicrobial Use in Conventional Versus Natural Cattle Feedlots on the Microbiome and Resistome

Antimicrobial use (AMU) in the livestock industry has been associated with increased levels of antimicrobial resistance. Recently, there has been an increase in the number of "natural" feedlots in the beef cattle sector that raise cattle without antibiotics. Shotgun metagenomics was employed to characterize the impact of AMU in feedlot cattle on the microbiome, resistome, and mobilome. Sequenced fecal samples identified a decline (q < 0.01) in the genera Methanobrevibacter and Treponema in the microbiome of naturally vs. conventionally raised feedlot cattle, but this difference was not (q > 0.05) observed in catch basin samples. No differences (q > 0.05) were found in the class-level resistome between feedlot practices. In fecal samples, decreases from conventional to natural (q < 0.05) were noted in reads for the antimicrobial-resistant genes (ARGs) mefA, tet40, tetO, tetQ, and tetW. Plasmid-associated ARGs were more common in feces from conventional than natural feedlot cattle. Interestingly, more chromosomal- than plasmid-associated macrolide resistance genes were observed in both natural and conventional feedlots, suggesting that they were more stably conserved than the predominately plasmid-associated tetracycline resistance genes. This study suggests that generationally selected resistomes through decades of AMU persist even after AMU ceases in natural production systems.

Genomic analysis of almost 8,000 Salmonella genomes reveals drivers and landscape of antimicrobial resistance in China

IMPORTANCE

We established the largest Salmonella genome database from China and presented the landscape and spatiotemporal dynamics of antimicrobial resistance genes. We also found that economic, climatic, and social factors can drive the rise of antimicrobial resistance. The Chinese local Salmonella genome database version 2 was released as an open-access database (https://nmdc.cn/clsgdbv2) and thus can assist surveillance studies across the globe. This database will help inform interventions for AMR, food safety, and public health.

The mobilome of Lactobacillus crispatus M247 includes two novel genetic elements: Tn7088 coding for a putative bacteriocin and the siphovirus prophage ΦM247

Lactobacillus crispatus is a member of the vaginal and gastrointestinal human microbiota. Here we determined the complete genome sequence of the probiotic strain M247 combining Nanopore and Illumina technologies. The M247 genome is organized in one circular chromosome of 2 336 109 bp, with a GC content of 37.04 % and 2303 ORFs, of which 1962 could be annotated. Analysis of the M247 mobilome, which accounts for 14 % of the whole genome, revealed the presence of: (i) Tn7088, a novel 14 105 bp long integrative and mobilizable element (IME) containing 16 ORFs; (ii) ΦM247, a novel 42 510 bp long siphovirus prophage containing 52 ORFs; (iii) three clustered regularly interspaced short palindromic repeats (CRISPRs); and (iv) 226 insertion sequences (ISs) belonging to 14 different families. Tn7088 has a modular organization including a mobilization module encoding FtsK homologous proteins and a relaxase, an integration/excision module coding for an integrase and an excisionase, and an adaptation module coding for a class I bacteriocin and homologous to the listeriolysin S (lls) locus of Listeria monocytogenes. Genome-wide homology search analysis showed the presence of Tn7088-like elements in 12 out of 23 L. crispatus complete public genomes. Mobilization and integration/excision modules are essentially conserved, while the adaptation module is variable since it is the target site for the integration of different ISs. Prophage ΦM247 contains genes for phage structural proteins, DNA replication and packaging, lysogenic and lytic cycles. ΦM247-like prophages are present in seven L. crispatus complete genomes, with sequence variability mainly due to the integration of ISs. PCR and sequencing showed that the Tn7088 IME excises from the M247 chromosome producing a circular form at a concentration of 4.32×10-5 copies per chromosome, and reconstitution of the Tn7088 chromosomal target site occurred at 6.65×10-4 copies per chromosome. The ΦM247 prophage produces an excised form and a reconstituted target site at a level of 3.90×10-5 and 2.48×10-5 copies per chromosome, respectively. This study identified two novel genetic elements in L. crispatus. Tn7088 represents the first example of an IME carrying a biosynthetic gene cluster for a class I bacteriocin in L. crispatus.

One Earth: The Equilibrium between the Human and the Bacterial Worlds

Misuse and abuse of antibiotics on humans, cattle, and crops have led to the selection of multi-resistant pathogenic bacteria, the most feared 'superbugs'. Infections caused by superbugs are progressively difficult to treat, with a subsequent increase in lethality: the toll on human lives is predicted to reach 10 million by 2050. Here we review three concepts linked to the growing resistance to antibiotics, namely (i) the Resistome, which refers to the collection of bacterial genes that confer resistance to antibiotics, (ii) the Mobilome, which includes all the mobile genetic elements that participate in the spreading of antibiotic resistance among bacteria by horizontal gene transfer processes, and (iii) the Nichome, which refers to the set of genes that are expressed when bacteria try to colonize new niches. We also discuss the strategies that can be used to tackle bacterial infections and propose an entente cordiale with the bacterial world so that instead of war and destruction of the 'fierce enemy' we can achieve a peaceful coexistence (the One Earth concept) between the human and the bacterial worlds. This, in turn, will contribute to microbial biodiversity, which is crucial in a globally changing climate due to anthropogenic activities.

The Mobilome-Enriched Genome of the Competence-Deficient Streptococcus pneumoniae BM6001, the Original Host of Integrative Conjugative Element Tn5253, Is Phylogenetically Distinct from Historical Pneumococcal Genomes

Streptococcus pneumoniae is an important human pathogen causing both mild and severe diseases. In this work, we determined the complete genome sequence of the S. pneumoniae clinical isolate BM6001, which is the original host of the ICE Tn5253. The BM6001 genome is organized in one circular chromosome of 2,293,748 base pairs (bp) in length, with an average GC content of 39.54%; the genome harbors a type 19F capsule locus, two tandem copies of pspC, the comC1-comD1 alleles and the type I restriction modification system SpnIII. The BM6001 mobilome accounts for 15.54% (356,521 bp) of the whole genome and includes (i) the ICE Tn5253 composite; (ii) the novel IME Tn7089; (iii) the novel transposon Tn7090; (iv) 3 prophages and 2 satellite prophages; (v) 5 genomic islands (GIs); (vi) 72 insertion sequences (ISs); (vii) 69 RUPs; (viii) 153 BOX elements; and (ix) 31 SPRITEs. All MGEs, except for the GIs, produce excised circular forms and attB site restoration. Tn7089 is 9089 bp long and contains 11 ORFs, of which 6 were annotated and code for three functions: integration/excision, mobilization and adaptation. Tn7090 is 9053 bp in size, flanked by two copies of ISSpn7, and contains seven ORFs organized as a single transcriptional unit, with genes encoding for proteins likely involved in the uptake and binding of Mg²⁺ cations in the adhesion to host cells and intracellular survival. BM6001 GIs, except for GI-BM6001.4, are variants of the pneumococcal TIGR4 RD5 region of diversity, pathogenicity island PPI1, R6 Cluster 4 and PTS island. Overall, prophages and satellite prophages contain genes predicted to encode proteins involved in DNA replication and lysogeny, in addition to genes encoding phage structural proteins and lytic enzymes carried only by prophages. ΦBM6001.3 has a mosaic structure that shares sequences with prophages IPP69 and MM1 and disrupts the competent comGC/cglC gene after chromosomal integration. Treatment with mitomycin C results in a 10-fold increase in the frequency of ΦBM6001.3 excised forms and comGC/cglC coding sequence restoration but does not restore competence for genetic transformation. In addition, phylogenetic analysis showed that BM6001 clusters in a small lineage with five other historical strains, but it is distantly related to the lineage due to its unique mobilome, suggesting that BM6001 has progressively accumulated many MGEs while losing competence for genetic transformation.

Metagenomics analysis of the neonatal intestinal resistome

Introduction

The intestinal microbiome forms a major reservoir for antibiotic resistance genes (ARGs). Little is known about the neonatal intestinal resistome.

Objective

The objective of this study was to investigate the intestinal resistome and factors that influence the abundance of ARGs in a large cohort of neonates.

Methods

Shotgun metagenomics was used to analyse the resistome in stool samples collected at 1 week of age from 390 healthy, term-born neonates who did not receive antibiotics.

Results

Overall, 913 ARGs belonging to 27 classes were identified. The most abundant ARGs were those conferring resistance to tetracyclines, quaternary ammonium compounds, and macrolide-lincosamide-streptogramin-B. Phylogenetic composition was strongly associated with the resistome composition. Other factors that were associated with the abundance of ARGs were delivery mode, gestational age, birth weight, feeding method, and antibiotics in the last trimester of pregnancy. Sex, ethnicity, probiotic use during pregnancy, and intrapartum antibiotics had little effect on the abundance of ARGs.

Conclusion

Even in the absence of direct antibiotic exposure, the neonatal intestine harbours a high abundance and a variety of ARGs.

Genome comparison reveals that Halobacterium salinarum 63-R2 is the origin of the twin laboratory strains NRC-1 and R1

The genome of Halobacterium strain 63-R2 was recently reported and provides the opportunity to resolve long-standing issues regarding the source of two widely used model strains of Halobacterium salinarum, NRC-1 and R1. Strain 63-R2 was isolated in 1934 from a salted buffalo hide (epithet "cutirubra"), along with another strain from a salted cow hide (91-R6^T , epithet "salinaria," the type strain of Hbt. salinarum). Both strains belong to the same species according to genome-based taxonomy analysis (TYGS), with chromosome sequences showing 99.64% identity over 1.85 Mb. The chromosome of strain 63-R2 is 99.99% identical to the two laboratory strains NRC-1 and R1, with only five indels, excluding the mobilome. The two reported plasmids of strain 63-R2 share their architecture with plasmids of strain R1 (pHcu43/pHS4, 99.89% identity; pHcu235/pHS3, 100.0% identity). We detected and assembled additional plasmids using PacBio reads deposited at the SRA database, further corroborating that strain differences are minimal. One plasmid, pHcu190 (190,816 bp) corresponds to pHS1 (strain R1) but is even more similar in architecture to pNRC100 (strain NRC-1). Another plasmid, pHcu229, assembled partially and completed in silico (229,124 bp), shares most of its architecture with pHS2 (strain R1). In deviating regions, it corresponds to pNRC200 (strain NRC-1). Further architectural differences between the laboratory strain plasmids are not unique, but are present in strain 63-R2, which contains characteristics from both of them. Based on these observations, it is proposed that the early twentieth-century isolate 63-R2 is the immediate ancestor of the twin laboratory strains NRC-1 and R1.

Staphylococcal diversity in atopic dermatitis from an individual to a global scale

Atopic dermatitis (AD) is a multifactorial, chronic relapsing disease associated with genetic and environmental factors. Among skin microbes, Staphylococcus aureus and Staphylococcus epidermidis are associated with AD, but how genetic variability and staphylococcal strains shape the disease remains unclear. We investigated the skin microbiome of an AD cohort (n = 54) as part of a prospective natural history study using shotgun metagenomic and whole genome sequencing, which we analyzed alongside publicly available data (n = 473). AD status and global geographical regions exhibited associations with strains and genomic loci of S. aureus and S. epidermidis. In addition, antibiotic prescribing patterns and within-household transmission between siblings shaped colonizing strains. Comparative genomics determined that S. aureus AD strains were enriched in virulence factors, whereas S. epidermidis AD strains varied in genes involved in interspecies interactions and metabolism. In both species, staphylococcal interspecies genetic transfer shaped gene content. These findings reflect the staphylococcal genomic diversity and dynamics associated with AD.

A survey on antimicrobial resistance genes of frequently used probiotic bacteria, 1901 to 2022

BackgroundAntimicrobial resistance (AMR) is caused by AMR determinants, mainly genes (ARGs) in the bacterial genome. Bacteriophages, integrative mobile genetic elements (iMGEs) or plasmids can allow ARGs to be exchanged among bacteria by horizontal gene transfer (HGT). Bacteria, including bacteria with ARGs, can be found in food. Thus, it is conceivable that in the gastrointestinal tract, bacteria from the gut flora could take up ARGs from food.AimThe study objective was to gain insight into the ARG set carried by commonly used probiotic bacteria that may enter the human body with non-fermented foods, fermented foods, or probiotic dietary supplements (FFPs) and to assess ARG mobility.MethodsNext generation sequencing whole genome data from 579 isolates of 12 commonly employed probiotic bacterial species were collected from a public repository. Using bioinformatical tools, ARGs were analysed and linkage with mobile genetic elements assessed.ResultsResistance genes were found in eight bacterial species. The ratios of ARG positive/negative samples per species were: Bifidobacterium animalis (65/0), Lactiplantibacillus plantarum (18/194), Lactobacillus delbrueckii (1/40), Lactobacillus helveticus (2/64), Lactococcus lactis (74/5), Leucoconstoc mesenteroides (4/8), Levilactobacillus brevis (1/46), Streptococcus thermophilus (4/19). In 66% (112/169) of the ARG-positive samples, at least one ARG could be linked to plasmids or iMGEs. No bacteriophage-linked ARGs were found.ConclusionThe finding of potentially mobile ARGs in probiotic strains for human consumption raises awareness of a possibility of ARG HGT in the gastrointestinal tract. In addition to existing recommendations, screening FFP bacterial strains for ARG content and mobility characteristics might be considered.

Avian strains of emerging pathogen Escherichia fergusonii are phylogenetically diverse and harbor the greatest AMR dissemination potential among different sources: Comparative genomic evidence

Introduction

Escherichia fergusonii is regarded as an emerging pathogen with zoonotic potential. In the current study, we undertook source-wise comparative genomic analyses (resistome, virulome, mobilome and pangenome) to understand the antimicrobial resistance, virulence, mobile genetic elements and phylogenetic diversity of E. fergusonii.

Methods

Six E. fergusonii strains (5 multidrug resistant strains and 1 biofilm former) were isolated from poultry (duck faeces and retail chicken samples). Following confirmation by phenotypic and molecular methods, the isolates were further characterized and their genomes were sequenced. Comparative resisto-virulo-mobilome analyses and pangenomics were performed for E. fergusonii genomes, while including 125 other E. fergusonii genomes available from NCBI database.

Results and discussion

Avian and porcine strains of E. fergusonii were found to carry significantly higher number of antimicrobial resistance genes (p < 0.05) and mobile genetic elements (plasmids, transposons and integrons) (p < 0.05), while the pathogenic potential of bovine strains was significantly higher compared to other strains (p < 0.05). Pan-genome development trends indicated open pan-genome for all strains (0 < γ < 1). Genomic diversity of avian strains was found to be greater than that from other sources. Phylogenetic analysis revealed close clustering among isolates of similar isolation source and geographical location. Indian isolates of E. fergusonii clustered closely with those from Chinese and a singleton Australian isolate. Overall, being the first pangenomic study on E. fergusonii, our analysis provided important cues on genomic features of the emerging pathogen E. fergusonii while highlighting the potential role of avian strains in dissemination of AMR.

New insight into the microbiome, resistome, and mobilome on the dental waste water in the context of heavy metal environment

Object

Hospital sewage have been associated with incorporation of antibiotic resistance genes (ARGs) and mobile genetic elements (MGEs) into microbes, which is considered as a key indicator for the spread of antimicrobial resistance (AMR). The compositions of dental waste water (DWW) contain heavy metals, the evolution of AMR and its effects on the water environment in the context of heavy metal environment have not been seriously investigated. Thus, our major aims were to elucidate the evolution of AMR in DWW.

Methods

DWW samples were collected from a major dental department. The presence of microbial communities, ARGs, and MGEs in untreated and treated (by filter membrane and ozone) samples were analyzed using metagenomics and bioinformatic methods.

Results

DWW-associated resistomes included 1,208 types of ARGs, belonging to 29 antibiotic types/subtypes. The most abundant types/subtypes were ARGs of multidrug resistance and of antibiotics that were frequently used in the clinical practice. Pseudomonas putida, Pseudomonas aeruginosa, Chryseobacterium indologenes, Sphingomonas laterariae were the main bacteria which hosted these ARGs. Mobilomes in DWW consisted of 93 MGE subtypes which belonged to 8 MGE types. Transposases were the most frequently detected MGEs which formed networks of communications. For example, ISCrsp1 and tnpA.5/4/11 were the main transposases located in the central hubs of a network. These significant associations between ARGs and MGEs revealed the strong potential of ARGs transmission towards development of antimicrobial-resistant (AMR) bacteria. On the other hand, treatment of DWW using membranes and ozone was only effective in removing minor species of bacteria and types of ARGs and MGEs.

Conclusion

DWW contained abundant ARGs, and MGEs, which contributed to the occurrence and spread of AMR bacteria. Consequently, DWW would seriously increase environmental health concerns which may be different but have been well-documented from hospital waste waters.

First comparative genomic characterization of the MSSA ST398 lineage detected in aquaculture and other reservoirs

Staphylococcus aureus ST398 can cause diseases in several different animals. In this study we analyzed ten S. aureus ST398 previously collected in three different reservoirs in Portugal (humans, gilthead seabream from aquaculture and dolphin from a zoo). Strains tested against sixteen antibiotics, by disk diffusion or minimum inhibitory concentration, showed decreased susceptibility to benzylpenicillin (all strains from gilthead seabream and dolphin) and to erythromycin with an iMLS_B phenotype (nine strains), and susceptibility to cefoxitin (methicillin-susceptible S. aureus, MSSA). All strains from aquaculture belonged to the same spa type, t2383, whereas strains from the dolphin and humans belonged to spa type t571. A more detailed analysis using single nucleotide polymorphisms (SNPs)-based tree and a heat map, showed that all strains from aquaculture origin were highly related with each other and the strains from dolphin and humans were more distinct, although they were very similar in ARG, VF and MGE content. Mutations F3I and A100V in glpT gene and D278E and E291D in murA gene were identified in nine fosfomycin susceptible strains. The blaZ gene was also detected in six of the seven animal strains. The study of the genetic environment of erm(T)-type (found in nine S. aureus strains) allowed the identification of MGE (rep13-type plasmids and IS431R-type), presumably involved in the mobilization of this gene. All strains showed genes encoding efflux pumps from major facilitator superfamily (e.g., arlR, lmrS-type and norA/B-type), ATP-binding cassettes (ABC; mgrA) and multidrug and toxic compound extrusion (MATE; mepA/R-type) families, all associated to decreased susceptibility to antibiotics/disinfectants. Moreover, genes related with tolerance to heavy metals (cadD), and several VF (e.g., scn, aur, hlgA/B/C and hlb) were also identified. Insertion sequences, prophages, and plasmids made up the mobilome, some of them associated with ARG, VF and genes related with tolerance to heavy metals. This study highlights that S. aureus ST398 can be a reservoir of several ARG, heavy metals resistance genes and VF, which are essential in the adaption and survival of the bacterium in the different environments and an active agent in its dissemination. It makes an important contribution to understanding the extent of the spread of antimicrobial resistance, as well as the virulome, mobilome and resistome of this dangerous lineage.

Comparative genome analysis of the genus Shewanella unravels the association of key genetic traits with known and potential pathogenic lineages

Shewanella spp. are Gram-negative rods widely disseminated in aquatic niches that can also be found in human-associated environments. In recent years, reports of infections caused by these bacteria have increased significantly. Mobilome and resistome analysis of a few species showed that they are versatile; however, comprehensive comparative studies in the genus are lacking. Here, we analyzed the genetic traits of 144 genomes from Shewanella spp. isolates focusing on the mobilome, resistome, and virulome to establish their evolutionary relationship and detect unique features based on their genome content and habitat. Shewanella spp. showed a great diversity of mobile genetic elements (MGEs), most of them associated with monophyletic lineages of clinical isolates. Furthermore, 79/144 genomes encoded at least one antimicrobial resistant gene with their highest occurrence in clinical-related lineages. CRISPR-Cas systems, which confer immunity against MGEs, were found in 41 genomes being I-E and I-F the more frequent ones. Virulome analysis showed that all Shewanella spp. encoded different virulence genes (motility, quorum sensing, biofilm, adherence, etc.) that may confer adaptive advantages for survival against hosts. Our data revealed that key accessory genes are frequently found in two major clinical-related groups, which encompass the opportunistic pathogens Shewanella algae and Shewanella xiamenensis together with several other species. This work highlights the evolutionary nature of Shewanella spp. genomes, capable of acquiring different key genetic traits that contribute to their adaptation to different niches and facilitate the emergence of more resistant and virulent isolates that impact directly on human and animal health.

mobileOG-db: a Manually Curated Database of Protein Families Mediating the Life Cycle of Bacterial Mobile Genetic Elements

Bacterial mobile genetic elements (MGEs) encode functional modules that perform both core and accessory functions for the element, the latter of which are often only transiently associated with the element. The presence of these accessory genes, which are often close homologs to primarily immobile genes, incur high rates of false positives and, therefore, limits the usability of these databases for MGE annotation. To overcome this limitation, we analyzed 10,776,849 protein sequences derived from eight MGE databases to compile a comprehensive set of 6,140 manually curated protein families that are linked to the “life cycle” (integration/excision, replication/recombination/repair, transfer, stability/transfer/defense, and phage-specific processes) of plasmids, phages, integrative, transposable, and conjugative elements. We overlay experimental information where available to create a tiered annotation scheme of high-quality annotations and annotations inferred exclusively through bioinformatic evidence. We additionally provide an MGE-class label for each entry (e.g., plasmid or integrative element), and assign to each entry a major and minor category. The resulting database, mobileOG-db (for mobile orthologous groups), comprises over 700,000 deduplicated sequences encompassing five major mobileOG categories and more than 50 minor categories, providing a structured language and interpretable basis for an array of MGE-centered analyses. mobileOG-db can be accessed at mobileogdb.flsi.cloud.vt.edu/, where users can select, refine, and analyze custom subsets of the dynamic mobilome.

IMPORTANCE The analysis of bacterial mobile genetic elements (MGEs) in genomic data is a critical step toward profiling the root causes of antibiotic resistance, phenotypic or metabolic diversity, and the evolution of bacterial genera. Existing methods for MGE annotation pose high barriers of biological and computational expertise to properly harness. To bridge this gap, we systematically analyzed 10,776,849 proteins derived from eight databases of MGEs to identify 6,140 MGE protein families that can serve as candidate hallmarks, i.e., proteins that can be used as “signatures” of MGEs to aid annotation. The resulting resource, mobileOG-db, provides a multilevel classification scheme that encompasses plasmid, phage, integrative, and transposable element protein families categorized into five major mobileOG categories and more than 50 minor categories. mobileOG-db thus provides a rich resource for simple and intuitive element annotation that can be integrated seamlessly into existing MGE detection pipelines and colocalization analyses.

Development of Resistance to Clarithromycin and Amoxicillin-Clavulanic Acid in Lactiplantibacillus plantarum In Vitro Is Followed by Genomic Rearrangements and Evolution of Virulence

Ensuring the safety of the use of probiotics is a top priority. Obviously, in addition to studying the beneficial properties of lactic acid bacteria, considerable attention should be directed to assessing the virulence of microorganisms as well as investigating the possibility of its evolution under conditions of selective pressure. To assess the virulence of probiotics, it is now recommended to analyze the genomes of bacteria in relation to the profiles of the virulome, resistome, and mobilome as well as the analysis of phenotypic resistance and virulence in vitro. However, the corresponding procedure has not yet been standardized, and virulence analysis of strains in vivo using model organisms has not been performed. Our study is devoted to testing the assumption that the development of antibiotic resistance in probiotic bacteria under conditions of selective pressure of antimicrobial drugs may be accompanied by the evolution of virulence. In this regard, special attention is required for the widespread in nature commensals and probiotic bacteria actively used in pharmacology and the food industry. As a result of step-by-step selection from the Lactiplantibacillus plantarum 8p-a3 strain isolated from the “Lactobacterin” probiotic (Biomed, Russia), the L. plantarum 8p-a3-Clr-Amx strain was obtained, showing increased resistance simultaneously to amoxicillin-clavulanic acid and clarithromycin (antibiotics, the combined use of which is widely used for Helicobacter pylori eradication) compared to the parent strain (MIC_{8p-a3-Clr-Amx} of 20 μg/mL and 10 μg/mL, and MIC_8p-a3 of 0.5 μg/mL and 0.05 μg/mL, respectively). The results of a comparative analysis of antibiotic-resistant and parental strains indicate that the development of resistance to the corresponding antimicrobial drugs in L. plantarum in vitro is accompanied by the following: (i) significant changes in the genomic profile (point mutations as well as deletions, insertions, duplications, and displacement of DNA sequences) associated in part with the resistome and mobilome; (ii) changes in phenotypic sensitivity to a number of antimicrobial drugs; and (iii) an increase in the level of virulence against Drosophila melanogaster, a model organism for which L. plantarum is considered to be a symbiont. The data obtained by us indicate that the mechanisms of adaptation to antimicrobial drugs in L. plantarum are not limited to those described earlier and determine the need for comprehensive studies of antibiotic resistance scenarios as well as the trajectories of virulence evolution in probiotic bacteria in vivo and in vitro to develop a standardized system for detecting virulent strains of the corresponding microorganisms.

IMPORTANCE Ensuring the safety of the use of probiotics is a top priority. We found that increased resistance to popular antimicrobial drugs in Lactiplantibacillus plantarum is accompanied by significant changes in the genomic profile and phenotypic sensitivity to a number of antimicrobial drugs as well as in the level of virulence of this bacterium against Drosophila. The data obtained in our work indicate that the mechanisms of antibiotic resistance in this bacterium are not limited to those described earlier and determine the need for comprehensive studies of the potential for the evolution of virulence in lactic acid bacteria in vivo and in vitro and to develop a reliable control system to detect virulent strains among probiotics.

Emergence of a VIM-2-producing XDR Pseudomonas aeruginosa ST309 in South America: a comparative genomic analysis

Pseudomonas aeruginosa has been considered a top priority pathogen associated with elevated morbidity and mortality. Worldwide outbreaks have been associated with few high-risk epidemic P. aeruginosa lineages. However, the biological features involved with the persistence and spread of such lineages among clinical settings remain to be unraveled. This study reports the emergence of an XDR ST309 P. aeruginosa in South America/Brazil, more precisely, in the Amazon region. Genomic analyses were performed with 42 complete and draft ST309 genomes, giving insights into its epidemiology, resistome, and mobilome. Heterogeneous distribution of acquired antimicrobial resistance genes among ST309 genomes was observed, which included bla_VIM-2, bla_IMP-15, and qnrVC1, all of them associated with class 1 integrons. The mobilome mining showed the presence of Integrative and Conjugative Elements, transposons, and genomic islands harbouring a huge arsenal of heavy metal resistance determinants, which probably provided adaptive advantages to ST309 lineage.

ecc_finder: A Robust and Accurate Tool for Detecting Extrachromosomal Circular DNA From Sequencing Data

Extrachromosomal circular DNA (eccDNA) has been observed in different species for decades, and more and more evidence shows that this specific type of DNA molecules may play an important role in rapid adaptation. Therefore, characterizing the full landscape of eccDNA has become critical, and there are several protocols for enriching eccDNAs and performing short-read or long-read sequencing. However, there is currently no available bioinformatic tool to identify eccDNAs from Nanopore reads. More importantly, the current tools based on Illumina short reads lack an efficient standardized pipeline notably to identify eccDNA originating from repeated loci and cannot be applied to very large genomes. Here, we introduce a comprehensive tool to solve both of these two issues. Applying ecc_finder to eccDNA-seq data (either mobilome-seq, Circle-Seq and CIDER-seq) from Arabidopsis, human, and wheat (with genome sizes ranging from 120Mb to 17 Gb), we document the improvement of computational time, sensitivity, and accuracy and demonstrate ecc_finder wide applicability and functionality.

Mobile Antimicrobial Resistance Genes in Probiotics

Even though people worldwide tend to consume probiotic products for their beneficial health effects on a daily basis, recently, concerns were outlined regarding the uptake and potential intestinal colonisation of the bacteria that they carry. These bacteria are capable of executing horizontal gene transfer (HGT) which facilitates the movement of various genes, including antimicrobial resistance genes (ARGs), among the donor and recipient bacterial populations. Within our study, 47 shotgun sequencing datasets deriving from various probiotic samples (isolated strains and metagenomes) were bioinformatically analysed. We detected more than 70 ARGs, out of which rpoB mutants conferring resistance to rifampicin, tet(W/N/W) and potentially extended-spectrum beta-lactamase (ESBL) coding TEM-116 were the most common. Numerous ARGs were associated with integrated mobile genetic elements, plasmids or phages promoting the HGT. Our findings raise clinical and public health concerns as the consumption of probiotic products may lead to the transfer of ARGs to human gut bacteria.

RNAi suppression of DNA methylation affects the drought stress response and genome integrity in transgenic poplar

Trees are long-lived organisms that continuously adapt to their environments, a process in which epigenetic mechanisms are likely to play a key role. Via downregulation of the chromatin remodeler DECREASED IN DNA METHYLATION 1 (DDM1) in poplar (Populus tremula × Populus alba) RNAi lines, we examined how DNA methylation coordinates genomic and physiological responses to moderate water deficit. We compared the growth and drought response of two RNAi-ddm1 lines to wild-type (WT) trees under well-watered and water deficit/rewatering conditions, and analyzed their methylomes, transcriptomes, mobilomes and phytohormone contents in the shoot apical meristem. The RNAi-ddm1 lines were more tolerant to drought-induced cavitation but did not differ in height or stem diameter growth. About 5000 differentially methylated regions were consistently detected in both RNAi-ddm1 lines, colocalizing with 910 genes and 89 active transposable elements. Under water deficit conditions, 136 differentially expressed genes were found, including many involved in phytohormone pathways; changes in phytohormone concentrations were also detected. Finally, the combination of hypomethylation and drought led to the mobility of two transposable elements. Our findings suggest major roles for DNA methylation in regulation of genes involved in hormone-related stress responses, and the maintenance of genome integrity through repression of transposable elements.

Exploration of the Diversity of Clustered Regularly Interspaced Short Palindromic Repeats-Cas Systems in Clostridium novyi sensu lato

Classified as the genospecies Clostridium novyi sensu lato and distributed into four lineages (I-IV), Clostridium botulinum (group III), Clostridium novyi, and Clostridium haemolyticum are clostridial pathogens that cause animal diseases. Clostridium novyi sensu lato contains a large mobilome consisting of plasmids and circular bacteriophages. Here, we explored clustered regularly interspaced short palindromic repeats (CRISPR) arrays and their associated proteins (Cas) to shed light on the link between evolution of CRISPR-Cas systems and the plasmid and phage composition in a study of 58 Clostridium novyi sensu lato genomes. In 55 of these genomes, types I-B (complete or partial), I-D, II-C, III-B, III-D, or V-U CRISPR-Cas systems were detected in chromosomes as well as in mobile genetic elements (MGEs). Type I-B predominated (67.2%) and was the only CRISPR type detected in the Ia, III, and IV genomic lineages. Putative type V-U CRISPR Cas14a genes were detected in two different cases: next to partial type-IB CRISPR loci on the phage encoding the botulinum neurotoxin (BoNT) in lineage Ia and in 12 lineage II genomes, as part of a putative integrative element related to a phage-inducible chromosomal island (PICI). In the putative PICI, Cas14a was associated with CRISPR arrays and restriction modification (RM) systems as part of an accessory locus. This is the first time a PICI containing such locus has been detected in C. botulinum. Mobilome composition and dynamics were also investigated based on the contents of the CRISPR arrays and the study of spacers. A large proportion of identified protospacers (20.2%) originated from Clostridium novyi sensu lato (p1_Cst, p4_BKT015925, p6_Cst, CWou-2020a, p1_BKT015925, and p2_BKT015925), confirming active exchanges within this genospecies and the key importance of specific MGEs in Clostridium novyi sensu lato.

Integrated Genomic Analyses From Low-Depth Sequencing Help Resolve Phylogenetic Incongruence in the Bamboos (Poaceae: Bambusoideae)

The bamboos (Bambusoideae, Poaceae) comprise a major grass lineage with a complex evolutionary history involving ancient hybridization and allopolyploidy. About 1700 described species are classified into three tribes, Olyreae (herbaceous bamboos), Bambuseae (tropical woody bamboos), and Arundinarieae (temperate woody bamboos). Nuclear analyses strongly support monophyly of the woody tribes, whereas plastome analyses strongly support paraphyly, with Bambuseae sister to Olyreae. Our objectives were to clarify the origin(s) of the woody bamboo tribes and resolve the nuclear vs. plastid conflict using genomic tools. For the first time, plastid and nuclear genomic information from the same bamboo species were combined in a single study. We sampled 51 species of bamboos representing the three tribes, estimated their genome sizes and generated low-depth sample sequence data, from which plastomes were assembled and nuclear repeats were analyzed. The distribution of repeat families was found to agree with nuclear gene phylogenies, but also provides novel insights into nuclear evolutionary history. We infer two early, independent hybridization events, one between an Olyreae ancestor and a woody ancestor giving rise to the two Bambuseae lineages, and another between two woody ancestors giving rise to the Arundinarieae. Retention of the Olyreae plastome associated with differential dominance of nuclear genomes and subsequent diploidization in some lineages explains the paraphyly observed in plastome phylogenetic estimations. We confirm ancient hybridization and allopolyploidy in the origins of the extant woody bamboo lineages and propose biased fractionation and diploidization as important factors in their evolution.

Distinct Retrotransposon Evolution Profile in the Genome of Rabbit (Oryctolagus cuniculus)

Although the rabbit genome has already been annotated, it is mobilome remains largely unknown. Here, multiple pipelines were used to de novo mine and annotate the mobilome in rabbit. Four families and 19 subfamilies of LINE1s, two families and nine subfamilies of SINEs, and 12 ERV families were defined in rabbit based on sequence identity, structural organization, and phylogenetic tree. The analysis of insertion age and polymerase chain reaction suggests that a number of families are very young and may remain active, such as L1B, L1D, OcuSINEA, and OcuERV1. RepeatMasker annotation revealed a distinct transposable element landscape within the genome, with approximately two million copies of SINEs, representing the greatest proportion of the genome (19.61%), followed by LINEs (15.44%), and LTRs (4.11%), respectively, considerably different from most other mammal mobilomes except hedgehog and tree shrew, in which LINEs have the highest proportion. Furthermore, a very high rate of insertion polymorphisms (>85%) for the youngest subfamily (OcuSINEA1) was identified by polymerase chain reaction. The majority of retrotransposon insertions overlapped with protein-coding regions (>80%) and lncRNA (90%) genes. Genomic distribution bias was observed for retrotransposons, with those immediately upstream (-1 kb) and downstream (1 kb) of genes significantly depleted. Local GC content in 50-kb widows had significantly negative correlations with LINE (rs=-0.996) and LTR (rs=-0.829) insertions. The current study revealed a distinct mobilome landscape in rabbit, which will assist in the elucidation of the evolution of the genome of lagomorphs, and even other mammals.

Diverse ATPase Proteins in Mobilomes Constitute a Large Potential Sink for Prokaryotic Host ATP

Prokaryote mobilome genomes rely on host machineries for survival and replication. Given that mobile genetic elements (MGEs) derive their energy from host cells, we investigated the diversity of ATP-utilizing proteins in MGE genomes to determine whether they might be associated with proteins that could suppress related host proteins that consume energy. A comprehensive search of 353 huge phage genomes revealed that up to 9% of the proteins have ATPase domains. For example, ATPase proteins constitute ∼3% of the genomes of Lak phages with ∼550 kbp genomes that occur in the microbiomes of humans and other animals. Statistical analysis shows the number of ATPase proteins increases linearly with genome length, consistent with a large sink for host ATP during replication of megaphages. Using metagenomic data from diverse environments, we found 505 mobilome proteins with ATPase domains fused to diverse functional domains. Among these composite ATPase proteins, 61.6% have known functional domains that could contribute to host energy diversion during the mobilome infection cycle. As many have domains that are known to interact with nucleic acids and proteins, we infer that numerous ATPase proteins are used during replication and for protection from host immune systems. We found a set of uncharacterized ATPase proteins with nuclease and protease activities, displaying unique domain architectures that are energy intensive based on the presence of multiple ATPase domains. In many cases, these composite ATPase proteins genomically co-localize with small proteins in genomic contexts that are reminiscent of toxin-antitoxin systems and phage helicase-antibacterial helicase systems. Small proteins that function as inhibitors may be a common strategy for control of cellular processes, thus could inspire future biochemical experiments for the development of new nucleic acid and protein manipulation tools, with diverse biotechnological applications.

Escherichia coli Causing Recurrent Urinary Tract Infections: Comparison to Non-Recurrent Isolates and Genomic Adaptation in Recurrent Infections

Recurrent urinary tract infection (rUTI) remains a major problem for many women and therefore the pursuit for genomic and phenotypic traits which could define rUTI has been ongoing. The present study applied a genomic approach to investigate recurrent urinary tract infections by comparative analyses of recurrent and non-recurrent Escherichia coli isolates from general practice. From whole-genome sequencing data, phylogenetic clustering and genomic traits were studied on a collection of isolates which caused recurrent infection compared to non-recurrent isolates. In addition, genomic variation between the 1st and following infection was studied on a subset of the isolates. Evidence of limited adaptation between the recurrent infections based on single nucleotide polymorphism analyses with a range of 0-13 non-synonymous single nucleotide polymorphisms (SNPs) between the paired isolates. This included an overrepresentation of SNPs in metabolism genes. We identified several genes which were more common in rUTI isolates, including nine fimbrial genes, however, not significantly after false-discovery rate. Finally, the results show that recurrent isolates of the present dataset are not distinctive by variation in the core genome, and thus, did not cluster distinct from non-rUTI isolates in a SNP phylogeny.

Combining Functional Genomics and Whole-Genome Sequencing to Detect Antibiotic Resistance Genes in Bacterial Strains Co-Occurring Simultaneously in a Brazilian Hospital

(1) Background: The rise of multi-antibiotic resistant bacteria represents an emergent threat to human health. Here, we investigate antibiotic resistance mechanisms in bacteria of several species isolated from an intensive care unit in Brazil. (2) Methods: We used whole-genome analysis to identify antibiotic resistance genes (ARGs) and plasmids in 34 strains of Gram-negative and Gram-positive bacteria, providing the first genomic description of Morganella morganii and Ralstonia mannitolilytica clinical isolates from South America. (3) Results: We identified a high abundance of beta-lactamase genes in resistant organisms, including seven extended-spectrum beta-lactamases (OXA-1, OXA-10, CTX-M-1, KPC, TEM, HYDRO, BLP) shared between organisms from different species. Additionally, we identified several ARG-carrying plasmids indicating the potential for a fast transmission of resistance mechanism between bacterial strains. Furthermore, we uncovered two pairs of (near) identical plasmids exhibiting multi-drug resistance. Finally, since many highly resistant strains carry several different ARGs, we used functional genomics to investigate which of them were indeed functional. In this sense, for three bacterial strains (Escherichia coli, Klebsiella pneumoniae, and M. morganii), we identified six beta-lactamase genes out of 15 predicted in silico as those mainly responsible for the resistance mechanisms observed, corroborating the existence of redundant resistance mechanisms in these organisms. (4) Conclusions: Systematic studies similar to the one presented here should help to prevent outbreaks of novel multidrug-resistant bacteria in healthcare facilities.

Oxidative and radiation stress induces transposable element transcription in Drosophila melanogaster

It has been shown that stressors are capable of activating transposable elements (TEs). Currently, there is a hypothesis that stress activation of TEs may be involved in adaptive evolution, favouring the increase in genetic variability when the population is under adverse conditions. However, TE activation under stress is still poorly understood. In the present study, we estimated the fraction of differentially expressed TEs (DETEs) under ionizing radiation (144, 360 and 864 Gy) and oxidative stress (dioxin, formaldehyde and toluene) treatments. The stress intensity of each treatment was estimated by measuring the number of differentially expressed genes, and we show that several TEs families are activated by stress whereas others are repressed. The proportion of DETEs was positively related to stress intensity. However, even under the strongest stress, only a small fraction of TE families were activated (9.28%) and 17.72% were repressed. Considering all treatments together, the activated proportion was 19.83%. Nevertheless, as several TEs are incomplete or degenerated, only 10.55% of D. melanogaster mobilome is, at same time, activated by the stressors and able to transpose or at least code a protein. Thus, our study points out that although stress activates TEs, it is not a generalized activation process, and for some families, the stress induces repression.

Mobilome Analysis of Achromobacter spp. Isolates from Chronic and Occasional Lung Infection in Cystic Fibrosis Patients

Achromobacter spp. is an opportunistic pathogen that can cause lung infections in patients with cystic fibrosis (CF). Although a variety of mobile genetic elements (MGEs) carrying antimicrobial resistance genes have been identified in clinical isolates, little is known about the contribution of Achromobacter spp. mobilome to its pathogenicity. To provide new insights, we performed bioinformatic analyses of 54 whole genome sequences and investigated the presence of phages, insertion sequences (ISs), and integrative and conjugative elements (ICEs). Most of the detected phages were previously described in other pathogens and carried type II toxin-antitoxin systems as well as other pathogenic genes. Interestingly, the partial sequence of phage Bcep176 was found in all the analyzed Achromobacter xylosoxidans genome sequences, suggesting the integration of this phage in an ancestor strain. A wide variety of IS was also identified either inside of or in proximity to pathogenicity islands. Finally, ICEs carrying pathogenic genes were found to be widespread among our isolates and seemed to be involved in transfer events within the CF lung. These results highlight the contribution of MGEs to the pathogenicity of Achromobacter species, their potential to become antimicrobial targets, and the need for further studies to better elucidate their clinical impact.

Genomic and Transcriptomic Survey Provides New Insight into the Organization and Transposition Activity of Highly Expressed LTR Retrotransposons of Sunflower (Helianthus annuus L.)

LTR retrotransposons (RTEs) play a crucial role in plant genome evolution and adaptation. Although RTEs are generally silenced in somatic plant tissues under non-stressed conditions, some expressed RTEs (exRTEs) escape genome defense mechanisms. As our understanding of exRTE organization in plants is rudimentary, we systematically surveyed the genomic and transcriptomic organization and mobilome (transposition) activity of sunflower (Helianthus annuus L.) exRTEs. We identified 44 transcribed RTEs in the sunflower genome and demonstrated their distinct genomic features: more recent insertion time, longer open reading frame (ORF) length, and smaller distance to neighboring genes. We showed that GAG-encoding ORFs are present at significantly higher frequencies in exRTEs, compared with non-expressed RTEs. Most exRTEs exhibit variation in copy number among sunflower cultivars and one exRTE Gagarin produces extrachromosomal circular DNA in seedling, demonstrating recent and ongoing transposition activity. Nanopore direct RNA sequencing of full-length RTE RNA revealed complex patterns of alternative splicing in RTE RNAs, resulting in isoforms that carry ORFs for distinct RTE proteins. Together, our study demonstrates that tens of expressed sunflower RTEs with specific genomic organization shape the hidden layer of the transcriptome, pointing to the evolution of specific strategies that circumvent existing genome defense mechanisms.

Genomic and Resistance Epidemiology of Gram-Negative Bacteria in Africa: a Systematic Review and Phylogenomic Analyses from a One Health Perspective

Antibiotic resistance (AR) is one of the major public health threats and challenges to effective containment and treatment of infectious bacterial diseases worldwide. Here, we used different methods to map out the geographical hot spots, sources, and evolutionary epidemiology of AR. Escherichia coli, Klebsiella pneumoniae, Salmonella enterica, Acinetobacter baumannii, Pseudomonas aeruginosa, Enterobacter spp., Neisseria meningitis/gonorrhoeae, Vibrio cholerae, Campylobacter jejuni, etc., were common pathogens shuttling AR genes in Africa. Transmission of the same clones/strains across countries and between animals, humans, plants, and the environment was observed. We recommend Enterobacter spp. or K. pneumoniae as better sentinel species for AR surveillance.

Antibiotic resistance (AR) remains a major threat to public and animal health globally. However, AR ramifications in developing countries are worsened by limited molecular diagnostics, expensive therapeutics, inadequate numbers of skilled clinicians and scientists, and unsanitary environments. The epidemiology of Gram-negative bacteria, their AR genes, and geographical distribution in Africa are described here. Data were extracted and analyzed from English-language articles published between 2015 and December 2019. The genomes and AR genes of the various species, obtained from the Pathosystems Resource Integration Center (PATRIC) and NCBI were analyzed phylogenetically using Randomized Axelerated Maximum Likelihood (RAxML) and annotated with Figtree. The geographic location of resistant clones/clades was mapped manually. Thirty species from 31 countries and 24 genera from 41 countries were analyzed from 146 articles and 3,028 genomes, respectively. Genes mediating resistance to β-lactams (including bla_TEM-1, bla_CTX-M, bla_NDM, bla_IMP, bla_VIM, and bla_OXA-48/181), fluoroquinolones (oqxAB, qnrA/B/D/S, gyrA/B, and parCE mutations, etc.), aminoglycosides (including armA and rmtC/F), sulfonamides (sul1/2/3), trimethoprim (dfrA), tetracycline [tet(A/B/C/D/G/O/M/39)], colistin (mcr-1), phenicols (catA/B, cmlA), and fosfomycin (fosA) were mostly found in Enterobacter spp. and Klebsiella pneumoniae, and also in Serratia marcescens, Escherichia coli, Salmonella enterica, Pseudomonas, Acinetobacter baumannii, etc., on mostly IncF-type, IncX_3/4, ColRNAI, and IncR plasmids, within IntI1 gene cassettes, insertion sequences, and transposons. Clonal and multiclonal outbreaks and dissemination of resistance genes across species and countries and between humans, animals, plants, and the environment were observed; Escherichia coli ST103, K. pneumoniae ST101, S. enterica ST1/2, and Vibrio cholerae ST69/515 were common strains. Most pathogens were of human origin, and zoonotic transmissions were relatively limited.

IMPORTANCE Antibiotic resistance (AR) is one of the major public health threats and challenges to effective containment and treatment of infectious bacterial diseases worldwide. Here, we used different methods to map out the geographical hot spots, sources, and evolutionary epidemiology of AR. Escherichia coli, Klebsiella pneumoniae, Salmonella enterica, Acinetobacter baumannii, Pseudomonas aeruginosa, Enterobacter spp., Neisseria meningitis/gonorrhoeae, Vibrio cholerae, Campylobacter jejuni, etc., were common pathogens shuttling AR genes in Africa. Transmission of the same clones/strains across countries and between animals, humans, plants, and the environment was observed. We recommend Enterobacter spp. or K. pneumoniae as better sentinel species for AR surveillance.

Whole Genome Sequencing Reveals Virulence Potentials of Helicobacter pylori Strain KE21 Isolated from a Kenyan Patient with Gastric Signet Ring Cell Carcinoma

Helicobacter pylori (H.pylori) infection is etiologically associated with severe diseases including gastric cancer; but its pathogenicity is deeply shaped by the exceptional genomic diversification and geographic variation of the species. The clinical relevance of strains colonizing Africa is still debated. This study aimed to explore genomic features and virulence potentials of H. pylori KE21, a typical African strain isolated from a native Kenyan patient diagnosed with a gastric cancer. A high-quality circular genome assembly of 1,648,327 bp (1590 genes) obtained as a hybrid of Illumina Miseq short reads and Oxford Nanopore MinION long reads, clustered within hpAfrica1 population. This genome revealed a virulome and a mobilome encoding more than hundred features potentiating a successful colonization, persistent infection, and enhanced disease pathogenesis. Furthermore, through an experimental infection of gastric epithelial cell lines, strain KE21 showed the ability to promote interleukin-8 production and to induce cellular alterations resulting from the injection of a functional CagA oncogene protein into the cells. This study shows that strain KE21 is potentially virulent and can trigger oncogenic pathways in gastric epithelial cells. Expended genomic and clinical explorations are required to evaluate the epidemiological importance of H. pylori infection and its putative complications in the study population.

Comparative genomics of the fish pathogens Edwardsiella ictaluri 93-146 and Edwardsiella piscicida C07-087

Edwardsiella ictaluri and Edwardsiella piscicida are important fish pathogens affecting cultured and wild fish worldwide. To investigate the genome-level differences and similarities between catfish-adapted strains in these two species, the complete E. ictaluri 93-146 and E. piscicida C07-087 genomes were evaluated by applying comparative genomics analysis. All available complete (10) and non-complete (19) genomes from five Edwardsiella species were also included in a systematic analysis. Average nucleotide identity and core-genome phylogenetic tree analyses indicated that the five Edwardsiella species were separated from each other. Pan-/core-genome analyses for the 29 strains from the five species showed that genus Edwardsiella members have 9474 genes in their pan genome, while the core genome consists of 1421 genes. Orthology cluster analysis showed that E. ictaluri and E. piscicida genomes have the greatest number of shared clusters. However, E. ictaluri and E. piscicida also have unique features; for example, the E. ictaluri genome encodes urease enzymes and cytochrome o ubiquinol oxidase subunits, whereas E. piscicida genomes encode tetrathionate reductase operons, capsular polysaccharide synthesis enzymes and vibrioferrin-related genes. Additionally, we report for what is believed to be the first time that E. ictaluri 93-146 and three other E. ictaluri genomes encode a type IV secretion system (T4SS), whereas none of the E. piscicida genomes encode this system. Additionally, the E. piscicida C07-087 genome encodes two different type VI secretion systems. E. ictaluri genomes tend to encode more insertion elements, phage regions and genomic islands than E. piscicida. We speculate that the T4SS could contribute to the increased number of mobilome elements in E. ictaluri compared to E. piscicida. Two of the E. piscicida genomes encode full CRISPR-Cas regions, whereas none of the E. ictaluri genomes encode Cas proteins. Overall, comparison of the E. ictaluri and E. piscicida genomes reveals unique features and provides new insights on pathogenicity that may reflect the host adaptation of the two species.

Genomic analysis of a multidrug-resistant clinical Providencia rettgeri (PR002) strain with the novel integron ln1483 and an A/C plasmid replicon

Whole-genome sequence analysis was performed on a multidrug-resistant Providencia rettgeri PR002 clinical strain isolated from the urine of a hospitalized patient in Pretoria, South Africa, in 2013. The resistome, mobilome, pathogenicity island(s), as well as virulence and heavy-metal resistance genes of the isolate, were characterized using whole-genome sequencing and bioinformatic analysis. PR002 had a genome assembly size of 4,832,624 bp with a GC content of 40.7%, an A/C₂ plasmid replicase gene, four integrons/gene cassettes, 17 resistance genes, and several virulence and heavy metal resistance genes, confirming PR002 as a human pathogen. A novel integron, In1483, harboring the gene bla_OXA-2 , was identified, with other uncharacterized class 1 integrons harboring aacA4cr and dfrA1. Aac(3')-IIa and bla_SCO-1 , as well as bla_PER-7 , sul2, and tet(B), were found bracketed by composite Tn3 transposons, and IS91, IS91, and IS4 family insertion sequences, respectively. PR002 was resistant to all antibiotics tested except amikacin, carbapenems, cefotaxime-clavulanate, ceftazidime-clavulanate, cefoxitin, and fosfomycin. PR002 was closely related to PR1 (USA), PRET_2032 (SPAIN), DSM_1131, and NCTC7477 clinical P. rettgeri strains, but not close enough to suggest it was imported into South Africa from other countries. Multidrug resistance in P. rettgeri is rare, particularly in clinical settings, making this case an important incident requiring urgent attention. This is also the first report of an A/C plasmid in P. rettgeri. The array, multiplicity, and diversity of resistance and virulence genes in this strain are concerning, necessitating stringent infection control, antibiotic stewardship, and periodic resistance surveillance/monitoring policies to preempt further horizontal and vertical spread of these resistance genes.

Genomic Analysis of Shewanella sp. O23S-The Natural Host of the pSheB Plasmid Carrying Genes for Arsenic Resistance and Dissimilatory Reduction

Shewanella sp. O23S is a dissimilatory arsenate reducing bacterial strain involved in arsenic transformations within the abandoned gold mine in Zloty Stok (SW Poland). Previous physiological studies revealed that O23S may not only release arsenic from minerals, but also facilitate its immobilization through co-precipitation with reduced sulfur species. Given these uncommon, complementary characteristics and the application potential of the strain in arsenic-removal technologies, its genome (~5.3 Mbp), consisting of a single chromosome, two large plasmids (pSheA and pSheB) and three small plasmid-like phages (pSheC-E) was sequenced and annotated. Genes encoding putative proteins involved in heavy metal transformations, antibiotic resistance and other phenotypic traits were identified. An in-depth comparative analysis of arsenic respiration (arr) and resistance (ars) genes and their genetic context was also performed, revealing that pSheB carries the only copy of the arr genes, and a complete ars operon. The plasmid pSheB is therefore a unique natural vector of these genes, providing the host cells arsenic respiration and resistance abilities. The functionality of the identified genes was determined based on the results of the previous and additional physiological studies, including: the assessment of heavy metal and antibiotic resistance under various conditions, adhesion-biofilm formation assay and Biolog^TM metabolic preferences test. This combined genetic and physiological approach shed a new light on the capabilities of O23S and their molecular basis, and helped to confirm the biosafety of the strain in relation to its application in bioremediation technologies.

Into the Thermus Mobilome: Presence, Diversity and Recent Activities of Insertion Sequences Across Thermus spp

A high level of transposon-mediated genome rearrangement is a common trait among microorganisms isolated from thermal environments, probably contributing to the extraordinary genomic plasticity and horizontal gene transfer (HGT) observed in these habitats. In this work, active and inactive insertion sequences (ISs) spanning the sequenced members of the genus Thermus were characterized, with special emphasis on three T. thermophilus strains: HB27, HB8, and NAR1. A large number of full ISs and fragments derived from different IS families were found, concentrating within megaplasmids present in most isolates. Potentially active ISs were identified through analysis of transposase integrity, and domestication-related transposition events of ISTth7 were identified in laboratory-adapted HB27 derivatives. Many partial copies of ISs appeared throughout the genome, which may serve as specific targets for homologous recombination contributing to genome rearrangement. Moreover, recruitment of IS1000 32 bp segments as spacers for CRISPR sequence was identified, pointing to the adaptability of these elements in the biology of these thermophiles. Further knowledge about the activity and functional diversity of ISs in this genus may contribute to the generation of engineered transposons as new genetic tools, and enrich our understanding of the outstanding plasticity shown by these thermophiles.

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.

Evolution of Sphingomonad Gene Clusters Related to Pesticide Catabolism Revealed by Genome Sequence and Mobilomics of Sphingobium herbicidovorans MH

Bacterial degraders of chlorophenoxy herbicides have been isolated from various ecosystems, including pristine environments. Among these degraders, the sphingomonads constitute a prominent group that displays versatile xenobiotic-degradation capabilities. Four separate sequencing strategies were required to provide the complete sequence of the complex and plastic genome of the canonical chlorophenoxy herbicide-degrading Sphingobium herbicidovorans MH. The genome has an intricate organization of the chlorophenoxy-herbicide catabolic genes sdpA, rdpA, and cadABCD that encode the (R)- and (S)-enantiomer-specific 2,4-dichlorophenoxypropionate dioxygenases and four subunits of a Rieske non-heme iron oxygenase involved in 2-methyl-chlorophenoxyacetic acid degradation, respectively. Several major genomic rearrangements are proposed to help understand the evolution and mobility of these important genes and their genetic context. Single-strain mobilomic sequence analysis uncovered plasmids and insertion sequence-associated circular intermediates in this environmentally important bacterium and enabled the description of evolutionary models for pesticide degradation in strain MH and related organisms. The mobilome presented a complex mosaic of mobile genetic elements including four plasmids and several circular intermediate DNA molecules of insertion-sequence elements and transposons that are central to the evolution of xenobiotics degradation. Furthermore, two individual chromosomally integrated prophages were shown to excise and form free circular DNA molecules. This approach holds great potential for improving the understanding of genome plasticity, evolution, and microbial ecology.

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.

The CRISPR Spacer Space Is Dominated by Sequences from Species-Specific Mobilomes

Clustered regularly interspaced short palindromic repeats and CRISPR-associated protein (CRISPR-Cas) systems store the memory of past encounters with foreign DNA in unique spacers that are inserted between direct repeats in CRISPR arrays. For only a small fraction of the spacers, homologous sequences, called protospacers, are detectable in viral, plasmid, and microbial genomes. The rest of the spacers remain the CRISPR “dark matter.” We performed a comprehensive analysis of the spacers from all CRISPR-cas loci identified in bacterial and archaeal genomes, and we found that, depending on the CRISPR-Cas subtype and the prokaryotic phylum, protospacers were detectable for 1% to about 19% of the spacers (~7% global average). Among the detected protospacers, the majority, typically 80 to 90%, originated from viral genomes, including proviruses, and among the rest, the most common source was genes that are integrated into microbial chromosomes but are involved in plasmid conjugation or replication. Thus, almost all spacers with identifiable protospacers target mobile genetic elements (MGE). The GC content, as well as dinucleotide and tetranucleotide compositions, of microbial genomes, their spacer complements, and the cognate viral genomes showed a nearly perfect correlation and were almost identical. Given the near absence of self-targeting spacers, these findings are most compatible with the possibility that the spacers, including the dark matter, are derived almost completely from the species-specific microbial mobilomes.

The principal function of CRISPR-Cas systems is thought to be protection of bacteria and archaea against viruses and other parasitic genetic elements. The CRISPR defense function is mediated by sequences from parasitic elements, known as spacers, that are inserted into CRISPR arrays and then transcribed and employed as guides to identify and inactivate the cognate parasitic genomes. However, only a small fraction of the CRISPR spacers match any sequences in the current databases, and of these, only a minority correspond to known parasitic elements. We show that nearly all spacers with matches originate from viral or plasmid genomes that are either free or have been integrated into the host genome. We further demonstrate that spacers with no matches have the same properties as those of identifiable origins, strongly suggesting that all spacers originate from mobile elements.

Phylogenomic, Pan-genomic, Pathogenomic and Evolutionary Genomic Insights into the Agronomically Relevant Enterobacteria Pantoea ananatis and Pantoea stewartii

Pantoea ananatis is ubiquitously found in the environment and causes disease on a wide range of plant hosts. By contrast, its sister species, Pantoea stewartii subsp. stewartii is the host-specific causative agent of the devastating maize disease Stewart's wilt. This pathogen has a restricted lifecycle, overwintering in an insect vector before being introduced into susceptible maize cultivars, causing disease and returning to overwinter in its vector. The other subspecies of P. stewartii subsp. indologenes, has been isolated from different plant hosts and is predicted to proliferate in different environmental niches. Here we have, by the use of comparative genomics and a comprehensive suite of bioinformatic tools, analyzed the genomes of ten P. stewartii and nineteen P. ananatis strains. Our phylogenomic analyses have revealed that there are two distinct clades within P. ananatis while far less phylogenetic diversity was observed among the P. stewartii subspecies. Pan-genome analyses revealed a large core genome comprising of 3,571 protein coding sequences is shared among the twenty-nine compared strains. Furthermore, we showed that an extensive accessory genome made up largely by a mobilome of plasmids, integrated prophages, integrative and conjugative elements and insertion elements has resulted in extensive diversification of P. stewartii and P. ananatis. While these organisms share many pathogenicity determinants, our comparative genomic analyses show that they differ in terms of the secretion systems they encode. The genomic differences identified in this study have allowed us to postulate on the divergent evolutionary histories of the analyzed P. ananatis and P. stewartii strains and on the molecular basis underlying their ecological success and host range.

Proteomics technique opens new frontiers in mobilome research

A large proportion of the genome of most eukaryotic organisms consists of highly repetitive mobile genetic elements. The sum of these elements is called the "mobilome," which in eukaryotes is made up mostly of transposons. Transposable elements contribute to disease, evolution, and normal physiology by mediating genetic rearrangement, and through the "domestication" of transposon proteins for cellular functions. Although 'omics studies of mobilome genomes and transcriptomes are common, technical challenges have hampered high-throughput global proteomics analyses of transposons. In a recent paper, we overcame these technical hurdles using a technique called "proteomics informed by transcriptomics" (PIT), and thus published the first unbiased global mobilome-derived proteome for any organism (using cell lines derived from the mosquito Aedes aegypti). In this commentary, we describe our methods in more detail, and summarise our major findings. We also use new genome sequencing data to show that, in many cases, the specific genomic element expressing a given protein can be identified using PIT. This proteomic technique therefore represents an important technological advance that will open new avenues of research into the role that proteins derived from transposons and other repetitive and sequence diverse genetic elements, such as endogenous retroviruses, play in health and disease.

Analysis of the Genome and Mobilome of a Dissimilatory Arsenate Reducing Aeromonas sp. O23A Reveals Multiple Mechanisms for Heavy Metal Resistance and Metabolism

Aeromonas spp. are among the most ubiquitous microorganisms, as they have been isolated from different environmental niches including waters, soil, as well as wounds and digestive tracts of poikilothermic animals and humans. Although much attention has been paid to the pathogenicity of Aeromonads, the role of these bacteria in environmentally important processes, such as transformation of heavy metals, remains to be discovered. Therefore, the aim of this study was a detailed genomic characterization of Aeromonas sp. O23A, the first representative of this genus capable of dissimilatory arsenate reduction. The strain was isolated from microbial mats from the Zloty Stok mine (SW Poland), an environment strongly contaminated with arsenic. Previous physiological studies indicated that O23A may be involved in both mobilization and immobilization of this metalloid in the environment. To discover the molecular basis of the mechanisms behind the observed abilities, the genome of O23A (∼5.0 Mbp) was sequenced and annotated, and genes for arsenic respiration, heavy metal resistance (hmr) and other phenotypic traits, including siderophore production, were identified. The functionality of the indicated gene modules was assessed in a series of minimal inhibitory concentration analyses for various metals and metalloids, as well as mineral dissolution experiments. Interestingly, comparative analyses revealed that O23A is related to a fish pathogen Aeromonas salmonicida subsp. salmonicida A449 which, however, does not carry genes for arsenic respiration. This indicates that the dissimilatory arsenate reduction ability may have been lost during genome reduction in pathogenic strains, or acquired through horizontal gene transfer. Therefore, particular emphasis was placed upon the mobilome of O23A, consisting of four plasmids, a phage, and numerous transposable elements, which may play a role in the dissemination of hmr and arsenic metabolism genes in the environment. The obtained results indicate that Aeromonas sp. O23A is well-adapted to the extreme environmental conditions occurring in the Zloty Stok mine. The analysis of genome encoded traits allowed for a better understanding of the mechanisms of adaptation of the strain, also with respect to its presumable role in colonization and remediation of arsenic-contaminated waters, which may never have been discovered based on physiological analyses alone.

A Pan-Genomic Approach to Understand the Basis of Host Adaptation in Achromobacter

Over the past decade, there has been a rising interest in Achromobacter sp., an emerging opportunistic pathogen responsible for nosocomial and cystic fibrosis lung infections. Species of this genus are ubiquitous in the environment, can outcompete resident microbiota, and are resistant to commonly used disinfectants as well as antibiotics. Nevertheless, the Achromobacter genus suffers from difficulties in diagnosis, unresolved taxonomy and limited understanding of how it adapts to the cystic fibrosis lung, not to mention other host environments. The goals of this first genus-wide comparative genomics study were to clarify the taxonomy of this genus and identify genomic features associated with pathogenicity and host adaptation. This was done with a widely applicable approach based on pan-genome analysis. First, using all publicly available genomes, a combination of phylogenetic analysis based on 1,780 conserved genes with average nucleotide identity and accessory genome composition allowed the identification of a largely clinical lineage composed of Achromobacter xylosoxidans, Achromobacter insuavis, Achromobacter dolens, and Achromobacter ruhlandii. Within this lineage, we identified 35 positively selected genes involved in metabolism, regulation and efflux-mediated antibiotic resistance. Second, resistome analysis showed that this clinical lineage carried additional antibiotic resistance genes compared with other isolates. Finally, we identified putative mobile elements that contribute 53% of the genus's resistome and support horizontal gene transfer between Achromobacter and other ecologically similar genera. This study provides strong phylogenetic and pan-genomic bases to motivate further research on Achromobacter, and contributes to the understanding of opportunistic pathogen evolution.

The Expanding Family of Virophages

Virophages replicate with giant viruses in the same eukaryotic cells. They are a major component of the specific mobilome of mimiviruses. Since their discovery in 2008, five other representatives have been isolated, 18 new genomes have been described, two of which being nearly completely sequenced, and they have been classified in a new viral family, Lavidaviridae. Virophages are small viruses with approximately 35–74 nm large icosahedral capsids and 17–29 kbp large double-stranded DNA genomes with 16–34 genes, among which a very small set is shared with giant viruses. Virophages have been isolated or detected in various locations and in a broad range of habitats worldwide, including the deep ocean and inland. Humans, therefore, could be commonly exposed to virophages, although currently limited evidence exists of their presence in humans based on serology and metagenomics. The distribution of virophages, the consequences of their infection and the interactions with their giant viral hosts within eukaryotic cells deserve further research.

Genomic and metagenomic technologies to explore the antibiotic resistance mobilome

Antibiotic resistance is a relevant problem for human health that requires global approaches to establish a deep understanding of the processes of acquisition, stabilization, and spread of resistance among human bacterial pathogens. Since natural (nonclinical) ecosystems are reservoirs of resistance genes, a health-integrated study of the epidemiology of antibiotic resistance requires the exploration of such ecosystems with the aim of determining the role they may play in the selection, evolution, and spread of antibiotic resistance genes, involving the so-called resistance mobilome. High-throughput sequencing techniques allow an unprecedented opportunity to describe the genetic composition of a given microbiome without the need to subculture the organisms present inside. However, bioinformatic methods for analyzing this bulk of data, mainly with respect to binning each resistance gene with the organism hosting it, are still in their infancy. Here, we discuss how current genomic methodologies can serve to analyze the resistance mobilome and its linkage with different bacterial genomes and metagenomes. In addition, we describe the drawbacks of current methodologies for analyzing the resistance mobilome, mainly in cases of complex microbiotas, and discuss the possibility of implementing novel tools to improve our current metagenomic toolbox.

Evolution and Diversity of the Antimicrobial Resistance Associated Mobilome in Streptococcus suis: A Probable Mobile Genetic Elements Reservoir for Other Streptococci

Streptococcus suis is a previously neglected, newly emerging multidrug-resistant zoonotic pathogen. Mobile genetic elements (MGEs) play a key role in intra- and interspecies horizontal transfer of antimicrobial resistance (AMR) determinants. Although, previous studies showed the presence of several MGEs, a comprehensive analysis of AMR-associated mobilome as well as their interaction and evolution has not been performed. In this study, we presented the AMR-associated mobilome and their insertion hotspots in S. suis. Integrative conjugative elements (ICEs), prophages and tandem MGEs were located at different insertion sites, while 86% of the AMR-associated MGEs were inserted at rplL and rum loci. Comprehensive analysis of insertions at rplL and rum loci among four pathogenic Streptococcus species (Streptococcus agalactiae, Streptococcus pneumoniae, Streptococcus pyogenes, and S. suis) revealed the existence of different groups of MGEs, including Tn5252, ICESp1108, and TnGBS2 groups ICEs, Φm46.1 group prophage, ICE_ICE and ICE_prophage tandem MGEs. Comparative ICE genomics of ICESa2603 family revealed that module exchange and acquisition/deletion were the main mechanisms in MGEs' expansion and evolution. Furthermore, the observation of tandem MGEs reflected a novel mechanism for MGE diversity. Moreover, an in vitro competition assay showed no visible fitness cost was observed between different MGE-carrying isolates and a conjugation assay revealed the transferability of ICESa2603 family of ICEs. Our statistics further indicated that the prevalence and diversity of MGEs in S. suis is much greater than in other three species which prompted our hypothesis that S. suis is probably a MGEs reservoir for other streptococci. In conclusion, our results showed that acquisition of MGEs confers S. suis not only its capability as a multidrug resistance pathogen, but also represents a paradigm to study the modular evolution and matryoshkas of MGEs.