Sex-specific associations between particulate matter exposure and gene expression in independent discovery and validation cohorts of middle-aged men and women

BACKGROUND: Particulate matter (PM) exposure leads to premature death, mainly due to respiratory and cardiovascular diseases. OBJECTIVES: Identification of transcriptomic biomarkers of air pollution exposure and effect in a healthy adult population. METHODS: Microarray analyses were performed in 98 healthy volunteers (48 men, 50 women). The expression of eight sex-specific candidate biomarker genes (significantly associated with PM(10) in the discovery cohort and with a reported link to air pollution-related disease) was measured with qPCR in an independent validation cohort (75 men, 94 women). Pathway analysis was performed using Gene Set Enrichment Analysis. Average daily PM(2.5) and PM(10) exposures over 2-years were estimated for each participant’s residential address using spatiotemporal interpolation in combination with a dispersion model. RESULTS: Average long-term PM(10) was 25.9 (± 5.4) and 23.7 (± 2.3) μg/m(3) in the discovery and validation cohorts, respectively. In discovery analysis, associations between PM(10) and the expression of individual genes differed by sex. In the validation cohort, long-term PM(10) was associated with the expression of DNAJB5 and EAPP in men and ARHGAP4 (p = 0.053) in women. AKAP6 and LIMK1 were significantly associated with PM(10) in women, although associations differed in direction between the discovery and validation cohorts. Expression of the eight candidate genes in the discovery cohort differentiated between validation cohort participants with high versus low PM(10) exposure (area under the receiver operating curve = 0.92; 95% CI: 0.85, 1.00; p = 0.0002 in men, 0.86; 95% CI: 0.76, 0.96; p = 0.004 in women). CONCLUSIONS: Expression of the sex-specific candidate genes identified in the discovery population predicted PM(10) exposure in an independent cohort of adults from the same area. Confirmation in other populations may further support this as a new approach for exposure assessment, and may contribute to the discovery of molecular mechanisms for PM-induced health effects. CITATION: Vrijens K, Winckelmans E, Tsamou M, Baeyens W, De Boever P, Jennen D, de Kok TM, Den Hond E, Lefebvre W, Plusquin M, Reynders H, Schoeters G, Van Larebeke N, Vanpoucke C, Kleinjans J, Nawrot TS. 2017. 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 125:660–669; http://dx.doi.org/10.1289/EHP37

Dose and Time Dependencies in Stress Pathway Responses during Chemical Exposure: Novel Insights from Gene Regulatory Networks

Perturbation of biological networks is often observed during exposure to xenobiotics, and the identification of disturbed processes, their dynamic traits, and dose–response relationships are some of the current challenges for elucidating the mechanisms determining adverse outcomes. In this scenario, reverse engineering of gene regulatory networks (GRNs) from expression data may provide a system-level snapshot embedded within accurate molecular events. Here, we investigate the composition of GRNs inferred from groups of chemicals with two distinct outcomes, namely carcinogenicity [azathioprine (AZA) and cyclophosphamide (CYC)] and drug-induced liver injury (DILI; diclofenac, nitrofurantoin, and propylthiouracil), and a non-carcinogenic/non-DILI group (aspirin, diazepam, and omeprazole). For this, we analyzed publicly available exposed in vitro human data, taking into account dose and time dependencies. Dose–Time Network Identification (DTNI) was applied to gene sets from exposed primary human hepatocytes using four stress pathways, namely endoplasmic reticulum (ER), NF-κB, NRF2, and TP53. Inferred GRNs suggested case specificity, varying in interactions, starting nodes, and target genes across groups. DILI and carcinogenic compounds were shown to directly affect all pathway-based GRNs, while non-DILI/non-carcinogenic chemicals only affected NF-κB. NF-κB-based GRNs clearly illustrated group-specific disturbances, with the cancer-related casein kinase CSNK2A1 being a target gene only in the carcinogenic group, and opposite regulation of NF-κB subunits being observed in DILI and non-DILI/non-carcinogenic groups. Target genes in NRF2-based GRNs shared by DILI and carcinogenic compounds suggested markers of hepatotoxicity. Finally, we indicate several of these group-specific interactions as potentially novel. In summary, our reversed-engineered GRNs are capable of revealing dose dependent, chemical-specific mechanisms of action in stress-related biological networks

MagiCMicroRna: a web implementation of AgiMicroRna using shiny

BACKGROUND: MicroRNA expression can be quantified using sequencing techniques or commercial microRNA-expression arrays. Recently, the AgiMicroRna R-package was published that enabled systematic preprocessing and statistical analysis for Agilent microRNA arrays. Here we describe MagiCMicroRna, which is a user-friendly web interface for this package, together with a new filtering approach. RESULTS: We used MagiCMicroRna to normalize and filter an Agilent miRNA microarray dataset of cancerous and normal tissues from 14 different patients. With the standard filtering procedure, 250 out of 817 microRNAs remained, whereas the new group-specific filtering approach resulted in broader datasets for further analysis in most groups (>279 microRNAs remaining). CONCLUSIONS: The user-friendly web interface of MagiCMicroRna enables researchers to normalize and filter Agilent microarrays by the click of one button. Furthermore, MagiCMicroRna provides flexibility in choosing the filtering method. The new group-specific filtering approach lead to an increased number and additional tissue-specific microRNAs remaining for subsequent analysis compared to the standard procedure. The MagiCMicroRna web interface and source code can be downloaded from https://bitbucket.org/mutgx/magicmicrorna.git. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s13029-015-0035-5) contains supplementary material, which is available to authorized users

Evaluating microRNA profiles reveals discriminative responses following genotoxic or non-genotoxic carcinogen exposure in primary mouse hepatocytes

Chemical carcinogenesis can be induced by genotoxic (GTX) or non-genotoxic (NGTX) carcinogens. GTX carcinogens have a well-described mode of action. However, the complex mechanisms by which NGTX carcinogens act are less clear and may result in conflicting results between species [e.g. Wy-14,643 (Wy)]. We hypothesise that common microRNA response pathways exist for each class of carcinogenic agents. Therefore, this study compares and integrates mRNA and microRNA expression profiles following short term acute exposure (24 and 48h) to three GTX [aflatoxin B1 (AFB1), benzo[a]pyrene (BaP) and cisplatin (CisPl)] or three NGTX (2,3,7,8-tetrachloordibenzodioxine (TCDD), cyclosporine A (CsA) and Wy) carcinogens in primary mouse hepatocytes. Discriminative gene sets, microRNAs (not for 24h) and processes were identified following 24 and 48h of exposure. From the three discriminative microRNAs found following 48h of exposure, mmu-miR-503-5p revealed to have an interaction with mRNA target gene cyclin D2 (Ccnd2 - 12444) which was involved in the discriminative process of p53 signalling and metabolism. Following exposure to NGTX carcinogens Mmu-miR-503-5p may have an oncogenic function by stimulating Ccnd2 possibly leading to a tumourigenic cell cycle progression. By contrast, after GTX carcinogen exposure it may have a tumour-suppressive function (repressing Ccnd2) leading to cell cycle arrest and to increased DNA repair activities. In addition, compound-specific microRNA-mRNA interactions [mmu-miR-301b-3p-Papss2 (for AFB1), as well as mmu-miR-29b-3p-Col4a2 and mmu-miR-24-3p-Flna (for BaP)] were found to contribute to a better understanding of microRNAs in cell cycle arrest and the impairment of the DNA damage repair, an important hallmark of GTX-induced carcinogenesis. Overall, our results indicate that microRNAs represent yet another relevant intracellular regulatory level in chemical carcinogenesis

Exploiting microRNA and mRNA profiles generated in vitro from carcinogen-exposed primary mouse hepatocytes for predicting in vivo genotoxicity and carcinogenicity

The well-defined battery of in vitro systems applied within chemical cancer risk assessment is often characterised by a high false-positive rate, thus repeatedly failing to correctly predict the in vivo genotoxic and carcinogenic properties of test compounds. Toxicogenomics, i.e. mRNA-profiling, has been proven successful in improving the prediction of genotoxicity in vivo and the understanding of underlying mechanisms. Recently, microRNAs have been discovered as post-transcriptional regulators of mRNAs. It is thus hypothesised that using microRNA response-patterns may further improve current prediction methods. This study aimed at predicting genotoxicity and non-genotoxic carcinogenicity in vivo, by comparing microRNA-and mRNA-based profiles, using a frequently applied in vitro liver model and exposing this to a range of well-chosen prototypical carcinogens. Primary mouse hepatocytes (PMH) were treated for 24 and 48 h with 21 chemical compounds [genotoxins (GTX) vs. non-genotoxins (NGTX) and non-genotoxic carcinogens (NGTX-C) versus non-carcinogens (NC)]. MicroRNA and mRNA expression changes were analysed by means of Exiqon and Affymetrix microarray-platforms, respectively. Classification was performed by using Prediction Analysis for Microarrays (PAM). Compounds were randomly assigned to training and validation sets (repeated 10 times). Before prediction analysis, pre-selection of microRNAs and mRNAs was performed by using a leave-one-out t-test. No microRNAs could be identified that accurately predicted genotoxicity or non-genotoxic carcinogenicity in vivo. However, mRNAs could be detected which appeared reliable in predicting genotoxicity in vivo after 24 h (7 genes) and 48 h (2 genes) of exposure (accuracy: 90% and 93%, sensitivity: 65% and 75%, specificity: 100% and 100%). Tributylinoxide and para-Cresidine were misclassified. Also, mRNAs were identified capable of classifying NGTX-C after 24 h (5 genes) as well as after 48 h (3 genes) of treatment (accuracy: 78% and 88%, sensitivity: 83% and 83%, specificity: 75% and 93%). Wy-14,643, phenobarbital and ampicillin trihydrate were misclassified. We conclude that genotoxicity and non-genotoxic carcinogenicity probably cannot be accurately predicted based on microRNA profiles. Overall, transcript-based prediction analyses appeared to clearly outperform microRNA-based analyses