Next Article in Journal
Molecular Detection of Multidrug Resistant Salmonella Species Isolated from Broiler Farm in Bangladesh
Next Article in Special Issue
Prevention of Intramammary Infections by Prepartum External Application of a Teat Dip Containing Lactic Acid Bacteria with Antimicrobial Properties in Dairy Heifers
Previous Article in Journal
Effects of Inactivated Mycobacterium bovis Vaccination on Molokai-Origin Wild Pigs Experimentally Infected with Virulent M. bovis
Previous Article in Special Issue
Witch Hazel Significantly Improves the Efficacy of Commercially Available Teat Dips
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptome Analysis of The Inflammatory Responses of Bovine Mammary Epithelial Cells: Exploring Immunomodulatory Target Genes for Bovine Mastitis

1
Food and Feed Immunology Group, Laboratory of Animal Products Chemistry, Graduate School of Agricultural Science, Tohoku University, Sendai 980-8572, Japan
2
Livestock Immunology Unit, International Education and Research Center for Food and Agricultural Immunology (CFAI), Graduate School of Agricultural Science, Tohoku University, Sendai 980-8572, Japan
3
Department of Medicine, Faculty of Veterinary Science, Bangladesh Agricultural University, Mymensingh-2202, Bangladesh
4
Laboratory of Immunobiotechnology, Reference Centre for Lactobacilli, (CERELA-CONICET), Tucuman 980-0845, Argentina
5
Scientific Computing Laboratory, Computer Science Department, Faculty of Exact Sciences and Technology, National University of Tucuman, Tucuman 980-0845, Argentina
6
Infection Immunity Unit, International Education and Research Center for Food and Agricultural Immunology (CFAI), Graduate School of Agricultural Science, Tohoku University, Sendai 980-8572, Japan
7
Cell Biology Laboratory, Graduate School of Agricultural Science, Tohoku University, Sendai 980-8572, Japan
8
Graduate School of Food, Agriculture and Environment, Miyagi University, Sendai 980-8572, Japan
9
Department of Infectious Diseases and Immunology, Faculty of Veterinary Medicine, Utrecht University, 3584 CL Utrecht, The Netherlands
10
Department of Veterinary Tropical Diseases, Faculty of Veterinary Science, University of Pretoria, Private bag X20, Hatfield 0028, South Africa
*
Authors to whom correspondence should be addressed.
These authors have contributed equally to this work.
Pathogens 2020, 9(3), 200; https://doi.org/10.3390/pathogens9030200
Submission received: 23 February 2020 / Revised: 6 March 2020 / Accepted: 7 March 2020 / Published: 9 March 2020
(This article belongs to the Collection Mastitis in Dairy Ruminants)

Abstract

:
Bovine mastitis is the inflammatory reaction of the mammary gland and is commonly caused by bacterial infections in high-yielding dairy cows. The detailed investigation of the immunotranscriptomic response of bovine mammary epithelial (BME) cells to pattern recognition receptors (PRRs) activation by microbial-associated molecular patterns (MAMPs) can be of great importance for understanding the innate immune defense mechanisms, and for exploring the immunomodulatory candidate genes. In this work, we investigated the transcriptome modifications of BME cells after the in vitro stimulation with Escherichia coli derived lipopolysaccharide (LPS) and heat-killed Staphylococcus aureus JE2 and S. aureus SA003. In addition, the effect of Pam3CSK4 (a synthetic triacylated lipopeptide that activates Toll-like receptor 2 (TLR2)), and the intracellular chemotactic protein cyclophilin A (CyPA), which is secreted by BME cells during mastitis, in the expression changes of selected cytokines and chemokines were evaluated by qPCR. Microarray analysis identified 447, 465 and 520 differentially expressed genes (DEGs) in the BME cells after LPS, S. aureus JE2 and S. aureus SA003 stimulation, respectively. A major differential response in the inflammatory gene expression was noticed between the stimulation of LPS and S. aureus strains. Unlike the S. aureus strains, LPS stimulation resulted in significant upregulation of CCL2, CXCL2, CXCL3, CXCL8, IL1α and IL1β, which were confirmed by qPCR analysis. Pam3CSK4 was not able to induce significant changes in the expression of cytokines and chemokines in challenged BME cells. The exogenous CyPA administration was able to upregulate CXCL2, CXCL3, CXCL8, IL1α and IL1β expression in BME cells indicating its ability to promote inflammation. The identification of transcriptional markers of mastitis specific for individual inflammatory factors such as LPS, Pam3CSK4 or CyPA, which can be evaluated in vitro in BME cells, may enable the development of novel diagnostics and/or immunomodulatory treatments, providing new tools for the effective management of mastitis in dairy cows. The results of this work are an advance in this regard.

1. Introduction

Bovine mastitis, defined as the inflammation of the mammary gland, is a frequent and contagious infectious disease of high-yielding cows. Alterations in milk composition, reduced milk quality and quantity, and increased costs for treatments, labor and culling account not only for huge economic losses but also create concerns for animal welfare and human health [1,2]. Mastitis is classified as clinical or subclinical based on the appearance of the inflammation of the udder [3]. Escherichia coli and Staphylococcus aureus are among the most prevalent Gram-negative and Gram-positive bacterial pathogens that cause mammary gland infection in dairy cows [4]. It has reported that E. coli infection results in clinical mastitis which is characterized by acute symptoms of inflammation in the milk collecting cistern and the teat by a reduced milk production and an elevated somatic cell count [5]. On the other hand, S. aureus is responsible for one-third of cases of clinical and subclinical mastitis in dairy cattle which is characterized by less severe inflammation and is sometimes asymptomatic [3]. The severity of mastitis largely depends on the patterns of interactions between invading pathogens and the bovine mammary epithelial (BME) cells [6]. Accumulated research revealed that Gram-negative bacteria provoke a strong inflammatory response through a vigorous stimulation of cytokine synthesis in the mammary gland, resulting in the activation of the local and systemic inflammatory response [5,7]. On the other hand, it was reported that Gram-positive bacteria elicit a much weaker immune reaction of the udder and generally no strong systemic immune response is detected [8,9]. Therefore, in-depth understanding of the pathogen-specific molecular mechanisms involved in the generation of mammary gland immune responses could be of great importance to explore and select effective control measures of specific pathogen-induced mastitis in dairy cows.
When pathogenic bacteria enter the udder lumen via the teat canal, they interact with BME cells in order to establish colonization. This pathogen-BME cells interaction results in the release of inflammatory mediators and chemo-attractants that recruit and stimulate immune cells which exert their antibacterial activities locally and amplify the inflammatory response [10,11]. Therefore, it is considered that BME cells stand at the frontline in the resistance against bacterial infections in mammary glands. A number of studies have shown that BME cells are able to sense bacteria or bacterial products, and that they react by up-regulating several sets of genes involved in the inflammatory response [12,13,14,15,16,17]. The innate immune response of mammary gland initiates through the recognition of microbes associated molecular patterns (MAMPs) by the patterns recognition receptor (PRRs), such as Toll-like receptors (TLRs) expressed in BME cells. The MAMPs-mediated activation of TLRs results in several downstream cell-signaling events that induce the expression of cytokine and chemokines and trigger inflammatory responses [12,13,14,15,16,17]. Although it has been demonstrated that the recognition of MAMPs by TLRs expressed in BME cells is a key event in the generation of mammary inflammation [12], detailed transcriptomic studies evaluating the response of those cells to TLRs activation has not been widely performed [18,19].
In vivo studies to uncover the mastitis-associated gene expression changes in BME cells of lactating mammary gland require the use of a large number of animals to obtain statistically robust results because these out-bred populations exhibit considerable genetic variation. Though short-term in vitro experiments using the primary cell cultures have some advantages of reflecting the appropriate mitogenic responses, isolation of primary epithelial cells from the mammary gland tissue of lactating cows is relatively difficult as compared to that of prepubertal animals [20]. On the other hand, the use of untransformed cell lines has the advantage over primary cultures in that they are able to replicate over several passages, keeping the same functional characteristics. In this work, we conducted transcriptome profiling of BME cells that were cloned from mammary tissue of a 200-day pregnant Holstein cow [21] after in vitro stimulation with either E. coli derived lipopolysaccharide (LPS) or heat-killed S. aureus JE2 or S. aureus SA003 strains. In addition, we compared the innate immune response triggered by LPS in BME cells with those induced by the stimulation with Pam3CSK4 (a synthetic triacylated lipopeptide that activates toll-like receptor 2 (TLR2)) and cyclophilin A (CyPA), an intracellular chemotaxis protein secreted by BME cells during inflammatory reactions which has being proposed as a marker for the early diagnosis of bovine mastitis [18,22,23,24,25]. The identification of transcriptional markers of mastitis specific for individual inflammatory factors such as LPS, Pam3CSK4 and CyPA, which can be evaluated in vitro in BME cells, may enable the development of novel diagnostics and/or immunomodulatory treatments, providing new tools for the effective management of mastitis in dairy cows.

