- Research article
- Open Access
DNA methylation analysis of porcine mammary epithelial cells reveals differentially methylated loci associated with immune response against Escherichia coli challenge
BMC Genomics volume 20, Article number: 623 (2019)
Epigenetic changes such as cytosine (CpG) DNA methylations regulate gene expression patterns in response to environmental cues including infections. Microbial infections induce DNA methylations that play a potential role in modulating host-immune response. In the present study, we sought to determine DNA methylation changes induced by the mastitis causing Escherichia coli (E. coli) in porcine mammary epithelial cells (PMEC). Two time points (3 h and 24 h) were selected based on specific transcriptomic changes during the early and late immune responses, respectively.
DNA methylation analysis revealed 561 and 898 significant (P < 0.01) differentially methylated CpG sites at 3 h and 24 h after E. coli challenge in PMEC respectively. These CpG sites mapped to genes that have functional roles in innate and adaptive immune responses. Significantly, hypomethylated CpG sites were found in the promoter regions of immune response genes such as SDF4, SRXN1, CSF1 and CXCL14. The quantitative transcript estimation indicated higher expression associated with the DNA CpG methylation observed in these immune response genes. Further, E. coli challenge significantly reduced the expression levels of DNMT3a, a subtype of de novo DNA methylation enzyme, in PMEC indicating the probable reason for the hypomethylation observed in the immune response genes.
Our study revealed E. coli infection induced DNA methylation loci in the porcine genome. The differentially methylated CpGs were identified in the regulatory regions of genes that play important role in immune response. These results will help to understand epigenetic mechanisms for immune regulation during coliform mastitis in pigs.
Epigenomic changes are dynamically regulated by environmental cues. Among the various epigenetic modifications, cytosine (CpG) methylation of genomic DNA is an important reversible gene regulatory mechanism. DNA methylation plays a crucial role in transcriptional regulation by affecting the recruitment of regulatory factors onto promoters and enhancers . Microbial infections, where induced DNA methylation may modulate host immune responses [2, 3]. It has been shown that human papillomavirus infection cause aberrant host DNA methylation through direct interaction of viral protein E7 with DNA methyltransferase 1 (DNMT1) . Escherichia coli infection (UPEC) of human uroepithelial cells was reported to regulate the expression of DNMT1, and CpG hypermethylation down-regulated the cell-cycle inhibitor CDKN2A likely inhibiting apoptosis and increasing proliferation of uroepithelial cells . Infection by the protozoan, Leishmania donovani also alters DNA methylation profiles in human macrophages suppressing the host immune response and enabling intracellular survival of the protozoan . Similar observations were made for Mycobacterium tuberculosis (TB) infection that rapidly methylates host DNA at distal enhancer elements and associated chromatin remodelling . Schistosoma parasite infection induces hypermethylation in transcription factors that inhibit IFN-γ signalling in CD4+ T cells of children, with a significant effect on the downstream TB-specific immune phenotypes . These epigenetic effects of bacterial and parasitic infections and subsequent regulations of immune response highlight the role of DNA methylations in host-pathogen interactions.
Coliform mastitis (CM) causes postpartum dysagalactia syndrome (PDS), an important disease in pigs. The affected animals show high fever, loss of appetite, pain and inflammation of teats. PDS is a disease of economic significance as it severely affects the health and milk production of sows leading to poor survival of piglets . Gram-negative bacteria such as Escherichia coli (E. coli) are the most prominent causative pathogens isolated from PDS- affected sows. Lipopolysaccharide (LPS), an outer membrane endotoxin component of E. coli, is the major pathogenic factor that can induce inflammatory responses in sows with PDS . In our previous studies, an E. coli bacterial challenge of porcine mammary epithelial cells (PMECs) was used as a model for porcine mastitis at two time points (3 h and 24 h) that represented early and late transcriptional responses . There was a distinct transcription pattern at 3 h and 24 h post challenge of E. coli in porcine mammary epithelial cells and found upregulated set of genes, including of cytokines, chemokines, and cell adhesion factors, which together coordinate the immune response of host cells . These time points broadly represented early and late transcriptional responses. These changes are mediated through epigenetic mechanism including microRNA . We hypothesized that the E. coli induced transcriptomic changes in PMEC also follow other epigenomic modifications. Therefore, our present study focuses on host-cell DNA methylation changes induced by the mastitis causing E. coli.
Genome-wide DNA methylation in PMEC
Genome-wide DNA methylation profiling of PMEC (unchallenged control, 3 h post-challenge (hpc) E. coli and 24hpc E. coli) generated 25–40 million Illumina sequencing reads for each of the nine reduced-representation bisulfite sequencing (RRBS) libraries. Mapping to the pig genome (Sscrofa 11.1) using Bismark (Bowtie 2) revealed that overall 50% of the generated reads uniquely mapped to the genome. Twenty percent of the CpGs analysed across all samples mapped to the 5′ regulatory region (promoter) of genes with functional annotations. Further, approximately 45% of mapped CpGs were enriched in known CpG islands of the pig genome. The average CpG methylation levels was approximately 45% and the non-CpG level was less than 10% in both the control and treated groups (Table 1, Additional file 1: Figure S1). The analysis of methylation levels of the CpG sites in different regions of the genome indicated that CpGs located in the upstream promoter regions of genes had lower methylation levels (< 10%); whereas, gene bodies and intragenic spaces had higher methylation levels (> 30%). In total, 49,921 CpG positions which for common across all the control and treated samples were identified for further study after quality checks and normalization with at least 10X coverage. These methylation levels were similar and no significant differences were observed across different samples. Further, the methylation levels of the CpGs presented a similar bimodal distribution in all the samples of control and treatment, which is consistent with the results of earlier studies on distribution of CpG methylation (Additional file 2: Figure S2).
Differentially methylated CpGs in E. coli challenged PMEC
After quality checks and normalization with at least 10X coverage by using pairwise analysis, 82,833, 132,624 and 116,632 CpG positions common between control and E. coli 3 hpc group, between control and E. coli 24 hpc group and between E. coli 3 hpc group and 24 hpc group, respectively, were further used for differential methylation analysis. There were a total of 561 differentially methylated CpGs when comparing the E. coli 3 hpc group and the unchallenged control group (P < 0.01, Additional file 6). A total of 898 differentially methylated CpGs when comparing the E. coli 24 hpc group and unchallenged control group (P < 0.01, Additional file 7). The analysis between 3 hpc and 24 hpc groups indicated 855 differentially methylated CpGs loci (P < 0.01, Additional file 8) and their similar distribution in different genomic features.
The Venn diagram shown in Fig. 1a indicates the numbers of differentially methylated loci at both comparing (3 vs control and 24 hpc vs control). The proportion of these loci present in different genomic regions (promoter, exon, intron or intergenic) is shown in Fig. 1b. Manhattan plots were generated to show the chromosome-wide distribution of differentially methylated CpG sites identified in this study (Fig. 2). Statistically significant CpG sites are indicated as being above the cut off negative P values (3) on the log 10 scale. The chromosome wide distributions of CpG sites for the comparisons groups (3 hpc vs. control, 24 hpc vs. control) were shown (Fig. 2a & b). Similarly, the volcano plots depicted for the significant CpG sites based on both P values and methylation differences of 20% between E. coli challenged and control PMEC the groups (3 hpc vs. control, 24 hpc vs. control) (Fig. 3a & b).
The three sample groups of PMEC (Control, E. coli 3 hpc & E. coli 24 hpc) were examined by cluster analysis using hierarchical clustering method. The heatmap shows differential methylation of the top 100 CpG sites between different groups (Fig. 4). The top 10 significant CpG sites possibly differentially methylated at the in early phase of infection were selected by pair-wise comparison between E. coli 3 hpc and control samples. These differentially methylated mapped to genes including CYC1 (cg605583), GDNF (cg22964888), PTHR20855_SF21 (cg3759161), MYO1B (cg95955381), GNG7 (cg76083105) and LARP1. Similarly, the top ten differentially methylated loci at the late phase of infection were identified by pairwise comparison of E. coli 24 hpc and control samples. These differentially methylated loci mapped to genes including CYB5R1 (cg24937951), COLA4 (cg77052099), SDF4 (cg63545568), NTN1 (cg54274409), DNAH17 (cg3518148), GSE1 (cg3323937), MMP20 (cg33280374), ABCA4 (cg123320105) and SF3B2 (cg6288374) (Table 2). Figure 5 displays representative differentially methylated loci and their associated genes that exhibited similar methylation patterns (increased or decreased methylation) at both treatment time points (3 hpc and 24 hpc).
Functional analysis of differentially methylated genes
The annotated genes associated with differentially methylated CpG sites were submitted to functional analysis. This revealed the molecular functions, cellular processes, and biological disorders associated with genes identified during the early (3 hpc) and late (24 hpc) phases of E. coli challenge. The top 4 molecular and cellular functions at 3 hpc were cellular function and maintenances, cell death & survival, cell morphology, cellular assembly and compromise. Similar functions were found at 24 hpc in addition to post-translation modification, cellular growth, movement and proliferation. The diseases most associated with these genes were inflammatory diseases, organismal injury & abnormalities, and immunological diseases (Table 3). By comparing 3 hpc and 24 hpc groups, the top 4 molecular and cellular functions were cellular compromise, cell morphology, cell cycle and
cell-to-cell signaling and interaction. No inflammatory diseases or immunological diseases genes was found by comparing 3 hpc and 24 hpc groups. Further, differentially methylated CpGs that are located within the promoter regions (TSS ± 2000 bp) were grouped by k-means clustering and the affected top biological functions were identified. The innate immune regulation and stress activated mechanisms were identified in cluster groups at 3 h and 24 h respectively (Additional file 3: Figure S3). Additional k-means clusters of the union of differentially methylated CpGs of all three groups the most important biological functions of the group were performed (Additional file 4: Figure S4).
The most enriched transcription factors were identified by considering genes with differentially methylated CpGs within their promoter regions and using DAVID with the Homo sapiens UCSC TFBS function, a function not yet supported in Sus scrofa spp. the results revealed the enriched transcription factors identified by comparing E. coli 3 hpc to control samples (Table 4). One such transcription factor, PAX5, is known to influence B-cell differentiation and trafficking by regulating large numbers of downstream genes involved in immune functions. Other enriched factors with interesting roles in immune response include MSX1 (innate immunity), CREB (cAMP signalling), and IRF2 (Interferon Regulatory factor-2). Similarly, E. coli 24 hpc vs unchallenged analysis revealed enrichment of transcription factors (Table 4). Among these PAX5, CREB, and AP4 were again found to be enriched along with other transcription factors such as XBP1 and E2F that are involved in immune response and other biological functions.
The genome wide analysis identified differentially methylated CpG sites mapped to different genes. Based on these results, SENP6 (cg90300054), SRXN1 (cg34587223), JAK2 (cg217001619), AQP2 (cg62156997) and ZMYM2 (651169) were selected for bisulfite PCR and pyrosequencing analysis. The results validated the methylation patterns of CpG sites for SRXN1, SENP6, and JAK2 that were hypomethylated in E. coli after 24 h challenged cells compared to the control cells. However, the pyrosequencing validations were not significant for CpG sites in AQP2, and ZMYM2 genes (Fig. 6 and Additional file 5: Figure S5).
Expression analysis using quantitative RT-PCR
The differentially methylated CpG loci within promoter or regulatory regions were deemed quite likely to influence the expression of their downstream genes. We selected genes with the significant hypomethylated CpG sites in the promoter (− 2 kb to + 1 kb from TSS) as identified in our study for quantitative expression analysis. Significantly higher expression was noticed for stromal derived factor 4 (SDF4), sulforedoxin (SRXN1), Colony stimulating factor 1 (CSF1), and CXCL14. The zinc finger transcription factor (ZMYM2) gene exhibited no change in expression between when comparing control and E. coli challenged PMEC (Fig. 7). The expression of DNA methyl transferases (DNMTs) was also examined. There was a significant reduction in the expression of DNMT3a (de novo methyl transferase) in 3 hpc cells, however no major changes were observed for other related genes (DNMT1 and DNMT3b).
E. coli infection is the main source of mastitis in domestic animals, and epigenetic mechanisms have a strongly influence on the expression of genes during E. coli and Staphylococcus aureus (S. aureus) associated mastitis [13,14,15]. The LPS endotoxin, a major E. coli-derived pathogenic factor for inflammatory response in mastitis, induces changes in innate immunity genes in dairy cows, and a role for DNA methylation has been demonstrated for this process . Exposure of bovine dermal fibroblasts to a demethylating agent, 5-aza-2′-deoxycitidine (AZA), altered their response to an LPS challenge, indicating the role of DNA methylation in immune response . Also, a recent genome wide methylation analysis revealed DNA methylation differences contributing to LPS-induced immune response in bovine fibroblasts .
Pigs are an important disease model of coliform mastitis, and the present study is the first to investigate the underlying epigenetic mechanisms mediating host-pathogen interaction in cultured porcine mammary cells. To achieve this, we have taken advantage of the clear differences in gene expression patterns occurring during the immediate (3 h) and late immune response (24 h) observed in our earlier studies [11, 12] where we demonstrated that pathogenic E. coli challenge of PMEC resulted in up-regulation of immune response genes, such as chemokines, cytokines, and cell adhesion molecules.
In the current study, we examined changes in DNA methylation induced by E. coli challenge. The genome wide methylation patterns indicated comparatively lower CpG methylation in upstream promoter regions compared to gene body regions, and a bimodal distribution of CpG methylation levels was observed. These findings are consistent with prior studies on genome wide CpG methylation in pig [19, 20]. Our analysis revealed 561 and 898 differentially methylated CpG loci associated with the immediate and late responses to an E. coli challenge, respectively at P < 0.01. Annotation of these CpG probes identified differentially methylated genes, and pathway analysis revealed that the highest enriched functional associations of these genes were inflammatory diseases and organismal injury. In terms of molecular and cellular functions, the top differentially methylated probes associated with cell death & survival, post-translation modifications, cell morphology, cell growth & proliferations. Further differentially methylated CpG’s were clustered into different groups of related genes that indicated immune response as one of the major biological process affected including other cellular processes. Promoter regions are known to have binding sites for different transcription factors that drive gene expression. The genes with differentially methylated CpG’s in their promoters were used to identify relevant transcription factors, and several such identified factors have known prominent roles in inflammation, including PAX5, AP4, CREB, IRF2 and XBP1. PAX5 is a transcription factor that regulates various B cell functions including activation of NF-kB . IRF2 (Interferon regulatory factor 2) regulates IFN-β expression and was found to inhibit LPS induced proinflammatory responses . CREB (cAMP-responsive element binding protein) transcriptionally activates pro-inflammatory genes, such as IL-2, IL-6, and TNF-α by binding to the cAMP response element (CRE) in their promoter sequences . XBP1 (X-box binding protein 1) is a positive regulator of TLR gene induction and plays a major role in interlinking LPS-associated TLR activation to Endoplasmic Reticulum (ER) stress .
In general, it is known that hypomethylation of CpG motifs in the promoter regions of genes enhances their transcription. As shown in earlier studies, infection induced hypomethylation of CpGs leads to higher expression of immune genes [7, 25]. In our study, we examined the expression of SDF4, SRXN1, CSF1, CXCL14 and ZMYM2, all genes which we found to have a significant challenge-dependent reduction in CpG methylation in their upstream regulatory regions around TSS. SDF4, SRXN1, and CXCL14 were selected from hypomethylation CpGs at 24 h after E. coli-challenged, while CSF1 were selected from hypomethylation CpGs at 3 h after E. coli-challenged. The quantitative transcript estimation showed that only the expression of SDF4 and CSF1 in E.coli-challenged cells at 24 h and 3 h was significantly higher than in control cells. However expression pattern and methylation level of SRXN1, and CXCL14 didn’t show a negative relationship. Many studies reported different methylation sites associated with expression, regardless of the directional change in expression and methylation level . This may because DNA methylation is not exclusively associated with repression of transcription initiation . The expression of ZMYM2 was also higher in the challenged cells, although this difference was not significant. These results confirmed that the CpG methylation in the promoter regions of these genes correlates with changes in gene expression.
The present study revealed promoter hypomethylation of SDF4 and significantly higher expression of SDF4 in E.coli-challenged PMEC. SDF (stromal derived factors) belong to the CXC subfamily of proteins, are expressed in different tissues, and function as chemokines to attract inflammatory cells . Increased levels of SDF1 (CXCL12) can be seen in inflammatory conditions, such as sub-acromial bursitis and acute liver injury [29, 30]. SDF4 expression was found to be significantly reduced in mammary tumours compared to normal tissue and low levels of SDF4 are linked to a poor prognosis . Our results indicate a possible epigenetic mechanism for regulating SDF expression and immune response in mammary tissues.
In this study, the upstream regulatory region of the SRXN1 gene had reduced CpG methylation and significantly increased expression in E. coli-challenged cells. Sulfiredoxin (SRXN1) is an antioxidant enzyme that prevents ROS injury to cells, can reduces oxidized cysteine residues of peroxiredoxin proteins (Prx I- IV), and facilitate removal of free radicals . LPS released form E. coli induces reactive oxygen species (ROS) formation leading to upregulation of antioxidant enzymes, including SRXN1, and SRXN1 production in macrophages was shown to protect mice from LPS induced endotoxic shock . Earlier studies indicate that the promoter region of the SRXN1 gene has putative NF-kB and AP1 transcription factor binding sites required for up-regulation of SRXN1 expression by LPS treatment . Hypomethylation may enhance the binding of these transcription factors at the promoter region to drive the SRXN1 expression.
In the present study, we found increased expression of CSF1 in the immediate early phase of E.coli challenge but not at the later phase. Colony stimulating factor-1 (CSF-1) has been shown to differentiate peripheral blood monocytes into tissue macrophages [35, 36]. CSF-1 expression is induced by uropathogenic E.coli infection and has a critical role in bacterial clearance during infection [37, 38]. The higher expression of CSF1 corresponds to hypomethylation in the upstream regulatory regions of the CSF1 gene. Curiously, expression of the CSF1 receptor (CSF1R) is also known to be regulated by DNA methylation of its promoter region . CXCL14 is an additional gene found to be regulated in this study. CXCL14 is a chemokine that plays inflammatory modulator and host defence roles in epithelial tissues [40, 41]. It was shown to have antimicrobial activity in the respiratory tract and contributed to the clearance of Streptococcus pneumoniae . There is evidence for epigenetic regulation of CXCL14 in prostate cancer cells. Treatment of these cells with the demethylating agent 5-aza-2-deoxycytidine affects a hypermethylated CpG island in the CXCL14 gene promoter resulting in the recovery of CXCL14 expression and restoration of chemotaxis . Methylation-mediated control of CXCL14 can also be inferred from our study where an E.coli challenge induced hypomethylation in the promoter region of the CXCL14 gene associated with higher levels of expression. These results indicate the potential involvement of epigenetic mechanisms in regulating host cell response to E. coli infection.
Some of the genomic regions showed different methylation between infection groups, including CpGs at SSC5 39.78–40.07 Mb (OSBPL8) and at SSC9 63.38–63.39 Mb (MROH9). These regions (9:63272406–63401079 bp and 5: 39774063–39828561 bp) contained large CpG islands (CGI) (genome order: Sscrofa11.1). These regions also showed heterogeneity in methylation and changes in the degree of methylation between infection groups. But both of these regions are far away from the promoter site of the transcripts.
DNA methyltransferases (DNMT) are the principle enzymes responsible for controlling epigenetic modifications, and DNMT3a and DNMT3b are responsible for de novo DNA methylation. In the present study, the quantitative expression results indicated significantly reduced expression of DNMT3a in the E. coli infected PMEC compared to the control group; however, no changes in the patterns of either DNMT1 or DNMT3B were observed. Hypomethylation in upstream regions of immunity genes may partially be explained by the reduced expression of DNMT3a. It was previously shown that UPEC infections increase DNMT1 expression, the enzyme responsible for maintenance DNA methylation . The level of various DNA methyl transferase enzymes is likely to modulate the expression of many genes during bacterial infection.
In conclusion, the present study identified first time genome wide differential CpG methylation patterns induced by E. coli challenge in PMEC. CpG methylation changes in the upstream regulatory regions were used to identify enriched transcription factors that regulate immune response pathways. Further, reduced DNA CpG methylation was observed in the immune response genes with corresponding increases in their expression. These results indicate potential epigenetic mechanisms that regulate inflammation during coliform mastitis in pigs.
Primary culture of PMEC
German Law of Animal Protection guidelines were followed for collecting the tissues. Animal Care Committee at Leibniz Institute of Farm Animal Biology (FBN), Dummerstorf 18196, Germany approved the experiments. The sows were weighed and slaughtered by electronarcosis followed by exsanguination in the experimental abattoir of the FBN. Tissues from mammary complexes cranial of the navel were collected aseptically immediately after slaughter from each individual. After tissue collection all animals underwent routine processes of the slaughterhouse. Veterinary inspection of the animals before slaughtering and of their carcasses and organs after slaughter proofed that they were without any impairment, disease symptoms and pathological signs. Primary cultures of PMEC were obtained as described in our earlier studies . Briefly, tissues were isolated from mammary glands of three lactating sows, minced and digested with collagenase enzyme. The cells were washed and suspended in complete growth medium. Primary cultures of PMECs were established by removing other cell types, such as fibroblasts and adipocytes, by selective trypsinization (Trypsin/EDTA-0.25/0.02%, Sigma-Aldrich). The fibroblasts and adipocytes detach more rapidly; whereas, epithelial cell islands remain adhered to the surface of the culture flasks. The procedure was repeated several times until a uniform and confluent monolayer of epithelial cells was obtained.
E. coli challenge to the cultured PMEC
The present study uses the same E. coli strain (gMEc240, sequence type 101, phylogroup B1, C+) isolated from milk of PDS-positive sows described in our previous study [11, 12]. Briefly, approximately 4.4 × 105 PMEC from each 3 sow (3 biological replicates) were seeded and cultured in collagen-coated 6-well plates (1:10 collagen R in distilled water, Menal, Emmendingen, Germany) in complete medium without APS (three technical replicates per individual and treatment condition). After 24 h, the medium was changed. Forty-eight hours after seeding, the cells reached ~ 90% confluency. Then, PMEC were challenged with 107/ml heat-inactivated E. coli for 3 h or for 24 h. The control PMEC cells were not challenged with E. coli. The medium was discarded, and the cells washed three times with phosphate buffered saline (PBS, PAA) to remove bacteria after incubation periods. The experiments included in triplicates for each three animals in three groups (control, 3 h and 24 h). In total 27 genomic DNA and total RNA samples were isolated from treated and control PMEC.
Preparation of RRBS libraries
Equivalent amount of genomic DNA from three technical replicates per individual animal were pooled. In total 9 pooled samples of DNA, three for each E. coli challenged PMEC at two different time points (3 h post-challenge (hpc) and 24 hpc) and unchallenged control. PMEC were used for libraries construction. RRBS libraries were prepared using 2 μg of pooled genomic DNA with a 1% spike-in control (unmethylated cl857 Sam Lambda DNA (Promega)). The genomic DNA was digested with Msp I and TaqαI. Double-enzyme (MspI and TaqαI) digestion RRBS with increased size-selected fragments will enhance genome-wide CpG coverage. The digested fragments were end repaired, A-tailed and ligated with the C-methylated adaptor sequences TruSeq Nano DNA Sample Preparation kit (Illumina) by following the manufacturer’s protocol (Illumina, San Diego, CA). The DNA fragments were later size selected for 40–200 bp with a 2.5% NuSieve 3:1 agarose gel and extracted using the Zyomclean™ Gel DNA Recovery Kit (Zymo Research). The purified DNA was treated with bisulfite using the EZ DNA Methylation-Gold Kit™ (Zymo Research). The preparative scale PCR was performed for 15 cycles and PCR products were purified with the DNA Clean and Concentrator Kit™ (Zymo Research). The qualities of the RRBS libraries were assessed using an Agilent DNA 1000 kit (Agilent Technologies). NGS of the 9 RRBS libraries were performed on an Illumina HiSeq2500 for single-reads of 114 bp at the FBN, Dummerstorf. The bcl2fastq2 conversion software v2.19 was used to convert base call (BCL) files from a sequencing run into FASTQ files that were used for further analysis.
The sequence reads were assessed for quality using FastQC and bases with Phred score greater than 20 were retained for further downstream analysis. RRBS introduces artificial CpG at the 3′ end of the fragments that are removed along with the adaptor sequences to avoid their inclusion in the methylation calling. The default settings for Trim Galore (v0.1.1.1, Babraham Bioinformatics, UK) were used for Illumina adaptor trimming as they specifically remove the first two bases from the 3′ end of the sequence such that the additional ‘C’ closest to the enzyme cut site is removed. The trimmed reads were mapped to the in-silico bisulfite converted porcine genome (11.1) using the Bismark alignment tool (v0.13.1, Babraham Bioinformatics, UK). Bisulfite treatment converts unmethylated cytosines to uracils whereas, methylated cytosine is not affected. The sequence reads were mapped to the pre-converted reference genome (Sscrofa 11.1), reads aligned to the multiple regions are removed and best uniquely mapped reads were used for methylation calling.
Differential methylation analysis and annotation
Methyl call files from the Bismark aligner with the percent methylation score per base were taken as input files for analysis. The reads which cover all the treatment and control samples with a minimum coverage of 10 were only considered to increase the power of statistical test. The reads from the sex chromosomes, mitochondria, unannotated genome segments and those showing no methylation variation across all samples were filtered out. Differential methylation analysis was done using the MethylKit . Logistic regression with “MN” basic overdispersion correction was applied for testing methylation proportion of each CpG between the treatment and control group samples. The methylation differences between the groups were considered to be statistically significant at P < 0.01. A heatmap was used to show the methylation differences between the groups using selected DNA methylation loci. The differentially methylated CpG were annotated to genomic features by using the genomation R/Bioconductor package.
Bisulfite PCR and pyrosequencing
Differentially methylated CpGs identified by genome-wide analysis were validated using bisulfite PCR and pyrosequencing methods. The same genomic DNA which was used for genome wide methylation analysis was treated with bisulfite using EZ DNA Methylation-Gold Kit™ (Zymo Research). Primers were designed using the Pyrosequencing Assay Design Software (Additional file 9) and target region was amplified with PCR using AmpliTaq Gold DNA Polymerase (Applied Biosystems, Cat. No. 4311814). As most of the differentially methylated loci are in CpGislands, in addition to target CpGs, the additional adjacent CpGs were also included for pyrosequencing. Pyrosequencing was performed with PSQ™96MA as per the manufacturer’s instructions. The CpG methylation percentages were calculated using PyroMark CpG software 2.0 based on the height of the T and C peaks at the target CpG sites.
Functional network analysis was done to gain biological insights into top differentially methylated loci between E. coli challenged and control PMEC. Genes annotated from the selected CpG with P < 0.01 were included in the gene function network analysis. Ingenuity pathway analysis (Ingenuity Systems, Inc., CA, USA) with its core analysis features was used. Differentially methylated CpGs present between − 2.5 kb and + 1 kb from an annotated transcription start site (TSS) were considered as being in the promoter region of a gene. Such promoter gene IDs were used for transcription factor enrichment, and transcription factor binding site (TFBS) analysis was done using the default parameters of the UCSC_TFBS track of Homo sapiens due to the lack of a Sus scrofa data track and considering that most of the TFBS are conserved.
Quantitative gene expression analysis
Expression analysis was done for selected genes with differentially methylated regulatory regions. RNA was isolated from 27 samples after the E. coli challenged and control PMEC using the TRI reagent (Sigma-Aldrich) following the manufacturer’s instructions. Isolated RNA was purified by the RNeasy Mini Kit (Qiagen) and DNase I treatment was done to remove the contaminating DNA. First strand cDNA was synthesized using SuperScript III MMLV reverse transcriptase (Invitrogen) with 1 μg of RNA, 500 ng oligo (dT) and 500 ng random hexamer primers (Promega). Quantitative real-time PCR was performed using the LightCycler® 480 Real-Time PCR System (Roche Diagnostics). GAPDH and RPL32 were used as internal housekeeping control genes. The sequences of primers for the selected test genes and internal control genes were designed using Primer3 (v.0.4.1) (Additional file 9). The reaction conditions were: 95 °C for 5 min (initial denaturation), 45 cycles consisting of 60 °C for 15 s (denaturation), 72 °C for 25 s (extension/fluorescence acquisition). Melting curve analysis and agarose gel electrophoresis indicated single amplicon without any primer dimers. Threshold cycle values were normalized to the expression of the internal control genes, and significant differences in expression were assessed with ANOVA and P < 0.05 was considered statistically significant.
Availability of data and materials
All RRBS sequencing data have been deposited in the ArrayExpress database at EMBL-EBI (www.ebi.ac.uk/arrayexpress) under the accession number E-MTAB-7363.
Postpartum dysagalactia syndrome
Porcine mammary epithelial cells
Reduced-representation bisulfite sequencing
Jones PA. Functions of DNA methylation: islands, start sites, gene bodies and beyond. Nat Rev Genet. 2012;13(7):484.
De Monerri NCS, Kim K. Pathogens hijack the epigenome: a new twist on host-pathogen interactions. Am J Pathol. 2014;184(4):897–911.
McGrath-Morrow S, Ndeh R, Helmin KA, Chen S-Y, Anekalla KR, Abdala-Valencia H, D'Alessio FR, Collaco JM, Singer BD. DNA methylation regulates the neonatal CD4+ T-cell response to pneumonia in mice. J Biol Chem. 2018;293(30):11772–83.
Laurson J, Khan S, Chung R, Cross K, Raj K. Epigenetic repression of E-cadherin by human papillomavirus 16 E7 protein. Carcinogenesis. 2010;31(5):918–26.
Tolg C, Sabha N, Cortese R, Panchal T, Ahsan A, Soliman A, Aitken KJ, Petronis A, Bägli DJ. Uropathogenic E. coli infection provokes epigenetic downregulation of CDKN2A (p16INK4A) in uroepithelial cells. Lab Invest. 2011;91(6):825.
Marr AK, MacIsaac JL, Jiang R, Airo AM, Kobor MS, McMaster WR. Leishmania donovani infection causes distinct epigenetic DNA methylation changes in host macrophages. PLoS Pathog. 2014;10(10):e1004419.
Pacis A, Tailleux L, Morin AM, Lambourne J, MacIsaac JL, Yotova V, Dumaine A, Danckaert A, Luca F, Grenier J-C. Bacterial infection remodels the DNA methylation landscape of human dendritic cells. Genome Res. 2015.
DiNardo AR, Nishiguchi T, Mace EM, Rajapakshe K, Mtetwa G, Kay A, Maphalala G, Secor WE, Mejia R, Orange JS. Schistosomiasis Induces Persistent DNA Methylation and Tuberculosis-Specific Immune Changes. J Immunol. 2018;201(1):124–33.
Niemi JK, Bergman P, Ovaska S, Sevón-Aimonen M-L, Heinonen M. Modeling the costs of postpartum Dysgalactia syndrome and locomotory disorders on sow productivity and replacement. Front vet sci. 2017;4:181.
Kaiser M, Jacobson M, Andersen PH, Bækbo P, Cerón JJ, Dahl J, Escribano D, Jacobsen S. Inflammatory markers before and after farrowing in healthy sows and in sows affected with postpartum dysgalactia syndrome. BMC Vet Res. 2018;14(1):83.
Jaeger A, Bardehle D, Oster M, Günther J, Muráni E, Ponsuksili S, Wimmers K, Kemper N. Gene expression profiling of porcine mammary epithelial cells after challenge with Escherichia coli and Staphylococcus aureus in vitro. Vet Res. 2015;46(1):50.
Jaeger A, Hadlich F, Kemper N, Lübke-Becker A, Muráni E, Wimmers K, Ponsuksili S. MicroRNA expression profiling of porcine mammary epithelial cells after challenge with Escherichia coli in vitro. BMC Genomics. 2017;18(1):660.
Vanselow J, Yang W, Herrmann J, Zerbe H, Schuberth H-J, Petzl W, Tomek W, Seyfert H-M. DNA-remethylation around a STAT5-binding enhancer in the αS1-casein promoter is associated with abrupt shutdown of αS1-casein synthesis during acute mastitis. J Mol Endocrinol. 2006;37(3):463–77.
Modak R, Mitra SD, Vasudevan M, Krishnamoorthy P, Kumar M, Bhat AV, Bhuvana M, Ghosh SK, Shome BR, Kundu TK. Epigenetic response in mice mastitis: role of histone H3 acetylation and microRNA (s) in the regulation of host inflammatory gene expression during Staphylococcus aureus infection. Clin Epigenetics. 2014;6(1):12.
Song M, He Y, Zhou H, Zhang Y, Li X, Yu Y. Combined analysis of DNA methylome and transcriptome reveal novel candidate genes with susceptibility to bovine Staphylococcus aureus subclinical mastitis. Sci Rep. 2016;6:29390.
Doherty R, O’Farrelly C, Meade K. Epigenetic regulation of the innate immune response to LPS in bovine peripheral blood mononuclear cells (PBMC). Vet Immunol Immunopathol. 2013;154(3–4):102–10.
Green BB, Kerr DE. Epigenetic contribution to individual variation in response to lipopolysaccharide in bovine dermal fibroblasts. Vet Immunol Immunopathol. 2014;157(1–2):49–58.
Korkmaz FT, Kerr DE. Genome-wide methylation analysis reveals differentially methylated loci that are associated with an age-dependent increase in bovine fibroblast response to LPS. BMC Genomics. 2017;18(1):405.
Choi M, Lee J, Le MT, Nguyen DT, Park S, Soundrarajan N, Schachtschneider KM, Kim J, Park J-K, Kim J-H. Genome-wide analysis of DNA methylation in pigs using reduced representation bisulfite sequencing. DNA Res. 2015;22(5):343–55.
Schachtschneider KM, Madsen O, Park C, Rund LA, Groenen MA, Schook LB. Adult porcine genome-wide DNA methylation patterns support pigs as a biomedical model. BMC Genomics. 2015;16(1):743.
Wang D, Chen J, Li R, Wu G, Sun Z, Wang Z, Zhai Z, Fang F, Guo Y, Zhong Y. PAX5 interacts with RIP2 to promote NF-κB activation and drug-resistance in B-lymphoproliferative disorders. J Cell Sci. 2016;129(11):2261–72.
Cui H, Banerjee S, Guo S, Xie N, Liu G. IFN Regulatory Factor 2 Inhibits Expression of Glycolytic Genes and Lipopolysaccharide-Induced Proinflammatory Responses in Macrophages. J Immunol. 2018;200(9):3218–30.
Wen AY, Sakamoto KM, Miller LS. The role of the transcription factor CREB in immune function. J Immunol. 2010;185(11):6413–9.
Martinon F, Chen X, Lee A-H, Glimcher LH. TLR activation of the transcription factor XBP1 regulates innate immune responses in macrophages. Nat Immunol. 2010;11(5):411.
Gomes AVBT, de Souza Morais SM, Menezes-Filho SL, de Almeida LGN, Rocha RP, Ferreira JMS, Dos Santos LL, Malaquias LCC, Coelho LFL. Demethylation profile of the TNF-α promoter gene is associated with high expression of this cytokine in dengue virus patients. J Med Virol. 2016;88(8):1297–302.
Schachtschneider KM, Liu Y, Rund LA, Madsen O, Johnson RW, Groenen MA, Schook LB. Impact of neonatal iron deficiency on hippocampal DNA methylation and gene transcription in a porcine biomedical model of cognitive development. BMC Genomics. 2016;17(1):856.
Spruijt CG, Vermeulen M. DNA methylation: old dog, new tricks? Nat Struct Mol Biol. 2014;21:949–54.
Hall JM, Korach KS. Stromal cell-derived factor 1, a novel target of estrogen receptor action, mediates the mitogenic effects of estradiol in ovarian and breast cancer cells. Mol Endocrinol. 2003;17(5):792–803.
Yang F, Sun W, Yang Y, Wang Y, Li C-L, Fu H, Wang X-L, Yang F, He T, Chen J. SDF1–CXCR4 signaling contributes to persistent pain and hypersensitivity via regulating excitability of primary nociceptive neurons: involvement of ERK-dependent Nav1. 8 up-regulation. J Neuroinflammation. 2015;12(1):219.
Chalin A, Lefevre B, Devisme C, Pronier C, Carrière V, Thibault V, Amiot L, Samson M. Serum CXCL10, CXCL11, CXCL12, and CXCL14 chemokine patterns in patients with acute liver injury. Cytokine. 2018;111:500–4.
Kang H, Escudero-Esparza A, Douglas-Jones A, Mansel RE, Jiang WG. Transcript analyses of stromal cell derived factors (SDFs): SDF-2, SDF-4 and SDF-5 reveal a different pattern of expression and prognostic association in human breast cancer. Int J Oncol. 2009;35(1):205–11.
Jeong W, Bae SH, Toledano MB, Rhee SG. Role of sulfiredoxin as a regulator of peroxiredoxin function and regulation of its expression. Free Radic Biol Med. 2012;53(3):447–56.
Planson A-G, Palais G, Abbas K, Gerard M, Couvelard L, Delaunay A, Baulande S, Drapier J-C, Toledano MB. Sulfiredoxin protects mice from lipopolysaccharide-induced endotoxic shock. Antioxid Redox Signal. 2011;14(11):2071–80.
Kim K-H, Jeong J-Y, Surh Y-J, Kim K-W. Expression of stress-response ATF3 is mediated by Nrf2 in astrocytes. Nucleic Acids Res. 2009;38(1):48–59.
Sordet O, Rébé C, Plenchette S, Zermati Y, Hermine O, Vainchenker W, Garrido C, Solary E, Dubrez-Daloz L. Specific involvement of caspases in the differentiation of monocytes into macrophages. Blood. 2002;100(13):4446–53.
Boulakirba S, Pfeifer A, Mhaidly R, Obba S, Goulard M, Schmitt T, Chaintreuil P, Calleja A, Furstoss N, Orange F. IL-34 and CSF-1 display an equivalent macrophage differentiation ability but a different polarization potential. Sci Rep. 2018;8(1):256.
Ingersoll MA, Kline KA, Nielsen HV, Hultgren SJ. G-CSF induction early in uropathogenic Escherichia coli infection of the urinary tract modulates host immunity. Cell Microbiol. 2008;10(12):2568–78.
Mavromatis C, Bokil NJ, Totsika M, Kakkanat A, Schaale K, Cannistraci CV, Ryu T, Beatson SA, Ulett GC, Schembri MA. The co-transcriptome of uropathogenic E scherichia coli-infected mouse macrophages reveals new insights into host–pathogen interactions. Cell Microbiol. 2015;17(5):730–46.
Zasłona Z, Scruggs AM, Peters-Golden M, Huang SK. Protein kinase a inhibition of macrophage maturation is accompanied by an increase in DNA methylation of the colony-stimulating factor 1 receptor gene. Immunology. 2016;149(2):225–37.
Lu J, Chatterjee M, Schmid H, Beck S, Gawaz M. CXCL14 as an emerging immune and inflammatory modulator. J Inflamm. 2016;13(1):1.
Collins PJ, McCully ML, Martínez-Muñoz L, Santiago C, Wheeldon J, Caucheteux S, Thelen S, Cecchinato V, Laufer JM, Purvanov V. Epithelial chemokine CXCL14 synergizes with CXCL12 via allosteric modulation of CXCR4. FASEB J. 2017;31(7):3084–97.
Dai C, Basilico P, Cremona TP, Collins P, Moser B, Benarafa C, Wolf M. CXCL14 displays antimicrobial activity against respiratory tract bacteria and contributes to clearance of Streptococcus pneumoniae pulmonary infection. J Immunol. 2015;1402634.
Song EY, Shurin MR, Tourkova IL, Gutkin DW, Shurin GV. Epigenetic mechanisms of promigratory chemokine CXCL14 regulation in human prostate cancer cells. Cancer Res. 2010;70(11):4394–401.
Akalin A, Kormaksson M, Li S, Garrett-Bakelman FE, Figueroa ME, Melnick A, Mason CE. methylKit: a comprehensive R package for the analysis of genome-wide DNA methylation profiles. Genome Biol. 2012;13(10):R87.
The authors thank Joana Bittner, Nicole Gentz and Annette Jugert for excellent technical assistance.
The Leibniz Institute for Farm Animal Biology (FBN) provided own funding. The funding body plays no role in the design, collection, analysis, interpretation of data and in writing this manuscript.
Ethics approval and consent to participate
Animal care and tissue collection procedures followed the guidelines of the German Law of Animal Protection and the experimental protocol was approved by the Animal Care Committee of the Leibniz Institute for Farm Animal Biology as well as by the State Mecklenburg-Western Pomerania (Landesamt für Landwirtschaft, Lebensmittelsicherheit und Fischerei; LALLF M-V/TSD/7221.3–2.1-020/09). All animals used in this study belong to the Leibniz Institute for Farm Animal Biology (FBN). The Animal Care Committee of the FBN has approved this procedure. The experimental protocol was carried out in accordance with the approved guidelines for safeguarding good scientific practice at the institutions in the Leibniz Association and the measures were taken to minimize pain and discomfort and accord with the guidelines laid down by the European Communities Council Directive of 24 November 1986 (86/609/EEC).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Genomic distribution of mapped CpG residues from RRBS libraries to the known functional annotations of porcine RefSeq genes in-terms of reads fraction (A) and percentages (B). The genomic distribution of RRBS reads to the porcine genomic CpG island/CpG shore regions in terms of reads fraction (C) and percentages (D). (TIF 855 kb)
Figure S2. Methylation levels of identified CpG sites. The bimodal distribution of CpG methylation was observed in all the samples (A). The methylation levels (%) at different genomic features such as CpG islands, CpG shores (B) and at Promotes, Exons, Introns and Intergenic regions (C) represented. (TIF 1017 kb)
Figure S3. k-means clustering of differentially methylated genes (CpG in TSS ± 2000) with k = 2 and scaled as Z-score across rows. Top five gene ontology (GO) biological processes derived from each k-means clusters ranked based on the fold enrichment. A) E coli 3 hpc vs control, B) E coli 24 hpc vs control. (TIF 2836 kb)
Figure S4. k-means clustering of differentially methylated common CpG that were present in both E. coli 3 hpc vs control and E. coli 24 hpc vs control and top five enriched biological biological process. (TIF 1258 kb)
Figure S5. Differentially methylated CpG sites identified between E. coli 3hpc or 24hpc compared to the unchallenged control group from NGS data compare to pyrosequencing including SENP6 (cg90300054), SDF4 (cg63545568), JAK2 (cg217001619), SRXN1 (cg34587223), ZMYM2 (c651169). The y-axis for both box plots represents methylation level. Genes associated with the CpG are given. Box plot represents the range of variation and median value. (TIF 2466 kb)
Differentially methylated CpG sites in E. coli 3hpc compared to the unchallenged control group. (XLSX 55 kb)
Differentially methylated CpG sites in E. coli 24hpc compared to the unchallenged control group. (XLSX 75 kb)
Differentially methylated CpG sites in E. coli 3hpc compared to E. coli 24hpc. (XLSX 71 kb)
List of primers sequences used for quantitative gene expression and pyrosequencing. (DOCX 15 kb)