Skip to main content

 

Main menu

  • Home
  • Articles
    • Current Issue
    • Archive
    • Editorials
    • Immunology Education
    • Most Read
    • On The Horizon Collection
  • Info
    • About the Journal
    • For Authors
    • Journal Policies
    • For Advertisers
  • Editors
  • Submit
    • Submit a Manuscript
    • Instructions for Authors
    • Journal Policies
  • Services
    • Email Alerts
    • RSS Feeds
  • More
    • AAI Disclaimer
    • Feedback
  • Other Publications
    • American Association of Immunologists
    • The Journal of Immunology

User menu

  • Log in

Search

  • Advanced search
ImmunoHorizons
  • Other Publications
    • American Association of Immunologists
    • The Journal of Immunology
  • Log in
ImmunoHorizons

Advanced Search

  • Home
  • Articles
    • Current Issue
    • Archive
    • Editorials
    • Immunology Education
    • Most Read
    • On The Horizon Collection
  • Info
    • About the Journal
    • For Authors
    • Journal Policies
    • For Advertisers
  • Editors
  • Submit
    • Submit a Manuscript
    • Instructions for Authors
    • Journal Policies
  • Services
    • Email Alerts
    • RSS Feeds
  • More
    • AAI Disclaimer
    • Feedback
  • Follow ImmunoHorizons on Twitter
  • Follow ImmunoHorizons on RSS
Open Access

DNA Methylation Changes in Regional Lung Macrophages Are Associated with Metabolic Differences

David A. Armstrong, Youdinghuan Chen, John A. Dessaint, Daniel S. Aridgides, Jacqueline Y. Channon, Diane L. Mellinger, Brock C. Christensen and Alix Ashare
ImmunoHorizons July 1, 2019, 3 (7) 274-281; DOI: https://doi.org/10.4049/immunohorizons.1900042
David A. Armstrong
*Department of Medicine, Dartmouth-Hitchcock Medical Center, Lebanon, NH 03756;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for David A. Armstrong
Youdinghuan Chen
†Department of Epidemiology, Geisel School of Medicine at Dartmouth, Lebanon, NH 03756;
‡Department of Molecular and Systems Biology, Geisel School of Medicine at Dartmouth, Lebanon, NH 03756;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
John A. Dessaint
*Department of Medicine, Dartmouth-Hitchcock Medical Center, Lebanon, NH 03756;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Daniel S. Aridgides
*Department of Medicine, Dartmouth-Hitchcock Medical Center, Lebanon, NH 03756;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jacqueline Y. Channon
§Department of Microbiology and Immunology, Geisel School of Medicine at Dartmouth, Lebanon, NH 03756; and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Jacqueline Y. Channon
Diane L. Mellinger
*Department of Medicine, Dartmouth-Hitchcock Medical Center, Lebanon, NH 03756;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Brock C. Christensen
†Department of Epidemiology, Geisel School of Medicine at Dartmouth, Lebanon, NH 03756;
‡Department of Molecular and Systems Biology, Geisel School of Medicine at Dartmouth, Lebanon, NH 03756;
¶Department of Community and Family Medicine, Geisel School of Medicine at Dartmouth, Lebanon, NH 03756
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Brock C. Christensen
Alix Ashare
*Department of Medicine, Dartmouth-Hitchcock Medical Center, Lebanon, NH 03756;
§Department of Microbiology and Immunology, Geisel School of Medicine at Dartmouth, Lebanon, NH 03756; and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF + SI
  • PDF
Loading

Abstract

A number of pulmonary diseases occur with upper lobe predominance, including cystic fibrosis and smoking-related chronic obstructive pulmonary disease. In the healthy lung, several physiologic and metabolic factors exhibit disparity when comparing the upper lobe of the lung to lower lobe, including differences in oxygenation, ventilation, lymphatic flow, pH, and blood flow. In this study, we asked whether these regional differences in the lung are associated with DNA methylation changes in lung macrophages that could potentially lead to altered cell responsiveness upon subsequent environmental challenge. All analyses were performed using primary lung macrophages collected via bronchoalveolar lavage from healthy human subjects with normal pulmonary function. Epigenome-wide DNA methylation was examined via Infinium MethylationEPIC (850K) array and validated by targeted next-generation bisulfite sequencing. We observed 95 CpG loci with significant differential methylation in lung macrophages, comparing upper lobe to lower lobe (all false discovery rate < 0.05). Several of these genes, including CLIP4, HSH2D, NR4A1, SNX10, and TYK2, have been implicated as participants in inflammatory/immune-related biological processes. Functionally, we identified phenotypic differences in oxygen use, comparing upper versus lower lung macrophages. Our results support a hypothesis that epigenetic changes, specifically DNA methylation, at a multitude of gene loci in lung macrophages are associated with metabolic differences regionally in lung.

Introduction

Numerous lung diseases display upper lobe predominance, including cystic fibrosis, sarcoidosis, tuberculosis, and smoking-related chronic obstructive pulmonary disease (1). There are well-described physiologic differences in oxygenation and ventilation between the upper and lower lobes of the lung in the upright position (2). The origin of these differences has been largely attributed to gravity’s effects on the ventilation/perfusion ratio (VA/Q) in the lung, which is defined as the amount of inspired air divided by the pulmonary blood flow in an alveolus (3). Metabolic factors, such as lymphatic flow, pH, and blood flow, differ because of these inequalities of the VA/Q in the lung. Additionally, an oxygen gradient occurs in the lung because of this VA/Q from the apex to the base of the lung, with the upper lobe of the lung being notably more oxygen rich than the lower lobe at end expiration in an upright position (1). Hypoxia is known to induce broad effects on various aspects of cell biology, including alteration of the transcription of hundreds of genes that allow cells to adapt to the changing oxygen availability in the local environment (4). However, how metabolic factors contribute to upper lobe predominant lung disease is unclear.