2. Results

2.1. Expression Patterns of Toll-like Receptors in BME Cells

We first aimed to evaluate the expression patterns of the members of TLR family in BME cells. The results demonstrated that the ten members of bovine TLR family are expressed in this bovine epithelial cell line (Figure 1). Among the TLR members, TLR1, TLR3, TLR4 and TLR6 were strongly expressed with a copy number greater than 103.

2.2. Differential Gene Expressions in BME Cells after Stimulation with LPS and Heat-killed S. auresus Strains

A total of 1432 transcripts were differentially expressed in BME cells after stimulation with LPS or S. aureus JE2 or S. aureus SA003 (Table 1). LPS stimulation resulted in the differential regulation of 447 transcripts, of which 222 were upregulated and 207 were downregulated. Likewise, S. aureus JE2 stimulation as compared to the control showed 465 transcripts, with 227 and 238 transcripts upregulated and downregulated, respectively. Furthermore, S. aureus SA003 stimulation as compared to control showed 520 transcripts, with 226 and 294 transcripts upregulated and downregulated, respectively.
In addition, there were 139 differentially expressed genes shared by LPS, S. aureus JE2 and S. aureus SA003 stimulations. Among them, 36 were up- and 103 were downregulated genes (Figure 2A,B). A total of 189 differentially expressed genes were common between LPS and S. aureus JE2 stimulation, where 60 and 129 were up- and downregulated, respectively. In addition, 100 differentially expressed genes were common between LPS and S. aureus SA003. A total of 202 genes were shared by S. aureus JE2 and S. aureus SA003 stimulations, of which 68 were up- and 134 were downregulated genes, respectively (Figure 2A,B). Hierarchical clustering of differentially expressed genes (DEGs) indicate that LPS and S. aureus stimulation are differently clustered, and the two S. aureus strains also differed between them (Figure 2C). The DEGs were distributed into ten distinct clusters (Figure 2C). The full list of DEGs are provided with the Supplementary files Tables S1–S3.

2.3. Gene Ontology (GO) Terms Enriched by DEGs

The LPS-induced DEGs were found to be involved in the significant enrichment of GO biological terms including “Cell chemotaxis,” “Inflammatory response,” “Positive regulation of nitric-oxide synthetase biosynthesis process,” “Response to molecule of bacterial origin” and “Negative regulation of extrinsic apoptotic signaling pathway” (Table 2). The DEGs after S. aureus JE2 stimulation were predicted to be involved in the significant enrichment of GO biological terms, namely “Cell differentiation,” “Cell adhesion,” “Inflammatory response,” “Blood coagulation” and “Chemotaxis” while S. aureus SA003-induced DEGs are involved with the enrichment of GO terms including “Innate immune response,” “Cell surface receptor signaling pathway,” “Blood coagulation,” “signal transduction” and “Positive regulation of NF-kB signaling” (Table 2).

2.4. Differential Responses in Inflammatory Gene Gxpression in BME Cells after Stimulation with LPS and Heat-Killed and S. aureus Strains

Next, we extracted the inflammatory response related genes from the list of DEGs and compared their expression patterns between LPS and S. aureus stimulations. Results demonstrated that expression of IL1α, IL1β, CCL2, CXCL2, CXCL3, CXCL5, CXCL8, CXCL9 and C6 were significantly upregulated in BME cells following LPS stimulation while their expressions were not significantly modified after stimulation with either S. aureus strains (Table 3). The expression of CRP was increased after S. aureus stimulation while not modified after LPS stimulation. The expression of PPIA (also called CyPA) expression remained unchanged irrespective of stimulation (Table 3).

2.5. Inflammatory Gene Expression in BME Cells after LPS or Pam3CSK4 Stimulation

In order to confirm the microarray expression results, we measured the mRNA expression dynamics of six selected inflammatory genes after the stimulation of BME cells with LPS or Pam3CSK4 by using RT-qPCR. As shown in Figure 3, the activation of TLR4 significantly increased the expression of IL1α and IL1β that was in accordance to microarray expression results. The levels of IL1β continuously increased until 24 hours after stimulation, while Pam3CSK4 stimulation did not induce significant variations in the expression levels of IL1β when compared to basal levels, showing a marked contrast with the results observed for the LPS stimulation (Figure 3). Pam3CSK4 was able to induce significant increases in the expressions of IL1α in BME cells at hour 12 post-stimulation. However, the increases of IL1α in Pam3CSK4-treated BME cells were lower than the observed for LPS-treated cells. In addition, the expressions of the chemokines CCL2 (MCP-1), CXCL8 (IL8), CXCL2 (MIP-2α) and CXCL3 (MIP-2β) were evaluated in BME cells after LPS or Pam3CSK4 stimulation (Figure 3). The mRNA levels of CCL2 and CXCL8 continuously increased until 24 hours after LPS stimulation while CXCL2 and CXCL3 reached their highest values on hour six and then they maintained similar expression levels. Pam3CSK4 did not induce significant variations in the expression levels of four chemokines evaluated when compared to basal levels (Figure 3).

2.6. Inflammatory Gene Expression in BME Cells after CyPA Stimulation

Cyclophilins, also called immunophilins, are a group of proteins having peptidyl-prolyl cis-trans isomerase activity that are secreted from mammalian cells in response to infection, cellular hypoxia or oxidative stress [22,23]. CyPA is able to promote cell signaling activation through its interaction with the intracellular receptor for the immunosuppressive drug cyclosporine A [24] or the extracellular matrix metalloproteinase inducer (EMMPRIN, also called CD147) [22]. In a recent study, we observed that extracellular CyPA possesses chemotactic activity for bovine peripheral blood cells [25], which indicate its ability to influence the inflammatory immune responses. Then, we next aimed to evaluate whether CyPA might have the ability to modulate the innate immune response of BME cells. For this purpose, BME cells were stimulated with bovine recombinant CyPA at three different doses and the expression of inflammatory cytokines and chemokines were evaluated at different time points by qPCR. As shown in Figure 4, all the doses of CyPA significantly increased the expression of IL1α at hours 3 and 6, as well as IL1β at hour 24 in BME cells when compared with unstimulated cells. In addition, all the doses of CyPA were able to significantly increase the expression of CXCL3 at hour 6, while no variations were observed in CCL2 and CXCL2 (Figure 4). Significant increases of CXCL8 were detected in BME cells at hours 6 and 24 when 1000 and 100 ng/mL of CyPA were used, respectively (Figure 4).

3. Discussion

