Skip to main content

Transcriptome-wide analyses indicate mitochondrial responses to particulate air pollution exposure



Due to their lack of repair capacity mitochondria are critical targets for environmental toxicants. We studied genes and pathways reflecting mitochondrial responses to short- and medium-term PM10 exposure.


Whole genome gene expression was measured in peripheral blood of 98 adults (49% women). We performed linear regression analyses stratified by sex and adjusted for individual and temporal characteristics to investigate alterations in gene expression induced by short-term (week before blood sampling) and medium-term (month before blood sampling) PM10 exposure. Overrepresentation analyses (ConsensusPathDB) were performed to identify enriched mitochondrial associated pathways and gene ontology sets. Thirteen Human MitoCarta genes were measured by means of quantitative real-time polymerase chain reaction (qPCR) along with mitochondrial DNA (mtDNA) content in an independent validation cohort (n = 169, 55.6% women).


Overrepresentation analyses revealed significant pathways (p-value <0.05) related to mitochondrial genome maintenance and apoptosis for short-term exposure and to the electron transport chain (ETC) for medium-term exposure in women. For men, medium-term PM10 exposure was associated with the Tri Carbonic Acid cycle. In an independent study population, we validated several ETC genes, including UQCRH and COX7C (q-value <0.05), and some genes crucial for the maintenance of the mitochondrial genome, including LONP1 (q-value: 0.07) and POLG (q-value: 0.04) in women.


In this exploratory study, we identified mitochondrial genes and pathways associated with particulate air pollution indicating upregulation of energy producing pathways as a potential mechanism to compensate for PM-induced mitochondrial damage.

Peer Review reports


Mitochondria are cellular organelles specialized in energy production and produce the majority of intracellular reactive oxygen species (ROS), which are continually generated as toxic by-products by the electron transport chain (ETC). There is a fine balance in ROS signalling maintained by the redox environment. ROS production may be altered as a consequence of exposure to particulate matter (PM). Oxidative stress can occur both when the intracellular and/or intramitochondrial environments are either highly reduced or highly oxidized [1]. Mitochondrial DNA (mtDNA) repairs DNA damage less efficiently compared to nuclear DNA, making it susceptible to ROS and environmental toxicants such as PM [2]. Accumulation of mtDNA damage can cause disturbed mtDNA replication and elimination of damaged mtDNA, and in turn lead to decreased levels of mtDNA [3, 4]. Furthermore, components of the mitochondrial membrane rich in unsaturated fatty acids, such as cardiolipin, are especially sensitive to peroxidation by ROS, resulting in reactive aldehydes which can further damage mitochondrial structures [3, 5]. Increased levels of mtROS or/and accumulation of mitochondrial DNA damage may ultimately lead to programmed cell death [6]. Moreover, oxidative stress and mitochondrial dysfunction are linked with several age-related diseases such as diabetes, cancer, cardiovascular, and neurodegenerative diseases [7,8,9,10].

Here, we explored sex-specific associations of PM exposure on expression of mitochondrial associated genes. Furthermore, a pathway analysis was performed on genome wide transcriptome data to investigate whether mitochondrial pathways are highly affected by air pollution exposure. This hypothesis-generating approach identified sex-specific mitochondrial related genes associated with short- and medium-term PM10 exposure that were analysed further in a validation study by means of real-time quantitative PCR (qRT-PCR).


Study design

As our study aim was to investigate the association of short- and medium-term PM10 exposure with mitochondrial-associated transcriptomic responses in peripheral blood, we performed sex-stratified microarray analyses in a discovery cohort of 98 adults. At gene level, we specifically investigated associations between PM10 exposure and expression of 1064 genes listed in the “Human MitoCarta2.0” inventory [11, 12] which are known to encode proteins with mitochondrial localization. Furthermore, we performed pathway analyses on all 15,589 measured transcripts. Based on the gene level and pathway level analyses, we selected 13 MitoCarta genes contributing to top ranked mitochondrial pathways for validation by means of qRT-PCR in an independent validation cohort (n = 169). To substantiate the mitochondrial response to air pollution exposure, we further investigate the link between PM10 exposure and mtDNA content in peripheral blood in the validation cohort.

Study population

Discovery cohort

The original study population is part of the first Flemish Environment and Health Survey (FLEHS I) and consisted of 398 individuals from eight different regions of residence in Flanders (Belgium) [13]. Inclusion criteria were living in the region of Flanders >5 years, age 50 till 65 years and being able to complete questionnaires in Dutch. Informed consent was obtained from all subjects. Sampling took place between September 2004 and June 2005. Participants donated a blood and urine sample, body height and weight were measured in a standardised way. Demographic data, life style factors and health parameters were provided through an extensive self-assessment questionnaire. A subset of 98 non-smokers was selected for whole genome microarray analysis. The selection procedure was previously described by Vrijens et al. [14].

Validation cohort

To validate a selection of MitoCarta genes identified as being associated to PM10 exposure in the discovery cohort, we measured whole blood gene expression levels using qRT-PCR in an independent study population of 169 subjects being part of the third Flemish Environment and Health Survey (FLEHS III). Additionally, mtDNA content was determined in peripheral blood of 150 individuals. Inclusion criteria and data collection were similar as for FLESH I. Informed consent was obtained from all participants. The sampling campaign lasted from May until November 2014.