Macrophages are the sentinel innate immune cells of the lung and possess remarkable immune plasticity, with the ability to sense and adapt to the local milieu (5). Lung macrophages play key roles in bacterial recognition and elimination as well as in polarization of innate and adaptive immunity. Depending on the local environment, macrophages sometimes play a role in anti-inflammatory responses, tissue repair, and homeostasis, whereas at other times, they promote inflammatory and phagocytic processes through the complex production of cytokines and cellular interaction (6, 7).

Epigenetics is the study of heritable changes in gene function caused by mechanisms other than changes in the underlying DNA sequence (8). Epigenetic mechanisms have emerged as modulators of host defenses that can lead to a more-prominent immune response and shape the course of inflammation in the host, both driving the production of specific inflammatory mediators and controlling the magnitude of the host response (9). It is now clear that both the innate and adaptive immune responses use epigenetic mechanisms of gene regulation to maintain long-term phenotypes (10). One important mode of epigenetic regulation is DNA methylation (11, 12).

In this study, we investigate whether DNA methylation changes in lung macrophages are associated with the metabolic factors that contribute to the regional phenotypic differences in lung. To examine this question, we initiated an epigenome-wide DNA methylation profiling study of lung macrophages, comparing upper lobe cells to lower lobe–derived cells.

Materials and Methods

This study was approved by the Committee for the Protection of Human Subjects at the Geisel School of Medicine at Dartmouth (approval no. 22781). All subjects were healthy adults and nonsmokers without any known underlying lung disease.

Bronchoalveolar lavage and macrophage isolation

Following local anesthesia with nebulized lidocaine and i.v. conscious sedation, subjects underwent flexible bronchoscopy. Briefly, a flexible fiberoptic bronchoscope was inserted transorally and advanced through the vocal cords. Bronchoalveolar lavage (BAL) fluid was obtained from tertiary airways in the right upper lobes (RUL) and right lower lobes (RLL). BAL was performed sequentially in the RUL and RLL with 20 ml of sterile saline followed by 10 ml of air, and this was repeated for a total of five times per airway. Lung macrophages were isolated, as previously described (13, 14). BAL fluid was filtered through two-layer gauze, centrifuged, and washed twice in 0.9% NaCl. Cells were counted with a T10 automated cell counter (Bio-Rad, Hercules, CA). Cytospins were performed using a Shandon Cytospin 3 centrifuge (Thermo Fisher Scientific, Waltham, MA). Briefly, 75,000 cells resuspended in 200 μl of 0.9% NaCl were loaded into a cytology funnel (Fisher Scientific, Pittsburgh, PA) and centrifuged for 10 min. Cells were allowed to air dry and processed for viewing via Hema 3 Stat pack (Fisher). Imaging was done on an Olympus (Waltham, MA) BX41 microscope with DP2-BSW software (version 2007).

DNA methylation array

Epigenome-wide DNA methylation profiling was performed via the Infinium Methylation EPIC Bead Chips (Illumina, San Diego, CA) for the determination of methylation levels of more than 850,000 CpG sites, as previously described (15, 16). DNA was extracted from lung macrophages via Qiagen (Germantown, MD) DNeasy Blood and Tissue Kit. DNA was quantitated on a Qubit 3.0 Fluorometer (Life Technologies, Carlsbad, CA). Bisulfite conversion of DNA was carried out with the Zymo EZ DNA methylation kit (Zymo Research, Irvine, CA), and EPIC array hybridization and scanning were performed at the University of Southern California Molecular Genomics Core.

DNA methylation array data processing

Raw intensity data files from the MethylationEPIC BeadChips were processed by the minfi R/Bioconductor analysis pipeline (version 1.21) (17) with annotation file version ilm10b3.hg19. All samples showed consistent and high-quality profiles. CpGs annotated as single nucleotide polymorphisms, technical probes, and sex chromosome associated, as well as those failing to meet a detection p value < 0.05 in ≥20% samples, were excluded. This preprocessing procedure left 813,096 CpGs with high-quality methylation data in the final data set. DNA methylation microarray data has been deposited into Gene Expression Omnibus under accession no. GSE132547.

Targeted, next-generation bisulfite sequencing of candidate genes

Targeted next-generation bisulfite sequencing (tNGBS) was performed on nine of the original 13 BAL specimens measured in the DNA EPICmethylation array by EpigenDx (Hopkinton, MA). DNA bisulfite modification was done using EZ-96 DNA methylation kit (Zymo), followed by multiplex PCR with Qiagen HotStar Taq and products purified with QIAquick PCR purification kit. Libraries were prepared using the KAPA Library Preparation Kit for Ion Torrent platforms and Ion Xpress Barcode Adapters (Thermo Fisher Scientific). Library products were purified using Agencourt AMPure XP beads (Beckman Coulter, Indianapolis, IN) and quantified using the Qiagen QIAxcel Advanced System. Barcoded samples were then pooled in an equimolar fashion before template preparation and enrichment were performed on the Ion Chef system (Thermo Fisher Scientific) using Ion 520 and Ion 530 Chef reagents. Enriched, template-positive library products were sequenced on the Ion S5 sequencer using Ion 530 sequencing chips (Thermo Fisher Scientific). FASTQ files from the Ion Torrent S5 server were aligned to the local reference database using open-source Bismark Bisulfite Read Mapper with the Bowtie2 alignment algorithm. Methylation levels were calculated in Bismark by dividing the number of methylated reads by the total number of reads.

Extracellular flux assay