The generation of inflammatory response is a vital part of the innate defenses against pathogens infecting the mammary gland mucosa. BME cells have an important role in the activation of the mammary gland innate immunity through their ability to recognize pathogens and release of inflammatory mediators that initiate the inflammatory response [10,11]. After the recognition of MAMPs by the host PRRs such as TLRs, mammary epithelial cells initiate a downstream signal transduction, which in turn lead to the activation of transcription factors including NF-kB (12). Those transcription factors modulate the expression of genes on mammary epithelial cells and induce the secretion of proinflammatory cytokines and chemokines, type I interferons, and antimicrobial peptides which coordinately induce local and systemic inflammatory responses that allows the elimination of the pathogens [12]. The results of the present study demonstrated the expression of TLR family members in BME cells, supporting the hypothesis that BME cells are capable of generating inflammatory responses against different types of pathogen-derived MAMPs through the activation of TLR pathways [12,14,26]. Moreover, by using a transcriptomic approach, we showed that BME cells are capable of modifying their gene expression patterns in response to TLR4 activation. GO enrichment of DEGs revealed that the stimulation of BME cells with LPS induced the upregulation of genes involved in cell chemotaxis, response to bacterial molecules, and inflammatory responses.
Notably, there was a major difference in expression of inflammatory mediators in BME cells between LPS and S. aureus stimulations. The 12 h post LPS stimulation resulted in remarkable changes in expression of chemokines but it was not the case for infection with both S. aureus strains. Several in vitro studies based on exposure of BME cells to killed bacteria or purified bacterial MAMPs have shown that BME cells respond differently to E. coli and S. aureus [4,14,15,17,27,28,29]. This may explain the considerable variations in the initial recognition and subsequent pathogenesis of mastitis caused by these two pathogens. A previous study reported that severe inflammatory symptoms may be caused by the induction of IL1 and TNF-α expression in E. coli-challenged bovine mammary gland epithelial cells, while S. aureus quickly triggered an increased expression of IL6 [29]. Stimulation of BME cells with LPS for 24 h induced a marked increase in IL1, CXCL6, CXCL8, TNF-α and β-defensin expression while lipoteichoic acid (LTA, MAPM of Gram-positive bacteria) stimulation did not show any significant changes of these genes [15]. Activation of type IFN I pathway along with higher upregulation of chemokines (CXCL10, CCL2, CCL5, and CCL20) were reported in the case of LPS stimulation but not for S. aureus stimulation in BME cells [14]. However, the ability of S. aureus infection to evade the host’s immunity and to continue life-long mild inflammation [9] indicate that chemokine could be a late event in S. aureus/BME cells interaction, on the opposite of LPS/BME cell interaction. In this regard, it was reported the over expression of CCL5 and CXCL10 only after 24 h exposure of BME cells to heat killed S. aureus [30]. It should be noted that the present work has the limitation of testing only a single time point; therefore, a time series investigation of transcriptome shifting would provide better insight of S. aureus mediated inflammatory response in BME cells.
Our transcriptomic analysis revealed that the most notable change in LPS-challenged BME cells was the significant increase in the expression of IL1α, IL1β, CCL2, CXCL8, CXCL2 and CXCL3 (Figure 5A). Those inflammatory cytokines and chemokines have been associated with protection against pathogens infecting the bovine mammary gland since bacterial growth can be efficiently inhibited through the activity of leucocytes recruited from the blood into the mammary tissue [31,32]. Earlier studies by Waller [33], evaluating the effect of a polyclonal antibody to recombinant bovine IL-1β during LPS-induced inflammation in the bovine teat cistern, demonstrated the chemotactic effects of this cytokine and its protective role in the recruitment of neutrophils into the mammary gland. Studies in primary cell cultures of BME cells obtained from high or low mastitis-resistant heifers challenged with heat-inactivated E. coli showed higher expression levels of IL1β and IL8 in resistant animals when compared with susceptible cows [34]. Interestingly, it was reported that mammary epithelial cells constitutively secrete the chemokine CXCL3 and that normal milk contains active concentrations of this neutrophil chemoattractant protein [35]. The expression of this chemokine is increased in the course of mammary gland infections to act in synergy with other chemokines or cytokines to recruit leucocytes [36]. On the other hand, in addition of being a powerful neutrophil chemoattractant, CXCL2 shows direct antimicrobial properties against E. coli [37].
The results of the present work indicate that our in vitro BME cells system could be of value for the evaluation of immunomodulatory treatments aimed to improve the TLR4-mediated innate immune response in the bovine mammary gland and the protection against infections by Gram-negative pathogens. It should be noted that the inflammatory response in the mammary epithelium is beneficial for the clearance of pathogens and for the protection against infection; however, uncontrolled inflammation can be harmful and even fatal to the host [38,39]. The appropriate regulation of the production of proinflammatory mediators is highly important to maintain immune homeostasis and to prevent excessive inflammatory injury to the mammary gland. Thus, our in vitro BME cells system could be also of value for the selection and characterization of immunomodulatory treatments aimed to modulate the excessive inflammatory response that could become harmful for the bovine host.
The findings of this study also indicate that our in vitro BME cells system would be efficient for the selection and study of immunomodulatory treatments aimed to improve the innate immune response against Gram-positive pathogens in the bovine mammary gland. This limitation of our system could be related to the biology of MAPMs-PRRs interactions in the bovine mammary mucosa. It has been described that the production of chemotactic factors and the infiltration of immune cells in the mammary gland differs between Gram-positive and Gram-negative pathogens [15]. Several in vitro studies based on the exposure of BME cells to viable or non-viable bacteria as well as purified bacterial MAMPs have shown that mammary epithelial cells respond differently to Gram-positive (mainly Staphylococcus aureus) and Gram-negative (mainly E. coli) pathogens [14,15,16,27,28,29]. It was reported that in experimentally induced mastitis by E. coli challenge, increased concentration of chemokines and cytokines were found in milk samples. However, those inflammatory mediators were undetectable or were present in low concentrations during S. aureus mastitis [40]. In addition, the challenge of BME cells with LPS significantly increased the expression of TNF-α, IL1β and CXCL6 while the stimulation of these cells with Pam3CSK4 did not alter the expression of these genes [15]. Our results are in line with those previous findings since significant differential expression patterns in chemotactic factors were observed between LPS and Pam3CSK4-stimulated BME cells. We confirmed the improved expression of IL1β in LPS-stimulated BME cells and the lack of variation in this cytokine in Pam3CSK4-stimulated BME cells. Moreover, we extended those previous findings by demonstrating the upregulation of IL1α, CCL2, CXCL8, CXCL2 and CXCL3 in LPS-stimulated but not in Pam3CSK4-stimulated BME cells (Figure 5B).
It should be noted that the experiments in this work were limited to the evaluation of Pam3CSK4 that signals mainly through TLR2 [41]. Several other PRRs are involved in the recognition of Gram-positive pathogens and in the generation of innate immune responses. The combination of several MAMPs from Gram-positive pathogens acting on different PRRs may explain the findings of Lee et al. [42], who described higher levels of IL-8 in cells from mammary gland quarters infected with S. aureus relative to those infected with E. coli. Therefore, more detailed transcriptomic studies evaluating the effects of different Gram-positive pathogens-derived MAMPs (lipoproteins, peptidoglycans, teichoic acids) alone or in combinations as well as viable bacteria such as S. aureus in BME cells would be of importance to definitively establish the value of our system for the study of the bovine mammary gland inflammatory response triggered by Gram-positive bacteria.
In addition to MAMPs from pathogens, several other factors are able to trigger and amplify inflammatory responses in the mammary gland including cytokines, chemokines or reactive oxygen species. As mentioned before, CyPA is a protein secreted from mammalian cells in response to infection or inflammation [22,23,24]. It was reported that reactive oxygen species induce CyPA secretion from vascular smooth muscle cells resulting in the expression of the endothelial cell adhesion molecule and the aberrant infiltration of inflammatory cells in the aortic wall [43,44]. In addition, our recent study with bovine peripheral blood cells in vitro treated with recombinant bovine CyPA revealed that the granulocytes migrate towards this chemotactic molecule and the migration could be inhibited by pre-treatment with an anti-bovine CyPA antibodies (25). Then, we were also interested in evaluating the ability of recombinant bovine CyPA to modulate the inflammatory response of BME cells, not only to increase our knowledge about the immunobiology of this molecule, but also to determine if our system can be used to study the inflammatory response triggered by an eukaryotic molecule not included in the MAMPs category. Our results showed that the stimulation of BME cells with CyPA resulted in an increase of inflammatory cytokines and chemokine expression (Figure 5C). In particular, the CyPA was able to significantly increase the expression of IL1β, IL1α, CXCL8 and CXCL3 in BME cells. As mentioned earlier, IL-1β has chemotactic effects inducing the recruitment of neutrophils [33]. On the other hand, it is well established that CXCL8 induces chemotaxis in target cells, primarily neutrophils but also other granulocytes, causing them to migrate toward the site of infection while CXCL3 and CXCL2 controls migration and adhesion of neutrophils and monocytes [32]. Thus, the results obtained here indicate that CyPA, in addition to having its own chemotactic properties, can enhance this function by increasing the expression of chemokines by BME cells.

