FastMLST: A Multi-core Tool for Multilocus Sequence Typing of Draft Genome Assemblies

Characterization of acquired β-lactamases in Pseudomonas aeruginosa and quantification of their contributions to resistance

Pseudomonas aeruginosa is a highly problematic opportunistic pathogen that causes a range of different infections. Infections are commonly treated with β-lactam antibiotics, including cephalosporins, monobactams, penicillins, and carbapenems, with carbapenems regarded as antibiotics of last resort. Isolates of P. aeruginosa can contain horizontally acquired bla genes encoding β-lactamase enzymes, but the extent to which these contribute to β-lactam resistance in this species has not been systematically quantified. The overall aim of this research was to address this knowledge gap by quantifying the frequency of β-lactamase-encoding genes in P. aeruginosa and by determining the effects of β-lactamases on susceptibility of P. aeruginosa to β-lactams. Genome analysis showed that β-lactamase-encoding genes are present in 3% of P. aeruginosa but are enriched in carbapenem-resistant isolates (35%). To determine the substrate antibiotics, 10 β-lactamases were expressed from an integrative plasmid in the chromosome of P. aeruginosa reference strain PAO1. The β-lactamases reduced susceptibility to a variety of clinically used antibiotics, including carbapenems (meropenem, imipenem), penicillins (ticarcillin, piperacillin), cephalosporins (ceftazidime, cefepime), and a monobactam (aztreonam). Different enzymes acted on different β-lactams. β-lactamases encoded by the genomes of P. aeruginosa clinical isolates had similar effects to the enzymes expressed in strain PAO1. Genome engineering was used to delete β-lactamase-encoding genes from three carbapenem-resistant clinical isolates and increased susceptibility to substrate β-lactams. Our findings demonstrate that acquired β-lactamases play an important role in β-lactam resistance in P. aeruginosa, identifying substrate antibiotics for a range of enzymes and quantifying their contributions to resistance.IMPORTANCEPseudomonas aeruginosa is an extremely problematic pathogen, with isolates that are resistant to the carbapenem class of β-lactam antibiotics being in critical need of new therapies. Genes encoding β-lactamase enzymes that degrade β-lactam antibiotics can be present in P. aeruginosa, including carbapenem-resistant isolates. Here, we show that β-lactamase genes are over-represented in carbapenem-resistant isolates, indicating their key role in resistance. We also show that different β-lactamases alter susceptibility of P. aeruginosa to different β-lactam antibiotics and quantify the effects of selected enzymes on β-lactam susceptibility. This research significantly advances the understanding of the contributions of acquired β-lactamases to antibiotic resistance, including carbapenem resistance, in P. aeruginosa and by implication in other species. It has potential to expedite development of methods that use whole genome sequencing of infecting bacteria to inform antibiotic treatment, allowing more effective use of antibiotics, and facilitate the development of new antibiotics.

ColiSeq: a multiplex amplicon assay that provides strain level resolution of Escherichia coli directly from clinical specimens

Escherichia coli is a diverse pathogen, causing a range of disease in humans, from self-limiting diarrhea to urinary tract infections (UTIs). Uropathogenic E. coli (UPEC) is the most frequently observed uropathogen in UTIs, a common disease in high-income countries, incurring billions of dollars yearly in treatment costs. Although E. coli is easily grown and identified in the clinical laboratory, genotyping the pathogen is more complicated, yet critical for reducing the incidence of disease. These goals can be achieved through whole-genome sequencing of E. coli isolates, but this approach is relatively slow and typically requires culturing the pathogen in the laboratory. To genotype E. coli rapidly and inexpensively directly from clinical samples, including but not limited to urine, we developed and validated a multiplex amplicon sequencing assay, called ColiSeq. The assay consists of targets designed for E. coli species confirmation, high resolution genotyping, and mixture deconvolution. To demonstrate its utility, we screened the ColiSeq assay against 230 clinical urine samples collected from a hospital system in Flagstaff, Arizona, USA. A limit of detection analysis demonstrated the ability of ColiSeq to identify E. coli at a concentration of ~2 genomic equivalent (GEs)/mL and to generate high-resolution genotyping at a concentration of 1 × 105 GEs/mL. The results of this study suggest that ColiSeq could be a valuable method to understand the source of UPEC strains and guide infection mitigation efforts. As sequence-based diagnostics become accepted in the clinical laboratory, workflows such as ColiSeq will provide actionable information to improve patient outcomes.IMPORTANCEUrinary tract infections (UTIs), caused primarily by Escherichia coli, create an enormous health care burden in the United States and other high-income countries. The early detection of E. coli from clinical samples, including urine, is important to target therapy and prevent further patient complications. Additionally, understanding the source of E. coli exposure will help with future mitigation efforts. In this study, we developed, tested, and validated an amplicon sequencing assay focused on direct detection of E. coli from urine. The resulting sequence data were demonstrated to provide strain level resolution of the pathogen, not only confirming the presence of E. coli, which can focus treatment efforts, but also providing data needed for source attribution and contact tracing. This assay will generate inexpensive, rapid, and reproducible data that can be deployed by public health agencies to track, diagnose, and potentially mitigate future UTIs caused by E. coli.