Oxygen consumption rate (OCR) in lung macrophages was measured by Seahorse assay (Agilent, Santa Clara, CA), according to manufacturer’s instructions. Briefly, cells were plated at 50,000 per well and allowed to equilibrate overnight. The following morning, cells were washed twice and placed in XF base medium. The Seahorse instrument sensor cartridge was lowered to 200 μm above the cells, creating the transient microchamber. The sensor contains two fluorophores, one measuring oxygen (O2), the other measuring pH. Changes in oxygen concentration and pH were analyzed via XF Mito Stress Test Report Generator software (version 2.0) and reported as OCR in real time, enabling time-resolved, kinetic data of cell metabolism. Seahorse assays were carried out in DartLab, the Immune Monitoring and Flow Cytometry Shared Resource at the Norris Cotton Cancer Center at Dartmouth-Hitchcock Medical Center.

Statistical analysis

To determine the data structure of methylation profiles with respect to sample characteristics, we performed an unsupervised clustering analysis using the recursively partitioned mixture model (RPMM) method (18), assuming two terminal clusters on 10,341 most variable CpGs (all intersample variance > 0.01). RPMM was implemented in the RPMM R package (version 1.25).

To investigate the putative cell-type proportions in all samples, the RefFreeEWAS algorithm (R package version 2.1) (19) was applied to 10,000 most variable CpGs across all samples. The proportion of putative cell types was calculated iteratively for the number of cell types K from 2 to 10. The optimal number of putative cell types K = 2 was selected, as it minimized the variance of the bootstrapped deviance matrix. Subsequently, the two putative cell types by sample were displayed as stratified boxplots for visualization purposes. Putative cell-type 1 proportion was included as a fixed-effect term in the differential methylation analysis linear model, as described below.

The distribution of intersample variance in methylation β value was examined prior to differential methylation analysis: 27,107 CpGs with intersample variance in β values exceeding 0.005 were selected for further investigation. To account for the upper and lower lobe measurements collected from the same subject, the correlation coefficient (0.935) for the 27,107 CpG β values among the unique subjects (n = 13) was first calculated by the duplicateCorrelation function in the limma R/Bioconductor package (version 3.38.3). Differential methylation analysis was implemented by passing logit-transformed β values (i.e., M values), the consensus correlation coefficient for matched subject pairs, into the lmFit and eBayes functions in limma, with adjustment of subject age, sex, and RefFree cell-type 1 as fixed effects and the subject as a random effect in the model, such that Y = β0 + XRUL βRUL + βage Xage + βsex Xsex + βcell type 1 Xcell type 1 + Random Effect (Subject Pair), where Y is the methylation β value, β0 is the intercept, X is a given covariate, β is the respective model coefficient, XRUL is an indicator variable of sample from the right upper lobe, and βRUL is the linear model coefficient associated with it. Data processing, statistical analysis and data visualization R code can be found at https://github.com/Christensen-Lab-Dartmouth/CF_Epigenetics.

Results

Healthy subjects (n = 13) underwent flexible bronchoscopy to obtain BAL fluid from tertiary airways in both the RUL and RLL, and lung macrophages were isolated as described above. Seven female subjects and six male subjects had a mean age of 27.4 y (SD = 5.3). Genome-scale DNA methylation was measured in lung macrophages from both RUL and RLL, and the most variably methylated 10,341 CpGs from the EPIC methylation array were selected for unsupervised model-based clustering with RPMM (Fig. 1). In nearly all patients, matched RUL and RLL samples cluster immediately adjacent to each other, indicating that intersample variance in DNA methylation is larger than intrasample variance. For all subjects, both the RUL and RLL sample share RPMM cluster membership. The genomic context for the most variable probes was evenly distributed across promoters/nonpromoters, enhancers/nonenhancers, and CpG islands/non-CpG islands.

FIGURE 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
FIGURE 1.

Unsupervised clustering analysis, upper lobe versus lower lobe macrophages.

RPMM of the 10,000 CpGs (rows) with greatest sample variance across subjects. Individual subjects 1–13 are shown in columns with sample status bar at top: black (upper) and gray (lower). Blue color represents increased sample methylation.

To determine the putative cell-type composition, we applied the RefFreeEWAS algorithm (R package v.2.1, Linux R version 3.5) to 10,000 most variable CpGs, with 500 bootstraps, 10 iterations, and nine presumptive cell types (i.e., K = 2–10), as previously described (16). We chose K = 2 as this value minimized the column variance of the deviance matrix. The distribution of the two putative cell types was very similar between matched RUL and RLL samples (Supplemental Fig. 1). We also explored phenotypic differences between lung macrophages of the upper and lower lobes. Visual evaluation of lung macrophage morphology via cytospin revealed cells similar in presentation from both lobes: typically, a monocytoid appearance with reniform nuclei and ample cytosol (Fig. 2A, 2B), comprising what appeared to be >85% of the total cell population.

FIGURE 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
FIGURE 2.

BAL cell cytospins.

BAL cell samples were obtained from tertiary airways from both the RUL (A) and RLL (B) of healthy subjects (n = 13). BAL cells from healthy subjects are visually a mostly homogeneous population of cells, primarily lung macrophages with oval to reniform nuclei and abundant cytosol. Images shown are representative of multiple subjects. BAL cells were stained with H&E. Original magnification ×400.

Next, we investigated the relationship between DNA methylation in upper lobe versus lower lobe lung macrophages. Differential analysis comparing RUL and RLL revealed a total of 95 CpGs at a false discovery rate (FDR) < 0.05 (Fig. 3) and 372 CpGs at an FDR < 0.10 (Supplemental Table I). Of the 95 CpGs with an FDR < 0.05, the top 15 CpGs associated with genes of known function are shown in Table I. CpG location varied across multiple University of California, Santa Cruz reference gene groups, including 5′ untranslated region (UTR), TSS1500, TSS200, and body. The observed change in methylation values ranged from −8.2 to 10.7% for upper lobe cells. Upper lobe cells exhibited increased methylation at CpGs tracking to several genes previously implicated in host–pathogen interaction or in inflammatory/immune-related biological processes, including CAP-Gly domain-containing linker protein family member 4 (CLIP4), HSH2D, NR4A1, SNX10, and TYK2.