4. Materials and Methods

4.1. Cell Line

The bovine mammary epithelial (BME) cell line used in this study was originally developed by our group and the detail protocols for establishing the cell clone are available in the previous publication [21]. BME cells were isolated from mammary gland tissue taken from a 200-days pregnant Holstein cow as per the methods described and cryopreserved in liquid nitrogen. It has been demonstrated that this BME cell clone is functionally unaffected by cryopreservation and it is responsive to mitogen and lactogenic hormones [21].

4.2. Growth and Maintenance of BME Cell Line

For the present study, the cryopreserved BME cells were thawed and cultured in Dulbecco’s modified Eagle medium (DMEM, Gibco, Paisley, Scotland, UK) supplemented with 20% fetal calf serum (FCS; Sigma-Aldrich, Tokyo, Japan), 100 U/mL penicillin and 100 μg/mL streptomycin (Gibco™ 15140122, Life Technologies), transferrin (5 mg/mL) and sodium acetate (5 mM) as a growth medium. For the passage, BME cells were seeded 1 × 104/cm2 in the cell culture flask (Sumitomo Bakelite Co., Ltd., Tokyo, Japan), and incubated at 37 °C in a humidified atmosphere of 5% CO2. The cell culture medium was changed every 48 h and at least three consecutive passages were performed before the challenge experiments were conducted.

4.3. TLRs Copy Number Analysis

For absolute quantification of TLRs family mRNA, four pools of BME cells were subjected to prepare the cDNA standards for bovine TLR 1-10. Total RNA was isolated from BME cells and cDNAs were synthesized by reverse transcription from 5 μg of total RNA using olio (dT) primers and ThermoScript RNase H-reverse transcriptase as previously described [45]. The cDNA standards were purified and quantified spectrophotometrically. TLR and β-actin cDNAs were amplified by PCR using the primers listed in Supplementary Table S4. The purified PCR products were inserted into the vector pGEM-Teasy DNA (Promega, Madison, WI, USA). We confirmed the homology of each insert with dideoxy chain termination method using a DNA sequencer (4000 L; Li-Cor, Lincoln, NE, USA) and the SequiTherm EXCEL™ II DNA Sequencing Kit-LC (Epicentre Biotechnologies, USA). The copy number for each mRNA was determined using the following formula: Copy number (copies) = 6.02 × 1023 (copies/mol) × measurement (g)/MW (g/mol). MW = size in bp × 660 (g/mol/bp). Aliquots of standard cDNA containing 107–102 pg/μl were created for TLRs and used as assay standards.

4.4. Preparations of Stimulants

Two staphylococcal strains: S. aureus JE2 and S. aureus SA003 were grown in the antibiotic-free α-minimum essential medium. The bacteria were propagated up to mid-logarithmic phase (OD620 = 0.4) and washed twice with phosphate buffered saline (PBS). The determination of bacterial density was made by limiting dilution of washed bacteria using the SLGC bacterial calculator (Sansho Co. Ltd., Tokyo, Japan) according to the manufacturer’s instructions. Bacteria were resuspended in PBS, heat-killed by incubation at 95 °C for 20 minutes and stored at −20 °C until use. For mimicking the E. coli infection, we used commercially available lipopolysaccharide derived from Escherichia coli 55: B5(L6529-1mg, Sigma-Aldrich, Tokyo, Japan). In addition, we used commercially available Pam3CSK4 (Nacalai Tesque, Inc., Kyoto, Japan) for demonstrating TLR2 mediated inflammation.

4.5. Experimental Challenge to BME Cells and Sampling for Gene Expression Study

Five days-confluent BME cells (5 × 105 cells/well) were seeded in a twelve-well cell culture plate. After settling down, three pool of adherent cells were stimulated with either LPS (1.0 µg/mL) or heat-killed bacteria of S. aureus JE2 (5 × 107 cells/well) or S. aureus SA003 (5 × 107 cells/well) strains or kept unstimulated as control. The cultures were maintained at 37 °C in a humidified atmosphere of 5% CO2 for 12 h and then harvested for RNA extraction. The microarray-based gene expression analyses was performed in a single replicate for each treatment, followed by the RT-qPCR validation in three biological replicate samples as described below.
In a second set of experiments, BME cells (5 × 105 cells/well) were plated in a twelve-well cell culture plate and stimulated with either LPS (1.0 µg/mL) or Pam3CSK4 (100 ng/mL) or recombinant bovine CyPA (10, 100 or 1000 ng/mL) or kept unstimulated as control. The cultures were maintained at 37 °C in a humidified atmosphere of 5% CO2 for 3, 6, 12 or 24 h and then harvested at every time point for RNA extraction and RT-qPCR analysis as described below.

4.6. RNA Isolation, Labeling and Microarray Hybridization

Total RNA was extracted from the treated and untreated control cells using PureLink RNA Mini Kit (Life Technology INC., USA) along with on-column DNase treatment. RNA integrity, quality and quantity were evaluated with microcapillary electrophoresis (2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) using the Agilent RNA 6000 Nano Kit (Agilent Technologies, Santa Clara, CA, USA). Samples with RNA integrity number (RIN) greater than 9 were used for microarray study.
Microarray analyses were performed using a one-color microarray-based Gene Expression Analysis Kit (Agilent, USA) and Bovine (V2) Gene Expression Microarray 4 × 44K oligonucleotide slides (Agilent, USA). Cyanine-3 (Cy3) labeled cRNA was prepared from 200 ng RNA using the One-Color Low Input Quick Amp Labeling Kit (Agilent) according to the manufacturer’s instructions. Dye incorporation and cRNA yield were checked with the NanoDrop ND-1000 Spectrophotometer. The Gene Expression Hybridization kit (Agilent Technologies) was used for hybridization. In brief, 1650 ng of Cy3-labelled cRNA (specific activity >10.0 pmol Cy3/ug cRNA) was fragmented at 60 °C for 30 minutes in a reaction volume of 250 mL containing 1x Agilent fragmentation buffer and 2x Agilent blocking agent following the manufacturer’s instructions. On completion of the fragmentation reaction, 250 mL of 2x Agilent hybridization buffer was added to the fragmentation mixture and hybridized to Agilent Bovine (V2) Gene Expression Microarray, 4 × 44K (G2519F) for 17 hours at 65 °C in a rotating Agilent hybridization oven. After hybridization, microarrays were washed 1 minute at room temperature with GE Wash Buffer 1 (Agilent) and 1 minute with 37 °C GE Wash buffer 2 (Agilent), then dried immediately by brief centrifugation.

4.7. Statistical Analysis of Microarray Data

The normalization and differential expression analysis of microarray data were performed with GeneSpring GX software (v13.1, Agilent Technologies, USA). The log2 transformed expression values of probes were normalized based on 75 percentile shifts. For identification of differentially expressed genes (DEGs) following each stimulation, we compared gene expression profiles between control and each treatment: LPS, S. aureus JE2 and S. aureus SA003. In order to determine the differential expression of genes, an unpaired t-test was performed between untreated control and stimulated samples. Benjamini and Hochberg (B-H) adjustment method was applied for multiple test correction. Significant differentially expressed genes were selected on the basis of two criteria: an adjusted p-value (FDR, false discover rate) of less than 0.05 and a cutoff in fold change of at least 1.5. The annotation and conversion of DEGs into human orthologous ensembl gene ID were performed by using the bioDBnet tool [46] and BovineMine [47], which were subjected for downstream functional analysis. The MIMAE (minimum information about a microarray experiment) standard raw microarray dataset have been submitted to the NCBI-GEO database under the access number GSE139612.

4.8. Gene Ontology Analyses

For biological interpretation of differential gene expressions, gene ontology (GO) enrichment analysis was performed using the InnateDB online tool (v5.4, [48]). The ensemble gene ID were uploaded to the InnateDB web portal and performed GO analysis. For each case, the over representation analysis (ORA) was performed using hypergeometric algorithm with Benjamini-Hochberg (B-H) multiple test correction method.

4.9. RT-qPCR Validation