Culture-independent detection of low-abundant Clostridioides difficile in environmental DNA via PCR

Clostridioides difficile represents a major burden to public health. As a well-known nosocomial pathogen whose occurrence is highly associated with antibiotic treatment, most examined C. difficile strains originated from clinical specimen and were isolated under selective conditions employing antibiotics. This suggests a significant bias among analyzed C. difficile strains, which impedes a holistic view on this pathogen. In order to support extensive isolation of C. difficile strains from environmental samples, we designed a detection PCR that targets the hpdBCA-operon and thereby identifies low abundances of C. difficile in environmental samples. This operon encodes the 4-hydroxyphenylacetate decarboxylase, which catalyzes the production of the antimicrobial compound para-cresol. Amplicon-based analyses of diverse environmental samples demonstrated that the designed PCR is highly specific for C. difficile and successfully detected C. difficile despite its absence in general 16S rRNA gene-based detection strategies. Further analyses revealed the potential of the hpdBCA detection PCR sequence for initial phylogenetic classification, which allows assessment of C. difficile diversity in environmental samples via amplicon sequencing. Our findings furthermore showed that C. difficile strains isolated under antibiotic treatment from environmental samples were originally dominated by other strains according to PCR amplicon results. This provided evidence for selective cultivation of under-represented but antibiotic-resistant isolates. Thereby, we revealed a substantial bias in C. difficile isolation and research.IMPORTANCEClostridioides difficile is a main cause of diarrheic infections after antibiotic treatment with serious morbidity and mortality worldwide. Research on this pathogen and its virulence has focused on bacterial isolation from clinical specimens under antibiotic treatment, which implies a substantial bias in isolated strains. Comprehensive studies, however, require an unbiased strain collection, which is accomplished by isolation of C. difficile from diverse environmental samples and avoidance of antibiotic-based enrichment strategies. Thus, isolation can significantly benefit from our C. difficile-specific detection PCR, which rapidly verifies C. difficile presence in environmental samples and further allows estimation of the C. difficile diversity by using next-generation sequencing.

Metagenomic profiling of viral and microbial communities from the pox lesions of lumpy skin disease virus and sheeppox virus-infected hosts

Introduction

It has been recognized that capripoxvirus infections have a strong cutaneous tropism with the manifestation of skin lesions in the form of nodules and scabs in the respective hosts, followed by necrosis and sloughing off. Considering that the skin microbiota is a complex community of commensal bacteria, fungi and viruses that are influenced by infections leading to pathological states, there is no evidence on how the skin microbiome is affected during capripoxvirus pathogenesis.

Methods

In this study, shotgun metagenomic sequencing was used to investigate the microbiome in pox lesions from hosts infected with lumpy skin disease virus and sheep pox virus.

Results

The analysis revealed a high degree of variability in bacterial community structures across affected skin samples, indicating the importance of specific commensal microorganisms colonizing individual hosts. The most common and abundant bacteria found in scab samples were Fusobacterium necrophorum, Streptococcus dysgalactiae, Helcococcus ovis and Trueperella pyogenes, irrespective of host. Bacterial reads belonging to the genera Moraxella, Mannheimia, Corynebacterium, Staphylococcus and Micrococcus were identified.

Discussion

This study is the first to investigate capripox virus-associated changes in the skin microbiome using whole-genome metagenomic profiling. The findings will provide a basis for further investigation into capripoxvirus pathogenesis. In addition, this study highlights the challenge of selecting an optimal bioinformatics approach for the analysis of metagenomic data in clinical and veterinary practice. For example, direct classification of reads using a kmer-based algorithm resulted in a significant number of systematic false positives, which may be attributed to the peculiarities of the algorithm and database selection. On the contrary, the process of de novo assembly requires a large number of target reads from the symbiotic microbial community. In this work, the obtained sequencing data were processed by three different approaches, including direct classification of reads based on k-mers, mapping of reads to a marker gene database, and de novo assembly and binning of metagenomic contigs. The advantages and disadvantages of these techniques and their practicality in veterinary settings are discussed in relation to the results obtained.