FIGURE 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
FIGURE 3.

Epigenome-wide differential methylation in RUL macrophages.

Comparative analysis of DNA methylation in upper versus lower lobe macrophages in healthy subjects identified 95 differentially methylated CpGs (FDR p < 0.05). CpGs hypermethylated in RUL (green), and hypomethylated CpGs (blue) are plotted as log2-fold increase or decrease in methylation M value (x-axis) versus log10 FDR-adjusted p value (y-axis). Statistically significant CpGs associated with specific genes are labeled, and unlabeled points represent CpGs associated with no known gene at that location.

View this table:
  • View inline
  • View popup
Table I. Top EPIC CpGs with differential methylation in upper lobe macrophages associated with genes of known function

Subsequently, we used tNGBS for validation of the results from Infinium MethylationEPIC array. Assays were designed and run for NR4A1 (cg17323256), A2M (cg00889217), RUNX2 (cg13427753), and CLIP4 (cg26118047). A number of the original Infinium EPIC array CpGs were situated within repeat elements (LINE or SINE), and tNGBS assay design was unsuccessful (i.e., HSH2D, TYK2, SNX10, and SCNN1A). Methylation values from tNGBS were consistent with the EPIC array. Percent methylation difference (Δβ) values for select gene-associated CpGs are presented in Table II. Additionally, an example of tNGBS assay design for CLIP4 is illustrated in Supplemental Fig. 2.

View this table:
  • View inline
  • View popup
Table II. Comparison of percent methylation changes (Δβ) between EPIC array and tNGBS platforms at gene-specific CpGs

At all targeted gene loci, the two methylation platforms, Illumina EPIC array versus tNGBS, showed significant correlations: NR4A1, CLIP4, RUNX2, A2M, PDE4D, and ZFHX3 (Fig. 4). Correlation and p values ranged from 0.769 and p = 1.93e−04 (ZFHX3) to 0.96 and p = 2.79e−10 (PDE4D).

FIGURE 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
FIGURE 4.

Correlation of DNA methylation analysis between Illumina EPIC methylation array and tNGBS platforms.

Correlation coefficients and linear regression p values were determined for CpGs associated with the genes NR4A1, CLIP4, RUNX2, A2M, PDE4D, and ZFHX3 as methodology validation, directly comparing the Illumina EPIC DNA methylation array and tNGBS assays.

Lung macrophages were also evaluated for differences in metabolic function by measuring basal respiration levels. OCR was assayed as a measure of basal mitochondrial function in cells comparing upper versus lower. The XF Cell Mito Stress Test performed on the Seahorse extracellular flux analyzer revealed a 21.4% difference (p < 0.001) in basal respiration between RUL versus RLL cells, with lower lobe cells exhibiting a higher OCR (Fig. 5).

FIGURE 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
FIGURE 5.

Oxygen consumption is higher in lower lobe macrophages.

Primary lung macrophages isolated from BAL fluid were plated and allowed to equilibrate in vitro for 16 h. OCR was measured via a two-fluorophore sensor cartridge via the Seahorse Mito Stress assay. Data are pooled from six to eight wells per experiment, with lung macrophages from three separate donors. Values were normalized to percent RUL OCR. ***p < 0.001 by unpaired Student t test.

Discussion

It has been recognized for more than a half century that differences in physiologic and metabolic factors, such as lymphatic flow, ventilation, pH, and oxygenation levels, exist regionally within the lung (2, 3, 20–22). However, the impact of these physiologic differences on disease pathology or the immune response in upper lobe predominant lung diseases is unclear. We hypothesized that regional physiologic differences may contribute to epigenetic variability in lung macrophages isolated from the upper and lower lobes of the lung, as these cells are long lived and would be prone to epigenetic modifications in response to environmental differences. The goal of this study was to determine if epigenetic changes in lung macrophages, specifically DNA methylation changes, exist regionally in the lung (i.e., upper versus lower lobe) that may ultimately predispose upper lobe lung macrophages to respond differently than lower lobe macrophages.

Limited studies to date have examined lung-related DNA methylation profiling or differences in regional innate immune cells of the lung. Many such DNA methylation studies are lung cancer related (23–27). Additionally, several groups have looked at the effects of cigarette smoking on lung-related DNA methylation (28, 29), and various studies have investigated changes in the DNA methylome in the single cell–type setting (30, 31). One of our recent studies has identified differential methylation in cystic fibrosis subjects at 109 gene-associated CpGs in lung macrophages, using the Illumina Infinium MethylationEPIC (850 K) array (16). In addition, we have demonstrated variant regional lung microRNA and cytokine expression associated with hypoxia in lung macrophages and BAL fluid (upper versus lower lobe) (32).

In this study, we report DNA methylation differences from regional lung macrophages, providing a direct comparison of the epigenetic molecular signature of upper lobe versus lower lobe innate immune cells from healthy individuals. Our cell collection methodology (flexible bronchoscopy) is a rare approach used in research studies and affords us a unique opportunity to analyze innate immune cells isolated directly from the airway and alveolar space.

