Transcriptomic Analysis in Human 3D Skin Model Injected with Resorbable Hyaluronic Acid Fillers Reveals Foreign Body Response.

Usage of injectable dermal fillers applied for aesthetic purposes has extensively increased over the years. As such, the number of related adverse reactions has increased, including patients showing severe complications such as product migration, topical swelling and inflammatory reactions of the skin. In order to understand the underlying molecular events of these adverse reactions we performed a genome-wide gene expression study on the multi-cell type human Phenion® Full-Thickness Skin Model exposed to five experimental hyaluronic acid (HA) preparations with increasing cross-linking degree, four commercial fillers from Perfectha®, and non-resorbable filler Bio-Alcamid®. In addition, we evaluated whether cross-linking degree or particle size of the HA-based fillers could be associated with the occurrence of adverse effects. In all cases, exposure to different HA fillers resulted in a clearly elevated gene expression of cytokines and chemokines related to acute inflammation as part of the foreign body response. Furthermore, for one experimental filler genes of OXPHOS complexes I-V were significantly down-regulated (adjusted p-value < 0.05), resulting in mitochondrial dysfunction which can be linked to over-expression of pro-inflammatory cytokines TNFα and IL-1ß and chemokine CCL2. Our hypothesis that cross-linking degree or particle size of the HA-based fillers is related to the biological responses induced by these fillers could only partially be confirmed for particle size. In conclusion, our innovative approach resulted in gene expression changes from a human 3D skin model exposed to dermal fillers that mechanistically substantiate aforementioned adverse reactions, and thereby adds to the weight of evidence that these fillers may induce inflammatory and fibrotic responses.

Gene expression responses reflecting 5-FU-induced toxicity: Comparison between patient colon tissue and 3D human colon organoids.

Capecitabine is a chemotherapeutic drug that is widely used as a monotherapy option in advanced cancer patients. After administration, it is converted into its active metabolite 5-fluorouracil (5-FU), a cytotoxic compound that may also induce adverse side effects in the gastrointestinal (GI) tract. Although these side effects can interfere with the continuation of the chemotherapy, diagnostic tools to detect early onset and prevention strategies are not available. In this explorative case study, we aim to identify differentially expressed genes (DEGs) that provide insight into the molecular mechanisms of toxicity induced by 5-FU in healthy colon tissue of breast cancer patients receiving capecitabine. Gene expression responses observed in patients were compared with those established in an in vitro model of healthy colon organoids. Colon biopsies from two patients with advanced breast cancer were collected before and after the treatment with capecitabine and used for RNA sequencing to determine transcriptomic responses. Differential expression analysis resulted in 31 affected genes, showing that the most affected pathways were transport of small molecules, cellular responses to stress, folate metabolism, NF-kB signalling pathway and immune system responses. The most biologically relevant genes were haemoglobin subunits encoding genes, involved in several processes; ATP12A, SLC26A3 and AQP8, involved in the transport of ions and water; TRIM31, a regulator of NF-kB signalling pathway; MST1P2 and MST1L, stimulators of macrophages. Comparison of human in vitro and in vivo responses showed that the gene expression of TRIM31 was similarly altered in the colon organoids exposed to 5-FU. Therefore, this gene constitutes a potential biomarker of colon toxicity that might be used in future in vitro drug safety design and screening.

Evaluation of Adverse Effects of Resorbable Hyaluronic Acid Fillers: Determination of Macrophage Responses.

Resorbable tissue fillers for aesthetic purposes can induce severe complications including product migration, late swelling, and inflammatory reactions. The relation between product characteristics and adverse effects is not well understood. We hypothesized that the degree of cross-linking hyaluronic acid (HA) fillers was associated with the occurrence of adverse effects. Five experimental HA preparations similar to HA fillers were synthesized with an increasing degree of cross-linking. Furthermore, a series of commercial fillers (Perfectha®) was obtained that differ in degradation time based on the size of their particulate HA components. Cytotoxic responses and cytokine production by human THP-1-derived macrophages exposed to extracts of the evaluated resorbable HA fillers were absent to minimal. Gene expression analysis of the HA-exposed macrophages revealed the responses related to cell cycle control and immune reactivity. Our results could not confirm the hypothesis that the level of cross-linking in our experimental HA fillers or the particulate size of commercial HA fillers is related to the induced biological responses. However, the evaluation of cytokine induction and gene expression in macrophages after biomaterial exposure presents promising opportunities for the development of methods to identify cellular processes that may be predictive for biomaterial-induced responses in patients.

