Comparative metagenomics of microbial communities and resistome in southern farming systems: implications for antimicrobial stewardship and public health

Agricultural practices significantly influence microbial diversity and the distribution of virulence and antimicrobial resistance (AMR) genes, with implications for ecosystem health and food safety. This study used metagenomic sequencing to analyze 60 samples (30 per state) including water, soil, and manure (10 each) from Alabama (a mix of cattle and poultry sources) and Tennessee (primarily from cattle). The results highlighted a rich microbial diversity, predominantly comprising Bacteria (67%) and Viruses (33%), with a total of over 1,950 microbial species identified. The dominant bacterial phyla were Proteobacteria, Cyanobacteria, Actinobacteria, Firmicutes, and Bacteroidetes, with the viral communities primarily represented by Phixviricota and Uroviricota. Distinct state-specific microbial profiles were evident, with Alabama demonstrating a higher prevalence of viral populations and unique bacterial phyla compared to Tennessee. The influence of environmental and agricultural practices was reflected in the microbial compositions: soil samples were notably rich in Actinobacteria, water samples were dominated by Proteobacteria and Cyanobacteria, and manure samples from Alabama showed a predominance of Actinobacteria. Further analyses, including diversity assessment and enterotype clustering, revealed complex microbial structures. Tennessee showed higher microbial diversity and phylogenetic complexity across most sample types compared to Alabama, with poultry-related samples displaying distinct diversity trends. Principal Coordinate Analysis (PCoA) highlighted notable state-specific variations, particularly in manure samples. Differential abundance analysis demonstrated elevated levels of Deinococcus and Ligilactobacillus in Alabama, indicating regional effects on microbial distributions. The virulome analysis revealed a significant presence of virulence genes in samples from Alabama. The community resistome was extensive, encompassing 109 AMR genes across 18 antibiotic classes, with manure samples displaying considerable diversity. Ecological analysis of the interactions between AMR gene subtypes and microbial taxa revealed a sophisticated network, often facilitated by bacteriophages. These findings underscore the critical role of agricultural practices in shaping microbial diversity and resistance patterns, highlighting the need for targeted AMR mitigation strategies in agricultural ecosystems to protect both public health and environmental integrity.

Multidrug-Resistant Extended-Spectrum β-Lactamase-Producing Escherichia coli Pathotypes in North Eastern Region of India: Backyard Small Ruminants-Human-Water Interface

A total of 648 diarrheagenic Escherichia coli (DEC) were isolated from calves (n = 219), lambs (n = 87), kids (n = 103), human (n = 193), and water (n = 46) samples. The presence of enteropathogenic E. coli (EPEC), enterotoxigenic E. coli (ETEC), and shigatoxigenic E. coli (STEC) was confirmed by PCR-based detection of the Shiga toxin, intimin, hemolysin, and enterotoxin genes. All the isolates were tested for antimicrobial resistance (AMR) by disc diffusion assay. Extended-spectrum β-lactamase (ESBL), carbapenemase, and metallo-beta-lactamase production were determined by double-disk synergy test, modified Hodge test, and combined disk test assays. AMR genes (bla_TEM, bla_SHV, bla_CTX-M, bla_CMY-2, bla_NDM, bla_KPC, bla_VIM, and bla_IMP) were detected by PCR using specific primers. Majority of the isolates from human and water exhibited resistance (>80%) against amoxicillin, ampicillin, aztreonam, cefotaxime, cefixime, gentamicin, ceftazidime, and cefalexin, and against imipenem (70.98%), doripenem (70.47%), and ertapenem (60.62%). Bovine isolates were sensitive to carbapenems. Many isolates (5.75-24.35%) from human, water, calves, kids, and lambs were multidrug resistant (MDR), with resistance against three or more classes of antimicrobials. A total of 170/648 (26.23%) isolates were classified as STEC (9.88%), EPEC (4.32%), and ETEC (12.04%). The AMR genes, including bla_TEM, bla_CMY2, bla_CTX-M, and bla_SHV were detected in the E. coli from all sources. but bla_NDM and bla_KPC were detected only in the isolates from human and water. Three STEC isolates from human origin possessed multiple ESBLs, carbapenemase and metallo-beta-lactamase genes reported for the first time. ESBLs producing EPEC and ETEC in lambs and kids are also reported under this study. Presence of MDR-DEC in domestic animals and common potable water poses public health concern in this region.

Investigation of Virulence Genes Detected in Antimicrobial-Resistance Pathogens Isolates for Five Countries across the World