RNA isolation

Total RNA was extracted from 2.5 ml of whole blood in vacutainers using the Paxgene Blood RNA system (PreAnalytiX, Qiagen, Hilden, Germany), according to the manufacturer’s guidelines. A globin reduction assay (GLOBINclear™ Kit by Ambion, Austin, USA) was performed to remove hemoglobin mRNA from samples assessed in microarray analyses. RNA purity was measured spectrophotometrically and RNA integrity was checked using the BioAnalyzer (Agilent, Palo Alto, USA). Labelled samples were assessed for specific activity and dye incorporation.

Microarray preparation, hybridization and preprocessing

From each sample of the discovery cohort, 0.2 μg total RNA was used to synthesize fluorescent cyanine-3-labeled cRNA following the Agilent one-color Quick-Amp labelling protocol (Agilent Technologies). Samples were hybridized on Agilent Human Whole Genome 4x44K microarrays (design ID 014850). Microarrays signals were detected with an Agilent G2505C DNA Microarray Scanner (Agilent Technologies). Raw data were entered in an in-house developed quality control pipeline in R software applying following preprocessing steps: local background correction, omission of controls, flagging of bad quality spots (based on the size of the spot, the number of pixels per spot, the mean vs. median ratio of the pixel intensity, intensity of spot is not above background, and/or saturation of the spot), and spots with too low intensity, log2-transformation and quantile normalization. Information about the flagging and the R-scripts of the pipeline are available at Further preprocessing included the omission of probes showing >30% flagged data, merging of replicate probes based on median, and the imputation of missing values using K-nearest neighbor imputation (K = 15). If multiple probes represent the same gene, the probe with the largest interquartile range was selected. The final dataset consisted of 15,589 unique Entrez Gene IDs.

Exposure assessment

PM10 and PM2.5 exposure (μg/m3) concentrations were modelled using a spatial temporal interpolation method (Kriging) [15] for each participants’ residential address in combination with a dispersion model. The interpolation method uses land-cover data obtained from satellite images (CORINE land-cover data set) and pollution data collected from a governmental stationary monitoring network. Overall model performance was evaluated by leave-one-out cross-validation including 58 and 34 monitoring points for PM10 and PM2.5 respectively. Validation statistics of the interpolation tool explained >70% of the temporal variability for hourly and annual PM10 and PM2.5 averages in the Flemish Region of Belgium [16]. Coupled with a dispersion model (Immission Frequency Distribution Model, IFDM) [16, 17] that uses emissions from point sources and line sources, this model chain provides PM values in a dense irregular receptor grid. Previous studies conducted a thorough intercomparison of different models currently in use for regulatory purposes in Europe including IFDM [18,19,20,21,22]. To explore potentially critical exposure windows, we averaged residential one week exposure as a proxy for recent exposure, one month exposure as a proxy for medium-term exposure, two-year exposure as a proxy for long-term exposure. Note that PM2.5 exposure estimates were only available for the validation cohort. The Belgian Royal Meteorological Institute provided meteorological data consisting of mean daily air temperature and relative humidity. Apparent temperature was calculated [23, 24] and averaged over the same exposure time window as PM10.

Real-time quantitative PCR (qRT-PCR)

For the validation cohort, total RNA was reverse transcribed to cDNA using the GoScript Reverse Transcription System (Promega, Madison, WI, USA). Gene expression was measured in a 10 μl PCR reaction consisting of 2 μL of a 5 ng/μL dilution of cDNA, TaqMan Fast Advanced Master Mix (Life Technologies, Foster City, CA, USA) and PrimeTimeTM assay (Integrated DNA Technologies, Coralville, IA, USA). Samples were analyzed in triplicate with a 7900HT Fast Real-Time PCR system (Life Technologies, Foster City, CA, USA) applying standard cycling conditions. SDS 2.3 provided threshold cycle (C p ) values which were further processed to normalized relative gene expression values with qBase plus (Biogazelle, Zwijnaarde, Belgium). Replicates were included if the difference in C p values was <0.5. HPRT, IPO8 and YWHAZ were used for data normalization.

DNA extraction and measurement of mtDNA content

For the validation cohort, DNA was isolated from peripheral blood using the QIAmp DNA Mini Kit (QIAGEN GmbH, Hilden, Germany), following the manufacturer’s guidelines. The quantity and purity of the extracted DNA were determined using a Nanodrop spectrophotometer (ND-1000; Isogen Life Science B.V., De Meern, the Netherlands). The DNA samples were diluted to 2.4 ng/μL. MtDNA was measured by calculating the relative ratio of two mitochondrial sequences [MT-ND1 and mitochondrial forward primer from nucleotide 3212 and reverse primer from nucleotide 3319 (MTF3212/R3319)] to a single housekeeping nuclear gene (RPLP0)] measured using a qPCR assay [25]. qPCR was performed using 2.5 μl extracted DNA and 7.5 μl master mix containing Fast SYBR Green dye 2× (Applied Biosystems, Inc., Foster City, California), forward and reverse primers diluted to 300 nM per well, and RNase-free water. Samples were run in triplicate. Each 384-well plate contained 6 interrun calibrators and 2 no-template controls. qPCR was performed using the 7900HT Fast Real-Time PCR System (Life Technologies, Foster City, CA, United States) with following thermal cycling profile: 20 s at 95 °C, followed by 40 cycles of 1 s at 95 °C and 20 s at 60 °C. A melting curve analysis was performed at the end of each run to confirm the absence of nonspecific products. Replicates were included if the difference in C p values was <0.5. qBase software (Biogazelle, Zwijnaarde, Belgium) was used to normalize C p values of the two mtDNA sequences relative to the nuclear gene and to correct for run-to-run differences [26].

