Unveiling the microbiome during post-partum uterine infection: a deep shotgun sequencing approach to characterize the dairy cow uterine microbiome

BACKGROUND

The goal of this study was to assess the microbial ecology and diversity present in the uterus of post-partum dairy cows with and without metritis from 24 commercial California dairy farms using shotgun metagenomics. A set subset of 95 intrauterine swab samples, taken from a larger selection of 307 individual cow samples previously collected, were examined for α and β diversity and differential abundance associated with metritis. Cows within 21 days post-partum were categorized into one of three clinical groups during sample collection: control (CT, n = 32), defined as cows with either no vaginal discharge or a clear, non-purulent mucus vaginal discharge; metritis (MET, n = 33), defined as a cow with watery, red or brown colored, and fetid vaginal discharge; and purulent discharge cows (PUS, n = 31), defined as a non-fetid purulent or mucopurulent vaginal discharge.

RESULTS

All three clinical groups (CT, MET, and PUS) were highly diverse, with the top 12 most abundant genera accounting for 10.3%, 8.8%, and 10.1% of mean relative abundance, respectively. The α diversity indices revealed a lower diversity from samples collected from MET and PUS when compared to CT cows. PERMANOVA statistical testing revealed a significant difference (P adjusted < 0.01) in the diversity of genera between CT and MET samples (R2 = 0.112, P = 0.003) and a non-significant difference between MET and PUS samples (R2 = 0.036, P = 0.046). ANCOM-BC analysis revealed that from the top 12 most abundant genera, seven genera were increased in the natural log fold change (LFC) of abundance in MET when compared to CT samples: Bacteroides, Clostridium, Fusobacterium, Phocaeicola, Porphyromonas, Prevotella, and Streptococcus. Two genera, Dietzia and Microbacterium, were decreased in natural LFC of abundance when comparing MET (regardless of treatment) and CT, while no changes in natural LFC of abundance were observed for Escherichia, Histophilus, and Trueperella.

CONCLUSIONS

The results presented here, are the current deepest shotgun metagenomic analyses conducted on the bovine uterine microbiome to date (mean of 256,425 genus-level reads per sample). Our findings support that uterine samples from cows without metritis (CT) had increased α-diversity but decreased β-diversity when compared to metritis or PUS cows, characteristic of dysbiosis. In summary, our findings highlight that MET cows have an increased abundance of Bacteroides, Porphyromonas, and Fusobacterium when compared to CT and PUS, and support the need for further studies to better understand their potential causal role in metritis pathogenesis.

Antimicrobial resistance of Vibrio spp. from the coastal California system: discordance between genotypic and phenotypic patterns

Antimicrobial resistance in Vibrio species poses risks to both human and marine mammal health. Whole genome sequencing of Vibrio spp. can be utilized to screen for antimicrobial resistance genes and allelic variants to provide mechanistic insights in ways that PCR screening and phenotypic interpretation cannot. Our goals were to (i) characterize antimicrobial resistance patterns of Vibrio spp. pathogens isolated from southern sea otters (Enhydra lutris nereis), northern sea otters (Enhydra lutris kenyoni), and environmental samples from the central California coast using whole genome sequencing, and (ii) compare the presence of antimicrobial resistance genes with phenotypic interpretation from antibiotic susceptibility testing. Unexpectedly, genomic classification identified an understudied species, Vibrio diabolicus, in sea otter and environmental isolates that were previously identified as Vibrio alginolyticus. A total of 489 Vibrio spp. isolates were sequenced, and frequently detected antimicrobial resistance genes included multidrug efflux pumps and genes associated with resistance to ß-lactams and tetracyclines. Genes associated with resistance to fluoroquinolones, aminoglycosides, chloramphenicol, and sulfonamides were uncommon. Sea otter isolates were phenotypically susceptible to tetracycline despite carrying genes tet34 and tet35. Both between- and within-species variations in ampicillin resistance were observed despite the ubiquitous presence of blaCARB genes in V. alginolyticus, V. diabolicus, and Vibrio parahaemolyticus. Discordance between phenotypic and genotypic ampicillin resistance was especially noted for V. parahaemolyticus and was partially attributed to the allelic variation of the blaCARB genes. Tetracyclines and fluoroquinolones, but not ß-lactams, are likely to be effective treatments for vibriosis in sea otters.

IMPORTANCE

Vibriosis (infection with non-cholera Vibrio spp.) is the most common seafood-borne illness globally, with major impacts on public health, food security, and wildlife health. Potential treatments of antimicrobial-resistant Vibrio spp. in humans, aquaculture, and marine wildlife rehabilitation are complicated by current diagnostic challenges regarding bacterial species identification and interpretation of antimicrobial resistance patterns. Unexpected detection of previously misidentified Vibrio diabolicus in sea otters suggests that a broader taxonomic group of Vibrio infect sea otters than previously described. We also determined that the presence of ß-lactamase genes alone in sea otter isolates does not necessarily correlate with an ampicillin-resistant phenotype, likely due to deleterious amino acid substitutions in certain blaCARB alleles. Continued monitoring of Vibrio spp. phenotypes and genotypes in sea otters is warranted to observe biologically relevant changes in antimicrobial resistance.

Antimicrobial Resistance Genes in Respiratory Bacteria from Weaned Dairy Heifers