A large portion of annual deaths worldwide are due to infections caused by disease-causing pathogens. These pathogens contain virulence genes, which encode mechanisms that facilitate infection and microbial survival in hosts. More recently, antimicrobial resistance (AMR) genes, also found in these pathogens, have become an increasingly large issue. While the National Center for Biotechnology Information (NCBI) Pathogen Detection Isolates Browser (NPDIB) database has been compiling genes involved in microbial virulence and antimicrobial resistance through isolate samples, few studies have identified the genes primarily responsible for virulence and compared them to those responsible for AMR. This study performed the first multivariate statistical analysis of the multidimensional NPDIB data to identify the major virulence genes from historical pathogen isolates for Australia, China, South Africa, UK, and US—the largely populated countries from five of the six major continents. The important virulence genes were then compared with the AMR genes to study whether there is correlation between their occurrences. Among the significant genes and pathogens associated with virulence, it was found that the genes fdeC, iha, iss, iutA, lpfA, sslE, ybtP, and ybtQ are shared amongst all five countries. The pathogens E. coli and Shigella, Salmonella enterica, and Klebsiella pneumoniae mostly contained these genes and were common among four of the five studied countries. Additionally, the trend of virulence was investigated by plotting historical occurrences of gene and pathogen frequency in the annual samples. These plots showed that the trends of E. coli and Shigella and Salmonella enterica were similar to the trends of certain virulence genes, confirming the two pathogens do indeed carry important virulence genes. While the virulence genes in the five countries are not significantly different, the US and the UK share the largest amount of important virulence genes. The plots from principal component analysis and hierarchical clustering show that the important virulence and AMR genes were not significantly correlated, with only few genes from both types of genes clustered into the same groups.

Coinfection with Human Norovirus and Escherichia coli O25:H4 Harboring Two Chromosomal blaCTX-M-14 Genes in a Foodborne Norovirus Outbreak in Shizuoka Prefecture, Japan

ABSTRACT

Hospital-acquired infections caused by extended-spectrum β-lactamase (ESBL)-producing Escherichia coli are a global problem. Healthy people can carry ESBL-producing E. coli in the intestines; thus, E. coli from healthy people can potentially cause hospital-acquired infections. Therefore, the transmission routes of ESBL-producing E. coli from healthy persons should be determined. A foodborne outbreak of human norovirus (HuNoV) GII occurred at a restaurant in Shizuoka, Japan, in 2018. E. coli O25:H4 was isolated from some of the HuNoV-infected customers. Pulsed-field gel electrophoresis showed that these E. coli O25:H4 strains originated from one clone. Because the only epidemiological link among the customers was eating food from this restaurant, the customers were concurrently infected with E. coli O25:H4 and HuNoV GII via the restaurant food. Whole genome analysis revealed that the E. coli O25:H4 strains possessed genes for regulating intracellular iron and expressing the flagellum and flagella. Extraintestinal pathogenic E. coli often express these genes on the chromosome. Additionally, the E. coli O25:H4 strains had plasmids harboring nine antimicrobial resistance genes. These strains harbored ESBL-encoding blaCTX-M-14 genes on two loci of the chromosome and had higher ESBL activity. Multilocus sequence typing and fimH subtyping revealed that the E. coli O25:H4 strains from the outbreak belonged to the subclonal group, ST131-fimH30R, which has been driving ESBL epidemics in Japan. Because the E. coli O25:H4 strains isolated in the outbreak belonged to a subclonal group spreading in Japan, foods contaminated with ESBL-producing E. coli might contribute to spreading these strains among healthy persons. The isolated E. coli O25:H4 strains produced ESBL and contained plasmids with multiple antimicrobial resistance genes, which may make it difficult to select antimicrobials for treating extraintestinal infections caused by these strains.

HIGHLIGHTS

Genomic Variability among Field Isolates and Laboratory-Adapted Strains of Leptospira borgpetersenii Serovar Hardjo

Leptospira borgpetersenii serovar Hardjo colonizes cattle kidneys and may occasionally infect humans and other mammals. Strains belonging to two clonal subtypes (types A and B) with marked differences in their pathogenicity in the hamster experimental model have been described for this serovar. Such differences have been attributed to point mutations in individual genes, although those genes have not yet been characterized. In order to better understand genetic variability among L. borgpetersenii serovar Hardjo isolates, we sequenced and compared the genomes of two laboratory-adapted strains and three abattoir-derived field isolates of L. borgpetersenii serovar Hardjo. Relatively low genetic variability was observed within isolates of the same subtype, with most of the mutations of moderate or high impact found in the laboratory-adapted isolates. In contrast, several differences regarding gene content and genetic variants were observed between the two subtypes. Putative type-specific genes appear to encode proteins associated with functions that are critical for infection. Some of these genes seem to be involved in transcriptional regulation, possibly leading to a distinct regulatory pattern in each type. These results show that changes in regulatory mechanisms, previously suggested to be critical during Leptospira speciation, may occur in L. borgpetersenii. In addition, the bioinformatics methodology used in this study for variant calling can be useful to other groups working with nonmodel prokaryotic organisms such as Leptospira species.

Genotypic Features of Clinical and Bovine Escherichia coli O157 Strains Isolated in Countries with Different Associated-Disease Incidences