Unravelling Mechanisms of Doxorubicin-Induced Toxicity in 3D Human Intestinal Organoids.

Doxorubicin is widely used in the treatment of different cancers, and its side effects can be severe in many tissues, including the intestines. Symptoms such as diarrhoea and abdominal pain caused by intestinal inflammation lead to the interruption of chemotherapy. Nevertheless, the molecular mechanisms associated with doxorubicin intestinal toxicity have been poorly explored. This study aims to investigate such mechanisms by exposing 3D small intestine and colon organoids to doxorubicin and to evaluate transcriptomic responses in relation to viability and apoptosis as physiological endpoints. The in vitro concentrations and dosing regimens of doxorubicin were selected based on physiologically based pharmacokinetic model simulations of treatment regimens recommended for cancer patients. Cytotoxicity and cell morphology were evaluated as well as gene expression and biological pathways affected by doxorubicin. In both types of organoids, cell cycle, the p53 signalling pathway, and oxidative stress were the most affected pathways. However, significant differences between colon and SI organoids were evident, particularly in essential metabolic pathways. Short time-series expression miner was used to further explore temporal changes in gene profiles, which identified distinct tissue responses. Finally, in silico proteomics revealed important proteins involved in doxorubicin metabolism and cellular processes that were in line with the transcriptomic responses, including cell cycle and senescence, transport of molecules, and mitochondria impairment. This study provides new insight into doxorubicin-induced effects on the gene expression levels in the intestines. Currently, we are exploring the potential use of these data in establishing quantitative systems toxicology models for the prediction of drug-induced gastrointestinal toxicity.

A Transcriptomic Approach to Elucidate the Mechanisms of Gefitinib-Induced Toxicity in Healthy Human Intestinal Organoids.

Gefitinib is a tyrosine kinase inhibitor (TKI) that selectively inhibits the epidermal growth factor receptor (EGFR), hampering cell growth and proliferation. Due to its action, gefitinib has been used in the treatment of cancers that present abnormally increased expression of EGFR. However, side effects from gefitinib therapy may occur, among which diarrhoea is most common, that can lead to interruption of the planned therapy in the more severe cases. The mechanisms underlying intestinal toxicity induced by gefitinib are not well understood. Therefore, this study aims at providing insight into these mechanisms based on transcriptomic responses induced in vitro. A 3D culture of healthy human colon and small intestine (SI) organoids was exposed to 0.1, 1, 10 and 30 µM of gefitinib, for a maximum of three days. These drug concentrations were selected using physiologically-based pharmacokinetic simulation considering patient dosing regimens. Samples were used for the analysis of viability and caspase 3/7 activation, image-based analysis of structural changes, as well as RNA isolation and sequencing via high-throughput techniques. Differential gene expression analysis showed that gefitinib perturbed signal transduction pathways, apoptosis, cell cycle, FOXO-mediated transcription, p53 signalling pathway, and metabolic pathways. Remarkably, opposite expression patterns of genes associated with metabolism of lipids and cholesterol biosynthesis were observed in colon versus SI organoids in response to gefitinib. These differences in the organoids' responses could be linked to increased activated protein kinase (AMPK) activity in colon, which can influence the sensitivity of the colon to the drug. Therefore, this study sheds light on how gefitinib induces toxicity in intestinal organoids and provides an avenue towards the development of a potential tool for drug screening and development.

Integrating SNPs-based genetic risk factor with blood epigenomic response of differentially arsenic-exposed rural subjects reveals disease-associated signaling pathways.