Quantitative two-step real-time PCR was performed to confirm the microarray results and to evaluate the effect of Pam3CSK4 and CyPA by quantifying the mRNA expression levels of six selected genes in BME cells. Primer sequences are presented in Supplementary Table S4. Total RNA was isolated from each sample using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) followed by treated with gDNA Wipeout Buffer (Qiagen, Tokyo, Japan). All cDNAs were synthesized using the Quantitect reverse transcription Kit (Qiagen, Tokyo, Japan) according to the manufacturer’s recommendations.
The RT-qPCR was performed using 7300 real-time PCR system (Applied Biosystems, Warrington, UK) using the TaqMan® gene expression assay kit (Life Technologies) and TaqMan® Universal Master Mix II, with UNG (Applied Biosystems, Warrington, UK). The PCR cycling conditions were 2 min at 50 °C, followed by 10 min at 95 °C, and then 40 cycles of 15 s at 95 °C, 1 min at 60 °C. The reaction mixtures contained 2.5 μL of sample cDNA, 1 μL gene expression assay and 10 μL TaqMan® Universal Master mix II, with UNG, and 6.5 μL distilled water. According to the minimum information for publication of quantitative real-time PCR experiments guidelines, β-actin was used as a house-keeping gene because of its high stability across various bovine tissues [49]. Relative index was calculated as the ratio of target mRNA expression to β-actin. Then, raw data were transferred from the mean Ct values of replicated samples to copy number of the established standard curve. The normality of the log2 transformed data distribution were checked by Kolmogorov-Smirnov test. Comparisons between mean values were carried out using one-way ANOVA and Fisher’s least significant difference test. For every case, p < 0.05 was considered statistically significant.

5. Conclusions

We herein presented the transcriptome profile associated with the inflammatory response of BME cells following in vitro stimulation with LPS or heat-killed S. aureus JE2 or S. aureusSA003. The challenge of BME cells with LPS resulted in the activation of the TLR4 signaling pathway that conducted to a strong inflammatory response compared with the responses induced by S. aureus strains or other inflammatory factors such as Pam3CSK4 or CyPA. The genes with marked variations in response to TLR4 activation in BME cells identified in the present study including IL1β, IL1α, CCL2, CXCL8, CXCL2 and CXCL3 could be used as potential biomarkers for the selection of effective immunomodulatory treatments for improving the outcome of mastitis through the use of the BME cells in vitro evaluation system. Precise studies of the transcriptional changes induced by the activation of PRRs by pathogens in BME cells are necessary to provide the scientific basis for the generation of immunointervention tools that improve the outcome of infections and protect against inflammatory damage in the mammary gland of the bovine host. The results of this work are an advance in this regard.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-0817/9/3/200/s1. There are four additional files provided within a zip folder. Table S1: List of DEGs after LPS stimulation to BME cells. Table S2: List of DEGs after Staphylococcus aureus JE2 stimulation to BME cells Table S3: List of DEGs after Staphylococcus aureusSA003 stimulation to BME cells. Table S4: The primer sequences used for qPCR study.

Author Contributions

Conceptualization, H.A., H.K., M.A.I. and M.T.; methodology, M.T., K.F., R.K. and M.A.I.; software, M.A.I. and L.A.; validation, M.T., R.K.; formal analysis, M.A.I. and J.V.; investigation, J.V. and H.K.; resources, H.A.; data curation, M.A.I. and M.T.; writing—original draft preparation, M.A.I.; writing—review and editing, J.V., Y.S., W.I.-O., T.N., V.R., and W.v.E.; visualization, M.A.I., L.A., and J.V.; supervision, H.K.; project administration, H.K.; funding acquisition, H.A. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the project of NARO Bio-oriented Technology Research Advancement Institution (Research Program on development of innovation technology, No. 29006A) to HA, and a Grant-in-Aid for Scientific Research (A) (19H00965), (B) (16H05019), Challenging Exploratory Research (26660216, 16K15028), Open Partnership Joint Projects of JSPS Bilateral Joint Research Projects from the Japan Society for the Promotion of Science (JSPS) to HK. MAI was supported by JSPS (Postdoctoral Fellowship for Foreign Researchers, Program No. 18F18081). This work was also supported by JSPS Core-to-Core Program, A. Advanced Research Networks entitled Establishment of international agricultural immunology research-core for a quantum improvement in food safety.

Conflicts of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationship that could be constructed as a potential conflict of interest.