There is great geographical variation in the frequency of Escherichia coli O157 infections that correlates with important differences in the bovine reservoir of each country. Our group carried out a broad molecular characterization of human and bovine E. coli O157 strains circulating in Argentina using different methodologies. Our data allows us to conclude that in Argentina, a high homogeneity is observed in both cattle and human strains, with almost exclusive circulation of strains belonging to the hypervirulent clade 8 described by Manning. The aim of this review was to compare the genetic background of E. coli O157 strains isolated in countries that have conducted similar studies, to try to correlate specific O157 genotypes with the incidence and severity of E. coli O157 associated diseases. The characteristics of the strains that cause disease in humans reflect the predominant genotypes in cattle in each of the countries analyzed. The main features clearly linked to high incidence or severity of E. coli O157 infections are lineage-specific polymorphism assay-6 lineage I/II, clade 8 strains and probably, clade 6 strains, the stx_2a/stx_2c genotype, the presence of q933 and q21 simultaneously, and putative virulence factor EC_3286. In countries with an absence of these features in O157 strains, the overall incidence of O157 disease is low. Argentina, where these characteristics are detected in most strains, shows the highest incidence of hemolytic uremic syndrome (HUS) worldwide.

Predicting Essential Metabolic Genome Content of Niche-Specific Enterobacterial Human Pathogens during Simulation of Host Environments

Microorganisms have evolved to occupy certain environmental niches, and the metabolic genes essential for growth in these locations are retained in the genomes. Many microorganisms inhabit niches located in the human body, sometimes causing disease, and may retain genes essential for growth in locations such as the bloodstream and urinary tract, or growth during intracellular invasion of the hosts' macrophage cells. Strains of Escherichia coli (E. coli) and Salmonella spp. are thought to have evolved over 100 million years from a common ancestor, and now cause disease in specific niches within humans. Here we have used a genome scale metabolic model representing the pangenome of E. coli which contains all metabolic reactions encoded by genes from 16 E. coli genomes, and have simulated environmental conditions found in the human bloodstream, urinary tract, and macrophage to determine essential metabolic genes needed for growth in each location. We compared the predicted essential genes for three E. coli strains and one Salmonella strain that cause disease in each host environment, and determined that essential gene retention could be accurately predicted using this approach. This project demonstrated that simulating human body environments such as the bloodstream can successfully lead to accurate computational predictions of essential/important genes.

Colonization of Beef Cattle by Shiga Toxin-Producing Escherichia coli during the First Year of Life: A Cohort Study

Each year Shiga toxin-producing Escherichia coli (STEC) are responsible for 2.8 million acute illnesses around the world and > 250,000 cases in the US. Lowering the prevalence of this pathogen in animal reservoirs has the potential to reduce STEC outbreaks in humans by controlling its entrance into the food chain. However, factors that modulate the colonization and persistence of STEC in beef cattle remain largely unidentified. This study evaluated if animal physiological factors such as age, breed, sex, and weight gain influenced the shedding of STEC in beef cattle. A cohort of beef calves (n = 260) from a multi-breed beef calf population was sampled every three months after birth to measure prevalence and concentration of STEC during the first year of life. Metagenomic analysis was also used to understand the association between the STEC colonization and the composition of gut microflora. This study identified that beef calves were more likely to shed STEC during the first 6 months and that STEC shedding decreased as the animal matured. Animal breed group, sex of the calf, and average weight gain were not significantly associated with STEC colonization. The metagenomic analysis revealed for the first time that STEC colonization was correlated with a lower diversity of gut microflora, which increases as the cattle matured. Given these findings, intervention strategies that segregate younger animals, more likely to be colonized by STEC from older animals that are ready to be harvested, could be investigated as a method to reduce zoonotic transmission of STEC from cattle to humans.

Genomic Comparative Study of Bovine Mastitis Escherichia coli

Escherichia coli, one of the main causative agents of bovine mastitis, is responsible for significant losses on dairy farms. In order to better understand the pathogenicity of E. coli mastitis, an accurate characterization of E. coli strains isolated from mastitis cases is required. By using phylogenetic analyses and whole genome comparison of 5 currently available mastitis E. coli genome sequences, we searched for genotypic traits specific for mastitis isolates. Our data confirm that there is a bias in the distribution of mastitis isolates in the different phylogenetic groups of the E. coli species, with the majority of strains belonging to phylogenetic groups A and B1. An interesting feature is that clustering of strains based on their accessory genome is very similar to that obtained using the core genome. This finding illustrates the fact that phenotypic properties of strains from different phylogroups are likely to be different. As a consequence, it is possible that different strategies could be used by mastitis isolates of different phylogroups to trigger mastitis. Our results indicate that mastitis E. coli isolates analyzed in this study carry very few of the virulence genes described in other pathogenic E. coli strains. A more detailed analysis of the presence/absence of genes involved in LPS synthesis, iron acquisition and type 6 secretion systems did not uncover specific properties of mastitis isolates. Altogether, these results indicate that mastitis E. coli isolates are rather characterized by a lack of bona fide currently described virulence genes.