Original Research ARTICLE
Blood Transcriptome Response to Environmental Metal Exposure Reveals Potential Biological Processes Related to Alzheimer's Disease
- 1Department of Toxicogenomics, Maastricht University, Maastricht, Netherlands
- 2Section of Sustainable Health, Department of Public Health and Clinical Medicine, Umeå University, Umeå, Sweden
- 3Odontology, Umeå University, Umeå, Sweden
- 4Cancer Risk Factors and Life-Style Epidemiology Unit, Institute for Cancer Research, Prevention and Clinical Network - ISPRO, Florence, Italy
- 5Division of Occupational and Environmental Medicine, Lund University Hospital, Lund, Sweden
- 6Institute of Chemical Biology, National Hellenic Research Foundation, Athens, Greece
Alzheimer's disease (AD) is a neurodegenerative disease which is manifested by a progressive and irreversible decline of cognition, memory loss, a shortened attention span, and changes in personality. Aging and genetic pre-dispositions, particularly the presence of a specific form of apolipoprotein E (APOE), are main risk factors of sporadic AD; however, a large body of evidence has shown that multiple environmental factors, including exposure to toxic metals, increase the risk for late onset AD. Lead (Pb) and cadmium (Cd) are ubiquitous toxic metals with a wide range of applications resulting in global distribution in the environment and exposure of all living organisms on earth. In addition to being classified as carcinogenic (Cd) and possibly carcinogenic (Pb) to humans by the International Agency for Research on Cancer, both compounds disrupt metal homeostasis and can cause toxic responses at the cellular and organismal levels. Pb toxicity targets the central nervous system and evidence for that has emerged also for Cd. Recent epidemiological studies show that both metals possibly are etiological factors of multiple neurodegenerative diseases, including Alzheimer's disease (AD). To further explore the association between metal exposure and AD risk we applied whole transcriptome gene expression analysis in peripheral blood leukocytes (PBLs) from 632 subjects of the general population, taken from the EnviroGenomarkers project. We used linear mixed effect models to associate metal exposure to gene expression after adjustment for gender, age, BMI, smoking, and alcohol consumption. For Pb exposure only few associations were identified, including a downregulation of the human eukaryotic translation initiation factor 5 (eIF5). In contrast, Cd exposure, particularly in males, revealed a much stronger transcriptomic response, featuring multiple pathways related to pathomolecular mechanisms of AD, such as endocytosis, neutrophil degranulation, and Interleukin−7 signaling. A gender stratified analysis revealed that the Cd responses were male-specific and included a downregulation of the APOE gene in men. This exploratory study revealed novel hypothetical findings which might contribute to the understanding of the neurotoxic effects of chronic Pb and Cd exposure and possibly improve our knowledge on the molecular mechanisms linking metal exposure to AD risk.
Alzheimer's disease (AD) is a progressive neurodegenerative disease and the most common form of dementia. Currently, close to 50 million people are afflicted with dementia across the globe, and this number is estimated to triple by 2050 (1). AD is manifested in a progressive and irreversible decline of cognition, memory loss, a shortened attention span and changes in personality (2).
Despite decades of research, the etiopathology of AD remains largely unknown, and no effective treatment is currently available. Pathological hallmarks include a degenerative process of damaging cholinergic neurons across brain regions that serve important functional roles in conscious awareness, attention, and memory; extracellular aggregates of amyloid beta (Aβ) peptides resulting in disrupted cell functioning; intracellular accumulation of hyperphosphorylated tau, causing synaptic dysfunction, along with glial activation, and inflammation (2, 3).
A rare form of AD is early-onset, or familial AD, which accounts for only 2–10% of total AD cases and occurs at ages below 65. Genome-wide association studies have revealed three high-penetrance mutations to cause autosomal dominant AD: Amyloid precursor protein (APP), Presenilin 1 (PSEN1), and Presenilin 2 (PSEN2) (4). However, the vast majority of patients suffer from late onset, or sporadic AD, with clinical symptoms only occurring at an age above 65 years. Sporadic AD is a multifactorial disease with aging and genetic pre-dispositions, particularly the presence of the apolipoprotein E ε4 (APOE ε4) allele, being the main risk factors. Moreover, epidemiological studies have also identified multiple environmental risk factors, including exposure to toxic metals, as drivers in late onset AD (4–7).
Humans are exposed to toxic metals through inhalation of contaminated air, dermal absorption of metals contained in soil, and ingestion of contaminated water and food (7). Lead (Pb) and Cadmium (Cd) are naturally occurring metals that are widely used in e.g., industrial and domestic applications. They also occur in some agricultural applications. This has resulted in a global distribution in the environment. Even at low levels of exposure, Pb, and Cd can induce multiple organ damage, and therefore present a major human health threat. Cd has been classified as carcinogenic and Pb as possibly carcinogenic to humans by the International Agency for Research on Cancer (8). Next to the carcinogenic effects, a large body of evidence has emerged that Pb and Cd toxicity also target the central nervous system (CNS), causing neurodegenerative effects (9), and increasing evidence has demonstrated that both Pb and Cd exposure are possible etiological factors for sporadic AD. Population-based research has shown that Pb, and possibly also Cd, exposure results in AD-like pathology such as memory loss and deficits in intelligence, attention, language, and emotion (10, 11).
Once in the blood, Pb is distributed throughout the body, thereby making it available to other tissues (12). In the brain vessels Pb substitutes for calcium ions, and can therefore rapidly cross the blood-brain barrier (13). Brain Pb levels alter neuronal differentiation and cause severe damage to the brain (14). Even low Pb levels have been shown to cause neurotransmitter alterations, possibly resulting in malfunctioning of the GABAergic, dopaminergic, and cholinergic systems (15). Intracellularly, Pb replaces calcium, and other essential metals, and disrupts the corresponding biometal-dependent mechanisms (12). Pb exposure was associated to increased levels of Aβ peptides in animal and cell culture studies (16, 17). In rodents, increased Pb levels were also associated to elevated expression of APP as well as increased production of Aβ peptides (18). A case study has reported that Pb poisoning in childhood has led to neurofibrillary tangles in the brain (19). A study on 55 young adults who had participated as newborns in a prospective cohort study showed inverse correlations of umbilical cord Pb concentration and expression of potential AD biomarkers e.g., ADAM9, RTN4, and LRPAP1 genes, hence demonstrating that early-life Pb exposure influences biological processes involved in AD pathogenesis (20).
Chronic Cd exposure is associated with multiple health outcomes including hypertension, kidney dysfunction, and also neurological diseases (12). Cd has been shown to cross the blood-brain barrier, and has even been shown to alter its permeability (21). The mechanisms underlying Cd neurotoxicity are not completely understood, however, they include oxidative damage, inflammation, neuronal apoptosis, and possibly interaction with other metals such as cobalt and zinc (12, 22). Other studies have shown that Cd also interacts with neurotransmitters, e.g., causing a decrease in exocytotic release of glutamate (23), possibly by blocking the influx of Ca2+ through membrane channels (24).
The potential role of Cd in sporadic AD arises from the presence of significantly higher Cd concentrations in brain tissue of AD patients as compared to age-matched control subjects (25). In vitro experiments have illustrated that Cd potentially causes self-aggregation of the tau peptides (26), and may directly interact with Aβ protein (27).
Evidently, humans are not exposed to single metals but rather to mixtures, which makes it even more complicated to identify molecular mechanisms triggered by individual compounds. Several studies have shown that chemicals, including heavy metals, can enhance the effect of other chemicals, so that as a mixture they can exert a larger effect (28). This synergistic effect for Pb and Cd has been shown to occur in model organisms (29). In addition it has been suggested that high Pb levels affect mental and psychomotor development in children exposed to high prenatal Cd levels (30).
Together, these studies suggest potential implications of Pb and Cd in the pathomolecular mechanisms which might increase sporadic AD risk, although further research is needed to better understand the neurotoxic mechanisms of chronic low-level Pb and Cd exposure, as well as its potential role in the etiopathology of AD. In the present study we assume that Pb and Cd induce responses in blood cells that are similar to responses in the CNS. Previous studies have shown that differentially expressed genes in the brain and blood revealed a significant overlap of gene expression patterns (31, 32). Therefore, we have applied whole transcriptome gene expression analysis in blood cells from 632 apparently healthy subjects chronically exposed to Pb and Cd from the EnviroGenomarkers project (http://www.envirogenomarkers.net).
Materials and Methods
Selection of the Population
To investigate the effect of metal exposure on disease risk we performed a global transcriptome analysis in archived peripheral blood leukocytes (PBLs) samples of 632 and exposure/health data were derived from -at that time- generally healthy subjects from the EnviroGenomarkers Project. The EnviroGenomarkers project was designed as a nested case-control study, including controls and “future disease” subjects developing either lymphoma or breast cancer after a minimum 2 years from sample collection. The project is based on two prospective cohort studies that have previously been described: The Northern Sweden Health and Disease Study (NSHDS) (33), in which its subcohort the Vasterbotten Intervention Programme (n = 406) (34) was used, and EPIC-Italy (n = 226) (35). For the present study we analyzed blood samples from 632 participants, including 393 females and 239 males (Table 1). Smoking status was assessed at the time of sampling. The EnviroGenomarkers project was approved by the Regional Ethical Review Board in Umea and the Florence Health Unit Local Ethical Committee and all methods were carried out in accordance with the approved guidelines. All participants gave written informed consent.
Table 1. Population characteristics. Smoking status was self-reported at the time of blood sampling.
Erythrocyte concentrations of Pb and Cd were determined by means of inductively coupled plasma–mass spectrometry (ICP–MS) (36) at Lund University Hospital, Lund, Sweden, for both the EPIC–Italy and NSHDS samples. The limit of detection, calculated as three times the standard deviation of a blank, was 0.03 μg/L for cadmium and 0.09 μg/L for lead in EPIC–Italy and NSHDS samples, as described earlier (37). In addition to the numerical variable for each metal we calculated a ZScore for each subject as a representative of the internal exposure to the mixture of both metals. This ZScore is centered on the population mean and therefore negative values correspond to the raw scores below the population mean, while positive values correspond to the raw scores above the population mean. The ZScore was defined as the mean of the ZScores for each metal: ZScore = (X–μ)/σ (where X represents the value of the subject, μ the mean and σ the standard deviation of the population) (38–40).
PBLs for RNA extraction were taken at the same time as samples for measurements of metal concentrations. Total RNA was extracted from PBLs with the RiboPure™ Blood kit (Ambion, Austin, TX, United States). Genome-wide analysis of gene expression (Agilent 4 × 44K human whole genome microarray platform) and the corresponding data quality assessment and preprocessing, were conducted as described previously (41, 42). Briefly: RNA samples were reverse transcribed into cDNA and labeled with cyanine 3 prior to hybridization. Subsequently, slides were washed and scanned using the Agilent Technologies G2565CA DNA Microarray Scanner. Quality control was established by visual evaluation of the scan images before and after within- and between-array normalization (using LOESS and A-quantile, respectively) (42). Probes with <75% of the maximum possible number of pixels were left out. Imputation of the missing values was done using k nearest neighbors approach (41). Resulting signal intensities were log2 transformed before proceeding.
Statistical Data Analysis
Linear mixed models (LMM) were used to identify associations between the gene expressions and the individual Pb and Cd levels as well as the cumulative ZScore. The analyses were performed using the R package “lme4,” correcting for the potential confounders gender, age, BMI, smoking and alcohol intake. Furthermore, the technical variations introduced by the different batches of RNA isolation, RNA labeling and hybridization were added as random effects. Covariates included in the model were selected by the following method: Using the R package “lmer/lm” we performed a model m0 explaining the variation in the transcriptome by exposure, and a model m1 explaining the variation in the transcriptome by exposure and controlling for a single covariate. Both models were subjected to an anova test and resulting p-values were plotted in a pqq plot. If the distribution of p-values from the respective anova deviated significantly from the normal distribution a covariate was included in the final model (Supplementary Figure 1). Scaled, numerical Pb, Cd and the cumulative ZScore were defined as the respective variables of interest. Resulting p-values were controlled for by the False Discovery Rate (FDR) set at 5 and 20%. Reactome and KEGG pathway analyses were performed using the R package “enrichPlot.” Furthermore, we have used Cytoscape (version 3.8.0) and the Reactome FI package to perform a network analysis, pathway enrichment analysis and gene set enrichment analysis (GSEA). The GeneScore used for the GSEA and network analysis was calculated by multiplying the β-coefficient and –log10 of the corresponding p-value. A full KEGG pathway analysis was done using the web interface of ConsensusPathDB. To test the stability of the model we performed modified k-means cross validation. Briefly: 10% of the subjects were randomly left out at a time; LMM was performed on the remaining 90% subjects. This step was repeated for 100 times. At the end, we compared all resulting 100 coefficient vectors with the coefficient vector of the full dataset set including all subjects.
Population and Exposure
To investigate the effect of metal exposure on disease risk we performed a global transcriptome analysis in PBL samples of 632 subjects from the EnviroGenomarkers study population of 393 females and 239 males (Table 1).
Next to the internal exposure to the individual pollutants Pb and Cd we also attempted to investigate the potential synergistic effect using a cumulative ZScore of both metals. Across the cohorts we found significantly higher exposure levels in EPIC-Italy for ZScore and Pb when compared to NHDS. For cadmium no significant differences were observed (Table 2). Males showed significantly higher levels of Pb exposure, while females had higher Cd levels. The ZScore, however, did not show a significant difference between genders (Table 3).
Transcriptome-Wide Metal Exposure Associations
The LMM identified the cumulative ZScore of metal exposure to be significantly associated to 16 Agilent Probe IDs (AgIDs) when controlling the FDR at 5 percent, and 954 AgIDs when controlling the FDR at 20 percent. These corresponded to 16 (14↑, 2↓) and 767 (734↑, 33↓) transcripts, respectively (Figure 1A). The individual pollutant Pb was significantly associated with only 1 AgID when controlling the FDR at 5 percent, and 5 AgIDs when controlling the FDR at 20 percent. These corresponded to 1 (↓) and 4 (2↑, 2↓) transcripts, respectively (Figure 1B). For Cd the LMM identified 20 AgIDs when controlling the FDR at 5 percent, and 1,039 AgIDs when controlling the FDR at 20 percent. These corresponded to 17 (17↑) and 792 (606↑, 186↓) transcripts, respectively (Figure 1C). The cumulative ZScore showed 3 overlapping transcripts with Pb and 389 with Cd (FDR at 20%). Consequently, 374 transcripts were only identified with the ZScore and are therefore potentially associated to the synergistic effect of Pb and Cd. Since relatively few associations were identified, even taken the low significance threshold of FDR <20%, we hypothesized that gender-differences might have an impact on the analysis. Consequently, we performed a subsequent analysis on the data stratified for males and females.
Figure 1. LMM results for the full dataset. These volcano plots show the resulting p-values of the LMM when associating the AgIDs to the cumulative metal ZScore (A), individual Pb (B) and Cd (C) exposure. The red color indicates transcripts significant with a FDR <20%; labeled transcripts remained significant at FDR < 0.05.
When stratifying the dataset by gender we observed the following associations: the cumulative ZScore was associated to 3 and 5 AgIDs in the 393 females when controlling the FDR at 5 and 20 percent, respectively. These corresponded to 3 (1↑, 2↓) and 4 (1↑, 3↓) transcripts, respectively (Figure 2A). For the 239 males we identified 9 and 235 associated AgiDs which corresponded to 9 (8↑, 1↓) and 199 (198↑, 1↓) transcripts, respectively (Figure 2B). The individual pollutant Pb was significantly associated with 7 and 12 AgIDs in females when controlling the FDR at 5 and 20 percent. These corresponded to 7 (2↑, 5↓) and 11 (3↑, 8↓) transcripts, respectively (Figure 2C). For males, however, we identified only 1 (1↓) AgID/transcript at both significance thresholds (Figure 2D). For the individual Cd we identified only 1 (1↑) AgID/transcript in females at both significance thresholds (Figure 2E). In contrast, we identified 2 and 4204 AgIDs in males when controlling the FDR at 5 and 20 percent. These corresponded to 2 (2↑ and 2862 (1,662↑, 1,200↓) transcripts, respectively (Figure 2F). The strong male-specific response to Cd exposure, observed at the lower significance threshold of FDR <20%, prompted us to focus the following analysis of potentially perturbed pathways in males only.
Figure 2. LMM results for the gender stratified dataset. These volcano plots show the resulting p-values of the LMM when associating the AgIDs to the cumulative ZScore for females (A) and males (B), individual Pb for females (C) and males (D) and Cd for females (E) and males (F) exposure. The red color indicates transcripts significant with a FDR <0.2; labeled transcripts remained significant at FDR <0.05.
To validate the stability of the results we performed a modified k-means cross validation (see methods). The cross validation revealed that the results of the cumulative ZScore and Cd showed high stability, while the results of Pb only resulted in moderate stability for full and stratified analysis (Supplementary Figure 1). A stratified LMM analysis for each individual cohort was implemented; however, only few significant associations were found as a consequence of the smaller sample size. Briefly; for the ZScore 27 and 1 AgIDs, and for Pb 3 and 3 AgIDs were identified for EPIC and NHDS, respectively (FDR at 20%). No associations for Cd were identified (FDR at 20%).
To test the impact of a future cancer on our findings we stratified the data to only controls. Despite the small sample size (n = 313) we identified 3, 2, and 763 AgIDs to be associated with the ZScore, Pb, and Cd, respectively. Of these findings in controls 2, 1, and 387 AgIDs were overlapping with the associations identified in the full dataset. Stratified analysis for future lymphoma (n = 239) and breast cancer (n = 80) cases returned no significant associations for Cd and were deemed underpowered.
Since smoking is a major source of cadmium exposure we also stratified our data into never-smokers, non-smokers, and smokers. For never- and non-smokers we did not find any significant association, while for the smokers we only observed 12 positively associated AgIDs (FDR at 20%).
Perturbed Disease Pathways and Network Analysis
The strong response to Cd in males prompted us to focus the downstream analysis of perturbed pathways solely on this gender-specific effect. Pathway analysis of all Cd-associated transcripts showed an overrepresentation of genes involved in Reactome pathways such as interleukin signaling, clathrin-mediated endocytosis, and neutrophil degranulation (Figure 3A); KEGG pathways included chemokine signaling (Figure 3B). Figures only include the top 20 significant pathways. Additionally we attached a full list of enriched Reactome and KEGG pathways, and a GSEA (Supplementary Tables 1–3). Interestingly, taken the known accumulation of Cd in kidney, the KEGG pathways analysis revealed a significant enrichment of genes involved in renal cell carcinoma (Figure 3, Supplementary Table 2). When feeding only positively associated transcripts into the analysis we also identified Reactome pathways among others related to hemostasis, rho GTPase, and membrane trafficking (Figure 3A); KEGG pathways also featured such as endocytosis, neurotrophin signaling pathways, and Toll-like receptor signaling (Figure 3B). Inversely associated transcripts identified Reactome pathways such as related to GPCR and Keratinization (Figure 3A); KEGG pathways also featured such as neuroactive ligand-receptor interaction (Figure 3B).
Figure 3. Reactome and KEGG overrepresentation analysis. Figure shows enriched Reactome (A) and KEGG (B) pathways for genes significantly associated to Cd in males (FDR <0.2). Each figure shows enriched pathways for all, only positive and only inversely associated transcripts. P-values of the enrichment analysis have been corrected by FDR and are indicated by color. Bullet sizes represent the gene ratio of transcripts from the associated genes. Numbers in brackets are transcripts that were mapped to the listed pathways.
The identified response on enriched pathways was not specific to AD, but included AD relevant pathways e.g., inflammatory response and endocytosis. To obtain more specific information on the cellular perturbations as a result of Cd exposure we mapped all Cd-associated transcripts onto the KEGG AD pathway (Figure 4) and identified in total 14 (4↑, 10↓) transcripts to be affected. Hence, we have identified these 14 transcripts to be potentially involved in the molecular modes of action by which Cd exposure influences AD risk in males. We assessed the overlap of the Cd associated transcripts with the reported 470 genes of a recent genetic meta-analysis of diagnosed AD patients (43). We identified 101 genes overlapping with our list of identified genes. A hypergeometric test returned a significant overlap (p = 0.02, Supplementary Figure 2A). Furthermore, we also screened for overlap with a recent single cell transcriptomic analysis of AD (44) and identified 100 genes overlapping with our list and 756 genes found in no-pathology vs. pathology in excitatory neurons. The hypergeometric test, however, showed that this was not a significant overlap (p = 0.99, Supplementary Figure 2B). Additionally, we explored the Genotype-Tissue Expression (GTEx) database to see if the identified genes are also expressed in blood. We tested the top 50 Cd associated genes ranked by significance which were also present in the GTEx database and revealed that the expression of these genes in whole blood is largley reassembled in brain and CNS tissue (Supplementary Figure 3).
Figure 4. KEGG AD pathway including Cd-associated transcripts in males. The figure shows the LMM results for transcripts associated to Cd exposure in males after mapping onto the AD pathway from KEGG. Red color indicates significant positive and green significant inverse associations with Cd exposure.
Additionally, we constructed a protein-protein interaction (PPI) network of the 2,862 Cd associated genes in males. The network was based on PPIs derived from Reactome FIs using Cytoscape. We clustered the data by using the integrated MCL Cluster function to construct networks with a minimum of 40 connected nodes. This resulted in 8 subnetworks (Supplementary Figures 4A–H). We colored the nodes by the previously calculated GeneScore derived from the multiplication of the β-coefficient and –log10(p-value). Consequently, red shades indicates a significant positive, a blue shade a significant negative association of the respective transcripts. The networks show that perturbations are centered around multiple hub genes related to cancer such as RAC1, MAPK1, TP53, UBA52 (Supplementary Figures 4A–E), as well as inflammation incl. TYK2, JAK3, and STAT3 (Supplementary Figure 4D).
The present study investigates the effect of metal exposure on molecular mechanisms involved in the etiology of AD by using the transcriptomic response of human PBLs as a surrogate model of the CNS. The whole transcriptome analysis of human subjects exposed to Pb and Cd revealed a pollutant- and gender-specific response.
It remains presently unknown whether single metals or rather, mixtures of environmental factors contribute to AD onset. We aimed to model synergistic effects of Pb and Cd by using a cumulative ZScore; however, since very few associations were found for Pb, and many associations were identified in the ZScore only, we consider these results suggestive and consequently, we have focused on the effect of individual Pb and Cd exposure.
Pb exposure was only associated with 4 (1↑, 3↓) transcripts (eIF5, PIM1, PTPN23, and SLC20A1), with the downregulation of the human eukaryotic translation initiation factor 5 (eIF5) being by far the most significant observation (Figure 1B). During initiation of cytoplasmic translation, eIF5 forms together with other proteins the 43S preinitiation complex. Translation initiation and elongation plays an important role in learning and memory, neuronal perturbations of translational control have been linked to numerous cognitive disorders (45, 46). Subtle disruptions in translational regulation can have dramatic consequences for survival of neurons, especially in the aging brain (47). Hypothetically, inappropriate reductions of translation initiation, e.g., through a suppression of eIF5 by chronic Pb exposure, might contribute to the pathogenesis of AD. When stratifying the data for gender we identified a slightly stronger response to Pb in females 11 (3↑, 8↓) (Figure 2C), thereby losing the association to eIF5. For males, however, only the inverse association of eIF5 remained significant (Figure 2D).
The LMM identified 792 (606↑, 186↓) transcripts to be associated with internal Cd exposure. Among the multiple neuroinflammatory genes with increased expression (e.g., PTEN, MAPK14, CAT, IL8RB, SASH1, TXN), the most significant association was with increased P2Y6 expression (Figure 1C). Interestingly, accumulating evidence reveals a prominent role for P2Y receptors in AD pathology. Subtle changes in P2Y receptor expression have been shown to influence molecular pathways relevant for AD pathology, including Aβ production and elimination, neuroinflammation, neuronal function, and cerebral blood flow [reviewed in (48, 49)]. A recent study has identified microglia-activated neuroinflammation through increased levels of P2Y6 in Parkinson disease patients when compared to controls (50).
For females we only observed the previously discussed association of Cd and P2Y6 (Figure 2E). Remarkably, for males we identified a very strong response to Cd exposure (Figure 2F) consisting of 2,862 (1,662↑, 1,200↓) transcripts. This is the first study showing a gender-specific effect of Cd exposure on relevant AD genes in human. In a very recent study AD mouse models, that expressed either the E4 or E3 alleles of the human APOE gene, were exposed to low-doses of Cd. Particularly male mice were vulnerable to the observed accelerated cognitive impairment induced by Cd (51). In line with the suggested gene-environment interaction between Cd and the APOE gene (51), our study identified an inverse association between Cd and APOE transcript (β-coefficient = −0.22; p = 0.0050; FDR = 0.1191) in males only. Multiple studies have reported neuronal degeneration and increased AD risk with increased expression of APOE transcript (52). In transgenic models, mouse APOE expression appears to increase Aβ deposition (53); however, other studies have also shown that human APOE expression decreases brain Aβ load (54, 55). A study on 114 patients with early and sporadic AD reported 60% decrease in APOE expression in lymphocytes from AD cases vs. controls, and also reported an inverse correlation with Aβ. The study concluded that reduced APOE expression in human may influence risk and constitute a determinant factor in Aβ loading in AD (54). Consequently, based on our lymphocyte transcriptome analysis, we suggest that Cd-dependent, chronic suppression of APOE increases Aβ deposition, suppresses Aβ clearance and/or perturbs other downstream molecular pathways, thus increasing AD risk. Furthermore, we observed an upregulation of the membrane metallo-endopeptidase (MME), an enzyme involved in Aβ metabolism (56), as well as FADD and CASP8, both suspected to mediate Aβ induced neuronal apoptosis (57). Moreover, we identified a perturbed expression of genes associated to decreased energy production and mitochondrial dysfunction with an upregulation of GAPDH and SDHA, and a suppression of COX7A1 (58), upregulation of transcripts related to neuronal cell death such as MAPK1 and GSK3B (59), and multiple deregulated genes involved in the calcium signaling (60) (Figure 4, Table 4).
We subjected the list of Cd-associated transcripts derived from males to a Reactome and KEGG overrepresentation analysis and identified multiple pathways associated to neuroinflammation (Figure 3). Thereby we identified an upregulation of neutrophil degranulation, the process by which the neutrophils induce inflammation. A large body of publications report on how activated neutrophils in the brain vessels play an important role in the pathogenesis of BBB damage, and subsequent injury of surrounding neurons and AD etiology [reviewed in (61)]. Furthermore, several other inflammatory pathways linked to AD etiology were increased in our data such as hemostasis (62), signaling by interleukins (63), and chemokines (62, 64) as well as TLR signaling (65).
We also identified an overrepresentation of upregulated genes involved in neurotrophin signaling and endocytosis. The signaling endosome hypothesis and the role of neurotrophic factors and endocytosis in AD pathogenesis have been reviewed (66). The upregulation of endocytosis is thought to increase the chances of proteins from the cell surface being packaged into endosomes. As a consequence of this process, increasing amounts of APP are brought into contact with the secretases of the endosomes and increase Aβ production (67). The importance of endocytosis in AD is further supported by the fact that enlarged endosomes are linked to age, and have also been reported in postmortem AD patients (68, 69).
Owing to the design of the EnviroGenoMarkers project the investigated population is below the age to develop AD. Furthermore, no follow up of the participants with respect to AD diagnosis has been made, and subjects were not enrolled in memory performance tests. Therefore, our study can be solely considered an exploratory effort with the aim to identify novel hypothetical association's between heavy metal exposure and molecular pathways which are potentially involved in the etiology of AD. Certainly, these findings have to be investigated in a more dedicated population.
The difference in Pb levels between both cohort and genders is most likely a result from different diets (70). In our statistical model we adjusted for the confounder cohort. All associations were only found when combining both cohorts. Certainly, we cannot exclude that there was any residual confounding; however, the stratified analyses per cohort did not return any significant associations.
The higher Cd levels in females could be a reason for a slightly higher ratio of smokers in the females (males 18% vs. females 25%). Even with the adjustment of smoking in the LMM of Cd exposure, we could not exclude any residual confounding of smoking. Therefore, we have performed a stratified analysis to never-smokers, non-smokers, and smokers. In smokers we observed only 12 positively associated AgIDs, which were also associated with Cd exposure (FDR at 20%). We also assessed the overlap of Cd associations with associations found in a previous study on the EnviroGenomarkers project assessing the influence of smoking on transcriptomics (41). The study identified 351 smoking associated AgIDs, of which only 7 overlapped with the identified Cd associated AgIDs in males. Together, these outcomes indicate that the residual confounding of smoking on our Cd analysis is limited.
To summarize, the data suggest that Pb exposure affects the transcription initiation possibly by neuronal perturbation of the 43S preinitiation complex through downregulation of eIF5, while Cd exposure, particularly in males, induces APOE suppression, neuroinflammation, and increased endocytosis. These findings possibly contribute to the understanding of the neuronal molecular response of chronic Pb and Cd exposure as well as improve our knowledge on the molecular mechanisms linking metal exposure and AD etiology.
Data Availability Statement
The datasets generated and analyzed during the current study are not publicly available due to restrictions imposed by Swedish legislation on the protection of personal data. Requests to access the datasets should be directed to Ingvar A. Bergdahl, firstname.lastname@example.org.
The studies involving human participants were reviewed and approved by Ethical Review Board in Umea and the Florence Health Unit Local Ethical Committee. The patients/participants provided their written informed consent to participate in this study.
JK, TdK, SK, and JCK designed the research. IB, AJ, and DP organized the epidemiologic part of the work. TL performed the experiments, and JK analyzed the data. TdK and JCK supervised the study. JK, TdK, and JCK co-wrote and all authors commented on the paper.
This work was supported by the EnviroGenomarkers project (Grant No. 226756) and Swedish Research Council (VR 2017-00650).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
NSHDS investigators thank the Biobank Research Unit at Umeå University, Västerbotten Intervention Programme, and the County Council of Västerbotten for providing data and samples and acknowledge the contribution from Biobank Sweden.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2020.557587/full#supplementary-material
3. Nordengen K, Kirsebom BE, Henjum K, Selnes P, Gísladóttir B, Wettergreen M, et al. Glial activation and inflammation along the Alzheimer's disease continuum. J Neuroinflammation. (2019) 16:46. doi: 10.1186/s12974-019-1399-2
6. Gatz M, Reynolds CA, Fratiglioni L, Johansson B, Mortimer JA, Berg S, et al. Role of genes and environments for explaining Alzheimer disease. Arch Gen Psychiatry. (2006) 63:168–74. doi: 10.1001/archpsyc.63.2.168
7. Yegambaram M, Manivannan B, Beach TG, Halden RU. Role of environmental contaminants in the etiology of Alzheimer's disease: a review. Curr Alzheimer Res. (2015) 12:116–46. doi: 10.2174/1567205012666150204121719
8. Tchounwou PB, Yedjou CG, Patlolla AK, Sutton DJ. Molecular, clinical and environmental toxicicology volume 3: environmental toxicology. Mol Clin Environ Toxicol. (2012) 101:133–64. doi: 10.1007/978-3-7643-8340-4_6
11. Viaene MK, Masschelein R, Leenders J, Swerts LJVC, de Groof M, Roels HA. Neurobehavioural effects of occupational exposure to cadmium: a cross sectional epidemiological study. Occup Environ Med. (2000) 57:19–27. doi: 10.1136/oem.57.1.19
12. Huat TJ, Camats-Perna J, Newcombe EA, Valmas N, Kitazawa M, Medeiros R. Metal toxicity links to Alzheimer's disease and neuroinflammation. J Mol Biol. (2019) 431:1843–68. doi: 10.1016/j.jmb.2019.01.018
14. Senut MC, Sen A, Cingolani P, Shaik A, Land SJ, Ruden DM. Lead exposure disrupts global DNA methylation in human embryonic stem cells and alters their neuronal differentiation. Toxicol Sci. (2014) 139:142–61. doi: 10.1093/toxsci/kfu028
16. Gu H, Robison G, Hong L, Barrea R, Wei X, Farlow MR, et al. Increased β-amyloid deposition in Tg-SWDI transgenic mouse brain following in vivo lead exposure. Toxicol Lett. (2012) 213:211–9. doi: 10.1016/j.toxlet.2012.07.002
17. Huang H, Bihaqi SW, Cui L, Zawia NH. In vitro Pb exposure disturbs the balance between Aβ production and elimination: the role of AβPP and neprilysin. Neurotoxicology. (2011) 32:300–6. doi: 10.1016/j.neuro.2011.02.001
18. Basha MR, Wei W, Bakheet SA, Benitez N, Siddiqi HK, Ge YW, et al. The fetal basis of amyloidogenesis: exposure to lead and latent overexpression of amyloid precursor protein and β-amyloid in the aging brain. J Neurosci. (2005) 25:823–9. doi: 10.1523/JNEUROSCI.4335-04.2005
20. Mazumdar M, Xia W, Hofmann O, Gregas M, Sui SH, Hide W, et al. Prenatal lead levels, plasma amyloid β levels, and gene expression in young adulthood. Environ Health Perspect. (2012) 120:702–7. doi: 10.1289/ehp.1104474
21. Shukla A, Shukla GS, Srimal RC. Cadmium-induced alterations in blood - Brain barrier permeability and its possible correlation with decreased microvessel antioxidant potential in rat. Hum Exp Toxicol. (1996) 15:400–5. doi: 10.1177/096032719601500507
23. Borisova T, Krisanova N, Sivko R, Kasatkina L, Borysov A, Griffin S, et al. Presynaptic malfunction: the neurotoxic effects of cadmium and lead on the proton gradient of synaptic vesicles and glutamate transport. Neurochem Int. (2011) 59:272–9. doi: 10.1016/j.neuint.2011.05.014
25. Panayi AE, Spyrou NM, Iversen BS, White MA, Part P. Determination of cadmium and zinc in Alzheimer's brain tissue using inductively coupled plasma mass spectrometry. J Neurol Sci. (2002) 195:1–10. doi: 10.1016/S0022-510X(01)00672-4
26. Jiang LF, Yao TM, Zhu ZL, Wang C, Ji LN. Impacts of Cd(II) on the conformation and self-aggregation of Alzheimer's tau fragment corresponding to the third repeat of microtubule-binding domain. Biochim Biophys Acta. (2007) 1774:1414–21. doi: 10.1016/j.bbapap.2007.08.014
27. Notarachille G, Arnesano F, Calò V, Meleleo D. Heavy metals toxicity: effect of cadmium ions on amyloid beta protein 1-42. Possible implications for Alzheimer's disease. BioMetals. (2014) 27:371–388. doi: 10.1007/s10534-014-9719-6
29. Yoo JW, Cho H, Lee KW, Won EJ, Lee YM. Combined effects of heavy metals (Cd, As, and Pb): comparative study using conceptual models and the antioxidant responses in the brackish water flea. Comp Biochem Physiol Part C Toxicol Pharmacol. (2021) 239:108863. doi: 10.1016/j.cbpc.2020.108863
30. Ciesielski T, Weuve J, Bellinger DC, Schwartz J, Lanphear B, Wright RO. Cadmium exposure and neurodevelopmental outcomes in US children. Environ Health Perspect. (2012) 120:758–63. doi: 10.1289/ehp.1104152
31. Ramsay L, Quillé M, Orset C, Grange P, Rousselet E, Férec C, et al. Blood transcriptomic biomarker as a surrogate of ischemic brain gene expression. Ann Clin Transl Neurol. (2019) 6:1681–95. doi: 10.1002/acn3.50861
32. Rollins B, Martin MV, Morgan L, Vawter MP. Analysis of whole genome biomarker expression in blood and brain. Am J Med Genet Part B Neuropsychiatr Genet. (2010) 153:919–36. doi: 10.1002/ajmg.b.31062
33. Hallmans G, Agren A, Johansson G, Johansson A, Stegmayr B, Jansson JH, et al. Cardiovascular disease and diabetes in the Northern Sweden health and disease study cohort- evaluation of risk factors and their interactions. Scand J Public Health. (2003) 31:18–24. doi: 10.1080/14034950310001432
35. Palli D, Berrino F, Vineis P, Tumino R, Panico S, Masala G, et al. A molecular epidemiology project on diet and cancer: the epic-italy prospective study. Design and baseline characteristics of participants. Tumori J. (2003) 89:586–93. doi: 10.1177/030089160308900602
36. Kelly RS, Lundh T, Porta M, Bergdahl IA, Palli D, Johansson AS, et al. Blood erythrocyte concentrations of cadmium and lead and the risk of B-cell non-Hodgkin's lymphoma and multiple myeloma: a nested case-control study. PLoS ONE. (2013) 8:e0081892. doi: 10.1371/journal.pone.0081892
37. Gaudet MM, Deubler EL, Kelly RS, Ryan Diver W, Teras LR, Hodge JM, et al. Blood levels of cadmium and lead in relation to breast cancer risk in three prospective cohorts. Int J Cancer. (2019) 144:1010–6. doi: 10.1002/ijc.31805
38. Krauskopf J, de Kok TM, Hebels DG, Bergdahl IA, Johansson A, Spaeth F, et al. MicroRNA profile for health risk assessment: environmental exposure to persistent organic pollutants strongly affects the human blood microRNA machinery. Sci Rep. (2017) 7:9262. doi: 10.1038/s41598-017-10167-7
39. Govarts E, Remy S, Bruckers L, Den Hond E, Sioen I, Nelen V, et al. Combined effects of prenatal exposures to environmental chemicals on birth weight. Int J Environ Res Public Health. (2016) 13:495. doi: 10.3390/ijerph13050495
40. Espín-Pérez A, de Kok TMCM, Jennen DGJ, Hendrickx DM, de Coster S, Schoeters G, et al. Distinct genotype-dependent differences in transcriptome responses in humans exposed to environmental carcinogens. Carcinogenesis. (2015) 36:1154–61. doi: 10.1093/carcin/bgv111
41. Georgiadis P, Hebels DG, Valavanis I, Liampa I, Bergdahl IA, Johansson A, et al. Omics for prediction of environmental health effects: blood leukocyte-based cross-omic profiling reliably predicts diseases associated with tobacco smoking. Sci Rep. (2016) 6:20544. doi: 10.1038/srep20544
42. Hebels DGAJ, Georgiadis P, Keun HC, Athersuch TJ, Vineis P, Vermeulen R, et al. Performance in omics analyses of blood samples in long-term storage: opportunities for the exploitation of existing biobanks in environmental health research. Environ Health Perspect. (2013) 121:480–7. doi: 10.1289/ehp.1205657
43. Kunkle BW, Grenier-Boley B, Sims R, Bis JC, Damotte V, Naj AC, et al. Genetic meta-analysis of diagnosed Alzheimer's disease identifies new risk loci and implicates Aβ, tau, immunity and lipid processing. Nat Genet. (2019) 51:414–30. doi: 10.1038/s41588-019-0358-2
49. Ajit D, Woods LT, Camden JM, Thebeau CN, El-Sayed FG, Greeson GW, et al. Loss of P2Y2 nucleotide receptors enhances early pathology in the TgCRND8 mouse model of Alzheimer's disease. Mol Neurobiol. (2014) 49:1031–42. doi: 10.1007/s12035-013-8577-5
50. Yang X, Lou Y, Liu G, Wang X, Qian Y, Ding J, et al. Microglia P2Y6 receptor is related to Parkinson's disease through neuroinflammatory process. J Neuroinflammation. (2017) 14:38. doi: 10.1186/s12974-017-0795-8
51. Zhang L, Wang H, Abel GM, Storm DR, Xia Z. The effects of gene-environment interactions between cadmium exposure and apolipoprotein E4 on memory in a mouse model of Alzheimer's disease. Toxicol Sci. (2020) 173:189–201. doi: 10.1093/toxsci/kfz218
52. Haasdijk ED, Vlug A, Mulder MT, Jaarsma D. Increased apolipoprotein E expression correlates with the onset of neuronal degeneration in the spinal cord of G93A-SOD1 mice. Neurosci Lett. (2002) 335:29–33. doi: 10.1016/S0304-3940(02)01159-X
54. Lambert JC, Mann D, Richard F, Tian J, Shi J, Thaker U, et al. Is there a relation between APOE expression and brain amyloid load in Alzheimer's disease? J Neurol Neurosurg Psychiatry. (2005) 76:928–33. doi: 10.1136/jnnp.2004.048983
55. Holtzman DM, Bales KR, Wu S, Bhat P, Parsadanian M, Fagan AM, et al. Expression of human apolipoprotein E reduces amyloid-β deposition in a mouse model of Alzheimer's disease. J Clin Invest. (1999) 103:R15–R21. doi: 10.1172/JCI6179
56. Tikhonova MA, Amstislavskaya TG, Belichenko VM, Fedoseeva LA, Kovalenko SP, Pisareva EE, et al. Modulation of the expression of genes related to the system of amyloid-beta metabolism in the brain as a novel mechanism of ceftriaxone neuroprotective properties. BMC Neurosci. (2018) 19:15–28. doi: 10.1186/s12868-018-0412-5
58. Liang WS, Reiman EM, Valla J, Dunckley T, Beach TG, Grover A, et al. Alzheimer's disease is associated with reduced expression of energy metabolism genes in posterior cingulate neurons. Proc Natl Acad Sci USA. (2008) 105:4441–6. doi: 10.1073/pnas.0709259105
59. del Pino J, Zeballos G, Anadón MJ, Moyano P, Díaz MJ, García JM, et al. Cadmium-induced cell death of basal forebrain cholinergic neurons mediated by muscarinic M1 receptor blockade, increase in GSK-3β enzyme, β-amyloid and tau protein levels. Arch Toxicol. (2016) 90:1081–92. doi: 10.1007/s00204-015-1540-7
66. Chen XQ, Mobley WC. Exploring the pathogenesis of Alzheimer disease in basal forebrain cholinergic neurons:Converging insights from alternative hypotheses. Front Neurosci. (2019) 13:446. doi: 10.3389/fnins.2019.00446
67. Ubelmann F, Burrinha T, Salavessa L, Gomes R, Ferreira C, Moreno N, et al. Bin1 and CD 2 AP polarise the endocytic generation of beta-amyloid. EMBO Rep. (2017) 18:102–22. doi: 10.15252/embr.201642738
68. Cataldo AM, Barnett JL, Pieroni C, Nixon RA. Increased neuronal endocytosis and protease delivery to early endosomes in sporadic Alzheimer's disease: neuropathologic evidence for a mechanism of increased β-amyloidogenesis. J Neurosci. (1997) 17:6142–51. doi: 10.1523/JNEUROSCI.17-16-06142.1997
69. Cataldo AM, Peterhoff CM, Troncoso JC, Gomez-Isla T, Hyman BT, Nixon RA. Endocytic pathway abnormalities precede amyloid β deposition in sporadic alzheimer's disease and down syndrome: differential effects of APOE genotype and presenilin mutations. Am J Pathol. (2000) 157:277–86. doi: 10.1016/S0002-9440(10)64538-5
Keywords: Alzheimer's disease, metals, gene expression, transcriptomics, microarray, APOE
Citation: Krauskopf J, Bergdahl IA, Johansson A, Palli D, Lundh T, Kyrtopoulos SA, de Kok TM and Kleinjans JC (2020) Blood Transcriptome Response to Environmental Metal Exposure Reveals Potential Biological Processes Related to Alzheimer's Disease. Front. Public Health 8:557587. doi: 10.3389/fpubh.2020.557587
Received: 30 April 2020; Accepted: 16 September 2020;
Published: 21 October 2020.
Edited by:Sonia Dagnino, Imperial College London, United Kingdom
Reviewed by:Periasamy Viayapuri Subbarayan, King Saud University, Saudi Arabia
Mariona Bustamante, Instituto Salud Global Barcelona (ISGlobal), Spain
Joseph H. Lee, Columbia University, United States
Copyright © 2020 Krauskopf, Bergdahl, Johansson, Palli, Lundh, Kyrtopoulos, de Kok and Kleinjans. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Julian Krauskopf, email@example.com