References

  1. Seegers, H.; Fourichon, C.; Beaudeau, F. Production effects related to mastitis and mastitis economics in dairy cattle herds. Vet. Res. 2003, 34, 475–491. [Google Scholar] [CrossRef] [Green Version]
  2. Dos Reis, C.B.M.; Barreiro, J.R.; Mestieri, L.; de Felício Porcionato, M.A.; dos Santos, M.V. Effect of somatic cell count and mastitis pathogens on milk composition in Gyr cows. BMC Vet. Res. 2013, 9, 67. [Google Scholar]
  3. Pumipuntu, N.; Kulpeanprasit, S.; Santajit, S.; Tunyong, W.; Kong-Ngoen, T.; Hinthong, W.; Indrawattana, N. Screening method for Staphylococcus aureus identification in subclinical bovine mastitis from dairy farms. Vet. World 2017, 10, 721–726. [Google Scholar] [CrossRef] [Green Version]
  4. Jensen, K.; Günther, J.; Talbot, R.; Petzl, W.; Zerbe, H.; Schuberth, H.J.; Seyfert, H.M.; Glass, E.J. Escherichia coli- and Staphylococcus aureus-induced mastitis differentially modulate transcriptional responses in neighbouring uninfected bovine mammary gland quarters. BMC Genom. 2013, 14, 36. [Google Scholar] [CrossRef] [Green Version]
  5. Burvenich, C.; Van Merrid, V.; Mehrzad, J.; ez-Fraile, A.; Duchateau, L. Severity of E. coli mastitis is mainly determined by cow factors. Vet. Res. 2003, 34, 521–564. [Google Scholar] [CrossRef] [Green Version]
  6. Schukken, Y.H.; Günther, J.; Fitzpatrick, J.; Fontaine, M.C.; Goetze, L.; Holst, O.; Leigh, J.; Petzl, W.; Schuberth, H.-J.; Sipka, A.; et al. Host-response patterns of intramammary infections in dairy cows. Vet. Immunol. Immunopathol. 2011, 144, 270–289. [Google Scholar] [CrossRef]
  7. Mitterhuemer, S.; Petzl, W.; Krebs, S.; Mehne, D.; Klanner, A.; Wolf, E.; Zerbe, H.; Blum, H. Escherichia coli infection induces distinct local and systemic transcriptome responses in the mammary gland. BMC Genom. 2010, 11, 138. [Google Scholar] [CrossRef] [Green Version]
  8. Bannerman, D.D. Pathogen-dependent induction of cytokines and other soluble inflammatory mediators during intramammary infection of dairy cows. J. Anim. Sci. 2009, 87, 10–25. [Google Scholar] [CrossRef] [Green Version]
  9. Gresham, H.D.; Lowrance, J.H.; Caver, T.E.; Wilson, B.S.; Cheung, A.L.; Lindberg, F.P. Survival of Staphylococcus aureus inside neutrophils contributes to infection. J. Immunol. 2000, 164, 3713–3722. [Google Scholar] [CrossRef] [Green Version]
  10. Paape, M.J.; Mehrzad, J.; Zhao, X.; Detilleux, J.; Burvenich, C. Defense of the bovine mammary gland by polymorphonuclear neutrophil leukocytes. J. Mammary Gland Biol. 2002, 7, 109–121. [Google Scholar] [CrossRef]
  11. Paape, M.J.; Bannerman, D.D.; Zhao, X.; Lee, J.W. The bovine neutrophil: Structure and function in blood and milk. Vet. Res. 2003, 34, 597–627. [Google Scholar] [PubMed] [Green Version]
  12. Kawai, T.; Akira, S. The role of pattern-recognition receptors in innate immunity: Update on Toll-like receptors. Nat. Immunol. 2010, 11, 373–384. [Google Scholar] [CrossRef] [PubMed]
  13. Aitken, S.L.; Corl, C.M.; Sordillo, L.M. Immunopathology of mastitis: Insights into disease recognition and resolution. J. Mammary Gland Biol. Neoplasia 2011, 16, 291–304. [Google Scholar] [CrossRef]
  14. Gilbert, F.B.; Cunha, P.; Jensen, K.; Glass, E.J.; Foucras, G.; Robert-Granié, C.; Rupp, R.; Rainard, P. Differential response of bovine mammary epithelial cells to Staphylococcus aureus or Escherichia coli agonists of the innate immune system. Vet. Res. 2013, 44, 40. [Google Scholar] [CrossRef] [Green Version]
  15. Strandberg, Y.; Gray, C.; Vuocolo, T.; Donaldson, L.; Broadway, M.; Tellam, R. Lipopolysaccharide and lipoteichoic acid induce different innate immune responses in bovine mammary epithelial cells. Cytokine 2005, 31, 72–86. [Google Scholar] [CrossRef]
  16. Günther, J.; Koczan, D.; Yang, W.; Nurnberg, G.; Repsilber, D.; Schuberth, H.J.; Park, Z.; Maqbool, N.; Molenaar, A.; Seyfert, H.M. Assessment of the immune capacity of mammary epithelial cells: Comparison with mammary tissue after challenge with Escherichia coli. Vet. Res. 2009, 40, 31. [Google Scholar] [CrossRef] [Green Version]
  17. Griesbeck-Zilch, B.; Meyer, H.H.; Kuhn, C.H.; Schwerin, M.; Wellnitz, O. Staphylococcus aureus and Escherichia coli cause deviating expression profiles of cytokines and lactoferrin messenger ribonucleic acid in mammary epithelial cells. J. Dairy Sci. 2008, 91, 2215–2224. [Google Scholar] [CrossRef] [Green Version]
  18. Genini, S.; Badaoui, B.; Sclep, G.; Bishop, S.C.; Waddington, D.; van der Pinard, L.M.H.; Klopp, C.; Cabau, C.; Seyfert, H.M.; Petzl, W.; et al. Strengthening insights into host responses to mastitis infection in ruminants by combining heterogeneous microarray data sources. BMC Genom. 2011, 12, 225. [Google Scholar] [CrossRef] [Green Version]
  19. Wang, X.; Xiu, L.; Hu, Q.; Cui, X.; Liu, B.; Tao, L.; Wang, T.; Wu, J.; Chen, Y.; Chen, Y. Deep Sequencing-Based Transcriptional Analysis of Bovine Mammary Epithelial Cells Gene Expression in Response to In Vitro Infection with Staphylococcus aureus Stains. PLoS ONE 2013, 8, e82117. [Google Scholar] [CrossRef]
  20. McGrath, M.F. A novel system for mammary epithelial cell culture. J. Dairy Sci. 1987, 70, 1967–1980. [Google Scholar] [CrossRef]
  21. Rose, T.M.; Aso, H.; Yonekura, S.; Komatsu, T.; Hagino, A.; Ozutsumi, K.; Obara, Y. In Vitro differentiation of a cloned bovine mammary epithelial cell. J. Dairy Res. 2002, 69, 345–355. [Google Scholar] [CrossRef] [PubMed]
  22. Yurchenko, V.; Constant, S.; Eisenmesser, E.; Bukrinsky, M. Cyclophilin-CD147 interactions: A new target for anti-inflammatory therapeutics. Clin. Exp. Immunol. 2010, 160, 305–317. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Nigro, P.; Pompilio, G.; Capogrossi, M.C. Cyclophilin A: A key player for human disease. Cell Death Dis. 2013, 4, e888. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Handschumacher, R.E.; Harding, M.W.; Rice, J.; Drugge, R.J.; Speicher, D.W. Cyclophilin: A specific cytosolic binding protein for cyclosporin A. Science 1984, 226, 544–547. [Google Scholar] [CrossRef]
  25. Takanashi, S.; Nochi, T.; Abe, M.; Itaya, N.; Urakawa, M.; Sato, K.; Zhuang, T.; Umemura, S.; Hayashi, T.; Kiku, Y.; et al. Extracellular cyclophilin A possesses chemotaxic activity in cattle. Vet. Res. 2015, 46, 80. [Google Scholar] [CrossRef] [Green Version]
  26. Zhang, W.Y.; Wang, H.; Qi, S.; Wang, X.; Li, X.; Zhou, K.; Zhang, Y.; Gao, M.Q. CYP1A1 Relieves Lipopolysaccharide-Induced Inflammatory Responses in Bovine Mammary Epithelial Cells. Mediat. Inflamm. 2018, 2018, 4093285. [Google Scholar] [CrossRef]
  27. McClenahan, D.J.; Sotos, J.P.; Czuprynski, C.J. Cytokine response of bovine mammary gland epithelial cells to Escherichia coli, coliform culture filtrate, or lipopolysaccharide. Am. J. Vet. Res. 2005, 66, 1590–1597. [Google Scholar] [CrossRef]
  28. Petzl, W.; Zerbe, H.; Günther, J.; Yang, W.; Seyfert, H.M.; Nürnberg, G.; Schuberth, H.J. Escherichia coli, but not Staphylococcus aureus, triggers an early increased expression of factors contributing to the innate immune defense in the udder of the cow. Vet. Res. 2008, 39, 18. [Google Scholar] [CrossRef] [Green Version]
  29. Günther, J.; Esch, K.; Poschadel, N.; Petzl, W.; Zerbe, H.; Mitterhuemer, S.; Blum, H.; Seyfert, H.M. Comparative kinetics of Escherichia coli-and Staphylococcus aureus-specific activation of key immune pathways in mammary epithelial cells demonstrates that S. aureus elicits a delayed response dominated by interleukin-6 (IL-6) but not by IL-1A or tumor necrosis factor alpha. Infect. Immun. 2011, 79, 695–707. [Google Scholar]
  30. Brand, B.; Hartmann, A.; Repsilber, D.; Griesbeck-Zilch, B.; Wellnitz, O.; Kuhn, C.; Ponsuksili, S.; Meyer, H.H.; Schwerin, M. Comparative expression profiling of E. coli and S. aureus inoculated primary mammary gland cells sampled from cows with different genetic predispositions for somatic cell score. Genet. Sel. Evol. 2011, 43, 24. [Google Scholar] [CrossRef] [Green Version]
  31. Rainard, P.; Riollet, C. Innate immunity of the bovine mammary gland. Vet. Res. 2006, 37, 369–400. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Wellnitz, O.; Bruckmaier, R.M. The innate immune response of the bovine mammary gland to bacterial infection. Vet. J. 2012, 192, 148–152. [Google Scholar] [CrossRef] [PubMed]
  33. Waller, K.P. Modulation of endotoxin-induced inflammation in the bovine teat using antagonists/inhibitors to leukotrienes, platelet activating factor and interleukin 1 beta. Vet. Immunol. Immunopathol. 1997, 57, 239–251. [Google Scholar] [CrossRef]
  34. Griesbeck-Zilch, B.; Osman, M.; Kühn, C.; Schwerin, M.; Bruckmaier, R.H.; Pfaffl, M.W.; Hammerle-Fickinger, A.; Meyer, H.H.; Wellnitz, O. Analysis of key molecules of the innate immune system in mammary epithelial cells isolated from marker-assisted and conventionally selected cattle. J. Dairy Sci. 2009, 92, 4621–4633. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Reinard, P.; Riollet, C.; Berthon, P.; Cunha, P.; Fromageau, A.; Rossingnol, C.; Gilbert, F.B. The chemokine CXCL3 is responsible for the constitutive chemotactic activity of bovine milk for neutrophils. Mol. Immunol. 2008, 45, 4020–4027. [Google Scholar] [CrossRef] [PubMed]
  36. Struyf, S.; Gouwy, M.; Dillen, C.; Proost, P.; Opdenakker, G.; Van Damme, J. Chemokines synergize in the recruitment of circulating neutrophils into inflamed tissue. Eur. J. Immunol. 2005, 35, 1583–1591. [Google Scholar] [CrossRef]
  37. Yang, D.; Chen, Q.; Hoover, D.M.; Staley, P.; Tucker, K.D.; Lubkowski, J.; Oppenheim, J.J. Many chemokines including CCL20/MIP-3alpha display antimicrobial activity. J. Leukoc. Biol. 2003, 74, 448–455. [Google Scholar] [CrossRef]
  38. Liew, F.Y.; Xu, D.; Brint, E.K.; O’Neill, L.A. Negative regulation of toll-like receptor-mediated immune responses. Nat. Rev. Immunol. 2005, 5, 446–458. [Google Scholar] [CrossRef]
  39. Liu, M.; Song, S.; Li, H.; Jiang, X.; Yin, P.; Wan, C.; Liu, X.; Liu, F.; Xu, J. The protective effect of caffeic acid against inflammation injury of primary bovine mammary epithelial cells induced by lipopolysaccharide. J. Dairy Sci. 2014, 97, 2856–2865. [Google Scholar] [CrossRef] [Green Version]
  40. Bannerman, D.D.; Paape, M.J.; Lee, J.W.; Zhao, X.; Hope, J.C.; Rainard, P. Escherichia coli and Staphylococcus aureus elicit differential innate immune responses following intramammary infection. Clin. Diagn. Lab. Immunol. 2004, 11, 463–472. [Google Scholar] [CrossRef] [Green Version]
  41. Schröder, N.W.; Morath, S.; Alexander, C.; Hamann, L.; Hartung, T.; Zahringer, U.; Gobel, U.B.; Weber, J.R.; Schumann, R.R. Lipoteichoic acid (Pam3CSK4) of Streptococcus pneumoniae and Staphylococcus aureus activates immune cells via Toll-like receptor (TLR)-2, lipopolysaccharide binding protein (LBP), and CD14, whereas TLR-4 and MD-2 are not involved. J. Biol. Chem. 2003, 278, 15587–15594. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Lee, J.W.; Bannerman, D.D.; Paape, M.J.; Huang, M.K.; Zhao, X. Characterization of cytokine expression in milk somatic cells during intramammary infections with Escherichia coli or Staphylococcus aureus by real-time PCR. Vet. Res. 2006, 37, 219–229. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Satoh, K.; Nigro, P.; Berk, B.C. Oxidative stress and vascular smooth muscle cell growth: A mechanistic linkage by cyclophilin A. Antioxid. Redox Signal. 2010, 12, 675–682. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Satoh, K.; Nigro, P.; Matoba, T.; O’Dell, M.R.; Cui, Z.; Shi, X.; Mohan, A.; Yan, C.; Abe, J.; Illig, K.A.; et al. Cyclophilin A enhances vascular oxidative stress and the development of angiotensin II-induced aortic aneurysms. Nat. Med. 2009, 15, 649–656. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Tohno, M.; Shimosato, T.; Kawai, Y.; Ikegami, S.; Taketomo, N.; Saito, T.; Kitazawa, H. Advanced molecular immunoassay system for immunobiotic lactic acid bacteria using a transfectant of Toll-like receptor 2. Anim. Sci. J. 2007, 78, 195–205. [Google Scholar] [CrossRef]
  46. Mudunuri, U.; Che, A.; Yi, M.; Stephens, R.M. bioDBnet: The biological database network. Bioinformatics 2009, 25, 555–556. [Google Scholar] [CrossRef] [Green Version]
  47. Elsik, C.G.; Unni, D.R.; Diesh, C.M.; Tayal, A.; Emery, M.L.; Nguyen, H.N.; Hagen, D.E. Bovine Genome Database: New tools for gleaning function from the Bos taurus genome. Nucleic Acids Res. 2016, 44, D834–D839. [Google Scholar] [CrossRef] [Green Version]
  48. Breuer, K.; Foroushani, A.K.; Laird, M.R.; Chen, C.; Sribnaia, A.; Lo, R.; Winsor, G.L.; Hancock, R.E.; Brinkman, F.S.; Lynn, D.J. InnateDB: Systems biology of innate immunity and beyond--recent updates and continuing curation. Nucleic Acids Res. 2013, 41, 24. [Google Scholar] [CrossRef]
  49. Bustin, S.A.; Benes, V.; Garson, J.A.; Hellemans, J.; Huggett, J.; Kubista, M.; Mueller, R.; Nolan, T.; Pfaffl, M.W.; Shipley, G.L.; et al. The MIQE guidelines: Minimum information for publication of quantitative real-time PCR experiments. Clin. Chem. 2009, 55, 611–622. [Google Scholar] [CrossRef] [Green Version]
