Blood as a route of transmission of uterine pathogens from the gut to the uterus in cows
© The Author(s). 2017
Received: 20 June 2017
Accepted: 20 August 2017
Published: 25 August 2017
Metritis is an inflammatory disease of the uterus caused by bacterial infection, particularly Bacteroides, Porphyromonas, and Fusobacterium. Bacteria from the environment, feces, or vagina are believed to be the only sources of uterine contamination. Blood seeps into the uterus after calving; therefore, we hypothesized that blood could also be a seeding source of uterine bacteria. Herein, we compared bacterial communities from blood, feces, and uterine samples from the same cows at 0 and 2 days postpartum using deep sequencing and qPCR. The vaginal microbiome 7 days before calving was also compared.
There was a unique structure of bacterial communities by sample type. Principal coordinate analysis revealed two distinct clusters for blood and feces, whereas vaginal and uterine bacterial communities were more scattered, indicating greater variability. Cluster analysis indicated that uterine bacterial communities were more similar to fecal bacterial communities than vaginal and blood bacterial communities. Nonetheless, there were core genera shared by all blood, feces, vaginal, and uterine samples. Major uterine pathogens such as Bacteroides, Porphyromonas, and Fusobacterium were part of the core genera in blood, feces, and vagina. Other uterine pathogens such as Prevotella and Helcococcus were not part of the core genera in vaginal samples. In addition, uterine pathogens showed a strong and significant interaction with each other in the network of blood microbiota, but not in feces or vagina. These microbial interactions in blood may be an important component of disease etiology. The copy number of total bacteria in blood and uterus was correlated; the same did not occur in other sites. Bacteroides heparinolyticus was more abundant in the uterus on day 0, and both B. heparinolyticus and Fusobacterium necrophorum were more abundant in the uterus than in the blood and feces on day 2. This indicates that B. heparinolyticus has a tropism for the uterus, whereas both pathogens thrive in the uterine environment early postpartum.
Blood harbored a unique microbiome that contained the main uterine pathogens such as Bacteroides, Porphyromonas, and Fusobacterium. The presence of these pathogens in blood shortly after calving shows the feasibility of hematogenous spread of uterine pathogens in cows.
KeywordsBlood microbiota Fecal microbiota Uterine microbiota Dairy cows Droplet digital PCR Bacteroides heparinolyticus Fusobacterium necrophorum Network analysis Metagenomics Uterine disease
According to the World Health Organization, maternal sepsis from infection of the uterus (metritis) postpartum is still prevalent (~ 5%) in developing countries and accounted for 77,000 maternal deaths worldwide in 2000 . The incidence of metritis in dairy cows is even greater (~ 20%) [2–4]; therefore, cows represent a good model for the study of uterine infections in large mammals. In our previous study, the uterine microbiota at calving was discriminated between cows that remained healthy and cows that later developed metritis, showing that the uterus begins to establish a bacterial community towards either health or disease shortly after calving . Management strategies to control the source of bacterial contamination or to control the proliferation of pathogenic bacteria in the uterus such as vaccination  or the use of probiotics  could help prevent metritis in dairy cows.
It is widely believed that uterine bacteria ascend from the vagina or through the vagina from the environment or feces, when the cervix, which serves as an anatomical and immunological barrier, opens during parturition [8, 9]. Pathogens associated with the development of metritis such as Bacteroides, Fusobacterium, and Porphyromonas [5, 10, 11] are part of the normal flora of the rumen in cows  and are shed in feces; therefore, ascending uterine contamination from the environment could contribute to the development of metritis. The vaginal microbiota of beef and dairy cows have also been shown to harbor the main uterine pathogens such as Bacteroides, Fusobacterium, and Porphyromonas [13, 14]; therefore, ascending uterine contamination from the vagina is also possible. Nonetheless, one specific uterine pathogen, Fusobacterium necrophorum (F. necrophorum), is known to gain access to the circulation, probably during episodes of rumen acidosis, and cause liver abscesses in cows . Interestingly, F. necrophorum is usually co-cultured with Trueperella (Arcanobacterium) pyogenes , another important uterine pathogen . Furthermore, Helcococcus ovis, an emerging uterine pathogen [5, 11], was first reported to cause valvular endocarditis in cattle . Therefore, hematogenous transmission must be considered as a possible route of dissemination of uterine pathogens in addition to ascending contamination from or through the vagina. Indeed, blood is a normal component of lochia; therefore, maternal blood is naturally transferred to the lumen of the uterus after birth. In cows, this mostly occurs because of degenerative vascular changes characterized by pyknosis followed by karyorrhexis of the endothelial cells and the cells of the tunica media of small blood vessels, changes that are observed within 24 h after calving . Indeed, the endometrium becomes highly edematous during the first 24 h after calving , which is probably the result of the aforementioned vascular changes, hence leading to leakage of blood components into the uterine lumen. Although bacteria could be free-floating in blood, recent studies in mice and cows showed that bacteria could be transported to extraintestinal sites by mononuclear cells [19, 20]. Interestingly, translocation to extraintestinal sites was more common in mice in late gestation and shortly after parturition than in mice that were not pregnant or were in early or mid-lactation . In cows as well as in other species, there is massive migration of leukocytes to the uterus with impending parturition and into the uterus after parturition [21, 22]; therefore, free-floating bacteria or bacteria engulfed by monocytes/macrophages could be readily transferred to the uterine lumen after calving.
The blood microbiome has not been investigated in peripartum dairy cows. Here, we hypothesized that cow’s blood would have a microbiome that would contain the main uterine pathogens such as Fusobacterium, Bacteroides, and Porphyromonas, therefore making the hematogenous route a feasible transmission route from the gut to the uterus. Given that microbiota in feces represents the distal portion of the gut microbiota , we sought to characterize bacterial communities in the blood, feces, and uterine samples to examine how these communities are related to each other. We collected samples from the same 12 individuals within 60 min after calving to minimize exogenous bacterial contamination and at 2 days after calving to examine the change in bacterial community. Blood was collected aseptically from the jugular vein. Because the vaginal microbiome also harbors uterine pathogens, we included the vaginal microbiome from a previous study  in our analysis for comparison. This study provides insight into the origin of uterine bacteria as well as the potential role of the gut and blood microbiota in uterine disease.
Characteristics of the study samples
Blood, fecal, and uterine samples were collected from the same 12 Holstein dairy cows on day 0 (the day of calving) to minimize exogenous bacterial contamination and on day 2 (2 days postpartum) to examine the change in bacterial community. The V4 region of the bacterial 16S ribosomal RNA (rRNA) gene was PCR-amplified from all 72 samples, and sequencing was performed on the Illumina MiSeq platform. Samples that failed quality control were excluded for taxonomic classification; therefore, 19 blood samples (9 on day 0 and 10 on day 2), 22 fecal samples (10 on day 0 and 12 on day 2), and 20 uterine samples (10 on day 0 and 10 on day 2) were analyzed. 16S rRNA sequencing resulted in 6,818,977 reads, with an average of 85,651 ± 5888 reads (standard error of the mean) from blood samples, 156,566 ± 62,878 reads (standard error of the mean) from fecal samples, and 87,358 ± 4956 reads (standard error of the mean) from uterine samples. Rarefaction curves of 61 samples at the minimum cutoff of 97% sequence identity nearly reached a plateau, which indicates that sampling depth is sufficient to characterize bacterial communities (Additional file 1: Figure S1).
For comparison, vaginal data were obtained from a previously published study , which were generated using the same sequencing technique and quality control as other samples. In our analysis, we used data from vaginal samples collected on day − 7 (7 days prepartum) from 105 Holstein dairy cows that were different from the cows used for collection of samples from the uterus, feces, and blood. Although samples collected from different individuals have been used to compare the microbiome from different body sites such as the oral cavity, the gut, and the human placenta , these samples do not allow for direct comparisons between body sites such as comparison of means or evaluation of correlations between individual taxon. Although vaginal samples from day 0 were also available, those would be contaminated with uterine discharge, therefore not being able to differentiate between the uterine and vaginal microbiome.
Dissimilarity of blood, fecal, and uterine microbiota
Similarity of bacterial community composition in blood, feces, and uterine samples
Association of uterine pathogens with blood and fecal microbiota
Comparison of blood, fecal, and uterine microbiota with vaginal microbiota
Absolute quantification of bacteria
To assess if the microbiota from one body site affected the microbiota from another site, we evaluated the correlations between the 16S rRNA gene copies in blood and uterus and between feces and uterus in cows that had samples from all three sites. We observed a tendency for a negative correlation (Spearman’s r s = −0.69, P = 0.06) in the 16S rRNA gene copies between the uterus and blood (Additional file 6: Figure S3), meaning that blood bacteria decreased as uterine bacteria increased. This result also supports the idea of bacterial transport from the blood to the uterus. The correlations between feces and uterus (Spearman’s r s = 0.60, P = 0.12) and between feces and blood (Spearman’s rs = 0.05, P = 0.91) were not significant. The lack of correlation is probably because of the immensely greater bacterial biomass in the gut compared with other body sites; therefore, transfer of bacteria from the gut to the uterus may be more related to the integrity of physical barriers such as the vulva, vagina, and cervix, and transfer of bacteria from the gut to blood may be more related to the integrity of the gastrointestinal epithelium.
Previous metagenomic studies identified Bacteroides heparinolyticus and Fusobacterium necrophorum as uterine pathogens because these bacteria were more prevalent in the uterus of dairy cows with metritis, compared with those in that of the healthy cows [5, 26]. Therefore, to determine if uterine pathogens are present in blood and feces and, if they are, how many of them occupy their habitats, we used ddPCR for copy numbers of B. heparinolyticus and F. necrophorum (Fig. 5b). The mean copy numbers of B. heparinolyticus in uterine samples were 2.46 ± 0.28 (log) on day 0 and 2.82 ± 0.27 (log) on day 2, which were significantly greater (P < 0.01) than those in the blood and feces, with mean copy numbers of 1.15 ± 0.10 (log) on day 0 and 1.14 ± 0.11 (log) on day 2 for blood and mean copy numbers of 1.23 ± 0.10 (log) on day 0 and 1.16 ± 0.11 (log) on day 2 for feces. With regard to F. necrophorum, there was no statistical difference (P > 0.30) on day 0 among sample types, where the mean copy number was 1.69 ± 0.17 (log) for blood, 1.88 ± 0.19 (log) for feces, and 1.55 ± 0.28 (log) for uterine samples. Meanwhile, on day 2, the abundance of F. necrophorum was elevated in the uterus at 2.46 ± 0.43 (log), which was significantly more abundant (P = 0.03) than in the blood at 1.53 ± 0.22 (log) and tended to be higher (P < 0.10) than in the feces at 1.81 ± 0.11 (log). Both B. heparinolyticus and F. necrophorum showed an increasing trend in the uterus, although not significant. On the other hand, B. heparinolyticus and F. necrophorum remained at low abundance in blood and feces. Taken together, uterine pathogens were detected in blood and feces, but the uterine environment in early postpartum seems to provide the ideal conditions for them to thrive.
The current paradigm of the origin of uterine bacteria is that physical barriers are compromised during parturition, which allows for bacteria to ascend the genital tract from the vagina or through the vagina from the environment as well as the animal’s skin and feces [8, 9]. Nonetheless, as stated before, uterine pathogens such as F. necrophorum and Trueperella pyogenes  cause liver abscess and Helcococcus ovis causes valvular endocarditis ; therefore, a hematogenous route of colonization of the uterus cannot be discounted. Previously, we reported that the uterus had an established microbiome within 60 min of calving , which indicates that colonization occurred before or shortly after calving either by ascending contamination from or through the vagina, via the bloodstream, or both. Herein, we showed that blood harbored a unique microbiome that contained the main uterine pathogens such as Bacteroides, Fusobacterium, and Porphyromonas. The presence of uterine pathogens in feces and blood indicate the feasibility of a hematogenous spread of bacteria from the gut to the uterus. Nevertheless, ascending contamination of the uterus cannot be discarded because the vagina also harbors the main uterine pathogens. Interestingly, other uterine pathogens such as Prevotella, Helcococcus, Filifactor, Campylobacter, and Arcanobacterium were not part of the core vaginal microbiome. In addition, vaginal microbiota was distinct from uterine microbiota on NMDS plot (Fig. 4a), although both microbiota showed a similarity with fecal microbiota. Thus, these results indicate that both vaginal and uterine microbiota are influenced by the gut microbiota, as suggested in a previous study . Nonetheless, it is worth noting that vaginal samples were collected from a different group of cows, which limits our ability to perform direct comparisons between the uterine and vaginal microbiomes. Therefore, sampling within the same animals is required in future studies to confirm our findings.
Previous studies have shown that network analysis is a powerful tool to investigate microbial interactions in complex environments such as soil and water [29, 30]. Thus, we applied network analysis to our samples to find bacterial genera that are important in the structure of their microbiota. Interestingly, uterine pathogens such as Bacteroides, Porphyromonas, and Fusobacterium formed similar networks in blood and uterus, despite the difference in microbiota abundance and composition (Additional file 7: Figure S4). The same was not observed in the vagina or feces. It is not clear at this point why this network is formed in blood, but it may be an important factor for their transmission to other body sites such as the uterus or liver. Of particular interest was the fact that Coxiella was found to be part of the blood network that included uterine pathogens. Coxiella is a bacterium that infects and multiplies inside of monocytes and has a tropism for the uterus [5, 31]; therefore, it is likely that the influx of leukocytes to the uterine lumen after calving contributes to the high prevalence of Coxiella on day 0. The correlation of Coxiella with uterine pathogens indicates that uterine pathogens may also be transported and transferred into the uterine lumen by monocytes. Although gut microbes were observed in blood leukocytes in cows in mid-lactation, uterine pathogens were not found in all cows ; therefore, presence of uterine pathogens in blood leukocytes shortly after calving should be investigated. Indeed, a study showed that bacterial translocation to extraintestinal sites was more common in late gestation and shortly after parturition than in the non-pregnant state or in early or mid-lactation . Nonetheless, free-floating bacteria in blood could also be transferred to the uterus because of the degenerative vascular changes that occur shortly after calving .
Herein, we used ddPCR to quantify uterine-specific pathogens and all bacteria. Droplet digital PCR is considered a third-generation PCR technology and has been shown to be more accurate than the real-time PCR [32, 33]. We found that the total number of bacteria was lower in the uterus than in blood and feces on day 0 (Fig. 5a), in spite of high species richness and diversity (Additional file 2: Figure S2). This shows that blood is a reasonably abundant source of bacteria for seeding of other tissues. The observation of higher abundance of B. heparinolyticus in the uterus than in blood or feces at calving indicates a tropism of this bacterium for the uterus. It is not clear how this bacterium concentrates in the uterus but the same was not observed for F. necrophorum, although F. necrophrum also showed an adaptation for the uterine environment early postpartum. The negative correlation in total bacteria between the uterus and blood (Additional file 6: Figure S3) supports the idea of bacterial transport from the blood to the uterus via blood leukocytes, but the synergism among uterine pathogens in blood and their specific mechanisms of invasion of the uterus warrants further investigation.
The blood, feces, vagina, and uterus have unique environments, and thereby, unique structures of bacterial communities were observed depending on body sites. Nonetheless, ordination and cluster analysis revealed that fecal bacterial communities are closely related to vaginal and uterine bacterial communities. Additionally, high abundance of core genera was shared by blood, fecal, vaginal, and uterine samples. More importantly, major uterine pathogens such as Bacteroides, Porphyromonas, and Fusobacterium and other uterine pathogens such as Prevotella and Helcococcus were part of the core genera in blood samples. Interestingly, although major uterine pathogens such as Bacteroides, Porphyromonas, and Fusobacterium were also part of the core genera in vaginal samples, other uterine pathogens such as Prevotella and Helcococcus were not part of the core genera in vaginal samples, which indicates that the blood may be the most important route of transmission of some uterine pathogens. Furthermore, uterine pathogens formed similar networks in blood and uterus, which may be an important factor for transmission, and warrants further investigation. The copy number of total bacteria in blood was correlated with the total bacteria in the uterus. On the other hand, the copy number of total bacteria in feces and uterus and feces and blood were not correlated. The copy number of total bacteria was higher in feces than in blood and uterus. In contrast, B. heparinolyticus was more abundant in the uterus on day 0, and both B. heparinolyticus and F. necrophorum were more abundant in the uterus than in the blood and in the feces on day 2. Our findings indicate that bacteria originating from the gut may be translocated to the uterus via the bloodstream. This study shows the feasibility of hematogenous spread of uterine pathogens in cows, although it does not exclude the possibility of direct fecal contamination or contamination from the vagina. In fact, direct fecal contamination or contamination from the vagina are likely to occur as well.
Animals and sampling
Holstein cows (n = 12) from a commercial dairy in Central Florida milking 5000 cows were used in this study. Blood, feces, and uterine samples were collected from the same individuals quickly after calving (within 60 min of calving; mean time = 20 min; SD = 14 min) to avoid or minimize the chance for contamination of the uterus from the environment ascending through the vagina. As part of the routine management, cows were changed from a high-fiber to a high-concentrate diet after calving, which leads to changes in the rumen and fecal microbiome , which could lead to changes in the blood microbiome. Therefore, we collected samples 2 days after calving to capture shifts in each separate microbiome. Cows were followed until 8 days postpartum for the diagnosis of metritis, characterized by fetid red-brownish watery uterine discharge, as previously reported . Because there are no major differences in the uterine microbiome between healthy and metritic cows at calving and at 2 days postpartum  and because we wanted to focus on the source of the uterine contamination rather than the differences between metritic and healthy cows, we included six cows that later developed metritis and six cows that remained healthy in the study. All the cows were healthy at the time of sampling because the clinical signs of metritis did not develop until 6 ± 2 days postpartum .
Uterine samples were collected using a sterile swab (Har-VetTM McCullough Double-Guarded Uterine Culture Swab) as previously described . Blood samples were collected from the jugular vein using vacutainer tubes with EDTA after surgically prepping a 150-cm2 area over the vein with iodine scrub and alcohol-soaked gauze pads. Fecal samples were collected from the rectum using sterile cotton-tipped swabs. Samples were stored at − 80 °C until DNA extraction.
Because this study aimed to confirm our hypothesis that gut bacteria could be transported to the uterus via the bloodstream, we did not sample the vagina. Nonetheless, for completeness, we compared our metagenomic data with vaginal data from a previous study as the reference . Vaginal samples were collected on day − 7 (7 days prepartum) from 105 Holstein dairy cows that were different from the cows used for collection of other samples, using sterile cotton-tipped swabs.
The gDNA was isolated from uterine swabs using the QIAamp DNA Mini kit (Qiagen), from 400 μl of blood using the QIAamp DNA Blood Mini kit (Qiagen), and from 200 mg of feces using the QIAamp DNA Stool Mini kit (Qiagen). The steps were performed as directed by the manufacturer with a modification; all samples were incubated with 400 μg of lysozyme for 1 h at 37 °C to maximize bacterial DNA extraction. The DNA concentrations of samples were measured using NanoDrop® ND-2000 (NanoDrop Technologies). The gDNA from vaginal swabs was extracted using the PowerSoil DNA Isolation Kit (MO BIO Laboratory Inc., Carlsbad, CA) after disruption of the sample using a bead beater homogenizer (Mini-Beadbeater-8, Biospec Products).
All the samples were sequenced, by our collaborators from Cornell University, using the same technique. The V4 hypervariable region of the 16S rRNA gene was amplified by PCR as previously described . PCR products were tagged with a 12 bp error-correcting Golay barcodes. The 5′ barcoded amplicons were prepared in triplicate using 10 μM of primer 515F and 806R, 1× GoTaq Green Master Mix (Promega), 1 mM MgCl2, and DNA template as follows: an initial denaturing step at 94 °C for 3 min, followed by 35 cycles of 94 °C for 45 s, 50 °C for 1 min, and 72 °C for 90 s, and a final elongation step at 72 °C for 10 min. Replicate amplicons were pooled and purified with a QIAquick PCR Purification Kit (Qiagen), followed by electrophoresis to visualize PCR products. Purified amplicon was quantified using the Quant-iT™ PicoGreen® dsDNA Assay Kit (Life Technologies Corporation) to normalize the concentration of all DNA libraries. Normalized libraries were pooled and sequenced using the MiSeq reagent kit V2-300 cycles on the MiSeq platform (Illumina Inc.). The reads were demultiplexed and filtered in each sample, allowing a single mismatch in index recognition, a quality score of 30, and a minimum length of 100 nt. Taxonomy was assigned using the Metagenomics workflow based on an Illumina-curated version of the Greengenes database.
Metagenomic and statistical analysis
Sequencing depth was evaluated using rarefaction curves in the Metagenomics RAST (https://metagenomics.anl.gov/) with the following parameters: annotation source Greengenes, maximum e-value cutoff 1e− 5, minimum identity % cutoff 97%, and minimum alignment length cutoff 100 bp. The Chao1 and Shannon indices were calculated in the R “fossil” and “vegan” packages, respectively (http://www.r-project.org). The relative abundance of bacterial phyla was compared between day 0 and day 2 using the Wilcoxon test. To represent the distance between samples, the NMDS of Bray-Curtis dissimilarity was carried out using PAST3 (http://folk.uio.no/ohammer/past/), in which non-parametric multivariate analysis of variance (PERMANOVA) was used to test significant difference among groups. To measure a difference in bacterial communities between groups and to identify which taxa are primarily responsible for the difference, the SIMPER analysis was conducted in the PAST3. To compare microbiota composition, the UPGMA clustering analysis was performed at the genus level based on Jaccard similarity using PAST3. The Jaccard similarity index was calculated in pairwise comparisons of the communities. Venn diagrams showing the number of core genera in blood, fecal, and uterine samples were created using the Bioinformatics & Evolutionary Genomics (http://bioinformatics.psb.ugent.be/webtools/Venn/). To understand the interrelationships of core genera within body habitats, co-occurrence patterns of core genera were evaluated in the network interface using pairwise Spearman’s rank correlations based on bacterial abundance. Strong (Spearman’s r s < − 0.7 or r s > 0.7) and significant (P < 0.01) correlations between core genera were considered a valid co-occurrence event. In the network, the nodes represent core genera and edges indicate relations among nodes. The topology of the network including average node connectivity, clustering coefficient, and modularity was calculated  and was visualized in the ForceAtlas2 algorithm using the Gephi (http://gephi.org) [37, 38].
Droplet digital PCR
Absolute quantification of bacteria was examined by ddPCR using a DNA binding dye (EvaGreen) according to the manufacturer’s instructions. Universal 16S primers which were designed by Clifford et al.  were used for the identification of all bacteria, and species-specific primers for B. heparinolyticus and F. necrophorum were designed in this study using Primer3 (Additional file 8: Table S4). The ddPCR reaction mixture contained 10 μl Supermix (Bio-Rad), 250 nM primers, and gDNA (~ 40 ng) in a final volume of 20 μl and combined with 20 μl of droplet generation oil (Bio-Rad), which partitioned into approximately 20,000 droplets in the QX200 droplet generator (Bio-Rad). The droplets generated from each sample were amplified by PCR on the PTC-100 (Bio-Rad) with the following condition: 95 °C for 10 min, 40 cycles of 94 °C for 30 s, 55 °C for 45 s, and 72 °C for 50 s, followed by 72 °C for 5 min and a hold at 4 °C. All samples were run in duplicate. After amplification, each droplet was read by the QX200 droplet reader (Bio-Rad) to count the number of positive and negative droplets, and target DNA molecules were presented as copies per microliter by the QuantaSoft™ software (Bio-Rad). The original concentration of the target DNA was log10-transformed in copies per 1 μL of gDNA and was analyzed by ANOVA using JMP Pro 13. For statistical analysis, differences with P ≤ 0.05 were considered significant, and differences with 0.05 < P ≤ 0.10 were considered to have a tendency towards statistical significance.
The authors thank the owners and staff of Alliance dairy for allowing the use of their cows in this experiment. We are grateful to Drs. ML Bicalho and RC Bicalho for providing the 16S sequencing data from the vaginal samples.
This project was supported by the USDA-NIFA-CRIS program (Accession Number: 1002880).
Availability of data and materials
The metagenome sequences analyzed during the current study are available from the MG-RAST under the ID numbers. The detailed information is described in metadata (Additional file 9: Table S5).
SJ analyzed and interpreted the metagenomics data and was a major contributor in writing the manuscript; SJ and KG contributed to the design of the experiment; AV and KG collected the samples; SL and RB performed the sequencing; KG and FC revised the manuscript. All authors approved the final version of the manuscript to be published.
Ethics approval and consent to participate
This study was approved by the University of Florida Institutional Animal Care and Use Committee (IACUC Protocol No.: 201207405).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Dolea C, Stein C. Global burden of maternal sepsis in the year 2000. Geneva: World Health Organization; 2003.Google Scholar
- Benzaquen ME, Risco CA, Archbald LF, Melendez P, Thatcher MJ, Thatcher WW. Rectal temperature, calving-related factors, and the incidence of puerperal metritis in postpartum dairy cows. J Dairy Sci. 2007;90(6):2804–14.View ArticlePubMedGoogle Scholar
- Goshen T, Shpigel NY. Evaluation of intrauterine antibiotic treatment of clinical metritis and retained fetal membranes in dairy cows. Theriogenology. 2006;66(9):2210–8.View ArticlePubMedGoogle Scholar
- Lima FS, Vieira-Neto A, Vasconcellos GS, Mingoti RD, Karakaya E, Sole E, Bisinotto RS, Martinez N, Risco CA, Galvao KN, Santos JE. Efficacy of ampicillin trihydrate or ceftiofur hydrochloride for treatment of metritis and subsequent fertility in dairy cows. J Dairy Sci. 2014;97(9):5401–14.View ArticlePubMedGoogle Scholar
- Jeon SJ, Vieira-Neto A, Gobikrushanth M, Daetz R, Mingoti RD, Parize AC, de Freitas SL, da Costa AN, Bicalho RC, Lima S, Jeong KC, Galvao KN. Uterine microbiota progression from calving until establishment of metritis in dairy cows. Appl Environ Microbiol. 2015;81(18):6324–32.View ArticlePubMedPubMed CentralGoogle Scholar
- Machado VS, Bicalho ML, Meira Junior EB, Rossi R, Ribeiro BL, Lima S, Santos T, Kussler A, Foditsch C, Ganda EK, Oikonomou G, Cheong SH, Gilbert RO, Bicalho RC. Subcutaneous immunization with inactivated bacterial components and purified protein of Escherichia coli, Fusobacterium necrophorum and Trueperella pyogenes prevents puerperal metritis in Holstein dairy cows. PLoS One. 2014;9(3):e91734.View ArticlePubMedPubMed CentralGoogle Scholar
- Deng Q, Odhiambo JF, Farooq U, Lam T, Dunn SM, Ametaj BN. Intravaginal lactic acid bacteria modulated local and systemic immune responses and lowered the incidence of uterine infections in periparturient dairy cows. PLoS One. 2015;10(4):e0124167.View ArticlePubMedPubMed CentralGoogle Scholar
- Sheldon IM, Dobson H. Postpartum uterine health in cattle. Anim Reprod Sci. 2004;82–83:295–306.View ArticlePubMedGoogle Scholar
- Tharpe N. Postpregnancy genital tract and wound infections. J Midwifery Womens Health. 2008;53(3):236–46.View ArticlePubMedGoogle Scholar
- Jeon SJ, Cunha F, Ma X, Martinez N, Vieira-Neto A, Daetz R, Bicalho RC, Lima S, Santos JE, Jeong KC, Galvao KN. Uterine microbiota and immune parameters associated with fever in dairy cows with metritis. PLoS One. 2016;11(11):e0165740.View ArticlePubMedPubMed CentralGoogle Scholar
- Machado VS, Oikonomou G, Bicalho ML, Knauer WA, Gilbert R, Bicalho RC. Investigation of postpartum dairy cows’ uterine microbial diversity using metagenomic pyrosequencing of the 16S rRNA gene. Vet Microbiol. 2012;159(3–4):460–9.View ArticlePubMedGoogle Scholar
- Lima FS, Oikonomou G, Lima SF, Bicalho ML, Ganda EK, Filho JC, Lorenzo G, Trojacanec P, Bicalhoa RC. Prepartum and postpartum rumen fluid microbiomes: characterization and correlation with production traits in dairy cows. Appl Environ Microbiol. 2015;81(4):1327–37.View ArticlePubMedPubMed CentralGoogle Scholar
- Swartz JD, Lachman M, Westveer K, O’Neill T, Geary T, Kott RW, Berardinelli JG, Hatfield PG, Thomson JM, Roberts A, Yeoman CJ. Characterization of the vaginal microbiota of ewes and cows reveals a unique microbiota with low levels of Lactobacilli and near-neutral pH. Front Vet Sci. 2014;1:19.View ArticlePubMedPubMed CentralGoogle Scholar
- Bicalho ML, Santin T, Rodrigues MX, Marques CE, Lima SF, Bicalho RC. Dynamics of the microbiota found in the vaginas of dairy cows during the transition period: associations with uterine diseases and reproductive outcome. J Dairy Sci. 2017;100(4):3043–58.View ArticlePubMedGoogle Scholar
- Tadepalli S, Narayanan SK, Stewart GC, Chengappa MM, Nagaraja TG. Fusobacterium necrophorum: a ruminal bacterium that invades liver to cause abscesses in cattle. Anaerobe. 2009;15(1–2):36–43.View ArticlePubMedGoogle Scholar
- Kutzer P, Schulze C, Engelhardt A, Wieler LH, Nordhoff M. Helcococcus ovis, an emerging pathogen in bovine valvular endocarditis. J Clin Microbiol. 2008;46(10):3291–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Archbald LF, Schultz RH, Fahning ML, Kurtz HJ, Zemjanis R. A sequential histological study of the post-partum bovine uterus. J Reprod Fertil. 1972;29(1):133–6.View ArticlePubMedGoogle Scholar
- Gier HT, Marion GB. Uterus of the cow after parturition: involutional changes. Am J Vet Res. 1968;29(1):83–96.PubMedGoogle Scholar
- Perez PF, Dore J, Leclerc M, Levenez F, Benyacoub J, Serrant P, Segura-Roggero I, Schiffrin EJ, Donnet-Hughes A. Bacterial imprinting of the neonatal immune system: lessons from maternal cells? Pediatrics. 2007;119(3):e724–32.View ArticlePubMedGoogle Scholar
- Young W, Hine BC, Wallace OA, Callaghan M, Bibiloni R. Transfer of intestinal bacterial components to mammary secretions in the cow. PeerJ. 2015;3:e888.View ArticlePubMedPubMed CentralGoogle Scholar
- Thomson AJ, Telfer JF, Young A, Campbell S, Stewart CJ, Cameron IT, Greer IA, Norman JE. Leukocytes infiltrate the myometrium during human parturition: further evidence that labour is an inflammatory process. Hum Reprod. 1999;14(1):229–36.View ArticlePubMedGoogle Scholar
- Hansen PJ. Physiology and Endocrinology Symposium: maternal immunological adjustments to pregnancy and parturition in ruminants and possible implications for postpartum uterine health: is there a prepartum-postpartum nexus? J Anim Sci. 2013;91(4):1639–49.View ArticlePubMedGoogle Scholar
- de Oliveira MN, Jewell KA, Freitas FS, Benjamin LA, Totola MR, Borges AC, Moraes CA, Suen G. Characterizing the microbiota across the gastrointestinal tract of a Brazilian Nelore steer. Vet Microbiol. 2013;164(3–4):307–14.View ArticlePubMedGoogle Scholar
- Aagaard K, Ma J, Antony KM, Ganu R, Petrosino J, Versalovic J. The placenta harbors a unique microbiome. Sci Transl Med. 2014;6(237):237ra65.View ArticlePubMedPubMed CentralGoogle Scholar
- Blondel VD, Guillaume J, Lambiotte R, Lefebvre E. Fast unfolding of communities in large networks. J Stat Mech. 2008;10:10008.View ArticleGoogle Scholar
- Peng Y, Wang Y, Hang S, Zhu W. Microbial diversity in uterus of healthy and metritic postpartum Holstein dairy cows. Folia Microbiol (Praha). 2013;58(6):593–600.View ArticleGoogle Scholar
- Nagaraja TG, Chengappa MM. Liver abscesses in feedlot cattle: a review. J Anim Sci. 1998;76(1):287–98.View ArticlePubMedGoogle Scholar
- Laguardia-Nascimento M, Branco KM, Gasparini MR, Giannattasio-Ferraz S, Leite LR, Araujo FM, Salim AC, Nicoli JR, de Oliveira GC, Barbosa-Stancioli EF. Vaginal microbiome characterization of Nellore cattle using metagenomic analysis. PLoS One. 2015;10(11):e0143294.View ArticlePubMedPubMed CentralGoogle Scholar
- Barberan A, Bates ST, Casamayor EO, Fierer N. Using network analysis to explore co-occurrence patterns in soil microbial communities. ISME J. 2012;6(2):343–51.View ArticlePubMedGoogle Scholar
- Li B, Yang Y, Ma L, Ju F, Guo F, Tiedje JM, Zhang T. Metagenomic and network analysis reveal wide distribution and co-occurrence of environmental antibiotic resistance genes. ISME J. 2015;9(11):2490–502.View ArticlePubMedPubMed CentralGoogle Scholar
- Bernard H, Brockmann SO, Kleinkauf N, Klinc C, Wagner-Wiening C, Stark K, Jansen A. High seroprevalence of Coxiella burnetii antibodies in veterinarians associated with cattle obstetrics, Bavaria, 2009. Vector Borne Zoonotic Dis. 2012;12(7):552–7.View ArticlePubMedGoogle Scholar
- Kim TG, Jeong SY, Cho KS. Comparison of droplet digital PCR and quantitative real-time PCR for examining population dynamics of bacteria in soil. Appl Microbiol Biotechnol. 2014;98(13):6105–13.View ArticlePubMedGoogle Scholar
- Dreo T, Pirc M, Ramsak Z, Pavsic J, Milavec M, Zel J, Gruden K. Optimising droplet digital PCR analysis approaches for detection and quantification of bacteria: a case study of fire blight and potato brown rot. Anal Bioanal Chem. 2014;406(26):6513–28.View ArticlePubMedGoogle Scholar
- Plaizier JC, Li S, Danscher AM, Derakshani H, Andersen PH, Khafipour E. Changes in microbiota in rumen digesta and feces due to a grain-based subacute ruminal acidosis (SARA) challenge. Microb Ecol. 2017;74(2):485–95.Google Scholar
- Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Huntley J, Fierer N, Owens SM, Betley J, Fraser L, Bauer M, Gormley N, Gilbert JA, Smith G, Knight R. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 2012;6(8):1621–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Newman MEJ. The structure and function of complex networks. SIAM Rev. 2003;45:167–256.View ArticleGoogle Scholar
- Bastian M, Heymann S, Jacomy M. Gephi: an open source software for exploring and manipulating networks. ICWSM. 2009;8:361–2.Google Scholar
- Jacomy M, Venturini T, Heymann S, Bastian M. ForceAtlas2, a continuous graph layout algorithm for handy network visualization designed for the Gephi software. PLoS One. 2014;9(6):e98679.View ArticlePubMedPubMed CentralGoogle Scholar
- Clifford RJ, Milillo M, Prestwood J, Quintero R, Zurawski DV, Kwak YI, Waterman PE, Lesho EP, Mc GP. Detection of bacterial 16S rRNA and identification of four clinically important bacteria by real-time PCR. PLoS One. 2012;7(11):e48558.View ArticlePubMedPubMed CentralGoogle Scholar