Our initial strategy was to collect lung macrophages of both the upper lobe and lower lobe in healthy subjects to analyze epigenome-wide DNA methylation patterns. DNA methylation profiling via unsupervised clustering analysis showed that matched upper lobe and lower lobe samples cluster immediately adjacent to each other in nearly all subjects. Additionally, we identified 95 differentially methylated CpGs (all FDR < 0.05), of which 80 were hypermethylated in upper lobe cells and 15 were hypomethylated in upper lobe cells as well as 372 differential CpGs (all FDR < 0.10). Differentially methylated CpGs included those associated with genes such as NR4A1, SNX10, SCNN1A, and A2M. The orphan nuclear receptor NR4A1, which has been shown to be rapidly induced in response to stimuli, such as LPS (33), has been linked to the regulation of glucose metabolism and inflammation as well as implicated in hypoxia-induced apoptosis (34, 35). SNX10, a member of the family of sorting nexins, is a protein involved in intracellular protein sorting, trafficking, and signal transduction (36). Moreover, SNX10 has been reported to regulate endosomal morphology, which might be crucial for macrophage function, including phagocytosis and digestion of pathogens, inflammatory response, and Ag presentation (37). The SCNN1A gene encodes for the α subunit of the epithelial sodium channel (ENaC), which plays a critical role in the ion and fluid regulation of the lung. Altered activity of ENaC, which regulates the airway surface liquid layer, might lead to airway dehydration, mucus stasis, and bacterial overgrowth, as can be seen in cystic fibrosis and chronic bronchitis (38); however, an exact role for ENaC in lung macrophages, specifically, has yet to be fully established.

In addition to epigenetic changes in lung macrophages, we have also identified altered oxygen use in upper versus lower lobe cells. Lower lobe macrophages consistently exhibit a higher OCR than upper lobe macrophages. This increased OCR may represent a compensation mechanism to counter the decreased oxygen present in the lower lobe. Mitochondrial and cellular metabolic processes are central to immune cell responses, hypoxia sensing, and apoptosis in addition to their well-known roles of substrate oxidation and ATP production (39–41). Additionally, glucose metabolism is pivotal during macrophage pathogen phagocytosis (42). Although we do not currently have a direct causal link, one gene we have identified as differentially methylated, NR4A1, has been shown to be a transcriptional regulator of multiple glycolytic enzymes (35) and hence could ultimately be connected to changes in oxygen use levels in lung macrophages.

This study has several major strengths. To the best of our knowledge, this is the first study to demonstrate gene-associated DNA methylation differences in lung macrophages, comparing upper versus lower lobe. We have used the Infinium MethylationEPIC array, spanning 850,000 sites across the human genome, and validated these observations via tNGBS. In addition, we provide evidence of phenotypic differences of upper versus lower lobe cells, namely, a divergence in basal metabolism. However, a limitation of this study is that although we have identified epigenetic and a metabolic phenotypic difference between upper lobe and lower lobe cells in healthy subjects, a causative relationship between the DNA methylation changes and cellular phenotype is yet to be ascertained. Future work with larger sample size and power should investigate such a relationship.

In conclusion, through epigenome-wide DNA methylation profiling, we have identified 95 differential methylated CpGs in lung macrophages from the RUL relative to the RLL. Additionally, we have demonstrated energy metabolome-related phenotypical difference (i.e., oxygen consumption) in regional macrophages, upper versus lower. These observations support a hypothesis that epigenetic changes, specifically DNA methylation, at a multitude of gene loci in lung macrophages are associated with the well-established metabolic regional differences in lung and may, at least in part, be contributing to a differential innate immune response in alternate regions of the lung.

Disclosures

The authors have no financial conflicts of interest.

Footnotes

  • ↵1 Contributed equally to the manuscript.

  • This work was supported by National Institutes of Health (NIH) Grants R01HL122372 (to A.A.) and R01CA216265 (to B.C.C.). Subject enrollment and clinical sample collection were achieved with the assistance of the Cystic Fibrosis Translational Research Core at Dartmouth, which is jointly funded by NIH Grant P30GM106394 (to Bruce A. Stanton, principal investigator) and Cystic Fibrosis Foundation Award STANTO11R0 (to Bruce A. Stanton, principal investigator). This work was also supported by a Burroughs Wellcome Big Data in Life Sciences training grant and DartLab–National Cancer Institute Cancer Center Support Grant 5P30 CA023108-37.

  • The DNA methylation microarray data presented in this article have been submitted to the Gene Expression Omnibus (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE132547) under accession number GSE132547.

  • The online version of this article contains supplemental material.

  • Abbreviations used in this article:

    BAL
    bronchoalveolar lavage
    CLIP4
    CAP-Gly domain-containing linker protein family member 4
    ENaC
    epithelial sodium channel
    FDR
    false discovery rate
    OCR
    oxygen consumption rate
    RLL
    right lower lobe
    RPMM
    recursively partitioned mixture model
    RUL
    right upper lobe
    tNGBS
    targeted next-generation bisulfite sequencing
    UTR
    untranslated region
    VA/Q
    ventilation/perfusion ratio.

  • Received June 3, 2019.
  • Accepted June 5, 2019.
  • Copyright © 2019 The Authors

This article is distributed under the terms of the CC BY-NC 4.0 Unported license.

