Search

Chlamydia are Gram-negative, obligate intracellular bacterial pathogens that infect eukaryotic cells and reside within a host-derived vacuole known as the inclusion. To facilitate intracellular replication, these bacteria must engage in host-pathogen interactions to obtain nutrients and membranes required for the growth of the inclusion, thereby sustaining prolonged bacterial colonization. Autophagy is a highly conserved process that delivers cytoplasmic substrates to the lysosome for degradation. Pathogens have developed strategies to manipulate and/or exploit autophagy to promote their replication and persistence. This review delineates recent advances in elucidating the interplay between Chlamydia trachomatis infection and autophagy in recent years, emphasizing the intricate strategies employed by both the Chlamydia pathogens and host cells. Gaining a deeper understanding of these interactions could unveil novel strategies for the prevention and treatment of Chlamydia infection.

Abstract Objectives Pseudomonas aeruginosa has intrinsic antibiotic resistance and the strong ability to acquire additional resistance genes. However, a limited number of investigations provide detailed modular structure dissection and evolutionary analysis of accessory genetic elements (AGEs) and associated resistance genes (ARGs) in P. aeruginosa isolates. The objective of this study is to reveal the prevalence and transmission characteristics of ARGs by epidemiological investigation and bioinformatics analysis of AGEs of P. aeruginosa isolates taken from a Chinese hospital. Methods Draft-genome sequencing was conducted for P. aeruginosa clinical isolates (n = 48) collected from a single Chinese hospital between 2019 and 2021. The clones of P. aeruginosa isolates, type 3 secretion system (T3SS)-related virulotypes, and the resistance spectrum were identified using multilocus sequence typing (MLST), polymerase chain reaction (PCR), and antimicrobial susceptibility tests. In addition, 17 of the 48 isolates were fully sequenced. An extensive modular structure dissection and genetic comparison was applied to AGEs of the 17 sequenced P. aeruginosa isolates. Results From the draft-genome sequencing, 13 STs were identified, showing high genetic diversity. BLAST search and PCR detection of T3SS genes (exoT, exoY, exoS, and exoU) revealed that the exoS+/exoU- virulotype dominated. At least 69 kinds of acquired ARGs, involved in resistance to 10 different categories of antimicrobials, were identified in the 48 P. aeruginosa isolates. Detailed genetic dissection and sequence comparisons were applied to 25 AGEs from the 17 isolates, together with five additional prototype AGEs from GenBank. These 30 AGEs were classified into five groups -- integrative and conjugative elements (ICEs), unit transposons, IncpPBL16 plasmids, Incp60512−IMP plasmids, and IncpPA7790 plasmids. Conclusion This study provides a broad-scale and deeper genomics understanding of P. aeruginosa isolates taken from a single Chinese hospital. The isolates collected are characterized by high genetic diversity, high virulence, and multiple drug resistance. The AGEs in P. aeruginosa chromosomes and plasmids, as important genetic platforms for the spread of ARGs, contribute to enhancing the adaptability of P. aeruginosa in hospital settings.

Providencia as an opportunistic pathogen can cause serious infection, and moreover the emergence of multi-drug-resistant Providencia strains poses a potentially life-threatening risk to public health. However, a comprehensive genomic study to reveal the population structure and dissemination of Providencia is still lacking. In this study, we conducted a genomic epidemiology analysis on the 580 global sequenced Providencia isolates, including 257 ones sequenced in this study (42 ones were fully sequenced). We established a genome sequence-based species classification scheme for Providencia, redefining the conventional 11 Providencia species into seven genocomplexes that were further divided into 18 genospecies, providing an extensively updated reference for Providencia species discrimination based on the largest Providencia genome dataset to date. We then dissected the profile of antimicrobial resistance genes and the prevalence of multi-drug-resistant Providencia strains among these genocomplexes/genospecies, disclosing the presence of diverse and abundant antimicrobial resistance genes and high resistance ratios against multiple classes of drugs in Providencia. We further dissected the genetic basis for the spread of blaNDM-1 in Providencia. blaNDM-1 genes were mainly carried by five incompatible (Inc) groups of plasmids: IncC, IncW, IncpPROV114-NR, IncpCHS4.1-3, and IncpPrY2001, and the last three were newly designated in this study. By tracking the spread of blaNDM-1-carrying plasmids, IncC, IncpPROV114-NR, IncpCHS4.1-3, and IncpPrY2001 plasmids were found to be highly involved in parallel horizontal transfer or vertical clonal expansion of blaNDM-1 among Providencia. Overall, our study provided a comprehensive genomic view of species differentiation, antimicrobial resistance prevalence, and plasmid-mediated blaNDM-1 dissemination in Providencia.

ABSTRACT Three novel imipenemase (IMP)-type metallo-β-lactamases (MBLs), referred to as IMP-89, IMP-91, and IMP-96, were detected in three clinical isolates from China. Antimicrobial susceptibility tests indicated these novel enzymes were resistant to most β-lactams, and IMP-96 with a Ser262Gly mutation had higher activity against meropenem than its point mutant. We then collected sequence data on all 91 available IMP variants for phylogenetic analysis. To further analyze the genetic environment of blaIMP, an extensive comparison was applied to nine accessory genetic elements (AGEs), including six sequenced blaIMP-carrying AGEs in this study and three others from GenBank. These nine AGEs were divided into three groups: three IncpJBCL41 plasmids, Tn6417 and its two derivatives, and three Tn6879-related integrative and conjugative elements (ICEs). All blaIMP genes in this study were captured by class 1 integrons. In the integrons, blaIMP genes usually coexisted with other resistance genes, which further impeded clinical antibacterial treatment. The emergence of new IMP variants and the diversity and complexity of their genetic environment make the prevention and control of drug-resistant strains critical, requiring serious attention from clinical and public health management departments. IMPORTANCE The spread of IMP-type MBLs has increased dramatically in recent years. We discovered three novel IMP variants from three clinical isolates in China. We summarized the classification and evolutionary relationship of all available IMP variants. Moreover, we detailed the genetic characteristics of blaIMP-carrying accessory genetic elements in five clinical isolates. Given the risk of rapid and extensive spread of blaIMP genes, we suggest that continuous surveillance is crucial to combat the acquisition and transmission of blaIMP genes by bacteria, which can impede clinical therapy effectiveness.

