Global DNA methylation changes and differential gene expression in Anaplasma phagocytophilum-infected human neutrophils
© Sinclair et al. 2015
Received: 26 March 2015
Accepted: 30 June 2015
Published: 29 July 2015
Anaplasma phagocytophilum is an obligate intracellular prokaryotic pathogen that both infects and replicates within human neutrophils. The bacterium represses multiple antimicrobial functions while simultaneously increasing proinflammatory functions by reprogramming the neutrophil genome. Previous reports show that many observed phenotypic changes are in part explained by altered gene transcription. We recently identified that large chromosomal regions of the neutrophil genome are differentially expressed during A. phagocytophilum infection. Because of this, we sought to determine whether gene expression programs altered by infection were the result of changes in the host neutrophil DNA methylome.
Within 24 h of infection, marked increases in DNA methylation were observed genome-wide as compared with mock-infected controls and pharmacologic inhibition of DNA methyltransferases resulted in decreased bacterial growth. New regions of DNA methylation were enriched at intron and exon junctions; however, intragenic methylation did not correlate with altered gene expression. In contrast, intergenic DNA methylation was associated with A. phagocytophilum-induced gene expression changes. Within the major histocompatibility complex locus on chromosome 6, a region with marked changes in infection-induced differential gene expression, new regions of methylation were localized to boundaries of active and inactive chromatin.
These data strongly suggest that A. phagocytophilum infection, in addition to altering histone structure, alters DNA methylation and the epigenome of its host cell to promote survival and replication, providing evidence that such bacterial infection can radically alter the epigenome of its host cell.
KeywordsDNA methylation Anaplasma phagocytophilum Hypermethylation Cellular reprogramming MHC locus
Anaplasma phagocytophilum is an obligate intracellular prokaryotic pathogen of mammalian neutrophils. Transmitted by Ixodes spp. ticks, A. phagocytophilum infection results in manifestations ranging from subclinical infection or mild self-limited fever to severe infection requiring hospitalization or causing death [1, 2]. The neutrophil is an unusual host cell for any microorganism, especially bacteria, due to its role in innate immunity where its primary function is to recognize and kill microbes. Despite this inhospitable environment, A. phagocytophilum reprograms its host cell so that the bacterium can survive long enough for continued replication. The major phenotypic alterations among infected human neutrophils include complex functions such as decreased antimicrobial activity and respiratory burst, reduced phagocytosis, reduced margination and emigration across the endothelium, delayed apoptosis, and an increased production and release of proinflammatory cytokines, chemokines, and proteases [3–8]. These functional alterations allow the bacterium to survive long enough to replicate and spread to newly recruited host cells and to ensure that its host cell remains viable within the intravascular network until it can be accessed by a feeding tick.
The observed changes in host cell function are partially explained by alterations in granulocyte gene transcription [9–11]. Decreased respiratory burst allows for prolonged inhibition of superoxide production favoring bacterial survival and likely results from the downregulation of two key components of the NADPH oxidase required for proper oxidase assembly: CYBB, which encodes gp91 phox , and RAC2, a GTPase [3, 12]. Likewise, transcription of BCL2 family member genes are maintained during A. phagocytophilum infection which results in delayed neutrophil apoptosis [7, 13] and increased transcript levels of cytokines and chemokines such as IL-1α and IL-8 contribute to exaggerated inflammatory responses that recruit new neutrophil hosts [14–16].
Precisely how A. phagocytophilum coordinates the transcriptional changes that allow it to reprogram host cell functions is not known but is unlikely to result from an accumulation of individual signaling pathways targeted by unique pathogen protein effectors. We described decreased histone acetylation and methylation at multiple defense gene promoters with infection  and enhanced intracellular bacterial growth with HDAC-1 and −2 expression and activity, as well as the dose-dependent decrease in bacterial load effects of HDAC1 silencing and HDAC pharmacologic inhibition . We also showed that transfection of the nucleomodulin and type IV system secretion substrate AnkA of A. phagocytophilum alone mimics many transcriptional changes observed with infection by directly binding host cell DNA in multiple genomic regions, including at the CYBB proximal promoter where it decreases histone H3 acetylation and gene expression . Previously, we showed that the A. phagocytophilum nucleomodulin AnkA recruits HDAC1 to the CYBB promoter, and the altered chromatin configuration excludes RNA polymerase II recruitment thereby silencing expression . These observations prove an epigenetic basis for silencing at CYBB and suggest that A. phagocytophilum infection alters transcriptional programs of its host cells via genome-wide epigenetic alterations.
HDACs target nucleosomal histone proteins to create modifications that alter gene promoter accessibility to transcriptional regulators. Interestingly, HDACs also complex with DNA methyl binding proteins (MBDs) and work to coordinately alter transcriptional programs. DNA methylation at CpG dinucleotides is an important epigenetic regulator of cellular differentiation, neoplasia, and metastasis [19–22]. Moreover, expression of DNA methyltransferase 3A (DNMT3A), a key enzyme for de novo DNA methylation, is upregulated in A. phagocytophilum-infected neutrophils ex vivo . Given the complex interplay between these two epigenetic modifiers and the induction of histone deacetylation by A. phagocytophilum, we hypothesized that A. phagocytophilum also alters the DNA methylome of human neutrophils that in part determines transcriptional reprogramming.
Precedence for bacterial influence on host DNA methylation patterns has been established in both Helicobacter pylori and uropathogenic Escherichia coli (UPEC) infections . Gastric mucosa infected with H. pylori shows increased DNA methylation patterns at CpG islands but is hypothesized to be a result of infection-related inflammation rather than bacterial-induced DNA methylation (Na and Woo, 2014). It is important to note that H. pylori is generally considered to be an extracellular organism so its effects on gastric mucosa are not due to direct invasion and replication inside the stomach tissue. UPEC, which does actively invade uroepithelial cells, was shown to directly influence the methylation of target genes, but its genome-wide effects have not yet been elucidated . Here, we found that A. phagocytophilum induces hypermethylation of the human neutrophil genome within 24 h of infection. Although several bacteria are implicated in alterations of host cell DNA methylation , this is the first demonstration that infection by an intracellular bacterium alters host DNA methylation patterns on a genome-wide scale. The data strongly suggest that this is part of a coordinated reprogramming of host cell functions which leads to improved microbial fitness by promoting intracellular survival, replication, and spread.
A. phagocytophilum induces genome-wide hypermethylation
Inhibition of DNMTs with 5-azacytidine slows A. phagocytophilum growth
Characterization of DNA methylation across common gene features
DNA methylation and transcriptional alterations in the A. phagocytophilum-infected neutrophil
We then sought to determine whether or not DNA methylation influenced transcription of the most differentially expressed genes. Genes that were up- or downregulated with a fold-change of 3–6 or >6 standard deviations (SD) from the mean fold-change were assessed for DNA methylation peaks. Since DNA methylation is found to be correlated across genomic regions at a scale of <2 kb , if a region of DNA methylation was located within the gene body or ±3 kb from TSS or TTS, it was considered associated with that gene. Using these criteria, we compared the fold-change in expression of 516 genes associated with DNA methylation to the remaining 404 genes not associated with DNA methylation. The average log2 differential gene expression (DGE) of 361 methylated genes with between 3–6 SD from the mean was −0.037, and no statistical association with DNA methylation could be identified when compared with 269 genes lacking DNA methylation (p = 0.565). At this level, 120 of 361 genes with DNA methylation were upregulated and 100 of 269 genes lacking DNA methylation were upregulated. In contrast, of the 155 methylated genes that were differentially expressed ≥6 SD above the mean, 86 were upregulated (log2 DGE 0.694) rather than downregulated, although the resulting association with DNA methylation was not significant compared to 135 non-methylated genes with the same differential gene expression. This is consistent with a previously published report by Borjesson et al. that the most differentially expressed genes induced by A. phagocytophilum infection were upregulated .
Although A. phagocytophilum infection induced genome hypermethylation, the native level of genome methylation in neutrophils was high. Paired with the relatively low fold-changes in gene expression and the small proportion of all genes that were differentially regulated, these data could mask correlations that would provide evidence for the effects that A. phagocytophilum has on DNA methylation that in turn lead to altered gene expression. Because of this, we set out to identify local effects at the chromosomal level especially given that large regions of chromosomes are regulated similarly during infection .
Correlation of individual new meDNA and differential gene transcription at the MHC locus
The major histocompatibility complex (MHC) locus on chromosome 6 was shown to have local clustering of increased gene expression across a large chromosomal area spanning nearly 4 Mb . This locus also has well-described epigenetic regulation; thus, we focused our attention on the association of gene expression and DNA methylation here. To examine whether increased DNA methylation at the MHC locus on chromosome 6 was related to the observed cluster of differentially regulated genes, we examined newly methylated DNA (meDNA) fold-enrichment at individual genes and created a sliding scale window for both the log2 fold-change in expression and newly methylated sites in those windows.
On chromosome 6, there were 687 unique newly methylated genomic sites for which differential gene expression data was also obtained, including 53 unique newly methylated sites in the MHC locus and 22 unique newly methylated sites in the position immediately downstream of MHC that we call “PROX”. When individual newly methylated sites across the entire chromosome were examined, no correlation between individual differential gene expression and meDNA fold-enrichment (ρ −0.0391, p = 0.306) or –log10 p value (ρ −0.0174 and p = 0.649) could be identified. Although correlation between individual differential gene expression with enrichment of individual newly methylated sites over the MHC locus was higher (ρ 0.2664, p = 0.053 and –log10 p value ρ 0.2311; p = 0.096), it did not reach a level of significance we defined. Likewise, there was no significant correlation at the downstream control PROX region (meDNA fold-enrichment ρ 0.0864, p = 0.702, −log10 p value ρ −0.0212 and p = 0.925).
Correlation of new meDNA enrichment with differential gene transcription over long contiguous linear regions at the MHC locus and chromosome 6
To deconstruct this further, we subdivided individual DNA methylation peaks and differentially expressed genes over chromosome 6 into their gene regions: intron, CDS, intergenic, and promoter. Here, no single gene region of chromosome 6 had new DNA methylation fold-enrichment that significantly correlated, either positively or negatively, with differential gene transcription (Spearman p ≥ 0.091), although the number of meDNA peaks in many categories was relatively small. Only 16 peaks for meDNA enrichment were found in 9 distinct gene promoters (BMP6, NOL7, FRS3, PTK7, RCAN2, GJA1, SYNJ2, TULP4, and PLG), yielding a Spearman ρ value of −0.438 (p = 0.091). Because those peaks were located outside of gene bodies, they were grouped with intergenic regions. This grouping increased the intergenic region Spearman coefficient to −0.017, suggesting that although these 16 regions could influence expression of their respective genes, the likely impact across all genes on chromosome 6 is not high. There was insufficient data information about whether and how DNA methylation at other genomic regions affects gene expression since similar findings were obtained when the windows corresponding to these regions were evaluated. Here, the greatest correlation was in the 5′ UTR regions, yet the ρ value was −0.223, with a p value of 0.247. The window regions corresponding to the MHC locus could not be evaluated in the same manner because of the paucity of data over windows in this region. Yet in specific genomic loci, the impact of newly meDNA regions is important, especially considering that newly meDNA was present in 47 % of intergenic and promoter regions across chromosome 6 but in 58 % of the intergenic regions in the MHC locus.
Alteration of DNA methylation patterns is well documented in stem cell differentiation, cellular reprogramming, cancers, and other disease states but not in bacterial infections [21, 22, 31–34]. Previously, several groups suggested that infection by microbial pathogens induces hypermethylation of candidate genes involved in oncogenesis [35, 36, 24]. Helicobacter pylori infection causes chronic inflammation of surrounding gastric tissue resulting in aberrant DNA methylation patterns which can contribute to neoplastic transformation . Uropathogenic E. coli can induce hypermethylation of CDKN2A, a cyclin-dependent kinase inhibitor and tumor suppressor gene, leading to increased proliferation of the urothelium . Most recently, methylomes of ocular adnexal extranodal marginal zone B-cell lymphoma (EMZL) tumors from patients with Chlamydophila psittaci infection were compared to EMZL C. psittaci-negative tumors . Differentially methylated regions of infected samples hierarchically clustered together suggesting changes in DNA methylation patterns were specific to C. psittaci infection . Though informative, two of these studies investigated tissue samples which inevitably contain heterogeneous cell populations, and the other only investigated several candidate genes. Here, we show a direct link between infection of A. phagocytophilum and genome-wide alterations of DNA methylation profiles of the human neutrophil.
A. phagocytophilum infection induces widespread hypermethylation of the neutrophil genome within 24 h of intracellular infection. Hypermethylation of the neutrophil genome benefits the bacterium and allows for increased colonization of the host cell since inhibition of DNMTs by 5-azacytidine resulted in decreased bacterial loads at 24-72 h. This further suggests that even small changes in DNA methylation affect bacterial growth. Previous studies that investigated neutrophil methylation patterns focused on changes induced by myeloid differentiation and regions of hypomethylation [33, 34]. Despite the hypermethylation of the host genome induced by A. phagocytophilum, many of the key regions identified as being hypomethylated remained unmethylated yet had reduced expression in infected cells, suggesting a level of regulation beyond DNA methylation in these regions.
Regions of methylation were found to localize within gene bodies and gene-associated DNA within 3 kb of the transcriptional start and termination sites. Intragenic DNA methylation is associated with regulation of alternative promoters, alternative splicing, and non-coding RNAs (ncRNAs) [38–42]. Genes with intragenic DNA methylation continue to be expressed or have increased expression. It is thought that CpG methylation in intragenic regions results in decreased transcription from alternate promoters, and thus, decreased transcription of regulating ncRNAs that in turn promotes increased target gene transcription [39, 42]. In contrast, Jjingo et al. suggest that increased transcription and gene body methylation could be a result of DNA accessibility to methylating complexes, rather than regulation by ncRNA or expression of alternate promoters . Previously published reports regarding neutrophil transcription in response to A. phagocytophilum infection show that more genes are upregulated than downregulated [9–11], and this observation is borne out by the demonstration that among neutrophil genes with newly methylated DNA after A. phagocytophilum infection, more with the most significant differential transcription are upregulated than downregulated; however, a demonstration of which genomic features, when DNA hypermethylated, influence differential gene expression could not be deduced by this approach. The chromosome 6 MHC locus was one of the regions with the most differentially expressed genes at 24 h in these datasets, and these findings correlate with new regions of new DNA methylation. This discovery, in light of genome-wide hypermethylation, further supports the hypothesis that epigenetic alterations coordinate infected granulocyte transcriptional reprogramming. Whether the result of ncRNAs, the effects of methylating complexes, or other mechanisms, hypermethylation of intragenic sequences could support marked changes in transcriptional programs, including those that benefit A. phagocytophilum survival in neutrophils. Owing to the suboptimal infection levels obtained in the Borjesson et al. studies, the proximate cause of these changes will require further investigation .
Closer examination of new methylation sites in infected samples shows higher than expected localization within exons and intron/exon junctions. Rönnerblad et al. independently confirmed a large number of differentially methylated regions in human neutrophils localized after the first exon in gene bodies . DNA methylation at intron/exon junctions is associated with increased H3K36 trimethylation resulting in alternative splice variants . Currently it is not known whether DNA methylation is the precipitating factor for histone H3K36 methylation or vice versa. Given these findings, we hypothesize that deep RNA-seq of A. phagocytophilum-infected neutrophils could identify splice variants of highly expressed genes. The localization of newly methylated regions in gene bodies suggests that transcription of these genes could be enhanced and that splice variants are likely to exist; these observations could shed light on some aspects of the transcriptional reprogramming that is key to A. phagocytophilum survival and transmission. Previous reports utilizing transcriptional profiling by microarrays [9, 10] are unlikely to detect or quantify the majority of potential variants.
Despite the marked increase in gene-associated DNA methylation induced by A. phagocytophilum, there is not a strong genome-wide correlation between infection-induced alterations in gene transcription and annotated regions of meDNA, especially at individual gene or DNA methylation marks. Intragenic DNA methylation seems to have no effect on expression on a genome-wide scale. This could suggest that most newly methylated regions were not transcriptionally active and that the potential for upregulated expression is abrogated by A. phagocytophilum infection. We now have evidence to suggest that intergenic DNA methylation influences transcriptional programs of nearby genes as observed across the genome as a whole. For example, within the MHC locus on chromosome 6, regions significantly enriched for meDNA appear to demarcate boundaries of chromatin with increased expression and localize within regions of differential expression. This particular arrangement suggests that DNA methylation, and potentially histone marks, co-contribute to alterations over large chromosomal territories that influence expression of surrounding genes and accompanying transcriptional programs. Based on these observations, DNA methylation or factors that influence DNA methylation and DGE across large genomic regions could be important aspects of A. phagocytophilum neutrophil infection, and thus, disease pathogenesis.
The conceptual framework regarding genomic regions of differential methylation has scarcely been examined. However, recent studies in colonic carcinoma reveal that focal regions, >100 kb in size, correlate with regions of attachment to the nuclear lamina, suggesting that factors that bring together distant coordinately regulated genomic regions, such as those tethered by matrix attachment region binding proteins, could guide DNA methylation, histone acetylation, and perhaps ncRNA expression to direct transcriptional reprogramming and pathogen survival . Much remains to be determined, and an important question includes whether A. phagocytophilum plays a passive or active role in the process. A compelling hypothesis could revolve around the known translocation of the A. phagocytophilum nucleomodulin AnkA that binds to DNA in multiple genomic locations and recruits histone-modifying complexes including those known to interact with DNA methyltransferases and methylated CpG-binding domain proteins [17, 18]. The specificity of AnkA-DNA binding could then dictate genome-wide epigenetic alterations that reprogram cells to support A. phagocytophilum survival instead of microbial killing. Indeed, whether such reprogramming events occur as a result of other prokaryotic nucleomodulins is an area of intense research by investigators of microbial pathogenesis. The epigenetic research tools pioneered by oncologists, stem cell biologists, and those who study reprogramming of somatic cells now find applications in microbial pathogenesis studies as well. With mounting evidence of the impact of nucleomodulins from Shigella, Listeria, Legionella, Chlamydia, Mycobacterium, Toxoplasma, and other pathogenic microbes comes the realization that their impact surmounts the classical appreciation of influences over signaling pathways and intracellular trafficking, to the control of entire cellular programs.
Within 24 h of infection, A. phagocytophilum induces hypermethylation of neutrophil DNA on a genome-wide scale. Hypermethylation of the genome is associated with enhanced bacterial growth and colonization of the host cell. DNA methylation was predominantly associated with gene bodies or within 3 kb of transcriptional start and termination sites and tended to localize to exons and intron/exon junctions. Intragenic DNA methylation did not correlate strongly with A. phagocytophilum-induced differential gene expression over all chromosomes. However, intergenic DNA methylation is associated with infection-related transcriptional alterations. Targeted investigation of the MHC locus on chromosome 6 demonstrated that newly methylated sites fall within inactive chromatin domains and line the boundaries of active chromatin and that DNA methylation is associated with differential gene expression. Together, these data suggest that A. phagocytophilum-induced hypermethylation of the neutrophil genome is important for initial survival and propagation and that its influence over cellular epigenetics can lead to cellular reprogramming of even terminally-differentiated neutrophils.
Cell lines and cell culture
Primary peripheral blood neutrophils were isolated from venous blood of three healthy adult donors (two females, one male) as approved by the Johns Hopkins Medicine IRB and as previously described . Briefly, EDTA-anticoagulated blood was dextran-sedimented and leukocyte-rich plasma centrifuged through a Ficoll-Paque gradient. Mononuclear cells were removed and the remaining erythrocytes were lysed in hypotonic saline. HL-60 promyelocytic cells (ATCC CCL-240) were purchased from American Type Tissue Culture (Manassas, VA, USA). Both neutrophils and HL-60 cells were maintained or propagated in RMPI 1640 (Hyclone, Thermo Fisher Scientific, Waltham, MA, USA) with 10 % fetal bovine serum (Thermo Fisher Scientific) and Glutamax (Life Technologies, Carlsbad, CA, USA). HL-60 cells were differentiated for 3 to 5 days prior to infection with 1 μM all-trans retinoic acid (ATRA). All cells were grown at 37 °C in a humidified incubator with 5 % CO2. Where stated, HL-60 cells were treated with 5-azacytidine (Sigma-Aldrich, St. Louis, MO, USA) in DMSO for 24 h and separately with RG108 (Sigma-Aldrich, St. Louis, MO, USA) in DMSO for 24, 48, and 72 h, or with DMSO alone as a vehicle control.
Anaplasma phagocytophilum propagation, neutrophil infection, and quantitation
A. phagocytophilum was maintained in HL-60 cells as previously described . Cell-free bacteria were obtained from HL-60 cells of which 90–95 % were infected. Infected cells were centrifuged at 500× g for 5 min and resuspended in 1× PBS. Infected cells were lysed by sonication using an output of four of a Branson Sonifier 250 (Branson Ultrasonics, Danbury, CT, USA) for 15 s or via syringe lysis using a 22 g syringe. Lysed cells were then centrifuged at 1000× g for 10 min, and the bacteria-enriched supernatant was centrifuged at 13,000× g for 30 min. The bacterial pellets were suspended in RMPI 1640 medium.
Neutrophils were infected ex vivo in triplicate using a multiplicity of infection (MOI) of 25:1 and 100:1 for experiments using HL-60 cells. An estimate of 10 infectious bacteria per infected HL-60 cell was used. The proportion of infected cells was determined at 24 h after Romanowsky staining (Protocol HEMA3, Thermo Fisher Scientific) of cytocentrifuged cells. In addition and for experiments to test the effect of DNA methyltransferase inhibitors on A. phagocytophilum propagation, DNA extracted from these cultures (DNeasy Blood and Tissue Kit, Qiagen, Valencia, CA, USA) was used in a quantitative real time 5′ nuclease PCR assay that targets A. phagocytophilum msp2/p44 normalized to human ACTB . A standard curve of cloned msp2 was used to accurately estimate the number of bacteria per cell, providing that the A. phagocytophilum genome encodes approximately 100 msp2/p44 genes . For experiments that examined the effects of DNA methyltransferase inhibition on A. phagocytophilum growth, 8-azacytosine and RG108 were initially diluted in DMSO, then diluted to working concentrations (200, 100, and 50 nM for 8-azacytosine and 100, 50, and 25 μg/mL for RG-108) in RPMI 1640 medium. DMSO diluted to equivalent concentrations as used in the 200 nM 8-azacytosine or 100 μg/mL RG-108 preparations was used as vehicle controls.
MBD-Seq library preparation
Enrichment of 5-methylcytosine modified DNA was performed using DNA extracted from ex vivo infected and uninfected neutrophils from there healthy donor subjects 24 h post ex vivo infection. DNA libraries isolated by MBD enrichment were analyzed by MBD-seq to determine sites of MBD protein binding and DNA methylation (SKCCC Next Generation Sequencing Core). MBD-seq was carried out as described previously . Briefly, 2 μg of DNA was sonicated to an average molecular size of ~150–250 bp and end-repaired using the NEBNext SOLiD DNA library preparation kit end-repair module following the manufacturer’s protocol (New England Biolabs, Ipswich, MA, USA). Fragments were purified using a Qiagen PCR purification kit (Qiagen). SOLiD P1 and P2 adapters lacking 5′ phosphate groups (Life Technologies) were ligated to the fragment using the NEBNext adapter ligation module, column-purified, and subjected to isothermal nick-translation by treating with Platinum Taq polymerase to remove the nick. The resulting library was divided into two fractions, a total input fraction, and an enriched methylated fraction. The enriched methylated fraction was subjected to affinity enrichment of methylated DNA fragments using recombinant C-terminal 6xHis-tagged MBD2-MBD polypeptides immobilized on magnetic beads, similar to previously described methods [47–49]. The resulting enriched methylated fraction and the total input fraction were then subjected to library amplification using the NEBNext amplification module according to the manufacturer’s protocols, using 4–6 cycles for the total input, and 10–12 cycles for the enriched methylated fractions. Library fragments between 200–300 bp were selected after agarose gel electrophoresis. The libraries were then subjected to emulsion PCR and bead enrichment following the SOLiD emulsion PCR protocol (Life Technologies). The resulting beads were then deposited on the SOLiD flow cell and subjected to massively parallel 50 bp single-read sequencing on a Life Technologies SOLiD sequencer. Reads were mapped to the human genome hg18 build using Bioscope software. The total number of reads per sample and the percentage alignment are summarized in Additional file 3: Table S1. Regions of methylated DNA enrichment were identified using MACS  and SICER  as described previously . Additional file 4: Table S2 summarizes the number of peaks called with both MACS and SICER as well as delineates the total number of peaks above various thresholds of −10*log10 (p value) >50, >100, and >200. These and all subsequent DNA methylation studies were conducted independently using the three individual donors’ neutrophils, and all data is shown as a single donor representative of all three, or as an average of all three donors’ results, as noted in the figures.
Luminometric Methylation Assay and 5-methylcytosine DNA ELISA
LUMA was performed by EpigenDx (Hopkinton, MA, USA) as described by Karimi et al. . Briefly, samples were digested with HpaII + EcoRI or MspI + EcoRI, and fragments were amplified by Pyrosequencing™ (Qiagen). The ratio of (dGTP + dCTP)/dATP was calculated for HpaII/EcoRI and MspI/EcoRI fragments, and the percentage of methylation was determined by the ratio of (HpaII/EcoRI)/(MspI/EcoRI). Samples were tested in duplicate. 5-Methylcytosine was quantified in DNA preparations from RG108-treated cultures at days 0–3 postinfection using a 5-mC DNA ELISA (Zymo Research Corp., Irvine, CA, USA). The quantity of 5-methyl cytosine DNA was quantified by comparison with a standard curve and normalized to the content of day 0 no drug control cultures.
Identifying gene regions enriched for DNA methylation
Publically available cis-regulatory element annotation (CEAS) software [52, 53] was used to characterize regions of methylated DNA enrichment by calculating the fraction present in different gene regions (introns, exons, 5′ UTR, 3′UTR, and distal intergenic regions), promoters, bi-directional promoters, and regions downstream of gene bodies; this was also used to create average profile analyses of the smoothed adjusted log2 (M/T) values across transcriptional start sites (TSS), transcriptional termination sites (TTS), long (2715–11,673 bp), medium (842–2715 bp), and short (158–842 bp) intronic sequences and long (164–483 bp), medium (109–164 bp), and short (66–109 bp) exonic sequences.
Identifying chromosomal locations enriched for DNA methylation and differential gene expression
Publically available transcription microarray data from A. phagocytophilum-infected neutrophils was re-analyzed using RMA . Only the 24 h time point was used for comparison with the same interval in the DNA methylation studies. Differential expression was determined based on the standard deviation of the fold-change from the mean. Infected and uninfected fold-change gene expression data were transformed to log2 for all analyses.
Regions enriched for DNA methylation and differential expression were organized linearly among their respective chromosomes and divided into bins of ~one million base pairs. The number of gene features (expression) and methylation peaks within each bin were counted and compared to the average number across the genome to determine locations of interest. Regions of interest were defined by those bins which had significantly more counts than average for that chromosome.
DNA methylation and transcription in A. phagocytophilum-infected neutrophils
Because A. phagocytophilum infection-induced genome-wide DNA hypermethylation, we investigated whether this affected gene transcription. Borjesson et al. previously investigated the effect of A. phagocytophilum infection on transcription of the neutrophil using microarrays . To access and compare these data to the DNA methylation analyses, the transcriptional profiling data were obtained from Gene Expression Omnibus (GEO accession number GSE2405) and re-normalized using RMA, as described above. Differential expression was determined for 24 h postinfection, and genes were sorted by chromosome and physical location based on annotations in human genome assembly hg18 used for both these data and the DNA methylation studies. To account for underestimations of significant discoveries using p values, we calculated the FDR adjusted p values, or q values, for each Spearman correlation between windows of DNA methylation and windows of differential gene expression on chromosome 6. We selected a FDR of a maximum of 10 % as a significant correlation.
Spearman correlation coefficient (rho)
B-cell lymphoma 2 gene
cyclin-dependent kinase inhibitor
cis-regulatory element annotation software
cytochrome b-245 beta chain gene
differential gene expression
extranodal marginal zone B-cell lymphoma
Gene Expression Omnibus
heme-binding subunit of NADPH oxidase
histone 3 lysine 26
a human promyelocytic cell line
Luminometric Methylation Assay
model-based analysis of ChIP-Seq software
DNA methyl binding protein
methylated DNA binding domain protein immunoprecipitation of bound DNA followed by sequence analysis
major histocompatibility complex
nicotinamide adenine dinucleotide phosphate
genomic region downstream of MHC locus on chromosome 6
ras-related C3 botulinum toxin substrate 2 gene
robust multi-array average normalization
Statistical model for Identification of ChIP-Enriched Regions software
transcriptional start site
transcription termination site
The authors wish to acknowledge the contributions and discussions with members of SHGS’s thesis committee, including Stephen Baylin, M.D. (Johns Hopkins University School of Medicine), Srinivasan Yegnasubramanian M.D., Ph.D. (Johns Hopkins University School of Medicine), and Frank DeLeo, Ph.D. (Rocky Mountain Laboratory, National Institute of Allergy and Infectious Disease). The authors would also like to thank Kristen Rennoll-Bankert, Ph.D. and Jose Carlos Garcia-Garcia, Ph.D. for their discussion and day-to-day support and Valeria Pappas-Brown, Ph.D. and Emily Clemens for help with technical work. Thank you to S. Yegnasubramanian, M.D., Ph.D. and Sarah Wheelan, M.D., Ph.D. and the Johns Hopkins University Sidney Kimmel Comprehensive Cancer Center Next Generation Sequencing Center (http://nextgenseq.onc.jhmi.edu) for performing the MBD-seq and Conover Talbot, B.S. from the Johns Hopkins Deep Sequencing and Microarray Core for his assistance re-analyzing microarray using current tools and for the thoughtful discussion. Additional thanks to Cosme Adrover-Pacheco, Ph.D. for help with biocomputing analyses. Funding to support this work was provided by NIAID grant R01AI044102 to JSD. SY was supported by NIH/NCI grant P30CA006973, the Cleveland Foundation Masenhimer Fellowship.
- Chen SM, Dumler JS, Bakken JS, Walker DH. Identification of a granulocytotropic Ehrlichia species as the etiologic agent of human disease. J Clin Microbiol. 1994;32:589–95.PubMed CentralPubMedGoogle Scholar
- Bakken JS, Dumler JS, Chen SM, Eckman MR, Van Etta LL, Walker DH. Human granulocytic ehrlichiosis in the upper Midwest United States. A new species emerging? JAMA. 1994;272:212–8.PubMedView ArticleGoogle Scholar
- Carlyon JA, Chan WT, Galan J, Roos D, Fikrig E. Repression of rac2 mRNA expression by Anaplasma phagocytophila is essential to the inhibition of superoxide production and bacterial proliferation. J Immunol. 2002;169:7009–18.PubMedView ArticleGoogle Scholar
- Carlyon JA, Fikrig E. Mechanisms of evasion of neutrophil killing by Anaplasma phagocytophilum. Curr Opin Hematol. 2006;13:28–33.PubMedView ArticleGoogle Scholar
- Carlyon JA, Fikrig E. Invasion and survival strategies of Anaplasma phagocytophilum. Cell Microbiol. 2003;5:743–54.PubMedView ArticleGoogle Scholar
- Dumler JS, Choi KS, Garcia-Garcia JC, Barat NS, Scorpio DG, Garyu JW, et al. Human granulocytic anaplasmosis and Anaplasma phagocytophilum. Emerg Infect Dis. 2005;11:1828–34.PubMed CentralPubMedView ArticleGoogle Scholar
- Choi KS, Park JT, Dumler JS. Anaplasma phagocytophilum delay of neutrophil apoptosis through the p38 mitogen-activated protein kinase signal pathway. Infect Immun. 2005;73:8209–18.PubMed CentralPubMedView ArticleGoogle Scholar
- Garyu JW, Choi KS, Grab DJ, Dumler JS. Defective phagocytosis in Anaplasma phagocytophilum-infected neutrophils. Infect Immun. 2005;73:1187–90.PubMed CentralPubMedView ArticleGoogle Scholar
- Borjesson DL, Kobayashi SD, Whitney AR, Voyich JM, Argue CM, Deleo FR. Insights into pathogen immune evasion mechanisms: Anaplasma phagocytophilum fails to induce an apoptosis differentiation program in human neutrophils. J Immunol. 2005;174:6364–72.PubMedView ArticleGoogle Scholar
- de la Fuente J, Ayoubi P, Blouin EF, Almazan C, Naranjo V, Kocan KM. Gene expression profiling of human promyelocytic cells in response to infection with Anaplasma phagocytophilum. Cell Microbiol. 2005;7:549–59.PubMedView ArticleGoogle Scholar
- Lee HC, Kioi M, Han J, Puri RK, Goodman JL. Anaplasma phagocytophilum-induced gene expression in both human neutrophils and HL-60 cells. Genomics. 2008;92:144–51. http://dx.doi.org/10.1016/j.ygeno.2008.05.005.PubMedView ArticleGoogle Scholar
- Banerjee R, Anguita J, Roos D, Fikrig E. Cutting edge: infection by the agent of human granulocytic ehrlichiosis prevents the respiratory burst by down-regulating gp91phox. J Immunol. 2000;164:3946–9.PubMedView ArticleGoogle Scholar
- Ge Y, Yoshiie K, Kuribayashi F, Lin M, Rikihisa Y. Anaplasma phagocytophilum inhibits human neutrophil apoptosis via upregulation of bfl-1, maintenance of mitochondrial membrane potential and prevention of caspase 3 activation. Cell Microbiol. 2005;7:29–38.PubMedView ArticleGoogle Scholar
- Akkoyunlu M, Malawista SE, Anguita J, Fikrig E. Exploitation of interleukin-8-induced neutrophil chemotaxis by the agent of human granulocytic ehrlichiosis. Infect Immun. 2001;69:5577–88.PubMed CentralPubMedView ArticleGoogle Scholar
- Klein MB, Hu S, Chao CC, Goodman JL. The agent of human granulocytic ehrlichiosis induces the production of myelosuppressing chemokines without induction of proinflammatory cytokines. J Infect Dis. 2000;182:200–5.PubMedView ArticleGoogle Scholar
- Scorpio DG, Akkoyunlu M, Fikrig E, Dumler JS. CXCR2 blockade influences Anaplasma phagocytophilum propagation but not histopathology in the mouse model of human granulocytic anaplasmosis. Clin Diagn Lab Immunol. 2004;11:963–8.PubMed CentralPubMedGoogle Scholar
- Garcia-Garcia JC, Barat NC, Trembley SJ, Dumler JS. Epigenetic silencing of host cell defense genes enhances intracellular survival of the rickettsial pathogen Anaplasma phagocytophilum. PLoS Pathog. 2009;5:e1000488.PubMed CentralPubMedView ArticleGoogle Scholar
- Garcia-Garcia JC, Rennoll-Bankert KE, Pelly S, Milstone AM, Dumler JS. Silencing of host cell CYBB gene expression by the nuclear effector AnkA of the intracellular pathogen Anaplasma phagocytophilum. Infect Immun. 2009;77:2385–91.PubMed CentralPubMedView ArticleGoogle Scholar
- Stirzaker C, Song JZ, Davidson B, Clark SJ. Transcriptional gene silencing promotes DNA hypermethylation through a sequential change in chromatin modifications in cancer cells. Cancer Res. 2004;64:3871–7. doi:10.1158/0008-5472.CAN-03-3690.PubMedView ArticleGoogle Scholar
- Srivastava S, Mishra RK, Dhawan J. Regulation of cellular chromatin state: insights from quiescence and differentiation. Organogenesis. 2010;6:37–47.PubMed CentralPubMedView ArticleGoogle Scholar
- Baylin SB, Herman JG. DNA hypermethylation in tumorigenesis: epigenetics joins genetics. Trends Genet. 2000;16:168–74.PubMedView ArticleGoogle Scholar
- Baylin SB, Esteller M, Rountree MR, Bachman KE, Schuebel K, Herman JG. Aberrant patterns of DNA methylation, chromatin formation and gene expression in cancer. Hum Mol Genet. 2001;10:687–92.PubMedView ArticleGoogle Scholar
- Bierne H, Hamon M, Cossart P. Epigenetics and bacterial infections. Cold Spring Harb Perspect Med. 2012;2:a010272. doi:10.1101/cshperspect.a010272.PubMed CentralPubMedView ArticleGoogle Scholar
- Tolg C, Bagli DJ. Uropathogenic Escherichia coli infection: potential importance of epigenetics. Epigenomics. 2012;4:229–35. doi:10.2217/epi.12.5.PubMedView ArticleGoogle Scholar
- Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008;9:R137. doi:10.1186/gb-2008-9-9-r137.PubMed CentralPubMedView ArticleGoogle Scholar
- Zang C, Schones DE, Zeng C, Cui K, Zhao K, Peng W. A clustering approach for identification of enriched domains from histone modification ChIP-Seq data. Bioinformatics. 2009;25:1952–8. doi:10.1093/bioinformatics/btp340.PubMed CentralPubMedView ArticleGoogle Scholar
- Bell JT, Pai AA, Pickrell JK, Gaffney DJ, Pique-Regi R, Degner JF, et al. DNA methylation patterns associate with genetic and gene expression variation in HapMap cell lines. Genome Biol. 2011;12:R10. doi:10.1186/gb-2011-12-1-r10.PubMed CentralPubMedView ArticleGoogle Scholar
- Sinclair SH, Rennoll-Bankert KE, Dumler JS. Effector bottleneck: Microbial reprogramming of parasitized host cell transcription by epigenetic remodeling of chromatin structure. Front Genet. 2014;5:274. doi:10.3389/fgene.2014.00274.PubMed CentralPubMedView ArticleGoogle Scholar
- Fullwood MJ, Liu MH, Pan YF, Liu J, Xu H, Mohamed YB, et al. An oestrogen-receptor-alpha-bound human chromatin interactome. Nature. 2009;462:58–64. doi:10.1038/nature08497.PubMed CentralPubMedView ArticleGoogle Scholar
- Berman BP, Weisenberger DJ, Aman JF, Hinoue T, Ramjan Z, Liu Y, et al. Regions of focal DNA hypermethylation and long-range hypomethylation in colorectal cancer coincide with nuclear lamina-associated domains. Nat Genet. 2012;44:40–6. doi:10.1038/ng.969.PubMed CentralView ArticleGoogle Scholar
- Ciechomska M, van Laar JM, O’Reilly S. Emerging role of epigenetics in systemic sclerosis pathogenesis. Genes Immun. 2014;15:433–9. doi:10.1038/gene.2014.44.PubMedView ArticleGoogle Scholar
- de Mello VD, Pulkkinen L, Lalli M, Kolehmainen M, Pihlajamaki J, Uusitupa M. DNA methylation in obesity and type 2 diabetes. Ann Med. 2014;46:103–13. doi:10.3109/07853890.2013.857259.PubMedView ArticleGoogle Scholar
- Ronnerblad M, Andersson R, Olofsson T, Douagi I, Karimi M, Lehmann S, et al. Analysis of the DNA methylome and transcriptome in granulopoiesis reveals timed changes and dynamic enhancer methylation. Blood. 2014;123:e79–89. doi:10.1182/blood-2013-02-482893.PubMedView ArticleGoogle Scholar
- Zilbauer M, Rayner TF, Clark C, Coffey AJ, Joyce CJ, Palta P, et al. Genome-wide methylation analyses of primary human leukocyte subsets identifies functionally important cell-type-specific hypomethylated regions. Blood. 2013;122:e52–60. doi:10.1182/blood-2013-05-503201.PubMed CentralPubMedView ArticleGoogle Scholar
- Kondo T, Oka T, Sato H, Shinnou Y, Washio K, Takano M, et al. Accumulation of aberrant CpG hypermethylation by Helicobacter pylori infection promotes development and progression of gastric MALT lymphoma. Int J Oncol. 2009;35:547–57.PubMedGoogle Scholar
- Matsusaka K, Funata S, Fukayama M, Kaneda A. DNA methylation in gastric cancer, related to Helicobacter pylori and Epstein-Barr virus. World J Gastroenterol. 2014;20:3916–26. doi:10.3748/wjg.v20.i14.3916.PubMed CentralPubMedView ArticleGoogle Scholar
- Lee MJ, Min BJ, Choung HK, Kim N, Kim YA, Khwarg SI. Genome-wide DNA methylation profiles according to Chlamydophila psittaci infection and the response to doxycycline treatment in ocular adnexal lymphoma. Mol Vis. 2014;20:1037–47.PubMed CentralPubMedGoogle Scholar
- Luco RF, Allo M, Schor IE, Kornblihtt AR, Misteli T. Epigenetics in alternative pre-mRNA splicing. Cell. 2011;144:16–26. doi:10.1016/j.cell.2010.11.056.PubMed CentralPubMedView ArticleGoogle Scholar
- Maunakea AK, Nagarajan RP, Bilenky M, Ballinger TJ, D’Souza C, Fouse SD, et al. Conserved role of intragenic DNA methylation in regulating alternative promoters. Nature. 2010;466:253–7. doi:10.1038/nature09165.PubMed CentralPubMedView ArticleGoogle Scholar
- Jjingo D, Conley AB, Yi SV, Lunyak VV, Jordan IK. On the presence and role of human gene-body DNA methylation. Oncotarget. 2012;3:462–74.PubMed CentralPubMedGoogle Scholar
- Shenker N, Flanagan JM. Intragenic DNA methylation: implications of this epigenetic mechanism for cancer research. Br J Cancer. 2012;106:248–53. doi:10.1038/bjc.2011.550.PubMed CentralPubMedView ArticleGoogle Scholar
- Kulis M, Queiros AC, Beekman R, Martin-Subero JI. Intragenic DNA methylation in transcriptional regulation, normal differentiation and cancer. Biochim Biophys Acta. 1829;2013:1161–74. doi:10.1016/j.bbagrm.2013.08.001.Google Scholar
- Rennoll-Bankert KE, Sinclair SH, Lichay MA, Dumler JS. Comparison and characterization of granulocyte cell models for Anaplasma phagocytophilum infection. Pathog Dis. 2014;71:55–64. doi:10.1111/2049-632X.12111.PubMed CentralPubMedView ArticleGoogle Scholar
- Goodman JL, Nelson C, Vitale B, Madigan JE, Dumler JS, Kurtti TJ, et al. Direct cultivation of the causative agent of human granulocytic ehrlichiosis. N Engl J Med. 1996;334:209–15.PubMedView ArticleGoogle Scholar
- Dunning Hotopp JC, Lin M, Madupu R, Crabtree J, Angiuoli SV, Eisen J, et al. Comparative genomics of emerging human ehrlichiosis agents. PLoS Genet. 2006;2:e21.PubMed CentralPubMedView ArticleGoogle Scholar
- Guerrero-Preston R, Michailidi C, Marchionni L, Pickering CR, Frederick MJ, Myers JN, et al. Key tumor suppressor genes inactivated by “greater promoter” methylation and somatic mutations in head and neck cancer. Epigenetics. 2014;9:1031–46. doi:10.4161/epi.29025.PubMed CentralPubMedView ArticleGoogle Scholar
- Aryee MJ, Liu W, Engelmann JC, Nuhn P, Gurel M, Haffner MC, et al. DNA methylation alterations exhibit intraindividual stability and interindividual heterogeneity in prostate cancer metastases. Sci Transl Med. 2013;5:169ra10. doi:10.1126/scitranslmed.3005211.PubMed CentralPubMedView ArticleGoogle Scholar
- Yegnasubramanian S, Lin X, Haffner MC, DeMarzo AM, Nelson WG. Combination of methylated-DNA precipitation and methylation-sensitive restriction enzymes (COMPARE-MS) for the rapid, sensitive and quantitative detection of DNA methylation. Nucleic Acids Res. 2006;34:e19. doi:10.1093/nar/gnj022.PubMed CentralPubMedView ArticleGoogle Scholar
- Yegnasubramanian S, Wu Z, Haffner MC, Esopi D, Aryee MJ, Badrinath R, et al. Chromosome-wide mapping of DNA methylation patterns in normal and malignant prostate cells reveals pervasive methylation of gene-associated and conserved intergenic sequences. BMC Genomics. 2011;12:313. doi:10.1186/1471-2164-12-313.PubMed CentralPubMedView ArticleGoogle Scholar
- Easwaran H, Johnstone SE, Van Neste L, Ohm J, Mosbruger T, Wang Q, et al. A DNA hypermethylation module for the stem/progenitor cell signature of cancer. Genome Res. 2012;22:837–49. doi:10.1101/gr.131169.111.PubMed CentralPubMedView ArticleGoogle Scholar
- Karimi M, Johansson S, Stach D, Corcoran M, Grander D, Schalling M, et al. LUMA (LUminometric Methylation Assay)—a high throughput method to the analysis of genomic DNA methylation. Exp Cell Res. 2006;312:1989–95. doi:10.1016/j.yexcr.2006.03.006.PubMedView ArticleGoogle Scholar
- Ji X, Li W, Song J, Wei L, Liu XS. CEAS: cis-regulatory element annotation system. Nucleic Acids Res. 2006;34:W551–4. doi:10.1093/nar/gkl322.PubMed CentralPubMedView ArticleGoogle Scholar
- Shin H, Liu T, Manrai AK, Liu XS. CEAS: cis-regulatory element annotation system. Bioinformatics. 2009;25:2605–6. doi:10.1093/bioinformatics/btp479.PubMedView ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.