Recombination-mediated dissemination of Methicillin-resistant S. aureus clonal complex 1 in the Egyptian health care settings

BACKGROUND

Methicillin-resistant Staphylococcus aureus (MRSA) is a rapidly evolving pathogen that is frequently associated with outbreaks and sustained epidemics. This study investigated the population structure, resistome, virulome, and the correlation between antimicrobial resistance determinants with phenotypic resistance profiles of 36 representative hospital-acquired MRSA isolates recovered from hospital settings in Egypt.

RESULTS

The community-acquired MRSA lineage, clonal complex 1 (CC1) was the most frequently detected clone, followed by three other globally disseminated clones, CC121, CC8, and CC22. Most isolates carried SCCmec type V and more than half of isolates demonstrated multi-drug resistant phenotypes. Resistance to linezolid, a last resort antibiotic for treating multidrug resistant MRSA, was observed in 11.11% of the isolates belonging to different genetic backgrounds. Virulome analysis indicated that most isolates harboured a large pool of virulence factors and toxins. Genes encoding aureolysin, gamma hemolysins, and serine proteases were the most frequently detected virulence encoding genes. CC1 was observed to have a high pool of AMR resistance determinants including cfr, qacA, and qacB genes, which are involved in linezolid and quaternary ammonium compounds resistance, as well as high content of virulence-related genes, including both of the PVL toxin genes. Molecular clock analysis revealed that CC1 had the greatest frequency of recombination (compared to mutation) among the four major clones, supporting the role of horizontal gene transfer in modulating AMR and hypervirulence in this clone.

CONCLUSIONS

This pilot study provided evidence on the dissemination success of CA-MRSA clone CC1 among Egyptian hospitals. Co-detection of multiple AMR and virulence genes in this lineage pose a broad public health risk, with implications for successful treatment. The results of this study, together with other surveillance studies in Egypt, should be used to develop strategies for controlling MRSA infections in Egyptian health-care settings.

Genome characteristics related to the virulence of Streptococcus suis in Swedish pigs

The impact of S. suis on Swedish pig production has increased in recent years, and characterization of the strains present in the pig population is needed to aid in surveillance and prevention. Therefore, the aim of this study was to identify and characterize differences in the genomes between Swedish S. suis isolates associated with disease and isolates from healthy animals. Isolates categorized as being pathogenic (n = 100) or non-pathogenic (n = 117) were whole-genome sequenced, serotyped in silico, and sequence-typed using traditional MLST and core-genome MLST, and a genome-wide association study was performed to identify virulence-associated genes. In decreasing order, serotypes 2, 1, and 7 were the most common in the pathogenic group, and serotypes 15 and 12 were the most common in the non-pathogenic group. Among the commonly disease-associated sequence types, ST28 and ST25 were identified, whereas ST1 was scarcely found. The majority of isolates belonged to novel sequence types, revealing differences between Swedish isolates and those reported from other countries. The genomes of the pathogenic isolates were on average smaller and less heterogenic as compared to those of the non-pathogenic isolates. Although a majority of the previously published virulence-associated genes included in the study were found in the genomes of both pathogenic and non-pathogenic isolates, several new, significantly virulence-associated genes were identified.

Genomic analysis of Enterococcus faecium strain RAOG174 associated with acute chorioamnionitis carried antibiotic resistance gene: is it time for precise microbiological identification for appropriate antibiotic use?

BACKGROUND

Preterm labor syndrome is associated with high perinatal morbidity and mortality, and intra-amniotic infection is a cause of preterm labor. The standard identification of causative microorganisms is based on the use of biochemical phenotypes, together with broth dilution-based antibiotic susceptibility from organisms grown in culture. However, such methods could not provide an accurate epidemiological aspect and a genetic basis of antimicrobial resistance leading to an inappropriate antibiotic administration. Hybrid genome assembly is a combination of short- and long-read sequencing, which provides better genomic resolution and completeness for genotypic identification and characterization. Herein, we performed a hybrid whole genome assembly sequencing of a pathogen associated with acute histologic chorioamnionitis in women presenting with PPROM.

RESULTS