Resistance and tolerance are two important strategies employed by the host immune response to defend against pathogens. Multidrug-resistant bacteria affect the resistance mechanisms involved in pathogen clearance. Disease tolerance, defined as the ability to reduce the negative impact of infection on the host, might be a new research direction for the treatment of infections. The lungs are highly susceptible to infections and thus are important for understanding host tolerance and its precise mechanisms. This review focuses on the factors that induce lung disease tolerance, cell and molecular mechanisms involved in tissue damage control, and the relationship between disease tolerance and sepsis immunoparalysis. Understanding the exact mechanism of lung disease tolerance could allow better assessment of the immune status of patients and provide new ideas for the treatment of infections.

Mature biofilm formation by Vibrio parahaemolyticus requires exopolysaccharide (EPS), type IV pili, and capsular polysaccharide (CPS). Production of each is strictly regulated by various control pathways including quorum sensing (QS) and bis-(3′–5′)-cyclic di-GMP (c-di-GMP). QsvR, an AraC-type regulator, integrates into the QS regulatory cascade via direct control of the transcription of the master QS regulators, AphA and OpaR. Deletion of qsvR in wild-type or opaR mutant backgrounds altered the biofilm formation by V. parahaemolyticus, suggesting that QsvR may coordinate with OpaR to control biofilm formation. Herein, we demonstrated both QsvR and OpaR repressed biofilm-associated phenotypes, c-di-GMP metabolism, and the formation of V. parahaemolyticus translucent (TR) colonies. QsvR restored the biofilm-associated phenotypic changes caused by opaR mutation, and vice versa. In addition, QsvR and OpaR worked coordinately to regulate the transcription of EPS-associated genes, type IV pili genes, CPS genes and c-di-GMP metabolism-related genes. These results demonstrated how QsvR works with the QS system to regulate biofilm formation by precisely controlling the transcription of multiple biofilm formation-associated genes in V. parahaemolyticus.

Abstract Background The carbapenem-resistance genes bla VIM are widely disseminated in Pseudomonas, and frequently harbored within class 1 integrons that reside within various mobile genetic elements (MGEs). However, there are few reports on detailed genetic dissection of bla VIM-carrying MGEs in Pseudomonas. Methods This study presented the complete sequences of five bla VIM-2/-4-carrying MGEs, including two plasmids, two chromosomal integrative and mobilizable elements (IMEs), and one chromosomal integrative and conjugative element (ICE) from five different Pseudomonas isolates. Results The two plasmids were assigned to a novel incompatibility (Inc) group IncpSTY, which included only seven available plasmids with determined complete sequences and could be further divided into three subgroups IncpSTY-1/2/3. A detailed sequence comparison was then applied to a collection of 15 MGEs belonging to four different groups: three representative IncpSTY plasmids, two Tn6916-related IMEs, two Tn6918-related IMEs, and eight Tn6417-related ICEs and ten of these 15 MGEs were first time identified. At least 22 genes involving resistance to seven different categories of antibiotics and heavy metals were identified within these 15 MGEs, and most of these resistance genes were located within the accessory modules integrated as exogenous DNA regions into these MGEs. Especially, eleven of these 15 MGEs carried the bla VIM genes, which were located within 11 different concise class 1 integrons. Conclusion These bla VIM-carrying integrons were further integrated into the above plasmids, IMEs/ICEs with intercellular mobility. These MGEs could transfer between Pseudomonas isolates, which resulted in the accumulation and spread of bla VIM among Pseudomonas and thus was helpful for the bacteria to survival from the stress of antibiotics. Data presented here provided a deeper insight into the genetic diversification and evolution of VIM-encoding MGEs in Pseudomonas.

ABSTRACT Vibrio parahaemolyticus, a causative agent of seafood-associated gastroenteritis, undergoes opaque-translucent (OP-TR) colony switching associated with capsular polysaccharide (CPS) production. Here, we showed that V. parahaemolyticus was also able to naturally and reversibly switch between wrinkly and smooth phenotypes. More than 1,000 genes were significantly differentially expressed during colony morphology switching, including the major virulence gene loci and key biofilm-related genes. The genes responsible for type III secretion system 1 (T3SS1), type VI secretion systems (T6SS1 and T6SS2), and flagellar synthesis were downregulated in the wrinkly spreader phenotype, whereas genes located on the pathogenicity island Vp-PAI and those responsible for chitin-regulated pili (ChiRP) and Syp exopolysaccharide synthesis were upregulated. In addition, we showed that the wrinkly spreader grew faster, had greater motility and biofilm capacities, and produced more c-di-GMP than the smooth type. A dozen genes potentially associated with c-di-GMP metabolism were shown to be significantly differentially expressed, which may account for the differences in c-di-GMP levels between the two phenotypes. Most importantly, dozens of putative regulators were significantly differentially expressed, and hundreds of noncoding RNAs were detected during colony morphology switching, indicating that phenotype switching is strictly regulated by a complex molecular regulatory network in V. parahaemolyticus. Taken together, the presented work highlighted the gene expression profiles related to wrinkly-smooth switching, showing that the significantly differentially expressed genes were involved in various biological behaviors, including virulence factor production, biofilm formation, metabolism, adaptation, and colonization. IMPORTANCE We showed that Vibrio parahaemolyticus was able to naturally and reversibly switch between wrinkly and smooth phenotypes and disclosed the gene expression profiles related to wrinkly-smooth switching, showing that the significantly differentially expressed genes between the two colony morphology phenotypes were involved in various biological behaviors, including virulence factor production, biofilm formation, metabolism, adaptation, and colonization.

