- Open Access
Induction of mastitis by cow-to-mouse fecal and milk microbiota transplantation causes microbiome dysbiosis and genomic functional perturbation in mice
Animal Microbiome volume 4, Article number: 43 (2022)
Mastitis pathogenesis involves a wide range of opportunistic and apparently resident microorganims including bacteria, viruses and archaea. In dairy animals, microbes reside in the host, interact with environment and evade the host immune system, providing a potential for host-tropism to favor mastitis pathogenesis. To understand the host-tropism phenomena of bovine-tropic mastitis microbiomes, we developed a cow-to-mouse mastitis model.
A cow-to-mouse mastitis model was established by fecal microbiota transplantation (FMT) and milk microbiota transplantation (MMT) to pregnant mice to assess microbiome dysbiosis and genomic functional perturbations through shotgun whole metagenome sequencing (WMS) along with histopathological changes in mice mammary gland and colon tissues.
The cow-to-mouse FMT and MMT from clinical mastitis (CM) cows induced mastitis syndromes in mice as evidenced by histopathological changes in mammary gland and colon tissues. The WMS of 24 samples including six milk (CM = 3, healthy; H = 3), six fecal (CM = 4, H = 2) samples from cows, and six fecal (CM = 4, H = 2) and six mammary tissue (CM = 3, H = 3) samples from mice generating 517.14 million reads (average: 21.55 million reads/sample) mapped to 2191 bacterial, 94 viral and 54 archaeal genomes. The Kruskal–Wallis test revealed significant differences (p = 0.009) in diversity, composition, and relative abundances in microbiomes between CM- and H-metagenomes. These differences in microbiome composition were mostly represented by Pseudomonas aeruginosa, Lactobacillus crispatus, Klebsiella oxytoca, Enterococcus faecalis, Pantoea dispersa in CM-cows (feces and milk), and Muribaculum spp., Duncaniella spp., Muribaculum intestinale, Bifidobacterium animalis, Escherichia coli, Staphylococcus aureus, Massilia oculi, Ralstonia pickettii in CM-mice (feces and mammary tissues). Different species of Clostridia, Bacteroida, Actinobacteria, Flavobacteriia and Betaproteobacteria had a strong co-occurrence and positive correlation as the indicator species of murine mastitis. However, both CM cows and mice shared few mastitis-associated microbial taxa (1.14%) and functional pathways regardless of conservation of mastitis syndromes, indicating the higher discrepancy in mastitis-associated microbiomes among lactating mammals.
We successfully induced mastitis by FMT and MMT that resulted in microbiome dysbiosis and genomic functional perturbations in mice. This study induced mastitis in a mouse model through FMT and MMT, which might be useful for further studies- focused on pathogen(s) involved in mastitis, their cross-talk among themselves and the host.
Mastitis is one of the most prevalent infectious diseases in the dairy animals worldwide. It has a negative impact on the agro-economy due to reduced milk production, early culling, and high therapeutic costs [1,2,3,4]. The disease is caused by a wide range of apparently resident and opportunistic microbes including bacteria, viruses, and archaea of variable origin, where the severity and outcome of the disease depends on the cross-talk between host and pathogen [5,6,7,8]. During mastitis, an imbalance of microbiomes occurs due to the inclusion of the pathogenic and opportunistic microorganisms in susceptible mammary glands. This invading encroachment is favoured by the compromised immune status of the host and the interactions between the opportunistic pathogens and the resident microbiota of the mammary gland [5, 7,8,9,10,11]. Depending on the host–pathogen interactions [12,13,14], bovine mastitis can manifest clinical features of apparent changes in the color of milk, swelling, redness, warmth and pain in the affected udder, systemic symptoms like fever and anorexia and sometimes death due to toxemia [15, 16]. Despite intensive research and implementation of various strategies over the last few decades to manage mastitis, the control of this dairy disease is still elusive because of the dynamic changes in the etiology . Moreover, microbial communities colonizing the mammary gland have evolved several novel mechanisms to facilitate their opportunistic proliferation manifesting the pathogenesis of CM [5, 7, 9]. The microbiomes of bovine CM comprise both contagious udder pathogens including Staphylococcus aureus, Streptococcus agalactiae, Streptococcus dysgalactiae, Mycoplasma spp., and Corynebacterium bovis [7, 11], and environmental pathogens such as Escherichia coli, Klebsiella pneumoniae, Klebsiella oxytoca, Enterobacter aerogenes, Streptococcus dysgalactiae, and Streptococcus uberis [8, 17]. However, microbiomes of the mammary glands and milk are considered to be highly similar, and the origin of microbes in the milk could be deep inside the upper parts of the mammary gland. Therefore, it is very likely that many of these microbes can migrate to the mammary gland from extra-mammary sites such as the gut, via the entero-mammary axis and/or the environment [17, 18].
In the last decade, the advent of shotgun metagenomics and downstream bioinformatics analyses allow for robust surveys of host-microbiome interactions in dairy populations [5, 7, 10]. A large body of literature suggests that disruption of gut microbiota homeostasis is linked to metabolic, immune, gastrointestinal, and mammary gland diseases [5, 19,20,21,22]. Previous studies demonstrated that changes in the gut microbiota can result in proliferation of specific pathogenic bacteria which can thereafter enter the mammary gland through the entero-mammary pathway [23,24,25]. The colonization of bovine gut microbiota into germ free (GF) mouse model demonstrates the influence of bovine gut microbes on the development of mastitis [5, 24, 25]. The usefulness of a mouse mastitis model as a complementary tool to investigate differences between bovine associated coagulase-negative Staphylococci (CNS) has already been established . Ma et al.  recently reported that cow-to-mouse fecal microbiota transplantation can modulate intestinal microbiome dysbiosis, which was one cause of mastitis. In this study, we found that in addition to fecal microbiota, CM milk microbiota transplantation is one of the important cuases of mastitis. Therefore, the cows-to-mouse mastitis milk microbiota transplantations (FMT and MMT) might be useful for further molecular biological studies on mastitis using the mouse model. Previous studies reported that humanized gnotobiotic mouse models have contributed to advancements in biomedicine, bridging the gap between human and animal gut pathophysiology [25, 27, 28].
In the pathophysiology of bovine mastitis, dysbiosis of the milk microbiomes occurs with an increase in opportunistic pathogens, and a reduction in healthy milk commensal microbes [8, 29]. The opportunistic pathogens of the gut use very efficient strategies to evade host defenses in order to colonize and invade mammary tissues. Thereafter, damage occurs to mammary epithelial cells and disruption of the cow immune system causes clinical episodes of CM [7, 30]. Therefore, to investigate the host-tropism of bovine mastitis microbiomes using mouse model, this study attempted to explore whether dysbiosis of fecal and milk microbiota can lead to clinical episodes of mastitis in pregnant mice. In the present study, we induced mastitis in mice by transplantation of bovine fecal or milk microbiota from CM cows into GF pregnant mice and investigated the changes in the fecal and mammary tissue microbiota using a high-throughput WMS followed by metagenomic and metabolic functional analyses along with histopathological changes in mice mammary gland and colon tissues (Graphical abstract). Our study revealed that both CM cows and mastitis induced mice shared a limited number of microbial taxa and related genomic functional potentials. We also observed a high degree of discrepancy in the microbiomes in milk, mammary tissue and feces of the mastitis affected animals. The developed mouse model would be useful for identification of the primary pathogen(s) of bovine mastitis, and their molecular interactions with hosts that are critical for sustainable management of this economically important disease in the dairy industry worldwide. Although, the baseline data presented here are promising, further studies are recommended using a larger sample size, and with the inclusion of gut/rumen microbiome sampling in addition to the milk samples for direct testing of microbial transfer across entero-mammary axis to confirm the dysbiosis of microbiomes and associated genomic functional features.
Fecal or milk microbiome transplantion from mastitis cows induced mastitis in GF mice
To elucidate whether CM-associated cow microbiota could develop mastitis syndrome in mice, fecal or milk microbiota from seven CM and five H cows were transplanted to 40 GF pregnant mice. The sample groups were- CCMF: cow clinical mastitis feces; HCF: healthy cow feces; CCMM: cow clinical mastitis milk; HCM: healthy cow milk; MCMF: mouse clinical mastitis feces; HMF: healthy mouse feces; MCMMT: mouse clinical mastitis mammary tissue; and HMMT: healthy mouse mammary tissue (Additional file 1). Among the 40 recipient mice, 20 underwent to FMT and MMT from CM samples while the remaining 20 were treated with H cow fecal and milk microbiota (controls). By comparing the post-intervention (FMT and MMT) inflammatory responses among the four groups, we found that fecal microbiota from CM cows induced a higher incidence of mastitis (90.0%) in mice compared to milk microbiota (80.0%). However, clinical signs of CM were not apparent before 10 days-post challenges (i.e., before 27 days of gestation). Remarkably, none of the mice in the control groups (either H-cows FMT or MMT) developed mastitis syndromes. Comparison of murine post-microbiota transplantation inflammatory responses revealed that gut and milk microbiota from CM cows induced a much greater inflammatory (of mammary gland and colon) response than those from healthy cows. On mammary gland surface, severe gross pathological changes (e.g., swelling and congestion of mammary glands) that corresponded to CM were observed in the mastits groups (CCMF and CCMM), however no pathological changes were visually apparent in healthy groups i.e., HCF or HCM (Fig. 1A). The gross syndrome of CM was further supported and confirmed by histopathological changes in the mammary and colon tissues. The CM-associated histopathological changes included damage of the mammary gland tissue (e.g., loss of architecture of alveoli, haemorrhages, involuted alveoli, degenerative changes in the epithelium of the alveoli and ducts, mammary alveolus thickening, hyperemia, and edema) and extent of inflammatory cell infiltration (i.e., infiltration of polymorphonuclear cells). For example, observable histopathological lesions under hematoxylin–eosin staining includes broken lobules of the mammary gland, damaged acinuses, and destroyed epithelial cells, with inflammatory cells including macrophages, neutrophils, and blood cells detected in the mammary lobule, supporting connective tissue and linning of the epithelium (Fig. 1B); on the other hand, in control groups (HCF or HCM), no pathological changes were apparent. In addition, predominant changes in the colon tissues of the CM mice (Fig. 1C) included moderate to severe inflammatory cell infiltration into mucosa and submucosa, disorder in mucosal structure (epithelial necrosis, extension of the subepithelial space, and structural damage of villi). Conversely, the mice receiving H cows fecal and milk microbiota did not show such pathological changes in the colon tissues and exhibited normal intestinal mucosa with well-arranged villi structure (Fig. 1C). Moreover, the CM mice had significantly higher (p = 0.0013, Kruskal–Wallis test) pathological grade of injury in their mammary gland and colon tissues than those of H mice (Fig. 1D). Generally, healthy mammary tissue possesses more fat vacuoles than mastitis udder , and thus, occurrence of mastitis was not affected by the lactation periods.
Microbiome structure and composition in gut, milk and mammary tissues differed in cows and mice
To test for differences in taxonomic diversity and structure of microbiomes, we analyzed the WMS data using an open-source cloud-based metagenomic mapping-based method of IDSeq . We found significant differences in microbial α-diversity (i.e., within-sample diversity) estimated through observed species in each sample (i.e., richness), Chao1, Shannon, and Simpson indices in the CM and H samples regardless of the host. The species level α-diversity remained significantly higher in H control samples of both cows (p = 0.0079, Kruskal–Wallis test) and mice (p = 0.0068, Kruskal–Wallis test) compared to CM samples (Fig. 2A). The PCoA plot based on Bray–Curtis distances (Fig. 2B) and NMDS plot based on weighted UniFrac distances (Fig. 2C), showed significant microbial disparity between CM and H samples (p = 0.001, Kruskal–Wallis test) and hosts (i.e., cows and mice) (p = 0.005, Kruskal–Wallis test). The phylum-level microbiome composition also showed distinct differences (p = 0.002, Kruskal–Wallis test) across the detected microbial domains (i.e., bacteria, archaea and viruses) (Fig. 2D).
At the domain level, bacteria were the most abundant (99.78%) microbial community followed by viruses (0.13%) and archaea (0.09%) (Additional file 2). In this study, we detected 1731 bacterial species including 1590 and 979 in CM and H cow fecal samples, respectively (Fig. 3A), of which 43.44% species had sole association with CM (Fig. 3A, Additional file 3). Similarly, 618 species of bacteria were identified in cows milk including 592 and 79 species in CM and H milk samples, respectively, and the CM samples had sole association of 87.22% species (Fig. 3B, Additional file 3). Conversely, 1065 species including 853 in CM and 721 in H fecal samples of mice were detected (Fig. 3C), and of them, the CM fecal samples had sole association of 32.30% species (Fig. 3C). Likewise, 115 species of bacteria were identified in mouse mammary tissue including 69 in CM and 77 in H tissues (Fig. 3D), and the CM mice mammary tissues had sole association of 33.04% species (Fig. 3D, Additional file 3). Comparing the bacterial taxa in CM and H sample in both hosts, we found that only 1.14% and 0.65% bacterial species shared across the CM (Fig. 3E) and H (Fig. 3F) samples, respectively.
In addition to bacterial fraction of the microbiomes, we detected 94 viral (Additional file 4), and 54 archaeal (Additional file 5) genera from both cows and mice (Additional file 3). By comparing these genera across the sample categories, we found 78 and 62 genera in cow CM and H samples, respectively, and majority of these genera (75.0%) were found to be shared between CM and H samples. The milk samples of both CM and H cows harboured 29 viral genera (CM = 13 and H = 26) and, only 10.44% genera had sole association with CM. Likewise, 74 viral genera were detected in mice fecal samples (CM = 69, H = 57), of which 70.27% genera were shared between the conditions. The mammary tissue of mice possessed only 24 viral genera (CM = 20, H = 21), and among these genera, 70.84% were shared between the conditions. In this study, we found that fecal samples from both cows and mice harboured 54 archaeal genera (CM = 54, H = 54, in each category), and of them, 100% of the genera were found to be shared between the conditions (CM and H) (Additional file 5). Moreover, 32 and 46 archaeal genera were detected in cow’s milk and mice mammary tissues, respectively. In this study, the facal and milk samples of both CM and H cows were found to share 17 archaeal genera (Additional file 5) while the facal and mammary tissue samples of both CM and H mice shared 16 archaeal genera (Additional file 5).
Gut and milk bacteria of mastitis cows are distinct from those of healthy cows
To test associations between microbiomes (gut and milk) and clinical mastitis (CM), fecal and milk microbiota from seven crossbred Holstein cows diagnosed with CM were compared to five physically similar, age-matched, crossbred Holstein cows that served as the control (Additional file 1; Materials and method). We found significant differences (p = 0.009, Kruskal–Wallis test) in the relative abundance of the bacterial species in fecal and milk samples of dairy cows. Among the detected bacterial species, Pseudomonas aeruginosa (21.0%), Lactobacillus crispatus (16.10%), Klebsiella oxytoca (10.34%), Enterococcus faecalis (10.0%), Nocardia pseudobrasiliensis (5.0%), Lactobacillus vaginalis (4.85%), Clostridioides difficile (4.23%), Ralstonia insidiosa (4.0%), Bifidobacterium pseudolongum (3.80%), Muribaculum sp. (2.35%), Duncaniella sp. (2.32%), and Duncaniella dubosii (2.0%) were the top abundant species in the CM-fecal samples of cows (Fig. 3, Additional file 6). Conversely, P. aeruginosa (32.7%), E. faecalis (30.05%), Lachnospiraceae bacterium (3.03%), Clostridiales bacterium (2.28%), and Phocaeicola dorei (2.15%) had higher relative abundances in H-cow fecal samples compared to those of CM cows (Fig. 3, Additional file 6). Similarly, the CM-cow milk samples were dominated by Pantoea dispersa (24.29%), K. oxytoca (20.68%), Actinoalloteichus sp. (11.63%), N. pseudobrasiliensis (8.87%), Staphylococcus aureus (4.7%), C. botulinum (4.19%), Acinetobacter baumannii (4.0%), Acinetobacter johnsonii (3.5%), K. pneumoniae (3.44%) and Escherichia coli (3.02%). On the contrary, P. aeruginosa had the highest relative abundance (40.35%) in H-cow milk samples followed by P. dispersa (13.86%), A. baumannii (11.0%), Prevotella melaninogenica (10.0%), Actinoalloteichus sp. (5.25%), K. oxytoca (4.0%) and N. pseudobrasiliensis (3.05%). The remaining bacterial species detected from the metagenomic data had relatively lower (< 3.0%) abundances and remained mostly abundant in CM-associated fecal and milk samples of cows (Fig. 3, Additional file 6).
Cow-to-mouse FMT and MMT reveal distinct gut and mammary gland bacteria between mastitis and healthy mice
We further investigated whether cow-to-mouse FMT and MMT treatment could produce distinct disease outcomes among the challenged mice. The WMS of both fecal and mammary tissues obtained from seven-CM and five-H mice at Day 27 of gestation (Methods) showed distinct changes in both composition and relative abundances of bacteria at the species level. The CM-mice fecal samples had a higher number of bacterial species (n = 853) than H-mice fecal samples (n = 561). The CM-related fecal metagenome of mice (MCMF) was dominated by Muribaculum sp. (38.30%), Duncaniella sp. (10.17%), Muribaculum intestinale (9.61%), Bifidobacterium animalis (8.36%), D. dubosii (7.14%), E. faecalis (5.0%), Akkermansia muciniphila (4.42%), L. crispatus (2.87%), B. pseudolongum (2.15%) and Lactobacillus murinus (2.10%) (Fig. 3, Additional file 6). In contrast, the H-mice fecal metagenome (HMF) was dominated by Muribaculum sp. (26.74%), A. muciniphila (20.61%), Duncaniella sp. (7.89%), M. intestinale (6.82%), D. dubosii (4.71%), K. pneumoniae (3.5%), P. aeruginosa (3.2%), and E. faecalis (2.0%), and rest of the species detected in both groups had comparatively lower relative abundances than that of CM-mice (Fig. 4, Additional file 6). The mammary tissues of the challenged mice had a significantly lower number of bacterial species compared to fecal samples (115 vs. 1065). In the mammary tissue of CM mice (MCMMT), E. coli (42.48%) was identified as the single most prevalent species followed by S. aureus (9.7%), Massilia oculi (5.90%), Ralstonia pickettii (4.13%), Curtobacterium flaccumfaciens (3.84%), P. aeruginosa (2.9%), A. johnsonii (2.37%), A. junii (2.36%), and Cutibacterium acnes (2.06%) (Fig. 4). Conversely, the H mice mammary tissue metagenome (HMMT) was mostly dominated by Ralstonia pickettii (9.84%), S. aureus (9.53%), A. johnsonii (8.85%), P. aeruginosa (8.74%), Helicobacter cinaedi (7.049%), M. oculi (5.74%), A. junii (4.51%), and E. coli (3.28%). The remaining species detected in both MCMMT and HMMT metagenomes had relatively lower (< 3.0%) abundances (Fig. 4, Additional file 6).
Indicator and shared bacterial taxa in the mastitis and healthy mice metagenomes
The indicator species analysis identified 46 differentially abundant (IndVal values ≥ 0.6, p < 0.01) bacterial species in mouse CM (MCMF and MCMMT) and H (HMF and HMMT) metagenomes (Fig. 5). In this study, mice CM samples showed the highest number of differentially enriched species (including 24 species in MCMF samples and eight species in MCMMT samples). Among the indicator species identified, Paenibacillus durus, Stenotrophomonas maltophilia, Pontibacter russatus, Caproiciproducens sp. NJN-50, Pseudobutyrivibrio xylanivorans, Treponema brennaborense, Capnocytophaga sputigena, Roseimicrobium sp. ORNL1, Christensenella massiliensis, and Blautia obeum were the top scoring (IndVal ≥ 0.8, p = 0.002) species in MCMF samples. Similarly, P. polymyxa, Niabella ginsenosidivorans, Hymenobacter sedentarius and several species that were classified at higher taxonomic ranking (IndVal ≥ 0.7, p = 0.03) were differentially abundant in the MCMMT samples. Remarkably, indicator species analysis confirmed that P. durus, S. maltophilia, P. russatus and Paenibacillus polymyxa were all good indicators of the murine CM (Fig. 5, Additional file 2). On the contrary, the H-mice samples had only seven differentially enriched bacterial species. The most abundant species (IndVal ≥ 0.78, p = 0.01) in the HMF samples were Bifidobacterium choerinum, Desulfovibrio sp. IOR2, and M. intestinale while Azoarcus olearius, Bacteroides sp. A1C1, B. uniformis and B. animalis were enriched (IndVal ≥ 0.79, p = 0.01) in HMMT samples (Fig. 5).
A microbial co-occurrence network was built based on correlations of relative abundance of the indicator species between mice-CM and -H samples. The network analysis (Fig. 6) presented 46 nodes and 1449 edges (significant positive correlations, r > 0.6 and Spearman’s corrected p = 0.001), all connected into one cluster with a clustering coefficient of 0.503. In the co-occurrence network, Clostridia had the highest number of edges (19.57%) followed by Bacteroida (13.04%), Actinobacteria (10.87%), Flavobacteriia (8.7%), Betaproteobacteria (8.7%), Gammaproteobacteria (6.52%), and Spirochaetia (6.52%), indicating strong co-occurrence among the species of these classes (Fig. 6). Pseudoclostridium thermosuccinogenes and Pseudobutyrivibrio xylanivorans (phylum: Firmicutes, class: Clostridia) had the highest co-occurrence (Correlation coefficient: 0.55, p = 0.005, in each) showing positive correlation with 45 indicator species and had only a negative correlation with Azoarcus species. Likewise, Labilibaculum antarcticum (phylum: Bacteroidetes, class: Bacteroida) was correlated with 25 bacterial species including P. russatus, Polynucleobacter necessaries, Treponema succinifaciens etc. and had no negative correlation with other indicator species (Fig. 6). In contrast, Azoarcus sp. and Lacrimispora sphenoides were mutually exclusive in the mouse habitats, and their abundance in mouse gut and mammary tissues was negatively correlated (Fig. 6).
CM-associated changes in viral and archaeal fraction (genus-level) of microbiomes in mice
Consistent with the variation in the bacterial component of the microbiomes, we concurrently found notable differences in the relative abundances of the viral (Fig. 7, Additional file 4) and archaeal (Fig. 8, Additional file 5) components in CM and H samples of both cows and mice (Additional file 2). For instance, the CM mice fecal samples (MCMF) were enriched with higher relative abundances of Siphovirus (61.98%), Mastadenovirus (14.17%) and Myovirus (12.21%), and Gammaretrovirus (37.77%), whereas the CM mice mammary tissues (MCMMT) were predominated by Ichnovirus (21.90%), Betaretrovirus (15.45%), and Macavirus (12.69%) (Fig. 7, Additional file 2). In contrast, Siphovirus (61.98%), Mastadenovirus (14.17%) and Myovirus (12.21%) in H mice fecal samples (HMF), and Gammaretrovirus (36.12%), Ichnovirus (23.0%), Betaretrovirus (14.43%), Macavirus (13.11%) and Rhadinovirus (5.04%) in H mice mammary tissues (HMMMT) had higher relative abundances compared to the CM samples. The remaining viral genera in the CM and H samples of mice had relatively lower abundances (< 5.0%) (Fig. 7, Additional file 2).
The present microbiome study demonstrated significant differences (p = 0.012, Kruskal–Wallis test) among the archaeal community in the sample categories of both hosts (cow and mice) with mastitis at the genus level. The CM cow fecal metagenome (CCMF) was dominated by Aeropyrum (11.12%), Ferroplasma (3.81%), Cenarchaeum (3.72%), Desulfurococcus (3.31%), Haloarcula (2.90%), Halobacterium (2.77%), Halogeometricum (2.33%) and Methanobrevibacter (1.06%) genera whereas Methanosphaerula (33.05%), Methanoculleus (5.08%), Methanoplanus (2.54%) and Pyrococcus (2.54%) were identified as the top abundant archaeal genera in CM cows milk samples (CCMM) (Fig. 8, Additional file 7). Among the identified archaeal genera in the mastitis induced mice metagenomes, Acidilobus (23.74%), Aciduliprofundum (14.47%), Archaeoglobus (8.34%), Aeropyrum (6.93%), Caldivirga (5.12%), Candidatus (3.04%), Cenarchaeum (2.65%), Desulfurococcus (2.54%), Ferroglobus (2.51%), Ferroplasma (2.26%), Halalkalicoccus (2.26%), Haloarcula (2.01%), Haloferax (1.98%), Halobacterium (1.79%), Halogeometricum (1.67%), Halomicrobium (1.63%) in MCMF metagenome, and Methanosaeta (46.23%), Methanoplanus (22.64%), Methanosphaerula (%11.32%), Methanosarcina (3.77%), Acidilobus (1.89%), Aciduliprofundum (1.89%), Methanobrevibacter (1.89%), Methanoculleus (1.89%), and Pyrococcus (1.89%) in MCMMT metagenome were the most dominant genera (Fig. 8, Additional files 2, 7). Conversely in healthy mice, Aciduliprofundum (15.02%) in HMF, and Archaeoglobus (11.03%), Halobacterium (4.26%), Halalkalicoccus (3.01%), Halogeometricum (2.76%), Haloquadratum (2.51%), Methanoculleus (2.51%), Halomicrobium (2.01%) and Ignicoccus (2.01%) in HMMT were the most abundant archaeal genera. Although, the remaining archaeal genera detected in this study had a relatively lower abundance (< 1.0%), their abundance always remained higher in CM associated murine samples (Fig. 8, Additional files 2, 7).
CM-associated genomic functional perturbations of microbiomes in murine mastitis
The WMS data were further analyzed using an assembly-based hybrid method of MG-RAST 4.0 (MR)  to compare the genomic functional potentials of the microbiomes. In MR analysis, 90.29 million reads (43.15% of cleaned reads) mapped to putative genes with known protein functions after filtering against host associated reads (Data S1). We identified 154 differentially abundant KEGG orthologues; KOs (including MCMF = 149, HMF = 144, MCMMT = 55, HMMT = 109, CCMF = 142, HCF = 131, CCMM = 110 and HCM = 129) and 61 SEED functions (including MCMF = 59, HMF = 59, MCMMT = 44, HMMT = 56, CCMF = 59, HCF = 61, CCMM = 56 and HCM = 60) at different subsystem levels across the bovine and murine metagenomes (CM and H). By comparing the composition and relative abundances of the different Kos or SEED functions in the same KEGG pathway or SEED subsystem between CM- and H-metagenomes, we found significant differences (p = 0.003, Kruskal Wallis test) in their relative abundances indicating positive correlations with CM in both hosts (Fig. 9, Additional file 2). Moreover, by measuring the relative abunadances of these functional pathways between CM- and H-mice, our analysis revealed that bacterial chemotaxis (67.18%), primary immunodeficiency; ADA (32.11%), methanogenesis (29.16%), phosphotransferase system (17.88%) in CM-mice feces (MCMF), and glycolysis and gluconeogenesis (61.90%), reactive oxygen species; ROS (43.90%), cell-to-cell communication (43.15%), ABC transporters (40.25%), recombination-activating proteins; RAG1/RAG2 (34.21%), chemotaxis protein; motB (23.76%), Jak-STAT signaling pathway (32.23%) and one-carbon metabolism (18.75%) in CM-mice mammary tissues (MCMMT) were the predominantly enriched metabolic functions compared to their H counterparts (Fig. 9, Additional file 8). On the contrary, genes coding for quorum sensing: autoinducer-2 synthesis (46.81%) and HIF-1 signalling pathway (21.16%) in H-mice feces (HMF), and flagellar assembly (46.08%), bacterial secretion system (39.25%), and citrate synthase; gltA (15.61%) in H-mice mammary tissues (HMMT) were the enriched metabolic functions compared to the CM associated microbiomes. Moreover, remaining KOs and SEED modules identified in this study also varied in their relative abundances and had relatively lower abundance in H metagenomes (Fig. 9, Additional file 8).
Relationships between predominant microbial species and their genomic functional potentials in murine mastitis
In this study, B. pseudolongum, Bacteroides thetaiotaomicron, D. dubosii, Duncaniella sp. B8, Faecalibaculum rodentium, L. bacterium, L. murinus and Muribaculum sp. TLL-A4 had strongest positive correlations (Spearman’s correlation; r > 0.5, p < 0.01, Fig. 10) with different SEED functions including membrane transport, quorum sensing and biofilm formation, oxidative stress, rubrerythrin mediated oxidative stress, regulation of oxidative stress response, protein YjgK cluster linked to biofilm formation, virulence, diseases and defenses, proteolytic pathways, GTP or GMP pathways, methanogenesis, glycolysis and gluconeogenesis, and one-carbon metabolism (Fig. 10). Conversely, C. botulinum, A. johnsonii and Actinoalloteichus sp. AHMU revealed significant negative correlations (Spearman’s correlation; r ≥ − 0.4, p < 0.05) with most of the SEED modules (Fig. 10). Simultaneously, A. muciniphila, B. pseudolongum, B. thetaiotaomicron, B. animalis, D. dubosii, Duncaniella sp. B8, F. rodentium, L. bacterium, L. murinus and Muribaculum sp. TLL-A4 showed significantly higher positive correlations (Spearman’s correlation; r ≥ 0.6, p < 0.01) with KOs like cheV, motB, ABC transporters, phosphotransferase and two-component systems, ADA, RAG1/RAG2 and mcp (Additional file 9). Consistent with SEED functions, C. botulinum, M. oculi, A. johnsonii, A. junii, Actinoalloteichus sp. AHMU and H. cinaedi displayed negative correlations (Spearman’s correlation; r ≥ − 0.4, p < 0.05) with most of the KOs identified in this study (Additional file 9).
Metagenomics is a powerful tool for shedding light on microbiome signature and concurrent genomic features associated with the process of animal disease. Our cow-to-mice induced mastitis model developed in this study clearly indicated that microbiome dysbiosis and concurrent genomic functional perturbations are associated with mammalian mastitis. In this study, more than 8.5 million metagenomic reads were assigned to taxonomic composition which is thought to be sufficient enough to capture maximum microbial richness and corroborated with our previous studies [7, 8, 12]. Several previous studies [34, 35] reported that in shotgun metagenome (WMS) study, more than 90% sequence reads may come from host DNA supporting our present findings. One of the hallmark findings of this study is the identification of several new bacterial species associated with both bovine (e.g., L. crispatus, N. pseudobrasiliensis, L. vaginalis, C. difficile, R. insidiosa, B. pseudolongum, Muribaculum sp., Duncaniella sp., D. dubosii., and Actinoalloteichus sp.) (Fig. 4) and murine (e.g., Muribaculum sp., Duncaniella sp., M. intestinale, B. animalis, D. dubosii, A. muciniphila, L. crispatus, B. pseudolongum, L. murinus, M. oculi, R. pickettii, C. flaccumfaciens, and C. acnes) mastitis (Fig. 5). So far, none of these species have been reported to be associated with bovine and murine CM and other lactating mammals.
Microbial dysbiosis and host-tropism occurs in mastitis pathogenesis
In the present study, the fecal microbiota of cows and mice differed significantly at the species level. The prevalence of these microbiotas significantly varied across the CM and H metagenomes of both the hosts. The association of P. aeruginosa, K. oxytoca, P. dispersa, E. faecalis, S. aureus, A. baumannii, A. johnsonii, K. pneumoniae and E. coli in bovine CM corroborated with our previous studies [7, 8, 12]., However, the relative abundance of these species remained substantially higher in this study than what we reported in earlier studies [7, 8, 12]. Consistent with these findings of taxonomic discrepencies in microbiome signature and abundances between two respective sample categories (CM vs. H). We also observed the distinct changes in the genomic functional features of the microbiomes in the two hosts (cow and mice) and sample categories between groups (CM and H) (Fig. 9).
The findings of the present study revealed that majority of the species of Lactobacillus were identified with higher abundances in the CM fecal samples of both cows and mice, which support their positive association with CM as also reported previously in a number of studies [24, 36, 37]. Lactic-acid bacteria of the genera Lactobacillus, Bifidobacterium, Bacillus and Enterococcus are commonly found in the human , avian  and animals [37, 40] gut. Different species and/or strains of these genera are good candidates to compete with pathogens for mammary gland colonization . However, no direct scientific evidence is currently available regarding the typical role of these microbes in the pathophysiology of mastitis. K. pneumoniae and N. pseudobrasiliensis are the emerging opportunistic environmental pathogens, and transmission of these pathogens to the mammary glands might occur from feces and bedding materials [7, 8, 42]. E. faecalis, one of the major pathogens of bovine mastitis, generally present in organic bedding materials and feces being opportunistic invaders of mammary glands [17, 30]. Additionally, a novel species of Actinobacteria phylum, Actinoalloteichus sp. was predominantly found only in CM cows milk samples (not in mouse samples), which argues against Actinoalloteichus sp. as the causative agent of murine mastitis. However, Actinoalloteichus has recently been detected as the most abundant genus in milk  supporting our present findings. P. dispersa can cause a variety of infections in immunocompromised dairy cows [8, 10], however, P. melaninogenica, a member of the normal microbiota of the human breast tissue , has never been described in cow feces and milk until now. Notably in this study, microbial communities originating from H-control samples have differed from those of CM samples. For instance, the detection of P. aeruginosa, E. faecalis, P. dispersa, A. baumannii and P. melaninogenica at higher levels within normal H fecal and milk is intriguing (Fig. 4). Pseudomonas is an animal skin microbe, which can colonize the udders and/or quarters from contaminated water and fecal sources [12, 45] and thus, potentially act as opportunistic pathogen by causing intramammary infections [8, 10]. This study therefore revealed a close association between the gut microbiome and milk microbes in the pathogenesis of bovine CM as also reported previously [7, 8, 46]. Additionally, M. oculi which was predominantly identified in the mammary tissue of both CM and H mice had significantly lower abundance in fecal and milk samples of cows. This soil bacterium has recently been detected in human clinical specimen , however, this species has not been reported in either cows or mice.
Notably, we detected some organisms that are indicative of the mastitis in GF mice (Fig. 5). Our results showed that CM mice sample groups had the highest number of indicator bacterial species. For example, P. durus, S. maltophilia, P. russatus and P. polymyxa were identified as good indicators of the murine CM (IndVal ≥ 0.8, p = 0.002). This high indicator value suggested that these species were found in most of the samples from CM mice and were comparatively less abundant in H mice samples (Fig. 5). Moreover, different species of Paenibacillus have potential contributions to maintaining mammary gland homeostasis , and we therefore suggest that association of Paenibacillus and Pontibacter in mammalian mastitis could be as opportunistic pathogens. In addition, the multidrug resistant bovine CM pathogen, S. maltophilia [17, 48] was found with higher relative abundances in CM mice fecal samples indicating its potential role in the induction of mastitis in GF mice. We furhter found that different species of Clostridia, Bacteroida, Actinobacteria, Flavobacteriia and Betaproteobacteria had strong co-occurrence and positive correlation as indicator species of murine mastitis. In this study, P. thermosuccinogenes, P. xylanivorans, M. intestinale, B. uniformis, P. polymyxa, and B. animalis showed more connections and overlap suggesting their cooperative or syntrophic interactions during the pathogenesis of murine mastitis (Fig. 6). This is further supported by the higher relative abundance of these microbial taxa in CM mice samples.
Multi-microbial pathogenesis favors mastitis through metabolomic cross-talk
Mastitis is a polymicrobial (including bacteria, viruses, archaea) disease where both archaea and viruses are considered as traditionally neglected microbes. Unlike the bacteria, the diversity, composition, and the relative abundances of viruses (Fig. 7) and archaea (Fig. 8) remained much lower in this metagenome study. Our previous studies [7, 8, 10, 12] also reported the concurrent occurrence of archaea and viruses in bovine mastitis highlighting the notion that mastitis is a multietiological disease [8, 10, 29, 49]. Archaea are major colonisers of rumen or the intestinal tracts of animals or humans , and thus, their abundance in milk always remains much lower. Recently, majority of the archaeal genera we detected, have been identified from a diverse group of samples including the fecal sample of cattle , canine  and swine , bovine milk [7, 8] and saliva of human  supporting our present findings. Mastitis is the result of direct interactions of bacteria with hosts (host–pathogen interactions) under immunosuppression or stressed conditions (when the cow suffers from a severe negative energy balance at the onset, in-and-around lactation, and other environmental stress e.g., heat stress). During the progression of bacterial mastitis, viruses jump into, and reach the site of inflammation in the mammary tissues, and triggers pathogenic pathways by the lysis of macrophages . This phagocytic macrophage storming further aggravates the pathogenesis and creates a micro-aerobic/anaerobic condition which ultimately favors the archaeal growth [56, 57]. Therefore, both viruses and archaea may act as a predisposing factor as well as a primary etiological agent for more severe and prolonged mastitis [7, 8, 10, 12].
In this study, we also found alteration in relative abundances of some important predicted genomic functions among different sample groups of CM and H microbiomes (Fig. 9). The metabolic features identified in the same KEGG pathway or SEED subsystem varied across mastitis samples in both cows and mice, suggesting their possible association in the early colonization and disease progression [8, 17]. Higher abundance of genes associated with bacterial chemotaxis, two-component system, GTP or GMP, NF-kappa-B and Jak-STAT signalling pathways in CM microbiomes of both hosts suggests their potential roles in mastitis through several complex biologic processes including immune disorders, cell differentiation, migration, proliferation, expression of many cytokines, quorum-sensing, microbial group behaviours and oxidative stress mediators, which likely accounts for the high systemic pathogen burden [7, 58,59,60]. Moreover, the higher abundance of ROS in CM-samples may contribute to the development of oxidative stress and inflammatory response  to further aggravate the pathogenesis of mastitis. Besides, quorum sensing, biofilm formation, glutathione non-redox reactions, and methyl-accepting chemotaxis genes that were predominantly identified in CM-pathogens play an important role in many opportunistic bacterial infections [8, 30]. Conversely, the PI3K-Akt pathway related genes remained highly expressed in H-microbes, and this pathway plays a critical role in the regulation of cell growth, proliferation, survival, motility, differentiation, angiogenesis, and metabolism .
The observed differences in microbiome composition and their corresponding genomic functional properties are considered to be the co-selection factor for mammary gland pathogenesis. Notably, the enriched consortia of B. pseudolongum, B. thetaiotaomicron, D. dubosii, Duncaniella sp. B8, F. rodentium, L. bacterium, L. murinus, Muribaculum sp. TLL-A4, A. muciniphila and B. animalis species had the strongest positive correlations with most of the KEGG and SEED functional pathways, indicating these microbial members could play potential roles in the pathophysiology of mastitis (Fig. 10).
Cow-to-mouse model may be a useful protocol for mastitis diagonosis, curative and preventive studies
Laboratory animals like mice provide an effective experimental model for understanding the underlying mechanisms of host–microbe interactions [5, 63]. The mouse mastitis model seems to be a good model to study bovine mastitis compared to other laboratory animals for ease of handling, controlled environments, and low maintenance cost . Despite restoring the common mastitis syndromes in both cows and mice, we observed marked differences in microbiome composition and relative abundances irrespective of the sample categories. Therefore, the findings of this study imply that mastitis is not solely caused by the resident microbiota of the mammary gland and/or its secretory product milk but can also be impacted by the alterations of gut microbiota and their genomic functional potentials. Our results highlight the paramount importance of existing entero-mammary pathways through which gut and milk microbiota (from CM host only) could transfer to induce mastitis. The results of the present study corroborated with the previous findings of Ma et al.  who reported that bovine mastitis is not necessarily a local infection of mammary glands, rather can be caused by a dysfunctional intestinal microbiota. Taken together, the core microbiota identified in the present study distinguished gut and mammary microbiota not only with different sample sources, but also with the health and/or disease state of mice (CM vs. H). Therefore, restoration of gut and mammary gland ecosystem function, for example the mastitis-associated pathways identified in this study, could possibly serve as an effective therapeutic target for bovine mastitis, which may deserve further validation in a larger representative cohort of dairy cows. However, further research is needed to understand the mechanism that allows gut and environmental microorganisms to invade the mammary glands and/or quarters and their relationships with mammalian hosts immune system.
The omics approach employed to study mastitis pathogenesis clearly showed microbiome dysbiosis. Species bias was dependent on the host and the predicted genomic functional features in CM hosts were significantly different from the H control counterparts. Our results show that few mastitis-associated microbial taxa and/or genomic functions were shared between diseased (CM) cows and mice regardless of conservation of mastitis syndromes. Cow-to-mouse microbiota transplantation protocols for induction of mastitis might be useful for further molecular studies of mastitis, which will ultimately improve the prevention and treatment strategies in both human and animal species. Taken together, a high-level abundance of the dominant and indicator microbial communities, associated genomic functional potentials, and their simultaneous correlations with the pathogenesis of mammary glands are considered to be driving factors for the mammalian mastitis.
Dairy cow selection and sampling
Twelve (n = 12) Holstein crossbred cows (including CM = 7 and H = 5) from seven dairy farms in Dhaka district (23.81 N, 90.41 E) of Bangladesh were used as the donors of fecal and milk samples. The age, parity and lactation of the cows ragend from 2.5 to 6 years, 1 to 5, and 7 to 45 days, respectively (Additional file 1). The cows gave birth randomly throughout the year (no particular control breeding), were milked once daily with their calves used for stimulating milk letdown. The cows were fed on rice straw, cut-and-carry grasses and milling by-products as concentrate (crashed rice and/or sometimes mustard oil cake) with limited grazing . California mastitis test (CMT) was employed initially to diagnose CM following a previously published protocol , and manufacturer’s instruction (CMT®, Original Schalm reagent, ThechniVet, USA). In brief, about 2 mL of milk sample was squirted in each cup of mastitis paddle, and an equal volume (2 mL) of CMT reagent was added to the cups. The reactions were developed within 20 s in positive samples and was categorized into five grades based on gel formation (Scandinavian scoring system) in the reaction mixture viz. 0 (negative), T (trace, possible infections), 1 (weak positive), 2 (distinct positive), and + 3 (strong positive) . The cows having a CMT score of ≥ 2 along with gross visible signs of mastitis were designated as CM cows . A total of 12 milk samples (including CM = 7 and H = 5) were collected from lactating cows. Approximately, 15–20 mL of milk from each cow was collected in a sterile falcon tube during the morning milking (9.0–11.0 am) with an emphasis on pre-sampling disinfection of teat-ends and hygiene during milk sampling . Simultaneously, we collected 12 fresh fecal samples (including CM = 7 and H = 5) from the selected cows under hygienic condition. Aprroximately, 5–10 g of fecal sample was collected through the rectum of each cow wearing a hygienic disposable plastic glove. No lubricants were used during sample collection. Collected fecal samples were then mixed-well and placed in a sterile falcon tube after proper labelling (15 mL) . For either the CM or H group, both milk and fecal samples were freshly collected, transported to the laboratory keeping in an ice box (at 4 °C temp), the content was thereafter processed and divided into aliquots. A portion of the processed samples proceeded to FMT and MMT immediately after collection, and rest of the aliquoted samples were stored at − 80 °C for metagenomic DNA extraction.
Microbiota transfer experiments
For the FMT and MMT procedures, forty (n = 40) timed-pregnant (at Day 17 of breeding) Swiss albino mice were procured from the ICDDR’B, Dhaka, Bangladesh. The mice were randomly divided into four groups: Group-I (FMT from CM cows, n = 10) Group-II (FMT from H cows, n = 10), Group-III (MMT from CM cows, n = 10) and Group-IV (MMT from H cows, n = 10). The mice were challenged on the same day of sampling. At the day of challenge (Day 17 after mating), 0.5 g fecal sample obtained from each of the CM and H cow was resuspended with twice the fecal volume of sterile physiological saline. After thorough mixing and resting (to minimize the number of bacteria lost), the supernatant was collected, and FMT was performed by a single oral administration of 1 g/kg fecal suspension to each mouse of Group-I and Group-II . Like wise, 10 mL whole milk tubes (from both CM and H cows) were centrifuged at room temperature for 10 min at 5000× g . After centrifugation, both the cream layer and supernatant liquid were removed. The pellet was resuspended to the initial sample volume of sterile physiological saline, and MMT was performed by a single oral administration of 0.5 mL of milk suspension to each mouse of Group-III and Group-IV . The mice were housed in GF environment on a 12 h light/dark cycle with unlimited access to food and water. In order to prevent cross-contamination of gut microbiota, the four groups of mice were physically separated into different GF isolators after challenge. Moreover, each mouse was housed in a separate cage with safe distance apart within each of the individual GF isolator, so as to prevent any island effects [66, 69]. At the end of Day 27 of mating (10 days of challenge), the mice were sacrificed, and fresh fecal samples, mammary and gut (colon) tissues were collected. The collected fecal samples were then mixed with freshly prepared phosphate buffered saline (PBS), and finally stored at − 80 °C until further processing and DNA extraction.
Collection of milk and subsequent milk somatic cell count (SCC) was not feasible in challenged mice, therefore histopathological examination was performed to assess the alterations and inflammatory changes of mammary gland and colon tissues during mastitis. The mammary gland and colon tissues were dissected from mice, kept in glass tubes, and fixed in 4% paraformaldehyde (PFA) (Sigma-Aldrich) for 24 h at 4 °C. Paraffin-embedded tissues were cut on a Leica Rotary Microtome (Leica Microsystems), placed onto SuperFrost Plus slides, and dried overnight at 37 °C. Sections were deparaffinized with xylene, and gradually rehydrated through graded alcohols for staining with standard hematoxylin and eosin (H & E staining) sectioned . To assess the degree of tissue injury to the mammary gland and colon, we used the Chiu Scoring System  in a blinded manner, where the number of inflammatory cells was counted in 12 randomly selected fields from each slide at a magnification of 400× . The slides were observed for severe, diffuse interstitial and/or alveolar infiltrate of inflammatory cells, focal to multifocal areas of tissue damage, epithelial abnormalities and extensive necrotic areas in mastitic mice . The degree of necrosis in mammary gland tissues was scored on a scale of 0–3 (normal 0, mild 1, moderate 2, severe 3). The degree of colon injury was scored as grades 0 (normal mucosa), 1 (development of subepithelial spaces at villus tips), 2 (extension of the subepithelial space with moderate lifting of the epithelial layer), 3 (massive epithelial lifting with a few denuded villi), 4 (denuded villi with exposed capillaries), and 5 (disintegration of the lamina propria, ulceration, and hemorrhage) [66, 70]. The slides were assessed under an Olympus BX51 upright microscope (40× objective), and finally images were collected using an Olympus DP73 camera through cellSens entry software (Olympus Corporation, Japan), and visualized using image J software (https://imagej.nih.gov/ij/).
Genomic DNA extraction and whole metagenome sequencing
Total genomic DNA from six milk (CM = 3, H = 3), six mammary tissue (CM = 3, H = 3), and 12 fecal (CM = 8, H = 4, from both cow and mice) samples (Data S1) was extracted using Maxwell 16 automated DNA extraction platform (Promega Corporation, Madison, WI 53711-5399, USA). In brief, DNA from fecal (400 µL fecal suspension) and milk (400 µL of whole milk) samples was extracted using Maxwell® 16 FFS Nucleic Acid Extraction Kit  and Maxwell® 16 FFPE Plus LEV DNA Purification Kit , respectively. For DNA extraction from fecal samples, 80 mg fecal material was added to 400 μL lysis buffer and mixed thoroughly by vortexing to make a suspension. The Maxwell® 16 Tissue DNA Purification Kit (Promega Corporation, Madison, WI 53711-5399, USA) was used for DNA extraction from mammary tissue (50 mg) samples following previously published protocol . DNA quantity and purity were determined using NanoDrop ND-2000 spectrophotometer (ThermoFisher, USA) by measuring 260/280 absorbance ratio. Libraries for shotgun WMS were prepared with Nextera XT DNA Library Preparation Kit  according to the manufacturer’s instructions, and paired-end (2 × 150 bp) sequencing was performed using a NovaSeq 6000 sequencer (Illumina Inc., USA). Our metagenomic DNA yielded 517.14 million reads with an average of 21.55 million (maximum = 22.96 million, minimum = 18.76 million) reads per sample (Data S1).
WMS data processing and microbiome analysis
The generated FASTQ files were concatenated and filtered through BBDuk (available from https://sourceforge.net/projects/bbmap/) with options k = 21, mink = 6, ktrim = r, ftm = 5, qtrim = rl, trimq = 20, minlen = 30, overwrite = true  to remove Illumina adapters, known Illumina artifacts, and phiX. Any sequence below these thresholds or reads containing more than one ‘N’ were discarded. The WMS data were analyzed using both open-source cloud-based metagenomic mapping-based and assembly-based hybrid methods of IDSeq  and MG-RAST 4.0 (MR) , respectively. IDseq is an open-source cloud-based pipeline for taxonomic assignments with NTL (NTL; nucleotide alignment length in bp) ≥ 50 and NT % identy ≥ 90. In IDSeq analysis, a ‘target’ genome library was constructed containing all prokaryotic sequences from the NCBI Database. The WMS reads were then aligned against the target libraries using the very sensitive Bowtie 2 algorithm , and filtered with Trimmomatic  to remove the reads aligned with the cattle (bosTau8), mouse (GRCm39) and human (hg38) genomes. The raw sequences were simultaneously uploaded to the MR server with properly embedded metadata, and were subjected to quality filtering with dereplication, host DNA removal (filtering against the set reference genome of both cattle, mouse and human), and finally screening for metabolic functional assignment. We used minimum identity of ≥ 90% for metabolic functional analysis through KEGG pathways and SEED subsystems in the MR pipeline. In IDSeq pipeline, 4.06% reads (of total cleaned reads) mapped to the target reference genomes whereas in MR pipeline 43.15% reads mapped to different known protein functions in KEGG pathways and SEED subsystems after filtering the cow, mouse and human genomes (Additional file 2).
Alpha diversity (diversity within samples) was estimated using the observed species, Chao1, Shannon and Simpson diversity indices  for IDSeq read assignments and counts. To visualize differences in microbiome diversity, a principal coordinate analysis (PCoA) based on the Bray–Curtis distance method , and non-metric multidimensional scaling (NMDS) measured by weighted-UniFrac distance on IDSeq data were performed through Phyloseq R package, version 4.1 . Indicator species specific to a given sample group in mouse model (having ≥ 1000 reads assigned to a taxon) were identified based on the normalized abundances of species using the R package, indicspecies , and the significant indicator value (IndVal) index was calculated by the 999-permutation test. Data were processed through the Phyloseq R package, visualized by using ggplot2 . The Venn diagrams representing taxonomic composition were generated through a stand-alone software tool; FunRich (http://www.funrich.org/). The differences in the microbiome abundances (Z-scores) across the study metagenomes was calculated using Spearman’s correlation test. Z-scores were calculated to construct heatmap for showing the relative abundance of the microbes in each sample group with the formula: Z = (x − μ)/σ, where x is the relative abundance of microbes in each sample category, μ is the mean value of relative abundances of microbes in all samples, and σ is the standard deviation of relative abundances . Finally, a heatmap of normalized Z-scores (from 0 to 40) of relative abundances of bacterial populations was produced. In addition, the Spearman’s correlation coefficient and significance tests for the indicator species were performed using the R package Hmisc (https://cran.r-project.org/web/packages/Hmisc/index.html). A correlation network was constructed and visualized with Gephi (ver. 0.9.2) (https://gephi.org/) to explore the co-occurrence patterns of the indicator species.
Metabolic functional potential analysis
The genomic functional profile of the microbiomes was annotated according to the Kyoto Encyclopedia of Genes and Genomes (KEGG) Orthology , and SEED subsystem  databases in MR pipeline using a “Best Hit Classification” method. The functional mapping was performed with the partially modified set parameters (e-value cutoff: 1 × 10−30, min. % identity cutoff: ≥ 90%, and min. alignment length cutoff: 20) of the MR server . We also investigated the relationships between the relative abundances of the predominating bacterial species and their genomic functional potentials in the pathogenesis of murine mastitis after renormalization and permutations.
The pair-wise non-parametric Kruskal–Wallis rank sum test was used to evaluate differences in the relative percent abundance of the microbial taxa and differentially abundant SEED or KEGG functions (at different levels) in CM and H animal (cow and mice) groups. Indicator species analysis calculated an IndVal which was the product of the relative frequency and relative abundance of a species in a cluster. To test the significance of the IndVal, p values were calculated with 100 iterations, where in each iteration, the sample groupings were randomly assigned and an IndVal determined. The p values for the IndVal calculation were corrected for multiple comparisons using the false discovery rate correction . To explore the relationship among the relative abundance of indicator species, we calculated Spearman’s rank correlation coefficients using the ‘multtest’ package in R , and the correction was made using Benjamini–Hochberg FDR. In the co-occurrence network, each node represents one species and each edge stands for correlation between the species abundances . Spearman’s correlation among the top abundant 30 bacterial species, and their genomic metabolic functions (SEED and KEGG functional pathways were performed in Hmisc and corrplot R packages as described above (Methods).
Availability of data and materials
The sequence data and related metadata reported in this paper are available in the SRA repository of the NCBI database under BioProject Accession ID of PRJNA753312. Supplementary information supporting the findings of the study are available in this article as Additional files 1, 2, 3, 4, 5, 6, 7, 8 and 9.
Fecal microbiota transplantation
Milk microbiota transplantation
Whole shotgun metagenomic sequencing
Cow clinical mastitis feces
Healthy cow feces
Cow clinical mastitis milk
Healthy cow milk
Mouse clinical mastitis feces
Healthy mouse feces
Mouse clinical mastitis mammary tissue
Healthy mouse mammary tissue
Thompson-Crispi K, Atalla H, Miglior F, Mallard BA. Bovine mastitis: frontiers in immunogenetics. Front Immunol. 2014;5:493.
Hoque MN, Das ZC, Talukder AK, Alam MS, Rahman ANMA. Different screening tests and milk somatic cell count for the prevalence of subclinical bovine mastitis in Bangladesh. Trop Anim Health Prod. 2015;47:79–86.
Hoque M, Das Z, Rahman A, Hoque M. Effect of administration of vitamin E, selenium and antimicrobial therapy on incidence of mastitis, productive and reproductive performances in dairy cows. Int J Vet Sci Med. 2016;4:63–70.
Al Amin M, Hoque MN, Siddiki AZ, Saha S, Kamal MM. Antimicrobial resistance situation in animal health of Bangladesh. Vet World. 2020;13:2713.
Ma C, Sun Z, Zeng B, Huang S, Zhao J, Zhang Y, Su X, Xu J, Wei H, Zhang H. Cow-to-mouse fecal transplantations suggest intestinal microbiome as one cause of mastitis. Microbiome. 2018;6:1–17.
Boix-Amorós A, Hernández-Aguilar MT, Artacho A, Collado MC, Mira A. Human milk microbiota in sub-acute lactational mastitis induces inflammation and undergoes changes in composition, diversity and load. Sci Rep. 2020;10:1–14.
Hoque MN, Istiaq A, Clement RA, Sultana M, Crandall KA, Siddiki AZ, Hossain MA. Metagenomic deep sequencing reveals association of microbiome signature with functional biases in bovine mastitis. Sci Rep. 2019;9:1–14.
Hoque MN, Istiaq A, Rahman MS, Islam MR, Anwar A, Siddiki AZ, Sultana M, Crandall KA, Hossain MA. Microbiome dynamics and genomic determinants of bovine mastitis. Genomics. 2020;112:5188–203.
Patel SH, Vaidya YH, Patel RJ, Pandit RJ, Joshi CG, Kunjadiya AP. Culture independent assessment of human milk microbial community in lactational mastitis. Sci Rep. 2017;7:1–11.
Hoque M, Sultana M, Hossain A. Dynamic changes in microbiome composition and genomic functional potentials in bovine mastitis. J Data Min Genomics Proteomics. 2021;12:232.
Oikonomou G, Bicalho ML, Meira E, Rossi RE, Foditsch C, Machado VS, Teixeira AGV, Santisteban C, Schukken YH, Bicalho RC. Microbiota of cow’s milk; distinguishing healthy, sub-clinically and clinically diseased quarters. PLoS ONE. 2014;9:e85904.
Hoque MN, Istiaq A, Clement RA, Gibson KM, Saha O, Islam OK, Abir RA, Sultana M, Siddiki A, Crandall KA. Insights into the resistome of bovine clinical mastitis microbiome, a key factor in disease complication. Front Microbiol. 2020;11:860.
Hoque M, Das Z, Rahman A, Haider M, Islam M. Molecular characterization of Staphylococcus aureus strains in bovine mastitis milk in Bangladesh. Int J Vet Sci Med. 2018;6:53–60.
Gryaznova MV, Syromyatnikov MY, Dvoretskaya YD, Solodskikh SA, Klimov NT, Mikhalev VI, Zimnikov VI, Mikhaylov EV, Popov VN. Microbiota of cow’s milk with udder pathologies. Microorganisms. 1974;2021:9.
Royster E, Wagner S. Treatment of mastitis in cattle. Vet Clin Food Anim Pract. 2015;31:17–46.
Gomes F, Saavedra MJ, Henriques M. Bovine mastitis disease/pathogenicity: evidence of the potential role of microbial biofilms. Pathogens Dis. 2016;74:ftw006.
Derakhshani H, Fehr KB, Sepehri S, Francoz D, De Buck J, Barkema HW, Plaizier JC, Khafipour E. Invited review: microbiota of the bovine udder: contributing factors and potential implications for udder health and mastitis susceptibility. J Dairy Sci. 2018;101:10605–25.
Porcellato D, Meisal R, Bombelli A, Narvhus JA. A core microbiota dominates a rich microbial diversity in the bovine udder and may indicate presence of dysbiosis. Sci Rep. 2020;10:1–14.
Kamada N, Seo S-U, Chen GY, Núñez G. Role of the gut microbiota in immunity and inflammatory disease. Nat Rev Immunol. 2013;13:321–35.
Cheng S, Ma X, Geng S, Jiang X, Li Y, Hu L, Li J, Wang Y, Han X. Fecal microbiota transplantation beneficially regulates intestinal mucosal autophagy and alleviates gut barrier injury. Msystems. 2018;3:e00137-e1118.
Amin N, Schwarzkopf S, Kinoshita A, Tröscher-Mußotter J, Dänicke S, Camarinha-Silva A, Huber K, Frahm J, Seifert J. Evolution of rumen and oral microbiota in calves is influenced by age and time of weaning. Anim Microbiome. 2021;3:1–15.
Chiu CY, Miller SA. Clinical metagenomics. Nat Rev Genet. 2019;20:341–55.
Zebeli Q, Ametaj B. Relationships between rumen lipopolysaccharide and mediators of inflammatory response with milk fat production and efficiency in dairy cows. J Dairy Sci. 2009;92:3800–9.
Ma C, Zhao J, Xi X, Ding J, Wang H, Zhang H, Kwok L. Bovine mastitis may be associated with the deprivation of gut Lactobacillus. Benef Microbes. 2016;7:95–102.
Hu X, Guo J, Zhao C, Jiang P, Maimai T, Yanyi L, Cao Y, Fu Y, Zhang N. The gut microbiota contributes to the development of Staphylococcus aureus-induced mastitis in mice. ISME J. 2020;14:1897–910.
Gogoi-Tiwari J, Williams V, Waryah CB, Costantino P, Al-Salami H, Mathavan S, Wells K, Tiwari HK, Hegde N, Isloor S. Mammary gland pathology subsequent to acute infection with strong versus weak biofilm forming Staphylococcus aureus bovine mastitis isolates: a pilot study using non-invasive mouse mastitis model. PLoS ONE. 2017;12:e0170668.
Weingarden AR, Vaughn BP. Intestinal microbiota, fecal microbiota transplantation, and inflammatory bowel disease. Gut Microbes. 2017;8:238–52.
Choi HH, Cho Y-S. Fecal microbiota transplantation: current applications, effectiveness, and future perspectives. Clin Endosc. 2016;49:257.
Jiménez E, de Andrés J, Manrique M, Pareja-Tobes P, Tobes R, Martínez-Blanch JF, Codoñer FM, Ramón D, Fernández L, Rodríguez JM. Metagenomic analysis of milk of healthy and mastitis-suffering women. J Hum Lact. 2015;31:406–15.
Gomes F, Saavedra MJ, Henriques M. Bovine mastitis disease/pathogenicity: evidence of the potential role of microbial biofilms. FEMS Pathogens Dis. 2016;74:ftw006.
Bianchi RM, Schwertz CI, de Cecco BS, Panziera W, De Lorenzo C, Heck LC, Snel GG, Lopes BC, da Silva FS, Pavarini SP. Pathological and microbiological characterization of mastitis in dairy cows. Trop Anim Health Prod. 2019;51:2057–66.
Kalantar KL, Carvalho T, de Bourcy CF, Dimitrov B, Dingle G, Egger R, Han J, Holmes OB, Juan Y-F, King R. IDseq—an open source cloud-based pipeline and analysis service for metagenomic pathogen detection and monitoring. Gigascience. 2020;9:giaa111.
Glass EM, Wilkening J, Wilke A, Antonopoulos D, Meyer F. Using the metagenomics RAST server (MG-RAST) for analyzing shotgun metagenomes. Cold Spring Harbor Protoc. 2010;2010:pdb. prot5368.
Pereira-Marques J, Hout A, Ferreira RM, Weber M, Pinto-Ribeiro I, Van Doorn L-J, Knetsch CW, Figueiredo C. Impact of host DNA and sequencing depth on the taxonomic resolution of whole metagenome sequencing for microbiome analysis. Front Microbiol. 2019;10:1277.
Ong CT, Ross EM, Boe-Hansen GB, Turni C, Hayes BJ, Tabor AE. overcoming host contamination in bovine vaginal metagenomic samples with nanopore adaptive sequencing. J Anim Sci. 2022;100:skab344.
Qiao J, Kwok L, Zhang J, Gao P, Zheng Y, Guo Z, Hou Q, Huo D, Huang W, Zhang H. Reduction of Lactobacillus in the milks of cows with subclinical mastitis. Benef Microbes. 2015;6:485–90.
Derakhshani H, Plaizier JC, De Buck J, Barkema HW, Khafipour E. Composition and co-occurrence patterns of the microbiota of different niches of the bovine mammary gland: potential associations with mastitis susceptibility, udder inflammation, and teat-end hyperkeratosis. Anim Microbiome. 2020;2:1–17.
Arokiyaraj S, Seo SS, Kwon M, Lee JK, Kim MK. Association of cervical microbial community with persistence, clearance and negativity of Human Papillomavirus in Korean women: a longitudinal study. Sci Rep. 2018;8:1–9.
Wei S, Morrison M, Yu Z. Bacterial census of poultry intestinal microbiome. Poult Sci. 2013;92:671–83.
Lin W-C, Ptak CP, Chang C-Y, Ian M-K, Chia M-Y, Chen T-H, Kuo C-J. Autochthonous lactic acid bacteria isolated from dairy cow feces exhibiting promising probiotic properties and in vitro antibacterial activity against foodborne pathogens in cattle. Front Vet Sci. 2020;7:239.
Bouchard DS, Seridan B, Saraoui T, Rault L, Germon P, Gonzalez-Moreno C, Nader-Macias FM, Baud D, François P, Chuat V. Lactic acid bacteria isolated from bovine mammary microbiota: potential allies against bovine mastitis. PLoS ONE. 2015;10:e0144831.
Chen W, Liu Y, Zhang L, Gu X, Liu G, Shahid M, Gao J, Ali T, Han B. Nocardia cyriacigeogica from bovine mastitis induced in vitro apoptosis of bovine mammary epithelial cells via activation of mitochondrial-caspase pathway. Front Cell Infect Microbiol. 2017;7:194.
Rubiola S, Chiesa F, Dalmasso A, Di Ciccio P, Civera T. Detection of antimicrobial resistance genes in the milk production environment: impact of host DNA and sequencing depth. Front Microbiol. 1983;2020:11.
Urbaniak C, Cummins J, Brackstone M, Macklaim JM, Gloor GB, Baban CK, Scott L, O’Hanlon DM, Burton JP, Francis KP. Microbiota of human breast tissue. Appl Environ Microbiol. 2014;80:3007–14.
Kuehn JS, Gorden PJ, Munro D, Rong R, Dong Q, Plummer PJ, Wang C, Phillips GJ. Bacterial community profiling of milk samples as a means to understand culture-negative bovine clinical mastitis. PLoS ONE. 2013;8:e61959.
Maga EA, Weimer BC, Murray JD. Dissecting the role of milk components on gut microbiota composition. Gut Microbes. 2013;4:136–9.
Song W, Wang S, Shen J, Zhu B. Complete genome sequence of Massilia oculi sp. Nov. CCUG 43427 T (= DSM 26321 T), the Type Strain of M. oculi, and comparison with genome sequences of other Massilia Strains. Curr Microbiol. 2019;76:1082–6.
Ohnishi M, Sawada T, Marumo K, Harada K, Hirose K, Shimizu A, Hayashimoto M, Sato R, Uchida N, Kato H. Antimicrobial susceptibility and genetic relatedness of bovine Stenotrophomonas maltophilia isolates from a mastitis outbreak. Lett Appl Microbiol. 2012;54:572–6.
Ojo-Okunola A, Nicol M, Du Toit E. Human breast milk bacteriome in health and disease. Nutrients. 2018;10:1643.
Seedorf H, Kittelmann S, Henderson G, Janssen PH. RIM-DB: a taxonomic framework for community structure analysis of methanogenic archaea from the rumen and other intestinal environments. PeerJ. 2014;2:e494.
Ngetich DK, Bett RC, Gachuiri CK, Kibegwa FM. Diversity of gut methanogens and functional enzymes associated with methane metabolism in smallholder dairy cattle. Res Square. 2022. https://doi.org/10.21203/rs.3.rs-1120235/v1.
Swanson KS, Dowd SE, Suchodolski JS, Middelbos IS, Vester BM, Barry KA, Nelson KE, Torralba M, Henrissat B, Coutinho PM. Phylogenetic and gene-centric metagenomics of the canine intestinal microbiome reveals similarities with humans and mice. ISME J. 2011;5:639–49.
Zeineldin M, Megahed A, Blair B, Aldridge B, Lowe J. Metagenomic analysis of the fecal archaeome in suckling piglets following perinatal tulathromycin metaphylaxis. Animals. 1825;2021:11.
Lettieri GM, Santiago LM, Lettieri GC, dos Anjos Borges LG, Marconatto L, De Oliveira LA, Damé-Teixeira N, Salles LP. Oral phenotype and salivary microbiome of individuals with papillon-lefèvre syndrome. Front Cell Infect Microbiol. 2021;11:720790.
Argov T, Sapir SR, Pasechnek A, Azulay G, Stadnyuk O, Rabinovich L, Sigal N, Borovok I, Herskovits AA. Coordination of cohabiting phage elements supports bacteria–phage cooperation. Nat Commun. 2019;10:1–14.
Ruiz L, García-Carral C, Rodriguez JM. Unfolding the human milk microbiome landscape in the omics era. Front Microbiol. 2019;10:1378.
Lurie-Weinberger MN, Gophna U. Archaea in and on the human body: health implications and future directions. PLoS Pathog. 2015;11:e1004833.
Hall CL, Lee VT. Cyclic-di-GMP regulation of virulence in bacterial pathogens. Wiley Interdiscip Rev RNA. 2018;9:e1454.
Khan MZ, Khan A, Xiao J, Ma J, Ma Y, Chen T, Shao D, Cao Z. Overview of research development on the role of NF-κB signaling in mastitis. Animals. 2020;10:1625.
Usman T, Yu Y, Liu C, Wang X, Zhang Q, Wang Y. Genetic effects of single nucleotide polymorphisms in JAK2 and STAT5A genes on susceptibility of Chinese Holsteins to mastitis. Mol Biol Rep. 2014;41:8293–301.
Turk R, Koledić M, Maćešić N, Benić M, Dobranić V, Đuričić D, Cvetnić L, Samardžija M. The role of oxidative stress and inflammatory response in the pathogenesis of mastitis in dairy cows. Mljekarstvo/Dairy. 2017;67:91–101.
Bamodu OA, Chang H-L, Ong J-R, Lee W-H, Yeh C-T, Tsai J-T. Elevated PDK1 expression drives PI3K/AKT/MTOR signaling promotes radiation-resistant and dedifferentiated phenotype of hepatocellular carcinoma. Cells. 2020;9:746.
Fritz JV, Desai MS, Shah P, Schneider JG, Wilmes P. From meta-omics to causality: experimental models for human microbiome research. Microbiome. 2013;1:1–15.
Nguyen TLA, Vieira-Silva S, Liston A, Raes J. How informative is the mouse for human gut microbiota research? Dis Model Mech. 2015;8:1–16.
Huque K, Sarker N. Feeds and feeding of livestock in Bangladesh: performance, constraints and options forward. Bangladesh J Anim Sci. 2014;43:1–10.
Jamali H, Barkema HW, Jacques M, Lavallée-Bourget E-M, Malouin F, Saini V, Stryhn H, Dufour S. Invited review: Incidence, risk factors, and effects of clinical mastitis recurrence in dairy cows. J Dairy Sci. 2018;101:4729–46.
Lima SF, Bicalho MLdS, Bicalho RC. Evaluation of milk sample fractions for characterization of milk microbiota from healthy and clinical mastitis cows. PLoS ONE. 2018;13:e0193671.
Wang Y, Lin L, Yin C, Othtani S, Aoyama K, Lu C, Sun X, Yoshikai Y. Oral administration of bovine milk from cows hyperimmunized with intestinal bacterin stimulates lamina propria T lymphocytes to produce Th1-biased cytokines in mice. Int J Mol Sci. 2014;15:5458–71.
Koren O, Goodrich JK, Cullender TC, Spor A, Laitinen K, Bäckhed HK, Gonzalez A, Werner JJ, Angenent LT, Knight R. Host remodeling of the gut microbiome and metabolic changes during pregnancy. Cell. 2012;150:470–80.
Bonfiglio R, Galli F, Varani M, Scimeca M, Borri F, Fazi S, Cicconi R, Mattei M, Campagna G, Schönberger T. Extensive histopathological characterization of inflamed bowel in the dextran sulfate sodium mouse model with emphasis on clinically relevant biomarkers and targets for drug development. Int J Mol Sci. 2028;2021:22.
Chiu C-J, McArdle AH, Brown R, Scott HJ, Gurd FN. Intestinal mucosal lesion in low-flow states: I. A morphological, hemodynamic, and metabolic reappraisal. Arch Surg. 1970;101:478–83.
Burbach K, Seifert J, Pieper DH, Camarinha-Silva A. Evaluation of DNA extraction kits and phylogenetic diversity of the porcine gastrointestinal tract based on Illumina sequencing of two hypervariable regions. Microbiologyopen. 2016;5:70–82.
Sisea C, Pamfil D. Comparison of DNA extraction methods for GMO analysis of food products. Bull USAMV-CN. 2007;63:64.
Sato MP, Ogura Y, Nakamura K, Nishida R, Gotoh Y, Hayashi M, Hisatsune J, Sugai M, Takehiko I, Hayashi T. Comparison of the sequencing bias of currently available library preparation kits for Illumina sequencing of bacterial genomes and metagenomes. DNA Res. 2019;26:391–8.
Rahman MS, Hoque MN, Puspo JA, Islam MR, Das N, Siddique MA, Hossain MA, Sultana M. Microbiome signature and diversity regulates the level of energy production under anaerobic condition. Sci Rep. 2021;11:1–23.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.
Koh H. An adaptive microbiome α-diversity-based association analysis method. Sci Rep. 2018;8:1–12.
Beck J, Holloway JD, Schwanghart W. Undersampling and the measurement of beta diversity. Methods Ecol Evol. 2013;4:370–82.
McMurdie PJ, Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE. 2013;8:e61217.
De Cáceres M, Legendre P, Wiser SK, Brotons L. Using species combinations in indicator value analyses. Methods Ecol Evol. 2012;3:973–82.
Wickham H. ggplot2: elegant graphics for data analysis. New York: Springer; 2016.
Yang H, Huang X, Fang S, Xin W, Huang L, Chen C. Uncovering the composition of microbial community structure and metagenomics among three gut locations in pigs with distinct fatness. Sci Rep. 2016;6:1–11.
Kanehisa M, Sato Y, Furumichi M, Morishima K, Tanabe M. New approach for understanding genome variations in KEGG. Nucleic Acids Res. 2019;47:D590–5.
Overbeek R, Begley T, Butler RM, Choudhuri JV, Chuang H-Y, Cohoon M, de Crécy-Lagard V, Diaz N, Disz T, Edwards R. The subsystems approach to genome annotation and its use in the project to annotate 1000 genomes. Nucleic Acids Res. 2005;33:5691–702.
Keegan KP, Glass EM, Meyer F. MG-RAST, a metagenomics service for analysis of microbial community structure and function. In: Microbial environmental genomics (MEG). Springer; 2016. pp. 207–233.
Techtmann SM, Fortney JL, Ayers KA, Joyner DC, Linley TD, Pfiffner SM, Hazen TC. The unique chemistry of Eastern Mediterranean water masses selects for distinct microbial communities by depth. PLoS ONE. 2015;10:e0120605.
Pollard KS, Dudoit S, van der Laan MJ. Multiple testing procedures: the multtest package and applications to genomics. In: Bioinformatics and computational biology solutions using R and bioconductor. Springer; 2005. pp. 249–271.
The first author of this manuscript, M.N.H. received his Ph.D. Fellowship from the Bangabandhu Science and Technology Fellowship Trust (BSTFT), Ministry of Science and Technology, Government of Bangladesh. The authors would like to thank Animal Research Farm (ARF), ICDDR’B for providing experimental mice and technical support in mouse model experiment. We also acknowledge the cooperation of the dairy farmers for allowing us to conduct the study in their farm.
This study was partially supported by the grants received from the Bangladesh Bureau of Educational Information and Statistics (BANBEIS), Ministry of Education, People’s Republic of Bangladesh (Grant No.: LS 2017314). The funding organization (BANBEIS) played no role in the design of the study; collection, analysis, and interpretation of data; or in writing the manuscript.
Ethics approval and consent to participate
The protocol for milk sample collection from dairy cows and in vivo mouse model challenge was approved by the Animal Experimentation Ethical Review Committee (AEERC), Faculty of Biological Sciences, University of Dhaka (under reference number: 79/Biol.Scs), and International Centre for Diarrhoeal Disease Research, Bangladesh (ICDDR’B). Permission was also obtained from the owners of sampled dairy farms. Animal care and treatment were conducted in accordance with the institutional guidelines and international laws and policies (Directive 2010/63/EU) on the protection of animals used for scientific purposes. Every effort was made to minimize animal suffering.
Consent for publication
Consent was obtained from participating farms prior to the start of the study.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1.
Study information (animal information, location, breeds, lactation, and parity).
Additional file 2.
Information on metagenomic data, microbial structure, composition and relative abundances in each metagenome.
Additional file 3.
Taxonomic structure of microbiomes identified in different sample groups.
Additional file 4.
Taxonomic composition of viruses.
Additional file 5.
Taxonomic composition of archaea.
Additional file 6.
Taxonomic information of top seventy bacterial species identified in different metagenomic groups.
Additional file 7.
Archaeal genera detected (with their relative abundances) in different metagenomes.
Additional file 8.
Changes in metabolic functional profile in the study metagenomes.
Additional file 9.
The correlation between predominantly abundant bacterial species and their genomic functional potentials in different KEGG orthologues (KOs) in CM and H control metagenomes of both cows and mice.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Hoque, M.N., Rahman, M.S., Islam, T. et al. Induction of mastitis by cow-to-mouse fecal and milk microbiota transplantation causes microbiome dysbiosis and genomic functional perturbation in mice. anim microbiome 4, 43 (2022). https://doi.org/10.1186/s42523-022-00193-w
- Mammary gland