MicroRNA expression profiles of bovine monocyte-derived macrophages infected in vitro with two strains of Streptococcus agalactiae
- Anna Monika Lewandowska-Sabat1Email author,
- Silje Furre Hansen2,
- Trygve Roger Solberg3,
- Olav Østerås4,
- Bjørg Heringstad3, 5,
- Preben Boysen†2 and
- Ingrid Olsaker†1
Received: 16 May 2017
Accepted: 9 March 2018
Published: 10 April 2018
Abstract
Background
MicroRNAs (miRNAs) are short, non-coding RNAs that regulate gene expression at the post-transcriptional level and play a key role in the control of innate and adaptive immune responses. For a subclinical infection such as bovine streptococcal mastitis, early detection is a great challenge, and miRNA profiling could potentially assist in the diagnosis and contribute to the understanding of the pathogenicity and defense mechanisms. We have examined the miRNA repertoire and the transcript level of six key immune genes [tumor necrosis factor alpha (TNFα), interleukin-1 beta (IL-1β), interleukin-6 (IL-6), interleukin-8 (IL-8), interleukin-10 (IL-10) and transforming growth factor beta 1 (TGFβ1)] during the early phase response of bovine immature macrophages to in vitro infection with live Streptococcus agalactiae. Next generation sequencing of small RNA libraries from 20 cultures of blood monocyte-derived macrophages exposed to either one of two sequence types of S. agalactiae (ST103 or ST12) for 6 h in vitro and unchallenged controls was performed.
Results
Analyzes of over 356 million high quality sequence reads, revealed differential expression of 17 and 44 miRNAs (P < 0.05) in macrophages infected with ST103 and ST12, respectively, versus unchallenged control cultures. We also identified the expression of 31 potentially novel bovine miRNAs. Pathway analysis of the differentially regulated miRNAs and their predicted target genes in the macrophages infected with ST12 revealed significant enrichment for inflammatory response and apoptosis, while significant enrichment for integrin and GABA signaling were found in ST103 infected macrophages. Furthermore, both bacterial strains regulated miRNAs involved in the alternative activation of macrophages. The transcript levels of TNF-α, IL-1β, IL-6, IL-8 and IL-10 were significantly up-regulated by both bacterial strains, however the expression of TGFβ1 was significantly down-regulated only by ST12.
Conclusions
Our study identified pathogen-induced differential regulation of miRNAs controlling inflammation and polarization in bovine macrophages. This implies that miRNAs have potential to serve as biomarkers for early detection of bacterial infection.
Keywords
Background
MicroRNAs (miRNAs) are an abundant class of non-coding, small RNA molecules (19–24 nt long), which bind to the 3’UTR of target mRNAs to repress the translation into protein or accelerate the decay of expressed transcripts (affecting 60% of all mRNA transcripts [1]). Microorganisms are known to alter the host transcriptome and proteome in several ways, affecting miRNA is one of several such mechanisms. For example, miRNA induced by herpesviridae and Hepatitis C virus appear to work as an escape mechanism by dampening the host’s immune system [2] and miRNA induced by the latter is currently being targeted for therapy [3]. Up-regulation of miRNA will mostly attenuate the immune response, but it is unclear if this should be understood as a pathogen escape mechanism or a survival mechanism on the part of the host to avoid immunopathology.
Only a few studies describe miRNA regulation in the context of bovine bacterial mastitis [4]. Modification of miRNA expression were reported in response to Streptococcus uberis infection of primary bovine mammary epithelial cells (BMEs) and circulating monocytes from blood and milk [5, 6]. In addition, profiling of the complete miRNA content (miRNome) of BMEs infected with Staphylococcus aureus and Escherichia coli [7], and milk exosomes from S. aureus infected cows [8], revealed several pathogen directed miRNAs with enriched role in immunity, infection and cellular processes. MiRNAs have also gained prominence as potential biomarkers for a range of infections and diseases, and more recent studies profiling serum miRNAs from a bovine Mycobacterium avium subsp. paratuberculosis infection model demonstrated high stability of circulating miRNAs [9].
Macrophages are critical effectors and regulators of inflammation serving as the first line of defense against invading pathogens. Intramammary infections will activate macrophages to produce pro-inflammatory cytokines [e.g. tumor necrosis factor alpha (TNFα), interleukin-1 beta (IL-1β), interleukin-6 (IL-6) and interleukin-8 (IL-8)] required to kill intracellular pathogens. However, the balance between pro- and anti-inflammatory signals [e.g. interleukin-10 (IL-10) and transforming growth factor beta 1 (TGFβ1)] is crucial for immune regulation of the inflammation and preventing chronic conditions. The presence of microbes together with the host microenvironment underlies the diversity of activation of macrophages [10], and in the bovine mammary gland, the status of differential macrophage activation may be pivotal for the defense and resolution of mastitis [11].
In this study, we examined the expression of selected genes and the miRNA profiles of immature macrophages (primary cells) infected in vitro with the ST103 or ST12 strain of Streptococcus agalactiae, a leading causative agent of subclinical bovine mastitis. ST103 is prevalent in farms with considerable environmental contamination, while ST12, a strain usually associated with colonization of pregnant women, has been found in cattle herds with no positive environmental samples [12, 13]. To the best of our knowledge, this study represents the first report of next generation sequencing used to profile the host miRNAs response to S. agalactiae. Furthermore, the identification of differentially expressed miRNAs in infected immature macrophages may provide a basis for development of biomarker assays for early detection of subclinical mastitis.
Methods
Animals and cell isolation
Six healthy Norwegian Red (NR) cows aged 4–7 years were used for miRNAs sequencing, and six healthy NR cows aged 2.5–7 years were used for reverse transcription-quantitative PCR (RT-qPCR). Five of the cows were included in both experiments (for details see Additional file 1: Table S1 and Additional file 2: Figure S1). All animals were maintained under uniform housing conditions and nutritional regimens at the Norwegian University of Life Sciences (NMBU) herd. Blood sampling was performed by certified personnel and conducted in agreement with the provisions enforced by the Norwegian Animal Research Authority. Five hundred ml of blood was collected from the neck of each animal in sterile glass bottles with sodium citrate as anticoagulant. Peripheral blood mononuclear cells (PBMC) were extracted by density gradient centrifugation (2210×g, 30 min) on lymphoprep (Axis-Shield, Norway). CD14+ cells were extracted by positive selection of monocyte differentiation antigen CD14 using anti-human CD14 MACS MicroBeads (Miltenyi Biotec GmbH, Bergisch Gladbach, Germany), according to the manufacturer’s instructions and as described earlier by Lewandowska-Sabat et al. [14]. Purity of selected cells was verified by flow cytometry by staining positively selected cells with PE conjugated anti-mouse IgG2a (Southern Biotech, Birmingham, AL, USA), analyzing in a Gallios flow cytometer (Beckman Coulter), and found to be in the range of 95–98%. The CD14+ cells were subsequently grown in 6-well dishes at a density of 1.5 × 106 cells per well in RPMI medium supplemented with 10% FCS (Invitrogen, Carlsbad, USA). Cells were left over night at 37 °C in an atmosphere with 5% CO2. The phenotypic morphology of cells, i.e. differentiation of monocytes into an early-stage adherent macrophage phenotype, was visualized and confirmed by phase contrast microscopy.
Bacterial infection
Two previously described S. agalactiae strains (ST103 and ST12) were obtained from The Norwegian Veterinary Institute [12]. These bovine adapted strains were originally isolated from milk samples. Bacteria were collected from blood agar plates and grown in “Todd Hewitt broth” (Sigma-Aldrich) until mid-log phase. Growth was measured by optical density (OD) at 600 nm. The cultures were further aliquoted and frozen in 20% glycerol stocks in − 70 °C, and the final number of colony-forming units (CFU) was determined by serial dilutions and plating on blood agar plates. Bacteria used in this study all came from aliquots of the same batch.
For each individual animal the wells with immature macrophages were grouped into four classes with as equal number of wells and cells per class as possible. Two classes were infected with ST103 or ST12, in a multiplicity of infection (MOI) of 1 (1 bacterium per cell, on average). The third cell class (positive control) was exposed to 1 mg/mL of lipopolysaccharides (LPS, rough strains) from Salmonella minnesota Re 595 (Re mutant, Sigma-Aldrich) and the last class was left uninfected (negative control). After 1 h of exposure, 1% of penicillin/streptomycin (60 pg/mL penicillin and 100 μg/mL streptomycin) were added to prevent growth of remaining extracellular bacteria. The controls and the infected cells were treated equally. Inhibition of bacterial growth by antibiotics was verified by microscopy. Incubation was continued for one additional hour for LPS-exposed cells and for 5 more hours for bacteria infected and negative control cell classes, a total of 6 h incubation (Additional file 2: Figure S1). Media was aspirated and the cells were collected using cell scraper. Cells were centrifuged (400×g, 5 min), the pellet was washed with cold PBS buffer, snap frozen in liquid nitrogen and stored at − 70 °C.
RNA extraction and reverse transcription-quantitative PCR
Twenty-four RNA samples from six animals were used (Additional file 1: Table S1, Additional file 2: Figure S1). Total RNA was isolated from control and infected cells using the MirVANA isolation kit (Ambion, Austin, TX) following the manufacturer’s instructions. All RNA samples were treated with amplification grade DNase I (Invitrogen) to remove any traces of genomic DNA. RNA concentration and quality was measured using NanoDrop 1000 (Thermo Fisher Scientific, Wilmington, USA) and 2100 BioAnalyzer (Agilent RNA 6000 Nano, Agilent Technologies, Palo Alto, USA), respectively. The RNA integrity numbers (RIN), concentrations and OD A260/A280 ratios are listed in Additional file 1: Table S1. For RT-qPCR a total of 200 ng RNA was used for cDNA synthesis using Tetro cDNA synthesis kit (Nordic BioSite, Norway), and cDNA equivalent to 5 ng of total RNA was used in qPCR reactions set up in triplicate for each sample using Express SYBR GreenER SuperMix with premixed ROX (Invitrogen) according to the manufacturer’s recommendations. Transcript levels were analyzed using a 7900HT Fast Real-Time PCR System (Applied Biosystems) and the standard program: 50 °C for 2 min, 95 °C for 2 min, 40 cycles of 95 °C for 15 s and 60 °C for 1 min, followed by melting curve analyses. Gene-specific primers were either derived from literature or designed using Primer3 ver. 0.4.0 [15]. The transcript levels of the following genes were analyzed: TNF-α, IL-1β, IL-6, IL-8, IL-10 and TGFβ1. The primer sequences are presented in Additional file 3: Table S2. The efficiencies of all primer pairs were tested by template dilution series using pooled cDNA from control and infected cells and were 100% (±10). Negative controls with no added template were included for all primer pairs (no template control), and no RT control reactions for each sample and each primer pair were run in qPCR in order to check for genomic DNA contamination (no RT control). The peptidylprolyl isomerase A (PPIA) housekeeping gene were used in the current study, as it has been shown to be one of the most stable genes for gene expression studies in cattle macrophages [14] and lymphocytes [16], and in human LPS-stimulated monocytes [17]. In the experiment, PPIA was expressed at the same level in the cells stimulated with both bacterial strains and LPS as in the negative controls. Initial analysis of the RT-qPCR data was performed using RQ Manager 1.2 (Applied Biosystems). Standard deviation of ≤0.3 per triplicate was accepted. The ΔCt method was used to calculate RT-qPCR data, i.e. ΔCt = Cttarget gene – Ctreference gene, and normalized gene expression was calculated as 2 (-ΔCt). Fold change was calculated relative to the negative, unexposed control. Reciprocal values of fold change were used for down-regulated gene expression (i.e. TGFβ1) in order to facilitate the interpretation of the results. The differences of normalized gene expression levels between control and infected cells for each gene were tested using Wilcoxon matched-pairs signed rank test using GraphPad Prism version 7.00 for Windows (GraphPad Software, La Jolla California USA, www.graphpad.com). The pairwise differences between responses to LPS, ST103 and ST12 were tested using log-transformed fold change values for each treatment and gene by RM one-way ANOVA and Tukey’s multiple comparison test using GraphPad Prism. The significance level was determined at P < 0.05.
MicroRNA sequencing and data analyses
Twenty RNA samples from six animals were used (control = 6, LPS = 4, ST12 = 5, ST103 = 5; see Additional file 1: Table S1 and Additional file 2: Figure S1 for details) for short non-coding RNA deep sequencing. Total RNA was isolated from control and infected cells using the MirVANA isolation kit, which retains short RNA fragments, and RNA concentration and quality were measured as described above for RT-qPCR (see Additional file 1: Table S1 for details). RNA-seq libraries were prepared with NEBNext® Multiplex Small RNA Library Prep Set for Illumina (New England Biolabs), according to the manufacturer’s protocol. Libraries were sequenced (75 bp single-end, 3–400 million reads) on one lane with an Illumina NextSeq 500 machine. Preparation of RNA-seq libraries and sequencing were performed by the Norwegian Sequencing Centre (Oslo, Norway; http://www.sequencing.uio.no/).
The miRNA sequence data was analyzed using Oasis ver. 2.0, a web application that allows for fast and flexible online analysis of small-RNA-seq (sRNA-seq) data ([18]; http://oasis.dzne.de). Preliminary quality control analysis of the 20 fastq files was carried out with FASTQC software ver. 0.11.2 (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/). Cutadapt ver. 1.2.1 (https://cutadapt.readthedocs.io/en/stable/) was then used to trim 3′ adaptor sequences. Reads, shorter than 17 and longer than 32 nucleotides after trimming were discarded. Preprocessed reads, which successfully passed filtering were aligned to the bovine reference genome assembly (bostau8) using STAR ver. 2.30.0e_r291 [19] in non splice-junction-aware mode, allowing 5% mismatch of the read length (0 mismatches for reads with length 17–19, 1 mismatch for reads with length 20–32). Reads that could only be aligned using soft trimming (trimming of the beginning or end of reads that is not counted as mismatches) were removed from STAR’s output using a custom script. Quantification of uniquely mapped reads was performed using featureCounts ver. 1.4.6 [20], and assignment to the short RNA species was performed using miRBase ver. 21 (release date: 2014–6-22) and Ensembl (v84).
In addition to profiling the expression of known miRNAs, the prediction of potentially novel miRNAs was performed by miRDeep2 ver. 2.0.0.5 software [21] on the combined set of all input fastq files (i.e., all the input samples were merged into a unique file containing the total reads). miRDeep2 analysis were performed with the default parameters except for the minimum read depth (−a). This parameter is automatically computed according to the total number of uniquely mapped reads, targeting a performance of at least 80% true positives and 80% recall. Further filters were applied in order to identify those novel miRNAs that have the highest likelihood of being true positives. These are based on significant randfold p-value (default mirdeep2 threshold), miRDeep2 score > 5, where both the mature and star read counts were expressed with a minimum of 5 reads each and rFAM alerts for other types of RNA (e.g. rRNAs and tRNAs). As Oasis maps to the mature miRNA sequence, miRDeep2 code was modified in order to show the mature miRNA position instead of the pre-mature miRNA position.
Differential expression analysis
In order to identify statistically significant differentially expressed (DE) miRNAs between control and treatment samples, DESeq2 ver. 1.4.5. [22] was applied. The analyses were performed pairwise, i.e. negative control and one of the treatment groups at the time (LPS, ST103 or ST12, respectively) were analyzed. In addition, a pairwise comparison between the two bacterial treatment groups (i.e. ST103 and ST12) was performed. The miRNA read counts identified by miRDeep2 were normalized across all samples using DESeq2 normalization method [22]. DE miRNAs were defined as having a Benjamini and Hochberg corrected P-value of < 0.05 and with the minimum number of miRNAs (mean expression values) with an average of 5 reads for either biological condition. Principal component analysis (PCA) and hierarchical clustering of the samples were performed in DESeq2, in order to determine how well samples cluster together based on the similarity of their sRNA expression. The analyses were performed using Oasis.
Target prediction and pathway analysis
DE miRNAs identified in ST103 (n = 17) vs. negative control and in ST12 (n = 44) vs. negative control were used as input lists in Ingenuity Pathway Analysis (IPA; https://www.qiagenbioinformatics.com/products/ingenuity-pathway-analysis/) in order to identify biological functions and networks that were overrepresented in the datasets.
Target genes that are potentially regulated by DE miRNAs were predicted using TargetScan ver.7.1 [23]. Given the high false positive rates for miRNA target prediction, we identified only those potential target genes that had cumulative weighted context++ score < − 0.5 [23].
Identification of statistically overrepresented biological pathways in the lists of target genes of the DE miRNAs were performed by pathway analysis using InnateDB, a curated database of innate immunity genes, pathways and molecular interactions (www.innatedb.com; [24]). Target genes predicted by TargetScan were converted to human homologs and further analyzed in InnateDB. The pathway overrepresentation analyses were performed using hypergeometric test and Benjamini and Hochberg corrected P-value ≤0.1 was defined as significant.
Results
Reverse transcription-quantitative PCR
mRNA gene expression levels in bovine monocyte-derived macrophages as determined by reverse transcription-quantitative PCR (RT-qPCR). Genes displaying a significant difference in expression between the Streptococcus agalactiae strains ST103 and ST12, and the lipopolysaccharides (LPS) stimulated cells, respectively, compared to uninfected controls are denoted with * (P ≤ 0.05). The significant differences between responses to LPS vs. ST103, and LPS vs. ST12, respectively, are denoted with # (P ≤ 0.05). Data are presented as fold change relative to negative control, and as a mean values with SD. Tumor necrosis factor alpha (TNFα), interleukin-1 beta (IL-1β), interleukin-6 (IL-6), interleukin-8 (IL-8), interleukin-10 (IL-10) and transforming growth factor beta 1 (TGFβ1). Peptidylprolyl isomerase A (PPIA) was used as a reference gene
Identification of small RNAs in bovine monocytes
Next, we conducted miRNA-seq on 20 sample libraries from the bovine monocyte-derived macrophages. Preliminary quality control analysis of the resulting 20 fastq files revealed that all libraries passed the quality criteria with a Phred score > 32. The sequencing resulted in a total of 356,741,900 high quality reads. Among them, 349,785,433 sequences (98%) remained after adaptor trimming. Approximately 43% of all reads were discarded due to length filtering after adaptor trimming. The proportion of reads uniquely aligning to different RNA biotypes demonstrated that miRNAs were the dominant class of small RNAs sequenced in our small RNA libraries. The majority (> 86%) of all reads in the 17 to 32 nucleotide fraction aligned to known miRNAs. Almost all the remaining reads mapped to small nucleolar RNAs (snoRNAs). The number of initial reads per sample, the percentage of trimmed reads and the percentage of uniquely mapped reads for each sample are presented in Additional file 4: Table S3.
The top ten most abundant miRNAs across all samples. Bovine monocyte-derived macrophages infected in vitro with Streptococcus agalactiae strains ST103 or ST12 or lipopolysaccharides (LPS), and uninfected negative controls. Mean miRNAs expression values (log transformed normalized mean number of reads) are presented on the x-axis, as calculated by DESeq2 software (at least 5 mean reads for either biological condition)
Novel microRNAs identified in bovine monocyte-derived macrophages with or without challenge with lipopolysaccharides (LPS), Streptococcus agalactiae strains ST103 and ST12
provisional ID | mature read count | star read count | miRDeep2 score | example miRBase miRNA with the same seed | consensus mature sequence | precursor coordinate |
|---|---|---|---|---|---|---|
chr25_2078 | 63,905 | 8,00 | 32,589,40 | mmu-miR-106b-3p | ccgcacuguggguacuugcugc | chr25:36892057..36892117:+ |
chr8_3368 | 48,256,00 | 44,00 | 24,641,00 | hsa-let-7d-3p | cuauacgaccugcugccuuucu | chr8:86887438..86887514:+ |
chrX_3864 | 40,373,00 | 6,00 | 20,590,60 | hsa-miR-223-5p | cguguauuugacaagcugaguug | chrX:99936328..99936387:- |
chr7_3248 | 14,956,00 | 13,00 | 7636,00 | mmu-miR-24-1-5p | gugccuacugagcugaaacacag | chr7:12981645..12981702:- |
chr26_2162 | 10,315,00 | 54,00 | 5290,60 | eca-miR-146b-3p | ugcccuagggacucaguucuggu | chr26:22930915..22930975:+ |
chr11_440 | 4397,00 | 11,00 | 2251,80 | ssc-miR-181d-3p | accaccgaccguugacuguacc | chr11:95709449..95709509:+ |
chr8_3359 | 3633,00 | 95,00 | 1990,00 | mmu-miR-27b-5p | agagcuuagcugauuggugaaca | chr8:83009841..83009902:+ |
chr25_2117 | 3826,00 | 7,00 | 1958,50 | gra-miR7486h | cagcaacuaaagaucccucagg | chr25:34409297..34409357:- |
chr22_1816 | 3325,00 | 78,00 | 1741,70 | bmo-miR-3000 | cugcgcuuggauuucguuccc | chr22:51543484..51543548:+ |
chr3_2544 | 1854,00 | 9,00 | 951,00 | bta-miR-2285f | aaaaaccugaaugacccuuuug | chr3:94548590..94548649:+ |
chr8_3366 | 1064,00 | 754,00 | 932,60 | hsa-let-7a-3p | cuauacaaucuauugccuuccc | chr8:86885231..86885309:+ |
chr7_3219 | 1392,00 | 74,00 | 748,70 | hsa-miR-378a-5p | cuccugacuccagguccugugu | chr7:63067305..63067362:+ |
chr3_2643 | 863,00 | 137,00 | 513,90 | hsa-miR-30c-2-3p | cugggagaggguuguuuacucc | chr3:106059376..106059437:- |
chr16_1041 | 937,00 | 11,00 | 486,50 | hsa-miR-181a-3p | accaucgaccguugauuguacc | chr16:79685955..79686018:- |
chr14_712 | 672,00 | 34,00 | 364,10 | hsa-miR-30a-3p | cuuucagucagauguuugcugcu | chr14:8080297..8080360:+ |
chr1_143 | 643,00 | 10,00 | 335,60 | mmu-miR-15b-3p | cgaaucauuauuugcugcucuag | chr1:107923396..107923457:- |
chr12_549 | 565,00 | 17,00 | 301,20 | hsa-miR-92a-1-5p | agguugggaucgguugcaaugcu | chr12:66227265..66227321:+ |
chr1_111 | 475,00 | 6,00 | 249,50 | mmu-miR-125b-2-3p | acaagucaggcucuugggacc | chr1:19881359..19881419:- |
chr16_1003 | 359,00 | 103,00 | 240,00 | efu-miR-9283 | uguggccucuggguguguacccuc | chr16:33022297..33022356:- |
chrX_3750 | 463,00 | 6,00 | 240,00 | ssc-miR-374a-3p | uuaucagguuguauuguaauu | chrX:81951234..81951286:+ |
chr18_1154 | 297,00 | 42,00 | 181,00 | hsa-let-7e-3p | cuauacggccuccuagcuuucc | chr18:58015043..58015110:+ |
chr4_2734 | 287,00 | 41,00 | 168,50 | – | acacgcguccuuggauccugacu | chr4:119142851..119142912:+ |
chrX_3767 | 228,00 | 37,00 | 139,50 | hsa-miR-222-5p | cucaguagccaguguagaucc | chrX:103538171..103538234:+ |
chrX_3861 | 197,00 | 16,00 | 113,20 | hsa-let-7a-3p | cuauacaacuuacuacuuuccc | chrX:96382645..96382725:- |
chr19_1372 | 157,00 | 5,00 | 82,60 | – | agggagucccugguaguucagu | chr19:46741763..46741851:- |
chr19_1308 | 103,00 | 22,00 | 65,00 | – | ccccggcuuuuccucccccagg | chr19:62308493..62308542:+ |
chr5_2875 | 86,00 | 8,00 | 52,80 | ssc-miR-7134-5p | auguccgcggguucccugucc | chr5:112083860..112083922:+ |
chr19_1281 | 70,00 | 5,00 | 43,00 | hsa-miR-152-5p | agguucugugauacacuccgacu | chr19:39081179..39081236:+ |
chr18_1135 | 13,528,00 | 95,00 | 5,60 | mmu-miR-140-5p | cagugguuuuacccuaugguag | chr18:37088153..37088219:+ |
chr19_1237 | 93,00 | 15,00 | 5,40 | ahy-miR3511-5p | accagggcuggaagcugcuucu | chr19:9534051..9534107:+ |
chr2_1498 | 1928,00 | 20,00 | 5,30 | hsa-miR-26b-3p | ccuguucuccauuacuuggcucg | chr2:107133408..107133466:+ |
Differentially expressed miRNAs in macrophages challenged with S. agalactiae
The difference between the negative control and the ST103 and ST12 infected samples, respectively, was confirmed by PCA of the small RNA sequence reads. For ST103 infected samples, PC1 explained 31% and PC2 explained 22% of the overall miRNA expression variability. While for ST12 infected samples, PC1 explained 39% and PC2 explained 20% of the overall miRNA expression variability. No clear difference between the negative control and the LPS challenged samples was revealed by PCA (Additional file 6: Fig. S2).
MicroRNAs significantly differentially expressed between bovine monocyte-derived macrophages infected with Streptococcus agalactiae strains ST103 or ST12 and the respective uninfected controls
ST103 | ST12 | ||||
|---|---|---|---|---|---|
miRNA ID | Log2FC | P-value | miRNA ID | Log2FC | P-value |
p-bta-miR-3 | 8,22 | 8,77E-45 | p-bta-miR-3 | 9,98 | 1,10E-53 |
bta-miR-2284i | 7,37 | 2,28E-29 | bta-miR-2284i | 8,62 | 9,53E-41 |
bta-miR-2285m | 4,97 | 5,87E-13 | bta-miR-2438 | 7,93 | 2,46E-25 |
bta-miR-7858 | 4,24 | 1,15E-06 | bta-miR-2285 m | 6,08 | 9,41E-17 |
bta-miR-222 | 1,75 | 1,05E-05 | bta-miR-2478 | −2,29 | 4,69E-15 |
bta-miR-146b | 2,34 | 1,88E-04 | bta-miR-223 | 1,27 | 3,63E-14 |
bta-miR-2427 | −2,53 | 1,60E-03 | bta-miR-1249 | −3,29 | 1,77E-13 |
bta-miR-2898 | −1,52 | 2,88E-03 | bta-miR-128 | −1,91 | 2,34E-11 |
bta-miR-2478 | −0,93 | 4,00E-03 | bta-miR-2427 | −4,19 | 1,55E-09 |
bta-miR-628 | 1,20 | 6,59E-03 | bta-miR-2898 | −2,92 | 2,98E-09 |
bta-miR−1306 | −1,73 | 8,43E-03 | bta-miR-500 | 1,39 | 4,76E-09 |
bta-miR-1249 | −1,54 | 1,14E-02 | bta-miR-92b | −2,27 | 2,08E-08 |
bta-miR-708 | 1,08 | 1,85E-02 | bta-miR-484 | −1,81 | 3,88E-08 |
bta-miR-221 | 0,69 | 3,71E-02 | bta-miR-365-3p | −2,36 | 4,13E-07 |
bta-miR-1246 | −0,94 | 4,33E-02 | bta-miR-1306 | −3,34 | 8,34E-07 |
bta-miR-2892 | − 1,57 | 4,97E-02 | bta-miR-374b | 1,57 | 6,66E-06 |
bta-miR-9-5p | 1,13 | 4,97E-02 | bta-miR-628 | 1,50 | 7,96E-06 |
bta-miR-197 | −2,21 | 9,36E-06 | |||
bta-miR-425-5p | 0,93 | 3,38E-05 | |||
bta-miR-340 | −1,39 | 1,28E-04 | |||
bta-miR-30d | −1,31 | 1,72E-04 | |||
bta-miR-425-3p | 1,04 | 1,82E-04 | |||
bta-miR-505 | −1,15 | 1,19E-03 | |||
bta-miR-125a | −2,01 | 1,22E-03 | |||
bta-miR-146b | 2,11 | 1,64E-03 | |||
bta-miR-1343-3p | −1,59 | 1,64E-03 | |||
bta-miR-2388-5p | −2,67 | 2,77E-03 | |||
bta-miR-423-3p | −0,98 | 2,85E-03 | |||
bta-miR-328 | −1,84 | 3,39E- 03 | |||
bta-miR-30b-5p | 1,26 | 3,63E-03 | |||
bta-miR-92a | −0,91 | 6,29E-03 | |||
bta-miR-1468 | −0,88 | 6,32E-03 | |||
bta-miR-30f | −1,34 | 6,58E-03 | |||
bta-miR-125b | −1,18 | 8,86E-03 | |||
bta-miR-10a | −0,76 | 9,31E-03 | |||
bta-miR-2431-3p | −2,25 | 9,76E-03 | |||
bta-miR-155 | 1,33 | 1,76E- 02 | |||
bta-miR-361 | −0,97 | 1,97E-02 | |||
bta-miR-122 | 2,12 | 2,12E-02 | |||
bta-miR-221 | 0,65 | 2,16E-02 | |||
bta-miR-2284ab | −1,06 | 3,02E-02 | |||
bta-miR-2284w | −0,74 | 3,02E-02 | |||
bta-miR-30c | −0,82 | 3,62E-02 | |||
bta-miR-669 | −0,85 | 4,39E-02 | |||
Heatmap of the top most significant differentially expressed microRNAs (P < 0.05) between bovine monocyte-derived macrophages infected in vitro by a) Streptococcus agalactiae strain ST103 or b) Streptococcus agalactiae strain ST12 and non-infected negative controls, as determined by DESeq2 analyses. Red color indicates high expression of the microRNA
Venn diagram of differentially expressed miRNAs (P < 0.05) between Streptococcus agalactiae strain ST103 or ST12 infected and respective non-infected (negative controls) bovine monocyte-derived macrophages. Red indicates miRNAs that were up-regulated, and blue indicates miRNAs that were down-regulated compared to the negative control samples
One novel miRNA (p-bta-miR-3) was found to be DE in response to both bacterial strains (Fig. 4, Additional file 7: Table S5). This 20 nt long miRNA was mapped to bovine chromosome 9 and had a significant homology to mmu-mir-7025-5p (E-value 6.2).
Enriched pathways of differentially regulated miRNAs and their predicted target genes
Functional network overrepresented in the list of differentially expressed miRNAs in response to infection of bovine macrophages with Streptococcus agalactiae strain (a) ST103 and (b) ST12. Ingenuity Pathway Analysis (IPA) identified 6 associated molecules with a network score of 16 for ST103 and 13 associated molecules with a network score of 31 for ST12
Top significant pathways overrepresented among target genes of differentially expressed miRNAs in response to exposure of bovine macrophages in vitro to Streptococcus agalactiae strains ST103 or ST12
Pathway name | P-value | P-value (B-H) | Genes | ||
|---|---|---|---|---|---|
ST103 | miRNAs up-regulated | Mucin type O-Glycan biosynthesis | 1.4E-4 | 0.01 | GALNT3; GALNT4; GALNT9; POC1B-GALNT4; |
GABA receptor activation | 1.0E-4 | 0.02 | GABRA1; GABRB2; GNG10; GNG5; KCNJ2; | ||
Regulation of gene expression in early pancreatic precursor cells | 0.001 | 0.06 | ONECUT1; ONECUT3; | ||
miRNAs down-regulated | Integrin signaling pathway | 0.004 | 0.10 | ACTN1; ACTN2; BCR; CAV1; FYN; ITGA1; MAPK8; PXN; SOS1; TNS1; VCL; | |
Sodium/Calcium exchangers | 0.003 | 0.11 | SLC24A1; SLC24A2; SLC24A3; SLC24A4; SLC8A1; SLC8A2; SLC8A3; | ||
Reduction of cytosolic Ca++ levels | 0.004 | 0.11 | ATP2A2; ATP2B1; ATP2B2; ATP2B3; ATP2B4; SLC8A1; SLC8A2; SLC8A3; | ||
ST12 | miRNAs up-regulated | Activation of G protein gated Potassium channels | 2.3E-4 | 0.04 | GNB3; GNG10; GNG2; GNG5; KCNJ6; |
G protein gated Potassium channels | 2.3E-4 | 0.04 | GNB3; GNG10; GNG2; GNG5; KCNJ6; | ||
Inhibition of voltage gated Ca2+ channels via Gbeta/gamma subunits | 2.3E-4 | 0.04 | GNB3; GNG10; GNG2; GNG5; KCNJ6; | ||
miRNAs down-regulated | IL4-mediated signaling events | 5.1E-4 | 0.15 | BCL2L1; CCL11; CCL26; FCER2; IL10; IRF4; IRS2; SOCS1; SOCS5; STAT6; | |
Intrinsic Pathway for Apoptosis | 3.6E-4 | 0.16 | BAK1; BAX; BCL2L1; BMF; CASP7; TFDP1; YWHAB; YWHAG; | ||
Role of parkin in ubiquitin-proteasomal pathway | 3.1E-4 | 0.28 | PARK2; UBE2E2; UBE2G1; UBE2L3; |
Discussion
MiRNAs are suggested to be fine-tuners of gene expression during inflammatory response [25, 26]. This indicates a great potential of the genome-wide miRNome profiling for investigation of inflammation. In this study, we have identified the early phase miRNome of bovine monocyte-derived macrophages infected in vitro with two strains of S. agalactiae derived from Norwegian dairy herds. S. agalactiae is a major causative agent of subclinical mastitis and an increasing problem in Norway. Furthermore, S. agalactiae is widely recognized as a main cause of life-threatening infections in human neonates, pregnant females, and elderly adults [27]. The two strains used here differ in their ability to survive in the environment and transmit within dairy herds. ST103 is the most predominant isolate in bovine mastitis found in farms with substantial environmental contamination, while ST12 was found in cattle herds with no positive environmental samples [12]. Using modern sequencing of miRNAs from macrophages infected with S. agalactiae and analyzing over 356 million reads, we found that several miRNAs were significantly differentially regulated in response to infection (Figs. 3 and 4). We observed that fewer miRNAs were differentially regulated in macrophages challenged with strain ST103 than with strain ST12 (17 and 44, respectively; Fig. 4). This suggests that ST103 induces relatively subtle changes in miRNA expression during the early stage of infection. However, the expression pattern of the examined pro- and anti-inflammatory genes was not significantly different between these two strains, except for TGFβ1, that was significantly down-regulated only by ST12 (Fig. 1).
It has been shown that different bacterial strains are often associated with a variety of virulence factors, which, in turn, are associated with the outcome of infection in the host [28, 29]. The latest study of molecular and virulence characteristics of the most prevalent bovine S. agalactiae isolates has shown that these factors in strain ST103 (but also in other of the most prevalent bovine isolates) are associated with growth ability in milk, biofilm formation and adhesion to bovine mammary epithelial cells. Furthermore, ST103 showed significantly higher hemolytic activity and cytotoxicity compared to most of the other bovine strains [30]. On the contrary, the characterization of ST12 revealed that this strain shows high rates of resistance to erythromycin and clindamycin [31] and can survive within macrophages [32]. The observed differences in the miRNA expression profiles during the early stage of macrophage infection may be explained by the different characteristics of these two S. agalactiae strains such as the differences in virulence factors, the different niches and modes of transmission.
Both strains increased the level of miR-146b and IL-10 transcripts. This may be associated with toll-like receptor (TLR) recognition of bacterial patterns, as shown in monocytes for other bacterial TLR-agonists. MiR-146b has been reported to reduce the expression of several pro-inflammatory cytokines and chemokines in human monocytes upon stimulation with LPS [33], which suggests that miR-146b has an anti-inflammatory activity in monocytes. In a previous transcriptomic study, we found indications that another gram-positive bacterial mastitis pathogen, S. aureus, promotes bovine macrophages to be activated in an alternative manner [14]. Here we observed several indications of alternative activation on the miRNA level as a result of streptococcal uptake in macrophages. Both strains induced miR-221, ST12 induced miR-30b, miR-223, miR-374b and miR-500 but down-regulated miR-125a and miR-125b, while ST103 induced miR-222, all associated with alternative macrophage activation [34–37]. While one study reported that miR-125a promotes classical macrophage activation [38], and several of the DE miRNA found in our study are not accounted for in the literature, the results suggest that S. agalactiae may induce an alternative-like macrophage activation similar to S. aureus. However, this result needs to be confirmed by more comprehensive functional studies.
One of the novel miRNA (p-bta-miR-3) was found to be DE in response to both bacterial strains (Fig. 4). This miRNA shares 100% sequence identity and covers 95% of the query (i.e. 19 out of 20 nucleotides match) when blasted against the genome sequences of several S. agalactiae strains (e.g. GenBank: CP025028.1; https://blast.ncbi.nlm.nih.gov/Blast.cgi). This suggests that it may be a bacterial-derived miRNA released by infection of the macrophages. Bacterial-derived miRNAs have recently been extensively studied [39, 40]. It is likely that these small RNA molecules may repress the ability of the host cell to resist bacterial invasion, as has been shown for viruses and fungi [41–43]. It is not clear if bacterial miRNAs require the host cell miRNA processing machinery for its biogenesis or if it is a bacterial RNA fragment released into the host cell cytoplasm during bacterial lysis. Whether p-bta-miR-3 observed in our study is an authentic miRNA of bacterial origin needs to be further tested.
Mapping of bovine miRNAs to their human homologs may reveal conserved pathways and biological functions associated with the regulation of gene expression during immune response to bacteria. To characterize the distinct DE miRNA profiles of the two bacterial strains we performed pathway analysis using IPA. As a result 2 genes, dicer 1 ribonuclease III (DICER1) and argonaute 2 (AGO2), were found to be associated with the ST103 DE miRNAs (e.g., miR-9-5p and miR-708; Fig. 5a). DICER and AGO2 are miRNA processing enzymes, and DICER1 knockdown cell lines and animals have been correlated with a global drop of mature miRNA expression levels, shifting macrophages towards a classically activated phenotype [44–46]. These results show that strain ST103 induces miRNAs that likely contribute to maintain an alternative macrophage phenotype.
For miRNAs down-regulated in ST103 infected macrophages, putative target genes with a role in the integrin signaling pathway may be noted (Table 3). Integrin activation is necessary for leukocyte arrest before migration across the endothelial cell barrier and many bacterial species, such as group A streptococci [47] S. aureus [48], and E. coli [49] use host integrins for adhering to and invading host cells. Moreover, the study of virulence factors of bovine S. agalactiae isolates has identified virulence factors associated with adhesion ability to bovine mammary epithelial cells in strain ST103 [30]. These findings, together with our results may indicate that the integrin signaling pathway is essential for S. agalactiae ST103 strain invasion of macrophages or mammary epithelial cells.
Furthermore, target genes of ST103 up-regulated miRNAs are involved in gamma-aminobutyric acid (GABA) receptor activation (Table 3). GABA plays an inhibitory role in autoimmune inflammation [50] and GABA receptor transcripts are present in different immune cells [50–53]. Inflammatory cytokine production in peripheral macrophages and T cell autoimmunity is decreased after GABA treatment [51, 52]. Moreover, GABA and GABA type A receptor agonists reduced cytotoxic immune responses of T cells [54]. Other target genes associated with up-regulated ST103 miRNAs belong to cancer, organismal injury, abnormalities, and reproductive system disease pathways, without obvious connection to the bacterial infection context (Additional file 8: Table S6). Taken together, ST103 induced miRNA regulation of integrin and GABA signaling, and maintenance of alternative macrophage activation suggest that miRNAs may be involved in the pathogenesis of subclinical mastitis. However, whether this is an evasion strategy developed by S. agalactiae ST103 to avoid host defense and promote its intracellular replication and persistence remains to be determined.
In contrast, IPA revealed that strain ST12 induced DE miRNAs were associated with inflammatory disease and inflammatory response (Fig. 5b, Additional file 8: Table S6). A number of genes associated with the ST12 induced miRNA networks are key players of immune response, i.e. IL-1, nuclear factor kappa β (NFKβ), interleukin-12 (IL-12), Janus kinase (JaK) and SMAD family member 2/3 (Smad2/3) (Fig. 5b). For three of the ST12-unique DE miRNAs (bta-miR-155, bta-miR-125b and bta-miR-223) evidence of induction by bacteria has been previously documented. It has been shown that miR-155 and miR-125b are involved in regulation of TNF production during mycobacterial infection [55, 56]. MiR-155 was reported to induce toll-like receptor 2 (TLR-2) recognition of Streptococcus pneumoniae, and promote effective bacterial clearance in the nasopharynx [57]. MiR-223 regulates inflammation (reviewed in [58]) and was found to be up-regulated in bovine mammary tissue infected with S. uberis [59]. Moreover, an earlier study shows that Salmonella can up-regulate intestinal epithelial miR-128 expression, which, in turn, decreases macrophage recruitment [60]. Finally, IPA revealed that miR-128 and miR-92a are associated with Smad2/3 proteins (Fig. 5b), critical downstream mediators of TGFβ-signaling (reviewed in [61]). This is consistent with the qPCR findings of down-regulation of TGFβ1 by ST12 (Fig. 1).
The predicted target genes of the down-regulated miRNAs in ST12 infected macrophages were significantly enriched for genes with roles in interleukin-4 (IL-4)-mediated signaling events and intrinsic pathway for apoptosis (Table 3). These genes have multiple roles in regulation of apoptosis [e.g. BCL-2 like 1 (BCL2L1), BCL-2 antagonist killer 1 (BAK1), BCL-2 associated X (BAX)] and anti-inflammation (e.g. IL-10), and display chemotactic activity for eosinophils [e.g. C-C motif chemokine ligand 11 (CCL11), C-C motif chemokine ligand 26 (CCL26)]. Pro-apoptotic BAK1 is a direct target of miR-125b [62] and in our study miR-125b was found to be down-regulated in response to ST12 (Fig. 4). Apoptosis plays a critical role in the pathogenesis of sepsis [63] and S. agalactiae strains are frequent agents of life-threatening sepsis and meningitis in human neonates and adults. This may indicate the putative mechanisms of pathogenesis, i.e. impaired immune responses due to extensive death of immune system cells [63]. However, further in vivo studies are needed in order to confirm this hypothesis.
The predicted target genes of ST12 up-regulated miRNAs have multiple roles in G protein signaling (Table 3). It has been shown that G-protein-coupled receptors (GPCRs) are linked to heterotrimeric G-proteins composed of α, β, and γ subunits (reviewed in [64]). Inflammatory cells such as macrophages express a large number of GPCRs for classic chemoattractants and chemokines, and these receptors are critical for both enhancement of inflammation and promotion of its resolution (reviewed in [65]). This may suggests that ST12-mediated up-regulation of miRNAs, which in turn regulate G-protein genes, may have impact on regulation of inflammatory gene expression, particularly those facilitating chemotaxis. Taken together, we found that strain ST12 induces a stronger inflammatory immune response in bovine macrophages than strain ST103, with miRNAs associated with down-regulation of the TGFβ/Smad signaling, chemotactic pathways and enhancement of apoptotic pathways and alternative macrophage activation.
Three out of thirteen DE miRNAs (i.e. bta-miR-2478, bta-miR-2898 and bta-miR-1249) detected among macrophages challenged with ST12 compared to macrophages challenged with ST103 were also found among DE miRNAs commonly regulated by both S. agalactiae strains (Fig. 4 and Additional file 7: Table S5). It was reported recently that bta-miR-2478 inhibits TGFβ1 expression during mammary gland development in goats and a strong negative correlation between miR-2478 and TGFβ1 expression was demonstrated [66]. We could not observe this negative correlation in our study, as the lowest transcript level of TGFβ1 (Fig. 1) was associated with the lowest expression of bta-miR-2478 (Additional file 5: Table S4) in macrophages challenged with ST12, compared to the other groups. However, it has also been reported that this miRNA is associated with energy metabolism and feeding-induced changes in muscles, and regulates different target genes in cattle [67, 68]. This suggests that bta-miR-2478 may target other gene(s) than TGFβ1 in infected macrophages. Bta-miR-2898 was found to be up-regulated in the mammary gland tissues of mastitis-infected cows, while bta-miR-1249 was identified as miRNA, which targets genome of H5N1 influenza A virus during in vitro infection [69, 70]. The significantly stronger down-regulation of these three miRNAs in macrophages challenged with ST12 compared to ST103 may be explained partly by the differences in virulence factors between these strains and suggest an important role of these miRNAs during the early stage of macrophage infection.
We were not able to identify any significant DE miRNAs in response to LPS in our study, however, several pro- and anti-inflammatory genes were up-regulated after exposing macrophages to LPS (Fig. 1). This may be due to higher individual variation in miRNA expression between LPS-treated samples combined with lower number of biological replicates (n = 4) compared to ST12 and ST103 (Additional file 6: Figure S2). However, we were not able to identify any outliers within the LPS-treated group based on the variation pattern (Additional file 6: Figure S2), in order to improve the statistical analysis. The number of biological replicates needed to ensure valid biological interpretation of high-throughput sequencing results is unclear. However, some studies attempt to answer this question and provide guidelines for experimental design. Earlier it has been proposed that the absolute minimum of three biological replicates should be used for next-generation sequencing studies, however, Schurch et al. [71] suggested that at least six biological replicates should be used in order to assure statistical power in the detection of differentially expressed genes. In the present study, we used four replicates for LPS-stimulated and five replicates for each of the bacterial strain-challenged macrophages. Furthermore, using five biological replicates for ST103- and ST12-treated samples may also have a negative impact on the number of detected significant miRNAs. Increasing the number of biological replicates would be beneficial to increase the statistical power and accuracy of the RNA and miRNA-seq analyses, but unfortunately, this was not possible within the scope of the present project. Originally, we had planned for six samples of each class of exposure, but some of the cultures failed to yield sufficient number of cells.
We observed approximately 250 miRNAs expressed in bovine macrophages (Additional file 5: Table S4). This is in agreement with the previous reports on miRNA sequencing of primary bovine mammary epithelial cells infected with S. uberis [6] and with S. aureus and E. coli [7], but it is approximately 3 times more than identified in bovine alveolar macrophages [72]. In our study, some miRNAs were expressed at very high levels (Fig. 2), while the majority were expressed at low levels. This is also observed in several previous studies of macrophages and epithelial cells miRNome [6–8, 72].
Conclusions
Our study have shown that S. agalactiae strain ST12 induces miRNAs associated with a stronger inflammatory response than S. agalactiae strain ST103. Both seem to drive the macrophages away from a classical (M1) type of activation, but further studies are needed to classify them into canonical activation phenotypes. Taken together these analyses suggest that the differentially expressed miRNAs identified in this study during in vitro infection of bovine macrophages with two different strains of S. agalactiae likely are crucial regulators of the innate immune response to this pathogen, and thus represent potential biomarkers of infection and inflammation. We have also identified several candidate pathways likely involved in the pathogenesis of subclinical mastitis. The study contributes to better understanding of the pathogenic mechanisms of different bacteria strains within one species.
Abbreviations
- AGO2:
-
Argonaute 2
- BAK1:
-
BCL-2 antagonist killer 1
- BAX:
-
BCL-2 associated X
- BCL2L1:
-
BCL-2 like 1
- BMEs:
-
Bovine mammary epithelial cells
- CCL11:
-
C-C motif chemokine ligand 11
- CCL26:
-
C-C motif chemokine ligand 26
- CFU:
-
Colony-forming units
- DE:
-
Differentially expressed
- DICER1:
-
Dicer 1, ribonuclease III
- GABA:
-
Gamma-aminobutyric acid
- GPCRs:
-
G-protein-coupled receptors
- IL-10:
-
Interleukin-10
- IL-12:
-
Interleukin-12
- IL-1β:
-
Interleukin-1 beta
- IL-4:
-
Interleukin-4
- IL-6:
-
Interleukin-6
- IL-8:
-
Interleukin-8
- IPA:
-
Ingenuity Pathway Analysis
- JaK:
-
Janus kinase
- LPS:
-
Lipopolysaccharides
- M1:
-
Classically activated macrophages
- miRNA:
-
microRNA
- miRNome:
-
Complete miRNA content
- MOI:
-
Multiplicity of infection
- NFKβ:
-
Nuclear factor kappa β
- NMBU:
-
Norwegian University of Life Sciences
- NR:
-
Norwegian Red
- OD:
-
Optical density
- PBMC:
-
Peripheral blood mononuclear cells
- PCA:
-
Principal component analysis
- PPIA:
-
Peptidylprolyl isomerase A
- RIN:
-
RNA integrity numbers
- RT-qPCR:
-
Reverse transcription-quantitative PCR
- SMAD2/3:
-
SMAD family member 2/3
- small nucleolar RNAs:
-
snoRNAs
- sRNA-seq:
-
Small-RNA-seq
- TGFβ1:
-
Transforming growth factor beta 1
- TLR:
-
Toll-like receptor
- TLR2:
-
Toll-like receptor 2
- TNFα:
-
Tumor necrosis factor alpha
Declarations
Acknowledgments
We would like to thank Anita Holm and Haakon Aaen from NMBU for collecting the blood samples, and Tormod Mørk, Bjørg Kvitle and Hanna Jørgensen from the Norwegian Veterinary Institute for providing the S. agalactiae isolates. We are also thankful to Rune Andreassen for technical comments on the study design.
Funding
This paper is a part of the “Multimast” project (project no. 233778) funded by the Research Council of Norway, Tine and Geno. The Research Council of Norway was not involved in the experimental design, data collection/interpretation, or preparation of this manuscript, however, representatives from Tine and Geno have been involved in the development of the grant application and discussion of the results.
Availability of data and materials
The datasets analyzed during the current study have been deposited in ArrayExpress database at EMBL-EBI (www.ebi.ac.uk/arrayexpress/) under accession number E-MTAB-5937.
Authors’ contributions
AMLS participated in the design of the study and the RT-qPCR experiment, carried out the infection experiments, analyzed the RT-qPCR and microRNA sequencing data, and drafted the manuscript. SFH participated in the design of the study, cultured the bacteria, carried out infection experiments and performed the RT-qPCR. PB, BH and IO participated in the study design, discussion and interpretation of the results and manuscript drafting. TRS and OØ participated in the study design and discussion of the results. All authors read and approved the final version of the manuscript.
Ethics approval and consent to participate
Blood sampling was conducted in accordance with the laws and regulations controlling experiments using live animals in Norway and the study was approved by the Norwegian Animal Research Authority (Norwegian Food Safety Authority; FOTS id: 8194).
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
Authors’ Affiliations
References
- Bi Y, Liu G, Yang R. MicroRNAs: novel regulators during the immune response. J Cell Physiol. 2009;218(3):467–72.View ArticlePubMedGoogle Scholar
- Flor TB, Blom B. Pathogens use and abuse MicroRNAs to deceive the immune system. Int J Mol Sci. 2016;17(4):538.View ArticlePubMedPubMed CentralGoogle Scholar
- Janssen HL, Reesink HW, Lawitz EJ, Zeuzem S, Rodriguez-Torres M, Patel K, et al. Treatment of HCV infection by targeting microRNA. N Engl J Med. 2013;368(18):1685–94.View ArticlePubMedGoogle Scholar
- Gigli I, Maizon DO. microRNAs and the mammary gland: a new understanding of gene expression. Genet Mol Biol. 2013;36(4):465–74.View ArticlePubMedPubMed CentralGoogle Scholar
- Lawless N, Reinhardt TA, Bryan K, Baker M, Pesch B, Zimmerman D, et al. MicroRNA regulation of bovine monocyte inflammatory and metabolic networks in an in vivo infection model. G3 (Bethesda). 2014;4(6):957–71.View ArticlePubMed CentralGoogle Scholar
- Lawless N, Foroushani AB, McCabe MS, O'Farrelly C, Lynn DJ. Next generation sequencing reveals the expression of a unique miRNA profile in response to a gram-positive bacterial infection. PLoS One. 2013;8(3):e57543.View ArticlePubMedPubMed CentralGoogle Scholar
- Jin W, Ibeagha-Awemu EM, Liang G, Beaudoin F, Zhao X, Guan le L. Transcriptome microRNA profiling of bovine mammary epithelial cells challenged with Escherichia coli or Staphylococcus aureus bacteria reveals pathogen directed microRNA expression profiles. BMC Genomics. 2014;15:181.View ArticlePubMedPubMed CentralGoogle Scholar
- Sun J, Aswath K, Schroeder SG, Lippolis JD, Reinhardt TA, Sonstegard TS. MicroRNA expression profiles of bovine milk exosomes in response to Staphylococcus aureus infection. BMC Genomics. 2015;16:806.View ArticlePubMedPubMed CentralGoogle Scholar
- Shaughnessy RG, Farrell D, Riepema K, Bakker D, Gordon SV. Analysis of Biobanked serum from a Mycobacterium avium subsp paratuberculosis bovine infection model confirms the remarkable stability of circulating miRNA profiles and defines a bovine serum miRNA repertoire. PLoS One. 2015;10(12):e0145089.View ArticlePubMedPubMed CentralGoogle Scholar
- Sica A, Erreni M, Allavena P, Porta C. Macrophage polarization in pathology. Cellular and molecular life sciences : CMLS. 2015;72(21):4111–26.View ArticlePubMedGoogle Scholar
- Duvel A, Frank C, Schnapper A, Schuberth HJ, Sipka A. Classically or alternatively activated bovine monocyte-derived macrophages in vitro do not resemble CD163/calprotectin biased macrophage populations in the teat. Innate immunity. 2012;18(6):886–96.View ArticlePubMedGoogle Scholar
- Jorgensen HJ, Nordstoga AB, Sviland S, Zadoks RN, Solverod L, Kvitle B, et al. Streptococcus agalactiae in the environment of bovine dairy herds--rewriting the textbooks? Vet Microbiol. 2016;184:64–72.View ArticlePubMedGoogle Scholar
- Manning SD, Springman AC, Lehotzky E, Lewis MA, Whittam TS, Davies HD. Multilocus sequence types associated with neonatal group B streptococcal sepsis and meningitis in Canada. J Clin Microbiol. 2009;47(4):1143–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Lewandowska-Sabat AM, Boman GM, Downing A, Talbot R, Storset AK, Olsaker I. The early phase transcriptome of bovine monocyte-derived macrophages infected with Staphylococcus aureus in vitro. BMC Genomics. 2013;14:891.View ArticlePubMedPubMed CentralGoogle Scholar
- Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, et al. Primer3--new capabilities and interfaces. Nucleic Acids Res. 2012;40(15):e115.View ArticlePubMedPubMed CentralGoogle Scholar
- Spalenza V, Girolami F, Bevilacqua C, Riondato F, Rasero R, Nebbia C, et al. Identification of internal control genes for quantitative expression analysis by real-time PCR in bovine peripheral lymphocytes. Vet J. 2011;189(3):278–83.View ArticlePubMedGoogle Scholar
- Piehler AP, Grimholt RM, Ovstebo R, Berg JP. Gene expression results in lipopolysaccharide-stimulated monocytes depend significantly on the choice of reference genes. BMC Immunol. 2010;11:21.View ArticlePubMedPubMed CentralGoogle Scholar
- Capece V, Garcia Vizcaino JC, Vidal R, Rahman RU, Pena Centeno T, Shomroni O, et al. Oasis: online analysis of small RNA deep sequencing data. Bioinformatics. 2015;31(13):2205–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29(1):15–21.View ArticlePubMedGoogle Scholar
- Liao Y, Smyth GK, featureCounts SW. An efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30(7):923–30.View ArticlePubMedGoogle Scholar
- Friedlander MR, Mackowiak SD, Li N, Chen W, Rajewsky N. miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 2012;40(1):37–52.View ArticlePubMedGoogle Scholar
- Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.View ArticlePubMedPubMed CentralGoogle Scholar
- Agarwal V, Bell GW, Nam JW, Bartel DP. Predicting effective microRNA target sites in mammalian mRNAs. elife. 2015;4Google Scholar
- Breuer K, Foroushani AK, Laird MR, Chen C, Sribnaia A, Lo R, et al. InnateDB: systems biology of innate immunity and beyond--recent updates and continuing curation. Nucleic Acids Res. 2013;41(Database issue):D1228–33.View ArticlePubMedGoogle Scholar
- O'Neill LA, Sheedy FJ, McCoy CE. MicroRNAs: the fine-tuners of toll-like receptor signalling. Nat Rev Immunol. 2011;11(3):163–75.View ArticlePubMedGoogle Scholar
- Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009;136(2):215–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Henneke P, Berner R. Interaction of neonatal phagocytes with group B streptococcus: recognition and response. Infect Immun. 2006;74(6):3085–95.View ArticlePubMedPubMed CentralGoogle Scholar
- Budd KE, Mitchell J, Keane OM. Lineage associated expression of virulence traits in bovine-adapted Staphylococcus aureus. Vet Microbiol. 2016;189:24–31.View ArticlePubMedGoogle Scholar
- Tassi R, McNeilly TN, Sipka A, Zadoks RN. Correlation of hypothetical virulence traits of two Streptococcus uberis strains with the clinical manifestation of bovine mastitis. Vet Res. 2015;46:123.View ArticlePubMedPubMed CentralGoogle Scholar
- Pang M, Sun L, He T, Bao H, Zhang L, Zhou Y, et al. Molecular and virulence characterization of highly prevalent Streptococcus agalactiae circulated in bovine dairy herds. Vet Res. 2017;48(1):65.View ArticlePubMedPubMed CentralGoogle Scholar
- Jiang H, Chen M, Li T, Liu H, Gong Y, Li M. Molecular characterization of Streptococcus agalactiae causing community- and hospital-acquired infections in shanghai, China. Front Microbiol. 2016;7:1308.PubMedPubMed CentralGoogle Scholar
- Korir ML, Laut C, Rogers LM, Plemmons JA, Aronoff DM, Manning SD. Differing mechanisms of surviving phagosomal stress among group B Streptococcus strains of varying genotypes. Virulence. 2017;8(6):924–37.View ArticlePubMedGoogle Scholar
- Curtale G, Mirolo M, Renzi TA, Rossato M, Bazzoni F, Locati M. Negative regulation of toll-like receptor 4 signaling by IL-10-dependent microRNA-146b. Proc Natl Acad Sci U S A. 2013;110(28):11499–504.View ArticlePubMedPubMed CentralGoogle Scholar
- Cobos Jimenez V, Willemsen AM, Bradley EJ, Baas F, van Kampen AH, Kootstra NA. Next-generation sequencing of microRNAs in primary human polarized macrophages. Genomics data. 2014;2:181–3.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang Y, Zhang M, Zhong M, Suo Q, Lv K. Expression profiles of miRNAs in polarized macrophages. Int J Mol Med. 2013;31(4):797–802.View ArticlePubMedGoogle Scholar
- Graff JW, Dickson AM, Clay G, McCaffrey AP, Wilson ME. Identifying functional microRNAs in macrophages with polarized phenotypes. J Biol Chem. 2012;287(26):21816–25.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhuang G, Meng C, Guo X, Cheruku PS, Shi L, Xu H, et al. A novel regulator of macrophage activation: miR-223 in obesity-associated adipose tissue inflammation. Circulation. 2012;125(23):2892–903.View ArticlePubMedGoogle Scholar
- Banerjee S, Cui H, Xie N, Tan Z, Yang S, Icyuz M, et al. miR-125a-5p regulates differential activation of macrophages and inflammation. J Biol Chem. 2013;288(49):35428–36.View ArticlePubMedPubMed CentralGoogle Scholar
- Furuse Y, Finethy R, Saka HA, Xet-Mull AM, Sisk DM, Smith KL, et al. Search for microRNAs expressed by intracellular bacterial pathogens in infected mammalian cells. PLoS One. 2014;9(9):e106434.View ArticlePubMedPubMed CentralGoogle Scholar
- Shmaryahu A, Carrasco M, Valenzuela PD. Prediction of bacterial microRNAs and possible targets in human cell transcriptome. J Microbiol. 2014;52(6):482–9.View ArticlePubMedGoogle Scholar
- Cullen BR. MicroRNAs as mediators of viral evasion of the immune system. Nat Immunol. 2013;14(3):205–10.View ArticlePubMedPubMed CentralGoogle Scholar
- Grundhoff A, Sullivan CS. Virus-encoded microRNAs. Virology. 2011;411(2):325–43.View ArticlePubMedPubMed CentralGoogle Scholar
- Weiberg A, Wang M, Lin FM, Zhao H, Zhang Z, Kaloshian I, et al. Fungal small RNAs suppress plant immunity by hijacking host RNA interference pathways. Science (New York, NY). 2013;342(6154):118–23.View ArticleGoogle Scholar
- Zhang C, Huys A, Thibault PA, Wilson JA. Requirements for human dicer and TRBP in microRNA-122 regulation of HCV translation and RNA abundance. Virology. 2012;433(2):479–88.View ArticlePubMedGoogle Scholar
- Schmitter D, Filkowski J, Sewer A, Pillai RS, Oakeley EJ, Zavolan M, et al. Effects of dicer and Argonaute down-regulation on mRNA levels in human HEK293 cells. Nucleic Acids Res. 2006;34(17):4801–15.View ArticlePubMedPubMed CentralGoogle Scholar
- Baer C, Squadrito ML, Laoui D, Thompson D, Hansen SK, Kiialainen A, et al. Suppression of microRNA activity amplifies IFN-gamma-induced macrophage activation and promotes anti-tumour immunity. Nat Cell Biol. 2016;18(7):790–802.View ArticlePubMedGoogle Scholar
- Ozeri V, Rosenshine I, Ben-Ze'Ev A, Bokoch GM, Jou TS, Hanski E. De novo formation of focal complex-like structures in host cells by invading streptococci. Mol Microbiol. 2001;41(3):561–73.View ArticlePubMedGoogle Scholar
- Agerer F, Lux S, Michel A, Rohde M, Ohlsen K, Hauck CR. Cellular invasion by Staphylococcus aureus reveals a functional link between focal adhesion kinase and cortactin in integrin-mediated internalisation. J Cell Sci. 2005;118(Pt 10):2189–200.View ArticlePubMedGoogle Scholar
- Plancon L, Du Merle L, Le Friec S, Gounon P, Jouve M, Guignot J, et al. Recognition of the cellular beta1-chain integrin by the bacterial AfaD invasin is implicated in the internalization of afa-expressing pathogenic Escherichia coli strains. Cell Microbiol. 2003;5(10):681–93.View ArticlePubMedGoogle Scholar
- Bhat R, Axtell R, Mitra A, Miranda M, Lock C, Tsien RW, et al. Inhibitory role for GABA in autoimmune inflammation. Proc Natl Acad Sci U S A. 2010;107(6):2580–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Reyes-Garcia MG, Hernandez-Hernandez F, Hernandez-Tellez B, Garcia-Tamayo F. GABA (a) receptor subunits RNA expression in mice peritoneal macrophages modulate their IL-6/IL-12 production. J Neuroimmunol. 2007;188(1–2):64–8.View ArticlePubMedGoogle Scholar
- Tian J, Lu Y, Zhang H, Chau CH, Dang HN, Kaufman DL. Gamma-aminobutyric acid inhibits T cell autoimmunity and the development of inflammatory responses in a mouse type 1 diabetes model. J Immunol. 2004;173(8):5298–304.View ArticlePubMedGoogle Scholar
- Alam S, Laughton DL, Walding A, Wolstenholme AJ. Human peripheral blood mononuclear cells express GABAA receptor subunits. Mol Immunol. 2006;43(9):1432–42.View ArticlePubMedGoogle Scholar
- Tian J, Chau C, Hales TG, Kaufman DL. GABA(a) receptors mediate inhibition of T cell responses. J Neuroimmunol. 1999;96(1):21–8.View ArticlePubMedGoogle Scholar
- Rajaram MV, Ni B, Morris JD, Brooks MN, Carlson TK, Bakthavachalu B, et al. Mycobacterium tuberculosis lipomannan blocks TNF biosynthesis by regulating macrophage MAPK-activated protein kinase 2 (MK2) and microRNA miR-125b. Proc Natl Acad Sci U S A. 2011;108(42):17408–13.View ArticlePubMedPubMed CentralGoogle Scholar
- Rothchild AC, Sissons JR, Shafiani S, Plaisier C, Min D, Mai D, et al. MiR-155-regulated molecular network orchestrates cell fate in the innate and adaptive immune response to Mycobacterium tuberculosis. Proc Natl Acad Sci U S A. 2016;113(41):E6172–e81.View ArticlePubMedPubMed CentralGoogle Scholar
- Verschoor CP, Dorrington MG, Novakowski KE, Kaiser J, Radford K, Nair P, et al. MicroRNA-155 is required for clearance of Streptococcus pneumoniae from the nasopharynx. Infect Immun. 2014;82(11):4824–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Aziz F. The emerging role of miR-223 as novel potential diagnostic and therapeutic target for inflammatory disorders. Cell Immunol. 2016;303:1–6.View ArticlePubMedGoogle Scholar
- Naeem A, Zhong K, Moisa SJ, Drackley JK, Moyes KM, Loor JJ. Bioinformatics analysis of microRNA and putative target genes in bovine mammary tissue infected with Streptococcus uberis. J Dairy Sci. 2012;95(11):6397–408.View ArticlePubMedGoogle Scholar
- Zhang T, Yu J, Zhang Y, Li L, Chen Y, Li D, et al. Salmonella enterica serovar enteritidis modulates intestinal epithelial miR-128 levels to decrease macrophage recruitment via macrophage colony-stimulating factor. J Infect Dis. 2014;209(12):2000–11.View ArticlePubMedPubMed CentralGoogle Scholar
- Luo K. Signaling cross talk between TGF-beta/Smad and other signaling pathways. Cold Spring Harb Perspect Biol. 2017;9:1.View ArticleGoogle Scholar
- Zhou M, Liu Z, Zhao Y, Ding Y, Liu H, Xi Y, et al. MicroRNA-125b confers the resistance of breast cancer cells to paclitaxel through suppression of pro-apoptotic Bcl-2 antagonist killer 1 (Bak1) expression. J Biol Chem. 2010;285(28):21496–507.View ArticlePubMedPubMed CentralGoogle Scholar
- Hotchkiss RS, Nicholson DW. Apoptosis and caspases regulate death and inflammation in sepsis. Nat Rev Immunol. 2006;6(11):813–22.View ArticlePubMedGoogle Scholar
- Newton K, Dixit VM. Signaling in innate immunity and inflammation. Cold Spring Harb Perspect Biol. 2012;4(3)Google Scholar
- Sun L, Ye RD. Role of G protein-coupled receptors in inflammation. Acta Pharmacol Sin. 2012;33(3):342–50.View ArticlePubMedPubMed CentralGoogle Scholar
- Li Z, Lan X, Han R, Wang J, Huang Y, Sun J, et al. miR-2478 inhibits TGFbeta1 expression by targeting the transcriptional activation region downstream of the TGFbeta1 promoter in dairy goats. Sci Rep. 2017;7:42627.View ArticlePubMedPubMed CentralGoogle Scholar
- Romao JM, Jin W, He M, McAllister T, Guan le L. MicroRNAs in bovine adipogenesis: genomic context, expression and function. BMC Genomics. 2014;15:137.View ArticlePubMedPubMed CentralGoogle Scholar
- Muroya S, Shibata M, Hayashi M, Oe M, Ojima K. Differences in circulating microRNAs between grazing and grain-fed wagyu cattle are associated with altered expression of intramuscular microRNA, the potential target PTEN, and Lipogenic genes. PLoS One. 2016;11(9):e0162496.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang R, Zhang YY, Lu JS, Xia BH, Yang ZX, Zhu XD, et al. The highly pathogenic H5N1 influenza a virus down-regulated several cellular MicroRNAs which target viral genome. J Cell Mol Med. 2017;21(11):3076–86.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang XG, Huang JM, Feng MY, Ju ZH, Wang CF, Yang GW, et al. Regulatory mutations in the A2M gene are involved in the mastitis susceptibility in dairy cows. Anim Genet. 2014;45(1):28–37.View ArticlePubMedGoogle Scholar
- Schurch NJ, Schofield P, Gierlinski M, Cole C, Sherstnev A, Singh V, et al. How many biological replicates are needed in an RNA-seq experiment and which differential expression tool should you use? RNA (New York, NY). 2016;22(6):839–51.View ArticleGoogle Scholar
- Vegh P, Foroushani AB, Magee DA, McCabe MS, Browne JA, Nalpas NC, et al. Profiling microRNA expression in bovine alveolar macrophages using RNA-seq. Vet Immunol Immunopathol. 2013;155(4):238–44.View ArticlePubMedGoogle Scholar