Data analysis

Statistical analyses were carried out using SAS software (version 9.3, SAS Institute Inc., Cary, NC, USA). Continuous data were presented as mean and 10th–90th percentiles and categorical data as percentages and frequencies.

Discovery cohort

For each gene, a multivariable linear regression was fitted to investigate the association between log2-transformed gene expression levels and PM10 exposure estimates (short- and medium-term exposure). Previous epidemiological studies reported that environmental stressors have sex-specific immunological responses, with women being more susceptible to smoking than men [27, 28]. Thereupon, in this study we performed sex-stratified analyses to explore both sex-specific and non-specific PM-induced effects and we adjusted for age, body mass index (BMI), socio-economic status (lower secondary or less, higher secondary, higher education), season (medium [April–May, September–November] or cold [December–March]), time of blood sampling (between 08.20 am and 02.30 pm), apparent temperature and microarray batch number (3 scan dates) to correct for batch effects. Of the 15,589 measured genes, 1064 were “Human MitoCarta2.0” genes [11, 12] which are known to encode proteins with mitochondrial localization. Firstly, because of the specific interest in mitochondria, we performed Human MitoCarta gene-wide association scan, with p-values adjusted for multiple testing (for the 1064 genes) by controlling the Benjamini-Hochberg (BH) false discovery rate at 5%. FDR adjusted p-values are referred to as q-values. Secondly, we performed pathway analyses. Of the 15,589 genes, genes with unadjusted p-value <0.05 were uploaded into the online overrepresentation analysis (ORA) tool ConsensusPathDB ( [29] developed at the Max Planck Institute for Molecular Genetics, to identify processes altered by PM10 exposure. Pathways with a p-value <0.05 were considered significant.

Validation cohort

Based on the results of the pathway analyses, 7 MitoCarta genes with a q-value <0.25 (4 in association with short-term exposure and 3 in association with medium-term exposure), contributing to the top 15 ranked pathways/GO terms and with a well-known functional role within mitochondria were selected in women for validation by qRT-PCR. For men, 6 MitoCarta genes (unadjusted p-value <0.05) in relation to medium-term exposure and contributing to the top 15 ranked mitochondrial pathways were selected for validation. We examined the association of expression levels measured by qRT-PCR of the 13 selected genes and of mtDNA content with short-, medium-, and long-term PM10 and PM2.5 exposure. We adjusted for age, BMI, socio-economic status, smoking, temperature and time as well as season of blood sampling, white blood cell (WBC) count and percentage of neutrophils. Residuals were plotted to check whether significance was driven by outlying gene expression values. Over different time windows, for PM10 and PM2.5 separately, the BH multiple testing method was applied to correct for the false discover rate (FDR). For the validation cohort, a q-value <0.05 was considered significant.


The characteristics of the discovery and validation cohort are listed in Table 1 for women and men separately. All participants were of European origin. Both cohorts did not differ in the distribution of age and BMI. For both cohorts, age ranged between 50 and 65 years. BMI averaged (range) 26.6 (20.9–38.5) kg/m2 in the discovery cohort and 25.8 (16.8–39.4) kg/m2 in the validation cohort. Overall, short- and medium-term PM10 exposure estimates were higher in the discovery cohort compared to the validation cohort. In the discovery cohort, more subjects were recruited during the cold period of the year compared to the validation cohort (81.6 vs 39.7%). The discovery cohort consisted only of non-smokers, whilst the validation cohort included 21 (12.4%) smokers. In the validation cohort a higher percentage of participants (53.3%) had a high socio-economic status compared to 28.6% in the discovery cohort.

Table 1 Descriptive characteristics for women and men of the discovery and validation cohort

Gene level analysis

For short- and medium-term exposure, volcano plots of all measured transcripts are presented for both sexes in (Additional file 1: Figure S1). Overall responses to PM10 exposure seem to differ between women and men.

Table 2 lists the top 10 Human MitoCarta genes and their corresponding fold changes for an increase in short-term PM10 exposure of 10 μg/m3 for women and men. For women, 8 genes were significantly associated with short-term PM10 exposure. The top significant gene for women, POLG, encoding the catalytic subunit of the mtDNA polymerase, was downregulated. For men, no significant genes after correction for multiple testing were found. The top ranked gene was IDI1 required in the mevalonate pathway. None of the 8 significant genes in women were in the top 100 of men.

Table 2 Top 10 most strongly associated Human MitoCarta genes with short-term PM10 exposure for women and men

Table 3 list the top 10 mitochondria-localized genes based on their p-value and there corresponding fold changes for an increase in medium-term PM10 exposure of 10 μg/m3 for women and 10 highest ranked genes for men. ALDH7A1 (q-value: 0.21) and MRPS15 (q-value: 0.47) were the top ranked genes for women and men respectively.

Table 3 Top 10 most strongly associated Human MitoCarta genes with medium-term PM10 exposure for women and men

Overrepresentation analysis

Sex-specific PM10 effects were further explored by overrepresentation analyses (ORA). Tables 4 and 5 represent the top 15 significant pathways, with at least 15 measured genes and a total gene size of at most 150 genes, related to, respectively, short- and medium-term PM10 exposure for both sexes. For pathways with the same contributing genes, only the most significant pathway is shown. Mitochondrial pathways, containing mainly MitoCarta genes, are marked with an asterisk. Human MitoCarta genes are indicated in bold font.

Table 4 Top 15 significant pathways associated with short-term exposure
Table 5 Top 15 significant pathways associated with medium-term exposure

In women, mitochondrial translation was the top ranked mitochondrial pathway associated with short-term PM10 exposure of one week before sampling (Table 4). In addition, PM10 exposure was associated with the reaction pathway of busulfan and other DNA damaging agents (p-value: 0.004), by deregulating the expression of pro-apoptotic (e.g. BNIP3, LTBR and BCL2L1, isoform Bcl-xS), anti-apoptotic (e.g. BCL2L1, isoform Bcl-xL), DNA repair (e.g. MLH1), and detoxifying genes (e.g. GSTP1 and GGT1). Of these, BNIP3, BCL2L1 and MLH1 encode proteins (partially) localized in the mitochondria [11, 12].

Mitochondrial GO terms associated with short-term PM10 exposure included “mitochondrial respiratory chain complex I biogenesis” (p-value: 0.001), of which most genes were downregulated, “regulation of mitochondrial membrane permeability” (p-value: 0.0007) playing a crucial role in apoptosis and “mitochondrial genome maintenance” (p-value: 0.026) including genes important for mitochondrial biogenesis and cardiolipin biosynthesis (e.g. STOML2), mtDNA replication (e.g. POLG), mitochondria-mediated apoptosis (e.g. DNAJA3), and unfolded protein response in the mitochondrial matrix (e.g. LONP1). These four Human Mitocarta genes (q-value <0.25) were selected for validation in an independent study population. All were downregulated by medium-term PM10 exposure, except for STOML2.

In men, ORA did not reveal any mitochondrial pathways/GO terms associated with short-term PM10 exposure.

The top 15 significant overrepresented pathways associated to medium-term exposure to PM10 of one month before sampling are listed for women and men in Table 5. In women, top significant mitochondrial processes altered by medium-term PM10 exposure included mitochondrial translation (p-value: 0.001) and the respiratory electron transport chain (p-value: 0.004). A more detailed overview of the respiratory chain is given in Fig. 1. All contributing genes in association to PM10 were upregulated except for ATP5L, a gene encoding a protein of the ATP synthase complex. NDUFA13, UQCRH, and COX7C (q-value <0.25) were selected for further validation.

Fig. 1

Schematic overview of the mitochondrial respiratory electron transport chain and the genes significantly associated with medium-term PM10 exposure per complex in women. Green and red boxes indicate significantly up- and down-regulated genes respectively. Cyt c: cytochrome C

For men, the Tri Carbonic Acid (TCA) cycle (p-value: 0.0004) was the top mitochondrial pathway associated with medium-term PM10 exposure and is represented in Fig. 2. Contributing genes in this pathway were all upregulated. MDH2, IDH2, PC, SUCLA2, SDHA, and ACO2 (p-value <0.05) were validated in an independent study population. Other significant pathways including contributing Human MitoCarta genes, were 3-phosphoinositide biosynthesis, IL6-signaling pathway, and histidine metabolism. Phosphoinositide 3-kinases (PI3K) are crucial for various general cellular processes, including cell survival and apoptosis. The IL6-signaling pathway contains the MitoCarta gene BCL2L1 of which the expression was upregulated in men exposed to relatively high medium-term PM10 exposure. Furthermore, the expression of some mitochondria-localized aldehyde dehydrogenases (ALDH7A1, ALDH1B1), participating in the histidine metabolism pathway, was upregulated by medium-term PM10 exposure in men.

Fig. 2

Schematic overview of the TCA cycle. Green boxes indicate significantly up-regulated genes and their corresponding protein in men in association to short-term PM10 exposure

ConsensuspathDB analyses revealed overrepresented GO terms regarding mitochondrial functioning consistent with the pathway analysis such as the electron transport chain (p-value: 0.008) and mitochondrial translational (p-value: 0.001) in women and the TCA cycle (p-value: 0.009) in men.


For the validation cohort, both PM10 and PM2.5 estimates were available. Results on long-term PM10 exposure in the discovery cohort were published in a previous paper [14]. However, analysis of short- and medium-term exposure with microarray data in the discovery cohort and qPCR validation of MitoCarta genes in an independent cohort is novel. Table 6 present the fold changes (95% CI) and p-values for the linear association between the 7 MitoCarta genes selected for women in the discovery cohort and short-, medium-, and long-term PM10 and PM2.5 exposure in women and men of the validation cohort. For women, several of the selected genes were associated with medium- and/or long-term PM exposure. Of the genes contributing to mitochondrial genome maintenance, expression levels of POLG and LONP1 were negatively associated with long-term PM2.5 exposure (q-value: 0.04 and 0.07 respectively) and DNAJA3 and LONP1 were downregulated by medium-term PM2.5 exposure (q-value: 0.05 and 0.07 respectively). ETC genes were upregulated by PM10 and PM2.5 for all time windows, however only significantly for the association between long-term PM2.5 exposure and COX7C and UQCRH gene expression (q-value <0.05). For men, none of the selected TCA contributing genes could be validated. However, consistent with the observations in women, LONP1 was negatively associated and UQCRH and NDUFA13 were positively associated with long-term PM10 and PM2.5 exposure in men of the validation cohort. However, after FDR correction associations in men were not significant.

Table 6 Association between expression levels of the selected genes and PM exposure in women (n = 94) and men (n = 75) of the validation cohort

Mitochondrial DNA content

Mitochondrial DNA content was negatively associated with short-, medium-, and long-term PM10 and PM2.5 exposure in women (Table 7). For men, mtDNA content was negatively associated with short-term PM10 and PM2.5 exposure whilst medium-term PM10 and PM2.5 exposure revealed a trend towards significance. No significant associations were observed between mtDNA content and the expression of the 13 selected Human MitoCarta genes.

Table 7 Association between mtDNA content and PM exposure


The current study identified several mitochondrial-related genes and pathways significantly associated with fine particle exposure at different exposure time windows: short- (one week before the blood sampling) and medium-term (one month before the blood sampling) PM10 exposure. For women, PM exposure affected, among others, pathways contributing to mitochondrial genome maintenance (short-term), electron transport chains (short-, medium-term) and mitochondrial translation (short- and medium-term). For men, the TCA cycle was positively associated with medium-term PM10 exposure. Furthermore, we were able to validate a selection of mitochondrial-linked genes in an independent study population.

Transcriptome-wide long-term (two-year averages) results of the discovery cohort were described in a previous paper identifying potential gene expression biomarkers of PM exposure [14]. In line with the results for medium-term exposure, the electron transport chain was significantly associated with long-term PM10 and PM2.5 exposure in women.

For women, we selected three genes (COX7C, UQCRH, and NDUFA13), associated with medium-term PM10 exposure in the discovery cohort, encoding proteins contributing to the electron transport chain complexes and four genes (POLG, STOML2, DNAJA3, and LONP1) significantly associated with short-term PM10 exposure in the discovery cohort, of which their corresponding proteins play a role in mitochondrial genome maintenance. For all genes, the direction of association, by PM exposure during the significant time window of the discovery cohort, were replicated in the validation cohort by qPCR. We validated POLG, LONP1, COX7C and UQCRH in relation to long-term PM exposure and DNAJA3 and LONP1, for medium-term exposure. For men, none of the selected TCA contributing genes could be validated. However, consistent with the results for women, UQCRH, and NDUFA13 were upregulated and LONP1 was down-regulated by long-term PM exposure in the validation cohort. Possibly, we could validate most genes only for long-term exposure because this exposure estimate is independent of season of blood sampling, which partly differs between the two study cohorts. For the validation cohort more pronounced effects were observed for PM2.5 compared to PM10, however, for each time window the correlation coefficient between PM2.5 and PM10 was >0.85.

In accordance with our study results, Hoffmann and colleagues observed significantly increased ETC protein levels (Complex II, III, V) in the human bronchial epithelial cell line BEAS-2B in response to cigarette smoke exposure for 6 months [30]. The upregulated expression of ETC genes (Fig. 1), as observed in the current study, and ETC proteins due to environmental toxicant exposure may indicate increased energy demand required to eliminate damage to cellular components.

To further explore mitochondrial responses to PM exposure, we investigated mtDNA content in the validation cohort. In line with the results at the gene expression level (downregulation of genes important for mitochondrial genome replication such as POLG, POLG2, and POLRMT), mtDNA content, measured in the validation cohort, was decreased among both women and men exposed to relatively high PM levels for short-, medium- and long-term. To date, several studies reported a deregulation of mtDNA content in response to environmental factors [3, 25, 31,32,33,34]. However, the direction of effect is not consistent over these studies. Differences in exposure levels and population characteristics often make it difficult to compare study findings. The findings of the current study are in line with previous evidence reporting the ability of ultrafine particles to induce oxidative stress and mitochondrial damage [35], and the selective elimination of damaged mtDNA in order to help maintain mtDNA integrity [4]. Cline hypothesised that poly aromatic hydrocarbons, toxic components of PM, can block mtDNA polymerase and topoisomerase activity and in turn reduce mtDNA replication [3]. In this study, we observe a reduction at the level of gene expression of POLG, POLG2 and POLRMT which may further explain the decreased mtDNA content in individuals exposed to relatively high PM exposure. Presumably, decreased mtDNA content and mitochondrial damage stimulate transcription factors regulating the expression of electron transport genes in order to provide the required energy to eliminate cellular damage. Mitochondrial dysfunction can augment ROS production which may in turn activate the mitochondrial apoptotic pathway. The altered expression of pro- and anti-apoptotic genes (such as gene members of the BCL-2 family and caspases) in response to short-term PM10 exposure in women and of members of the PI3K/AKT (busulfan) pathway, which delivers an anti-apoptotic signal, in men further supports the theory that PM-induced formation of ROS can influence mitochondrial function and regulate cell fate [36]. Fig. 3 shows a schematic overview of the potential effect of PM on mitochondria.

Fig. 3

Schematic representation of the hypothesised pathway by which PM exposure alters mitochondrial functioning and genome maintenance. Within the mitochondria, PM can interact with the electron transport chain inducing increased levels of ROS production. ROS can damage mtDNA leading to further mitochondrial dysfunctioning and ROS production. To repair or eliminate damaged cellular components, the electron transport genes are upregulated to provide the required energy. Elimination of damaged mtDNA and perturbation of mtDNA replication results in a reduction of mtDNA content. Eventually, accumulation of mitochondrial damage can lead to mitochondrial apoptotic signalling. TCA: Tri Carbonic Acid cycle; ETC: electron transport chain; mtDNA: mitochondrial DNA

Overall, differences were observed between the response to PM exposure in men and women. However, for both sexes interacting pathways are altered by PM2.5 exposure; the mitochondrial apoptotic pathway is tightly regulated by several factors such as Bcl-2 family members, altered in women, and the PI3K/AKT (busulfan) pathway modulated in men [37, 38]. Moreover, the TCA cycle, deregulated in men, donates high-energy molecules to the ETC of which genes were differently expressed by PM exposure in women of both cohorts and in men of the validation cohort. In the validation cohort, the effects of air pollution on mtDNA content and expression levels of respiratory electron chain genes and genes contributing to the mitochondrial genome maintenance were more pronounced in women compared to men. Possibly, men are more effectively protected against environmental toxicants and ROS as implied by previous studies [39, 40]. Both studies reported more oxidative damage in female smokers compared to male smokers [39, 40]. In accordance, we observed in men exposed to relatively high PM levels, augmented expression levels of some aldehyde dehydrogenases (ALDH7A1↑, ALDH1B1↑), which convert reactive aldehydes (produced by oxidation of unsaturated fatty acids by ROS) to less toxic products, whilst in women ALDH7A1 was down-regulated.

A strength of our study is that we validated genes in an independent validation cohort by means of qPCR. Moreover, in addition to gene expression, we analyzed mtDNA content in regard to PM exposure in the validation cohort. Our study has some limitations. First, observational studies do not allow to establish causality. Second, PM2.5 estimates were only available for the validation cohort. Third, the large number of tests in combination with the observational study design reduces the power of the transcriptome-wide study. However, this was addressed through focused analyses on MitoCarta genes and mitochondrial pathways using the ORA approach.


Peripheral blood mtDNA content and expression of several genes related to mitochondrial genome maintenance, apoptosis and energy production were altered by PM exposure in a population of healthy middle-aged men and women, potentially reflecting mitochondrial and cellular damage. Future studies at different omics level may further clarify the effect of air pollution on mitochondria functioning and biogenesis.



Apparent temperature


Body mass index


Confidence interval

C p :

Threshold cycle


Electron transport chain


False discovery rate


Flemish environment and health survey


Gene Ontology


Hypoxanthine phosphoribosyltransferase 1


Immission frequency distribution model


Importine 8


Belgian interregional environment agency


Mitochondrial DNA


Mitochondrially encoded NADH dehydrogenase 1


Overrepresentation analysis


Particulate matter

PM10 :

Particulate matter with a diameter < 10 μm

PM2.5 :

Particulate matter with a diameter < 2.5 μm


Quantitative real-time polymerase chain reaction


Real-time quantitative PCR


Reactive oxygen species


Tri carbonic acid


White blood cell


World Health Organization


Tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta


  1. 1.

    Aon MA, Cortassa S, O'Rourke B. Redox-optimized ROS balance: a unifying hypothesis. Biochim Biophys Acta. 2010;1797:865–77.

    CAS  Article  Google Scholar 

  2. 2.

    Lee HC, Wei YH. Mitochondrial biogenesis and mitochondrial DNA maintenance of mammalian cells under oxidative stress. Int J Biochem Cell Biol. 2005;37:822–34.

    CAS  Article  Google Scholar 

  3. 3.

    Cline SD. Mitochondrial DNA damage and its consequences for mitochondrial gene expression. Biochim Biophys Acta. 1819;2012:979–91.

    Google Scholar 

  4. 4.

    Liu P, Demple B. DNA repair in mammalian mitochondria: much more than we thought? Environ Mol Mutagen. 2010;51:417–26.

    CAS  Google Scholar 

  5. 5.

    Paradies G, Petrosillo G, Pistolese M, Ruggiero FM. Reactive oxygen species affect mitochondrial electron transport complex I activity through oxidative cardiolipin damage. Gene. 2002;286:135–41.

    CAS  Article  Google Scholar 

  6. 6.

    Li X, Fang P, Mai J, Choi ET, Wang H, Yang XF. Targeting mitochondrial reactive oxygen species as novel therapy for inflammatory diseases and cancers. J Hematol Oncol. 2013;6:19.

    CAS  Article  Google Scholar 

  7. 7.

    Wallace DC. Colloquium paper: bioenergetics, the origins of complexity, and the ascent of man. Proc Natl Acad Sci U S A. 2010;107(Suppl 2):8947–53.

    CAS  Article  Google Scholar 

  8. 8.

    Sahin E, Colla S, Liesa M, Moslehi J, Muller FL, Guo M, et al. Telomere dysfunction induces metabolic and mitochondrial compromise. Nature. 2011;470:359–65.

    CAS  Article  Google Scholar 

  9. 9.

    Lin MT, Beal MF. Mitochondrial dysfunction and oxidative stress in neurodegenerative diseases. Nature. 2006;443:787–95.

    CAS  Article  Google Scholar 

  10. 10.

    Boland ML, Chourasia AH, Macleod KF. Mitochondrial dysfunction in cancer. Front Oncol. 2013;3:292.

    Article  Google Scholar 

  11. 11.

    Pagliarini DJ, Calvo SE, Chang B, Sheth SA, Vafai SB, Ong SE, et al. A mitochondrial protein compendium elucidates complex I disease biology. Cell. 2008;134:112–23.

    CAS  Article  Google Scholar 

  12. 12.

    Calvo SE, Clauser KR, Mootha VK. MitoCarta2.0: an updated inventory of mammalian mitochondrial proteins. Nucleic Acids Res. 2016;44:D1251–7.

    CAS  Article  Google Scholar 

  13. 13.

    van Leeuwen DM, Gottschalk RW, Schoeters G, van Larebeke NA, Nelen V, Baeyens WF, et al. Transcriptome analysis in peripheral blood of humans exposed to environmental carcinogens: a promising new biomarker in environmental health studies. Environ Health Perspect. 2008;116:1519–25.

    Article  Google Scholar 

  14. 14.

    Vrijens K, Winckelmans E, Tsamou M, Baeyens W, De Boever P, Jennen D et al. Sex-specific Associations between Particulate Matter Exposure and Gene Expression in Independent Discovery and Validation Cohorts of Middle-aged Men and Women. Environ Health Perspect. In Press. DOI: 10.1289/EHP370.

  15. 15.

    Janssen S, Dumont G, Fierens F, Mensink C. Spatial interpolation of air pollution measurements using CORINE land cover data. Atmos Environ. 2008;42:4884–903.

    CAS  Article  Google Scholar 

  16. 16.

    Maiheu B, Veldeman B, Viaene P, De Ridder K, Lauwaet D, Smeets N et al. Identifying the best available large-scale concentration maps for air quality in Belgium. 2012. Available from:

  17. 17.

    Lefebvre W, Degrawe B, Beckx C, Vanhulsel M, Kochan B, Bellemans T, et al. Presentation and evaluation of an integrated model chain to respond to traffic- and health-related policy questions. Environ Model Softw. 2013;40:160–70.

    Article  Google Scholar 

  18. 18.

    Olesen H. The model validation exercise at Mol: overview of results. Int J Environ Pollut. 1995;5:761–84.

    CAS  Google Scholar 

  19. 19.

    Maes G, Cosemans G, Kretzschmar J, Janssen L, Van Tongerloo J. Comparison of six Gaussian dispersion models used for regulatory purposes in different countries of the EU. Int J Environ Pollut. 1995;5:734–47.

    CAS  Google Scholar 

  20. 20.

    Cosemans G, Kretzschmar J, Janssen L, Maes G. The third workshops environmental impact assessment model intercomparison exercise. Int J Environ Pollut. 1995;5:785–98.

    CAS  Google Scholar 

  21. 21.

    Mensink C, Maes G. Comparative sensitivity study for operational short-range atmospheric dispersion models. Int J Environ Pollut. 1996;8:356–66.

    Google Scholar 

  22. 22.

    Cosemans G, Ruts R, Kretzschmar JG. Impact assessment with the Belgian dispersion model IFDM and the New Dutch National Model. Belgitrate: 7th Int Conf on Harmonisation within Atmospheric Dispersion Modelling for Regulatory Purposes; 2001:125–9.

  23. 23.

    Steadman RG. The assessment of sultriness. Part II: effects of wind, extra radiation and barometric pressure on apparent temperature. J Appl Meteorol. 1979;18:874–85.

    Article  Google Scholar 

  24. 24.

    Kalkstein LS, Valimont KM. An Evaluation of Summer Discomfort in the United-States Using a Relative Climatological Index. Bulletin of the American Meteorological Society.  1986;67:842–8.

  25. 25.

    Janssen BG, Munters E, Pieters N, Smeets K, Cox B, Cuypers A, et al. Placental mitochondrial DNA content and particulate air pollution during in utero life. Environ Health Perspect. 2012;120:1346–52.

    CAS  Article  Google Scholar 

  26. 26.

    Hellemans J, Mortier G, De Paepe A, Speleman F, Vandesompele J. qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biol. 2007;8:R19.

    Article  Google Scholar 

  27. 27.

    Faner R, Gonzalez N, Cruz T, Kalko SG, Agusti A. Systemic inflammatory response to smoking in chronic obstructive pulmonary disease: evidence of a gender effect. PLoS One. 2014;9:e97491.

    Article  Google Scholar 

  28. 28.

    Paul S, Amundson SA. Differential Effect of Active Smoking on Gene Expression in Male and Female Smokers. J Carcinog Mutagen. 2014;5:198.

  29. 29.

    Kamburov A, Pentchev K, Galicka H, Wierling C, Lehrach H, Herwig R. ConsensusPathDB: toward a more complete picture of cell biology. Nucleic Acids Res. 2011;39:D712–7.

    CAS  Article  Google Scholar 

  30. 30.

    Hoffmann RF, Zarrintan S, Brandenburg SM, Kol A, de Bruin HG, Jafari S, et al. Prolonged cigarette smoke exposure alters mitochondrial structure and function in airway epithelial cells. Respir Res. 2013;14:97.

    Article  Google Scholar 

  31. 31.

    Bonner MR, Shen M, Liu CS, Divita M, He X, Lan Q. Mitochondrial DNA content and lung cancer risk in Xuan Wei. China Lung Cancer. 2009;63:331–4.

    Article  Google Scholar 

  32. 32.

    Hou L, Zhu ZZ, Zhang X, Nordio F, Bonzini M, Schwartz J, et al. Airborne particulate matter and mitochondrial damage: a cross-sectional study. Environ Health. 2010;9:48.

    Article  Google Scholar 

  33. 33.

    Zhong J, Cayir A, Trevisi L, Sanchez-Guerra M, Lin X, Peng C, et al. Traffic-related air pollution, blood pressure, and adaptive response of mitochondrial abundance. Circulation. 2016;133:378–87.

    CAS  Article  Google Scholar 

  34. 34.

    Pieters N, Koppen G, Smeets K, Napierska D, Plusquin M, De Prins S, et al. Decreased mitochondrial DNA content in association with exposure to polycyclic aromatic hydrocarbons in house dust during wintertime: from a population enquiry to cell culture. PLoS One. 2013;8:e63208.

    CAS  Article  Google Scholar 

  35. 35.

    Li N, Sioutas C, Cho A, Schmitz D, Misra C, Sempf J, et al. Ultrafine particulate pollutants induce oxidative stress and mitochondrial damage. Environ Health Perspect. 2003;111:455–60.

    CAS  Article  Google Scholar 

  36. 36.

    Dagher Z, Garcon G, Billet S, Gosset P, Ledoux F, Courcot D, et al. Activation of different pathways of apoptosis by air pollution particulate matter (PM2.5) in human epithelial lung cells (L132) in culture. Toxicology. 2006;225:12–24.

    CAS  Article  Google Scholar 

  37. 37.

    Engelman JA. Targeting PI3K signalling in cancer: opportunities, challenges and limitations. Nat Rev Cancer. 2009;9:550–62.

    CAS  Article  Google Scholar 

  38. 38.

    Fulda S. Modulation of mitochondrial apoptosis by PI3K inhibitors. Mitochondrion. 2013;13:195–8.

    CAS  Article  Google Scholar 

  39. 39.

    Hakim IA, Harris R, Garland L, Cordova CA, Mikhael DM, Sherry Chow HH. Gender difference in systemic oxidative stress and antioxidant capacity in current and former heavy smokers. Cancer Epidemiol Biomark Prev. 2012;21:2193–200.

    CAS  Article  Google Scholar 

  40. 40.

    Mooney LA, Perera FP, Van Bennekum AM, Blaner WS, Karkoszka J, Covey L, et al. Gender differences in autoantibodies to oxidative DNA base damage in cigarette smokers. Cancer Epidemiol Biomark Prev. 2001;10:641–8.

    CAS  Google Scholar 

Download references


The authors thank all institutions for supporting this project. In particular we would like to acknowledge Environment, Nature and Energy Department of the Flemish Government and Steunpunt Milieu- en Gezondheid.


The project was funded by the Environment, Nature and Energy Department of the Flemish government (LNE/OL201100023/13034/M&G), Steunpunt Milieu- en Gezondheid and European Research Council (ERC-2012-StG 310,898). Ellen Winckelmans has a PhD. fellowship of Hasselt University (BOF program). Karen Vrijens has a Postdoctoral Fellowship of the Research Foundation Flanders (12D7714N).

Availability of data and materials

The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.

Author information




EW, TSN and KV designed the current study within cohorts of the environmental health surveillance program with the help of WB, EDH, MP, HR, GS and NVL. TMDK and JK were responsible for the transcriptome analysis. EW performed the statistical analysis and, with contribution of MT, the bioinformatical analysis. CV and WL did the air pollution modelling. MP performed qPCR. EW wrote the first draft of the manuscript with the help of KV and TSN. All authors were involved in data interpretation and critical revision of the manuscript.

Corresponding author

Correspondence to Tim S Nawrot.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Ethics committee of Antwerp University Hospital (A02–045–46-47; B300201318591), and complies with the Helsinki declaration. All participants provided written informed consent.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional file

Additional file 1: Figure S1.

Volcano plots for short-term exposure (A: women, B: men) and for medium-term PM10 exposure (C: women, D: men). (DOCX 492 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Winckelmans, E., Nawrot, T.S., Tsamou, M. et al. Transcriptome-wide analyses indicate mitochondrial responses to particulate air pollution exposure. Environ Health 16, 87 (2017).

Download citation


  • Ambient air pollution
  • Particulate matter
  • Transcriptome-wide analyses
  • Sex-specific, mitochondria