Particulate matter-induced hypomethylation of Alu and LINE1 in normal human bronchial epithelial cells and epidermal keratinocytes
Genes and Environment volume 44, Article number: 8 (2022)
Airborne particulate matter (PM), a complex mixture of organic and inorganic compounds, is a major public health concern due to its adverse health effects. Understanding the biological action of PM is of particular importance in the improvement of public health. Differential methylation of repetitive elements (RE) by PM might have severe consequences for the structural integrity of the genome and on transcriptional activity, thereby affecting human health. This study aimed to evaluate the effect of inhaled and non-inhaled PM (PM2.5, PM10, and PM10-PAH) exposure on DNA methylation. We quantitatively measured the methylation content of Alu and LINE1 in PM-treated normal human bronchial epithelial cells (NHBE) and normal human epidermal keratinocytes (NHEK) by using whole-genome bisulfite sequencing and pyrosequencing.
All PMs exposure significantly lowered Alu and LINE1 methylation in both cells than in mock-treated controls. Hypomethylation was more prominent in PM10-PAH exposed-NHBE and PM10 exposed-NHEK. Alu and LINE1 methylation change exhibited different sensitivity according to the subfamily evolutionary ages, with stronger effects on the oldest L1-M and Alu J in NHBE, and oldest L1-M and youngest Alu S in NHEK.
These results demonstrate that the differential susceptibility of PM-induced hypomethylation of Alu and LINE1 depends upon RE evolutionary age and PM type.
Air pollution has become an important health concern and ranks as the sixth leading risk factors for premature death globally [1, 2]. Exposure to air pollution is ubiquitous and typically beyond the control of individuals, and the resulting health burden on the population can be high. Particulate matter (PM), one of the most toxic forms of air pollution, is recognized as a major health hazard worldwide, and is associated with respiratory, cardiovascular, and skin diseases [3,4,5]. However, the mechanisms linking PM exposure to adverse health outcomes have not been completely clarified. The size and composition determine the toxicity of the particle . PM consists of a mixture of volatile organic compounds, polycyclic aromatic hydrocarbons (PAH) and inorganic chemicals such as heavy metals that, both individually and together, cause adverse health effects. PM constitutes of microscopic particles of solid or liquid matter suspended in the air. Airborne PM is usually classified as coarse PM10 (<10 mm) and fine PM2.5 (<2.5 mm), depending on the aerodynamic diameter of the particles. PM10 is composed of inhalable particles from dusts, industrial emissions, and traffic emissions. PM2.5 is primarily composed of organic carbon compounds, nitrates, and sulfates. Understanding the biological action of PM is of particular importance in improvement of public health.
Alterations in DNA methylation are associated with various health outcomes, representing an interface between the environment and human disease . Emerging data indicate that PM exposure modulates DNA methylation, a major genomic mechanism of gene expression control, and that these changes might in turn influence inflammation, disease development, and exacerbation risk . However, whether such effects are targeted to specific sites or scattered across the genome globally remain challenging. Alu and long interspersed nucleotide element 1 (LINE1, L1) are significant components of repetitive transposable DNA elements, constituting approximately 17% and 11% of the human genome, respectively , representing as a surrogate marker for genome-wide global methylation levels. Interestingly, transposable repeats are considered as a responser to environmental stressors [10, 11] and their reactivation through hypomethylation can increase genome instability, reactivate lowly expressed genes, or disrupted gene function, thereby potentially contributing to disease-related pathological consequences [12, 13] and provide promising candidate biomarkers for human disease including cancer . Unfortunately, PM-induced DNA methylation of repetitive elements (RE) reported in most previous studies is measured in blood cells, representing overall results from the body organs [15,16,17,18]. However, methylation in skin and lung, major targets of air pollution, has not been well studied. In the present study, we treated the normal human bronchial epithelial cells (NHBE) and normal human epidermal keratinocytes (NHEK) with inhaled and non-inhaled PM (PM2.5, PM10, and PM10-PAH) and then determined the changes in global DNA methylation using whole-genome bisulfite sequencing (WGBS) and pyrosequencing.
Materials and methods
PM2.5, which is a standard diesel PM (SRM1650b) issued by the National Institute of Standard and Technology (Gaithersburg, MD, USA), was bought from Sigma-Aldrich (St. Louis, MO, USA). It was dissolved in dimethyl sulfoxide (DMSO) at 50-mg/ml concentration. PM10-like fine dusts (ERM-CZ100 and ERM-CZ120), which are issued by the European Reference Materials (ERM, Belgium), were brought from Sigma-Aldrich. The former (PM10-PAH) includes several PAHs (benzoanthracene, benzopyrene, benzofluoranthene, and dibenzoanthracene, etc.) in ambient PM10, the latter (PM10) contains heavy metals (arsenic, cadmium, lead, and nickel). They were suspended in phosphate buffered saline (PBS) at 5-mg/ml concentration. PM was prepared just before cell application and sonicated in an ultrasonic bath for 10 min to avoid variability in PM composition and aggregation of particles.
Cell culture and PM treatment
NHBE and NHEK were obtained from the American Type Culture Collection (ATCC, Manassas, VA, USA). NHEK (ATCC PCS-200-011) was grown in Dermal Cell Basal Media (ATCC PCS-200-030) supplemented with Keratinocyte Growth Kit (ATCC PCS-200-040) to propagate in serum-free conditions. NHBE (ATCC PCS-300-011) was cultured in serum-free Airway Epithelial Cell Basal Media (ATCC PCS-300-030) supplemented with Bronchial Epithelial Cell Growth Kit (ATCC PCS-300-040). Both cells were grown at 50% confluency and were treated with PM2.5, PM10, and PM10-PAH for 3 days at a final 50-µg/ml concentration of without a medium change. Cells maintained in culture medium with vehicle (0.1% DMSO or 1% PBS) were used as untreated control groups.
WGBS library preparation and sequencing
The cells were washed with PBS, and genomic DNA was extracted using a QIAamp DNA Mini Kit (Qiagen, Valencia, CA, USA) according to the manufacturer’s instruction. The concentration and quality of the DNA were determined using an Agilent Bioanalyzer 2100 (Agilent Technologies, Santa Clara CA, USA) and agarose gel electrophoresis. DNAs were fragmented using a Bioruptor (Diagenode, Liege, Belgium) to an average size of approximately 250 bp, followed by the blunt ending, 3’-end addition of dA, and adaptor ligation (in this case of methylated adaptors to protect from bisulfite conversion). Ligated DNA was bisulfite-converted using the EZ DNA Methylation-Gold kit (Zymo Research Corp, Irvine, CA, USA). Fragments pf length 200–250 bp were excised from a 2% TAE agarose gel, purified using a QIAquick Gel Extraction Kit (Qiagen), and then amplified via PCR. Libraries were constructed from PCR products with BGI’s DNA nanoball (DNB) technology. The qualified libraries were sequenced using the DNBSEQ®-platform (BGI, Shenzhen, China). Base-calling was performed using the BGISEQ-500 software (v 0.3.8.1111).
Data filtering was conducted using the elimination of contaminating DNA and low-quality reads. Low-quality reads include three types and the reading that accord with one of them will be removed: (1) contain adaptor sequence; (2) N-base number >10%; (3) the number of bases with a quality of <20% and >10% was trimmed. Only clean data was used for further analyses.
Reads mapping and differentially methylated level analysis
Clean reads of each sample were mapped to human UCSC hg19 reference genome using BSMAP software (v2.90) to obtain BAM file. BAM files were sorted and indexed using Samtools software (v0.1.18). The parameters in the mapping and the results of mapping were shown in Supplementary Table 1. Methylation level was determined by dividing the number of reads covering each methylcytosine by the total reads covering the cytosine. MOABS software (v1.3.2) was used to calculate the methylation level of every cytosine in every sample, and to determine differentially methylated cytosine (DMC). After calculating the methylation level, cytosine was considered as “hypomethylated” when the methylation level of cytosine ≤0.2, and Fisher’s Exact Test p-value was <0.05.
Alu and LINE1 methylation analysis
To analyze the methylation level of six evolutionary subfamilies in Alu (Alu Y, Alu S, and Alu J) and LINE1 (L1-H, L1-P, and L1-M), genomic coordinates of all repeats based on hg19 were extracted and obtained from UCSC genome browser using RepeatMasker track as Sae-Lee et al. . All subtypes of 6 target repeats (for example, AluJb, AluJo, AluJr, etc. for AluJ subtypes) were extracted, and methyl-cytosine data within 6 repeats were collected from DMC data of MOABS results in each sample using Python scripts. Data analysis workflow was shown in Supplementary Fig. 1. All CG sites in hg19 and DMCs within six repeats were counted through the samples.
Pyrosequencing of Alu and LINE1
Bisulfite-converted DNA was amplified with PCR primers under previously described conditions . The PCR products were then assayed on the PyroMark Q24 with PyroMark Gold Q24 Reagents (Qiagen) and then analyzed with accompanying software. The degree of methylation was expressed for each DNA locus as a percentage of methylated cytosines over the sum of methylated and unmethylated cytosines. We used non-CpG cytosine residues as built-in controls to verify the bisulfite conversion. In every pyrosequencing run, three controls were included. One well was filled with water to ensure no contamination, and two wells were filled with CpGenome universal methylated and unmethylated DNA (Chemicon, Temecula, CA, USA) as positive and negative control to weigh the repeatability of the assay.
Data are presented as the means ± standard error (SE) of three independent experiments. One-way ANOVA was used for the mean difference test between the groups, and Bonferroni’s correction p-value was used for post-hoc comparison the two groups when the ANOVA was significant. The statistical analyses were conducted using SAS 9.4 (SAS Inc., Cary, NC, USA), and the plots were constructed using R version 4.1.0 (The R Foundation for Statistical Computing, Vienna, Austria). The significance level of the statistical test was set to 5%.
We comprehensively examined the effect of inhaled and non-inhaled PM exposure on RE Alu and LINE1 by WGBS. Based on the reference genome and the UCSC RepeatMasker, approximately 34.6% of all 28 million CpG sites are in Alu (25.0%) and LINE1 (9.6%). The RepeatMasker library mapped 7,040,695 Alu and 2,651,373 LINE1 loci in the UCSC hg19 reference genome assembly, corresponding to 10.1% and 17.1% of the human genome respectively (Table 1). Because human Alu and LINE1 are heavily methylated in normal tissues, all PMs exposure significantly increased the unmethylated CpGs of Alu and LINE1 in both NHBE and NHEK comparing with mock-treated (0.1% DMSO or 1% PBS) cells (Table 1), indicating a potentially defective functionality of these RE. Moreover, PM-induced hypomethylation was prominent in NHBE compared with NHEK as well as in LINE1 compared with Alu (Fig. 1). Although there was a narrow margin between each PM, PM10-PAH and PM10 exhibited the strongest effect on RE hypomethylation in NHBE and NHEK, respectively (Fig. 1). Furthermore, because CpG content and DNA methylation levels dramatically differ across subfamilies, we evaluated the sensitivity of DNA methylation in differentially-evolved Alu and LINE1 subfamilies to different types of airborne PM. We subdivided Alu and LINE1 into three evolutionary subfamilies; oldest Alu J and L1-M, intermediate Alu S and L1-P, and youngest Alu Y and L1-H. Interestingly, Alu J and L1-M showed the strongest hypomethylation in NHBE following treatment with three PMs, whereas Alu Y and L1-M exhibited the strongest hypomethylation in the NHEK (Fig. 2), indicating the association of differential susceptibility of the RE hypomethylation with evolutionary ages of subfamilies. In addition, although pyrosequencing is not expected to comprehensively reflect DNA methylation patterns within individual subfamilies, we have validated the WGBS results by pyrosequencing for accuracy and reproducibility of methylation levels. Significantly reduced methylation of Alu and LINE1 was detectable in both NHBE and NHEK following all PMs exposure (Fig. 3). Likewise, PM10-PAH and PM10 exhibited the strongest effect on RE hypomethylation in NHBE and NHEK, respectively (Fig. 3), showing that a similar trend might be present between WGBS and pyrosequencing methods.
The first major finding of this study was that Alu and LINE1 methylation was significantly lower after PM exposures, providing for the first direct experimental evidence that PM exposure induces DNA hypomethylation in NHBE and NHEK. Recent in vitro experiments have shown that oxidative DNA damage by PM can interfere with the ability of DNA methyltransferase, resulting in RE hypomethylation . DNA methylation is a common feature of eukaryotic genomes and is a core epigenetic process that influences numerous biological processes, such as gene repression, control of cellular development and differentiation, RE silencing, and maintenance of genome stability . DNA methylation mainly changed at locus-specific and genome-wide levels. A number of methods are available for the analysis of global DNA methylation levels . Recently, WGBS has revolutionized the way of interrogating the methylome to realize genome-wide methylation analysis at a single-base resolution . The current study determined uniquely mappable WGBS data to be the most reproducible and accurate measurement of global DNA methylation levels by comparing with pyrosequencing assays of RE, providing WGBS as the gold standard method in methylomics for its unsurpassed resolution and coverage.
Alu is the largest family of short interspersed nuclear elements in the human genome and LINE1 is a predominant member of LINEs . The former is a non-autonomous, transposable element (TE) to be mobilized in trans by LINE1, but the latter is only autonomous TE. Moreover, they constitute the critical regulators of genetic information expression by providing regulatory sequences or introducing alternative start or stop codons into functional genes , providing to act as global modifiers of gene expression through changes in their own methylation state. Accordingly, growing evidence has shown that the altered methylation states of TEs might be associated with aging, autoimmune diseases, cardiovascular disease, or cancer development and progression , suggesting that these changes are not the simple consequences of the disease, but may often drive the pathogenesis. Interestingly, Alu and LINE1 initiate the spread of CpG island (CGI) methylation and the CGI length is associates with their distribution , indicating the potential centers for de novo methylation events. Unfortunately, current investigations have focused only on analyzing a single common sequence for Alu and LINE1 through pyrosequencing assays, which is easier to do than previous methods to quantify total genomic 5-methylcytosine . Moreover, recent reports have shown that the methylation of a common sequence is not correlated with global methylation content in normal tissues and that CpG content is a primary determinant of changes over time in DNA methylation at individual CpG sites [29, 30].
The second novel finding of the current study was that the evolutionary age of RE subfamilies determined differential susceptibility of DNA hypomethylation to ambient PM. Sparse data are available on the effects of environmental exposures across different subfamilies of TEs. Based on the peak period of amplification and the level of nucleotide substitutions, Alu and LINE1 are subdivided into each subfamily with different evolutionary ages; young (Alu Y and L1-H), intermediate (Alu S and L1-P), and old (Alu J and L1-M) subfamilies [30, 31]. Because of frequent deamination of methylated cytosines in CpG dinucleotides, older subfamilies remain less rich in CpG sites and show weaker or no transposon activity, whereas young subfamilies are richer in CpGs and still transcriptionally active in the human genome. Moreover, CpG content and DNA methylation levels vary dramatically across subfamilies. Interestingly, Byun et al.  have demonstrated that the effect of PM10 exposure on DNA methylation depends on the subfamily evolutionary age, with a stronger negative effect on older LINE1 and younger Alu. Recently, older subfamilies of Alu and LINE1 elements (Alu J and L1-M) exhibit great hypomethylation in chronic lymphocytic leukemia . Furthermore, the Alu Y sequence shows remarkable differences in DNA methylation state across colorectal cancer drug resistance . Taken together, these results suggest that the evolutionary age of TE subfamilies might determine differential vulnerability of DNA methylation to environmental exposures.
The present study showed PM-induced hypomethylation of Alu and LINE1 elements with differential susceptibility of the evolutionary subfamily, suggesting that RE hypomethylation might be a vital mechanism underlying the harmful effects of airborne PM and that monitoring of the methylation status for a specific subset of RE could serve as interface sensors between PM and DNA methylation. Furthermore, these results could provide a better understanding of the effects of PM exposure on RE subfamilies and the role of RE in response to environmental risk factors related to human health and disease. This study is the first to utilize the WGBS platform in the analysis of the subfamily-specific methylation of RE in PM-exposed human skin and lung tissues. However, further work to analyze locus-specific hypomethylation of Alu and LINE1 remains challenging.
Availability of data and materials
All data generated or analyzed during this study are included in this published article.
Differentially methylated cytosine
Long interspersed nucleotide element 1
Normal human bronchial epithelial cells
Normal human epidermal keratinocytes
Polycyclic aromatic hydrocarbons
Whole-genome bisulfite sequencing
Brunekreef B, Holgate ST. Air pollution and health. Lancet. 2002;360:1233–42.
Schraufnagel DE, Balmes JR, Cowl CT, De Matteis S, Jung SH, Mortimer K, et al. Air pollution and noncommunicable diseases: A review by the forum of International Respiratory Societies’ Environmental Committee, Part 2: Air pollution and organ systems. Chest. 2019;155:417–26.
Samet J, Krewski D. Health effects associated with exposure to ambient air pollution. J Toxicol Environ Health A. 2007;70:227–42.
Ruckerl R, Schneider A, Breitner S, Peters A. Health effects of particulate air pollution: A review of epidemiological evidence. Inhal Toxicol. 2011;23:555–92.
Kim KE, Cho D, Park HJ. Air pollution and skin disease: Adverse effects of airborne particulate matter on various skin diseases. Life Sci. 2016;152:126–34.
Valavanidis A, Fiotakis K, Vlachogianni T. Airborne particulate matter and human health: Toxicological assessment and importance of size and composition of particles for oxidative damage and carcinogenic mechanisms. J Environ Sci Health C Environ Carcinog Ecotoxicol Rev. 2008;26:339–62.
Ozanne SE, Constancia M. Mechanisms of disease: the development origins of disease and the role of the epigenotype. Nat Clin Prat Endocrinol Metab. 2007;3:539–46.
Rider CF, Carlsten C. Air pollution and DNA methylation: effects of exposure in humans. Clin Epigenetics. 2019;11:131.
Lander ES, Linton LM, Birren B, Nusbaum C, Zody MC, Baldwin J, et al. Initial sequencing and analysis of the human genome. Nature. 2001;409:860–921.
Mioussea IR, Chalbota MG, Lumenb A, Fergusona A, Kavourasa IG, Koturbasha I. Transposable elements in response to environmental stressors. Mutat Res Rev Mutat Res. 2015;765:19–39.
Mirouze M, Lanciano S. Transposable elements: All mobile, all different, some stress responsive, some adaptive. Curr Opin Genet Dev. 2018;49:106–14.
Kazazian HH, Moran JV. Mobile DNA in health and disease. N Engl J Med. 2017;377:361–70.
Burns KH. Our conflict with transposable elements and its implications for human disease. Annu Rev Pathol. 2020;15:51–70.
Ponomaryova AA, Rykova EY, Gervas PA, Cherdyntseva NV, Mamedov IZ, Azhikina TL. Aberrant methylation of LINE-1 transposable elements: A search for cancer biomarker. Cells. 2020;9:2017.
Tarantini L, Bonzini M, Apostoli P, Pegoraro V, Bollati V, Marinelli B, et al. Effects of particulate matter on genomic DNA methylation content and iNOS promoter methylation. Environ Health Perspect. 2009;117:217–22.
Madrigano J, Baccarelli A, Mittleman MA, Wright RO, Sparrow D, Vokonas PS, et al. Prolonged exposure to particulate pollution, genes associated with glutathione pathways, and DNA methylation in a cohort of older men. Environ Health Perspect. 2011;119:977–82.
De Prins S, Koppen G, Jacobs G, Dons E, Van de Mieroop E, Nelen V, et al. Influence of ambient air pollution on global DNA methylation in healthy adults: A seasonal follow-up. Environ Int. 2013;59:418–24.
Ding R, Jin Y, Liu X, Zhu Z, Zhang Y, Wang T, et al. Characteristics of DNA methylation changes induced traffic-related air pollution. Mutat Res Genet Toxicol Environ Mutagen. 2016;796:46–53.
Sae-Lee C, Biasi JD, Robinson N, Barrow TM, Mathers JC, Koutsidis G, et al. DNA methylation patterns LINE-1 and Alu for pre-symptomatic dementia in type 2 diabetes. PLoS One. 2020;15:e0234578.
Kim KY, Kim DS, Lee SK, Lee IK, Kang JH, Chang YS, et al. Association of low-dose exposure to persistent organic pollutants with global DNA hypomethylation in healthy Koreans. Environ Health Perspect. 2010;118:370–4.
Valinluck V, Tsai HH, Rogstad DK, Burdzy A, Bird A, Sowers LC. Oxidative damage to methyl-CpG sequence inhibits the binding of the methyl-CpG binding domain (MBD) of methyl-CpG binding protein 2 (MeCP2). Nucleic Acids Res. 2004;32:4100–8.
Li E, Zhang Y. DNA methylation in mammals. Cold Spring Harb Perspect Biol. 2014;6:a019133.
Terry MB, Delgado-Cruzatak L, Vin-Raviv N, Wu HC, Santella RM. DNA methylation in white blood cells: Association with risk factors in epidemiologic studies. Epigenetics. 2011;6:828–37.
Li Q, Hermanson PJ, Springer NM. Detection of DNA methylation by whole-genome bisulfite sequencing. Methods Mol Biol. 2018;1676:185–96.
Ali A, Han K, Liang P. Role of transposable elements in gene regulation in the human genome. Life (Basel). 2021;11:118.
Wilson AS, Power BE, Molloy PL. DNA hypomethylation and human diseases. Biochim Biophys Acta. 2007;1775:138–62.
Kang MI, Rhyu MG, Kim YH, Jung YC, Hong SJ, Cho CS, et al. The length of CpG islands is associated with the distribution of Alu and L1 retroelements. Genomics. 2006;87:580–90.
Yang AS, Estecio MR, Doshi K, Kondo Y, Tajara EH, Issa JP. A simple method for estimating global DNA methylation using bisulfite PCR of repetitive DNA elements. Nucleic Acids Res. 2004;32:e38.
Wolff EM, Byun HM, Han HF, Sharma S, Nichols PW, Siegmund KD, et al. Hypomethylation of a LINE-1 promoter activates an alternate transcript of the MET oncogene in bladders with cancer. PLoS Genet. 2010;6:e1000917.
Byun HM, Nordio F, Coull BA, Tarantini L, Hou L, Bonzini M, et al. Temporal stability of epigenetic markers: Sequence characteristics and predictors of short-term DNA methylation variations. PLoS One. 2012;7:e39220.
Ohshima K, Hattori M, Yada T, Gojobori T, Sakaki Y, Okada N. Whole-genome screening indicates a possible burst of formation of processed pseudogenes and Alu repeats by particular L1 subfamilies in ancestral primates. Molecular reconstitution of extinct LINE-1 elements. Genome Biol. 2003;4:R74.
Giordano J, Ge Y, Gelfand Y, Abrusan G, Benson G, Warburton PE. Evolutionary history of mammalian transposons determined by genome-wide defragmentation. PLoS Comput Biol. 2007;3:e137.
Byun HM, Motta V, Panni T, Bertazzi PA, Apostoli P, Hou L, et al. Evolutionary age of repetitive element subfamilies and sensitivity of DNA methylation to airborne pollutants. Part Fibre Toxicol. 2013;10:28.
Barrow TM, Doo NW, Milne RL, Giles GG, Willmore E, Strathdee G, et al. Analysis of retrotransposon subfamily DNA methylation reveals novel early epigenetic changes in chronic lymphocytic leukemia. Haematologica. 2021;106:98–110.
Lin X, Stenvang J, Rasmussen MH, Zhu S, Jensen NF, Tarpgaard LS, et al. The potential role of Alu Y in the development of resistance to SN38 (Irinotecan) or oxaliplatin in colorectal cancer. BMC Genom. 2015;16:404.
Authors thank Dr. Yong Chul Boo for his critical reading of this manuscript.
This research was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (No. NRF-2020R1I1A3067217) and by a grant of the Korea Health Technology R&D Project through the Korea Health Industry Development Institute (KHIDI) funded by the Ministry of Health & Welfare, Republic of Korea (grant number: HP20C0004).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no conflict of interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Fig. 1. Graphical overview and representative of Alu and LINE1 data analysis. (A) Workflow chart. (B) Analysis of methylation index in every CG-cytosine using MOABS. (C) Extraction of repetitive elements including Alu and LINE1. (D) Six repeats (AluJ, AluS, AluY, LiH, L1M, and L1P) extraction using python script. (E) CGs extraction within 6 repeats.
Overview of WGBS data analysis.
About this article
Cite this article
Lee, J.Y., Lee, W.K. & Kim, D.S. Particulate matter-induced hypomethylation of Alu and LINE1 in normal human bronchial epithelial cells and epidermal keratinocytes. Genes and Environ 44, 8 (2022). https://doi.org/10.1186/s41021-022-00235-4
- Particulate matter
- DNA methylation