References

  1. ↵
    1. Nemec S. F.,
    2. A. A. Bankier,
    3. R. L. Eisenberg
    . 2013. Upper lobe-predominant diseases of the lung. AJR Am. J. Roentgenol. 200: W222–W237.
    OpenUrlCrossRefPubMed
  2. ↵
    1. West J. B.
    1962. Regional differences in gas exchange in the lung of erect man. J. Appl. Physiol. 17: 893–898.
    OpenUrlCrossRefPubMed
  3. ↵
    1. West J. B.
    1963. Blood-flow, ventilation, and gas exchange in the lung. Lancet 2: 1055–1058.
    OpenUrlPubMed
  4. ↵
    1. Cummins E. P.,
    2. C. T. Taylor
    . 2005. Hypoxia-responsive transcription factors. Pflugers Arch. 450: 363–371.
    OpenUrlCrossRefPubMed
  5. ↵
    1. Bruscia E. M.,
    2. T. L. Bonfield
    . 2016. Cystic fibrosis lung immunity: the role of the macrophage. J. Innate Immun. 8: 550–563.
    OpenUrl
  6. ↵
    1. Hirayama D.,
    2. T. Iida,
    3. H. Nakase
    . 2017. The phagocytic function of macrophage-enforcing innate immunity and tissue homeostasis. Int. J. Mol. Sci. 19: 92–105.
    OpenUrl
  7. ↵
    Rossi, G., A. Cavazza, P. Spagnolo, S. Bellafiore, E. Kuhn, P. Carassai, L. Caramanico, G. Montanari, G. Cappiello, A. Andreani, et al. 2017. The role of macrophages in interstitial lung diseases: Number 3 in the series “pathology for the clinician.” P. Dorfmüller and A. Cavazza, eds. Eur. Respir. Rev. 26: 1–16.
  8. ↵
    1. Wu C. T.,
    2. J. R. Morris
    . 2001. Genes, genetics, and epigenetics: a correspondence. Science 293: 1103–1105.
    OpenUrlAbstract/FREE Full Text
  9. ↵
    1. Morandini A. C.,
    2. C. F. Santos,
    3. Ö. Yilmaz
    . 2016. Role of epigenetics in modulation of immune response at the junction of host-pathogen interaction and danger molecule signaling. Pathog. Dis. 74: 1–8.
    OpenUrl
  10. ↵
    1. DiNardo A. R.,
    2. T. Nishiguchi,
    3. E. M. Mace,
    4. K. Rajapakshe,
    5. G. Mtetwa,
    6. A. Kay,
    7. G. Maphalala,
    8. W. E. Secor,
    9. R. Mejia,
    10. J. S. Orange, et al
    . 2018. Schistosomiasis induces persistent DNA methylation and tuberculosis-specific immune changes. J. Immunol. 201: 124–133.
    OpenUrlAbstract/FREE Full Text
  11. ↵
    1. Marsit C. J.,
    2. S. S. Brummel,
    3. D. Kacanek,
    4. G. R. Seage III,
    5. S. A. Spector,
    6. D. A. Armstrong,
    7. B. M. Lester,
    8. K. Rich, Pediatric HIV/AIDS Cohort Studies Network
    . 2015. Infant peripheral blood repetitive element hypomethylation associated with antiretroviral therapy in utero. Epigenetics 10: 708–716.
    OpenUrl
  12. ↵
    1. Selman M.,
    2. F. J. Martinez,
    3. A. Pardo
    . 2019. Why does an aging smoker’s lung develop idiopathic pulmonary fibrosis and not chronic obstructive pulmonary disease? Am. J. Respir. Crit. Care Med. 199: 279–285.
    OpenUrl
  13. ↵
    1. Bessich J. L.,
    2. A. B. Nymon,
    3. L. A. Moulton,
    4. D. Dorman,
    5. A. Ashare
    . 2013. Low levels of insulin-like growth factor-1 contribute to alveolar macrophage dysfunction in cystic fibrosis. J. Immunol. 191: 378–385.
    OpenUrlAbstract/FREE Full Text
  14. ↵
    1. Monick M. M.,
    2. L. S. Powers,
    3. C. W. Barrett,
    4. S. Hinde,
    5. A. Ashare,
    6. D. J. Groskreutz,
    7. T. Nyunoya,
    8. M. Coleman,
    9. D. R. Spitz,
    10. G. W. Hunninghake
    . 2008. Constitutive ERK MAPK activity regulates macrophage ATP production and mitochondrial integrity. J. Immunol. 180: 7485–7496.
    OpenUrlAbstract/FREE Full Text
  15. ↵
    1. Kling T.,
    2. A. Wenger,
    3. S. Beck,
    4. H. Carén
    . 2017. Validation of the MethylationEPIC BeadChip for fresh-frozen and formalin-fixed paraffin-embedded tumours. Clin. Epigenetics 9: 33–38.
    OpenUrl
  16. ↵
    1. Chen Y.,
    2. D. A. Armstrong,
    3. L. A. Salas,
    4. H. F. Hazlett,
    5. A. B. Nymon,
    6. J. A. Dessaint,
    7. D. S. Aridgides,
    8. D. L. Mellinger,
    9. X. Liu,
    10. B. C. Christensen,
    11. A. Ashare
    . 2018. Genome-wide DNA methylation profiling shows a distinct epigenetic signature associated with lung macrophages in cystic fibrosis. Clin. Epigenetics 10: 152–161.
    OpenUrl
  17. ↵
    1. Aryee M. J.,
    2. A. E. Jaffe,
    3. H. Corrada-Bravo,
    4. C. Ladd-Acosta,
    5. A. P. Feinberg,
    6. K. D. Hansen,
    7. R. A. Irizarry
    . 2014. Minfi: a flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics 30: 1363–1369.
    OpenUrlCrossRefPubMed
  18. ↵
    1. Houseman E. A.,
    2. B. C. Christensen,
    3. R. F. Yeh,
    4. C. J. Marsit,
    5. M. R. Karagas,
    6. M. Wrensch,
    7. H. H. Nelson,
    8. J. Wiemels,
    9. S. Zheng,
    10. J. K. Wiencke,
    11. K. T. Kelsey
    . 2008. Model-based clustering of DNA methylation array data: a recursive-partitioning algorithm for high-dimensional data arising as a mixture of beta distributions. BMC Bioinformatics 9: 365–379.
    OpenUrlCrossRefPubMed
  19. ↵
    1. Houseman E. A.,
    2. M. L. Kile,
    3. D. C. Christiani,
    4. T. A. Ince,
    5. K. T. Kelsey,
    6. C. J. Marsit
    . 2016. Reference-free deconvolution of DNA methylation data and mediation by cell composition effects. BMC Bioinformatics 17: 259.
    OpenUrlCrossRefPubMed
  20. ↵
    1. West J. B.
    1962. The effect of uneven blood flow in the lung on regional gas exchange. Med. Thorac. 19: 392–398.
    OpenUrlPubMed
    1. West J. B.
    1963. Distribution of gas and blood in the normal lungs. Br. Med. Bull. 19: 53–58.
    OpenUrlPubMed
  21. ↵
    1. West J. B.,
    2. R. A. Holland,
    3. C. T. Dollery,
    4. C. M. Matthews
    . 1962. Interpretation of radioactive gas clearance rates in the lung. J. Appl. Physiol. 17: 14–20.
    OpenUrlPubMed
  22. ↵
    1. Shi Y. X.,
    2. Y. Wang,
    3. X. Li,
    4. W. Zhang,
    5. H. H. Zhou,
    6. J. Y. Yin,
    7. Z. Q. Liu
    . 2017. Genome-wide DNA methylation profiling reveals novel epigenetic signatures in squamous cell lung cancer. BMC Genomics 18: 901–912.
    OpenUrl
    1. Wang X.,
    2. W. Shang,
    3. Y. Chang,
    4. X. Li
    . 2018. Methylation signature genes identification of the lung squamous cell carcinoma occurrence and recognition research. J. Comput. Biol. 25: 1161–1169.
    OpenUrl
    1. Um S. W.,
    2. Y. Kim,
    3. B. B. Lee,
    4. D. Kim,
    5. K. J. Lee,
    6. H. K. Kim,
    7. J. Han,
    8. H. Kim,
    9. Y. M. Shim,
    10. D. H. Kim
    . 2018. Genome-wide analysis of DNA methylation in bronchial washings. Clin. Epigenetics 10: 65–73.
    OpenUrl
    1. Stueve T. R.,
    2. W. Q. Li,
    3. J. Shi,
    4. C. N. Marconett,
    5. T. Zhang,
    6. C. Yang,
    7. D. Mullen,
    8. C. Yan,
    9. W. Wheeler,
    10. X. Hua, et al
    . 2017. Epigenome-wide analysis of DNA methylation in lung tissue shows concordance with blood studies and identifies tobacco smoke-inducible enhancers. Hum. Mol. Genet. 26: 3014–3027.
    OpenUrl
  23. ↵
    1. Kettunen E.,
    2. H. Hernandez-Vargas,
    3. M. P. Cros,
    4. G. Durand,
    5. F. Le Calvez-Kelm,
    6. K. Stuopelyte,
    7. S. Jarmalaite,
    8. K. Salmenkivi,
    9. S. Anttila,
    10. H. Wolff, et al
    . 2017. Asbestos-associated genome-wide DNA methylation changes in lung cancer. Int. J. Cancer 141: 2014–2029.
    OpenUrl
  24. ↵
    1. de Vries M.,
    2. D. A. van der Plaat,
    3. I. Nedeljkovic,
    4. R. N. Verkaik-Schakel,
    5. W. Kooistra,
    6. N. Amin,
    7. C. M. van Duijn,
    8. C. A. Brandsma,
    9. C. C. van Diemen,
    10. J. M. Vonk,
    11. H. Marike Boezen
    . 2018. From blood to lung tissue: effect of cigarette smoke on DNA methylation and lung function. Respir. Res. 19: 212–220.
    OpenUrl
  25. ↵
    1. Choukrallah M. A.,
    2. N. Sierro,
    3. F. Martin,
    4. K. Baumer,
    5. J. Thomas,
    6. S. Ouadi,
    7. J. Hoeng,
    8. M. C. Peitsch,
    9. N. V. Ivanov
    . 2019. Tobacco heating system 2.2 has a limited impact on DNA methylation of candidate enhancers in mouse lung compared with cigarette smoke. Food Chem. Toxicol. 123: 501–510.
    OpenUrl
  26. ↵
    1. Shi Y.,
    2. T. Zhao,
    3. X. Yang,
    4. B. Sun,
    5. Y. Li,
    6. J. Duan,
    7. Z. Sun
    . 2019. PM2.5-induced alteration of DNA methylation and RNA-transcription are associated with inflammatory response and lung injury. Sci. Total Environ. 650: 908–921.
    OpenUrl
  27. ↵
    1. Clifford R. L.,
    2. N. Fishbane,
    3. J. Patel,
    4. J. L. MacIsaac,
    5. L. M. McEwen,
    6. A. J. Fisher,
    7. C. A. Brandsma,
    8. P. Nair,
    9. M. S. Kobor,
    10. T. L. Hackett,
    11. A. J. Knox
    . 2018. Altered DNA methylation is associated with aberrant gene expression in parenchymal but not airway fibroblasts isolated from individuals with COPD. Clin. Epigenetics 10: 32–45.
    OpenUrl
  28. ↵
    1. Armstrong D. A.,
    2. A. B. Nymon,
    3. C. S. Ringelberg,
    4. C. Lesseur,
    5. H. F. Hazlett,
    6. L. Howard,
    7. C. J. Marsit,
    8. A. Ashare
    . 2017. Pulmonary microRNA profiling: implications in upper lobe predominant lung disease. Clin. Epigenetics 9: 56–64.
    OpenUrl
  29. ↵
    1. Ipseiz N.,
    2. S. Uderhardt,
    3. C. Scholtysek,
    4. M. Steffen,
    5. G. Schabbauer,
    6. A. Bozec,
    7. G. Schett,
    8. G. Krönke
    . 2014. The nuclear receptor Nr4a1 mediates anti-inflammatory effects of apoptotic cells. J. Immunol. 192: 4852–4858.
    OpenUrlAbstract/FREE Full Text
  30. ↵
    1. Wohlkoenig C.,
    2. K. Leithner,
    3. A. Olschewski,
    4. H. Olschewski,
    5. A. Hrzenjak
    . 2017. TR3 is involved in hypoxia-induced apoptosis resistance in lung cancer cells downstream of HIF-1α. Lung Cancer 111: 15–22.
    OpenUrl
  31. ↵
    1. Chao L. C.,
    2. K. Wroblewski,
    3. Z. Zhang,
    4. L. Pei,
    5. L. Vergnes,
    6. O. R. Ilkayeva,
    7. S. Y. Ding,
    8. K. Reue,
    9. M. J. Watt,
    10. C. B. Newgard, et al
    . 2009. Insulin resistance and altered systemic glucose metabolism in mice lacking Nur77. Diabetes 58: 2788–2796.
    OpenUrlAbstract/FREE Full Text
  32. ↵
    1. Worby C. A.,
    2. J. E. Dixon
    . 2002. Sorting out the cellular functions of sorting nexins. [Published erratum appears in 2003 Nat. Rev. Mol. Cell Biol. 4: 156.] Nat. Rev. Mol. Cell Biol. 3: 919–931.
    OpenUrlCrossRefPubMed
  33. ↵
    1. Lou J.,
    2. X. Li,
    3. W. Huang,
    4. J. Liang,
    5. M. Zheng,
    6. T. Xu,
    7. J. Lyu,
    8. D. Li,
    9. Q. Xu,
    10. X. Jin, et al
    . 2017. SNX10 promotes phagosome maturation in macrophages and protects mice against Listeria monocytogenes infection. Oncotarget 8: 53935–53947.
    OpenUrlCrossRefPubMed
  34. ↵
    1. Hamacher J.,
    2. Y. Hadizamani,
    3. M. Borgmann,
    4. M. Mohaupt,
    5. D. N. Männel,
    6. U. Moehrlen,
    7. R. Lucas,
    8. U. Stammberger
    . 2018. Cytokine-ion channel interactions in pulmonary inflammation. Front. Immunol. 8: 1644–1688.
    OpenUrl
  35. ↵
    1. Dimeloe S.,
    2. A. V. Burgener,
    3. J. Grählert,
    4. C. Hess
    . 2017. T-cell metabolism governing activation, proliferation and differentiation; a modular view. Immunology 150: 35–44.
    OpenUrl
    1. Ochocki J. D.,
    2. M. C. Simon
    . 2013. Nutrient-sensing pathways and metabolic regulation in stem cells. J. Cell Biol. 203: 23–33.
    OpenUrlAbstract/FREE Full Text
  36. ↵
    1. Smith R. A.,
    2. R. C. Hartley,
    3. H. M. Cochemé,
    4. M. P. Murphy
    . 2012. Mitochondrial pharmacology. Trends Pharmacol. Sci. 33: 341–352.
    OpenUrlCrossRefPubMed
  37. ↵
    1. Barghouthi S.,
    2. K. D. Everett,
    3. D. P. Speert
    . 1995. Nonopsonic phagocytosis of Pseudomonas aeruginosa requires facilitated transport of D-glucose by macrophages. J. Immunol. 154: 3420–3428.
    OpenUrlAbstract