Arsenic (As) contamination in groundwater is responsible for numerous adverse health outcomes among millions of people. Epigenetic alterations are among the most widely studied mechanisms of As toxicity. To understand how As exposure alters gene expression through epigenetic modifications, a systematic genome-wide study was designed to address the impact of multiple important single nucleotide polymorphisms (SNPs) related to As exposure on the methylome of drinking water As-exposed rural subjects from Pakistan. Urinary As levels were used to stratify subjects into low, medium and high exposure groups. Genome-wide DNA methylation was investigated using MeDIP in combination with NimbleGen 2.1 M Deluxe Promotor arrays. Transcriptome levels were measured using Agilent 8 × 60 K expression arrays. Genotyping of selected SNPs (As3MT, DNMT1a, ERCC2, EGFR and MTHFR) was measured and an integrated genetic risk factor for each respondent was calculated by assigning a specific value to the measured genotypes based on known risk allele numbers. To select a representative model related to As exposure we compared 9 linear mixed models comprising of model 1 (including the genetic risk factor), model 2 (without the genetic risk factor) and models with individual SNPs incorporated into the methylome data. Pathway analysis was performed using ConsensusPathDB. Model 1 comprising the integrated genetic risk factor disclosed biochemical pathways including muscle contraction, cardio-vascular diseases, ATR signaling, GPCR signaling, methionine metabolism and chromatin modification in association with hypo- and hyper-methylated gene targets. A unique pathway (direct P53 effector) was found associated with the individual DNMT1a polymorphism due to hyper-methylation of CSE1L and TRRAP. Most importantly, we provide here the first evidence of As-associated DNA methylation in relation with gene expression of ATR, ATF7IP, TPM3, UBE2J2. We report the first evidence that integrating SNPs data with methylome data generates a more representative epigenome profile and discloses a better insight in disease risks of As-exposed individuals.

New insights into the mechanisms underlying 5-fluorouracil-induced intestinal toxicity based on transcriptomic and metabolomic responses in human intestinal organoids.

5-Fluorouracil (5-FU) is a widely used chemotherapeutical that induces acute toxicity in the small and large intestine of patients. Symptoms can be severe and lead to the interruption of cancer treatments. However, there is limited understanding of the molecular mechanisms underlying 5-FU-induced intestinal toxicity. In this study, well-established 3D organoid models of human colon and small intestine (SI) were used to characterize 5-FU transcriptomic and metabolomic responses. Clinically relevant 5-FU concentrations for in vitro testing in organoids were established using physiologically based pharmacokinetic simulation of dosing regimens recommended for cancer patients, resulting in exposures to 10, 100 and 1000 µM. After treatment, different measurements were performed: cell viability and apoptosis; image analysis of cell morphological changes; RNA sequencing; and metabolome analysis of supernatant from organoids cultures. Based on analysis of the differentially expressed genes, the most prominent molecular pathways affected by 5-FU included cell cycle, p53 signalling, mitochondrial ATP synthesis and apoptosis. Short time-series expression miner demonstrated tissue-specific mechanisms affected by 5-FU, namely biosynthesis and transport of small molecules, and mRNA translation for colon; cell signalling mediated by Rho GTPases and fork-head box transcription factors for SI. Metabolomic analysis showed that in addition to the effects on TCA cycle and oxidative stress in both organoids, tissue-specific metabolic alterations were also induced by 5-FU. Multi-omics integration identified transcription factor E2F1, a regulator of cell cycle and apoptosis, as the best key node across all samples. These results provide new insights into 5-FU toxicity mechanisms and underline the relevance of human organoid models in the safety assessment in drug development.

Small RNA Sequencing of Aqueous Humor and Plasma in Patients With Primary Open-Angle Glaucoma.

Purpose: Identify differentially expressed microRNAs (miRNAs) in aqueous humor (AH) and blood of primary open-angle glaucoma (POAG) patients by using small RNA sequencing. These may provide insight into POAG pathophysiology or serve as diagnostic biomarker. Methods: AH and plasma of nine POAG patients and 10 cataract control patients were small RNA sequenced on Illumina NovaSeq 6000. Identification of gene transcripts targeted by differentially expressed miRNAs was done with miRWalk and MirPath. These targets were used for pathway analysis and Gene Ontology enrichment. Diagnostic potential was evaluated by receiver operating characteristics analysis. Results: We identified 715 miRNAs in plasma and 62 miRNAs in AH. Plasma miRNA profile did not differ between POAG and control. In contrast, in AH, seven miRNAs were differentially expressed. Hsa-miR-30a-3p, hsa-miR-143-3p, hsa-miR-211-5p, and hsa-miR-221-3p were upregulated, whereas hsa-miR-92a-3p, hsa-miR-451a, and hsa-miR-486-5p were downregulated in POAG. Compared to previous studies, hsa-mir-143-3p, hsa-miR-211-5p, and hsa-miR-221-3p were reported previously, strengthening their involvement in POAG whereas hsa-miR-30a-3p, hsa-miR-92a-3p, and hsa-miR-486-5p are implicated in POAG for the first time. Identified gene transcripts were involved in several pathways, some implicated in glaucoma before (e.g., TGF-ß and neurotrophin signaling), whereas others are new (e.g., prolactin and apelin signaling). In respect to diagnostics, AH concentration of hsa-mir-143-3p had an area under the curve (AUC) of 0.889. Combined with hsa-miR-221-3p, AUC improved to 0.96. Conclusions: Small RNA sequencing identified seven differentially expressed miRNAs in AH of POAG patients. The differentially expressed miRNAs may be useful as POAG biomarkers or could become targets for new therapeutic strategies.

Generation and characterization of iPSC-derived renal proximal tubule-like cells with extended stability.

The renal proximal tubule is responsible for re-absorption of the majority of the glomerular filtrate and its proper function is necessary for whole-body homeostasis. Aging, certain diseases and chemical-induced toxicity are factors that contribute to proximal tubule injury and chronic kidney disease progression. To better understand these processes, it would be advantageous to generate renal tissues from human induced pluripotent stem cells (iPSC). Here, we report the differentiation and characterization of iPSC lines into proximal tubular-like cells (PTL). The protocol is a step wise exposure of small molecules and growth factors, including the GSK3 inhibitor (CHIR99021), the retinoic acid receptor activator (TTNPB), FGF9 and EGF, to drive iPSC to PTL via cell stages representing characteristics of early stages of renal development. Genome-wide RNA sequencing showed that PTL clustered within a kidney phenotype. PTL expressed proximal tubular-specific markers, including megalin (LRP2), showed a polarized phenotype, and were responsive to parathyroid hormone. PTL could take up albumin and exhibited ABCB1 transport activity. The phenotype was stable for up to 7 days and was maintained after passaging. This protocol will form the basis of an optimized strategy for molecular investigations using iPSC derived PTL.

microRNA expression profiles and personal monitoring of exposure to particulate matter.

An increasing number of findings from epidemiological studies support associations between exposure to air pollution and the onset of several diseases, including pulmonary, cardiovascular and neurodegenerative diseases, and malignancies. However, intermediate, and potentially mediating, biological mechanisms associated with exposure to air pollutants are largely unknown. Previous studies on the human exposome have shown that the expression of certain circulating microRNAs (miRNAs), regulators of gene expression, are altered upon exposure to traffic-related air pollutants. In the present study, we investigated the relationship between particulate matter (PM) smaller than 2.5 µm (PM2.5), PM2.5 absorbance (as a proxy of black carbon and soot), and ultrafine-particles (UFP, smaller than 0.1 µm), measured in healthy volunteers by 24 h personal monitoring (PEM) sessions and global expression levels of peripheral blood miRNAs. The PEM sessions were conducted in four European countries, namely Switzerland (Basel), United Kingdom (Norwich), Italy (Turin), and The Netherlands (Utrecht). miRNAs expression levels were analysed using microarray technology on blood samples from 143 participants. Seven miRNAs, hsa-miR-24-3p, hsa-miR-4454, hsa-miR-4763-3p, hsa-miR-425-5p, hsa-let-7d-5p, hsa-miR-502-5p, and hsa-miR-505-3p were significantly (FDR corrected) expressed in association with PM2.5 personal exposure, while no significant association was found between miRNA expression and the other pollutants. The results obtained from this investigation suggest that personal exposure to PM2.5 is associated with miRNA expression levels, showing the potential for these circulating miRNAs as novel biomarkers for air pollution health risk assessment.

Transcriptome responses in blood reveal distinct biological pathways associated with arsenic exposure through drinking water in rural settings of Punjab, Pakistan.

BACKGROUND: Groundwater Arsenic (As) contamination is a global public health concern responsible for various health implications and a neglected area of environmental health research in Pakistan. Because of interindividual differences in genetic predisposition, As-related health issues may not be equally distributed among the As-exposed population. However, till date, no studies have been conducted including multiple SNPs involved in As metabolism and disease risk using a linear mixed effect model approach to analyze peripheral blood transcriptomics results. OBJECTIVES: In order to detect early responses on the gene expression level and to evaluate the impact of selected SNPs inferring disease risks associated with As exposure, we designed a systematic study to investigate blood transcriptomics profiles of 57 differentially exposed rural subjects living in drinking water As-contaminated settings of Lahore and Kasur districts in Punjab Province in southeast Pakistan. Exposure among the subjects was correlated with individual transcriptome responses applying urinary As profiles as the main biomarker for risk stratification. METHODS: We performed whole genome gene expression analysis in blood of subjects using microarrays. Linear effect mixed models were applied for evaluating the combined impact of SNPs hypothetically increasing the risk for As exposure-induced health effects (GSTM1, GSTT1, As3MT, DNMT1, MTHFR, ERCC2 and EGFR). RESULTS: Our findings confirmed important signaling, growth factor, cancer and other disease related pathways known to be associated with increased As exposure levels. In addition, upon implementing our integrative SNPs-based genetic risk factor, pathways associated with an increased risk of NAFLD and diabetes appeared significantly enhanced by down-regulation of genes NDUFV3, IKBKB, IL6R, ADIPOR1, PPARA, OGT and FOXO1. CONCLUSION: We report the first comprehensive study applying state-of-the-art bioinformatics approaches to address multiple SNP-based inter-individual variability in adverse molecular responses among subjects exposed to drinking water As contamination in Pakistan thereby providing strong evidence of various gene expression targets associated with development of known As-related diseases.

Acetaminophen Overdose as a Potential Risk Factor for Parkinson's Disease.

Four complementary approaches were used to investigate acetaminophen overdose as a risk factor for Parkinson's disease (PD). Circulating microRNAs (miRNAs) serum profiles from acetaminophen-overdosed patients were compared with patients with terminal PD, revealing four shared miRNAs. Similarities were found among molecular structures of dopamine (DA), acetaminophen, and two known PD inducers indicating affinity for dopaminergic transport. Potential interactions between acetaminophen and the human DA transporter were confirmed by molecular docking modeling and binding free energy calculations. Thus, it is plausible that acetaminophen is taken up by the dopaminergic transport system into the substantia nigra (SN). A ChEMBL query identified proteins that are similarly targeted by DA and acetaminophen. Here, we highlight CYP3A4, present in the SN, a predominant metabolizer of acetaminophen into its toxic metabolite N-acetyl-p-benzoquinone imine and shown to be regulated in PD. Overall, based on our results, we hypothesize that overdosing of acetaminophen is a potential risk factor for parkinsonism.

Drug-induced gene expression profile changes in relation to intestinal toxicity: State-of-the-art and new approaches.

One of the major complications that patients experience during pharmacological treatment is the occurrence of adverse drug reactions (ADRs). The most affected organs are the liver, kidney, heart and the gastrointestinal-immune system. In comparison to the other organs, less progress has been made on human-relevant prediction of drug-induced intestinal toxicity, evidencing current large data gaps. The most widely used drugs that are associated with intestinal damage include chemotherapeutics, such as 5-Fluorouracil or Tyrosine Kinase Inhibitors (TKIs), as well as non-steroidal anti-inflammatory drugs (NSAIDs). Chemotherapeutics are regarded as inducers of acute intestinal toxicity whereas NSAIDs are associated with chronic inflammation of the intestine. In view of the fact that only a few studies have been dedicated to studying cellular and genomic responses in relation to drug-induced intestinal ADRs, little is known about how intestinal toxicity develops after exposure to such drugs or which molecular mechanisms are involved. Therefore, new models and experiments are required to establish transcriptomic responses and alterations of molecular markers induced by different medicines. This review summarizes the available information about transcriptomic responses and biomarkers of toxicity induced by 5-FU, NSAIDS or TKIs in different experimental models. Future investigation should address the challenges in predicting intestinal toxicity induced by drugs and unveil specific gene expression profiles that can be applied in the development of safer drugs.

DNA methylation profiling implicates exposure to PCBs in the pathogenesis of B-cell chronic lymphocytic leukemia.

OBJECTIVES: To characterize the impact of PCB exposure on DNA methylation in peripheral blood leucocytes and to evaluate the corresponding changes in relation to possible health effects, with a focus on B-cell lymphoma. METHODS: We conducted an epigenome-wide association study on 611 adults free of diagnosed disease, living in Italy and Sweden, in whom we also measured plasma concentrations of 6 PCB congeners, DDE and hexachlorobenzene. RESULTS: We identified 650 CpG sites whose methylation correlates strongly (FDR < 0.01) with plasma concentrations of at least one PCB congener. Stronger effects were observed in males and in Sweden. This epigenetic exposure profile shows extensive and highly statistically significant overlaps with published profiles associated with the risk of future B-cell chronic lymphocytic leukemia (CLL) as well as with clinical CLL (38 and 28 CpG sites, respectively). For all these sites, the methylation changes were in the same direction for increasing exposure and for higher disease risk or clinical disease status, suggesting an etiological link between exposure and CLL. Mediation analysis reinforced the suggestion of a causal link between exposure, changes in DNA methylation and disease. Disease connectivity analysis identified multiple additional diseases associated with differentially methylated genes, including melanoma for which an etiological link with PCB exposure is established, as well as developmental and neurological diseases for which there is corresponding epidemiological evidence. Differentially methylated genes include many homeobox genes, suggesting that PCBs target stem cells. Furthermore, numerous polycomb protein target genes were hypermethylated with increasing exposure, an effect known to constitute an early marker of carcinogenesis. CONCLUSIONS: This study provides mechanistic evidence in support of a link between exposure to PCBs and the etiology of CLL and underlines the utility of omic profiling in the evaluation of the potential toxicity of environmental chemicals.

Identification of Sex-Specific Transcriptome Responses to Polychlorinated Biphenyls (PCBs).

PCBs are classified as xenoestrogens and carcinogens and their health risks may be sex-specific. To identify potential sex-specific responses to PCB-exposure we established gene expression profiles in a population study subdivided into females and males. Gene expression profiles were determined in a study population consisting of 512 subjects from the EnviroGenomarkers project, 217 subjects who developed lymphoma and 295 controls were selected in later life. We ran linear mixed models in order to find associations between gene expression and exposure to PCBs, while correcting for confounders, in particular distribution of white blood cells (WBC), as well as random effects. The analysis was subdivided according to sex and development of lymphoma in later life. The changes in gene expression as a result of exposure to the six studied PCB congeners were sex- and WBC type specific. The relatively large number of genes that are significantly associated with PCB-exposure in the female subpopulation already indicates different biological response mechanisms to PCBs between the two sexes. The interaction analysis between different PCBs and WBCs provides only a small overlap between sexes. In males, cancer-related pathways and in females immune system-related pathways are identified in association with PCBs and WBCs. Future lymphoma cases and controls for both sexes show different responses to the interaction of PCBs with WBCs, suggesting a role of the immune system in PCB-related cancer development.

Phosphorylation of eIF2α promotes cell survival in response to benzo[a]pyrene exposure.

Cellular adaptation is important to cope with various stresses induced by altered environmental conditions. By controlling mRNA translation rates cells may adapt to stress to promote survival. Phosphorylation of eIF2α at serine 51 is one of the pathways controlling mRNA translation. Here we investigated the role of phosphorylated eIF2α during exposure to the environmental carcinogen benzo(a)pyrene (BaP). For our study we used mouse embryonic fibroblasts with a wild type eIF2α (MEF WT) and mouse embryonic fibroblasts with an eIF2α S51A knock-in mutation that cannot be phosphorylated. Here, we show that eIF2α phosphorylation occurs in MEF WT cells but not in MEF S51A cells. Survival of MEF S51A cells is profoundly reduced compared to MEF WT controls after BaP exposure. No differences in DNA damage or ROS production were observed between MEF WT and S51A cells. Disruption of eIF2α phosphorylation caused increased levels of apoptosis in response to BaP. This work demonstrates that eIF2α phosphorylation is important for reducing apoptosis and promoting cell survival in order to adapt to BaP exposure.

Persistence of Epigenomic Effects After Recovery From Repeated Treatment With Two Nephrocarcinogens.

The discovery of the epigenetic regulation of transcription has provided a new source of mechanistic understanding to long lasting effects of chemicals. However, this information is still seldom exploited in a toxicological context and studies of chemical effect after washout remain rare. Here we studied the effects of two nephrocarcinogens on the human proximal tubule cell line RPTEC/TERT1 using high-content mRNA microarrays coupled with miRNA, histone acetylation (HA) and DNA methylation (DM) arrays and metabolomics during a 5-day repeat-dose exposure and 3 days after washout. The mycotoxin ochratoxin A (OTA) was chosen as a model compound for its known impact on HA and DM. The foremost effect observed was the modulation of thousands of mRNAs and histones by OTA during and after exposure. In comparison, the oxidant potassium bromate (KBrO3) had a milder impact on gene expression and epigenetics. However, there was no strong correlation between epigenetic modifications and mRNA changes with OTA while with KBrO3 the gene expression data correlated better with HA for both up- and down-regulated genes. Even when focusing on the genes with persistent epigenetic modifications after washout, only half were coupled to matching changes in gene expression induced by OTA, suggesting that while OTA causes a major effect on the two epigenetic mechanisms studied, these alone cannot explain its impact on gene expression. Mechanistic analysis confirmed the known activation of Nrf2 and p53 by KBrO3, while OTA inhibited most of the same genes, and genes involved in the unfolded protein response. A few miRNAs could be linked to these effects of OTA, albeit without clear contribution of epigenetics to the modulation of the pathways at large. Metabolomics revealed disturbances in amino acid balance, energy catabolism, nucleotide metabolism and polyamine metabolism with both chemicals. In conclusion, the large impact of OTA on transcription was confirmed at the mRNA level but also with two high-content epigenomic methodologies. Transcriptomic data confirmed the previously reported activation (by KBrO3) and inhibition (by OTA) of protective pathways. However, the integration of omic datasets suggested that HA and DM were not driving forces in the gene expression changes induced by either chemical.

Comparison of statistical methods and the use of quality control samples for batch effect correction in human transcriptome data.

Batch effects are technical sources of variation introduced by the necessity of conducting gene expression analyses on different dates due to the large number of biological samples in population-based studies. The aim of this study is to evaluate the performances of linear mixed models (LMM) and Combat in batch effect removal. We also assessed the utility of adding quality control samples in the study design as technical replicates. In order to do so, we simulated gene expression data by adding "treatment" and batch effects to a real gene expression dataset. The performances of LMM and Combat, with and without quality control samples, are assessed in terms of sensitivity and specificity while correcting for the batch effect using a wide range of effect sizes, statistical noise, sample sizes and level of balanced/unbalanced designs. The simulations showed small differences among LMM and Combat. LMM identifies stronger relationships between big effect sizes and gene expression than Combat, while Combat identifies in general more true and false positives than LMM. However, these small differences can still be relevant depending on the research goal. When any of these methods are applied, quality control samples did not reduce the batch effect, showing no added value for including them in the study design.

Short-term transcriptome and microRNAs responses to exposure to different air pollutants in two population studies.

Diesel vehicle emissions are the major source of genotoxic compounds in ambient air from urban areas. These pollutants are linked to risks of cardiovascular diseases, lung cancer, respiratory infections and adverse neurological effects. Biological events associated with exposure to some air pollutants are widely unknown but applying omics techniques may help to identify the molecular processes that link exposure to disease risk. Most data on health risks are related to long-term exposure, so the aim of this study is to investigate the impact of short-term exposure (two hours) to air pollutants on the blood transcriptome and microRNA expression levels. We analyzed transcriptomics and microRNA expression using microarray technology on blood samples from volunteers participating in studies in London, the Oxford Street cohort, and, in Barcelona, the TAPAS cohort. Personal exposure levels measurements of particulate matter (PM10, PM2.5), ultrafine particles (UFPC), nitrogen oxides (NO2, NO and NOx), black carbon (BC) and carbon oxides (CO and CO2) were registered for each volunteer. Associations between air pollutant levels and gene/microRNA expression were evaluated using multivariate normal models (MVN). MVN-models identified compound-specific expression of blood cell genes and microRNAs associated with air pollution despite the low exposure levels, the short exposure periods and the relatively small-sized cohorts. Hsa-miR-197-3p, hsa-miR-29a-3p, hsa-miR-15a-5p, hsa-miR-16-5p and hsa-miR-92a-3p are found significantly expressed in association with exposures. These microRNAs target also relevant transcripts, indicating their potential relevance in the research of omics-biomarkers responding to air pollution. Furthermore, these microRNAs are also known to be associated with diseases previously linked to air pollution exposure including several cancers such lung cancer and Alzheimer's disease. In conclusion, we identified in this study promising compound-specific mRNA and microRNA biomarkers after two hours of exposure to low levels of air pollutants during two hours that suggest increased cancer risks.