We identified Enterococcus faecium, namely E. faecium strain RAOG174, with several antibiotic resistance genes, including vancomycin and aminoglycoside. Virulence-associated genes and potential bacteriophage were also identified in this genome.

CONCLUSION

We report herein the first study demonstrating the use of hybrid genome assembly and genomic analysis to identify E. faecium ST17 as a pathogen associated with acute histologic chorioamnionitis. The analysis provided several antibiotic resistance-associated genes/mutations and mobile genetic elements. The occurrence of E. faecium ST17 raised the awareness of the colonization of clinically relevant E. faecium and the carrying of antibiotic resistance. This finding has brought the advantages of genomic approach in the identification of the bacterial species and antibiotic resistance gene for E. faecium for appropriate antibiotic use to improve maternal and neonatal care.

A local-scale One Health genomic surveillance of Clostridioides difficile demonstrates highly related strains from humans, canines, and the environment

Although infections caused by Clostridioides difficile have historically been attributed to hospital acquisition, growing evidence supports the role of community acquisition in C. difficile infection (CDI). Symptoms of CDI can range from mild, self-resolving diarrhoea to toxic megacolon, pseudomembranous colitis, and death. In this study, we sampled C. difficile from clinical, environmental, and canine reservoirs in Flagstaff, Arizona, USA, to understand the distribution and transmission of the pathogen in a One Health framework; Flagstaff is a medium-sized, geographically isolated city with a single hospital system, making it an ideal site to characterize genomic overlap between sequenced C. difficile isolates across reservoirs. An analysis of 562 genomes from Flagstaff isolates identified 65 sequence types (STs), with eight STs being found across all three reservoirs and another nine found across two reservoirs. A screen of toxin genes in the pathogenicity locus identified nine STs where all isolates lost the toxin genes needed for CDI manifestation (tcdB, tcdA), demonstrating the widespread distribution of non-toxigenic C. difficile (NTCD) isolates in all three reservoirs; 15 NTCD genomes were sequenced from symptomatic, clinical samples, including two from mixed infections that contained both tcdB+ and tcdB- isolates. A comparative single nucleotide polymorphism (SNP) analysis of clinically derived isolates identified 78 genomes falling within clusters separated by ≤2 SNPs, indicating that ~19 % of clinical isolates are associated with potential healthcare-associated transmission clusters; only symptomatic cases were sampled in this study, and we did not sample asymptomatic transmission. Using this same SNP threshold, we identified genomic overlap between canine and soil isolates, as well as putative transmission between environmental and human reservoirs. The core genome of isolates sequenced in this study plus a representative set of public C. difficile genomes (n=136), was 2690 coding region sequences, which constitutes ~70 % of an individual C. difficile genome; this number is significantly higher than has been published in some other studies, suggesting that genome data quality is important in understanding the minimal number of genes needed by C. difficile. This study demonstrates the close genomic overlap among isolates sampled across reservoirs, which was facilitated by maximizing the genomic search space used for comprehensive identification of potential transmission events. Understanding the distribution of toxigenic and non-toxigenic C. difficile across reservoirs has implications for surveillance sampling strategies, characterizing routes of infections, and implementing mitigation measures to limit human infection.

Metagenomic surveillance and comparative genomic analysis of Chlamydia psittaci in patients with pneumonia

Chlamydia psittaci, a strictly intracellular bacterium, is an underestimated etiologic agent leading to infections in a broad range of animals and mild illness or pneumonia in humans. In this study, the metagenomes of bronchoalveolar lavage fluids from the patients with pneumonia were sequenced and highly abundant C. psittaci was found. The target-enriched metagenomic reads were recruited to reconstruct draft genomes with more than 99% completeness. Two C. psittaci strains from novel sequence types were detected and these were closely related to the animal-borne isolates derived from the lineages of ST43 and ST28, indicating the zoonotic transmissions of C. psittaci would benefit its prevalence worldwide. Comparative genomic analysis combined with public isolate genomes revealed that the pan-genome of C. psittaci possessed a more stable gene repertoire than those of other extracellular bacteria, with ~90% of the genes per genome being conserved core genes. Furthermore, the evidence for significantly positive selection was identified in 20 virulence-associated gene products, particularly bacterial membrane-embedded proteins and type three secretion machines, which may play important roles in the pathogen-host interactions. This survey uncovered novel strains of C. psittaci causing pneumonia and the evolutionary analysis characterized prominent gene candidates involved in bacterial adaptation to immune pressures. The metagenomic approach is of significance to the surveillance of difficult-to-culture intracellular pathogens and the research into molecular epidemiology and evolutionary biology of C. psittaci.