Bovine respiratory disease (BRD) is the leading cause of mortality and antimicrobial drug (AMD) use in weaned dairy heifers. Limited information is available regarding antimicrobial resistance (AMR) in respiratory bacteria in this population. This study determined AMR gene presence in 326 respiratory isolates (Pasteurella multocida, Mannheimia haemolytica, and Histophilus somni) from weaned dairy heifers using whole genome sequencing. Concordance between AMR genotype and phenotype was determined. Twenty-six AMR genes for 8 broad classes of AMD were identified. The most prevalent, medically important AMD classes used in calf rearing, to which these genes predict AMR among study isolates were tetracycline (95%), aminoglycoside (94%), sulfonamide (94%), beta-lactam (77%), phenicol (50%), and macrolide (44%). The co-occurrence of AMR genes within an isolate was common; the largest cluster of gene co-occurrence encodes AMR to phenicol, macrolide, elfamycin, β-lactam (cephalosporin, penam cephamycin), aminoglycoside, tetracycline, and sulfonamide class AMD. Concordance between genotype and phenotype varied (Matthew's Correlation Coefficient ranged from -0.57 to 1) by bacterial species, gene, and AMD tested, and was particularly poor for fluoroquinolones (no AMR genes detected) and ceftiofur (no phenotypic AMR classified while AMR genes present). These findings suggest a high genetic potential for AMR in weaned dairy heifers; preventing BRD and decreasing AMD reliance may be important in this population.

WGS of intrauterine E. coli from cows with early postpartum uterine infection reveals a non-uterine specific genotype and virulence factors

Escherichia coli has been attributed to playing a major role in a cascade of events that affect the prevalence and severity of uterine disease in cattle. The objectives of this project were to (i) define the association between the prevalence of specific antimicrobial resistance and virulence factor genes in E. coli with the clinical status related to uterine infection, (ii) identify the genetic relationship between E. coli isolates from cows with diarrhea, with mastitis, and with and without metritis, and (iii) determine the association between the phenotypic and genotypic antimicrobial resistance identified on the E. coli isolated from postpartum cattle. Bacterial isolates (n = 148) were obtained from a larger cross-sectional study. Cows were categorized into one of three clinical groups before enrollment: metritis, cows with purulent discharge, and control cows. For genomic comparison, public genomes (n = 130) from cows with diarrhea, mastitis, and metritis were included in a genome-wide association study, to evaluate differences between the drug classes or the virulence factor category among clinical groups. A distinct E. coli genotype associated with metritis could not be identified. Instead, a high genetic diversity among the isolates from uterine sources was present. A virulence factor previously associated with metritis (fimH) using PCR was not associated with metritis. There was moderate accuracy for whole-genome sequencing to predict phenotypic resistance, which varied depending on the antimicrobial tested. Findings from this study contradict the traditional pathotype classification and the unique intrauterine E. coli genotype associated with metritis in dairy cows.IMPORTANCEMetritis is a common infectious disease in dairy cattle and the second most common reason for treating a cow with antimicrobials. The pathophysiology of the disease is complex and is not completely understood. Specific endometrial pathogenic Escherichia coli have been reported to be adapted to the endometrium and sometimes lead to uterine disease. Unfortunately, the specific genomic details of the endometrial-adapted isolates have not been investigated using enough genomes to represent the genomic diversity of this organism to identify specific virulence genes that are consistently associated with disease development and severity. Results from this study provide key microbial ecological advances by elucidating and challenging accepted concepts for the role of Intrauterine E. coli in metritis in dairy cattle, especially contradicting the existence of a unique intrauterine E. coli genotype associated with metritis in dairy cows, which was not found in our study.

Whole genome sequence analysis reveals high genomic diversity and potential host-driven adaptations among multidrug-resistant Escherichia coli from pre-weaned dairy calves

Food-producing animals such as dairy cattle are potential reservoirs of antimicrobial resistance (AMR), with multidrug-resistant (MDR) organisms such as Escherichia coli observed in higher frequency in young calves compared to older cattle. In this study, we characterized the genomes of enteric MDR E. coli from pre-weaned dairy calves with and without diarrhea and evaluated the influence of host-level factors on genomic composition. Whole genome sequence comparative analysis of E. coli (n = 43) revealed substantial genomic diversity that primarily clustered by sequence type and was minimally driven by calf diarrheal disease status (healthy, diarrheic, or recovered), antimicrobial exposure, and dietary zinc supplementation. Diverse AMR genes (ARGs)-including extended-spectrum beta-lactamase genes and quinolone resistance determinants-were identified (n = 40), with unique sets of ARGs co-occurring in gene clusters with large AMR plasmids IncA/C2 and IncFIB(AP001918). Zinc supplementation was not significantly associated with the selection of individual ARGs in E. coli, however analysis of ARG and metal resistance gene pairs identified positive associations between certain aminoglycoside, beta-lactam, sulfonamide, and trimethoprim ARGs with acid, tellurium and mercury resistance genes. Although E. coli in this study lacked the typical virulence factors of diarrheagenic strains, virulence genes overlapping with those in major pathotypes were identified. Among the 103 virulence genes detected, the highest abundance and diversity of genes corresponded to iron acquisition (siderophores and heme uptake). Our findings indicate that the host-level factors evaluated in this study were not key drivers of genomic variability, but that certain accessory genes in enteric MDR E. coli may be enriched. Collectively, this work provides insight into the genomic diversity and host-microbe interface of MDR E. coli from pre-weaned dairy calves.