As a newly emerging Klebsiella pathogen, more and more Klebsiella michiganensis drug resistant strains have been reported in recent years, which posed serious threats to public health. Here we first reported a multidrug-resistant K. michiganensis strain 12084 with two blaSIM-1 and one mcr-9.2 genes isolated from the sputum specimen of a patient in the Second Affiliated Hospital of Zhejiang University School of Medicine and analyzed its genetic basis and drug-resistance phenotypes. Genetic analysis showed that this strain harbored three different incompatibility groups (IncHI2, IncHI5, and IncFIIpKPHS2:IncFIB-4.1) of plasmids (p12084-HI2, p12084-HI5, and p12084-FII). A total of 26 drug-resistance genes belonging to 12 classes of antibiotics were identified, most of which (24) were located on two plasmids (p12084-HI2 and p12084-HI5). Interestingly, two blaSIM-1 genes were identified to locate on p12084-HI2 and p12084-HI5, respectively, both of which were embedded in In630, indicating their genetic homogeny. It was noting that one blaSIM-1 gene was situated in a novel unit transposon (referred to as Tn6733) on the p12084-HI5 plasmid. We also discovered an mcr-9.2 gene on the p12084-HI2 plasmid. To the best of our knowledge, this is the first report of a blaSIM-1 and mcr-9.2 harboring K. michiganensis strain. We then investigated the population structure/classification, and antibiotic resistance for all 275 availably global K. michiganensis genomes. Population structure revealed that K. michiganensis could be divided into two main clades (Clade 1 and Clade 2); the most popular ST29 was located in Clade 1, while other common STs (such as ST50, ST27, and ST43) were located in Clade 2. Drug-resistance analysis showed 25.5% of the K. michiganensis strains (70/275) harboring at least one carbapenemase gene, indicating severe drug resistance of K. michiganensis beyond our imagination; this is a dangerous trend and should be closely monitored, especially for ST27 K. michiganensis with the most drug-resistant genes among all the STs. Overall, we reported a blaSIM-1 and mcr-9.2 harboring K. michiganensis strain, and further revealed the population structure/classification, and drug-resistance of K. michiganensis, which provided an important framework, reference, and improved understanding of K. michiganensis.

USA300, a dominant clone of community-acquired methicillin-resistant Staphylococcus aureus (CA-MRSA), is circulating globally and can cause necrotizing pneumonia with high morbidity and mortality. To further reveal the host anti-MRSA infection immune response, we established a mouse model of acute primary MRSA pneumonia challenged with aerosols of the USA300 clone. A time-course transcriptome analysis of the lungs collected at 0, 12, 24, 48 and 96 h post-infection (hpi) was conducted using RNA sequencing (RNA-seq) and multiple bioinformatic analysis methods. The change trend of histopathology and five innate immune cell (neutrophils, mononuclear cells, eosinophils, macrophages, DC cells) proportions in the lungs after infection was also examined. We observed a distinct acute pulmonary recovery process. A rapid initiation period of inflammation was present at 12 hpi, during which the IL-17 pathway dominantly mediated inflammation and immune defense. The main stages of host inflammatory response occurred at 24 and 48 hpi, and the regulation of interferon activation and macrophage polarization played an important role in the control of inflammatory balance at this stage. At 96 hpi, cellular proliferation processes associated with host repair were observed, as well as adaptive immunity and complement system responses involving C1q molecules. More importantly, the data provide new insight into and identify potential functional genes involved in the checks and balances occurring between host anti-inflammatory and proinflammatory responses. To the best of our knowledge, this is the first study to investigate transcriptional responses throughout the inflammatory recovery process in the lungs after MRSA infection. Our study uncovers valuable research targets for key regulatory mechanisms underlying the pathogenesis of MRSA lung infections, which may help to develop novel treatment strategies for MRSA pneumonia.

Hypervirulent Klebsiella pneumoniae (hvKp) can cause life-threatening community-acquired infections among healthy young individuals and is thus of concern for global dissemination. In this study, a mouse model of acute primary hvKp pneumonia was established via aerosolized intratracheal (i.t.) inoculation, laying the foundation for conducting extensive studies related to hvKp. Subsequently, a time-course transcriptional profile was created of the lungs from the mouse model at 0, 12, 24, 48 and 60 hours post-infection (hpi) using RNA Sequencing (RNA-Seq). RNA-Seq data were analyzed with the use of Mfuzz time clustering, weighted gene co-expression network analysis (WGCNA) and Immune Cell Abundance Identifier for mouse (ImmuCellAI-mouse). A gradual change in the transcriptional profile of the lungs was observed that reflected expected disease progression. At 12 hpi, genes related to acute phase inflammatory response increased in expression and lipid metabolism appeared to have a pro-inflammatory effect. At 24 hpi, exacerbation of inflammation was observed and active IFN-γ suggested that signaling promoted activation and recruitment of macrophages occurred. Genes related to maintaining the structural integrity of lung tissues showed a sustained decrease in expression after infection and the decrease was especially marked at 48 hpi. TNF, IL-17, MAPK and NF-kB signaling pathways may play key roles in the immunopathogenesis mechanism at all stages of infection. Natural killer (NK) cells consistently decreased in abundance after infection, which has rarely been reported in hvKp infection and could provide a new target for treatment. Genes Saa1 and Slpi were significantly upregulated during infection. Both Saa1, which is associated with lipopolysaccharide (LPS) that elicits host inflammatory response, and Slpi, which encodes an antimicrobial protein, have not previously been reported in hvKp infections and could be important targets for subsequent studies. To t our knowledge, this paper represents the first study to investigate the pulmonary transcriptional response to hvKp infection. The results provide new insights into the molecular mechanisms underlying the pathogenesis of hvKp pulmonary infection that can contribute to the development of therapies to reduce hvKp pneumonia.