Non-Toxigenic Clostridioides difficile Strain E4 (NTCD-E4) Prevents Establishment of Primary C. difficile Infection by Epidemic PCR Ribotype 027 in an In Vitro Human Gut Model

Clostridioides difficile infection (CDI) remains a significant healthcare burden. Non-toxigenic C. difficile (NTCD) strains have shown a benefit in preventing porcine enteritis and in human recurrent CDI. In this study, we evaluated the efficacy of metronidazole-resistant NTCD-E4 in preventing CDI facilitated by a range of antimicrobials in an in vitro human gut model. NTCD-E4 spores (at a dose of 107) were instilled 7 days before a clinical ribotype (RT) 027 (at the same dose) strain (210). In separate experiments, four different antimicrobials were used to perturb gut microbiotas; bacterial populations and cytotoxin production were determined using viable counting and Vero cell cytotoxicity, respectively. RT027 and NTCD-E4 proliferated in the in vitro model when inoculated singly, with RT027 demonstrating high-level cytotoxin (3-5-log10-relative units) production. In experiments where the gut model was pre-inoculated with NTCD-E4, RT027 was remained quiescent and failed to produce cytotoxins. NTCD-E4 showed mutations in hsmA and a gene homologous to CD196-1331, previously linked to medium-dependent metronidazole resistance, but lacked other metronidazole resistance determinants. This study showed that RT027 was unable to elicit simulated infection in the presence of NTCD-E4 following stimulation by four different antimicrobials. These data complement animal and clinical studies in suggesting NTCD offer prophylactic potential in the management of human CDI.

Strong pathogen competition in neonatal gut colonisation

Opportunistic bacterial pathogen species and their strains that colonise the human gut are generally understood to compete against both each other and the commensal species colonising this ecosystem. Currently we are lacking a population-wide quantification of strain-level colonisation dynamics and the relationship of colonisation potential to prevalence in disease, and how ecological factors might be modulating these. Here, using a combination of latest high-resolution metagenomics and strain-level genomic epidemiology methods we performed a characterisation of the competition and colonisation dynamics for a longitudinal cohort of neonatal gut microbiomes. We found strong inter- and intra-species competition dynamics in the gut colonisation process, but also a number of synergistic relationships among several species belonging to genus Klebsiella, which includes the prominent human pathogen Klebsiella pneumoniae. No evidence of preferential colonisation by hospital-adapted pathogen lineages in either vaginal or caesarean section birth groups was detected. Our analysis further enabled unbiased assessment of strain-level colonisation potential of extra-intestinal pathogenic Escherichia coli (ExPEC) in comparison with their propensity to cause bloodstream infections. Our study highlights the importance of systematic surveillance of bacterial gut pathogens, not only from disease but also from carriage state, to better inform therapies and preventive medicine in the future.

Metagenomic Insights into Pathogenic Characterization of ST410 Acinetobacter nosocomialis Prevalent in China

Acinetobacter nosocomialis is a prevalent opportunistic pathogen that causes hospital-acquired infections. The increasing threats from A. nosocomialis infections have led to attention from the scientific and medical communities. Metagenomic next-generation sequencing (mNGS) was performed for an exudate specimen collected from an ICU patient with wound infection, followed by sepsis, in Tongji Hospital. Three assembly strategies were employed to recover the genome of A. nosocomialis in the metagenomic sample. Together with publicly available genomes of A. nosocomialis, the features of population genetics and molecular epidemiology were deeply analyzed. A draft genome was reconstructed for the metagenomic strain WHM01, derived from the ST410 A. nosocomialis dominating the microbial community, thereby prompting its highly pathogenic risk, which is associated with infection and persistence. The structure of the bacterial pangenome was characterized, including the 1862 core and 11,815 accessory genes present in the 157 strains. The genetic diversity of the genes coding for the 128 virulence factors assigned to 14 functional categories was uncovered in this nosocomial pathogen, such as the lipooligosaccharide, capsule, type IV pilus, and outer membrane proteins. Our work revealed genomic properties of ST410 A. nosocomialis, which is prevalent in China, and further highlighted that metagenomic surveillance may be a prospective application for evaluating the pathogenic characteristics of the nosocomial opportunistic pathogens.

Intra-species diversity of Clostridium perfringens: A diverse genetic repertoire reveals its pathogenic potential