PreviousNext
Back to top

In this issue

ImmunoHorizons: 3 (7)
ImmunoHorizons
Vol. 3, Issue 7
1 Jul 2019
  • Table of Contents
  • About the Cover
  • Editorial Board (PDF)
Print
Download PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for your interest in spreading the word about ImmunoHorizons.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
DNA Methylation Changes in Regional Lung Macrophages Are Associated with Metabolic Differences
(Your Name) has forwarded a page to you from ImmunoHorizons
(Your Name) thought you would like to see this page from the ImmunoHorizons web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
DNA Methylation Changes in Regional Lung Macrophages Are Associated with Metabolic Differences
David A. Armstrong, Youdinghuan Chen, John A. Dessaint, Daniel S. Aridgides, Jacqueline Y. Channon, Diane L. Mellinger, Brock C. Christensen, Alix Ashare
ImmunoHorizons July 1, 2019, 3 (7) 274-281; DOI: 10.4049/immunohorizons.1900042

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Share
DNA Methylation Changes in Regional Lung Macrophages Are Associated with Metabolic Differences
David A. Armstrong, Youdinghuan Chen, John A. Dessaint, Daniel S. Aridgides, Jacqueline Y. Channon, Diane L. Mellinger, Brock C. Christensen, Alix Ashare
ImmunoHorizons July 1, 2019, 3 (7) 274-281; DOI: 10.4049/immunohorizons.1900042
del.icio.us logo Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Disclosures
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • PDF + SI
  • PDF

Related Articles

Cited By...

More in this TOC Section

  • Caspase-8 Blocks Receptor-Interacting Protein Kinase-1 Kinase–Independent Necroptosis during Embryogenesis
  • Systems Immunology Analyses of STAT1 Gain-of-Function Immune Phenotypes Reveal Heterogeneous Response to IL-6 and Broad Immunometabolic Roles for STAT1
  • Transcriptional and Cytotoxic Responses of Human Intestinal Organoids to IFN Types I, II, and III
Show more INNATE IMMUNITY

Similar Articles

Navigate

  • Home
  • Current Issue
  • Archive

For Authors

  • Submit a Manuscript
  • Instructions for Authors
  • About the Journal
  • Journal Policies
  • Editors

General Information

  • FAR 889
  • Privacy Policy
  • Disclaimer

Journal Services

  • Email Alerts
  • RSS Feeds

Copyright© 2022 by The American Association of Immunologists, Inc

Online ISSN 2573-7732