Figure 1. RT-qPCR-based absolute quantification of Toll-like receptors (TLRs) mRNA expression levels in a 25 ng of cDNA prepared from untreated bovine mammary epithelial (BME) cells. Results presented are as mean ± SD of four independent experiments. Different letters (a, b, c) indicates statistically significant (p < 0.05) difference between them.
Figure 1. RT-qPCR-based absolute quantification of Toll-like receptors (TLRs) mRNA expression levels in a 25 ng of cDNA prepared from untreated bovine mammary epithelial (BME) cells. Results presented are as mean ± SD of four independent experiments. Different letters (a, b, c) indicates statistically significant (p < 0.05) difference between them.
Pathogens 09 00200 g001
Figure 2. Overview of the differentially expressed genes (DEGs) in BME cells. A. Number of upregulated genes, B. Number of downregulated genes. C. Hierarchical clustering of DEGs.
Figure 2. Overview of the differentially expressed genes (DEGs) in BME cells. A. Number of upregulated genes, B. Number of downregulated genes. C. Hierarchical clustering of DEGs.
Pathogens 09 00200 g002
Figure 3. Expression patterns of inflammatory cytokines and chemokines genes in bovine mammary epithelial (BME) cells after lipopolysaccharide (LPS) and Pam3CSK4 stimulation. The expression of interleukin 1 alpha (IL1α), IL1β, C-C motif chemokine ligand 2 (CCL2), C-X-C motif chemokine ligand 2 (CXCL2), CXCL3 and CXCL8 were determined in BME cells after 3, 6, 12 and 24 hours of LPS and Pam3CSK4 stimulation by qPCR. BME cells with no LPS treatment were used as basal controls (hour 0). The results represent data (mean ± SD) from three independent experiments. Significantly different compared to basal control: *, p < 0.05, **p < 0.01.
Figure 3. Expression patterns of inflammatory cytokines and chemokines genes in bovine mammary epithelial (BME) cells after lipopolysaccharide (LPS) and Pam3CSK4 stimulation. The expression of interleukin 1 alpha (IL1α), IL1β, C-C motif chemokine ligand 2 (CCL2), C-X-C motif chemokine ligand 2 (CXCL2), CXCL3 and CXCL8 were determined in BME cells after 3, 6, 12 and 24 hours of LPS and Pam3CSK4 stimulation by qPCR. BME cells with no LPS treatment were used as basal controls (hour 0). The results represent data (mean ± SD) from three independent experiments. Significantly different compared to basal control: *, p < 0.05, **p < 0.01.
Pathogens 09 00200 g003
Figure 4. Expression patterns of inflammatory cytokines and chemokines in bovine mammary epithelial (BME) cells after recombinant bovine cyclophilin A (CyPA) stimulation. The expression of interleukin 1 alpha (IL1α), IL1β, C-C motif chemokine ligand 2 (CCL2), C-X-C motif chemokine ligand 2 (CXCL2), CXCL3 and CXCL8 were determined by qPCR in BME cells after 3, 6, 12 and 24 hours of stimulation with different concentrations of CyPA (10, 100 or 1000 ng/mL). BME cells with no CyPA treatment were used as controls (0 ng/mL). The results represent data (mean ± SD) from three independent experiments. Significantly different compared to control in the same time point: *, p < 0.05.
Figure 4. Expression patterns of inflammatory cytokines and chemokines in bovine mammary epithelial (BME) cells after recombinant bovine cyclophilin A (CyPA) stimulation. The expression of interleukin 1 alpha (IL1α), IL1β, C-C motif chemokine ligand 2 (CCL2), C-X-C motif chemokine ligand 2 (CXCL2), CXCL3 and CXCL8 were determined by qPCR in BME cells after 3, 6, 12 and 24 hours of stimulation with different concentrations of CyPA (10, 100 or 1000 ng/mL). BME cells with no CyPA treatment were used as controls (0 ng/mL). The results represent data (mean ± SD) from three independent experiments. Significantly different compared to control in the same time point: *, p < 0.05.
Pathogens 09 00200 g004
Figure 5. Comparative profiles of selected inflammatory gene expressions in BME cells after LPS (A), Pam3CSK4 (B) and CyPA (C) stimulation. The mRNA of IL1α, IL1β, CCl2, CXCL8, CXCL2, CXCL3, were upregulated after LPS stimulation. Among them only IL-1α, TLR2 and NOD2 were upregulated but others have no changes after Pam3CSK4 stimulation. While CXCL3, CXCL8, IL1α, and IL1β were upregulated and CXCL2 and CCL2 had no change after CyPA stimulation. FI, Fold increase.
Figure 5. Comparative profiles of selected inflammatory gene expressions in BME cells after LPS (A), Pam3CSK4 (B) and CyPA (C) stimulation. The mRNA of IL1α, IL1β, CCl2, CXCL8, CXCL2, CXCL3, were upregulated after LPS stimulation. Among them only IL-1α, TLR2 and NOD2 were upregulated but others have no changes after Pam3CSK4 stimulation. While CXCL3, CXCL8, IL1α, and IL1β were upregulated and CXCL2 and CCL2 had no change after CyPA stimulation. FI, Fold increase.
Pathogens 09 00200 g005
Table 1. The number of differentially expressed genes (DEGs) in BME cells after stimulation with lipopolysaccharide (LPS) or heat-killed S. aureus strains. Statistically significant DEGs were detected based on the criteria of p-value < 0.05, false discovery rate (FDR) of 0.1 and fold change (FC) 1.5, as compared to the untreated control.
Table 1. The number of differentially expressed genes (DEGs) in BME cells after stimulation with lipopolysaccharide (LPS) or heat-killed S. aureus strains. Statistically significant DEGs were detected based on the criteria of p-value < 0.05, false discovery rate (FDR) of 0.1 and fold change (FC) 1.5, as compared to the untreated control.
Control vs. LPSControl vs. S. aureus JE2Control vs. S. aureus SA003
Upregulated222227226
Downregulated207238294
Total 447465520
Table 2. Gene ontology (GO) terms enriched after stimulation with LPS or heat-killed S. aureus strains to BME cells.
Table 2. Gene ontology (GO) terms enriched after stimulation with LPS or heat-killed S. aureus strains to BME cells.
GO Terms and IDsGenes Involved p-Value
LPS stimulation:
Cell chemotaxis [GO:0060326]AGTR1, CCL2, CXCL2, CXCL3, CXCL5, CXCL8, CXCL90.0001
Inflammatory response [GO:0006954]CCL20, CCL2, CXCL2, CXCL8, CXCL9, GGT5, IL1α, IL1β0.002
Positive regulation of nitric-oxide synthetase production [GO:0051770]CCL20, CCL2, NOD20.0002
Response molecules of bacterial origin [GO:0002237]CXCL2, CXCL80.003
Negative regulation of extrinsic apoptotic signaling pathway [GO:2001240]IL1α, IL1β, IL17, UNC5B0.006
S. aureus JE2 stimulation:
Cell differentiation [GO:0030154]CCDC135, DLL1, HEMGN, MGP, NEUROD6, NKX2-2, PDX1, PIWIL3, SEMA4G, TCF230.001
Cell adhesion [GO:0007155]ABL1, CLDN11, ITGBL1, MAG, MPDZ, NCAN, PRKCE0.02
Inflammatory response [GO:0006954]CCL26, CNR2, CRP, CXCR2, IL220.04
Blood coagulation [GO:0007596]ABL1, EFEMP2, ESAM, KIF2B, MAG, PRKCE, PRKG2, SERPINB20.08
Chemotaxis [GO:0006935]CCL26, CXCR2, CXCR50.19
S. aureus SA003 stimulation:
Innate immune response [GO:0045087]AHSG, APOA1, C8A, F11, FREM1, KLRD1, MAFB, MAPK10, OTUB1, PRKCE, SFTPA10.01
Cell surface receptor signaling pathway [GO:0007166]GPR97, KLRD1, LEPR, LIFR0.03
Blood coagulation [GO:0007596]APOA1, ESAM, F11, KIF3C, PRKCE, TRPC30.03
Signal transduction [GO:0007165]CAPN3, ECM1, GDF3, MAPK10, NPAS2, PDE4B, PRKCE, SYNGAP10.04
Positive regulation of NF-kB signaling [GO:0043123]ECM1, HTR2B, PRKCE0.06
Table 3. Expression patterns of selected inflammation related genes in BME cells after stimulation with LPS or heat-killed S. aureus strains.
Table 3. Expression patterns of selected inflammation related genes in BME cells after stimulation with LPS or heat-killed S. aureus strains.
Gene SymbolGene Name LPS S. aureus JE2S. aureus SA003
IL1αInterleukin 1 Alpha1.463N.MN.M
IL1βInterleukin 1 Beta2.94N.MN.M
CCL2C-C Motif Chemokine Ligand 2 (also called monocyte chemoattractant protein -1, MCP-1)1.682N.MN.M
CX3CL1C-X3-C Motif Chemokine Ligand 11.146N.MN.M
CXCL2C-X-C Motif Chemokine Ligand 22.953N.MN.M
CXCL3C-X-C Motif Chemokine Ligand 31.194N.MN.M
CXCL5C-X-C Motif Chemokine Ligand 52.007N.MN.M
CXCL8C-X-C Motif Chemokine Ligand 8 (also called Interleukin 8, IL8)3.013N.MN.M
CXCL9C-X-C Motif Chemokine Ligand 94.379N.MN.M
C6Complement component 61.313N.MN.M
CRPC-Reactive ProteinN.M3.012N.M
M-SAA3.2Mammary Serum Amyloid A3.23.413N.MN.M
SAA3Serum Amyloid A33.869N.MN.M
PPIAPeptidylprolyl Isomerase A (also called Cyclophilin A, CyPA)N.MN.MN.M
Values shown as the log2 fold-change of corresponding genes, N.M., Not modified (<−1.5 to > 1.5).