ABSTRACT In this study, draft-genome sequencing was conducted for 60 Chinese Morganella isolates, and furthermore, 12 of them were fully sequenced. Then, a total of 166 global sequenced Morganella isolates, including the above 60, were collected to perform average nucleotide identity-based genomic classification and core single nucleotide polymorphism-based phylogenomic analysis. A genome sequence-based species classification scheme for Morganella was established, and accordingly, the two conventional Morganella species were redefined as two complexes and further divided into four and two genospecies, respectively. At least 88 acquired antimicrobial resistance genes (ARGs) were disseminated in these 166 isolates and were prevalent mostly in the isolates from hospital settings. IS26/IS15DI, IS10 and IS1R, and Tn3-, Tn21-, and Tn7-subfamily unit transposons were frequently presented in these 166 isolates. Furthermore, a detailed sequence comparison was applied to 18 Morganella chromosomal accessory genetic elements (AGEs) from the fully sequenced 12 isolates, together with 5 prototype AGEs from GenBank. These 23 AGEs were divided into eight different groups belonging to composite/unit transposons, transposable prophages, integrative and mobilizable elements, and integrative and conjugative elements, and they harbored at least 52 ARGs involved in resistance to 15 categories of antimicrobials. Eleven of these 23 AGEs acquired large accessory modules, which exhibited complex mosaic structures and contained many antimicrobial resistance loci and associated ARGs. Integration of ARG-containing AGEs into Morganella chromosomes would contribute to the accumulation and dissemination of ARGs in Morganella and enhance the adaption and survival of Morganella under complex and diverse antimicrobial selection pressures. IMPORTANCE This study presents a comprehensive genomic epidemiology analysis on global sequenced Morganella isolates. First, a genome sequence-based species classification scheme for Morganella is established with a higher resolution and accuracy than those of the conventional scheme. Second, the prevalence of accessory genetic elements (AGEs) and associated antimicrobial resistance genes (ARGs) among Morganella isolates is disclosed based on genome sequences. Finally, a detailed sequence comparison of eight groups of 23 AGEs (including 19 Morganella chromosomal AGEs) reveals that Morganella chromosomes have evolved to acquire diverse AGEs harboring different profiles of ARGs and that some of these AGEs harbor large accessory modules that exhibit complex mosaic structures and contain a large number of ARGs. Data presented here provide a deeper understanding of the classification and evolution of Morganella species and also those of ARG-containing AGEs in Morganella at the genomic scale.

Abstract Background Coxiella burnetii is an obligate intracellular Gram-negative bacterium that causes a zoonotic disease commonly called Q fever globally. In this study, an up-converting phosphor technology-based lateral flow (UPT-LF) assay was established for the rapid and specific detection of phase I strains of C. burnetii. Results Specific monoclonal antibodies (10B5 and 10G7) against C. burnetii phase I strains were prepared and selected for use in the UPT-LF assay by the double-antibody-sandwich method. The detection sensitivity of the Coxiella-UPT-LF was 5 × 104 GE/ml for a purified C. burnetii phase I strain and 10 ng/ml for LPS of C. burnetii Nine Mile phase I (NMI). Good linearity was observed for C. burnetii phase I and NMI LPS quantification (R2 ≥ 0.989). The UPT-LF assay also exhibited a high specificity to C. burnetii, without false-positive results even at 108 GE/ml of non-specific bacteria, and good inclusivity for detecting different phase I strains of C. burnetii. Moreover, the performance of the Coxiella-UPT-LF assay was further confirmed using experimentally and naturally infected samples. Conclusions Our results indicate that Coxiella-UPT-LF is a sensitive and reliable method for rapid screening of C. burnetii, suitable for on-site detection in the field.

ABSTRACTKlebsiella pneumoniae carbapenemase (KPC)-producing K. pneumoniae (KPC-KP) have disseminated worldwide and emerged as major threats to public health. Of epidemiological significance, the international pandemic of KPC-KP is primarily associated with CG258 isolates and blaKPC-IncF plasmids. CRISPR-Cas system is an adaptive immune system that can hinder gene expansion driven by horizontal gene transfer. Because of blaKPC-IncF plasmids are favored by CG258 K. pneumoniae, it was of interest to examine the co-distribution of CRISPR and blaKPC-IncF plasmids in such isolates. We collected 459 clinical K. pneumoniae isolates in China and collected 203 global whole-genome sequences in GenBank to determine the prevalence of CRISPR-Cas systems. We observed that CRISPR-Cas system was significantly scarce in the CG258 lineage and blaKPC-positive isolates. Furthermore, the results of conjugation and plasmid stability assay fully demonstrated the CRIPSR-Cas system in K. pneumoniae could effectively hindered blaKPC-IncF plasmids invasion and existence. Notably, most blaKPC-IncF plasmids were also proved to be good targets of CRISPR owing to carry matched and functional protospacers and PAMs. Overall, our work suggests that type I-E CRISPR-Cas systems could impact the spread of blaKPC in K. pneumoniae populations, and the scarcity of CRISPR-Cas system was one of potential factors leading to the propagation of blaKPC-IncF plasmids in CG258 K. pneumoniae.

Vibrio parahaemolyticus, the leading cause of seafood-associated gastroenteritis worldwide, has a strong ability to form biofilms on surfaces. Quorum sensing (QS) is a process widely used by bacteria to communicate with each other and control gene expression via the secretion and detection of autoinducers. OpaR is the master QS regulator of V. parahaemolyticus operating under high cell density (HCD). OpaR regulation of V. parahaemolyticus biofilm formation has been reported, but the regulatory mechanisms are still not fully understood. bis-(3′-5′)-cyclic di-GMP (c-di-GMP) is an omnipresent intracellular second messenger that regulates diverse behaviors of bacteria including activation of biofilm formation. In this work, we showed that OpaR repressed biofilm formation and decreased the intracellular concentration of c-di-GMP in V. parahaemolyticus RIMD2210633. The OpaR box-like sequences were detected within the regulatory DNA regions of scrA, scrG, VP0117, VPA0198, VPA1176, VP0699, and VP2979, encoding a group of GGDEF and/or EAL-type proteins. The results of qPCR, LacZ fusion, EMSA, and DNase I footprinting assays demonstrated that OpaR bound to the upstream DNA regions of scrA, VP0117, VPA0198, VPA1176, and VP0699 to repress their transcription, whereas it positively and directly regulated the transcription of scrG and VP2979. Thus, transcriptional regulation of these genes by OpaR led directly to changes in the intracellular concentration of c-di-GMP. The direct association between QS and c-di-GMP metabolism in V. parahaemolyticus RIMD2210633 would be conducive to precise control of gene transcription and bacterial behaviors such as biofilm formation.

In this study, a detailed genetic dissection of the huge and complex blaNDM-carrying genetic elements and their related mobile genetic elements was performed in Enterobacteriaceae. An extensive comparison was applied to 12 chromosomal genetic elements, including six sequenced in this study and the other six from GenBank. These 12 genetic elements were divided into five groups: a novel IME Tn6588; two related IMEs Tn6523 (SGI1) and Tn6589; four related ICEs Tn6512 (R391), Tn6575 (ICEPvuChnBC22), Tn6576, and Tn6577; Tn7 and its derivatives Tn6726 and 40.7-kb Tn7-related element; and two related IMEs Tn6591 (GIsul2) and Tn6590. At least 51 resistance genes, involved in resistance to 18 different categories of antibiotics and heavy metals, were found in these 12 genetic elements. Notably, Tn6576 carried another ICE Tn6582. In particular, the six blaNDM-carrying genetic elements Tn6588, Tn6589, Tn6575, Tn6576, Tn6726, and 40.7-kb Tn7-related element contained large accessory multidrug resistance (MDR) regions, each of which had a very complex mosaic structure that comprised intact or residual mobile genetic elements including insertion sequences, unit or composite transposons, integrons, and putative resistance units. Core blaNDM genetic environments manifested as four different Tn125 derivatives and, notably, two or more copies of relevant Tn125 derivatives were found in each of Tn6576, Tn6588, Tn6589, and 40.7-kb Tn7-related element. The huge and complex blaNDM-carrying genetic elements were assembled from complex transposition and homolog recombination. Firstly identified were eight novel mobile elements, including three ICEs Tn6576, Tn6577, and Tn6582, two IMEs, Tn6588 and Tn6589, two composite transposons Tn6580a and Tn6580b, and one integron In1718.

More than a decade of gut microbiome studies have a common goal of improving human health. However, while most of the disease studies have focused on the elderly or the middle-aged, a reference cohort for the gut microbiome in young individuals has been lacking. It is also not clear what other omics data need to be measured to better understand the gut microbiome. Here, we present a cohort including 2183 adults with high-depth metagenomic shotgun sequencing data for the fecal microbiome and other omics data. In this multiomic cohort, we observe a number of vitamins, hormones, amino acids, and trace elements that correlated with the gut microbiome. Many of the associations are validated in an additional cohort consisting of 1404 individuals. Our comprehensive data are poised to provide advice to future populations and mechanistic study designs to better understand and manage our gut microbiome.

Pseudomonas aeruginosa is noted for its intrinsic antibiotic resistance and capacity of acquiring additional resistance genes. In this study, the genomes of nine clinical P. aeruginosa isolates were fully sequenced. An extensive genetic comparison was applied to 18 P. aeruginosa accessory genetic elements (AGEs; 13 of them were sequenced in this study and located within P. aeruginosa chromosomes) that were divided into four groups: five related integrative and conjugative elements (ICEs), four related integrative and mobilizable elements (IMEs), five related unit transposons, and two related IMEs and their two derivatives. At least 45 resistance genes, involved in resistance to 10 different categories of antibiotics and heavy metals, were identified from these 18 AGEs. A total of 10 β-lactamase genes were identified from 10 AGEs sequenced herein, and nine of them were captured within class 1 integrons, which were further integrated into ICEs and IMEs with intercellular mobility, and also unit transposons with intracellular mobility. Through this study, we identified for the first time 20 novel MGEs, including four ICEs Tn6584, Tn6585, Tn6586, and Tn6587; three IMEs Tn6853, Tn6854, and Tn6878; five unit transposons Tn6846, Tn6847, Tn6848, Tn6849, and Tn6883; and eight integrons In1795, In1778, In1820, In1784, In1775, In1774, In1789, and In1799. This was also the first report of two resistance gene variants blaCARB-53 and catB3s, and a novel ST3405 isolate of P. aeruginosa. The data presented here denoted that complex transposition and homologous recombination promoted the assembly and integration of AGEs with mosaic structures into P. aeruginosa chromosomes.

Plague, which is caused by Yersinia pestis, is one of the most dangerous infectious diseases. No FDA-approved vaccine against plague is available for human use at present. To improve the immune safety of Y. pestis EV76 based live attenuated vaccine and to explore the feasibility of aerosolized intratracheal inoculation (i.t.) route for vaccine delivery, a plasminogen activator protease (pla) gene deletion mutant of the attenuated Y. pestis strain EV76-B-SHU was constructed, and its residual virulence and protective efficacy were evaluated in a mouse model via aerosolized intratracheal inoculation (i.t.) or via subcutaneous injection (s.c.). The residual virulence of EV76-B-SHUΔpla was significantly reduced compared to that of the parental strain EV76-B-SHU following i.t. and s.c. infection. The EV76-B-SHUΔpla induced higher levels of mucosal antibody sIgA in the bronchoalveolar lavage fluid of mice immunized by i.t. but not by s.c.. Moreover, after lethal challenge with Y. pestis biovar Microtus strain 201 (avirulent in humans), the protective efficacy and bacterial clearance ability of the EV76-B-SHUΔpla-i.t. group were comparable to those of the EV76-B-SHUΔpla-s.c. and EV76-B-SHU immunized groups. Thus, the EV76-B-SHUΔpla represents an excellent live-attenuated vaccine candidate against pneumonic plague and aerosolized i.t. represents a promising immunization route in mouse model.

Mobile gene elements play an important role in the continuous evolution of the prophage DNA of bacteria, promoting the emergence of new gene structures. This study explored the evolution of four strains of Klebsiella pneumoniae harboring prophages, 19051, 721005, 911021, and 675920, and 16 genomes of K. pneumoniae from GenBank. The results revealed a wide range of genetic variation in the prophage DNA inserted into the sap sites of K. pneumoniae chromosomes. From analysis and comparison of the sequences of the 20 prophage DNAs determined from the four strains and the 16 GenBank genomes of K. pneumoniae using high-throughput sequencing and antimicrobial susceptibility tests, we identified a novel transposon, Tn6556. We also identified at least nine large genetic structures with massive genetic acquisitions or losses and five hotspot sites showing a tendency to undergo insertion of gene elements such as IS1T, IS1R, IS26, ISKpn26, ISKpn28, Tn6556, MDR, and In27-related regions as variable regions; however, the only highly conserved core genes were int and umuCD among the 20 prophage DNAs. These findings provide important insights into the evolutionary diversity of bacteriophage DNA contained in K. pneumoniae.

Three different MDR plasmids p16005813A, p16005813B, and p16005813C, which carried a total of 18 non-redundant resistance genes or gene loci, were identified in a single clinical isolate of Leclercia adecarboxylata. The p16005813A backbone showed very low levels of identity to all DNA sequences available in public databases and carried a repA gene that could not assigned into any of known incompatibility groups. The IncFII-family p16005813B and pECAZ161_KPC had essentially identical backbones. p16005813C belonged to an IncR single-replicon plasmid. p16005813A, p16005813B, and p16005813C harbored three different novel MDR regions as their sole accessory modules. The MDR region of p16005813B manifested as Tn6505, which was generated from insertion of blaIMP–8-carrying In655 instead of In4 into the Tn1696 backbone. Other key antibiotic resistance elements included Tn2, IS26–mph(A)–mrx–mphR(A)–IS6100 unit, chrA region, In27, and aacC2–tmrB region in the MDR region of p16005813A, and ΔTn9 carrying catA1, In609, and IS26–tetA(C)–tetR(C)–IS26 unit in the MDR region of p16005813C. This was the first report of coexistence of three different MDR plasmids, and that of occurrence of IMP-encoding plasmid and blaIMP–8 gene in L. adecarboxylata.

A collection of 11 IncC plasmids from China were fully sequenced herein and compared with reference plasmids pR148 and pR55. These 13 plasmids could be assigned into three different subgroups: type 1, type 2, and type 1/2 hybrid. Type 1/2-hybrid plasmids most likely emerged from homologous recombination between type 1 and type 2 plasmids. Different IncC plasmids had evolved to acquire quite different profiles of accessory modules and thus different collections of resistance genes. The accessory resistance modules included not only the blaCMY-carrying region, the ARI-A island, and the ARI-B island, but also various additional kinds of resistance islands such as the blaCTX–M-carrying regions and the MDR regions. Insertion of accessory modules was sometimes accompanied by deletion, inversion, and translocation of surrounding backbone regions. pR148 and pR55 were confirmed to have the most complete backbones for type 1 and type 2, respectively. This was the first report of a blaIMP–8-carrying IncC plasmid, and that of three novel mobile elements: a Tn1696-derived unit transposon Tn6395, a class 2 integron In2-76, and an insertion sequence ISEcl10.

This study presents three novel integrons In1394, In1395, and In1443, three novel unit transposons Tn6392, Tn6393, and Tn6403, one novel conjugative element (ICE) Tn6413, and the first sequenced IncP-7 resistance plasmid p1160-VIM from clinical Pseudomonas aeruginosa. Detailed sequence comparison of p1160-VIM (carrying Tn6392 and Tn6393) and Tn6413 (carrying Tn6403) with related elements were performed. Tn6392, Tn6393, and Tn6403 were generated from integration of In1394 (carrying blaVIM–24), In1395 and In1443 (carrying blaVIM–4) into prototype Tn3-family unit transposons Tn5563, Tn1403, and Tn6346, respectively. To the best of our knowledge, this is the first report of a blaVIM–24-carrying P. aeruginosa isolate.

IncI-complex plasmids can be divided into seven subgroups IncI1, IncI2, IncIγ, IncB/O, IncK1, IncK2, and IncZ. In this study, a replicon-based scheme was proposed for typing IncI-complex plasmids into four separately clustering subgroups IncI2, IncI1/B/O, IncIγ/K1 and IncK2/Z, the last three of which were combined from IncI1 and IncB/O, IncIγ and IncK1, and IncK2 and IncZ, respectively. Four IncIγ/K1 plasmids p205880-NR2, p14E509-CTXM, p11011-CTXM and p61806-CTXM were fully sequenced and compared with IncIγ/K1 reference pCT, IncI2 reference R721, IncI1/B/O reference R64 and IncK2/Z reference pO26-CRL-125. These plasmids shared conserved gene organization in the replication and conjugal transfer regions, but displaying considerable sequence diversity among different subgroups. Remarkable modular differences were observed in the maintenance and transfer leading regions. p205880-NR2 contained no resistance genes or accessory modules, while the other seven plasmids acquired one or more accessory modules, which harbored mobile elements [including unit transposons, insertion sequence (IS)-based transposition units and individual IS elements] and associated resistance markers (especially including those involved in resistance to β-lactams, aminoglycosides, tetracyclins, phenicols, streptomycins, trimethoprims, sulphonamides, tunicamycins and erythromycins). Data presented here provided a deeper insight into diversification and evolution of IncI-complex plasmids.

IncHI plasmids could be divided into five different subgroups IncHI1–5. In this study, the complete nucleotide sequences of seven blaIMP- or blaVIM-carrying IncHI5 plasmids from Klebsiella pneumoniae, K. quasipneumoniae, and K. variicola were determined and compared in detail with all the other four available sequenced IncHI5 plasmids. These plasmids carried conserved IncHI5 backbones composed of repHI5B and a repFIB-like gene (replication), parABC (partition), and tra1 (conjugal transfer). Integration of a number of accessory modules, through horizontal gene transfer, at various sites of IncHI5 backbones resulted in various deletions of surrounding backbone regions and thus considerable diversification of IncHI5 backbones. Among the accessory modules were three kinds of resistance accessory modules, namely Tn10 and two antibiotic resistance islands designated ARI-A and ARI-B. These two islands, inserted at two different fixed sites (one island was at one site and the other was at a different site) of IncHI5 backbones, were derived from the prototype Tn3-family transposons Tn1696 and Tn6535, respectively, and could be further discriminated as various intact transposons and transposon-like structures. The ARI-A or ARI-B islands from different IncHI5 plasmids carried distinct profiles of antimicrobial resistance markers and associated mobile elements, and complex events of transposition and homologous recombination accounted for assembly of these islands. The carbapenemase genes blaIMP-4, blaIMP-38 and blaVIM-1 were identified within various class 1 integrons from ARI-A or ARI-B of the seven plasmids sequenced in this study. Data presented here would provide a deeper insight into diversification and evolution history of IncHI5 plasmids.

Biofilm formation is critical for blocking flea foregut and hence for transmission of Y. pestis by flea biting. In this study, we identified the regulatory role of the AraC-family transcriptional regulator BfvR (YPO1737 in strain CO92) in biofilm formation and virulence of Yersinia pestis biovar Microtus. Crystal violet staining, Caenorhabditis elegans biofilm assay, colony morphology assay, intracellular c-di-GMP concentration determination, and BALB/c mice challenge were employed to reveal that BfvR enhanced Y. pestis biofilm formation while repressed its virulence in mice. Further molecular biological assays demonstrated that BfvR directly stimulated the expression of hmsHFRS, waaAE-coaD, and hmsCDE, which, in turn, affected the production of exopolysaccharide, LPS, and c-di-GMP, respectively. In addition, BfvR directly and indirectly repressed psaABC and psaEF transcription, respectively. We concluded that the modulation of biofilm- and virulence-related genes by BfvR led to increased biofilm formation and reduced virulence of Y. pestis biovar Microtus.

Vibrio parahaemolyticus, the leading causative agent of seafood-associated gastroenteritis, harbors two major virulence gene loci T3SS1 and Vp-PAI (T3SS2 and tdh2). ToxR is a virulence regulator of vibrios. Cell density-dependent transcriptional pattern of toxR and its regulatory actions on T3SS1 and Vp-PAI have been previously reported, but the detailed regulatory mechanisms are still obscure. In the present work, we showed that the highest transcription level of toxR occurs at an OD600 = 0.2–0.4, which may be due to the subtle repression of ToxR and the quorum-sensing (QS) master regulator AphA. We also showed that ToxR is involved in regulating the mouse lethality, enterotoxicity, cytotoxicity, and hemolytic activity of V. parahaemolyticus. ToxR binds to the multiple promoter-proximal DNA regions within the T3SS1 locus to repress their transcription. In addition, ToxR occupies the multiple promoter-proximal DNA regions of Vp-PAI locus to activate their transcription. Thus, ToxR regulates the multiple virulence phenotypes via directly acting on the T3SS1 and Vp-PAI genes. Data presented here provide a deeper understanding of the regulatory patterns of ToxR in V. parahaemolyticus.

Forty-five KPC-producing Enterobacteriaceae strains were isolated from multiple departments in a Chinese public hospital from 2014 to 2015. Genome sequencing of four representative strains, namely Proteus mirabilis GN2, Serratia marcescens GN26, Morganella morganii GN28, and Klebsiella aerogenes E20, indicated the presence of blaKPC-2-carrying IncX6 plasmids pGN2-KPC, pGN26-KPC, pGN28-KPC, and pE20-KPC in the four strains, respectively. These plasmids were genetically closely related to one another and to the only previously sequenced IncX6 plasmid, pKPC3_SZ. Each of the plasmids carried a single accessory module containing the blaKPC-2/3-carrying ΔTn6296 derivatives. The ΔTn6292 element from pGN26-KPC also contained qnrS, which was absent from all other plasmids. Overall, pKPC3_SZ-like blaKPC-carrying IncX6 plasmids were detected by PCR in 44.4% of the KPC-producing isolates, which included K. aerogenes, P. mirabilis, S. marcescens, M. morganii, Escherichia coli, and Klebsiella pneumoniae, and were obtained from six different departments of the hospital. Data presented herein provided insights into the genomic diversity and evolution of IncX6 plasmids, as well as the dissemination and epidemiology of blaKPC-carrying IncX6 plasmids among Enterobacteriaceae in a hospital setting.

Vibrio parahaemolyticus, the leading cause of seafood-associated gastroenteritis, harbors two separate T6SSs on chromosomes 1 and 2, i.e., T6SS1 (VP1386-1420) and T6SS2 (VPA1025-1046). T6SS1 contains at least 7 putative operons: VP1386-1387, VP1388-1390, VP1392-1391, VP1393-1406, VP1400-1406, VP1409-1407, and VP1410-1420. V. parahaemolyticus AphA and OpaR are the two master regulators of quorum sensing (QS) system that are highly expressed at low cell density and high cell density, respectively. ToxR is a membrane-bound virulence regulatory protein conserved across the Vibrio family. In the present work, we show that ToxR coordinates with AphA and OpaR to repress T6SS1 expression in V. parahaemolyticus. OpaR binds to the promoters of VP1388-1390, VP1400-1406, and VP1409-1407 to repress their transcription, but it appears to negatively regulate VP1393-1406 transcription in an indirect manner. By contrast, AphA negatively regulated the above four T6SS1 operons in an indirect manner. In addition, ToxR binds to the promoters of VP1400-1406 and VP1409-1407 to inhibit their transcription, but it presents an indirect interaction with VP1388-1390 and VP1393-1406 promoters. Notably, the expression of ToxR also manifested in a QS-dependent manner and the highest expression occurred at LCD. Meanwhile, the highest expression of T6SS1 occurred at an OD600 value of 0.6 to 0.8 due to the tight regulation of ToxR and QS, suggesting T6SS1 functions only during the mid-logarithmic growth phase. These observations provide significant insight into the molecular mechanism of T6SS1 gene regulation by QS and ToxR in V. parahaemolyticus.

Novel class 1 integrons In1085 and In1086, containing the class D β-lactamase -encoding gene blaOXA, were identified in clinical enterobacterial strains. In this study, we aimed to characterize the genetic contexts of In1085 and In1086, with the goal of identifying putative mechanisms of integron mobilization. Four plasmids, approximately 5.3, 5.3, 5.7, and 6.6 kb, from 71 clinical Enterobacteriaceae strains were found to contain class 1 integrons (In37, In62, In1085, and In1086, respectively). Two of these plasmids, pEco336 and pNsa292, containing In1085 and In1086, respectively, were further characterized by antibiotic susceptibility testing, conjugation experiments, PCR, sequencing, and gene mapping. The OXA-type carbapenemase activities of the parental strains were also assessed. The results revealed that the novel integrons had different genetic environments, and therefore demonstrated diverse biochemical characteristics. Using evolutionary inferences based on the recombination of gene cassettes, we also identified a role for attCaadA16 structural features during attI1 × attC insertion reactions. Our analysis showed that gene cassette insertions in the bottom strand of attCaadA16 in the correct orientation lead to the expression the encoded genes from the Pc promoter. Our study suggests that the genetic features harbored within the integrons are inserted in a discernable pattern, involving the stepwise and parallel evolution of class 1 integron variations under antibiotic selection pressures in a clinical setting.

The purpose of this study was to characterize mechanisms of plasmid-mediated antimicrobial resistance in Shigella boydii. S. boydii strain 2246 with resistance to ciprofloxacin, ceftriaxone and azithromycin was isolated from a human case of watery diarrhea in a Chinese public hospital. Resistance in strain 2246 to ceftriaxone and azithromycin was attributable to the presence of blaCTX-M-14, and erm(B) and mph(A), respectively, which were co-located on a multidrug-resistant (MDR) plasmid p2246-CTXM. p2246-CTXM represented a novel IncFII-type MDR plasmid with a very complex chimera structure. Its master backbone was genetically closely related to the R100 plasmid, but p2246-CTXM had evolved to integrate additional R100-unrelated backbone regions as well as massive exogenous mobile elements that carried multiple resistance determinants. In p2246-CTXM, erm(B) together with its leading peptide gene erm(C), mph(A) together with its regulatory genes mrx and mphR(A), and blaCTX-M-14 were captured by three different mobile elements Tn6295, the IS26-mph(A)-mrx-mphR(A)-IS6100 unit, and a truncated ISEcp1-blaCTX-M-14-IS903D-iroN transposition unit, respectively, all of which were harbored in a large Tn3-family transposon Tn6285. p2246-CTXM still carried additional resistance determinants mer (mercury resistance), aacA4 (aminoglycoside resistance), cmlA1 (chloramphenicol resistance) and qacED1 (quaternary ammonium compound resistance). This is the first report of identifying a clinical S. boydii strain simultaneously resistant to ciprofloxacin, ceftriaxone and azithromycin, and determining the complete sequence of a resistance plasmid from S. boydii.