Clostridium perfringens is the causative agent of many enterotoxic diseases in humans and animals, and it is present in diverse environments (soil, food, sewage, and water). Multilocus Sequence Typing (MLST) and Whole Genome Sequencing (WGS) have provided a general approach about genetic diversity of C. perfringens; however, those studies are limited to specific locations and often include a reduced number of genomes. In this study, 372 C. perfringens genomes from multiple locations and sources were used to assess the genetic diversity and phylogenetic relatedness of this pathogen. In silico MLST was used for typing the isolates, and the resulting sequence types (ST) were assigned to clonal complexes (CC) based on allelic profiles that differ from its founder by up to double-locus variants. A pangenome analysis was conducted, and a core genome-based phylogenetic tree was created to define phylogenetic groups. Additionally, key virulence factors, toxinotypes, and antibiotic resistance genes were identified using ABRicate against Virulence Factor Database (VFDB), TOXiper, and Resfinder, respectively. The majority of the C. perfringens genomes found in publicly available databases were derived from food (n = 85) and bird (n = 85) isolates. A total of 195 STs, some of them shared between sources such as food and human, horses and dogs, and environment and birds, were grouped in 25 CC and distributed along five phylogenetic groups. Fifty-three percent of the genomes were allocated to toxinotype A, followed by F (32%) and G (7%). The most frequently found virulence factors based on > 70% coverage and 99.95% identity were plc (100%), nanH (99%), ccp (99%), and colA (98%), which encode an alpha-toxin, a sialidase, an alpha-clostripain, and a collagenase, respectively, while tetA (39.5%) and tetB (36.2%), which mediate tetracycline resistance determinants, were the most common antibiotic resistance genes detected. The analyses conducted here showed a better view of the presence of this pathogen across several host species. They also confirm that the genetic diversity of C. perfringens is based on a large number of virulence factors that vary among phylogroups, and antibiotic resistance markers, especially to tetracyclines, aminoglycosides, and macrolides. Those characteristics highlight the importance of C. perfringens as a one of the most common causes of foodborne illness.

Metagenomic Assessment of the Pathogenic Risk of Microorganisms in Sputum of Postoperative Patients With Pulmonary Infection

Respiratory infections are complicated biological processes associated with an unbalanced microbial community and a wide range of pathogens. To date, robust approaches are still required for distinguishing the pathogenic microorganisms from the colonizing ones in the clinical specimens with complex infection. In this study, we retrospectively analyzed the data of conventional culture testing and metagenomic next-generation sequencing (mNGS) of the sputum samples collected from 50 pulmonary infected patients after cardiac surgery from December 2020 and June 2021 in Ruijin Hospital. Taxonomic classification of the sputum metagenomes showed that the numbers of species belonging to bacteria, fungi, and viruses were 682, 58, and 21, respectively. The full spectrum of microorganisms present in the sputum microbiome covered all the species identified by culture, including 12 bacterial species and two fungal species. Based on species-level microbiome profiling, a reference catalog of microbial abundance detection limits was constructed to assess the pathogenic risks of individual microorganisms in the specimens. The proposed screening procedure detected 64 bacterial pathogens, 10 fungal pathogens, and three viruses. In particular, certain opportunistic pathogenic strains can be distinguished from the colonizing ones in the individual specimens. Strain-level identification and phylogenetic analysis were further performed to decipher molecular epidemiological characteristics of four opportunistic etiologic agents, including Klebsiella pneumoniae, Corynebacterium striatum, Staphylococcus aureus, and Candida albicans. Our findings provide a novel metagenomic insight into precision diagnosis for clinically relevant microbes, especially for opportunistic pathogens in the clinical setting.

Antimicrobial resistance and population genomics of multidrug-resistant Escherichia coli in pig farms in mainland China

The expanding use of antimicrobials in livestock is an important contributor to the worldwide rapid increase in antimicrobial resistance (AMR). However, large-scale studies on AMR in livestock remain scarce. Here, we report findings from surveillance of E. coli AMR in pig farms in China in 2018-2019. We isolated E. coli in 1,871 samples from pigs and their breeding environments, and found AMR in E. coli in all provinces in mainland China. We detected multidrug-resistance in 91% isolates and found resistance to last-resort drugs including colistin, carbapenems and tigecycline. We also identified a heterogeneous group of O-serogroups and sequence types among the multidrug-resistant isolates. These isolates harbored multiple resistance genes, virulence factor-encoding genes, and putative plasmids. Our data will help to understand the current AMR profiles of pigs and provide a reference for AMR control policy formulation for livestock in China.