Share and Cite

MDPI and ACS Style

Islam, M.A.; Takagi, M.; Fukuyama, K.; Komatsu, R.; Albarracin, L.; Nochi, T.; Suda, Y.; Ikeda-Ohtsubo, W.; Rutten, V.; Eden, W.v.; et al. Transcriptome Analysis of The Inflammatory Responses of Bovine Mammary Epithelial Cells: Exploring Immunomodulatory Target Genes for Bovine Mastitis. Pathogens 2020, 9, 200. https://doi.org/10.3390/pathogens9030200

AMA Style

Islam MA, Takagi M, Fukuyama K, Komatsu R, Albarracin L, Nochi T, Suda Y, Ikeda-Ohtsubo W, Rutten V, Eden Wv, et al. Transcriptome Analysis of The Inflammatory Responses of Bovine Mammary Epithelial Cells: Exploring Immunomodulatory Target Genes for Bovine Mastitis. Pathogens. 2020; 9(3):200. https://doi.org/10.3390/pathogens9030200

Chicago/Turabian Style

Islam, Md. Aminul, Michihiro Takagi, Kohtaro Fukuyama, Ryoya Komatsu, Leonardo Albarracin, Tomonori Nochi, Yoshihito Suda, Wakako Ikeda-Ohtsubo, Victor Rutten, Willem van Eden, and et al. 2020. "Transcriptome Analysis of The Inflammatory Responses of Bovine Mammary Epithelial Cells: Exploring Immunomodulatory Target Genes for Bovine Mastitis" Pathogens 9, no. 3: 200. https://doi.org/10.3390/pathogens9030200

APA Style

Islam, M. A., Takagi, M., Fukuyama, K., Komatsu, R., Albarracin, L., Nochi, T., Suda, Y., Ikeda-Ohtsubo, W., Rutten, V., Eden, W. v., Villena, J., Aso, H., & Kitazawa, H. (2020). Transcriptome Analysis of The Inflammatory Responses of Bovine Mammary Epithelial Cells: Exploring Immunomodulatory Target Genes for Bovine Mastitis. Pathogens, 9(3), 200. https://doi.org/10.3390/pathogens9030200

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop